线性系统分析-上机实验.docxVIP

线性系统分析-----上机实验报告姓名：戚煜华学号：3120140024实验一编写一个程序，判读一个线性定常系统的能控性，对于完全能控的系统变换成能控标准型，对于不完全能控的系统进行能控性分解。通过算例进行验证。程序：A,b,c为输入的线性定常系统。function [A_2,b_2,c_2]=zuoye_1(A,b,c)n=length(A); %%%确定系统状态变量的维数M=zeros(n,n); %%%能控性矩阵 Mfor i=1:n M(:,i)=A^(i-1)*b;endrank_M=rank(M);if rank_M==n %%若完全能控则转换成能控标准1型 disp(这个线性定常系统状态完全能控); a=poly(A); %%求A的特征多项式系数 for i=1:n X(i:n,i)=a(1:end-i); end Tc1=zeros(n,n); %%变换矩阵Tc1 for i=1:n Tc1(:,i)=A^(n-i)*b; end Tc1=Tc1*X; A_2=inv(Tc1)*A*Tc1; b_2=inv(Tc1)*b; c_2=c*Tc1;else