| 查看: 386 | 回复: 1 | ||
[求助]
烦请@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 自己选参数值，选了很久了，帮帮我吧！ |
回复此楼» 本帖@通知
» 猜你喜欢
售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急
已经有7人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急
已经有8人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急
已经有8人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急
已经有11人回复
-
为什么中国大学工科教授们水了那么多所谓的顶会顶刊，但还是做不出宇树机器人？
已经有8人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急
已经有5人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急
已经有6人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急
已经有9人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急
已经有4人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急
已经有5人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 急求！！！求用1stOpt软件微分求解方程参数。
已经有21人回复
- 刚开始用1stopt拟合曲线，但执行不了计算，求大神帮看一下，感谢感谢~
已经有4人回复
- 求助大神帮忙用1stopt拟合下模型参数
已经有9人回复
- 1stopt自变量超多时是否有简单语法？AND 自变量多 数据量大的话1stopt都会算的特慢么
已经有7人回复
- 谁能帮忙用1Stopt运行下程序呀？谢谢，我的1.5版本的运行不了啊
已经有6人回复
- 求大神跑下1stOpt，一个sharemodel，谢谢！
已经有7人回复
- 1stopt拟合的动力学参数相关系数不高，有哪些原因
已经有3人回复
- 求助1stopt拟合动力学参数
已经有4人回复
- 1stopt求最小值
已经有37人回复
- 一个1stOpt的初级问题：复数方程组计算问题（有代码）
已经有5人回复
- 1stopt做最小二乘法回归参数，请求大家帮忙指导下
已经有3人回复
- 用1stopt拟合非线性方程结果与其他软件拟合结果差异大
已经有5人回复
- 大家在使用1stOpt进行数据拟合的时候有没有遇到过这种情况？
已经有8人回复
- 【求助】用matlab最优化方法进行参数拟合
已经有17人回复
2楼2015-05-19 16:46:27
