线性系统结构与分解及标准型_new重点解读.doc

实验题目：线性系统结构分析与分解及标准型 【实验时间】 【实验地点】 【实验目的】 掌握线性系统状态空间标准型、解及其模型转换 【实验设备与软件】 MATLAB 【实验原理】 标准型变换、矩阵Jordan型变换、特征值 (1)标准型变换 命令格式 csys=canon(sys,’type’) (2) 矩阵Jordan规范型 命令格式 [V J]=Jordan(A) (3)求矩阵特征值和特征向量 命令格式 [V J]=eig(A) cv=eig(A) 2、状态模型的相似变换: 命令格式 sysb=ss2ss(sys,T) 传递函数模型与状态空间模型之间的相互转换: 命令格式 [A,B,C,D]=tf2ss(num,den) [num,den]=ss2tf(A,B,C,D,iu) zpk模型与空间状态模型之间的相互转换： 命令格式 [A,B,C,D]=zp2ss(z,p,k) [z,p,k]=ss2zp(A,B,C,D,iu) 线性定常系统的可控性与可观性及结构分解 (1)状态可控性 Qc=ctrb(A,B) r=rank(Qc) l=length(A) (2)状态可观性 Qo=obsv(A,C) r=rank(Qo) l=size(A,1) (3)输出可控性 Co=ctrb(A,B) m=size(C,1) Qyc=[C*Co,D] Tm=rank(Qyc) 4、定常线性系统的标准型（转换限于SISO系统） Jordan标准型 能控标准型，限于SISO系统 能控标准型，限于SISO系统 【实验内容、方法、过程与分析】 1、已知线性系统 在MATLAB中输入该系统的状态空间描述： 判断其状态可控性、可观性和传递函数的关系，并加以说明分析。 ?A=[-6?-0.625?0.75;8?0?0;0?2?0];B=[1;0;0];C=[1?-0.25?0.0625];??Qc=ctrb(A,B)Qc?=?????1????-6????31?????0?????8???-48?????0?????0????16?Rc=rank(Qc)Rc?=?????3?Co=ctrb(A,B);?m=size(C,1);?Qyc=[C*Co,1];?Tm=rank(Qyc)Tm?=?????1??m=size(C,1)m?=?????1??Qo=obsv(A,C)Qo?=????1.0000???-0.2500????0.0625???-8.0000???-0.5000????0.7500???44.0000????6.5000???-6.0000??Ro=rank(Qo)Ro?=?????2 可知系统是完全能控但不完全能观的 对系统分别按能控性分解、能观性分解以及能控能观性分解，因为系统能控所以不用分解了。 按能观性分解： [Abar ,Bbar Cbar,T,K]=obsvf(A,B,C) Abar = 1.0000 -2.2772 2.7134 -0.0000 0.3407 -5.5238 0.0000 2.6251 -7.3407 Bbar = 0.0558 -0.2432 0.9684 Cbar = 0.0000 -0.0000 1.0327 T = 0.0558 0.4465 0.8930 -0.2432 -0.8614 0.4459 0.9684 -0.2421 0.0605 K = 1 1 0 2、在Matlab中建立并运行如下的.m代码，回答下面的问题。 num=[1 2 3]; den=conv([1 6 25],[1 12 35]); G=tf(num,den) Gs=ss(G) [V J]=Jordan(Gs.a) %求特征向量和Gs.a的Jordan标准型 Gss=ss2ss(Gs,inv(V)) %Jordan型系统 Gsm=canon(Gs,’model’) %模态型系统 Gsf=canon(Gs,’companion’)