PID的数学模型
' B s( l& y0 u: ^+ O, Y0 X1 B) T5 P
在工业应用中PID及其衍生算法是应用最广泛的算法之一,是当之无愧的万能算法,如果能够熟练掌握PID算法的设计与实现过程,对于一般的研发人员来讲,应该是足够应对一般研发问题了,而难能可贵的是,在很多控制算法当中,PID控制算法又是最简单,最能体现反馈思想的控制算法,可谓经典中的经典。经典的未必是复杂的,经典的东西常常是简单的,而且是最简单的。PID算法的一般形式:9 S9 ^6 H5 D3 g. V- r' |
# Z4 ^: |2 A7 g* e5 a/ x% d
PID算法通过误差信号控制被控量,而控制器本身就是比例、积分、微分三个环节的加和。这里我们规定(在t时刻):
4 P2 a5 |/ u0 b* U, v* B+ z
% r% `. E% h: L: |1.输入量为:O(t)
/ x( D9 `& Y- [. }3 V% I
# s; }" w/ ?/ m1 W. |2.输出量为:i(t)/ F* _/ B& \5 s1 e
; p2 N, y. Y) Z! X' a
3.偏差量为& f; v) M3 L- I0 \. T- {
& S' D# s: M- K* M
9 X6 I/ y0 ]4 S
4 W' U @! u. J7 }/ h( }, u$ X2 q6 q- I1 T) u$ P `" V
- ]* d3 E5 H3 q( o, L! W# V
% Y$ G& L; l) B3 T2 @9 }) ^# A
◎ PID算法的数字离散化" w, }' O& K- k5 O3 ]% l; i5 S
4 W9 i' c$ U5 S- r+ q
假设采样间隔为T,则在第K个T时刻:
# d, @, [( V J6 s" f9 m- N, T9 r3 c* f' M+ y2 V
偏差=# C( ]$ K3 |2 U6 X U# z
& t4 ?: v& e+ g" J1 a' K1 m% Z5 O( ]9 r) U
* Z! J/ \" n% a- a; ?积分环节用加和的形式表示,即% J. J: {" e; x( O+ R
1 g6 X2 [1 f. p/ Z4 N7 n6 P2 [
/ g x2 M0 h9 F
5 O( S0 ^" @% C" x1 g5 A
微分环节用斜率的形式表示,即
, S1 W/ ~ M( U, O; d
% e8 G( V0 X7 ^ t j: @$ \1 l( L4 k/ {! m
' b3 ]) F* m$ |
PID算法离散化后的式子:! E" ?% _3 z# x, c1 `; A
t& D, K( V1 X& p
; r& ]0 X5 H ]: }3 s5 n- ~, U
( ^# A2 [# k2 @2 z: R0 F+ ~则可表示成为:! ]' H0 @3 ~, _, G, S6 {$ w
! {) W; t! d4 w+ g( D5 T; {
! _! P9 v6 d+ Q
8 W% \& T' _& m. j* ^
7 H: N& B0 `$ w; b$ {% D& l( _; f
其中式中:# E. H! N$ ?4 X0 c- Q. t" l# W4 B
: D2 A1 j6 `2 t4 R* X
比例参数kp:控制器的输出与输入偏差值成比例关系。系统一旦出现偏差,比例调节立即产生调节作用以减少偏差。特点:过程简单快速、比例作用大,可以加快调节,减小误差;但是使系统稳定性下降,造成不稳定,有余差。- @+ r2 }$ a9 ^" b; k4 i' U' O
1 O2 b5 u# Y* i# Y+ T4 V6 {9 k, H积分参数ki:积分环节主要是用来消除静差,所谓静差,就是系统稳定后输出值和设定值之间的差值,积分环节实际上就是偏差累计的过程,把累计的误差加到原有系统上以抵消系统造成的静差。) H2 L2 l/ ^" D" F$ C
, X8 X6 [( H; |# O6 v( }* Y微分参数kd:微分信号则反应了偏差信号的变化规律,或者说是变化趋势,根据偏差信号的变化趋势来进行超前调节,从而增加了系统的快速性。& L# F! y. b/ f0 v0 `! Y; W
1 z( O9 \1 X6 O: q8 r* r
PID的基本离散表示形式如上。目前的这种表述形式属于位置型PID,另外一种表述方式为增量式PID,由上述表达式可以轻易得到:
* D- N" f% Q' n$ V' v
& p* ~3 |" a1 s3 m0 u0 q
! G( @, z# k5 m9 s$ o# R0 @* i2 D* C% a5 U6 g$ t9 V! u
' ?5 i& I) m8 j( M* ^那么:
5 i! M0 i+ H$ J( D3 y! \, k+ c7 z3 c, R1 f8 F, L, G$ s
: v z" y8 b+ b6 x8 O- u O! s8 g$ M1 ^* ]
; M9 @3 B( }8 s# v5 J; P
上式就是离散化PID的增量式表示方式,由公式可以看出,增量式的表达结果和最近三次的偏差有关,这样就大大提高了系统的稳定性。需要注意的是最终的输出结果应该为:& e- x: d7 d% M/ A9 n7 s# C
- R; p" D, t; l$ [输出量 = u(k)+ 增量调节值
7 h+ ~: z1 } G8 j) r% P+ B; q
◎ 目的
2 M7 `8 Y3 }* P% D# F( Q; N$ S7 Q" h% D# \: _. Z
PID 的重要性应该无需多说了,这个控制领域的应用最广泛的算法了。
: Q, y. e, C2 B3 O6 V( F
0 t0 ]$ {" ^8 q; D$ K, U本篇文章的目的是希望通过一个例子展示算法过程,并解释以下概念:
- ~2 ?! ^: c3 k0 U# o
* g8 N9 L6 G, M# R l) Q(1)简单描述何为PID, 为何需要PID,PID 能达到什么作用。0 x6 Z7 u7 L5 n- v- a" o$ i- u
8 |2 ~7 L7 ~" i: @ T+ u(2)理解P(比例环节)作用:基础比例环节。
- x) U9 {' ?2 U% Y( |8 I
2 [$ v: P, M# B2 b: w' ~ L缺点: 产生稳态误差.
2 z; x, y$ U1 F) x
1 h) w/ }( ?) n$ z7 T, i7 V疑问: 何为稳态误差 为什么会产生稳态误差.) ^- x* _9 f1 k& h
( W1 X- |+ c5 B& `7 ^8 l& k/ z(3)理解I(积分环节)作用:消除稳态误差./ n3 q. ]- C8 R/ k& ]/ {+ f2 X
# j" d. h7 c L- L5 h缺点: 增加超调: D7 ~, C, U# Y' T
1 y S! c2 u6 C
疑问: 积分为何能消除稳态误差?
5 n. V- j% g \7 |1 P! ^, E% P5 \8 W
& f/ Z5 ], ?3 M(4) 理解D(微分环节)作用:加大惯性响应速度,减弱超调趋势
9 J. ^, a( d' i7 m$ \
6 v6 T& ]% @, H& r疑问: 为何能减弱超调
7 }9 d6 J2 B6 R; R" f$ t
0 g2 p3 h) y7 S B(5)理解各个比例系数的作用
1 y( w" V1 L0 |9 ?& G" b" u& `
- z: g3 O3 C! y7 {7 B! o! o& x
h3 x2 k) d5 S4 L8 K4 D g4 `( w
& G4 a4 Y2 e7 b9 S+ a$ r
& ~2 I) e+ N! _5 A
何为PID以及为何需要PID?3 a; k( A$ k3 P
; @% `/ `/ z0 V/ c" D* h以下即PID 控制的整体框图,过程描述为:3 ?, r7 }$ T( h6 `! r
+ R0 F, ^0 t' s7 M; G
设定一个输出目标,反馈系统传回输出值,如与目标不一致,则存在一个误差,PID 根据此误差调整输入值,直至输出达到设定值.
! `! ^/ X6 c# U0 t: P1 K; @, d
1 z3 ^0 b/ [5 P3 X( r
8 V% i p+ M1 y8 i1 e' w1 c! h5 r! b0 h, g& N
$ h6 U2 n' B2 l& d8 m5 i
疑问:$ ]& l/ b+ Q( r8 Y7 s( q
) _9 f8 x/ O5 {
那么我们为什么需要PID 呢,比如我控制温度,我不能监控温度值,温度值一到就停止吗?
( R! X' z3 c* ^/ m/ g% L0 f5 W7 K
5 a; N3 x( U, n. G1 u- f6 A W, a+ j ~这里必须要先说下我们的目标,因为我们所有的控制无非就是想输出能够达到我们的设定,即如果我们设定了一个目标温度值,那么我们想要一个什么样的温度变化呢.
4 C3 ]6 K l6 f0 z6 _
# ?; v; x8 g4 O' L7 L8 b比如设定目标温度为30度, 目标无非是希望达到图1 希望其能够快速而且没有抖动的达到30度.
: {9 m+ a0 U# g( U
! q& G& _) F8 V6 l那这样大家应该就明白,如果使用温度一到就停止的办法,当然如果要求不高可能也行,当肯定达不到图1 这样的要求,因为温度到了后余温也会让温度继续升高.而且温度自身也会通过空气散热的.
8 {/ b* n1 ?* z- J: b4 N) L) B7 m1 A
% X% d6 j- ^; ]1 U" S
& z3 E( q7 T6 A& l3 V
图 系统输出的响应目标) Y: P$ e% S. K
5 W0 h8 ]8 q* L! P- ?综上所述,我们需要PID的原因无非就是普通控制手段没有办法使输出快速稳定的到达设定值。
2 @ ]% J$ u# g
! W5 @4 H3 F/ K+ |: G◎ 控制器的P,I,D项选择1 Y2 h. I* {1 ^& P
3 X5 N7 J) W' f: x6 s
下面将常用的各种控制规律的控制特点简单归纳一下:' q8 M5 g4 M" p. o# {5 Y! b! v
* {+ q% \8 z% ~# p3 o! H [$ B(1)、比例控制规律P:采用P控制规律能较快地克服扰动的影响,它的作用于输出值较快,但不能很好稳定在一个理想的数值,不良的结果是虽较能有效的克服扰动的影响,但有余差出现。它适用于控制通道滞后较小、负荷变化不大、控制要求不高、被控参数允许在一定范围内有余差的场合。如:金彪公用工程部下设的水泵房冷、热水池水位控制;油泵房中间油罐油位控制等。! L5 M, t6 ~" @! Y! A: H4 b
; ^0 a4 s c" t! {1 y
(2)、比例积分控制规律(PI):在工程中比例积分控制规律是应用最广泛的一种控制规律。积分能在比例的基础上消除余差,它适用于控制通道滞后较小、负荷变化不大、被控参数不允许有余差的场合。如:在主线窑头重油换向室中F1401到F1419号枪的重油流量控制系统;油泵房供油管流量控制系统;退火窑各区温度调节系统等。0 P3 X1 N, v. V) V! J! v; W
, U# ~4 f# _' n) i5 J0 c(3)、比例微分控制规律(PD):微分具有超前作用,对于具有容量滞后的控制通道,引入微分参与控制,在微分项设置得当的情况下,对于提高系统的动态性能指标,有着显著效果。因此,对于控制通道的时间常数或容量滞后较大的场合,为了提高系统的稳定性,减小动态偏差等可选用比例微分控制规律。如:加热型温度控制、成分控制。需要说明一点,对于那些纯滞后较大的区域里,微分项是无能为力,而在测量信号有噪声或周期性振动的系统,则也不宜采用微分控制。如:大窑玻璃液位的控制。/ k i T+ ?8 ?2 ]
$ F F* k! S; Q$ _5 [) b) z(4)、例积分微分控制规律(PID):PID控制规律是一种较理想的控制规律,它在比例的基础上引入积分,可以消除余差,再加入微分作用,又能提高系统的稳定性。它适用于控制通道时间常数或容量滞后较大、控制要求较高的场合。如温度控制、成分控制等。
2 U1 I; e, F% Z$ m4 }
' V; k v& [7 O; w' s8 W( F鉴于D规律的作用,我们还必须了解时间滞后的概念,时间滞后包括容量滞后与纯滞后。其中容量滞后通常又包括:测量滞后和传送滞后。测量滞后是检测元件在检测时需要建立一种平衡,如热电偶、热电阻、压力等响应较慢产生的一种滞后。而传送滞后则是在传感器、变送器、执行机构等设备产生的一种控制滞后。纯滞后是相对与测量滞后的,在工业上,大多的纯滞后是由于物料传输所致,如:大窑玻璃液位,在投料机动作到核子液位仪检测需要很长的一段时间。
$ L4 Z6 X0 i' W) a$ t, b( G, U
5 b! M6 [0 {5 X- ?: ~- z4 I K总之,控制规律的选用要根据过程特性和工艺要求来选取,决不是说PID控制规律在任何情况下都具有较好的控制性能,不分场合都采用是不明智的。如果这样做,只会给其它工作增加复杂性,并给参数整定带来困难。当采用PID控制器还达不到工艺要求,则需要考虑其它的控制方案。如串级控制、前馈控制、大滞后控制等。' p! D' w: t$ O c
, n- n4 P: E0 f3 ~$ j' w1 i* p
Kp,Ti,Td三个参数的设定是PID控制算法的关键问题。一般说来编程时只能设定他们的大概数值,并在系统运行时通过反复调试来确定最佳值。因此调试阶段程序须得能随时修改和记忆这三个参数。
4 E' q z; t+ y- ?8 d; X5 X
" Y g9 b5 F$ R) L数字PID控制器3 ~5 W3 H9 \8 K$ z, E
" M$ O5 K2 m# `$ n l
(1)模拟PID控制规律的离散化
+ x" b! N0 _* V5 n2 ^7 H& x0 y. d# @ Q+ o3 `
0 L' U! B) C9 o y4 E
0 x0 g% ?7 F. U/ `" ^
- J# D" O; O9 V8 p* w/ B
(2)数字PID控制器的差分方程
7 X' y' d2 z1 @
0 y8 N8 g" q/ |# s! e
q3 k$ L; J; s, X4 P
6 J v* W+ o% k/ D& n
6 \1 o) e" V7 s4 R& [- c& F参数的自整定
3 P" i/ e. v; t- ~. R7 r( z/ C" H
, [9 C* G( r2 m" j( X在某些应用场合,比如通用仪表行业,系统的工作对象是不确定的,不同的对象就得采用不同的参数值,没法为用户设定参数,就引入参数自整定的概念。实质就是在首次使用时,通过N次测量为新的工作对象寻找一套参数,并记忆下来作为以后工作的依据。具体的整定方法有三种:临界比例度法、衰减曲线法、经验法。
/ P. [, ~5 r% ]3 v _$ _0 r( w6 a5 E2 z6 y
1、临界比例度法(Ziegler-Nichols)7 B7 i4 w4 [, s8 @! H- _. x
( d2 B9 l' Z! v1.1 在纯比例作用下,逐渐增加增益至产生等副震荡,根据临界增益和临界周期参数得出PID控制器参数,步骤如下:
& Y/ k% O5 D1 {; }/ y# h
/ \/ P% ?) k0 o9 Z" d% v' r% L(1)将纯比例控制器接入到闭环控制系统中(设置控制器参数积分时间常数Ti =∞,实际微分时间常数Td =0)。
' t/ y% f1 p0 o" G2 V/ P
0 g; [# X5 o# O- [# T' b(2)控制器比例增益K设置为最小,加入阶跃扰动(一般是改变控制器的给定值),观察被调量的阶跃响应曲线。
5 J3 T. v; r( D
1 k$ ^9 t" x2 v7 `(3)由小到大改变比例增益K,直到闭环系统出现振荡。" S" q; P6 e1 b* A$ F4 Y& l
0 Q. e( m; o: Z
(4)系统出现持续等幅振荡时,此时的增益为临界增益(Ku),振荡周期(波峰间的时间)为临界周期(Tu)。) l Y+ p( {" J8 ]$ t" W
1 u5 z6 K8 s0 x9 B
(5) 由表1得出PID控制器参数。
: a9 w N& l( S9 }, H, T* {( O
; A c% P' y2 `; R) o, t表1
0 A" c4 v R' S7 j) l7 I
9 g* p; j8 F. J4 x1 ^. Y9 a( y7 x2 L
* N& r h; G& v) P$ ]9 C% y" r, t
, s. U9 V" Y, w5 C3 z4 T. m8 }5 w
1.2 采用临界比例度法整定时应注意以下几点:& b9 @ T; z3 o8 R; v- H
, Z* e/ L! ^% ~(1)在采用这种方法获取等幅振荡曲线时,应使控制系统工作在线性区,不要使控制阀出现开、关的极端状态,否则得到的持续振荡曲线可能是“极限循环”,从线性系统概念上说系统早已处于发散振荡了。5 k) @$ Q8 ?) [ K! }
, {! ^ K; P7 n; u(2)由于被控对象特性的不同,按上表求得的控制器参数不一定都能获得满意的结果。对于无自平衡特性的对象,用临界比例度法求得的控制器参数往住使系统响应的衰减率偏大(ψ>0.75 )。而对于有自平衡特性的高阶等容对象,用此法整定控制器参数时系统响应衰减率大多偏小(ψ<0.75 )。为此,上述求得的控制器参数,应针对具体系统在实际运行过程中进行在线校正。
3 `, u: C& [# F% K5 l
) C7 j T0 [3 q0 Y0 X(3) 临界比例度法适用于临界振幅不大、振荡周期较长的过程控制系统,但有些系统从安全性考虑不允许进行稳定边界试验,如锅炉汽包水位控制系统。还有某些时间常数较大的单容对象,用纯比例控制时系统始终是稳定的,对于这些系统也是无法用临界比例度法来进行参数整定的。
0 X' C- W) Y3 X9 ^3 Q& s7 m* k% n2 U; B; H( e
(4)只适用于二阶以上的高阶对象,或一阶加纯滞后的对象,否则,在纯比例控制情况下,系统不会出现等幅振荡。
. E$ X$ [& H3 q. }, I1 t! ?8 y- n7 ?4 w8 z, o6 P( X
1.3 若求出被控对象的静态放大倍数KP=△y/△u ,则增益乘积KpKu可视为系统的最大开环增益。通常认为Ziegler-Nichols闭环试验整定法的适用范围为:2<KpKu<20+ y4 h3 Y) @4 b; y8 P+ v2 U9 N
! X; j# t! O* K x; s+ r(1) 当KpKu > 20时,应采用更为复杂的控制算法,以求较好的调节效果。 n9 x0 ]6 I r7 ?" Y+ W- T4 p
, r9 T9 t1 J% d Z# a8 |(2)当KpKu < 2时,应使用一些能补偿传输迟延的控制策略。" c- u/ o, p! g* A5 N
& U I4 L T# _2 _6 |; ]! k6 g(3)当1.5<kpku< 2时,在对控制精度要求不高的场合仍可使用pid控制器,但需要对表1进行修正。在这种情况下,建议采用smith预估控制和imc控制策略。< span>% z: E5 h! b+ X" E+ Z0 p
& B/ ^% y& d2 ^# Z6 U
(4)当KpKu< 1.5时,在对控制精度要求不高的场合仍可使用PI控制器,在这种情况下,微分作用已意义不大。
( P) G- O% ~; k: R% A5 i
5 ~4 @4 N" \& n% x( ?/ }2、衰减曲线法
: v8 ` [- j5 U' o4 X7 Y
. ^. O; Y' B* D& a' ^3 M1 R衰减曲线法与临界比例度法不同的是,闭环设定值扰动试验采用衰减振荡(通常为4:1或10:l),然后利用衰减振荡的试验数据,根据经验公式求取控制器的整定参数。整定步骤如下:
) u f+ |$ Q5 {* Y! |# }3 T$ s
' o3 y& I! m& y- M5 W% `(1)在纯比例控制器下,置比例增益K为较小值,并将系统投入运行。
" K, b1 F. ]+ Z; L5 l1 e/ e9 x# w3 @# V& O S7 P! K) ?; V
(2)系统稳定后,作设定值阶跃扰动,观察系统的响应,若系统响应衰减太快,则减小比例增益K;反之,应增大比例增益K。直到系统出现如下图(a)所示的4:1衰减振荡过程,记下此时的比例增益Ks及和振荡周期Ts数值。" o) _& ~4 G; I& J, x6 r$ i* m
( R& s. f% r; I) U9 J
. b7 b( y& ]5 h+ y; ~. I- C. R& |) d
: ~' n$ X; C2 ^* \& _(3)利用Ks和Ts值,按下表给出的经验公式,计算出控制器的参数整定值。
- c8 E, y. ?2 R W. E! v5 ]3 } R* J- f' k* ~
3 t5 u* G" `( x
, c* A7 @0 e. ]) Y( Z. w8 A# W3 g
' m# A1 [( ?& F/ K(4)10:1衰减曲线法类似,只是用Tr带入计算。
% H& B, b" c& A w& g! W$ r& Q- B5 q0 F' j% H" a) T3 q# J
采用衰减曲线法必须注意几点:% U& ?( f( r4 p3 ~0 z
; X9 w6 B+ Y+ M# ]; x) W(1)加给定干扰不能太大,要根据生产操作要求来定,一般在5%左右,也有例外的情况。4 t! u; E. ]* ?) g$ A) Z3 {
' O$ ^/ _0 M, Y; {/ ?; m$ u
(2)必须在工艺参数稳定的情况下才能加给定干扰,否则得不到正确得 整定参数。9 d* v6 @, _8 k. T
4 P4 I7 G: p- R2 L! y+ C3 p8 g
(3)对于反应快的系统,如流量、管道压力和小容量的液位调节等,要得到严格的4:1衰减曲线较困难,一般以被调参数来回波动两次达到稳定,就近似地认为达到4:1衰减过程了。
$ E* U4 X0 z, k7 z8 s( R. B4 C( [! l+ n
7 i* ^" y- V+ Y3 P6 S3 L(4)投运时,先将K放在较小的数值,把Ti减少到整定值,把Td逐步放大到整定值,然后把K拉到整定值(如果在K=整定值的条件下很快地把Td放到整定值,控制器的输出会剧烈变化)。
5 I4 k1 f1 q `
# a6 q2 i( t; ~* x6 m, C# Z' P" ] |
© 声明:本文仅表作者或发布者个人观点,与环保之家[2TECH.CN]无关。其原创性及陈述文字、内容、数据及图片均未经证实,对本文及其全部或部分内容、图片、文字的真实性、完整性、及时性本站不作任何保证或承诺,仅做参考并自行核实。如有侵权,请联系我们处理,在此深表歉意。
|