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

24小时热门版块排行榜    

CyRhmU.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » MATLAB应用四阶龙格库塔法求解反应动力学问题?
221/3返回列表123下一页
查看: 3108  |  回复: 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下一页
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭