基于代理模型的麥克風相位陣列設(shè)計技術(shù)研究

丁存?zhèn)? 李周復(fù), 張 雪, 周國成

(中國航空工業(yè)空氣動力研究院, 哈爾濱 150001)

0 引 言

麥克風相位陣列測量技術(shù)的應(yīng)用[1-5]，使得聲學風洞試驗?zāi)軌驅(qū)υ囼災(zāi)Ｐ偷闹饕肼曉催M行定位，獲取聲源強度、頻率等豐富的聲學特性信息，成為當前飛行器聲學試驗研究中不可或缺的測量手段。

麥克風相位陣列測量性能主要通過陣列分辨率(Resolution)和最大旁瓣水平(Maximum Sidelobe Level, MSL)進行衡量。陣列分辨率和最大旁瓣水平會隨著目標聲源頻率的改變而發(fā)生明顯變化。同一期聲學風洞試驗中模型構(gòu)型、來流風速等試驗狀態(tài)的不同，可能導致主要聲源位置、尺寸和噪聲頻率的顯著改變，對麥克風陣列分辨率和最大旁瓣水平的需求也隨之變化，從而超出了原有陣列的最佳適用范圍，導致某些狀態(tài)下測量結(jié)果的精準度受到影響。為了實現(xiàn)最佳測量效果，麥克風相位陣列也應(yīng)針對聲源特性的變化而實時改進，但是，這一改進過程必須滿足風洞試驗對現(xiàn)場操作的高效率要求。目前聲學試驗中普遍采用格柵式快速插拔陣列支架，能夠支持麥克風陣列的快速拆裝，因此，需要解決的問題是：根據(jù)試驗狀態(tài)對陣列分辨率和最大旁瓣水平的需求，如何在較短的時間內(nèi)完成麥克風相位陣列的優(yōu)化設(shè)計。

在麥克風相位陣列設(shè)計方面，國內(nèi)外研究人員研究了陣列尺寸和麥克風布置方式等因素對陣列分辨率和最大旁瓣水平的影響，并且將粒子群算法、模擬退火法等優(yōu)化設(shè)計方法引入了麥克風陣列設(shè)計中[4,6-13]。但是以上研究方法并未考慮陣列設(shè)計的耗時問題。盡管陣列設(shè)計中通常采用計算量較小的點擴散函數(shù)進行陣列性能分析，但是當聲源定位面內(nèi)掃描點數(shù)較多，優(yōu)化搜索空間較大時，計算量將會顯著增加，完成整個陣列優(yōu)化工作可能需要耗時數(shù)小時至數(shù)十個小時不等。這種計算耗時無法滿足風洞試驗過程中及時改進麥克風陣列的應(yīng)用需求，必須研究更加高效的麥克風相位陣列設(shè)計方法。

基于Kriging模型的優(yōu)化方法很好地解決了計算量與計算精度的問題，顯著提升了優(yōu)化效率，已經(jīng)在飛行器氣動設(shè)計領(lǐng)域得到了大量的應(yīng)用[14-21]。文獻[14]開展了旋翼槳尖氣動與噪聲綜合優(yōu)化設(shè)計研究，其中旋翼的氣動性能和噪聲特性需要應(yīng)用非定常流場求解方法進行計算，單點計算耗時較高。如優(yōu)化過程中全部采用非定常方法進行目標函數(shù)計算，則其計算量將急劇增加，無法滿足工程應(yīng)用需求。文章采用Kriging模型代替復(fù)雜費時的非定常計算，采用遺傳算法進行優(yōu)化，有效解決了計算耗時問題，建立了一種工程實用的旋翼槳尖優(yōu)化設(shè)計方法。此外，為了解決計算量過大的問題，文獻[15]應(yīng)用Kriging模型進行了機翼氣動/結(jié)構(gòu)剛度耦合優(yōu)化設(shè)計，文獻[16]將Kriging模型應(yīng)用于高升力翼型的優(yōu)化設(shè)計研究。可見，Kriging模型的應(yīng)用是提升優(yōu)化效率的有效方法。

為了最大限度提升麥克風相位陣列設(shè)計效率，本文借鑒飛行器氣動優(yōu)化設(shè)計領(lǐng)域的研究方法，創(chuàng)新性地將Kriging代理模型方法引入陣列優(yōu)化設(shè)計研究中。采用點擴散函數(shù)進行樣本點性能分析，建立Kriging模型，應(yīng)用遺傳算法進行優(yōu)化，形成了一種效率極高的麥克風陣列設(shè)計方法。該方法可有效改善陣列對目標聲源的測量性能，滿足聲學風洞試驗中針對特定聲源及時改進麥克風相位陣列的需求，有助于改善試驗測量精準度，提升試驗效率。

1 基于代理模型的陣列優(yōu)化設(shè)計方法

1.1 陣列性能分析方法

陣列性能分析是通過給定模擬點聲源，計算陣列對該聲源激勵的波束形成響應(yīng)來評估陣列性能的一種技術(shù)手段。陣列設(shè)計過程中，進行陣列性能量化分析時，波束形成響應(yīng)通常表示為點擴散函數(shù)(Point Spread Function, PSF)[22]。采用基于點擴散函數(shù)的計算方法分析陣列對模擬聲源的定位效果，是陣列設(shè)計過程中較為常用且可靠的性能分析方法，可以定量地評價麥克風陣列定位性能。計算方法如下：

陣列測量位于yj處的單位聲源，獲得互譜矩陣(Cross Spectrum Matrix, CSM)：

(1)

(2)

則位于yj的單位聲源在掃描點k處的波束形成輸出可表示為：

(3)

上式為PSF的定義式，Wk為對應(yīng)麥克風的權(quán)向量，是關(guān)于引導向量g的函數(shù)，取決于麥克風陣列與掃描平面的位置、測量環(huán)境與聲源頻率。

采用點擴散函數(shù)完成全部掃描點上的波束形成輸出計算后，陣列性能分析結(jié)果通常以聲壓級云圖的形式給出，從中可以獲取陣列分辨率和最大旁瓣水平兩個主要性能參數(shù)。其中，陣列分辨率描述了陣列定位聲源和分辨臨近聲源的能力；陣列的最大旁瓣水平表明陣列放大目標聲源、抑制背景噪聲的能力。

1.2 設(shè)計變量與約束條件

與其他陣列形式相比，多臂螺旋陣列能夠在較寬的頻率范圍內(nèi)保持良好的旁瓣抑制性能，因此在聲學試驗中得到了廣泛應(yīng)用，本文主要針對這種陣列進行優(yōu)化設(shè)計舉例。已有的多臂螺旋陣列設(shè)計方法中，螺旋線上麥克風的位置主要采用等弧長分布，或者等環(huán)形區(qū)域面積分布的方式確定[22]，無需進行優(yōu)化計算，但是設(shè)計空間較為單一。本文將多臂螺旋陣列的螺旋線軌跡、麥克風分布位置、旋臂數(shù)等信息進行參數(shù)化，拓展了陣列設(shè)計空間，通過對設(shè)計參數(shù)的優(yōu)化，搜索得到性能優(yōu)良的陣列方案。多臂螺旋陣列設(shè)計參數(shù)包括：懸臂數(shù)目n，螺旋線系數(shù)a和k，麥克風沿螺旋線分布參數(shù)q。

對數(shù)螺旋線可用極坐標表示如下：

r=a·ekθ

(4)

其中，r為極徑，θ為極角，a、k為常系數(shù)。

當系數(shù)a、k確定之后，螺旋線的軌跡即為確定值。麥克風在螺旋線上的分布通過參數(shù)q進行確定，具體方法如下：

相鄰兩麥克風間螺旋線弧長成等比增長，比值為q。對于最大極徑為R的對數(shù)螺旋線，當q=1時，螺旋線上麥克風為均勻分布；q≠1時，第m個麥克風的極坐標為：

(5)

(6)

其中，j為單條懸臂上麥克風數(shù)量，由整個陣列可用的麥克風總數(shù)和懸臂條數(shù)決定。

對于有n條旋臂的螺旋陣列，將陣列面板按2π/n弧度等分為n份，并通過麥克風坐標點的旋轉(zhuǎn)復(fù)制得到整個陣列坐標。

試驗環(huán)境、設(shè)備以及試驗對象等因素決定了麥克風陣列設(shè)計需要考慮的約束條件。試驗場地的空間限制和聲源定位范圍約束了陣列最大直徑；試驗中可供使用的麥克風總數(shù)量決定了陣列單元數(shù)；麥克風尺寸和安裝方式?jīng)Q定了兩麥克風最小間距。因此，優(yōu)化過程中的約束條件包括：陣列直徑、麥克風數(shù)量、麥克風最小間距。

1.3 基于Kriging模型的遺傳算法

遺傳算法不僅具有全局性優(yōu)化的特點，而且具有良好的魯棒性、可靠性和移植性，在工程優(yōu)化設(shè)計領(lǐng)域應(yīng)用較為廣泛。但是遺傳算法通常需要的種群規(guī)模龐大，容易導致計算量過大。本文采用基于Kriging模型的遺傳算法開展陣列優(yōu)化設(shè)計。Kriging模型在保證陣列性能預(yù)測精度的同時，極大地提升了優(yōu)化效率，同時使遺傳算法的種群規(guī)模不再受到計算量的限制，有效拓展了優(yōu)化搜索空間。

Kriging代理模型是一種插值模型，由一個參數(shù)模型和一個非參數(shù)隨機過程構(gòu)成，通過選擇適當?shù)膬?yōu)化算法計算得到模型參數(shù)最優(yōu)值，進而建立Kriging模型[18-19]。圖1給出了基于代理模型的陣列優(yōu)化方法流程圖。其基本思想為：

(1) 采用試驗設(shè)計方法生成設(shè)計空間內(nèi)的樣本點；

(2) 通過陣列性能分析程序計算樣本點的響應(yīng)值；

(3) 由樣本信息及其響應(yīng)值構(gòu)建代理模型；

(4) 采用遺傳算法進行優(yōu)化，種群適應(yīng)度計算由代理模型完成；

(5) 按照選定的加點方法，搜索性能更優(yōu)的樣本點；

(6) 采用陣列性能分析程序?qū)ψ顑?yōu)樣本點進行分析；

(7) 判斷是否達到設(shè)計目標，若是，則輸出結(jié)果，停止優(yōu)化；若否，則對原樣本集進行更新，以建立更加精確的代理模型，返回步驟(3)，重復(fù)優(yōu)化過程直至滿足停止條件。

圖1 優(yōu)化設(shè)計方法流程圖

2 麥克風陣列優(yōu)化設(shè)計算例

針對分辨率和最大旁瓣水平兩個指標分別進行陣列優(yōu)化，以驗證該優(yōu)化設(shè)計方法的實用性。當分辨率為目標函數(shù)時，設(shè)置滿足試驗測量要求的最大旁瓣水平下限指標作為約束條件；以最大旁瓣水平為目標函數(shù)時，則給定分辨率指標作為約束。

2.1 多臂螺旋陣列分辨率優(yōu)化

針對1m直徑量級63通道多臂螺旋陣列進行優(yōu)化設(shè)計，優(yōu)化問題總結(jié)如下：

設(shè)計狀態(tài)：2kHz單極子模擬聲源，定位區(qū)域1m×1m；

目標函數(shù)：分辨率最小，最大旁瓣水平≤-14dB；

設(shè)計變量：懸臂數(shù)目n，螺旋線系數(shù)a，k，分布參數(shù)q；

設(shè)計約束：陣列直徑D=1m，麥克風最小間距dmin≥0.02m，麥克風總數(shù)Nm=63。

共調(diào)用60次陣列性能分析程序進行樣本點計算，建立了代理模型。表1給出了Kriging代理模型的預(yù)測結(jié)果與數(shù)值模擬分析結(jié)果對比，代理模型預(yù)測誤差小于1%，預(yù)測精度完全滿足優(yōu)化設(shè)計需求。遺傳算法優(yōu)化過程中，以代理模型進行種群適應(yīng)度計算，種群規(guī)模為200，交叉概率0.8，變異概率0.05，共進行了100代遺傳優(yōu)化。優(yōu)化程序采用單核運行，整個優(yōu)化過程耗時約10min，獲得了性能較優(yōu)異的陣列方案。可見，Kriging代理模型的應(yīng)用顯著減少了計算時間，極大地拓展了優(yōu)化搜索的范圍。

表1 Kriging模型預(yù)測誤差Table 1 Prediction error of Kriging model

取每個設(shè)計變量上、下限的平均值為初始陣列，作為設(shè)計參考，優(yōu)化前后設(shè)計變量取值如表2所示。優(yōu)化前后陣列單元分布情況如圖2所示，可見優(yōu)化后陣列從中心向邊緣方向，單條螺旋線上麥克風間距逐漸減小，陣列單元的分布有向陣列邊緣移動的趨勢。圖3給出了2kHz模擬聲源下陣列響應(yīng)云圖，可見優(yōu)化陣列的分辨率有明顯改善。

表2 分辨率優(yōu)化陣列設(shè)計變量Table 2 Design variables of resolution optimized array

(a) 初始陣列

(b) 分辨率優(yōu)化陣列

Fig.2Schematicofinitialarrayandresolutionoptimizedmulti-armspiralmicrophonearray

與全尺寸飛機相比，航空聲學風洞試驗中的試驗?zāi)Ｐ屯ǔＪ墙?jīng)過大幅縮比的，主要噪聲的頻率被相應(yīng)放大。因此，本文主要進行1～6kHz頻率范圍內(nèi)麥克風陣列測量性能分析，能夠滿足航空聲學風洞試驗的大部分測量需求。表3給出了優(yōu)化后陣列與初始陣列性能對比數(shù)據(jù)，在1～6kHz頻率范圍內(nèi)，優(yōu)化陣列的分辨率均優(yōu)于初始陣列，以設(shè)計點2kHz頻率性能為例，優(yōu)化陣列分辨率從0.187m改善為0.164m，旁瓣聲壓級(SPL)略有升高，但最大旁瓣水平(MSL)滿足-14dB的約束要求。

(a) 初始陣列結(jié)果

(b) 分辨率優(yōu)化陣列結(jié)果

Fig.3Responsesofinitialarrayandresolutionoptimizedmicrophonearray

表3 初始陣列與分辨率優(yōu)化麥克風陣列性能Table 3 Performance of initial array and resolution optimized microphone array

2.2 多臂螺旋陣列最大旁瓣水平優(yōu)化

麥克風陣列最大旁瓣水平優(yōu)化算例的設(shè)計狀態(tài)、設(shè)計變量、設(shè)計約束等均保持不變，目標函數(shù)變化為：最大旁瓣水平降低、分辨率≤0.2m。

優(yōu)化后陣列單元分布如圖4所示，可見優(yōu)化后陣列單元向中心區(qū)域集中。圖5給出了2kHz模擬聲源下的陣列響應(yīng)云圖，可見優(yōu)化陣列的最大旁瓣水平有明顯改善。表4給出了優(yōu)化后陣列與初始陣列性能對比數(shù)據(jù)，可以看出在1～6kHz頻率范圍內(nèi)，優(yōu)化陣列的最大旁瓣水平明顯優(yōu)于初始陣列，以設(shè)計點2kHz頻率下性能為例，優(yōu)化陣列最大旁瓣水平有5dB的改善，分辨率滿足0.2m的優(yōu)化約束條件。

圖4 最大旁瓣水平優(yōu)化陣列示意圖

圖5 最大旁瓣水平優(yōu)化陣列響應(yīng)圖

初始陣列優(yōu)化陣列頻率/kHz最大旁瓣水平/dB分辨率/m最大旁瓣水平/dB分辨率/m1-17.4120.388-20.5370.4122-16.8820.187-21.9220.1983-15.5350.123-19.3460.1294-15.2520.089-16.1050.0945-14.7710.070-15.6980.0746-14.1880.057-14.0090.061

以上陣列優(yōu)化算例表明：本文發(fā)展的麥克風相位陣列設(shè)計方法能夠根據(jù)試驗具體需求，優(yōu)化設(shè)計出性能優(yōu)良的陣列方案，且優(yōu)化效率極高，滿足聲學試驗過程中隨時改進麥克風陣列的特殊需求。對最大旁瓣水平進行優(yōu)化時，陣列單元沿螺旋線分布有向陣列中心聚集的趨勢，同時分辨率會略變差。相反，對陣列分辨率進行優(yōu)化時，陣列單元沿螺旋線分布有向陣列邊緣移動的趨勢，最大旁瓣水平則會略變差。可見，最大旁瓣水平和分辨率這兩個性能指標在一定條件下是此消彼長的關(guān)系，無法同時優(yōu)化到最優(yōu)值。陣列設(shè)計工作需要考慮聲學試驗的聲源強度、特征尺寸等具體因素，在一個指標滿足基本約束的情況下，盡量提升最重要的核心指標。作為約束的指標可能會有所下降，但仍滿足基本需求，而核心指標的大幅提升則能夠明顯改善試驗測量結(jié)果的精準度。

2.3 矩形陣列設(shè)計

該陣列優(yōu)化設(shè)計方法具有較強的可拓展性，只需針對參數(shù)化部分進行修改，即可應(yīng)用于其他形式的陣列設(shè)計。針對81通道的1m邊長矩形陣列，將陣列單元的橫、縱坐標作為設(shè)計參數(shù)。以陣列單元在平面內(nèi)的均勻分布作為初始方案，通過對陣元分布趨勢的優(yōu)化，提升陣列最大旁瓣水平。圖6給出了優(yōu)化前后矩形陣列單元分布示意圖。

(a) 初始矩形陣列

(b) 優(yōu)化矩形陣列

Fig.6ComparisonofinitialarrayandMSLoptimizedrectangularmicrophonearray

圖7給出了2kHz模擬聲源下陣列響應(yīng)云圖，可見旁瓣水平有明顯改善。表5給出了1～6kHz頻率下陣列性能對比數(shù)據(jù)，優(yōu)化陣列的最大旁瓣水平改善顯著，分辨率略變差。其中，聲源頻率為2kHz時，優(yōu)化后矩形陣列最大旁瓣水平改善了4.2dB，提升約35.85%；分辨率下降約5.62%。整體而言，陣列性能有明顯提升。

(a) 初始矩形陣列結(jié)果

(b) 優(yōu)化矩形陣列結(jié)果

Fig.7ResponsesofinitialarrayandMSLoptimizedrectangularmicrophonearrayresponses

表5 初始陣列與優(yōu)化矩形麥克風陣列性能Table 5 Performance of initial array and optimized rectangular microphone array

3 結(jié) 論

將基于Kriging代理模型的遺傳算法引入麥克風相位陣列設(shè)計研究中，形成了高效率麥克風相位陣列設(shè)計方法，得到研究結(jié)論如下：

(1) 陣列優(yōu)化結(jié)果表明：該優(yōu)化設(shè)計方法可以根據(jù)試驗研究的具體需求進行麥克風相位陣列設(shè)計，有效提升麥克風相位陣列的分辨率或最大旁瓣水平性能，是一種可靠實用的麥克風相位陣列設(shè)計方法。

(2) Kriging代理模型的引入使得陣列優(yōu)化設(shè)計過程耗時極短，可以滿足聲學風洞試驗中隨時針對聲源特點而改進麥克風相位陣列的特殊應(yīng)用需求，有利于改善聲學風洞試驗測量數(shù)據(jù)的精準度，顯著提升試驗效率。

(3) 麥克風相位陣列的兩個主要性能參數(shù)(分辨率和最大旁瓣水平)通常無法同時達到極值，陣列設(shè)計時需要根據(jù)試驗測量的實際需求，搜尋陣列性能的最佳平衡點。

(4) 該陣列優(yōu)化設(shè)計方法具有較強的可拓展性，針對參數(shù)化部分進行修改后，即可應(yīng)用于其他類型陣列的設(shè)計工作。