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);