%《MATLAB在電子信息課程中的應用》第八章例8.16程序q816
% 化為可控階梯形
% 電子工業出版社出版  陳懷琛 吳大正 高西全合著 2001年10月 

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