matlab用SS函数求解2元微分方程图像!

你老师的这个EXAMPLE是将解线性微分方程转换为控制系统中的状态方程,并利用matlab中的控制系统工具箱中的系列函数进行求解。我想既然你老师用该工具箱求解,你应该也就学过控制系统的相关理论知识了。ABCD定义了状态方程x'=Ax+Bu;y=Cx+Du,在这个EXAMPLE里就是x1'=x2x2'=-10x1+10x3x3'=x4x4'=x1-11x3y1=x1y2=x3sys=ss(A,B,C,D)就定义了上述状态方程sys它等价于下述微分方程组y1''=-y1+10y2y2''=y1-11y2然后initial(sys,x0,tspan)即得到状态方程x'=Ax;y=Cx在初值条件x(0)=x0时的解曲线y(t),t∈tspan在该EXAMPLE中即是方程x1'=x2x2'=-10x1+10x3x3'=x4x4'=x1-11x3y1=x1y2=x3x1(0)=1x2(0)=0x3(0)=0。05x4(0)=0的解曲线y1和y2,即y1''=-y1+10y2y2''=y1-11y2y1(0)=1y1'(0)=0y2(0)=0。05y2'(0)=0的解曲线y1和y2plot(t,y);即画出了y1和y2所以你的题目可以照样画葫芦,我们以x1''=-20x1+10x2x2''=10x1+20x2response when x1(0)=2, x2(0)=4; dx1(0)/dt=0, dx2(0)/dt=0;为例下面是matlab代码A=[0 1 0 0;-20 0 10 0; 0 0 0 1; 10 0 20 0]; B=[0; 0 ;0 ;0];C=[1 0 0 0; 0 0 1 0];D=0;sys=ss(A,B,C,D); %状态方程t=0:0。01:1;[y,t1,x]=initial(sys,[2 0 4 0],t); %初始条件仿真plot(t,y);