精品课设经典PID与专家PID控制的对比与分析

Posted 2023-02-14 昔时扬尘处

tags:

篇首语：本文由小常识网(cha138.com)小编为大家整理，主要介绍了精品课设经典PID与专家PID控制的对比与分析相关的知识，希望对你有一定的参考价值。

【精品课设】经典PID与专家PID控制的对比与分析（二）

1.本代码基于MATLAB2019a版本，低版本或者不同版本可能会报错，mdl文件或slx文件打开可能会失败；
2.如果运行时间过长，请观察设置参数是否一致。
3.本博客附上算法运行图并详细介绍，如果转载请注明出处；
4.如果本博客恰巧与您的研究有所关联，欢迎您的咨询qq1366196286 。

1经典PID控制的设计与仿真

经典PID控制的设计与仿真可参照上一篇博客，【精品课设】经典PID与模糊PID控制的对比与分析（一）

1.1 被控对象传递函数的设计

被控对象的传递函数为：

1.2 经典PID控制的仿真

Traditional_PID.m



clc,
close all;
clear all; %#ok<CLALL>
warning off;

rin = 1.0; % 设定系统的参考输入量-单位阶跃响应
ts=0.001; %设置系统采样的时间--运行时间1.0s

sys=tf(5.235e005,[1, 87.35, 1.047e004, 0]); %设置被控对象的传递函数

dsys=c2d(sys,ts,'tustin');    %离散化处理-Z域
[num,den]=tfdata(dsys,'v');


u_1=0.0;u_2=0.0;u_3=0.0; % 系统初始值的0初始化
y_1=0;y_2=0;y_3=0;

x=[0,0,0]';  %中间变量的初始化
e_1=0; %对误差的初始化
ec_1=0;

%% 未进行整定
kp0=0.90; % PID控制参数的大小
ki0=0.05;
kd0=0.2;

% 对系统进行建模
for k=1:1:1000
    
time(k)=k*ts;
% r(k)=sign(sin(2*pi*k*ts)); %%脉冲输入
r(k)=rin; %%单位阶跃输入

u(k)=kp0*x(1)+ki0*x(2)+kd0*x(3);

% if k==300     % Adding disturbance(1.0v at time 0.3s)施加干扰量
%    u(k)=u(k)+1.0;
% end

y(k)=-den(2)*y_1-den(3)*y_2-den(4)*y_3+num(1)*u(k)+num(2)*u_1+num(3)*u_2+num(4)*u_3;
e(k)=r(k)-y(k);
%%%%%%%%%%%%%%Return of PID parameters%%%%%%%%%%%%%%%
   u_3=u_2;
   u_2=u_1;
   u_1=u(k);
   
   y_3=y_2;
   y_2=y_1;
   y_1=y(k);
   
   x(1)=e(k);            % Calculating P
   x(2)=e(k)-e_1;        % Calculating D
   x(3)=x(3)+e(k)*ts;    % Calculating I

   ec_1=x(2);
   e_2=e_1;
   e_1=e(k);
end
%% 绘制响应输出
clf,
figure(1);
plot(time,r,'b',time,y,'r');
xlabel('time(s)');ylabel('rin,yout');

程序运行结果图如下：