楊 垚 李 鋆 何元安

(1 哈爾濱工程大學水聲工程學院 哈爾濱 150001)

(2 中國船舶工業(yè)系統(tǒng)工程研究院 北京 100036)

0 引言

水聲材料作為水下聲波產生、傳播的重要介質，在我國海洋戰(zhàn)略不斷推進的今天，其重要性日益上升，因此，對水聲材料聲學參數(shù)的測量是水聲材料研究的基礎和關鍵。近年來，隨著對水聲材料的測試頻率逐漸降低，模擬海洋環(huán)境條件下的各種水聲材料及構件樣品的聲學性能測量需求越來越迫切。在試驗室條件下，對聲學材料小樣的聲學參數(shù)測量常在聲管中進行，常見的聲管包括脈沖管、行波管、駐波管[1]。

現(xiàn)有的行波場建立方法分為美國學者提出的傳遞矩陣法[2]與俄羅斯學者提出的迭代法[3]。中船重工第七一五研究所在“十五”期間建立了我國第一套行波管聲管，采用迭代法進行主動消聲，并通過對聲速的校準提高了測量效果[4]。2010年，七一五研究所采用虛擬行波法進一步提高了測量的效果，據(jù)此實現(xiàn)行波管的單頻及寬帶測量[5]。2012年，中國船舶工業(yè)系統(tǒng)工程研究院也建成了行波管，采用迭代法進行行波場建立，在原來國內行波管的基礎上提高了管內水聽器的相幅一致性、信號采集的同步能力及信噪比，采用稀土型縱向換能器作為耐高靜壓的低頻發(fā)射換能器，增強了發(fā)射換能器的低頻聲場穩(wěn)定性。

現(xiàn)有行波場建立方法尚存在一定問題，傳遞矩陣法使用時需要對傳遞矩陣中各參數(shù)進行測量計算，測量精度將影響行波場建立效果；迭代法使用中需經過較長迭代時間后才可實現(xiàn)行波場建立。為了避免現(xiàn)有方法存在的問題，改善行波場建立效果，本文基于中國船舶工業(yè)系統(tǒng)工程研究院的行波管，設計了一種基于黃金分割搜索算法的行波環(huán)境建立方法，從而滿足測量材料性能所需的環(huán)境，達到提升聲管測試性能的目的。

1 算法說明

行波場建立通過有源吸聲的方法，使用次發(fā)換能器發(fā)射抵消信號，令次發(fā)波形與主發(fā)波形經次發(fā)換能器反射面產生的反射波反相抵消，從而實現(xiàn)行波管的“半無限長”特性。

若輸入至主發(fā)換能器的信號為x，聲信號輸入主發(fā)換能器發(fā)射，經過待測聲學材料到達次發(fā)換能器反射面的傳遞函數(shù)為H主，入射至次發(fā)換能器反射面的聲信號為P主，有

其中，?表示卷積。假設P主到達次發(fā)換能器反射面產生反射聲信號P反的傳遞函數(shù)為H反，有

為了滿足行波場建立的條件，次發(fā)換能器發(fā)射的聲信號P次應滿足P反=?P次的關系，綜合以上各式可以得到P次為

對式(3)取傅里葉變換，則有

令輸入至次發(fā)換能器的抵消聲信號為x次，考慮信號輸入至次發(fā)換能器發(fā)射聲信號P次的傳遞函數(shù)為H次，有

可解得

其中，若令x次(ω)=x(ω)Aejωτ時延，并令M(ω)ejφ(ω)可得

其中，M(ω)和φ(ω)分別為建立行波場所需的最佳幅度譜與相位譜，而Aopt與τopt為次發(fā)換能器發(fā)射時輸入的最佳幅度與最佳時延。為了獲得最佳的Aopt與τopt從而逼近M(ω)與φ(ω)以提高抵消效果，本文采用了一維黃金分割搜索算法用于搜索每一個測量頻點的幅值A與時延τ時延，從而實現(xiàn)行波場的建立。

黃金分割算法通過以固定比率(即黃金分割率)減小搜索的范圍，進而可以逼近區(qū)間內的極值點。在輸入搜索時延或者幅值區(qū)間[a,b]后，由雙水聽器法得到此時的次發(fā)反射系數(shù)r1與r2，并在比較兩者的大小后，令其中較大的一方按0.618 的黃金分割倍率縮小區(qū)間范圍，然后開始新一輪的循環(huán)。當區(qū)間[a,b]的長度低于門限c之時，算法搜索完畢，具體流程如圖1所示。

圖1 黃金分割搜索算法原理框圖Fig.1 Golden section search algorithm functional block diagram

令黃金分割率為Φ，則算法的搜索次數(shù)n與區(qū)間[a,b]及門限c關系如下：

對于時延搜索而言，令fs為采樣率，f為當前頻點頻率，且令搜索區(qū)間為當前頻點下的半周期長的采樣點數(shù)，門限寬度設置為1個采樣點(即s)，則時延搜索次數(shù)n時延有

同理可得幅值搜索所需的次數(shù)。設置幅值搜索范圍時需考慮功放輸出許可范圍，且由于行波場建立對時延更為敏感，故實踐中可以略微放寬幅值搜索門限從而加快搜索次數(shù)?？紤]到換能器的頻響以及待測材料的聲學特性，實踐中應對每一測量頻點進行幅值與時延搜索。

2 數(shù)值仿真

下面對使用黃金分割搜索算法建立行波場的方法進行數(shù)值仿真。建立如圖2所示的行波管模型，使用主發(fā)換能器發(fā)射單頻連續(xù)(Continuous wave,CW)正弦脈沖信號，該信號經由次發(fā)換能器界面反射，并使用次發(fā)換能器發(fā)射抵消信號對回波進行抵消，使用雙水聽器法對回波進行分離。圖2 為數(shù)值仿真位置關系示意圖，行波管長L=4.5 m，仿真中假設聲信號在聲管中的信道為純時延信道，水聽器的靈敏度為1，并假設兩端換能器與外圍電路的頻響是由隨機數(shù)組成的32 階有限脈沖響應(Finite impulse response,FIR)濾波器。

圖2 行波管仿真模型Fig.2 Traveling wave tube simulation model

圖3 為主發(fā)換能器發(fā)射1.5 kHz 的CW 正弦脈沖信號，此時次發(fā)換能器未發(fā)射抵消信號，聲信號由次發(fā)換能器界面反射。通過雙水聽器法解算得到的次發(fā)換能器反射界面處的入射信號與反射信號。圖4(a)為未考慮時延與幅度影響，次發(fā)換能器直接使用τ理論抵消反射波從而進行行波場建立的結果?？梢娪捎趽Q能器與外圍電路的頻響影響導致的時延偏移以及幅值偏移使得行波場建立失敗。使用一維黃金分割搜索算法對時延與幅值進行搜索，行波場建立結果如圖4(b)所示，搜索過程中的次發(fā)表面反射系數(shù)如圖5 所示?？梢娛褂命S金分割搜索算法成功地搜索得到正確的時延與幅值，從而完成了行波場的建立。

圖3 1.5 kHz 時主發(fā)信號時域波形Fig.3 Time domain waveform of main signal at 1.5 kHz

圖4 搜索前后行波場建立結果Fig.4 Establishment of traveling wave field before and after search

圖5 反射系數(shù)搜索曲線Fig.5 Reflection coefficient search curve

3 試驗驗證

為了驗證算法和系統(tǒng)性能，作者在中國船舶工業(yè)系統(tǒng)工程研究院的行波管系統(tǒng)中進行了試驗，使用LabView 編寫了試驗系統(tǒng)。該試驗系統(tǒng)如圖6所示。

以水柱為對象，使用本文算法在1.5 kHz 頻率下進行時延與幅值搜索，以驗證使用黃金分割算法建立行波場的效果。首先在固定幅值時進行時延搜索，得到最佳時延，并使用該時延進行幅值搜索，得到最佳幅值。1.5 kHz 頻率下搜索時延與幅值時次發(fā)換能器處反射系數(shù)r隨搜索次數(shù)變化的情況如圖7 所示，其中橫軸表示搜索次數(shù)，縱軸表示次發(fā)換能器處反射系數(shù)r。在給定初始幅值的情況下，以τ理論的正負半周期范圍內進行最佳時延的搜索，隨著搜索的進行，不斷縮小搜索區(qū)間，當搜索區(qū)間寬度小于門限值后，時延搜索停止，得到最佳時延。然后使用得到的最佳時延，在允許輸出范圍內進行幅值搜索，以得到最佳幅值?？梢娫谒阉鲿r延的情況下，次發(fā)表面的反射系數(shù)r迅速下降，而搜索幅值的情況下反射系數(shù)r下降緩慢，說明此時行波場建立已經趨于收斂。當搜索結束后，次發(fā)反射系數(shù)r約等于0.05，可以認為此時行波場已經建立。

圖6 行波管測量系統(tǒng)框圖Fig.6 Traveling wave tube measurement system block diagram

圖7 1.5 kHz 下次發(fā)表面反射系數(shù)搜索結果Fig.7 Search result of interface reflection coefficient of secondary transducer at 1.5 kHz

對水柱進行1 kHz～3 kHz 范圍內以100 Hz 為間隔的聲學參數(shù)測量，結果如圖8所示，其中橫軸表示頻率點，縱軸表示無量綱系數(shù)?？梢娡干湎禂?shù)與水柱的理論值符合較好，但是反射系數(shù)測量誤差較大，其原因可能與雙水聽器法誤差有關，產生雙水聽器法誤差的原因包括傳遞函數(shù)的測量誤差與反射、透射系數(shù)的計算誤差與聲速測量誤差[6]。

圖8 水柱聲學參數(shù)測量結果Fig.8 Measurement result of acoustic parameters of water cylinder

對30 mm 厚度的鋼板進行1 kHz～3 kHz 范圍內以100 Hz為間隔的聲學參數(shù)測量，并與理論值進行對比，結果如圖9 所示，其中橫軸表示頻率點，縱軸表示無量綱系數(shù)。可見，透射系數(shù)t測量結果與理論值符合得較好，說明本系統(tǒng)測量透射系數(shù)較為準確，但是反射系數(shù)r比理論值偏小，且曲線不光滑，測量反射系數(shù)誤差較大。將t與r的測量值與理論值的相對誤差做柱狀圖如圖10 所示，其中橫軸表示頻率點，縱軸表示無量綱相對誤差?？梢?，除了少數(shù)頻點，透射系數(shù)的相對誤差都低于10%，但是反射系數(shù)的誤差較大。反射系數(shù)誤差較大、透射系數(shù)測量誤差較小，這與水柱聲學參數(shù)測量結果趨勢一致。

圖9 30 mm 鋼板聲學參數(shù)理論值與測量值對比Fig.9 Comparison between theoretical and measured acoustic parameters of 30 mm steel plate

圖10 30 mm 鋼板測量相對誤差Fig.10 Measurement relative error of 30 mm steel plate

4 結論

本文通過一維黃金分割搜索算法，對幅度A與時延τ時延進行搜索，實現(xiàn)CW 脈沖信號的行波場建立。該方法能夠快速建立行波場，從而實現(xiàn)聲學參數(shù)測量，但是該方法并不能對發(fā)射波形進行精細的修改，從而更加匹配建立最佳行波場所需的最佳幅度譜M與相位譜φ，進一步的工作可以考慮對發(fā)射波形進行均衡，從而提高行波場建立效果。本文使用該算法進行了行波場建立，并使用測量系統(tǒng)對水柱以及鋼板的聲學參數(shù)進行了測量，驗證了該系統(tǒng)的可行性。具體如下：

(1)使用黃金分割算法能夠得到抵消反射波所需的時延與幅值，實現(xiàn)有源吸聲，完成行波場的建立。

(2)對鋼板和水柱的測量結果表明，本系統(tǒng)聲學參數(shù)測量與理論值符合較好。