| 查看: 3349 | 回复: 21 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
anxiaoxin2金虫 (小有名气)
|
[求助]
MATLAB应用四阶龙格库塔法求解反应动力学问题? 已有2人参与
|
||
|
各位英雄大侠，你们好。本人在做铁水喷吹脱硫的优化，遇到一个问题，文献中是应用MATLAB四阶龙格库塔法分析的，本人对MATLAB不熟，希望大家能帮帮忙！ 谢谢。。。 具体问题见附件。 |
回复此楼» 本帖附件资源列表
-
欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com - 附件 1 : 脱硫动力学微分方程.docx
2015-03-05 14:00:32, 46.03 K
» 猜你喜欢
B区学生调剂-兰州交通大学材料科学与工程学院
已经有12人回复
-
山西大同大学物理学专业还有调剂名额，欢迎调剂！
已经有13人回复
-
物理学I论文润色/翻译怎么收费?
已经有211人回复
-
桂林理工大学物理学专业招收调剂
已经有18人回复
-
VASP 的一组 GPU / CPU 基准测试记录
已经有0人回复
-
津理工大学晶体材料全国重点实验室刘红军教授课题组招收博士生一名
已经有0人回复
-
【原创讨论】从电子约束到物质编辑：一套可迭代的环形磁场科技树
已经有0人回复
-
【方案分享】单环磁场+轴心控制+偏转导出电子束约束系统（可行性实验）
已经有6人回复
-
【修正版】单环用磁约束低速电子实验方案（简化版）
已经有0人回复
-
桂林理工大学物理学专业招收调剂，还有三个名额！！！
已经有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
- 金币: 68132.7
- 散金: 1938
- 红花: 443
- 沙发: 4
- 帖子: 4373
- 在线: 3291.4小时
- 虫号: 1122189
- 注册: 2010-10-14
- 专业: 宇宙学
- 管辖: 计算模拟区
4楼2015-03-05 21:00:42
20