求高版本1stopt帮忙解决dingd老师帮忙啊

郑发 · 发表于 2013-5-18 15:09

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

您需要登录才可以下载或查看，没有账号？我要加入

x

Parameters a, b, c, d;
Variable x, y;
ODEFunction x'=a*x-b*x*y;    //Bristow & Campbell (1984)
                  y'=-c*y+d*x*y;
Data;
45.79    41.40
53.03    38.90
64.05    36.78
75.40    36.04
90.36    33.78
107.14    35.40
127.79    34.68
150.77    36.61
179.65    37.71
211.82    41.98
249.91    45.72
291.31    53.10
334.95    65.44
380.67    83.00
420.28    108.74
445.56    150.01
447.63    205.61
414.04    281.60
347.04    364.56
265.33    440.30
187.57    489.68
128.00    512.95
85.25    510.01
57.17    491.06
39.96    462.22
29.22    430.15
22.30    396.95
16.52    364.87
14.41    333.16
11.58    304.97
10.41    277.73
10.17    253.16
7.86    229.66
9.23    209.53
8.22    190.07
8.76    173.58
7.90    156.40
8.38    143.05
9.53    130.75
9.33    117.49
9.72    108.16
10.55    98.08
13.05    88.91
13.58    82.28
16.31    75.42
17.75    69.58
20.11    62.58
23.98    59.22
28.51    54.91
31.61    49.79
37.13    45.94
45.06    43.41
53.40    41.30
62.39    40.28
72.89    37.71
86.92    36.58
103.32    36.98
121.70    36.65
144.86    37.87
171.92    39.63
202.51    42.97
237.69    46.95
276.77    54.93
319.76    64.61
362.05    81.28
400.11    105.50
427.79    143.03
434.56    192.45
410.31    260.84
354.18    339.39
278.49    413.79
203.72    466.94
141.06    494.72
95.08    499.37
66.76    484.58
45.41    460.63
33.13    429.79
25.89    398.77
20.51    366.49
17.11    336.56
12.69    306.39
11.76    279.53
11.22    254.95
10.29    233.50
8.82    212.74
9.51    193.61
8.69    175.01
9.53    160.59
8.68    146.12
10.82    131.85

dingd · 发表于 2013-5-22 15:31

数据少了一列积分区间的数据吧，如t：

Parameters a, b, c, d;
Variable t, x, y;
ODEFunction x'=a*x-b*x*y; //Bristow & Campbell (1984)
y'=-c*y+d*x*y;
Data;
1 45.79 41.40
2 53.03 38.90
3 64.05 36.78
4 75.40 36.04
5 90.36 33.78
6 107.14 35.40
7 127.79 34.68
8 150.77 36.61
9 179.65 37.71
10 211.82 41.98
11 249.91 45.72
12 291.31 53.10
13 334.95 65.44
14 380.67 83.00
15 420.28 108.74
16 445.56 150.01
17 447.63 205.61
18 414.04 281.60
19 347.04 364.56
20 265.33 440.30
21 187.57 489.68
22 128.00 512.95
23 85.25 510.01
24 57.17 491.06
25 39.96 462.22
26 29.22 430.15
27 22.30 396.95
28 16.52 364.87
29 14.41 333.16
30 11.58 304.97
31 10.41 277.73
32 10.17 253.16
33 7.86 229.66
34 9.23 209.53
35 8.22 190.07
36 8.76 173.58
37 7.90 156.40
38 8.38 143.05
39 9.53 130.75
40 9.33 117.49
41 9.72 108.16
42 10.55 98.08
43 13.05 88.91
44 13.58 82.28
45 16.31 75.42
46 17.75 69.58
47 20.11 62.58
48 23.98 59.22
49 28.51 54.91
50 31.61 49.79
51 37.13 45.94
52 45.06 43.41
53 53.40 41.30
54 62.39 40.28
55 72.89 37.71
56 86.92 36.58
57 103.32 36.98
58 121.70 36.65
59 144.86 37.87
60 171.92 39.63
61 202.51 42.97
62 237.69 46.95
63 276.77 54.93
64 319.76 64.61
65 362.05 81.28
66 400.11 105.50
67 427.79 143.03
68 434.56 192.45
69 410.31 260.84
70 354.18 339.39
71 278.49 413.79
72 203.72 466.94
73 141.06 494.72
74 95.08 499.37
75 66.76 484.58
76 45.41 460.63
77 33.13 429.79
78 25.89 398.77
79 20.51 366.49
80 17.11 336.56
81 12.69 306.39
82 11.76 279.53
83 11.22 254.95
84 10.29 233.50
85 8.82 212.74
86 9.51 193.61
87 8.69 175.01
88 9.53 160.59
89 8.68 146.12
90 10.82 131.85

复制代码

均方差(RMSE): 3.01974154808004
残差平方和(SSE): 1623.15334506175
相关系数(R): 0.999906678609375
相关系数之平方(R^2): 0.999813365927633
确定系数(DC): 0.999810731430232
F统计(F-Statistic): 71454.6315934645

参数最佳估算
-------------------- -------------
a 0.213767217206914
b 0.00121118818676084
c 0.104331027261455
d 0.000962784941586087

账号		自动登录	找回密码
密码			我要加入

[1stopt] 求高版本1stopt帮忙解决dingd老师帮忙啊

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。