PID控制的MATLAB仿真(1)

PID控制的MATLAB仿真(1)PID控制应用广泛,说起来PID特别的简单,在比例负反馈的基础上加入微分项实现快速调节,加入积分项实现无静差。MATLAB强大的功能让学习变得非常简单。在学习智能控制这门课的时候,老师推荐了刘金琨老师的先进PID控制一书,边看边学边调程序。下文程序选自此书。首先从模拟PID开始,被控系统一般是下图这种结构:其中PID控制器:一种用模拟PID控制的磁悬浮:电路图左侧是霍尔传

大家好,又见面了,我是你们的朋友全栈君。

1. 简介

PID控制应用广泛,说起来PID特别的简单,在比例负反馈的基础上加入微分项实现快速调节,加入积分项实现无静差。MATLAB强大的功能让学习变得非常简单。在学习智能控制这门课的时候,老师推荐了刘金琨老师的先进PID控制一书,边看边学边调程序。下文程序选自此书。

2. 模拟PID

首先从模拟PID开始,被控系统一般是下图这种结构:
这里写图片描述

其中PID控制器:
这里写图片描述

一种用模拟PID控制的磁悬浮:
这里写图片描述
这里写图片描述
电路图
这里写图片描述
左侧是霍尔传感器获得磁体的位置,中间是控制器,右侧是驱动部分。U2D是比例放大器,后面的网络引入了微分。

对系统的仿真,可以使用SIMULINK,还可以通过S函数实现:

%S-function for continuous state equation
function [sys,x0,str,ts]=s_function(t,x,u,flag)
switch flag,
%Initialization
  case 0,
    [sys,x0,str,ts]=mdlInitializeSizes;
%Outputs
  case 3,
    sys=mdlOutputs(t,x,u);
%Unhandled flags
  case { 
   2, 4, 9 }
    sys = [];
%Unexpected flags
  otherwise
    error(['Unhandled flag = ',num2str(flag)]);
end
%mdlInitializeSizes
function [sys,x0,str,ts]=mdlInitializeSizes
sizes = simsizes;
sizes.NumContStates  = 0;
sizes.NumDiscStates  = 0;
sizes.NumOutputs     = 1;
sizes.NumInputs      = 3;
sizes.DirFeedthrough = 1;
sizes.NumSampleTimes = 0;
sys=simsizes(sizes);
x0=[];
str=[];
ts=[];
function sys=mdlOutputs(t,x,u)
error=u(1);
derror=u(2);
errori=u(3);

kp=60;
ki=1;
kd=3;
ut=kp*error+kd*derror+ki*errori;
sys(1)=ut;

3. 离散系统PID

有了方便的计算机,很多时候数字控制器变得特别方便。
对一个系统

1/(1+s)

%Discrete PID control for continuous plant
clear all;
close all;

ts=0.001;  %Sampling time
xk=zeros(2,1);
e_1=0;
u_1=0;

for k=1:1:2000
time(k) = k*ts;

yd(k)=0.50*sin(1*2*pi*k*ts);

para=u_1;
tSpan=[0 ts];
[tt,xx]=ode45('chap1_6plant',tSpan,xk,[],para);
xk = xx(length(xx),:);
y(k)=xk(1); 

e(k)=yd(k)-y(k);
de(k)=(e(k)-e_1)/ts; 

u(k)=20.0*e(k)+0.50*de(k);
%Control limit
if u(k)>10.0
   u(k)=10.0;
end
if u(k)<-10.0
   u(k)=-10.0;
end

u_1=u(k);
e_1=e(k);

D=1;
if D==1  %Dynamic Simulation Display
    plot(time,yd,'b',time,y,'r');
    pause(0.00000000000000000);
end
end
figure(1);
plot(time,yd,'r',time,y,'k:','linewidth',2);
xlabel('time(s)');ylabel('yd,y');
legend('Ideal position signal','Position tracking');
figure(2);
plot(time,yd-y,'r','linewidth',2);
xlabel('time(s)'),ylabel('error');
版权声明:本文内容由互联网用户自发贡献,该文观点仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请联系我们举报,一经查实,本站将立刻删除。

发布者:全栈程序员-站长,转载请注明出处:https://javaforall.net/142051.html原文链接:https://javaforall.net

(0)
上一篇 2022年5月24日 下午3:20
下一篇 2022年5月24日 下午3:20


相关推荐

  • 基于opencv人脸识别

    基于opencv人脸识别一、准备所有代码上传至https://pan.baidu.com/s/1ETUetRnU83iytb6ZcE5GUQ提取码:0jei环境:win10,vs2013,opencv3.2注:每个代码都可以在vs中单独运行该项目所有文件如下:其中只需要使用的文件如下:详细代码://这个是生成照片的代码#include<opencv2…

    2022年6月4日
    46
  • PyCharm使用lxml库报错(ModuleNotFoundError: No module named ‘lxml’)解决办法

    PyCharm使用lxml库报错(ModuleNotFoundError: No module named ‘lxml’)解决办法1 在学习爬虫时遇到 ModuleNotFou Nomodulename lxml 错误 2 起初考虑是自己的 lxml 库有问题 卸载以后重新安装 pipinstalllx ihttp pypi douban com simple trusted hostpypi douban com3 重新安装 lxml 库时 错误没有解决 百度各种方式后无果 经

    2026年3月18日
    3
  • 上帝掷骰子吗–量子物理史话

    上帝掷骰子吗–量子物理史话  上帝掷骰子吗–量子物理史话    第一章黄金时代    一    我们的故事要从1887年的德国开始。位于莱茵河边的卡尔斯鲁厄是一座风景秀丽的城市,在它的城中心,矗立着著名的18世纪的宫殿。郁郁葱葱的森林和温暖的气候也使得这座小城成为了欧洲的一个旅游名胜。然而这些怡人的景色似乎没有分散海因里希?鲁道夫?赫兹(HeinrichRudolfHertz)的注意力:现在他正

    2022年6月7日
    37
  • Android自定义控件传送,使用 Leanback 传输控件

    Android自定义控件传送,使用 Leanback 传输控件Leanbackandr 库中具有新的播放控件 可改善用户体验 对于视频应用 传输控件支持使用前进 后退控件进行视频拖动 拖动时 屏幕上会显示缩略图来帮助浏览视频 该库中包含抽象类和开箱即用的实现 可为开发者提供更精细的控制 使用预构建实现时 您无需编写大量代码即可快速构建功能丰富的应用 如果需要进行更多自定义设置 您可以扩展该库的任何预构建组件 控件和播放器 Leanback 库将带有传

    2026年3月17日
    3
  • PHP加密技术 附源码

    PHP加密技术 附源码

    2022年2月13日
    56
  • 解决: java.net.ConnectException: Connection refused: connect

    解决: java.net.ConnectException: Connection refused: connectjava.net.ConnectException:Connectionrefused:connect1.报错如上。2.原因和解决:1)端口号被占用:杀死占用端口号(一般为8080)的进程,释放端口;彻底关闭Tomcat服务器,或者重新启动项目。2)使用以下命令检查服务器中允许的最大进程数:ulimit-u。如果将此值设置为较低值,例如1024,则使用以下…

    2022年4月30日
    44

发表回复

您的邮箱地址不会被公开。 必填项已用 * 标注

关注全栈程序员社区公众号