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
傳遞函數爲:
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
則離散的傳遞函數爲:
4.求矩陣特徵值
已知矩陣A,
[x,y]=eig(A);
返回之中,x爲特徵向量組成的矩陣,y的對角元素代表特徵值。
5.求矩陣的能觀性與能觀性
CONT=ctrb(A,B) //能控性
rank(CONT);
CONT爲觀測矩陣,返回的秩爲滿秩則可觀
CONT=obsv(A,B) //能觀性
rank(CONT);