摘 要：以地下爆炸產(chǎn)生的震動(dòng)信號(hào)為研究對(duì)象，基于LabVIEW軟件，主要設(shè)計(jì)軟件模塊完成信號(hào)時(shí)域、頻域、時(shí)頻域分析。整個(gè)軟件采用自頂向下集成的模塊化編程思路，大大提高了軟件的執(zhí)行效率和可擴(kuò)展性。測(cè)試表明，該軟件系統(tǒng)可以很好的應(yīng)用于震動(dòng)信號(hào)分析方面。

關(guān)鍵詞：震動(dòng)信號(hào)；時(shí)頻域分析；系統(tǒng)設(shè)計(jì)

在實(shí)際震動(dòng)測(cè)試試驗(yàn)中，獲取的震動(dòng)信號(hào)屬于非平穩(wěn)信號(hào)。直接些信號(hào)進(jìn)行觀察分析往往不能滿足們的需求，很多信號(hào)的特征也不能直觀的被顯示出來(lái)。因此，對(duì)信號(hào)進(jìn)行時(shí)頻域分析則成了研究的重點(diǎn)課題。時(shí)域分析可以直觀的給出信號(hào)在時(shí)域的特征，頻域分析可以給出信號(hào)在頻域的特征，時(shí)頻域聯(lián)合分析可以從圖中分析信號(hào)時(shí)域與頻域之間的聯(lián)系。這些往往對(duì)信號(hào)的特征提取分析有著重要的作用。在地下采集到的震動(dòng)信號(hào)，往往由于傳播介質(zhì)的復(fù)雜性，使其含有大量的噪聲信號(hào)或者和其他不明成分，對(duì)后續(xù)震源位置定位有著重要的影響。因此，文章以LabVIEW為軟件開(kāi)發(fā)平臺(tái)，設(shè)計(jì)時(shí)頻域分析模塊，對(duì)震動(dòng)信號(hào)進(jìn)行相關(guān)的分析，為后續(xù)的震源定位提供良好的分析條件。

1 軟件整體結(jié)構(gòu)設(shè)計(jì)

設(shè)計(jì)的軟件要良好遵從的界面交互性、軟件可移植性、開(kāi)放性等原則。文章基于LabVIEW開(kāi)放平臺(tái)，其具有計(jì)算機(jī)強(qiáng)大的圖形環(huán)境，采用可視化的圖形編程語(yǔ)言和平臺(tái)，跟其他編程語(yǔ)言相比，其具有如下優(yōu)勢(shì)：（1）虛擬儀器具有良好的擴(kuò)展性能，能實(shí)時(shí)根據(jù)需要完成軟件的擴(kuò)展。（2）由于虛擬儀器測(cè)控測(cè)試技術(shù)，其能與硬件很好連接，工程上經(jīng)常用硬件開(kāi)發(fā)平臺(tái)，利用虛擬儀器作為上位機(jī)完成相應(yīng)的項(xiàng)目。（3）虛擬儀器是強(qiáng)大的圖形化語(yǔ)言，可以直觀方面的供開(kāi)發(fā)人員開(kāi)發(fā)程序。設(shè)計(jì)軟件從功能劃分為時(shí)域分析、頻域分析、時(shí)頻域分析。時(shí)域分析又包括濾波功能，頻域分析包括FFT功能，時(shí)頻域分析包括小波去噪功能。

2 軟件詳細(xì)功能設(shè)計(jì)

2.1 時(shí)域分析

時(shí)域分析是主要是根據(jù)分析信號(hào)的特點(diǎn)，有些特點(diǎn)我們不能直觀的看出來(lái)，因此時(shí)域分析要根據(jù)信號(hào)時(shí)間信息來(lái)獲得相應(yīng)的特征 [1]。本系統(tǒng)設(shè)計(jì)的時(shí)域分析包括濾波功能。

濾波功能模塊設(shè)計(jì)：軟件設(shè)計(jì)的濾波有Butterworth濾波器、Chebshev濾波器、橢圓濾波器、貝賽爾濾波器。

（1）巴特沃斯濾波器的特點(diǎn)是其時(shí)頻域曲線較為平臺(tái)，沒(méi)有較大的起伏，在阻帶內(nèi)從最高點(diǎn)慢慢下降至零。（2）切比雪夫?yàn)V波和理想濾波的頻率幅值響應(yīng)曲線相互誤差系數(shù)為最小，但是在其頻率的通頻帶內(nèi)存在幅度上下波動(dòng)。（3）橢圓濾波器相比其他類型的濾波器，在系數(shù)階數(shù)相同的條件下有著很小通頻帶和阻帶輕微波動(dòng)。（4）貝賽爾（Bessel）濾波器是具有最大平坦值的群延遲（線性相位響應(yīng)）的線性過(guò)濾器。貝賽爾濾波器通常應(yīng)用在音頻聲響系統(tǒng)中。貝賽爾濾波容易造成真?zhèn)€區(qū)域內(nèi)的群延遲。

上述濾波器的實(shí)現(xiàn)，可以通過(guò)調(diào)用LabVIEW中的函數(shù)。用戶則通過(guò)設(shè)置高低截止頻率，濾波器的類型完成對(duì)信號(hào)的濾波操作。

2.2 頻域分析

FFT功能模塊設(shè)計(jì)：FFT是一種DFT的高效算法，稱為快速傅立葉變換（fast Fourier transform），它根據(jù)離散傅氏變換的奇、偶、虛、實(shí)等特性，對(duì)離散傅立葉變換的算法進(jìn)行改進(jìn)獲得的。在編寫(xiě)該功能模塊時(shí)，可以調(diào)用LabVIEW中的FFT功能函數(shù)。用戶則可以通過(guò)設(shè)置FFT的頻率來(lái)分析信號(hào)。

2.3 時(shí)頻域分析

小波閾值去噪：由于爆炸產(chǎn)生的震動(dòng)信號(hào)是非平穩(wěn)信號(hào)，通常工程上所用的信號(hào)去噪方法包括傅里葉變換、維納濾波、卡爾曼濾波、自適應(yīng)濾波等，其都要求信號(hào)是線性的，并不適于震動(dòng)信號(hào)這樣的非平穩(wěn)信號(hào)[2]。小波去噪的思想：震動(dòng)信號(hào)經(jīng)過(guò)若干層小波變換后，在幅值較大、數(shù)量較少數(shù)幾個(gè)小波系數(shù)上集中其大部分的能量；而噪聲經(jīng)過(guò)小波變換后，在其高頻信息則主要集中了噪聲的特征，其能量值主要分布于幅值很小、小波系數(shù)較多的幾個(gè)小波系數(shù)上。根據(jù)上述特點(diǎn)，得出噪聲的小波系數(shù)值必然小于信號(hào)的小波系數(shù)值。因此，根據(jù)以上分析，軟件使用利用小波閾值變換能夠有效地實(shí)現(xiàn)信號(hào)與噪聲的分離[3]。

在軟件設(shè)計(jì)中，該小波去噪模塊主要通過(guò)LabVIEW和MathScript結(jié)合完成對(duì)信號(hào)的消噪設(shè)計(jì)。

3 軟件界面效果圖

根據(jù)上述功能的實(shí)現(xiàn)，編寫(xiě)了軟件界面圖如圖1所示。用戶通過(guò)設(shè)定采樣頻率、小波分解層數(shù)、小波的類型、低通截止頻率、高通截止頻率以及濾波的類型，得出相應(yīng)的波形信號(hào)圖。

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

本系統(tǒng)主要針對(duì)震動(dòng)信號(hào)分析處理為目的，設(shè)計(jì)功能模塊完成對(duì)信號(hào)分析處理。為后續(xù)震動(dòng)信號(hào)的定位提供了良好的前提。

參考文獻(xiàn)

[1]馬瑞恒，等.時(shí)頻分布在爆破震動(dòng)信號(hào)處理中的應(yīng)用[J].工程爆破，2004（10）：8-12.

[2]何軍，于亞倫，梁文基.爆破震動(dòng)信號(hào)的小波分析[J].巖土工程學(xué)報(bào)，1998，20（1）：47-50.

[3]唐進(jìn)元，等.一種新的小波閾值函數(shù)及其在振動(dòng)信號(hào)去噪分析中的應(yīng)用[J].振動(dòng)與沖擊，2009（28）：118-121.

作者簡(jiǎn)介：宋璐（1990-），女，山西太原人，碩士研究生，主要從事通信與信息系統(tǒng)方向的研究。