鲁棒

  • 1 鲁棒
  • 2 H ∞ H_\infty H∞​
  • 3 鲁棒 H ∞ H∞ H∞ 控制
  • 一、 H 2 H_2 H2​ 范数的计算
  • 二、 H ∞ H_\infty H∞​ 范数的计算
  • 绪论
    • 1. H 2 H_2 H2​ 控制
    • 2. 鲁棒控制
    • 3. H ∞ H_\infty H∞​ 控制
  • 第一章 数学基础知识
  • 第二章 线性二次型最优调节器
  • 第三章 参数化控制器设计
  • 第四章 H ∞ H_\infty H∞​ 标准控制
  • 第五章 时滞系统的 H ∞ H_\infty H∞​ 控制
  • 第六章 H 2 H_2 H2​ 标准控制
  • 第七章 H 2 / H ∞ H_2/H_\infty H2​/H∞​ 混合控制

首先要明白什么是鲁棒,什么是 H ∞ H_\infty H∞​。这两个概念是独立的。

1 鲁棒

鲁棒指系统在不确定参数(parameter uncertainties)摄动下,系统维持其自身某种性能(稳定、某指标)的特性。鲁棒性的主要问题是高频的未建模动态。

应该注意的是,我们在进行鲁棒稳定性分析时,要将干扰与控制输入置0。

2 H ∞ H_\infty H∞​

H ∞ H_\infty H∞​ 是对传递函数增益大小的一个度量指标,简单说就是一个系统输入输出的放大倍数。

H H H 无穷的标准问题实际是:求解已真有理(解析有界)的控制器 K K K,使外输入 W W W 到加权输入 Z Z Z 的传递函数阵的 H H H 无穷范数最小。翻译过来就是,由于我们一般考虑线性系统,直接运用闭环系统稳定性分析的 LMI 进一步公式推导,直至导出矩阵元素中含有某一参数(性能指标)的 LMI,运用 MATLAB 求解即可。

3 鲁棒 H ∞ H∞ H∞ 控制

鲁棒 H ∞ H∞ H∞ 控制就是抑制噪声到期望输出之间的传递函数集的最大增益,从而达到抗扰的目的。其实顾名思义,鲁棒 H H H 无穷就是研究系统既包含不确定性,又要求达到某一性能指标。鲁棒 H ∞ H∞ H∞ 控制具有较强的保守性。

总结:研究系统鲁棒性,就是研究包含不确定参数的系统。研究H无穷,就是研究满足某性能指标的系统。

From: 如何理解鲁棒H无穷控制理论?


现代鲁棒控制理论中,特别是多入多出系统(MIMO), H 2 H_2 H2​ 问题和 H ∞ H_\infty H∞​ 问题占据着中心地位,其中 H 2 H_2 H2​ 问题(也是一种特殊情况的 LQG问题 )研究的最多,最成熟,但是应用的却最少,相反却很大的促进了对线性系统的研究,也为 H ∞ H_\infty H∞​ 问题的提出和解决铺平了道路。

H 2 H_2 H2​ 问题以最小化输出能量为控制其设计目标,上世纪六七十年代,人们把主要的研究投入到这个问题,由 H 2 H_2 H2​ 问题引出的 H 2 H_2 H2​ 控制,优化目标中没有考虑扰动(噪声),把只针对白噪声或者某一个噪声的信号有效果,此外现实中广为使用的乘法不确定性也很难考虑,而且 H 2 H_2 H2​ 范数不是诱导范数,数学运算受到很大限制,虽然 H 2 H_2 H2​ 控制可以把目标值优化为0,但实际应用中,不确定噪声很大衰减了性能及 其稳定性。

在80年代,传递函数的 H ∞ H_\infty H∞​ 范数才被提出并被Doyle成功应用到鲁棒控制中, 线性系统的鲁棒控制得以完成, H ∞ H_\infty H∞​ 控制最小化传递函数(阵)的 H ∞ H_\infty H∞​ 范数,也是传递函数矩阵最大奇异值的最大值,与 H 2 H_2 H2​ 控制对比,它以最小化最大奇异值的最大值(也可以形象地说是传递函数的峰值),只要求扰动和不确定性有界,而且 H ∞ H_\infty H∞​ 范数是 L 2 L_2 L2​诱导范数,数学运算更加方便。九十年代只是针 对结构不确定性的 H ∞ H_\infty H∞​ 问题提出了结构奇异值综合法,并针对 H ∞ H_\infty H∞​ 问题提出了更好的算法,分别针对解析解( H ∞ H_\infty H∞​ 次优解等)和数值解 (LMI-线形矩阵不等式)两个方向发展,其中的LMI得以广泛的发展和推崇。

当然非线性系统的鲁棒控制, H ∞ H_\infty H∞​ 的思想也有很好的应用。

第一次看鲁棒控制感觉扑朔迷离,理论艰深,特别是本科刚毕业时,需要很多的知识需要补充,矩阵分析是最基本的,复变函数,线性泛函,线性系统理论都是必修课。这些完成之后就可以学习鲁棒控制了,个人认为 H 2 H_2 H2​ 问题虽然应用少,但是理解了 H 2 H_2 H2​, H ∞ H_\infty H∞​ 就顺理成章,有的牛人讲课时,只讲 H 2 H_2 H2​ 问题, H ∞ H_\infty H∞​ 一笔带过。

Ref: 鲁棒控制中的H2问题和H-inf问题-转


一、 H 2 H_2 H2​ 范数的计算

考虑线性时不变系统(多变量)
x ˙ = A x + B u y = C x \dot{x} = A x + B u \\ y = Cx x˙=Ax+Buy=Cx

求解下列 Lyapunov 方程
A L c + L c A T = − B B T A L_c + L_c A^T = -B B^T ALc​+Lc​AT=−BBT

代入
∥ G ( s ) ∥ 2 = t r a c e ( C L c C T ) \|G(s)\|_2 = \sqrt{trace( C L_c C^T )} ∥G(s)∥2​=trace(CLc​CT) ​

或者求解
A T L o + L o A = − C T C A^T L_o + L_o A = -C^T C ATLo​+Lo​A=−CTC

代入
∥ G ( s ) ∥ 2 = t r a c e ( B T L o B ) \|G(s)\|_2 = \sqrt{trace( B^T L_o B )} ∥G(s)∥2​=trace(BTLo​B) ​

Matlab 求解 Lyapunov 方程
A X + X A T + C = 0 X = l y a p ( A , C ) AX+XA^T+C=0\\ X = lyap(A,C) AX+XAT+C=0X=lyap(A,C)

对于单变量系统的传递函数
∥ G ( s ) ∥ 2 = 1 j 2 π ∮ G ( − s ) G ( s ) d s \|G(s)\|_2 = \sqrt{\frac{1}{j2\pi}\oint G(-s) G(s) ds} ∥G(s)∥2​=j2π1​∮G(−s)G(s)ds ​

等于 G ( − s ) G ( s ) G(-s)G(s) G(−s)G(s) 在左半平面上的留数之和的平方根。

例子:
G ( s ) = 1 s + 1 G(s) = \frac{1}{s+1} G(s)=s+11​

G ( s ) G ( − s ) G(s)G(-s) G(s)G(−s) 在左半平面上的极点为 s = − 1 s=-1 s=−1,那么在该极点上的留数为
lim ⁡ s → − 1 ( s + 1 ) 1 − s + 1 1 s + 1 = 1 2 \lim_{s\rightarrow -1} (s+1) \frac{1}{-s+1} \frac{1}{s+1} = \frac{1}{2} s→−1lim​(s+1)−s+11​s+11​=21​

所以
∥ G ( s ) ∥ 2 = 1 2 \|G(s)\|_2 = \frac{1}{\sqrt{2}} ∥G(s)∥2​=2 ​1​

function [gamma]=H2()
%{
程序功能:
1、使用迭代法计算传递函数G为矩阵时的H2范数。
2、求解ALc+LcA'+BB'=0,代入sqrt(trace(CLcC'))
3、求解A'Lo+LoA+C'C=0,代入sqrt(trace(B'LoB))
4、Matlab求解Lyapunov方程
AX+XA^T+C=0
X=lyap(A,C)
%}clear,clcA=[2,3 ;5,3];B=[3;5];C=[1, 0];D=0; %给出系统矩阵n=max(size(A));I=eye(2,2);Lc=lyap(A,B*B');gamma=sqrt(trace(C*Lc*C')) %方法1Lo=lyap(A',C'*C) ;gamma=sqrt(trace(B'*Lo*B)) %方法2end

二、 H ∞ H_\infty H∞​ 范数的计算

Ref: 鲁棒控制理论(十)计算系统H2范数和H∞范数


《 H 2 H_2 H2​ 和 H ∞ H_\infty H∞​ 优化控制理论》

绪论

1. H 2 H_2 H2​ 控制

最优控制理论中的线性二次型性能指标
J = ∫ 0 ∞ [ x T ( t ) Q x ( t ) + u T ( t ) R u ( t ) ] d t J = \int_0^\infty [x^T(t) Q x(t) + u^T(t) R u(t)]dt J=∫0∞​[xT(t)Qx(t)+uT(t)Ru(t)]dt

如果被控过程的数学模型用随机过程模型描述,且干扰信号是高斯白噪声(有限谱信号),则相应的线性二次型最优控制问题称为 LQG 问题,其性能指标为
J = lim ⁡ T → ∞ 1 T E { ∫ 0 ∞ [ x T ( t ) Q x ( t ) + u T ( t ) R u ( t ) ] d t } J = \lim_{T\rightarrow\infty} \frac{1}{T} E\{ \int_0^\infty [x^T(t) Q x(t) + u^T(t) R u(t)]dt \} J=T→∞lim​T1​E{∫0∞​[xT(t)Qx(t)+uT(t)Ru(t)]dt}

其中 E E E 为数学期望。

可以证明,上述 LQG 问题的性能指标等价于从随机干扰信号到系统被控输出信号的传函(阵)的 H 2 H_2 H2​ 范数。

上述 LQG 问题可以等价地表述为:设计反馈控制器 K K K,使闭环系统稳定,同时使从 w w w 到 z z z 的闭环传函阵 G z w ( s ) G_{zw}(s) Gzw​(s) 的 H 2 H_2 H2​ 范数达到极小,即
min ⁡ K ∥ G z w ( s ) ∥ 2 = γ 0 \min_{K} \| G_{zw}(s) \|_2 = \gamma_0 Kmin​∥Gzw​(s)∥2​=γ0​

称为 H 2 H_2 H2​ 最优控制问题。若给定 γ > γ 0 \gamma>\gamma_0 γ>γ0​,求反馈控制器 K K K,使
∥ G z w ( s ) ∥ 2 < γ \| G_{zw}(s) \|_2 < \gamma ∥Gzw​(s)∥2​<γ

则称为 H 2 H_2 H2​ 次优控制问题。统称为 H 2 H_2 H2​ 标准控制问题。

2. 鲁棒控制

一种典型的不确定系统模型的基本形式是
y = ( P + Δ ) u + n y = (P + \Delta) u + n y=(P+Δ)u+n

这里 y y y 是输出, u u u 是输入, P P P 是标称对象的传函,模型的不确定性以两种形式出现: n n n 是未知噪声或干扰, Δ \Delta Δ 是未知的对象摄动(未建模动态)。

一般地说,鲁棒性概念是指反馈控制系统的某项性能对于某集合中的每一个对象都是成立的。具体地说,我们有如下鲁棒性定义。

(1)鲁棒稳定性。设计一个控制器,使对每一个摄动后的对象,都能保证闭环系统的稳定性。

(2)鲁棒性能。设计一个控制器,使对每一个摄动后的对象,闭环系统都满足稳定性和某种特定的系统性能。

可见,上述鲁棒性概念是指给定一个控制器,如果某集合中的每一个对象都能保持某种特性成立,则称控制器对此特性是鲁棒的。因此,谈及鲁棒性必要求有一个控制器,有一个对象的集合某些系统特性

3. H ∞ H_\infty H∞​ 控制

考虑从 r r r 到 e e e 的误差传函(灵敏度函数)
S ( s ) = E ( s ) R ( s ) = 1 1 + P ( s ) K ( s ) S(s) = \frac{E(s)}{R(s)} = \frac{1}{1+ P(s)K(s)} S(s)=R(s)E(s)​=1+P(s)K(s)1​

我们的任务是设计控制器 K ( s ) K(s) K(s),使得
∣ S ( j ω ) ∣ < γ , ∀ ω |S(j\omega)| < \gamma, \forall \omega ∣S(jω)∣<γ,∀ω

其中 0 < γ ≤ 1 0<\gamma\le1 0<γ≤1。即将控制误差(偏差)控制在工程允许的范围内。上式两端对 ω \omega ω 取上确界,得到
sup ⁡ ω ∣ S ( j ω ) ∣ < γ \sup_{\omega} |S(j\omega)| < \gamma ωsup​∣S(jω)∣<γ

将上式左端定义为 S ( s ) S(s) S(s) 的 ∞ \infty ∞-范数(norm),即
∥ S ( s ) ∥ ∞ = sup ⁡ ω ∣ S ( j ω ) ∣ \|S(s)\|_\infty = \sup_\omega |S(j\omega)| ∥S(s)∥∞​=ωsup​∣S(jω)∣

H ∞ H_\infty H∞​ 最优控制问题:设计控制器 K ( s ) K(s) K(s),使闭环系统保持稳定,且使 ∥ S ( s ) ∥ ∞ \|S(s)\|_\infty ∥S(s)∥∞​ 达极小。即
inf ⁡ K { sup ⁡ ω ∣ S ( j ω ) ∣ } = inf ⁡ K ∥ S ( j ω ) ∥ ∞ = γ 0 \inf_{K} \{\sup_{\omega} |S(j\omega)|\} = \inf_{K} \|S(j\omega)\|_\infty = \gamma_0 Kinf​{ωsup​∣S(jω)∣}=Kinf​∥S(jω)∥∞​=γ0​

亦称为极大极小化问题。

H ∞ H_\infty H∞​ 次优控制问题:给定 γ > γ 0 \gamma>\gamma_0 γ>γ0​,设计控制器 K ( s ) K(s) K(s),使闭环系统保持稳定,且使
∥ S ( s ) ∥ ∞ < γ \|S(s)\|_{\infty} < \gamma ∥S(s)∥∞​<γ

与 H 2 H_2 H2​ 控制问题不同, H ∞ H_\infty H∞​ 最优控制问题难于求解,通常是以 H ∞ H_\infty H∞​ 次优控制问题的解去逼近 H ∞ H_\infty H∞​ 最优控制问题的解,并将 H ∞ H_\infty H∞​ 次优控制问题称为 H ∞ H_\infty H∞​ 标准控制问题。

对于多变量系统,传函阵 G z w G_{zw} Gzw​ 的 H ∞ H_\infty H∞​ 范数定义为
∥ G z w ∥ ∞ = sup ⁡ ω σ ˉ { G z w ( j ω ) } \|G_{zw}\|_\infty = \sup_\omega \bar{\sigma} \{ G_{zw}(j\omega) \} ∥Gzw​∥∞​=ωsup​σˉ{Gzw​(jω)}

其中 σ ˉ = λ max ⁡ ( G ∗ G ) 1 2 \bar{\sigma} = {\lambda_{\max} (G^* G) }^{\frac12} σˉ=λmax​(G∗G)21​,为 G G G 的最大奇异值, G ∗ G^* G∗ 为 G G G 的共轭转置阵, λ max ⁡ \lambda_{\max} λmax​ 为最大特征值。

在图1所示标准控制问题中,如果评价信号 z z z 为向量信号,并将其分块为 z = { z 1 T , z 2 T , ⋯ , z m T } T z=\{z_1^T, z_2^T, \cdots, z_m^T \}^T z={z1T​,z2T​,⋯,zmT​}T,记从 w w w 到 z i ( i = 1 , 2 , ⋯ , m ) z_i(i=1,2,\cdots,m) zi​(i=1,2,⋯,m) 的闭环传函阵为 G z i w ( s ) G_{z_i w}(s) Gzi​w​(s),取 G z i w ( s ) G_{z_i w}(s) Gzi​w​(s) 的 H 2 H_2 H2​ 范数或 H ∞ H_{\infty} H∞​ 范数,就得到相应的 H 2 H_2 H2​ 或 H ∞ H_{\infty} H∞​ 多目标控制问题。

如果将 z z z 分块为 z = { z 1 T , z 2 T } T z=\{z_1^T, z_2^T\}^T z={z1T​,z2T​}T,且取
∥ G z 1 w ( s ) ∥ 2 < γ 1 , ∥ G z 2 w ( s ) ∥ ∞ < γ 2 \|G_{z_1w}(s)\|_2 < \gamma_1, \|G_{z_2w}(s)\|_\infty < \gamma_2 ∥Gz1​w​(s)∥2​<γ1​,∥Gz2​w​(s)∥∞​<γ2​

则相应的标准控制问题称为 H 2 / H ∞ H_2/H_\infty H2​/H∞​ 混合控制问题。

第一章 数学基础知识

第二章 线性二次型最优调节器

第三章 参数化控制器设计

第四章 H ∞ H_\infty H∞​ 标准控制

第五章 时滞系统的 H ∞ H_\infty H∞​ 控制

第六章 H 2 H_2 H2​ 标准控制

第七章 H 2 / H ∞ H_2/H_\infty H2​/H∞​ 混合控制

【控制】鲁棒性 H2 H无穷相关推荐

  1. H2/H∞半车悬架控制仿真分析

    目录 前言 1.悬架模型 2.LMI求解 3.simulink仿真分析 3.1结论 前言 对于H2/H∞控制的鲁棒项相比不用多说,之前也写过两篇关于1/4车的H2/H∞控制文章,链接如下: 基于LMI ...

  2. H无穷控制理论与应用案例分析

    0.知识背景 概念.定义 内稳定: BIBO稳定: 镇定:对于一个控制系统来说,如果通过某种反馈可以使系统实现渐近稳定,即闭环系统极点具有负实部,则称该系统是能镇定的. 信号的H2范数: 信号的H无穷 ...

  3. matlab使用LMI对二阶系统进行H无穷控制

    MATLAB使用LMI对二阶系统进行H无穷控制 1. 二阶系统建立 2. 二阶系统实现 3. MATLAB代码 本文需要对LMI理论和H无穷控制理论具有一定的了解. 1. 二阶系统建立 建立一个二阶系 ...

  4. 【控制】《多智能体系统一致性与复杂网络同步控制》郭凌老师-第7章-Lurie 型动态网络的鲁棒 H 无穷簇同步

    第6章 回到目录 第8章 第7章-Lurie 型动态网络的鲁棒 H 无穷簇同步 7.1 引言 7.2 问题描述 7.3 主要结果 7.3.1 无时滞 Lurie 型子系统 7.3.2 时滞 Lurie ...

  5. 1.2 引入不定度量空间的动机——有限时域H2/H∞滤波

    在本节中,通过回顾H2/H∞状态空间理论的一些结果来引入Krein空间(不定度量空间) 目录 1.2.1 状态空间估计 H2 方法:随机扰动过程 问题1.2.1(H2估计问题) 定理1.2.1(卡尔曼 ...

  6. H无穷滤波(H∞ filtering )

    从今年年初开始,由于项目的需要接触到H无穷滤波,无奈中文互联网上相关内容太少,这部分知识点又过于复杂.学的相当痛苦,花费了大量时间还是一知半解.后来找到一本系统讲述H无穷滤波的书(Indefinite ...

  7. matlab中h无穷状态反馈控制,求教 :用LMI方法求倒立摆H无穷状态反馈控制器程序...

    求教 :用LMI方法求倒立摆H无穷状态反馈控制器程序用LMI方法求倒立摆H无穷状态反馈控制器程序 但运行出现没有合适的gama值 显示"could not establish feasibi ...

  8. 【MATLAB教程案例69】基于MATLAB编程实现H无穷控制器与性能分析

    欢迎订阅<FPGA学习入门100例教程>.<MATLAB学习入门100例教程> 目录 1.软件版本 2.H无穷控制器简介 3.MATLAB编程实现H无穷控制器

  9. H无穷控制学习笔记——系统性能分析

    一.连续时间系统 1.系统增益指标 (1)系统性能指标定义 对于线性时不变的连续时间系统:x˙(t)=Ax(t)+Bw(t)z(t)=Cx(t)+Dw(t)\dot{x}(t)=Ax(t)+Bw(t) ...

最新文章

  1. opencv reshape函数详解
  2. 如何提升你的能力?给年轻程序员的几条建议
  3. hook 驱动 截屏
  4. RHEL 5.4 安装Oracle 11gR2, 准备篇...
  5. HarmonyOS之公共事件的发布、订阅与退订
  6. Netty实战 IM即时通讯系统(五)客户端启动流程
  7. java.io.IOException: No FileSystem for scheme: hdfs转载加自己笔记
  8. 动态表单实现客户端二次过滤及字段汇总统计
  9. jooq sql_用jOOQ用Java编写SQL
  10. 三星电子时隔近3年再次成为全球最大半导体厂商
  11. 致敬那些用键盘改变世界的劳动者
  12. 什么是JavaSE/JavaEE/JavaME?
  13. linux下开发51单片机
  14. 图像压缩-《Learned Image Compression with Discretized Gaussian Mixture Likelihoods and Attention Modules》
  15. 【FFmpeg+Qt开发】转码流程 H.264 转(mov、mp4、avi、flv)等视频格式 示例详解
  16. win10系统做游戏服务器,Win10专业版如何提升游戏流畅度?Win10游戏流畅度的三种提升方法...
  17. Flask框架——flask-caching缓存
  18. shell判断所输整数是否为质数
  19. 关注疫情调查实践报告
  20. xbox360游戏下载_如何将游戏文件移动到另一个Xbox

热门文章

  1. python的flask框架
  2. setTimeout 实现 setInterval
  3. 队列 (Queue)
  4. 《嵌入式 - 嵌入式大杂烩》 处理器架构与指令集
  5. RevitAPI: 修改视图View裁剪区域Cropbox的大小
  6. 使用64位任务管理器转储的32位.NET进程dump,WinDbg如何调试
  7. 【HackTheBox】 meow
  8. 词袋模型和词向量模型
  9. cae属于计算机技术吗,CAE、CAD、CAPP与CAM的区别——扫盲贴
  10. 最全最详细的PHP面试题(带有答案)