| 查看: 490 | 回复: 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 |
回复此楼» 猜你喜欢
大豆异黄酮分离
已经有0人回复
-
湖南大学材料学院急招2026年博士生，临时增加一名博士联培指标
已经有10人回复
-
物理学I论文润色/翻译怎么收费?
已经有141人回复
-
天津理工大学晶体材料全国重点实验室刘红军教授课题组招收博士生1-2名
已经有1人回复
-
中国科学院物理研究所谌志国研究员团队招收2027年博士研究生
已经有5人回复
-
2026年中德博士后交流项目 - 新型量子和磁性材料：材料制备表征和中子散射研究
已经有12人回复
-
26申博推荐：南京航空航天大学国际前沿院光学方向招收博士生！
已经有1人回复
-
如何从铁电相到顺电相。
已经有1人回复
2楼2016-03-04 15:05:18