Ziegler-Nichols对PID的调参方法

PID控制器是工业过程控制中广泛采用的一种控制器,其中,P、I、D分别为比例(Proportion)、积分(Integral)、微分(Differential)的简写;将偏差的比例、积分和微分通过线性组合构成控制量,用该控制量对受控对象进行控制,称为PID算法 。

其传递关系为:

控制器的传递函数可写为:

差分方程为:

其中KP、KI、KD分别为比例系数、积分系数、微分系数。

比例系数KP :反应系统当前最基本的误差,系数大,可以加快调节,减小误差,但是过大的比例使系统稳定性下降,甚至造成系统的不稳定。

积分系数KI  :反应系统的累计误差,使系统消除稳态误差,提高无差度,只要有误差,积分调节就会起作用。

微分系数KD :反应系统误差的变化率,具有预见性,们可以预见偏差的变化趋势,产生超前的控制效果。因此可以改善系统的动态性能。但是微分对噪声有放大作用,会减弱系统的抗干扰性。

(一)PID控制器参数整定的方法有很多,概括起来主要有两大类:一是理论计算整定法,二是通过在线实验的工程整定法。

  1. 理论计算整定法

它主要是依据被控对象准确的数学模型,经过理论计算确定控制器参数。这种方法一般较难做到,同时,得到的计算数据未必可以直接使用,还必须通过工程实际进行调整和修改。

  2.工程整定法

它不需要得到过程模型,主要依赖工程经验,在控制系统的试验中直接进行参数整定。方法简单实用,计算简便且易于掌握,可以解决一般的实际问题,在工程实际中被广泛采用。PID控制器参数的工程整定法,主要有临界比例度法(又称稳定边界法)、反应曲线法和4:1衰减法。其共同点都是通过实验,然后按照工程经验公式对控制器参数进行整定。然而,无论采用哪一种方法整定所得到的控制参数,都需要在实际运行中进行最后的调整与完善。
2.1试凑法就是通过仿真或实际运行,观察系统对典型输入作用的响应曲线,根据各控制参数对系统的影响,反复调节试凑,直到满意为止,从而确定PID参数。
  参数整定找最佳,从小到大顺序查,
  先是比例后积分,最后再把微分加,
  曲线振荡很频繁,比例度盘要放大,
  曲线漂浮绕大湾,比例度盘往小扳,
  曲线偏离回复慢,积分时间往下降,
  曲线波动周期长,积分时间再加长,
  曲线振荡频率快,先把微分降下来,
  动差大来波动慢,微分时间应加长,
  理想曲线两个波,前高后低四比一,
  一看二调多分析,调节质量不会低


2.2 Ziegler-Nichols方法 已经发展成为一种在参数设定中,处于经验和计算法之间的中间方法。这种方法可以为控制器确定非常精确的参数,在此之后也可进行微调。
Ziegler-Nichols基于时域响应曲线的整定
      1反应曲线法:用阶跃响应曲线来整定控制器的参数。
设想对被控对象(开环系统)施加一个阶跃信号,通过实验方法,测出其响应信号,根据这条阶跃响应曲线定出一些能反映控制对象动态特性的参数。如图所示,以曲线的拐点作一条切线得到三个参数:K是控制对象的增益,L是等效滞后时间,T是等效滞后时间常数。则输出信号可由图中的形状近似确定参数K,L和T(或α) ,其中 。如果获得了参数K,L和T(或α) 后,则可根据表A-1确定PID控制器的有关参数。

    2稳定边界法:用系统的等幅振荡曲线来整定控制器的参数。

先测出系统处于闭环状态下的对象的等幅振荡曲线,根据等幅振荡曲线定出一些能反映控制对象动态特性的参数。设系统为只有比例控制的闭环系统,则当KP增大时,闭环系统若能产生等幅振荡,如测出其振幅kp'和振荡周期p',然后由表A-2整定PID参数。

上述二法亦适用于系统模型已知的系统。但是此二法在应用中也有约束,因为许多系统并不与上述系统匹配,例如第一法无法应于开环传递中含积分项的系统,第二法就无法直接应用于二阶系统。如就无法利用Ziegler-Nichols法进行整定。

Ziegler-Nichols基于频域法的整定

如果实验数据是由频率响应得到的,则可先画出其对应的Nyquist图,从图中可以容易得到系统的剪切频率 与系统的极限增益 ,若令 ,同样我们从表3-3给出的经验公式可以得到PID控制器对应的参数。事实上,此法即时域法的第二法。

(二)选用Ziegler-Nichols基于时域响应曲线的整定反应曲线法
已知受控对象为传递函数 
传递函数在阶跃信号激励下测出其响应信号,根据这条阶跃响应曲线定出一些能反映控制对象动态特性的参数。如图 4所示,以曲线的拐点作一条切线得到三个参数

图 4开环阶跃响应曲线确定PID 参数

由图可知,取L=0.6 T=3.186 k =0.8000。于读图存在误差,因此参数仍需整定。已知对象的K、L 和T 值后,根据Ziegler — Nichols整定公式编写一个MATLAB函数

clc;clear

s=tf('s');clc;clear;

 num=[8];

den=conv([1 5],conv([1 2],[2 1]));

 G=tf(num,den)

 step(G,15);;hold on

 k=dcgain(num,den)

 L=0.6; T=3.186 ;K =0.8;

 a=K*L/T

PIDKp=1.2*T/(K*L); 

PIDTi=2*L;

PIDTd=0.5*L;

PIDG=PIDKp*(1+1/(PIDTi*s)+PIDTd*s)

step(feedback(PIDG*G,1)),hold on     

 

  • 0
    点赞
  • 12
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值