24小时热门版块排行榜

返回列表

你的CEO

铁虫 (小有名气)

应助: 2 (幼儿园)
金币: 175.6
帖子: 71
在线: 18小时
虫号: 3071920
注册: 2014-03-19
性别: GG
专业: 能源化工

[求助] 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

回复此楼