吳思勵,孔 玲,衣向舉,孔全存
(1.93277部隊,沈陽 110000;2.93303部隊,沈陽 110000;3.北京信息科技大學,北京 100000)
在探測裝備使用過程中,以雷達裝備為例,由于外界環(huán)境的作用以及元器件頻譜相互影響,都能使實際探測威力發(fā)生改變,所以有必要對探測威力圖定期修正。本文在上述概率計算方法基礎上,提出基于海量云數(shù)據(jù)的探測概率計算方法,對飛行軌跡進行網格分割,在已具備探測裝備實際威力圖區(qū)間參數(shù)情況下,根據(jù)塔臺飛行自動化系統(tǒng)的情報數(shù)據(jù)(云數(shù)據(jù))進行數(shù)據(jù)比對,利用數(shù)理分析的方法計算出任意空間航測點實際探測區(qū)間,實現(xiàn)探測威力范圍的再次修正。
目標反射截面積(RCS),它是一個等效面積(以雷達為例,用σ表示目標截面積,σ是基本方程中的一個因子)。由波長、極化方式、飛行器參數(shù)及相對視向角等多個變量組成了相關符合函數(shù),表述最大發(fā)現(xiàn)距離的公式為:
式中,G為雷達天線增益;Pt為雷達發(fā)射脈沖功率;λ為雷達波長;k為波爾茲曼常數(shù);T0為工作噪聲溫度;Fn為噪聲系數(shù);Bn為噪聲帶寬;(S/N)min為檢測器輸入端所需最小信噪比;L為總損耗系數(shù)。
并且當同一雷達對不同飛行物進行探測時,在其他參數(shù)相同且等信噪比條件下,不同飛行物的最大發(fā)現(xiàn)距離與目標反射截面積有如下關系:
式中,R1,R2分別為不同飛行物的最大發(fā)現(xiàn)距離;σ1,σ2分別為反射目標截面積。在已知各種型號飛行器的雷達目標反射面積的條件下,可以算出最大發(fā)現(xiàn)距離。
3.1.1 計算實際發(fā)現(xiàn)概率
根據(jù)雷達掃描周期和云數(shù)據(jù)航跡參數(shù),統(tǒng)計發(fā)現(xiàn)點和未發(fā)現(xiàn)點,然后分段計算發(fā)現(xiàn)概率。把距離取樣間隔△D分段,在航跡不同高度區(qū)間內,各距離取樣間隔發(fā)現(xiàn)概率P表述為:
式中,M為距離取樣間隔內的發(fā)現(xiàn)點數(shù);N為觀測量。需要說明的是,將距離取樣間隔交點處作為較近距離間隔內的觀測點,并且只統(tǒng)計計入一次。
最后將計算結果即雷達探測概率值合并填入空間分割數(shù)據(jù)單元中。圖1為Matlab仿真結果。
圖1 Matlab仿真結果
3.1.2 修正雷達最大探測距離
(1)在上述得到的發(fā)現(xiàn)概率曲線仿真結果中,取發(fā)現(xiàn)概率P=0.5的兩個航跡點。還是以雷達作為探測裝備為例,一個航跡點是波瓣在某高度發(fā)現(xiàn)目標的探測距離,另一個航跡點是波瓣在某高度丟失目標的探測距離,分別用rd,rv(km)表示,同理探測裝備對飛行物發(fā)現(xiàn)點和消失點相對高度分別用Hd,Hv(m)表示,探測裝備對飛行物發(fā)現(xiàn)和消失時的仰角用θd,θv(度)表示,Re為地球半徑。
兩種情況下其仰角公式如下:
(2)設定探測距離為R0、飛行物視在仰角θ0,F(xiàn)(θ0)為該仰角時傳播因子。探測裝備最大探測距離Rmax公式如下:
3.1.3 繪制概率曲線
建立坐標系,距離R為橫坐標,發(fā)現(xiàn)概率P為縱坐標,隨著距離增加發(fā)現(xiàn)概率值發(fā)生變化,按已得到的發(fā)現(xiàn)概率可以畫出不同高度區(qū)間P與R的關系曲線。關系曲線逐漸平順后,可以查看比對出與飛行目標發(fā)現(xiàn)概率P0所對應探測裝備的最大探測距離R0。
針對起伏動目標,垂直波瓣圖是雷達探測區(qū)和盲區(qū)的曲線圖。是在給定檢驗飛行器平均截面積、虛警概率和探測概率條件下通過經常性飛行數(shù)據(jù)檢驗獲得的。對比可知,在3000米高度,雷達的最大作用距離從300公里減小到了290公里。同時用Matlab對垂直波瓣進行修正,通過最小二乘多項式曲線擬合求出擬合曲線方程,得出垂直波瓣擬合曲線,以實現(xiàn)探測范圍直觀可視。
不同于以往探測概率計算方法,本文是基于塔臺飛行自動化系統(tǒng)中日常積累的海量云數(shù)據(jù)的背景下,云數(shù)據(jù)能準確的反映雷達長期以來的實測概率區(qū)間。在今后陸續(xù)的實測工作中,可通過數(shù)據(jù)挖掘從這些海量云數(shù)據(jù)中找出更多可借鑒的數(shù)據(jù)參數(shù),在諸如定位精度、氣象環(huán)境、電磁干擾等相關領域繼續(xù)深入研究。