曹 敏,林 聰,李 波,劉清蟬,梅二召,王 行,謝 濤,李 川

(1.云南電網(wǎng)有限責(zé)任公司電力科學(xué)研究院，云南 昆明 650217；2.中國(guó)南方電網(wǎng)公司電能計(jì)量重點(diǎn)實(shí)驗(yàn)室，云南 昆明 650217；3.昆明理工大學(xué)信息工程與自動(dòng)化學(xué)院，云南 昆明 650500)

小波包結(jié)合奇異值分解在通信信道降噪中的研究

曹 敏1,2,林 聰1,2,李 波1,2,劉清蟬1,2,梅二召3,王 行3,謝 濤3,李 川3

(1.云南電網(wǎng)有限責(zé)任公司電力科學(xué)研究院，云南 昆明 650217；2.中國(guó)南方電網(wǎng)公司電能計(jì)量重點(diǎn)實(shí)驗(yàn)室，云南 昆明 650217；3.昆明理工大學(xué)信息工程與自動(dòng)化學(xué)院，云南 昆明 650500)

針對(duì)營(yíng)配一體化通信信號(hào)受噪聲干擾的問(wèn)題，為了提高營(yíng)配一體化通信信號(hào)的降噪效果，提出一種小波包結(jié)合奇異值分解的降噪方法。對(duì)電力線高速寬帶載波信號(hào)進(jìn)行小波包分解，有選擇地對(duì)小波包分解系數(shù)進(jìn)行奇異值分解和變換，并與小波包的多尺度時(shí)頻分析及重構(gòu)算法相結(jié)合，將其應(yīng)用于通信信道降噪的研究中。針對(duì)不同的背景環(huán)境，挑選三個(gè)采集地點(diǎn)，采用TPS-2024數(shù)字信號(hào)示波器在現(xiàn)場(chǎng)采集數(shù)據(jù)，并對(duì)采集到的通信信號(hào)分別進(jìn)行小波包降噪和小波包結(jié)合奇異值降噪。測(cè)試結(jié)果顯示，小波包結(jié)合奇異值降噪較小波降噪的信噪比(SNR)提高10 dB左右，且均方誤差(MSE)都有所降低。通過(guò)對(duì)比測(cè)試數(shù)據(jù)結(jié)果可知，采用小波包結(jié)合奇異值分解的降噪方法可以有效地保留信號(hào)的有效成分，降噪效果明顯優(yōu)于小波包降噪。

電力； 智能電網(wǎng)； 通信； 網(wǎng)絡(luò)； 載波； 降噪； 數(shù)據(jù)

0 引言

營(yíng)配一體化是實(shí)現(xiàn)智能電網(wǎng)、推進(jìn)企業(yè)資源計(jì)劃系統(tǒng)漸進(jìn)的過(guò)程中，在技術(shù)層面上不可或缺的主體功能[1-3]。營(yíng)配一體化系統(tǒng)是以統(tǒng)一的電力信息模型與電網(wǎng)設(shè)備、基礎(chǔ)資料以及拓?fù)潢P(guān)系為基礎(chǔ)的系統(tǒng)，是利用高新信息技術(shù)，應(yīng)用于線損四分管理、配網(wǎng)建設(shè)規(guī)劃、客戶停電管理等領(lǐng)域的一體化信息系統(tǒng)[4-5]。電力線高速寬帶通信(high-speed power line communication，H-PLC)作為電力通信網(wǎng)絡(luò)的重要組成部分，不僅能夠滿足當(dāng)前智能電網(wǎng)中各種智能終端的數(shù)據(jù)接入，還能避免由于使用GPRS、ADSL等公共網(wǎng)絡(luò)資源帶來(lái)的通信效率低、穩(wěn)定性差和安全性問(wèn)題。實(shí)現(xiàn)H-PLC通信的技術(shù)難點(diǎn)之一是解決噪聲干擾問(wèn)題[6-7]。

小波包具有良好的時(shí)域局域化功能和優(yōu)異的帶通濾波功能，因此，小波包分解被廣泛地應(yīng)用到信號(hào)降噪的研究中。2006年，史賢俊等提出了一種基于 Shannon 熵準(zhǔn)則的最優(yōu)小波包基信號(hào)降噪算法，對(duì)某型導(dǎo)彈慣導(dǎo)系統(tǒng)陀螺儀信號(hào)進(jìn)行降噪[8]；2010年，于文新等采用類 Sigmoid 函數(shù)作為改進(jìn)的軟閾值函數(shù)，通過(guò)改變閾值來(lái)提高信號(hào)的信噪比和減小均方根誤差[9]；2011年，Sheybani等將自適應(yīng)小波包變換應(yīng)用于電力線通信中，能夠有效地去除通信信號(hào)中的背景噪聲和脈沖噪聲[10]。為進(jìn)一步提高小波包降噪精度，陳桂明等將小波包和奇異值結(jié)合，并分別進(jìn)行仿真試驗(yàn)、聲音信號(hào)處理以及飛行數(shù)據(jù)的降噪，均取得顯著的降噪效果[11-13]。

本文應(yīng)用小波包變換，對(duì)中壓信道信號(hào)進(jìn)行分解，根據(jù)分解信號(hào)所在的頻帶范圍，針對(duì)性地對(duì)不同頻段的小波包系數(shù)進(jìn)行奇異值分解(singular value decompositions,SVD)和變換，然后重構(gòu)信號(hào)。通過(guò)分析去噪后信號(hào)的信噪比(signal noise ratio,SNR)以及均方誤差(mean square error,MSE)，來(lái)比較小波包結(jié)合奇異值分解對(duì)H-PLC通信信號(hào)的降噪效果。

1 小波包降噪和奇異值分解算法

1.1 小波包降噪算法

當(dāng)進(jìn)行小波分解時(shí)，所得到的高頻信號(hào)和低頻信號(hào)在這個(gè)頻域中分別占據(jù)一半的頻帶寬度；再一次用小波分解時(shí)，又把低頻部分分為兩個(gè)同樣寬的頻帶。而小波包分解需要對(duì)小波分解過(guò)程中沒(méi)被分解的高頻信號(hào)進(jìn)行再分解，其目的是根據(jù)代價(jià)函數(shù)選取合適的基波函數(shù)，然后對(duì)小波包基函數(shù)分解后的信號(hào)進(jìn)行時(shí)頻分析[14-15]。因此，小波包分解不僅具有小波的優(yōu)點(diǎn)，還擴(kuò)展了小波分析的應(yīng)用。

對(duì)于給定的正交尺度函數(shù)φ(x)及其對(duì)應(yīng)的小波函數(shù)Ψ(x)，有雙尺度方程：

(1)

式中:{gk}和{hk}為多尺度分析(multiscalanalysis,MSA)中定義的共軛濾波器。

定義μ0=φ(x)、μ1=Ψ(x)，推導(dǎo)出μn(x)。

(2)

式中：函數(shù)簇{μn(x)|n∈Z+}為相對(duì)于正交尺度函數(shù)φ(x)的正交小波包。

(3)

綜上所述，小波包分解可以逐步將頻率分解成若干頻段的信號(hào)。一個(gè)含噪聲的中壓載波信道信號(hào)的模型可以表示為：

y(n)=x(n)+e(n)n=1,2,…,N

(4)

式中:y(n)為有噪聲的混合信號(hào)，即取樣的時(shí)域信號(hào)；x(n)為真實(shí)的無(wú)噪信號(hào)；e(n)為噪聲。

小波包降噪就是根據(jù)噪聲與信號(hào)在不同頻帶上的不同特性，去除各尺度上的小波分量，尤其是那些噪聲分量占主導(dǎo)地位的小波分量，則保留下來(lái)的小波包系數(shù)即為原始小波包系數(shù)。小波包降噪的關(guān)鍵是濾除由噪聲產(chǎn)生的小波包分解系數(shù)[16]。采用系統(tǒng)默認(rèn)閾值對(duì)小波包分解的系數(shù)進(jìn)行處理，就是將大于閾值的小波包系數(shù)保留下來(lái)重構(gòu)信號(hào)，而對(duì)小于或等于默認(rèn)閾值的小波包系數(shù)進(jìn)行置零處理。

1.2 奇異值分解

SVD是一種非線性濾波。根據(jù)文獻(xiàn)[17]～文獻(xiàn)[19]可知，在SVD理論中，任何m×n階矩陣A的SVD為：

A=USVT

(5)

由于Λ是一對(duì)角矩陣，因此SVD可以將一個(gè)秩為k的m×n階矩陣A表示為k個(gè)秩為1的m×n階子矩陣的和。其中，每個(gè)子矩陣由2個(gè)特征矢量和權(quán)值相乘得到：

(6)

在處理信號(hào)的過(guò)程中，根據(jù)信號(hào)特征合理地將對(duì)角矩陣Λ中的較小元素置零，然后代入式(7)：

A0=US0VT

(7)

式中:S0為去除小奇異值后的對(duì)角矩陣；A0為處理信號(hào)后的去噪矩陣。在中壓信道信號(hào)中，噪聲主要集中在高頻帶中，在總信號(hào)能量分布中所占比重較小，因此可以利用奇異值分頻段地去除信號(hào)噪聲[18]。

1.3 評(píng)價(jià)指標(biāo)

(8)

(9)

2 白邑線營(yíng)配一體化通信平臺(tái)

本試驗(yàn)數(shù)據(jù)采集依托35 kV白邑變電站的營(yíng)配一體化智能信息通信平臺(tái)。白邑變電站擁有5條10 kV線路，該項(xiàng)目選擇從該變電站輸出的10 kV白邑線，該線路是5條線路中最長(zhǎng)的一條。主干線路全長(zhǎng)15.487 km，主干分為3段，支線全長(zhǎng)25.854 km，安裝有變壓器56臺(tái)，變壓器容量為3 691 kVA，供電半徑為17.438 km。線路途徑眾多村落，有平地、丘陵、山地等多種地形；并且最長(zhǎng)支線距離為3 km。從以上數(shù)據(jù)可知，該10 kV線路具有較強(qiáng)的代表性。

該項(xiàng)目設(shè)計(jì)的營(yíng)配一體化智能信息通信平臺(tái)，利用H-PLC通信設(shè)備，通過(guò)現(xiàn)有的電纜線路實(shí)現(xiàn)通信；耦合器在變壓器和中壓斷路器處作耦合改造；利用濾波器濾除雜信，建立起通信通道。營(yíng)配一體化智能信息通信平臺(tái)組網(wǎng)示意圖如圖1所示。

圖1 智能信息通信平臺(tái)組網(wǎng)示意圖

營(yíng)配一體化智能信息通信系統(tǒng)由主站(服務(wù)器)、H-PLC頭端、中壓信號(hào)耦合器、H-PLC中繼以及配電自動(dòng)化終端(FTU、TTU等)和低壓集抄系統(tǒng)組成。H-PLC頭端設(shè)備通過(guò)變電站內(nèi)部的光纖或局域網(wǎng)接入到電力系統(tǒng)中，并與營(yíng)配一體化服務(wù)器通信。H-PLC信號(hào)通過(guò)中壓信號(hào)耦合器耦合到10 kV的線路上，遠(yuǎn)端的H-PLC中繼設(shè)備將信號(hào)放大轉(zhuǎn)發(fā)，送至各個(gè)節(jié)點(diǎn)的末端設(shè)備上；末端的H-PLC設(shè)備通過(guò)RJ45以太網(wǎng)接口或RS-232、RS-485等接口與各種配電自動(dòng)化設(shè)備和低壓集抄設(shè)備進(jìn)行通信，并將數(shù)據(jù)傳回營(yíng)配一體化服務(wù)器或向設(shè)備發(fā)送來(lái)自服務(wù)器端的操作指令。該線路覆蓋主干線需要使用H-PLC頭端4套，H-PLC中繼12套，智能開(kāi)關(guān)采集終端3套，單相10 kV電壓互感器10臺(tái)以及支線的TTU、DTU、FTU接口設(shè)備和寬帶通信模塊等設(shè)備15套。

3 測(cè)試數(shù)據(jù)

在白邑線營(yíng)配一體化通信平臺(tái)的H-PLC線路中，針對(duì)不同的背景環(huán)境，挑選3個(gè)采集地點(diǎn)。根據(jù)采集地點(diǎn)的實(shí)際情況采取相應(yīng)的采集方式，實(shí)際采集點(diǎn)及采集方式如下。

①第一采集點(diǎn)距離終端基站680 m，周邊基本沒(méi)有其他設(shè)施的干擾，采用單耦合方式，信號(hào)較強(qiáng)。

②第二采集點(diǎn)距離第一變電站1.6 km，由于該點(diǎn)附近有500 kV的高壓輸電線路通過(guò)，故采用負(fù)控終端以及集差雙耦合技術(shù)來(lái)弱化干擾。因此，針對(duì)此環(huán)境，采集該處的信號(hào)具有相對(duì)的代表性。

③第三采集點(diǎn)距離第二變電站1.4 km，附近有800 kV的高壓輸電線路通過(guò)，采用的也是負(fù)控終端以及集差雙耦合技術(shù)來(lái)弱化干擾。因此，針對(duì)此環(huán)境，采集該處的信號(hào)具有相對(duì)的代表性。

采用TPS-2024數(shù)字信號(hào)示波器在現(xiàn)場(chǎng)采集數(shù)據(jù)。在不同環(huán)境下均采集2組數(shù)據(jù)，按照相關(guān)規(guī)范要求，專業(yè)技術(shù)人員將示波器連接在PLC系統(tǒng)回路上，采集到混合信號(hào)后，將該混合信號(hào)的時(shí)域波形圖以及數(shù)據(jù)分別存儲(chǔ)起來(lái)?，F(xiàn)場(chǎng)分別采集多組數(shù)據(jù)，以便統(tǒng)計(jì)并分析比較。

4 測(cè)試數(shù)據(jù)分析

根據(jù)TPS-2024數(shù)字信號(hào)示波器保存數(shù)據(jù)的特點(diǎn)，對(duì)采集到的通信信號(hào)分別進(jìn)行小波包降噪和小波包結(jié)合奇異值降噪。圖2給出了采集點(diǎn)的采集信號(hào)、小波包降噪、小波包結(jié)合奇異值降噪的測(cè)試結(jié)果示意圖。

圖2 采集點(diǎn)測(cè)試圖

對(duì)比圖2(b)和圖2(c)可知，單獨(dú)由小波包降噪得到的曲線毛刺較多，整體上明顯沒(méi)有小波包結(jié)合奇異值降噪后的曲線光滑。對(duì)以上處理結(jié)果進(jìn)行量化指標(biāo)評(píng)判濾波的優(yōu)劣，如表1所示。

表1 降噪結(jié)果對(duì)比表

從表1可以看出，三個(gè)采集點(diǎn)的小波包結(jié)合奇異值降噪結(jié)果較小波降噪的信噪比提高10 dB左右，且均方誤差都有所降低。同時(shí)，對(duì)比第一采集點(diǎn)和第二采集點(diǎn)可以發(fā)現(xiàn)，受環(huán)境噪聲的影響，第二采集點(diǎn)的信噪比較小。但考慮到第一采集點(diǎn)和第二采集點(diǎn)信號(hào)采集方式的差異，對(duì)比第三采集點(diǎn)和第二采集點(diǎn)可以得出，第三采集點(diǎn)的信噪比較小。從整體特征分析可知，環(huán)境噪聲使得H-PLC通信信號(hào)的信噪比有所降低。采用小波包和小波包結(jié)合奇異值降噪都可以提高通信信號(hào)的信噪比，并且降低通信信號(hào)的均方誤差。綜合分析可以得出，小波包結(jié)合奇異值處理中壓信道信號(hào)的降噪效果優(yōu)于單獨(dú)采用小波包處理的效果。

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

為了提高營(yíng)配一體化通信信號(hào)的降噪效果，提出了小波包結(jié)合奇異值降噪算法。在小波包分解得到各頻帶信號(hào)的基礎(chǔ)上，對(duì)前兩個(gè)頻帶保留小波包系數(shù)，后幾個(gè)頻帶的小波包系數(shù)采用奇異值分解、變換和重構(gòu)；然后對(duì)奇異值分解降噪后的頻帶與剩余頻帶信號(hào)用小波包重構(gòu)，得出去噪后的通信信號(hào)。

通過(guò)對(duì)比直接采用默認(rèn)閾值的小波包降噪與采用小波包結(jié)合奇異值降噪的結(jié)果可知，采用小波包結(jié)合奇異值分解降噪算法，能夠提高通信信號(hào)的信噪比，降低通信信號(hào)的均方誤差。因此，小波包結(jié)合奇異值分解是提高營(yíng)配一體化通信信號(hào)質(zhì)量的一種有效降噪方法。

[1] 何劍輝,麥楚良.基于GIS營(yíng)配一體化中信息共享的研究[J].科技資訊,2013(5):24.

[2] 張東霞,苗新,劉麗平,等.智能電網(wǎng)大數(shù)據(jù)技術(shù)發(fā)展研究[J].中國(guó)電機(jī)工程學(xué)報(bào),2015(1):2-12.

[3] 黃少雄,黃太貴,程栩.智能變電站信息一體化平臺(tái)建設(shè)方案研究[J].東北電力技術(shù),2012(2):21-24.

[4] 畢建軍,蔡明明.電網(wǎng)營(yíng)銷(xiāo)與配電一體化信息通信技術(shù)平臺(tái)研究[J].中國(guó)電力教育,2010(S2):550-554.

[5] 李瀚.探析營(yíng)配一體化系統(tǒng)及其信息集成平臺(tái)[J].科技創(chuàng)新與應(yīng)用,2012(21):46-47.

[6] 謝志遠(yuǎn),聶立賢.基于PLC通信的噪聲特性及多載波調(diào)制技術(shù)研究[J].電子技術(shù)應(yīng)用,2009,35(11):93-96.

[7] 林虹,林東.PLC信道噪聲建模及噪聲發(fā)生器的設(shè)計(jì)[J].電力系統(tǒng)通信,2011,32(7):61-65.

[8] 史賢俊,林颯,李瑞亮.基于最優(yōu)小波包基的信號(hào)去噪算法及其應(yīng)用[J].海軍航空工程學(xué)院學(xué)報(bào),2006,21(5):506-509.

[9] 于文新,張謙.基于改進(jìn)閾值消噪算法的小波包降噪[J].通信技術(shù),2010,43(6):7-9.

[10]SHEYBANI E,RASHIDI N.Application of wavelet packets in power line communications[A].The IASTED International Conference on Power and Energy Systems,EuroPES,2011:211-215.

[11]陳桂明,楊斌,何慶飛.一種基于奇異值的小波包改進(jìn)算法[J].噪聲與振動(dòng)控制,2013(1):63-67,71.

[12]盧德林,郭興明.基于奇異譜分析的心音信號(hào)小波包去噪算法研究[J].振動(dòng)與沖擊,2013,32(18):63-69.

[13]葉艷,王少云,孫永帥.基于奇異值分解和小波分析的飛行數(shù)據(jù)野值綜合處理方法[J].光學(xué)儀器,2014,36(5):394-398.

[14]CHING P C,SO H C,WU S Q.On wavelet denoising and its applications to time delay estimation[J].IEEE Transaction Signal Processing,1999,47(10):2879-2882.

[15]VIDAKOVIC B,LOZOYA C B.On time-dependent wavelet denoising[J].IEEE Trans Signal Processing,1998,46(9):2549-2551.

[16]李陳靈.基于小波馬爾科夫鏈的電力線噪聲建模研究[D].南京：南京理工大學(xué),2014.

[17]何慶飛,姚春江,陳桂明,等.基于改進(jìn)小波包奇異值法的齒輪泵振動(dòng)信號(hào)去噪[J].機(jī)械科學(xué)與技術(shù),2012,31(9):1445-1448.

[18]郭昊坤.低壓電力線通信信道噪聲建模研究[D].南京：南京理工大學(xué),2013.

[19]吳軍基,郭昊坤,孟紹良,等.電力線通信信道背景噪聲建模研究[J].電力系統(tǒng)保護(hù)與控制,2011,39(23):6-10.

[20]左飛飛,王海彬,馬捷,等.基于小波變換的改進(jìn)去噪閾值函數(shù)[J].探測(cè)與控制學(xué)報(bào),2015(1):80-85.

Study on the Application of Wavelet Packet+Singular Value Decomposition in Noise Reduction for Communication Channel

CAO Min1,2,LIN Cong1,2,LI Bo1,2,LIU Qingchan1,2,MEI Erzhao3,WANG Xing3,XIE Tao3,LI Chuan3

(1.Electric Power Research Institute,Yunnan Power Grid Co.,Ltd,Kunming 650217，China; 2.Key Laboratory of CSG for Electric Power Measurement，Kunming 650217，China; 3.Faculty of Information Engineer and Automation,Kunming University of Science and Technology,Kunming 650500，China)

To solve the problem of noise interference in communication of integrated marketing and distribution,and in order to improve the denoising effect of communication signals in integration of marketing and distribution,the method which is combining wavelet packet and singular value decomposition is proposed to reduce the noise.The high speed broadband carrier signals of the power line is decomposed by wavelet packet,and the singular value decomposition and transformation are selectively conducted for wavelet packet decomposition coefficients;and combined with the wavelet packet multi-scale time-frequency analysis and reconstruction algorithm,then applied in the research on noise reduction of communication channel.According to the different background environment,three of the collection sites are selected,using TPS-2024 digital signal oscilloscope to collect data in the field,and the communication signals collected are analyzed by denoising of wavelet packet and wavelet packet combining with singular value.The results show that the SNR of wavelet packet + singular value de-noising can be improved about 10 dB,and the mean square error (MSE) is reduced.The results of comparing the test show that using wavelet packet + singular value decomposition method can effectively keep the effective components of the signal,the noise reduction result is better than wavelet packet denoising.

Power； Smart grid； Communication； Network; Carrier； Noise reduction; Data

資金項(xiàng)目:國(guó)家自然科學(xué)基金地區(qū)科學(xué)基金項(xiàng)目(51567013)

曹敏(1961—)，男，學(xué)士,教授級(jí)高級(jí)工程師，主要從事電能表、互感器、在線監(jiān)測(cè)裝置、光纖光柵儀器儀表檢測(cè)技術(shù)方向的研究。E-mail：547365190@qq.com。 李川(通信作者)，男，博士，教授,博士生導(dǎo)師，主要從事傳感器的研制與檢測(cè)方向的研究。E-mail：1625677252@qq.com。

TH6;TP273

A

10.16086/j.cnki.issn 1000-0380.201702014

修改稿收到日期:2016-04-08