朱榮釗，劉 晶，劉文超，周艷玲，曾張帆

(1.湖北大學(xué) 計(jì)算機(jī)與信息工程學(xué)院，武漢 430062; 2.武漢晴川學(xué)院 計(jì)算機(jī)學(xué)院，武漢 430204)

0 引言

隨著科技的發(fā)展，我們身邊的噪聲源越來(lái)越多。廣場(chǎng)舞音響、汽車(chē)?guó)Q笛、店鋪促銷(xiāo)喇叭等產(chǎn)生的噪聲在不同程度上影響了人們的正常生活[1]。聲頻定向技術(shù)是一種可以實(shí)現(xiàn)可聽(tīng)聲在一定區(qū)域內(nèi)傳播的新技術(shù)，從而可以構(gòu)造出獨(dú)立的語(yǔ)音空間。基于此現(xiàn)狀，本文設(shè)計(jì)出一種聲頻定向揚(yáng)聲器系統(tǒng)，相比于已有的聲頻定向系統(tǒng)，本系統(tǒng)的創(chuàng)新點(diǎn)在于基于FPGA的信號(hào)處理模塊的搭建以及功率放大電路的設(shè)計(jì)，經(jīng)過(guò)一定環(huán)境下的測(cè)試，實(shí)驗(yàn)結(jié)果表明本系統(tǒng)具有較強(qiáng)的指向性，可以在一定程度上緩解噪聲污染。

1 系統(tǒng)結(jié)構(gòu)與原理

聲頻定向系統(tǒng)可以產(chǎn)生高指向性可聽(tīng)聲主要是利用了超聲波在空氣中的非線(xiàn)性傳播效應(yīng)，Helmholtz提出的非線(xiàn)性聲學(xué)理論可以對(duì)這種非線(xiàn)性過(guò)程進(jìn)行初步解釋[2]。該理論指出，若向超聲波換能器陣列輸入兩列頻率分別為f1、f2的超聲波電信號(hào)時(shí)，在空氣的非線(xiàn)性自解調(diào)作用下，將會(huì)產(chǎn)生頻率分別為f1、f2、f1+f2、f1-f2以及各階諧波在內(nèi)的復(fù)雜聲波，由于聲波在空氣中的衰減系數(shù)與聲波自身頻率的平方成正比，所以頻率較高的f1+f2、f1、f2以及各階次諧波將會(huì)很快被空氣吸收掉，只有處于可聽(tīng)聲范圍內(nèi)的差頻信號(hào)繼續(xù)傳播。若超聲波頻率選擇的合適，則可產(chǎn)生與原輸入可聽(tīng)聲信號(hào)頻率一致的高指向性可聽(tīng)聲，這便是實(shí)現(xiàn)聲頻定向揚(yáng)聲器系統(tǒng)的基本原理，如圖1所示。

圖1 非線(xiàn)性傳播效應(yīng)

聲頻定向揚(yáng)聲器系統(tǒng)基于參量聲學(xué)陣?yán)碚摱O(shè)計(jì)。目前，主要存在兩種理論，分別為“Westervelt方程”和“Berktay遠(yuǎn)場(chǎng)解”[4]?！癢estervelt方程”可解釋如何通過(guò)兩列超聲波產(chǎn)生高指向性差頻波，而“Berktay遠(yuǎn)場(chǎng)解”則可解釋如何通過(guò)參量聲學(xué)陣原理產(chǎn)生寬帶高指向性聲頻波。

本文中涉及到的參量陣系統(tǒng)的輸入信號(hào)是寬帶信號(hào)，所以下面簡(jiǎn)要介紹一下“Berktay遠(yuǎn)程解”。Berktay在“Westervelt方程”的基礎(chǔ)上，推導(dǎo)了寬帶信號(hào)作為參量聲學(xué)陣輸入時(shí)通過(guò)非線(xiàn)性產(chǎn)生的聲壓為：

(1)

其中:β為非線(xiàn)性系數(shù)，S為聲源面積，p0為原波聲壓幅值，E(τ)為調(diào)制包絡(luò)函數(shù)，z為聲波傳播距離，τ=t-z/c0為延遲時(shí)間。

式(1)明確表明寬帶參量聲學(xué)陣最后在空氣中自解調(diào)得到的聲壓p(t)與包絡(luò)E(τ)平方的兩次時(shí)間導(dǎo)成正比，即為“Berkaty遠(yuǎn)程解”[5]。本系統(tǒng)就是基于“Berktay遠(yuǎn)場(chǎng)解”理論，將輸入聲頻信號(hào)合理融入到包絡(luò)平方函數(shù)的E2(τ)中去，利用寬帶參量陣實(shí)現(xiàn)該聲頻信號(hào)的定向傳播?，F(xiàn)有聲頻定向系統(tǒng)信號(hào)處理方法皆是建立在“Berktay遠(yuǎn)場(chǎng)解”基礎(chǔ)上的。

2 信號(hào)處理系統(tǒng)的設(shè)計(jì)

該聲頻定向系統(tǒng)在結(jié)構(gòu)上包括前端音頻放大模塊、FPGA中信號(hào)處理模塊、ADDA模塊、功率放大模塊等，整體結(jié)構(gòu)如圖2所示。

圖2 系統(tǒng)組成框圖

需要對(duì)該系統(tǒng)組成進(jìn)行說(shuō)明的有：

1)超聲波載波信號(hào)的產(chǎn)生和信號(hào)的調(diào)制均在FPGA中進(jìn)行。

2)為了得到最好的指向性效果，在探究了超聲波換能器陣列數(shù)目及形狀的基礎(chǔ)上，最終確定該系統(tǒng)的超聲波換能器陣列由中心頻率為40kHz[6],直徑為16mm，數(shù)目為91的超聲換能器組成，形狀為正六邊形。

3)由于驅(qū)動(dòng)超聲波換能器需要較大的電壓，所以本系統(tǒng)選用OPA548芯片，自主設(shè)計(jì)功率放大電路，經(jīng)過(guò)測(cè)試分析，結(jié)果表明該功率放大電路可以較好地驅(qū)動(dòng)超聲波換能器模塊。

本系統(tǒng)軟件方面主要是基于DE1-SOC提供的硬件開(kāi)發(fā)環(huán)境和Quartus13.1提供的軟件開(kāi)發(fā)環(huán)境，采用SoC核心思想，通過(guò)ARM和FPGA之間的交互訪(fǎng)問(wèn)，實(shí)現(xiàn)對(duì)聲頻定向系統(tǒng)的靈活控制，包括LED燈、SDCard和聲束的偏轉(zhuǎn)角度等外圍設(shè)備。同時(shí)，實(shí)現(xiàn)了多種方式的輸入，包括麥克風(fēng)輸入、模擬音頻信號(hào)輸入以及SDCard數(shù)字音頻信號(hào)輸入等。

下面分別介紹關(guān)鍵模塊的設(shè)計(jì)方法。

2.1 前端音頻放大模塊的設(shè)計(jì)

若采用手機(jī)、MP3等輸入方式，則輸出的模擬音頻信號(hào)的電壓幅值較小，在A(yíng)D轉(zhuǎn)換模塊中，無(wú)法較好的進(jìn)行模/數(shù)轉(zhuǎn)換，所以需要將前端音頻信號(hào)進(jìn)行放大。OP07芯片是一種低噪聲，非斬波穩(wěn)零的雙極性運(yùn)算放大器集成電路，具有輸入偏置電流低和開(kāi)環(huán)增益高的特點(diǎn)，很適合作為本系統(tǒng)的前端小信號(hào)放大器。采用OP07芯片實(shí)現(xiàn)前端小信號(hào)的放大，放大倍數(shù)可以由滑動(dòng)變阻器進(jìn)行調(diào)節(jié)，具體電路圖如圖3，其中R1=R3=1 kΩ，滑動(dòng)變阻器R2選用3 296 W 103系列，阻值變化范圍為0～10 kΩ。

圖3 前端電壓放大電路

若采用SD卡輸入，其實(shí)現(xiàn)方法為：利用Quartus13.1里面的Qsys軟件快速生成讀取SD卡里面.WAV文件音頻數(shù)字信號(hào)所需的硬件配置，待系統(tǒng)生成成功之后，再到Quartus13.1里面對(duì)不同的硬件進(jìn)行相應(yīng)的引腳分配，同時(shí)可以在Eclipse軟件里面編寫(xiě)C語(yǔ)言程序?qū)崿F(xiàn)所需要的邏輯，最終成功捕捉到SD卡里面純正的數(shù)字音頻信號(hào)。

2.2 FPGA中信號(hào)處理模塊的設(shè)計(jì)

2.2.1 麥克勞林預(yù)處理

作為一種新概念聲源,聲頻定向系統(tǒng)與傳統(tǒng)揚(yáng)聲器的發(fā)聲機(jī)理有本質(zhì)區(qū)別。傳統(tǒng)揚(yáng)聲器可近似看作線(xiàn)性系統(tǒng),一般不需要進(jìn)行信號(hào)預(yù)處理就可以直接產(chǎn)生無(wú)指向性可聽(tīng)聲,而聲頻定向系統(tǒng)則是一種典型的非線(xiàn)性系統(tǒng),必須對(duì)信號(hào)進(jìn)行處理才能產(chǎn)生失真度小、具有高指向性的可聽(tīng)聲。

其具體實(shí)現(xiàn)方法為：根據(jù)非線(xiàn)性聲學(xué)理論，將音頻信號(hào)送入FPGA中進(jìn)行麥克勞林N階近似平方根處理[7]，本系統(tǒng)直接采用了Quartus13.1軟件里面的宏模塊，包括加法器、乘法器以及除法器等，通過(guò)DE1-SoC核心板進(jìn)行高速數(shù)學(xué)運(yùn)算，最終完成了對(duì)可聽(tīng)聲信號(hào)的預(yù)處理。

2.2.2 頻率為40 kHz的載波信號(hào)的生成

在信號(hào)調(diào)制部分的一個(gè)關(guān)鍵問(wèn)題是如何產(chǎn)生調(diào)制所需要的載波[7]。在本系統(tǒng)中，采用Altera公司的DE1-SoC開(kāi)發(fā)板作為信號(hào)處理核心器件，憑借FPGA的自身優(yōu)勢(shì)，自主制作40 kHz的載波信號(hào)的過(guò)程為：先用VS2013軟件編寫(xiě)C語(yǔ)言程序，程序的功能為輸出一個(gè)位寬為8位，深度為500的正弦表，然后將程序格式制作成MIF格式文件，再調(diào)用FPGA中宏模塊ROM，將制作的MIF格式文件固化到單口ROM里面，然后用Verilog語(yǔ)言寫(xiě)一個(gè)計(jì)算器對(duì)ROM進(jìn)行地址讀取，從而實(shí)現(xiàn)對(duì)ROM里面的正弦表數(shù)據(jù)進(jìn)行讀取，完成正弦信號(hào)的制作。

2.2.2 信號(hào)調(diào)制

利用FPGA對(duì)信號(hào)進(jìn)行超聲調(diào)制，使調(diào)制后的超聲波信號(hào)在空氣的非線(xiàn)性自解調(diào)作用下還原出可聽(tīng)聲。目前常用的調(diào)制算法有DSB法、平方根法、SSB法等，但是DSB法不滿(mǎn)足無(wú)失真包絡(luò)平方條件，平方根法的實(shí)現(xiàn)需要一個(gè)具有無(wú)限帶寬的物理系統(tǒng)，而實(shí)際物理系統(tǒng)無(wú)法提供這種無(wú)限帶寬，在輸入寬帶信號(hào)的條件下，SSB法存在較為嚴(yán)重的互調(diào)失真?；诖爽F(xiàn)狀，本作品采用N階近似平方根法[8]進(jìn)行信號(hào)調(diào)制，解決了平方根法的無(wú)限帶寬及信號(hào)失真問(wèn)題。采用Signal Tap采集到的調(diào)制波形如圖5，第一行為輸入1 kHz的正弦波，第二行為從ROM中調(diào)用的40 kHz載波，第三行為調(diào)制后信號(hào)的波形。

圖5 SignalTap采集的調(diào)制波形

2.3 信號(hào)的轉(zhuǎn)換與放大

經(jīng)過(guò)FPGA核心板高速運(yùn)算的調(diào)制信號(hào)從FPGA的GPIO口以8位的數(shù)據(jù)位寬并行輸出與外圍設(shè)備高速ADA模塊進(jìn)行對(duì)接，這里使用了FPGA提供的PLL IP Core，將產(chǎn)生的高速率時(shí)鐘送到高速ADA的DA 時(shí)鐘線(xiàn)上，通過(guò)高速率的精確的DA轉(zhuǎn)換，實(shí)現(xiàn)了調(diào)制信號(hào)從數(shù)字信號(hào)到模擬信號(hào)的跨變。

實(shí)現(xiàn)調(diào)制信號(hào)從數(shù)字信號(hào)轉(zhuǎn)變?yōu)槟M信號(hào)后，還存在從FPGA核心板輸出的調(diào)制信號(hào)強(qiáng)度很小的問(wèn)題。即經(jīng)由數(shù)模轉(zhuǎn)換后，信號(hào)強(qiáng)度還是很難驅(qū)動(dòng)超聲波換能器陣列。OPA548作為高電壓、高電流、寬輸出電壓擺幅功率運(yùn)算放大器，可以用來(lái)實(shí)現(xiàn)調(diào)制后輸出信號(hào)的功率放大。因此本系統(tǒng)獨(dú)立設(shè)計(jì)了對(duì)超聲波調(diào)制信號(hào)進(jìn)行放大的功率放大電路，實(shí)現(xiàn)了將7 V左右的超聲波調(diào)制信號(hào)放大到35 V左右，可以較好驅(qū)動(dòng)超聲波換能器。功率放大電路的設(shè)計(jì)如圖6所示。

圖6 功率放大電路的設(shè)計(jì)

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

聲頻定向揚(yáng)聲器系統(tǒng)整體實(shí)物圖如7和8所示。

圖7 聲頻定向揚(yáng)聲器系統(tǒng) 圖8 超聲波換能器 實(shí)物側(cè)視圖

圖7為聲頻定向揚(yáng)聲器系統(tǒng)側(cè)視圖，從上到下分為三層，第一層為Altera的DE1-SoC開(kāi)發(fā)板，在FPGA中實(shí)現(xiàn)信號(hào)的預(yù)處理和調(diào)制，第二層包括自主設(shè)計(jì)研發(fā)并完成的前端音頻放大模塊，AD/DA模塊、功率放大模塊、散熱模塊以及DC-DC降壓模塊，第三層為開(kāi)關(guān)電源NES-200-24，通過(guò)降壓等模塊，給整個(gè)系統(tǒng)提供合適的電壓。

圖8為自主設(shè)計(jì)的超聲波換能器，由91個(gè)直徑為15 mm的超聲波陣列頭組成，整體呈正六邊形。

除此之外，為了驗(yàn)證此聲頻定向系統(tǒng)設(shè)計(jì)方案的正確性和可行性，并檢驗(yàn)系統(tǒng)的定向傳聲性能，本章對(duì)聲頻定向進(jìn)行系統(tǒng)測(cè)試。通過(guò)研究相關(guān)的實(shí)驗(yàn)數(shù)據(jù)，從而得出聲場(chǎng)指向性、反射性、虛擬聲源效果等系統(tǒng)的性能參數(shù)，同時(shí)也為優(yōu)化系統(tǒng)設(shè)計(jì)提供依據(jù)。

在測(cè)試聲頻定向系統(tǒng)性能的過(guò)程中，必須選定以下測(cè)試條件：

1)大面積空曠場(chǎng)地；

2)弱噪聲干擾環(huán)境。

根據(jù)上述條件，按照以下步驟進(jìn)行具體測(cè)試：

1)選擇大面積空曠場(chǎng)地?cái)[放聲頻定向系統(tǒng)；

2)給系統(tǒng)供電，啟動(dòng)系統(tǒng)的各個(gè)模塊；

3)分別進(jìn)行定向傳聲測(cè)試、反射性測(cè)試、遠(yuǎn)距離傳播測(cè)試等；

4)記錄聲音強(qiáng)度，距離等參數(shù)；

5)將記錄的數(shù)據(jù)與理論結(jié)果對(duì)比，分析系統(tǒng)的優(yōu)勢(shì)及不足。

在測(cè)試系統(tǒng)的硬件性能過(guò)程中，測(cè)試主要針對(duì)系統(tǒng)邏輯電路在FPGA端的邏輯利用率、相關(guān)運(yùn)算速度等。測(cè)試在Quartus II 13.1中進(jìn)行。

3.1 定向性測(cè)試

根據(jù)以上測(cè)試條件，本系統(tǒng)在無(wú)噪聲或者小噪聲的環(huán)境下進(jìn)行了測(cè)試，采用單一變量的測(cè)試原則，分別用普通揚(yáng)聲器和聲頻定向系統(tǒng)播放相同的音頻，以播放器為中心，在半徑為5米的圓周上，分貝儀進(jìn)行180度旋轉(zhuǎn)采集，記錄分貝儀上的波形變化，測(cè)試結(jié)果如圖9所示。圖9(a)為設(shè)計(jì)的定向系統(tǒng)測(cè)試結(jié)果，圖9(b)為普通非定向揚(yáng)聲器的測(cè)試結(jié)果。

圖9 測(cè)試結(jié)果圖

從圖9(a)可以看出，在13～17 s的時(shí)間內(nèi)，分貝儀位于聲頻定向系統(tǒng)的正前方區(qū)域，聲音可以被聽(tīng)到，分貝儀測(cè)試曲線(xiàn)出現(xiàn)了升高現(xiàn)象，驗(yàn)證了聲頻定向系統(tǒng)具有定向性的特點(diǎn)。

由圖9(b)可以看出，對(duì)于普通揚(yáng)聲器，在以5 m為半徑的圓周上旋轉(zhuǎn)采集信號(hào)強(qiáng)度時(shí)，分貝儀曲線(xiàn)基本不變，大約為50 dB，驗(yàn)證了普通揚(yáng)聲器的聲音是360°全向傳播的。

3.2 遠(yuǎn)距離傳聲測(cè)試

根據(jù)理論知識(shí)可知，聲頻定向系統(tǒng)發(fā)出的可聽(tīng)聲具有能量集中，可遠(yuǎn)距離傳播的特點(diǎn)。下面將分別在一定距離處測(cè)試普通揚(yáng)聲器聲音強(qiáng)度與聲頻定向系統(tǒng)聲音強(qiáng)度，數(shù)據(jù)匯總表如表1所示。

表1 遠(yuǎn)距離傳播測(cè)試數(shù)據(jù)匯總表

測(cè)試標(biāo)號(hào)表示4次測(cè)試樣本，測(cè)試均以相同的測(cè)試方法、測(cè)試設(shè)備、測(cè)試環(huán)境下完成。分析數(shù)據(jù)可以得到以下結(jié)論：

1)在相同的傳播距離下，聲頻定向系統(tǒng)的聲音強(qiáng)度大于普通揚(yáng)聲器的聲音強(qiáng)度，符合聲頻定向系統(tǒng)具有能量集中的特點(diǎn)。

2)隨著傳播距離的增加，聲頻定向系統(tǒng)聲音強(qiáng)度的衰減速度小于揚(yáng)聲器聲音強(qiáng)度的衰減速度，符合聲頻定向系統(tǒng)的聲音可以遠(yuǎn)距離傳播的特點(diǎn)。

4 結(jié)語(yǔ)

本設(shè)計(jì)基于DE1-SoC硬件集成開(kāi)發(fā)環(huán)境，借助SoC嵌入式軟件編程技術(shù)，實(shí)現(xiàn)了聲頻定向傳播系統(tǒng)，打破了傳統(tǒng)揚(yáng)聲器全方位音頻傳播的特征，構(gòu)建了獨(dú)立的語(yǔ)音空間，可以緩解噪聲污染，提升通話(huà)保密性，在民用和軍用場(chǎng)合都具有極大的應(yīng)用潛力。除此之外，通過(guò)使用HPS，在大大提升系統(tǒng)靈活性的同時(shí)，降低了系統(tǒng)的設(shè)計(jì)難度，設(shè)計(jì)平臺(tái)的搭建也更為迅速。

從實(shí)際出發(fā)，該系統(tǒng)還具有以下設(shè)計(jì)特點(diǎn)：

1)將頻率為20～20 K的可聽(tīng)聲與頻率為40 K的高頻信號(hào)進(jìn)行調(diào)制，改變了可聽(tīng)聲全方位傳播的特性，使可聽(tīng)聲只能在指定區(qū)域內(nèi)傳播。最終創(chuàng)造出一種新概念聲源[11]，使可聽(tīng)聲具有定向性，反射性，能量集中等特點(diǎn)。

2)支持麥克風(fēng)、手機(jī)、Mp3等模擬音頻信號(hào)的輸入的同時(shí)，還支持SD卡等數(shù)字音頻信號(hào)的輸入，擴(kuò)大系統(tǒng)的使用范圍。

3)通過(guò)對(duì)換能器陣列形式、陣列數(shù)目、陣列間距的研究，采用正六邊形陣列，使系統(tǒng)具有更好的寬帶特性及平坦的頻率響應(yīng)特性。

4)在實(shí)現(xiàn)了機(jī)械式聲波偏轉(zhuǎn)的基礎(chǔ)上，我們后續(xù)的研究方向主要為相位控制聲束偏轉(zhuǎn)技術(shù)。相比于機(jī)械式聲束偏轉(zhuǎn)成本低、易實(shí)現(xiàn)的特點(diǎn)，相控陣聲束偏轉(zhuǎn)具有操作靈活、變化速率快、全方位偏轉(zhuǎn)等優(yōu)點(diǎn)。

但是，初步設(shè)計(jì)的聲頻定向揚(yáng)聲器系統(tǒng)的音質(zhì)無(wú)法與傳統(tǒng)揚(yáng)聲器媲美，若向大規(guī)模推向市場(chǎng)，還需進(jìn)一步研究如何提升音質(zhì)以及如何減少功耗的問(wèn)題。

[1] 倪曉鳴. 廣場(chǎng)舞噪音擾民多視角分析[J]. 體育文化導(dǎo)刊，2015(2)：13-16.

[2] Jamea J, Norris J O.HSS white paper[Z].USA: American Technology Corporation,2005.

[3] 陳 敏.聲頻定向系統(tǒng)理論與關(guān)鍵技術(shù)研究[D].成都：電子科技大學(xué)，2008.

[4] Berktay H O. Possible exploitation of non-linear acoustics in underwater transmitting applications[J]. Journal of Sound & Vibration, 1965, 2(4):435-461.

[5] Westervelt P J. Parametric Acoustic Array[J]. Journal of the Acoustical Society of America, 1963, 35(4):535-537.

[6] James J, Norris J O. HSS white paper[Z]. American Technology Corporation, San Diego, California, USA, 2005.

[7] Shi C, Kajikawa Y. Automatic gain control for parametric array loudspeakers[A]. International Conference on Acoustics, Speech, and Signal Processing[C]. 2016.

[8] Yoneyama M, Fujimoto J, Yu K, et al. The audio spotlight: an application of nonlinear interaction of sound waves to a new type of loudspeaker design[J]. Journal of the Acoustical Society of America, 1983, 73(57):1972-1986.

[9] 姬培鋒. 指向性聲源的理論與實(shí)驗(yàn)研究[D]. 濟(jì)南：山東科技大學(xué), 2005.

[10] 賈寶賢, 邊文鳳, 趙萬(wàn)生,等. 壓電超聲換能器的應(yīng)用與發(fā)展[J]. 壓電與聲光, 2005, 27(2):131-135.

[11] 董麗杰, 熊傳溪, 陳 娟,等. 聚偏氟乙烯-鈦酸鋇復(fù)合材料的介電性能[J]. 復(fù)合材料學(xué)報(bào), 2003, 20(3):122-126.

[12] H.O.Berktay, 錢(qián)祖文. 非線(xiàn)性聲學(xué)[J]. 聲學(xué)技術(shù), 1975(1).

[13] 王永兵, 侯 靜, 李志國(guó). 超聲相控陣檢測(cè)技術(shù)研究進(jìn)展[J]. 中國(guó)化工裝備, 2012, 14(6):15-18.

[14] 韓 捷, 張 坡, 尤國(guó)雷. 國(guó)家標(biāo)準(zhǔn)《揚(yáng)聲器主要性能測(cè)試方法》解讀[J]. 電聲技術(shù), 2012, 36(11):86-90.

[15] Hua H, Wang Y, Yan D. A low-cost dynamic range-finding device based on amplitude-modulated continuous ultrasonic wave[J]. Instrumentation & Measurement， IEEE, 2002, 51(2):362-367.

[16] Smith P R, Cowell D M J, Raiton B, et al. A PLL-based phased array method to minimize phase quantization errors and reducephasing-lobes[J]. 2010:1837-1840.