| 查看: 570 | 回复: 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] |
回复此楼» 猜你喜欢
基元I理论下三大核心空间现象精准推导与细节解析
已经有0人回复
-
基于基元 I 统一理论的反重力理论推导
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有58人回复
-
基于基元I统一理论的量子力学本源推导
已经有1人回复
-
推荐一款可以AI辅助写作的Latex编辑器SmartLatexEditor，超级好用，AI润色，全免费
已经有20人回复
-
【EI|Scopus 双检索】第六届智能机器人系统国际会议（ISoIRS 2026）
已经有1人回复
-
2026年第四届电动车与车辆工程国际会议（CEVVE 2026）
已经有0人回复
2楼2015-08-11 09:08:23
月只蓝
主管区长 (职业作家)
-
专家经验: +1059 - 计算强帖: 8
- 应助: 1712 (讲师)
- 贵宾: 8.888
- 金币: 68126.7
- 散金: 1938
- 红花: 443
- 沙发: 4
- 帖子: 4373
- 在线: 3291.4小时
- 虫号: 1122189
- 注册: 2010-10-14
- 专业: 宇宙学
- 管辖: 计算模拟区
3楼2015-08-11 09:27:15