求助 ~~语法有错误运行不了~~帮忙看看...

bolzanow

谢谢...明天就要交论文了
搞不出一个结果还是满郁闷的

my test programm;
Parameter a(1:3),x(1:3);
Parameter b(1:3),y(1:3);
Parameter c(1:3),z(1:3);
Variable t,s,p,u;
function u=(sum(i=1:3)a*e^(x*s))* (sum(j=1:3)b[j]*e^(y[j]*t))* (sum(k=1:3)c[k]*e^(z[k]*p));
data;
// t,s,p,u
-2.0005 0 0 999.67
-1.7894 0.1005 0.50251 1000
-1.5783 0.20101 1.005 1000.4
-1.3672 0.30151 1.5075 1000.8
-1.1561 0.40201 2.0101 1001.1
-0.94495 0.50251 2.5126 1001.5
-0.73384 0.60302 3.0151 1001.8
-0.52274 0.70352 3.5176 1002.2
-0.31163 0.80402 4.0201 1002.5
-0.10053 0.90452 4.5226 1002.9
0.11058 1.005 5.0251 1003.2
0.32168 1.1055 5.5276 1003.5
0.53279 1.206 6.0302 1003.9
0.74389 1.3065 6.5327 1004.2
0.955 1.407 7.0352 1004.6
1.1661 1.5075 7.5377 1004.9
1.3772 1.608 8.0402 1005.2
1.5883 1.7085 8.5427 1005.6
1.7994 1.809 9.0452 1005.9
2.0105 1.9095 9.5477 1006.2
2.2216 2.0101 10.05 1006.5
2.4327 2.1106 10.553 1006.9
2.6438 2.2111 11.055 1007.2
2.855 2.3116 11.558 1007.5
3.0661 2.4121 12.06 1007.8
3.2772 2.5126 12.563 1008.1
3.4883 2.6131 13.065 1008.4
3.6994 2.7136 13.568 1008.7
3.9105 2.8141 14.07 1009.1
4.1216 2.9146 14.573 1009.4
4.3327 3.0151 15.075 1009.7
4.5438 3.1156 15.578 1010
4.7549 3.2161 16.08 1010.3
4.966 3.3166 16.583 1010.6
5.1771 3.4171 17.085  1010.9
5.3882 3.5176 17.588 1011.2
5.5993 3.6181 18.09 1011.5
5.8104 3.7186 18.593 1011.8
6.0215 3.8191 19.095 1012.1
6.2327 3.9196 19.598 1012.3
6.4438 4.0201 20.101 1012.6
6.6549 4.1206 20.603 1012.9
6.866 4.2211 21.106 1013.2
7.0771 4.3216 21.608 1013.5
7.2882 4.4221 22.111 1013.8
7.4993 4.5226 22.613 1014.1
7.7104 4.6231 23.116 1014.3
7.9215 4.7236 23.618 1014.6
8.1326 4.8241 24.121 1014.9
8.3437 4.9246 24.623 1015.2
8.5548 5.0251 25.126 1015.4
8.7659 5.1256 25.628 1015.7
8.977 5.2261 26.131 1016
9.1882 5.3266 26.633 1016.3
9.3993 5.4271 27.136 1016.5
9.6104 5.5276 27.638 1016.8
9.8215 5.6281 28.141 1017.1
10.033 5.7286 28.643 1017.3
10.244 5.8291 29.146 1017.6
10.455 5.9296 29.648 1017.9
10.666 6.0302 30.151 1018.1
10.877 6.1307 30.653 1018.4
11.088 6.2312 31.156 1018.6
11.299 6.3317 31.658 1018.9
11.51 6.4322 32.161 1019.2
11.721 6.5327 32.663 1019.4
11.933 6.6332 33.166 1019.7
12.144 6.7337 33.668 1019.9
12.355 6.8342 34.171 1020.2
12.566 6.9347 34.673 1020.4
12.777 7.0352 35.176 1020.7
12.988 7.1357 35.678 1020.9
13.199 7.2362 36.181 1021.2
13.41 7.3367 36.683 1021.4
13.621 7.4372 37.186 1021.7
13.833 7.5377 37.688 1021.9
14.044 7.6382 38.191 1022.1
14.255 7.7387 38.693 1022.4
14.466 7.8392 39.196 1022.6
14.677 7.9397 39.698 1022.9
14.888 8.0402 40.201 1023.1
15.099 8.1407 40.704 1023.3
15.31 8.2412 41.206 1023.6
15.521 8.3417 41.709 1023.8
15.733 8.4422 42.211 1024
15.944 8.5427 42.714 1024.3
16.155 8.6432 43.216 1024.5
16.366 8.7437 43.719 1024.7
16.577 8.8442 44.221 1025
16.788 8.9447 44.724 1025.2
16.999 9.0452 45.226 1025.4
17.21 9.1457 45.729 1025.7
17.421 9.2462 46.231 1025.9
17.632 9.3467 46.734 1026.1
17.844 9.4472 47.236 1026.3
18.055 9.5477 47.739 1026.6
18.266 9.6482 48.241 1026.8
18.477 9.7487 48.744 1027
18.688 9.8492 49.246 1027.2
18.899 9.9497 49.749 1027.4
19.11 10.05 50.251 1027.7
19.321 10.151 50.754 1027.9
19.532 10.251 51.256 1028.1
19.744 10.352 51.759 1028.3
19.955 10.452 52.261 1028.5
20.166 10.553 52.764 1028.7
20.377 10.653 53.266 1029
20.588 10.754 53.769 1029.2
20.799 10.854 54.271 1029.4
21.01 10.955 54.774 1029.6
21.221 11.055 55.276 1029.8
21.432 11.156 55.779 1030
21.644 11.256 56.281 1030.2
21.855 11.357 56.784 1030.4
22.066 11.457 57.286 1030.6
22.277 11.558 57.789 1030.8
22.488 11.658 58.291 1031
22.699 11.759 58.794 1031.3
22.91 11.859 59.296 1031.5
23.121 11.96 59.799 1031.7
23.332 12.06 60.302 1031.9
23.544 12.161 60.804 1032.1
23.755 12.261 61.307 1032.3
23.966 12.362 61.809 1032.5
24.177 12.462 62.312 1032.7
24.388 12.563 62.814 1032.9
24.599 12.663 63.317 1033.1
24.81 12.764 63.819 1033.3
25.021 12.864 64.322 1033.5
25.232 12.965 64.824 1033.6
25.444 13.065 65.327 1033.8
25.655 13.166 65.829 1034
25.866 13.266 66.332 1034.2
26.077 13.367 66.834 1034.4
26.288 13.467 67.337 1034.6
26.499 13.568 67.839 1034.8
26.71 13.668 68.342 1035
26.921 13.769 68.844 1035.2
27.132 13.869 69.347 1035.4
27.344 13.97 69.849 1035.6
27.555 14.07 70.352 1035.8
27.766 14.171 70.854 1035.9
27.977 14.271 71.357 1036.1
28.188 14.372 71.859 1036.3
28.399 14.472 72.362 1036.5
28.61 14.573 72.864 1036.7
28.821 14.673 73.367 1036.9
29.032 14.774 73.869 1037.1
29.244 14.874 74.372 1037.2
29.455 14.975 74.874 1037.4
29.666 15.075 75.377 1037.6
29.877 15.176 75.879 1037.8
30.088 15.276 76.382 1038
30.299 15.377 76.884 1038.1
30.51 15.477 77.387 1038.3
30.721 15.578 77.889 1038.5
30.932 15.678 78.392 1038.7
31.144 15.779 78.894 1038.9
31.355 15.879 79.397 1039
31.566 15.98 79.899 1039.2
31.777 16.08 80.402 1039.4
31.988 16.181 80.905 1039.6
32.199 16.281 81.407 1039.7
32.41 16.382 81.91 1039.9
32.621 16.482 82.412 1040.1
32.832 16.583 82.915 1040.3
33.044 16.683 83.417 1040.4
33.255 16.784 83.92 1040.6
33.466 16.884 84.422 1040.8
33.677 16.985 84.925 1040.9
33.888 17.085 85.427 1041.1
34.099 17.186 85.93 1041.3
34.31 17.286 86.432 1041.4
34.521 17.387 86.935 1041.6
34.732 17.487 87.437 1041.8
34.944 17.588 87.94 1042
35.155 17.688 88.442 1042.1
35.366 17.789 88.945 1042.3
35.577 17.889 89.447 1042.5
35.788 17.99 89.95 1042.6
35.999 18.09 90.452 1042.8
36.21 18.191 90.955 1042.9
36.421 18.291 91.457 1043.1
36.632 18.392 91.96 1043.3
36.844 18.492 92.462 1043.4
37.055 18.593 92.965 1043.6
37.266 18.693 93.467 1043.8
37.477 18.794 93.97 1043.9
37.688 18.894 94.472 1044.1
37.899 18.995 94.975 1044.2
38.11 19.095 95.477 1044.4
38.321 19.196 95.98 1044.6
38.532 19.296 96.482 1044.7
38.744 19.397 96.985 1044.9
38.955 19.497 97.487 1045
39.166 19.598 97.99 1045.2
39.377 19.698 98.492 1045.4
39.588 19.799 98.995 1045.5
39.799 19.899 99.497 1045.7
40.01 20 100 1045.8

bolzanow

p.s.

我也不知道如果只给定数据的话怎么开始搜索公式,反正我的1.0死活不认我说的话...
不知道哪里写错了....

shamohu

公式表达形式有错，如果不熟悉sum的用法，自己展开写就行了。
function u=(sum(j=1:3)(a[j]*exp(x[j]*s)))* (sum(j=1:3)(b[j]*exp(y[j]*t)))*(sum(j=1:3)(c[j]*exp(z[j]*p)));

结果还是不错的，但估计不是唯一的：

均方差(RMSE)        3.88327979791036E-02
残差平方和(SSE)        .301597239777175
相关系数(R)        .999995649102781
相关系数之平方(R^2)        .999991298224493

参数名        参数值
a1        8.08747033297878E-04
a2        -8.11147853084382E-05
a3        38.2730207917541
x1        2.33046115315307E-02
x2        2.46392056298638E-02
x3        -4.40996432816221E-02
b1        2.12294910640009E-02
b2        8.50730959694048E-19
b3        4.12921809808961
y1        -7.59729044632777E-02
y2        .22795811657472
y3        -9.15967468031746E-02
c1        -2.17007801098777E-02
c2        -.273176694188079
c3        5.53519254478785
z1        1.54676401075122E-02
z2        3.75219235195406E-02
z3        4.73459468609683E-02

[ 本帖最后由 shamohu 于 2007-6-6 09:42 编辑 ]

bolzanow

是不是3个变量的就不能自动搜索公式了呢？

wanglei5201118

是的，不过你可以自己设计公式在1stopt，也挺方便的