基于繼電反饋的組合積分系統(tǒng)參數(shù)辨識(shí)與整定

閆子豪,任正云,陳安鋼

(東華大學(xué) 信息科學(xué)與技術(shù)學(xué)院，上海 201620)

基于繼電反饋的組合積分系統(tǒng)參數(shù)辨識(shí)與整定

閆子豪,任正云,陳安鋼

(東華大學(xué) 信息科學(xué)與技術(shù)學(xué)院，上海 201620)

將一種改進(jìn)型的繼電辨識(shí)方法推廣到組合積分系統(tǒng)的參數(shù)辨識(shí)中,取代了以一階或二階加純滯后模型近似描述組合積分系統(tǒng)進(jìn)行參數(shù)估計(jì)的方法，在得到估計(jì)參數(shù)的基礎(chǔ)上設(shè)計(jì)了組合積分控制器。利用相角偏移的概念，補(bǔ)償繼電試驗(yàn)中振蕩點(diǎn)與被控對(duì)象的極限環(huán)理論臨界點(diǎn)之間的偏差，辨識(shí)得到頻域范圍內(nèi)的組合積分對(duì)象的參數(shù)，并且辨識(shí)精度較高；通過(guò)相角偏移的在線計(jì)算，避免了多次迭代試驗(yàn)獲取先驗(yàn)信息的過(guò)程，縮短了調(diào)節(jié)時(shí)間。仿真結(jié)果驗(yàn)證了該算法的實(shí)用性和有效性。

組合積分 參數(shù)辨識(shí) 參數(shù)整定 相角偏移

工業(yè)過(guò)程中，大多數(shù)過(guò)程對(duì)象通常被近似為一階加純滯后(FOPDT)或二階加純滯后(SOPDT)模型。組合積分系統(tǒng)[1-2]作為近年來(lái)被提出的一種新型的流程工業(yè)系統(tǒng)，在過(guò)去很長(zhǎng)的一段時(shí)間內(nèi)被近似為上述兩種過(guò)程對(duì)象之一。雖然，這種近似處理使得PID[3]和Smith預(yù)估器[4]以及預(yù)測(cè)PI[5]等算法能在該種對(duì)象的控制上取得一定的效果，但是模型近似導(dǎo)致的控制精度、魯棒性無(wú)法進(jìn)一步提高，限制了資源的高效利用。同時(shí)，工業(yè)設(shè)備結(jié)構(gòu)老化和作業(yè)環(huán)境改變會(huì)引起系統(tǒng)參數(shù)漂移，從而使得控制效果惡化。為了解決控制的局限性，應(yīng)對(duì)參數(shù)變化，改善控制效果，針對(duì)組合積分系統(tǒng)的參數(shù)辨識(shí)與自整定控制器的研究是十分有必要的。

20世紀(jì)90年代Astrom等[6]提出利用繼電反饋的方法提取系統(tǒng)的臨界信息，進(jìn)而設(shè)計(jì)了PID自整定控制器。Luyben在此基礎(chǔ)之上，提出了ATV(auto tuning variation)方法[7]，該方法將繼電試驗(yàn)得到的臨界信息與理論上的臨界信息近似，推導(dǎo)得到模型參數(shù)并設(shè)計(jì)控制器參數(shù)。由于簡(jiǎn)單直接，該方法也成為了工業(yè)應(yīng)用中標(biāo)準(zhǔn)的繼電辨識(shí)方法。但是，試驗(yàn)得到的臨界信息與理論上的臨界信息之間的誤差積累導(dǎo)致了辨識(shí)精度始終不高。因此，許多學(xué)者提出了改進(jìn)的方法來(lái)減少這種近似誤差，但改進(jìn)方法大都以增加實(shí)驗(yàn)迭代次數(shù)或者辨識(shí)結(jié)構(gòu)復(fù)雜度來(lái)提高辨識(shí)精度。

2006年，Wang等[8]在標(biāo)準(zhǔn)繼電辨識(shí)方法的基礎(chǔ)上，提出了“相角偏移”的概念，找到了繼電試驗(yàn)中實(shí)際振蕩信息與臨界信息的關(guān)系，消除了臨界值近似誤差，避免了多次迭代實(shí)驗(yàn)，縮短了調(diào)節(jié)時(shí)間，提高了辨識(shí)精度。對(duì)于一般時(shí)滯系統(tǒng)而言，以一階或二階加純滯后模型描述系統(tǒng)，并利用該方法進(jìn)行參數(shù)辨識(shí)，可得到非常準(zhǔn)確的估計(jì)參數(shù)。但對(duì)于組合積分系統(tǒng)而言，利用該方法始終無(wú)法求得以一階或二階加純滯后模型描述系統(tǒng)的估計(jì)參數(shù)。

因此，筆者提出了基于繼電反饋的組合積分系統(tǒng)的參數(shù)辨識(shí)，運(yùn)用“相角偏移”的概念，實(shí)現(xiàn)了組合積分系統(tǒng)多個(gè)參數(shù)的單次辨識(shí)，并且辨識(shí)精度較高。同時(shí)，設(shè)計(jì)了控制器并整定了控制器參數(shù)。

1 繼電反饋的基本原理

繼電器反饋基本原理如圖1所示，圖1中過(guò)程對(duì)象即被控對(duì)象的輸入信號(hào)為u(t)和輸出信號(hào)為y(t)，系統(tǒng)運(yùn)行分為辨識(shí)和控制兩個(gè)部分。辨識(shí)過(guò)程中，繼電器的輸出作為過(guò)程對(duì)象輸入信號(hào)作用于過(guò)程對(duì)象，繼電器接收反饋信號(hào)，實(shí)現(xiàn)開關(guān)切換。

圖1 繼電反饋基本原理示意

辨識(shí)過(guò)程中，系統(tǒng)的輸出在工作點(diǎn)附近呈現(xiàn)頻率為ωosc的周期振蕩，即對(duì)應(yīng)描述函數(shù)中的振蕩點(diǎn)。傳統(tǒng)繼電辨識(shí)根據(jù)描述函數(shù)近似的方法，將該點(diǎn)近似為系統(tǒng)臨界點(diǎn)，即ωosc≈ωu，ωu為臨界振蕩頻率。理論情況下，系統(tǒng)在臨界點(diǎn)會(huì)產(chǎn)生頻率為ωu的周期振蕩，當(dāng)頻率ω=ωu時(shí)，∠G(jω)=-π,振蕩頻率ωosc與臨界頻率ωu存在的偏差以相角關(guān)系表示為∠G(jnωosc)=-π-δpn,n=1, 3, 5…。系統(tǒng)的描述函數(shù)的分析如圖2所示，系統(tǒng)振蕩點(diǎn)與臨界點(diǎn)的相角偏移為δpn。圖2中Im為虛軸，Re為實(shí)軸。

圖2 描述函數(shù)分析示意

2 相角偏移

根據(jù)文獻(xiàn)[3]中的定理，設(shè)繼電反饋系統(tǒng)滿足極限環(huán)條件，令

(1)

(2)

式中：Tosc——振蕩周期；n——諧波階數(shù)。

兩者均為非零常量，則系統(tǒng)開環(huán)Nyquist曲線上振蕩點(diǎn)和臨界點(diǎn)存在相角偏移，且

(3)

式中：n=(2k+1),k=0, 1, 2, 3…。

3 組合積分系統(tǒng)

2) 時(shí)滯τ1i和τ2i必須滿足方程τ2i=τ2(i-1)-τ1(i-1), 1

這類對(duì)象同時(shí)含有非右半平面的極點(diǎn)和零點(diǎn)。從對(duì)象的時(shí)域特性可知，其同時(shí)具有積分特性和非最小相位的特性。

不失一般性，有如下1～5種典型的組合積分對(duì)象：

(Ⅰ)

(Ⅱ)

(Ⅲ)

(Ⅳ)

(Ⅴ)

這些傳遞函數(shù)是根據(jù)測(cè)試得到的輸入、輸出數(shù)據(jù)建模而得到的；或者根據(jù)過(guò)程對(duì)象的機(jī)理特性分析得到的。

4 理論振蕩周期

對(duì)于第二種組合積分對(duì)象(Ⅱ)，假設(shè)系統(tǒng)在理想繼電器的作用下建立了周期振蕩，那么根據(jù)文獻(xiàn)[2]的定理可知：

(4)

式中：Z——變換符號(hào)；z——z域算子。

理想繼電器作用下即d=1且ε=0，有H(T/2, -1)=0，可得理論振蕩周期，即：

T0=τ1+2τ2

(5)

5 組合積分系統(tǒng)的參數(shù)辨識(shí)

組合積分系統(tǒng)的特性不同于一階或二階加純滯后系統(tǒng)特性，在滿足建立極限環(huán)的條件下，振蕩波形輸出表明了兩者的特性差異，如圖3所示。第二種組合積分系統(tǒng)呈現(xiàn)典型的組合積分特性，其振蕩輸出為梯形波。解釋了對(duì)組合積分系統(tǒng)，以一階或二階加滯后模型近似描述系統(tǒng)時(shí)，利用相角偏移的辨識(shí)方法難以解得估計(jì)參數(shù)的情況。

圖3 組合積分系統(tǒng)振蕩輸出波形示意

在系統(tǒng)建立極限環(huán)的條件下，可通過(guò)分析振蕩波形得到振蕩周期Tosc和幅值A(chǔ)m并利用數(shù)值計(jì)算得到相角偏移δpn。最后，根據(jù)建立的頻域信息關(guān)系，推導(dǎo)系統(tǒng)參數(shù)。

下面主要介紹第四種組合積分對(duì)象(Ⅳ)辨識(shí)過(guò)程，其頻率特性為

(6)

通過(guò)相角偏移的補(bǔ)償，建立新的系統(tǒng)頻域關(guān)系，可表示如下3個(gè)式子：

(7)

(8)

(9)

由式(7)～式(9)幅值關(guān)系，可得：

(10)

(11)

(12)

由相角關(guān)系得到：

(13)

式(7)～式(9)中的kn,n=1,3,5可由估計(jì)對(duì)象頻率響應(yīng)利用數(shù)值計(jì)算得到，即

(14)

式中：h——理想繼電器輸出電平。

聯(lián)立式(10)～(14)，由于上述方程均為超越方程沒有解析解，利用1STOPT最優(yōu)化工具箱進(jìn)行計(jì)算，可得到各個(gè)參數(shù)即k,τ1,τ2,T。對(duì)于第二種組合積分對(duì)象，推導(dǎo)方式同上，可增加1組約束條件即式(5)，使得方程組不必分次求解，單次求解即可得到符合要求的解。

6 仿真結(jié)果與分析

為了驗(yàn)證算法的有效性，選幾組不同的組合積分對(duì)象模型試驗(yàn)比較。

假設(shè)ATV先驗(yàn)信息已知，對(duì)第二種組合積分對(duì)象進(jìn)行參數(shù)辨識(shí)結(jié)果，見表1所列，由圖4的Nyquist曲線對(duì)比可知： 利用相角偏移方法辨識(shí)得到系統(tǒng)模型幾乎與實(shí)際系統(tǒng)重合，ATV方法辨識(shí)結(jié)果存在一定偏差。

對(duì)第四種組合積分對(duì)象，直接利用Pade近似模型與實(shí)際對(duì)象對(duì)比，參數(shù)辨識(shí)結(jié)果見表2所列。由圖5的Nyquist曲線對(duì)比可知： 相角偏移方法辨識(shí)得到的模型幾乎與實(shí)際系統(tǒng)重合，效果優(yōu)于以Pade近似模型描述第四種組合積分對(duì)象。

表1 第二種組合積分對(duì)象的辨識(shí)結(jié)果比較

圖4 估計(jì)模型的Nyquist曲線比較示意

Tτ1τ2實(shí)際對(duì)象Pade近似模型相角偏移方法23213s(2s+1)(1-e-3s)e-2s12s2+3．5s+1e-2s0．98493．001s(1．9497s+1)(1-e-3．001s)e-1．9877s34514s(3s+1)(1-e-4s)e-5s16s2+5s+1e-5s1．00173．9986s(3．0111s+1)(1-e-3．9986s)e-4．7201s

圖5 Pade近似、估計(jì)模型Nyquist曲線比較

7 組合積分控制器的設(shè)計(jì)

考慮第四種組合積分對(duì)象，由于組合積分對(duì)象具有良好的開環(huán)特性，當(dāng)給定1個(gè)階躍時(shí)，系統(tǒng)能夠快速響應(yīng)達(dá)到給定值。二階組合積分系統(tǒng)的階躍響應(yīng)曲線如圖6所示。

圖6 二階組合積分系統(tǒng)的階躍響應(yīng)曲線示意

為了使得整個(gè)閉環(huán)系統(tǒng)具有如圖6所示的良好特性，選擇所期望的閉環(huán)傳遞函數(shù)具有以下的結(jié)構(gòu)形式：

(15)

式中：λ,τ10,τ20——整定參數(shù)。當(dāng)λ=1，開環(huán)的響應(yīng)時(shí)間與閉環(huán)響應(yīng)時(shí)間相同；當(dāng)λ>1，開環(huán)的響應(yīng)時(shí)間比閉環(huán)響應(yīng)時(shí)間快；當(dāng)λ<1，開環(huán)的響應(yīng)時(shí)間比閉環(huán)的響應(yīng)時(shí)間慢。這樣可以推導(dǎo)控制器傳遞函數(shù)：

(16)

假定λ=1,τ10=τ1,τ20=τ2,k0=k,有

(17)

控制器在時(shí)間域的輸入、輸出關(guān)系為

(18)

式(18)第一部分為比例項(xiàng)，第二部分可以解釋為控制器在t時(shí)刻的輸出是由控制器在過(guò)去時(shí)間[t-(τ1+τ2),t-τ2]的輸出預(yù)測(cè)而得到。

考慮1個(gè)二階組合積分對(duì)象，利用IMC-PID[9]，預(yù)測(cè)PI[5]，PID[10]，組合積分控制器進(jìn)行仿真比較：

(19)

根據(jù)辨識(shí)得到的參數(shù)，可知系統(tǒng)的傳遞函數(shù)，即：

(20)

以上控制算法得到含有階躍干擾情況下的階躍響應(yīng)曲線比較，如圖7所示。

圖7 標(biāo)稱系統(tǒng)響應(yīng)曲線示意

由圖7可見，組合積分控制器快速性優(yōu)良，幾乎沒有超調(diào)，而IMC-PID與PID都或多或少存在一定的超調(diào)，并且調(diào)節(jié)時(shí)間較長(zhǎng)。預(yù)測(cè)PI雖然無(wú)超調(diào)，但速度稍慢。總體而言，組合積分控制器的控制效果快速平穩(wěn)。

8 結(jié)束語(yǔ)

本文將相角偏移的概念運(yùn)用到組合積分系統(tǒng)的繼電反饋辨識(shí)中，避免了以一階或二階加純滯后模型描述組合積分系統(tǒng)進(jìn)行參數(shù)辨識(shí)無(wú)法解得估計(jì)參數(shù)的問題。利用相角偏移補(bǔ)償繼電實(shí)驗(yàn)振蕩

點(diǎn)與理論臨界值之間的偏差，從而消除了近似誤差。在不需要先驗(yàn)信息的情況下，可辨識(shí)得到組合積分系統(tǒng)的多個(gè)參數(shù)，并在此基礎(chǔ)之上設(shè)計(jì)得到了組合積分控制器，且控制效果快速平穩(wěn)。

[1] 任正云,鄭達(dá).流程工業(yè)的組合積分系統(tǒng)及其先進(jìn)控制[J].控制理論與應(yīng)用,2009,26(09)： 1009-1013.

[2] 鄭達(dá).組合積分系統(tǒng)的控制理論與應(yīng)用研究[D].上海： 東華大學(xué),2009.

[3] 邵惠鶴.工業(yè)過(guò)程高級(jí)控制[M].2版.上海： 上海交通大學(xué)出版社,2003.

[4] Smith O J M. A Controller to Overcome Dead Time[J]. Isa Transactions, 1959, 6(02)： 28-33.

[5] 任正云,邵惠鶴,張立群.幾種特殊動(dòng)態(tài)特性對(duì)象的預(yù)測(cè)PI控制[J].儀器儀表學(xué)報(bào),2004,25(05)： 615-619.

[6] ?str?m K J, H?gglund T. Automatic Tuning of Simple Regulators With Specifications on Phase and Amplitude Margins [J]. Automatica, 1984, 20(05)： 645-651.

[7] Luyben W L. Derivation of Transfer Functions for Highly Nonlinear Distillation Columns[J]. Ind. Eng. Chem. Res, 1987, 26(12)： 2490-2495.

[8] 王萍,顧誕英,張衛(wèi)東.一種新的繼電自整定控制器設(shè)計(jì)方法[J].上海交通大學(xué)學(xué)報(bào),2007,41(08)： 1314-1319

[9] 孫功武,聶紅偉,蘇義鑫,等.時(shí)滯系統(tǒng)的二自由度IMC-PID控制研究[J].計(jì)算機(jī)應(yīng)用研究,2014,31(08)： 2357-2360.

IdentificationandParameterTuningofCombinedIntegratingSystemBasedonRelayFeedback

Yan Zihao, Ren Zhengyun, Chen Angang

(College of Information Science and Technology, Donghua University, Shanghai,201620, China)

An improved relay identification method is extended to parameter identification of combined integrating system. The parameter estimation method approximately describing combined integrating system with the first order or the second order plus dead-time model is replaced. A combined integrating controller is designed based on obtaining estimated parameters. Using the concept of phase angle offset, the deviation between the oscillation point in the relay test and the critical value of the limit cycles of the controlled object is compensated,as well as the parameters of the combined integrating object in the frequency domain are identified and the identification accuracy is higher than before. The process of obtaining a priori information by multiple iterations is avoided through online calculation of phase angle offset.The adjustment time is shortened. The simulation results verify the practicability and validity of the algorithm.

combined integrating; parameter identification; parameter tuning; phase angle offset

TP273

B

1007-7324(2017)05-0033-06

稿件收到日期： 2017-06-29，修改稿收到日期2017-07-20。

閆子豪(1992—)，男，湖北宜昌人，在讀碩士研究生，主要從事先進(jìn)過(guò)程控制方向的研究。