24小时热门版块排行榜

返回列表

1064126551

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 916
帖子: 15
在线: 25.4小时
虫号: 3099587
注册: 2014-03-29
专业: 分离过程

[求助] 烦请@dingd大神用1stOpt软件运行下程序-参数拟合

自己编的程序如下：
VarConstant N0=[0.479,0.314,0.281,0.314,0.314,0.281,0.281,0.314,0.314,0.314,0.314,0.314,0.314],
                  P=[0.1,0.1,0.1,0.2,0.3,0.2,0.3,0.2,0.1,0.3,0.1,0.2,0.3],
                  T=[433.15,433.15,433.15,433.15,433.15,433.15,433.15,423.15,423.15,423.15,443.15,443.15,438.15],
                  x1=[0.526,0.519,0.528,0.519,0.519,0.528,0.528,0.52,0.52,0.52,0.52,0.52, 0.52],
                  x2=[0.335,0.341,0.334,0.341,0.341,0.334 ,0.334,0.34,0.34,0.34,0.34,0.34,0.34],
                  x3=[0.139,0.141,0.137,0.141,0.141 ,0.137,0.137,0.14,0.14,0.14,0.14,0.14,0.14];
ConstStr r1=k1*exp(-k2/T)*(P^1.21588)*(y1*(1.5*x1-2)+0.5*y2*x1+x1)^(-0.82016)* (y1*(1.5*x2-1)+y2*(0.5*x2-1)+x2)^1.01727*(y1*(1.5*x3-0.5)+y2*(0.5*x3-0.5)+x3)^1.01877,
         r2=k3*exp(-k4/T)*P^(1.15102)*(y1*(1.5*x2-1)+y2*(0.5*x2-1)+x2)^0.52552*(y1*(1.5*x3-0.5)+y2*(0.5*x3-0.5)+x3)^0.6255;
ParameterDomain = [0,];
Variable  w,y1,y2;
ODEFunction    y1'=(1+1.5*y1+0.5*y2)*((1.5*y1+1)*r1 +0.5*y1*r2 )/N0;
               y2'=(1+1.5*y1+0.5*y2)*((1+0.5*y2)*r2+1.5*y2*r1 )/N0;
Data;  //1
0       0    0
0.42    0.021361816       0.007788162
Data; //2
0       0       0
0.42    0.030774673       0.019455253
Data;    //3
0       0          0
0.42    0.023121387       0.008477842
Data;    //4
0       0          0
0.42    0.023424294       0.02777107
Data;    //5
0       0          0
0.42    0.023082793       0.049217856
Data;    //6
0       0          0
0.42    0.024403619       0.026323005
Data;    //7
0       0          0
0.42    0.023542354       0.053025303
Data;    //8
0       0          0
0.42    0.007462037       0.023348255
Data;    //9
0       0          0
0.42    0.017986431       0.025487851
Data;    //10
0       0             0
0.42    0.014846134       0.102384543
Data;    //11
0          0          0
0.42       0.00879869       0.118733604
Data;       //12
0       0             0
0.42    0.014725875       0.373568961
Data;       //13
0       0             0
0.42    0.01665493       0.017564107
自己选参数值，选了很久了，帮帮我吧！

回复此楼