MATLAB 控制系统设计与仿真 - 41

鲁棒控制的其他函数 - 回路成型函数 loopsyn

灵敏度问题由鲁棒控制工具箱中的loopsyn就可以直接求解,该函数采用H无穷回路成型算法设计控制器,函数的调用格式为:

Matlab 复制代码
[K,CL,gamma,info] = loopsyn(G,Gd) % G为受控对象模型
                                  % Gd为期望的回路传递函数
                                  % K为回路成型控制器模型
                                  % CL为在该控制器下的闭环系统模型
                                  % gamma为成型精度,当gamma=1时表示设计出精确的成型控制器
[K,CL,gamma,info] = loopsyn(G,Gd,alpha) % alpha为[0 1]之间的数,当alpha趋于0时,
                                        % 设计出的控制器具有较好的闭环控制效果
                                        % 当alpha趋于1时,控制器具有较好的抗干扰能力
[K,CL,gamma,info] = loopsyn(G,Gd,alpha,ord) %orde可以用来设定控制器的阶数

给定一下系统

选择两个回路的模型均为

试设计回路成型控制器。

实现的MATLAB代码如下:

Matlab 复制代码
clear all;clc;
g11=tf([0.806 0.264],[1 1.15 0.202]);
g12=tf([-15 -1.42],[1 12.8 13.6 2.36]);
g21=tf([1.95 2.12 0.49],[1 9.15 9.39 1.62]);
g22=tf([7.15 25.8 9.35],[1 20.8 116.4 111.6 18.8]);
G=[g11 g12; g21 g22];
s=tf('s');
Gd=500/(s+1); 
[K,CL,gamma]=loopsyn(G,Gd);
figure(1)
step(CL,0.1);
hold on
step(feedback(G*K,eye(2)),0.1)
hold off
legend('step reponse of CL','step response of G*K/(1+G*K)')
figure(2)
bodemag(G*K,[Gd 0; 0 Gd])

程序运行结果如下:

有结果可知,设计的效果一般,因为耦合度还是比较高的。

如果我们设计控制器时要求控制系统的系统优先于模型不确定性。可以做如下修改

Matlab 复制代码
[K,CL,gamma]=loopsyn(G,Gd,0.0001);
figure(1)
step(CL,0.1);
hold on
step(feedback(G*K,eye(2)),0.1)
hold off
legend('step reponse of CL','step response of G*K/(1+G*K)')
figure(2)
bodemag(G*K,[Gd 0; 0 Gd])

程序运行结果如下:

由上图可知,系统的耦合得到了大大改善,小伙伴们也可以自己画一下闭环系统的传递函数比较一下。但是在实际应用中,要考虑到系统的不确定性,否则系统会不稳定。

最后,欢迎大家有问题给我留言。

非常感谢小伙伴们的-点赞-收藏-加关注。

相关推荐
qingfeng15415几秒前
企业微信 API 可以做什么?
java·开发语言·python·自动化·企业微信
梧桐和风几秒前
2026 年 Java 趋势:AI 浪潮下,Java 会过时吗?
java·开发语言·人工智能
lsx2024061 分钟前
React 组件详解
开发语言
Hesionberger1 分钟前
LeetCode105:前序中序构建二叉树(三解法)
java·数据结构·python·算法·leetcode·深度优先
嗨嗨的迷子1 分钟前
JDK 17 远程调试连不上 5005:从 attach timeout 到 JDWP 监听地址变更
java·开发语言
@小柯555m3 分钟前
算法(移动零)
数据结构·算法·leetcode
Chase_______3 分钟前
【Java杂项】为什么 long 可以自动转 float?宽化基本类型转换与精度丢失详解
java·开发语言·python
listhi5204 分钟前
基于QT的串口心电波形实时显示系统
开发语言·qt
之歆4 分钟前
DAY_12JavaScript DOM 完全指南(三):高级工程篇
开发语言·前端·javascript·ecmascript
invicinble4 分钟前
java数组相关的信息量
java·开发语言·python