冷 宇 李 棟

（1.上海科學技術(shù)職業(yè)學院，上海201800；2.華域電動系統(tǒng)有限公司，上海201203）

0 引言

電動機是現(xiàn)代化工業(yè)時代必不可少的生產(chǎn)動力提供設備，其在工業(yè)、農(nóng)業(yè)、商業(yè)領域都有廣泛應用。對于電動機故障的檢測和預防可以為生產(chǎn)減少損失，并避免不必要的安全事故。在電機故障檢測中，最重要也是最常用的數(shù)據(jù)就是電機軸承的振動數(shù)據(jù)，為了實現(xiàn)電機軸承的振動數(shù)據(jù)（加速度、速度、位移）檢測，本文采用了ICP傳感器（加速度傳感器）。為了實現(xiàn)本文設計的良好拓展性，將處理后的數(shù)據(jù)使用RS485通信協(xié)議與上位機進行數(shù)據(jù)交換，從而進一步實現(xiàn)更多功能的數(shù)據(jù)分析，以進一步完善振動檢測。

1 方案介紹

本文所設計的檢測裝置硬件電路主要分為六大塊，分別為輔助電源電路、DSP控制電路、RS485通信電路、信號調(diào)理電路、報警電路、LCD顯示電路，其功能框圖如圖1所示。

圖1 功能框圖

2 硬件電路設計與仿真

本小節(jié)主要介紹硬件電路的設計。由圖1可以看出來硬件電路主要分為五大塊，即輔助電源、DSP控制電路、RS485通信電路、信號調(diào)理電路和LCD顯示與報警電路。

輔助電源主要實現(xiàn)對各個電路模塊的供電，整個系統(tǒng)需要的電源有+5 V、-5 V、4 V、恒流源和隔離+5 V。由于本文的ICP傳感器是使用恒流源供電，選擇LM317芯片組成恒流源。本文所選擇的ICP傳感器正常工作時所需要的外接恒流源需提供2～20 mA電流，這里設從1號管腳ADJ所流出的電流為Iadj，由于Iadj?Icc，所以忽略Iadj，有：

所以，調(diào)節(jié)電阻R23即可得到想輸出的恒流值。

主控芯片選擇TMS320F2812。通信電路設計，本文選擇了RS485通信協(xié)議。使用MAX491芯片作為收發(fā)器實現(xiàn)通信，均采用光耦進行隔離。

壓電式ICP傳感器內(nèi)置電荷放大器電路模塊，直接輸出電壓，無需設計專門的電荷放大器，使用簡單，抗干擾性能好，可長線使用，其型號為JM411-22。

本文使用二階低通巴特沃斯濾波器，截止頻率理論計算：

U2A運放實現(xiàn)了隔直和電壓抬升的作用，設輸入電壓為U2A，輸出電壓為U2Ao，則有公式：

U1A運放實現(xiàn)了比例縮小，設輸入電壓為U2B，輸出電壓為U2Bo，則有公式：

為了驗證電路設計的可靠性和穩(wěn)定性，本文使用Multism對設計電路進行了仿真，仿真結(jié)果如圖2所示。

圖2 仿真結(jié)果

3 軟件算法設計

軟件控制主要分為五大塊：小波消噪、頻域濾波、去除趨勢項、數(shù)字積分和信號顯示，用于實現(xiàn)對采樣信號的處理。軟件算法框圖如圖3所示。

下面詳細介紹各個模塊的作用與功能。

3.1 小波消噪

直接采集的加速度信號具有較多的噪聲信號，這些信號的存在會給信號處理帶來較大誤差，為使速度信號和位移信號具有較高的準確度，則需進行消噪，本文選用小波消噪法。

3.1.1 小波消噪的數(shù)學模型及基本方法

首先，為了方便建立數(shù)學模型，必須對噪聲信號作出三個假設：（1）噪聲經(jīng)小波變換后大多數(shù)小波系數(shù)為零或近似為零；（2）噪聲均勻地分布在所有系數(shù)中；（3）噪聲水平不是太高。

圖3 軟件算法框圖

設觀測得到的數(shù)據(jù)：

由公式可知，觀測數(shù)據(jù)是由真實信號g和噪聲ε疊加而成，f=g+ε為向量形式。在理想情況下，εi需要滿足以下三點：（1）正態(tài)噪聲；（2）不相關(guān)噪聲；（3）方差為常量。

由于現(xiàn)實中較難達到要求，所以必須將假設條件放寬，否則不能滿足實際測量的需要。最終目的是使用觀測到的f，得到估計的g。設觀測數(shù)據(jù)由小波變換后得：

設W（·）和W-1（·）分別表示小波變換和逆變換算子，則小波消噪的過程可定義為：

（1）w=W（f）——小波變換；

（2）wt=D（w，t）——對小波系數(shù)進行非線性處理，以濾除噪聲；

（3）g?=W-1（wt）——小波逆變換。

非線性濾波算子D（·，·）為去噪問題的關(guān)鍵，通過選擇不同的W（·）和D（·，·），則有多種去噪的辦法。小波消噪的關(guān)鍵在于第（2）步——按照一定的準則修改小波系數(shù)，從而在損失較小信號的前提下達到去噪效果。目前小波消噪有貝葉斯法和非貝葉斯法，一般采用非貝葉斯法，非貝葉斯法大致可分為三種：（1）模極大值重構(gòu)濾波法消噪；（2）空域相關(guān)濾波消噪；（3）小波閾值濾波消噪。

對比以上三種濾波算法可知，小波閾值濾波消噪效果相對較好，實現(xiàn)也較簡單。本文采用該方法，僅需改變閾值即可得出不同的濾波效果。

3.1.2 小波閾值濾波消噪

有用信號經(jīng)過小波變換后的小波系數(shù)具有幅值大、數(shù)量少的特點，并且包含信號的重要信息，而噪聲的小波系數(shù)則幅值較小。正是基于以上特點，在不同尺度上選取不同閾值，將小于該值的小波系數(shù)歸零，保留大于閾值的系數(shù)，就可以實現(xiàn)去噪，最后進行逆變換，實現(xiàn)信號復原。

Donoho對閾值濾波提出了兩個前提條件：

（1）光滑性——在大概率情況下，g?至少跟g有同樣的光滑度；

（2）適應性——g?是最小均方差估計。

由條件（1）可得，當N→∞時，下式近似為1的概率成立：

式中，C1為常數(shù)，小波域則有式（8）的成立。

式中，1≤i≤N為位置，j為尺度；θj，i為真實信號在尺度j上的第i個小波系數(shù)。

對于條件（2），可以理解為對式（9）求最小值，與對式（10）求最小值互為等價，Donoho證明，此時g?、θ?必須滿足式（11）。

由式（11）可知，對于任何，取θj，i=0，將滿足式（10）。這相當于認為當，θj，i由噪聲所產(chǎn)生，因此，可取閾值為

閾值濾波有硬閾值濾波和軟閾值濾波之分，下面分別介紹兩種濾波的方法和步驟。

硬閾值濾波消噪的步驟為：

（1）對信號求小波變換；

（2）除了最粗尺度信號外，對細節(jié)信號作閾值處理，閾值當某位置小波變換值大于閾值時，保留原值，否則置零，用公式表示為：

（3）利用小波變換重構(gòu)，求出信號的濾波值。

軟閾值濾波消噪的步驟為：

（1）對信號求小波變換；

（2）除了最粗尺度信號外，將各細節(jié)信號作閾值處理，閥一旦出現(xiàn)閾值小于小波變換值時，向著減小系數(shù)幅值的方向做一個收縮t，否則置零，用公式表示為：

式中，sgn（x）為符號函數(shù)。

（3）進行小波變換重構(gòu)，求出信號的濾波值。

綜上所述，如果要實現(xiàn)整個算法，則需知道噪聲的方差，噪聲標準差σ的估計已經(jīng)由Donoho給出。

閾值和閾值函數(shù)是小波閾值濾波算法中的兩個基本要素，閾值對濾波效果具有較大影響，因此閾值的確定較為關(guān)鍵，而目前的文獻中給出了很多閾值的確定方法。

使用Matlab仿真，結(jié)果如圖4所示。通過仿真可以看出，小波降噪具有良好的去噪效果。

3.2 FFT濾波

消噪后的信號仍存在干擾信號，所以需要做濾波處理，本設計采用FFT濾波（頻域濾波）。FFT濾波是基于FFT快速傅里葉算法，將采樣信號進行離散傅里葉變換，分析頻譜，并根據(jù)要求將濾除頻率成分設置為零，最后采用IFFT（傅里葉逆變換）快速算法恢復出時域信號。該方法的頻域表達式為：

圖4 仿真結(jié)果

式中，X為輸入信號X的離散傅里葉變換；H為濾波響應函數(shù)。

設fu為上限截止頻率，Δf為頻率分辨率，本文選擇低通濾波器，即FFT濾波中采用低通濾波器濾波，其頻率響應函數(shù)如下：

3.3 去除趨勢項

在實際測量中，由于現(xiàn)實環(huán)境的不斷變化，會產(chǎn)生零點漂移、傳感器測量特性改變等等不良因素。而這些因素會導致測量結(jié)果偏離基線，并隨時間而改變，這一現(xiàn)象被稱為趨勢項。而數(shù)據(jù)一旦進行積分，趨勢項所帶來的影響則不可忽視，若是再二次積分，則有可能使積分的結(jié)果完全失真。本文選用多項式最小二乘法，原理如下：

將｛xk｝（k=1，2，3，…，n）設為實際所獲得的振動信號的采樣，因為在現(xiàn)實中是固定步長進行采樣，令Δt=1為采樣步長時間，并設如下多項式：

由x?k與xk的差值平方最小可求aj（j=0，1，2，…，m），即：

由式（16）可得：

通過對式（18）求解，則可以求出m+1個待定系數(shù)aj（j=0，1，2，…，m）。若m=0，則可解得：

若m=1，則可解得：

當m≥2時，也可以同理算出。值得注意的是，實際采樣中一般取m=1～3來去除趨勢項。

使用Matlab進行仿真，仿真結(jié)果如圖5所示。

3.4 數(shù)字積分

在經(jīng)過一系列的軟件信號調(diào)理之后，接下來就需要進行數(shù)字積分，以求出速度與位移信號，并將位移信號與設定的保護值實時比較，一旦發(fā)現(xiàn)比設定值大，那么將會進入中斷，執(zhí)行報警程序。

圖5 仿真結(jié)果

本文所設計的數(shù)字積分采取簡單的梯形法。梯形法是將一個被積函數(shù)的曲線分成很多個小區(qū)域，并將分割的小區(qū)域看為一個個梯形。令x（i-1）和x（i）為相鄰的兩個點，Δt為積分步長，則可得到數(shù)學表達式：

其中，k=1，2，3，…，N。

使用Matlab對梯形法積分算法進行仿真，得到如圖6所示仿真結(jié)果。

圖6 梯形法積分仿真結(jié)果

4 結(jié)語

本文總結(jié)并給出了電機軸承處振動在線監(jiān)測設備的硬件和軟件設計，并給出小波消噪的基本算法，實現(xiàn)對振動數(shù)據(jù)的監(jiān)測，從而最終實現(xiàn)對電機軸承的保護，并將相關(guān)數(shù)據(jù)上傳至上位機，以便上位機實現(xiàn)對大數(shù)據(jù)的分析。

[1]羊初發(fā).基于EMD的時頻分析與濾波研究[D].成都：電子科技大學，2009.

[2]楊君，呂鏡清.FFT濾波誤差分析[J].通信技術(shù)，2010，43（6）：235-237，240.

[3]朱可恒.滾動軸承振動信號特征提取及診斷方法研究[D].大連：大連理工大學，2013.

[4]HAMMAN J,VAN DER MERWE F S.Voltage harmonics generated by voltage-fed inverters using PWM natural sampling[J].IEEE Transactions on Power Electronics,1988,3(3):297-302.

[5]王賽男，邢冬梅.基于MATLAB的小波閾值濾波與應用[J].軟件，2015，36（10）：102-104.