q816.m

matlab在电子信息中的应用和实例中的源代码

%《MATLAB及在电子信息课程中的应用》第八章例8.16程序q816
% 化为可控阶梯形
% 电子工业出版社出版  陈怀琛 吴大正 高西全合著 2001年10月初版，2003年7月第二版
%
clear,
A=[-2,2,-1;0,-2,0;1,-4,3];B=[0;0;1];C=[1,-1,1];D=0;	% 系数矩阵赋值
s1=ss(A,B,C,D);							            % 构成LTI模型
[Abar,Bbar,Cbar,T,k]=ctrbf(s1.a,s1.b,s1.c)	        % 化为可控阶梯型
rA=rank(A)						                    % 状态方程系数矩阵的秩
rc=sum(k)								            % 判断可控矩阵的秩
[Abaro,Bbaro,Cbaro,To,ko]=obsvf(s1.a,s1.b,s1.c);	% 化为可观阶梯型
ro=sum(ko)								            % 判断可观矩阵的秩