專利名稱:一種基于面元散射的合成孔徑聲納面目標(biāo)快速仿真方法
技術(shù)領(lǐng)域:
本發(fā)明屬于聲納信號(hào)處理領(lǐng)域,特別涉及一種基于面元散射的合成孔徑聲納面目標(biāo)快速仿真方法。
背景技術(shù):
合成孔徑聲納(Synthetic Aperture SonarSAS)是一種先進(jìn)的高分辨率水聲成像聲納,其基本原理是利用小孔徑基陣在方位向的移動(dòng)形成虛擬大孔徑,通過(guò)對(duì)不同位置的聲納回波進(jìn)行相干處理,從而獲得方位向的高分辨率。
面目標(biāo)的回波數(shù)據(jù)仿真是進(jìn)行合成孔徑聲納信號(hào)處理研究的重要途徑。利用仿真的面目標(biāo)回波數(shù)據(jù)可以驗(yàn)證各種SAS成像算法、運(yùn)動(dòng)補(bǔ)償及自聚焦算法的有效性,同時(shí)還可以分析系統(tǒng)性能,對(duì)系統(tǒng)參數(shù)的選擇具有指導(dǎo)作用,此外水下目標(biāo)的識(shí)別更是需要大量的數(shù)據(jù)以建立目標(biāo)匹配模型。固然SAS的回波數(shù)據(jù)可以通過(guò)湖海試獲得,但其成本過(guò)高且特殊的場(chǎng)景不易實(shí)施,面目標(biāo)的回波仿真被國(guó)內(nèi)外研究機(jī)構(gòu)所關(guān)注。
面目標(biāo)的表征主要包括點(diǎn)模型和面元模型兩種。基于點(diǎn)模型表征的面目標(biāo),其回波數(shù)據(jù)仿真原理簡(jiǎn)單、易于操作,其缺點(diǎn)是仿真運(yùn)算量龐大。基于面元模型表征的面目標(biāo),利用赫姆霍茲-基爾霍夫積分公式可計(jì)算出單個(gè)小面元的聲散射場(chǎng),對(duì)所有小面元聲散射場(chǎng)進(jìn)行疊加處理結(jié)合聲源信號(hào)可計(jì)算出面目標(biāo)對(duì)應(yīng)的回波數(shù)據(jù)。針對(duì)三角面元的赫姆霍茲-基爾霍夫積分公式,一般通過(guò)通過(guò)兩種技術(shù)途徑解決第一,利用空間變換將任一三角面元變換至二維平面上,在此基礎(chǔ)上計(jì)算其面積分;第二,利用仿射變換將任一三角面元的散射場(chǎng)通過(guò)對(duì)單位直角三角面元散射場(chǎng)的插值處理獲得。第二種方法相對(duì)于第一種方法提高了計(jì)算效率,然而由于插值處理的存在,仿真精度有待進(jìn)一步提高。
發(fā)明內(nèi)容
本發(fā)明的目的在于,為克服現(xiàn)有的插值信號(hào)處理以及面元形狀選擇造成的仿真精度低等缺點(diǎn),從而提出一種基于面元散射的合成孔徑聲納面目標(biāo)快速仿真方法。
為了解決上述問(wèn)題,本發(fā)明的目的在于提供一種基于面元散射的合成孔徑聲納面目標(biāo)快速仿真方法。該方法利用格林公式結(jié)合三角面元的三條邊參量方程,將赫姆霍茲-基爾霍夫積分公式進(jìn)一步化簡(jiǎn),在此基礎(chǔ)上實(shí)現(xiàn)合成孔徑聲納面目標(biāo)的快速仿真。
一種基于面元散射的合成孔徑聲納面目標(biāo)快速仿真方法,該方法用格林公式結(jié)合三角面元的三條邊參量方程,化簡(jiǎn)赫姆霍茲-基爾霍夫積分公式,,實(shí)現(xiàn)合成孔徑聲納面目標(biāo)的快速仿真,所述的方法包含如下步驟 (1)設(shè)定合成孔徑聲納系統(tǒng)的仿真參數(shù),構(gòu)造面目標(biāo)場(chǎng)景模型; 其中,所述的合成孔徑聲納系統(tǒng)的仿真參數(shù)主要包括中心頻率、信號(hào)帶寬、脈沖寬度、脈沖重復(fù)周期、采樣頻率、發(fā)射陣孔徑、接收陣孔徑、陣元個(gè)數(shù)、采樣點(diǎn)數(shù)、最小距離、聲速或拖體速度;所述的構(gòu)造面目標(biāo)場(chǎng)景模型,通過(guò)ANSYS仿真分析創(chuàng)建面目標(biāo)模型,并對(duì)其進(jìn)行三角面元的網(wǎng)格劃分。
(2)依據(jù)聲納基陣在方位向上的時(shí)間順序,確定基陣位置,針對(duì)確定的基陣位置,確定表征面目標(biāo)的三角面元被照射狀態(tài),所述的三角面元的狀態(tài)分為完全被照射、完全未被照射和部分被照射三種狀態(tài); (3)對(duì)所有完全被照射的三角面元,用格林公式結(jié)合每個(gè)三角面元的三條邊參量方程化簡(jiǎn)赫姆霍茲-基爾霍夫積分公式,再用化簡(jiǎn)得到的公式計(jì)算所有被照射三角面元的散射場(chǎng)的頻域信號(hào),對(duì)得到的所有三角面元的散射場(chǎng)的頻域信號(hào)進(jìn)行疊加,最后進(jìn)行逆傅里葉變換獲得當(dāng)前位置對(duì)應(yīng)的面目標(biāo)回波數(shù)據(jù); (4)按照方位向的順序重復(fù)步驟(2)~(4),從而獲得當(dāng)前基陣位置對(duì)應(yīng)的回波信號(hào)的頻域表示,再對(duì)頻域信號(hào)進(jìn)行逆傅里葉變換獲得面目標(biāo)對(duì)應(yīng)的完整回波數(shù)據(jù),實(shí)現(xiàn)合成孔徑聲納面目標(biāo)的快速仿真。
上述技術(shù)方案采用格林公式結(jié)合三角面元的三條邊參量方程,得到化簡(jiǎn)的赫姆霍茲-基爾霍夫積分公式如下 其中,
為三角面元的單位法線向量,P為發(fā)射陣位置,H為接收陣位置,
為單位入射向量,
為單位散射向量,記
為
在三角面元上的投影,
對(duì)上式求其逆傅里葉變換可以得到該三角面元對(duì)應(yīng)聲納基陣位置的時(shí)域回波數(shù)據(jù)。
作為本發(fā)明的一個(gè)改進(jìn),所述的三角面元的三種照射狀態(tài),通過(guò)下式判斷
作為本發(fā)明的又一改進(jìn),所述的部分被照射的三角面元,通過(guò)進(jìn)一步劃分為若干完全被照射和完全未被照射的小三角面元進(jìn)行處理,具體劃分方法步驟如下 對(duì)于部分被照射三角面元,按照其處于波束照射內(nèi)的頂點(diǎn)個(gè)數(shù)分為兩類,其中,對(duì)于一個(gè)頂點(diǎn)在被照射區(qū)域的三角面元Tk,假設(shè)為Vk2處于波束內(nèi)時(shí),A1和A2為波束多棱錐和三角面元Tk的兩個(gè)交點(diǎn),則Tk可進(jìn)一步劃分成一個(gè)完全被照射小三角面元Tk1和兩個(gè)完全未被照射的小三角面元Tk2、Tk3,其中Tk1由頂點(diǎn)Vk2、A1、A2構(gòu)成,Tk2由頂點(diǎn)Vk1、A1、A2構(gòu)成,Tk3由頂點(diǎn)Vk1、Vk3、A2構(gòu)成;對(duì)于兩個(gè)頂點(diǎn)在被照射區(qū)域的三角面元Tk,假設(shè)有兩個(gè)頂點(diǎn)Vk2、Vk3處于波束內(nèi)時(shí),A1和A2為波束多棱錐和三角面元Tk的兩個(gè)交點(diǎn),則Tk可進(jìn)一步劃分成一個(gè)完全為被照射小三角面元Tk1和兩個(gè)完全被照射的小三角面元Tk2、Tk3,其中Tk1由頂點(diǎn)Vk1、A1、A2構(gòu)成,Tk2由頂點(diǎn)Vk2、Vk3、A1構(gòu)成,Tk3由頂點(diǎn)Vk3、A1、A2構(gòu)成。所述的劃分得到的若干完全被照射小三角面元依據(jù)完全被照射三角面元處理方法得到回波數(shù)據(jù)。
所述的基陣位置,依據(jù)方位向上對(duì)應(yīng)的時(shí)間順序和拖體速度確定聲納基陣的位置,其中當(dāng)需要驗(yàn)證運(yùn)動(dòng)補(bǔ)償和自聚焦算法時(shí),拖體的速度可以根據(jù)實(shí)際需要進(jìn)行設(shè)置。
依據(jù)方位向上對(duì)應(yīng)的時(shí)間順序和拖體速度可以確定聲納基陣的位置。合成孔徑聲納成像的理想航跡為拖體沿方位向方向做勻速直線運(yùn)動(dòng),當(dāng)需要驗(yàn)證運(yùn)動(dòng)補(bǔ)償和自聚焦算法時(shí),拖體的速度可以根據(jù)實(shí)際需要進(jìn)行設(shè)置。上述技術(shù)方案所述的聲納基陣的發(fā)射波束由多棱錐近似模型表示,依據(jù)三角面元質(zhì)心指向基陣的向量和三角面元的法線向量的空間關(guān)系,結(jié)合發(fā)射波束,可以將表征面目標(biāo)的三角面元分為完全被照射、部分被照射和完全未被照射三種情況。針對(duì)部分被照射的三角面元依據(jù)發(fā)射波束在面目標(biāo)的投影可進(jìn)一步拆分成完全被照射和完全未被照射的小三角面元。
針對(duì)赫姆霍茲-基爾霍夫積分公式,利用格林公式和三角面元的三條表的參量方程可以將其進(jìn)一步化簡(jiǎn)。對(duì)所有完全被照射三角面元依據(jù)化簡(jiǎn)后的赫姆霍茲-基爾霍夫積分公式,逐一求出其散射場(chǎng),然后進(jìn)行疊加處理,從而獲得當(dāng)前基陣位置對(duì)應(yīng)的回波信號(hào)的頻域表示,再對(duì)頻域信號(hào)進(jìn)行逆傅里葉變換獲得回波數(shù)據(jù)。
本發(fā)明的優(yōu)點(diǎn)主要有 (1)計(jì)算單個(gè)三角面元回波時(shí),不存在插值處理,結(jié)果更為精確,計(jì)算效率更高; (2)針對(duì)面目標(biāo)的被照射區(qū)域的準(zhǔn)確獲得,進(jìn)一步提高仿真精度; (3)利用已有的計(jì)算機(jī)輔助設(shè)計(jì)軟件(ANSYS等),場(chǎng)景生成便捷。
圖1是本發(fā)明的合成孔徑聲納面目標(biāo)快速仿真流程圖; 圖2是本發(fā)明的合成孔徑聲納面目標(biāo)回波模型; 圖3是三角面元被照射狀態(tài)示意圖; 圖4是本發(fā)明的一個(gè)頂點(diǎn)處于波束內(nèi)的三角面元示意圖; 圖5是本發(fā)明的兩個(gè)頂點(diǎn)處于波束內(nèi)的三角面元示意圖; 圖6是本發(fā)明的三角面元散射計(jì)算模型。
具體實(shí)施例方式 下面結(jié)合附圖1~4對(duì)本發(fā)明即“一種基于面元散射的合成孔徑聲納面目標(biāo)快速仿真方法”進(jìn)行詳細(xì)的說(shuō)明。
圖1是本發(fā)明即“一種基于面元散射的合成孔徑聲納面目標(biāo)快速仿真方法”的流程圖,該方法由設(shè)置仿真參數(shù)、生成面目標(biāo)場(chǎng)景、確定聲納基陣、確定目標(biāo)面元被照射狀態(tài)、對(duì)完全被照射面元計(jì)算其頻域響應(yīng)及對(duì)應(yīng)的回波信號(hào)等步驟組成。
如圖2所示,設(shè)y方向?yàn)楹铣煽讖铰暭{成像的方位向,x方向?yàn)楹铣煽讖铰暭{成像的距離向,z方向?yàn)榫嗟咨疃?。設(shè)置合成孔徑聲納面目標(biāo)的仿真參數(shù),主要包括中心頻率、信號(hào)帶寬、脈沖寬度、脈沖重復(fù)周期、采樣頻率、發(fā)射陣孔徑、接收陣孔徑、陣元個(gè)數(shù)、采樣點(diǎn)數(shù)、最小距離、聲速、拖體速度等。利用ANSYS、AutoCAD、3dMAX等計(jì)算機(jī)輔助設(shè)計(jì)軟件創(chuàng)建面目標(biāo)模型,并對(duì)其進(jìn)行三角面元網(wǎng)格劃分,保存面目標(biāo)對(duì)應(yīng)的三角面元網(wǎng)格信息,生成面目標(biāo)場(chǎng)景文件。
合成孔徑聲納理想航跡為拖體沿方位向做勻速直線運(yùn)動(dòng),當(dāng)需要仿真存在運(yùn)動(dòng)誤差的面目標(biāo)回波數(shù)據(jù)時(shí)可以根據(jù)實(shí)際需要設(shè)定拖體速度v(vx,vy,vz)。設(shè)發(fā)射陣發(fā)射脈沖的重復(fù)周期為prt,拖體上裝載的聲納基陣起始坐標(biāo)為(0,0,H),其中H為距底深度,則第n個(gè)采樣時(shí)刻聲納基陣坐標(biāo)為(n·vx·prt,n·vy·prt,H+n·vz·prt)。
對(duì)于面目標(biāo)場(chǎng)景中由點(diǎn)Vk1、Vk2和Vk3構(gòu)成的三角面元Tk,其單位法向量為
質(zhì)心為Ck,設(shè)
表示由質(zhì)心點(diǎn)Ck指向發(fā)射陣的單位向量,Bi為發(fā)射陣所處位置對(duì)應(yīng)的發(fā)射波束照射范圍(發(fā)射波束引用多棱錐近似模型表示),對(duì)于三角面元Tk則有
如圖3所示,圖中分別示意了完全被照射、完全未被照射以及部分被照射的三個(gè)三角面元,針對(duì)部分被照射的三角面元,其處于波束內(nèi)的頂點(diǎn)個(gè)數(shù)只能為1和2,分別是圖4和5所示。如圖4所示,對(duì)于三角面元Tk,當(dāng)僅有一個(gè)頂點(diǎn)假設(shè)為Vk2處于波束內(nèi)時(shí),A1和A2為波束多棱錐和三角面元Tk的兩個(gè)交點(diǎn),則Tk可進(jìn)一步劃分成一個(gè)完全被照射三角面元Tk1和兩個(gè)完全未被照射的三角面元Tk2、Tk3,其中Tk1由頂點(diǎn)Vk2、A1、A2構(gòu)成,Tk2由頂點(diǎn)Vk1、A1、A2構(gòu)成,Tk3由頂點(diǎn)Vk1、Vk3、A2構(gòu)成。如圖5所示,對(duì)于三角面元Tk,當(dāng)有兩個(gè)頂點(diǎn)假設(shè)為Vk2、Vk3處于波束內(nèi)時(shí),A1和A2為波束多棱錐和三角面元Tk的兩個(gè)交點(diǎn),則Tk可進(jìn)一步劃分成一個(gè)完全為被照射三角面元Tk1和兩個(gè)完全被照射的三角面元Tk2、Tk3,其中Tk1由頂點(diǎn)Vk1、A1、A2構(gòu)成,Tk2由頂點(diǎn)Vk2、Vk3、A1構(gòu)成,Tk3由頂點(diǎn)Vk3、A1、A2構(gòu)成。
如圖6所示,對(duì)于任一三角面元,其散射場(chǎng)可以由赫姆霍茲-基爾霍夫積分公式計(jì)算如下 其中R反射系數(shù),G(H|M)為自由空間的格林函數(shù),定義如下 圖6中,
為三角面元的單位法線向量,P為發(fā)射陣位置,H為接收陣位置,
為單位入射向量,
為單位散射向量,記
為
在三角面元上的投影,
利用格林公式將赫姆霍茲-基爾霍夫積分公式由面積分化簡(jiǎn)為線積分如下
對(duì)于上式中的線積分項(xiàng),針對(duì)三角面元的第i條邊可用參量方程表示如下 li=(1-t)ai+tai+1 t∈
則有 針對(duì)第i條邊的線積分可表示為 可以進(jìn)一步推出三角面元對(duì)應(yīng)的回波的頻域表示為 對(duì)上式求其逆傅里葉變換可以得到對(duì)應(yīng)聲納基陣位置的面目標(biāo)時(shí)域回波數(shù)據(jù)。依據(jù)方位向的時(shí)間順序依次重復(fù)以上步驟可以得出面目標(biāo)對(duì)應(yīng)的完整回波數(shù)據(jù),實(shí)現(xiàn)合成孔徑聲納面目標(biāo)的快速仿真。
需要說(shuō)明的是,以上介紹的本發(fā)明的實(shí)施方案而并非限制。本領(lǐng)域的技術(shù)人員應(yīng)當(dāng)理解,任何對(duì)本發(fā)明技術(shù)方案的修改或者等同替代都不脫離本發(fā)明技術(shù)方案的精神和范圍,其均應(yīng)涵蓋在本發(fā)明的權(quán)利要求范圍內(nèi)。
權(quán)利要求
1.一種基于面元散射的合成孔徑聲納面目標(biāo)快速仿真方法,該方法用格林公式結(jié)合三角面元的三條邊參量方程,化簡(jiǎn)赫姆霍茲-基爾霍夫積分公式,實(shí)現(xiàn)合成孔徑聲納面目標(biāo)的快速仿真,所述的方法包含如下步驟
(1)設(shè)定合成孔徑聲納系統(tǒng)的仿真參數(shù),構(gòu)造面目標(biāo)場(chǎng)景模型;
(2)依據(jù)聲納基陣在方位向上的時(shí)間順序,確定基陣位置,針對(duì)確定的基陣位置,確定表征面目標(biāo)的三角面元被照射狀態(tài),所述的三角面元被照射狀態(tài)包括完全被照射、完全未被照射或部分被照射三種狀態(tài);
(3)對(duì)所有完全被照射的三角面元,用格林公式結(jié)合每個(gè)三角面元的三條邊參量方程化簡(jiǎn)赫姆霍茲-基爾霍夫積分公式,再用化簡(jiǎn)得到的公式計(jì)算所有被照射三角面元的散射場(chǎng)的頻域信號(hào),對(duì)得到的頻域信號(hào)進(jìn)行疊加,最后進(jìn)行逆傅里葉變換獲得當(dāng)前位置對(duì)應(yīng)的面目標(biāo)回波數(shù)據(jù);
(4)按照方位向的順序重復(fù)步驟(2)~(4),從而獲得當(dāng)前基陣位置對(duì)應(yīng)的回波信號(hào)的頻域表示,再對(duì)頻域信號(hào)進(jìn)行逆傅里葉變換獲得面目標(biāo)對(duì)應(yīng)的完整回波數(shù)據(jù),實(shí)現(xiàn)合成孔徑聲納面目標(biāo)的快速仿真。
2.根據(jù)權(quán)利要求1所述的基于面元散射的合成孔徑聲納面目標(biāo)快速仿真方法,其特征在于,所述的步驟(3)中的化簡(jiǎn)的赫姆霍茲-基爾霍夫積分公式如下
其中,
為三角面元的單位法線向量,P為發(fā)射陣位置,H為接收陣位置,
為單位入射向量,
為單位散射向量,記
為
在三角面元上的投影,
對(duì)上式求其逆傅里葉變換可以得到該三角面元對(duì)應(yīng)聲納基陣位置的時(shí)域回波數(shù)據(jù)。
3.根據(jù)權(quán)利要求1所述的基于面元散射的合成孔徑聲納面目標(biāo)快速仿真方法,其特征在于,所述的步驟(2)中三角面元的三種照射狀態(tài),通過(guò)下式判斷
4.根據(jù)權(quán)利要求1或3所述的基于面元散射的合成孔徑聲納面目標(biāo)快速仿真方法,其特征在于,所述的部分被照射的三角面元,通過(guò)進(jìn)一步劃分為若干完全被照射和完全未被照射的小三角面元進(jìn)行處理,具體劃分方法步驟如下
對(duì)于部分被照射三角面元,按照其處于波束照射內(nèi)的頂點(diǎn)個(gè)數(shù)分為兩類,其中,對(duì)于一個(gè)頂點(diǎn)在被照射區(qū)域的三角面元Tk,假設(shè)為Vk2處于波束內(nèi)時(shí),A1和A2為波束多棱錐和三角面元Tk的兩個(gè)交點(diǎn),則Tk可進(jìn)一步劃分成一個(gè)完全被照射小三角面元Tk1和兩個(gè)完全未被照射的小三角面元Tk2、Tk3,其中Tk1由頂點(diǎn)Vk2、A1、A2構(gòu)成,Tk2由頂點(diǎn)Vk1、A1、A2構(gòu)成,Tk3由頂點(diǎn)Vk1、Vk3、A2構(gòu)成;對(duì)于兩個(gè)頂點(diǎn)在被照射區(qū)域的三角面元Tk,假設(shè)有兩個(gè)頂點(diǎn)Vk2、Vk3處于波束內(nèi)時(shí),A1和A2為波束多棱錐和三角面元Tk的兩個(gè)交點(diǎn),則Tk可進(jìn)一步劃分成一個(gè)完全為被照射小三角面元Tk1和兩個(gè)完全被照射的小三角面元Tk2、Tk3,其中Tk1由頂點(diǎn)Vk1、A1、A2構(gòu)成,Tk2由頂點(diǎn)Vk2、Vk3、A1構(gòu)成,Tk3由頂點(diǎn)Vk3、A1、A2構(gòu)成。
5.根據(jù)權(quán)利要求1或4所述的基于面元散射的合成孔徑聲納面目標(biāo)快速仿真方法,其特征在于,所述的劃分得到的若干完全被照射小三角面元依據(jù)完全被照射三角面元處理方法得到回波數(shù)據(jù)。
6.根據(jù)權(quán)利要求1所述的基于面元散射的合成孔徑聲納面目標(biāo)快速仿真方法,其特征在于,所述的步驟(1)中的合成孔徑聲納系統(tǒng)的仿真參數(shù)主要包括中心頻率、信號(hào)帶寬、脈沖寬度、脈沖重復(fù)周期、采樣頻率、發(fā)射陣孔徑、接收陣孔徑、陣元個(gè)數(shù)、采樣點(diǎn)數(shù)、最小距離、聲速或拖體速度。
7.根據(jù)權(quán)利要求1所述的基于面元散射的合成孔徑聲納面目標(biāo)快速仿真方法,其特征在于,所述的步驟(1)中的構(gòu)造面目標(biāo)場(chǎng)景模型,通過(guò)ANSYS仿真分析創(chuàng)建面目標(biāo)模型,并對(duì)其進(jìn)行三角面元的網(wǎng)格劃分。
8.根據(jù)權(quán)利要求1所述的基于面元散射的合成孔徑聲納面目標(biāo)快速仿真方法,其特征在于,所述的步驟(2)中的基陣位置,依據(jù)方位向上對(duì)應(yīng)的時(shí)間順序和拖體速度確定聲納基陣的位置,其中當(dāng)需要驗(yàn)證運(yùn)動(dòng)補(bǔ)償和自聚焦算法時(shí),拖體的速度可以根據(jù)實(shí)際需要進(jìn)行設(shè)置。
全文摘要
本發(fā)明提出一種基于面元散射的合成孔徑聲納面目標(biāo)快速仿真方法,所述的方法包含如下步驟(1)設(shè)定合成孔徑聲納系統(tǒng)的仿真參數(shù),構(gòu)造面目標(biāo)場(chǎng)景模型;(2)依據(jù)聲納基陣在方位向上的時(shí)間順序,確定基陣位置,確定表征面目標(biāo)的三角面元被照射狀態(tài)(3)對(duì)所有完全被照射的三角面元,用化簡(jiǎn)得到的赫姆霍茲-基爾霍夫積分公式計(jì)算所有被照射三角面元的散射場(chǎng)的頻域信號(hào),再對(duì)得到的所有三角面元的散射場(chǎng)的頻域信號(hào)進(jìn)行疊加和逆傅里葉變換獲得當(dāng)前位置對(duì)應(yīng)的面目標(biāo)回波數(shù)據(jù);(4)按照方位向的順序重復(fù)步驟(2)~(4),從而獲得當(dāng)前基陣位置對(duì)應(yīng)的回波信號(hào)的頻域表示,再對(duì)頻域信號(hào)進(jìn)行逆傅里葉變換獲得面目標(biāo)對(duì)應(yīng)的完整回波數(shù)據(jù)。
文檔編號(hào)G01S7/52GK101762811SQ20101011224
公開日2010年6月30日 申請(qǐng)日期2010年2月8日 優(yōu)先權(quán)日2009年12月8日
發(fā)明者李保利, 劉維, 劉紀(jì)元, 張春華 申請(qǐng)人:中國(guó)科學(xué)院聲學(xué)研究所