matlab在线性系统理论中的应用

《matlab在线性系统理论中的应用》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、个人收集整理仅供参考学习MATLAB在线性系统理论中地应用第一章传递函数与状态空间表达式1.1传递函数与状态空间表达式之间地转换用ss命令来建立状态空间模型.对于连续系统，其格式为sys=ss(A,B,C,D),其中a,b,c,d为描述线性连续系统地矩阵.b5E2RGbCAP当sys1是一个用传递函数表示地线性定常系统时，可以用命令sys=ss(sys1)将其转换成为状态空间形式，也可以用命令sys=ss(sys1,’min’)计算出系统sys地最小实现.p1EanqFDPwexample1：系数传递函数到状态空间表达式>>num=[172424];den=[110355024];

2、g=tf(num,den);sys=ss(g)theansweris：a=x1x2x3x4x1-10-4.375-3.125-1.5x28000x30200x40010b=u1x12x20x30x40c=x1x2x3x4y10.50.43750.750.75d=u1y10Continuous-timemodel.12/12个人收集整理仅供参考学习example2：由传递函数系数，将离散系统脉冲传递函数模型转换成状态空间表达式>>num=[0.310.570.380.89];den=[13.233.982.220.47];DXDiTa9E3dgyu=tf(num,den,'ts',0

3、.1)theansweris:Transferfunction:0.31z^3+0.57z^2+0.38z+0.89-----------------------------------------z^4+3.23z^3+2.98z^2+2.22z+0.47Samplingtime:0.1Pzmap(gyu)%绘制零极点分布图sys=ss(gyu)%将离散系统脉冲传递函数模型转换成状态空间表达式.Theansweris:a=x1x2x3x4x1-3.23-1.49-1.11-0.235x22000x30100x40010b=u1x11x20x30x40c=12/12个人收集整理仅供

4、参考学习x1x2x3x4y10.310.2850.190.445d=u1y10Samplingtime:0.1Discrete-timemodel.Example3：用s求逆矩阵法从系统矩阵a,b,c,d求得传递函数>>symss;a=[01;-2-3];b=[10;11];c=[21;11;-2-1];d=[30;00;01];i=[10;01];f=inv(s*i-a)g=simple(simple(c*f*b)+d)Theansweris:f=[(s+3)/(s^2+3*s+2),1/(s^2+3*s+2)][-2/(s^2+3*s+2),s/(s^2+3*s+2)]g=[3

5、/(s+1)+3,1/(s+1)][2/(s+2),1/(s+2)][-3/(s+1),-1/(s+1)+1]Example4eig()指令，求特征根矩阵和特征向量矩阵函数eig()Example5约旦标准型函数jordan()>>a=[010;001;2-54]a=0100012-54>>[q,j]=jordan(a)q=12/12个人收集整理仅供参考学习1-202-2-24-2-4j=200011001Example6从状态转移矩阵到传递函数地转化举例：cleara=[010;001;-6-11-6];b=[0;0;1];c=[111];d=[0];v=ss(a,b,c,d)[

6、num,den]=ss2tf(a,b,c,d);printsys(num,den)[z,p,k]=ss2zp(a,b,c,d);zpk(z,p,k)x0=[2;0;1];figure(1)step(v)figure(2)initial(v,x0)t=0:0.1:60;u=t;figure(3)lsim(v,u,t);%figure(3)第二章状态转移矩阵与状态方程地解Example1Collect函数地作用是合并同类项，ilaplace（）函数地作用地求取laplace反变换，函数det（）地作用是求方阵地行列式.RTCrpUDGiTsymsstx0xtaophiphi0;%声明变

7、量a=[01;-2-3];I=[10;01];e=s*I-a;c=det(e);12/12个人收集整理仅供参考学习d=collect(inv(e));phi0=ilaplace(d)x0=[1;0];x=phi0*x0%公式与关系：sinh是双曲正弦函数.cosh是双曲余弦函数.　　带h地都是双曲函数.　　sinh(x)=(exp(x)-exp(-x))/2.0;　　cosh(x)=(exp(x)+exp(-x))/2.0;　　tanh(x)=sinh(x)/cosh