| 查看: 469 | 回复: 1 | ||
[求助]
1stopt代码求有正版软件的大神帮运行下，求计算结果，小弟不胜感激
|
|
注y’=dy/dl Variable k（1）、k（2）、k（3）、k（4）、A、B、C、D、E、F、G ODEFunction1 yO2=0.0769*(1+0.5*y(1)-0.5*y(3))-(0.5*y(1)+y(2)+0.5*y(3)); yMEOH=0.3077*(1+0.5*y(1)-0.5*y(3))-(3*y(1)+2*y(2)+y(3)+2*y(4)); y(1)’=0.0432*(2+y(1)-y(3))/(4*0.154284)*((2+y(1))*k(1)*6^(A+B)*yMEOH^A*yO2^B-y(1)*k(3)*6^(E+F)*yMEOH^E*yO2^F); y(2)’=0.0432*(2+y(1)-y(3))/(4*0.154284)*(2*k(2)*6^(C+D)*yMEOH^C*yO2^D+ y(2)*k(1)*6^(A+B)*yMEOH^A*yO2^B-y(2)*k(3)*6^(E+F)*yMEOH^E*yO2^F); y(3)’=0.0432*(2+y(1)-y(3))/(4*0.154284)*((2-y(3))*k(3)*6^(E+F)*yMEOH^E*yO2^F+ y(3)*k(1)*6^(A+B)*yMEOH^A*yO2^B); y(4)’=0.0432*(2+y(1)y(3))/(4*0.154284)*(2*k(4)*6^G*yMEOH^G+y(4)*k(1)*6^(A+B)*yMEOH^A*yO2^B-y(4)*k(3)*6^(E+F)*yMEOH^E*yO2^F); ]; Data1 L y（1） y（2） y（3） y（4） 0 0 0 0 0 0.466 0.0325 0.0016 0.0002 0.0002 ODEFunction2 yO2=0.0882*(1+0.5*y(1)-0.5*y(3))-(0.5*y(1)+y(2)+0.5*y(3)); yMEOH=0.2941*(1+0.5*y(1)-0.5*y(3))-(3*y(1)+2*y(2)+y(3)+2*y(4)); y(1)’=0.0432*(2+y(1)-y(3))/(4*0.161404)*((2+y(1))*k(1)*8^(A+B)*yMEOH^A*yO2^B-y(1)*k(3)*8^(E+F)*yMEOH^E*yO2^F); y(2)’=0.0432*(2+y(1)-y(3))/(4*0.161404)*(2*k(2)*8^(C+D)*yMEOH^C*yO2^D+ y(2)*k(1)*8^(A+B)*yMEOH^A*yO2^B-y(2)*k(3)*8^(E+F)*yMEOH^E*yO2^F); y(3)’=0.0432*(2+y(1)-y(3))/(4*0.161404)*((2-y(3))*k(3)*8^(E+F)*yMEOH^E*yO2^F+ y(3)*k(1)*8^(A+B)*yMEOH^A*yO2^B); y(4)’=0.0432*(2+y(1)-y(3))/(4*0.161404)*(2*k(4)*8^G*yMEOH^G+y(4)*k(1)*8^(A+B)*yMEOH^A*yO2^B-y(4)*k(3)*8^(E+F)*yMEOH^E*yO2^F); ]; Data2 L y（1） y（2） y（3） y（4） 0 0 0 0 0 0.466 0.0288 0.0014 0.0005 0.0002 ODEFunction3 yO2=0.1176*(1+0.5*y(1)-0.5*y(3))-(0.5*y(1)+y(2)+0.5*y(3)); yMEOH=0.2941*(1+0.5*y(1)-0.5*y(3))-(3*y(1)+2*y(2)+y(3)+2*y(4)); y(1)’=[0.0432*(2+y(1)-y(3))/(4*0.161404)*((2+y(1))*k(1)*5^(A+B)*yMEOH^A*yO2^B-y(1)*k(3)*5^(E+F)*yMEOH^E*yO2^F); y(2)’=0.0432*(2+y(1)-y(3))/(4*0.161404)*(2*k(2)*5^(C+D)*yMEOH^C*yO2^D+ y(2)*k(1)*5^(A+B)*yMEOH^A*yO2^B-y(2)*k(3)*5^(E+F)*yMEOH^E*yO2^F); y(3)’=0.0432*(2+y(1)-y(3))/(4*0.161404)*((2-y(3))*k(3)*5^(E+F)*yMEOH^E*yO2^F+ y(3)*k(1)*5^(A+B)*yMEOH^A*yO2^B); y(4)’=0.0432*(2+y(1)-y(3))/(4*0.161404)*(2*k(4)*5^G*yMEOH^G+y(4)*k(1)*5^(A+B)*yMEOH^A*yO2^B-y(4)*k(3)*5^(E+F)*yMEOH^E*yO2^F); ]; Data3 L y（1） y（2） y（3） y（4） 0 0 0 0 0 0.466 0.0255 0.0012 0.0003 0.0003 ODEFunction4 yO2=0.1429*(1+0.5*y(1)-0.5*y(3))-(0.5*y(1)+y(2)+0.5*y(3)); yMEOH=0.2857*(1+0.5*y(1)-0.5*y(3))-(3*y(1)+2*y(2)+y(3)+2*y(4)); y(1)’=[0.0432*(2+y(1)-y(3))/(4*0.166152)*((2+y(1))*k(1)*7^(A+B)*yMEOH^A*yO2^B-y(1)*k(3)*7^(E+F)*yMEOH^E*yO2^F); y(2)’=0.0432*(2+y(1)-y(3))/(4*0.166152)*(2*k(2)*7^(C+D)*yMEOH^C*yO2^D+ y(2)*k(1)*7^(A+B)*yMEOH^A*yO2^B-y(2)*k(3)*7^(E+F)*yMEOH^E*yO2^F); y(3)’=0.0432*(2+y(1)-y(3))/(4*0.166152)*((2-y(3))*k(3)*7^(E+F)*yMEOH^E*yO2^F+ y(3)*k(1)*7^(A+B)*yMEOH^A*yO2^B); y(4)’=0.0432*(2+y(1)-y(3))/(4*0.166152)*(2*k(4)*7^G*yMEOH^G+y(4)*k(1)*7^(A+B)*yMEOH^A*yO2^B-y(4)*k(3)*7^(E+F)*yMEOH^E*yO2^F); ]; Data4 L y（1） y（2） y（3） y（4） 0 0 0 0 0 0.466 0.0264 0.0015 0.0004 0.0004; ODEFunction5 yO2=0.0857*(1+0.5*y(1)-0.5*y(3))-(0.5*y(1)+y(2)+0.5*y(3)); yMEOH=0.1429*(1+0.5*y(1)-0.5*y(3))-(3*y(1)+2*y(2)+y(3)+2*y(4)); y(1)’=[0.0432*(2+y(1)-y(3))/(4*0.332303)*((2+y(1))*k(1)*9^(A+B)*yMEOH^A*yO2^B-y(1)*k(3)*9^(E+F)*yMEOH^E*yO2^F); y(2)’=0.0432*(2+y(1)-y(3))/(4*0.332303)*(2*k(2)*9^(C+D)*yMEOH^C*yO2^D+ y(2)*k(1)*9^(A+B)*yMEOH^A*yO2^B-y(2)*k(3)*9^(E+F)*yMEOH^E*yO2^F); y(3)’=0.0432*(2+y(1)-y(3))/(4*0.332303)*((2-y(3))*k(3)*9^(E+F)*yMEOH^E*yO2^F+ y(3)*k(1)*9^(A+B)*yMEOH^A*yO2^B); y(4)’=0.0432*(2+y(1)-y(3))/(4*0.332303)*(2*k(4)*9^G*yMEOH^G+y(4)*k(1)*9^(A+B)*yMEOH^A*yO2^B-y(4)*k(3)*9^(E+F)*yMEOH^E*yO2^F); ]; Data5.. L y（1） y（2） y（3） y（4） 0 0 0 0 0 0.466 0.0109 0.0006 0.0002 0.0003 五个方程联立（不能分开来求解）同时求解出参数k（1）、k（2）、k（3）、k（4）、A、B、C、D、E、F、G |
回复此楼» 猜你喜欢
光学工程学硕调剂信息
已经有32人回复
-
欢迎加入课题组
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有112人回复
-
散金币，求好运，祝面上顺利！
已经有79人回复
-
华东师范大学芯片设计徐珑真实水平咋样
已经有1人回复
-
调剂
已经有0人回复
-
【新加坡】纳米电子器件项目组有“联合培养博士生”名额
已经有0人回复
-
中德博士后2026项目5月30日截止- 国内博士后入站、德国两年全职
已经有0人回复
-
面向全球申请人的中德博士后2026项目5月30日截止- 国内博后入站德国两年全职
已经有0人回复
-
【新加坡】纳米电子器件项目组有“联合培养博士生”名额
已经有2人回复
-
重庆交大材料学院26年硕士招生复试实施细则已出。欢迎加入光子学微结构与器件课题组。
已经有0人回复
2楼2016-03-04 15:05:18
