![]() |
急!!!新手求助,请各位高手帮帮忙
a=0.35;A(1)=-29;A(2)=138;A(3)=129;alpher(1)=0.03;alpher(2)=0.13;alpher(3)=0.625;b=0.5;beta=33;
dlambdadt=1.6380*(1-exp(-beta*t))./b; Y2=dsolve('DQ1=A(1)*1.6380*(1-exp(-beta*t))/b-alpher(1)*Q1','Q1(0)=1'); Y3=dsolve('DQ2=A(2)*1.6380*(1-exp(-beta*t))/b-alpher(2)*Q2','Q2(0)=1'); %Q2(T0)=0 Y4=dsolve('DQ3=A(3)*1.6380*(1-exp(-beta*t))/b-alpher(3)*Q3','Q3(0)=1'); Q=((819/500*A(1)/b*(1/alpher(1)*exp(alpher(1)*t)-1/(alpher(1)-beta)*exp(t*(alpher(1)-beta)))+1-819/500*A(1)/b*(1/alpher(1)-1/(alpher(1)-beta)))*exp(-alpher(1)*t))+((-819/500*A(2)/b*(-1/alpher(2)*exp(alpher(2)*t)+1/(alpher(2)-beta)*exp(t*(alpher(2)-beta)))+1+819/500*A(2)/b*(-1/alpher(2)+1/(alpher(2)-beta)))*exp(-alpher(2)*t))+((A(3)/b*(819/500/alpher(3)*exp(alpher(3)*t)-819/500/(alpher(3)-beta)*exp(t*(alpher(3)-beta)))+1-A(3)/b*(819/500/alpher(3)-819/500/(alpher(3)-beta)))*exp(-alpher(3)*t)) ??? Error using ==> mtimes Inner matrix dimensions must agree. QQ=subs(Q); % 在不同t时刻Q的值,赋值给QQ(向量) for q=QQ if q<0 T=T0*(a*q+1)/(1-q); else T=T0*(1+(2+a)*q/(1-q)); end end t=1:max(size(T)); plot(t,T) T0在其它程序已经算出,为一列数据向量~请各位高手帮帮忙,非常感谢~~~~~~ |
回复: 急!!!新手求助,请各位高手帮帮忙
其中Q=Y2+Y3+Y4,因为我把上式的Y2,Y3,Y4微分方程解出,Y2,Y3,Y4为函数式,直接代入Q,所以Q式特别长~请各位高手帮帮忙,非常感谢~
|
所有时间均为北京时间。现在的时间是 06:24。 |
Powered by vBulletin
版权所有 ©2000 - 2025,Jelsoft Enterprises Ltd.