布萊克窗函數(shù)頻域插值電力諧波計(jì)算

唐軼，陳奎，韓智強(qiáng)，陳雷，楊洛

（中國(guó)礦業(yè)大學(xué)信息與電氣工程學(xué)院，徐州 221008）

布萊克窗函數(shù)頻域插值電力諧波計(jì)算

唐軼，陳奎，韓智強(qiáng)，陳雷，楊洛

（中國(guó)礦業(yè)大學(xué)信息與電氣工程學(xué)院，徐州 221008）

提出了一種電力諧波高精度估計(jì)的新方法。首先使用CZT（Chirp Z變換）從含有諧波的電力信號(hào)中高精度提取基波相量，并從中減去基波信號(hào)；然后加布萊克（Blackman）窗截?cái)嚯娏π盘?hào)，用快速傅里葉變換（FFT）計(jì)算出剩余信號(hào)的頻譜；再依據(jù)基波頻率精確估計(jì)出各諧波的頻率；最后依據(jù)各諧波頻率對(duì)Blackman窗函數(shù)在頻域插值，精確計(jì)算出各電力諧波的參數(shù)。仿真試驗(yàn)表明，該方法與加窗雙譜線插值FFT分析電力諧波方法有基本同等的估值精度，但計(jì)算量約為1/2。

線性調(diào)頻Z變換；快速傅里葉變換；諧波分析；布萊克窗；頻域內(nèi)插

電力電子技術(shù)的發(fā)展，非線性負(fù)荷在電力系統(tǒng)中的應(yīng)用的增加，電力系統(tǒng)諧波污染日益嚴(yán)重，諧波已成為影響電能質(zhì)量的主要問(wèn)題[1]。快速傅里葉變換FFT（fast Fourier transform）是諧波分析最快捷的工具[2]。FFT精確分析頻譜的前提是信號(hào)的同步采樣和整周期截?cái)啵鴮?shí)際電網(wǎng)頻率通常在工頻附近波動(dòng)的，雖然從理論上說(shuō)通過(guò)鎖相環(huán)技術(shù)可解決同步采樣和整周期截?cái)鄦?wèn)題[3，4]。電網(wǎng)頻率并非恒定值，鎖相環(huán)響應(yīng)需要時(shí)間，因而難以保證完全同步采樣。普遍采用的思路是通過(guò)選擇譜能量主要集中在主瓣的窗函數(shù)，以減小譜間干擾[5～12]；通過(guò)雙譜線間插值擬合修正，以減小柵欄效應(yīng)，進(jìn)而提高諧波估計(jì)精度[6～13]。加窗插值擬合法有效地提高了諧波估計(jì)的精度，但隨著插值修正曲線擬合函數(shù)的階次增高及諧波含有次數(shù)的增多，諧波估計(jì)精度提高的同時(shí)計(jì)算量大量增加。本文將提出諧波高精度估計(jì)的另一思路，直接對(duì)窗函數(shù)連續(xù)頻譜內(nèi)插精確計(jì)算諧波相量。

1 CZT簡(jiǎn)介

設(shè)x（n）為信號(hào)采樣序列，其長(zhǎng)度為N，n= 0，1，…，n，…N-1，其Z變換為

式中，zk為z平面抽樣點(diǎn)，k=0，1，…，M-1。

如果沿z平面上一段螺線作等分角抽樣，如圖1所示，抽樣點(diǎn)為

式中：A0是起始抽樣半徑；θ為起始抽樣點(diǎn)的角頻率；W0是螺線的伸展率，W0〉1向外延伸，W0〈1向內(nèi)收縮；φ為相鄰兩抽樣點(diǎn)之間的角頻率差，φ〉0表示抽樣路徑逆時(shí)針旋轉(zhuǎn)，φ〈0表示抽樣路徑順時(shí)針旋轉(zhuǎn)。

圖1 CZT在z平面上的螺旋抽樣Fig.1CZT screw slag sampling in z-plane

式中，V（k）=y（k）h（k）可以看成是一個(gè)數(shù)字網(wǎng)絡(luò)的單位脈沖響應(yīng)的輸出，可用FFT快速求得：V（k）= IFFT[FFT[y（k）]FFT[h（k）]]。由于式（4）的單位脈沖響應(yīng)h（n）的頻率隨時(shí)間成線性增加的線性調(diào)頻信號(hào)相似，因此這種變換就稱為線性調(diào)Z變換（Chirp-Z transform）或CZT[14]。

考察式（4）可見(jiàn)：其一，可以計(jì)算單位圓上任意一段線段上的頻率樣本，而且可以任意給定起始頻率θ，終止頻率θ+（M-1）φ；其二，頻率樣本的點(diǎn)數(shù)M可以任意給定，而頻率抽樣間隔φ隨著M的增大而縮小，這一特點(diǎn)可達(dá)到細(xì)化頻率、密集譜線、減小柵欄效應(yīng)的目的。這兩個(gè)特點(diǎn)顯然具有精確地從含有諧波的電力信號(hào)中提取基波參數(shù)的優(yōu)勢(shì)。根據(jù)標(biāo)準(zhǔn)[15]，電網(wǎng)最大允許頻率偏差為±0.5Hz，因此，可設(shè)待分析頻段的起始頻率為θ=2π×49.5，終止頻率為θ+（M-1）φ=2π×50.5，即分析帶寬為1 Hz。如果信號(hào)采樣頻率為fs=6 400 Hz，若取M=1 280（10個(gè)周波），則頻率抽樣間隔φ=2π/ 127 9=2π×7.818 61×10-4。若取模最大值所對(duì)應(yīng)的頻率作為基波頻率估計(jì)值f1，頻率測(cè)量的誤差為|Δf|=（φ/2π）/2≤3.909 30×10-4Hz。因此，按圖2計(jì)算流程可計(jì)算得基波信號(hào)的頻譜X（k），進(jìn)而分別計(jì)算得基波的幅值A(chǔ)1、頻率f1和相位φ1分別為

圖2 快速CZT算法流程Fig.2Fast CZT algorithm process

2 窗函數(shù)頻域插值諧波相量的計(jì)算

單頻電力諧波信號(hào)為

式中：p∈[0，P]，p為諧波次數(shù)，是正整數(shù)，P為最高諧波次數(shù)；Ap為第p次諧波幅值；f1為電力信號(hào)基波角頻率；φp為第p次諧波初相角。

以采樣頻率fs=1/Ts對(duì)式（8）離散化，得數(shù)字序列為

式中：Ts為采樣周期；ωp=2πpf1Ts為p次諧波數(shù)字角頻率。x（n）經(jīng)Fourier變換后的連續(xù)頻譜為

Blackman窗函數(shù)[13]為

可預(yù)定以電網(wǎng)的基波額定頻率f1′=50 Hz的2i倍（i=1，2，…）采樣，即采樣頻率fs=1/Ts=50× 2i。這樣，電網(wǎng)基波和諧波的額定頻率是與Blackman窗函數(shù)同步采樣并能整周期截?cái)嗟摹tωp′=2πpf1′Ts是p次電力諧波信號(hào)與Blackman窗函數(shù)同步采樣并整周期截?cái)嗟臄?shù)字角頻率，以角頻率ωp′在Blackman窗函數(shù)連續(xù)頻譜上在頻域內(nèi)抽樣，必然有WBla（ωp-ωp′）=1，如圖3所示（僅示出正頻段）。

圖3 同步采樣和正周期截?cái)嗪蟮念l譜Fig.3Spectrum with synchronous sampling and integral period truncation

同時(shí)p次諧波信號(hào)的實(shí)際數(shù)字角頻率ωp= 2πpf1Ts也在同一Blackman窗函數(shù)連續(xù)頻譜上在頻域內(nèi)抽樣求得WBla（ωp-ωp′），如圖4所示（僅示出正頻段）。

由此可求得p次電力諧波。

幅值校正系數(shù)為

圖4 非同步采樣和非整周期截?cái)嗪蟮念l譜Fig.4Spectrum with nonsynch ronous sampling or non-integral period truncating

對(duì)xBla（n）的FFT實(shí)際上就是以Δω=2π/N等間隔對(duì)XBla（ω）進(jìn)行抽樣的結(jié)果，即

-0.5Δω≤δωp≤0.5Δω，將式（18）代入式（13）和式（14）得p次諧波，則p次諧波的幅值和相位分別為

式中，φp′為電力p次諧波信號(hào)的FFT離散頻譜的主譜XBla（kpΔω）的相位。

3 仿真實(shí)驗(yàn)

本文采用的仿真軟件為Matlab7.7.0（R2008b）。仿真采用的信號(hào)模型為

式中：基波頻率f1為50.4Hz；采樣頻率fs為6400Hz，截?cái)嘈盘?hào)的數(shù)據(jù)長(zhǎng)度N為1 280點(diǎn)；電力信號(hào)的參數(shù)如表1所示。

表1 仿真信號(hào)的基波及諧波成分Tab.1Components of the simulated harmonic signal

利用CZT，取M=N=1 280，高精度估計(jì)出基波頻率、幅值和相位分別為f1=50.400 000 Hz，A1= 220.004 3 V，φ1=0.499 3°。

如圖5所示，本文提出的分析方法與基于Blackman窗的雙譜線插值擬合法具有基本同等的分析精度。用同一臺(tái)計(jì)算機(jī)和同一個(gè)仿真軟件進(jìn)行計(jì)算，兩種方法計(jì)算耗時(shí)分別為7.889 ms和15.830 ms，本文所提方法計(jì)算時(shí)間減少了一半。

圖5 幅值與相位絕對(duì)誤差對(duì)比Fig.5Comparisons of absolute errors in calculating amplitude and phase

4 結(jié)語(yǔ)

本文提出了一種電力諧波分析的新思路，即先高精度地估計(jì)出含有諧波的電力信號(hào)中的基波相量，進(jìn)而，通過(guò)窗函數(shù)頻域內(nèi)插值的方法精確計(jì)算出各次諧波的校正系數(shù)，分別計(jì)算出各次諧波的參數(shù)。通過(guò)同一電力信號(hào)的仿真試驗(yàn)證明，它與加Blackman窗雙譜線插值擬合估計(jì)電力諧波的分析方法有基本等同的估計(jì)精度，而這種算法在計(jì)算速度上有明顯的優(yōu)勢(shì)，是一種很有實(shí)用價(jià)值的電力諧波進(jìn)行高精度估計(jì)的方法。

[1]Mahmoud Aly A.Effects of power system harmonics on power system equipment and loads[J].IEEE Trans on Power Apparatus and Systems，1985，104（9）：2555-2563.

[2]Heydt G T，F(xiàn)jeld P S，Liu C C，et al.Applications of the windowed FFT to electric power quality assessment[J]. IEEE Trans on Power Delivery，1999，14（4）：1411-1416.

[3]Hidalgo R M，F(xiàn)ernandez J G，Rivera R R，et al.A simple adjustable window algorithm to improve FFT measurements[J].IEEE Trans on Instrumentation and Measurement，2002，51（1）：31-36.

[4]門長(zhǎng)有，王榮華，譚年熊（Men Changyou，Wang Rong hua，Tan Nianxiong）.一種用于諧波測(cè)量的全數(shù)字同步采樣算法（An entire digital synchronous sampling algorithm for harmonic measurements）[J].電力系統(tǒng)自動(dòng)化（Automation of Electric Power Systems），2008，32（22）：83-86.

[5]吳靜，趙偉（Wu Jing，Zhao Wei）.一種用于分析電網(wǎng)諧波的多譜線插值算法（An algorithm of MICA for analyzing harmonics in power system）[J].中國(guó)電機(jī)工程學(xué)報(bào)（Proceedings of the CSEE），2006，26（8）：55-60.

[6]Offelli C，Petri D.Interpolation techniques for real-time multifrequency waveform analysis[J].IEEE Trans on Instrumentation and Measurement，1990，39（1）：106-111.

[7]Agrez D.Weighted multipoint interpolated DFT to improve amplitude estimation of multifrequency signal[J].IEEE Trans on Instrumentation and Measurement，2002，51（2）：287-292.

[8]Xue H，Yang R.Optimal interpolating windowed discrete Fourier transform algorithms for harmonic analysis in power systems[J].IEE Proceedings of Generation，Transmission and Distribution，2003，150（5）：583-587.

[9]曾博，滕召勝，高云鵬，等（Zeng Bo，Teng Zhaosheng，Gao Yunpeng，et al）.基于Rife-Vincent窗的高準(zhǔn)確度電力諧波相量計(jì)算方法（An accurate approach for power harmonic phasor calculation based on Rife-Vincent window）[J].電工技術(shù)學(xué)報(bào)（Transactions of China Electrotechnical Society），2009，24（8）：154-158.

[10]曾博，滕召勝，溫和，等（Zeng Bo，Teng Zhaosheng，Wen He，et al）.萊夫-文森特窗插值FFT諧波分析方法（An approach for harmonic analysis based on Rife-Vincent window interpolation FFT）[J].中國(guó)電機(jī)工程學(xué)報(bào)（Proceedings of the CSEE），2009，29（10）：115-120.

[11]卿柏元，滕召勝，高云鵬，等（Qing Baiyuan，Teng Zhaosheng，Gao Yunpeng，et al）.基于Nuttall窗雙譜線插值FFT的電力諧波分析方法（An approach for electrical harmonic analysis based on Nuttall window doublespectrum-line interpolation FFT）[J].中國(guó)電機(jī)工程學(xué)報(bào)（Proceedings of the CSEE），2008，28（25）：153-158.

[12]Reljin I S，Reljin B D，Papic V D.Extremely flat-top windows for harmonic analysis[J].IEEE Trans on Instrumentation and Measurement，2007，56（3）：1025-1041.

[13]龐浩，李東霞，俎云霄，等（Pang Hao，Li Dongxia，Zu Yunxiao，et al）.應(yīng)用FFT進(jìn)行電力系統(tǒng)諧波分析的改進(jìn)算法（An improved algorithm for harmonic analysis of power system using FFT technique）[J].中國(guó)電機(jī)工程學(xué)報(bào)（Proceedings of the CSEE），2003，23（6）：50-54.

[14]Rabiner L R，Schafer R W，Rader C M.The Chirp Z-transform algorithm[J].IEEE Trans on Audio and Electro-acoustics，1969，17（2）：86-92.

[15]GB/T 15945 2008，電能質(zhì)量電力系統(tǒng)頻率允許偏差[S].

New Calculation Method of Electric Power Harmonic in Frequency Domain Interpolation of Blackman Windows

TANG Yi，CHEN Kui，HAN Zhi-qiang，CHEN Lei，YANG Luo
（School of Information and Electrical Engineering，China University of Mining and Technology，Xuzhou 221008，China）

In this paper，a new method with high accuracy for electric harmonic estimation is presented.First，the CZT（Chirp Z-transform）is used to extract the parameters（amplitude，frequency and phase）of fundamental frequecy signals precisely from the multifrequency signals of power system.Second，after subtracting the fundamental frequency signal from the multifrequency signals，the harmonic signals are truncated with the Blackman window and their spectrums are calculated by fast Fourier transform（FFT）.Then，according to the fundamental frequency，each frequency harmonic is calculated.At last，interpolating on the frequency domain of the Blackman window，the parameters（amplitude and phase）of each frequency harmonic signals can be calculated precisely.The simulation results show that the proposed method almost has the same estimation accuracy with the method which is based on double-spectrum-line interpolation FFT，but the new method only costs about half computation times.

Chirp Z transform；fast Fourier transform（FFT）；harmonic analysis；Blackman window；interpolation on frequency domain

TM935-470.40

A

1003-8930（2013）03-0142-05

唐軼（1957—），男，博士，教授，博士生導(dǎo)師，從事電力系統(tǒng)自動(dòng)化和電能質(zhì)量控制領(lǐng)域的教學(xué)、研究工作。Email：tangyi@cumt.edu.cn

2011-06-14；

2011-08-17

江蘇省科技支撐計(jì)劃（工業(yè)）資助項(xiàng)目（BE2012042）

陳奎（1973—），男，博士，副教授，從事電力系統(tǒng)自動(dòng)化方面的教學(xué)和研究工作。Email：chenkui73@cumt.edu.cn

韓智強(qiáng)（1987—），男，碩士研究生，從事電能質(zhì)量的監(jiān)測(cè)與控制的研究。E-mail：hzqiang87@163.com