| 查看: 513 | 回复: 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 |
回复此楼» 猜你喜欢
PVA溶解
已经有2人回复
-
基金申请
已经有44人回复
-
物理学I论文润色/翻译怎么收费?
已经有91人回复
-
CSC与新西兰维多利亚大学PhD奖学金项目
已经有0人回复
-
新西兰Robinson研究所 招聘CSC公派访问人员
已经有0人回复
-
帮我的英语口语老师找学生
已经有0人回复
-
什么时候开奖？
已经有13人回复
-
散金币祈福
已经有94人回复
-
青基已中
已经有1人回复
-
散金币祈福
已经有22人回复
2楼2016-03-04 15:05:18
30