深入浅出PID算法

前言

博主是工业互联网行码农一枚,虽然不是算法工程师和自动化方向的,但经常参加同事介绍控制算法原理的培训,慢慢的对小部分控制算法有一定了解,其中使用频率最高的控制算法非PID莫属。很多同学在学习PID的时候,会被繁杂的数学公式吓倒,今天我们就抛开数学公式,用逻辑和例子给大家讲明白“到底什么是PID”?

PID算法简介

PID是一种控制算法,是 Proportional(比例)、Integral(积分)、Differential(微分)的缩写。它是连续系统中技术最为成熟、应用最为广泛的一种控制算法。

PID控制器主要适用于基本上线性,且动态特性不随时间变化的系统。简单来说就是:类似于需要将某一个物理量“保持稳定”的场合,PID基本都能派上用场。

工业中PID典型的应用场景有:温度控制、流量控制、液位控制等。

image.png
生活中PID的应用也很常见:自来水的压力控制、空调的温度控制、平衡车的平衡控制、汽车的定速巡航控制、无人机的悬停控制、火箭飞机的姿态调整等。

image.png

实际运行经验和理论分析都表明,运用PID控制算法对许多工业过程进行控制都能得到比较满意的效果。

PID的概念及公式

PID的基本思路是根据偏差量的大小,运用比例、积分、微分计算出一个控制量,将这个控制量输入被控制的系统,系统接收到该输入量后会输出一个相应的输出量,PID控制器再检测该输出量,并再计算偏差,然后再循环以上过程。
image.png
原始公式如下:
u ( t ) = K p ( e ( t ) + 1 T i ∫ 0 t e ( t ) d t + T d d e ( t ) d t ) u(t)=K_p \Bigg( e(t)+\frac{1}{T_i}\int^t_0e(t)dt+T_d\frac{de(t)}{dt}\Bigg) u(t)=Kp(e(t)+Ti10te(t)dt+Tddtde(t))

  • K p K_p Kp —— 控制器的比例系数;
  • T i T_i Ti —— 积分时间常数;
  • T d T_d Td —— 微分时间常数;
  • u ( t ) u(t) u(t) —— PID控制器的输出信号;
  • e ( t ) e(t) e(t) —— 给定值 r ( t ) r(t) r(t)与测量值之差。

公式可简化为:
u ( t ) = K p ∗ e ( t ) + K i ∗ ∑ n = 0 t e ( t ) + K d ∗ ( e ( t ) − e ( t − 1 ) ) u(t)=K_p * e(t) + K_i * \sum_{n=0}^t e(t) + K_d * (e(t) - e(t-1)) u(t)=Kpe(t)+Kin=0te(t)+Kd(e(t)e(t1))

  • K p K_p Kp —— 控制器的比例系数;
  • K i K_i Ki —— 控制器的积分系数;
  • K d K_d Kd —— 控制器的微分系数;
  • u ( t ) u(t) u(t) —— PID控制器的输出信号;
  • e ( t ) e(t) e(t) —— 给定值 r ( t ) r(t) r(t)与测量值之差。
  • e ( t − 1 ) e(t-1) e(t1) —— 上一次给定值 r ( t ) r(t) r(t)与测量值之差。

很多同学可能会被上面数学公式整懵了,没关系,我们先看完下面例子,再回头看来公式,相信你能对以上公式有新的认识。

在学习PID之前,我们需要先了解什么是 开环控制闭环控制,这2个概念能够帮助我们更好的理解PID。

开放回路控制系统

开环控制(Open Loop Control System):不将控制的结果反馈回来影响当前控制的系统。例如:

  • 开关——按下开关后的一瞬间,控制活动已经结束,灯是否亮起已对按开关的这个活动没有影响;
  • 投篮——篮球出手后就无法再继续对其控制,无论球进与否,球出手的一瞬间控制活动即结束。

闭环回路控制系统

闭环控制(Closed Loop Control System):需要将控制的结果反馈回来与希望值比较,并根据它们的误差调整控制作用的系统。例如:

  • 调节水龙头——首先在头脑中对水流有一个期望的流量,水龙头打开后由眼睛观察现有的流量大小与期望值进行比较,并不断的用手进行调节形成一个反馈闭环控制;
  • 骑自行车——同理,需要不断的修正行进的方向与速度形成闭环控制。

闭环控制系统通常会由以下6个环节组成。
image.png
我们拿“维持水缸水位高度在1米”为例,来详解闭环控制中的每个环节:

  • 传感器:人工测量当前水位高度
  • 目标量:维持水缸水位高度在1米
  • 偏差量目标量 - 当前水位高度
  • 控制器:根据偏差量计算出执行量
  • 执行量:传给执行器的入参
  • 执行器:人工用水桶向水缸中加水

其中传感器目标量偏差量执行量执行器 这5个环节都比较简单,一眼就能看明白意思。最关键的环节控制器相对复杂一些,在控制器环节选择不同的控制算法,根据偏差量计算出执行量也不同,今天我们就具体看看PID的控制效果。

PID之比例P

前面有说到PID的基本思路是根据偏差量的大小,运用比例P积分I微分D计算出一个控制量。但PID的3个参数,并不是非要一起使用,可以单用比例P来控制,也可以两两联合用比例P+积分I比例P+微分D来控制。

我们先看看仅用比例P能带来什么样的控制效果?继续上面的“水缸”例子:

  • 传感器:人工测量水缸初始水位高度是0.2米
  • 目标量:维持水缸水位高度在1米
  • 偏差量:1米 - 0.2米 = 0.8米
  • 控制器:使用PID算法的比例控制,根据偏差量计算出执行量
  • 执行量:传给执行器的入参
  • 执行器:人工用水桶向水缸中加水

通常情况偏差量执行量之间的单位不同,在当前例子里偏差量的单位是水缸,而执行量的单位是水桶,他们的容积不一样,所以我们需要一个系数,来放大或者缩小两者的关系。具体公式可抽象为:

执行量 = 比例P = 偏差量 * 比例P系数

这个时候,假设水缸旁边站着一个人,用水桶往水缸里加水来控制水位的高度。如果单纯的用比例控制算法,假设比例P系数是0.5(相当于两桶水的容量等于一缸水的容量),我们开始模拟加水实验:
image.png
如上图结果所示,人工用水桶向水缸中加水8次,就能把水缸水位加到1米高度。把当前水位转换成曲线,大概效果如下。

上述例子比较简单,单靠比例P就能完成任务,但现实情况往往只有比例P是不够的,我们再来看一个更复杂的情况,假设比例P系数仍然是0.5,但在每次加水的间隔,水缸都会漏掉0.1米高度的水。我们再从头模拟加水实验,看看具体结果会怎样?
image.png
如上图结果所示,当水位达到0.8之后,水位就不会继续增加了。因为,当水位等于0.8时,偏差量是0.2,每次往水缸中加水的量为执行量= 比例P= 0.2 * 0.5 = 0.1,而每次加水的间隔,水缸都会漏掉0.1米高度的水,所以加入的水和流出的水相抵消。

虽然,此时控制系统已达稳定状态,但实际值目标量之间的会存在一个稳定差值,这个差值叫稳态误差稳态误差非常常见,比如:控制空调温度会因为空气温差而降温、控制无人机固定高度会受重力影响往下掉、控制汽车定速巡航会有空气阻力和摩擦力的影响而减速,这些场景都会产生稳态误差

单单只用比例P控制闭环回路,是无法避免稳态误差问题的,因为比例P系数无法根据时间或次数累加,当存在外界干扰因素时,比例P系数无法动态调整大小,那么稳态误差就会一直存在。

假如我们引入时间的维度,就能获得2个神器 增幅器-积分I抑制器-微分D,他们分别解决比例P过小和过大的情况,比例P过小的话由增幅器-积分I 补充即可解决稳态误差问题,比例P过大由抑制器-微分D消减来防止过度震荡。

PID之积分I(增幅器)

我们先来看看增幅器积分I,它也可以理解为累加经验,当比例P过小,可以由积分I补充,它的原理是利用过去的时间不断累加。具体公式可抽象为:

比例P = 偏差量 * 比例P系数
积分I = 上一次积分I + 偏差量 * 积分I系数
执行量 = 比例P + 积分I

想象一下,如果用只比例P控制存在稳态误差,说明比例P过小,这时积分I就会不断累加到一个很大的值,来补充比例P,从而影响执行量

如果最终控制效果在目标量附近抖动,我们就能得到一个正负交替的偏差量,会在目标附近不断产生正负数累加到积分I积分I就会不断趋近于零,最终使控制效果趋于稳定。

我们继续用上面“水缸加水”的例子,假设目标量是1,每次漏水0.1,比例P系数为0.5,积分I系数为0.2,我们再从头模拟加水实验,具体数值如下。
image.png
把当前水位转换成曲线,大概效果如下。
image.png
上一轮实验只用比例P控制存在稳态误差。本轮实验加入了积分I之后就有了累加效果,在未达到控制效果之前积分I会持续累加,在达到目标量之后积分I因为惯性会继续过量控制 ,同时偏差量会由正转负再转正,积分I也会由正转负再转正,最终积分I会持续抵消掉每次漏水0.1,控制效果趋于稳定 。

增幅器有他的危险性,如果系统出现意外或错误,增幅器可能会被累加到无限大,导致系统不可用,所以增幅器需要有一定的限制。

  • 1、限制幅度,在任意时刻都给积分I设定最大值和最小值。
  • 2、不运行时清零,当系统判断没有运行时,主动将积分I清零。

PID之微分D(抑制器)

我继续再看看抑制器微分D,它也可以理解为预测未来,用当前的偏差量上一次偏差量,得到的结果就可能是下一次偏差量,用下一次偏差量提前参与到计算中,就可以防止执行量过大,产生超出目标量的问题。具体公式可抽象为:

比例P = 偏差量 * 比例P系数
积分I = 上一次积分I + 偏差量 * 积分I系数
微分D = (偏差量 - 上一次偏差量) * 微分D系数
执行量 = 比例P + 积分I + 微分D

换个“汽车刹车”的例子,平稳驾驶的车辆,当发现前面有红灯时,为了使得行车平稳,基本上提前几十米就松油门踩刹车。当车辆离停车线非常近的时候,则使劲踩刹车使车辆停下来,整个过程可以看做一个加入微分D的控制策略。

可以看到,在刹车过程中,因为偏差量是越来越小的,所以 微分D= (偏差量-上一次偏差量)*微分D系数一定是负数,在控制中加入一个负数项,他存在的作用就是为了防止汽车由于刹车不及时而闯过停车线。

从常识上理解,越靠近停车线,就越应该踩深刹车,不能让车过线,所以这个 微分D的作用,可以理解为刹车。当车离停车线很近,并且车速还很快时,这个 微分D的绝对值(实际上是一个负数)就会很大,表示应该大力踩刹车尽快让车停下来。

再回到上面“水缸加水”的例子,当发现水缸里的水快要接近目标量时,加入 微分D可以减少过量加水的幅度,说白了就是减少控制过程中的震荡。假设目标量是1,每次漏水0.1,比例P系数为0.5,积分I系数为0.2,微分D系数为0.2,我们再从头模拟加水实验,具体数值如下。
image.png
上一轮实验使用比例P积分I联合控制,最高水位达到1.3,超过目标量之后,水位最低回落至0.83。本轮实验加入了抑制器微分D之后,最高水位仅达到1.23,超过目标量之后,水位最低回落至0.88,相比上一轮实验,本轮震荡幅度明显减小,这就是微分D的抑制作用。

分享一个动图,很好的展示了比例P积分I微分D的控制效果,其中红色虚线是目标量,曲线是用当前值的变化趋势。结合这个动图,大家再回想一下“水缸加水”的例子,控制效果是不是很相似呀~

PID原理总结

我们用了“水缸加水”的例子,详细解释了PID的比例P积分I微分D三个参数的控制原理,相信小伙伴们已经明白什么是PID了。大家再回过头看这个公式,是不是也觉得挺亲切的。
u ( t ) = K p ∗ e ( t ) + K i ∗ ∑ n = 0 t e ( t ) + K d ∗ ( e ( t ) − e ( t − 1 ) ) u(t)=K_p * e(t) + K_i * \sum_{n=0}^t e(t) + K_d * (e(t) - e(t-1)) u(t)=Kpe(t)+Kin=0te(t)+Kd(e(t)e(t1))

  • K p K_p Kp —— 控制器的比例系数;
  • K i K_i Ki —— 控制器的积分系数;
  • K d K_d Kd —— 控制器的微分系数;
  • u ( t ) u(t) u(t) —— PID控制器的输出信号;
  • e ( t ) e(t) e(t) —— 给定值 r ( t ) r(t) r(t)与测量值之差(偏差量);
  • e ( t − 1 ) e(t-1) e(t1) —— 上一次给定值 r ( t ) r(t) r(t)与测量值之差。

最后再总结一下:

  • 比例P积分I微分D都跟偏差量有关
  • 比例P取决于当前的偏差量
  • 积分I累计过去所有偏差量之和
  • 微分D预测下一时刻偏差量

所以,经常有人说比例P是现在,积分I是过去,微分D是未来,是不是有一种哲学的感觉,哈哈~

PID调参口诀

最后分享个PID调参口诀,写得挺好还挺押韵的~

参数整定找最佳,从小到大顺序查;
先是比例后积分,最后再把微分加;
曲线振荡很频繁,比例度盘要放大;
曲线漂浮绕大湾,比例度盘往小扳;
曲线偏离回复慢,积分时间往下降;
曲线波动周期长,积分时间再加长;
曲线振荡频率快,先把微分降下来;
动差大来波动慢。微分时间应加长;
理想曲线两个波,前高后低4比1;
一看二调多分析,调节质量不会低;
若要反应增快,增大P减小I;
若要反应减慢,减小P增大I;
如果比例太大,会引起系统振荡;
如果积分太大,会引起系统迟钝。

参考引用

如有任何问题或建议,欢迎前往公众号【劼哥舍】留言~

  • 15
    点赞
  • 64
    收藏
    觉得还不错? 一键收藏
  • 2
    评论
总所周知,PID算法是个很经典的东西。而做自平衡小车,飞行器PID是一个必须翻过的坎。因此本节我们来好好讲解一下PID,根据我在学习中的体会,力求通俗易懂。并举出PID的形象例子来帮助理解PID。 一、首先介绍一下PID名字的由来: P:Proportion(比例),就是输入偏差乘以一个常数。 I :Integral(积分),就是对输入偏差进行积分运算。 D:Derivative(微分),对输入偏差进行微分运算。 注:输入偏差=读出的被控制对象的值-设定值。比如说我要把温度控制在26度,但是现在我从温度传感器上读出温度为28度。则这个26度就是”设定值“,28度就是“读出的被控制对象的值”。然后来看一下,这三个元素对PID算法的作用,了解一下即可,不懂不用勉强。 P,打个比方,如果现在的输出是1,目标输出是100,那么P的作用是以最快的速度达到100,把P理解为一个系数即可;而I呢?大家学过高数的,0的积分才能是一个常数,I就是使误差为0而起调和作用;D呢?大家都知道微分是求导数,导数代表切线是吧,切线的方向就是最快到至高点的方向。这样理解,最快获得最优解,那么微分就是加快调节过程的作用了。 二、然后要知道PID算法具体分两种:一种是位置式的 ,一种是增量式的。在小车里一般用增量式,为什么呢?位置式PID的输出与过去的所有状态有关,计算时要对e(每一次的控制误差)进行累加,这个计算量非常大,而明显没有必要。而且小车的PID控制器的输出并不是绝对数值,而是一个△,代表增多少,减多少。换句话说,通过增量PID算法,每次输出是PWM要增加多少或者减小多少,而不是PWM的实际值。所以明白增量式PID就行了。 三、接着讲PID参数的整定,也就是PID公式中,那几个常数系数Kp,Ti,Td等是怎么被确定下来然后带入PID算法中的。如果要运用PID,则PID参数是必须由自己调出来适合自己的项目的。通常四旋翼,自平衡车的参数都是由自己一个调节出来的,这是一个繁琐的过程。本次我们可以不管,关于PID参数怎么确定的,网上有很多经验可以借鉴。比如那个经典的经验试凑口诀: 参数整定找最佳, 从小到大顺序查。 先是比例后积分, 最后再把微分加。 曲线振荡很频繁, 比例度盘要放大。 曲线漂浮绕大弯, 比例度盘往小扳。 曲线偏离回复慢, 积分时间往下降。 曲线波动周期长, 积分时间再加长。 曲线振荡频率快, 先把微分降下来。 动差大来波动慢, 微分时间应加长。 理想曲线两个波, 前高后低四比一。 一看二调多分析, 调节质量不会低。 四、接下来我们用例子来辅助我们把常用的PID模型讲解了。(PID控制并不一定要三者都出现,也可以只是PI、PD控制,关键决定于控制的对象。)(下面的内容只是介绍一下PID模型,可以不看,对理解PID没什么用) 例子:我们要控制一个人,让他一PID的控制方式来行走110步后停下来。 1)P比例控制,就是让他按照一定的比例走,然后停下。比如比例系数为108,则走一次就走了108步,然后就不走了。 说明:P比例控制是一种最简单的控制方式,控制器的输出与输入误差信号成比例关系。但是仅有比例控制时系统输出存在稳态误差。比如上面的只能走到108,无论怎样都走不到110。 2)PI积分控制,就是按照一定的步伐走到112步然后回头接着走,走到108步位置时,然后又回头向110步位置走。在110位置处来回晃荡几次,最后停在110步的位置。说明:在积分I控制中,控制器的输出与输入误差信号的积分成正比关系。对一个自动控制系统来说,如果在进入稳态后存在稳态误差,则称这个控制系统是有稳态误差的或简称有差系统。为了消除稳态误差,在控制器中必须引入“积分项”。积分项对误差的影响取决于时间的积分,随着时间的增加,积分项会增大。这样,即便误差很小,积分项也会随着时间的增加而加大,它推动控制器的输出增大,从而使稳态误差进一步减小,直到等于0。因此,比例+积分(PI)控制器可以使系统在进入稳态后无稳态误差。 3)PD微分控制,就是按照一定的步伐走到一百零几步后,再慢慢地走向110步的位置靠近,如果最后能精确停在110步的位置,就是无静差控制;如果停在110步附近(如109步或111步位置),就是有静差控制。 说明:在微分控制D中,控制器的输出与输入误差信号的微分(即误差的变化率)成正比关系。 自动控制系统在克服误差的调节过程中可能会出现振荡甚至失稳,原因是存在较大惯性组件(环节)或滞后组件,具有抑制误差的作用,其变化总是落后于误差的变化。解决的办法是使抑制误差作用的变化“超前”,即在误差接近于零时,抑制误差的作用就应该是零。这就是说,在控制器中仅引入“比例P”项往往是不够的,比例项的作用仅是放大误差的幅值,而目前需要增加的是“微分项”,它能预测误差变化的趋势。这样,具有比例+微分的控制器就能够提前使抑制误差的控制作用等于零,甚至为负值,从而避免了被控量的严重超调。所以对有较大惯性或滞后的被控对象,比例P+微分D(PD)控制器能改善系统在调节过程中的动态特性。 五、用小明来说明PID: 小明接到这样一个任务:有一个水缸有点漏水(而且漏水的速度还不一定固定不变),要求水面高度维持在某个位置,一旦发现水面高度低于要求位置,就要往水缸里加水。 小明接到任务后就一直守在水缸旁边,时间长就觉得无聊,就跑到房里看小说了,每30分钟来检查一次水面高度。水漏得太快,每次小明来检查时,水都快漏完了,离要求的高度相差很远,小明改为每3分钟来检查一次,结果每次来水都没怎么漏,不需要加水,来得太频繁做的是无用功。几次试验后,确定每10分钟来检查一次。这个检查时间就称为采样周期。 开始小明用瓢加水,水龙头离水缸有十几米的距离,经常要跑好几趟才加够水,于是小明又改为用桶加,一加就是一桶,跑的次数少了,加水的速度也快了,但好几次将缸给加溢出了,不小心弄湿了几次鞋,小明又动脑筋,我不用瓢也不用桶,老子用盆,几次下来,发现刚刚好,不用跑太多次,也不会让水溢出。这个加水工具的大小就称为比例系数。 小明又发现水虽然不会加过量溢出了,有时会高过要求位置比较多,还是有打湿鞋的危险。他又想了个办法,在水缸上装一个漏斗,每次加水不直接倒进水缸,而是倒进漏斗让它慢慢加。这样溢出的问题解决了,但加水的速度又慢了,有时还赶不上漏水的速度。于是他试着变换不同大小口径的漏斗来控制加水的速度,最后终于找到了满意的漏斗。漏斗的时间就称为积分时间。 小明终于喘了一口,但任务的要求突然严了,水位控制的及时性要求大大提高,一旦水位过低,必须立即将水加到要求位置,而且不能高出太多,否则不给工钱。小明又为难了!于是他又开努脑筋,终于让它想到一个办法,常放一盆备用水在旁边,一发现水位低了,不经过漏斗就是一盆水下去,这样及时性是保证了,但水位有时会高多了。他又在要求水面位置上面一点将水缸要求的水平面处凿一孔,再接一根管子到下面的备用桶里这样多出的水会从上面的孔里漏出来。这个水漏出的快慢就称为微分时间。 六、在代码中理解PID:(好好看注释,很好理解的。注意结合下面PID的公式) 首先看PID的增量型公式: PID=Uk+KP*【E(k)-E(k-1)】+KI*E(k)+KD*【E(k)-2E(k-1)+E(k-2)】 在单片机中运用PID,出于速度和RAM的考虑,一般不用浮点数,这里以整型变量为例来讲述PID在单片机中的运用。由于是用整型来做的,所以不是很精确。但是对于一般的场合来说,这个精度也够了,关于系数和温度在程序中都放大了10倍,所以精度不是很高,但是大部分的场合都够了,若不够,可以再放大10倍或者100倍处理,不超出整个数据类型的范围就可以了。一下程序包括PID计算和输出两部分。当偏差>10度时全速加热,偏差在10度以内时为PID计算输出。 程序说明:下面的程序,先看main函数。可知在对定时器0初始化后就一直在执行PID_Output()函数。在PID_Output()函数中先用iTemp变量来得到PID运算的结果,来决定是启动加热丝加热还是不启动加热丝。下面的if语句结合定时器来决定PID算法多久执行一次。PID_Operation()函数看似很复杂,其实就一直在做一件事:根据提供的数据,用PID公式把最终的PID值算出来。 [C] 纯文本查看 复制代码 ? 001 002 003 004 005 006 007 008 009 010 011 012 013 014 015 016 017 018 019 020 021 022 023 024 025 026 027 028 029 030 031 032 033 034 035 036 037 038 039 040 041 042 043 044 045 046 047 048 049 050 051 052 053 054 055 056 057 058 059 060 061 062 063 064 065 066 067 068 069 070 071 072 073 074 075 076 077 078 079 080 081 082 083 084 085 086 087 088 089 090 091 092 093 094 095 096 097 098 099 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 149 150 151 152 153 154 155 156 157 158 159 160 161 162 163 164 165 166 #include typedef unsigned char uChar8; typedef unsigned int uInt16; typedef unsigned long int uInt32; sbit ConOut = P1^1; //加热丝接到P1.1口 typedef struct PID_Value { uInt32 liEkVal[3]; //差值保存,给定和反馈的差值 uChar8 uEkFlag[3]; //符号,1则对应的为负数,0为对应的为正数 uChar8 uKP_Coe; //比例系数 uChar8 uKI_Coe; //积分常数 uChar8 uKD_Coe; //微分常数 uInt16 iPriVal; //上一时刻值 uInt16 iSetVal; //设定值 uInt16 iCurVal; //实际值 }PID_ValueStr; PID_ValueStr PID; //定义一个结构体,这个结构体用来存算法中要用到的各种数据 bit g_bPIDRunFlag = 0; //PID运行标志位,PID算法不是一直在运算。而是每隔一定时间,算一次。 /* ******************************************************** /* 函数名称:PID_Operation() /* 函数功能:PID运算 /* 入口参数:无(隐形输入,系数、设定值等) /* 出口参数:无(隐形输出,U(k)) /* 函数说明:U(k)+KP*[E(k)-E(k-1)]+KI*E(k)+KD*[E(k)-2E(k-1)+E(k-2)] ******************************************************** */ void PID_Operation(void) { uInt32 Temp[3] = {0}; //中间临时变量 uInt32 PostSum = 0; //正数和 uInt32 NegSum = 0; //负数和 if(PID.iSetVal > PID.iCurVal) //设定值大于实际值否? { if(PID.iSetVal - PID.iCurVal > 10) //偏差大于10否? PID.iPriVal = 100; //偏差大于10为上限幅值输出(全速加热) else //否则慢慢来 { Temp[0] = PID.iSetVal - PID.iCurVal; //偏差 PID.liEkVal[1]) //E(k)>E(k-1)否? { Temp[0] = PID.liEkVal[0] - PID.liEkVal[1]; //E(k)>E(k-1) PID.uEkFlag[0] = 0; //E(k)-E(k-1)为正数 } else { Temp[0] = PID.liEkVal[1] - PID.liEkVal[0]; //E(k) Temp[2]) //E(k-2)+E(k)>2E(k-1)否? { Temp[2] = (PID.liEkVal[0] + PID.liEkVal[2]) - Temp[2]; PID.uEkFlag[2]=0; //E(k-2)+E(k)-2E(k-1)为正数 } else //E(k-2)+E(k) PID.iCurVal(即E(K)>0)才进入if的, 那么就没可能为负,所以打个转回去就是了 */ /* ========= 计算KD*[E(k-2)+E(k)-2E(k-1)]的值 ========= */ if(PID.uEkFlag[2]==0) PostSum += Temp[2]; //正数和 else NegSum += Temp[2]; //负数和 /* ========= 计算U(k) ========= */ PostSum += (uInt32)PID.iPriVal; if(PostSum > NegSum) //是否控制量为正数 { Temp[0] = PostSum - NegSum; if(Temp[0] 0,才有必要减“1” uCounter++; if(100 == uCounter) { PID_Operation(); //每过0.1*100S调用一次PID运算。 uCounter = 0; } } } /* ******************************************************** /* 函数名称:PID_Output() /* 函数功能:PID输出控制 /* 入口参数:无(隐形输入,U(k)) /* 出口参数:无(控制端) ******************************************************** */ void Timer0Init(void) { TMOD |= 0x01; // 设置定时器0工作在模式1下 TH0 = 0xDC; TL0 = 0x00; // 赋初始值 TR0 = 1; // 开定时器0 EA = 1; // 开总中断 ET0 = 1; // 开定时器中断 } void main(void) { Timer0Init(); while(1) { PID_Output(); } } void Timer0_ISR(void) interrupt 1 { static uInt16 uiCounter = 0; TH0 = 0xDC; TL0 = 0x00; uiCounter++; if(100 == uiCounter) { g_bPIDRunFlag = 1; } } #include typedef unsigned char uChar8; typedef unsigned int uInt16; typedef unsigned long int uInt32; sbit ConOut = P1^1; //加热丝接到P1.1口 typedef struct PID_Value { uInt32 liEkVal[3]; //差值保存,给定和反馈的差值 uChar8 uEkFlag[3]; //符号,1则对应的为负数,0为对应的为正数 uChar8 uKP_Coe; //比例系数 uChar8 uKI_Coe; //积分常数 uChar8 uKD_Coe; //微分常数 uInt16 iPriVal; //上一时刻值 uInt16 iSetVal; //设定值 uInt16 iCurVal; //实际值 }PID_ValueStr; PID_ValueStr PID; //定义一个结构体,这个结构体用来存算法中要用到的各种数据 bit g_bPIDRunFlag = 0; //PID运行标志位,PID算法不是一直在运算。而是每隔一定时间,算一次。 /* ******************************************************** /* 函数名称:PID_Operation() /* 函数功能:PID运算 /* 入口参数:无(隐形输入,系数、设定值等) /* 出口参数:无(隐形输出,U(k)) /* 函数说明:U(k)+KP*[E(k)-E(k-1)]+KI*E(k)+KD*[E(k)-2E(k-1)+E(k-2)] ******************************************************** */ void PID_Operation(void) { uInt32 Temp[3] = {0}; //中间临时变量 uInt32 PostSum = 0; //正数和 uInt32 NegSum = 0; //负数和 if(PID.iSetVal > PID.iCurVal) //设定值大于实际值否? { if(PID.iSetVal - PID.iCurVal > 10) //偏差大于10否? PID.iPriVal = 100; //偏差大于10为上限幅值输出(全速加热) else //否则慢慢来 { Temp[0] = PID.iSetVal - PID.iCurVal; //偏差 PID.liEkVal[1]) //E(k)>E(k-1)否? { Temp[0] = PID.liEkVal[0] - PID.liEkVal[1]; //E(k)>E(k-1) PID.uEkFlag[0] = 0; //E(k)-E(k-1)为正数 } else { Temp[0] = PID.liEkVal[1] - PID.liEkVal[0]; //E(k) Temp[2]) //E(k-2)+E(k)>2E(k-1)否? { Temp[2] = (PID.liEkVal[0] + PID.liEkVal[2]) - Temp[2]; PID.uEkFlag[2]=0; //E(k-2)+E(k)-2E(k-1)为正数 } else //E(k-2)+E(k) PID.iCurVal(即E(K)>0)才进入if的, 那么就没可能为负,所以打个转回去就是了 */ /* ========= 计算KD*[E(k-2)+E(k)-2E(k-1)]的值 ========= */ if(PID.uEkFlag[2]==0) PostSum += Temp[2]; //正数和 else NegSum += Temp[2]; //负数和 /* ========= 计算U(k) ========= */ PostSum += (uInt32)PID.iPriVal; if(PostSum > NegSum) //是否控制量为正数 { Temp[0] = PostSum - NegSum; if(Temp[0] 0,才有必要减“1” uCounter++; if(100 == uCounter) { PID_Operation(); //每过0.1*100S调用一次PID运算。 uCounter = 0; } } } /* ******************************************************** /* 函数名称:PID_Output() /* 函数功能:PID输出控制 /* 入口参数:无(隐形输入,U(k)) /* 出口参数:无(控制端) ******************************************************** */ void Timer0Init(void) { TMOD |= 0x01; // 设置定时器0工作在模式1下 TH0 = 0xDC; TL0 = 0x00; // 赋初始值 TR0 = 1; // 开定时器0 EA = 1; // 开总中断 ET0 = 1; // 开定时器中断 } void main(void) { Timer0Init(); while(1) { PID_Output(); } } void Timer0_ISR(void) interrupt 1 { static uInt16 uiCounter = 0; TH0 = 0xDC; TL0 = 0x00; uiCounter++; if(100 == uiCounter) { g_bPIDRunFlag = 1; } }
2.1 机器人控制 机器人控制是PID控制在前沿科技领域的应用之一。通过PID控制技术,可以实现机器人在运动、定位、姿态控制等方面的精准控制。例如,在工业生产线上,机器人通常需要完成复杂的动作和操作,如对零部件进行装配、焊接等。这些动作需要机器人具有高精度的运动控制能力,以完成高质量的生产任务。 2.2 无人驾驶车辆控制 无人驾驶车辆是近年来快速发展的一项技术。PID控制在无人驾驶车辆的控制中起到关键作用。无人驾驶车辆需要实时地控制车辆的速度、转向角度以及刹车等操作,以实现安全高效的行驶。PID控制可以通过对车辆的传感器数据进行实时分析,从而对车辆的运动状态进行精确控制。 2.3 航空航天控制 PID控制在航空航天领域的应用也非常广泛。例如,在飞机的自动驾驶系统中,PID控制可以实现对飞机的空速、高度、航向和俯仰角等方面的控制。通过PID控制技术,可以实现飞机在空中的稳定飞行和精确导航。 2.4 人工智能控制 人工智能控制是当前前沿科技领域的热门话题之一。PID控制在人工智能控制中也有着广泛的应用。例如,在机器学习领域,PID控制可以用来优化神经网络的训练过程,从而提高机器学习的效率和准确性。此外,在机器人控制、自动化控制等领域,PID控制也可以用来实现人工智能系统的精确控制。
评论 2
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值