现代控制理论第六章分析

合集下载

现代控制理论 6-1 概念 6-2 李雅普诺夫第一法(间接法)

现代控制理论 6-1 概念  6-2 李雅普诺夫第一法(间接法)
无差别
渐近稳定
收敛至 平衡状态
y 一致稳定
对定常系统
tc 与初始时刻
无差别
李雅普诺夫稳定 (稳定)
无关
前页 返回
cae 渐近稳定 tcy 小球
李雅普诺夫 意义下稳定
/ 稳定
cn = 2 x − xe = (x1 − x1e )2 + (x2 ) − x2e 2 = c 表示状态空间中,以xe为圆心,半径为c的圆
y n = 3 x − xe = (x1 − x1e )2 + (x2 ) − x2e 2 + (x3 ) − x3e 2 = c
tc 表示状态空间中,以xe为圆心,半径为c的球
前页
返回
12
例:⎩⎨⎧xx&&12
= =
x2 −x1
平衡状态
xe
=
⎡0⎤ ⎢⎣0⎥⎦
cae tcy前页
返回
设系统初始状态位于以平衡状态 xe 为球心,δ 为半径的闭球域 S(δ)内,即
ex0 − xe ≤ δ t = t0
a若系统的平衡状态 xe不仅具有李雅普诺夫意义
下的稳定性,且有
c lim t→∞
定,不一定大范围渐进稳定。
δ → ∞ S(δ ) → ∞
x2
x0
xe x0
x1
前页 返回
例:机械位移系统
aex(t), x&(t) cm
k
tcy前页 返回 18
内部稳定/状态稳定
初始状态 任意
大范围一致渐近稳定
大范围渐近稳定
e 对线性系统 无差别
对定常系统 无差别
对线性系统 无差别
a一致渐近稳定 c 对定常系统

江苏大学线性系统理论(现代控制理论)考试必备--第6章.答案

江苏大学线性系统理论(现代控制理论)考试必备--第6章.答案

=
C R
P1
CP1
RP
1
I qq 0
0 I ( n q )( n q )
再来讨论(n-q)维状态观测器的构建,用线性变换 x = Px,
将方程(1)变换成
x = PAP-1x + PBu y = CP-1x = CP-1x = Iqq 0 x
记 : A=PAP-1 B=PB
C CP1
以足够快的速度趋近于零,也就是说,不管状态观测器的
初始状态如何,状态观测器所重构的状态变量 xˆ 终将逐渐
趋近于实际状态 x ,所以,这样的状态观测器也称之为渐 进状态观测器。该性质也使其在实际使用中毋需设置初始 状态。
第6章 状态观测器
江苏大学电气学院
值得一提的是,虽然 (A-MC) 特征值的负实部离虚
i (A C M ) i , i =1,2, , n
求出M后,即可构成闭环状态观测器:
xˆ = (A - MC)xˆ + My + Bu
(8)
第6章 状态观测器
江苏大学电气学院
全维状态观测器的另一种设计方法是,先对被观测系
统进行非奇异变换 z=T,x 再从形式上列出类似于式(8)
的观测器方程。
B
x
x C
y
A
xˆ 0
B

xˆ C

A
第6章 状态观测器
江苏大学电气学院
这样的观测器称为开环状态观测器,从开环状态观测
器中取出 xˆ 可作为 x 的估计值近似替代,当然希望 xˆ 与x 是相等的。用 x 来表示 x 和 xˆ 的偏差,即 x x xˆ , 下面来简单分析估计偏差 x的特性。式(1)和式(2)相减得

现代控制理论基础第六章书上第三章(1)PPT课件

现代控制理论基础第六章书上第三章(1)PPT课件
是系统输出对输入的稳定性
两种稳定性既有区别,又有内在的联系
2
⑶ 本章内容

稳定性:内部稳定性与外部稳定性 本章重点是内部稳定性
•李雅普诺夫稳定性理论和方法
适用范围:线性系统、非线性系统和离散系统 常用的判据:李雅普诺夫函数法稳定性判据
李雅普诺夫方程稳定性判据
3
3.1 线性系统的外部稳定性
线性系统的外部稳定性或零状态响应的稳定性,是对应于系 统输入输出描述的稳定性 。是有界输入有界输出稳定性,简 称为BIBO 稳定性。
g (s)的一个极点2.5与零点对消,剩下一个负实极点 -1,所以系 统是 BIBO稳定的。
10
3.2 系统的内部稳定性
系统的内部稳定性是研究系统的零输入响应的稳定性。因
此只要讨论齐次状态方程
x f( x ,t)
x ( t0 ) x 0 ,t t0
(3-4)
由初始状态 x(t0)x0引起的响应的稳定性,是状态稳定性问题。
•对渐近稳定系统, A 总是非奇异的,零状态(原点)是系统的
唯一平衡状态。
12
例3-2 倒立摆系统
系统的齐次状态方程为
y(t)Cx(t)Du(t)
则系统的传递函数阵为
G (s ) C (s I A ) 1 B D 1 C a(s d I-A j)B ds e I tA )(
G (s)的极点必是 A的特征值。
(3-3)
如果 A的所有特征值具有负实部,则G (s)的所有极点必定具 有负实部,则系统是 BIBO稳定的。
4
3.1.1 单变量线性系统的 BIBO稳定性判据
⑴ 脉冲响应函数判据
定理3-1 线性系统的输入输出描述是
y(t)tt0g(t,)u(t)d

现代控制理论基础 第6章 线性系统的最优控制

现代控制理论基础 第6章 线性系统的最优控制

7
方法的比较
总的来说,当控制量无约束时,‘采用“变分法” ;当控制量有 约束时,采用“极小值原理” 或“动态规划”;如果系统是线性的, 采用“线性二次型”方法最好,因为,一方面,二次型指标反映了大 量实际的工程性能指标的要求;另方面,理论上的分析及求解较简单、 方便、规范,而且还有标准的计算机程序可供使用;得到的控制器易 于通过状态反馈实现闭环最优控制,工程实现方便。在实际的工程控 制中,目前线性二次型最优控制己得到了广泛的成功应用。
J 值为极值 J (最大值或最小值),这种泛函求极值的方法,实际上 就是数学上的“变分”问题,须采用数学中的“变分法” 。
5
采用直接变分法求解最优控制率,难于甚至“无法解决容许控 制属于闭集”的最优控制问题,所以受到实际工程应用上的限制, 例如,每台电动机都有最大功率的限制;船舶或飞机的操纵舵面 也有最大偏转角的限制。况且采用直接变分法设计出的系统,其 抗参数变化的能力,即系统的鲁棒性也不强。因此,工程应用上 有较小的实用价值。
线性系统二次型的最化控制,因为其性能指标具有明确的物理 意义,在大量的工程实际中具有代表性,而且最优控制率的求解 较简单,并具有统一的解析表达式,构成的最优控制系统具有简 单的线性状态反馈的型式,易于工程实现,所以在国内外实际的 工程中目前己得到广泛应用。本章主要介绍其基本概念、基本原 理和设计方法。
下面只介绍线性二次型最优控制的基本概念、求解原理及设 计中的一些主要结论。
8
第三节 线性二次型最优控制
一、控制对象数学模型
线性系统的状态空间表达式
x(t) A(t)x(t) B(t)u(t)
y(t) C(t)x(t)
式中,
n x(t) 为 维状态向量;
(6-4)

现代控制理论课件-第六章 极小值原理

现代控制理论课件-第六章  极小值原理

⑴ 满足正则方程
x*
k
1
H
x*
k
,u* k ,* k 1
k
1,k
f x* k ,u* k ,k
*
k
H
x*
k
,u* k ,* xk
k
1,k
⑵ 相对于最优控制,哈密尔顿函数达极小值,即
H x* k ,u* k ,* k 1,k H x* k ,uk ,* k 1,k
⑶ 及满足以下边界条件及横截条件
x*
0
x0,*
N
x* N ,N x N
同理,对不同的边界情况,只需选取相应的边界条 件及横截条件,条件1、2不变。当控制变量不受限 制时,则条件2与控制方程
等效。
H
x*
k ,u* k ,* uk
k
1,kபைடு நூலகம்
0
§ 6.3 极小值原理解最短时间控制问题
一般情况下,非线性受控系统的最短时间控制问题的 解析解是很困难的,本节只讨论线性定常受控系统的 最短时间控制问题。
比较上述极小值原理与变分法所得的结果,可以发现 两者的差别仅在⑵。 极小值原理的严格证明很复杂,下面的证明将重于物 理概念的阐述,尽量避免烦琐的数学推导。 设系统动态方程为:
xt f xt,ut,t
边界条件为:xt0 x0 ,为简单起见,假设终端时刻 t f
及终端状态 x t f 均为自由。控制变量 ut 受有界闭集 约束,即 utU
性能指标为:
J x
tf
,t f
tf t0
F
xt,ut,t dt
则使性能指标 J 达到极小的最优控制 u* t 及最优状态 轨线 x* t 必须满足以下条件:

第六章 最优控制(2) 现代控制理论

第六章 最优控制(2) 现代控制理论

x1(t)
x10
x20t
1t2 2
消去时间变量 t , 可得相应的最优轨线方程为
x1(t)
1 2
x22
(t)
C
(6-256)
在图6-16中用实线表示。
14
由于x2(t)=x20+t随t增大, 故最优轨线行进的方向自 下而上, 如曲线上箭头所示。
15
当 u= -1 时, 状态方程的解为
x2 (t) x20 t
在R-上 在+上
在R+上 到达原点
u 1,1 u 1 u 1, 1 u 0
19
进一步, 可综合为
u 1 当(x1, x2 ) R u 1 当(x1, x2 ) R
u 0 当(x1, x2 ) 0
若将开关曲线方程写成
h (x1, x2 )
x1
1 2
x2
x2
0
则最优控制律可表示成
x2 (t) u(t)
或写成矩阵形式
x(t)
0 0
1 0
x(t)
10u(t
)
初始条件 x(t0) x0
(6-248)
终端条件 x(t f ) 0
控制约束 1 u(t) 1, (t0 t t f )
性能指标
J
t f
t0
1 dt
求 最 优 控 制 u*(t) , 把 系 统 从 初 态 转 移 到 终 态 , 使
x1(t)
x10
x20t
1 2t2Fra bibliotek相应的最优轨线方程为
x1(t)
1 2
x22
(t)
C
在图6-16中用虚线表示。由于x2(t)随t减小, 故 曲线箭头方向自上而下。

现代控制理论 第6章


x&
(
A
BHC
)
x
Bv
1 1
2 3
x
0 1
v
其能控性矩阵和能观性矩阵的秩分别为
rank[ B
(
A
BHC
)
B
]
rank
0 1
2 3
定理 3:输出至 x&的反馈不改变系统的能观性但可能改变原系统的能
控性。
6.1线性反馈控制系统的基本结构
例1
设线性定常系统的状态空间模型为:
x&
1 3
2
1
x
0 1
u
y
1
2 x
并设状态反馈阵 K=[3 1] 和输出反馈 H=2。
试分析该系统的状态反馈闭环系统和输出反馈闭环系统的状态能控/能观性。
第6章 控制系统的状态空间综合
线性定常系统的综合
控制系统的分析和综合是研究控制系统的两大问题。
线性定常系统分析:在建立的数学模型的基础上分析系统的各种性能。 如: 能控性、能观性、稳定性等和定量运动规律分析 如 :系统运动轨迹、系统的性能品质指标等。
线性定常系统的综合
系统综合是系统分析的逆问题。 系统综合首先需要确定关于系统运动形式,或关于系统运动动态过程和 目标的某些特征的性能指标函数,然后据此确定控制规律。
线性定常系统综合: 给定被控对象,通过设计控制器的结构和参数,使系统满足性能指标要求。
分为: 常规综合和最优综合。
线性定常系统的综合
常规综合的性能指标是一类由不等式及等式约束的性能指标凸空间, 一般只要求解的控制规律对应的性能指标到达该凸空间即可。 而对非优化型性能指标一般存在解析方法求解控制规律,如极点配置方法。

现代控制理论基础课件第六章书上第三章(1)


⑴ 系统的平衡状态
对定常系统,齐次状态方程为
A(t ) x x x (t0 ) x0 , t t0
(3-5) (3-7)
如果系统所处的状态
xe
满足
e 0 x
这个状态称为平衡状态。 由平衡状态的定义,x 不会使系统产生运动,即
e
φ(t; t0 , xe , 0) xe
(3-8)
例 3-1 设系统的状态空间描述为
1 0 1 x x u 0 2.5 0 y 1 1 x
的特征值为 -1与2.5,不全为负实部。而其传递函数为
g (s) c(sI A) 1 b (s 2.5) 1 (s 1)(s 2.5) (s 1)
10
例3-3 如图所示的单摆, 当取状态变量为 x1 , x2 ,状态方程
0 g sin x l 1 x 0
l

图3-1
m
这是一个非线性系统,对其在 处进行线性化,可得线性化 方程
0
0 g cos 0 x l
1 x 0
i
i
1 1 1 , , , s pi (s pi ) 2 ( s p i ) mi
它们的拉氏反变换,或系统的单位脉冲响应相应地包含有下 列因子
e pi t , te pi t , , t m 1 e pi t
上列因子绝对可积的充分必要条件是 p 具有负实部,即系统 是 BIBO稳定的。 证毕
y(t ) tt0 g (t , )u(t )d
(3-1) (3-2)
则系统是 BIBO 稳定的充分必要条件是 式中, M 是一个有限常数。
证明 充分性:由式(3-1),有

现代控制理论第六章


式中,δx(t) 为宗量函数x(t)的变分, L[x(t), δx(t)] 是 δx(t) 的线性连续泛函,o[ x(t), δx(t)] 是关于 δx(t) 的高阶无穷 小,则定义泛函增量的线性主部
δJ = L[ x(t), δ x(t)]
(6-19)
为泛函 J[ x(t)] 的变分,记作 δJ 。若泛函有变分,则 称该泛函可微。
物体的升降速度,则上式可写成状态方程
& x1 (t) = x2 (t)
& x2 (t) = u(t) − mg
x 其初始条件是 x1 (t0 ) = x10 , 2 (t0 ) = x20 。现需寻找 一个能使物体以最短时间从初态 ( x10,x20 ) 到达终态 (0,0)的控制u(t)。定义系统的性能指标为
1. 始端时刻和终端时刻固定时的泛函极值问题
首先讨论不仅初始时刻 t0 、终端时刻 t f 固定,而 且初始状态 x(t 0 ) = x0 、终端状态 x(tf ) = xf固定这一最 简单情况下无约束条件的泛函极值问题(最优控制的 最优控制的 基本问题)。 基本问题
J = ∫ dt = t f − t0
tf t0
t 式中, t0为起始时刻, f 为终止时刻。要求时间最短, 即使性能指标J最小,这样求得的控制即为最优控制 u *(t) 。
2. 搅拌槽问题 设有一盛放液体的连续 搅拌槽,如图6-2所示。槽内 装有不停转动着的搅拌器S, 使液体经常处于完全混合状 态,槽中原放 0o C 的液体。 现需将其温度升高,为此在 入口处送进一定量的液体, 其温度为u(t),出口处流出 等量的液体,以保持槽内液
由式(6-20)得
∂ (J[x(t) + εδx(t)]) = ∂ ∫tt0f [x(t) + εδx(t)]2 dt ∂ε ∂ε ε =0

现代控制理论 6-3 李雅普诺夫第二法(直接法)

c et c y前页返回2 2c aet c y()xV&22xμ−=返回前页求出系统的李雅普诺夫第二法的基本思想ce tcy 1x 返回前页定理3渐近稳定cae tcy ()00≠=但x V&返回前页定理3⎪⎩⎪⎨⎧−−==21221x m x m k x x xμ&&⎥⎦⎤⎢⎣⎡=00e x 渐近稳定et c y返回前页1 xca e tcy ⎪⎩⎪⎨⎧−==1221x m k x x x &&⎥⎦⎤⎢⎣⎡=00e x 返回前页定理4李雅普诺夫意义下稳定cet c y1x返回前页定理3不稳定ca e tcy ()00≡=但x V&返回前页定理3⎪⎩⎪⎨⎧+−==21221x m x m k x x x μ&&⎥⎦⎤⎢⎣⎡=00e x 不稳定状态平面图状态仿真曲线注意tcy 前页返回前图?李氏函数选择不当!cet c y返回前页定理3et c y返回前页e虚构atcae tcy ()=V x ()02221>+=x x V x ()0 ≡x V &()0222≤−=x V x &ec ayt c etcy 返回前页定理4cae tcy ⎥⎦⎤⎢⎣⎡=0n πe x ⎪⎩⎪⎨⎧−== x L gx x x1221sin &&状态仿真曲线李雅普诺夫意义下稳定返回前页tcy 0≡返回前页定理3cae tcy 状态平面图状态仿真曲线()00≡=但x V&⎪⎩⎪⎨⎧−== x L g x x x1221sin &&2Dx −()L ,,,nn πe 2100±±=⎥⎦⎤⎢⎣⎡=x 垂直向下渐近稳定前页返回cae tcy 相平面图θL。

  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

系统
1 2 0 & x u : x 3 1 1
y [1 2]x
完全能控能观,引入反馈
u [3 1]x V
则闭环系统 K的状态空间表达式为
1 2 0 & K : x x v 0 0 1
y [1 2]x
1
1
1
定理 6.1.1 对于任何实常量矩阵 K,系统 K
完全能控的充要条件是系统 完全能控。 证 注意到系统 和 K 的能控性矩阵分别为
uc [B AB A2 B
An1B] ( A BK )n1 B]
uc ' [ B ( A BK ) B ( A BK )2 B
A TAT
1
Ac 0
A12 Ac
bc b Tb 0
且对任意 k [k1 , k 2 ] ,有
det(sI A bk) det(sI A bkT 1 ) det(sI A bk )
例6.2.1
0 0 0 1 & y 0 x 1 1 0 x 0 u 0 1 1 0
%:

& % % % % x Ax Bu % % % y Cx Du
1 a n 1 1 a1
P 为非奇异的实常量等价变换矩阵,且有 这里,
0 ~ A PAP1 a n
0 M ~ cP1 % Pb c b n n1 0 1 对式(6.2.2)引入状态反馈
(s 2)( s 1 j 3 )( s 1 j 3 ) s 3 4 s 2 8s 8
* a1 4, a2 8, a3 8.
k a3 a3 , a2 a2 , a1 a1 8,7,2
3)


4)
a2 a1 1 1 0 0 1 2 1 1 2 1 2 Q b Ab A b a1 1 0 0 1 1 2 1 0 1 1 0 1 0 0 0 0 1 1 0 0 1 0 0
1 式(6.3.2)可写为 y(s) G(s)u(s) C (sI A) Bu (s)
y1 ( s ) g11 ( s )u1 ( s ) g12 ( s )u2 ( s ) L L g1 p ( s )u p ( s ) y2 ( s ) g 21 ( s )u1 ( s ) g 22 ( s )u2 ( s ) L L g 2 p ( s )u p ( s ) M M yq ( s ) g q1 ( s )u1 ( s ) g q 2 ( s )u2 ( s ) L L g qp ( s )u p ( s )


比较两多项式同次幂的系数,有 :
2 k1 4 , 2 k1 k 2 1 8 , k1 k 2 k 3 8
得:k1 2, k2 3,k 3 3 即得状态反馈增益矩阵为:
k 2 3 3
与例6.2.1的结果相同 6.2.3 讨论 (1) 状态反馈不改变系统的维数,但是闭环传递 函数的阶次可能会降低,这是由分子分母的 公因子被对消所致。
% a * a 由此即有 k 1 n n
% * k2 an1 an1 M * % kn a1 a1
% % v Kx % u v Kx v KP1 x
又因为
所以
% K KP
必要性:采用反证法,设 不完全能控,则必
非奇异变换阵 T 使系统结构分解
6.3 应用状态反馈实现解耦控制
6.3.1 问题的提出 考虑MIMO系统
& Ax Bu :x y Cx
(6.3.1)
在 x(0) 0 的条件下,输出与输入之间的关系, 可用传递函数 G ( s) 描述:
y(s) G(s)u(s) C(sI A)1 Bu(s)
(6.3.2)
% % u v Kx
1
~ d d
% k % L % k K 1 2
% k n
% 的状态空间表达式为 则闭环系统 K
% % & % bK % % % x (A )x bv % : K % % % % % y (c dK )x dv
其中,显然有
1) 由
0 0 s s 3 2s 2 s det( sI A) det 1 s 1 0 0 1 s 1
a1 2, a2 1, a3 0.
* (s 1 ) (s *2 ) ( s * 3)
得 2) 由 得

K HC ,则 Kx Hy ,状态反馈就等价于输
出反馈 H 。
(2) D=0时,可以求得闭环系统 K 的传递函数阵
G(s;K,L) C[sI ( A BK )]1 BL
①利用矩阵运算直接可推出(见书)
G(s;K,L) G(s)[I K (sI A)1 B]1 L
第六章 状态反馈和状态观测器
6.1 状态反馈的定义及其性质 6.2 极点配置 6.3 应用状态反馈实现解耦控制
6.4 状态观测器
6.1 状态反馈的定义及其性质
给定系统
& Ax Bu :x y Cx Du
在系统中引入反馈控制律
u Lv Kx
则闭环系统 K 的结构如图 6.1.1 所示。
D
v
L
-
u
B
+

x
+
C
y
A K
图 6.1.1
K 的状态空间表达式为:
& K :x ( A BK ) x BLv y (C DK ) x DLv
若 D 0 ,则
& K :x ( A BK ) x BLv y Cx
状态反馈性质 (1) L I 时,为单纯的状态变量反馈。若


5)
P Q 1 1 2 1 1 1 0 1 0 0
1
0 0 1
0 1 2
1 1 1
6)
1 0 0 ~ 2 2 1 算法2:直接配置
n
同时,由指定的任意 n 个期望闭环极点*1 , * 2 ,, * n
可求得期望的闭环特征方程
(s )(s 2 )(s
* 1 * * n
) s a s
n
* n 1 1
a
* n 1
s an 0
*
通过比较系数,可知 ~ a1 k n a1* ~ * a 2 k n 1 a 2 ~ * a k a 1 n n
②在图 6.1.1 中令D 0 并改用图6.1.2 表示
a
v
L
-
u
B
+

x
I
b
C
y
A K
图 6.1.2
图中a和 b 之间的部分,可以看成是由系统
& x Ax Bu
% y Ix ( 为单位矩阵)
和输出反馈
% u v Ky
Gab (s) 不难用 所组成从到 b 的传递函数矩阵。
* 1 2 * 2, 3 1 j 3
给定系统的状态空间表达式为
1 1x
求状态反馈增益阵 K ,使反馈后闭环特征值为
1 0 0 3 A 2 b rank 0 1 1 0 0 1
解:因为
rank b

Ab

系统是状态完全能控,通过状态反馈控制律 能 任意 配置闭环特征值。

( A BK ) B 的列向量可以由 ( B AB)
( A BK ) B AB B( KB) ,可知
的列向量的线性组合表示。
( A BK ) B 的列向量可以由( B AB A B )的
2
2
列向量的线性组合表示。 依此类推,不难看出
[ B ( A BK ) B ( A BK )2 B ( A BK )n1 B] A B] 的列向量
输出反馈传递函数阵的公式求出,
Gab (s) (sI A)1 B[ I K (sI A)1 B]1
于是,从 即为
v
到 y 的传递函数矩阵 G(s; K , L)
G( s;K,L) C (sI A) B[ I K (sI A) B] L G(s)[ I K (sI A)1 B]1 L
n * n 1 * * f ( x) (s 1 )(s ) ( s ) s a s a s a 2 n 1 n 1 n
3) 列方程组 ai (k ) a , i 1,, n. 并求解 。
* i
其解 k [k1 ,, k n ] ,即为所求 例6.2.2 同例6.2.1。
s 3 2 k1 s 2 2k1 k 2 1s k1 k 2 k 3
根据要求的闭环期望极点,可求得闭环期望特征 多项式为
f s s 2 s 1 j 3 s 1 j 3 s 3 4s 2 8s 8

1 0 % % bK % (A ) % % an k1 an 1 k2 O L 1 % a1 k n
% 的闭环特征方程为 系统 K
~ n1 ~ ~ n 2 s (a1 kn )s (a2 kn1 )s (an k1 ) 0
解:设所需的状态反馈增益矩阵k为 k k1 k 2 k3 因为经过状态反馈 u v kx 后,闭环系统 的
相关文档
最新文档