解一个矩阵
1视图（过去30天）
显示旧的评论
pnh3=.046；%%托尔
po2=0.068；%%托尔
R1=1.987；%%cal/gmol*K
R=62.364；%%托尔升/克摩尔
vo=100/3600；%%信用证
P=.115；%%托尔
Per=40*300；%%微米
zf=1；%%厘米
vnh3=1；
vo2=5/2；
zspan=[01]；%%厘米
T=零（3,1）；
T（1）=473；
T（2）=673；
T（3）=800；
fio=Ljallo_f2_q1（T，R，vo，pnh3，po2）；
ic1=fio（1，：）；
ic2=fio（2，：）；
ic3=fio（3，：）；
符号Y
DFDT=Ljallo_f1_q1（T，R1，Per，vnh3，vo2，y，vo，R）；
diffy1=DFDT（1，：）；
diffy2=DFDT（2，：）；
diffy3=DFDT（3，：）；
f1=@（z，y）diffy1；
f2=@（z，y）diffy2；
f3=@（z，y）diffy3；
[z，y]=ode45（@（z，y）f1，zspan，ic1）；
[z，y]=ode45（@（z，y）f2，zspan，ic2）；
[z，y]=ode45（@（z，y）f3，zspan，ic3）；
作用dfdt=Ljallo_f1_q1（T，R1，Per，vnh3，vo2，y，vo，R）
A=零（3,1）；
B=零（3,1）；
C=零（3,1）；
rn=零（3,1）；
dfdt=零（3,2）；
对于i=1:3
A（i）=3.4*10^-8*exp（21700./（R1.*T（i））*（y（1）。/vo*R*T（i））*（y（2）。/vo*R*T（i））。^.5；
B（i）=1+8*10^-2*exp（4400/（R1.*T（i）））*（y（2）/vo*R*T（i））^.5；
C（i）=1+1.6*10^-3*exp（25500/（R1.*T（i））*（y（1）/vo*T（i））；
rn（i）=A（i）/（B（i）*C（i））；
dfdt（i，1）=rn（i）。*每*vnh3；
dfdt（i，2）=rn（i）。*每*vo2；
终止
作用iFlows=Ljallo_f2_q1（T、R、vo、pnh3、po2）
cnh3=零（3,1）；
finh3o=零（3,1）；
co2=零（3,2）；
fio2o=零（3,2）；
iFlows=零（3,2）；
对于i=1:3
cnh3（i，1）=pnh3/（R.*T（i））；
finh3o（i，1）=cnh3（i，1）。*vo；
iFlows（i，1）=finh3o（i，1）；
co2（i，2）=po2/（R.*T（i））；
fio2o（i，2）=co2（i，2）。*vo；
iFlows（i，2）=fio2o（i，2）；
终止
>>Ljallo_s_q1
指数超过数组元素数（1）。
错误在sym/subsref中（第890行）
R_tilde=内置(“subsref”，L_tilde，Idx）；
错误在Ljallo_f1_q1（第8行）
A（i）=3.4*10^-8*exp（21700./（R1.*T（i））*（y（1）。/vo*R*T（i））*（y（2）。/vo*R*T（i））。^.5；
错误在Ljallo_s_q1（第23行）
DFDT=Ljallo_f1_q1（T，R1，Per，vnh3，vo2，y，vo，R）；
在解决此问题时遇到问题，请使用矩阵3x2。
我把方程和初始条件分解成3个1x2矩阵。
我想同时解每一行。
0评论
答案（0）
您还可以从以下列表中选择网站：