%《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)								            % 判斷可觀矩陣的秩