请大神帮忙代跑下1stOpt,sharemodel,不胜感激!
请大神帮帮忙,看看结果。之前尝试过用matlab,没有做出来,麻烦大家了,谢谢!
Title "data";
SharedModel;
//Parameters ;
//Variable ;
//Function ;
Variable x1,x2,x3,x4,x5, y1,y2,y3,y4;
Parameters a(1:5)[0,1,0],b(1:5)[0,1,0],c(1:5)[0,1,0],d(1:5)[0,1,0];
Function y1=((a1*x1+a2*x2+a3*x3+a4*x4)*2+a5*x5*1.1)/((2*a1+2*b1+c1+d1)*x1+(2*a2+2*b2+c2+d2)*x2+(2*a3+2*b3+c3+d3)*x3+(2*a4+2*b4+c4+d4)*x4+0.55*(2*a5+2*b5+c5+d5)*x5 ) ;
y2=((b1*x1+b2*x2+b3*x3+b4*x4)*2+b5*x5*1.1)/((2*a1+2*b1+c1+d1)*x1+(2*a2+2*b2+c2+d2)*x2+(2*a3+2*b3+c3+d3)*x3+(2*a4+2*b4+c4+d4)*x4+0.55*(2*a5+2*b5+c5+d5)*x5 ) ;
y3=((c1*x1+c2*x2+c3*x3+c4*x4)+c5*x5*0.55)/((2*a1+2*b1+c1+d1)*x1+(2*a2+2*b2+c2+d2)*x2+(2*a3+2*b3+c3+d3)*x3+(2*a4+2*b4+c4+d4)*x4+0.55*(2*a5+2*b5+c5+d5)*x5 ) ;
y4=((d1*x1+d2*x2+d3*x3+d4*x4)+d5*x5*0.55)/((2*a1+2*b1+c1+d1)*x1+(2*a2+2*b2+c2+d2)*x2+(2*a3+2*b3+c3+d3)*x3+(2*a4+2*b4+c4+d4)*x4+0.55*(2*a5+2*b5+c5+d5)*x5 ) ;
a1+a2+a3+a4+a5=1;
b1+b2+b3+b4+b5=1;
c1+c2+c3+c4+c5=1;
d1+d2+d3+d4+d5=1;
data;
11.858 13.995 3.817 8.024 8.287 68.868 16.960 12.068 2.104
12.694 14.321 5.299 9.842 8.463 68.296 17.570 12.071 2.063
14.330 14.318 6.004 10.373 8.407 67.854 17.740 11.962 2.445
14.871 13.948 7.240 13.372 8.318 67.872 16.775 11.893 3.459
13.209 17.904 2.960 4.105 9.691 59.824 19.010 11.581 9.585
11.151 14.427 5.522 8.565 11.458 63.223 24.273 12.505 0.000
8.105 19.490 5.536 6.619 12.331 62.000 25.000 13.000 0.000
9.359 15.399 7.143 10.140 13.148 66.860 20.550 11.810 0.780
13.892 13.384 6.674 10.057 13.753 69.856 20.693 8.147 1.303
11.746 13.221 4.635 8.750 10.203 68.304 21.053 10.912 0.000
12.950 13.088 4.593 8.672 10.118 70.333 20.281 7.479 1.906
8.125 12.140 9.263 12.310 10.472 65.793 16.534 13.378 4.295
8.251 12.171 5.057 8.925 10.457 66.670 16.560 12.398 4.371
7.712 12.614 7.554 11.831 11.708 65.912 16.733 13.008 4.346
8.900 14.505 5.199 10.113 9.306 67.250 19.590 10.230 2.930
8.056 12.077 5.594 10.502 12.525 65.830 21.080 9.760 3.330
8.184 12.415 5.180 10.269 10.904 66.750 20.550 9.510 3.190
8.190 12.940 5.052 9.315 11.566 66.820 20.860 9.430 2.890
9.050 15.764 10.574 15.587 14.220 64.200 23.900 8.010 3.890
8.860 19.440 7.930 10.954 13.203 63.766 21.869 11.028 3.336
9.509 20.092 6.880 8.417 12.552 65.200 20.900 10.500 3.400
9.009 19.119 6.601 8.081 12.068 65.185 21.019 10.926 2.870
9.217 19.525 6.718 8.222 12.271 62.692 23.077 10.288 3.942
15.894 20.004 4.047 6.365 18.113 57.200 23.300 13.600 5.900
16.069 14.672 4.377 6.649 17.069 59.500 22.500 12.600 5.400
6.548 17.484 5.659 8.404 6.249 61.785 15.780 13.300 9.135
5.091 10.091 4.163 5.277 4.237 60.600 19.600 14.900 4.900
5.091 10.091 4.163 5.277 4.237 58.900 20.800 14.800 5.500
11.129 8.734 3.047 5.915 6.011 64.400 19.000 13.200 3.400
11.761 10.213 2.508 5.969 7.374 69.600 17.800 10.200 2.400
10.344 9.079 3.449 5.578 7.116 64.300 19.500 13.000 3.200
12.937 13.609 4.677 7.099 9.902 65.400 19.300 12.100 3.200
4.999 10.907 3.983 6.079 7.896 61.722 26.919 11.349 0.011
4.999 10.907 3.983 6.079 7.896 61.249 24.938 13.803 0.010
4.999 10.907 3.983 6.079 7.896 61.774 26.814 11.412 0.000
9.013 11.339 6.278 9.036 13.708 64.191 20.155 8.209 7.445
9.198 11.515 6.441 9.315 13.897 67.353 19.908 8.727 4.011
8.943 11.385 5.402 9.544 14.172 64.064 20.749 8.343 6.844
9.255 16.012 6.683 7.346 13.164 64.191 20.155 8.209 7.445
9.843 14.223 6.595 7.571 13.281 64.353 19.908 8.727 7.011
11.801 13.682 8.236 11.741 18.655 64.064 20.749 8.343 6.844
11.892 10.770 3.207 8.018 8.472 66.786 22.251 4.130 6.834
11.892 10.770 3.207 8.018 8.472 68.271 15.870 3.254 12.604
11.892 10.770 3.207 8.018 8.472 68.095 16.975 4.444 10.487
14.547 18.773 5.404 7.481 13.191 63.738 21.869 11.028 3.364
14.547 18.773 5.404 7.481 13.191 65.200 20.900 10.500 3.400
14.547 18.773 5.404 7.481 13.191 65.185 21.019 10.926 2.870
14.547 18.773 5.404 7.481 13.191 62.692 23.077 10.288 3.942
13.346 14.612 4.528 6.034 9.997 65.187 18.353 11.905 4.555
13.346 14.612 4.528 6.034 9.997 65.264 17.243 12.454 5.040
13.346 14.612 4.528 6.034 9.997 65.902 16.480 12.498 5.121
13.378 12.611 4.200 6.184 12.564 59.375 19.479 11.526 9.620
3.795 8.539 3.028 6.394 7.759 59.505 22.404 10.429 7.662
3.506 7.933 2.827 5.968 7.233 60.187 21.966 10.304 7.543
2.832 6.494 2.341 4.940 5.967 60.572 21.658 10.247 7.523
13.209 17.904 2.960 4.105 9.691 63.751 18.305 10.189 7.755
5.684 10.487 10.946 13.998 10.210 69.011 20.416 10.573 0.000
8.750 9.841 6.784 10.620 9.554 66.759 22.813 10.427 0.000
5.373 15.726 4.812 8.413 12.130 65.190 22.870 8.610 3.330
7.144 15.969 5.698 8.084 12.882 65.160 22.540 8.860 3.440
7.833 15.128 6.467 8.246 12.696 64.850 22.690 9.060 3.400
6.974 18.770 6.177 8.317 12.525 64.300 25.390 8.200 2.110
10.013 12.033 4.842 8.356 12.988 69.130 17.765 9.925 3.180
9.897 11.934 4.821 8.323 12.919 68.822 18.184 9.787 3.207
9.670 11.303 4.809 8.347 12.746 68.983 17.698 10.115 3.204
9.379 10.993 4.691 8.145 12.437 68.957 18.239 9.656 3.148
11.077 10.683 3.539 7.341 12.879 67.336 18.729 10.268 3.668
11.703 11.236 3.695 7.668 13.454 72.235 20.465 11.615 0.000
11.014 10.628 3.523 7.308 12.820 68.670 19.635 10.697 0.998
11.577 11.126 3.664 7.604 13.340 66.911 18.429 9.853 4.806
10.352 6.911 6.708 7.805 8.548 69.400 19.500 9.400 1.700
14.899 10.382 4.866 6.628 15.907 65.100 22.300 11.300 1.300
14.160 5.723 4.904 6.694 14.567 68.800 21.300 10.800 0.000
6.071 14.286 5.388 8.955 10.314 62.058 22.380 10.923 4.639
4.573 13.793 5.335 8.831 7.745 62.067 21.776 10.888 5.268
5.999 14.253 5.387 8.949 10.142 62.225 21.541 11.135 5.099
10.760 7.998 5.388 6.951 8.852 68.857 20.128 9.371 1.645
10.730 6.746 5.261 6.891 8.933 70.068 19.645 9.247 1.040
10.730 6.746 5.261 6.891 8.933 70.036 20.372 8.764 0.829
10.730 6.746 5.261 6.891 8.933 69.494 20.514 9.102 0.890
16.459 17.613 3.407 4.448 11.853 63.859 18.331 10.141 7.669 返回小木虫查看更多
ConstStr
a1=1-(a2+a3+a4+a5),
b1=1-(b2+b2+b4+b5),
c1=1-(c2+c3+c4+c5),
d1=1-(d2+d3+d4+d5);
Parameters a(2:5)[0,1],b(2:5)[0,1],c(2:5)[0,1],d(2:5)[0,1];
SharedModel;
Variable x1,x2,x3,x4,x5, y1,y2,y3,y4;
Function
y1=((a1*x1+a2*x2+a3*x3+a4*x4)*2+a5*x5*1.1)/((2*a1+2*b1+c1+d1)*x1+(2*a2+2*b2+c2+d2)*x2+(2*a3+2*b3+c3+d3)*x3+(2*a4+2*b4+c4+d4)*x4+0.55*(2*a5+2*b5+c5+d5)*x5);
y2=((b1*x1+b2*x2+b3*x3+b4*x4)*2+b5*x5*1.1)/((2*a1+2*b1+c1+d1)*x1+(2*a2+2*b2+c2+d2)*x2+(2*a3+2*b3+c3+d3)*x3+(2*a4+2*b4+c4+d4)*x4+0.55*(2*a5+2*b5+c5+d5)*x5);
y3=((c1*x1+c2*x2+c3*x3+c4*x4)+c5*x5*0.55)/((2*a1+2*b1+c1+d1)*x1+(2*a2+2*b2+c2+d2)*x2+(2*a3+2*b3+c3+d3)*x3+(2*a4+2*b4+c4+d4)*x4+0.55*(2*a5+2*b5+c5+d5)*x5);
y4=((d1*x1+d2*x2+d3*x3+d4*x4)+d5*x5*0.55)/((2*a1+2*b1+c1+d1)*x1+(2*a2+2*b2+c2+d2)*x2+(2*a3+2*b3+c3+d3)*x3+(2*a4+2*b4+c4+d4)*x4+0.55*(2*a5+2*b5+c5+d5)*x5);
data;
11.858 13.995 3.817 8.024 8.287 68.868 16.960 12.068 2.104
12.694 14.321 5.299 9.842 8.463 68.296 17.570 12.071 2.063
,
谢谢您,我按着您修改的再运行下试试。
网上下载的一般不能允许共享参数拟合。
是的,运行不出结果,看看能不能用matlab解决了
大神,能否帮我跑一下程序呢?
Title "NOLINER";
Parameters A,N,B;
Variable x1,y1,x2,y2,x3,y3;
SharedModel;
Function
y1=3*A*x1^(-N)+B;
y2=4*A*x2^(-N)+B;
y3=6*A*x3^(-N)+B;
Data;
x1,y1,x2,y2,x3,y3;
0.096 2400 0.1240 2760 0.1900 3180
0.110 2360 0.1380 2633 0.1970 3090
0.132 2150 0.1420 2590 0.2100 3020
0.145 2130 0.1550 2560 0.2460 2870
0.151 2100 0.1583 2520 0.2540 2840
0.168 1990 0.1980 2330 0.2549 2800
0.172 1980 0.2020 2320 0.2864 2660
0.193 1900 0.2100 2280 0.2895 2690
0.203 1890 0.2300 2160 0.2940 2660
0.212 1860 0.2520 2100 0.3700 2440
0.224 1840 0.2600 2080 0.3760 2460
0.240 1790 0.2940 2020 0.3852 2400
0.265 1750 0.3000 1990 0.4580 2295
0.279 1700 0.3300 1940 0.4610 2240
0.290 1715 0.3450 1935 0.4880 2220
金币不够五百,信息发不出去!