鲁棒控制与鲁棒控制器设计-课件

第7章

鲁棒控制与鲁棒控制器设计2023/7/202控制系统计算机辅助设计——MATLAB语言与应用主要内容线性二次型Gauss控制鲁棒控制问题的一般描述

鲁棒控制器的计算机辅助设计新鲁棒控制工具箱及应用分数阶控制系统分析与设计2023/7/203控制系统计算机辅助设计——MATLAB语言与应用7.1线性二次型Gauss控制7.1.1线性二次型Gauss问题假设对象模型的状态方程表示为为白噪声信号，分别表示模型的不确定性与输出信号的量测噪声。2023/7/204控制系统计算机辅助设计——MATLAB语言与应用定义最优控制的指标函数为2023/7/205控制系统计算机辅助设计——MATLAB语言与应用7.1.2使用MATLAB求解LQG问题带有Kalman滤波器的LQG结构2023/7/206控制系统计算机辅助设计——MATLAB语言与应用Kalman滤波器的增益矩阵式中,满足下面的Riccati代数方程2023/7/207控制系统计算机辅助设计——MATLAB语言与应用【例7-1】2023/7/208控制系统计算机辅助设计——MATLAB语言与应用LQG控制器设计的分离原理2023/7/209控制系统计算机辅助设计——MATLAB语言与应用基于观测器的LQG调节器设计2023/7/2010控制系统计算机辅助设计——MATLAB语言与应用由Kalman滤波器方程可以写出基于观测器的LQG调节器为2023/7/2011控制系统计算机辅助设计——MATLAB语言与应用2023/7/2012控制系统计算机辅助设计——MATLAB语言与应用【例7-2】2023/7/2013控制系统计算机辅助设计——MATLAB语言与应用7.1.3带有回路传输恢复的LQG控制LQG/LTR控制器设计算法使用LQG控制器，系统的开环传递函数表示为直接状态反馈系统的开环传递函数为2023/7/2014控制系统计算机辅助设计——MATLAB语言与应用【例7-3】2023/7/2015控制系统计算机辅助设计——MATLAB语言与应用2023/7/2016控制系统计算机辅助设计——MATLAB语言与应用回路传输恢复技术(looptransferrecovery，LTR)加权函数的选择2023/7/2017控制系统计算机辅助设计——MATLAB语言与应用

先求解标准的LQ问题，然后应用LTR技术2023/7/2018控制系统计算机辅助设计——MATLAB语言与应用2023/7/2019控制系统计算机辅助设计——MATLAB语言与应用【例7-4】对【例7-3】不同的q值应用LTR技术2023/7/2020控制系统计算机辅助设计——MATLAB语言与应用2023/7/2021控制系统计算机辅助设计——MATLAB语言与应用应用MATLAB求解LQG/LTR问题若想使得系统在输入端恢复环路传递函数，则若想在对象模型的输出端恢复环路传递函数，则2023/7/2022控制系统计算机辅助设计——MATLAB语言与应用【例7-5】对【例7-3】选定一个q向量，设计

LTR控制器，并绘制出不同q值下环路传递函数的Nyquist图。2023/7/2023控制系统计算机辅助设计——MATLAB语言与应用7.2

鲁棒控制问题的一般描述小增益定理鲁棒控制器的结构鲁棒控制系统的MATLAB描述2023/7/2024控制系统计算机辅助设计——MATLAB语言与应用7.2.1小增益定理(a)标准反馈控制结构(b)小增益定理示意图2023/7/2025控制系统计算机辅助设计——MATLAB语言与应用假设

为稳定的，则当且仅当小增益条件满足时图(b)中所示的系统对所有稳定的都是良定的，且是内部稳定的。小增益定理即如果系统的回路传递函数的范数小于1，则闭环系统将总是稳定的。2023/7/2026控制系统计算机辅助设计——MATLAB语言与应用7.2.2鲁棒控制器的结构闭环系统中引入的增广对象模型其对应的增广状态方程为2023/7/2027控制系统计算机辅助设计——MATLAB语言与应用闭环系统传递函数为2023/7/2028控制系统计算机辅助设计——MATLAB语言与应用最优控制问题

其中需求解；

最优控制问题

其中需求解；

控制问题

需要得出一个控制器满足鲁棒控制问题的三种形式：鲁棒控制的目的是设计出一个镇定控制器使得闭环系统的范数取一个小于1的值，亦即2023/7/2029控制系统计算机辅助设计——MATLAB语言与应用加权灵敏度问题的控制结构框图

加权函数，使得均正则。即传递函数在时均应该是有界的。2023/7/2030控制系统计算机辅助设计——MATLAB语言与应用式中假定系统对象模型的状态方程为，加权函数的状态方程模型为的状态方程模型为，而非正则的的模型表示为2023/7/2031控制系统计算机辅助设计——MATLAB语言与应用这时鲁棒控制问题可以集中成下面三种形式：灵敏度问题

并不指定稳定性与品质的混合鲁棒问题

假定为空一般的混合灵敏度问题

要求三个加权函数都存在。2023/7/2032控制系统计算机辅助设计——MATLAB语言与应用7.2.3鲁棒控制系统的

MATLAB描述

鲁棒控制工具箱中的系统描述方法建立鲁棒控制工具箱可以使用的系统模型2023/7/2033控制系统计算机辅助设计——MATLAB语言与应用2023/7/2034控制系统计算机辅助设计——MATLAB语言与应用2023/7/2035控制系统计算机辅助设计——MATLAB语言与应用【例7-6】2023/7/2036控制系统计算机辅助设计——MATLAB语言与应用分析与综合工具箱和LMI工具箱的模型描述2023/7/2037控制系统计算机辅助设计——MATLAB语言与应用变换出系统矩阵P2023/7/2038控制系统计算机辅助设计——MATLAB语言与应用【例7-7】用【例7-6】中的对象模型和加权函数，

得出其系统矩阵模型P2023/7/2039控制系统计算机辅助设计——MATLAB语言与应用7.3鲁棒控制器的

计算机辅助设计鲁棒控制工具箱的设计方法基于线性矩阵不等式工具箱的设计方法基于分析与综合工具箱的控制器设计基于回路成型技术的鲁棒控制器设计2023/7/2040控制系统计算机辅助设计——MATLAB语言与应用7.3.1鲁棒控制工具箱的

设计方法鲁棒控制器的状态方程表示其中X与Y由下面的两个代数Riccati方程求解2023/7/2041控制系统计算机辅助设计——MATLAB语言与应用控制器存在的前提条件为足够小,且满足；控制器Riccati方程的解为正定矩阵；观测器Riccati方程的解为正定矩阵；。该式说明两个Riccati方程的积矩阵的所有特征值均小于。

2023/7/2042控制系统计算机辅助设计——MATLAB语言与应用【例7-8】对【例7-6】中的增广的系统模型，分别

设计2023/7/2043控制系统计算机辅助设计——MATLAB语言与应用绘制在控制器作用下系统的开环Bode图和闭环阶跃响应曲线2023/7/2044控制系统计算机辅助设计——MATLAB语言与应用【例7-9】设计最优控制器，并绘制出该控制器作用下的阶跃响应曲线和开环系统的奇异值曲线。并设置加权矩阵2023/7/2045控制系统计算机辅助设计——MATLAB语言与应用2023/7/2046控制系统计算机辅助设计——MATLAB语言与应用【例7-10】带有双积分器的非最小相位受控对象，选择加权函数并选择极点漂移为设计系统的最优控制器。2023/7/2047控制系统计算机辅助设计——MATLAB语言与应用2023/7/2048控制系统计算机辅助设计——MATLAB语言与应用7.3.2基于线性矩阵不等式

工具箱的设计方法问题转换成线性矩阵不等式的最优化问题2023/7/2049控制系统计算机辅助设计——MATLAB语言与应用【例7-11】采用【例7-6】中增广的系统模型，用LMI

工具箱的相关函数设计最优控制器2023/7/2050控制系统计算机辅助设计——MATLAB语言与应用7.3.3基于分析与综合工具箱

的控制器设计【例7-12】采用【例7-6】中增广的系统模型，用

分析与综合工具箱的相关函数设计最优控制器2023/7/2051控制系统计算机辅助设计——MATLAB语言与应用2023/7/2052控制系统计算机辅助设计——MATLAB语言与应用7.3.4基于回路成型技术的

鲁棒控制器设计假设前向回路的数学模型为，由典型反馈系统有，则系统的灵敏度控制传递函数，灵敏度函数2023/7/2053控制系统计算机辅助设计——MATLAB语言与应用加权和数与回路成型示意图2023/7/2054控制系统计算机辅助设计——MATLAB语言与应用2023/7/2055控制系统计算机辅助设计——MATLAB语言与应用7.4新鲁棒控制工具箱

及应用7.4.1不确定系统的描述2023/7/2056控制系统计算机辅助设计——MATLAB语言与应用【例7-13】典型二阶开环传函选定标称值为构造不确定系统模型。2023/7/2057控制系统计算机辅助设计——MATLAB语言与应用对叠加型不确定性对乘积型的不确定性2023/7/2058控制系统计算机辅助设计——MATLAB语言与应用7.4.2灵敏度问题的鲁棒控制器设计一般情况下，受控对象G的D矩阵为非满秩矩阵时，不能得出精确的成型控制器，这时回路奇异值的上下限满足式子当时，控制器作用下实际回路奇异值介于之间。2023/7/2059控制系统计算机辅助设计——MATLAB语言与应用【例7-14】2023/7/2060控制系统计算机辅助设计——MATLAB语言与应用绘制在此控制器下的回路奇异值及闭环系统的阶跃响应曲线2023/7/2061控制系统计算机辅助设计——MATLAB语言与应用7.4.3混合灵敏度问题的鲁棒

控制器设计2023/7/2062控制系统计算机辅助设计——MATLAB语言与应用【例7-15】2023/7/2063控制系统计算机辅助设计——MATLAB语言与应用假设系统的不确定部分为乘积型的，且已知

，并已知不确定参数的变化范围为,设计固定的控制器2023/7/2064控制系统计算机辅助设计——MATLAB语言与应用7.5分数阶控制系统分析与设计7.5.1分数阶微积分学与数值计算分数阶微积分的定义2023/7/2065控制系统计算机辅助设计——MATLAB语言与应用当系数简单表示2023/7/2066控制系统计算机辅助设计——MATLAB语言与应用编写求取给定函数的分数阶微分函数2023/7/2067控制系统计算机辅助设计——MATLAB语言与应用Riemann-Liouville定义为目前最常用的分数阶微积分定义2023/7/2068控制系统计算机辅助设计——MATLAB语言与应用Caputo分数阶微分定义为Caputo分数阶积分定义为2023/7/2069控制系统计算机辅助设计——MATLAB语言与应用

分数阶微积分的性质①解析函数的分数阶导数对都是解析的。②为整数时，分数阶微分与整数阶微分的值完全一致，且。③分数阶微积分算子为线性的，即对任意常数，

有2023/7/2070控制系统计算机辅助设计——MATLAB语言与应用④分数阶微积分算子满足交换律，并满足叠加关系⑤函数分数阶微分的Laplace变换为特别地，若函数及其各阶导数的初值均为0，则2023/7/2071控制系统计算机辅助设计——MATLAB语言与应用7.5.2分数阶线性系统频域

与时域分析单变量线性系统的分数阶传递函数一般形式为2023/7/2072控制系统计算机辅助设计——MATLAB语言与应用分数阶系统的频域分析2023/7/2073控制系统计算机辅助设计——MATLAB语言与应用分数阶系统的时域分析2023/7/2074控制系统计算机辅助设计——MATLAB语言与应用7.5.3分数阶微分的滤波器近似及应用分数阶微分的滤波器近似Oustaloup算法滤波器零极点和增益为假设选定的拟合频率段为，则可以构造出连续滤波器的传递函数模型为2023/7/2075控制系统计算机辅助设计——MATLAB语言与应用编写设计连续滤波器的函数。2023/7/2076控制系统计算机辅助设计——MATLAB语言与应用【例7-16】2023/7/2077控制系统计算机辅助设计——MATLAB语言与应用2023/7/2078控制系统计算机辅助设计——MATLAB语言与应用

非线性分数阶系统的Simulink仿真2023/7/2079控制系统计算机辅助设计——MATLAB语言与应用【例7-17】用近似方法求解分数阶非线性微分方程2023/7/2080控制系统计算机辅助设计——MATLAB语言与应用7.5.4分数阶系统的模型降阶技术2023/7/2081控制系统计算机辅助设计——MATLAB语言与应用【例7-18】利用最优降阶函数opt_app()对其进行降阶处理，并绘制出高阶近似与最优降阶近似模型的阶跃响应曲线。2023/7/2082控制系统计算机辅助设计——MATLAB语言与应用7.5.5分数阶系统的控制器设计分数阶PID控制器的数学模型为2023/7/2083控制系统计算机辅助设计——MATLAB语言与应用【例7-19】根据Wang-Juang-Chan算法设计最优ITAE准则的PID控制器2023/7/2084控制系统计算机辅助设计——MATLAB语言与应用2023/7/2085控制系统计算机辅助设计——MATLAB语言与应用【例7-20】已知分数阶受控对象为，

其中分数阶次变化范围为，且标称

，选择滤波器近似的值