| 查看: 649 | 回复: 2 | |||
zhang66163银虫 (正式写手)
|
[求助]
matlab程序编写 已有2人参与
|
|
按照教程编写了一个程序：如下 function Batchreactor clear all;clc T=224.6+273 R=8.314 k0=[5.78052E+10, 3.92317E+12, 1.64254E+4, 6.264E+8] Ea=[124670, 150386, 77954, 111528] C0=[1, 0, 0, 0] tspan=[0, 1E+4] [t,C]=ode45(@Massequations,tspan,C0,[],k0,Ea,T,R); plot(t,C(:,1),'r-',t,C(:,2),'O-',t,C(:,3),'b-',t,C(:,4),'g-'); function dCdt=Massequations(t,C,k0,Ea,T,R); k=k0*exp(-Ea/(R*T)); k(5)=2.1667E-4; rA=-(k(1)+k(2))*C(1); rB=k(1)*C(1)-k(3)*C(2); rC=k(2)*C(1)-k(4)*C(3); rD=k(3)*C(2)-k(5)*C(4); rE=k(4)*C(3)+k(5)*C(4); 可是运行后，出现了错误；哪位大侠可以帮忙解释一下。谢谢 Error using * Inner matrix dimensions must agree. Error in concentration>Massequations (line 16) k=k0*exp(-Ea/(R*T)); k(5)=2.1667E-4; Error in odearguments (line 88) f0 = feval(ode,t0,y0,args{:}); % ODE15I sets args{1} to yp0. Error in ode45 (line 114) [neq, tspan, ntspan, next, t0, tfinal, tdir, y0, f0, odeArgs, odeFcn, ... Error in concentration (line 10) [t,C]=ode45(@Massequations,tspan,C0,[],k0,Ea,T,R); dCdt=[rA;rB;rC;rD;rE] |
回复此楼» 猜你喜欢
求助中医院临床中医中标的青年标书学习，谢谢
已经有6人回复
-
浙大物院表面物理课题组长期招收普博及直博学生
已经有259人回复
-
物理学I论文润色/翻译怎么收费?
已经有75人回复
-
【征稿】Springer Nature热电专题 Discover Materials / Discover Chemistry
已经有47人回复
-
25届硕士申请，材料计算方向
已经有2人回复
-
课题组招1-2名硕士
已经有1人回复
-
机时越来越贵之后，我把一部分 DFT 流程搬到了 AI 计算平台上
已经有0人回复
-
斯坦福大学与Elsevier世界前2%学者新排名查询
已经有1人回复
-
27申博推荐：南京航空航天大学国际前沿院光学方向招收博士生！
已经有19人回复
-
华南师范大学+双一流学科物理学+广东省杰青团队2027年博士招生+机器学习驱动材料设计
已经有1人回复
-
apl投稿状态
已经有2人回复
2楼2015-08-11 09:08:23
月只蓝
主管区长 (职业作家)
-
专家经验: +1059 - 计算强帖: 8
- 应助: 1712 (讲师)
- 贵宾: 8.888
- 金币: 68181.7
- 散金: 1938
- 红花: 443
- 沙发: 4
- 帖子: 4373
- 在线: 3292.4小时
- 虫号: 1122189
- 注册: 2010-10-14
- 专业: 宇宙学
- 管辖: 计算模拟区
3楼2015-08-11 09:27:15