一種聲吶應(yīng)答式干擾的剔除方法與流程

本發(fā)明屬于浮空平臺(tái)外輻射源雷達(dá)
技術(shù)領(lǐng)域：
，具體涉及一種聲吶應(yīng)答式干擾的剔除方法。
背景技術(shù)：
：在水下聲對(duì)抗中，如圖1所示，應(yīng)答式干擾由于與目標(biāo)回聲相似，常常被誤認(rèn)為是目標(biāo)，造成檢測虛警，大量的虛警使得目標(biāo)的檢測與航跡均無法完成，因此，該類干擾是聲吶面臨的額難題。且從圖1中可以看到，對(duì)于目標(biāo)而言，對(duì)于任意脈沖，距離沒有發(fā)生改變(改變很小，沒有反應(yīng)出來)，而應(yīng)答式干擾由于與聲吶發(fā)射信號(hào)存在異步效應(yīng)，因而每次回聲均不在同一時(shí)延，故常規(guī)的自適應(yīng)波束形成等方法對(duì)該類干擾難以起效，尤其是當(dāng)該類干擾從主瓣進(jìn)入聲吶系統(tǒng)。技術(shù)實(shí)現(xiàn)要素：針對(duì)上述問題，本發(fā)明提出一種聲吶應(yīng)答式干擾的剔除方法，采用線性調(diào)頻信號(hào)(LinearFrequencyModulation，LFM)發(fā)射，利用分?jǐn)?shù)階傅里葉變換技術(shù)(FractionFourierTransform，F(xiàn)RFT)，并結(jié)合自回歸(AutoRegressive，AR)模型超分辨譜估計(jì)技術(shù)，實(shí)現(xiàn)了多發(fā)射周期下的異步應(yīng)答式干擾的去除。實(shí)現(xiàn)上述技術(shù)目的，達(dá)到上述技術(shù)效果，本發(fā)明通過以下技術(shù)方案實(shí)現(xiàn)：一種聲吶應(yīng)答式干擾的剔除方法，包括以下步驟：步驟1、聲吶發(fā)射出多脈沖線性調(diào)頻信號(hào)，并接收包括目標(biāo)回聲信號(hào)和應(yīng)答式干擾信號(hào)的多脈沖回波信號(hào)；步驟2、針對(duì)其中的單個(gè)脈沖數(shù)據(jù)接收周期，根據(jù)系統(tǒng)設(shè)定的接收時(shí)間接收回波信號(hào)，對(duì)該接收到的回波信號(hào)采用分?jǐn)?shù)階傅里葉變換的方法進(jìn)行變換；步驟3、采用自回歸模型譜估計(jì)方法對(duì)分?jǐn)?shù)階傅里葉變換所得的輸出信號(hào)進(jìn)行頻率估計(jì)，得到譜估計(jì)結(jié)果；步驟4、根據(jù)譜估計(jì)結(jié)果，采用序列統(tǒng)計(jì)恒虛警算法獲得距離單元結(jié)果，實(shí)現(xiàn)目標(biāo)所處距離單元判斷與信號(hào)檢測；步驟5、針對(duì)所有脈沖的數(shù)據(jù)接收周期內(nèi)接收的回波信號(hào)均采用步驟2-4的方法進(jìn)行處理，然后設(shè)定門限閾值ε，對(duì)各個(gè)脈沖的檢測所得信號(hào)所處的距離單元結(jié)果進(jìn)行脈沖間比對(duì)，當(dāng)小于設(shè)定的門限閾值ε，則認(rèn)為檢測所得信號(hào)為目標(biāo)；當(dāng)大于設(shè)定的門限閾值ε，則舍棄，從而完成聲吶應(yīng)答式干擾的剔除。所述步驟1中，線性調(diào)頻信號(hào)為：s(t)=ej2π(f0-0.5μt)t---(1)]]>式中：t為時(shí)間變量，f0為載波頻率，μ為調(diào)頻斜率，ej[·]為復(fù)數(shù)；回波信號(hào)為：y(t)=Σi=0N-1βis(t-τi)ej2πfi(t-τi)+n(t)---(2)]]>式中：τi表示第i路時(shí)延，fi是第i個(gè)回波信號(hào)的多普勒頻移，βi是第i個(gè)回波信號(hào)的幅度系數(shù)；n(t)是接收機(jī)噪聲，服從高斯分布；N代表檢測到的回波信號(hào)個(gè)數(shù)；考慮到線性調(diào)頻信號(hào)的多普勒頻移不變性，因此僅關(guān)注時(shí)延方向的信號(hào)變化，則可在回波信號(hào)的表示式中省去多普勒頻移，得到：y(t)=Σi=0N-1βis(t-τi)+n(t)---(3)]]>將式(1)帶入式(3)得到：y(t)=Σi=0N-1βiej[2π(f0+μτi)t-πμt2+2πτi2-2πf0τi]+n(t)---(4)]]>所述步驟2中，對(duì)該接收到的回波信號(hào)采用分?jǐn)?shù)階傅里葉變換的方法進(jìn)行變換，具體為：分?jǐn)?shù)階傅里葉變換的核函數(shù)為:η(t,ω)=1-jcotα2πexp{jt2+ω22cotα-jωtcscα)α≠nπδ(t-ω)α=2nπδ(t+ω)α=(2n±1)π---(5)]]>式中：ω為角頻率，單位為弧度/秒；t為時(shí)間，δ(t)表示沖激響應(yīng)函數(shù)，α表示旋轉(zhuǎn)因子，cotα表示余切函數(shù)，cscα表示余割函數(shù)，exp{j·}表示復(fù)數(shù)；回波信號(hào)y(t)的分?jǐn)?shù)階傅里葉變換可表示為：ξ(ω)=∫t=0t=∞y(t)η(t,ω)dt---(6)]]>將式(5)代入式(6)可得到：ξ(ω)=1-jcotα2πexp{jω22cotα}·∫t=0t=∞y(t)exp{jcotα2t2}exp{-j·cscα·ωt}dt---(7)]]>采用分?jǐn)?shù)階傅里葉變換對(duì)信號(hào)進(jìn)行旋轉(zhuǎn)搜索處理，得到經(jīng)過分?jǐn)?shù)階傅里葉變換后的信號(hào)：u(t)=y(t)exp{jcotα2t2}---(8)]]>搜索α即可實(shí)現(xiàn)分?jǐn)?shù)階傅里葉變換的信號(hào)搜索處理，得到要檢測的信號(hào)；對(duì)比式(1)和式(7)t2項(xiàng)的系數(shù)，如果下式成立：cotα2=πμ---(9)]]>則式(9)將完成式(1)中線性調(diào)頻信號(hào)的調(diào)頻斜率μ的搜索，完成信號(hào)的匹配，用于恒虛警檢測，u(t)變?yōu)槎鄠€(gè)單頻信號(hào)的組合，如下所示：式中，為變換后的噪聲項(xiàng)，第一項(xiàng)為多個(gè)單頻信號(hào)的組合，用于后續(xù)采用自回歸模型譜估計(jì)方法進(jìn)行頻率估計(jì)，該多頻信號(hào)即對(duì)不同時(shí)延的目標(biāo)信號(hào)與應(yīng)答干擾信號(hào)。所述步驟3中，采用自回歸模型譜估計(jì)方法對(duì)經(jīng)過分?jǐn)?shù)階傅里葉變換的信號(hào)u(t)進(jìn)行頻率估計(jì)，將譜估計(jì)結(jié)果輸出作為下一步檢測所用數(shù)據(jù)，具體為：定義n為離散時(shí)間，自回歸模型譜估計(jì)方法進(jìn)行功率譜估計(jì)采用白噪聲序列ζ(n)激勵(lì)自回歸模型系統(tǒng)H(z)，從而得到信號(hào)u(n)，由u(n)的自相關(guān)函數(shù)來估計(jì)H(z)的參數(shù)，從而估計(jì)u(n)的功率譜，u(n)為式(10)中u(t)的離散時(shí)間表示，當(dāng)采用方差為σ2的白噪聲序列ζ(n)，u(n)經(jīng)過一個(gè)L階的自回歸模型系統(tǒng)：H(z)=1A(z)=11+Σk=1Lakz-k---(11)]]>式中ak為A(z)的時(shí)域系數(shù)，L表示模型階數(shù)，可獲得如下的信號(hào)：u(n)=Σk=1Laku(n-k)+ζ(n)---(12)]]>則u(n)的功率譜P(ejω)為：P(ejω)=σ2|1+Σk=1Lake-jωk|---(13)]]>當(dāng)確定了自回歸模型的參數(shù)ak，即可實(shí)現(xiàn)對(duì)回波信號(hào)頻率的估計(jì)，實(shí)現(xiàn)對(duì)不同時(shí)延的信號(hào)與應(yīng)答式干擾的估計(jì)，自回歸模型的正則方程即Yule-Walker方程具有如下形式：式中的r(m)為自相關(guān)函數(shù)，具有如下形式：r(m)=Σn=0Ns-1u(n)u*(n-m)---(15)]]>式中Ns是計(jì)算所需采樣點(diǎn)數(shù)，自回歸模型的參數(shù)ak通過如下表達(dá)式計(jì)算得到：am(m)=-[r(m)+Σk=1m-1am-1(k)r(m-k)]/σm-12---(16)]]>am(k)＝am-1(k)+am(m)am-1(m-k)(17)式(16)、(17)中k＝1,2,…,m-1，括號(hào)中的序號(hào)為循環(huán)遞推計(jì)算次數(shù)；σk2=(1-|am(m)|2)σk-12,σ02=r(0)---(18)]]>先計(jì)算階次m＝1時(shí)的預(yù)測系數(shù)am(k)＝a1(1)和直到完成計(jì)算m＝L時(shí)的預(yù)測系數(shù)aL(1),aL(2),…,aL(L)與計(jì)算完成后除去括號(hào)中的遞推次數(shù)，將系數(shù)a1(1)，a2(2)…,aL(L)記為a1，a2…,aL，代入式(13)即可完成自回歸模型譜估計(jì)的功率譜計(jì)算。所述步驟4中，采用序列統(tǒng)計(jì)恒虛警算法(OS-CFAR)實(shí)現(xiàn)目標(biāo)距離單元判斷，具體為：針對(duì)自回歸功率譜的估計(jì)結(jié)果P(ejω)，記其模值為Ql，采用l表示序列的序號(hào)，采用序列檢測恒虛警實(shí)現(xiàn)目標(biāo)檢測，其實(shí)現(xiàn)過程如下：1)對(duì)自回歸模型進(jìn)行譜估計(jì)，針對(duì)包絡(luò)輸出進(jìn)行檢測；2)獲取恒虛警門限：對(duì)于其中的第l點(diǎn)，其幅度值為Ql，利用其兩側(cè)的幅度值數(shù)據(jù)進(jìn)行門限選取，具體為：l點(diǎn)左側(cè)幅度值Ql-M～Ql-1由大到小排列，M是左側(cè)滑窗長度，M為偶數(shù)；Y1(k)＝Max-Min{Ql-M,…,Ql-1}(19)這里Max-Min{·}表示從大到小排列，選擇第K1個(gè)幅度值作為門限，記為Y1(K1)，K1的取值取M/2。l點(diǎn)右側(cè)幅度值Ql+1～Ql+M由從大到小排列，M是左側(cè)滑窗長度，取M為偶數(shù)，的序列如下Y2(K)＝Max-Min{Ql+1,…，Ql+M}(20)選擇第K2個(gè)幅度值作為門限，記為Y2(K2)，K2的取值為M/2；3)選擇Y1(K1)和Y2(K2)中較小的一個(gè)作為門限：u0＝Min{Y1(K1),Y2(K2)}(21)Min{·}表示求取最小值；4)當(dāng)滿足下式時(shí)：Ql≥u0(22)則單個(gè)脈沖數(shù)據(jù)接收周期內(nèi)接收的回波信號(hào)的檢測序列為：d(n)＝[0,0,0,…,1,1,0,0,…,0,0]T(23)式中[·]T表示為共軛轉(zhuǎn)置，標(biāo)記為1的點(diǎn)為檢測出目標(biāo)的位置，未檢測出目標(biāo)的位置標(biāo)記為0。所述各個(gè)脈沖的距離單元結(jié)果進(jìn)行脈沖間比對(duì)，具體比對(duì)方法為：設(shè)定閾值為ε，當(dāng)3個(gè)脈沖檢測結(jié)果中，如果為1的位置位差小于ε，則認(rèn)為為目標(biāo)，否則認(rèn)為為干擾。本發(fā)明的有益效果：本發(fā)明通過發(fā)射多脈沖線性調(diào)頻信號(hào)(LFM)，采用分?jǐn)?shù)階傅里葉變換(FrFT)的方法，將接收數(shù)據(jù)中的線性調(diào)頻脈沖信號(hào)進(jìn)行旋轉(zhuǎn)，將時(shí)延信號(hào)轉(zhuǎn)變?yōu)轭l率，并結(jié)合自回歸(AR)模型超分辨譜估計(jì)技術(shù)和序列統(tǒng)計(jì)恒虛警算法，實(shí)現(xiàn)了多發(fā)射周期下的異步應(yīng)答式干擾的去除。同時(shí)，相對(duì)于傳統(tǒng)的脈沖壓縮方法，本發(fā)明提高了目標(biāo)在時(shí)延維度的分辨率。附圖說明圖1為脈沖串工作時(shí)應(yīng)答式干擾的工作時(shí)序與目標(biāo)的差異示意圖。圖2為本發(fā)明的聲吶應(yīng)答式干擾的剔除方法流程示意圖。圖3(a)和(b)為回波脈沖的AR譜的幅度示意圖。圖4(a)和4(b)為采用本發(fā)明的剔除方法的結(jié)果示意圖。圖5為采用本發(fā)明的剔除方法對(duì)回波脈沖串的檢測結(jié)果示意圖。具體實(shí)施方式為了使本發(fā)明的目的、技術(shù)方案及優(yōu)點(diǎn)更加清楚明白，以下結(jié)合實(shí)施例，對(duì)本發(fā)明進(jìn)行進(jìn)一步詳細(xì)說明。應(yīng)當(dāng)理解，此處所描述的具體實(shí)施例僅僅用以解釋本發(fā)明，并不用于限定本發(fā)明。下面結(jié)合附圖對(duì)本發(fā)明的應(yīng)用原理作詳細(xì)的描述。如圖2所示，一種聲吶應(yīng)答式干擾的剔除方法，包括以下步驟：步驟1、聲吶發(fā)射出多脈沖線性調(diào)頻信號(hào)，并接收包括目標(biāo)回聲信號(hào)和應(yīng)答式干擾信號(hào)的多脈沖回波信號(hào)；步驟2、針對(duì)其中的單個(gè)脈沖數(shù)據(jù)接收周期，根據(jù)系統(tǒng)設(shè)定的接收時(shí)間接收回波信號(hào)，對(duì)該接收到的回波信號(hào)采用分?jǐn)?shù)階傅里葉變換的方法進(jìn)行變換；步驟3、采用自回歸模型譜估計(jì)方法對(duì)分?jǐn)?shù)階傅里葉變換所得的輸出信號(hào)進(jìn)行頻率估計(jì)，得到譜估計(jì)結(jié)果；步驟4、根據(jù)譜估計(jì)結(jié)果，采用序列統(tǒng)計(jì)恒虛警算法獲得距離單元結(jié)果，實(shí)現(xiàn)目標(biāo)所處距離單元判斷與信號(hào)檢測；步驟5、針對(duì)所有脈沖的數(shù)據(jù)接收周期內(nèi)接收的回波信號(hào)均采用步驟2-4的方法進(jìn)行處理，然后設(shè)定門限閾值ε，對(duì)各個(gè)脈沖的檢測所得信號(hào)所處的距離單元結(jié)果進(jìn)行脈沖間比對(duì)，當(dāng)小于設(shè)定的門限閾值ε，則認(rèn)為檢測所得信號(hào)為目標(biāo)；當(dāng)大于設(shè)定的門限閾值ε，則舍棄，從而完成聲吶應(yīng)答式干擾的剔除。所述步驟1中，線性調(diào)頻信號(hào)為：s(t)=ej2π(f0-0.5μt)t---(1)]]>式中：t為時(shí)間變量，f0為載波頻率，μ為調(diào)頻斜率，ej[·]為復(fù)數(shù)；回波信號(hào)為：y(t)=Σi=0N-1βis(t-τi)ej2πfi(t-τi)+n(t)---(2)]]>式中：τi表示第i路時(shí)延，fi是第i個(gè)回波信號(hào)的多普勒頻移，βi是第i個(gè)回波信號(hào)的幅度系數(shù)；n(t)是接收機(jī)噪聲，服從高斯分布；N代表檢測到的回波信號(hào)個(gè)數(shù)；考慮到線性調(diào)頻信號(hào)的多普勒頻移不變性，因此僅關(guān)注時(shí)延方向的信號(hào)變化，則可在回波信號(hào)的表示式中省去多普勒頻移，得到：y(t)=Σi=0N-1βis(t-τi)+n(t)---(3)]]>將式(1)帶入式(3)得到：y(t)=Σi=0N-1βiej[2π(f0+μτi)t-πμt2+2πτi2-2πf0τi]+n(t)---(4)]]>所述步驟2中，對(duì)該接收到的回波信號(hào)采用分?jǐn)?shù)階傅里葉變換的方法進(jìn)行變換，具體為：分?jǐn)?shù)階傅里葉變換的核函數(shù)為:η(t,ω)=1-jcotα2πexp{jt2+ω22cotα-jωtcscα)α≠nπδ(t-ω)α=2nπδ(t+ω)α=(2n±1)π---(5)]]>式中：ω為角頻率，單位為弧度/秒；t為時(shí)間，δ(t)表示沖激響應(yīng)函數(shù)，α表示旋轉(zhuǎn)因子，cotα表示余切函數(shù)，cscα表示余割函數(shù)，exp{j·}表示復(fù)數(shù)；回波信號(hào)y(t)的分?jǐn)?shù)階傅里葉變換可表示為：ξ(ω)=∫t=0t=∞y(t)η(t,ω)dt---(6)]]>將式(5)代入式(6)可得到：ξ(ω)=1-jcotα2πexp{jω22cotα}·∫t=0t=∞y(t)exp{jcotα2t2}exp{-j·cscα·ωt}dt---(7)]]>采用分?jǐn)?shù)階傅里葉變換對(duì)信號(hào)進(jìn)行旋轉(zhuǎn)搜索處理，得到經(jīng)過分?jǐn)?shù)階傅里葉變換后的信號(hào)：u(t)=y(t)exp{jcotα2t2}---(8)]]>搜索α即可實(shí)現(xiàn)分?jǐn)?shù)階傅里葉變換的信號(hào)搜索處理，得到要檢測的信號(hào)；對(duì)比式(1)和式(7)t2項(xiàng)的系數(shù)，如果下式成立：cotα2=πμ---(9)]]>則式(9)將完成式(1)中線性調(diào)頻信號(hào)的調(diào)頻斜率μ的搜索，完成信號(hào)的匹配，用于恒虛警檢測，u(t)變?yōu)槎鄠€(gè)單頻信號(hào)的組合，如下所示：式中，為變換后的噪聲項(xiàng)，第一項(xiàng)為多個(gè)單頻信號(hào)的組合，用于后續(xù)采用自回歸模型譜估計(jì)方法進(jìn)行頻率估計(jì)，該多頻信號(hào)即對(duì)不同時(shí)延的目標(biāo)信號(hào)與應(yīng)答干擾信號(hào)。所述步驟3中，采用自回歸模型譜估計(jì)方法對(duì)經(jīng)過分?jǐn)?shù)階傅里葉變換的信號(hào)u(t)進(jìn)行頻率估計(jì)，將譜估計(jì)結(jié)果輸出作為下一步檢測所用數(shù)據(jù)，具體為：定義n為離散時(shí)間，自回歸模型譜估計(jì)方法進(jìn)行功率譜估計(jì)采用白噪聲序列ζ(n)激勵(lì)自回歸模型系統(tǒng)H(z)，從而得到信號(hào)u(n)，由u(n)的自相關(guān)函數(shù)來估計(jì)H(z)的參數(shù)，從而估計(jì)u(n)的功率譜，u(n)為式(10)中u(t)的離散時(shí)間表示，當(dāng)采用方差為σ2的白噪聲序列ζ(n)，u(n)經(jīng)過一個(gè)L階的自回歸模型系統(tǒng)：H(z)=1A(z)=11+Σk=1Lakz-k---(11)]]>式中ak為A(z)的時(shí)域系數(shù)，L表示模型階數(shù)，可獲得如下的信號(hào)：u(n)=Σk=1Laku(n-k)+ζ(n)---(12)]]>則u(n)的功率譜P(ejω)為：P(ejω)=σ2|1+Σk=1Lake-jωk|---(13)]]>當(dāng)確定了自回歸模型的參數(shù)ak，即可實(shí)現(xiàn)對(duì)回波信號(hào)頻率的估計(jì)，實(shí)現(xiàn)對(duì)不同時(shí)延的信號(hào)與應(yīng)答式干擾的估計(jì)，自回歸模型的正則方程即Yule-Walker方程具有如下形式：式中的r(m)為自相關(guān)函數(shù)，具有如下形式：r(m)=Σn=0Ns-1u(n)u*(n-m)---(15)]]>式中Ns是計(jì)算所需采樣點(diǎn)數(shù)，自回歸模型的參數(shù)ak通過如下表達(dá)式計(jì)算得到：am(m)=-[r(m)+Σk=1m-1am-1(k)r(m-k)]/σm-12---(16)]]>am(k)＝am-1(k)+am(m)am-1(m-k)(17)式(16)、(17)中k＝1,2,…,m-1，括號(hào)中的序號(hào)為循環(huán)遞推計(jì)算次數(shù)；σk2=(1-|am(m)|2)σk-12,σ02=r(0)---(18)]]>先計(jì)算階次m＝1時(shí)的預(yù)測系數(shù)am(k)＝a1(1)和直到完成計(jì)算m＝L時(shí)的預(yù)測系數(shù)aL(1),aL(2),…,aL(L)與計(jì)算完成后除去括號(hào)中的遞推次數(shù)，將系數(shù)a1(1)，a2(2)…,aL(L)記為a1，a2…,aL，代入式(13)即可完成自回歸模型譜估計(jì)的功率譜計(jì)算。所述步驟4中，采用序列統(tǒng)計(jì)恒虛警算法(OS-CFAR)實(shí)現(xiàn)目標(biāo)距離單元判斷，具體為：針對(duì)自回歸功率譜的估計(jì)結(jié)果P(ejω)，記其模值為Ql，采用l表示序列的序號(hào)，采用序列檢測恒虛警實(shí)現(xiàn)目標(biāo)檢測，其實(shí)現(xiàn)過程如下：1)對(duì)自回歸模型進(jìn)行譜估計(jì)，針對(duì)包絡(luò)輸出進(jìn)行檢測；2)獲取恒虛警門限：對(duì)于其中的第l點(diǎn)，其幅度值為Ql，利用其兩側(cè)的幅度值數(shù)據(jù)進(jìn)行門限選取，具體為：l點(diǎn)左側(cè)幅度值Ql-M～Ql-1由大到小排列，M是左側(cè)滑窗長度，M為偶數(shù)；Y1(k)＝Max-Min{Ql-M,…,Ql-1}(19)這里Max-Min{·}表示從大到小排列，選擇第K1個(gè)幅度值作為門限，記為Y1(K1)，K1的取值取M/2。l點(diǎn)右側(cè)幅度值Ql+1～Ql+M由從大到小排列，M是左側(cè)滑窗長度，取M為偶數(shù)，的序列如下Y2(K)＝Max-Min{Ql+1,…，Ql+M}(20)選擇第K2個(gè)幅度值作為門限，記為Y2(K2)，K2的取值為M/2；3)選擇Y1(K1)和Y2(K2)中較小的一個(gè)作為門限：u0＝Min{Y1(K1),Y2(K2)}(21)Min{·}表示求取最小值；4)當(dāng)滿足下式時(shí)：Ql≥u0(22)則單個(gè)脈沖數(shù)據(jù)接收周期內(nèi)接收的回波信號(hào)的檢測序列為：d(n)＝[0,0,0,…,1,1,0,0,…,0,0]T(23)式中[·]T表示為共軛轉(zhuǎn)置，標(biāo)記為1的點(diǎn)為檢測出目標(biāo)的位置，未檢測出目標(biāo)的位置標(biāo)記為0。所述各個(gè)脈沖的距離單元結(jié)果進(jìn)行脈沖間比對(duì)，具體比對(duì)方法為：設(shè)定閾值為ε，當(dāng)3個(gè)脈沖檢測結(jié)果中，如果為1的位置位差小于ε，則認(rèn)為為目標(biāo)，否則認(rèn)為為干擾。實(shí)施例一采用計(jì)算機(jī)仿真驗(yàn)證本發(fā)明方法的性能。在本實(shí)施例中，采用3個(gè)發(fā)射脈沖，脈沖寬度為100ms，信號(hào)帶寬100Hz，發(fā)射間隔為1s。信號(hào)載波f0＝500HZ，μ＝107，F(xiàn)rFT旋轉(zhuǎn)因子α的搜索范圍為106-108，搜索步長為105，對(duì)于AR算法μ0取0.001，取ε為2，即2個(gè)采樣點(diǎn)，采樣頻率為2kHz，脈壓前信噪比為0dB，真目標(biāo)所處距離門為60，其它位置干擾隨機(jī)出現(xiàn)。首先采用第一個(gè)發(fā)射脈沖接收兩個(gè)采樣點(diǎn)期間觀察本發(fā)明的效果。圖3(a)和圖3(b)為脈沖壓縮方法的處理結(jié)果。圖4(a)和圖4(b)為本發(fā)明的方法的處理結(jié)果。本發(fā)明檢測到了包括干擾的6個(gè)目標(biāo)，而傳統(tǒng)脈沖壓縮方法檢測到4個(gè)目標(biāo)，由此可見，傳統(tǒng)的脈沖壓縮方法分辨率低于本發(fā)明方法。采用本發(fā)明方法實(shí)現(xiàn)3脈沖檢測，檢測結(jié)果見圖5。由于目標(biāo)的距離不變，因此可以通過距離對(duì)比實(shí)現(xiàn)信號(hào)檢測，見圖3的沿著箭頭所指方向的檢測目標(biāo)位置相同，距離位置誤差小于2。因此，該檢測結(jié)果即為目標(biāo)。以上顯示和描述了本發(fā)明的基本原理和主要特征和本發(fā)明的優(yōu)點(diǎn)。本行業(yè)的技術(shù)人員應(yīng)該了解，本發(fā)明不受上述實(shí)施例的限制，上述實(shí)施例和說明書中描述的只是說明本發(fā)明的原理，在不脫離本發(fā)明精神和范圍的前提下，本發(fā)明還會(huì)有各種變化和改進(jìn)，這些變化和改進(jìn)都落入要求保護(hù)的本發(fā)明范圍內(nèi)。本發(fā)明要求保護(hù)范圍由所附的權(quán)利要求書及其等效物界定。當(dāng)前第1頁1 2 3