控制系统计算机辅助设计.ppt

《控制系统计算机辅助设计.ppt》由会员分享，可在线阅读，更多相关《控制系统计算机辅助设计.ppt（124页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、第第 6 6 章章控制系统计算机辅助设计控制系统计算机辅助设计2022/11/281主要内容主要内容n n基于传递函数的控制器基于传递函数的控制器设计设计方法方法n n状态反馈控制状态反馈控制n n基于状态空间模型的控制器设计方法基于状态空间模型的控制器设计方法2022/11/2826.16.1基于传递函数的控制器基于传递函数的控制器设计方法设计方法6.1.1 串联超前滞后校正器串联超前滞后校正器2022/11/283 超前校正器超前校正器2022/11/284 滞后校正器滞后校正器2022/11/285 超前滞后校正器超前滞后校正器2022/11/2866.1.2 超前滞后校正器的设计方法超

2、前滞后校正器的设计方法 基于剪切频率和相位裕度的设计方法基于剪切频率和相位裕度的设计方法基于剪切频率和相位裕度的设计方法基于剪切频率和相位裕度的设计方法2022/11/287超前滞后校正器的设计规则：超前滞后校正器的设计规则：超前滞后校正器的设计规则：超前滞后校正器的设计规则：且且系统静态误差系数为系统静态误差系数为2022/11/2882022/11/289【例例6-16-1】2022/11/2810超前滞后校正器超前滞后校正器超前滞后校正器超前滞后校正器超前校正器超前校正器超前校正器超前校正器2022/11/28112022/11/28122022/11/2813 基于模型匹配算法的设计方

3、法基于模型匹配算法的设计方法基于模型匹配算法的设计方法基于模型匹配算法的设计方法假设受控对象的传递函数为假设受控对象的传递函数为 ，期望闭环系统的频域响应为期望闭环系统的频域响应为 ，超前滞后校正器的一般形式为超前滞后校正器的一般形式为使得在频率段使得在频率段 内闭环模型对期望闭环模型内闭环模型对期望闭环模型 匹配指标匹配指标为最小为最小2022/11/2814提出了下面的设计算法提出了下面的设计算法其中其中2022/11/2815其中，其中，gp 和和 f 分别为受控对象和期望闭环分别为受控对象和期望闭环系统的传递函数模型，系统的传递函数模型，w1 和和 w2 为需要拟为需要拟合的频率段上下

4、限。合的频率段上下限。2022/11/2816【例例例例6-26-26-26-2】受控对象模型为受控对象模型为受控对象模型为受控对象模型为2022/11/28176.1.3 控制系统工具箱中的设计界面控制系统工具箱中的设计界面控制器设计界面控制器设计界面 界面允许选择和修改控制器的结构，允许添界面允许选择和修改控制器的结构，允许添界面允许选择和修改控制器的结构，允许添界面允许选择和修改控制器的结构，允许添加零极点，调整增益，从而设计出控制器模型。加零极点，调整增益，从而设计出控制器模型。加零极点，调整增益，从而设计出控制器模型。加零极点，调整增益，从而设计出控制器模型。2022/11/2818

5、【例例例例6-36-36-36-3】受控对象和控制器的传递函数模型分别为受控对象和控制器的传递函数模型分别为受控对象和控制器的传递函数模型分别为受控对象和控制器的传递函数模型分别为2022/11/28196.2 6.2 基于状态空间模型的基于状态空间模型的 控制器设计方法控制器设计方法6.2.1 状态反馈控制状态反馈控制2022/11/2820 将将 代入开环系统的状态方代入开环系统的状态方程模型，则在状态反馈矩阵程模型，则在状态反馈矩阵 下，系统的闭环状下，系统的闭环状态方程模型可以写成态方程模型可以写成 如果系统如果系统 完全可控，则选择合适的完全可控，则选择合适的 矩矩阵，可以将闭环系统

6、矩阵阵，可以将闭环系统矩阵 的特征值配置的特征值配置到任意地方。到任意地方。2022/11/28216.2.2 线性二次型指标最优调节器线性二次型指标最优调节器假设线性时不变系统的状态方程模型为假设线性时不变系统的状态方程模型为设计一个输入量设计一个输入量 ,使得最优控制性能指标使得最优控制性能指标最小最小2022/11/2822则控制信号应该为则控制信号应该为由简化的由简化的 Riccati 微分方程微分方程 求出求出假设假设 ，其中其中 ，则则可以得出在状态反馈下的闭环系统的状态方程为可以得出在状态反馈下的闭环系统的状态方程为依照给定加权矩阵设计的依照给定加权矩阵设计的依照给定加权矩阵设计

7、的依照给定加权矩阵设计的 LQ LQ 最优控制器最优控制器最优控制器最优控制器2022/11/2823离散系统二次型性能指标离散系统二次型性能指标离散离散 Riccati 代数方程代数方程这时控制律为这时控制律为2022/11/2824【例例例例6-46-46-46-4】2022/11/28256.2.3 极点配置控制器设计极点配置控制器设计系统的状态方程为系统的状态方程为则系统的闭环状态方程为则系统的闭环状态方程为2022/11/28262022/11/2827n n Bass-Gura Bass-Gura 算法算法算法算法2022/11/2828基于此算法编写的基于此算法编写的 MATLA

8、B 函数函数2022/11/2829n nAckermann Ackermann 算法算法算法算法其中其中 为将为将 代入代入 得出的矩阵多项式的值得出的矩阵多项式的值n鲁棒极点配置算法鲁棒极点配置算法place()函数不适用于含有多重期望极点的问题函数不适用于含有多重期望极点的问题acker()函数可以求解配置多重极点的问题函数可以求解配置多重极点的问题2022/11/2830【例例例例6-56-56-56-5】2022/11/2831【例例6-66-6】2022/11/28326.2.4 观测器设计及基于观测器的观测器设计及基于观测器的 调节器设计调节器设计2022/11/28332022

9、/11/28342022/11/2835【例例例例6-76-76-76-7】2022/11/28362022/11/2837带有观测器的状态反馈控制结构图带有观测器的状态反馈控制结构图2022/11/28382022/11/28392022/11/2840如果参考输入信号如果参考输入信号 ，则控制结构，则控制结构 化简为化简为2022/11/2841【例例例例6-86-86-86-8】2022/11/28422022/11/28436.3 6.3 过程控制系统的过程控制系统的 PID 控制器设计控制器设计6.3.1 PID 控制器概述控制器概述n连续连续 PID 控制器控制器2022/11/2

10、844连续连续 PID 控制器控制器Laplace 变换形式变换形式2022/11/2845n离散离散 PID 控制器控制器2022/11/2846离散形式的离散形式的 PID 控制器控制器Z 变换得到的离散变换得到的离散 PID 控制器的传递函数控制器的传递函数2022/11/2847nPID 控制器的变形控制器的变形 积分分离式积分分离式 PID 控制器控制器在启动过程中，如果静态误差很大时，可以关闭积在启动过程中，如果静态误差很大时，可以关闭积分部分的作用，稳态误差很小时再开启积分作用，分部分的作用，稳态误差很小时再开启积分作用，消除静态误差消除静态误差2022/11/2848 离散增量

11、式离散增量式 PID 控制器控制器2022/11/2849 抗积分饱和抗积分饱和(anti-windup)PID 控制器控制器 2022/11/28506.3.2 过程系统的一阶延迟模型近似过程系统的一阶延迟模型近似带有时间延迟一阶模型带有时间延迟一阶模型(first-order lag plus delay，FOLPD)一阶延迟模型一阶延迟模型(FOLPD)的数学表示为的数学表示为2022/11/2851n n 由响应曲线识别一阶模型由响应曲线识别一阶模型 阶跃响应近似阶跃响应近似 Nyquist 图近似图近似 编写编写 MATLAB 函数函数 getfolpd(),key=12022/11

12、/2852n 基于频域响应的近似方法基于频域响应的近似方法调用编写的调用编写的 MATLAB 函数函数 getfolpd(),key=22022/11/2853n 基于传递函数的辨识方法基于传递函数的辨识方法调用编写的调用编写的 MATLAB 函数函数 getfolpd(),key=32022/11/2854n 最优降阶方法最优降阶方法调用编写的调用编写的 MATLAB 函数函数 getfolpd(),key=4【例例6-96-9】2022/11/28556.3.3 Ziegler-Nichols 参数整定方法参数整定方法n Ziegler-Nichols 经验公式经验公式编写编写 MATLA

13、B 函数函数 ziegler()2022/11/2856【例例6-106-10】2022/11/28572022/11/2858n改进的改进的 Ziegler-Nichols 算法算法2022/11/2859初始点初始点 A 增益增益期望点期望点 A1 增益增益PID 控制器控制器2022/11/2860n n PI 控制器2022/11/2861n n PID 控制器2022/11/2862【例例6-116-11】2022/11/28632022/11/2864n 改进改进 PID 控制结构与算法控制结构与算法n n微分动作在反馈回路的微分动作在反馈回路的 PID PID 控制器控制器202

14、2/11/2865n n精调的精调的 Ziegler-Nichols Ziegler-Nichols 控制器及算法控制器及算法2022/11/28662022/11/2867 若若 则保留则保留 Ziegler-Nichols 参数参数,同时为使超调量分别小于同时为使超调量分别小于 10%或或 20%，则则 若若 ，Ziegler-Nichols 控制器的控制器的 参数精调为参数精调为 若若 ，为使系统的超调量小于为使系统的超调量小于 10%，则，则 PID参数调为参数调为：2022/11/2868【例例6-126-12】用自编的用自编的 MATLAB 函数设计精调的函数设计精调的 Ziegl

15、er-Nichols PID 控制器控制器2022/11/2869n n改进的改进的 PID PID 结构结构一种一种 PID 控制器结构及整定算法的控制器模型为控制器结构及整定算法的控制器模型为:2022/11/28706.3.4 最优最优 PID 整定算法整定算法 最优化指标最优化指标 时间加权的指标时间加权的指标 IAE 和和 ITAE 指标指标2022/11/2871 庄敏霞与庄敏霞与 Atherton 教授提出了基于时间加权指标教授提出了基于时间加权指标的最优控制的最优控制 PID 控制器参数整定经验公式控制器参数整定经验公式适用范围适用范围 ，不适合于大时间延迟系统不适合于大时间延

16、迟系统2022/11/2872 Murrill 提出了使得提出了使得 IAE 准则最小的准则最小的 PID 控制器算法控制器算法2022/11/2873 对对 ITAE 指标进行最优化，得出的指标进行最优化，得出的 PID 控制器控制器设计经验公式设计经验公式 在在 范围内设计的范围内设计的 ITAE 最优最优 PID 控制器的经验公式控制器的经验公式 2022/11/2874【例例6-136-13】2022/11/28752022/11/28766.3.5 其他模型的其他模型的 PID 控制器参数控制器参数 整定算法整定算法 IPD 模型的模型的 PD 和和 PID 参数整定参数整定(int

17、egrator plus delay)2022/11/2877各种指标下的各种指标下的 PD 和和 PID 参数整定公式参数整定公式若选择若选择 ISE 指标指标,则则若选择若选择 ITSE 指标指标,则则若选择若选择 ISTSE 指标指标,则则2022/11/2878编写设计控制器的编写设计控制器的 MATLAB 函数函数2022/11/2879 FOLIPD 模型的模型的 PD 和和 PID 参数整定参数整定(first order lag and integrator plus delay)PID 控制器的整定算法控制器的整定算法 PD 控制器的设计算法控制器的设计算法2022/11/2

18、880编写设计控制器的编写设计控制器的 MATLAB 函数函数2022/11/2881【例例6-146-14】2022/11/2882 不稳定不稳定 FOLPD 模型的模型的 PID 参数整定参数整定 设计的设计的PID 控制器控制器若使若使 ISE 指标最小指标最小,则则若使若使 ITSE 指标最小指标最小,则则若使若使 ISTSE 指标最小，则指标最小，则2022/11/2883不稳定不稳定 FOLPD 模型的模型的 PID 控制器参数整定函数控制器参数整定函数2022/11/28846.3.6 基于基于 FOLPD 的的 PID 控制器控制器 设计程序设计程序 在在 MATLAB 提示符

19、下输入提示符下输入 pid_tuner。单击单击 Plant model 按钮，打开一个允许用户输入受按钮，打开一个允许用户输入受控对象模型参数的对话框。控对象模型参数的对话框。输入了受控对象模型后，单击输入了受控对象模型后，单击 Get FOLPD parameters 按钮获得按钮获得 FOLPD 模型，模型，亦即获得并显亦即获得并显示示 K,L,T 参数。参数。2022/11/2885 通过得出的通过得出的K,L,T 参数参数,设计所需的控制器。设计所需的控制器。单击单击 Design controller 按钮按钮,将自动设计出所需将自动设计出所需的的 PID 控制器模型，并将其显示出

20、来。控制器模型，并将其显示出来。单击单击 Closed-loop Simulation 按钮，则可以构造出按钮，则可以构造出 PID 控制器控制下的系统仿真模型，并在图形界面控制器控制下的系统仿真模型，并在图形界面上显示系统的阶跃响应曲线。上显示系统的阶跃响应曲线。2022/11/28866.4 6.4 最优控制器设计最优控制器设计6.4.1 最优控制的概念最优控制的概念在一定的具体条件下，要完成某个控制任务，使在一定的具体条件下，要完成某个控制任务，使得选定指标最小或增大的控制得选定指标最小或增大的控制.积分型误差指标、时间最短指标、能量最省指标等积分型误差指标、时间最短指标、能量最省指标等

21、2022/11/2887【例例6-166-16】设计最优控制器设计最优控制器2022/11/2888 为使得为使得 ITAE 准则最小化，可以编写如下的准则最小化，可以编写如下的 MATLAB 函数函数2022/11/28892022/11/2890为了降低超调量，改进的仿真框图为了降低超调量，改进的仿真框图2022/11/28912022/11/2892【例例6-176-17】考虑前面的例子，假设可以接受的控制信号考虑前面的例子，假设可以接受的控制信号 限幅值为限幅值为 20 202022/11/28932022/11/28946.4.2 基于基于 MATLAB/Simulink 的最优的最

22、优 控制程序及其应用控制程序及其应用最优控制器设计程序最优控制器设计程序(Optimal Controller Designer，OCD)的调用过程为：的调用过程为：在在 MATLAB 提示符下输入提示符下输入 ocd。建立一个建立一个 Simulink 仿真模型，该模型至少包含待仿真模型，该模型至少包含待优化的参数变量和误差信号的准则。优化的参数变量和误差信号的准则。将对应的将对应的 Simulink 模型名填写到界面的模型名填写到界面的Select a Simulink model 编辑框中。编辑框中。2022/11/2895 将待优化变量名填写到将待优化变量名填写到 Select var

23、iables to be optimized 编辑框中，且各个变量名之间用逗号分隔。编辑框中，且各个变量名之间用逗号分隔。估计指标收敛的时间段作为终止仿真时间估计指标收敛的时间段作为终止仿真时间,填写到填写到Simulation terminate time 栏目中去栏目中去。单击单击 Create File 按钮自动生成描述目标函数的按钮自动生成描述目标函数的 MATLAB 文件文件 opt_*.m。单击单击 Optimize 按钮将启动优化过程。按钮将启动优化过程。本程序允许用户指定优化变量的上下界，选择优化参本程序允许用户指定优化变量的上下界，选择优化参数的初值，选择不同的寻优算法，选择

24、离散仿真算法等。数的初值，选择不同的寻优算法，选择离散仿真算法等。2022/11/2896【例例例例6-186-186-186-18】受控对象的模型为受控对象的模型为受控对象的模型为受控对象的模型为 用最优控制器设计程序选择用最优控制器设计程序选择用最优控制器设计程序选择用最优控制器设计程序选择 PID PID PID PID 控制器参数。控制器参数。控制器参数。控制器参数。2022/11/2897自动生成目标函数的自动生成目标函数的 MATLAB：2022/11/2898【例例例例6-196-196-196-19】用用用用 OCD OCD OCD OCD 同时设计串级控制器同时设计串级控制器

25、同时设计串级控制器同时设计串级控制器2022/11/2899Simulink 仿真模型仿真模型2022/11/28100【例例例例6-206-206-206-20】对模型对模型对模型对模型 采用采用采用采用 ISE ISE ISE ISE 准准准准 则设计最优控制器。则设计最优控制器。则设计最优控制器。则设计最优控制器。2022/11/281016.4.3 最优控制程序的其他应用最优控制程序的其他应用【例例6-216-21】对模型对模型 采用采用 ITAE ITAE 准准 则，用则，用 OCD OCD 来进行最优降阶研究。来进行最优降阶研究。2022/11/28102 6.5 6.5 多变量系

26、统的频域多变量系统的频域 设计方法设计方法n n逆逆 Nyquist Nyquist 阵列方法阵列方法n n特征轨迹法特征轨迹法(characteristic locus method)characteristic locus method)n n反标架坐标法反标架坐标法(reversed-frame normalisationreversed-frame normalisation，RFN)RFN)n n序贯回路闭合方法序贯回路闭合方法(sequential loop closing)sequential loop closing)n n参数最优化方法参数最优化方法(parameters o

27、ptimisation method)parameters optimisation method)2022/11/281036.5.1 对角占优系统与伪对角化对角占优系统与伪对角化 为预补偿矩阵，它使得为预补偿矩阵，它使得 为对角占为对角占优矩阵。优矩阵。对所得对角占优矩阵作动态的补偿对所得对角占优矩阵作动态的补偿。2022/11/28104由以下步骤求取最优的补偿矩阵由以下步骤求取最优的补偿矩阵 ：选择一个函数的频率点选择一个函数的频率点 ，求出系统的逆求出系统的逆 Nyquist 阵阵 列列 。对各个对各个 值值 ，构成一个矩阵，构成一个矩阵 ，其中其中假设在假设在 频率处的系统传递函数

28、矩阵的逆频率处的系统传递函数矩阵的逆 Nyquist 阵列表示为阵列表示为2022/11/28105 求取求取 矩阵的特征值与特征向量，并将最小特征值的特矩阵的特征值与特征向量，并将最小特征值的特 征向量记作征向量记作 。由上面的各个由上面的各个 值得出的最小特征向量可以构成补偿矩阵值得出的最小特征向量可以构成补偿矩阵 选择选择 个频率点个频率点 ，并假设对第，并假设对第 个频率个频率点引入加权系数点引入加权系数 ，按照如下的方法构造矩阵按照如下的方法构造矩阵2022/11/28106由由 MATLAB 编写出为对角化函数编写出为对角化函数 pseudiag()2022/11/28107【例例

29、例例6-226-226-226-22】2022/11/281082022/11/28109【例例例例6-236-236-236-23】2022/11/28110引入动态补偿矩阵引入动态补偿矩阵2022/11/28111利用利用 Simulink 模型，绘制系统的阶跃响应曲线模型，绘制系统的阶跃响应曲线2022/11/281122022/11/281136.5.2 多变量系统的参数最优化设计多变量系统的参数最优化设计系统的闭环传递函数矩阵系统的闭环传递函数矩阵2022/11/28114控制器参数的最小二乘解控制器参数的最小二乘解2022/11/28115【例例例例6-246-246-246-24

30、】2022/11/28116系统选择闭环目标传递函数为系统选择闭环目标传递函数为2022/11/28117求目标控制器求目标控制器 ，并绘制并绘制 Bode 图图2022/11/28118按下面方式设置控制器的结构按下面方式设置控制器的结构,并建立分母矩阵并建立分母矩阵2022/11/28119绘制在控制器作用下，系统的阶跃响应输出曲线绘制在控制器作用下，系统的阶跃响应输出曲线2022/11/281206.5.3 基于基于 OCD 的多变量系统最优设计的多变量系统最优设计【例例6-256-25】采用加权采用加权 ITAE 准则下的最优准则下的最优 PI 控制器控制器 设计设计2022/11/2

31、81212022/11/281226.6 6.6 本章要点小结本章要点小结 超前、滞后与超前滞后串联校正器及其在系统控制超前、滞后与超前滞后串联校正器及其在系统控制 中的原理与意义，基于剪切频率与相位裕度配置的中的原理与意义，基于剪切频率与相位裕度配置的 校正器设计算法及其校正器设计算法及其 MATLAB MATLAB 实现，实现，MATLAB MATLAB 提提 供的基于根轨迹和供的基于根轨迹和 Bode Bode 图的控制器设计界面及其图的控制器设计界面及其 应用。应用。状态反馈的基本概念，及两种有影响的状态反馈控状态反馈的基本概念，及两种有影响的状态反馈控 制结构：基于二次型指标的最优控

32、制器设计及极点制结构：基于二次型指标的最优控制器设计及极点 配置控制器设计方法。观测器的概念，观测器的设配置控制器设计方法。观测器的概念，观测器的设 计方法，及基于观测器的控制结构及其应用。计方法，及基于观测器的控制结构及其应用。2022/11/28123各种常用的各种常用的 PID PID 控制器结构，侧重于受控对象模型控制器结构，侧重于受控对象模型 带有时间延迟的一阶模型带有时间延迟的一阶模型 FOLPD FOLPD 的近似和基于这类的近似和基于这类 受控对象模型的受控对象模型的 PID PID 控制器设计算法。控制器设计算法。最优控制的基本概念，基于数值最优化方法的最优最优控制的基本概念，基于数值最优化方法的最优 控制器设计，及作者编写的最优控制器设计程序控制器设计，及作者编写的最优控制器设计程序 OCD OCD，演示了它在最优控制器设计及模型最优拟合，演示了它在最优控制器设计及模型最优拟合 中的应用。中的应用。基于逆基于逆 NyquistNyquist 阵列的对角占优及伪对角化设计算阵列的对角占优及伪对角化设计算 法，参数最优化设计算法。法，参数最优化设计算法。2022/11/28124