线性时变系统状态方程的解
线性离散系统状态方程的解
Z变换法(2/7)
在Z反变换中对标量函数存在下述公式和性质:
Z 1 1 /(1 az 1 ) a k Z {W1 ( z )W2 ( z )} w1 (k i ) w2 (i )
1 i 0 k
Βιβλιοθήκη 其中W1(z)和W2(z)分别为w1(k)和w2(k)的Z变换。 将上述公式推广到向量函数和矩阵函数,则可得
离散卷积
Z变换法(3/7)—例3-14
因此,离散系统的状态方程的解为:
x(k ) G k x(0) G k j 1Hu( j )
j 0
k 1
该表达式与前面递推法求解结果一致。 例 已知某系统的状态方程和初始状态分别为 1 0 1 1 x(k 1) x( k ) u ( k ) x(0) 0.16 1 1 1
( k 1 , k0 ) G ( k ) ( k , k 0 ) ( k0 , k0 ) I
其解为
(k , k0 ) G(k 1)G(k 2)...G(k0 ) , k k0
线性时变离散系统状态方程的解(3/6)
与线性定常离散系统类似,线性时变离散系统的状态求解公 式可用迭代法证明。 对线性时变离散系统的状态方程，依次令k= k0, k0+1, k0+2, …,从而有
2. 引入状态转移矩阵概念和表示之后,线性连续系统和线 性离散系统的状态方程的求解公式在形式上一致,都由 零输入响应和零状态响应叠加组成, 只是相应的零状态响应在形式上略有不同,一为求 积分(卷积),一为求和(离散卷积),但本质是一致的。 3. 在由输入所引起的状态响应中,第k个时刻的状态只取决 于此采样时刻以前的输入采样值,而与该时刻的输入采 样值u(k)无关。
Matlab程序解现代控制理论与工程中的状态方程
Matlab程序解现代控制理论与工程中的状态方程用矩阵指数法解状态方程的matlab函数vslove1：函数vslove1:求解线性定常连续系统状态方程的解function[phit,phitbu]=vsolves1(a,b,ut)%vsolves1谋线性已连续系统状态方程x’=ax+bu的求解%[phit,phitbu]=vsolves1(a,b,ut)%a,b系数矩阵%ut控制输入，必须为时域信号的符号表达式，符号变量为t%phit――输出phi(t)%phitbu――输入phi(t-tao)*b*u(tao)在区间（0，t）的分数symsttao%定义符号变量t,taophit=expm(a*t);%求矩阵指数exp(at)if(b==0)b=zeros(size(a,l),l);%重构系数矩阵bendphi=sub(phit,’t’,’t-tao’);%求exp[a(t-tao)]phitbu=int(phi*b*ut,’tao’,’0’,’t’);%谋exp[a(t-tao)]*b*u(tao)在0~t区间的分数用拉氏变换法解状态方程的matlab函数vslove2：函数vslove2:解线性定常已连续系统状态方程的求解function[sl_a,sl_abu]=vsolves1(a,b,us)%vsolves2求线性连续系统状态方程x’=ax+bu的解%[sl_a,sl_abu]=vsolves1(a,b,ut)%a,b系数矩阵%us掌控输出，必须为拉氏转换后的符号表达式，符号变量为s%sl_a――输入矩阵(sl-a)^(-1)拉式反华转换的结果%sl_abu――输出(sl-a)^(-1)*b*u(s)拉式反变换后的结果symss%定义符号变量t,taoaa=s*eye(size(a))-a;%谋si-ainvaa=inv(aa);%求（si-a）矩阵的逆intaataa=ilaplace(intaa);%求intaa的拉氏反变换si_a=simplify;%简化拉式反变换的结果if(b==0)b=zeros(size(a,l),l);%重构系数矩阵bendtab=ilaplace(intaa*b*us);%求intaa*b*us的拉氏反变换si_abu=simplify(tab);%化简拉式反变换的结果解时变系统状态方程的matlab函数tslove：函数tslove:求解线性时变连续系统状态方程的解function[phi,phibu]=tsolves(a,b,u,x,a,n)%tsolves求时变系统状态方程%[phi,phibu]=vsolves1(a,b,u,x,a,n)%a,b时变系数矩阵%phi――状态迁移矩阵计算结果%phibu――THF1求解分量%u――控制输入向量，时域形式%x――符号变量，阐明矩阵a中的时变参数，通常为时间t%a――积分下限%n――时变状态转移矩阵中计算重积分的最大项数，n=0时无重积分项%n=1时包含二重积分项，.....phi=transmtx(a,x,a,n);%排序状态迁移矩阵phitao=subs(phi,x,’tao’);%谋phi(tao)if(b==0)btao=zeros(size(a,l),l);%谋b(tao)endutao=subs(u,x,’tao’);%求u(tao)phibu=simple(int(phitao*btao*utao,’tao’,a,x));%排序THF1分量求解时变系统转移矩阵的matlab函数transmtx：函数transmtx:解线性时变系统状态迁移矩阵functionphi=transmtx(a,x,a,n)%transmtx计算时变系统状态转移矩阵%phi=transmtx(a,x,a,n)%phi――状态迁移矩阵计算结果%a时变系数矩阵%x――符号变量，指明矩阵a中的时变参数，通常为时间t%a――积分下限%n――时变状态迁移矩阵中排序轻分数的最小项数，n=0时并无轻分数项%n=1时涵盖二重积分项，.....phi=eye(size(a));%初始化phi=iforlop=0:naa=a;fori=1:lopif(aa==0)break;endatemp=subs(aa,x,’taoi’);aa=simplify(a*int(atemp,’tao’,a,x));endif(aa==0)break;endatemp=subs(aa,x,’taoi’);aa=simplify(a*int(atemp,’tao’,a,x));%计算重积分phi=simplify(phi+aa);%修正phiend解线性定常离散系统状态方程的matlab函数disolve：函数disolve:求解线性定常离散系统状态方程的解function[ak,akbu]=disolve(a,b,uz)%disolve谋线性离散系统状态方程x(k+1)=ax(k)+bu(k)的求解%[ak,akbu]=disolve(a,b,uz)%a,b系数矩阵%uz控制输入，必须为z变换后的符号表达式，符号变量为z%ak――输出矩阵[((zi-a)^(-1)z]z反变换后的结果%akbu――输入矩阵[((zi-a)^(-1)*b*u(z)]z反华转换后的结果symsz%定义符号变量zaa=z*eye(size(a))-a;%求zi-ainvaa=inv(aa);%谋（zi-a）矩阵的逆intaataa=iztrans(intaa*z);%谋intaa*z的z 反华转换ak=simple(taa);%精简z反华转换的结果if(b==0)b=zeros(size(a,l),l);%重构系数矩阵bendtab=iztrans(intaa*b*uz);%谋intaa*b*uz的z反华转换akbu=simple(tab);%化简z 反华转换的结果求解线性时变离散系统状态方程的matlab函数tdsolve：函数tdsolve:解线性时变离散系统状态方程的求解functionxk=tsolve(ak,bk,uk,x0,kstart,kend)%tdsolve求线性时变离散系统状态方程x(k+1)=a(k)x(k)+b(k)u(k)的解%xk=tsolve(ak,bk,uk,x0,kstart,kend)%ak,bk系数矩阵%uk掌控输出，必须为时域符号表达式，符号变量为k%x0初始状态%kstart――起始时刻%kend――中止时刻%xk――输出结果，矩阵每一列分别对应x(k0+1),x(k0+2)....symsk%定义符号变量kif(bk==0)bk=zeros(size(a,l),l);%重构系数矩阵bendxk=[];forkk=kstart+1:kendaa=eye(size(k));fori=kstart:kk-1%排序a(k-1)a(k-2)....a(k0+1)a(k0)a=subs(ak,’k’,i);aa=a*aa;endaab=eye(size(ak));bb=zeros(size(bk));fori=kk-1:-1:kk+1%排序a(k-1)a(k-2)....a(j+1)b(j)u(j)的递增和a=subs(ak,’k’,i);aab=aab*a;b=subs(bk,’k’,kk-1+i+kstart);u=subs(uk,’k’,kk-1+i+kstart);bb=bb+aab*b*u;endb=subs(bk,’k’,kk-1);u=subs(uk,’k’,kk-1);bb=bb+b*u;xk=[xkaa*x0+bb];%计算x(k)end已连续系统状态方程线性化后的matlab符号函数sc2d：函数sc2d:线性连续系统状态方程的离散化function[ak,bk]=sc2d(a,b)%sc2d线性化线性已连续系统状态方程x’=ax+bu%sysd=sc2d(a,b)%a,b――连续系统的系数矩阵%ak,bk――离散系统系数符号矩阵%线性状态方程为：x(k+1)=ak*x(k)+bk*u(k)%ak,bk中变量t为取样周期symstt%定义符号变量ttphit=expm(a*t);%求矩阵指数exp(at)if(b==0)b=zeros(size(a,l),l);%重构系数矩阵bendphitb=int(phit*b,’t’,0,’t’);%求exp(at)*b在0~t区间的积分ak=simple(subs(phit,’t’,’t’));bk=simple(phitb);线性时变系统线性化后的matlab函数tc2d：函数tc2d:线性时变系统的离散化function[ak,bk]=tc2d(a,b,x,n)%tc2d线性时变系统的离散化%[ak,bk]=tc2d(a,b,x,n)%a,b――已连续系统的系数矩阵%ak――离散化后的系数矩阵a(kt)%bk――离散化后的系数矩阵b(kt)%x――符号变量，阐明矩阵a\\b中的时变参数，通常为时间t%n――时变状态转移矩阵中计算重积分的最大项数，n=0时无重积分项，%n=1时包含二重积分项，.....symsttphit=transmtx(a,x,k*t,n);%计算时变系统的状态转移矩阵ak=simplify(subs(phi,x,(k+1)*t));%计算离散化后的系数矩阵a(kt)phitao=subs(phi,x,’tao’);%谋phi(tao)if(b==0)btao=zeros(size(a,l),l);elsebtao=subs(b,x,’tao’);%谋b(tao)endphitb=simple(int(phitao*btao,’tao,k*t,x’));%计算受控分量bk=simplify(subs(phib,x,(k+1)*t));%排序线性化后的系数矩阵b(kt)定常系统可控规范i型变换函数ccanonl：函数ccanonl:谋线性定常系统的受控规范i型形式function[abar,bbar,cbar,t]=ccanonl(a,b,c)ìanonl求系统x’=ax+bu，y=cx的可控规范i型系数矩阵?ar,bbar,cbar,――变换后的可控规范i型系数矩阵%t――相似变换矩阵n=length(a);co=ctrb(a,b);if(rank(co)~=n),%判断系统可控性error(‘系统不可控！’);。
控制工程技术基础 第7章现代控制理论简介
7.2控制系统的状态空间表达式
7.2.1状态、状态变量
状态：系统运动信息的集合。 状态变量：可以完全确定系统的运动状态且数目最小的一组变量。所 谓完全确定，是指只要给定t0时刻的这组变量的值和系统在t ≥t0时系 统的输入函数，则系统在t > t0的任意时刻的状态就可完全确定。所谓 数目最小是指：如果变量数目大于该值，则必有不独立的变量；小于 该值，又不足以描述系统的运动状态。 状态向量：n个状态变量x1 (t)，x2 (t)，…， xn (t)所构成的向量X(t)就 是系统的状态向量，记作X(t)=[x1 (t)，x2 (t)，…， xn (t)]T
上一页
ห้องสมุดไป่ตู้下一页
返回
7.4最优控制
以上可见，邦特略京极小值原理实际上是把一个求性能指标J的 最小值问题，转化成一个求哈密顿函数H的最小值问题。 当系统的状态方程为
第7章现代控制理论简介
7.1概述 7.2控制系统的状态空间表达式 7.3状态反馈与输出反馈 7.4最优控制
7.1概述
现代控制理论的基本内容包括五个方面，简单说明如下。 1.最优控制 在图7-1所示系统中，有一组输入函数u (t)作用在受控系统上，其 相应状态变量是x (t)，通过量测系统可得到这些状态的某种组合y (t)， 此即系统输出。根据实际需要，可为受控系统指定一些目标（性能指 标）。 2.最优估计 图7-1所示系统中，输出量y (t)是通过量测系统由状态转换过来 的。但实际的量测系统常受到噪声v (t)的干扰，如图7-2所示。如果将 整个系统看成是一个信息传递系统，用输入噪声w( t)表示这个系统的 模型误差，也称动态噪声，则从y (t)中，克服w( t)和v (t)的影响估计 出状态x (t)来，称为最优状态估计问题。
现代控制理论-状态方程的解
3、复频域上
非齐次状态方程的解
2、说明
e At 状态转移矩阵
一般用 t 表示，即 t e At
考虑初始条件拉氏变换
sX ( s ) X (0 ) AX ( s ) BU ( s ) 有 ( sI A) 1 X ( s ) X ( 0 ) BU ( s ) 即 1 X ( s ) ( sI A) X (0) ( sI A) 1 BU ( s ) 则
e
d At e Ae At e At A dt
At 1
e At
[5]、对于 n n的方阵 A、 B 当且仅当 AB BA时 有 e At e Bt e( A B)t , 而当AB BA, e At e Bt e( A B)t。
电气工程学院
几个特殊的矩阵指数eAt
设单变量系统的差分方程为：
y(k n) an1 y(k n 1) a0 y(k ) bnu(k n) bn1u(k n 1) b0u(k )
相应的系统脉冲传递函数为
bn z n bn 1 z n 1 b1 z b0 G( z ) n z an 1 z n 1 a1 z a0
有
d At At AX ] e X e [X dt e At Bu(t )
考虑初始条件 拉氏变换得 sX ( s ) X ( 0 ) AX ( s )
将上式积分有 t t X (t ) 1 ( sI A) 1 X (0) A d A e Bu( ) d d e X ( ) 0 0 d 1 显然 e At 1 t ( sI A) At A X ( 0 ) e X ( t ) e 可得 At Bu( )d
《自动控制原理》线性系统的可控性与可观测性
将状态 x(t0 ) = 0 转移到 x(t f ) =x f 的控制作用，则称状态 x f 是 t0 时刻 可达的。若x f 对所有时刻都是可达的，则称状态x f 为完全可达或 一致可达。若系统对于状态空间中的每一个状态都是时刻 t0 可达的， 则称该系统是 t0 时刻状态完全可达的，或简称该系统是 t0 时刻可达
可观测性问题: 相应地，如果系统所有状态变量的任意形式 的运动均可由输出完全反映，则称系统是状态可观测的，简称为系 统可观测。反之，则称系统是不完全可观测的，或简称为系统不可 观测。
可控性与可观测性概念，是卡尔曼于20世纪60年代首先提出 来的，是用状态空间描述系统引伸出来的新概念，在现代控制理论 中起着重要的作用。它不仅是研究线性系统控制问题必不可少的重 要概念，而且对于许多最优控制、最优估计和自适应控制问题，也 是常用到的概念之一。
在研究可观测性问题时，输出 y 和输入 u 均假定为已知，只有初始
状态 x0 是未知的。因此，若定义
t
y(t) = y(t) − C(t) (t, )B( )u( )d − D(t)u(t) t0
则式（9-79）可写为
y(t) = C(t)(t,t0 )x0
（9-80）
这表明可观测性即x0 可由 y 完全估计的性能，由于 y 和 x0 可任意取
y = −6x2
这表明状态变量 x1 和 x2 都可通过选择控制量 u 而由始点达到原
点，因而系统完全可控。 如何判别？
但是，输出 y 只能反映状态变量 x2 ，而与状态变量 x1 既无直
接关系也无间接关系，所以系统是不完全可观测的。如何判别？
变化：（1）b1=0 ? （2）a12≠0 ？ (3) a21≠0 ?
值，所
现代控制理论基础第二章习题答案
第二章 状态空间表达式的解3-2-1 试求下列矩阵A 对应的状态转移矩阵φ（t ）。
（1） ⎥⎦⎤⎢⎣⎡-=2010A （2） ⎥⎦⎤⎢⎣⎡-=0410A （3） ⎥⎦⎤⎢⎣⎡--=2110A （4） ⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-=452100010A （5）⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡=0000100001000010A （6）⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡=λλλλ000100010000A 【解】：（1） （2） （3） （4）特征值为：2,1321===λλλ。
由习题3-1-7(3)得将A 阵化成约当标准型的变换阵P 为⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=421211101P ，⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡----=-1211321201P线性变换后的系统矩阵为：（5）为结构四重根的约旦标准型。
（6）虽然特征值相同，但对应着两个约当块。
或}0100010000{])[()(1111----⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡------=-=Φλλλλs s s s L A sI L t 3-2-2 已知系统的状态方程和初始条件 （1）用laplace 法求状态转移矩阵； （2）用化标准型法求状态转移矩阵； （3）用化有限项法求状态转移矩阵； （4）求齐次状态方程的解。
【解】：（1） （2）特征方程为： 特征值为：2,1321===λλλ。
由于112==n n ，所以1λ对应的广义特征向量的阶数为1。
求满足0)(11=-P A I λ的解1P ，得：0110000000312111=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡--P P P ，⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=0011P 再根据0)(22=-P A I λ，且保证1P 、2P 线性无关，解得：对于当23=λ的特征向量，由0)(33=-P A I λ容易求得： 所以变换阵为：[]⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-==110010001321P P P P ，⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=-1100100011P 线性变换后的系统矩阵为：（3）特征值为：2,1321===λλλ。
线性系统状态空间分析和运动解
线性系统状态空间分析和运动解状态空间分析方法是一种用来描述线性系统的分析方法。
它将系统的动态特性用一组状态变量来表示，并通过矩阵形式的状态方程进行分析和求解。
状态空间方法是目前广泛应用于自动控制系统设计与分析的一种方法，它可以对系统的稳定性、可控性、可观性以及性能等进行定量分析。
在状态空间分析方法中，首先需要将系统的微分方程表示为矩阵形式的状态方程。
状态方程描述了各个状态变量和它们的变化率之间的关系。
假设系统有n个状态变量x1, x2, ..., xn和m个输入变量u1, u2, ..., um，状态方程可以表示为:dx/dt = Ax + Bu其中，dx/dt是状态变量的变化率，A是状态矩阵，描述状态变量之间的耦合关系，B是输入矩阵，描述输入变量对状态变量的影响。
状态空间分析方法的基本思想是将系统转化为状态空间表达式，然后通过对状态方程进行分析和求解来得到系统的特性和响应。
常见的分析方法包括对系统的稳定性、可控性和可观性进行评估。
稳定性是系统的基本性质之一，用来描述系统在受到扰动时是否能够恢复到平衡状态。
在状态空间方法中，通过研究系统的特征根（或特征值）可以判断系统的稳定性。
特征根是状态方程的解的根，系统的稳定性与特征根的实部有关。
如果特征根的实部都小于零，则系统是稳定的；如果特征根存在实部大于零的情况，则系统是不稳定的。
可控性是指系统是否可以通过输入变量来控制系统的状态变量。
在状态空间方法中，通过可控性矩阵来判断系统的可控性。
如果可控性矩阵的秩等于系统的状态变量个数，则系统是可控的；如果可控性矩阵的秩小于系统的状态变量个数，则系统是不可控的。
可观性是指系统的状态变量是否可以通过观测变量来测量得到。
在状态空间方法中，通过可观性矩阵来判断系统的可观性。
如果可观性矩阵的秩等于系统的状态变量个数，则系统是可观的；如果可观性矩阵的秩小于系统的状态变量个数，则系统是不可观的。
除了稳定性、可控性和可观性外，状态空间分析方法还可以用来分析系统的性能指标，如系统的响应时间、稳态误差和系统的最大误差等。
线性系统理论-郑大钟(3-4章)
1
2 n
n 1 n
t e n
1
0 1
21
n 1 2
(n 1)1 (n 1)(n 2) n 3 1 2! n2 (n 1)1 n 1 1 1
矩阵指数函数的算法 1：定义法
e At I At
1 2 2 A t 2!
只能得到eAt的数值结果，难以获得eAt解析表达式，但用计算机计算，具 有编程简单和算法迭代的优点。 2：特征值法
A P 1 AP
A PA P 1
e At Pe A t P 1
P为变换A为约当规范型的变换矩阵 1）若A的特征值为两两互异
如果系统矩阵A(t),B(t)的所有元在时间定义区间[t0,tα]上为时间t的连续实函数，输 入u(t)的所有元为时间t的连续实函数，那么状态方程的解x(t)存在且唯一。 从数学观点，上述条件可减弱为： ①系统矩阵A(t)的各个元aij(t)在时间区间[t0,tα]上为绝对可积，即：
t
t0
| aij (t ) | dt ,
－1
te1t 1t e e3t
0 2tet e 2t 1 3tet 2et 2e 2t 2 tet et e 2t
e At 0 I 1 A 2 A2 (2tet e 2t ) I (3tet 2et 2e 2t ) A (tet et e 2t ) A2 2et e 2t 0 e t e 2t 0 et 0 2et 2e 2t 0 et 2e 2t
s3 ( s 1)( s 2) 2 ( s 1)( s 2)
线性系统的运动分析
A 有二种标准形式： 对角线矩阵、约当矩阵
（1）当A的特征值 1, 2,, n 为两两相异时：对角线标准型
e1t
0
e At Te AtT1 T
T 1
0
ent
其中： T为使A化为对角线标准型的非奇异变换矩阵。
求状态转移矩阵的步骤：
1） 先求得A阵的特征值 i 。
2） 求对应于 i 的特征向量 pi ，并得到T阵及T的逆阵。
j0
其中： a0(t), a1(t),, an1(t) 为t的标量函数，可按A的特 征值确定。
1） A的特征值 1, 2 ,, n 两两相异时，
注意求逆
a0(t)
1 1
12
n1 1
1
e 1t
a1 ( t )
1
1
22
n1 2
e
2t
an1 ( t )
1
1
2n
n1 n
e
nt
三、几个特殊的矩阵指数函数
（1）设 A diag 1 n ，即A为对角阵且具有互异元素时，
有
e1t
Φt
e2t
0
0
ent
（2）若A能通过非奇异变换为对角阵时，即
T -1AT Λ
e1t
ΦΦ(t)t
T
e2t
0
0
T
-1
(2-9)
ent
1
0
（3）设A为 (n n) 约当阵，即 A
a0
a1
(t (t
) )
1 1
1
1
e
1t
2
e
2t
在第3种方法中已经求得特征根，所以得：
状态空间表达式解
2.1.1 齐次状态方程的解 u=0
X·=AX
1、直接求解 设 n=1
X(t0) =X0
x·=ax 解为x(t)=eatx0
t0=0 且eat=1+at+a2t2/2！+…
对于 n阶， 解为X(t)=eAtX0 eAt=I+At+A2t2/2！+…
证明：设X(t)解的形式为
=(I+At1+A2t12/2！+…) (I+At2+A2t22/2！+…) = (t1) (t2)
1、状态转移矩阵的性质：设t0=0 （4）[(t)]–1= (–t)
证明：由 （1）(0)=I （3） (t1+t2)= (t1) (t2) 得 (t–t)= (t) (–t)=I
(–t +t)= (–t) (t) =I 所以 [(t)]–1= (–t) （5） (t2– t1) (t1– t0) = (t2– t0)
0 0 0… 0 0 0…
e1t tm–1/(m–1)！ e1t tm–2/(m–2)！
Q–1 te1t e1t
以A有三重特征值为例进行证明
1 1 0 J= Q–1AQ= 0 1 1
0 0 1
证明 eAt=I+At+A2t2/2！+… 则 Q–1eAtQ=Q–1IQ+ Q–1 AtQ+ Q–1 A2t2/2！Q+… =I+ Jt+ J2t2/2！+… eAt=Q(I+ Jt+ J2t2/2！+…) Q–1
1 k!
Akb0
2.1.1 齐次状态方程的解
3.5线性时变系统状态方程的解
=G( kT) x( kT) + H( kT) u( kT) .
其中： 其中：
G( kT)
H( kT) TB( kT) .
I +TA( kT) ,
第三章 状态方程的解 3.6.2 线性时不变系统状态方程的离散化 考虑系统： 考虑系统： & x( t) = A ( t) + B ( t) , x u 其状态方程的解为： 其状态方程的解为：
第三章 状态方程的解 第一项是由初始状态引起的响应； 第一项是由初始状态引起的响应； 第二项是由控制输入引起的响应。 第二项是由控制输入引起的响应。
连续系统的时间离散化 3.6 连续系统的时间离散化
3.6.1 近似离散化 考虑系统
& x( t) = A( t) x( t) + B( t) u( t) ,
t t t0 t0
t A(τ ) d x( t ) x( t) =exp ∫ τ 0 t0
0 =exp 0
( t −t0 ) ( t +1) ( t0 +1) x( t0 )
0
1 = 0
( t −t0 ) ( t +1) ( t0 +1) x t ( 0)
1
第三章 状态方程的解
t 1 & 例3.5.2 x( t) = x( t) 初始值为 x( 0) .求 x( t)。 1 t
解：
t2 2 t ∫t0 A(τ ) dτ = t t2 2,
t
A( t) ∫ A(τ ) dτ = ∫ A(τ ) dτA( t) ,
线性时变系统状态方程的解
2.4
线性时变系统状态方程的解
线性时变系统的状态方程为 & x = A ( t ) x + B( t )u （1）不一定有解 （2）当A（t)、B(t)在定义区间绝对可积时，对 每一初始态存在唯一解。 0 1 & 时变系统中A、B随t变化，如 x= x 0 t 一、齐次矩阵微分方程的解 & x(t ) = A(t ) x(t ) x ( t ) t = t = x ( t o )
0
其解为： x ( t ) = φ( t, t 0 ) x ( t 0 )
φ( t , t 0 ) 是n*n阶非奇异方阵，且满足
& φ(t,t0) = Aφ(t,t0) φ(t0,t0) = I
证明：将解带入齐次方程
d [φ ( t , t 0 ) x ( t 0 ) ] = A ( t )φ ( t , t 0 ) x ( t 0 ) dt
x ( t ) = φ( t, t 0 ) x (t 0 ) +
∫
t
t0
φ( t , τ) B ( τ) u ( τ)dτ
证明（略）
三，状态转移矩阵 φ( t , t 0 ) 1. φ( t , t 0 ) 与 φ( t ) φ( t − t 0 ) 对比 共同：形式和性质类似 区别：本质不同，但Φ(t,t0)既是t的函数，也 是t0的函数 Φ(t)是t的函数 e
& φ (t , t0 ) = A(t )φ (t , t0 ) 即 将t=t0代入解中 x(t ) =φ(t ,t )x(t )
∴ x ( t ) = φ( t , t 0 ) x ( t 0 ) 是齐次矩阵的解
现代控制理论-状态方程的解
e
At
Te JtT1
1
0
2
0
et
0
2
3
1 1 4 0 0 e 2 t 1 2
Iinv(A)
2
1
1
作业：p.536;9-9;9-10;9-11
2tete2t 3tet2et2e2t tetete2t
0 1 0
2(e2ttetet) 3tet5et4e2t tet2et2e2tA
iPi1APi1Pi
..................................
参见 现代控制论
教材p.490 刘豹 p.28
Matla 中，矩阵求逆 b 命令为
1 1 0
et tet 0
J
T1AT=
0
1
0
0 0 2
e Jt
0
0
et
0
0 e 2 t
1 1 1 e t te t 0 2 5
T
2 2
1 2
特征 向量 问题
Api i pi
T 1
1
1 2
1 1
T1AT
1 0
0 2
e At
TetT1
2 2
1 et
2
0
2et e2t 2et 2e2t
e
0
2
t
1
1
1 2 1
et e2t et 2e2t
利用拉氏变换的方法参见书 p.458 例9-4，9-5
证明→
1
0
0
2
1
，
n
2
...
n
e At
矩阵指数
的求法
现代控制理论习题解答(第二章)
第二章 状态空间表达式的解3-2-1 试求下列矩阵A 对应的状态转移矩阵φ（t ）。
（1） ⎥⎦⎤⎢⎣⎡-=2010A（2） ⎥⎦⎤⎢⎣⎡-=0410A（3）⎥⎦⎤⎢⎣⎡--=2110A （4）⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-=452100010A（5）⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡=00100001000010A （6）⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡=λλλλ0100010000A【解】： （1）⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡++=⎥⎦⎤⎢⎣⎡+-=-=Φ-----)2(10)2(11}201{])[()(11111s s s s L s sL A sI L t⎥⎥⎦⎤⎢⎢⎣⎡-=⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡++-=---ttees s s s L 22105.05.01)2(10)2(5.05.01（2）⎥⎦⎤⎢⎣⎡-=⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡+++-+=⎥⎦⎤⎢⎣⎡-=-=Φ-----t tt ts s s s s s L s sL A sI L t 2cos 2sin 22sin 5.02cos 444414}41{])[()(222211111（3）⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡++-+++=⎥⎦⎤⎢⎣⎡+-=-=Φ-----222211111)1()1(1)1(1)1(2}211{])[()(s s s s s s L s s L A sI L t⎥⎥⎦⎤⎢⎢⎣⎡--+=Φ------tttttt teetete e te t )(（4）特征值为：2,1321===λλλ。
由习题3-1-7(3)得将A 阵化成约当标准型的变换阵P 为⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=421211101P ，⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡----=-1211321201P线性变换后的系统矩阵为：⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡==-20010011~1AP P A⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=t tt ttA e ete e e2~0000⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡----⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡===Φ-121132120000421211101)(21~t t tttA Ate te eePPeet⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡--++-----++-----++--=Φt t ttt tt t t t t t t tt tt t ttt t tt t t e te eete ee te e e te e e te e ete e ete e e te e tee t 34838424225342222322)(222222222（5）为结构四重根的约旦标准型。
现代控制理论_哈尔滨工程大学中国大学mooc课后章节答案期末考试题库2023年
现代控制理论_哈尔滨工程大学中国大学mooc课后章节答案期末考试题库2023年1.已知线性定常系统如下所示，下面说法错误的是（）【图片】参考答案:引入状态反馈后，不改变系统的能观测性。
2.串联组合系统的传递函数矩阵为各串联子系统的传递函数矩阵之和。
参考答案:错误3.在最优控制问题中，如果系统的性能指标是状态变量和控制变量的二次型函数，则称为线性二次型最优控制问题，简称LQ（Linear Quadratic）问题。
参考答案:错误4.用不大的控制能量，使系统输出尽可能保持在零值附近，这类问题称为输出调节器问题。
参考答案:正确5.研究系统控制的一个首要前提是建立系统的数学模型，线性系统的数学模型主要有两种形式，即时间域模型和频率域模型。
参考答案:正确6.现代控制理论以多变量线性系统和非线性系统作为研究对象，以时域法，特别是状态空间方法作为主要的研究方法。
参考答案:正确7.1892年俄国数学家李亚普诺夫发表了论文《运动稳定性的一般问题》，用严格的数学分析方法全面地论述了稳定性问题。
参考答案:正确8.经典控制理论以单变量线性定常系统作为主要的研究对象，以时域法作为研究控制系统动态特性的主要方法。
参考答案:错误9.下述描述中哪些作为现代控制理论形成的标志（）参考答案:用于系统的整个描述、分析和设计过程的状态空间方法._最优控制中的Pontriagin极大值原理和Bellman动态规划。
_随机系统理论中的Kalman 滤波技术。
10.内部稳定性表现为系统的零初态响应，即在初始状态恒为零时，系统的状态演变的趋势。
参考答案:错误11.系统矩阵A所有特征值均具有负实部是线性时不变系统渐近稳定的充要条件。
参考答案:正确12.从物理直观性看，能观测性研究系统内部状态“是否可由输入影响的问题”。
参考答案:错误13.由系统结构的规范分解所揭示，传递函数矩阵一般而言只是对系统结构的不完全描述，只能反映系统中的能控能观测部分.参考答案:正确14.下面论述正确的是（）参考答案:李亚普诺夫意义下渐近稳定等同于工程意义下稳定。
第六章线性系统的状态方程
状态变量分析法的优点：1. 便于观察系统内部某些物理量的变化过程；2. 与系统的复杂程度无关，复杂系统和简单系统的数学模型相似，适于多输入多输出系统;3. 适于研究非线性或时变系统。
因为一阶微分方程或差分方程是研究非线性和时变系统的有效方法。
4. 便于研究系统的稳定性、可控性、可观测性及系统内部参数变化对系统特性的影响；5. 状态方程都是一阶微分方程或差分方程，便于采用数值解法在计算机上实现系统分析。
系数矩阵由系统的参数决定，非时变系统为常数，时变系统为时间的函数。
,A B 四、输出方程(output equation))(,),(),(21t y t y t y r Λ输出方程是由状态变量和激励信号的线性方程，因此对线性系统而言，输出方程是一组线性方程。
例如，假设系统有个输出，r mrm r r n rn r r r mm n n mm n n e d e d e d x c x c x c t y e d e d e d x c x c x c t y e d e d e d x c x c x c t y +++++++=+++++++=+++++++=ΛΛMΛΛΛΛ22112211222212122221212121211112121111)()()(则,A B矩阵形式为：)(10081910120010321'3'2'1t e x x x x x x ⎥⎥⎦⎤⎢⎢⎣⎡+⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡⎥⎥⎦⎤⎢⎢⎣⎡---=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡01000112198⎡⎤⎢⎥∴=⎢⎥⎢⎥---⎣⎦A 依此方法选择的状态变量常称为相变量状态变量，状态方程叫相变量状态方程。
状态方程和输出方程中的系数矩阵与输入输出方程有关。
[]⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=3210410)(x x x t y 001⎡⎤⎢⎥=⎢⎥⎢⎥⎣⎦B []1040=C 0=D矩阵形式为：1211012110''13'22'1)()(+--+++=+----====m m n n n nn x b x b x b t y t e x a x a x a x x xx x x x ΛΛM )(1000100010211210''2'1t e x x x a a a a x x x n n n ⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡+⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡----=⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎣⎡-M M ΛM ΛΛM[]001111n n n n n nb b a b b a b b a b --∴=---=C D L 当时，矩阵不再为0。
3.3 线性时变连续系统状态方程的解 (PPTminimizer)
线性时变连续系统齐次状态方程的解(3/3) 线性时变连续系统齐次状态方程的解
证明 对解表达式x(t)=Φ(t,t0)x(t0)求导,则有 且
ɺ ɺ x (t ) = Φ(t , t0 ) x (t0 ) = A(t )Φ(t , t0 ) x (t0 ) = A(t ) x (t )
ɺ Φ(t , t0 ) = A(t )Φ (t , t0 ) Φ(t0 , t0 ) = I
x(t0)=Φ(t0,t0)x(t0)=x(t0) 说明式x(t)=Φ(t,t0)x(t0)满足齐次状态方程及其初始条件。 根据微分方程解的唯一性，所以它是齐次状态方程的解。 时变系统齐次状态方程的解表示了系统自由运动的特性,也代 表了初始状态x(t0)的转移,其转移特性完全由状态转移矩阵 Φ(t,t0)决定。
t0
τ1
然后按此法继续迭代下去,并将各展开式代入式(3-59),可 得
I + τ 1 A(τ )Φ (τ , t )dτ dτ Φ (t , t0 ) = I + ∫ A(τ 1 ) 2 2 0 2 1 1 )dτ 1 + ∫ A(τ 1 ) ∫
状态转移矩阵的求解(7/7) 状态转移矩阵的求解
上述A(t)和∫A(τ)dτ可交换条件一般较难以检验是否成立。 事实上,根据该可交换条件有
∫ [ A(t ) A(τ ) − A(τ ) A(t )] dτ ≡ 0
t t0
上式对于任意时间变量t和t0都成立的充分必要条件是:对 t t : 于任意的t1和t2,下式成立 A(t1)A(t2)=A(t2)A(t1) 所以,实际上较易于检验的条件可取代A(t)和∫A(τ)dτ可交 换条件,成为时变系统的状态转移矩阵的解可表示为指数 矩阵形式的充分必要条件。
现代控制理论复习
G ( s ) = G1 ( s )[ I + G2 ( s )G1 ( s )]−1 或 = [ I + G1 ( s )G2 ( s )]−1 G1 ( s )
2 0 1 1 1 2 9 8
[ 第二章总结] 1 .线性定常齐次状态方程的解 2 .矩阵指数函数 e At 3 .状态转移矩阵 Φ ( t − t0 ) Φ ( t , t ) 0 4 .线性定常非齐次状态方程的解 5 .线性时变系统状态方程的解 1 、线性定常系统运动分析 1 ）齐次状态方程的解：
2 0 1 1 1 2 9 2 2
4 、李氏第二法判稳 李氏第二法判稳思路：寻找李氏函数 李氏第二法判稳思路 李氏第二法稳定性定理
G11 G12 G G Y(s) G(s) = = C(sI − A)−1 B + D = 21 22 M M U(s) Gm1 Gm2 L G1r L G2r L M L Gmr
G( s ) 的每个元素的含义：
Yi ( s ) 表示第i 个输出中，由第j 个输入变量所引 Gij ( s ) = 个输入变量间的传递关系 U j ( s ) 起的输出和第j
e At
2 0 1 1 1 2 9
e λ1t = P 0
0 −1 O P e λnt
1 0
约当标准型法：当A 的特征值为 λ1（n 重根）
λ1t e = Q M 0 0 te λ1t L O O L 1 t n−1e λ1t ( n − 1)! −1 O M Q O te λ1t 0 e λ1t
x ( t ) = Φ ( t − t0 ) x ( t0 )
2 ）非齐次状态方程的解：
x( t ) = Φ ( t ) x( 0) + ∫ Φ ( t − τ ) Bu(τ )dτ
江苏大学线性系统理论(现代控制理论)考试必备--第2章
t0 , t 上为时间
为存在且ቤተ መጻሕፍቲ ባይዱ一。 x (t )
第2章 线性系统的运动分析
江苏大学电气学院
从数学的观点上看，上述条件可能显得过强而可减弱 为如下3个条件： （1）系统矩阵A(t)的各个元 aij (t )在时间区间 t0 , t 上为 绝对可积，即有：
t aij (t ) dt , i, j 1, 2,, n
性质9 e 的相似变换 如果矩阵A的特征值互不相同，并且存在非奇异变换 矩阵P，使得 A PAP 1，由 e At定义可得到：
At
e At Pe At P 1
第2章 线性系统的运动分析
At e 三. 几种典型的矩阵指数函数
江苏大学电气学院
由于矩阵指数函数 e At在计算线性系统的响应起着十分
第2章 线性系统的运动分析
江苏大学电气学院
二. 状态方程解的存在性和惟一性条件
当所选的状态变量不同时，所得状态方程不同，故状 态方程不是唯一的。对任意的初始状态，只有当线性系统 的状态方程的解存在且唯一时，对系统的分析才有意义。 从数学上看，这就要求状态方程中的系数矩阵和输入作用 满足一定的假设，它们是保证状态方程的解存在且唯一所
上式表明，条件（2）和（3）还可进一步合并为要求
B(t)u(t)的各元在时间区间上绝对可积。 对于连续时间线性定常系统，系数矩阵A和B为常数矩 阵且各元为有限值，条件（1）和（2）自然满足，存在惟
一性条件只归结为条件（3）。 在本章随后各节的讨论中，总是假定系统满足上述存 在性和惟一性的条件，并在这一前提下分析系统状态运动 的演化规律。
At e 性质7 积的关系式
d At e Ae At eAt A dt
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
2.4
线性时变系统状态方程的解
线性时变系统的状态方程为 & x = A ( t ) x + B( t )u （1）不一定有解 （2）当A（t)、B(t)在定义区间绝对可积时，对 每一初始态存在唯一解。 0 1 & 时变系统中A、B随t变化，如 x= x 0 t 一、齐次矩阵微分方程的解 & x(t ) = A(t ) x(t ) x ( t ) t = t = x ( t o )
& φ (t , t0 ) = A(t )φ (t , t0 ) 即 将t=t0代入解中 x(t ) =φ(t ,t )x(t )
∴ x ( t ) = φ( t , t 0 ) x ( t 0 ) 是齐次矩阵的解
∴ (t0,t0) =I φ
0
00
0
证毕
讨论： 1 齐次解与定常系统一样，也是初始状态 的转移,φ(t,t0)称为时变系统的状态转移 矩阵。 2 将定常系统中的φ(t),φ(t-t0)改为 φ(t,t0)定常系统可推广到时变系统。 二、非齐次矩阵微分方程的解 & x (t) = A (t)x (t) + B(t)u (t) 解：
x ( t ) = φ( t, t 0 ) x (t 0 ) +
∫
t
t0
φ( t , τ) B ( τ) u ( τ)dτ
证明（略）
三，状态转移矩阵 φ( t , t 0 ) 1. φ( t , t 0 ) 与 φ( t ) φ( t − t 0 ) 对比 共同：形式和性质类似 区别：本质不同，但Φ(t,t0)既是t的函数，也 是t0的函数 Φ(t)是t的函数 e
0
其解为： x ( t ) = φ( t, t 0 ) x ( t 0 )
φ( t , t 0 ) 是n*n阶非奇异方阵，且满足
& φ(t,t0) = Aφ(t,t0) φ(t0,t0) = I
证明：将解带入齐次方程
d [φ ( t , t 0 ) x ( t 0 ) ] = A ( t )φ ( t , t 0 ) x ( t 0 ) dt
At
Φ(t-t0)是（t-t0）的函数 e
A( t −t0 )
2.Ф（t,t0）的计算式－展开成级数，可用
数值方法计算
τ1 A(τ )dτ dτ + L φ（t , t0 ) = I + ∫ A(τ )dτ + ∫ A(τ 1 ) ∫ 2 2 t0 t0 t0 1
t t
ቤተ መጻሕፍቲ ባይዱ
证明：两边各项求导，可得 & φ（t , t0 ) = A(t )φ (t , t0 ) φ (t0 , t0 ) = I + 0 + 0 + L = I 满足齐次解中的条件