标题: PID学习笔记(先翻译大神教程) [打印本页]

作者: 51黑ww    时间: 2016-4-11 16:24
标题: PID学习笔记(先翻译大神教程)
提高初学者PID原文地址:
http://brettbeauregard.com/blog/ ... s-pid-introduction/

提高初学者的PID:
这里是第一次接触PID需要学习的公式:

它能引导大多数人写出如下的PID控制器代码:
  1.     /*working variables*/
  2.     unsigned long lastTime;
  3.     double Input, Output, Setpoint;
  4.     double errSum, lastErr;
  5.     double kp, ki, kd;
  6.     void Compute()
  7.     {
  8.        /*How long since we last calculated*/
  9.        unsigned long now = millis();
  10.        double timeChange = (double)(now - lastTime);

  11.        /*Compute all the working error variables*/
  12.        double error = Setpoint - Input;
  13.        errSum += (error * timeChange);
  14.        double dErr = (error - lastErr) / timeChange;

  15.        /*Compute PID Output*/
  16.        Output = kp * error + ki * errSum + kd * dErr;

  17.        /*Remember some variables for next time*/
  18.        lastErr = error;
  19.        lastTime = now;
  20.     }

  21.     void SetTunings(double Kp, double Ki, double Kd)
  22.     {
  23.        kp = Kp;
  24.        ki = Ki;
  25.        kd = Kd;
  26.     }
复制代码

Compute()被称作定期或不定期的,它工作非常正常。虽然这个系列不是“工作的最好的”。如果我们想做出和工业PID控制器相近的驱动器,我们需要解决几个问题:
Sample Time(采样时间)——如果这是一个固定的时间间隔,PID算法的功能实现将是非常好的。如果已知了这个间隔时间,代码中也可以简化一些内部的数学运算。
Derivative Kick(微分的过冲)——不是最大的问题,但是很容易解决,所以我们也将处理这个问题。
On-The-Fly Tuning Changes——好的PID函数是当调整参数的时候不会干扰内部运算的。
Reset Windup Mitigation(缓解积分饱和)——我们将会了解什么是积分饱和,并且在有利的方向上进行解决方案的实施。
On/Off (Auto/Manual)(开关-自动或手动)——在大多数应用中,有时候我们希望关闭PID控制器手动调节输出而不受控制器的干涉。
Initialization(初始化)——当控制器打开的时候我们希望是“无扰切换”,即我们不希望输出值忽然变成一个新的值。
Controller Direction(控制器的方向)——这是最后一个不是在鲁棒本身名称下的变化。它是为了确保用户能输入正确的调优参数而设计的。

一旦我们解决了这些问题,我们将有一个对PID算法深刻的了解。我们还会拥有最新的Arduino PID控制库。所以不管你是想自己写出自己的PID算法还是想去了解PID算法里到底发生了什么,我希望这些都能帮上你。现在我们开始旅程吧。

提高初学者的PID——采样时间
初学者的PID被称作不规则的,这就有了以下两个问题:
》你没有从PID中得到一致的状态特性,因为有时它是非常快的变化,有时却没有。
》你需要额外的数学运算解决微分和积分,它们都同时依赖于时间的变化。

解决方法:
确保PID定义在一个固定的时间间隔里。我这样做的原因是让compute指令每个周期都被调用一次。根据之前设定好的采样周期,PID决定是该计算还是立刻返回值。

一旦我们知道了PID是在一个恒定的时间内运算,微分和积分也就变得简单了。

代码:
  1.     /*working variables*/
  2.     unsigned long lastTime;
  3.     double Input, Output, Setpoint;
  4.     double errSum, lastErr;
  5.     double kp, ki, kd;
  6.     int SampleTime = 1000; //1 sec
  7.     void Compute()
  8.     {
  9.        unsigned long now = millis();
  10.        int timeChange = (now - lastTime);
  11.        if(timeChange>=SampleTime)
  12.        {
  13.           /*Compute all the working error variables*/
  14.           double error = Setpoint - Input;
  15.           errSum += error;
  16.           double dErr = (error - lastErr);

  17.           /*Compute PID Output*/
  18.           Output = kp * error + ki * errSum + kd * dErr;

  19.           /*Remember some variables for next time*/
  20.           lastErr = error;
  21.           lastTime = now;
  22.        }
  23.     }

  24.     void SetTunings(double Kp, double Ki, double Kd)
  25.     {
  26.       double SampleTimeInSec = ((double)SampleTime)/1000;
  27.        kp = Kp;
  28.        ki = Ki * SampleTimeInSec;
  29.        kd = Kd / SampleTimeInSec;
  30.     }

  31.     void SetSampleTime(int NewSampleTime)
  32.     {
  33.        if (NewSampleTime > 0)
  34.        {
  35.           double ratio  = (double)NewSampleTime
  36.                           / (double)SampleTime;
  37.           ki *= ratio;
  38.           kd /= ratio;
  39.           SampleTime = (unsigned long)NewSampleTime;
  40.        }
  41.     }
复制代码

在第10和第11行,如果它的时间可以计算出来,那就将由算法本身决定。因为我们现在知道样本之间的时间是相同的,我们并不需要不断乘以时间的变化。我们只需要适当调整Ki和Kd(30和31行),虽然在数学上的结果是等价的,但更有效。

虽然这样做又一个小小的波动。如果用户在操作过程中决定改变采样时间,Ki和Kd将需要重新调整来以对这一新的变化做出反应。这就第39-42行代码所处理的问题。

另外请注意我在第29行将采样时间转换成秒s了。严格的的说这是没有必要的,只是允许用户输入的Ki和Kd是uint型的s而不是ms。








欢迎光临 (http://www.51hei.com/bbs/) Powered by Discuz! X3.1