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

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » MATLAB应用四阶龙格库塔法求解反应动力学问题?
221/3返回列表123下一页
查看: 3348  |  回复: 21
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

anxiaoxin2

金虫 (小有名气)

[求助] MATLAB应用四阶龙格库塔法求解反应动力学问题? 已有2人参与

各位英雄大侠，你们好。本人在做铁水喷吹脱硫的优化，遇到一个问题，文献中是应用MATLAB四阶龙格库塔法分析的，本人对MATLAB不熟，希望大家能帮帮忙！
谢谢。。。
具体问题见附件。
回复此楼

» 本帖附件资源列表

  • 欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
    本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com
  • 附件 1 : 脱硫动力学微分方程.docx
    • 2015-03-05 14:00:32, 46.03 K

» 收录本帖的淘帖专辑推荐

动力学拟合

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2015-03-05 14:01:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

月只蓝

主管区长 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ...
感谢参与，应助指数 +1
anxiaoxin2: 金币+300, ★★★很有帮助, 非常有用。。。 2015-03-06 08:37:41
CODE:
function solve_odes
clear all;clc
y0=[ ];    % [%s] [%Mg]在t=0时刻的数值，即其初值，待填入
tspan=; %时间t的范围，待填入
[t,y]=ode45(@my_odes],[2 0])
figure(1),plot(t,y(:,1));     % [%s]对时间t作图
figure(2),plot(t,y(:,2);      % [%Mg]对时间t作图

function dydt=my_odes
beta= ;    %β的数值
Ip= ;    Ip的数值
......   %以下依次把已经参数的数值填入。

dydt(1)=        ;  % d[%s]/dt的表达式
dydt(2)=        ;  % d[%Mg]/dt的表达式
赞一下(1人) 回复此楼
MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。
4楼2015-03-05 21:00:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

CelestialCYJ

木虫 (小有名气)

【答案】应助回帖

感谢参与，应助指数 +1
你不给数值怎么算
赞一下 回复此楼

» 本帖已获得的红花（最新10朵）

anxiaoxin2
2楼2015-03-05 16:11:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

anxiaoxin2

金虫 (小有名气)
送红花一朵
引用回帖:
2楼: Originally posted by CelestialCYJ at 2015-03-05 16:11:26
你不给数值怎么算

这些都是参数，可以设定，方便以后调整，你有编写的程序吗？
赞一下 回复此楼
3楼2015-03-05 17:00:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

anxiaoxin2

金虫 (小有名气)
谢谢版主。。。这个不需要m文件之类的吗？
赞一下 回复此楼
5楼2015-03-06 08:37:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)
引用回帖:
5楼: Originally posted by anxiaoxin2 at 2015-03-06 08:37:00
谢谢版主。。。这个不需要m文件之类的吗？

你好，因为你没有各个变量或参数的具体数据，我给出的是一个m文件的框架，你把你的参数和方程填入，就能计算的。
赞一下 回复此楼
MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。
6楼2015-03-06 08:40:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

anxiaoxin2

金虫 (小有名气)
好的，我先试试，不懂的再来请教。。。
赞一下 回复此楼
7楼2015-03-06 08:46:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

anxiaoxin2

金虫 (小有名气)
引用回帖:
6楼: Originally posted by 月只蓝 at 2015-03-06 08:40:39
你好，因为你没有各个变量或参数的具体数据，我给出的是一个m文件的框架，你把你的参数和方程填入，就能计算的。...

你好，按你提供的输入到matlab中了，显示
dydt(1)=(6*b*Ip*U1*ip*Pm*Km*([S]*exp(-Tmix/T)-[S]CaOe))/(Ds*Wm*Ps)+(6*Q*irt*Pm*Kb*([S]*exp(-Tmix/T)-[S]Mge))/Wm/Db+0.4167*k*Pm*[S]*exp(-Tmix/T)*[Mg]*exp(-Tmix/T)+A*Pz*Ks*Ls*([S]*exp(-Tmix/T)-(S)/Ls)/Wm;
dydt(2)=100*B*U3*Ip/Wm-0.75*(0.4167*k*Pm*[S]*exp(-Tmix/T)*[Mg]*exp(-Tmix/T));
??? function solve_odes
    |
Error: Function definitions are not permitted in this
context.
不知道是哪里出错了？该怎么修改？另外输入这些就可以作出[S]对t的图形吗？
赞一下 回复此楼
8楼2015-03-17 17:06:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)
引用回帖:
8楼: Originally posted by anxiaoxin2 at 2015-03-17 17:06:06
你好，按你提供的输入到matlab中了，显示
dydt(1)=(6*b*Ip*U1*ip*Pm*Km*(*exp(-Tmix/T)-CaOe))/(Ds*Wm*Ps)+(6*Q*irt*Pm*Kb*(*exp(-Tmix/T)-Mge))/Wm/Db+0.4167*k*Pm**exp(-Tmix/T)**exp(-Tmix/T)+A*Pz*Ks*Ls*(*ex ...

在MATLAB中，[S]这种代码写法不对的，用单个字母代替，如果字母有重复，可以用字母+数字的属性，比如S1 S2 S_1 S_2
赞一下 回复此楼
MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。
9楼2015-03-17 18:05:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

anxiaoxin2

金虫 (小有名气)
引用回帖:
9楼: Originally posted by 月只蓝 at 2015-03-17 18:05:45
在MATLAB中，这种代码写法不对的，用单个字母代替，如果字母有重复，可以用字母+数字的属性，比如S1 S2 S_1 S_2...

谢谢。。。
你说的问题我都已修改。。。
function solve_odes
clear all;clc
y0=[ ];    % [%S] [%Mg]在t=0时刻的数值，即其初值，待填入
tspan=1200; %时间t的范围，待填入
[t,y]=ode45(@my_odes],[2 0])
figure(1),plot(t,y(:,1));     % [%s]对时间t作图
figure(2),plot(t,y(:,2);      % [%Mg]对时间t作图
这段代码中y0=[];此处有两个初始值，怎么填入？
后边的微分方程中是不是所有的这两个[%S] [%Mg]的初始值都要改成y0?
赞一下 回复此楼
10楼2015-03-18 09:10:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 anxiaoxin2 的主题更新
221/3返回列表123下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿中科大材料与化工，353分还有调剂学校吗 +10 否极泰来2026 2026-04-15 12/600 2026-04-17 17:54 by mapenggao
[教师之家] 山东双非院校考核超级无底线，领导幸灾乐祸，教师遭殃恐 +4 qut2026 2026-04-11 8/400 2026-04-17 16:10 by 会飞的猪157
[考研] 一志愿沪9，326求生物学调剂 +11 刘墨墨 2026-04-13 11/550 2026-04-17 15:35 by wanganpin
[考研] 材料相关专业344求调剂双非工科学校或课题组 +23 hualkop 2026-04-12 25/1250 2026-04-16 22:12 by SUSE_CL
[考研] 294求调剂 +14 淡然654321 2026-04-15 14/700 2026-04-16 21:01 by lpl364211
[考研] 307中医考研调剂 +6 于以采蘩 2026-04-14 6/300 2026-04-16 16:20 by qingfeng258
[考研] 297，工科调剂?河南农业大学本科 +14 河南农业大学-能 2026-04-14 14/700 2026-04-16 14:41 by dingyanbo1
[考研] 322求调剂 +8 123安康 2026-04-12 15/750 2026-04-16 11:07 by Espannnnnol
[考研] 327求调剂 +26 Xxjc1107. 2026-04-13 29/1450 2026-04-16 10:52 by Espannnnnol
[考研] 求调剂 +11 小聂爱学习 2026-04-11 15/750 2026-04-15 21:57 by noqvsozv
[考研] 0854调剂 +13 长弓傲 2026-04-12 16/800 2026-04-15 13:45 by fenglj492
[教师之家] 转长聘了 +7 简单化xn 2026-04-13 7/350 2026-04-14 14:50 by xindong
[考研] 考研调剂 +13 长弓傲 2026-04-13 14/700 2026-04-14 14:44 by zs92450
[考研] B区0809 ，数一英一，290 求调剂 +3 泠潍1111 2026-04-12 4/200 2026-04-13 20:35 by 学员JpLReM
[考研] 一志愿085802 323分求调剂 +13 drizzle_9 2026-04-12 14/700 2026-04-13 10:26 by Faiz5552
[考研] 346分，工科0854求调剂，专硕 +6 moser233 2026-04-12 7/350 2026-04-12 22:11 by fqwang
[考研] 331求调剂 +5 王国帅 2026-04-11 5/250 2026-04-11 22:56 by 溪涧流水
[考研] 359求调剂 +5 胃痉挛累了 2026-04-11 5/250 2026-04-11 19:55 by lbsjt
[考研] 求调剂，一志愿大连理工大学354分 +5 雨声余生 2026-04-11 6/300 2026-04-11 16:12 by 雨声余生
[考研] 调剂 +5 文道星台 2026-04-11 5/250 2026-04-11 15:01 by 凯凯要变帅
信息提示
关闭
请填处理意见
关闭