急急急！matlab大神请赐教解微分方程组

显示全部楼层 · 2010-5-27 22:17:49

小弟对matlab编程不甚了解，请各位不吝赐教！感激不尽！
x=[x(1);x(2);x(3);x(4);x(5)]
满足
dx(1)/dt=a1+b1*(x(2)-x(1))
dx(2)/dt=c1*(x(1)-x(2))+c2*(x(3)-x(2))
dx(3)/dt=a2+b2*(x(2)+x(4)-2*x(3))
dx(4)/dt=c1*(x(5)-x(4))+c2*(x(3)-x(4))
dx(5)/dt=a1+b1*(x(4)-x(5))
其中
a1=-22.5000;
a2=18
b1=0.1588
b2=0.0623
c1=7.0285
c2=2.8463
要求输出x与时间t的关系曲线。

千问 · 2010-5-27 22:17:49

function dxdt=weifenequation(t,x)a1=-22.5000;a2=18;b1=0.1588;b2=0.0623;c1=7.0285;c2=2.8463;dxdt=zeros(5,1);dxdt(1)=a1+b1*(x(2)-x(1));dxdt(2)=c1*(x(1)-x(2))+c2*(x(3)-x(2));dxdt(3)=a2+b2*(x(2)+x(4)-2*x(3));dxdt(4)=c1*(x(5)-x(4))+c2*(x(3)-x(4));dxdt(5)=a1+b1*(x(4)-x(5));%以上保存在weifenequation.m文件[t,x]=ode4

千问 · 2010-5-27 22:17:49

要学会自己写

急急急！matlab大神请赐教 解微分方程组

急急急！matlab大神请赐教解微分方程组