当前位置: 首页 > 计算模拟 >请大神帮忙代跑下1stOpt,sharemodel,不胜感激!

请大神帮忙代跑下1stOpt,sharemodel,不胜感激!

作者 夏夜听雨love
来源: 小木虫 300 6 举报帖子
+关注

请大神帮帮忙,看看结果。之前尝试过用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 返回小木虫查看更多

今日热帖
  • 精华评论
  • lipenggg

    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

  • 夏夜听雨love

    引用回帖:
    2楼: Originally posted by lipenggg at 2018-03-19 10:36:04
    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),b(2:5),c(2:5),d(2:5);
    SharedModel;
    Variable x1,x2,x3,x4,x5, y1,y2,y3,y4;
    Function
    y ...

    谢谢您,我按着您修改的再运行下试试。

  • lipenggg

    网上下载的一般不能允许共享参数拟合。

  • 夏夜听雨love

    引用回帖:
    4楼: Originally posted by lipenggg at 2018-03-20 11:01:55
    网上下载的一般不能允许共享参数拟合。

    是的,运行不出结果,看看能不能用matlab解决了

  • Powerhot11

    引用回帖:
    4楼: Originally posted by lipenggg at 2018-03-20 11:01:55
    网上下载的一般不能允许共享参数拟合。

    大神,能否帮我跑一下程序呢?
    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

  • 680416

    金币不够五百,信息发不出去!

猜你喜欢
下载小木虫APP
与700万科研达人随时交流
  • 二维码
  • IOS
  • 安卓