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

：木文主要讨论线性系统解集的几何结构与系统能观性、 能控性和稳定性之间的关系。这一关系从两个方面来说明，第一部分 讲述系统解集几何结构与特征值和特征向量之间的关系，通过Matlab 仿真例子说明这一关系；第二部分分别讲述特征值和特征向量与系统 能观性、能控性和稳定性之间的关系，并讲述了能观性，能控性以及 稳定性的定义和判据，通过以约旦标准型为例来讲述相同特征值和不 同特征值情况下的能观性，能控性，最后在Simulink中仿真一定特 征值条件下系统的稳定性。从以上两个方面来说明解集的几何结构与 系统能观性、能控性和稳定性之间的关系。 1.零输入响应解集与特征值和特征向量之间的关系 一 x = Ax^Bu 线性定常系统状态方程『七+加宀(0)’八0的解为 x(t) = eAfx0 + J eA{l~r} Bu(T)dr.t 0。 为了研究线性泄常系统状态方程解集的几何结构与线性系统的特征之I可的关系，将 x = Ax 系统简化，只考虑系统为零输入的状态响应，即{ ° , x(O) = xo,r 0的解为 [y = Cx x(t) = eAfx0 o 所有的零输入状态响应组成了一个线性空I可，且该线性空I可屮有n个独立的元素, 它们的线性组合决定了所有零输入响应。所以可以通过选择一组线性独立的初始条件得 到一组零输入响应集中的基底。 下而先考虑最简单的零输入状态响应集的基底。 若九儿，…人是A的两两互异的特征值，且片宀，??必是相应的单位特征向量， 即 Av. = i = o 选兀0=片丿=12则 x(t) — eAtx0 — eAtv. TOC \o 1-5 \h \z = (/ +Az+ —A2/2 + — A3t3 +???)％ 2! 3! 1 =v + Av.t + ——A2v.t2 + + … ，，2! 1 3! 1 =岭+?叩 + ￡?# +￡人5” + … = 所以取兀o=0片+吋2+???+ 0必 时，相应的零输入响应为 x(t) = a1eAltv1 + a2e^fv2 + …+ ocne^vn 由此可以看岀线性定常系统的零输入响应解集的儿何结构可以由系统矩阵A的特 征值和特征向量来表征。即其解集由vI?v2?...vz,构成的n维坐标空间的线性组合。 上述结论的Matlab仿真程序和结果如下: 系统的状态方程为|一0—2丿(°丿取初始状态x(0)=a1 系统的状态方程为|一 0 —2丿 (°丿 取初始状态x(0)=a1*V1+a2 *V2,其 y = (2 3)x 打开MATLAB编辑器，编写如下程序， clear all; close all; clc; A=[-l 0;2 -2];B=[l;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(l ,length(t));% 输入为零 xO=alpha(l)*V(:,l )+alpha(2)* V(:,2);% 设定初始状态 subplot(211); [y,t,x]=lsim(sys,u,t,xO); t=0:0.01:5; plot(t,x);%绘制系统的零输入响应 xlabef时间t’titlcC原系统的零输入响应状态)； grid on; %绘制以特征值对应的特征向量为基底的零输入响应的状态 xz=alpha(l)* V(:, 1 )*exp(D( 1 )*t)+alpha(2)* V(:,2)*exp(D(4)*t); subplot(212);plot(t,xz); xlabclC时间f);titlcC以特征值对应的特征向量为基底的零输入响应的状态)； grid on; 输出结果如下图：由图可知，当系统的初始状态为系统特征值对应的特征向量的线 性组合时，系统的零输入响应也是相应特征向量的对应线性组合。即当x(O)=oti*V]+a2 *V2 时，eAt*x(0)= a1*eXI*t*V1+a2*eX2*t*V2o 2特征值、特征向量与系统特征之间的关系 2.1系统能控性与系统特征值和特征向量之间的关系 2.1.1能控性的定义 状态空I可中的任意两点Xo，Xi即X(to)=Xo，X(tj)=Xi，若存在控制信号u(t)能将状态 x(t)从to时刻的Xo在［to.切中驱动到XI，则系统能控。 2.1.2能控性判据 对于系统 x(Z)= Ax(t)+ Bu(t), 其响应为X(r) = % 「eA{thT}Bu(T)dT J% ⑴Wc矩阵为