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

24小时热门版块排行榜    

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

月只蓝

主管区长 (职业作家)
内容已删除
回复此楼
MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。
11楼2015-03-18 10:32:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

anxiaoxin2

金虫 (小有名气)
引用回帖:
11楼: Originally posted by 月只蓝 at 2015-03-18 10:32:46
比如t=0时，初值为1 2；那么  y0=;...

时间这样填可以不？
tspan=(0:1:1200); %时间t的范围，待填入
赞一下 回复此楼
12楼2015-03-18 11:27:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)
内容已删除
回复此楼
MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。
13楼2015-03-18 12:23:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

anxiaoxin2

金虫 (小有名气)
引用回帖:
4楼: Originally posted by 月只蓝 at 2015-03-05 21:00:42
function solve_odes
clear all;clc
y0=;    %  在t=0时刻的数值，即其初值，待填入
tspan=; %时间t的范围，待填入
=ode45(@my_odes],)
figure(1),plot(t,y(:,1));     % 对时间t作图
figure(2),plot(t,y(:, ...

你好，我在以下表达式中[%S],[%Mg]该如何表示？直接写y1,y2吗？dydt(1)=        ;  % d[%s]/dt的表达式
dydt(2)=        ;  % d[%Mg]/dt的表达式
赞一下 回复此楼
14楼2015-03-18 14:09:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)
引用回帖:
14楼: Originally posted by anxiaoxin2 at 2015-03-18 14:09:03
你好，我在以下表达式中,该如何表示？直接写y1,y2吗？dydt(1)=        ;  % d/dt的表达式
dydt(2)=        ;  % d/dt的表达式...

[%S],[%Mg]可直接写成y1 y2
代码按如下

function dydt=my_odes(t,y)
y1=y(1);
y2=y(2);
dydt(1)=        ;  % d[%s]/dt的表达式
dydt(2)=        ;  % d[%Mg]/dt的表达式
.......
赞一下 回复此楼

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

anxiaoxin2
MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。
15楼2015-03-18 14:35:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

anxiaoxin2

金虫 (小有名气)
送红花一朵
引用回帖:
15楼: Originally posted by 月只蓝 at 2015-03-18 14:35:28
,可直接写成y1 y2
代码按如下

function dydt=my_odes(t,y)
y1=y(1);
y2=y(2);
dydt(1)=        ;  % d/dt的表达式
dydt(2)=        ;  % d/dt的表达式
..........

版主，这个还是运行不了啊。。。
这是代码，你帮我看看呗
function solve_odes
clear all;clc
y0=[0.03,0.0001];    % [%S] [%Mg]在t=0时刻的数值，即其初值，待填入
tspan=[0:1:1200]; %时间t的范围，待填入
[t,y]=ode45(‘my_odes’, tspan, y0);
figure(1),plot(t,y(:,1));     % [%s]对时间t作图
figure(2),plot(t,y(:,2));      % [%Mg]对时间t作图

function dydt=my_odes(t)
b= 0.36392;   %β脱硫剂穿透比
Ip=0.233;  %脱硫剂喷粉速率
U1=0.4;  %CaO在脱硫剂中的质量分数
U3=0.6;  %Mg在脱硫剂中的质量分数
ip=25.85;  %粉剂在铁水中停留时间
irt=12.45;  %Mg气泡在铁水中的停留时间
Pm=7138;  %铁水密度
Ps=1380;  %脱硫剂密度       
Pz=3500;  %顶渣密度
Wm=210000;  %铁水重量
Wso=5000;  %初始渣量
Ds=0.001;  %脱硫剂径粒
Db=9.75E-05;  %Mg气泡直径
Km=1.42313E-05;  %铁水中硫的传质系数
Kb=0.001031861;  %硫在Mg气泡上的传质系数
Ks=4.19217E-09;  %硫在渣中的传质系数
k=5.20879E-05;  %表观反应常数
Tmix=240;  %混匀时间
Q=0.02777778;  %喷吹气量
A=8.24;  %渣金反应面积
Ls=169.6743082;  %硫的分配系数
B=0.9;  %Mg在铁液中的溶解比例
n=0.4;  %粉剂进入渣层的有效系数
SCaOe=5.71697E-05;  %CaO脱硫平衡时硫含量
SMge=1.42026E-05;  %Mg脱硫平衡时硫含量
S4=0.5;  %初始渣中硫含量
y1=y(1);
y2=y(2);
S2=(S4*Wso+(S3-y1*exp(-Tmix/t))*Wm)/(Wso+Ip*t*n);

dydt(1)=(6*b*Ip*U1*ip*Pm*Km*(y1*exp(-Tmix/t)-SCaOe))/(Ds*Wm*Ps)+(6*Q*irt*Pm*Kb*(y1*exp(-Tmix/t)-SMge))/Wm/Db+0.4167*k*Pm*y1*exp(-Tmix/t)*y2*exp(-Tmix/t)+A*Pz*Ks*Ls*(y1*exp(-Tmix/t)-S2/Ls)/Wm;
dydt(2)=100*B*U3*Ip/Wm-0.75*(0.4167*k*Pm*y1*exp(-Tmix/t)*y2*exp(-Tmix/t));
谢了。。。
赞一下 回复此楼
16楼2015-03-18 15:26:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)
引用回帖:
16楼: Originally posted by anxiaoxin2 at 2015-03-18 15:26:43
版主，这个还是运行不了啊。。。
这是代码，你帮我看看呗
function solve_odes
clear all;clc
y0=;    %  在t=0时刻的数值，即其初值，待填入
tspan=; %时间t的范围，待填入
=ode45(‘my_odes’, tspan, y0) ...

S3没有定义，由于不知道S3的数值，以下代码中S3=0.3；程序已经可以跑通
CODE:
function solve_odes
clear all;clc
y0=[0.03,0.0001];    % [%S] [%Mg]在t=0时刻的数值，即其初值，待填入
tspan=[0:1:100]; %时间t的范围，待填入
[t,y]=ode15s(@my_odes, tspan, y0);
figure(1),plot(t,y(:,1));     % [%s]对时间t作图
figure(2),plot(t,y(:,2));      % [%Mg]对时间t作图

function dydt=my_odes(t,y)
b= 0.36392;   %β脱硫剂穿透比
Ip=0.233;  %脱硫剂喷粉速率
U1=0.4;  %CaO在脱硫剂中的质量分数
U3=0.6;  %Mg在脱硫剂中的质量分数
ip=25.85;  %粉剂在铁水中停留时间
irt=12.45;  %Mg气泡在铁水中的停留时间
Pm=7138;  %铁水密度
Ps=1380;  %脱硫剂密度        
Pz=3500;  %顶渣密度
Wm=210000;  %铁水重量
Wso=5000;  %初始渣量
Ds=0.001;  %脱硫剂径粒
Db=9.75E-05;  %Mg气泡直径
Km=1.42313E-05;  %铁水中硫的传质系数
Kb=0.001031861;  %硫在Mg气泡上的传质系数
Ks=4.19217E-09;  %硫在渣中的传质系数
k=5.20879E-05;  %表观反应常数
Tmix=240;  %混匀时间
Q=0.02777778;  %喷吹气量
A=8.24;  %渣金反应面积
Ls=169.6743082;  %硫的分配系数
B=0.9;  %Mg在铁液中的溶解比例
n=0.4;  %粉剂进入渣层的有效系数
SCaOe=5.71697E-05;  %CaO脱硫平衡时硫含量
SMge=1.42026E-05;  %Mg脱硫平衡时硫含量
S4=0.5;  %初始渣中硫含量
y1=y(1);
y2=y(2);

S3=0.1;

S2=(S4*Wso+(S3-y1*exp(-Tmix/t))*Wm)/(Wso+Ip*t*n);



dydt(1)=(6*b*Ip*U1*ip*Pm*Km*(y1*exp(-Tmix/t)-SCaOe))/(Ds*Wm*Ps)+(6*Q*irt*Pm*Kb*(y1*exp(-Tmix/t)-SMge))/Wm/Db+0.4167*k*Pm*y1*exp(-Tmix/t)*y2*exp(-Tmix/t)+A*Pz*Ks*Ls*(y1*exp(-Tmix/t)-S2/Ls)/Wm;
dydt(2)=100*B*U3*Ip/Wm-0.75*(0.4167*k*Pm*y1*exp(-Tmix/t)*y2*exp(-Tmix/t));

dydt=dydt';

MATLAB应用四阶龙格库塔法求解反应动力学问题?
附图1.png
赞一下 回复此楼

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

anxiaoxin2
MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。
17楼2015-03-18 16:29:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)
引用回帖:
17楼: Originally posted by 月只蓝 at 2015-03-18 16:29:34
S3没有定义，由于不知道S3的数值，以下代码中S3=0.3；程序已经可以跑通

function solve_odes
clear all;clc
y0=;    %  在t=0时刻的数值，即其初值，待填入
tspan=; %时间t的范围，待填入
=ode15s(@my_odes ...

S3=0.1；不是S3=0.3哈。
赞一下 回复此楼
MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。
18楼2015-03-18 16:30:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

anxiaoxin2

金虫 (小有名气)
送红花一朵
引用回帖:
17楼: Originally posted by 月只蓝 at 2015-03-18 16:29:34
S3没有定义，由于不知道S3的数值，以下代码中S3=0.3；程序已经可以跑通

function solve_odes
clear all;clc
y0=;    %  在t=0时刻的数值，即其初值，待填入
tspan=; %时间t的范围，待填入
=ode15s(@my_odes ...

版主，早上好啊。。。以上代码是怎么运行的？具体操作步骤能说下吗？我这还是不会运行。。。
赞一下 回复此楼
19楼2015-03-19 08:09:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)
引用回帖:
19楼: Originally posted by anxiaoxin2 at 2015-03-19 08:09:30
版主，早上好啊。。。以上代码是怎么运行的？具体操作步骤能说下吗？我这还是不会运行。。。...

早上好。
代码框里面的代码完全复制，然后在MATLAB中 新建一个m文件 把m文件自带内容完全删除 粘贴代码 运行即可。
赞一下 回复此楼
MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。
20楼2015-03-19 08:36:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 anxiaoxin2 的主题更新
222/3返回列表上一页123下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 283求调剂 +18 那个噜子 2026-04-09 18/900 2026-04-09 18:57 by haironglove
[考研] 一志愿中科大070300化学，314分求调剂 +11 wakeluofu 2026-04-09 11/550 2026-04-09 18:03 by lijunpoly
[考研] 297求调剂 +5 Kwgyz 2026-04-09 5/250 2026-04-09 17:18 by wp06
[考研] 285求调剂 +7 恶法大二的气味� 2026-04-05 10/500 2026-04-08 14:34 by zhq0425
[考研] 调剂求助（生物与医药） +6 @6952 2026-04-06 6/300 2026-04-07 23:52 by lys0704
[考研] 081700学硕，323分，一志愿中国海洋大学求调剂学校 +19 披星河 2026-04-04 19/950 2026-04-07 15:00 by 上岸快快
[考研] 334分机械专硕求调剂 +3 蛋花紫菜汤 2026-04-03 3/150 2026-04-07 14:49 by 逍遥cocoa
[考研] 081200-11408-367学硕求调剂 +4 1_2_3111 2026-04-06 4/200 2026-04-07 08:13 by jp9609
[考研] 277工科求调剂 +10 1915668 2026-04-04 11/550 2026-04-06 23:53 by 心媛意北
[考研] 生物学求调剂 一志愿沪9，326分 +6 刘墨墨 2026-04-06 6/300 2026-04-06 19:36 by lijunpoly
[考研] 285求调剂 +5 mapmath 2026-04-06 6/300 2026-04-06 17:18 by 蓝云思雨
[考研] 一志愿 江南大学 085602 化工专硕 338分求调剂 +15 路痴小琪 2026-04-05 15/750 2026-04-06 09:27 by cql1109
[考研] 322求调剂 +3 嗯哼哼恒 2026-04-05 3/150 2026-04-05 19:52 by nepu_uu
[考研] 0860 求调剂 一志愿国科大 348 分 +3 WiiiP 2026-04-03 3/150 2026-04-05 17:43 by Ecowxq666！
[考研] 一志愿江南大学085501机械工程专硕326分，本科佳木斯大学 +5 顾若浮生 2026-04-03 9/450 2026-04-05 09:57 by 1753564080
[考研] 282电子信息0854专硕调剂 +4 202451007219 2026-04-02 6/300 2026-04-04 21:55 by laoshidan
[考研] 材料专业383求调剂 +8 郭阳阳阳成 2026-04-03 8/400 2026-04-04 10:29 by Rednal.
[考研] 281求调剂 +10 aaawhy 2026-04-03 10/500 2026-04-03 21:42 by lbsjt
[考研] 求材料调剂 一志愿南昌大学 328分 +5 yyy..... 2026-04-03 5/250 2026-04-03 13:46 by 百灵童888
[考研] 320求调剂 +3 农业工程与信息� 2026-04-03 3/150 2026-04-03 11:40 by 土木硕士招生
信息提示
关闭
请填处理意见
关闭