鲁棒控制与鲁棒控制器设计87205课件.ppt

1、2022-6-71控制系统计算机辅助设计MATLAB语言与应用东北大学信息学院第7 章鲁棒控制与鲁棒控制器设计薛定宇薛定宇 著著控制系统计算机辅助设计控制系统计算机辅助设计 MATLAB语言与应用语言与应用 清华大学出版社清华大学出版社2022-6-72控制系统计算机辅助设计MATLAB语言与应用东北大学信息学院主要内容 线性二次型线性二次型 Gauss 控制控制 鲁棒控制问题的一般描述鲁棒控制问题的一般描述 鲁棒控制器的计算机辅助设计鲁棒控制器的计算机辅助设计 新鲁棒控制工具箱及应用新鲁棒控制工具箱及应用 分数阶控制系统分析与设计分数阶控制系统分析与设计2022-6-73控制系统计算机辅助设

2、计MATLAB语言与应用东北大学信息学院7.1 线性二次型 Gauss 控制7.1.1 线性二次型 Gauss 问题假设对象模型的状态方程表示为 为白噪声信号，分别表示模型的不确定性与输出信号的量测噪声。2022-6-74控制系统计算机辅助设计MATLAB语言与应用东北大学信息学院定义最优控制的指标函数为2022-6-75控制系统计算机辅助设计MATLAB语言与应用东北大学信息学院7.1.2 使用 MATLAB 求解 LQG 问题 带有 Kalman 滤波器的 LQG 结构2022-6-76控制系统计算机辅助设计MATLAB语言与应用东北大学信息学院Kalman 滤波器的增益矩阵式中 , 满足

3、下面的 Riccati 代数方程2022-6-77控制系统计算机辅助设计MATLAB语言与应用东北大学信息学院【例7-1】2022-6-78控制系统计算机辅助设计MATLAB语言与应用东北大学信息学院 LQG 控制器设计的分离原理2022-6-79控制系统计算机辅助设计MATLAB语言与应用东北大学信息学院 基于观测器的 LQG 调节器设计2022-6-710控制系统计算机辅助设计MATLAB语言与应用东北大学信息学院由 Kalman 滤波器方程可以写出基于观测器的 LQG 调节器为2022-6-711控制系统计算机辅助设计MATLAB语言与应用东北大学信息学院2022-6-712控制系统计算

4、机辅助设计MATLAB语言与应用东北大学信息学院【例7-2】2022-6-713控制系统计算机辅助设计MATLAB语言与应用东北大学信息学院7.1.3 带有回路传输恢复的 LQG 控制 LQG/LTR 控制器设计算法使用 LQG 控制器，系统的开环传递函数表示为直接状态反馈系统的开环传递函数为2022-6-714控制系统计算机辅助设计MATLAB语言与应用东北大学信息学院【例7-3】2022-6-715控制系统计算机辅助设计MATLAB语言与应用东北大学信息学院2022-6-716控制系统计算机辅助设计MATLAB语言与应用东北大学信息学院 回路传输恢复技术回路传输恢复技术(loop tran

5、sfer recovery，LTR)加权函数的选择2022-6-717控制系统计算机辅助设计MATLAB语言与应用东北大学信息学院 2022-6-718控制系统计算机辅助设计MATLAB语言与应用东北大学信息学院2022-6-719控制系统计算机辅助设计MATLAB语言与应用东北大学信息学院【例7-4】 对【例7-3】不同的 q 值应用 LTR 技术2022-6-720控制系统计算机辅助设计MATLAB语言与应用东北大学信息学院2022-6-721控制系统计算机辅助设计MATLAB语言与应用东北大学信息学院 应用 MATLAB 求解 LQG/LTR 问题若想使得系统在输入端恢复环路传递函数，则

6、若想在对象模型的输出端恢复环路传递函数，则2022-6-722控制系统计算机辅助设计MATLAB语言与应用东北大学信息学院【例7-5】 对【例7-3】选定一个 q 向量，设计 LTR 控制器，并绘制出不同 q 值下环路传递函数 的 Nyquist 图。2022-6-723控制系统计算机辅助设计MATLAB语言与应用东北大学信息学院7.2鲁棒控制问题的一般描述 小增益定理 鲁棒控制器的结构 鲁棒控制系统的 MATLAB 描述2022-6-724控制系统计算机辅助设计MATLAB语言与应用东北大学信息学院7.2.1 小增益定理(a) 标准反馈控制结构(b) 小增益定理示意图2022-6-725控制

7、系统计算机辅助设计MATLAB语言与应用东北大学信息学院假设 为稳定的，则当且仅当小增益条件满足时图 (b) 中所示的系统对所有稳定的 都是良定的，且是内部稳定的。即如果系统的回路传递函数的范数小于 1，则闭环系统将总是稳定的。2022-6-726控制系统计算机辅助设计MATLAB语言与应用东北大学信息学院7.2.2 鲁棒控制器的结构闭环系统中引入的增广对象模型其对应的增广状态方程为2022-6-727控制系统计算机辅助设计MATLAB语言与应用东北大学信息学院闭环系统传递函数为2022-6-728控制系统计算机辅助设计MATLAB语言与应用东北大学信息学院 最优控制问题 其中需求解 ； 最优

8、控制问题 其中需求解 ； 控制问题 需要得出一个控制器满足鲁棒控制问题的三种形式： 鲁棒控制的目的是设计出一个镇定控制器 使得闭环系统 的范数取一个小于 1 的值，亦即2022-6-729控制系统计算机辅助设计MATLAB语言与应用东北大学信息学院加权灵敏度问题的控制结构框图 加权函数 ，使得 均正则。 即传递函数在 时均应该是有界的。2022-6-730控制系统计算机辅助设计MATLAB语言与应用东北大学信息学院式中假定系统对象模型的状态方程为 ，加权函数 的状态方程模型为 的状态方程模型为 ，而非正则的 的模型表示为 2022-6-731控制系统计算机辅助设计MATLAB语言与应用东北大学

9、信息学院这时鲁棒控制问题可以集中成下面三种形式： 灵敏度问题灵敏度问题 并不指定 稳定性与品质的混合鲁棒问题稳定性与品质的混合鲁棒问题 假定 为空 一般的混合灵敏度问题一般的混合灵敏度问题 要求三个加权函数都存在。2022-6-732控制系统计算机辅助设计MATLAB语言与应用东北大学信息学院7.2.3 鲁棒控制系统的 MATLAB 描述 鲁棒控制工具箱中的系统描述方法2022-6-733控制系统计算机辅助设计MATLAB语言与应用东北大学信息学院2022-6-734控制系统计算机辅助设计MATLAB语言与应用东北大学信息学院2022-6-735控制系统计算机辅助设计MATLAB语言与应用东北

10、大学信息学院【例7-6】2022-6-736控制系统计算机辅助设计MATLAB语言与应用东北大学信息学院 分析与综合工具箱和 LMI 工具箱的 模型描述2022-6-737控制系统计算机辅助设计MATLAB语言与应用东北大学信息学院变换出系统矩阵 P2022-6-738控制系统计算机辅助设计MATLAB语言与应用东北大学信息学院【例7-7】用【例7-6】中的对象模型和加权函数， 得出其系统矩阵模型 P 2022-6-739控制系统计算机辅助设计MATLAB语言与应用东北大学信息学院7.3 鲁棒控制器的 计算机辅助设计 鲁棒控制工具箱的设计方法鲁棒控制工具箱的设计方法 基于线性矩阵不等式工具箱的

11、设计方法基于线性矩阵不等式工具箱的设计方法 基于基于 分析与综合工具箱的分析与综合工具箱的 控制器设计控制器设计 基于回路成型技术的鲁棒控制器设计基于回路成型技术的鲁棒控制器设计2022-6-740控制系统计算机辅助设计MATLAB语言与应用东北大学信息学院7.3.1 鲁棒控制工具箱的 设计方法鲁棒控制器的状态方程表示其中X 与 Y 由下面的两个代数 Riccati 方程求解2022-6-741控制系统计算机辅助设计MATLAB语言与应用东北大学信息学院控制器存在的前提条件为 足够小, 且满足 ； 控制器 Riccati 方程的解为 正定矩阵； 观测器 Riccati 方程的解为 正定矩阵；

12、。该式说明两个 Riccati 方程的积矩阵的所有特征值均小于 。 2022-6-742控制系统计算机辅助设计MATLAB语言与应用东北大学信息学院【例7-8】对【例7-6】中的增广的系统模型，分别 设计2022-6-743控制系统计算机辅助设计MATLAB语言与应用东北大学信息学院绘制在控制器作用下系统的开环 Bode 图和闭环阶跃响应曲线2022-6-744控制系统计算机辅助设计MATLAB语言与应用东北大学信息学院【例7-9】设计最优 控制器，并绘制出该控制器作用下的阶跃响应曲线和开环系统的奇异值曲线。并设置加权矩阵2022-6-745控制系统计算机辅助设计MATLAB语言与应用东北大学

13、信息学院2022-6-746控制系统计算机辅助设计MATLAB语言与应用东北大学信息学院【例7-10】带有双积分器的非最小相位受控对象，选择加权函数并选择极点漂移为设计系统的最优 控制器。2022-6-747控制系统计算机辅助设计MATLAB语言与应用东北大学信息学院2022-6-748控制系统计算机辅助设计MATLAB语言与应用东北大学信息学院7.3.2 基于线性矩阵不等式 工具箱的设计方法 问题转换成线性矩阵不等式的最优化问题2022-6-749控制系统计算机辅助设计MATLAB语言与应用东北大学信息学院【例7-11】采用【例7-6】中增广的系统模型，用 LMI 工具箱的相关函数设计最优

14、控制器2022-6-750控制系统计算机辅助设计MATLAB语言与应用东北大学信息学院7.3.3 基于 分析与综合工具箱 的 控制器设计【例7-12】采用【例7-6】中增广的系统模型，用 分析与综合工具箱的相关函数设计最优 控制器2022-6-751控制系统计算机辅助设计MATLAB语言与应用东北大学信息学院2022-6-752控制系统计算机辅助设计MATLAB语言与应用东北大学信息学院7.3.4 基于回路成型技术的 鲁棒控制器设计假设前向回路的数学模型为 ，由典型反馈系统有 ，则系统的灵敏度控制传递函数 ，灵敏度函数2022-6-753控制系统计算机辅助设计MATLAB语言与应用东北大学信息

15、学院加权和数与回路成型示意图2022-6-754控制系统计算机辅助设计MATLAB语言与应用东北大学信息学院2022-6-755控制系统计算机辅助设计MATLAB语言与应用东北大学信息学院7.4 新鲁棒控制工具箱及应用7.4.1 不确定系统的描述2022-6-756控制系统计算机辅助设计MATLAB语言与应用东北大学信息学院【例7-13】典型二阶开环传函选定标称值为构造不确定系统模型。2022-6-757控制系统计算机辅助设计MATLAB语言与应用东北大学信息学院对叠加型不确定性对乘积型的不确定性2022-6-758控制系统计算机辅助设计MATLAB语言与应用东北大学信息学院7.4.2 灵敏度

16、问题的鲁棒控制器设计一般情况下，受控对象 G 的 D 矩阵为非满秩矩阵时，不能得出精确的成型控制器，这时回路奇异值的上下限满足式子当 时，控制器作用下实际回路奇异值介于 之间。2022-6-759控制系统计算机辅助设计MATLAB语言与应用东北大学信息学院【例7-14】2022-6-760控制系统计算机辅助设计MATLAB语言与应用东北大学信息学院绘制在此控制器下的回路奇异值及闭环系统的阶跃响应曲线2022-6-761控制系统计算机辅助设计MATLAB语言与应用东北大学信息学院7.4.3 混合灵敏度问题的鲁棒 控制器设计2022-6-762控制系统计算机辅助设计MATLAB语言与应用东北大学信

17、息学院【例7-15】2022-6-763控制系统计算机辅助设计MATLAB语言与应用东北大学信息学院假设系统的不确定部分为乘积型的，且已知 ，并已知不确定参数的变化范围为,设计固定的 控制器2022-6-764控制系统计算机辅助设计MATLAB语言与应用东北大学信息学院7.5 分数阶控制系统分析与设计7.5.1 分数阶微积分学与数值计算 分数阶微积分的定义2022-6-765控制系统计算机辅助设计MATLAB语言与应用东北大学信息学院当系数简单表示2022-6-766控制系统计算机辅助设计MATLAB语言与应用东北大学信息学院编写求取给定函数的分数阶微分函数2022-6-767控制系统计算机辅

18、助设计MATLAB语言与应用东北大学信息学院 Riemann-Liouville 定义为目前最常用的分数阶微积分定义2022-6-768控制系统计算机辅助设计MATLAB语言与应用东北大学信息学院Caputo 分数阶微分定义为Caputo 分数阶积分定义为2022-6-769控制系统计算机辅助设计MATLAB语言与应用东北大学信息学院 分数阶微积分的性质分数阶微积分的性质 解析函数 的分数阶导数 对 都是解析的。 为整数时，分数阶微分与整数阶微分的 值完全一致，且 。 分数阶微积分算子为线性的，即对任意常数 ， 有2022-6-770控制系统计算机辅助设计MATLAB语言与应用东北大学信息学院

19、 分数阶微积分算子满足交换律，并满足叠加关系 函数分数阶微分的 Laplace 变换为特别地，若函数 及其各阶导数的初值均为 0，则2022-6-771控制系统计算机辅助设计MATLAB语言与应用东北大学信息学院7.5.2 分数阶线性系统频域 与时域分析单变量线性系统的分数阶传递函数一般形式为2022-6-772控制系统计算机辅助设计MATLAB语言与应用东北大学信息学院 分数阶系统的频域分析2022-6-773控制系统计算机辅助设计MATLAB语言与应用东北大学信息学院 分数阶系统的时域分析2022-6-774控制系统计算机辅助设计MATLAB语言与应用东北大学信息学院7.5.3 分数阶微分

20、的滤波器近似及应用 分数阶微分的滤波器近似Oustaloup 算法滤波器零极点和增益为假设选定的拟合频率段为 ，则可以构造出连续滤波器的传递函数模型为2022-6-775控制系统计算机辅助设计MATLAB语言与应用东北大学信息学院编写设计连续滤波器的函数。2022-6-776控制系统计算机辅助设计MATLAB语言与应用东北大学信息学院【例7-16】2022-6-777控制系统计算机辅助设计MATLAB语言与应用东北大学信息学院2022-6-778控制系统计算机辅助设计MATLAB语言与应用东北大学信息学院 非线性分数阶系统的 Simulink 仿真2022-6-779控制系统计算机辅助设计MA

21、TLAB语言与应用东北大学信息学院【例7-17】用近似方法求解分数阶非线性微分方程2022-6-780控制系统计算机辅助设计MATLAB语言与应用东北大学信息学院7.5.4 分数阶系统的模型降阶技术2022-6-781控制系统计算机辅助设计MATLAB语言与应用东北大学信息学院【例7-18】利用最优降阶函数opt_app( )对其进行降阶处理，并绘制出高阶近似与最优降阶近似模型的阶跃响应曲线。2022-6-782控制系统计算机辅助设计MATLAB语言与应用东北大学信息学院7.5.5 分数阶系统的控制器设计分数阶 PID 控制器的数学模型为2022-6-783控制系统计算机辅助设计MATLAB语

22、言与应用东北大学信息学院【例7-19】根据 Wang-Juang-Chan 算法设计最优 ITAE 准则的PID 控制器2022-6-784控制系统计算机辅助设计MATLAB语言与应用东北大学信息学院2022-6-785控制系统计算机辅助设计MATLAB语言与应用东北大学信息学院【例7-20】 已知分数阶受控对象为 ， 其中分数阶次变化范围为 ，且标称 ，选择滤波器近似的值 ,选择加权函 数 ,设计最优 控制器2022-6-786控制系统计算机辅助设计MATLAB语言与应用东北大学信息学院2022-6-787控制系统计算机辅助设计MATLAB语言与应用东北大学信息学院【例7-21】 对【例7-

23、19】中的分数阶受控对象模型， 用 PID 控制的仿真、优化框图求解2022-6-788控制系统计算机辅助设计MATLAB语言与应用东北大学信息学院 分数阶系统的最优 PID 控制器设计框图2022-6-789控制系统计算机辅助设计MATLAB语言与应用东北大学信息学院 选择控制器参数为 若受控对象模型变化为2022-6-790控制系统计算机辅助设计MATLAB语言与应用东北大学信息学院7.6 本章要点简介 线性二次型线性二次型 Gauss 问题，即问题，即 LQG 问题，带有问题，带有 Kalman 滤波器滤波器 的的 LQG 控制器结构，并指出了该控制器结构存控制器结构，并指出了该控制器结

24、构存在的问题，引出了在的问题，引出了 LQG/LTR 控制器结构控制器结构 ，探讨了探讨了 q 值值选择对控制效果的影像。选择对控制效果的影像。 小增益定理以及基于范数的鲁棒控制三种形式：小增益定理以及基于范数的鲁棒控制三种形式： 控制控制、 控制及最优控制及最优 控制器，三种鲁棒控制问题，即灵控制器，三种鲁棒控制问题，即灵敏度问题、稳定性与品质的混合鲁棒问题及一般混合灵敏度问题、稳定性与品质的混合鲁棒问题及一般混合灵敏度问题，并介绍了在加权函数下受控对象模型的增广敏度问题，并介绍了在加权函数下受控对象模型的增广方法。方法。2022-6-791控制系统计算机辅助设计MATLAB语言与应用东北大

25、学信息学院 基于范数的鲁棒控制问题的基于范数的鲁棒控制问题的 MATLAB 描述方法和鲁描述方法和鲁棒控制器的计算机辅助设计的理论与求解方法，初步棒控制器的计算机辅助设计的理论与求解方法，初步探讨了加权函数的选取问题，并介绍了受控对象虚轴探讨了加权函数的选取问题，并介绍了受控对象虚轴上零点的双线性变换方法。上零点的双线性变换方法。 新版本的鲁棒控制工具箱将三种著名的方法，统一到新版本的鲁棒控制工具箱将三种著名的方法，统一到一个框架下，给出了统一的模型描述与设计函数。一个框架下，给出了统一的模型描述与设计函数。 分数阶微积分的定义与数值求解方法，分数阶系统的分数阶微积分的定义与数值求解方法，分数阶系统的时域与频域分析方法、分数阶系统的整数阶近似算法时域与频域分析方法、分数阶系统的整数阶近似算法、分数阶受控对象模型的整数阶、分数阶受控对象模型的整数阶 PID 控制器以及控制器以及 控制器的设计方法。控制器的设计方法。