新手碼的代碼，有很多內置函數都不清楚，歡迎大佬指路！

目錄

一、測試程序

二、能控能觀性判別函數

三、一般控制矩陣轉能觀標準型函數

---

一、測試程序

A=[-2 0 0 00 -5 1 00 0 -5 10 0 0 -5]; B=[2001]; C=[1 1 0 0]; [AA,BB,CC]=ss2obv(A,B,C);

二、能控能觀性判別函數

function[ks,kv]=con_obv_Rank(A,B,C) %% 能控判別 n=size(A,1); Qc=ctrb(A,B); rankQc=rank(Qc); if rankQc==nks=1; elseks=0; end %% 能觀判別 Q0=obsv(A,C); rankQc=rank(Q0); if rankQc==nkv=1; elsekv=0; end end

三、一般控制矩陣轉能觀標準型函數

function[A,B,C]=ss2obv(a,b,c) [~,kg]=con_obv_Rank(a,b,c); if kg==0error('the sys is unobservable')%如果系統不能觀，則函數直接報錯返回elsedisp('the sys is observable') end %% 求P矩陣 n=size(a,1);A=zeros(n,n); p=poly(a); W=zeros(n,n); w=zeros(1,n); for i=1:nw(n-i+1)=p(i); end for i=1:nfor j=1:n-i+1W(i,j)=w(i+j-1);end end Qc=obsv(a,c); iP=W*Qc; %% 求A,B,C C=c*inv(iP); B=zeros(n,1); B(n,:)=1; A=iP*a*inv(iP);

總結

以上是生活随笔為你收集整理的一般控制矩阵转能观标准型的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。