申春龍,郭 鋒,朱成銀

(1.西南科技大學(xué)信息工程學(xué)院,四川 綿陽 621010；2.中國工程物理研究院核物理與化學(xué)研究所,四川 綿陽 621900)

瞬態(tài)脈沖信號測量與處理系統(tǒng)設(shè)計

申春龍1,郭 鋒1,朱成銀2

(1.西南科技大學(xué)信息工程學(xué)院,四川 綿陽 621010；2.中國工程物理研究院核物理與化學(xué)研究所,四川 綿陽 621900)

針對高強度脈沖中子控制系統(tǒng)的測試，需要對通過LabVIEW控制的數(shù)字示波器所捕獲到的大量相似瞬態(tài)脈沖信號進(jìn)一步測量與處理，以便篩選出符合系統(tǒng)控制要求的瞬態(tài)脈沖信號。考慮到瞬態(tài)脈沖信號的控制與采集設(shè)計也是通過LabVIEW進(jìn)行控制，然而單獨采用Matlab對信號進(jìn)行分析處理，難以同時滿足對海量數(shù)據(jù)進(jìn)行測量與處理的設(shè)計需求。由于LabVIEW具有自動多線程技術(shù)、圖形化設(shè)計語言和豐富信號處理模塊，并且Matlab軟件在對信號進(jìn)行小波降噪方面具有優(yōu)勢，利用LabVIEW軟件與Matlab軟件混合編程，實現(xiàn)集時域測量、頻域測量和小波降噪處理于一體的海量信號分析處理系統(tǒng)。實測結(jié)果驗證了該系統(tǒng)應(yīng)用于處理海量脈沖數(shù)據(jù)并篩選出有用信號的正確性及可行性，同時也表明了該系統(tǒng)具有操作簡單、響應(yīng)速度快、測量精度高的優(yōu)點，具有較高的實用價值。

數(shù)據(jù)采集； 脈沖信號； LabVIEW； Matlab； 時域測量； 頻域測量； 有限長序列； 小波降噪

0 引言

在許多實際工程應(yīng)用中，經(jīng)常需要對瞬態(tài)脈沖信號進(jìn)行分析和處理，例如：高強度脈沖中子控制系統(tǒng)、機械故障診斷、機器部件的瞬態(tài)沖擊等[1]。由于瞬態(tài)脈沖信號具有短時、頻帶寬和能量集中等特點，然而在高強度脈沖中子控制系統(tǒng)的測試中會在短時間內(nèi)產(chǎn)生大量的相似瞬態(tài)脈沖信號，憑借傳統(tǒng)的示波器與頻譜分析儀難以快速和準(zhǔn)確地對大量數(shù)據(jù)進(jìn)行測量與處理。因此，能夠快速準(zhǔn)確地對瞬態(tài)脈沖信號進(jìn)行測量與處理是快速尋找到符合控制探測信號要求的重要工作。

本文在充分研究了現(xiàn)有信號測量與處理技術(shù)的基礎(chǔ)上，提出利用LabVIEW軟件與Matlab軟件混合編程[2-3]對信號進(jìn)行測量與處理，從而從大量信號數(shù)據(jù)中快速篩選出符合設(shè)計要求的脈沖信號。

1 系統(tǒng)信號來源

在瞬態(tài)脈沖信號測量與處理中，需要產(chǎn)生脈寬為微秒至毫秒級的脈沖信號來驗證信號測量與處理的精準(zhǔn)性，所以需要對脈沖信號的來源有簡單的了解，以便對測試結(jié)果作進(jìn)一步的對比。在此簡單介紹信號來源以及存儲方式，并通過一個簡單的硬件連接框圖來表達(dá)信號的產(chǎn)生過程。硬件結(jié)構(gòu)如圖1所示。

圖1 硬件結(jié)構(gòu)圖

硬件結(jié)構(gòu)主要包括安裝有LabVIEW、Matlab和DAQmx軟件的工業(yè)控制計算機和與工業(yè)控制計算機相連的顯示器、輸入設(shè)備、數(shù)據(jù)采集卡和屏蔽接線盒，以及需要控制的泰克DPO7000系列示波器。

首先，需要通過輸入設(shè)備以及顯示器向工控機內(nèi)的LabVIEW軟件設(shè)置脈沖寬度。其次，通過LabVIEW對數(shù)據(jù)采集卡的控制產(chǎn)生脈沖信號，通過屏蔽接線盒將脈沖信號引出，并使用特殊探頭將信號接入示波器；同時需要通過LabVIEW對示波器進(jìn)行控制，等待滿足觸發(fā)條件的脈沖信號。最后，將滿足基本條件的脈沖信號通過LabVIEW存儲到工控機內(nèi)，以便后續(xù)分析[4]。

2 軟件設(shè)計

LabVIEW軟件作為人機交互的紐帶，是整個設(shè)計系統(tǒng)的核心[5]。脈沖信號的分析主要包括時域測量和頻域分析，對脈沖信號的處理包括濾波和降噪。為了能夠簡單高效地實現(xiàn)對脈沖信號的分析和處理，采用對程序進(jìn)行分析和處理的模塊化設(shè)計。由于通過LabVIEW對示波器的控制將示波器顯示的波形存儲為有限長的數(shù)字信號序列，因此為后續(xù)的分析和處理提供了便利[6]。因為LabVIEW不僅可以通過“條件結(jié)構(gòu)”VI設(shè)置并同時選擇需要分析處理的任意多個數(shù)據(jù)文件，而且LabVIEW具有高效性和并行執(zhí)行的能力，所以能夠?qū)崿F(xiàn)對多個不同脈沖信號進(jìn)行時域測量、頻域測量與信號處理。因為需要快速地對海量數(shù)據(jù)進(jìn)行挖掘，所以通過對需要測量的某一參數(shù)與其參考值進(jìn)行比較。如果超出閥值，那么該組數(shù)據(jù)則無效；反之則有效，并且對有效數(shù)據(jù)進(jìn)行標(biāo)注，從而高效實現(xiàn)了海量數(shù)據(jù)的篩選工作。系統(tǒng)工作流程如圖2所示。

圖2 系統(tǒng)工作流程圖

2.1 時域測量

通過LabVIEW中“從文件讀取波形”VI的選擇，并通過“文件路徑”控件指定波形文件的位置，將存儲在計算機內(nèi)的脈沖數(shù)據(jù)進(jìn)行波形還原(波形回放)，還原脈寬為1 ms的脈沖波形如圖3所示。

圖3 脈沖波形還原圖

脈沖信號的時域測量根據(jù)信號時域的表達(dá)式或者觀察的數(shù)據(jù)進(jìn)行計算，以實現(xiàn)對信號的穩(wěn)定性、瞬態(tài)特性和穩(wěn)態(tài)性能進(jìn)行分析。對瞬態(tài)脈沖的時域測量包括以下幾個方面：脈沖寬度、上升沿時間、下降沿時間、脈沖幅度[7]。

時域測量如表1所示。

表1 時域測量參數(shù)表

2.2 頻域測量

因為硬件產(chǎn)生的瞬態(tài)信號非常相似，直接從時域?qū)π盘栠M(jìn)行觀測難以區(qū)分，所以需要對頻域進(jìn)行進(jìn)一步的觀測。一方面，頻譜分析方法能夠有效地抑制干擾信號對頻率測量的影響；另一方面，頻譜分析方法在頻率測量時不存在電平閾值的問題，所以測量精度較高[8]。由于針對瞬態(tài)脈沖信號的采集所選擇的高性能的示波器，能提供最低為10 GS/s、最高達(dá)40 GS/s的實時采樣率，因此所存儲的采樣數(shù)據(jù)點數(shù)非常多。在頻域測量中，快速傅里葉變換(fast fourier transformation,FFT)方法是一種行之有效的方法。

軟件使用LabVIEW調(diào)用原始的采集數(shù)據(jù)，通過FFT方法將信號從時域轉(zhuǎn)換至頻域，從而更加方便直接地從頻域?qū)π盘栠M(jìn)行觀測。為了驗證設(shè)計的正確性，選擇RIGOL公司的信號發(fā)生器產(chǎn)生頻率為1 kHz、峰-峰值為5 V的正弦波，將其作為與脈沖信號對比和驗證信號。通過對正弦波信號和圖3所示的1 ms脈沖信號進(jìn)行FFT得到的頻譜圖如圖4所示。

圖4 信號FFT頻譜圖

為了能與正弦波頻譜進(jìn)行對比，將正弦波頻率與脈沖信號頻譜進(jìn)行放大觀察。通過對比驗證，對脈沖信號的FFT符合頻譜測量要求。

3 脈沖信號處理

從圖3可以直觀看出，脈沖信號中存在噪聲干擾，為了降低脈沖信號中的噪聲，以此獲得更加清晰的脈沖信號，需要對脈沖信號降噪處理。從圖4可以看出，脈沖信號主要由低頻部分組成，可以采用低通濾波器對其進(jìn)行處理，這樣噪聲集中的高頻部分就被去除掉了。但是由于低通濾波器在通帶和阻帶之間有不連續(xù)性，所以選擇比低通濾波器更具優(yōu)勢的巴特沃斯低通濾波器。它不會像理想低通濾波器一樣將一些有用的高頻信息也濾除。如果想要得到好的濾除效果，就要提高巴特沃斯濾波器的濾波階數(shù)，但隨著高巴特沃斯濾波器的濾波階數(shù)增加會出現(xiàn)明顯的振鈴現(xiàn)象[9]。為了更好地對脈沖信號進(jìn)行噪聲處理，引入小波變換降噪。

3.1 小波降噪處理

通過圖4的頻譜分析可以看出，利用傳統(tǒng)濾波器的設(shè)計思想難以對信號進(jìn)行降噪處理，以獲取有用信息。然而小波變換具有低熵性、多分辨率特性和選基的靈活性等特點，所以小波變換在信號降噪方面得到了廣泛使用[10]。

由于波形信號存儲的格式為離散時間信號，所以對于小波變換主要討論的是離散小波變換。下文先從連續(xù)小波變換分析入手[11-13]。

3.1.1 連續(xù)小波變換

(1)

式中：a為尺度因子；b為平移因子。

(2)

對于任意函數(shù)在某一尺度a、平移點b上的小波變換系數(shù)，實質(zhì)上表征的是在b位置處，時間段aΔt內(nèi)，包含在中心頻率為ω0/a、帶寬為Δω/a的頻窗內(nèi)的頻率分量大小。隨著尺度a的變化，對應(yīng)窗口中心頻率ω*/a及窗口寬度Δω/a也發(fā)生變化[13]。

3.1.2 離散小波變換

在實際情況下，由于f(t)的連續(xù)小波變換系數(shù)Wψf(a,b)的信息量是冗余的，所以一般會對信號f(t)進(jìn)行離散化。具體做法就是將尺度因子a和平移因子b離散化，只取一些特定的值，便可以減少小波變換系數(shù)的冗余度，這樣用計算機分析和處理將會很方便。

(3)

從而式(2)的連續(xù)小波變換變?yōu)槿缡?4)所示的離散小波變換：

(4)

式中：m,n分別為頻率范圍指數(shù)和時間步長變化指數(shù)。

3.2 LabVIEW與Matlab聯(lián)合編程

Matlab軟件更加清晰直觀地表達(dá)對信號的降噪要求，同時結(jié)合LabVIEW圖形化的編輯語言的優(yōu)點實現(xiàn)對脈沖信號的降噪處理。Matlab程序通過LabVIEW調(diào)用Matlab-Script的節(jié)點，實現(xiàn)對脈沖信號的小波降噪處理。

為了選取更適合系統(tǒng)信號來源的降噪處理方式通過巴特沃斯濾波器和切比雪夫濾波器，對Daubenchies小波降噪、Symlets小波降噪、Coiflets小波降噪進(jìn)行對比。切比雪夫濾波如圖5所示。Daubenchies小波降噪如圖6所示。

圖5 切比雪夫濾波圖

圖6 Daubenchies小波降噪圖

對比圖5、圖6，可以看出：小波降噪的效果比濾波降噪效果明顯，因此小波降噪更加適合受噪聲干擾情況下的瞬態(tài)脈沖信號的降噪處理。

4 結(jié)束語

瞬態(tài)脈沖信號的測量和處理在許多領(lǐng)域中具有重要實用價值。本研究不僅通過LabVIEW實現(xiàn)了存儲在工控機內(nèi)的海量脈沖信號時域的測量、頻域的分析，而且通過LabVIEW與Matlab聯(lián)合編程實現(xiàn)了對信號的小波降噪處理。

試驗結(jié)果表明，該信號測量與處理系統(tǒng)具有較高的精確性和實用性，對信號測量的各項參數(shù)和對信號處理的性能均滿足對信號篩選的要求。該信號測量與處理系統(tǒng)的設(shè)計與實現(xiàn)可以為今后的實際工程和教學(xué)研究提供新方法。

[1] 田熙俏.脈沖中子煤質(zhì)分析中脈沖寬度、頻率的優(yōu)化研究[D].長春：東北師范大學(xué),2013.

[2] 張云亮,馮平法,鮑晟.LabVIEW在大數(shù)據(jù)量采集與處理軟件中的應(yīng)用[J].自動化儀表，2012,33(7):19-20.

[3] 陳錫輝,張銀鴻.LabVIEW8.20程序設(shè)計從入門到精通[M].北京:清華大學(xué)出版社，2007:284-306.

[4] ZHAI X W,ZHANG M L,BAI X N,et al.A signal analysis design based on LabVIEW[C]// 2015 International Conference on Cyber-Enabled Distributed Computing and Knowledge Discovery,2015:473-476.[5] 晏建洋,吳建星.基于LabVIEW和MATLAB的礦山微震信號小波分析與研究[J].安全與環(huán)境工程,2016,23(3)：125-128.

[6] KUMAR V,CHANDRASEKHAR N,ALBER S K,et al.Analysis of arc welding process using digital storage oscilloscope[J].Measurement,2016,3(81):1-12.

[7] 黃坤超,張靖悉,周燁.實時寬帶示波器在快沿脈沖測量中的應(yīng)用[J].電訊技術(shù).2013,53(11)：1532-1536.

[8] 杜玉環(huán),郭迎清,薛海東,等.基于LabVIEW的FFT光纖渦輪流量測量系統(tǒng)研究[J].儀表技術(shù)與傳感器,2015(11)：56-58.

[9] 楊小靜.基于LabVIEW和Matlab的圖像去噪研究[D].長沙：湖南師范大學(xué),2014.

[10]張德豐.MATLAB小波分析[M].北京：機械工業(yè)出版社，2009：49-97.

[11]羅淼,姚運萍.信號降噪中小波基與閾值的選取研究[J].機械制造,2015,53(12):42-45.

[12]ZHANG X P,TIAN L S,PENG Y N,et al.From the wavelet series to the discrete wavelet transform-the initialization [J].IEEE Transactions on Signal Processing,1996,44(1):129-133.

[13]趙秀玉,喬新勇.基于LabVIEW和Matlab混合編程的小波去噪方法[J].現(xiàn)代電子技術(shù)，2013,36(19) :46-48.

DesignoftheMeasurementandProcessingSystemofTransientPulseSignal

SHEN Chunlong1,GUO Feng1,ZHU Chengyin2

(1.School of Information Engineering,Southwest University of Science and Technology,Mianyang 621010,China；2.Institute of Nuclear Physics and Chemistry,China Academy of Engineering Physics,Mianyang 621900,China)

For the test of high intensity pulsed neutron control system,in order to sieve out the transient pulse signals that are fit the requirement of the control of the system,it is necessary to make further measurement and processing for a large number of similar transient pulse signals captured by the digital oscilloscope which is controlled by LabVIEW.Considering that the design of control and collection for transient pulse signals are also controlled by LabVIEW software,while using Matlab alone to analyze and process the signal is difficult to simultaneously satisfy the design demands for measurement and processing of massive amounts of data.Due to LabVIEW software has automatic multithreading technology,graphical design language and rich signal processing module,and Matlab software has advantage in wavelet noise reduction of signals,so the mixed programming of LabVIEW and Matlab is proposed to integrate time domain measurement,frequency domain measurement,and wavelet noise reduction processing,for implementing the design of massive signal analysis and processing system.The results of practical tests verify the correctness and feasibility of the system for processing massive pulsed data,and sieving out useful signals; these also indicate that the system has advantages of simple operation,fast response speed,and high measurement accuracy; thus possesses higher practical value.

Data acquisition; Pulse signal; LabVIEW; Matlab; Time domain measurement; Frequency domain measurement; Finite length sequence; Wavelet noise reduction

TH865；TP274

： A

10.16086/j.cnki.issn1000-0380.201709013

修改稿收到日期：2017-02-07

國家自然科學(xué)基金資助項目(51475453)

申春龍(1990—)，男，在讀碩士研究生，主要從事信號檢測與信號處理方向的研究。E-mail：1415725868@qq.com。 郭鋒(通信作者)，男，教授，主要從事智能儀器等方向的研究。E-mail：guofen9932@163.com。