嚴 旭 劉若晨 張?zhí)m春 孫見忠 周子元
(①江蘇理工學院汽車與交通工程學院,江蘇 常州 213001;②南京航空航天大學民航學院,南京 211106;③南京英田光學工程股份有限公司,江蘇 南京 210046)
滾動軸承在各種旋轉(zhuǎn)機械中的應用最為廣泛,也是整個設備中極易出現(xiàn)故障的通用部件,因此對滾動軸承進行故障診斷是旋轉(zhuǎn)機械設備中不可或缺的環(huán)節(jié)[1-2]。靜電監(jiān)測技術(shù)利用磨損區(qū)域荷電機理,感應探測缺陷磨損處的荷電水平,且不受運動部件的工況干擾[3]。相較于傳統(tǒng)振動監(jiān)測,靜電監(jiān)測避免了直接接觸耦合其他振動源,能更直觀地顯示故障位置的磨損情況。靜電監(jiān)測最早應用于航空航天發(fā)動機的燒蝕、碰磨等氣路故障監(jiān)測,隨后Harvey T J 等[4]驗證了靜電技術(shù)在磨損區(qū)域的可行性。張營[5]通過銷盤磨損試驗臺對軸承鋼初期故障的靜電監(jiān)測展開研究。滾動軸承故障診斷信號處理方法常見的有小波包分解算法[6]和EEMD 算法[7],但小波包分解算法受分解程度影響較大,而EEMD 算法的抗干擾能力受白噪聲幅值系數(shù)和總體平均次數(shù)等影響較大。
稀疏表示(sparse representation,SR)[8]是一種新型信號處理方法,利用稀疏分解將原始信號用少量原子的線性組合表示,主要分為字典構(gòu)造方法和稀疏恢復算法兩部分。對比其他信號處理方法,SR 能夠從原信號中提取并重構(gòu)故障特征成分,從而有效地將故障信號從噪聲中提取出來。近幾年,稀疏表示理論引起了學者們的關(guān)注,鄭勝等[9]提出了將級聯(lián)過完備字典與特征符號搜索算法相結(jié)合的軸承復合故障診斷方法。Sun R B 等[10]采用結(jié)構(gòu)化稀疏時頻分析方法進行齒輪故障診斷。鞏曉赟等[11]基于OMP 的稀疏表示方法對軸承耦合故障振動特性進行特征提取。He G L 等[12]完成了基于SR 的滾動軸承故障特征提取。
針對軸承靜電信號是一種微弱非平穩(wěn)信號,容易受到噪聲干擾,有效靜電信號易被噪聲淹沒。本文選取Laplace 原子構(gòu)建解析字典,采用聚類收縮分段正交匹配追蹤算法[13]求解稀疏表示系數(shù),通過StOMP 和CcStOMP 兩種算法的軸承實驗靜電信號故障特征提取效果對比,驗證CcStOMP 算法在滾動軸承靜電監(jiān)測信號的故障特征提取性能。
靜電監(jiān)測技術(shù)是針對滾動軸承磨損區(qū)域的帶電顆粒進行的一種狀態(tài)監(jiān)測技術(shù)。當帶電粒子通過靜電傳感器探極感應面的摩擦副間時,電場線將終止于探極感應面。由于靜電感應作用,摩擦副表面的靜電電荷源被吸引出相斥極性的電荷,探極內(nèi)大量自由移動的電子將聚集或離開探極感應面,同時正電子向相反方向運動形成感應電流,從而形成靜電感應。
在靜電感應過程中,傳感器末端與信號調(diào)理單元連接,將感應電荷量實時轉(zhuǎn)化為可供測量研究的靜電感應電壓信號,從而實現(xiàn)磨損區(qū)域的靜電監(jiān)測。當機械零件磨損加劇時,摩擦副周圍的帶電磨粒電荷量會大幅增加。因此,可通過靜電傳感器實時監(jiān)測機械零件的磨損情況。靜電監(jiān)測原理如圖1 所示。
圖1 靜電監(jiān)測原理
滾動軸承磨損區(qū)域靜電監(jiān)測過程中,記靜電傳感器探極面感應的電荷量為Q(t),采集系統(tǒng)顯示的電壓信號記為U(t),默認系統(tǒng)初始狀態(tài)為0。電荷量Q(t)與電壓量U(t)經(jīng)過拉普拉斯變換分別為Q(s)、U(s),兩者滿足的關(guān)系式為
式中:R為探測回路的等效電阻;s為t時刻對應復頻率;C為回路等效電容。
此外,油氣勘查開采監(jiān)管與服務也不斷加強。按照“雙隨機一公開”原則,在2017年5月隨機抽取了87個油氣勘查開采項目,在2017年底之前由各省(區(qū)、市)完成實地核查任務。根據(jù)核查結(jié)果,按照相關(guān)規(guī)定,將其中32個勘查項目列入異常名單,1個侵權(quán)開采項目按民事侵權(quán)糾紛案處理。
當R·C·s?1時,靜電信號監(jiān)測模型可等效為
因此,軸承采集系統(tǒng)采集的原始電壓信號U(t)和靜電傳感器感應面電荷量Q(t)的一階導數(shù)成正比,表明了靜電水平能直觀反映滾動軸承的故障程度。
稀疏表示的數(shù)學性質(zhì)是通過最佳原子加權(quán)和實現(xiàn)最優(yōu)化地表示信號,因此稀疏表示主要分為構(gòu)造字典和稀疏恢復算法兩部分。假設原始靜電信號為向量組y∈πA,其向量展開式為y=[y1,y2,y3,···,yA],若向量組存在絕大多數(shù)元素滿足yi=0,則向量y可為稀疏向量。若存在一個轉(zhuǎn)變矩陣B∈πA×K,其向量展開式B=[b1,b2,b3,···,bK],滿足 π矩陣列因子大于行因子,且 πA為滿秩矩陣,則可將原始靜電信號定義為欠定線性方程:
式中:d=[d1,d2,d3,···,dK]T∈πK為稀疏向量;γ為殘差信號。
1966年3月,父親與十幾萬南京軍區(qū)的官兵復員來到新疆兵團支邊。一來就分配在農(nóng)一師運輸公司汽車二連,單位駐地在阿拉爾,父親他們的任務就是拉運各種開發(fā)建設塔里木的物資,常年累月奔波往返于烏魯木齊、葉魯番、大河沿和阿拉爾這條運輸線上。那時都是簡易土公路,路況極差,從烏魯木齊往返阿拉爾一趟需要十一、二天的時間,車在路上顛簸,人在車里顛簸。出一趟車,不知遭多少罪。可父親從不抱怨,他們那代人心里裝著祖國,裝著大家,再苦再累也無怨無悔。
式中:δ為小正數(shù),用于衡量噪聲水平。
式中:‖d‖q為lq的范數(shù)。
實驗參數(shù)設置如下:軸承轉(zhuǎn)速為3 800 r/min,采樣頻率為12 kHz,采集時間1 s,根據(jù)軸承故障理論計算得外圈故障頻率和滾子故障頻率分別為193.17 Hz 和126.03 Hz。用圖4 搭建的滾動軸承靜電監(jiān)測實驗平臺分別采集軸承外圈和滾子重度的故障靜電信號,滾動軸承主要參數(shù)見表1。
因稀疏恢復算法主要從l0和l1兩種范數(shù)的稀疏問題條件出發(fā)。對于l0范數(shù)稀疏求解中,l0范數(shù)的最小化在數(shù)學上是一個非確定性多項式;對于l1稀疏問題中,常使用凸最優(yōu)算法,但l1范數(shù)最小化可能會極大降低求解的稀疏性,故選擇基于l0范數(shù)約束的稀疏表示,目標函數(shù)為
人大代表退出機制的建立,打破了代表身份終屆制,體現(xiàn)了“人民選我當代表,我當代表為人民”的宗旨,辭退代表不來虛的來實的,更加激發(fā)了代表履職的責任感。本屆以來,已對15名履職不認真、選民評價不高的市、鎮(zhèn)人大代表進行約談,其中已有9名主動辭去代表職務,勸辭2人,并將代表辭職情況及時向社會公開。
由2.1 可知,通過可轉(zhuǎn)變的矩陣B,密集的靜電信號可以用稀疏向量進行表示,且字典B是定義稀疏域的基礎要素。滾動軸承的故障信號為離散非平穩(wěn)信號,將信號分為故障瞬時沖擊特征成分x和隨機噪聲f,則滾動軸承沖擊響應模型y表示為
式中:η為極小的正數(shù)。
第二步:初始化。殘差r(0)=y、支撐集 Λ(0)=?、最大迭代次數(shù)maxiter。
為了對稀疏信號故障特征能更加凸顯,需要對小波字典中參數(shù)進行約束。由信號瞬時沖擊形態(tài)參數(shù)估值和,參數(shù)設置為
黨員網(wǎng)格化,建有詳細檔案,接受黨支部的領導,按時參加黨支部的組織生活,為黨支部工作獻計獻策,爭做“樹良好形象,做文明老人”的典范,主動“學雷鋒精神、傳遞正能量”,積極參加利民活動,弘揚了志愿者精神,倡導了社會新風。
構(gòu)造字典主要有解析字典和學習字典兩大類,解析字典求解速度快,由于稀疏表示對象具有較高的指向性,故選用解析字典方式構(gòu)造解析字典。軸承故障診斷要求對故障特征清晰識別,滾動軸承工作過程存在沖擊響應,具有震蕩衰減特性,Laplace小波可以有效提取故障特征分量[14],故本文采用Laplace 為原子構(gòu)造解析字典,Laplace 小波字典數(shù)學表達式為
2.3.1 StOMP 算法
2012 年,Donoho D 等[15-17]將正交匹配跟蹤算法與最小角度回歸(least angle regression,LAR)相結(jié)合[18],隨后提出分段正交匹配跟蹤算法StOMP?;趌0范數(shù)約束的稀疏優(yōu)化問題,由于原軸承靜電信號存在信號中背景噪聲、隨機脈沖等影響,將目標函數(shù)(6)改寫為
若稀疏向量d的解向量有e個非零元素并且滿足e?K時,則向量d的稀疏度為e。由式(3)可得,解向量d具有若干組解,由此引入稀疏約束條件,則稀疏求解的目標函數(shù)d為
StOMP 算法主要改進支撐集的更新方法,該算法通過設置閾值準則來替代最大化的準則,其具體的定義為:
式中:r(t-1)為t-1次迭代后的殘差值;ts為比例系數(shù);ts∈(0,1)。
通過設置閾值可以使支撐集里的原子單次更新匹配更多原子,繼而提高搜索匹配原子的速度,繼而提高StOMP 算法求解最優(yōu)解的速度。
StOMP 算法的核心為將解析字典中匹配原子與殘差做內(nèi)積,求解相關(guān)度,然后預設閾值,把相關(guān)度大于該閾值的所有匹配原子全部篩選出來,更新到支撐集 Λt中,再次由最小二乘法求出權(quán)值 θt,最后進行迭代算出,求出該組靜電信號重構(gòu)的最優(yōu)解。StOMP 算法求解過程如圖2 所示。
從專利數(shù)量對比情況來看,美國桔皮書所列專利中,非活性成分專利數(shù)量遠超過活性成分專利,即使新分子實體藥物也是如此。對1994-2006年新批準上市的214個新分子實體藥物進行分析發(fā)現(xiàn),在這些藥物中,共有716個專利列入桔皮書中,平均每個新分子實體藥物擁有3.34個專利,其中活性成分專利1.13個,非活性成分專利2.21個。將非活性成分進一步分類為使用方法專利和藥物產(chǎn)品專利,發(fā)現(xiàn)每個新分子實體藥物平均擁有1.34個使用方法專利和0.87個藥物產(chǎn)品專利,使用方法專利數(shù)量明顯高于藥物產(chǎn)品專利。④ 同注釋③。
圖2 StOMP 算法求解過程
靜電監(jiān)測技術(shù)的優(yōu)勢就是能夠在滾動軸承磨損區(qū)域監(jiān)測到軸承發(fā)生損壞的早期特征成分,而在實際實驗過程中,往往存在著工頻、背頻以及隨機脈沖等干擾,軸承故障并沒有凸顯出來,因此StOMP算法就會使軸承故障特征成分的識別精度下降。
2.3.2 CcStOMP 算法
CcStOMP 算法是在StOMP 基礎上進行優(yōu)化,在StOMP 算法中增加聚類收縮機制,在原子搜索過程中對待選原子進行兩次濾波得到的,使得CcStOMP 在故障特征識別精度上更有優(yōu)勢。為提高稀疏精度,提出聚類收縮機制,運用k近鄰算法對支撐集里面的冗余原子進行剔除,其核心思想是將原子映射到參數(shù)空間獲得 ψi,繼而將所用信號原子使用k近鄰聚類,保留最匹配的原子同時剔除冗余原子。聚類收縮機制原理如圖3 所示。
圖3 聚類收縮機制原理圖
在聚類收縮算法的基礎上對搜索原子進行二次過濾,最終可以完全剔除支撐集中的冗余原子。由式(9)知,小波字典中匹配原子通過各參數(shù)估值進行對鄰域的擴進,這些原子組成參數(shù)空間即span{ζ,f,τ},考慮到軸承故障特征成分是離散性的,根據(jù)閾值準則進行原子搜索時,針對進行原子搜索時可能存在若干個相似原子共同被篩選的問題,因此考慮迭代中間過程支撐集 Λt,函數(shù)為
式中:aki和akj是同一領域的兩個近似的原子,σ是極小的正數(shù)。由此可得支撐集滿足下列約束條件:
如果將上述計算公式中的Qi,t定義為第t期消費結(jié)構(gòu)中某一類消費性支出的占比,Qi,t-1定義為第t-1期該類消費性支出的占比時,就可根據(jù)Moore值算得不同時期消費結(jié)構(gòu)的變動情況。
完整的CcStOMP 算法求解過程如下。
第一步:輸入原始靜電信號y、字典B和最大距離尺度maxdis。
式中:f為震蕩頻率;ζ為阻尼比;τ為時間尺度;Wm為字典支撐區(qū)間。
為進一步驗證CcStOMP 算法對實際滾動軸承故障診斷的有效性,基于靜電監(jiān)測原理,搭建圖4所示的滾動軸承靜電監(jiān)測實驗平臺。實驗平臺主要由電機、聯(lián)軸器、傳動軸、傳感器夾具、磨損區(qū)域靜電傳感器、軸承座、磁粉制動器、振動傳感器、傳動帶、滾動軸承和電機控制儀等構(gòu)成。電機和傳動軸通過聯(lián)軸器1 接觸式連接,電機旋轉(zhuǎn)帶動旋轉(zhuǎn)軸轉(zhuǎn)動,進而使安裝在軸承座內(nèi)軸承旋轉(zhuǎn),兩者速度相同,軸承內(nèi)圈旋轉(zhuǎn),外圈不動,靜電傳感器由夾具固定在實驗平臺靠近磨損監(jiān)測區(qū)域處。本實驗主要監(jiān)測指標為摩擦副的荷電水平,所采用的磨損區(qū)域靜電傳感器主要結(jié)構(gòu)如圖5 所示,靜電傳感器主要由探極、絕緣層、屏蔽層、后端蓋和信號輸出端組成。實驗選用SKF-6204-2Z 型深溝球軸承為研究對象,通過在滾動軸承上預置外圈故障和滾子故障,選取缺口寬度為2.5 mm 作為軸承重度故障,滾動軸承靜電監(jiān)測實驗中所用到的滾動軸承實驗件分別如圖6a 和圖6b 所示。
扇三角洲沉積主要發(fā)育于沙二段和沙三段,可進一步將其劃分為扇三角洲平原亞相、扇三角洲前緣亞相和前扇三角洲亞相。
圖4 滾動軸承靜電監(jiān)測實驗平臺
圖5 靜電傳感器
圖6 滾動軸承實驗件
數(shù)據(jù)采集實驗步驟如下:電機控制儀預設轉(zhuǎn)速,電機帶動旋轉(zhuǎn)軸上的滾動軸承轉(zhuǎn)動發(fā)生磨損,并在磨損區(qū)域產(chǎn)生大量帶電顆粒。靜電傳感器探極近端產(chǎn)生的反電位電荷驅(qū)使同電位電荷到探極表面遠端,并通過信號調(diào)理電路使電荷信號轉(zhuǎn)換為電壓信號,經(jīng)過電荷放大器對靜電信號放大,然后通過數(shù)據(jù)采集卡轉(zhuǎn)為數(shù)字信號并傳輸?shù)綌?shù)據(jù)采集界面系統(tǒng)。靜電信號采集步驟如圖7 所示。
為此,筆者結(jié)合當前水稻高速插秧機的工作特點,開發(fā)了基于插秧機插植部傾角檢測的可橫向仿形的秧苗插深自適應調(diào)節(jié)系統(tǒng),并在虛擬仿真實驗平臺上開展了該系統(tǒng)靜態(tài)和動態(tài)實驗研究,力圖解決插秧盤隨水田表面橫向起伏引起的秧苗插植深度不一致現(xiàn)象,確保機插質(zhì)量穩(wěn)定,實現(xiàn)淺栽、穩(wěn)定、高產(chǎn)。
圖7 靜電信號采集步驟
獨立石油公司在資金和規(guī)模實力上的差距,使得其有選擇地重點進入某一些國家,其開展業(yè)務的國家數(shù)量遠低于國際大石油公司和國家石油公司。2014年國際油價大幅下降以來,獨立石油公司實施聚焦策略,選擇的國家目標更加精準。安納達科(Anadarko)除本土外,勘探資產(chǎn)主要分布在哥倫比亞;阿帕奇退出澳大利亞和加拿大,勘探區(qū)塊主要分布在埃及、北海、蘇里南。獨立石油公司勘探重心非常清晰,例如,康菲重點打造本土非常規(guī)以及阿拉斯加常規(guī)勘探,阿帕奇除了本土非常規(guī)外,還將埃及作為其重要勘探領域。
通過召開座談會、個別談話、走訪等形式的調(diào)研知道,改革轉(zhuǎn)型時期的企業(yè)員工的思想動態(tài)主要表現(xiàn)為以下幾種形式:
表1 實驗滾動軸承主要參數(shù)
軸承外圈實驗數(shù)據(jù)原始靜電信號時域圖、頻域圖分別如圖8a 和圖8b 所示,fBPFO為軸承外圈故障頻率;軸承滾子實驗數(shù)據(jù)原始靜電信號時域圖、頻域圖分別如圖9a 和圖9b 所示,fBSF為軸承滾子故障特征頻率。由圖8a 和圖9a 可知,外圈故障幅值峰值強度基本維持在10 mV 內(nèi),而滾子故障幅值變化較小,幅值峰值強度維持在5 mV,兩者靜電波動較大,說明靜電感應明顯。觀察圖8b 外圈故障頻域圖,在頻率191 Hz 處有明顯峰值,但該頻率附近也存在其他明顯峰值,從而對故障特征造成干擾,其2 倍頻不夠突出,易被周圍噪聲干擾;觀察圖9b 滾子故障頻域圖,能找到故障頻率127 Hz,但故障頻率峰值較低,不能快速定位其故障特征頻率,且故障頻率周圍也有其他突出的噪聲干擾,不利于故障特征提取。由此可得,靜電監(jiān)測技術(shù)能有效實現(xiàn)對滾動軸承的故障特征提取,但容易受到噪聲干擾,因此需要使用一種性能優(yōu)異的故障特征提取方法,完成滾動軸承故障診斷任務。
圖9 滾子故障實驗原始靜電信號
為驗證CcStOMP 算法在滾動軸承靜電監(jiān)測信號有較好去噪效果,在3.1 和3.2 兩小節(jié)中以軸承外圈和滾子故障為代表進行實驗結(jié)果分析,使用StOMP 算法和CcStOMP 算法對同一段靜電信號進行處理,得到重構(gòu)信號的時域圖、稀疏域圖和包絡譜圖,分析其信號重構(gòu)能力和去噪能力。
活性炭基脫氯劑的飽和氯容和穿透氯容采用自制的動態(tài)吸附裝置進行測定,如圖2所示。將活化后的活性炭基脫氯劑破碎至0.425~0.850 mm,取3 mL裝填入樣品管,按照圖2連接好裝置。將純氮氣以空速1000 h-1的模擬工況通入鹽酸溶液中,然后通過濃硫酸脫除水蒸氣,經(jīng)樣品管吸附后,用裝有500.0 mL預先加入溴甲酚橙指示劑的去離子水的吸收瓶吸收尾氣,待吸收瓶中水的顏色由藍色變?yōu)辄S色時,視為穿透。再取相同量的活性炭基脫氯劑,重復實驗,至穿透后黃色不再加深,視為飽和。穿透氯容和飽和氯容可由式(1)、(2)分別計算得到。
StOMP 處理后的外圈故障時域圖、稀疏域圖和包絡譜圖分別如圖10a、10b 和10c 所示,CcStOMP處理后的外圈故障時域圖、稀疏域圖和包絡譜圖分別如圖11a、11b 和11c 所示。
圖10 外圈故障StOMP 算法重構(gòu)結(jié)果
圖11 外圈故障CcStOMP 算法重構(gòu)結(jié)果
由圖10a 和圖11a 對比分析可知,在信號重構(gòu)的時域波形圖中,CcStOMP 算法故障特征頻率更加突出。觀察圖10b 和圖11b 稀疏系數(shù)分布圖,CcStOMP 算法較StOMP 算法稀疏性能更加優(yōu)越,同時支撐集中的匹配原子分布更加均勻密集,權(quán)值θt也較平緩。對比圖10c 和圖11c 可以看出,StOMP算法能識別出故障特征頻率,但其故障特征不夠突出,容易被附近噪聲干擾,且不能有效識別其故障特征倍頻;CcStOMP 算法能清晰識別出外圈故障的特征頻率fBPFO=191 Hz,同時能清晰識別出故障特征頻率的2 倍頻和3 倍頻,且故障特征頻率周圍噪聲干擾小,能快速定位故障特征頻率,兩種算法對軸承外圈故障特征提取效果見表2。
表2 外圈故障特征提取效果對比
StOMP 處理后的滾子故障時域圖、稀疏域圖和包絡譜圖分別如圖12a、12b 和12c 所示,CcStOMP處理后的滾子故障時域圖、稀疏域圖和包絡譜圖分別如圖13a、13b 和13c 所示。
工匠們喜歡不斷雕琢自己的產(chǎn)品,不斷改善自己的工藝,享受著產(chǎn)品在雙手中升華的過程。黃國平與其帶領的蘇印總廠,恰以工匠精神、創(chuàng)新的思想,悉心打磨每一件產(chǎn)品,讓每一個從該企業(yè)出去的作品都是讓客戶滿意的精品。
圖12 滾子故障StOMP 算法重構(gòu)結(jié)果
圖13 滾子故障CcStOMP 算法重構(gòu)結(jié)果
觀察圖12a 和圖13a 可知,在信號重構(gòu)時域波形上來看CcStOMP 算法故障特征頻率更加明顯。觀察圖12b 和圖13b 可知,StOMP 算法的稀疏性較差,支撐集中的匹配原子少,且分布零散;CcStOMP算法的稀疏性好,支撐集中的匹配原子密集,分布均勻,權(quán)值也較平緩。觀察圖12c 和圖13c 可知,CcStOMP 算法能夠比較清楚地識別出滾子故障特征頻率,同時也能清晰地識別出故障特征頻率的2倍頻而不易周圍噪聲干擾。由滾子故障理論故障頻率可知,實驗數(shù)據(jù)靜電信號的故障特征頻率與理論故障特征頻率之間存在0.02%的微弱誤差,分析存在誤差的原因可能是實際轉(zhuǎn)動頻率和理論轉(zhuǎn)動頻率存在偏差,繼而導致了該誤差存在,因誤差較小可以忽略該誤差的存在,具體兩種算法對軸承滾子故障特征提取效果見表3。
表3 滾子故障特征提取效果對比
兩類稀疏算法,StOMP 在算法中能識別故障特征頻率成分,但相比于CcStOMP 算法,識別精度存在明顯不足,且容易被周圍噪聲干擾。這說明CcStOMP 算法對滾動軸承故障信號重構(gòu)和故障特征提取有較好的性能。
針對滾動軸承故障靜電監(jiān)測信號故障特征不易提取,容易被噪聲淹沒的問題,將稀疏表示理論引入到靜電信號處理,采用基于CcStOMP 的稀疏表示方法。通過實驗數(shù)據(jù)對比分析,得出以下結(jié)論:
(1)針對滾動軸承靜電監(jiān)測技術(shù)特點,利用稀疏恢復算法對原瞬態(tài)故障特征進行信號重構(gòu),提取故障特征成分,再進行包絡分析,可以實現(xiàn)對滾動軸承故障靜電監(jiān)測信號的故障特征提取。
(2)將CcStOMP 與StOMP 兩種算法通過實驗靜電信號進行去噪效果對比分析。實驗結(jié)果分析表明,StOMP 能夠去除靜電信號中的噪聲,找到滾動軸承的故障特征頻率,但StOMP 算法精度低,準確性不高;CcStOMP 相比于StOMP 能夠有效消除噪聲干擾并保留有效靜電信號,驗證了CcStOMP算法對靜電信號故障特征提取的有效性。