MALAB语言及其应用课件

2024/1/261可编辑第6章

控制系统计算机辅助设计薛定宇著《控制系统计算机辅助设计—MATLAB语言与应用》第二版，清华大学出版社2006CAI课件开发：鄂大志、薛定宇2024/1/262可编辑主要内容超前滞后校正器设计方法基于状态空间模型的控制器设计方法过程控制系统的PID控制器设计最优控制器设计多变量系统的频域设计方法2024/1/263可编辑6.1超前滞后校正器

设计方法6.1.1串联超前滞后校正器2024/1/264可编辑超前校正器2024/1/265可编辑滞后校正器2024/1/266可编辑超前滞后校正器2024/1/267可编辑6.1.2超前滞后校正器的设计方法基于剪切频率和相位裕度的设计方法2024/1/268可编辑超前滞后校正器的设计规则：且系统静态误差系数为2024/1/269可编辑2024/1/2610可编辑【例6-1】2024/1/2611可编辑超前滞后校正器超前校正器2024/1/2612可编辑2024/1/2613可编辑2024/1/2614可编辑基于模型匹配算法的设计方法假设受控对象的传递函数为，期望闭环系统的频域响应为，超前滞后校正器的一般形式为使得在频率段内闭环模型对期望闭环模型匹配指标为最小2024/1/2615可编辑提出了下面的设计算法其中2024/1/2616可编辑其中，gp和f分别为受控对象和期望闭环系统的传递函数模型，w1和w2为需要拟合的频率段上下限。2024/1/2617可编辑【例6-2】受控对象模型为2024/1/2618可编辑6.1.3控制系统工具箱中的设计界面控制器设计界面界面允许选择和修改控制器的结构，允许添加零极点，调整增益，从而设计出控制器模型。2024/1/2619可编辑【例6-3】受控对象和控制器的传递函数模型分别为2024/1/2620可编辑6.2基于状态空间模型的控制器设计方法6.2.1状态反馈控制2024/1/2621可编辑将代入开环系统的状态方程模型，则在状态反馈矩阵下，系统的闭环状态方程模型可以写成如果系统完全可控，则选择合适的矩阵，可以将闭环系统矩阵的特征值配置到任意地方。2024/1/2622可编辑6.2.2线性二次型指标最优调节器假设线性时不变系统的状态方程模型为设计一个输入量,使得最优控制性能指标最小2024/1/2623可编辑则控制信号应该为由简化的Riccati微分方程求出假设，其中，则可以得出在状态反馈下的闭环系统的状态方程为依照给定加权矩阵设计的LQ最优控制器2024/1/2624可编辑离散系统二次型性能指标离散Riccati代数方程这时控制律为2024/1/2625可编辑【例6-4】2024/1/2626可编辑6.2.3极点配置控制器设计系统的状态方程为则系统的闭环状态方程为2024/1/2627可编辑2024/1/2628可编辑

Bass-Gura算法2024/1/2629可编辑基于此算法编写的MATLAB函数2024/1/2630可编辑Ackermann算法其中为将代入得出的矩阵多项式的值鲁棒极点配置算法place()函数不适用于含有多重期望极点的问题acker()函数可以求解配置多重极点的问题2024/1/2631可编辑【例6-5】2024/1/2632可编辑【例6-6】2024/1/2633可编辑6.2.4观测器设计及基于观测器的调节器设计2024/1/2634可编辑2024/1/2635可编辑2024/1/2636可编辑【例6-7】2024/1/2637可编辑2024/1/2638可编辑带有观测器的状态反馈控制结构图2024/1/2639可编辑2024/1/2640可编辑2024/1/2641可编辑如果参考输入信号，则控制结构化简为2024/1/2642可编辑【例6-8】2024/1/2643可编辑2024/1/2644可编辑6.3过程控制系统的PID

控制器设计6.3.1PID控制器概述连续PID控制器2024/1/2645可编辑连续PID控制器Laplace变换形式2024/1/2646可编辑离散PID控制器2024/1/2647可编辑离散形式的PID控制器Z变换得到的离散PID控制器的传递函数2024/1/2648可编辑PID控制器的变形积分分离式PID控制器在启动过程中，如果静态误差很大时，可以关闭积分部分的作用，稳态误差很小时再开启积分作用，消除静态误差2024/1/2649可编辑离散增量式PID控制器2024/1/2650可编辑抗积分饱和(anti-windup)PID控制器2024/1/2651可编辑6.3.2过程系统的一阶延迟模型近似带有时间延迟一阶模型(first-orderlagplusdelay，FOLPD)一阶延迟模型(FOLPD)的数学表示为2024/1/2652可编辑由响应曲线识别一阶模型阶跃响应近似Nyquist图近似编写MATLAB函数getfolpd(),key=12024/1/2653可编辑基于频域响应的近似方法调用编写的MATLAB函数getfolpd(),key=22024/1/2654可编辑基于传递函数的辨识方法调用编写的MATLAB函数getfolpd(),key=32024/1/2655可编辑最优降阶方法调用编写的MATLAB函数getfolpd(),key=4【例6-9】2024/1/2656可编辑6.3.3Ziegler-Nichols参数整定方法

Ziegler-Nichols经验公式编写MATLAB函数ziegler()2024/1/2657可编辑【例6-10】2024/1/2658可编辑2024/1/2659可编辑改进的Ziegler-Nichols算法2024/1/2660可编辑初始点A增益期望点A1增益PID控制器2024/1/2661可编辑PI控制器2024/1/2662可编辑PID控制器2024/1/26632024/1/2664可编辑【例6-11】2024/1/2665可编辑2024/1/2666可编辑

改进PID控制结构与算法微分动作在反馈回路的PID控制器2024/1/2667可编辑精调的Ziegler-Nichols控制器及算法2024/1/2668可编辑2024/1/2669可编辑若则保留Ziegler-Nichols参数,同时为使超调量分别小于10%或20%，则若，

Ziegler-Nichols控制器的参数精调为若，为使系统的超调量小于10%，则PID参数调为：2024/1/2670可编辑【例6-12】用自编的MATLAB函数设计精调的Ziegler-NicholsPID控制器2024/1/2671可编辑改进的PID结构一种PID控制器结构及整定算法的控制器模型为:2024/1/2672可编辑6.3.4最优PID整定算法最优化指标时间加权的指标IAE和ITAE指标2024/1/2673可编辑庄敏霞与Atherton教授提出了基于时间加权指标的最优控制PID控制器参数整定经验公式适用范围，不适合于大时间延迟系统2024/1/2674可编辑Murrill提出了使得IAE准则最小的PID控制器算法2024/1/2675可编辑对ITAE指标进行最优化，得出的PID控制器设计经验公式在范围内设计的ITAE最优PID控制器的经验公式2024/1/2676可编辑【例6-13】2024/1/2677可编辑2024/1/2678可编辑6.3.5其他模型的PID控制器参数整定算法IPD模型的PD和PID参数整定(integratorplusdelay)2024/1/2679可编辑各种指标下的PD和PID参数整定公式若选择ISE指标,则若选择ITSE指标,则若选择ISTSE指标,则2024/1/2680可编辑编写设计控制器的MATLAB函数2024/1/2681可编辑FOLIPD模型的PD和PID参数整定(firstorderlagandintegratorplusdelay)

PID控制器的整定算法

PD控制器的设计算法2024/1/2682可编辑编写设计控制器的MATLAB函数2024/1/2683可编辑【例6-14】2024/1/2684可编辑不稳定FOLPD模型的PID参数整定

设计的PID控制器若使ISE指标最小,则若使ITSE指标最小,则若使ISTSE指标最小，则2024/1/2685可编辑不稳定FOLPD模型的PID控制器参数整定函数2024/1/2686可编辑6.3.6基于FOLPD的PID控制器设计程序在MATLAB提示符下输入pid_tuner。单击Plantmodel按钮，打开一个允许用户输入受控对象模型参数的对话框。输入了受控对象模型后，单击GetFOLPDparameters按钮获得FOLPD模型，亦即获得并显示K,L,T参数。2024/1/2687可编辑通过得出的K,L,T参数,设计所需的控制器。单击Designcontroller按钮,将自动设计出所需的PID控制器模型，并将其显示出来。单击Closed-loopSimulation按钮，则可以构造出PID控制器控制下的系统仿真模型，并在图形界面上显示系统的阶跃响应曲线。2024/1/2688可编辑6.4最优控制器设计6.4.1最优控制的概念在一定的具体条件下，要完成某个控制任务，使得选定指标最小或增大的控制.积分型误差指标、时间最短指标、能量最省指标等2024/1/2689可编辑【例6-16】设计最优控制器2024/1/2690可编辑为使得ITAE准则最小化，可以编写如下的MATLAB函数2024/1/2691可编辑2024/1/2692可编辑为了降低超调量，改进的仿真框图2024/1/2693可编辑2024/1/2694可编辑【例6-17】考虑前面的例子，假设可以接受的控制信号限幅值为202024/1/2695可编辑2024/1/2696可编辑6.4.2基于MATLAB/Simulink的最优控制程序及其应用最优控制器设计程序(OptimalControllerDesigner，OCD)的调用过程为：在MATLAB提示符下输入ocd。建立一个Simulink仿真模型，该模型至少包含待优化的参数变量和误差信号的准则。将对应的Simulink模型名填写到界面的SelectaSimulinkmodel编辑框中。2024/1/2697可编辑将待优化变量名填写到Selectvariablestobeoptimized编辑框中，且各个变量名之间用逗号分隔。估计指标收敛的时间段作为终止仿真时间,填写到Simulationterminatetime栏目中去。单击CreateFile按钮自动生成描述目标函数的MATLAB文件opt_*.m。单击Optimize按钮将启动优化过程。本程序允许用户指定优化变量的上下界，选择优化参数的初值，选择不同的寻优算法，选择离散仿真算法等。2024/1/2698可编辑【例6-18】受控对象的模型为

用最优控制器设计程序选择PID控制器参数。2024/1/2699可编辑自动生成目标函数的MATLAB：2024/1/26100可编辑【例6-19】用OCD同时设计串级控制器2024/1/26101可编辑Simulink仿真模型2024/1/26102可编辑【例6-20】对模型采用ISE准

则设计最优控制器。2024/1/26103可编辑6.4.3最优控制程序的其他应用【例6-21】对模型采用ITAE准则，用OCD来进行最优降阶研究。2024/1/26104可编辑6.5多变量系统的频域

设计方法逆Nyquist阵列方法特征轨迹法(characteristiclocusmethod)反标架坐标法(reversed-framenormalisation，RFN)序贯回路闭合方法(sequentialloopclosing)参数最优化方法(parametersoptimisationmethod)2024/1/26105可编辑6.5.1对角占优系统与伪对角化为预补偿矩阵，它使得为对角占优矩阵。对所得对角占优矩阵作动态的补偿。2024/1/26106可编辑由以下步骤求取最优的补偿矩阵：选择一个函数的频率点，求出系统的逆Nyquist阵列。对各个

值，构成一个矩阵，其中假设在

频率处的系统传递函数矩阵的逆Nyquist阵列表示为2024/1/26107可编辑求取矩阵的特征值与特征向量，并将最小特征值的特征向量记作。由上面的各个值得出的最小特征向量可以构成补偿矩阵选择个频率点，并假设对第个频率点引入加权系数，按照如下的方法构造矩阵2024/1/26108可编辑由MATLAB编写出为对角化函数pseudiag()2024/1/26109可编辑【例6-22】2024/1/26110可编辑2024/1/26111可编辑【例6-23】2024/1/26112可编辑引入动态补偿矩阵2024/1/26113可编辑利用Simulink模型，绘制系统的阶跃响应曲线2024/1/26114可编辑2024/1/26115可编辑6.5.2多变量系统的参数最优化设计系统的闭环传递函数矩阵2024/1/26116可编辑控制器参数的最小二乘解202