找回密码
 注册
关于网站域名变更的通知
查看: 515|回复: 3
打印 上一主题 下一主题

PID算法~~比较形象的描述

[复制链接]

该用户从未签到

跳转到指定楼层
1#
发表于 2019-10-18 16:35 | 只看该作者 回帖奖励 |倒序浏览 |阅读模式

EDA365欢迎您登录!

您需要 登录 才可以下载或查看,没有帐号?注册

x
PID控制器% q* d6 }$ a- v  z

' [) c& @& r/ W7 X3 i) D7 N首先以一个通俗而又经典的小故事来介绍它。; c1 p6 Q+ ]$ o5 s5 \8 G$ n* |; Y
/************************小故事分界线*****************************/ - K0 S" F  b4 w/ L1 D: S
小明接到这样一个任务:
) }. p( a' ~" L+ E( }! \( O+ @; J& A: u有一个水缸点漏水(而且漏水的速度还不一定固定不变),
& l" Q2 L# J0 O  Q. X( o$ }3 R) e要求水面高度维持在某个位置,
4 }+ G% V8 t" `" h; H一旦发现水面高度低于要求位置,就要往水缸里加水。
4 S5 [. \  z5 J8 z& |5 v! \/ A& b9 R; B4 v
小明接到任务后就一直守在水缸旁边,2 r5 l; e2 a+ S& Z
时间长就觉得无聊,就跑到房里看小说了,
  I' Y/ D0 |. W$ x每30分钟来检查一次水面高度。水漏得太快,
; `% q3 V5 b8 I8 M每次小明来检查时,水都快漏完了,离要求的高度相差很远
5 ^' Q) _, w+ N: d+ h3 u  x( ^,小明改为每3分钟来检查一次,结果每次来水都没怎么漏6 W& @$ R5 k. P% p! O
,不需要加水,来得太频繁做的是无用功。几次试验后,! H1 ]! H; f4 y7 y8 H4 Q( P$ a
确定每10分钟来检查一次。这个检查时间就称为采样周期。+ p7 Y  x% c7 q7 l) N) H
$ X+ B" S7 q- O- S1 o" @
开始小明用瓢加水,水龙头离水缸有十几米的距离,% \2 }; t) t6 H9 M4 Y# A9 K
经常要跑好几趟才加够水,于是小明又改为用桶加,
4 f( Q4 Y( a8 b) Z, i. m9 d$ A" f, A一加就是一桶,跑的次数少了,加水的速度也快了,
& S. ]2 d# }7 k' t! S6 c4 h3 @( L但好几次将缸给加溢出了,不小心弄湿了几次鞋,小明又动脑筋,5 s1 q6 W- s* z1 N% I7 _5 ?
我不用瓢也不用桶,老子用盆,几次下来,发现刚刚好,不用跑太多次,; D" O0 t1 w9 @8 e( Y# [
也不会让水溢出。这个加水工具的大小就称为比例系数。 7 t7 E# b  P; I4 i4 H0 ?% o1 r

* s* e4 G" w6 e4 O) n: ^小明又发现水虽然不会加过量溢出了,有时会高过要求位置比较多5 t5 v4 r2 J$ r: p7 L. E; J8 @
,还是有打湿鞋的危险(还是说加一定次数后,水位一直有个误差??因为用加水过程中水有损耗,参见后面一篇文章)。他又想了个办法,在水缸上装一个漏斗,9 ~. l5 E/ {' c% I1 ?
每次加水不直接倒进水缸,而是倒进漏斗让它慢慢加。这样溢出的问题解决了,/ y3 H' v7 `4 R4 F1 S6 R
但加水的速度又慢了,有时还赶不上漏水的速度。
5 O- @; h& Y$ J  ^9 l2 d9 Z0 y于是他试着变换不同大小口径的漏斗来控制加水的速度# x. d* O9 V6 p5 m  M
,最后终于找到了满意的漏斗。漏斗的时间就称为积分时间 。; ~. B6 o- D* L: I! K9 u. i9 w% k
* W4 o* L4 E# u# d) f9 U
小明终于喘了一口,但任务的要求突然严了,0 K, z+ @4 N( `( v
水位控制的及时性要求大大提高,一旦水位过低,
4 a  |2 y% x# G) Q必须立即将水加到要求位置,而且不能高出太多,否则不给工钱。
1 t: }$ z3 A; x6 z& V小明又为难了!于是他又开努脑筋,终于让它想到一个办法,常放一盆备用水在旁边,9 ~& ^5 |: b4 e& z0 h
一发现水位低了,不经过漏斗就是一盆水下去,这样及时性是保证了,但水位有时会高多了。$ L5 M1 ~; F7 e" _' M5 F
他又在要求水面位置上面一点将水凿一孔,再接一根管子到下面的备用桶里这样多出的水会从上面的孔里漏出来。
& S; K4 I% T1 r# u' A这个水漏出的快慢就称为微分时间。   x2 v# _( z% |) c. R( F/ [
/************************小故事分界线*****************************/4 k( f, y- @7 N0 m, b# I. q6 C& M6 L
好,故事看完了,相信读者对PID算法有了初步认识,微分的比喻一点牵强,不过能帮助理解就行,接下来就深入一点进行讲解:9 f- ]) R4 |7 K% F' e+ ]6 s
PID控制器由比例单元(P)单元(I)单元(D),' A! ^: _" ~/ `7 i2 o, S3 L, n6 p
( W/ ^$ Q: o) q
①比例P控制& Y  m3 ?5 d3 Q: T' @3 Y; s
理解:比例控制中控制器的输出u与输入误差信号err成比例关系。假如某一刻t你的预期值是E,而实际值为A,那么误差值就为err(t)=E-A,这时候控制器输出就为u(t)=Kp*err(t)。
' ?$ h& V+ h% u) u7 r* k" T- n( ?*问:假设你有一个加减计算器,由于某种原因你每次只能加或减u(t)并且比例系数Kp=0.4,让你从0开始去逼近100。  L" r6 h9 k$ x; q8 r! O4 m- o
*解:显然,这时候你的预期值E为100,而实际值为0,那么误差值就为err(t1)=100-0,进行第一次加减运算(你按计算器加的值为u(t1)=0.4*100=40),好了,第一次计算器从0加到了40;
( x1 E1 D+ h/ g4 K7 t# C第二次,预期100,实际40,误差60,进行加减运算(你按计算机加的值为u(t2)=0.4*60=24),第二次计算器从40加到了64;
" m" G2 @+ o2 g+ U) N8 Z第三次,预期100,实际64,误差36,进行加减运算(你按计算器加的值为u(t3)=0.4*36=14.4),第三次计算器从64加到了78.4;
4 A" j4 x+ d" R2 ~6 t: h第四次从78.4加到87.04;0 w- \% ?4 G0 T3 `! ~
第五次87.04加到92.2249 F1 Z) S" N7 m, v4 i! E
.
, |$ T7 a! O6 L.
8 J7 H1 s2 B0 k$ v 1 V8 J  w0 {$ d2 c9 o( ?
由于比例系数较低,大家可以看到数据是慢慢地逼近100,实际上加减运算到第23次时,实际值A=99.999,显然比例控制已经到达一个比较可观的精度了。如果比例系数Kp设为1则可以一次运算到达100。! x, I+ G8 o' D' @. g1 B% f  n
! t$ p' F" C% y' T) q# O
②积分I控制
; ^& @; T" L1 [, c理解:在积分控制中控制器的输出与输入误差信号的积分成正比关系。如果系统进入稳态后存在稳态误差,必须引入积分环节以消除误差。积分项为误差对时间的积分,即便误差很小,随着时间的增加积分项也会增大,引导控制器增大输出使稳态误差减小,直至稳态误差为0。
9 _$ ?- c  j- v2 O
5 g) v9 v1 v! o  Y; b! J& C6 ^* g" x0 B③微分D控制
! r- E  X0 h6 o! b8 X8 L理解:在微分控制中,控制器的输出与输入误差信号的微分(即误差的变化率)成正比关系。误差变化越快,其微分绝对值越大。误差增大时,其微分为正;误差减小时,其微分为负;只有误差为常数时,微分才为0。总体来讲微分环节能预测变化的趋势来提前抑制误差直至微分为0。# r; X  m7 x) {+ _

' K1 E5 W& m" M% E: e# q$ f2 T/*****************************总结********************************/
$ z* h* Z4 S# l6 y5 u9 F% t6 E对于上述3种控制环节:7 {8 |+ P% j/ [  `0 E  S+ |
比例(P)+积分(I)控制器可以使系统在进入稳态后无稳态误差;9 s0 h7 |  q* d( U0 r5 ~0 Q
比例(P)+微分(D)控制器能改善系统在调节过程中的动态特性;- B- f% U9 G" L; a3 ?& ~3 D5 g
总的来说,PID控制其实是对反馈的误差进行控制的过程,如果误差为0,则比例环节不起作用,只有误差存在才进行比例调节;积分环节是用来消除稳态误差;微分环节根据误差信号的变化趋势进行超前调节。
( K: o: c' d6 ]# s( @5 x/ I/**************************************************************************/
' K+ F2 C. N7 p8 G1 y9 |; tPID基本原理讲解到此结束,接下来讲讲编程实现。) s3 S8 I9 H, g4 i; ?
原始公式是连续状态的,为了方便在处理器上实现,通常对PID连续系统离散化。
" e( D' v* W6 [( x6 W7 a' N假设采样周期间隔为T,则在KT时刻:. y1 ]7 }7 W5 ~$ f5 {: y- ?. c) ~( r6 U
偏差err(k)=E(k)-A(k);
: V0 i& y# s% z' I$ D积分环节用加和形式表示:all_err(k)=err(k)+err(K-1)+…err(0)(离散情况下积分其实就是做累加);' b' V4 |1 u$ Y3 {# w4 ^9 g
微分环节用斜率形式表示:[err(k)-err(k-1)]/T,在离散情况下,微分就是error的差值,就是t时刻和t-1时刻error的差,即u=kd*(error(t)-error(t-1)),其中的kd是一个系数项;
" B3 s  g7 N, B9 p比例系数:Kp;9 ^+ C/ k/ _/ r
积分系数:Ki=Kp*T/Ti;
, b7 L. Q# ]' Q% j微分系数:Kd=Kp*Td/T;/ S: r0 D0 B! m. h7 s* f2 q8 x5 S
从而得到PID离散表示形式:
2 X; y6 G$ p/ a5 [0 Hu(k)=Kp*err(k)+Ki*(all_err(k))+kd*(err(k)-err(k-1))
* C# k( A9 b; Y
& n; U9 b) q$ x% k+ v上式也是位置式PID的表达形式。5 N+ Q# U$ L3 f1 N2 P5 k! D5 U2 R
我们求出u(k-1),然后求出△u(k)=u(k)-u(k-1),就得到增量式PID的表达形式:
+ o9 f9 G9 T. G/ T) `" ^& O$ C, G% }1 w1 M0 q$ o# X
△u(k)=kp*(err(k)-err(k-1)) +ki*err(k)+kd*(err(k)-2err(k-1)+err(k-2))
0 M& n/ s, M$ S6 {' a/ b" L: @  g) s, u" T5 ~. ~; i

* s4 A. I5 |! P0 `4 _
- V5 a# f+ s1 i7 d0 u* Z/ I2 @; i" X+ m, R; {) L
既然表达式都有了,那现在就用C语言来实例化PID控制器:
. B- [$ `/ e6 M4 J: O9 k* Q" Mstruct _pid 6 S: a4 K, k! i- O2 B. M
{+ I: |7 y- E: m- j' U( B
        float ExpectedValue;//定义设定值
9 T. M- A* A: o) T3 I* p        float ActualValue;//定义实际值1 c& ~% W3 R- ]- t, d) o  [
        float err;//定义偏差值
. V& x1 A1 [& G: L9 @- i        float err_last;//定义上一个偏差值2 m) M2 z! }1 Y- N& a! H
        float err_prev;//定义前一个的偏差值7 A1 e& v, a( n2 I
        float Kp, Ki, Kd;//定义比例、积分、微分系数9 P4 f) F( v% ^1 D1 v
}pid;
/ G5 F$ c4 r) d& p4 \
1 E. ]' [, P# L5 G- Wvoid PID_Init()- f. r3 }% X) H! ^; u' x5 ^$ |
{
3 R2 L' i$ P. a' O3 X) }8 Q. X* q        pid.ExpectedValue = 0.0;, N2 @- V. G! l; F" z0 d3 G/ j
        pid.ActualValue = 0.0;
$ V3 B" t% g% P* M/ k        pid.err = 0.0;7 Y0 D( i2 \' J6 W
        pid.err_prev = 0.0;
, {- H* Y0 w. ^6 C        pid.err_last = 0.0;& I! C4 }$ W  k* L" ~+ x6 L
        pid.Kp = 0.4;
6 j" p: s5 B  u  ^  ^* J# \        pid.Ki = 0.2;
0 L$ x) ]0 E; F6 ~        pid.Kd = 0.0;
. u# j9 \! o, S% c; l! z# O7 q}
: j7 W% _$ O6 q 1 I) [  n4 x, f# Z6 L+ v4 S5 X
float PID_Realize(float speed) {- J  b( d6 a$ B" X
        float index;5 }1 a; S; w% _$ _
        pid.ExpectedValue = speed;8 ~- C! [+ z9 ~
        pid.err = pid.ExpectedValue - pid.ActualValue;- C/ c9 d9 w- q9 ~
        //增量式pid公式6 l; p3 x% l9 v7 l! t$ o* `
        float incrementValue = pid.Kp*(pid.err - pid.err_last) + pid.Ki*pid.err + pid.Kd*(pid.err - 2 * pid.err_last + pid.err_prev);
7 u5 D  R& U% V& r4 z        pid.ActualValue += incrementValue;* `9 R* M7 T: V" I) d6 W
        pid.err_prev = pid.err_last;
" }* p) a# Z4 e8 Z& \. p! s" ]3 a# f        pid.err_last = pid.err;
+ W) \% K6 q8 n, y6 Q% p8 S        return pid.ActualValue;
. g8 G, c) Q0 y  C, l}
% t) Q% a8 p& g- k  `2 _+ Y" r
$ D; A3 w& c$ q' E9 }- z( k5 \8 Y曲线1:设初始值为0,期望值ExpectedValue为100,比例系数Kp为0.5,积分系数Ki为0.2,微分系数kd为0.1,采样周期为20ms,让PID进行200次调节,得到如下曲线(data1):3 n" r! \; i+ x- Y2 v6 S

: B7 q/ H6 b5 T- i: N  v* i" X& f4 b  Z( B3 U" ]
设当实际值ActualValue达到199.999时认为稳态,通过查看数据得第84次运算时达到条件,所用时间为84*20(ms)=1680(ms)=1.68(s)
5 \4 w, {- M- `* L  B* O* ~3 p( C0 C  U0 e7 W3 k

) z) H0 {- S6 ~6 `  ~曲线2:设初始值为0,期望值ExpectedValue为100,比例系数Kp为0.5,积分系数Ki为0.3,微分系数kd为0.0,采样周期为20ms,让PID进行200次调节
) @# R( @# @3 E5 J( p& a7 p0 Z2 M+ F
; R2 P& S6 ?, w- C& @# a+ q: f# [设当ActualValue达到199.999时认为稳态,通过查看数据得第56次运算时达到条件,所用时间为56*20(ms)=1120(ms)=1.12(s),可见调参之后系统达到稳态时间缩减33.3%,并且调节过程更为稳定。# [% j* i4 Q( w9 j2 l5 c  ?. e4 `

9 I/ k1 i3 z( T6 ^5 Q*注:上面用反应曲线法进行粗略地参数整定,但PID控制器的参数整定往往需要实际工程操作进行调整修改才能得出最优解,这就要求大家要多动手实操了。+ h$ i3 W* f; M  D

' l% B9 f: l/ I- I: P; p. D! C
# q" M7 k$ x; ^* z
" b$ Q8 x9 n! w  p9 s- g2 s7 ]% ^

4 P8 q! W& I; e/ l, A, o$ z
8 j" o: y4 S; i
: S0 X  m6 I2 t& n
1 _. s+ x. E: e5 h$ I
  I) Q  l8 [* i

该用户从未签到

2#
 楼主| 发表于 2019-10-18 16:36 | 只看该作者
第一次发贴,自己顶一下

该用户从未签到

3#
发表于 2019-11-4 18:07 | 只看该作者
大牛真是太厉害了,通俗易懂

点评

别人写的,觉得很容易理解,就传上来了  详情 回复 发表于 2019-11-5 17:28

该用户从未签到

4#
 楼主| 发表于 2019-11-5 17:28 | 只看该作者
3290700544 发表于 2019-11-4 18:07% V9 H) I) z) m* C0 Y2 D! |3 R0 Z
大牛真是太厉害了,通俗易懂

2 R% {5 {- z- u4 I) I$ d. p别人写的,觉得很容易理解,就传上来了
您需要登录后才可以回帖 登录 | 注册

本版积分规则

关闭

推荐内容上一条 /1 下一条

EDA365公众号

关于我们|手机版|EDA365电子论坛网 ( 粤ICP备18020198号-1 )

GMT+8, 2025-10-7 05:50 , Processed in 0.140625 second(s), 24 queries , Gzip On.

深圳市墨知创新科技有限公司

地址:深圳市南山区科技生态园2栋A座805 电话:19926409050

快速回复 返回顶部 返回列表