利用小波能量特征的增長型自組織神經(jīng)網(wǎng)絡(luò)同調(diào)機(jī)組分群方法*

楊越，王濤，顧雪平，岳賢龍，徐振華，邱麗君

（1.華北電力大學(xué)新能源電力系統(tǒng)國家重點(diǎn)實(shí)驗(yàn)室，河北保定071003；2.國網(wǎng)福建省電力有限公司電力科學(xué)研究院，福州350007；3.中國電力科學(xué)研究院，北京100192）

0 引 言

電力系統(tǒng)中發(fā)電機(jī)功角不是一成不變的，而是隨著外界擾動(dòng)和運(yùn)行狀態(tài)的變化而變化，發(fā)電機(jī)功角在變化過程中，不同的機(jī)組可能具有一致的動(dòng)態(tài)變化特性，即機(jī)組之間存在著同調(diào)現(xiàn)象［1］?？梢砸罁?jù)電力系統(tǒng)中機(jī)組之間的同調(diào)性將整個(gè)系統(tǒng)分成各個(gè)互相解耦的子系統(tǒng)，進(jìn)而簡化系統(tǒng)模型。比如：當(dāng)系統(tǒng)發(fā)生低頻振蕩時(shí)，若存在著弱互聯(lián)區(qū)域，通常應(yīng)用機(jī)組之間的同調(diào)現(xiàn)象將系統(tǒng)分區(qū)域研究［2］；在EEAC法中，也是先進(jìn)行發(fā)電機(jī)的兩群假設(shè)，將發(fā)電機(jī)分為臨界機(jī)群和剩余機(jī)群［3］。除此之外，當(dāng)電網(wǎng)受到大沖擊，為了避免發(fā)生連鎖故障［4］、防止大停電而進(jìn)行的解列控制［5］也多以同調(diào)機(jī)組分群為基礎(chǔ)。所以精準(zhǔn)識(shí)別電力系統(tǒng)的同調(diào)特性具有重要意義。

目前，機(jī)組分群方法主要有兩類，第一類是建立系統(tǒng)發(fā)電機(jī)模型，然后選出一個(gè)反映系統(tǒng)發(fā)電機(jī)同調(diào)性的參數(shù)，根據(jù)參數(shù)的相似程度進(jìn)行分群；第二類是仿真擾動(dòng)后發(fā)電機(jī)功角的變化趨勢，通過辨識(shí)功角搖擺曲線相似程度進(jìn)行同調(diào)分群。其中第二種方法應(yīng)用越來越廣泛，很多學(xué)者采用時(shí)域分析、頻域分析、時(shí)頻分析來確定發(fā)電機(jī)功角搖擺曲線的相近程度。時(shí)域分析方法主要包括提取時(shí)域特征［6］、曲線相關(guān)分析［7］等。為了獲得更精確的反映發(fā)電機(jī)功角搖擺曲線趨勢的特征，一些學(xué)者利用頻域分析方法，采用傅里葉變換來提取頻譜分布特征［8］。時(shí)域分析能夠提取分辨率很高的時(shí)域特征，但是頻域分辨率卻很低，頻域分析能夠提取分辨率很高的頻域特征，但是時(shí)域分辨率卻很低，都不能全面準(zhǔn)確提取功角搖擺曲線的特征。為此，很多學(xué)者利用時(shí)頻分析來提取功角搖擺曲線特征，小波變換就是其中的典型代表。小波變換結(jié)合了時(shí)域分析和頻域分析的優(yōu)勢，在高頻部分具有很高的時(shí)間分辨率，在低頻部分具有很高的頻率分辨率，提取的特征能夠全面準(zhǔn)確地反映功角搖擺曲線特性［9］。文獻(xiàn)［10］采用小波分解方法將功角搖擺曲線分解為高頻信號(hào)和低頻信號(hào)，并且將低頻信息進(jìn)行重構(gòu)，將重構(gòu)的信息進(jìn)行相關(guān)性分析，這種方法既縮減了信息維數(shù)，又去除了信息噪聲，通過搖擺曲線基本升降規(guī)律識(shí)別同調(diào)機(jī)群，但是此方法計(jì)算量大。

同調(diào)機(jī)群的識(shí)別其本質(zhì)是一個(gè)聚類問題，即將同調(diào)的機(jī)組聚集在一起，不同調(diào)的機(jī)組分開。同調(diào)機(jī)組的聚類方法主要包括相關(guān)系數(shù)法［11］、層次分析法［12］、模糊聚類法［13］等。對(duì)于大系統(tǒng)而言，這些方法一般比較復(fù)雜，計(jì)算量大，不適合于大數(shù)據(jù)的處理。文獻(xiàn)［14］提出把Kohonen自組織神經(jīng)網(wǎng)絡(luò)運(yùn)用到同調(diào)識(shí)別，適合大數(shù)據(jù)的處理，但是需要預(yù)先指定競爭層網(wǎng)絡(luò)單元數(shù)目及其機(jī)構(gòu)形狀。

考慮到同調(diào)機(jī)群識(shí)別中存在的上述問題，首先使用小波變換的多尺度空間能量分布特征的提取方法提取機(jī)組功角搖擺曲線特征向量，然后將小波能量特征和時(shí)頻域特征組合，作為增長型自組織神經(jīng)網(wǎng)絡(luò)的輸入向量，在保證信號(hào)特征全面的情況下降低了特征的維度，通過控制閾值λ，最終在輸出層得到不同精度的分群結(jié)果。最后通過10機(jī)IEEE-39節(jié)點(diǎn)系統(tǒng)驗(yàn)證了此方法的計(jì)算量小，計(jì)算速度快，聚類結(jié)果具有更高的穩(wěn)定性和準(zhǔn)確性。

1 同調(diào)機(jī)組分群理論

發(fā)電機(jī)的同調(diào)性是指當(dāng)系統(tǒng)運(yùn)行方式改變或受擾后，系統(tǒng)中發(fā)電機(jī)功角變化呈現(xiàn)出相似的規(guī)律［15］。關(guān)于同調(diào)機(jī)群的定義多種多樣，最常用的定義如下，在計(jì)算時(shí)間t∈［0，T］內(nèi)滿足：

式中 Δδi=δi-δi0；Δδj=δj-δj0；δi、δj為發(fā)電機(jī)i、j的實(shí)時(shí)功角；δi0、δj0為發(fā)電機(jī)i、j的功角初值；Δδi、Δδj為發(fā)電機(jī)功角的實(shí)時(shí)變化量；ε是發(fā)電機(jī)之間同調(diào)程度的度量值。

但考慮到當(dāng)兩臺(tái)機(jī)組i，j擺幅都很小，并且一臺(tái)機(jī)組搖擺曲線上升，另一臺(tái)機(jī)組搖擺曲線下降的情況下也很有可能滿足式（1），而此時(shí)兩臺(tái)機(jī)組明顯不同調(diào)。因此利用相關(guān)性來定義機(jī)組同調(diào)：

k可以取任意常數(shù)使得Sij最小，令：

可求出：

綜合式（1）和式（4），定義機(jī)組同調(diào)識(shí)別判據(jù)為：

式中α為最大最小原則比重系數(shù)；β為均方根準(zhǔn)則比重系數(shù)（本文取α=0.6、β=0.4）。當(dāng)ε=0時(shí)，即發(fā)電機(jī)功角搖擺曲線走向完全一致，稱發(fā)電機(jī)完全同調(diào)；當(dāng)ε＞0，即發(fā)電機(jī)功角搖擺曲線最大差為ε，稱發(fā)電機(jī)ε-同調(diào)。在實(shí)際應(yīng)用中，沒有完全同調(diào)的發(fā)電機(jī)，所以一般使用ε-同調(diào)。

將滿足式（5）的發(fā)電機(jī)視為同調(diào)機(jī)組，不滿足式（5）的發(fā)電機(jī)視為不同調(diào)機(jī)組，這樣通過劃分外在功角搖擺曲線實(shí)質(zhì)劃分發(fā)電機(jī)內(nèi)在特性的過程就是機(jī)組同調(diào)識(shí)別。

2 小波能量特征提取

以短時(shí)傅里葉變換為過渡發(fā)展起來的小波變換，同時(shí)具有良好的時(shí)域趨勢描述、頻域局部特征刻畫功能，克服了傅里葉變換在時(shí)域上不能定位的缺點(diǎn)，并且小波的選取非常靈活，適合分析瞬時(shí)非平穩(wěn)信號(hào)，被廣泛應(yīng)用于特征提?。?6］。

對(duì)功角搖擺曲線在多尺度空間進(jìn)行小波分解，信號(hào)被分解為多個(gè)頻段，低頻段之前的系數(shù)為逼近系數(shù)，高頻段之前的系數(shù)為細(xì)節(jié)系數(shù)，合稱為小波系數(shù)，它們蘊(yùn)含著大量反映搖擺曲線特征的信息。通過小波系數(shù)計(jì)算各個(gè)頻段的小波能量，并且按頻次由大到小的順序排列，構(gòu)成小波能量特征向量，這就是小波變換多尺度空間能量分布特征的提取方法［17］。

假設(shè)f∈L2（R）是發(fā)電機(jī)的功角搖擺曲線，測量信號(hào)fj是f在尺度空間Vj的近似，則fj的多尺度分解可表示為：式中fj-k和wj-k分別為fj在空間Vj-k和Wj-k的逼近和細(xì)節(jié)。

為了計(jì)算快速方便，采用二進(jìn)小波變換中的Mallat算法實(shí)現(xiàn)［18］。利用Mallat算法對(duì)功角搖擺曲線逐層分解，每一層的逼近系數(shù)通過一組高通濾波器和低通濾波器得到下一層的逼近系數(shù)和細(xì)節(jié)系數(shù)［19］。功角搖擺曲線f（k），k=1，2，...，M，其二進(jìn)小波分解Mallat算法如下：

式中M是功角搖擺曲線f（k）的采樣點(diǎn)數(shù)；j是分解的層數(shù)；hn-2k，gn-2k分別是低通濾波系數(shù)和高通濾波系數(shù)；分別是功角搖擺曲線f（k）在分辨率2j-1下的逼近系數(shù)和細(xì)節(jié)系數(shù)；

對(duì)功角搖擺曲線f（k）進(jìn)行兩層小波分解，被分解為三個(gè)頻段，其中逼近系數(shù)為Ca2，細(xì)節(jié)系數(shù)為Cd1、Cd2。

提取小波變換多尺度空間能量特征的步驟如下：

首先給出小波能量定義，功角搖擺曲線f（k）在子空間Wj中小波能量定義為：

（1）確定信號(hào)各頻段小波系數(shù)的總能量，如下：

（2）構(gòu)造小波能量特征向量。依次計(jì)算各個(gè)頻段的小波能量值，并且按照頻段順序排列組合為特征向量：

利用上述流程可以對(duì)機(jī)組功角搖擺曲線信號(hào)進(jìn)行特征提取，提取出可以代表搖擺曲線的特征向量數(shù)據(jù)。即完成從搖擺曲線信號(hào)到特征向量數(shù)據(jù)的轉(zhuǎn)換。

3 自組織神經(jīng)網(wǎng)絡(luò)模型學(xué)習(xí)算法

增長型自組織神經(jīng)網(wǎng)絡(luò)不需要事先設(shè)定競爭層神經(jīng)元的個(gè)數(shù)，而是隨著閾值λ的改變，自動(dòng)生成不同數(shù)量的競爭層神經(jīng)元。該方法不受樣本輸入順序的影響，聚類結(jié)果準(zhǔn)確穩(wěn)定，被廣泛應(yīng)用于聚類分析［20］。

3.1 增長型自組織神經(jīng)網(wǎng)絡(luò)結(jié)構(gòu)

自組織神經(jīng)網(wǎng)絡(luò)和一般神經(jīng)網(wǎng)絡(luò)相比，沒有隱含層，只有由權(quán)值連接的輸入層和競爭層。如圖1所示。輸入層輸入樣本的特征向量，其中輸入層神經(jīng)元的個(gè)數(shù)等于樣本特征向量的維數(shù)。競爭層輸出樣本的分類結(jié)果，競爭層神經(jīng)元的個(gè)數(shù)等于最終的分類組數(shù)。競爭層神經(jīng)元個(gè)數(shù)是逐漸變化的，隨著輸入樣本的增加，樣本屬于不同類別的數(shù)量就會(huì)增加，競爭層神經(jīng)元的個(gè)數(shù)也就跟著增加。

圖1 自組織特征映射神經(jīng)網(wǎng)絡(luò)結(jié)構(gòu)圖Fig.1 Self-organizing neural network structure diagram

式中Xi為輸入向量；Wij為第j類組別模式的聚類中心；Dj為歐氏距離。

3.2 增長型自組織神經(jīng)網(wǎng)絡(luò)學(xué)習(xí)算法

增長型自組織神經(jīng)網(wǎng)絡(luò)屬于無導(dǎo)師學(xué)習(xí)型網(wǎng)絡(luò)，不需要先驗(yàn)知識(shí)的引導(dǎo)，自動(dòng)形成組群分類。該算法主要包括以下步驟：

（1）設(shè)置閾值λ。λ可以控制聚類結(jié)果的精細(xì)程度，聚類結(jié)果隨著λ的增大而變粗。在聚類前，通過試算λ值以應(yīng)對(duì)不同聚類情況；

（2）初始化權(quán)值。樣本輸入特征向量為：Xi=［xi1，xi2，...，xiN］T；N是樣本輸入特征向量維數(shù)；i為樣本數(shù)量。初始令競爭層神經(jīng)元個(gè)數(shù)為j=1，將任意樣本輸入特征向量作為輸入層與競爭層的連接權(quán)值Wij；

（3）競爭分類。輸入一個(gè)新的樣本，計(jì)算其與競爭層神經(jīng)元的歐氏距離Dj，歐式距離最小的競爭神經(jīng)元p獲勝：Dp=min｛Dj｝；

（4）更新權(quán)值。如果Dp＜λ，則將該樣本納入競爭神經(jīng)元k代表的組別模式，同時(shí)更新連接權(quán)值：Wij（n+1）=Wij（n）+η（n）（Xi-Wij（n）），η（n）是納入第j個(gè)動(dòng)態(tài)類型納入樣本數(shù)的倒數(shù)。然后轉(zhuǎn)入第（3）步；

（5）更新競爭層。如果Dp＞λ，則該樣本不納入任何一個(gè)競爭神經(jīng)元代表的組別模式，并且更新競爭層神經(jīng)元數(shù)量j=j+1，該樣本代表新的組別模式，

輸入層與競爭層中的連接權(quán)值起著中心指引的作用，相當(dāng)于整個(gè)自組織神經(jīng)網(wǎng)絡(luò)的分析系統(tǒng)，能夠記憶每個(gè)競爭神經(jīng)元代表的組別模式，隨著樣本劃分到不同的競爭神經(jīng)元，新增樣本的競爭神經(jīng)元的連接權(quán)值也要更新。當(dāng)一個(gè)樣本的特征向量與一個(gè)競爭神經(jīng)元的連接權(quán)值相似到一定程度時(shí)，那么該競爭神經(jīng)元就將其納入本組，否則，該樣本不屬于任何現(xiàn)有競爭神經(jīng)元，需要新增一個(gè)競爭神經(jīng)元，并用其連接權(quán)值記憶新的樣本模式［21］。

增長型自組織神經(jīng)網(wǎng)絡(luò)樣本輸入向量為N個(gè)功角搖擺曲線的特征，競爭層神經(jīng)元不斷更新的連接權(quán)值稱為該模式組別的聚類中心，根據(jù)輸入特征相似程度將樣本劃分到不同的組別模式中，使得功角搖擺曲線變化趨勢相似的發(fā)電機(jī)聚集在一起，不同的分開，輸入向量和組群中心的相似程度用距離來衡量，本文采用歐氏距離［22］：新組別模式與輸入層的連接權(quán)值Wij（n+1）為該樣本輸入特征向量；

（6）依次循環(huán)。當(dāng)全部樣本被納入競爭神經(jīng)元結(jié)束。競爭神經(jīng)元個(gè)數(shù)j就是聚類組數(shù)，每個(gè)競爭神經(jīng)元納入的樣本為具體分類結(jié)果。

上述增長型自組織神經(jīng)網(wǎng)絡(luò)學(xué)習(xí)的過程就是神經(jīng)網(wǎng)絡(luò)的構(gòu)建模型的過程，當(dāng)有新的實(shí)時(shí)樣本時(shí)，只需將新樣本輸入模型，在競爭層即可得到最終的分類結(jié)果，每個(gè)競爭神經(jīng)元代表不同的組別模式。

3.3 最優(yōu)閾值λ的確定

當(dāng)給定一個(gè)λ值時(shí)，系統(tǒng)機(jī)組被劃分為K類，即將系統(tǒng)中時(shí)序功角數(shù)據(jù)集 D=｛d1，d2，…，di，…，dN｝劃分為K個(gè)簇｛C1，…，Cj，…，CK｝，其中di∈Rm；N為系統(tǒng)發(fā)電機(jī)數(shù)量；m為時(shí)序功角向量維度。每個(gè)簇都有各自的聚類中心cj，計(jì)算所有簇內(nèi)各個(gè)數(shù)據(jù)點(diǎn)到聚類中心cj的歐式距離的平方和J（C），其函數(shù)值的大小反映了聚類結(jié)果的優(yōu)劣。

顯然，每個(gè)K值都對(duì)應(yīng)一個(gè)函數(shù)值J（C）。當(dāng)K值較小時(shí)，每個(gè)簇的歐氏距離很大，隨著K值的增大，函數(shù)J（C）減小的非常快。此時(shí)樣本快速分堆，具有相似特性的樣本聚集在一起，聚類效果迅速優(yōu)化；當(dāng)K值較大時(shí)，每個(gè)簇的歐式距離已經(jīng)很小，繼續(xù)增大K值，函數(shù)J（C）的變化將不再明顯，此時(shí)聚類結(jié)果已達(dá)到最優(yōu)。J-K關(guān)系曲線如圖2所示，圖中的拐點(diǎn)可確定最優(yōu)的分類K值，進(jìn)而根據(jù)該分類K值就可以得到最優(yōu)的閾值λ。

圖2 J-K關(guān)系曲線Fig.2 J-K relationship curve

4 算例分析

4.1 同調(diào)機(jī)群識(shí)別流程

選取發(fā)電機(jī)功角曲線的局部超調(diào)量d1ij和全局能量d2ij作為功角搖擺曲線時(shí)域特征［6］；選取發(fā)電機(jī)功角曲線的主頻成分ωmax作為功角搖擺曲線頻域特征；根據(jù)上述小波能量特征提取方法，構(gòu)建小波能量特征向量 F=［Ei0，Ei1，Ei2］。綜合時(shí)頻域特征、小波能量特征，構(gòu)成神經(jīng)網(wǎng)絡(luò)發(fā)電機(jī)樣本輸入向量為：G*i=［d1ij，d2ij，ωmax，Ei0，Ei1，Ei2］，并且將其歸一化：

以PSD-BPA為電力系統(tǒng)仿真平臺(tái)，模擬線路出現(xiàn)故障模式，對(duì)電力系統(tǒng)進(jìn)行故障模擬仿真。根據(jù)上文所訴內(nèi)容編寫計(jì)算機(jī)程序。程序流程如圖3所示。

圖3 程序流程圖Fig.3 Program flow chart

4.2 IEEE-39節(jié)點(diǎn)仿真分析算例

為驗(yàn)證本文方法的有效性，以IEEE-39節(jié)點(diǎn)系統(tǒng)為例進(jìn)行聚類分析，該系統(tǒng)包括39個(gè)節(jié)點(diǎn)，10臺(tái)發(fā)電機(jī)，其中31號(hào)發(fā)電機(jī)節(jié)點(diǎn)為平衡節(jié)點(diǎn)，30、32、33、34、35、36、37、38、39號(hào)發(fā)電機(jī)節(jié)點(diǎn)為 PV節(jié)點(diǎn)，負(fù)荷為恒定阻抗。IEEE-39節(jié)點(diǎn)系統(tǒng)如圖4所示。

圖4 IEEE 39節(jié)點(diǎn)系統(tǒng)接線圖Fig.4 Wiring diagram of IEEE 39-bus system

分別以發(fā)電機(jī)G30、G32為參考機(jī)，設(shè)BUS-3處發(fā)生三相斷路故障，故障持續(xù)時(shí)間0.12 s。切除故障前后網(wǎng)絡(luò)的拓?fù)浣Y(jié)構(gòu)保持不變。在BPA中得到每個(gè)發(fā)電機(jī)的功角變化數(shù)據(jù)，并經(jīng)過Matlab處理后分別如圖5、圖6所示。

圖5 IEEE 39發(fā)電機(jī)功角搖擺曲線Fig.5 IEEE 39 generators power angle of the rocking curve

圖6 IEEE 39發(fā)電機(jī)功角搖擺曲線Fig.6 IEEE 39 generators power angle of the rocking curve

樣本參數(shù)為S=｛G30，G31，G32，...，G39｝。首先只考慮時(shí)頻域特征作為特征向量，構(gòu)成10×3特征參數(shù)矩陣，將其輸入自組織神經(jīng)網(wǎng)絡(luò)進(jìn)行訓(xùn)練，得到的識(shí)別結(jié)果如表1所示；然后同時(shí)考慮小波能量特征和時(shí)頻域特征作為特征向量，構(gòu)成10×6特征參數(shù)矩陣，將其輸入增長型自組織神經(jīng)網(wǎng)絡(luò)進(jìn)行訓(xùn)練，得到的識(shí)別結(jié)果如表2所示。

表1 只考慮時(shí)頻域特征的同調(diào)機(jī)群劃分結(jié)果Tab.1 Coherent cluster division results when considering only the frequency domain characteristics of the coherent cluster division results

表2 同時(shí)考慮小波能量特征和時(shí)頻域特征的同調(diào)機(jī)群劃分結(jié)果Tab.2 Coherent fleet division results when taking into account the characteristics of wavelet energy and time frequency domain features coherent fleet division results

從表1可以看出，當(dāng)分類數(shù)為6時(shí)，只考慮時(shí)頻域特征的分類結(jié)果為：｛30｝｛31｝｛32｝｛37｝｛33，34，35，36｝｛38，39｝。從表2知，同時(shí)考慮小波特征和時(shí)頻域特征且分類數(shù)為6時(shí)，分類結(jié)果為：｛30，37｝｛31｝｛32｝｛33，34，35，36｝｛38｝｛39｝。按式（5）計(jì)算的同調(diào)判據(jù)值進(jìn)行分析，1.5 s之內(nèi)G38、G39同調(diào)判據(jù)值為0.075，G30與G37的同調(diào)判據(jù)值為0.053 5，G30與G37的同調(diào)性要強(qiáng)于G38與G39，從圖6也能直觀看出來，G30、G37的同調(diào)性明顯比G38、G39更強(qiáng)。因此，同時(shí)考慮小波特征和時(shí)頻域特征的分群結(jié)果更為合理。

當(dāng)分類數(shù)為4時(shí)，只考慮時(shí)頻域特征，分類結(jié)果為｛30，37｝｛31，32｝｛33，34，35，36｝｛38，39｝。同時(shí)考慮小波特征和時(shí)頻域特征的分類結(jié)果為｛30，37，38，39｝｛31｝｛32｝｛33，34，35，36｝。在1.5 s之內(nèi)G31、G32按式（5）計(jì)算同調(diào)判據(jù)值為0.270 6。而G30、G37、G38、G39相互之間按式（5）計(jì)算的同調(diào)判據(jù)最大值為0.148 0?？梢?，將｛30，37，38，39｝聚為一類是合理的。因此，同時(shí)考慮小波特征和時(shí)頻域特征的分群結(jié)果明顯要優(yōu)于只考慮時(shí)頻域特征的分群結(jié)果。

綜上所述，同時(shí)考慮小波能量特征和時(shí)頻域特征比只考慮時(shí)頻域特征作為增長型自組織神經(jīng)網(wǎng)絡(luò)輸入向量具有更高的準(zhǔn)確性，驗(yàn)證了所提出的方法更為有效。

5 結(jié)束語

提出了一種新的基于小波變換多尺度空間能量分布特征和增長型自組織神經(jīng)網(wǎng)絡(luò)的同調(diào)機(jī)組分群方法。改進(jìn)了同調(diào)機(jī)群識(shí)別的判據(jù)，利用小波變換多尺度空間能量分布特征的提取方法提取機(jī)組功角搖擺曲線信號(hào)的特征，并結(jié)合時(shí)頻域特征構(gòu)成輸入特征向量，減少了特征維數(shù)。利用融合自適應(yīng)共振理論和Kohonen網(wǎng)絡(luò)基本思想改進(jìn)了自組織神經(jīng)網(wǎng)絡(luò)，克服了自組織神經(jīng)網(wǎng)絡(luò)需要預(yù)先給定競爭層網(wǎng)絡(luò)單元數(shù)目及結(jié)構(gòu)形狀的問題。仿真算例可以看出，同時(shí)考慮小波能量特征和時(shí)頻域特征的聚類方法比只考慮時(shí)頻域特征具有更高的聚類準(zhǔn)確性。