求高版本1stopt帮忙解决dingd老师帮忙啊
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
数据少了一列积分区间的数据吧，如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
页:
[1]