m個信號采樣點，找出（Vi，a」)參數(shù)域中最大值對應(yīng)的參數(shù)，
[0042]
[0043] 令目標(biāo)的徑向速度估計\ =4，令目標(biāo)的徑向速度估計< =4 ;如果A vl> Δ心或 Δ ai> Δ af，則跳回第2)步，否則進入下一步；
[0044] 7)將參數(shù)(?)對應(yīng)的頻域分布V i, j(n，1)作為該信號采樣點的多普勒分布；
[0045] 本發(fā)明的有益效果是：
[0046] 對比現(xiàn)有技術(shù)，本技術(shù)方案所述的多項式拉東-多項式傅里葉變換的高超聲速目 標(biāo)檢測方法，有益效果在于：
[0047] 1)可以解決高超聲速隱身機動目標(biāo)帶來的距離走動和多普勒走動使得相參積累 能量不能聚焦的問題，實現(xiàn)了低信噪比環(huán)境下臨近空間高超聲速機動目標(biāo)的相參積累，提 尚了雷達的臨近空間尚超聲速目標(biāo)發(fā)現(xiàn)能力。
[0048] 2)該方法除了能估計目標(biāo)的徑向速度外，還能估計目標(biāo)的徑向加速度；
[0049] 3)該方法具有同時檢測臨近空間高超聲速目標(biāo)和傳統(tǒng)目標(biāo)的多目標(biāo)檢測能力；
[0050] 4)該方法中的多分辨率搜索方法為計算復(fù)雜度和檢測概率之間折中選擇提供了 可能，相當(dāng)于提供了一種提高算法實時性的選擇方法；
[0051] 5)該方法能依據(jù)相參積累時間、目標(biāo)的最大可能徑向速度、最大可能徑向加速度、 雷達相關(guān)參數(shù)自適應(yīng)的確定參數(shù)空間搜索的維度。
【附圖說明】
[0052]附圖1本發(fā)明的方法步驟流程圖；
[0053]附圖2是利用MTD方法相參積累后的距離-多普勒分布 [0054]附圖3是利用RFT方法相參積累后的距離-多普勒分布；
[0055] 附圖4是利用PRPFT方法相參積累后的距離-多普勒分布；
【具體實施方式】
[0056] 本發(fā)明的高超聲速目標(biāo)檢測方法，首先對待積累的Ν個周期的信號分別進行采樣， 提取慢時間-快時間二維平面中的目標(biāo)觀測值，然后確定多項式拉東-多項式傅里葉變換相 參積累檢測的初始化參數(shù)，利用多項式拉東-多項式傅里葉變換，在參數(shù)空間進行搜索補償 積累，遍歷所有信號采樣點，得到相參積累后的距離-多普勒分布圖，最后對距離-多普勒分 布圖進行恒虛警檢測，判斷目標(biāo)有無從而實現(xiàn)目標(biāo)的檢測。
[0057] 下面結(jié)合附圖，詳細(xì)描述本發(fā)明的技術(shù)方案，參照附圖2,本發(fā)明的具體步驟包括：
[0058] 步驟(一）、對待積累的Ν個周期的信號分別進行采樣，對采樣數(shù)據(jù)進行離散化處 理，提取慢時間-快時間二維平面中的目標(biāo)觀測值，然后對Ν個脈沖重復(fù)周期內(nèi)的采樣信號 分別進行脈沖壓縮，得到脈沖壓縮后的二維信號矩陣so(n，m)，其中η代表回波信號的個數(shù) 標(biāo)號，11 = 1，2，...少4為回波信號的總數(shù)，!11代表信號采樣點數(shù)的標(biāo)號，111=1，2，~，，為信 號采樣點總數(shù)；
[0059] 步驟（二）、多項式拉東-多項式傅里葉變換（Polynomial Radon-Polynomial Fourier Transform，PRPFT)高超聲速目標(biāo)相參積累檢測的參數(shù)初始化:根據(jù)雷達搜索范圍 和角速度，確定波束駐留時間，根據(jù)雷達脈沖重復(fù)頻率和波束駐留時間，確定相參積累的時 間和脈沖數(shù)，根據(jù)相參積累時間確定參數(shù)空間搜索的維數(shù)、搜索的分辨率、搜索步進的大 小;具體的
[0060] 1)波束駐留時間??Θβ/ο^，其中θβ代表波束寬度，ωΑ表雷達轉(zhuǎn)速；
[0061] 2)相參積累的時間T = N/fp，脈沖數(shù)N=Ti · fp，其中fp表示脈沖重復(fù)頻率；
[0062] 3)參數(shù)空間搜索的維數(shù)確定方法是：
[0063] 假設(shè)積累時間為T，高超聲速目標(biāo)的最大速度Vmax，最大加速度為amax，最大加加速 度為，光速為C，雷達信號帶寬為B、波長為λ，忽略目標(biāo)二階以上的運動對回波信號的影 響，如果對應(yīng)參數(shù)滿足
[0064]
[0065]則搜索參數(shù)維數(shù)為0，此時不用進行參數(shù)空間搜索，PRPFT退化為傳統(tǒng)的MTD方法； 如果上式不滿足，再判斷對應(yīng)參數(shù)是否滿足
[0066]
[0067]如果上式滿足，則搜索參數(shù)維數(shù)為1，即只對目標(biāo)速度參數(shù)進行搜索，此時PRPFT退 化為RFT(Radon-Fourier Transform)方法;如果上式不滿足，再判斷對應(yīng)參數(shù)是否滿足
[0068]
[0069]如果上式滿足，則搜索參數(shù)維數(shù)為2,此時利用PRPFT方法同時對目標(biāo)徑向速度和 徑向加速度參數(shù)進行搜索；
[0070] 4)搜索的分辨率nv取值范圍為[0，NV];
[0071] 5)搜索步進的大小確定方法是：令徑向速度搜索步進％ ，其中， 現(xiàn)+ 1，Δ Δ r/T，Δ r為距離采樣的間隔，r〇und( ·)表示取括號內(nèi)實數(shù)最近 的整數(shù)，I表示循環(huán)次數(shù);徑向加速度搜索步進~ = ~·其中，+ A r 為距離采樣的間隔，其中，l = maxhK、！，l表示目標(biāo)徑向飛行最大可能徑向加速度， 表示目標(biāo)切向飛行帶來的最大可能徑向加速度，'=i，maX{ ·，· }表示取大括號 內(nèi)兩個實數(shù)中較大的一個。
[0072] 步驟(三）、以第1個脈沖的每個信號采樣點為起點，利用多項式拉東-多項式傅里 葉變換，在參數(shù)空間進行搜索補償積累，得到每個信號采樣點對應(yīng)的參數(shù)空間檢測單元圖， 找出該檢測單元圖中信號幅度的最大值，并找出該最大值對應(yīng)的參數(shù)，將該參數(shù)對應(yīng)的補 償積累后頻域分布作為該距離采樣點的頻域分布。具體的步驟(三）中多項式拉東-多項式 傅里葉變換又包括
[0073] 1)令目標(biāo)徑向速度估計初始值化.=()__，_徑向加速度估計初始值1 = 〇_，1 = 1，Δ VI = ΝνΔ vr, Δ af = A/(4Τ2)；
[0074] 2)令I(lǐng) -1 + 1，如果Δ ν〗> Δ vr，
，否則Δ ν〗=Δ vr，nv表示根據(jù) 搜索實時性要求設(shè)置的搜索步進分辨率，nv取值范圍為[0，NV];
[0075] 4)如果 Δ ai> Δ af，令
，否則 Δ ai= Δ af，na取值范圍為[0，Na];
[0076] 5)F0R對任一i，i = l，2, · · ·，2nv_l，遍歷所有i
[0077] 令《，<- + (―. 6?)
[0078] FOR對任一k，k = 1，2，· · ·，2na_l，遍歷所有k
[0079] 令.& +. + .. Δ%.)
[0080] FOR對任一n，n = l，2, · · ·，N，遍歷所有η
[0081]
[0082] 其中Tn=l/PRF，PRF表示脈沖重復(fù)頻率;令
，進一步分別計算
[0083]
[0084] ENDF0R
[0085]利用FFT算法，得到第m個距離采樣點在參數(shù)(Vl，aJ下的距離-多普勒分布 [0086] G[rm，（Vi，aj)]=max( |FFT(P'）| )
[0087]其中，rm表示第m個距離采樣點對應(yīng)的距離，F(xiàn)FT( ·)表示對括號內(nèi)的向量進行FFT 變換；
[0088] ENDF0R
[0089] ENDF0R
[0090] 6)對第m個信號采樣點，找出（Vi，a」)參數(shù)域中最大值對應(yīng)的參數(shù)，
[0091]
[0092] 令目標(biāo)的徑向速度估計兔=包，令目標(biāo)的徑向速度估計^ ^如果Λ νι> Δ心或Δ ai> Δ af，則跳回第2)步，否則進入下一步；
[0093] 7)將參數(shù)(&Λ)對應(yīng)的多普勒分布P7 i,j(η, 1)作為該信號米樣點的多普勒分布； [0094]步驟（四）、重復(fù)步驟(三），遍歷所有信號采樣點，得到距離-PRPFT域的檢測時頻分 布圖 Glr^f)，
[0095]步驟(五）、對時頻分布圖進行恒虛警檢測，判斷目標(biāo)有無，具體方法為：
[0096] 對C(rm，f)與恒虛警門限比較
[0097]
[0098] 式