何睿 胡定玉 于淼 李曉龍
摘 要: 聲源可視化技術(shù)對傳聲器數(shù)量和采集系統(tǒng)的性能具有一定的要求,傳統(tǒng)駐極體傳聲器陣列采集系統(tǒng)存在功耗高、體積大、采集系統(tǒng)可移植性差和拓展性差等問題。本文設(shè)計了36陣元數(shù)字MEMS傳聲器陣列成像系統(tǒng),有效克服了以上問題。采集模塊選用FPGA為硬件平臺,完成數(shù)據(jù)采集和信號處理模塊設(shè)計,分析了I2S音頻傳輸協(xié)議的實現(xiàn)、數(shù)據(jù)傳輸機制的建立和上位機數(shù)據(jù)的交互方式。并在普通室內(nèi)環(huán)境中對此系統(tǒng)進行實驗驗證,系統(tǒng)的聲源成像效果良好、成像分辨率高,在實際的噪聲源定位與噪聲設(shè)備監(jiān)測中有廣泛的應(yīng)用前景。
關(guān)鍵詞: 聲源可視化; 波束形成; FPGA; MEMS聲傳感器
文章編號: 2095-2163(2021)07-0222-06中圖分類號:TP212,TB52+5文獻(xiàn)標(biāo)志碼: A
Design of microphone array imaging system based on FPGA
HE Rui1,? HU Dingyu1,2, YU Miao3, LI Xiaolong3
(1 School of Urban Railway Transportation, Shanghai University of Engineering Science, Shanghai 201620, China; 2? Shanghai
Engineering Research Center of Vibration and Noise Control Technologies for Rail Transit, Shanghai University of Engineering
Science, Shanghai 201620, China; 3 China Railway Harbin Bureau Group Co., Ltd. Harbin EMU, Harbin 150000, China)
【Abstract】Sound source visualization techniques have some requirements on the number of microphones and the performance of the data acquisition system. Traditional electret microphone array data acquisition systems have problems such as high-power consumption, large size, poor portability and poor expandability of the data acquisition system. In this paper, a 36-array digital MEMS microphone array imaging system is designed to effectively overcome the above problems. The FPGA is selected as the hardware platform for the data acquisition module, and the data acquisition and signal processing modules are designed. The implementation of I2S audio transmission protocol, the establishment of data transmission mechanism and the interaction mode of the upper computer data are analyzed. The system has been experimentally verified in a common environment, and the system has good sound source imaging effect and high imaging resolution, which has a wide application prospect in practical noise source localization and noise equipment monitoring.
【Key words】sound source visualization; beamforming; FPGA; MEMS microphone
0 引 言
聲源成像可視化技術(shù)[1]利用傳聲器陣列采集多通道聲壓數(shù)據(jù),通過信號處理算法,將聲信號轉(zhuǎn)換成圖像信號使得空間聲源分布可視化,可快速定位聲源位置、分布情況以及瞬時動態(tài)變化特征,以用于設(shè)備狀態(tài)評估和監(jiān)測。目前,傳聲器陣列采集系統(tǒng)可分為2種類型。一種是由數(shù)據(jù)采集卡、標(biāo)準(zhǔn)總線和計算機構(gòu)成,主要用于工業(yè)測試環(huán)境。其中,以美國國家儀器(NI)、泰克(Tektronix)、安捷倫(Agilent)為首的公司生產(chǎn)的儀器采集性能優(yōu)越,工作狀態(tài)穩(wěn)定,但該類型采集儀器具有高功耗、大體積的缺點[2],且需要配合相應(yīng)的數(shù)據(jù)采集卡使用,具有專一性強、可移植性差等局限性。另一種是由微處理器為采集核心的嵌入式采集系統(tǒng)構(gòu)成,主要用于實驗測試和部分工業(yè)測試。根據(jù)核心控制器類型不同,可選用單片機(MCU)、DSP(Digital Signal Processing)、ARM(Advanced RISC Machine)和FPGA(Field Programmable Gate Array)四種方式設(shè)計數(shù)據(jù)采集系統(tǒng)。其中,MCU一般為低速時鐘,采集效率不高,難以實現(xiàn)多通道數(shù)據(jù)采集控制。DSP適合于數(shù)字信號處理,難以處理多通道數(shù)據(jù)同步采集邏輯控制[3]。ARM雖然功能比較豐富,但其運算資源有限,同時也不適合對時鐘要求較高的場合。而基于FPGA的數(shù)據(jù)采集系統(tǒng)能并行處理和同步采集多通道信號,邏輯功能強大,接口豐富且使用靈活[4],便于大幅度拓展傳聲器采集通道數(shù)量,滿足其聲源成像可視化技術(shù)采集系統(tǒng)需求。
MEMS(Micro Electro-Mechanical Systems)技術(shù)是從半導(dǎo)體制造工藝中衍生出來的新型加工技術(shù),將傳感器、信號調(diào)理電路以及模數(shù)轉(zhuǎn)換器集成封裝[5]。自2006年Knowles Acoustics推出MEMS傳聲器以來,MEMS傳聲器逐漸取代傳統(tǒng)駐極體傳感器成為傳聲器陣列的默認(rèn)選擇,在噪聲源定位[6]、運動車輛狀態(tài)檢測[7]、語音信號處理[8]和聲學(xué)成像等場合廣泛應(yīng)用。Jelmer等人[9]利用數(shù)字MEMS傳聲器設(shè)計了一款52陣元圓形陣列,用于開闊地帶噪聲源定位,在25 m2消聲室測試中,定位精度達(dá)到厘米級,而在模擬的10 000 m2開闊場地中,其定位精度可達(dá)到5個寬帶聲源。Izquierdo等人[10]指出聲學(xué)成像系統(tǒng)設(shè)計取決于傳聲器數(shù)量的選擇和處理算法的集成度。對于聲學(xué)成像系統(tǒng)需要多通道信號采集的需求,該系統(tǒng)必須是模塊化和可拓展性的,便于構(gòu)建多維度陣列和增加陣列的孔徑。研究中通過選取MyRIO作為硬件平臺(帶有FPGA和ARM處理器的嵌入式硬件),提出了一種可拓展的聲學(xué)成像系統(tǒng)。通過研究生物識別領(lǐng)域的案例,驗證了該系統(tǒng)的可行性。陳松林等人[11]以260個陣元設(shè)計的嵌套式數(shù)字MEMS傳聲器陣列成像系統(tǒng),在普通車間環(huán)境下對機械噪聲和氣體泄漏噪聲進行成像測試,其實驗結(jié)果表明,該系統(tǒng)具有抗干擾性強,成像分辨率高的特點。
根據(jù)前述的系統(tǒng)開發(fā)探討研究,本文選用數(shù)字式MEMS傳聲器設(shè)計傳聲器陣列采集系統(tǒng),選用FPGA作為硬件采集平臺,并設(shè)計了36陣元Annular陣列,可實現(xiàn)多路數(shù)據(jù)的同步采編存儲。
1 系統(tǒng)總體架構(gòu)
數(shù)字MEMS傳聲器陣列并行采集處理系統(tǒng)如圖1所示。由圖1可知,基于FPGA的傳聲器陣列成像系統(tǒng)由數(shù)字MEMS傳聲器陣列、Zynq-7000和上位機組成。其中,傳聲器選用Knowles型號為SPH0645LM4H-B的數(shù)字MEMS傳聲器,最高采樣率可達(dá)到64 kHz,采樣位數(shù)為24 bit,數(shù)據(jù)精度為18 bit。硬件系統(tǒng)處理核心選用的FPGA芯片型號為Zynq-7000,該芯片由可編程邏輯(Programmable Logic, PL)和雙核ARM Cortex-9為核心構(gòu)成的處理系統(tǒng)(Processing System, PS)組成。
數(shù)字MEMS傳聲器陣列并行采集系統(tǒng)流程參見圖1。首先,聲壓信號通過MEMS傳聲器陣列采集,經(jīng)過通用I/O管腳送入Zynq-7000的PL端。接下來,I2S傳輸協(xié)議將音頻信號分成左右兩個通道進行采集,I2S傳輸協(xié)議可控制2個MEMS傳感器。本采集系統(tǒng)利用此傳輸特性,將其設(shè)為一組,分別將數(shù)據(jù)在一級FIFO中進行緩存。然后,通過設(shè)定數(shù)據(jù)輪詢機制,將各一級緩存FIFO中存儲的數(shù)據(jù)通過串并轉(zhuǎn)換緩存至二級FIFO中。PS端與PL端經(jīng)AXI總線進行傳輸通信,選用AXI-HP (High Performance Ports)高速通信接口,帶有FIFO緩沖來提供“批量”讀寫操作,以支持PL端和PS端中存儲器單元的高速率通信。在ARM中的PS端通過BRAM控制器讀取PL端緩存的數(shù)據(jù)。利用PL端的按鍵控制,實現(xiàn)采集的開始和數(shù)據(jù)的發(fā)送。最后,數(shù)據(jù)發(fā)送通過UART串口,將數(shù)據(jù)傳輸給上位機,以供波束形成等信號后處理算法進行分析。
2 系統(tǒng)模塊設(shè)計
2.1 傳聲器選型
MEMS傳聲器根據(jù)信號輸出類型可分為模擬型輸出和數(shù)字型輸出,模擬MEMS傳聲器是將聲壓信號轉(zhuǎn)換為相應(yīng)電壓信號輸出,數(shù)字MEMS傳聲器則是利用PDM(Pluse Density Modulation)和I2S傳輸協(xié)議輸出數(shù)字信號。與傳統(tǒng)駐極體式傳聲器相比,數(shù)字MEMS傳聲器頻響一致性好,其內(nèi)部使用的Sigma-Delta ADC 采用1 bit變換技術(shù),克服了采用多比特模數(shù)轉(zhuǎn)換時帶來的線性誤差、糾錯困難等缺點,可以提供更好的信噪比、抗射頻和抗電磁干擾能力[12]。此外,在性能方面,MEMS傳聲器采樣分辨率高,對環(huán)境溫度變化具有良好的魯棒性,頻率響應(yīng)平坦,且同批次的MEMS傳聲器具有良好的相位一致性。數(shù)字MEMS傳聲器的關(guān)鍵參數(shù)如下[13]:
(1)等效輸入噪聲(EIN),其值要在最小聲壓級和最大聲壓級之間,以適應(yīng)測量范圍。
(2)聲學(xué)過載點(AOP)或滿量程(max SPL),需滿足動態(tài)范圍大于100 dB的要求。
(3)平坦頻率響應(yīng)段,需選取靈敏度變化范圍在3 dB內(nèi)的頻率段。
SPH0645LM4H-B MEMS傳聲器如圖2所示,本文選擇Knowles公司的SPH0645LM4H-B傳聲器,其具體性能參數(shù)見表1。
2.2 傳聲器幾何陣列設(shè)計
聲成像性能主要體現(xiàn)在空間分辨率、主瓣寬度、有效動態(tài)范圍三方面[14]。而合適的傳聲器陣列幾何布置可以防止空間混疊、抑制旁瓣和消除鬼影,提高傳聲器陣列成像性能。關(guān)于傳聲器陣列的設(shè)計,相關(guān)學(xué)者已經(jīng)做過大量的研究和對比。Johnson等人[15]指出,將傳聲器陣列設(shè)計成非冗余間距可以有效地消除空間混疊現(xiàn)象的產(chǎn)生,此后,基于非冗余排列的螺旋陣列被廣泛應(yīng)用于波束形成。Prime等人[16]分別對單臂螺旋陣列和多臂螺旋陣列的傳聲器陣列性能進行比較研究。結(jié)果表明,等區(qū)域均勻分布的多臂螺旋陣列可在一定的動態(tài)范圍內(nèi)提供最佳的陣列分辨率,而在陣列中心區(qū)域增添傳聲器擺放的密度,可以獲得更好的動態(tài)范圍,但會增加主瓣寬度。Amaral等人[17]在Underbrink陣列基礎(chǔ)上提出一種Annular陣列,在設(shè)計上更加靈活,極大提高了動態(tài)范圍。Annular陣列設(shè)計參數(shù)圖如圖3所示。
Annular陣列設(shè)計的關(guān)鍵參數(shù):傳聲器陣列孔徑D,最內(nèi)圈傳聲器環(huán)直徑d,傳聲器陣列的環(huán)數(shù)Q,以及每個環(huán)的傳聲器數(shù)量Cq,M表示傳聲器的總數(shù),即:
2.3 采集驅(qū)動模塊設(shè)計
根據(jù)2.1節(jié)中采用的數(shù)字MEMS傳聲器的采樣特點、接口時序和輸出模式選擇,PL端需要設(shè)計串并轉(zhuǎn)換模塊、延時模塊、I2S數(shù)據(jù)輸出接口協(xié)議模塊以及數(shù)據(jù)解析模塊。其中,I2S接口輸出時序如圖5所示。
圖5中,WS信號控制左、右聲道的切換。當(dāng)WS信號為低電平時,輸出數(shù)據(jù)為左聲道數(shù)據(jù);當(dāng)WS信號為高電平時,輸出數(shù)據(jù)為右通道數(shù)據(jù)。本采集驅(qū)動模塊設(shè)計將左、右通道采集數(shù)據(jù)緩存至對應(yīng)通道(CH)寄存器,通過一根DATA傳輸數(shù)據(jù),減少其數(shù)字I/O的使用個數(shù)。CLK信號為MIC工作時鐘,其時序規(guī)定在時鐘的上升沿開始采集數(shù)據(jù),在時鐘的下降沿,將采集數(shù)據(jù)發(fā)送至DATA信號線上,采樣時鐘周期T與FPGA內(nèi)置時鐘clkFPGA的換算關(guān)系為:
其中,clk為采樣率,Sbit為數(shù)據(jù)采集位數(shù) 。
PL端在接收數(shù)據(jù)時,每路數(shù)據(jù)采集都對應(yīng)一個FIFO寄存器,F(xiàn)IFO的存儲深度設(shè)為2 048, 數(shù)據(jù)位寬設(shè)為32 bit,為了防止FIFO存儲溢出,每次采集的數(shù)據(jù)大小為4 KB。其I2S接口數(shù)據(jù)采集模塊通道選擇部分關(guān)鍵偽代碼具體如下。
為了驗證采集驅(qū)動模塊的正確性,在設(shè)計仿真時,利用[WT5HZ]$readmemb系統(tǒng)函數(shù)構(gòu)造32 bit激勵數(shù)據(jù),并設(shè)置循環(huán)次數(shù)為8。仿真后得到的時序圖如圖6所示。
由圖6的仿真結(jié)果可知,在mic_ws信號為低電平時,采集CH1數(shù)據(jù),在時鐘上升沿將數(shù)據(jù)緩存至left_mic_data寄存器中;CH2通道與之相反。仿真結(jié)果表明,PL端的I2S時序符合I2S輸出格式傳輸協(xié)議,數(shù)據(jù)傳輸與設(shè)定值一致。采樣周期為25 μs,滿足其設(shè)定的40 k采樣率,并且能按設(shè)計的數(shù)據(jù)存儲大小依次按通道存儲數(shù)據(jù)。此外,為了滿足采集有效控制的需求,在PL端增添了簡易人機交互功能,即通過按鍵和LED來控制采集的進行和完成。其中,分別設(shè)置采集進行LED1和采集完成LED2,當(dāng)采集進行時,LED1會以500 ms的頻率閃爍,表示采集正在進行;當(dāng)采集完成時,LED2點亮,表明采集完成。
2.4 傳輸機制模塊設(shè)計
為實現(xiàn)36通道數(shù)據(jù)與PS端DDR存儲數(shù)據(jù)之間無失真的交互,本文傳輸機制模塊設(shè)計了一套Ping-Pang數(shù)據(jù)輪詢交互機制。通過Ping-Pang操作,可以完成多通道數(shù)據(jù)并行無縫緩存和處理。傳輸機制流程圖如圖7所示。
具體過程為:根據(jù)I2S數(shù)據(jù)傳輸協(xié)議特性,一組I2S數(shù)據(jù)傳輸分別控制2個通道MEMS傳聲器采集,設(shè)為一組(GROUP),共記18組。每一組都對應(yīng)設(shè)置一個獨立的FIFO緩存空間,各通道每次固定向FIFO緩存器發(fā)送4 KB數(shù)據(jù)量,通過有限狀態(tài)機,控制各通道數(shù)據(jù)依次有序傳輸至DDR_FIFO緩存器中,完成Ping-Pang切換。DDR_FIFO是帶有AXI接口的自定義模塊,可將PL端采集到的數(shù)據(jù)通過AXI4高速總線,傳輸至PS端。AXI驅(qū)動模塊設(shè)置了PL端DDR_FIFO中數(shù)據(jù)向PS端DDR中數(shù)據(jù)傳輸協(xié)議。首先設(shè)置數(shù)據(jù)傳輸?shù)氖椎刂窞?x10000000,將DDR_FIFO中每傳輸4 KB數(shù)據(jù)打包成一個BLOCK塊,每個BLOCK塊對應(yīng)一組MEMS傳聲器數(shù)據(jù),通過上述數(shù)據(jù)輪詢機制,可實現(xiàn)多通道數(shù)據(jù)依次傳輸至DDR中。
本文利用36個FIFO塊來緩存PL端數(shù)據(jù),緩解PS端的處理壓力。若采用單緩存塊,頻繁的讀寫操作會導(dǎo)致PS端輪詢標(biāo)志位檢測頻率激增,并且讀寫處理需要固定周期執(zhí)行,易使系統(tǒng)整體處理速度下降。而采用多組FIFO塊緩存的方式,利用存儲資源換取存儲速率,可提高系統(tǒng)整體的處理速度。
3 實驗驗證
為了驗證本文設(shè)計采集系統(tǒng)的正確性以及數(shù)字MEMS傳聲器的聲源成像效果,在室內(nèi)普通環(huán)境進行實際測試。實驗現(xiàn)場布置圖如圖8(a) 所示,以陣列中心為原點,圖中所標(biāo)注的(X,Y,Z) 建立坐標(biāo)系。2個揚聲器由三腳架支撐,垂直高度為1.15 m(傳聲器陣列高度與其一致),間距為1 m,距離傳聲器2 m。聲源1坐標(biāo) (-0.5,0,2) m、聲源2坐標(biāo)(0.5,0,2) m,2個揚聲器同時播放高斯白噪聲。
設(shè)置采樣率為40 kHz,采樣時長為0.1 s,采用傳統(tǒng)頻域波束形成算法,選取3.8~3.9 kHz頻率段進行成像。成像效果如圖8(b) 所示,利用本文所設(shè)計的FPGA成像系統(tǒng)可以準(zhǔn)確分辨出2個揚聲器位置,聲成像圖與聲源實際分布一致,驗證了所設(shè)計成像系統(tǒng)的正確性。同時,成像主瓣較窄,有效抑制旁瓣,驗證了所設(shè)計傳聲器陣列的優(yōu)異性能。
4 結(jié)束語
本文設(shè)計的基于FPGA傳聲器陣列的成像系統(tǒng),采用了數(shù)字MEMS傳聲器研制的36陣元Annular陣列。其中,數(shù)字MEMS傳聲器,體積小,功耗低,具有較好的相位一致性;36陣元Annular陣列,具有較好的成像分辨率,有效抑制旁瓣、衰減鬼影,防止混疊。此外,本文以FPGA為硬件系統(tǒng)核心,實現(xiàn)了多通道數(shù)據(jù)同步采集和數(shù)據(jù)完整編存。與傳統(tǒng)的采用駐極體式傳聲器陣列采集系統(tǒng)相比,實現(xiàn)了數(shù)據(jù)采集系統(tǒng)的小型化,具有低功耗以及低成本的優(yōu)勢??偨Y(jié)可得到如下結(jié)論:
(1)數(shù)字MEMS傳聲器,體積小,功耗低,具有較好的相位一致性,不需要設(shè)計額外的外圍電路,極大降低了硬件系統(tǒng)板上資源。在多通道同步采集中具有優(yōu)勢,有利于聲學(xué)可視化技術(shù)的應(yīng)用推廣。
(2)FPGA作為硬性系統(tǒng)核心,其靈活的邏輯設(shè)計、可控制多通道數(shù)據(jù)同時采集和豐富的I/O引腳便于通道數(shù)量拓展。其低功耗、體積小以及魯棒性好的特點也為后信號處理算法的嵌入式設(shè)計提供穩(wěn)定的前端。
(3)所研制的36陣元Annular陣列,具有較好的成像分辨率,有效抑制旁瓣、衰減鬼影,防止混疊??蛇m用于普通室內(nèi)環(huán)境,具有良好的抗噪性能,成像效果穩(wěn)定,在實際工程環(huán)境中具有良好的應(yīng)用前景。
參考文獻(xiàn)
[1]褚志剛,楊洋,倪計民,等. 波束形成聲源識別技術(shù)研究進展[J]. 聲學(xué)技術(shù),2013,32(5):430-435.
[2]SHIMON R, ASAY B, DASCHER D, et al.InP IC technology powers agilent's infiniium 90000 X-Series Real time oscilloscope[C]//Compound Semiconductor Integrated Circuit Symposium.Monterey, CA, USA :IEEE, 2010:1-4.
[3]ZDEMIR A T, ATCI A. A TMS320C6416 DSP-based high-speed data acquisition system[C]// 2015 23nd Signal Processing and Communications Applications Conference (SIU) .Malatya, Turkey :IEEE, 2015:1-4.
[4]BAOShuang, YAN Hairong, CHI Qingping, et al. FPGA-based reconfigurable data acquisition system for industrial sensors[J]. IEEE Transactions on Industrial Informatics, 2017,13(4):1503-1512.
[5]SILVA B D, BRAEKEN A, TOUHAFI A. FPGA-based architectures for acousticbeamforming with microphone arrays: Trends, challenges and research opportunities[J]. Computers, 2018, 7(3):41.
[6]ZHANG Xin, SONG Enliang, HUANG Jingchang, et al. Acoustic source localization via subspace based method using small aperture MEMS arrays[J]. Journal of Sensors, 2014(11):1-14.
[7]ZHANG Xin, HUANG Jingchang, SONG Enliang, et al. Design of small MEMS microphone array systems for direction finding of outdoors moving vehicles[J]. Sensors, 2014, 14(3):4384-4398.
[8]HAFIZOVIC I, NILSEN C C, Kj[WT6BZ]φ[WT6B1]LERBAKKEN M, et al. Design and implementation of a MEMS microphone array system for real-time speech acquisition[J]. Applied Acoustics, 2012, 73(2):132-143.
[9]JELMER T, DOMINGUEZ F, SILVA B D, et al.SoundCompass: A distributed MEMS microphone array-based sensor for sound source localization[J]. Sensors (Basel, Switzerland), 2014, 14(2):1918-1949.
[10]IZQUIERDO A, VILLACORTA J, LARA D, et al. Design and evaluation of a scalable and reconfigurable multi-platform system for acoustic imaging[J]. Sensors, 2016, 16(10):1671.
[11]陳林松,曹躍云,郭文勇,等. 基于數(shù)字式MEMS傳感器的聲陣列成像系統(tǒng)設(shè)計[J]. 應(yīng)用聲學(xué),2015,34(5):439-444.
[12]CHRISTEN T. A 15-bit 140-[WT6BZ]μ[WT6B1]W scalable-bandwidth inverter-based [WT6BZ]ΔΣ[WT6B1]modulator for a MEMS microphone with digital output[J]. IEEE Journal of Solid State Circuits, 2013, 48(7):1605-1614.
[13]PFLUG P. Aspects of the use of MEMS microphones in phased array systems[C]// INTER-NOISE and NOISE-CON Congress and Conference Proceedings,.HongKong, China:the Hong Kong Institute of Acoustics and the Hong Kong Polytechnic University, 2017:5093-5103.
[14]褚志剛,楊洋,蔣忠翰. 波束形成傳聲器陣列性能研究[J]. 傳感技術(shù)學(xué)報,2011,24(5):665-670.
[15]JOHNSON D H, DUDGEON D E. Array signal processing: Concepts and techniques[M].New York, NY,United States:Simon & Schuster, Inc., 1992.
[16]PRIME Z, DOOLAN C. A comparison of popularbeamforming arrays[C]//Proceedings of Acoustics 2013 Victor Harbor: Science Technology and Amenity, Annual Conference of the Australian Acoustical Society .Victor Harbor, Australia:Australian Acoustical Society, 2013:1-7.
[17]AMARAL D, CARLOS S, DE M. Design of microphone phased arrays for acoustic beamforming[J]. Journal of the Brazilian Society of Mechanical Sciences & Engineering, 2018, 40(7):354.
基金項目: 國家自然科學(xué)基金(51605274); 上海市地方院校能力建設(shè)項目(20030501000)。
作者簡介: 何 睿(1998-),男,碩士研究生,主要研究方向:聲學(xué)成像、FPGA系統(tǒng)開發(fā); 胡定玉(1987-) ,男,博士,副教授,主要研究方向:聲陣列信號處理、軌道車輛故障診斷技術(shù)。
通訊作者: 胡定玉 Email:dyhu1987@sues.com
收稿日期: 2021-04-16