线性系统理论大作业2.doc

摘要：本文主要讨论线性系统解集的几何结构与系统能观性、能控性和稳定性之间的关系。这一关系从两个方面来说明，第一部分讲述系统解集几何结构与特征值和特征向量之间的关系，通过Matlab仿真例子说明这一关系；第二部分分别讲述特征值和特征向量与系统能观性、能控性和稳定性之间的关系，并讲述了能观性，能控性以及稳定性的定义和判据，通过以约旦标准型为例来讲述相同特征值和不同特征值情况下的能观性，能控性，最后在Simulink中仿真一定特征值条件下系统的稳定性。从以上两个方面来说明解集的几何结构与系统能观性、能控性和稳定性之间的关系。

1. 零输入响应解集与特征值和特征向量之间的关系
线性定常系统状态方程，的解为
。
为了研究线性定常系统状态方程解集的几何结构与线性系统的特征之间的关系，将系统简化，只考虑系统为零输入的状态响应，即，的解为。
所有的零输入状态响应组成了一个线性空间，且该线性空间中有n个独立的元素，它们的线性组合决定了所有零输入响应。所以可以通过选择一组线性独立的初始条件得到一组零输入响应集中的基底。
下面先考虑最简单的零输入状态响应集的基底。
若是A的两两互异的特征值，且是相应的单位特征向量，即。选，则

所以取时，相应的零输入响应为

由此可以看出线性定常系统的零输入响应解集的几何结构可以由系统矩阵A的特征值和特征向量来表征。即其解集由构成的n维坐标空间的线性组合。
上述结论的Matlab仿真程序和结果如下：
系统的状态方程为，取初始状态x(0)=α1*V1+α2 *V2，其中V1、V2为特征值对应的特征向量。 取[α1 α2]=[1 4]。
打开MATLAB编辑器，编写如下程序，
clear all;
close all;
clc;
A=[-1 0;2 -2];B=[1;0];C=[2 3];D=[0];%设定系统的状态方程参数
sys=ss(A,B,C,D);
[b,a]=ss2tf(A,B,C,D);%状态方程转换成传递函数
alpha =input('input alpha:');%输入alpha参数[1 4]
[V,D]=eig(A); %利用eig函数求系统的特征值D与特征向量V
t=0:0.01:5;u=zeros(1,length(t));%输入为零
x0=alpha(1)*V(:,1)+alpha(2)*V(:,2);%设定初始状态
subplot(211);
[y,t,x]=lsim(sys,u,t,x0);
t=0:0.01:5;
plot(t,x);%绘制系统的零输入响应
xlabel('时间t');title('原系统的零输入响应状态');
grid on;

%绘制以特征值对应的特征向量为基底的零输入响应的状态
xz=alpha(1)* V(:,1)*exp(D(1)*t)+alpha(2)* V(:,2)*exp(D(4)*t);
subplot(212);plot(t,xz);
xlabel('时间t');title('以特征值对应的特征向量为基底的零输入响应的状态');
grid on;
输出结果如下图：由图可知，当系统的初始状态为系统特征值对应的特征向量的线
性组合时，系统的零输入响应也是相应特征向量的对应线性组合。即当x(0)=α1*V1+α2 *V2时，eAt*x(0)= α1*eλ1*t*V1+α2*eλ2*t*V2。

2 特征值特征向量与系统特征之间的关系
2.1　系统能控性与系统特征值和特征向量之间的关系
2.1.1 能控性的定义
状态空间中的任意两点x0，x1即x(t0)=x0，x(t1)=x1，若存在控制信号u(t)能将状态x(t)从t0时刻的x0在[t0，t1]中驱动到x1，则系统能控。
2.1.2 能控性判据
对于系统,
其响应为
(1) Wc矩阵为非奇异矩阵系统能控
在t1时刻，，
设计控制器，
则
=
假设
则
若能保证Wc为非奇异矩阵，则能保证系统的能控性。这一命题的逆命题也成立，证明过程略。
(2)矩阵C=[B AB A2B…AN-1B]满秩系统能控
(3)PHB判据:矩阵[A-λI B]满秩系统能控
(4)AP+PAT=-BBT有唯一正定解系统能控
(5)矩阵B中所有列张成的子空间不属于A的任意一个不变子空间系统能控
不变子空间：若V是线性空间X的子空间，是V中任意一个元素，若对所有，（A为矩阵）则V称为A的不变子空间。
2.1.3 特征值和特征向量与能控性的关系
对于状态方程，由第一节中讨论可知，其零输入响应的解集是在由ξ1，ξ2，ξ3…组成的线性空间中。以下讨论以三阶系统为例说明能控性与特征值特征向量的关系。
设初始条件为x（0）=α1*ξ1+α2*ξ2+α3*ξ3，若B=α1*ξ1+α2*ξ2，则定义：

由以上推导可知，不满秩，系统不能控。可以得出结论：B矩阵的方向决定了能控性。即能控性判据（5），（A,B）能控需要B中所有列张成子空间不属于A的任意一个不变子空间（以特征向量为基底的线性空间）。
例如：对于系统，A矩阵的特征值为λ1=-1λ2=-2,
λ3=-3。对应的特征向量为ξ1=（1 0 0）T,ξ2=(0 1 0)T,ξ3=(0 -1 1)T。B=（ξ1+ξ2 ξ1+2ξ2），B的所有列张成子空间属于A的由ξ1=（1 0 0）T和ξ2=(0 1 0)T组成的不变子空间。根据以上结论可知，系统不能控。以下根据判据（2）来判断系统是否能控：
B AB A2B）=不满秩，系统不能控，与以上结论一致。

若B所有列张成的子空间属于A的某一个不变子空间，以下讨论对（A，B）进行能控分解。
若（A，B）不能控，系统可分解为

（1）当输入u=0 ，初始状态时，
则，从而。
说明输入为零，初始状态落在能控子空间里，则状态响应在能控子空间里。
（2）当出入u≠0，初始状态时，
则，。
（3）当出入u≠0，初始状态，都不为零时，
则，

说明初始条件不在特征向量张成子空间里，则状态响应不完全在能控子空间里。综上所述：当初始条件落在能控子空间，则一定能设计控制器使状态点到达预定点。
2.1.4 A的约旦标准型判断能控性
要判断一个系统基于约旦标准型的能控性，必须考虑其约旦标准型和输入矩阵的形式，这里通过举例说明这种方法。如果A不是约旦标准型，则可以通过线性变换将A化为约旦标准型。
对于系统
变换为约旦标准型
我们可以求出系统矩阵A的特征值，来直接写出系统的约旦标准型矩阵J。
当特征值无重根时，

当特征值有q个重根时，

对于已经通过线性变换转化成约旦标准型的系统：

当0时，是系统的能控模态；
当与相互线性独立时，是系统的能控模态。
2.2 系统能观性与系统特征值和特征向量之间的关系
2.2.1 能观性定义
若在给定时间段[t0,t1]内，根据系统的输入输出信号可唯一确定系统的初始
状态，则系统能观。
2.2.2 系统能观性判据
对于系统，其响应为：
系统输出为
(1) Wc矩阵为非奇异矩阵系统能观
令
两边乘以得

两边积分得

记 则

以上推导说明：当W0为非奇异矩阵时系统能观。该命题的逆命题也成立，证明略。
（2）能观性矩阵满秩系统能观。
（3）矩阵满秩系统能观。
（4）方程有唯一的正定解系统能观。
2.2.3 特征值和特征向量与能观性的关系
对于状态方程，由第一节中讨论可知，其零输入响应的解集是在由
ξ1，ξ2，ξ3…组成的线性空间中。以下说明能观性与特征值特征向量的关系。
根据A的特征值来判断系统的能观性，如果A不是约旦标准型，则可以通过线性变换将A化为约旦标准型。如下例子可以说明，对于已经通过线性变换转化成约旦标准型的系统，

当0时，是系统的能观模态；
当与相互线性独立时，是系统的能观模态。
不能观子空间：若x0是系统的初始状态从t0至t1，x0所产生的零输入响应为零，所有上述x0组成了一个子空间称为不能观子空间。
例如：，初始状态为，则，得，，系统不能观。以不能观模态对应向量为初始条件的输出为零。
2.3 系统稳定性与系统特征值和特征向量之间的关系
2.3.1 系统稳定性
定义：系统在有界输入条件下，输出有界。
系统的输入有界是指对u(t)存在一个常数满足，。对所有上述输入信号，存在一个常数M，使系统输出y(t)满足，。则称系统是有界输入有界输出稳定。
2.3.2 稳定性判据
系统脉冲响应y(t)绝对可积，即系统稳定。
2.3.3 约旦标准型的特征根与系统稳定性关系
研究线性定常系统的稳定性的方法很多，因，的解为
，所以可以根据的形态即可判别系统的稳定性。由于的形态是由决定的，因此可以用来研究系统的稳定性。除此之外，还可以根据A的特征值来研究系统的稳定性。
下面给出线性定常系统渐近稳定的判据：
定理1 [特征值判据] 线性定常系统为渐近稳定的充要条件是系统矩阵A的全部特征值都具有负实部。

定理2线性定常连续系统为渐近稳定的充要条件是对正定矩阵Q，方程有唯一的正定解。
证明：先证充分性。即证明线性定常系统为渐近稳定，则对任意的正定矩阵Q, 方程有唯一的正定解。
即证明若对任意的正定矩阵Q, 方程有唯一的正定解，则 是渐近稳定的。
令，则

所以P是上述方程的解。
若均为上述方程的解，即
, 则
等式两边同乘、得

积分得

再证必要性。
即证明若对任意的正定矩阵Q, 方程有唯一的正定解，则是渐近稳定的。
取A的特征值和其相应的特征向量，将方程两边同乘以、得

所以由定理1知A的全部特征值都具有负实部，则系统渐近稳定。
现在以一个例子来说明。线性定常系统为，其中A不是约旦标准型，由于线性变换不改变系统的特征值，进而不改变系统的稳定性，所以可通过线性变换将其化为约旦标准型，直接根据特征值的情况来讨论系统的稳定性。
于是
, 设，则。
若
则
若要稳定，则要求有负实部。
若,则
有界但不趋近于0，所以系统稳定但不是渐近稳定。
由上可得结论：
（1）若系统化为约旦标准型后，如有某个约旦标准型的特征值在s的左半平面（有负实部），则该约旦块渐近稳定。
（2）若某一个约旦块有一个纯虚根，则系统稳定但不是渐近稳定。

以下用simulink仿真例子说明系统稳定性：
以系统为例，在Simulink中仿真，设置初始状态为x1=[1;1];x2=[1;1]

打开示波器显示如下图形：

由图可以看出，当A中有一个不稳定的特征值时，系统的输出不稳定，对应的状态x1不稳定，而状态x2稳定。

13