%《MATLAB及在電子信息課程中的應(yīng)用》第八章例8.16程序q816
% 化為可控階梯形
% 電子工業(yè)出版社出版  陳懷琛 吳大正 高西全合著 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;	% 系數(shù)矩陣賦值
s1=ss(A,B,C,D);							            % 構(gòu)成LTI模型
[Abar,Bbar,Cbar,T,k]=ctrbf(s1.a,s1.b,s1.c)	        % 化為可控階梯型
rA=rank(A)						                    % 狀態(tài)方程系數(shù)矩陣的秩
rc=sum(k)								            % 判斷可控矩陣的秩
[Abaro,Bbaro,Cbaro,To,ko]=obsvf(s1.a,s1.b,s1.c);	% 化為可觀階梯型
ro=sum(ko)								            % 判斷可觀矩陣的秩