24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

zhang66163

银虫 (正式写手)

应助: 0 (幼儿园)
金币: 212.4
散金: 120
帖子: 322
在线: 322小时
虫号: 684233
注册: 2008-12-29
性别: GG
专业: 催化化学

[求助] 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]

回复此楼

» 猜你喜欢

求标准粉末衍射卡号 ICDD 01-076-1802 已经有0人回复
新西兰Robinson研究所招收全奖PhD 已经有0人回复
物理学I论文润色/翻译怎么收费? 已经有152人回复
石墨烯转移--二氧化硅衬底石墨烯已经有0人回复
笼目材料中量子自旋液体基态的证据已经有0人回复
数学教学论硕士可以读数学物理博士吗？已经有0人回复
德国亥姆霍兹Hereon中心汉堡分部招镁合金腐蚀裂变SCC课题方向2026公派博士生已经有4人回复
澳门大学应用物理及材料工程研究院潘晖教授课题组诚招博士后已经有11人回复
求助NH4V4O10晶体的CIF文件已经有0人回复
英国全奖博士招聘-深度学习与量子物理已经有0人回复
间接带隙半导体有效质量求助已经有0人回复

1楼 2015-08-10 21:31:35

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)

专家经验: +1059
计算强帖: 8
应助: 1712 (讲师)
贵宾: 8.888
金币: 68123.7
散金: 1938
红花: 443
沙发: 4
帖子: 4373
在线: 3291.4小时
虫号: 1122189
注册: 2010-10-14
专业: 宇宙学
管辖: 计算模拟区

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
zhang66163: 金币+40, ★★★很有帮助 2015-08-11 11:02:11

少一个初值条件，由C0=[1, 0, 0, 0]改为：C0=[1, 0, 0, 0,0]，即默认组分E初始浓度为0；如果初值条件不是这样的，请自行需要改，须填满5个，因为有5个方程。

CODE:

function Batchreactor

clear all;clc

global T R

T=224.6+273;

R=8.314;

C0=[1, 0, 0, 0,0];

tspan=[0, 1E+4];

[t,C]=ode45(@Massequations,tspan,C0);

[t C]

figure(1)

plot(t,C(:,1),'r-',t,C(:,2),'O-',t,C(:,3),'b-',t,C(:,4),'g-')

function dCdt=Massequations(t,C)

global T R

k0=[5.78052E+10, 3.92317E+12, 1.64254E+4, 6.264E+8];

Ea=[124670, 150386, 77954, 111528];

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);

dCdt=[rA;rB;rC;rD;rE];