張世弘，劉振興

(武漢科技大學(xué)信息科學(xué)與工程學(xué)院，湖北武漢 430081)

基于MATLAB－GUI的信號采集與分析系統(tǒng)設(shè)計

張世弘，劉振興

(武漢科技大學(xué)信息科學(xué)與工程學(xué)院，湖北武漢 430081)

在測控以及故障診斷等諸多場合，需要根據(jù)不同的需求對信號采取不同的時頻處理方法，現(xiàn)有的系統(tǒng)信號處理方法比較單一，難以滿足不同場合的需求。從時頻信號處理方法集成的角度出發(fā)，介紹了以MATLAB的圖形用戶界面(GUI)開發(fā)環(huán)境為基礎(chǔ)的信號采集和處理系統(tǒng)的整體設(shè)計方案。進而闡述了從信號采集與生成、分析處理以及數(shù)據(jù)管理3個方面實現(xiàn)該系統(tǒng)的方法。實驗結(jié)果表明:該信號采集與處理系統(tǒng)使用方便、穩(wěn)定性好、具有高精確性。該系統(tǒng)提供良好的交互性和擴展性，能夠用于實際信號分析與仿真實驗。

MATLAB－GUI;串口通信;信號分析集成;可擴展性

0 前言

信號采集和信號分析是測控系統(tǒng)與故障診斷系統(tǒng)的重要組成部分，目前幾乎所有的工程技術(shù)領(lǐng)域都涉及到信號處理。研究信號的時頻結(jié)構(gòu)，是信號處理過程中非常重要的分析手段，現(xiàn)有的產(chǎn)品提供的信號處理方法相對比較單一，有的只限于常規(guī)的頻譜分析，有的提供了少量的幾種方法，不具備擴展功能，只適合特定的專業(yè)應(yīng)用領(lǐng)域。適用范圍。而且普遍價格昂貴，加工工藝復(fù)雜，體積龐大，且生產(chǎn)技術(shù)要求較高［1］。近年來，隨著計算機總線技術(shù)、網(wǎng)絡(luò)技術(shù)、高速數(shù)據(jù)傳輸技術(shù)的快速發(fā)展提高和應(yīng)用，虛擬儀器在測量與分析等眾多領(lǐng)域發(fā)揮了越來越重要的作用［2］。文中充分利用MATLAB在處理數(shù)據(jù)和構(gòu)建人機交互界面上的優(yōu)勢，構(gòu)建虛擬儀器，將數(shù)據(jù)記錄、新的時頻信號處理方法、串口采集上位機等的功能進行集成，搭載多種時頻信號處理算法，實現(xiàn)了集信號生成與通信、信號多角度分析以及數(shù)據(jù)管理等功能于一體的信號采集與分析系統(tǒng)，既能滿足科學(xué)研究的需要，也能作為信號處理方法學(xué)習(xí)的示教系統(tǒng)。

1 MATLAB－GUI開發(fā)環(huán)境

MATLAB－GUI提供Windows基本控件的支持，并且具有良好的事件驅(qū)動機制，同時提供了MATLAB函數(shù)庫的接口，進行系統(tǒng)開發(fā)方便直接。

其設(shè)計圖形界面時有兩種方法:(1)直接通過編寫m文件來創(chuàng)建面向?qū)ο蟮膽?yīng)用程序;(2)使用MATLAB自帶的圖形界面開發(fā)環(huán)境GUIDE(Graph User Interface Developing Environment)實現(xiàn)。開發(fā)者可在GUIDE上直接創(chuàng)建按鈕、文本編輯器等，組建系統(tǒng)的用戶界面，并且通過直接設(shè)置控件的屬性及相應(yīng)代碼，實現(xiàn)各個控件的函數(shù)回調(diào)功能。文中選擇第二種方式。

在該系統(tǒng)的設(shè)計中，將通過模塊化思想利用MATLAB函數(shù)，搭建整體系統(tǒng)框架，繼而建立多個界面分別實現(xiàn)信號采集、分析與數(shù)據(jù)管理功能。最終，使用MATLAB系統(tǒng)LCC編譯器，將本MATLAB－GUI程序轉(zhuǎn)換為exe文件，使系統(tǒng)可以在不安裝MATLAB的環(huán)境中正常運行。

2 系統(tǒng)總體方案

在MATLAB－GUI開發(fā)環(huán)境下設(shè)計信號采集與分析系統(tǒng)需要能夠:

(1)既可以作為實際信號分析系統(tǒng)，也可以進行仿真實驗，用于理論分析;

(2)搭載多種信號處理方法，可以從多個角度對信號進行處理;

(3)提供良好的人機交互仿真環(huán)境，便于用戶針對不同信號選用不同分析方法進行信號檢測與試驗;

(4)提供一定的可擴展性，通過給定接口，方便用戶擴展新算法，新功能。

據(jù)此，將系統(tǒng)分為如下3個模塊:

(1)信號管理模塊。實現(xiàn)對采集、生成以及分析后數(shù)據(jù)的管理，包括數(shù)據(jù)傳輸、存儲等工作

(2)信號采集模塊。實現(xiàn)單通道、多通道信號的實時采集，而且設(shè)計信號發(fā)生器，實現(xiàn)自定義信號的生成;

(3)信號分析處理模塊。實現(xiàn)在線、離線信號分析處理功能，從多個角度對信號的時域、頻域和時頻域特征進行分析與處理，并且提供良好接口，實現(xiàn)新算法 (M文件或C文件)的動態(tài)加載。

信號采集與處理系統(tǒng)整體框架如圖1所示。

圖1 信號采集與處理系統(tǒng)整體框架

3 信號采集與管理設(shè)計

信號采集模塊，分為實時采集部分與信號發(fā)生器部分。實時采集部分負(fù)責(zé)采集實際信號并用作實際分析;信號發(fā)生器部分負(fù)責(zé)產(chǎn)生仿真實驗信號，作為理論分析的實驗信號。信號管理模塊為信號的調(diào)用、保存、新方法的開發(fā)提供保障。

3.1 信號采集功能實現(xiàn)

MATLAB的設(shè)備控制箱 (Instrument Control Toolbox，ICT)提供了對RS232與RS485串口通信的正式支持。MATLAB通過調(diào)用設(shè)備控制工具箱的serial類以及相關(guān)函數(shù)，來創(chuàng)建串口設(shè)備對象，得到設(shè)備的文件句柄，從而以操作文件的方式實現(xiàn)對PC串行口的讀寫操作。其基本步驟包括:創(chuàng)建串口對象并配置串口屬性、連接串口對象和外圍設(shè)備、串口讀寫數(shù)據(jù)、打開關(guān)閉串口以及清除串口對象。

(1)創(chuàng)建串口對象并配置串口屬性

在創(chuàng)建一個串口對象時，利用set函數(shù)來設(shè)置串口對象的屬性如:波特率、數(shù)據(jù)位、奇偶校驗位、停止位和結(jié)束符等。

(2)連接串口對象和外圍設(shè)備

在串口讀寫數(shù)據(jù)前，利用fopen函數(shù)打開串口，與相應(yīng)的設(shè)備連接起來。當(dāng)串口連接到設(shè)備，Status(串口是否連接到設(shè)備)屬性更新為open，輸入緩沖區(qū)和輸出緩沖區(qū)數(shù)據(jù)將清空。

(3)串口讀寫數(shù)據(jù)

串口在進行數(shù)據(jù)讀寫時，通過fread、fwrite完成讀寫操作。串口對象遵循配置的屬性值進行通信，通信數(shù)據(jù)可為二進制或者十六進制表示。

(4)關(guān)閉串口和清除串口對象

利用fclose函數(shù)斷開串口和設(shè)備的連接，使用delete函數(shù)將其從串口對象中清除，后再利用clear函數(shù)從MATLAB工作空間中將其清除。

據(jù)此設(shè)計信號通信流程如圖2所示。

圖2 數(shù)據(jù)通信流程

信號采集完成后，系統(tǒng)針對信號數(shù)據(jù)搭建了分頻系數(shù)——1/n分頻，便于滿足不同頻率分辨率的分析要求。其分頻原理如下:假設(shè)n=2，即1/2分頻，在信號數(shù)據(jù)中按順序從取出其中次序為2、4、6……的數(shù)據(jù)，另保存為新的信號。據(jù)此設(shè)計信號采集界面，實時采集信號過程如圖3所示。

圖3 信號采集界面

3.2 信號發(fā)生器功能實現(xiàn)

為了實現(xiàn)本系統(tǒng)可以進行信號處理仿真實驗的目的，作者設(shè)計了信號發(fā)生器功能。

MATLAB－GUI中的文本編輯器控件可以對編輯框中的內(nèi)容進行編輯、刪除和替換。用戶可以通過其編輯自定義函數(shù)，采樣頻率及信號長度等參數(shù)。進而，系統(tǒng)通過回調(diào)函數(shù)獲取這些信息，生成自定義信號，實現(xiàn)信號發(fā)生器功能。通過函數(shù)的選擇，可以實現(xiàn)任意信號的生成。如圖4所示，使用信號發(fā)生器生成:

其中采樣頻率為1 024 Hz，采樣點數(shù)為1 001的信號。

圖4 信號發(fā)生器界面

3.3 信號管理部分

該系統(tǒng)集成信號采集與分析等多種功能于一體，故而信號的存儲以及加載是系統(tǒng)必不可少的一部分。系統(tǒng)通過調(diào)用MATLAB－GUI開發(fā)環(huán)境中uiput-file和uigetfile函數(shù)，選擇文件名和保存/打開的路徑確定目標(biāo)位置，然后利用save和load函數(shù)完成信號數(shù)據(jù)的保存或者加載。

據(jù)此，在該系統(tǒng)的設(shè)計中，系統(tǒng)將采集、生成的信號和信號分析處理后的每一組數(shù)據(jù)進行保存，且按n×m矩陣形式保存為文本文件或者Excel文件，并且，可以在系統(tǒng)中加載符合以上格式的信號數(shù)據(jù)進行分析。

4 信號分析設(shè)計

信號采集與管理完成后，信號分析同樣是本系統(tǒng)的一個重要環(huán)節(jié)。系統(tǒng)設(shè)計對在線、離線信號從多個角度進行分析，包括時域、頻域、時頻域分析中多種信號處理方法。另外，設(shè)計具備友好的人機界面，集成MATLAB工具欄，方便使用者根據(jù)實際情況，在有限范圍內(nèi)進行適當(dāng)調(diào)整和提取所需的分析數(shù)據(jù)。

4.1 頻域、時域分析功能實現(xiàn)

系統(tǒng)中的信號頻譜分析主要涵蓋時域信號顯示、FFT變換、功率譜變換和IFFT變換等功能，并且集成不同的窗函數(shù)供用戶選擇，可以對信號進行截短分析，以減少頻譜能量泄漏對分析結(jié)果的影響。另外，該系統(tǒng)針對兩組信號間的信號分析，搭載相關(guān)性計算、軸心軌跡計算和相位差計算等功能。用戶通過加載信號數(shù)據(jù)文件，適當(dāng)調(diào)整采樣頻率和采樣點數(shù)等相關(guān)參數(shù)，即可對數(shù)據(jù)進行分析，進而獲得所需要的信息。

據(jù)此設(shè)計頻域、時域分析界面。在此通過信號采集模塊采集到軸承振動信號數(shù)據(jù)，其中采集頻率為10 000 Hz，采集點數(shù)為16 384。對其進行實時時域、頻域特征分析，如圖5、6所示。

圖5 頻域分析界面a

圖6 時域分析界面b

4.2 時頻分析功能實現(xiàn)

(1)極大重疊離散小波包變換分析功能

為了處理多分量的復(fù)雜信號，系統(tǒng)采用極大重疊離散小波包變換 (MODWPT)方法。

通過MODWPT可以將多分量的復(fù)雜信號分解為若干個瞬時頻率和瞬時幅值都具有經(jīng)典物理意義的分量之和，并且求出各個單分量信號的瞬時頻率和瞬時幅值，再進行組合得到原始復(fù)雜信號完整的時頻分布。

(2)Hilbert－Huang變換功能

為了準(zhǔn)確反映信號的固有特性，實現(xiàn)復(fù)雜的非平穩(wěn)信號的合理分析，系統(tǒng)采用了Hilbert－Huang變換(HHT)方法。

Hilbert－Huang變換由經(jīng)驗?zāi)B(tài)分解 (EMD)及Hilbert變換兩部分組成。EMD本質(zhì)是把1個信號進行平穩(wěn)化處理，其結(jié)果是將信號中不同尺度的波動或趨勢逐級分解開來，產(chǎn)生一系列具有不同特征尺度的數(shù)據(jù)序列—本征模態(tài)函數(shù)，再通過希爾伯特變換求得的瞬時頻率及幅值。得到信號振幅—瞬時頻率—時間的分布。

(3)廣義解調(diào)功能

針對分析非平穩(wěn)信號過程中可能產(chǎn)生的人為時頻分割、虛假分量和模態(tài)混淆等問題，并且達到盡可能還原原始信號的時頻分析曲線的目的，本系統(tǒng)還集成了廣義解調(diào)時頻算法。

廣義解調(diào)時頻算法通過選取一種信號變換的方法，將時頻分布為傾斜或非線性的信號變換為時頻分布為線性和平行于時間軸的信號。

據(jù)此，設(shè)計了時頻分析界面，在其中分別設(shè)置了MODWPT窗口、HHT變換窗口和廣義解調(diào)窗口，廣義解調(diào)信號分量提取窗口。并且設(shè)計多分量的復(fù)雜非平穩(wěn)信號:

其中采樣頻率為1 024 Hz，采樣點數(shù)為1 024，對其進行時頻分析，其結(jié)果如圖7所示。

圖7 時頻分析界面

在該界面用戶根據(jù)信號特征設(shè)置采樣頻率，MODWPT分解層數(shù)以及觀測窗口參數(shù)—最低頻率、最高頻率和頻率分辨率，即可對信號進行MODWPT、HHT和廣義解調(diào)分析。

其中，對HHT變換窗口，用戶可以通過設(shè)置觀測本征模態(tài)函數(shù)的層數(shù)，來觀測信號每個分量的特征。由于HHT變換中本征模態(tài)函數(shù)的層數(shù)是根據(jù)信號特征生成的，故而當(dāng)輸入觀測層數(shù)大于實際層數(shù)時，系統(tǒng)自動彈出提示，而且當(dāng)輸入層數(shù)為0時，系統(tǒng)將信號完整的HHT變換結(jié)果呈現(xiàn)出來，在此我們感測第一層本征模態(tài)函數(shù)如圖7左下窗口所示。

另外，系統(tǒng)可以通過設(shè)置信號廣義解調(diào)分量序號，提取廣義解調(diào)后的信號分量進行分析。如圖7中，信號采樣頻率為1 024 Hz，且分解層數(shù)5，故而信號被分解成25個分量，每個分量頻率的跨度=1 024/2×25=16 Hz。而通過逆廣義解調(diào)可以清楚的觀測到信號的兩個分量分別存在于70 Hz和90 Hz附近，即在分量5(64～80 Hz)與分量6(80～96 Hz)中。在此觀測分量5得到了滿意的結(jié)果，如圖7右下窗口所示。

4.3 可擴展功能實現(xiàn)

由于該系統(tǒng)是基于MATLAB－GUI開發(fā)的，提供了強大的數(shù)據(jù)處理功能以及良好拓展功能。據(jù)此作者設(shè)計了如下M文件拓展接口:

［x，y，z］ =functionname(data，con，other1，other2)

x，y，z為系統(tǒng)輸出;data為系統(tǒng)對象，con為約束條件，other1，other2為預(yù)留參數(shù)。用戶通過設(shè)置合理的接口并且通過加載連接算法文件即可實現(xiàn)擴展功能，在使用過程中應(yīng)能夠動態(tài)加載新算法，大大提升了系統(tǒng)的可擴展性。

5 結(jié)束語

利用MATLAB－GUI編程環(huán)境，建立了集信號生成、實時采集、多種信號分析方法和信號管理于一體的信號采集和處理系統(tǒng)，可以便捷地從多種角度對信號進行分析。且隨著信號處理方法的不斷進步，可以方便地擴展新的功能。具有良好的實用性、交互性和可擴展性。在智能化程度、處理能力、性能價格比和可操作性等方面有著明顯的優(yōu)勢。

［1］高新聞.智能儀表軟件開發(fā)的自動化技術(shù)研究［D］.上海:上海大學(xué).2007:95 －109.

［2］秦輝，席裕庚.基于Matlab－GUI的預(yù)測控制仿真平臺設(shè)計［J］.系統(tǒng)仿真學(xué)報，2006，18(10):2278 －2281.

［3］李頎，張建軍，李沙沙，等.基于串行通信的MATLAB GUI在壓力傳感器中的標(biāo)定軟件設(shè)計［J］.計算機測量與控制，2010，18(7):1680 －1683.

［4］王世香.精通MATLAB接口與編程［M］.北京:電子工業(yè)出版社，2007.

［5］張曉菲.廣義解調(diào)在復(fù)合時變信號處理中的應(yīng)用［D］.武漢:武漢科技大學(xué)，2012:88－95.

Design of Signal Acquisition and Analysis System Based on MATLAB－GUI

ZHANG Shihong，LIU Zhenxing
(School of Information Science and Engineering，Wuhan University of Science and Technology，Wuhan Hubei430081，China)

In manuy cases of test，monitor and fault dignosis，it is in need of using different time frequency processingmethod to signals according to different demands.The present processingmethod of system signal is rather simple，which is difficult tomeet demands in different conditions.Starting from the time frequency signal acquisition and processingmethod，the overall design schemewas introduced of the system of signal acquisition and processing based on the development environment of graphical user interface(GUI)program ofMATLAB －GUI.Furthermore，the system realizingmethod was presented in three aspects of signal acquisition and generation，processing and analysis and datamanagement.The experimental results indicate that the system runs rapidly and stably，aswell as the signal data are collected and analyzed promptly and accurately.The system provides a good interactive and extensibility，which can be used in the actual signal analysis and simulation experiments.

MATLAB－GUI;Serial communication;Signal acquisition and analysis;Scalability

TP274.2

A

1001－3881(2015)21－113－5

10.3969/j.issn.1001 －3881.2015.21.027

2014－09－18

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

張世弘 (1990—)，男，碩士研究生，研究方向為工業(yè)過程控制、電氣傳動及電機故障診斷。E－mail:595932627@qq.com。