王華慶 劉澤源 盧威 宋瀏陽(yáng) 韓長(zhǎng)坤

摘要：針對(duì)時(shí)域非平穩(wěn)振動(dòng)信號(hào)模式混疊、信噪比低，以及傳統(tǒng)稀疏表示算法模型復(fù)雜、優(yōu)化求解算法難以確定，導(dǎo)致故障特征提取難的問(wèn)題，提出了頻域組稀疏和群橋約束改進(jìn)迭代收縮閾值優(yōu)化的故障特征提取方法（Group Sparse Representation in Frequency Domain，GSRF）。將振動(dòng)信號(hào)轉(zhuǎn)換至頻域并對(duì)變量分組，構(gòu)造施加群橋約束的最小二乘回歸模型，準(zhǔn)確篩選沖擊相關(guān)變量；引入迭代重加權(quán)系數(shù)簡(jiǎn)化方程，以軟閾值收縮優(yōu)化求解頻域稀疏信號(hào)；對(duì)重構(gòu)的時(shí)域稀疏信號(hào)進(jìn)行包絡(luò)頻譜分析提取故障特征。試驗(yàn)結(jié)果表明，提出的頻域組稀疏算法優(yōu)于傳統(tǒng)的結(jié)合 L21范數(shù)約束的組稀疏索套方法，可有效提取微弱故障特征，實(shí)現(xiàn)稀疏域下的軸承故障診斷。

關(guān)鍵詞：故障診斷；滾動(dòng)軸承；微弱故障；稀疏表示；特征提取

中圖分類號(hào)： TH165+.3；TH133.33??? 文獻(xiàn)標(biāo)志碼： A??? 文章編號(hào)：1004-4523（2022）05-1242-08

DOI：10.16385/j .cnki .issn .1004-4523.2022.05.022

引言

滾動(dòng)軸承廣泛應(yīng)用在石化、軍工等領(lǐng)域，作為機(jī)械設(shè)備的重要傳動(dòng)零件，軸承在高速重載的工況下容易發(fā)生損傷失效，影響設(shè)備正常運(yùn)轉(zhuǎn)，因此對(duì)軸承早期故障的監(jiān)測(cè)和診斷至關(guān)重要[1]。分析振動(dòng)信號(hào)以提取故障特征是常用的診斷方法[2] ，但惡劣工況和復(fù)雜傳動(dòng)系統(tǒng)往往使測(cè)得的振動(dòng)信號(hào)呈現(xiàn)非高斯非平穩(wěn)、模式混疊、信噪比低的特性，故障特征被噪聲淹沒(méi)[3]。因此從早期微弱軸承故障信號(hào)中提取特征信息具有重要意義。

稀疏表示在故障診斷方向應(yīng)用廣泛，求解過(guò)程可以大致分為確定數(shù)學(xué)模型和優(yōu)化求解兩部分[4]。其數(shù)學(xué)模型是范數(shù)約束的最小二乘逼近函數(shù)，Hou 等[5]對(duì)希爾伯特解調(diào)信號(hào)施加 L1約束增強(qiáng)了分解結(jié)果的稀疏性，實(shí)現(xiàn)了未知轉(zhuǎn)速的軸承故障診斷。 Zheng 等[6]引入 L1、組2范數(shù)及無(wú)窮范數(shù)約束目標(biāo)函數(shù)，在保證重構(gòu)信號(hào)的稀疏性的同時(shí)提升了降噪效果，增強(qiáng)了早期微弱故障特征。優(yōu)化求解算法可以大致分為貪婪算法、優(yōu)化約束法、逼近算法、同倫算法。經(jīng)典的貪婪算法包含由Mallat等提出的匹配追蹤（Matching Pursuit，MP）及由 Pati 等改進(jìn)的正交匹配追蹤算法（Orthogonal? Matching? Pursuit， OMP），其重構(gòu)信號(hào)往往是近似解[7?8]。優(yōu)化約束法將原目標(biāo)函數(shù)（不光滑或不可導(dǎo)等不易快速求解的問(wèn)題）轉(zhuǎn)換為光滑的或可求導(dǎo)問(wèn)題，優(yōu)化后的目標(biāo)函數(shù)和原目標(biāo)函數(shù)不同，簡(jiǎn)化了求解過(guò)程。其代表算法為梯度下降和交替方向乘子法（Alternating Direc? tion of Method of Multipliers，ADMM）[9]。逼近算法將非光滑約束問(wèn)題拆解為多個(gè)簡(jiǎn)單子問(wèn)題，利用逼近算子分別求解這些子問(wèn)題以提高求解效率。經(jīng)典算法包括收縮算子、增廣拉格朗日算法、受控極小化（Majorization Minimization，MM）、迭代收縮閾值算法（Iterative? Shrinkage ? Thresholding? Algorithm ， ISTA）等[10?12]。

組稀疏回歸模型在求解中協(xié)同約束組內(nèi)、組間變量，促進(jìn)分解的稀疏性。An 等[13]提出多尺度周期組索套，利用多級(jí)小波子帶構(gòu)建稀疏模型，通過(guò) MM 優(yōu)化非凸罰函數(shù)，在小波域提取周期脈沖，實(shí)現(xiàn)強(qiáng)噪干擾下的特征提取。Zhao 等[14]構(gòu)建稀疏多周期組套索模型，使用 MM 解耦多重故障，實(shí)現(xiàn)了低信噪比復(fù)合故障信號(hào)的診斷。Zhao 等[15]基于增強(qiáng)稀疏組 Lasso 懲罰，自適應(yīng)選擇確定參數(shù)，并在 MM 優(yōu)化中嵌入周期先驗(yàn)，直接提取時(shí)域脈沖，實(shí)現(xiàn)了微弱故障信號(hào)的診斷。

雖然上述方法成功提取了故障特征，但仍存在如下問(wèn)題。1）模型過(guò)于復(fù)雜，多數(shù)模型為提升分解結(jié)果會(huì)引入多項(xiàng)罰參量，容易降低組內(nèi)、組間變量選擇的一致性。2）優(yōu)化保凸能力弱，復(fù)雜非凸模型求解對(duì)優(yōu)化算法要求高，當(dāng)懲罰項(xiàng)不易平滑時(shí)將顯著增加計(jì)算復(fù)雜度。

為解決上述問(wèn)題，提出了基于頻域的組稀疏特征提取方法，以組稀疏保證分解稀疏性的同時(shí)保留脈沖相關(guān)信息。軸承故障信號(hào)是故障特征信號(hào)與噪聲的耦合信號(hào)，時(shí)域降噪難以實(shí)現(xiàn)解耦，將信號(hào)變換到傅里葉正交空間，有利于實(shí)現(xiàn)信號(hào)解耦，降低噪聲，提取沖擊信號(hào)。因此，先建立沖擊敏感的數(shù)學(xué)模型并轉(zhuǎn)換至頻域的振動(dòng)信號(hào)按組劃分，再對(duì)分組后的頻域信號(hào)施加群橋約束以選取沖擊相關(guān)特征變量，引入迭代重加權(quán)系數(shù)，簡(jiǎn)化目標(biāo)函數(shù)方程便于 ISTA 優(yōu)化求解，最后重構(gòu)稀疏時(shí)域信號(hào)，在包絡(luò)譜上提取出信號(hào)故障特征。 GSRF 對(duì)同一懲罰項(xiàng)施加多種約束來(lái)保證變量篩選的一致性和稀疏性，利用權(quán)重系數(shù)凸優(yōu)化回歸模型，降低運(yùn)算難度，通過(guò)對(duì)軸承早期故障信號(hào)的分析驗(yàn)證了所提方法的有效性。

1 基本原理

1.1 稀疏表示

用于生物信息、數(shù)理統(tǒng)計(jì)、信號(hào)處理等領(lǐng)域的數(shù)據(jù)維數(shù)通常很大，對(duì)信號(hào)做稀疏表示有利于緩解數(shù)據(jù)的存儲(chǔ)、傳輸、提取的壓力。稀疏表示簡(jiǎn)單理解為原信號(hào)是用少量基原子線性組合而成的[16]，即：

式中振動(dòng)信號(hào) Y ∈ Rn，觀測(cè)矩陣 D∈ Rn × m 可為字典集或變換矩陣，稀疏系數(shù) a ∈ Rm，噪聲ε∈ Rn。

常用的稀疏表示模型可以總結(jié)為：

式中 p >0，懲罰系數(shù)λ>0。式（2）表示常用的范數(shù)約束求解，當(dāng) p=0時(shí)，是多項(xiàng)式復(fù)雜程度的非確定性問(wèn)題，難以直接求解，Donoho與Candes等證明 p=1時(shí)，式（2）轉(zhuǎn)變?yōu)樗魈姿惴ǎ↙east absolute shrinkage and selection operator，Lasso ），在滿足一定稀疏性條件下，與p=0有同樣的稀疏解。

1.2 橋組稀疏

傳統(tǒng) Lasso 算法僅能解決組內(nèi)變量選擇問(wèn)題，稀疏性不強(qiáng)，Yuan 等在 Lasso 基礎(chǔ)上引入組間變量選擇，Huang 等又加入了橋懲罰項(xiàng)，進(jìn)一步提高了運(yùn)算速度與稀疏性，廣泛應(yīng)用在數(shù)據(jù)分析統(tǒng)計(jì)中[17?18]。

其原理為，對(duì) n 維1列的信號(hào) y =（y1，y2，…，yn ）′，有 n 維 d 列設(shè)計(jì)向量 Dk =（ D 1k，D 2k，…，Dnk ）′（k=1，…，d），A1，…，AJ 表示設(shè)計(jì)向量 D 的已知分組（1≤ J ≤ d），以aAj =（ak，k ∈Aj ）′表示第 j 組的回歸系數(shù)，有橋組稀疏目標(biāo)函數(shù)為：

式中懲罰系數(shù)λ>0，橋懲罰系數(shù)γ>0。式（3）是非凸問(wèn)題，需轉(zhuǎn)化為便于求解的凸函數(shù)。

1.3 迭代重加權(quán)

任意向量 e ∈ Rn，其p 范數(shù)可定義為：

對(duì)式（4）變形：

利用式（6）及權(quán)重 w 即可以將式（3）中的非凸問(wèn)題轉(zhuǎn)變?yōu)橥箖?yōu)化求解[19]。

1.4 近端算子

有函數(shù)f（x ）=g（ x ）+h（x ），其中 g（ x ）是可導(dǎo)凸函數(shù)，h（x ）不可導(dǎo)，其最小值可由迭代循環(huán)求解：

式中? t 為步長(zhǎng)，常取李普希茲常數(shù)的倒數(shù)；?g （ x ）為梯度。式（7）可寫為不動(dòng)點(diǎn)迭代形式：

式中prox稱為近端算子，在廣義梯度下降更新求解中廣泛應(yīng)用，下文將結(jié)合迭代重加權(quán)的思想進(jìn)一步優(yōu)化約束。綜上所述，確定了數(shù)學(xué)回歸模型和優(yōu)化求解算法，便可通過(guò)迭代運(yùn)算求近似解。

2 基于頻域組稀疏的特征提取

2.1 組稀疏

傳統(tǒng)稀疏約束對(duì)向量?jī)?nèi)全部變量進(jìn)行篩選，沒(méi)有考慮變量間關(guān)系，組稀疏將 n 維向量 x 內(nèi)變量分組，示意圖如圖1所示。

圖1中，xgi表示第i組變量，每組變量數(shù)量可以根據(jù)信號(hào)特點(diǎn)劃分，若不同組間無(wú)重復(fù)的變量可認(rèn)為是非重疊組，反之認(rèn)為是重疊組劃分，即當(dāng)i≠j 時(shí)，非重疊組中xgi與xgj無(wú)相同元素，重疊組中xgi與xgj有相同元素。重疊組稀疏不僅約束組內(nèi)變量，還影響組間變量的選擇，顯著增強(qiáng)每次迭代對(duì)沖擊相關(guān)變量的篩選能力，使重構(gòu)時(shí)域信號(hào)具有更強(qiáng)的稀疏性[20]，沖擊特征更明顯。

為解決時(shí)域非平穩(wěn)振動(dòng)信號(hào)模式混疊、信噪比低及分解結(jié)果不稀疏的問(wèn)題，將待處理信號(hào)轉(zhuǎn)換到頻域。引入重疊組稀疏后，數(shù)學(xué)回歸模型寫為：

式中 x 為待求解的稀疏信號(hào)，y 為轉(zhuǎn)換至頻域的原始振動(dòng)信號(hào)，懲罰系數(shù)λ>0。2范數(shù)平方項(xiàng)保證重構(gòu)信號(hào)和原信號(hào)的一致性，對(duì)組變量xg施加1范數(shù)約束以增強(qiáng)沖擊特征選取能力[21]。

2.2 組橋約束

為進(jìn)一步提高分解結(jié)果的信噪比，引入群橋約束組的變量[22?23]，式（9）寫為：

式（3）中橋懲罰系數(shù)取0<γ<1時(shí)可以強(qiáng)化沖擊相關(guān)變量的選擇能力，此處依據(jù)經(jīng)驗(yàn)選取群橋系數(shù)γ=0.5，但式（10）是非凸優(yōu)化問(wèn)題，結(jié)合1.3節(jié)介紹的迭代重加權(quán)，式（10）轉(zhuǎn)化為：

上式形似式（2）的基礎(chǔ) Lasso 問(wèn)題。據(jù)1.4節(jié)介紹，式（12）的前一項(xiàng)可視為 g（ x ），后一項(xiàng)可視為 h（x ），對(duì) g（ x ）二階泰勒展開(kāi)并配方后，目標(biāo)函數(shù)式（12）存在近端算子：

式（13）有數(shù)值解，更新可得稀疏信號(hào) x。

2.3 基于頻域橋組稀疏的故障診斷方法

針對(duì)軸承的早期微弱故障，提出了以重疊組、群橋約束為基礎(chǔ)的稀疏模型，以迭代重加權(quán)和迭代收縮閾值為優(yōu)化算法的微弱故障特征提取方法，其具體實(shí)現(xiàn)步驟如圖2所示流程圖。

所提方法為減弱非平穩(wěn)時(shí)域噪聲干擾在頻域進(jìn)行特征提取，主要步驟如下：

（1）對(duì)轉(zhuǎn)換到頻域的振動(dòng)信號(hào)內(nèi)的變量按照重疊組劃分，確定最小二乘結(jié)合 L1范數(shù)約束組變量的基礎(chǔ)數(shù)學(xué)模型。

（2）施加群橋約束，進(jìn)一步強(qiáng)化篩選沖擊敏感的組間及組內(nèi)變量，并引入迭代重加權(quán)系數(shù)簡(jiǎn)化目標(biāo)函數(shù)，將非凸問(wèn)題轉(zhuǎn)化為形如加權(quán) Lasso 的凸函數(shù)。

（3）對(duì)優(yōu)化后的數(shù)學(xué)模型采用 ISTA 迭代更新求解頻域稀疏信號(hào)，再轉(zhuǎn)換回時(shí)域，最后對(duì)時(shí)域稀疏信號(hào)進(jìn)行包絡(luò)譜分析，提取故障特征頻率，實(shí)現(xiàn)故障診斷。

3 實(shí)驗(yàn)分析與驗(yàn)證

為驗(yàn)證 GSRF 的有效性，分別采用公開(kāi)軸承故障數(shù)據(jù)及課題組試驗(yàn)臺(tái)數(shù)據(jù)，與傳統(tǒng)非重疊組結(jié)合 L21范數(shù)約束的 Lasso 算法做了對(duì)比分析研究。

3.1 實(shí)驗(yàn)驗(yàn)證一

采用凱斯西儲(chǔ)大學(xué)公開(kāi)數(shù)據(jù)集中的 SKF 深溝球軸承，故障類型為內(nèi)圈及外圈故障，采樣頻率12×103 Hz，內(nèi)圈故障直徑為0.014 in（0.356 mm），電機(jī)轉(zhuǎn)速為1772 r/min，計(jì)算內(nèi)圈故障特征頻率約為159.9 Hz；外圈故障直徑0.021 in（0.533 mm），電機(jī)轉(zhuǎn)速1750 r/min，外圈故障頻率約為104.55 Hz。

如圖3所示，內(nèi)圈故障特征頻率完全被噪聲淹沒(méi)，分別以譜峭度、非重疊組結(jié)合 L21范數(shù)約束的 Lasso （傳統(tǒng)方法）和頻域群橋組稀疏（GSRF）做分析，圖4為時(shí)域重構(gòu)信號(hào)對(duì)比圖，圖5為包絡(luò)譜對(duì)比。

所提方法綜合考慮軸承轉(zhuǎn)速、采樣頻率及頻域分辨率的影響，選取頻域信號(hào)的重疊組變量系數(shù)為20。圖4中（ a ）為譜峭度，（b）為傳統(tǒng)方法，（ c ）為 GSRF 。由圖3和4對(duì)比可知，譜峭度和傳統(tǒng)方法恢復(fù)的時(shí)域信號(hào)更稀疏，GSRF 處理后的時(shí)域信號(hào)與原始振動(dòng)信號(hào)相比無(wú)明顯變化。圖 5（ a ）為譜峭度稀疏信號(hào)包絡(luò)譜，（b）為傳統(tǒng)方法稀疏信號(hào)包絡(luò)譜， GSRF 包絡(luò)譜如圖（ c ）所示。對(duì)比包絡(luò)譜，譜峭度和傳統(tǒng)方法實(shí)現(xiàn)了故障特征增強(qiáng)，但對(duì)噪聲成分的抑制效果較弱，GSRF 在增強(qiáng)故障特征的同時(shí)，較好地抑制了故障無(wú)關(guān)頻率成分。

沖擊明顯，包絡(luò)譜上特征頻率受噪聲干擾。圖 7，8中（ a ），（b），（ c ）分別是譜峭度、傳統(tǒng)方法及 GSRF 的分析結(jié)果。對(duì)比可知，傳統(tǒng)方法處理后的信號(hào)在時(shí)域上更稀疏，GSRF 則在頻域上消除了更多噪聲干擾。綜合分析認(rèn)為 GSRF 在稀疏分解中盡可能多的保留了沖擊信息，有效消除了故障無(wú)關(guān)成分的

圖9中，計(jì)算復(fù)雜的 GSRF 收斂相對(duì)更慢，也能在10步內(nèi)收斂。表 1記錄計(jì)算時(shí)間，分析認(rèn)為 GSRF 對(duì)變量的細(xì)致篩選會(huì)降低運(yùn)算速度，但卻能顯著提升診斷結(jié)果的準(zhǔn)確性。

3.2? 實(shí)驗(yàn)驗(yàn)證二

在如圖10所示的試驗(yàn)臺(tái)上采集振動(dòng)信號(hào)，對(duì) N205EM 圓柱滾子軸承做水平和豎直方向的數(shù)據(jù)采樣，采樣頻率為105 Hz，電機(jī)轉(zhuǎn)速為1300 r/min，線切割加工內(nèi)、外圈故障寬度均為0.3 mm，深度為0.05 mm，計(jì)算內(nèi)、外圈故障特征頻率分別為145.84，86.32 Hz。

內(nèi)圈故障信號(hào)時(shí)域頻譜如圖11所示，故障特征頻率已被噪聲淹沒(méi)。綜合考慮軸承轉(zhuǎn)速、采樣頻率及頻域分辨率的影響，對(duì)頻域信號(hào)按每組20個(gè)變量劃分重疊組，處理結(jié)果如圖12所示。譜峭度和傳統(tǒng)方法的包絡(luò)譜分析在圖12（ a ），（b）中故障特征頻率較原信號(hào)有一定增強(qiáng)，但仍無(wú)法消除其他噪聲成分的干擾，GSRF 處理的結(jié)果在圖12（ c ）中則有效抑制了非平穩(wěn)噪聲。

圖13為外圈故障信號(hào)的波形頻譜圖，故障特征頻率86.32 Hz 明顯被淹沒(méi)。圖 14（ a ），（b）是分別應(yīng)用譜峭度和傳統(tǒng)方法處理后的結(jié)果，沖擊特征頻率及其高次諧波明顯被增強(qiáng)，但無(wú)法消除旁瓣噪聲的影響，在圖14（ c ）中故障特征頻率及其高次諧波明顯且消除了干擾。表2為傳統(tǒng)方法與 GSRF 計(jì)算時(shí)間對(duì)比結(jié)果，GSRF 在沒(méi)有顯著增加運(yùn)算時(shí)間的情況下獲得了更好的故障識(shí)別效果，實(shí)現(xiàn)了微弱故障的診斷。

圖15和16討論了域變換和重疊組因素對(duì)分解結(jié)果的影響。圖 15和16中：（ a ）是時(shí)域重疊組橋稀疏算法分解結(jié)果；（b）是頻域無(wú)重疊組的橋稀疏算法結(jié)果，（ c ）是 GSRF 算法處理結(jié)果。

圖15（ a ），16（ a ）顯示，時(shí)域運(yùn)算分解結(jié)果稀疏性較強(qiáng)，但故障信息丟失較多，對(duì)比圖3，15（b）和16（b），無(wú)重疊組的橋稀疏算法增強(qiáng)了沖擊特征，但無(wú)法實(shí)現(xiàn)信號(hào)與噪聲的解耦，不能提取故障特征頻率，而 GSRF 在增強(qiáng)故障特征的同時(shí)，抑制了故障無(wú)關(guān)頻率成分，有效提取了故障特征頻率。

4 結(jié)論

確定以最小二乘結(jié)合 L1范數(shù)約束的基礎(chǔ)數(shù)學(xué)回歸模型，根據(jù)故障信號(hào)沖擊特點(diǎn)引入組稀疏及群橋約束，提出了基于頻域組稀疏結(jié)合群橋約束的故障特征提取方法。實(shí)驗(yàn)分析結(jié)果表明：GSRF 將低信噪比信號(hào)轉(zhuǎn)換至頻域進(jìn)行稀疏表示，可以減弱時(shí)域非平穩(wěn)噪聲擾動(dòng)帶來(lái)的模式混疊影響，提高了特征提取精度。組稀疏和群橋正則化按頻率差別篩選組內(nèi)、組間變量，在完整保留故障沖擊信息的同時(shí)濾除了無(wú)關(guān)頻率成分，進(jìn)一步增強(qiáng)了故障特征。通過(guò)引入迭代重加權(quán)系數(shù)優(yōu)化非凸目標(biāo)函數(shù)，顯著降低運(yùn)算復(fù)雜度，便于使用 ISTA 求解。與傳統(tǒng)非重疊組結(jié)合 L21范數(shù)約束的 Lasso 方法對(duì)比，驗(yàn)證了所提方法具有較強(qiáng)的特征提取能力和噪聲魯棒性。

參考文獻(xiàn)：

[1] 陳是扦，彭志科，周鵬.信號(hào)分解及其在機(jī)械故障診斷中的應(yīng)用研究綜述[ J ].機(jī)械工程學(xué)報(bào)，2020，56（17）：91-107.

Chen Shiqian，Peng Zhike，Zhou Peng . Review of sig? nal decomposition theory and its applications in machine fault diagnosis [ J ]. Journal of Mechanical Engineering，2020，56（17）：91-107.

[2] 魏永合，聶晨，李宏林. LMD 與優(yōu)化 OMP 算法的滾動(dòng)軸承故障診斷方法研究[ J ].沈陽(yáng)理工大學(xué)學(xué)報(bào)，2020，39（3）：61-66.

Wei Yonghe，Nie Chen，Li Honglin . Study on fault di? agnosis method of rolling bearing based on LMD and op? timized? OMP? algorithm [ J ]. Journal? of? Shenyang? Li? gong University，2020，39（3）：61-66.

[3] 李華，劉韜，伍星，等. EEMD 和優(yōu)化的頻帶熵應(yīng)用于軸承故障特征提取[ J ].振動(dòng)工程學(xué)報(bào)，2020，33（2）：414-423.

Li? Hua ，Liu? Tao ，Wu? Xing ，et? al . EEMD? andopti? mized frequency band entropy for fault feature extraction on? bearings [ J ]. Journal? of? Vibration? Engineering，2020，33（2）：414-423.

[4] 李清泉，王歡.基于稀疏表示理論的優(yōu)化算法綜述[ J ].測(cè)繪地理信息，2019，44（4）：1-9.

Li Qingquan，Wang Huan . Sparse representation-based optimization：a? survey [ J ]. Journal? of? Geomatics，2019，44（4）：1-9.

[5]? Hou F，Selesnick I，Chen J，et al . Fault diagnosis forrolling bearings under unknown time-varying speed con? ditions with sparse representation[ J ]. Journal of Sound and Vibration，2021，494：115854.

[6]? Zheng? K ，Bai? Y ，Xiong? J ，et? al . Simultaneously? lowrank and group sparse decomposition for rolling bearing fault diagnosis[ J ]. Sensors，2020，20（19）：1-27.

[7]? Liang K，Zhao M，Lin J，et al . An? information-basedK-singular-value? decomposition? method? for? rolling? ele? ment? bearing? diagnosis [ J ]. ISA? Transactions ，2020，96：444-456.

[8]? Zhang X，Liu Z，Wang L，et al . Bearing fault diagnosisbased? on? sparse? representations? using? an? improved OMP? with? adaptive? Gabor? sub-dictionaries [ J ]. ISA Transactions，2020，106：355-366.

[9]? Wang B，Liao Y，Ding C，et al . Periodical sparse low-rank? matrix? estimation? algorithm? for? fault? detection? of rolling? bearings [ J ]. ISA? Transactions ， 2020， 101：366-378.

[10] Zheng? Kai ，Yang? Dewei，Zhang? Bin ，et? al . A? groupsparse representation method in? frequency domain with adaptive? parameters? optimization? of&nbsp; detecting? incipient rolling? bearing? fault [ J ]. Journal? of? Sound? and? Vibra ? tion，2019，462：114931.

[11] Zhang W，Yu D，Yan X，et al . Weak multiple fault de?tection? based? on? weighted? morlet? wavelet-overlapping group sparse for rolling bearing fault diagnosis[ J ]. Ap? plied Sciences，2020，10（6）：2057.

[12] DiwuZhenkun，Cao Hongrui，Wang Lei，et al . Collab ?orative? double? sparse? period-group? Lasso? for? bearing fault diagnosis[ J ]. IEEE? Transactions on Instrumenta ? tion and Measurement，2021，70：1-10.

[13] An B，Zhao Z，Wang S，et al . Sparsity-assisted bearingfault diagnosis using multiscale period group Lasso[ J ]. ISA Transactions，2019，98：338-348.

[14] Zhao? Z ，Wang? S ，Sun? C ，et? al . Sparse? multiperiodgroup? Lasso? for bearing multifault diagnosis[ J ]. IEEE Transactions? on? Instrumentation? and? Measurement，2019，69（2）：419-431.

[15] Zhao Z，Wu S，Qiao B，et al . Enhanced sparse period-group Lasso for bearing fault diagnosis[ J ]. IEEE Trans ?actions on Industrial Electronics ，2019，66（3）：2143-2153.

[16]王華慶，任幫月，宋瀏陽(yáng)，等.基于終止準(zhǔn)則改進(jìn) K-SVD 字典學(xué)習(xí)的稀疏表示特征增強(qiáng)方法[ J ].機(jī)械工程學(xué)報(bào)，2019，55（7）：35-43.

Wang? Huaqing，Ren? Bangyue，Song? Liuyang，et? al .Sparse representation method based on termination crite? ria improved K-SVD dictionary learning? for feature en? hancement [ J ]. Journal? of? Mechanical? Engineering，2019，55（7）：35-43.

[17] Yuan M，Lin Y . Model selection? and? estimation in re ?gression with grouped variables[ J ]. Journal of the Royal Statistical Society：Series B，2006，68（1）：49-67.

[18] Huang J，Ma S，Xie H，et al . A group bridge approachfor? variable? selection [ J ]. Biometrika，2009，96（2）：339-355.

[19] Wang? X ，Wang? M . Adaptive? group? bridge? estimationfor high-dimensional partially linear models[ J ]. Journal of Inequalities and Applications，2017，2017（1）：1-18.

[20] Li? G ，Liu? X ，Chen? K . Integrative? multi-view? regres?sion：bridging? group-sparse? and? low-rank? models [ J ]. Biometrics，2019，75（2）：593-602.

[21] Park? C ，Yoon? Y? J . Bridge? regression ：adaptivity? andgroup selection[ J ]. Journal of Statistical Planning & In? ference，2011，141（11）：3506-3519.

[22] Burrus C? S，Barreto J A，Selesnick I W . Iterative re ?weighted? least-squares? design? of FIR? filters [ J ]. IEEE Transaction on Signal Processing，1994，42（11）：2926-2936.

[23]劉建偉，崔立鵬，羅雄麟.組稀疏模型及其算法綜述[ J ].電子學(xué)報(bào)，2015，43（4）：776-782.

Liu? Jianwei， Cui? Lipeng， Luo? Xionglin . Survey? on group sparse models and algorithms[ J ]. Acta Electroni? ca Sinica，2015，43（4）：776-782.

Fault feature extraction using group sparse representation in frequency domain

WANG Hua-qing1，LIU Ze-yuan1，LU Wei2，SONG Liu-yang1，HAN Chang-kun1????????????????? （1.College of Mechanical and Electrical Engineering，Beijing University of Chemical Technology，Beijing 100029，China；

2.Institute of Engineering Technology，Sinopec Catalyst Company Limited，Beijing 101100，China）

Abstract： A fault feature extraction method based on group sparsity and improved iterative shrinkage threshold optimization in fre ? quency domain （GSRF） is proposed to solve the issues in rolling bearing diagnosis about difficulty in mathematical model determi? nation，sparse constraints and optimization algorithm selection . The vibration signals are converted into the frequency domain and the variables? are divided by? overlapping rules . The? least square regression model with? group bridge? constraint is? constructed to screen impact related variables accurately . The iterative reweighting coefficient is introduced to simplify the equation，so that the sparse signal in frequency domain can be solved by iterative shrinkage-thresholding algorithm . The envelope spectrum analysis of reconstructed sparse signal in time domain is carried out to extract the fault features . The experimental results show that the pro? posed algorithm is superior to the traditional group sparse LASSO combined with L21 norm constraint . GSRF can effectively ex? tract weak fault features and achieve bearing fault diagnosis in the sparse domain .

Key words ： fault diagnosis；rolling bearing；weak fault；sparse representation；feature extraction

作者簡(jiǎn)介：王華慶（1973—），男，教授。電話：13801023830；E-mail：hqwang@mail .buct .edu .cn。

通訊作者：宋瀏陽(yáng)（1988—），女，副教授。E-mail：xq _0703@163.com。