实验三 系统的可控性与可观测性分析解析Word文档下载推荐.docx
《实验三 系统的可控性与可观测性分析解析Word文档下载推荐.docx》由会员分享,可在线阅读,更多相关《实验三 系统的可控性与可观测性分析解析Word文档下载推荐.docx(9页珍藏版)》请在冰豆网上搜索。
定义:
若对状态空间中任一非零初态x(t0),存在一个有限时刻t1>
t0,使得由输入u[t0,t1]和输出y[t0,t1]能够唯一确定初始状态x(t0),则称动态方程
在t0时刻是可观测的。
反之称为是不可观测的。
(二)、可控性和可观测性判据
1、可控性
构造一个相似变换矩阵
公式中,是系统的阶次;
矩阵称为系统的可控性变换矩阵。
矩阵可以由控制系统工具箱中提供的函数来产生。
其调用格式为
公式中,的秩,即称为系统的可控性指数,它的值表示系统中可控制的状态的数目。
如果,则系统是完全可控制的。
【例题1】考虑系统的状态方程模型为
分析系统的可控性。
A=[0,1,0,0;
0,0,-1,0;
0,0,0,1;
0,0,5,0]
B=[0;
1;
0;
-2]
Tc=ctrb(A,B)
rank(Tc)
结果如下:
>
rank(Tc)
ans=
4
可见,系统完全能控。
2、可观测性
构造一个相似变换矩阵如下
公式中,是系统的阶次。
矩阵称为系统的可观测变换矩阵。
公式中,的秩,即,称为系统的可观测性指数,它实际上是系统中可观测状态的数目。
如果,则系统是完全能观测的。
【例题2】考虑系统的状态方程模型为
分析系统的可观测性。
C=[1,0,0,0]
运行结果如下:
rank(To)
可见,系统是完全可观测的。
(二)、可控性和可观性的标准型实现
先来看什么是实现,所谓实现,就是根据描述系统输入输出动态关系的传递函数建立系统的状态空间表达式,所求得的状态空间表达式保持原来传递函数的输入输出关系不变,同时反映内部动态变化。
实现不是唯一的。
下面看一个实现的例子
【例题3】有以下状态空间模型
A=[1.25,-4,-1.25,0.5;
1.25,-0.5,-1.25,0.25;
0.25,-4.25,-1.25,0.5;
2.25,1.75,0.25,1]
B=[4,6;
3,4;
2,2;
1,0]
C=[0,0,0,1;
0,2,0,2]
D=zeros(2,2)
Gss1=ss(A,B,C,D)
Gtf=tf(Gss)%可以这样认为:
Gss就是Gtf的一个实现,有4个状态变量
现在我们继续用Gtf来完成一个实现Gss2,命令如下
Gss2=ss(Gtf)
看一下实现结果,这个实现有8个状态变量,它当然没有前面的4个状态变量的实现要好,虽然它们表示同一个系统。
大家不禁要问,到底那个实现好,还有没有标准了?
最小实现就是回答了这个问题。
所谓最小实现就是实现的阶次最低,或最低阶次的实现。
matlab最小实现函数为
Gmin=minreal(G)
其中,G为原系统的LTI对象,G1为最小实现后的LTI对象。
【例题4】对上例中的Gss2,求出最小实现命令为:
Gmin=minreal(Gss2)
从结果可以看是,系统阶次回到了4阶。
1、可控标准型I
设单输入系统的状态方程为:
设A的特征多项式
如果系统状态完全能控性
则可以通过线性变换,可以将其变成如下形式的能控标准形。
相似变换
变换矩阵为
此可控标准型也称为可控标准型I。
【例题5】已知能控的线性定常系统
求其能控标准型。
A=[101;
010;
100]
1]
C=[110]
(1)、能控性矩阵
Tc=ctrb(A,B)
系统完全可控。
(2)、A的特征多项式
symss
det(s*eye(3)-A)
结果:
s^3-2*s^2+1
相应系数为
(3)、计算变换矩阵
(4)、计算标准型
Abar=inv(Tc1)*A*Tc1
Bbar=inv(Tc1)*B
Cbar=C*Tc1
得:
即:
2、可控标准型
如果取相似变换矩阵为可控性矩阵,即
则,原状态空间表达式可变换为如下的可控标准型。
【例题6】将上面的例题变换为可控标准型。
D=0
Tc=ctrb(A,B)
[Ac,Bc,Cc,Dc]=ss2ss(A,B,C,D,inv(Tc))
Ac=
00-1
100
012
Bc=
1
0
Cc=
122
Dc=
(3)、能观测标准形
系统的能观测性矩阵为
,若,则系统能观测,通过线性变换可以将其变成如下形式的能观标准形。
变换矩阵可取为
同上,这个也可以称为能观测标准型
【例题7】将上面例题变换为能观测标准型
和前面一样,
取变换矩阵如下:
[Ao,Bo,Co,Do]=ss2ss(A,B,C,D,P)
最后结果为:
Ao=
Bo=
-2
Co=
001
Do=
还可以直接把可观测性矩阵取为变换矩阵,这样得到所谓的可观测标准型,matlab命令如下:
P=obsv(A,C)
010
-102
2
三、实验方法及步骤
打开计算机,运行MATLAB软件。
将上述内容写入程序编辑窗口并运行。
分析结果,写出实验报告
四、【作业】
1.从可控标准型和可观测标准型的A矩阵,可以得到什么结论?
2.已知系统的系数矩阵为,试判断它的可控性。
如果完全可控,将其转化为可控标准型。