Labfans是一个针对大学生、工程师和科研工作者的技术社区。 论坛首页 | 联系我们(Contact Us)
MATLAB爱好者论坛-LabFans.com
返回   MATLAB爱好者论坛-LabFans.com > 工程数学软件 > MATLAB论坛
MATLAB论坛 一切MATLAB相关问题在此讨论。
回复
 
主题工具 显示模式
旧 2009-04-27, 20:46   #1
M新手
初级会员
 
注册日期: 2009-04-27
帖子: 2
声望力: 0
M新手 正向着好的方向发展
默认 急!!!新手求助,请各位高手帮帮忙

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在其它程序已经算出,为一列数据向量~请各位高手帮帮忙,非常感谢~~~~~~
M新手 当前离线   回复时引用此帖
旧 2009-04-28, 12:44   #2
M新手
初级会员
 
注册日期: 2009-04-27
帖子: 2
声望力: 0
M新手 正向着好的方向发展
默认 回复: 急!!!新手求助,请各位高手帮帮忙

其中Q=Y2+Y3+Y4,因为我把上式的Y2,Y3,Y4微分方程解出,Y2,Y3,Y4为函数式,直接代入Q,所以Q式特别长~请各位高手帮帮忙,非常感谢~
M新手 当前离线   回复时引用此帖
回复


发帖规则
不可以发表新主题
不可以发表回复
不可以上传附件
不可以编辑自己的帖子

启用 BB 代码
论坛启用 表情符号
论坛启用 [IMG] 代码
论坛禁用 HTML 代码



所有时间均为北京时间。现在的时间是 17:00


Powered by vBulletin
版权所有 ©2000 - 2025,Jelsoft Enterprises Ltd.