一種基于改進模擬退火算法的極限稀疏陣列優(yōu)化方法
【技術(shù)領(lǐng)域】
[0001] 本發(fā)明涉及聲納換能器陣列設(shè)計領(lǐng)域，尤其涉及一種基于改進模擬退火算法的極 限稀疏陣列優(yōu)化方法。
【背景技術(shù)】
[0002] 全布陣換能器陣列是一種采用均勻、周期分布的換能器陣列。隨著海洋資源探測、 開發(fā)等領(lǐng)域的迅速發(fā)展，其對水下三維聲學(xué)成像系統(tǒng)在分辨率、系統(tǒng)復(fù)雜度、功耗以及成本 等方面的要求也越來越高，傳統(tǒng)的全布陣換能器陣列已經(jīng)難以滿足其日益增長的性能需 求。為解決上述問題，眾多學(xué)者采用一種非等間距陣列（Unequally spaced arrays)，即稀 疏陣列（Sparse arrays)，代替全布陣換能器陣列進行三維聲學(xué)成像，從而降低系統(tǒng)的復(fù)雜 度、功耗以及成本。稀疏陣列通過去除換能器陣列中某些陣元，并對保留陣元的位置和權(quán)重 系數(shù)進行優(yōu)化，從而大幅減小了換能器陣列陣元數(shù)量，同時也將其波束方向圖的主瓣寬度 和旁瓣峰值控制在預(yù)期范圍內(nèi)。
[0003] 根據(jù)優(yōu)化方法的不同，稀疏陣列的優(yōu)化算法主要可以分為：確定優(yōu)化算法和隨機 優(yōu)化算法。確定優(yōu)化算法相比于隨機優(yōu)化算法，具有效率高、迭代少、收斂快等優(yōu)勢，但其在 稀疏優(yōu)化過程中，對陣列的形狀有較高的要求，因此通常未必能夠?qū)さ萌肿顑?yōu)解。隨機 優(yōu)化算法首先基于陣列的主瓣寬度、旁瓣峰值以及波束方向圖形狀等指標構(gòu)建目標函數(shù)， 然后按照一定的隨機準則對目標函數(shù)進行多次迭代優(yōu)化，使其能夠逐步收斂到全局最優(yōu)狀 態(tài)，因此與確定優(yōu)化算法相比得到了更為廣泛的應(yīng)用。隨機優(yōu)化算法主要包含遺傳算法、粒 子群優(yōu)化算法和模擬退火算法三大類。其中，遺傳算法和粒子群優(yōu)化算法對于小規(guī)模的陣 列稀疏問題，能夠有效地進行優(yōu)化，但對于當(dāng)前應(yīng)用中越來越大的陣列規(guī)模，這兩種算法已 難以滿足其需求。模擬退火算法能夠同時對稀疏陣列的陣元位置和權(quán)重系數(shù)進行優(yōu)化，并 且換能器陣列的規(guī)模尺寸對其全局尋優(yōu)性能影響不大，因此在陣列稀疏優(yōu)化問題中，其應(yīng) 用更加廣泛。
[0004] 模擬退火算法可以用于解決各種優(yōu)化問題，尋求全局最優(yōu)解。該算法是根據(jù)固體 的退火原理衍生而來的，當(dāng)固體從充分高的溫度緩慢冷卻時，其內(nèi)部粒子會從無序狀態(tài)慢 慢達到有序排列的"結(jié)晶"狀態(tài)，其內(nèi)能也隨之逐漸減小，并最終達到最低。模擬退火算法 具備全局尋優(yōu)能力，其利用玻爾茲曼概率分布（Boltzmann Probability Distribution)跳 出局部極值區(qū)域，增大了尋求全局最優(yōu)解的概率。
[0005] 模擬退火算法作為一種啟發(fā)式的方法，適用于各種優(yōu)化問題。其在優(yōu)化過程中，能 夠?qū)さ萌肿顑?yōu)解，且收斂速度快，因此被廣泛應(yīng)用。采用模擬退火算法進行稀疏陣列優(yōu)化 時，需要定義一個目標能量函數(shù)（Energy Function)，通過多次的隨機迭代過程，尋求使能 量函數(shù)達到最小值的狀態(tài)，該過程屬于一種馬爾可夫鏈（Markov Chain)方法。

【發(fā)明內(nèi)容】

[0006] 本發(fā)明針對當(dāng)前稀疏陣列陣元數(shù)依舊偏大，難以滿足水下實時三維聲學(xué)成像系統(tǒng) 低功耗、小型化需求的問題，提供了一種基于改進模擬退火算法的極限稀疏陣列優(yōu)化方法。 該方法基于多頻率發(fā)射算法對模擬退火算法進行了改進，重新定義了其能量函數(shù)，以十字 型陣列為目標，對其進一步稀疏優(yōu)化，獲得了陣元數(shù)量最小化的極限稀疏陣列。該稀疏陣列 能夠以更少的陣元數(shù)量，獲得與其他稀疏陣列相同的波束方向圖指標（主瓣寬度、旁瓣峰 值等），最大程度降低了實時三維聲學(xué)成像系統(tǒng)的硬件復(fù)雜度。
[0007] -種基于改進模擬退火算法的極限稀疏陣列優(yōu)化方法，包括步驟：
[0008] (1)對十字型陣列的發(fā)射波束形成進行優(yōu)化，并利用十字型陣列的波束方向圖重 新定義模擬退火算法中的能量函數(shù)E(W，A);
[0010] 式中，W為十字型陣列陣元的權(quán)重系數(shù)矩陣，Β(Ι，ιι'，ν'，Α)為波束方向圖的表 達式，max (B (W，u'，ν'，f ρ)為B (W，u'，ν'，f J的最大值，h，k2,匕表示三項優(yōu)化目標參數(shù) 的權(quán)重，SLPd為預(yù)期的旁瓣峰值，A為十字型陣列陣元數(shù)量（發(fā)射和接收陣元數(shù)量總和）， 私為當(dāng)前狀態(tài)的陣元權(quán)重系數(shù)比，Rd為預(yù)期的陣元權(quán)重系數(shù)比，u'，v'分別定義為u' = sin a p-sin a Q，ν' = sin β q-sin β。，其中（α ρ, β q)為水平和垂直波束方向角，（α。，β Q)分 別為指向聲源的單位矢量的水平和垂直方向角，Ω代表（u'，v'）的取值情況集合；
[0011] 優(yōu)選地，Ω的取值滿足如下表達式：
[0013] (2)設(shè)置預(yù)期的旁瓣峰值SLPd(通常設(shè)置為-22dB)和陣元權(quán)重系數(shù)比R d(通常設(shè) 置為3)，對十字型陣列進行稀疏優(yōu)化，使能量函數(shù)達到最低值，獲得陣元數(shù)量最小化的極限 稀疏陣列。
[0014] 優(yōu)選的，對十字型陣列的發(fā)射波束形成進行優(yōu)化的步驟包括：
[0015] (1. 1)將十字型陣列中發(fā)射陣列的波束方向分割為多個扇面，在每個扇面內(nèi)依次 發(fā)射一系列不同頻率的扇形聲納波束信號，每個頻率的扇形聲納波束信號指向?qū)?yīng)扇面內(nèi) 的一個波束方向；
[0016] (1. 2)每個扇面內(nèi)所有頻率的扇形聲納波束信號發(fā)射結(jié)束后，利用十字型陣列中 的接收陣列接收聲納回波信號，通過離散傅里葉變換抽取各扇面內(nèi)所有扇形聲納波束信號 對應(yīng)的頻率信息，并在頻率信息對應(yīng)的頻域內(nèi)進行波束形成計算。
[0017] 本發(fā)明采用多頻率發(fā)射波束形成算法，將十字型陣列的Q個預(yù)設(shè)的垂直波束方向 分割為K個扇面，在每個扇面內(nèi)，發(fā)射換能器陣列通過各陣元間的相移補償，依次向預(yù)設(shè)的 J個垂直波束方向（Q = KXJ)發(fā)射不同頻率的扇形聲納波束信號，每個頻率的信號對應(yīng)一 個垂直波束方向；然后，當(dāng)該扇面內(nèi)所有頻率的聲納波束信號發(fā)射結(jié)束后，接收陣列收到聲 納回波信號，通過DFT運算，同時抽取回波中J個發(fā)射聲納波束信號對應(yīng)的頻率信息，并行 地在J個頻域上進行波束形成計算，生成P (水平波束方向數(shù)）XJ個波束強度結(jié)果；之后， 對其余扇面進行類似的處理。當(dāng)所有扇面完成上述操作后，則可得到完整的PXQ個方向的 波束強度結(jié)果。
[0018] 基于上述多頻率發(fā)射波束形成算法，縮短了十字型陣列的發(fā)射時間，提高了其實 時性；綜合考慮稀疏陣列波束方向圖的旁瓣峰值以及陣元的權(quán)重系數(shù)比，來重新定義模擬 退火算法中的能量函數(shù)E (W, A)。
[0019] 優(yōu)選的，所述十字型陣列的波束方向圖的表達式為：
[0021] 式中，BTrMFT(v'，fj)和BReMFT(u'，fj)分別為十字型陣列的發(fā)射波束方向圖表達式和 接收波束方向圖表達式，U'，V'分別定義為u' =sinap-sinaQ，v' =sin0q-sin0。。
[0022] 優(yōu)選的，對十字型陣列進行稀疏優(yōu)化的具體過程包括：
[0023] (2. 1)初始化所述的陣元權(quán)重系數(shù)矩陣W，并設(shè)置一個初始系統(tǒng)"溫度";作為優(yōu)選， 系統(tǒng)初始"溫度"通常設(shè)置在1〇〇〇以上；
[0024] (2. 2)開始迭代過程，每次對迭代變量引入一個隨機的微擾，若該微擾狀態(tài)使得能 量函數(shù)值降低，則接受該狀態(tài)，并更新參數(shù)配置；若該微擾狀態(tài)使得能量函數(shù)值增大，則根 據(jù)波爾茲曼概率分布決定該狀態(tài)是否被接受，概率分布由系統(tǒng)"溫度"決定，系統(tǒng)"溫度"越 尚，接受該微擾狀態(tài)的概率也越大；
[0025] (2.3)在當(dāng)前系統(tǒng)"溫度"下，所有迭代變量均被訪問過一次后，降低"溫度"，并重 復(fù)上述迭代過程；
[0026] (2. 4)隨著系統(tǒng)"溫度"不斷降低，能量函數(shù)值不再隨迭代次數(shù)的增加而下降，當(dāng)達 到結(jié)束條件后，停止迭代過程，此時能量函數(shù)達到最低值，當(dāng)前的陣列參數(shù)配置即為最終的 稀疏陣列優(yōu)化結(jié)果。
[0027] 優(yōu)選的，在第1次迭代中，隨機選取一個十字型陣列的陣元，根據(jù)陣元的狀態(tài)進行 如下操作：
[0028] (2. 2. 1)若選取的陣元權(quán)重系數(shù)為0,即其處于關(guān)閉狀態(tài)，則根據(jù)一個固定的重生 概率開啟該陣元，并賦予其一個隨機的權(quán)重系數(shù)ω ，同時更新陣元權(quán)重系數(shù)矩陣W和 陣元數(shù)量A ;若該陣元