ode15s解时出错：Index exceeds matrix dimensions

amygod · 发表于 2006-6-6 19:16

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

您需要登录才可以下载或查看，没有账号？我要加入

x

function dydt=dyncstr(t,y) global T0 UA Vr=10; v0=1000; CA0=0.3; CB0=0; Ta=273.15+60; CpA=200; CpB=CpA; CpC=CpA; E1=9900; E2=27000; R=1.987; DHr1=-55000; DHr2=-71500; T=y(1); CA=y(2); CB=y(3); k1=3.3*exp((E1/R)*(1/300.-1/T)); k2=4.58*exp((E2/R)*(1/500.-1/T)); dTdt=(UA*(Ta-T)+v0*CA0*CpA*(T0-T)-(Vr*(k1*CA)*DHr1+Vr*(k2*CB)*DHr2))/Vr/CA0/CpA; dCAdt=(v0*CA0-v0*CA-(k1*CA)*Vr)/Vr; dCBdt=(v0*CB0-v0*CB-(k2*CB-k1*CA)*Vr)/Vr; dydt(1)=dTdt; dydt(2)=dCAdt; dydt(3)=dCBdt; dydt=dydt'; *********************************************************************** clear; >> global T0 UA >> UA=60000;Tr=678.4069;Ca=8.5938e-004;Cb=0.0051; >> for i=1:200, T0=400-i; y0=[Tr Ca Cb];tspan=[0:0.01:0.49]; [t,y]=ode15s('dyncstr',tspan,y0); Tr=y(50,1);Ca=y(50,2);Cb=y0(50,3); for j=1:50, kk=(i-1)*50+j;tt(kk)=((i-1)*50+j)*0.01; yt(kk)=T0;yy(kk)=y(j,1); end plot(tt,yy,'-',tt,yt,'-') axis([0 100 100 800]);pause(0.01) end

amygod · 发表于 2006-6-6 19:19

??? Index exceeds matrix dimensions.怎么回事？

无水1324 · 发表于 2006-6-6 21:55

<DIV class=quote>以下是引用amygod在2006-6-6 19:16:23的发言： 
y0=[Tr Ca Cb];tspan=[0:0.01:0.49]; [t,y]=ode15s('dyncstr',tspan,y0); Tr=y(50,1);Ca=y(50,2);Cb=y0(50,3); for j=1:50,</DIV>
你上面定义的tspan是50个点， 改为 ............................ Tr=y(49,1);Ca=y(49,2);Cb=y0(49,3); for j=1:49, ................................. 应该就好了 我的matlab有问题，没有试，你试一下先

账号		自动登录	找回密码
密码			我要加入

[近似分析] ode15s解时出错：Index exceeds matrix dimensions

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

回复：（amygod）ode15s解时出错：Index exceeds matr...