1、Ch.3 线性系统的时域分析,目录(1/1),目 录概述 3.1 线性定常连续系统状态方程的解 3.2 状态转移矩阵及其计算 3.3 线性时变连续系统状态方程的解 3.4 线性定常连续系统的离散化 3.5 线性定常离散系统状态方程的解 3.6 Matlab问题 本章小结,线性时变连续系统状态方程的解(1/2),3.3 线性时变连续系统状态方程的解严格说来,实际控制对象都是时变系统,其系统结构或参数随时间变化。 如电机的温升导致电阻以及系统的数学模型变化;电子器件的老化使其特性也发生变化; 火箭燃料的消耗导致其质量以及运动方程的参数的变化等。 但是,由于时变系统的数学模型较复杂,且不易于系统分析
2、、优化和控制,因此只要实际工程允许,都可将慢时变系统在一定范围内近似地作为定常系统处理。 但对控制目标要求较高的高精度控制系统,需作为时变系统处理。,线性时变连续系统状态方程的解(2/2),下面将讨论线性时变连续系统状态方程的求解问题,依次讨论: 线性时变连续系统齐次状态方程的解 线性时变连续系统的状态转移矩阵 非齐次状态方程的解,线性时变连续系统齐次状态方程的解(1/3),3.3.1 线性时变连续系统齐次状态方程的解当系统没有外部输入作用时,线性时变连续系统的状态方程为齐次状态方程,可表示为 x(t)=A(t) x(t) 这里讨论其满足初始状态,的解,也就是由初始时刻t0的初始状态x(t0)
3、所引起的无输入强迫项(无外力)时的自由运动。 为保证该齐次状态方程解的存在性和唯一性,假设在系统的时间定义域t0,tf内, A(t)的各元素为时间t的分段连续函数。,线性时变连续系统齐次状态方程的解(2/3),下面证明时变系统齐次状态方程的解为 x(t)=(t,t0) x(t0) 式中,(t,t0)为时变系统的状态转移矩阵,它定义为如下矩阵微分方程的解。,对比 线性定常系统(t)=A(t), (t)|t=0=I,线性时变连续系统齐次状态方程的解(3/3),证明 对解表达式x(t)=(t,t0)x(t0)求导,则有且 x(t0)=(t0,t0)x(t0)=x(t0) 说明式x(t)=(t,t0)
4、x(t0)满足齐次状态方程及其初始条件。 根据微分方程解的唯一性，所以它是齐次状态方程的解。时变系统齐次状态方程的解表示了系统自由运动的特性,也代表了初始状态x(t0)的转移,其转移特性完全由状态转移矩阵(t,t0)决定。,x(t)= (t,t0) x(t0),x(t)=(t)x0=(t-t0)x(t0),对比线性定常系统,线性时变连续系统的状态转移矩阵(1/1),3.3.2 线性时变连续系统的状态转移矩阵下面进一步讨论前面引入的状态转移矩阵,主要内容为: 状态转移矩阵 的 求解 状态转移矩阵 的 性质,状态转移矩阵的求解(1/7),1. 状态转移矩阵的求解对于线性时变连续系统,状态转移矩阵(
5、t,t0)是如下矩阵微分方程和初始条件(t)=A(t) (t), (t)|t=0=I 的解,它是一个nn维的关于时间变量t和t0的矩阵函数。(t)=A (t), (t)|t=0=I 为了求得状态转移矩阵(t,t0)的表达式,可在时间域内对该矩阵微分方程积分,即有,对比线性定常系统,状态转移矩阵的求解(2/7),如果将上式中积分号内的(1,t0)再按上式展开,则有然后按此法继续迭代下去,并将各展开式代入式(3-59),可得,状态转移矩阵的求解(3/7),于是,可得一个由无穷项之和组成的状态转移矩阵(t,t0),即上式就是线性时变连续系统的状态转移矩阵的计算公式。在一般情况下,它不能写成封闭的解析
6、形式。 在实际应用此公式时,可按一定的精度要求,用数值积分计算方法去近似计算t1时刻的(t1,t0)的值。,状态转移矩阵的求解(4/7),当时变的系统矩阵A(t)满足如下条件(矩阵乘法可交换条件)时,时变系统的状态转移矩阵的解可以表示为的指数形式。 也就是说,只有 A(t) 与 A()d 满足矩阵乘法的可交换条件时,上述指数表达形式的解才成立。 下面对这个条件给予证明。,状态转移矩阵的求解(5/7),将该指数表达形式的右边展开成级数形式,有如果上式是系统的状态转移矩阵,它必须满足状态转移矩阵的定义式。 于是,将上式的两边对时间取导数,根据状态转移矩阵的解表达式,状态转移矩阵(t,t0)的导数可
7、表示为,状态转移矩阵的求解(6/7),比较上述两式可知,只有A(t)和A()d满足乘法可交换条件时,时变系统的状态转移矩阵可以表示为指数形式。 因此, A(t)满足矩阵乘法可交换条件下，线性 时变 连续 系统 齐次状态方程 的解 也可表示为指数形式, 即,状态转移矩阵的求解(7/7),上述A(t)和A()d可交换条件一般较难以检验是否成立。 事实上,根据该可交换条件有上式对于任意时间变量t和t0都成立的充分必要条件是:对于任意的t1和t2,下式成立 A(t1) A(t2)=A(t2) A(t1) 所以,实际上较易于检验的条件可取代A(t)和A()d可交换条件,成为时变系统的状态转移矩阵的解可表
8、示为指数矩阵形式的充分必要条件。,状态转移矩阵的性质(1/8),2. 状态转移矩阵的性质时变系统的状态转移矩阵的性质如下。 1) (t,t)=I 对比 (0)=eA0=I2) 传递性 (t2,t1) (t1,t0)=(t2,t0)对比 (t2-t1)(t1-t0)=(t2-t0),状态转移矩阵的性质(2/8),证明: 由于 x(t2)=(t2,t0)x(t0) 且 x(t2)=(t2,t1)x(t1)=(t2,t1)(t1,t0)x(t0) 故有 (t2,t0)x(t0)=(t2,t1)(t1,t0)x(t0) 由于上式对任意初始状态x(t0)都成立,所以有 (t2,t0)=(t2,t1)(t
9、1,t0),状态转移矩阵的性质(3/8),3) 可逆性 -1(t, t0)=(t0, t) 对比 (t2-t1)-1=(t1-t2) 证明 由性质1)和2),有 (t,t0)(t0,t)=(t,t)=I 故-1(t,t0)=(t0,t)成立。,状态转移矩阵的性质(4/8),4) 对角线矩阵的状态转移矩阵。 如果时变的系统矩阵A(t)如下表示的对角线矩阵。 A(t)=diaga11(t) a22(t) ann(t) 式中,aii(t)(i=1,2,n)为标量函数,则A(t)的状态转移矩阵(t,t0)为如下对角线矩阵： (t,t0)=diag11(t,t0) 22(t,t0) nn(t,t0) 式
10、中,ii(t,t0)(i=1,2,n)为满足如下标量微分方程的状态转移函数即,状态转移矩阵的性质(5/8),5) 块对角矩阵的状态转移矩阵。 如果时变的系统矩阵A(t)如下表示的块对角矩阵。 A(t)=block-diagA1(t) A2(t) Al(t) 式中,Ai(t)(i=1,2,l)为mimi维的分块矩阵函数,则A(t)的状态转移矩阵(t,t0)为如下块对角矩阵： (t,t0)=block-diag1(t,t0) 2(t,t0) l(t,t0) 式中,i(t,t0)(i=1,2,l)为满足如下矩阵微分方程的状态转移矩阵,状态转移矩阵的性质(6/8)例3-9,例3-9 求如下时变系统的状
11、态转移矩阵(t,t0)。解 首先检验矩阵A(t)和A()d与是否可交换。 为此计算,状态转移矩阵的性质(7/8)例3-9,因此 A(t1)A(t2)=A(t2)A(t1) t1，t2 即矩阵A(t)和A()d与满足可交换条件,可由指数展开式方法计算状态转移矩阵,即,状态转移矩阵的性质(8/8)例3-9,由于而于是,非齐次状态方程的解(1/8),3.3.3 非齐次状态方程的解当具有外加输入作用时,其状态方程为如下非齐次状态方程: x(t) = A(t) x(t) + B(t) u(t)该状态方程在初始状态为,下的解, 也就是由初始状态x(t0)和输入作用u(t)所引起的系统状态的运动轨迹。,非齐
12、次状态方程的解(2/8),下面将证明当输入u(t)为分段连续时,该非齐次状态方程的解为：证明 先设该非齐次状态方程的解为显然,有式中,(t)为待定函数。,非齐次状态方程的解(3/8),将所设的解代入该状态方程的左边,有将所设的解代入该非齐次状态方程的右边,有因此有即,非齐次状态方程的解(4/8),对上式两端积分,可得故该非齐次状态方程的解为当系统的状态空间模型中输出方程为 y(t)=C(t)x(t)+D(t)u(t) 时,系统的输出为：,非齐次状态方程的解(5/8),比较线性定常连续系统与线性时变连续系统状态方程的解的表示形式: 定常系统,时变系统,初始状态的影响,初始时刻后输入的影响,为脉冲
13、响应函数与输入的卷积,非齐次状态方程的解(6/8),与线性定常连续系统的状态方程和输出方程的解比较可知,线性时变连续系统与线性定常连续系统的解的结构和形式相同,都为状态的零输入响应和零状态响应的和。 线性定常连续系统的状态方程和输出方程的解可视为线性时变连续系统相应的解的一种特殊形式。 在A(t)为时不变A时,时变系统的状态转移矩阵(t,t0)即变为定常系统的状态转移矩阵(t-t0)。 由此可以看出引入状态转移矩阵的重要性。 只有引入状态转移矩阵,才能使时变系统和定常系统的求解公式建立统一的形式。,非齐次状态方程的解(7/8),例3-10 求如下时变系统在阶跃输入时的状态变量的值。解 由例3-9有,非齐次状态方程的解(8/8),由时变系统的状态方程的解表达式,有,
