第二章线性定常连续系统状态方程的解
2_状态方程求解

4) (t2 t1 ) (t1 t0 ) (t2 t0 )
e At e Bt e( A B )t 5)当且仅当AB=BA时
状态转移矩阵的基本性质
1) (t ) A (t ) (t ) A
2) (0) I
3)
(t )1 1 (t ) (t )
5 线性系统的脉冲响应矩阵
6 线性连续系统方程的离散化
7 线性离散系统的运动分析 8 用MATLAB求解系统方程
2.1
线性定常系统齐次状态方程的解
x (t ) Ax (t )
(1)
线性定常系统齐次状态方程为 这时系统的输入为零 先考察标量齐次微分方程的幂级数解法
x ax
假设其解为一幂级数
如果 [ sI A]为非奇异
[ sI A] x ( s) x (0)
1
x(s) [sI A]1 x(0) [sI A]1 x(0)
(10)
x (t ) L
1
{[sI A]1 x(0)} L
1
由微分方程解的唯一性
(t ) e At L
[sI A]1
而 b0 x(0) (4)
则解为 x(t ) (1 at 因为
1 2 2 1 a t a k t k ) x(0) e at x(0) 2! k!
1 2 2 1 a t akt k 2! k!
e at 1 at
模仿标量齐次微分方程的解法,假设线性定常系统齐次状态方程 (1)的解为 (5) x b0 b1t b2t 2 b3t 3 bk t k 将(5)式代入(1)式
e
Jt
现代控制理论-状态方程的解

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
现代控制理论第二章

= α n −1 (t ) An −1 + α n − 2 (t ) An − 2 + ⋯ + α1 (t ) A + α 0 (t ) I
【例2-5】见板书
(3)α i (t )的计算公式 A的特征值互异时 α 0 (t ) 1 λ1 α1 (t ) 1 λ2 ⋮ = ⋮ ⋮ α (t ) 1 λ n −1 n
λ λ λ
பைடு நூலகம்
2 1 2 2
⋮
2 n
⋯ λ e λ1t λ2 t ⋯ λ e ⋮ ⋮ λn t n −1 ⋯ λn e
At
2.变换A为约旦标准型 (1)A特征根互异 Λ = T −1 AT 有
例2-2 ,同例2-1
e At = Te ΛtT −1
(2)A特征值有重根
J = T AT e At = Te JtT −1
0 1 0 [例2 - 3]已知A = 0 0 1 , 求e At 2 - 5 4
若
σ ω A= −ω σ
则
cos ωt sin ωt σt e = Φ(t ) = e − sin ωt cos ωt
At
2.2.4 计算
1.根据 e At 或 Φ (t ) 的定义直接计算
1 2 2 1 33 1 n n e = I + At + A t + A t ⋯ A t + ⋯ 2! 3! k! 1 0 [例2 - 1]已知A = , 求e At − 2 − 3
线性控制理论 第2章 状态空间表达式的求解

12t 2 0 2 2 2 t 1 2! 0 2 2 n t
1 2 2 1 t t 0 1 1 2! 1 2 2 1 2 t 2 t 2! 1 2 2 0 1 n t n t 2!
1
1 2 1 m 1 t t 2! (m 1)! t (2-21) 1 2 1 t 2! t 1 mm
证明 因
12 1 1 0 1 2 ,A A 0 1 1 1 mm 21
x(t ) Φ(t ) x(0),t 0
上式表明齐次状态方程的解,在初始状 态确定情况下,由状态转移矩阵惟一确定,
即状态转移矩阵包含了系统自由运动的全部
信息,完全表征了系统的动态特性。
定义2.1
线性定常系统状态转移矩阵 Φ(t t0 ) 是
满足矩阵微分方程和初始条件
(t t ) AΦ (t t ), t t Φ 0 0 0 Φ (t0 t0 ) I
(2-3)
(t ) b1 2b2t kbk t x
( k 1)
k
Ax (t ) A(b0 b1t b2t bk t )
2ቤተ መጻሕፍቲ ባይዱ
比较上式两边t的同次幂可得
现代控制系统课件第2章

2021/1/4
5
2.2 矩阵指数函数——状态转移矩阵 2.2.1 状态转移矩阵
齐次微分方程的自由解为: x(t) eAt x0
或
x(t) e A(tt0 ) x0
从这个解的表达式可知,初始时刻的状态矢量x0, 到任意t>0或t>t0时刻的状态矢量x(t)的一种矢量变换 关系,变换矩阵就是矩阵指数函数 eAt 。
2021/1/4
27
例
x1 x2
0 2
1 3
x1 x2
0 1u
求 u(t) 为单位阶跃函数时,系统状态方程的解 (设
初始状态为零).
解:
(t)
e At
2et 2et
e2t 2e2t
et e2t
et
2e2t
x(t) e At x(0) 0te A(t ) Bu( )d
0 1
例:已知 A 2 3 ,求eAt
解: s 1
sI A 2 s 3
2021/1/4
19
(sI A) 1
1 sI A
adj (sI
A)
(s
1 1)(s
2)
s 3
2
1 s
s3
(
s
1)( s 2
2)
(s 1)( s 2)
1
(s
1)( s s
2)
(s 1)( s 2)
x(t) eAt x0 , t 0
2021/1/4
2
证明: 和标量微分方程求解类似,先假设式齐次状 态方程的解x(t)为t的矢量幂级数形式,即:
x(t) 0 1t 2t 2 iti
**
代入齐次状态方程中, 得
1 22t iit i1
连续系统的状态变量方程求解

连续系统的状态变量方程求解连续系统的状态变量方程求解通常采用数值方法,例如龙格-库塔法(Runge-Kutta)等。
在这个过程中,需要将连续系统的状态方程离散化,即将连续时间步长的微分方程转化为离散时间步长的离散方程。
求解离散方程可采用递推的方式,根据系统的初始条件和上一时刻的状态变量值,计算出当前时刻的状态变量值。
以下是一个求解连续系统状态变量方程的步骤:1. 确定连续系统的状态变量方程。
例如,给定线性定常系统dx/dt = Ax + Bu,其中x为状态变量,A和B为系统矩阵。
2. 离散化。
将状态变量方程转化为离散方程。
常见的离散化方法有前项差分变换、后项差分变换和Tustin变换。
具体变换方法取决于系统的特性以及所需的数值稳定性和精度。
例如,使用Tustin变换将连续系统离散化,得到离散状态方程x[k+1] = A*x[k] + B*u[k]。
3. 初始化。
给定初始条件,如x[0] 和u[0],初始化状态变量值。
4. 数值求解。
使用数值方法(如龙格-库塔法)递推计算离散方程,得到一系列状态变量值x[1], x[2], ...,以及对应的输出值y[1], y[2], ...。
5. 分析结果。
根据求解得到的状态变量值和输出值,分析系统的性能,如稳定性、收敛速度等。
在MATLAB中,可以使用ode45等函数求解连续系统的状态变量方程。
以下是一个简单的示例:```MATLAB定义系统矩阵A、B和输入信号uA = [1 0; -1 1];B = [0 1];u = [1; 0.5];定义初始条件x0 = [1; 2];设置求解参数tspan = [0, 10];options = odeset('RelTol', 1e-6, 'AbsTol', 1e-6);求解状态变量方程[x, u] = ode45(@(t, x) A*x + B*u, tspan, x0, options);绘制状态变量曲线figure;plot(t, x(:, 1), 'b', 'LineWidth', 2);hold on;plot(t, x(:, 2), 'r', 'LineWidth', 2);xlabel('Time');ylabel('State Variables');legend('x1', 'x2');```这个示例中,我们使用ode45函数求解了一个线性定常系统在给定输入信号下的状态变量演化。
现代控制理论第二章

第二章 控制系统状态空间表达式的解建立了控制系统状态空间表达式之后,就是讨论求解的问题,本章重点讨论状态转移矩阵的定义,性质和计算方法,从而导出状态方程的求解公式并讨论连续时间系统状态方程的离散化的问题。
§2-1线性定常齐次状态方程的解(自由解)所谓自由解是指系统输入为零时,由初始状态引起的自由运动。
状态方程为齐次矩阵微分方程:AX X= (2-1)若初始时刻0t 时的状态给定为00)(x t x =,则式(2-1)有唯一确定解。
0)(0)(x e t x t t A -=,0t t ≥(2-2)若初始时刻从0=t 开始,即0)0(x x =,则其解为:0)(x e t x At =, 0t t ≥(2-3)证:先假设式(2-1)的解)(t x 为t 的矢量幂级数形式,即:+++++=k k t b t b t b b t x 2210)((2-4)对上式求导: ++++=-1232132)(k k t kb t b t b b t x代人式(2-1)得:A x= ( +++++kk t b t b t b b 2210) (2-5)既然式(2-4)是(2-1)的解,则式(2-5)对任意时刻t 都成立,故t 的同次幂项的系数应相等,有:01Ab b =,0212!2121b A Ab b ==,0323!3131b A Ab b ==,… 01!11b A k Ab kb k k k ==-,… 在式(2-4)中,令0=t ,可得:00)0(x x b == 将以上结果代人式(2-4),故得:022)!1!211()(x t A k t A At t x k k +++++= (2-6)括号内的展开式是n n ⨯矩阵,它是一个矩阵指数函数,记为At e221112!!At k ke At A t A t K =+++++ (2-7)式(2-6)可表示为:0()At x t e x =再用)(0t t -代替)0(-t ,即在代替t 的情况下,同样证明0)(0)(x e t x t t A -=的正确性。
第2章 状态空间表达式求解

1 T 2. 若A能通过非奇异变换予以对角线化,即 AT
则
e1t e At (t ) T 0
e2t
0 T 1 n t e
证明:根据定义式
A2t 2 A3t 3 Ak t k e I At 2! 3! k 0 k! At
A2t 2 A3t 3 ( I At ) A e At A 信息与控制工程学院 2! 3!
5. 性质五
设有nxn矩阵A和B,当且仅当AB=BA 时,有eAteBt
= e(A+B)t ,而当AB≠BA 时,则eAteBt ≠ e(A+B)t 。
证明:根据定义式
e ( A B ) t ( A B ) 2 2 ( A B )3 3 I ( A B )t t t 2! 3! A2t 2 ABt2 BAt2 B 2t 2 I ( A B )t ( ) 2! 2! 2! 2! A3t 3 A2 Bt3 ABAt3 AB2t 3 BA2t 3 BABt3 ( 3! 3! 3! 3! 3! 3! B 2 At3 B 3t 3 ) 3! 3!
2 2 1 t 2! 1 1t 1 k k 2t At e At k 0 k! nt 1 0 0 k k 1 t k! k 0 0 2 2 2t 2!
(t )( ) (t ) (t )( t ) (t t ) I ( )(t ) ( t )
( t )(t ) ( t t ) I
从而证明了(t)与(-t)互为逆
信息与控制工程学院
4. 性质四
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
• 定义:一个系统,当且仅当对于任何两个容许 对 1 1 1 1
x (t ), u [t , ) x [t , ), y [t , ) x (t ), u [t , ) x [t , ), y [t , )
0 0 0 0 2 2 2 2 0 0 0 0
2 所构成的输入---状态--和任何实数1 和 输出对.
0, 0 0[t0 , ), 0[t0 , )
• 从而,如果系统是线性系统的话,则必有
当
x(t0 ) 0, u[t0 , ) 0
时
系统响应亦为零—这也是线性系统的一 个必要条件.
(2) 式(1)中, 若 1 2 1 称式(1)的关
系为可加性。 若 2 0 则称式(1)的关系为齐次性。
(3) 式(1)中,若设 1 2 1 ,及假定
x1 (t0 ) x(t0 ), u1 (0), x 2 (t0 ) 0, u 2 [t0 , ) u[t0 , )
则 x(t0 ) 0, 0 u[t0 , ) xt [t0 , ), yt [t0 , ) 或 x(t0 ), u[t0 , ) x(t0 ), 0 0, u[t0 , ) xt [t0 , ), yt [t0 , )
k k A t 1 22 At At e 或e I At A t 2! k 0 k !
1 k k At k!
并对于有限时间是绝对收敛的.
பைடு நூலகம்
• 结论:[零输入响应]线性定常连续系统的零
输入响应,即系统齐次方程的解,并具有如下
形式:
xou t e x 0
为初始状态的自由运动即零输入响应.
(2).零输入响应的形态.
• 对线性定常连续系统,零输入响应即自由运 动轨迹的形态,当且仅当由系统的矩阵指数函 数 e At 唯一地决定.不同的系统矩阵A,导致不 同形态的零输入响应,即自由运动轨道. 表明 e At 即A系统矩阵,包含了零输入响应即 自由运动形态的全部信息.
也是容许的,则称该系统是线性的,否则
该系统是非线性的.简而言之,满足迭加
原理的系统为线性系统.
2.对定义的讨论 (1)若设 1 2 1 并有
x (t0 ) x (t0 ), u [t0 , ) u [t0 , )
1 2 1 2
• 则如果是线性系统的话,按定义, • 则 x(t0 ) 0, u[t0 , ) 0 .
At
性质:
(1).e A0 I , t 0 (2).e e e e e At 1 At (3).(e ) e (4).若A,F R nn , 且AF FA, 则e( A F )t e At e Ft e Ft e At 若AF FA, 则e( A F )t e At e Ft
现代控制理论
第二章 线性定常连续系 统状态方程的解
线性定常连续系统状态方程的解
• 准备知识A1 • 准备知识A2 • 一.线性定常连续系统齐次方程的解 (零输入响应) • 二.状态转移矩阵 • 三.线性定常系统非齐次方程的解
• 准备知识A1
1.利用状态和状态方程来定义系统的线性性 质. • 用符号 u[t0 , ), x(t0 ) x[t0 , ), y[t0 , ) • 表示状态 x(t0 ) 和输入 u[t0 , ) 激励出输出 y (t ) 和状态 x(t ),t t0 ,并称其为输入-状态-输出对.
At
t0
• 推论: (1).零输入响应的运动特性. • 对于线性定常连续系统,其零输入响应 是由其齐次方程解的属性决定的,状态
空间中x(t)随时间演化轨道(几何表征),
属于由偏离系统平衡状态的初始状态 x0
引起的自由运动.
•一个典型的例子是:人造卫星在末级火箭 脱落后的运行轨道,以脱落时刻的运行状态
x(t ) e
A ( t t0 )
x (t0 )
t t0
(6).零输入响应的几何表征.
• 对线性定常连续系统,齐次方程解的表达 式表明:在时刻 ti 状态点 x(ti ) ,几何上对应 于状态空间中由初始状态点 x0 ,经线性变换
e x0 导出一个变换点.基于此,可推知,零输 入响应x(t ) 随时间t的演化过程,几何上即为 状态空间中由初始状态点出发和由各个时 刻变换点构成的一条轨迹.
时,系统的响应.
• 注意: 数学上,零状态响应 xox (t ) 即为零初始
状态下的强迫方程
x Ax Bu, x(t0 ) 0, t [t0 , )
的状态解. • 物理上,零状态响应 xox (t ) 代表系统状态由
输入u所激励的强迫运动
准备知识A2
• 不加证明地给出以下定理和定义.
(3)定理2.每一个基本矩阵 ,对(-∞,∞)中所有的t
而言,是非奇的.
(4)定义2.设 () 是 x A(t ) x(t ) 的任一基本矩
阵,对所有(-∞,∞)中的 (t , t0 )
1 ( t , t ) ( t ) (t0 ) 称 0
是 x A(t ) x(t ) 的状态转移矩阵.
• 从而 x(t ) L1[( sI A) 1 ]x(0)
2 I A A 由于 ( sI A)1 2 3 s s s
1 22 所以 L [( sI A) ] I At A t 2!
• 显然 b0 z (0) ,从而方程的解 z (t ) 可写为
1 22 1 kk z (t ) (1 at a t a t 2! k!
• 其中指数函数 1 22 at e at a t 2!
1 k k at k!
z (0) e z (0)
一.线性定常连续系统齐次方程的解(零 输入响应)
1.讨论
显然 x Ax 是矩阵微分方程,在解该方程之
前先观察纯量微分方程 z az 的解,其中 dz z dt
• 在解 z
az 时,先假定解
2
z (t ) b0 b1t b2t
• 代入方程得到
b1 2b2t 3b3t 2 a (b0 b1t b2t 2
零输入响应和零状态响应,而系统的全部响 应,则是它们的和. • 根据线性系统的性质:若
x(t0 ), 0 x1[t0 , ), y1[t0 , ) 0, u[t0 , ) x2 [t0 , ), y2 [t0 , )
• 则 x(t0 ), u[t0 , ) x(t0 ), 0 0, u[t0 , )
At At A ( t ) A A
d (5). (e At ) Ae At e At A dt d (6). (e At ) Ae At e At A dt (7).(e At ) m e A ( mt ) , m 0,1, 2,
3.齐次方程的拉普拉斯解法. • 同样先考虑纯量微分方程 z az 将方程两端作拉氏变换 sz ( s ) z (0) az ( s )
x Ax, x(t0 ) x0 , t [t0 , )
• 物理上,零输入响应代表系统状态的自由运动,特
点是响应形态只由系统矩阵所决定,不受外部输
入的影响.
• 定义: [零状态响应]:
线性系统的零状态响应 xox (t ) 定义为只有输
入作用,即 u (t ) 0 而无初始状态作用,即 x0 0
则 z ( s) ( s a) z (0)
1
• 将这种方法推广到矩阵微分方程的解
对 x Ax, x R nn , A R nn
两边取拉氏变换,则有
sx s x 0 Ax s 或 sI A x s x 0
即
x( s ) ( sI A) 1 x(0)
(4).零输入响应的计算.
• 根据解,则零输入响应计算的核心是计算矩 阵指数函数 e At。
(5).零输入响应表达式的更一般形式.
• •
对线性定常连续系统,通常习惯地取初始 时间 t0 0 。 由于线性时不变系统的分析只与相对时间 有关,这种处理也不失一般性.但若因某种需 要,将初始时间取为 t0 0 ,此时,零输入响 应更有一般的形式:
at
1 k k a t k 0 k !
• 仿上述纯量微分方程的解法,对于矩阵微分
方程 x Ax, x(0) x0 , t 0 其中 x R n1 , A R nn,
则
x e x(0), e R
At At
nn
称 e At 是按矩阵A定义的矩阵指数函数,并可 证明,若A是nn 的方阵时,则有:
(1)
1 2 1 2 x ( t ) x ( t ), u [ t , ) u [t0 , ) 1 0 2 0 1 0 2
1 2 1 2 x [ t , ) x [ t , ), y [ t , ) y [t0 , ) 1 0 2 0 1 0 2
(1)定理1. x A(t ) x (t ) 的全体解的集合,形成在实 数域上的n维向量空间.
(2)定义1. m n 矩阵函数 中,当且仅当n个列分 别是 x Ax 的n个线性无关解时,称 为 x A(t ) x (t ) 的基本矩阵,即 A ,且 (t0 ) 非奇.
Ati
2.解的性质(矩阵指数函数的性质) • 矩阵指数函数 e At 在线性系统分析中具有重
要意义,为此可基于定义,给出e At 的性质.
1 22 1 k k e I At A t A t 2! k! 1 k k A t A R nn , e At R nn k 0 k !