激光多普勒差頻條紋處理電路設(shè)計(jì)

王　偉,季園媛,張雄星,劉光海

(1.西安工業(yè)大學(xué) 光電工程學(xué)院,西安 710021;2.四川物科光學(xué)精密機(jī)械有限公司,綿陽 621900)

?

激光多普勒差頻條紋處理電路設(shè)計(jì)

王偉1,季園媛1,張雄星1,劉光海2

(1.西安工業(yè)大學(xué) 光電工程學(xué)院,西安 710021;2.四川物科光學(xué)精密機(jī)械有限公司,綿陽 621900)

摘要：為了解決全光纖激光速度干涉系統(tǒng)產(chǎn)生的多普勒信號(hào)的實(shí)時(shí)采集和處理問題,提出了一種以現(xiàn)場(chǎng)可編程邏輯門陣列(FPGA)為核心的信號(hào)采集和處理系統(tǒng)的實(shí)現(xiàn)方案.系統(tǒng)以Cyclone IV的FPGA為核心控制模塊,結(jié)合信號(hào)調(diào)理電路、A/D電路及D/A轉(zhuǎn)換電路構(gòu)成多普勒信號(hào)的采集和處理系統(tǒng),完成針對(duì)振動(dòng)信號(hào)的速度和位移信息的解調(diào).采用音叉搭建了振動(dòng)測(cè)試平臺(tái),通過示波器可同時(shí)觀測(cè)目標(biāo)的速度和位移變化.研究結(jié)果表明：所設(shè)計(jì)的測(cè)振系統(tǒng)速度測(cè)量范圍為0.1～30 m·s-1,頻率測(cè)量范圍為DC-1 MHz,頻率測(cè)量誤差為1%.

關(guān)鍵詞：激光干涉；現(xiàn)場(chǎng)可編程門陣列；信號(hào)采集；頻率測(cè)量

frequency measurement

測(cè)振技術(shù)目前分為兩大類接觸式測(cè)量和非接觸測(cè)量.接觸式測(cè)量是通過在測(cè)量物體上安裝壓電傳感器,必須依附在被測(cè)物體表面.在特殊的振動(dòng)測(cè)試中,由于傳感器大小的局限,會(huì)破壞原有的振動(dòng)狀態(tài),空間較小可能無法安裝,從而增加測(cè)試難度.激光多普勒測(cè)振[1]是利用光散射原理和多普勒效應(yīng),根據(jù)激光的相干性,通過獲得目標(biāo)物體速度變化引起的干涉條紋變化.光波比無線電波的波長(zhǎng)短、頻率高,以其頻移和波長(zhǎng)作為尺子去度量被測(cè)目標(biāo)的速度和位移,精度高.但光的多普勒頻移一般在幾十MHz量級(jí),光電探測(cè)器和數(shù)據(jù)采集電路的帶寬相應(yīng)的就需達(dá)到百M(fèi)Hz量級(jí),這對(duì)后續(xù)的光電探測(cè)和信號(hào)采集電路提出了較高要求,因此,激光多普勒信號(hào)的采集和處理是激光多普勒測(cè)量系統(tǒng)工程實(shí)現(xiàn)上的難點(diǎn)之一.文獻(xiàn)[2]提出在信號(hào)處理過程中,通過高速數(shù)據(jù)采集卡把多普勒信號(hào)采集到計(jì)算機(jī)中,借助Matlab平臺(tái),選取時(shí)頻分析中的Wigner_Ville分布對(duì)采集進(jìn)入的多普勒信號(hào)進(jìn)行分析.這種方式是在下位機(jī)中進(jìn)行采集,在上位機(jī)中解算,實(shí)現(xiàn)數(shù)據(jù)的離線處理；文獻(xiàn)[3]設(shè)計(jì)了基于數(shù)字信號(hào)處理器(Digitial Signal Processor,DSP)的數(shù)據(jù)采集和處理系統(tǒng),系統(tǒng)采用AD芯片采樣,片外隨機(jī)存取存儲(chǔ)器(Random Access Memory,RAM)緩存,液晶屏顯示的方式采用快速傅氏變換(Fast Fourier Transformation,FFT)算法進(jìn)行數(shù)據(jù)處理,這種處理算法復(fù)雜,計(jì)算量大；文獻(xiàn)[4]設(shè)計(jì)了以現(xiàn)場(chǎng)可編程邏輯門陣列(Field Programmable Gate Array,FPGA)為核心的處理模塊,在光路部分通過正交分解的方式產(chǎn)生兩路正交信號(hào)采用模數(shù)轉(zhuǎn)換芯片采集,經(jīng)過FPGA預(yù)處理后通過RS-232的方式傳輸上位機(jī)；文獻(xiàn)[5]研究了基于FPGA的直接測(cè)頻法的測(cè)頻精度分析,提出了采用直接測(cè)頻法誤差來源于開始階段采樣的量化誤差,以及終止采樣階段的量化誤差；文獻(xiàn)[6-7]提出了目前常見的方式主要是通過高速通信接口將數(shù)據(jù)傳遞到上位機(jī)中進(jìn)行離線數(shù)據(jù)處理.

基于此,本文從優(yōu)化數(shù)據(jù)處理算法實(shí)現(xiàn)在下位機(jī)中數(shù)據(jù)處理,實(shí)現(xiàn)實(shí)時(shí)觀測(cè)的角度進(jìn)行了研究.設(shè)計(jì)了一種針對(duì)多普勒信號(hào)的數(shù)據(jù)采集和處理電路.系統(tǒng)以FPGA為核心處理芯片,采用多周期同步測(cè)頻算法,實(shí)現(xiàn)對(duì)采集的多普勒信號(hào)解調(diào)處理,從而得到物體表面的速度和位移,通過模擬信號(hào)輸出方式.該方案將數(shù)據(jù)采集、信號(hào)處理算法以及輸出集成在一片電路中,實(shí)現(xiàn)了振動(dòng)信號(hào)的實(shí)時(shí)連續(xù)觀測(cè).搭建了實(shí)驗(yàn)測(cè)試平臺(tái),通過對(duì)不同頻率的音叉進(jìn)行測(cè)量,驗(yàn)證了本系統(tǒng)的可靠性和精度.

1硬件電路設(shè)計(jì)

全光纖多普勒測(cè)量系統(tǒng)輸出的光差頻信號(hào)經(jīng)光電探測(cè)器實(shí)現(xiàn)光電轉(zhuǎn)換,轉(zhuǎn)換后的電信號(hào)通過所設(shè)計(jì)的數(shù)據(jù)采集和處理板,實(shí)現(xiàn)信號(hào)的采集處理.所設(shè)計(jì)的系統(tǒng)框圖如圖1所示.電信號(hào)經(jīng)電流電壓轉(zhuǎn)換電路、經(jīng)濾波及放大電路預(yù)處理后送給A/D轉(zhuǎn)換成數(shù)字信號(hào),送入FPGA中.由FPGA中的硬件邏輯電路對(duì)數(shù)據(jù)進(jìn)行緩存,在FPGA的邏輯控制下,將解算出的速度及方向變化量通過模擬量的形式輸出.本設(shè)計(jì)的系統(tǒng)主要模塊包含了FPGA的最小系統(tǒng)、電源模塊、信號(hào)調(diào)理電路、A/D轉(zhuǎn)換電路及D/A轉(zhuǎn)換電路等部分.

圖1　數(shù)據(jù)采集和處理系統(tǒng)原理框圖

1.1信號(hào)調(diào)理電路設(shè)計(jì)

在本設(shè)計(jì)中,光路測(cè)量部分采用1 550 nm的分布反饋式激光器(Distributed Feedback Laser,DFB-LD),采用光外差檢測(cè)[8]輸出時(shí),光電探測(cè)器接收到的光信號(hào)功率僅為微瓦(μW)的數(shù)量級(jí),以至于輸出的多普勒信號(hào)非常微弱,噪聲成分很大,且具有較大的直流偏置.該信號(hào)很容易被其他噪聲干擾.為了提高信噪比,保證后續(xù)信號(hào)處理的準(zhǔn)確性,必須對(duì)光電探測(cè)器輸出的信號(hào)進(jìn)行去偏置、放大和濾波.本文采用帶通濾波器的設(shè)計(jì)方案,可以實(shí)現(xiàn)在特定頻段的頻率傳遞信號(hào),衰減這一頻段以外的所有信號(hào).由于多普勒信號(hào)的頻率較高,可達(dá)40 MHz,所以本文選擇了亞德諾半導(dǎo)體(Analog Devices,Inc.,ADI)公司的高速、高頻的運(yùn)放芯片AD8045,其具有低噪聲、電壓反饋和高速放大等特性,滿足本設(shè)計(jì)的需求.

如圖2所示,輸入信號(hào)為光電探測(cè)器輸出的包含光差頻信息的多普勒信號(hào).R1、C2構(gòu)成低通濾波器將高頻噪聲阻斷,R3、C1構(gòu)成高通濾波器將低頻噪聲阻斷,兩者串聯(lián)接在運(yùn)放的同相輸入端,通過反饋電阻R5來改善濾波特性.通過Pin6管腳將信號(hào)輸出到后續(xù)電路.

1.2A/D轉(zhuǎn)換電路設(shè)計(jì)

多普勒信號(hào)經(jīng)調(diào)理電路后,送入A/D轉(zhuǎn)換電路轉(zhuǎn)換成數(shù)字信號(hào).根據(jù)系統(tǒng)設(shè)定的測(cè)速范圍在0.1～30 m·s-1,多普勒的頻移fd表達(dá)式為

fd=2v/λ

(1)

式中：λ為激光波長(zhǎng)；v為振動(dòng)目標(biāo)的瞬時(shí)速度.

計(jì)算可得信號(hào)頻率最高不超過40 MHz,依據(jù)采樣定理知信號(hào)采樣速率至少為80 MHz.本系統(tǒng)選擇德州儀器公司的8 位高速模數(shù)轉(zhuǎn)換芯片ADC08200,該芯片的最高采樣速率可達(dá)200 MHz.電路原路圖如圖3所示.通過外接參考源ADR441ARZ芯片來產(chǎn)生2.5 V的電壓作為電壓參考源接入芯片Pin3管腳,其輸出的8位數(shù)據(jù)總線D0-D7和FPGA的I/O相連接,將數(shù)據(jù)傳輸?shù)紽PGA中.

圖2　信號(hào)調(diào)理電路設(shè)計(jì)

圖3　A/D轉(zhuǎn)換電路設(shè)計(jì)

1.3FPGA設(shè)計(jì)

多普勒信號(hào)經(jīng)A/D轉(zhuǎn)換后變?yōu)閿?shù)字信號(hào)傳送給FPGA,在FPGA中進(jìn)行數(shù)據(jù)的解算.本系統(tǒng)采用Altera公司的EP4CE15F17C8N,它具有15 408個(gè)邏輯單元.FPGA是一種多電源需求的芯片,根據(jù)本系統(tǒng)設(shè)計(jì)要求,系統(tǒng)采用±5 V外接電源供電,滿足前端信號(hào)調(diào)理電路和A/D轉(zhuǎn)換電路的需求,選用多塊電壓轉(zhuǎn)換芯片提供FPGA所需要的3.3 V、2.5 V和1.2 V的供電需求.所設(shè)計(jì)的FPGA數(shù)字邏輯模塊如圖4所示.

圖4　FPGA數(shù)字邏輯框圖

圖4包含了時(shí)鐘模塊、A/D邏輯模塊、數(shù)據(jù)處理模塊、FIFO緩存模塊、D/A邏輯模塊以及聯(lián)合測(cè)試工作組(Joint Test Action Group,JTAG)下載模塊.在工作過程中,通過時(shí)鐘模塊提供系統(tǒng)所需的工作時(shí)鐘.FPGA首先通過A/D控制模塊檢測(cè)是否接受到多普勒數(shù)據(jù)采集完畢的使能信號(hào),如果收到就讀取數(shù)據(jù)并采用異步FIFO緩存.當(dāng)FIFO接近全滿的狀態(tài)時(shí)啟動(dòng),在信號(hào)處理邏輯模塊進(jìn)行頻率測(cè)量.將解算出的速度信號(hào)和方向信號(hào)通過D/A控制模塊輸出,從而完成整個(gè)采集處理邏輯控制.

1.4D/A轉(zhuǎn)換電路設(shè)計(jì)

通過FPGA解算得到被測(cè)振動(dòng)物體表面的速度變換和方向變換,通過兩路D/A轉(zhuǎn)換模塊以在模擬信號(hào)的輸出方式,通過示波器可實(shí)時(shí)觀測(cè).由于系統(tǒng)是高速電路,在D/A的選型時(shí)需要考慮D/A芯片的轉(zhuǎn)換速率.本系統(tǒng)選用的是8位TLC5602C模數(shù)轉(zhuǎn)換芯片,常用于高速信號(hào)處理中,取樣速度最高可達(dá)20 MHz,可以滿足系統(tǒng)對(duì)轉(zhuǎn)換速率的需求,其設(shè)計(jì)電路如圖5所示.兩路D/A的數(shù)據(jù)總線D0～D7和時(shí)鐘控制口分別和FPGA的I/O相連,轉(zhuǎn)換結(jié)果通過芯片的Pin6輸出.

圖5　D/A轉(zhuǎn)換電路設(shè)計(jì)

2多周期同步頻率測(cè)量算法

A/D轉(zhuǎn)換模塊轉(zhuǎn)換得到數(shù)字量,傳輸給FIFO模塊,進(jìn)行數(shù)據(jù)的緩存.采用異步FIFO,在數(shù)據(jù)暫存的同時(shí),讀取數(shù)據(jù)進(jìn)行信號(hào)頻率測(cè)量.采用的多周期同步頻率測(cè)量[9-10]方法以直接測(cè)頻法[6]為基礎(chǔ),該方法保證了閘門信號(hào)與被測(cè)信號(hào)的同步,實(shí)現(xiàn)了全頻段等精度測(cè)量,測(cè)量原理如圖6所示.

圖6　多周期同步頻率測(cè)量原理圖

在系統(tǒng)的運(yùn)行過程中,首先開啟預(yù)設(shè)閘門信號(hào),所設(shè)的兩個(gè)計(jì)數(shù)器未開始計(jì)數(shù),需要等待被測(cè)信號(hào)fx的上升沿來作為觸發(fā)信號(hào),這時(shí)系統(tǒng)控制對(duì)應(yīng)的D觸發(fā)器電路,觸發(fā)器發(fā)出允許系統(tǒng)開始計(jì)數(shù)的信號(hào),計(jì)數(shù)器1開始對(duì)系統(tǒng)頻率信號(hào)f0的脈沖開始計(jì)數(shù),計(jì)數(shù)器2開始對(duì)被測(cè)信號(hào)fx的脈沖開始計(jì)數(shù).在關(guān)閉閘門信號(hào)時(shí),系統(tǒng)中的兩個(gè)計(jì)數(shù)器不停止計(jì)數(shù),需等到fx的上升沿作為終止信號(hào),才停止計(jì)數(shù).完成整個(gè)測(cè)頻過程.

圖6中計(jì)數(shù)器1所計(jì)得的脈沖個(gè)數(shù)為N0,計(jì)數(shù)器2所記得脈沖個(gè)數(shù)為Nx,由公式

N0·T0=Nx·Tx

(2)

可得被測(cè)信號(hào)的頻率為

fx=Nx/N0·f0

(3)

3實(shí)驗(yàn)結(jié)果及分析

3.1系統(tǒng)性能驗(yàn)證

采用音叉搭建振動(dòng)測(cè)試平臺(tái).實(shí)驗(yàn)所用儀器設(shè)備包括激光器、分束器、線性電源、示波器、音叉、數(shù)據(jù)采集和處理電路板等.采用Tektronix 數(shù)字示波器實(shí)行波形監(jiān)測(cè),可以實(shí)時(shí)觀測(cè)到被測(cè)對(duì)象的速度和位移變化.實(shí)驗(yàn)結(jié)果如圖7所示.

圖7中,通道1是光電探測(cè)器輸出的原始多普勒信號(hào),通道2為D/A輸出的被測(cè)音叉的瞬時(shí)速度變化,通道3是D/A輸出的被測(cè)音叉的位移方向變化.實(shí)驗(yàn)采用的音叉本征頻率為256 Hz,測(cè)振系統(tǒng)輸出的速度信號(hào)的頻率為256.1 Hz,實(shí)驗(yàn)結(jié)果表明系統(tǒng)誤差為0.04%.

圖7　音叉測(cè)試效果圖

實(shí)驗(yàn)采用了9組不同頻率(128～1 024 Hz)的本征音叉進(jìn)行驗(yàn)證,如圖8所示,由此結(jié)果可以看出,實(shí)際測(cè)量頻率和理論測(cè)量頻率誤差在1 Hz的范圍內(nèi),測(cè)量結(jié)果顯示該裝置對(duì)低頻段的振動(dòng)頻率測(cè)量準(zhǔn)確度很高,誤差小。

圖8　實(shí)測(cè)信號(hào)和本征音叉的對(duì)比

3.2誤差分析

算法誤差[13-15]表達(dá)式為fd=2v/λ,可以看出誤差主要來自fd和λ,通過誤差傳遞公式ΔV=λ/2×Δfd+fd/2×Δλ產(chǎn)生實(shí)驗(yàn)數(shù)據(jù)誤差,其中Δfd為真實(shí)頻率與測(cè)的頻率的差值,Δλ為標(biāo)稱波長(zhǎng)與實(shí)際波長(zhǎng)的差值.

Δfd的差值主要來源于算法誤差.本設(shè)計(jì)中采用的多周期同步測(cè)頻法,通過信號(hào)發(fā)生器產(chǎn)生不同頻率的方波,通過QuartusⅡ軟件的SignalTap在線監(jiān)控功能,所觀測(cè)到的測(cè)量頻率和實(shí)際頻率值對(duì)比見表1.

表1　頻率測(cè)量實(shí)際值和測(cè)量值對(duì)比

Δλ的誤差主要來源于DFB激光器頻率漂移產(chǎn)生的誤差,DFB激光器最大特點(diǎn)是具有較好的單色性(即光譜純度),其線寬普遍可以做到1 MHz以內(nèi),且具有較高的邊模抑制比(Side-Mode Suppression Ratio,SMSR),因此可認(rèn)為頻率漂移產(chǎn)生的誤差可以忽略不計(jì).

從上述分析可以看出,該裝置對(duì)振動(dòng)頻率測(cè)量準(zhǔn)確度很高,線性度較好.實(shí)驗(yàn)結(jié)果表明,所設(shè)計(jì)的信號(hào)采集和處理系統(tǒng)的精度較高,整體誤差不超過1%,可應(yīng)用于實(shí)際的工程中.

4結(jié) 論

1) 本文設(shè)計(jì)了全光纖激光多普勒測(cè)振儀的信號(hào)采集與實(shí)時(shí)處理系統(tǒng).解決了針對(duì)于多普勒信號(hào)本身具有較大偏置和噪聲的非穩(wěn)態(tài)信號(hào)的采集和在FPGA中進(jìn)行數(shù)據(jù)處理,對(duì)振動(dòng)信號(hào)的瞬時(shí)速率進(jìn)行實(shí)時(shí)輸出的問題.

2) 采用多周期同步測(cè)頻法,采用Verilog HDL編寫時(shí)序代碼,及各模塊之間的邏輯控制.通過仿真驗(yàn)證了整體誤差不超過1%,具有較高的精度.

3) 通過音叉實(shí)驗(yàn)驗(yàn)證整個(gè)系統(tǒng)的正確性,可以滿足設(shè)計(jì)需求,對(duì)振動(dòng)信號(hào)的實(shí)時(shí)解算和處理有實(shí)際的應(yīng)用價(jià)值.

參 考 文 獻(xiàn)：

[1]呂宏詩,劉彬.激光多普勒測(cè)振技術(shù)的最新進(jìn)展[J].激光技術(shù),2005,2(4):176.

LYU Hongshi,LIU Bin.Latest Development of Laser Doppler Technique in Vibration Measurement[J].Laser Technology,2005,2(4):176.(in Chinese)

[2]鐘瑩,張國(guó)雄,張濤,等.時(shí)頻分析在激光多普勒信號(hào)處理中的應(yīng)用[J].光學(xué)技術(shù),2007,33(2):281.

ZHONG Ying,ZHANG Guoxiong,ZHANG Tao,et al.Application of Time -frequency Analysis in Laser Doppler Signal Processing[J].Optical Technique,2007,33(2):281.(in Chinese)

[3]岳邦強(qiáng),張?jiān)?王鳴.基于DSP的激光自混合干涉測(cè)振儀信號(hào)采集與處理系統(tǒng)[J].光學(xué)與光電技術(shù),2009(2):44.

YUE Bangqiang,ZHANG Yuan,WANG Ming.Based on DSP of the Laser Self Mixing Interference Vibrometer Signal Acquisition and Processing System[J].Optical & Optoelectronic Technology,2009(2):44.(in Chinese)

[4]陳勇,錢劍敏,尚建華,等.基于FPGA的多普勒測(cè)振計(jì)信號(hào)采集與處理系統(tǒng)設(shè)計(jì)[J].現(xiàn)代電子技術(shù),2013,36(17):125.

CHEN Yong,QIAN Jianmin,SHANG Jianhua,et al.Design of FPGA-based Signal Acquisition and Processing System for Doppler Vibrometer[J].Modern Electronics Technique,2013,36(17):125.

(in Chinese)

[5]韓艷峰.直接計(jì)數(shù)法瞬時(shí)測(cè)頻的誤差分析[J].測(cè)控技術(shù),2014,33(9):154.

HAN Yanfeng.Error Analysis of Direct Counting Method of Instantaneous Frequency Measurement[J].Measurement and Control Technology,2014,33(9):154.

(in Chinese)

[6]高銳.基于FPGA的多周期同步頻率計(jì)設(shè)計(jì)[D].長(zhǎng)春：吉林大學(xué),2013.

GAO Rui.Design of Synchronous Multi-cycle Frequency Meter Based on FPGA[D].Changchun：Jilin University,2013.(in Chinese)

[7]范亞斌.基于FPGA的頻率計(jì)的設(shè)計(jì)[D].石家莊：河北科技大學(xué),2009.

FAN Yabin.Design on Digital Frequency Meter Based on FPGA[D].Shijiazhuang：Hebei University of Science and Technology,2009.(in Chinese)

[8]王德田.光纖多普勒差拍干涉測(cè)速技術(shù)研究[D].綿陽：中國(guó)工程物理研究院,2008.

WANG Detian.Research of the Fiber Doppler Velocimeter[D].Mianyang:China Academy of Engnieering Physics,2008.(in Chinese)

[9]楊世錫,梁文軍,于保華.振動(dòng)信號(hào)多通道同步整周期數(shù)據(jù)采集卡設(shè)計(jì)[J].振動(dòng)、測(cè)試與診斷,2013,33(1):23.

YANG Shixi,LIANG Wenjun,YU Baohua.Design of Vibration Signal’s Multi-channel Integer Period Synchronous Acquisition Card[J].Journal of Vibration,Measurement & Diagnosis,2013,33(1):23.

(in Chinese)

[10]魏西峰.全同步數(shù)字頻率測(cè)量方法的研究[J].現(xiàn)代電子技術(shù),2005,28(12):101.

WEI Xifeng.Complete Synchronization Digital Frequency Measure Method’s Research[J].Modern Electronics Technique,2005,28(12):101.

(in Chinese)

(責(zé)任編輯、校對(duì)潘秋岑)

Circuit of Laser Doppler Difference Frequency Fringe Design

WANGWei1,JIYuanyuan1,ZHANGXiongxing1,LIUGuanghai2

(1.School of Optoelectronic Engineering,Xi’an Technological University,Xi’an 710021,China;2.Sichuan P&S Optical and Fine Mechanics Co.,Ltd.,Mianyang 621900,China)

Abstract:In order to solve the problem of real time acquisition and processing of the Doppler signal generated by the whole fiber laser velocity interferometer system,a method of signal acquisition and processing system based on FPGA is proposed.The system is constituted with FPGA of cyclone IV as the core control module,combined with the signal conditioning circuit,circuit and A/D conversion circuit of A/D a Doppler signal acquisition and processing system,according to the velocity and displacement of vibration signal demodulation.The fork is used to build the vibration test platform,and the velocity and displacement changes are observed by the oscilloscope at the same time.The results show that the the range of velocity measurement of vibration system is in the 0.1～30 m·s-1,frequency measurement range is DC to 1 MHz,and the frequency measurement error is 1%.

Key words:laser interference; field programmable gate array(FPGA);signal acquisition;

DOI:10.16185/j.jxatu.edu.cn.2016.04.004

收稿日期：2015-09-28

作者簡(jiǎn)介：王偉(1973-),男,西安工業(yè)大學(xué)副教授,主要研究方向?yàn)槲C(jī)控制技術(shù).E-mail:15005804@qq.com.

文獻(xiàn)標(biāo)志碼:中圖號(hào)：TM932A

文章編號(hào)：1673-9965(2016)04-0277-06