资源描述：
《连续系统状态方程的数值解》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、实验26连续系统状态方程的数值解1学习使用Matlab的各种系统模型转换函数，加深对系统模型几种形式的理解。学习用Matlab计算连续系统状态方程的数值方法。加深对连续系统状态方程的的理解，对系统零输入响应、零状态响应的理解。实验目的2实验原理与说明（一）系统模型的相互转换线性非时变系统的系统模型有:状态空间型系统函数的多项式型系统函数的零极点型极点留数型它们都能描述系统的特性，但各有不同的应用场合。对于线性非时变系统，这几种模型是可以互相转换的。用Matlab就可以实现这一转换。3实验原理与说明1．状态空间型
2、与系统函数的多项式型互相转换Matlab提供的函数：[b,a]=ss2tf(A,B,C,D)将状态空间型转换成的多项式型，其中，b,a为的分子、分母多项式系统，A,B,C,D为状态空间型的系数矩阵。[A,B,C,D]=tf2ss(b,a)表示将的多项式型转换成状态空间型。2．状态空间型与系统函数的零极点型互相转换Matlab提供的函数：[z,p,k]=ss2zp(A,B,C,D)将状态空间型转换成的零极点型，其中，z,p,k为的零点、极点、增益，A,B,C,D为状态空间型的系数矩阵。[A,B,C,D]=zp2s
3、s(z,p,k)表示将的零极点型转换成状态空间型。4实验原理与说明3．系统函数的零极点型与多项式型互相转换将多项式型转换成零极点型：[z,p,k]=tf2zp(b,a)将的零极点型转换成多项式型：[b,a]=zp2tf(z,p,k)4．系统函数的极点留数型与多项式型互相转换将多项式型转换成极点留数型：[r,p,k]=residue(b,a)将的零极点型转换成多项式型：[b,a]=residue(r,p,k)（二）用lsim求系统响应的数值解在前面曾用到过这个函数，它的功能特别强，能对系统函数模型和状态空间模型对
4、线性非时变系统仿真，对状态空间模型可以求系统全响应、零输入响应、零状态响应的数值解。5计算示例11、已知描述系统的微分方程为求出它的四种模型。解用Matlab计算的程序如下：%系统模型相互转换的程序exp26_1.mformatcompactb=input('系统函数分子系数数组b=');a=input('系统函数分母系数数组a=');printsys(b,a,'s')disp('零极点型模型')[z,p,k]=tf2zp(b,a)disp('极点留数型模型')[r,p,k]=residue(b,a)disp(
5、'状态空间型模型')[A,B,C,D]=tf2ss(b,a)运行程序后，结果显示如下。系统函数分子系数数组b=[2-53]系统函数分母系数数组a=[2359]num/den=2s^2-5s+3-----------------------2s^3+3s^2+5s+9零极点型模型z=1.50001.0000p=-1.64410.0721+1.6528i0.0721-1.6528ik=1极点留数型模型r=-0.2322+0.4716i-0.2322-0.4716i1.4644p=0.0721+1.6528i0.07
6、21-1.6528i-1.6441k=[]状态空间型模型A=-1.5000-2.5000-4.50001.00000001.00000B=100C=1.0000-2.50001.5000D=06计算示例2设某系统的状态方程和输出方程为系统的初始状态为，输入信号，试求状态变量和输出的零输入响应、零状态响应和全响应。解：用Matlab计算的程序如下：7计算示例2%计算状态方程和输出方程的数值解exp26_2.mt=0:0.01:3;A=[10;1-3];B=[10]';C=[-0.251];D=[0];zi=[12
7、];%初始条件f=15*sin(2*pi*t);%输入信号sys=ss(A,B,C,D)[y,t,x]=lsim(sys,f,t,zi)%计算全响应f=zeros(1,length(t));%令输入为零yzi=lsim(sys,f,t,zi);%计算零输入响应f=15*sin(2*pi*t);zi=[00];%令初始条件为零yzs=lsim(sys,f,t,zi);%计算零状态响应8计算示例2figure(1)plot(t,x(:,1),'-',t,x(:,2),'-.','linewidth',2)legen
8、d('x(1)','x(2)')%显示图例title('状态变量波形')xlabel('t(sec)')figure(2)plot(t,y,'-',t,yzi,'-.',t,yzs,':','linewidth',2)legend('y','yzi','yzs')%显示图例title('系统响应,零输入响应,零状态响应')xlabel('t(sec)')运行程序后，系统全响应、零输入响
