新手碼的代碼，有很多內(nèi)置函數(shù)都不清楚，歡迎大佬指路！

目錄

一、測試程序

二、能控性判別函數(shù)

三、一般控制矩陣轉(zhuǎn)能控標(biāo)準(zhǔn)型函數(shù)

---

一、測試程序

A=[0 1 3 42 3 1 53 4 5 14 5 3 7]; C=[1 2 3 4]; B=[1;2;3;4]; [AA,BB,CC]=ss2con(A,B,C);

二、能控性判別函數(shù)

function[ks,kv]=controllability_RankCriterion(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

三、一般控制矩陣轉(zhuǎn)能控標(biāo)準(zhǔn)型函數(shù)

function[A,B,C]=ss2con(a,b,c) [ks,~]=controllability_RankCriterion(a,b,c); if ks==0error('the sys is uncontrollable') elsedisp('the sys is controllable') %% 求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=ctrb(a,b); P=Qc*W; %% 求A,B,C C=c*P; B=zeros(n,1); B(n,1)=1; den=zeros(1,n); for i=1:nden(i)=p(i+1); end for i=1:n-1A(i,i+1)=1; end dxn=fliplr(den); dxn=-1*dxn; A(n,:)=dxn; end

創(chuàng)作挑戰(zhàn)賽新人創(chuàng)作獎勵來咯，堅持創(chuàng)作打卡瓜分現(xiàn)金大獎

總結(jié)

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

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯，歡迎將生活随笔推薦給好友。