基于頻域建模與遺傳算法的電力電子電路參數(shù)辨識(shí)方法

孫鳳艷 王友仁 林 華 崔 江 姜媛媛

（南京航空航天大學(xué)自動(dòng)化學(xué)院 南京 210016）

1 引言

隨著航空與航天事業(yè)的不斷發(fā)展，航空航天飛行器的安全性和可靠性要求日益提高，由此飛機(jī)健康預(yù)報(bào)與管理技術(shù)得到越來(lái)越多的關(guān)注[1-2]。電力電子電路是飛機(jī)二次電源系統(tǒng)的重要組成部分[3]，對(duì)電力電子電路故障診斷和故障預(yù)測(cè)技術(shù)的研究具有重要意義與應(yīng)用前景。

電力電子電路故障分為結(jié)構(gòu)型故障和參數(shù)型故障。目前，電力電子電路故障診斷技術(shù)的研究多數(shù)采用智能診斷方法，如支持向量機(jī)[4]、模糊推理[5]和神經(jīng)網(wǎng)絡(luò)等[6-9]，這類(lèi)方法需要預(yù)先劃分元器件參數(shù)偏差等級(jí)來(lái)設(shè)置故障模式，利用不同故障模式下的樣本數(shù)據(jù)來(lái)訓(xùn)練分類(lèi)器，由此實(shí)現(xiàn)故障診斷。當(dāng)元器件參數(shù)偏差較小時(shí)，由于元器件容差的存在，會(huì)導(dǎo)致故障診斷正確率很低，甚至難以進(jìn)行故障診斷。為了有效診斷電力電子電路參數(shù)型故障，出現(xiàn)了基于參數(shù)辨識(shí)的電力電子電路故障診斷方法[10-12]。通過(guò)參數(shù)辨識(shí)能準(zhǔn)確獲知電路中故障元器件參數(shù)值，從而實(shí)現(xiàn)參數(shù)型故障診斷，這種方法一般需要建立電路數(shù)學(xué)模型，且參數(shù)估計(jì)方法直接決定了參數(shù)辨識(shí)的準(zhǔn)確性。

面向電力電子電路參數(shù)型故障診斷，研究一種電力電子電路參數(shù)辨識(shí)新方法。本文建立了電路頻域傳遞函數(shù)模型，給出了基于傅里葉變換的頻域特性求解方法，利用遺傳算法進(jìn)行電力電子電路參數(shù)辨識(shí)，以Buck電路為例，驗(yàn)證了該方法的有效性。

2 電力電子電路參數(shù)辨識(shí)方法

2.1 頻域傳遞函數(shù)模型

針對(duì)電力電子電路參數(shù)辨識(shí)，必須建立其電路數(shù)學(xué)模型。對(duì)于不同的電路，應(yīng)根據(jù)電路的功能和結(jié)構(gòu)特點(diǎn)，確定反映電路工作過(guò)程的頻域模型。

Buck電路是電力電子電路中實(shí)現(xiàn)降壓轉(zhuǎn)換的電路，是電力電子電路中最基本的電路。以圖1所示的 Buck電路為例，分析電力電子電路頻域傳遞函數(shù)模型的建立過(guò)程。

圖1 Buck電路圖Fig.1 Circuit diagram of Buck converter

圖1中，開(kāi)關(guān)MOSFET為IRF151，其開(kāi)關(guān)頻率 f =50kHz，占空比D=0.22；濾波電感L為43μH，VD1為續(xù)流二極管；輸出濾波電容C為220μF，此處采用理想電容C和等效串聯(lián)電阻RESR來(lái)表示實(shí)際電容；負(fù)載性質(zhì)為阻性負(fù)載1.25Ω。

通過(guò)電路分析可知，開(kāi)關(guān)管導(dǎo)通、截止及二極管的截止與導(dǎo)通對(duì)于電感、電容器及負(fù)載支路來(lái)講，電路的拓?fù)浣Y(jié)構(gòu)并沒(méi)有發(fā)生改變。Buck電路工作在連續(xù)工作方式，當(dāng)開(kāi)關(guān)管導(dǎo)通、二極管截止時(shí)，電感輸入端電壓值為直流輸入電壓與開(kāi)關(guān)導(dǎo)通壓降之差；當(dāng)開(kāi)關(guān)管截止、二極管導(dǎo)通時(shí)，電感輸入電壓值為二極管正向?qū)▔航?。在開(kāi)關(guān)管和二極管為理想器件情況下，該電感輸入端電壓信號(hào)可以等效為一個(gè)脈沖波形uin。其高電平等于輸入直流電壓、低電平為零。此時(shí) Buck電路的拓?fù)浣Y(jié)構(gòu)等效電路如圖2所示。

圖2 MOSFET開(kāi)關(guān)之后等效電路圖Fig.2 Equivalent circuit of Buck converter

根據(jù)頻域傳遞函數(shù)定義，圖2所示電路的傳遞函數(shù)H(s)為

式中 Uin(s)—電感輸入端脈沖電壓 uin的拉普拉斯變換；

Uo(s)—輸出電壓uo的拉普拉斯變換。

將s=jω代入式（1），可得

式中 Uin(jω)—電感輸入端脈沖電壓 uin的傅里葉變換；

Uo(jω)—輸出電壓uo的傅里葉變換；

H(jω)—電路頻域特性。

2.2 頻域特性求解

對(duì)于電力電子電路頻域特性或頻域模型參數(shù)的求解，可通過(guò)監(jiān)測(cè)有關(guān)電路狀態(tài)信號(hào)時(shí)域波形，再對(duì)其進(jìn)行傅里葉變換獲得頻域傳遞函數(shù)值。

設(shè)x(nTs)是狀態(tài)信號(hào)x(t)的采樣，Ts為采樣間隔，采樣數(shù)據(jù)長(zhǎng)度為 N。則離散傅里葉變換（Discrete Fourier Transform，DFT）算法如下：

利用基2時(shí)析型快速傅里葉變換（Fast Fourier Transform，F(xiàn)FT）來(lái)實(shí)現(xiàn) DFT。設(shè) N=2M，M 為正整數(shù)，如果序列長(zhǎng)度不滿足2的整數(shù)次冪可在序列后補(bǔ)零?；?時(shí)析型FFT將輸入時(shí)間序列x(n)按奇偶對(duì)分的原則分割序列得到奇、偶數(shù)號(hào)序列，即令n=2r和n=2r+1，r=0，1，…，N/2-1，于是

令

則輸入序列x(n)的離散傅里葉變換為

2.3 電路參數(shù)辨識(shí)

電力電子電路參數(shù)辨識(shí)的基本思路為：首先，根據(jù)電路原理建立電路頻域模型；其次，獲取電路狀態(tài)監(jiān)測(cè)信號(hào)，利用傅里葉分析得到電路頻率響應(yīng)特性；基于電路頻域模型及頻響特性，利用參數(shù)辨識(shí)算法進(jìn)行模型參數(shù)估計(jì)。考慮到模型參數(shù)估計(jì)過(guò)程中，涉及多變量高次方程求解十分困難，為此，把電路參數(shù)估計(jì)問(wèn)題轉(zhuǎn)換成優(yōu)化問(wèn)題求解，利用遺傳算法來(lái)實(shí)現(xiàn)電路參數(shù)辨識(shí)。

遺傳算法（Genetic Algorithm，GA）是一類(lèi)借鑒生物界的進(jìn)化規(guī)律演化而來(lái)的隨機(jī)搜索方法，具有簡(jiǎn)單通用、魯棒性強(qiáng)，且不易陷入局部最優(yōu)等優(yōu)點(diǎn)，已被廣泛地應(yīng)用于組合優(yōu)化、機(jī)器學(xué)習(xí)、信號(hào)處理、自適應(yīng)控制及曲線擬合等領(lǐng)域。利用遺傳算法對(duì)電力電子電路進(jìn)行參數(shù)辨識(shí)，電路參數(shù)辨識(shí)算法實(shí)現(xiàn)流程如圖3所示。算法實(shí)現(xiàn)主要過(guò)程如下：

圖3 電路參數(shù)辨識(shí)算法流程圖Fig.3 Flow chart of circuit parameter identification method

（1）編碼方案：根據(jù)電路中待求解元器件參數(shù)的變化范圍，采用二進(jìn)制編碼及解碼方案。

（2）適應(yīng)度函數(shù)與適應(yīng)度評(píng)估：把電路參數(shù)估計(jì)轉(zhuǎn)換成目標(biāo)優(yōu)化問(wèn)題。選擇電路頻域模型作為適應(yīng)度函數(shù)，而實(shí)際電路頻率響應(yīng)特性作為理想適應(yīng)度值（即對(duì)應(yīng)的優(yōu)化目標(biāo)最佳值）。通過(guò)計(jì)算個(gè)體適應(yīng)度值與最優(yōu)適應(yīng)度值之間的誤差大小來(lái)判斷個(gè)體的適應(yīng)度優(yōu)劣，利用遺傳算法進(jìn)化搜索得到最優(yōu)的個(gè)體。最優(yōu)個(gè)體解碼就對(duì)應(yīng)電路元器件參數(shù)，從而實(shí)現(xiàn)電力電子電路元器件參數(shù)辨識(shí)。

（3）遺傳操作：交叉操作是把兩個(gè)隨機(jī)選擇配對(duì)的染色體個(gè)體，按單點(diǎn)交叉方式相互交換其部分基因（交叉點(diǎn)按照概率選擇），以形成兩個(gè)新個(gè)體；變異操作是對(duì)染色體基因位進(jìn)行概率變異，每個(gè)染色體基因位的變異概率相同。

（4）種群生成：初始種群是隨機(jī)生成的，種群規(guī)模預(yù)先確定。在進(jìn)化過(guò)程中，種群規(guī)模保持不變；采用最優(yōu)個(gè)體保存策略，當(dāng)前種群中適應(yīng)度最高的個(gè)體不參與交叉和變異運(yùn)算，而是直接復(fù)制到下一代；利用“輪盤(pán)賭比例選擇法”來(lái)選擇染色體，某個(gè)體被選擇的概率取決于該個(gè)體的適應(yīng)度值與全部個(gè)體適應(yīng)度值之和的比值。隨機(jī)選擇兩個(gè)個(gè)體進(jìn)行遺傳操作，以生成兩個(gè)新個(gè)體。某個(gè)體的選擇概率大，則能被多次選中，其遺傳基因就會(huì)在種群中擴(kuò)大，而某個(gè)體的選擇概率小，則可能被淘汰。

（5）進(jìn)化終止條件：當(dāng)達(dá)到預(yù)定的進(jìn)化代數(shù)，或者最優(yōu)個(gè)體的適應(yīng)度值與理想適應(yīng)度值之差為某一閾值時(shí)，則進(jìn)化過(guò)程終止。

（6）最優(yōu)染色體解碼，獲得電路參數(shù)估計(jì)值。

3 實(shí)例分析

以圖 1所示 Buck電路為例，采用數(shù)字仿真與實(shí)驗(yàn)方案來(lái)驗(yàn)證評(píng)估電路參數(shù)辨識(shí)方法的可行性與有效性。

3.1 電路頻域模型響應(yīng)特性結(jié)果分析

利用FFT方法獲取Buck電路中電感輸入端電壓uin和電路輸出電壓uo的頻域特性。電路正常時(shí)，uin、uo的頻域響應(yīng)特性分別如圖4與圖5所示。由圖可知，頻譜分量主要集中在開(kāi)關(guān)頻率f0=50kHz及f0的整數(shù)倍處，且在f0時(shí)最大。

圖4 電感輸入端電壓FFT頻譜Fig.4 FFT spectrum of the inductor’s input voltage

圖5 輸出電壓FFT頻譜Fig.5 FFT spectrum of the output voltage

3.2 電路參數(shù)辨識(shí)算法實(shí)現(xiàn)及結(jié)果分析

3.2.1 Buck電路參數(shù)辨識(shí)遺傳算法編程實(shí)現(xiàn)

（1）編碼。選用二進(jìn)制編碼方式對(duì)待估計(jì)參數(shù)進(jìn)行編碼。電路模型有四個(gè)未知參數(shù)，則設(shè)置每個(gè)染色體個(gè)體有4個(gè)基因段，每個(gè)參數(shù)對(duì)應(yīng)的基因段有8個(gè)基因位，則染色體共有32位基因組成，各位所表示的含義見(jiàn)表1。

表1 染色體基因位含義表Tab.1 The present meaning of chromosome

設(shè)電容 C、串聯(lián)電阻 RESR、電感 L和電阻 RL的變化范圍分別為 220～130μF、0.45～1.35Ω、43～30μH、1.25～0.5Ω，將各參數(shù)的變化范圍作為其編碼與解碼區(qū)間。

（2）選擇算子。采用“輪盤(pán)賭”選擇方法，每個(gè)個(gè)體進(jìn)入下一代的概率為其適應(yīng)度值與整個(gè)種群中個(gè)體適應(yīng)度之和的比值。

（3）適應(yīng)度函數(shù)及適應(yīng)度評(píng)估。選擇式（2）作為適應(yīng)度函數(shù)，函數(shù)的模為適應(yīng)度值。由式（2）可知，染色體適應(yīng)度值是頻率ω 和電路參數(shù)的函數(shù)。選擇頻率ω 作為電路頻域傳遞函數(shù)的輸入，實(shí)際測(cè)量計(jì)算得到的|H(jω)|作為理想適應(yīng)度值。

適應(yīng)度評(píng)估是在N點(diǎn)頻率下計(jì)算染色體個(gè)體適應(yīng)度值與理想適應(yīng)度值之間的誤差平方和來(lái)進(jìn)行評(píng)估，誤差平方和越小表明染色體越好。選擇每代中誤差平方和最小的個(gè)體作為本代的最優(yōu)個(gè)體，并將本代的最優(yōu)個(gè)體與進(jìn)化過(guò)程中總體最優(yōu)個(gè)體做比較，誤差平方和最小的個(gè)體就是總體最優(yōu)個(gè)體。

（4）交叉操作。在每一代個(gè)體中，若產(chǎn)生的隨機(jī)數(shù)小于設(shè)定的交叉概率Pc=0.6，則隨機(jī)選擇兩個(gè)染色體，且隨機(jī)選擇交叉點(diǎn)位置，在交叉點(diǎn)處交換基因?qū)崿F(xiàn)父代的交叉。

（5）變異操作。變異過(guò)程是通過(guò)設(shè)置初始變異概率Pm=0.1，并隨著進(jìn)化代數(shù)的增加逐漸減小變異概率。

（6）終止條件。以終止代數(shù)、適應(yīng)度誤差作為進(jìn)化過(guò)程結(jié)束條件。設(shè)定終止代數(shù)為5000代，適應(yīng)度誤差為0.135。

3.2.2 參數(shù)辨識(shí)結(jié)果及分析

對(duì)電路正常時(shí)及電解電容器、電感、負(fù)載電阻參數(shù)、所有器件參數(shù)發(fā)生變化五種典型情況分別進(jìn)行電路仿真，監(jiān)測(cè)電感輸入電壓和電路輸出電壓，通過(guò)FFT分析得到電路模型傳遞函數(shù)幅值，從中選擇50～500kHz之間50kHz的整數(shù)倍頻作為理想適應(yīng)度計(jì)算和適應(yīng)度評(píng)估的頻率點(diǎn)，利用遺傳算法進(jìn)行電路參數(shù)辨識(shí)，有關(guān)結(jié)果見(jiàn)表2。

表2 仿真情況下Buck電路參數(shù)辨識(shí)結(jié)果Tab.2 The parameter identification results of Buck converter by simulation method

由表2可知，利用本文方法對(duì)電路元器件參數(shù)進(jìn)行辨識(shí)，無(wú)故障時(shí)相對(duì)誤差最??；當(dāng)負(fù)載電阻為0.8Ω時(shí)，辨識(shí)誤差最大，達(dá)到 14.5%。對(duì)電阻減小到1Ω以下的情況，做了多次實(shí)驗(yàn)，所得結(jié)果的誤差都比較大。而對(duì)于電阻大于1Ω情況，則電阻辨識(shí)結(jié)果的相對(duì)誤差都較??；經(jīng)分析，辨識(shí)誤差主要來(lái)源于FFT分析和遺傳算法參數(shù)辨識(shí)過(guò)程中。不同元器件參數(shù)變化情況下的辨識(shí)結(jié)果相對(duì)誤差都在15%之內(nèi)，可認(rèn)為本文方法能夠?qū)崿F(xiàn)Buck電路參數(shù)辨識(shí)。

3.3 電路參數(shù)辨識(shí)方法實(shí)驗(yàn)驗(yàn)證

為了驗(yàn)證所提出參數(shù)辨識(shí)方法對(duì)實(shí)際電路的參數(shù)辨識(shí)效果，搭建了 Buck實(shí)驗(yàn)電路。實(shí)驗(yàn)中開(kāi)關(guān)頻率為30kHz，電感L的標(biāo)稱值為100μH，電容C的標(biāo)稱值為2200μF，輸出負(fù)載RL為2Ω。為便于測(cè)試數(shù)據(jù)分析，采用PCI—1714UL板卡來(lái)實(shí)時(shí)采集實(shí)驗(yàn)電路監(jiān)測(cè)信號(hào)，A-D分辨率為12位，采樣頻率為600kHz。實(shí)驗(yàn)電路參數(shù)辨識(shí)結(jié)果見(jiàn)表 3，其中各元器件參數(shù)實(shí)際值是利用Agilent 4263B LCR分析儀測(cè)量得到。

表3 Buck實(shí)驗(yàn)電路參數(shù)辨識(shí)結(jié)果Tab.3 The parameter identification results of experimental Buck converter

由表 3可知，本文方法能夠進(jìn)行電路參數(shù)辨識(shí)，但是與仿真數(shù)據(jù)辨識(shí)結(jié)果相比，實(shí)驗(yàn)電路參數(shù)辨識(shí)結(jié)果誤差比較大。分析誤差產(chǎn)生的主要原因有：

（1）噪聲信號(hào)的影響。噪聲是影響參數(shù)辨識(shí)結(jié)果的重要原因。實(shí)驗(yàn) Buck電路是在面包板上連線搭建的，監(jiān)測(cè)信號(hào)容易受到噪聲影響，而電壓信號(hào)中紋波包含許多有用的頻譜分量，因紋波電壓幅值較小受到噪聲的影響也較大。

（2）監(jiān)測(cè)信號(hào)采樣精度的影響。因 A-D轉(zhuǎn)換位數(shù)有限，而輸出電壓中直流分量大導(dǎo)致采樣電壓范圍大，使得輸出電壓紋波信號(hào)部分的采樣分辨率較小，則電路頻域響應(yīng)特性計(jì)算誤差大，引起參數(shù)辨識(shí)結(jié)果誤差大。

（3）FFT頻譜分析和遺傳算法參數(shù)辨識(shí)算法對(duì)電路元器件參數(shù)估計(jì)也有一定的計(jì)算誤差影響。

4 結(jié)論

提出了一種基于電路頻域建模和遺傳算法參數(shù)求解的電力電子電路參數(shù)辨識(shí)方法。首先建立電路頻域模型，再對(duì)電路監(jiān)測(cè)信號(hào)進(jìn)行FFT分析獲得電路頻域響應(yīng)特性。然后，根據(jù)電路頻域特性，通過(guò)把電路參數(shù)辨識(shí)轉(zhuǎn)換成目標(biāo)優(yōu)化問(wèn)題，利用遺傳算法優(yōu)化搜索實(shí)現(xiàn)電路參數(shù)辨識(shí)，避免了多變量高次方程的求解，通用性較好。以 Buck電路為例，通過(guò)仿真與實(shí)驗(yàn)驗(yàn)證了新方法是可行有效的。新的電路參數(shù)辨識(shí)方法能用于電力電子電路參數(shù)型故障診斷，但實(shí)際電路參數(shù)辨識(shí)結(jié)果精度還有待進(jìn)一步提高，后續(xù)工作將深入研究采取措施以提高電路參數(shù)辨識(shí)精度，例如，監(jiān)測(cè)信號(hào)噪聲濾波、紋波電壓精密采樣等。

[1]Lall P, Hande M, Bhat C, et al. Prognostics health monitoring (PHM)for prior-damage assessment in electronics equipment under thermo-mechanical loads[C]. Proceedings of the Electronic Components and Technology Conference. New York, USA, 2007:1097-1111.

[2]Ma Zhangshan. A new life system approach to the prognostic and health management (PHM)with survival analysis, dynamic hybrid fault models,evolutionary game theory, and three-layer survivability analysis[C]. Proceedings of IEEE Aerospace Conference, Big Sky, Montana, USA, 2009: 1-20.

[3]Orsagh R, Brown D, Roemer M, et al. Prognostic health management for avionics system power supplies[C]. Proceedings of the IEEE Aerospace Conference, Big Sky, Montana, USA, 2005: 3585-3591.

[4]崔江, 王友仁, 劉權(quán). 基于高階譜與支持向量機(jī)的電力電子電路故障診斷技術(shù)[J]. 中國(guó)電機(jī)工程學(xué)報(bào),2007, 27(10): 62-66.Cui Jiang, Wang Youren, Liu Quan. The technique of power electronic circuit fault diagnosis based on higher-order spectrum analysis and support vector machines[J]. Proceedings of the CSEE, 2007, 27(10):62-66.

[5]崔江, 王友仁. 采用基于模糊推理的分類(lèi)器融合方法診斷電力電子電路參數(shù)故障[J]. 中國(guó)電機(jī)工程學(xué)報(bào), 2009, 29(18): 54-59.Cui Jiang, Wang Youren. Testing the parametric faults of power electronic circuits applying classifiers fusion method based on fuzzy inference[J].Proceedings of the CSEE, 2009, 29(18): 54-59.

[6]肖嵐, 李睿. 逆變器并聯(lián)系統(tǒng)功率管開(kāi)路故障診斷研究[J]. 中國(guó)電機(jī)工程學(xué)報(bào), 2006, 26(4): 99-104.Xiao Lan, Li Rui. Research on open-circuit fault diagnosis of transistor in inverter paralleling system[J].Proceedings of the CSEE, 2006, 26(4): 99-104.

[7]湯清泉, 顏世超, 盧松升, 等. 三電平逆變器的功率管開(kāi)路故障診斷[J]. 中國(guó)電機(jī)工程學(xué)報(bào), 2008,28(21): 26-32.Tang Qingquan, Yan Shichao, Lu Songsheng, et al.Opencircuit fault diagnosis of transistor in three-level inverter[J]. Proceedings of the CSEE, 2008, 28(21):26-32.

[8]陳如清. 采用新型粒子群算法的電力電子裝置在線故障診斷方法[J]. 中國(guó)電機(jī)工程學(xué)報(bào), 2008, 28(24):70-74.Chen Ruqing. A novel PSO based on-line fault diagnosis approach for power electronic system[J].Proceedings of the CSEE, 2008, 28(24): 70-74.

[9]羅慧, 王友仁, 崔江, 等. 電力電子電路多源特征層融合智能故障診斷方法[J]. 電機(jī)與控制學(xué)報(bào),2010, 14(4): 92-96.Luo Hui, Wang Youren, Cui Jiang, et al. Intelligent fault diagnosis for power electronic circuits based on multi-source feature-level fusion[J]. Electric Machines and Control, 2010, 14(4): 92-96.

[10]Shetty P, Mylaraswamy D, Ekambaram T, et al. A hybrid prognostic model formulation system identification and health estimation of auxiliary power units[C]. Proceedings of the IEEE Aerospace Conference, Big Sky, Montana, USA, 2006: 10.

[11]馬皓, 毛興云, 徐德鴻. 基于混雜系統(tǒng)模型的DC/DC電力電子電路參數(shù)辨識(shí)[J]. 中國(guó)電機(jī)工程學(xué)報(bào), 2005, 25(10): 50-54.Ma Hao, Mao Xingyun, Xu Dehong. Parameter identification of DC-DC power electronic circuit based on hybrid system model[J]. Proceedings of the CSEE, 2005, 25(10): 50-54.

[12]張志學(xué), 馬皓, 毛興云. 基于混雜模型系統(tǒng)理論的電力電子電路故障診斷[J]. 中國(guó)電機(jī)工程學(xué)報(bào),2005, 25(3): 49-53.Zhang Zhixue, Ma Hao, Mao Xingyun. Fault diagnosis for power electronic circuit based on hybrid system theory and event identification[J]. Proceedings of the CSEE, 2005, 25(3): 49-53.