q816.m

来自「dsp相关的matlab程序」· M 代码 · 共 12 行

%《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)								% 判断可观矩阵的秩