Matlab 常用指令整理

1.由傳遞函數到狀態空間

>> num=[10];
>> den=[1 0.1 3];
>> [A,B,C,D]=tf2ss(num,den)

A =

   -0.1000   -3.0000
    1.0000         0


B =

     1
     0


C =

     0    10


D =

     0

傳遞函數爲:

x(k+1)=\begin{vmatrix} -0.1&-3 \\1& 0\end{vmatrix}x(k)+\begin{bmatrix} 1\0 \end{bmatrix}u(k)

y=\begin{bmatrix} 0\\ 10 \end{bmatrix}x

2.由狀態空間到傳遞函數

[num den]=ss2tf(A,B,C,D)

num =

     0     0    10


den =

    1.0000    0.1000    3.0000

與上面的傳遞函數一致

3.連續空間的離散化

 狀態方程爲:(˙x)=Ax+Bu

[G,H]=c2d(A,B,Ts)

其中x(k+1)=Gx(k)+Hu(k)

仍舊以上面的狀態方程爲例

Ts=0.1;   //採樣時間爲0.1s

[G,H]=c2d(A,B,Ts)

G =

    0.9752   -0.2970
    0.0990    0.9851


H =

    0.9900
    0.0497

則離散的傳遞函數爲:
x(k+1)=\begin{vmatrix} 0.9752&-0.2970 \\0.0990 & 0.9851 \end{vmatrix}x(k)+\begin{bmatrix} 0.9900\\0.0497 \end{bmatrix}u(k)

4.求矩陣特徵值

已知矩陣A,

[x,y]=eig(A);

返回之中,x爲特徵向量組成的矩陣,y的對角元素代表特徵值。

5.求矩陣的能觀性與能觀性

CONT=ctrb(A,B)    //能控性

rank(CONT);

CONT爲觀測矩陣,返回的秩爲滿秩則可觀

CONT=obsv(A,B)    //能觀性

rank(CONT);

 

發表評論
所有評論
還沒有人評論,想成為第一個評論的人麼? 請在上方評論欄輸入並且點擊發布.
相關文章