，其中，。為第一門限值。
[0045] S104、重復S103,篩選出過第一門限值的另外N2個頻點G =[%，％，…，?/Λ:]。
[0046] S105、兩兩對應比較W與U中各元素，若二者差值的絕對值小于或等于第二門限 值，且S1ZiP的比值大于第三門限值，則判斷為物體回聲，其余值判斷為雜亂聲脈沖；其中， S1表示當前頻點對應的脈沖的各采樣點的能量平均值，P表示低于第一門限值的頻點對應 的脈沖的各采樣點值的平均值。
[0047] 下面對于本發(fā)明實施例進行仿真驗證，驗證條件做如下設定：干擾噪聲比20dB ; 信號噪聲比20dB(StretCh處理前）；發(fā)射的線性調(diào)頻連續(xù)波信號周期30ms，信號帶寬 5kHz。發(fā)射30個周期的調(diào)頻連續(xù)波信號。仿真中物體與雜亂脈沖的加入位置見表1。取 第一門限值σ電平取值為40,第二門限值δ為5,第三門限值μ為6。仿真中頻域脈沖2 為經(jīng)目標反射后的回波信號，即目標信號，其余為雜亂脈沖。取兩段信號分析頻域特性，進 行傅里葉變換的數(shù)據(jù)長度取為1. 5倍的連續(xù)波周期，即45ms。兩個分析段取相鄰的數(shù)據(jù)。
[0048] 圖4和圖5分別給出了本發(fā)明實施例方法處理后所獲得的物體回聲特性。各回聲 頻譜S1ZiP的值均大于6。由圖4和圖5可見，兩個對比周期內(nèi)的回聲接收期間，物體在頻 域上是較接近的。而其余的雜亂脈沖在兩次回聲處理周期中沒有對應出現(xiàn)。觀測表2所獲 得頻域峰值位置，通過比對過門限信號的峰值位置，確定位置誤差絕對值小于門限值的是 物體，其余過門限脈沖為雜亂聲脈沖。因此證明了本發(fā)明的正確性。
[0049] 表1仿真加入的雜亂脈沖與物體在距離門的位置
[0050]

[0051] 表2所取兩段數(shù)據(jù)獲得的信號頻率
[0052]
[0053] 本發(fā)明采用周期線性調(diào)頻連續(xù)波信號，利用雜亂聲脈沖在頻域峰值不確定，而經(jīng) 物體反射的回波脈沖在兩個比較周期期間頻域變化小的特點，通過門限位置的判斷，確定 回聲物體，從而達到去除雜亂脈沖的目的，本發(fā)明去除雜亂脈沖的效果好。
[0054] 需要說明的是，在本文中，術(shù)語"包括"、"包含"或者其任何其他變體意在涵蓋非排 他性的包含，從而使得包括一系列要素的過程、方法、物品或者裝置不僅包括那些要素，而 且還包括沒有明確列出的其他要素，或者是還包括為這種過程、方法、物品或者裝置所固有 的要素。在沒有更多限制的情況下，由語句"包括一個……"限定的要素，并不排除在包括 該要素的過程、方法、物品或者裝置中還存在另外的相同要素。
[0055] 專業(yè)人員還可以進一步意識到，結(jié)合本文中所公開的實施例描述的各示例的單元 及算法步驟，能夠以電子硬件、計算機軟件或者二者的結(jié)合來實現(xiàn)，為了清楚地說明硬件和 軟件的可互換性，在上述說明中已經(jīng)按照功能一般性地描述了各示例的組成及步驟。這些 功能究竟以硬件還是軟件方式來執(zhí)行，取決于技術(shù)方案的特定應用和設計約束條件。專業(yè) 技術(shù)人員可以對每個特定的應用來使用不同方法來實現(xiàn)所描述的功能，但是這種實現(xiàn)不應 認為超出本發(fā)明的范圍。
[0056] 結(jié)合本文中所公開的實施例描述的方法或算法的步驟可以直接用硬件、處理器執(zhí) 行的軟件模塊，或者二者的結(jié)合來實施。軟件模塊可以置于隨機存儲器（RAM)、內(nèi)存、只讀存 儲器（ROM)、電可編程R0M、電可擦除可編程R0M、寄存器、硬盤、可移動磁盤、CD-ROM、或技術(shù) 領域內(nèi)所公知的任意其它形式的存儲介質(zhì)中。
[0057] 對所公開的實施例的上述說明，使本領域?qū)I(yè)技術(shù)人員能夠?qū)崿F(xiàn)或使用本發(fā)明。 對這些實施例的多種修改對本領域的專業(yè)技術(shù)人員來說將是顯而易見的，本文中所定義的 一般原理可以在不脫離本發(fā)明的精神或范圍的情況下，在其它實施例中實現(xiàn)。因此，本發(fā)明 將不會被限制于本文所示的這些實施例，而是要符合與本文所公開的原理和新穎特點相一 致的最寬的范圍。
【主權(quán)項】
1. 一種適用于聲吶一維距離像穩(wěn)定特征提取的雜亂脈沖去除方法，其特征在于，包 括： 5101、 聲吶發(fā)射周期線性調(diào)頻連續(xù)波信號，并接收回聲信號，其中，所述回聲信號包括 雜亂脈沖和所述線性調(diào)頻脈沖信號經(jīng)物體反射后形成的回聲脈沖； 5102、 對所述回聲信號進行Stretch處理，得到頻域回聲信號； 5103、 對所述頻域回聲信號進行門限檢測，篩選出過第一門限值的&個頻點 5104、 重復S103,篩選出過第一門限值的另外隊個頻點U 5105、 兩兩對應比較W與U中各元素，若二者差值的絕對值小于或等于第二門限值，且 Si/P的比值大于第三門限值，則判斷為物體回聲，其余值判斷為雜亂聲脈沖；其中，Si表示 當前頻點對應的脈沖的各采樣點的能量平均值，P表示低于第一門限值的頻點對應的脈沖 的各采樣點值的平均值。2. 如權(quán)利要求1所述的適用于聲吶一維距離像穩(wěn)定特征提取的雜亂脈沖去除方法，其 特征在于，所述聲吶發(fā)射周期線性調(diào)頻連續(xù)波信號，并接收回聲信號，具體包括： 聲吶發(fā)射周期線性調(diào)頻連續(xù)波信號，其中，所述線性調(diào)頻脈沖信號為式中，t為時間變量，f。為載頻，|為調(diào)頻斜率，T為調(diào)頻連續(xù) 波周期，rect[?]表示周期方波的一個周期； 聲吶接收多個物體的回聲信號，其中，所述多個物體回聲信號為中，ai為第i個回聲脈沖的幅度復系數(shù)，ti為第i個回聲 脈沖的時延，n(t)為聲吶的接收通道噪聲，M為回聲物體的個數(shù)。3. 如權(quán)利要求2所述的適用于聲吶一維距離像穩(wěn)定特征提取的雜亂脈沖去除方法，其 特征在于，對所述回聲信號進行Stretch處理，得到頻域回聲信號，具體包括： 令周期連續(xù)波變換函數(shù)，對所述回聲信號進行Stretch處理， 得到信號z(t)，即式中0〈t〈N?T，N表示發(fā)射調(diào)頻連續(xù)波的所有周期個數(shù)； 選取所述信號z(t)其中的一段進行傅里葉變換，得到頻域回聲信號，其中，所述頻域 回聲信號為a為整數(shù)或分數(shù)，即選取的傅里葉變換的數(shù)據(jù)長度是 周期T的整數(shù)倍或者是分數(shù)倍。4.如權(quán)利要求3所述的適用于聲吶一維距離像穩(wěn)定特征提取的雜亂脈沖去除方 法，其特征在于，對所述頻域回聲信號進行門限檢測，篩選出過第一門限值的&個頻點 ,=[>%，叫，…，^ ]，具體包括： 求取所述頻域回聲信號的包絡.式中，yi(?)與yQ(?)為頻 域回波信號y(?)的實部與虛部； 篩選出滿足A〇) > 〇時《的&個取值，作為過第一門限值的Ni個頻點Fu'、]，其中，。為第一門限值。
【專利摘要】本發(fā)明公開了一種適用于聲吶一維距離像穩(wěn)定特征提取的雜亂脈沖去除方法，包括：S101、聲吶發(fā)射周期線性調(diào)頻連續(xù)波信號，并接收回聲信號；S102、對所述回波信號進行Stretch處理，得到頻域回波信號；S103、對所述頻域波信號進行門限檢測，篩選出過第一門限值的N1個頻點；S104、重復S103，篩選出過第一門限值的另外N2個頻點；S105、兩兩對應比較W與U中各元素，若二者差值的絕對值小于或等于第二門限值，且Si/ρ的比值大于第三門限值，則判斷為物體回聲，其余值判斷為雜亂脈沖。本發(fā)明可以從聲吶接收的信號中去除雜亂脈沖。
【IPC分類】G01S7/52, G01S7/527
【公開號】CN105093209
【申請?zhí)枴緾N201510502805
【發(fā)明人】王峰, 蔣德富
【申請人】河海大學
【公開日】2015年11月25日
【申請日】2015年8月14日