基于自適應(yīng)LMS的電路諧波電流檢測

徐鵬輝, 胥 飛

(上海電機學(xué)院 電氣學(xué)院，上海 201306)

基于自適應(yīng)LMS的電路諧波電流檢測

徐鵬輝, 胥 飛

(上海電機學(xué)院 電氣學(xué)院，上海 201306)

為濾除電路中的諧波，需要精確地檢測出電路中所含的諧波量。研究了自適應(yīng)最小均方誤差(LMS)諧波電流檢測算法。用自適應(yīng)濾波器取代低通濾波器，解決了諧波檢測實時性問題。自適應(yīng)LMS諧波電流檢測算法對電網(wǎng)電壓畸變、頻率波動等電網(wǎng)參數(shù)變化具有良好的自適應(yīng)調(diào)節(jié)能力，可運用于三相或單相電路，改變算法中的收斂因子，加快了諧波的檢測速度。Matlab仿真結(jié)果驗證了新算法的精確性和實時性。

諧波檢測； 自適應(yīng)濾波器； 最小均方誤差； 收斂因子

隨著電力電子技術(shù)廣泛應(yīng)用于工業(yè)領(lǐng)域，電網(wǎng)中的諧波污染問題也越來越嚴(yán)重。諧波的存在會使電氣設(shè)備出現(xiàn)過熱情況，加速了設(shè)備老化，降低了使用壽命，甚至引發(fā)設(shè)備故障或燒毀；另外，諧波也可能引起繼電保護和自動裝置產(chǎn)生誤動作，造成電力系統(tǒng)故障[1]。

為了消除電路中的諧波含量，就要對電路中的諧波量進行精確測量，需要采用一定的測量方法計算出電路中所含的諧波總量。目前，常用的諧波檢測方法是瞬時無功功率檢測法，由 Akagi 等[2]于1983年提出，即p-q理論，定義了瞬時有功功率p、瞬時無功功率q，該方法計算簡單，但誤差較大，諧波含量較高。文獻[3]中提出了ip-iq理論，即定義瞬時有功電流ip、瞬時無功電流iq，間接求出p、q,該方法計算簡單，但只適用于三相電路，不適合于單相電路，且無法檢測出諧波中的任意次諧波。文獻[4]中研究的d-q諧波和無功功率檢測法，諧波檢測量較高，但在三相不對稱電路中，因檢測電路復(fù)雜，降低了檢測精度。文獻[5]中研究的d-q-0諧波電流檢測法運算速度快，但是對電流、電壓的同步采樣要求較高，無法滿足實時性要求。文獻[6-7]中提出的自適應(yīng)諧波檢測法，檢測精度高，但是結(jié)構(gòu)中含有的濾波器器件較多，結(jié)構(gòu)復(fù)雜，計算繁瑣。文獻[8-9]中提出的基于神經(jīng)元的自適應(yīng)算法，控制和計算精度較高，但是計算過程復(fù)雜，運算速度較慢，不適用于實時變化的諧波電路。

本文提出的自適應(yīng)最小均方誤差(Least Mean Square, LMS)諧波電流檢測法，是在自適應(yīng)噪聲對消技術(shù)[10]及自適應(yīng)濾波器[11]基礎(chǔ)上發(fā)展起來的電流檢測算法，主要針對不確定對象，尤其是時變、非線性和隨機系統(tǒng),不僅可以精確檢測三相電路、單相電路中的諧波量，而且可以檢測非對稱電路中的諧波量。該方法自適應(yīng)能力強，所需模塊少，計算簡單，動態(tài)性高。

1 自適應(yīng)LMS濾波器

自適應(yīng)LMS諧波檢測法是在自適應(yīng)噪聲對消技術(shù)的基礎(chǔ)上，由含有1個90°移相因子的二階陷波[12](交流輸入電流從一相切換到另一相時產(chǎn)生的周期性電壓擾動)濾波器的閉環(huán)連續(xù)系統(tǒng)組成。

圖1 自適應(yīng)LMS濾波器結(jié)構(gòu)圖

假設(shè)電源電壓為標(biāo)準(zhǔn)的正弦信號us(t)=usin(lt)，其中，u為電壓幅值，l為周期，記流過非線性負(fù)載周期性、非線性電流為iL(t)，其傅氏級數(shù)展開式[14]為

iL(t) =I1cos(ψ1)sin(lt)+I1sin(ψ1)cos(lt)+

i1p(t)+i1q(t)+ih(t)=i1p(t)+if(t),

k=2,3,…，n,n∈N

(1)

同理，以采樣間隔Ts對iL(t)進行采樣，則

(2)

2 自適應(yīng)LMS電流諧波檢測算法

本文研究的自適應(yīng)LMS電流諧波檢測算法，用代價函數(shù)中對權(quán)值系數(shù)的瞬時梯度來代替最速下降算法[15]中的統(tǒng)計梯度。

由圖1知，權(quán)矢量為w1(n)和w2(n),輸入向量為x1(n)和x2(n),可得

y(n)=x1(n)w1(n)+x2(n)w2(n)

(3)

式(3)滿足輸入量與輸出量之間的線性組合關(guān)系。

設(shè)輸入向量X=(x1,x2,…,xk)T，則有k個權(quán)矢量W=(w1,w2,…,wk)T，由式(3)可得

y(n)=XT(n)W(n)=WT(n)X(n)

(4)

圖2給出了自適應(yīng)LMS算法的結(jié)構(gòu)圖。其中，I為單位矩陣，Z-1I為延時一段時間的矩陣，它屬于一個閉環(huán)控制系統(tǒng),由兩部分組成:一部分是濾波過程，由輸入信號u(n)和權(quán)矢量w(n)計算出輸出信號y(n)，與期望信號d(n)相比較后得到誤差值e(n)；另一部分是自適應(yīng)校正過程，將輸出e(n)與梯度矢量自動更新的權(quán)矢量w(n)反饋到輸入側(cè)得到w(n+1),重新計算輸出值，使之接近于期望信號，得到最佳權(quán)矢量。

圖2 自適應(yīng)LMS算法結(jié)構(gòu)圖

由圖2可知，誤差

e(n)=d(n)-y(n)

將式(4)代入后得

e(n)=d(n)-WT(n)X(n)

(5)

則

e2(n)=d2(n)-2d(n)WT(n)X(n)+
WT(n)X(n)XT(n)W(n)

(6)

對其取數(shù)學(xué)期望后得

E[e2(n)]=E[d2(n)]-2E[d(n)XT(n)]WT(n)+
WT(n)E[X(n)XT(n)]W(n)

(7)

式中，E[d2(n)]、E[X(n)XT(n)]均為常數(shù)。

令WT(n)為變量x，則滿足

E[e2(n)] =E[d2(n)]-2E[d(n)XT(n)]x+
E[X(n)XT(n)]x2(8)

式(8)為二階函數(shù)，且為開口向上的拋物線曲線，函數(shù)存在最小值，故只需求出最小誤差值就可以得到近乎所有的諧波含量。

只需對式(8)中的變量x進行一階求導(dǎo)，得

時，最小誤差值為

(9)

令ε(n)=E[e2(n)],由最速下降公式[15]

w(n+1)=w(n)-Kz(n)

其中，K為誤差系數(shù)，z(n)為第n次時的誤差，可得自適應(yīng)LMS諧波檢測算法的權(quán)值為

(10)

(11)

將式(5)代入式(11)，得

(12)

為了計算簡便，可令ε(n)=e2(n)，則得到瞬時梯度的估計量為

故有

w(n+1)=w(n)+2μe(n)X(n)

(13)

式中,μ控制著自適應(yīng)LMS諧波檢測算法的收斂速度和穩(wěn)態(tài)誤差，μ越大，收斂速度越大，且誤差越大。

由圖2可知，自適應(yīng)LMS算法的權(quán)矢量w(n)只與輸入數(shù)據(jù)u(n)有關(guān)，故當(dāng)u(n)獨立時，w(n)也是獨立的。為獲取μ的取值范圍，計算LMS算法梯度估計的數(shù)學(xué)期望為

-2E[d(n)u(n)-u(n)uT(n)w(n)]=

(14)

式中，R為濾波器輸入向量的相關(guān)矩陣；P為濾波器的輸入與期望響應(yīng)互相關(guān)向量

E[w(n)]-μ{RE[w(n)]-P}=

(I-μR)E[w(n)]-μP

(15)

由矩陣定義得，收斂因子的收斂范圍滿足

0<μ<1/tr(R)

(16)

式中,tr(R)為矩陣R的跡，即等于矩陣中主對角線元素之和。

3 仿真分析

為了驗證本文方法具有更高的精確性與更好的實時性，本文利用MATLAB軟件對自適應(yīng)LMS諧波電流檢測法與傳統(tǒng)的瞬時無功功率p-q諧波電流檢測法進行仿真比較。

(1)p-q諧波電流檢測。瞬時無功功率中的p-q諧波檢測法因其結(jié)構(gòu)簡單，動態(tài)響應(yīng)快的特點而被廣泛應(yīng)用[16]。該方法可以將電路中的電流分解為有功和無功電流，便于對諧波電流進行計算。圖3給出了利用p-q諧波電流檢測法檢測出的結(jié)果。

圖3 p-q法諧波電流檢測結(jié)果

由圖3可見，利用p-q法濾除電流中的諧波含量后，與標(biāo)準(zhǔn)正弦電流信號相比，電流中仍然含有較多的諧波量，誤差較大。可見，p-q法對誤差無法進行自行調(diào)整，而且隨著電路負(fù)載的變化，諧波量也會相應(yīng)變化，無法精確地檢測出電路中的實際諧波量。

(2) 自適應(yīng)LMS諧波檢測法。為驗證本文中自適應(yīng)LMS諧波檢測法的效果，在標(biāo)準(zhǔn)正弦信號中加入隨機噪聲，形成含諧波的待檢測電流信號；再使用本文方法檢測并濾除該合成電流信號中的諧波,將濾除諧波后的信號與標(biāo)準(zhǔn)正弦信號對比較，即可得知檢測的諧波是否完整。自適應(yīng)LMS諧波檢測算法流程如圖4所示。

圖4 自適應(yīng)LMS諧波電流檢測算法流程圖

圖5所示為用MATLAB軟件產(chǎn)生的100s隨機噪聲信號，用來模擬電網(wǎng)或其他電力系統(tǒng)中的電流諧波成分，然后將隨機噪聲信號與標(biāo)準(zhǔn)正弦周期信號疊加，生成電路中包含諧波成分的待檢測電流。

圖5 噪聲波形

圖6給出了利用自適應(yīng)LMS諧波電流檢測法對待檢測電流檢測的結(jié)果。由圖6可見，經(jīng)過自適應(yīng)濾波器濾除諧波后的電流信號不斷地逼近濾除諧波后的電流波形；兩波形之間的誤差也在不斷減小，誤差逐漸地接近于一條直線，這間接地說明，檢測出的諧波量也更接近于實際電路中的諧波量，諧波檢測率高。可見，與傳統(tǒng)的瞬時無功功率p-q諧波電流檢測法相比，自適應(yīng)LMS諧波電流檢測法諧波檢測的精確性更高。

圖6 自適應(yīng)LMS諧波電流檢測結(jié)果

為分析收斂因子μ對濾波器性能的影響。圖7給出了μ=10-3和μ=10-4時自適應(yīng)LMS電流信號輸出與濾除諧波后電流信號輸出的誤差。由圖可見，信號的收斂速度與收斂因子μ的取值有關(guān)，不同的收斂因子對濾波器的性能有很大影響。圖7(a)中，開始時濾波器的初始參數(shù)具有任意性，自適應(yīng)LMS輸出信號與濾除諧波后輸出信號相差較大，但經(jīng)過不斷修正，當(dāng)n=1 000時，自適應(yīng)LMS濾波器電流信號輸出已經(jīng)實現(xiàn)了對濾除諧波后電流信號輸出的精確逼近；圖7(b)中,當(dāng)n=8 000時，自適應(yīng)LMS輸出信號才接近諧波電流信號，自適應(yīng)諧波檢測速度慢，無法快速地在短時間內(nèi)完成收斂。

圖8給出了μ=10-3時自適應(yīng)LMS諧波電流檢測結(jié)果。與圖6比較后發(fā)現(xiàn)，保持誤差量e(n)和輸入矢量X(n)不變，僅增大收斂因子μ的值，可使自適應(yīng)LMS的諧波檢測速度明顯提高，也快速地達到了諧波的最佳狀態(tài)，并在此狀態(tài)下繼續(xù)進行檢測，最終使誤差值不斷趨向于零。

綜上所述，分別對瞬時無功功率理論中p-q法和自適應(yīng)LMS諧波檢測法進行仿真，兩者比較后發(fā)現(xiàn)，p-q法對于不斷變化的負(fù)載電流的檢測能力較差，無法對誤差進行不斷自我調(diào)整，誤差大；自適應(yīng)LMS諧波算法的檢測結(jié)果更加逼近于電路的實際諧波電流量，而且增加收斂因子值后，檢測速度更快，精度更高。

圖7 μ=10-3和μ=10-4時，自適應(yīng)LMS電流信號輸出與濾除諧波后電流信號輸出的誤差

圖8 μ=10-3時，自適應(yīng)LMS諧波電流檢測結(jié)果

4 結(jié) 語

諧波檢測是有源濾波器中必不可少的一部分，諧波檢測方法決定著電力系統(tǒng)和電力電子設(shè)備的濾波效果。本文研究了一種自適應(yīng)LMS諧波檢測法，它由程序中的循環(huán)部分不斷地對參數(shù)進行修正，使之逼近主電路中實際所含的諧波量；同時，由于收斂因子μ決定著信號收斂速度和穩(wěn)定性，通過調(diào)整收斂因子加快了算法的收斂速度，縮短計算時間，及時消除諧波。最后用MATLAB仿真驗證了自適應(yīng)LMS諧波檢測算法的精確性和實時性。

[1] 羅皓文.有源電力濾波器指定次諧波檢測與抑制技術(shù)研究 [D].武漢： 湖北工業(yè)大學(xué)，2013： 2-5.

[2] AKAGI H,KANAZAWA Y,NABAE A.Generalized theroy of the instantaneous reactive power in three-phase circuits [J].IPEC,1983： 1375-1386.

[3] 何英杰,劉進軍,王兆安,等.一種基于瞬時無功功率理論的數(shù)字諧波檢測 [J].電工技術(shù)學(xué)報，2010，25(8)： 185-192.

[4] 楊浩,葉明佳.d-q坐標(biāo)變換的改進型諧波和無功電流檢測法 [J].電源技術(shù)，2014，38(1)： 103-105，115.

[5] 吉曉青.有源電力濾波器的諧波檢測及控制算法的研究 [D].保定： 河北大學(xué)，2015： 13-20.

[6] 盤洪斌,羅安,唐杰，等.一種改進的基于最小二乘法的自適應(yīng)諧波檢測方法 [J].中國電機工程學(xué)報,2008，28(13)： 144-151.

[7] 牛俊萍.基于自適應(yīng)算法的電力系統(tǒng)諧波檢測分析[D].徐州： 中國礦業(yè)大學(xué)，2014： 25-30.

[8] 王群,吳寧,謝品芳.一種基于神經(jīng)元的自適應(yīng)諧波電流檢測法 [J].電力系統(tǒng)自動化，1997,21(10)： 13-16.

[9] PEREIA R R,da SILVA C H,Da SILVA L E B,et al.Improving the convergence time of adaptive notch filters to harmonic detection [C]∥36th Annual Conference on IEEE Industrial Electronics Society.[S.l.]： IEEE,2010： 521-525.

[10] 曹斌芳.自適應(yīng)噪聲抵消技術(shù)的研究 [D].長沙： 湖南大學(xué)，2007： 22-41.

[11] 孟小猛.自適應(yīng)濾波算法研究及應(yīng)用 [D].北京： 北京郵電大學(xué),2010.

[12] 殷桂梁,郭磊,李相男.基于自適應(yīng)陷波器單相電路無功功率測量新方法 [J].南方電網(wǎng)技術(shù)，2013,7(3):85-88.

[13] 陳艷,朱勇,易克初.一種改進的自適應(yīng)濾波器的原理與FPGA實現(xiàn) [J].電子工程師，2004,30(7):47-49，60.

[14] 喻翌,趙海全,何正友.基于ANCT和Adaline兩種自適應(yīng)諧波電流檢測模型的分析 [J].電力系統(tǒng)保護和控制,2013，41(16)： 71-77.

[15] 宗榮芳,李豐林.最速下降與共軛梯度在數(shù)字波束形成中的研究 [J].計算機工程與應(yīng)用，2010,46(19)： 151-153.

[16] 王兆安,李民,卓放.三相電路瞬時無功功率理論的研究 [J].電工技術(shù)學(xué)報,1992(3)： 55-59.

Harmonic Current Detection for Active Power Filter

XUPenghui,XUFei

(School of Electrical Engineering, Shanghai Dianji University, Shanghai 201306)

To filter out harmonics in a circuit, harmonic components in the circuit should be detected accurately. This paper proposes an adaptive least mean square (LMS) harmonic current detection algorithm. A low-pass filter is replaced with an adaptive filter to solve the problem of real-time harmonic detection. The adaptive LMS harmonic current detection algorithm is used for grid voltage distortion. As frequency fluctuation of power parameters has good adaptive ability, it can be applied to three-phase or single-phase circuits. To change convergence factor of the algorithm, detection speed of harmonic can be accelerated. With MATLAB simulation, accuracy and real-time performance of the proposed algorithm is verified.

harmonic detection; self-adapting filter; least mean square; convergence factor

2016 -09 -26

徐鵬輝(1991-)，男，碩士生，主要研究方向為電力電子技術(shù)及其應(yīng)用，E-mail： 1311697578@qq.com

2095 - 0020(2017)01 -0046 - 06

TM 933.1

A