2.能观测性判定 %获得系统维数 %构造能观测性判别阵Qo Qo=Q0: CA OI end if rank(Qo)=n%如果Qo的秩等于n,则完全能观 disp(the system is observable") else disp(the system is unobservable); end2.能观测性判定 n=length(A); %获得系统维数 Qo=[C]; %构造能观测性判别阵Qo for i=1:n-1 Qo=[Qo;C*A^(i)]; end if rank(Qo)==n %如果Qo的秩等于n,则完全能观 disp('the system is observable'); else disp('the system is unobservable'); end