matlab 微分方程组

来源:百度知道 编辑:UC知道 时间:2024/06/19 09:14:09
dx/dt=
方程为:
dx/dt=A0*A1*s*x-A2*x
N*(ds/dt)=-[N*A0*A1*s+A4*s*(A5+A1)]
其中A0、A1、A2、A3、A4、A5、N都是常数,

用MATLAB实现求解功能,并且实现画图。请高手指点一下呀!谢谢。

你的常数让人不知所措,好歹给组数据吧
数据自己编的

a0=1;
a1=1;
a2=1;a3=1;a4=1;a5=1;n=1;
key=dsolve('Dx=a0*a1*s*x-a2*x','n*Ds=-n*a0*a1*s-a4*s*(a5+a1)');
x=key.x
s=key.s
t=0:0.01:1;
C1=3;C2=4;
x=eval(x);
s=eval(s);

plot(t,x,'r--',t,s,'b:')

祝君好运哈