应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 1stopt代码求有正版软件的大神帮运行下，求计算结果，小弟不胜感激
21/1返回列表
查看: 468  |  回复: 1
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

你的CEO

铁虫 (小有名气)

[求助] 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
回复此楼

» 猜你喜欢
未来为你而来！
1楼 2016-03-04 15:04:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

你的CEO

铁虫 (小有名气)
求大神帮看下，顺便用软件算下
赞一下 回复此楼
未来为你而来！
2楼2016-03-04 15:05:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 你的CEO 的主题更新
21/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 26申博自荐 +3 whh869393 2026-03-24 3/150 2026-03-24 09:55 by 21018060
[考研] 生物学学硕求调剂 +5 小羊睡着了? 2026-03-23 5/250 2026-03-24 06:42 by ilovexiaobin
[考研] 一志愿华东理工大学081700，初试分数271 +4 kotoko_ik 2026-03-23 5/250 2026-03-23 23:23 by 呆呆师姐
[考研] 材料专业求调剂 +11 hanamiko 2026-03-18 11/550 2026-03-23 23:12 by peike
[考研] 0854 考研调剂 招生了！AI 方向 +4 pk3725069 2026-03-19 16/800 2026-03-23 23:09 by 汪！？！
[考研] 一志愿陕师大生物学071000，298分，求调剂 +3 SYA！ 2026-03-23 3/150 2026-03-23 19:09 by macy2011
[考研] 316求调剂 +7 梁茜雯 2026-03-19 7/350 2026-03-23 16:21 by lingjue
[考研] 材料与化工考研调剂 +4 孅華 2026-03-22 4/200 2026-03-23 16:13 by 一休哥FU
[考研] 招08考数学 +6 laoshidan 2026-03-20 14/700 2026-03-23 14:37 by 15614359529
[考研] 求调剂材料学硕080500，总分289分 5+3 @taotao 2026-03-19 21/1050 2026-03-23 10:17 by 冠c哥
[考研] 一志愿070300浙大化学358分，求调剂！ +4 酥酥鱼.. 2026-03-21 4/200 2026-03-23 08:12 by Iveryant
[考研] 一志愿华中科技大学071000，求调剂 +4 沿岸有贝壳6 2026-03-21 4/200 2026-03-22 07:21 by ilovexiaobin
[考研] 南昌大学材料专硕311分求调剂 +6 77chaselx 2026-03-20 6/300 2026-03-21 07:24 by JourneyLucky
[考研] 一志愿武理材料305分求调剂 +6 想上岸的鲤鱼 2026-03-18 7/350 2026-03-21 01:03 by JourneyLucky
[考研] 一志愿 西北大学 ，070300化学学硕，总分287，双非一本，求调剂。 +3 晨昏线与星海 2026-03-18 3/150 2026-03-21 00:46 by JourneyLucky
[考研] 317求调剂 +5 申子申申 2026-03-19 9/450 2026-03-20 22:26 by JourneyLucky
[考研] 一志愿吉林大学材料学硕321求调剂 +11 Ymlll 2026-03-18 15/750 2026-03-20 19:40 by 丁丁*
[论文投稿] 申请回稿延期一个月，编辑同意了。但系统上的时间没变，给编辑又写邮件了，没回复 10+3 wangf9518 2026-03-17 4/200 2026-03-19 23:55 by babero
[考研] 320求调剂0856 +3 不想起名字112 2026-03-19 3/150 2026-03-19 22:53 by 学员8dgXkO
[考研] 收复试调剂生 +4 雨后秋荷 2026-03-18 4/200 2026-03-18 14:16 by elevennnne
信息提示
关闭
请填处理意见
关闭