| 查看: 3154 | 回复: 21 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
anxiaoxin2金虫 (小有名气)
|
[求助]
MATLAB应用四阶龙格库塔法求解反应动力学问题? 已有2人参与
|
||
|
各位英雄大侠，你们好。本人在做铁水喷吹脱硫的优化，遇到一个问题，文献中是应用MATLAB四阶龙格库塔法分析的，本人对MATLAB不熟，希望大家能帮帮忙！ 谢谢。。。 具体问题见附件。 |
回复此楼» 本帖附件资源列表
-
欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com - 附件 1 : 脱硫动力学微分方程.docx
2015-03-05 14:00:32, 46.03 K
» 猜你喜欢
反铁磁体中的磁性切换：两种不同的机制已成功可视化
已经有0人回复
-
求标准粉末衍射卡号 ICDD 01-076-1802
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有86人回复
-
新西兰Robinson研究所招收全奖PhD
已经有0人回复
-
石墨烯转移--二氧化硅衬底石墨烯
已经有0人回复
-
笼目材料中量子自旋液体基态的证据
已经有0人回复
-
数学教学论硕士可以读数学物理博士吗？
已经有0人回复
-
德国亥姆霍兹Hereon中心汉堡分部招镁合金腐蚀裂变SCC课题方向2026公派博士生
已经有4人回复
-
澳门大学 应用物理及材料工程研究院 潘晖教授课题组诚招博士后
已经有11人回复
anxiaoxin2
金虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 2572.8
- 帖子: 178
- 在线: 86小时
- 虫号: 1042931
- 注册: 2010-06-17
- 性别: GG
- 专业: 无机非金属类光电信息与功
送红花一朵 |
版主，这个还是运行不了啊。。。 这是代码，你帮我看看呗 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
CelestialCYJ
木虫 (小有名气)
- 应助: 33 (小学生)
- 金币: 2482.6
- 红花: 12
- 帖子: 144
- 在线: 43.3小时
- 虫号: 3277506
- 注册: 2014-06-16
- 专业: 结构工程
2楼2015-03-05 16:11:26
anxiaoxin2
金虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 2572.8
- 帖子: 178
- 在线: 86小时
- 虫号: 1042931
- 注册: 2010-06-17
- 性别: GG
- 专业: 无机非金属类光电信息与功
3楼2015-03-05 17:00:56
月只蓝
主管区长 (职业作家)
-
专家经验: +1059 - 计算强帖: 8
- 应助: 1712 (讲师)
- 贵宾: 8.888
- 金币: 68123.7
- 散金: 1938
- 红花: 443
- 沙发: 4
- 帖子: 4373
- 在线: 3291.4小时
- 虫号: 1122189
- 注册: 2010-10-14
- 专业: 宇宙学
- 管辖: 计算模拟区
4楼2015-03-05 21:00:42