宋文瑰 沈安慰
(1.空軍工程大學(xué)航空航天工程學(xué)院,陜西 西安 710043;2.石家莊海山實業(yè)發(fā)展總公司,河北 石家莊 050208)
基于MATLAB GUI專利技術(shù)的雷達探測概率軟件平臺設(shè)計
宋文瑰1,2沈安慰1
(1.空軍工程大學(xué)航空航天工程學(xué)院,陜西 西安 710043;2.石家莊海山實業(yè)發(fā)展總公司,河北 石家莊 050208)
飛機的雷達隱身性能是飛機隱身性能的重要組成部分,是現(xiàn)階段研究的重點領(lǐng)域,當(dāng)今世界各國都在致力于這方面的研究,這些專利研究成果將對空戰(zhàn)的勝負產(chǎn)生決定性的影響。文中飛機雷達隱身性能的計算,采用從簡到繁的思路,首先對雷達基本探測原理進行分析計算,然后針對RCS不變時的檢測概率和可變時檢測概率分別進行計算,對于檢測概率與探測距離之間難以建立定量關(guān)系的問題,我們引入了中間變量信噪比(SNR),從而建立了間接的函數(shù)關(guān)系,使問題的分析得到簡化。
飛機雷達隱身性能;Swerling;SNR;MATLAB;軟件編程
從海灣戰(zhàn)爭到1999年科索沃戰(zhàn)爭,再從科索沃戰(zhàn)爭到2003年伊拉克戰(zhàn)爭,F(xiàn)-117A和B-2A隱身飛機在現(xiàn)代高科技戰(zhàn)場上展示了它出眾的突防能力和作戰(zhàn)效能[1-2]。隱身飛機的出現(xiàn)給現(xiàn)代雷達防空系統(tǒng)提出了巨大的挑戰(zhàn),對作戰(zhàn)樣式和防御系統(tǒng)建設(shè)產(chǎn)生了重大影響。由于隱身飛機對預(yù)警探測系統(tǒng)具有隱身功能,使敵方戰(zhàn)斗準(zhǔn)備和反應(yīng)時間減少到來不及應(yīng)對的程度,使戰(zhàn)役、戰(zhàn)斗進攻的主動性和突然性得到充分地發(fā)揮,在空襲時機、空襲目標(biāo)、空襲手段和方式等方面都有較大的選擇余地[3-4]。
但據(jù)目前統(tǒng)計分析,雷達依然是飛行器最致命的殺手,對飛機的威脅占到60%以上,包括機載火控雷達、記載預(yù)警偵察雷達、地面制導(dǎo)火控雷達和地面警戒或監(jiān)視雷達等[5-7]。從1999年科索沃戰(zhàn)場上F-117A隱身戰(zhàn)斗機看,隱身飛機并不是完全不可探測,它存在一定的弱點和局限性。隨著反隱身技術(shù)的發(fā)展,尤其是全數(shù)字化的有源相控陣技術(shù)在地空警戒或監(jiān)視雷達中的應(yīng)用,雷達的角度和距離精度都有了大幅度的提高,利用飛機在VHF 和UHF波段的RCS增加的特點,使得VHF和UHF波段的雷達似乎成為了對抗隱身飛機利器。所以在面對現(xiàn)代高科技中行之有效的防空預(yù)警探測系統(tǒng),特別是面對作為整個探測系統(tǒng)的核心裝備的雷達系統(tǒng),隱身性能成為現(xiàn)階段戰(zhàn)斗機設(shè)計制造關(guān)鍵性能,而研究飛機雷達隱身性能,就可以更好提高飛機生存力,可以使飛機高效地規(guī)避雷達的探測,就像是一只黑色幽靈突然降臨敵人的上空,實現(xiàn)隱蔽突防,奪取戰(zhàn)爭的主動權(quán)[8]。
1.1 機RCS的仿真建模計算
圖1
由于空中目標(biāo)對雷達波反向特性極其復(fù)雜,準(zhǔn)確地推導(dǎo)出地面(或機載)雷達在各種條件下對其探測概率是不可能的,所以必須進行近似計算。
理論上講,任意空中目標(biāo)的RCS值都是可以通過試驗測定的,但要做大量測試,其費用、安全性等因素又決定這是不現(xiàn)實的。所以,許多飛機往往只有一個或幾個典型狀態(tài)下的RCS值。進行仿真研究時,就要通過若干個數(shù)據(jù),建立擬合模型,最終得到通過的計算模型。
雷達相對空中目標(biāo)(飛機)的位置是以目標(biāo)的機體坐標(biāo)系為基準(zhǔn)的,如圖1所示。
從改點探測目標(biāo)時的角度定義如下:
φ——探測方位角,雷達與目標(biāo)連線(OA)在機體坐標(biāo)系oxyz平面內(nèi)投影OA’與ox軸的夾角,由ox軸逆時針旋轉(zhuǎn)至OA’時,φ為正,反之為負。
φ——探測俯仰角,雷達與目標(biāo)連線(OA)在機體坐標(biāo)系oxyz平面內(nèi)投影OA’之間的夾角,
因飛機一般是左右對稱的,左右探測時RCS值相同,或變化很小;但從前后方、上下方探測時是有很大區(qū)別的。
設(shè)從目標(biāo)正前方、正后方、正上方、正下方、正側(cè)方探測時RCS分別為σx,σy,σz,則給出RCS全向計算公式:
一般來講,通過查詢可以獲得目標(biāo)的典型值σx,σ-x,σy,σ-y,σz即可求得任意方向探測目標(biāo)時的RCS 值σ(φ,φ)
1.2 虛警概率
虛警概率Pfa定義為當(dāng)雷達中只有噪聲信號時,信號r(t)的樣本R超過門限電壓的概率。
虛警時間Tfa與虛警概率的關(guān)系為
其中,tint表示雷達的累積時間,或包絡(luò)檢波器的輸出超過門限電壓的平均時間。因此雷達的工作帶寬B是tint的逆,所以將式(2)代入式(3),可以將Tfa寫為
使虛警時間最小意味著增加門限值,結(jié)果雷達的最大檢測距離會減少。因此,Tfa可接受數(shù)值的選擇變成依賴雷達工作模式的折中[5-6]。
1.3 檢測概率
我們首先研究目標(biāo)截面積不變且不考慮累計損失的雷達探測概率。檢測概率是r(t)的樣本R在噪聲加信號的情況下超過門限電壓的概率[5-7]。
如果假設(shè)雷達信號是幅值為A的正弦波形,那么它的功率是/2?,F(xiàn)在,使用(單個脈沖的SNR)和,那么上式可以重寫為:
Q稱為MarcumQ函數(shù)。當(dāng)Pfa較小,PD相對較大,從而門限也較大時,式(6)可以近似寫為:
2.1 飛機生存力評估模型
飛機生存力是指飛機躲避或承受人為敵對環(huán)境的能力,即飛機不被擊中或擊中后能繼續(xù)執(zhí)行戰(zhàn)斗任務(wù)、控制飛行、完成強迫著陸等能力。飛機生存力由飛機的敏感性和易損性兩個基本要素構(gòu)成,可表示為
敏感性是指作戰(zhàn)飛機在完成任務(wù)過程中被殺傷機理擊中的可能性,用飛機被殺傷機理命中的概率PH來表示
式中PD為飛機被敵方威脅發(fā)現(xiàn)的概率,PT為被敵方威脅跟蹤的概率,PLGH為飛機被敵方武器系統(tǒng)成功跟蹤、制導(dǎo)并擊中的概率,該參數(shù)主要取決于發(fā)射裝置的效能、制導(dǎo)精度(或瞄準(zhǔn)精度)。在本軟件中PLGH可作為敵方威脅(主要為導(dǎo)彈)的屬性,直接作為輸入?yún)?shù)設(shè)定即可。
敵方威脅探測器主要分為三類,雷達探測器、紅外探測器和射頻探測器。而飛機面對的威脅情況可能是其中的一種或兩種或三種同時面對。因此飛機被敵方威脅發(fā)現(xiàn)并跟蹤的概率PDT=PD×PT可由下式計算
飛機在執(zhí)行一個完整的任務(wù)時,可能會遇到若干個形形色色的探測器組合。因此,在一個具體任務(wù)中,有若干個飛機被威脅探測并跟蹤的概率PDT。由上述公式即可計算出多次遭遇下的飛機生存力PS。假設(shè)一個任務(wù)下飛機經(jīng)歷了k次遭遇,因此對于一個完整任務(wù)下的生存力P
本軟件采用仿真方法來研究,即當(dāng)飛機遭遇威脅時,根據(jù)探測概率、跟蹤概率、擊中概率和易損性概率的不同,隨機確定飛機在某次仿真中到底能不能成功執(zhí)行一個完整的任務(wù)流程。整個仿真程序可用以下步驟完成:
步驟1:設(shè)定好紅藍雙方的兵力部署。暫考慮紅方只有一架飛機,藍方有若干陸基、?;涂栈綔y器和地空導(dǎo)彈,探測器包括雷達、紅外和射頻探測器。陸基和海基探測器默認是不動的,而空基(預(yù)警機)會以一定的概率隨機出現(xiàn)。只有飛機被上述預(yù)警探測設(shè)備發(fā)現(xiàn)并跟蹤之后,然后進入地空導(dǎo)彈的打擊范圍時,就有一定概率被擊中。(此時暫不考慮藍方有戰(zhàn)斗機出現(xiàn))。
步驟2:飛機按照指定路線飛行(飛機的任務(wù)一般是從基地起飛到達指定目標(biāo)點投彈后返回基地的一個過程),當(dāng)飛機第一次遭遇敵方威脅時,根據(jù)概率判斷飛機是否被發(fā)現(xiàn),若發(fā)現(xiàn)判斷飛機是否被擊中,若擊中判斷飛機是否毀傷,若飛機毀傷則終止本次仿真;否則繼續(xù)進行飛機的下一次遭遇,直到飛機完成整個飛行任務(wù)。若飛機完成飛行任務(wù),則計任務(wù)成功一次,否則計任務(wù)失敗一次。若任務(wù)失敗了,需記錄飛機毀傷時的時刻,以及被雷達/紅外/射頻中的哪種探測設(shè)備發(fā)現(xiàn)的。
步驟3:重復(fù)循環(huán)上述步驟N次,累積求和得到成功的次數(shù),除以總循環(huán)次數(shù)N即可得到飛機生存力。
2.2 軟件總體設(shè)計
雷達探測概率軟件平臺基于MATLAB GUI進行設(shè)計主要是為了依托MATLAB環(huán)境提供的基本科學(xué)計算減少開發(fā)周期,更為重要的是MATLAB開發(fā)環(huán)境提供了創(chuàng)建用戶界面的捷徑GUIDE,有Windows基本控件的支持并且有良好事件的驅(qū)動機制同時提供MATLAB數(shù)學(xué)庫的接口,還可以方便地創(chuàng)建各種圖形句柄對象,實現(xiàn)仿真平臺的用戶界面GUIDE可以根據(jù)用戶GUI的版面設(shè)計過程直接自動生成M文件框架,這樣就簡化了GUI應(yīng)用程序的創(chuàng)建工作,用戶可以直接使用這個框架來編寫自己的函數(shù)代碼【9】。
本軟件的界面布局設(shè)計采用自頂向下的設(shè)計方法即先設(shè)計主界面再設(shè)計各個子界面,界面的模塊主要是根據(jù)文章的內(nèi)容結(jié)構(gòu)來進行劃分各模塊,所包含的主要內(nèi)容雖不能涵蓋每個知識點,但是包含了其中的大部分重難點并且在人機界面中進行了可視化實現(xiàn)。
本文對Swerling模型進行比較深入的分析研究,并把該模型用于單雷達發(fā)現(xiàn)概率的仿真計算,對Swerling模型的理論進行實踐應(yīng)用,用MATLAB語言又編制了相應(yīng)的計算軟件,軟件的功能集成了典型的雷達探測概率模型,對隱身飛機雷達探測概率的計算提供了方便。
[1]師俊朋,胡國平,朱蘇北,羅亞宗.雷達反隱身技術(shù)分析及進展[J].現(xiàn)代防御技術(shù),2015,43(6):124-130.
[2]蔣忠進,趙書敏,耿江東,李希同.基于雷達成像的地面停留隱身飛機探測[J].雷達科學(xué)與技術(shù),2016(1):54-58.
[3]陶星.隱身通信系統(tǒng)中部分關(guān)鍵技術(shù)的研究與實現(xiàn)[D].武漢:華中師范大學(xué),2016.
[4]葉圣天,劉朝輝,成聲月,班國東.國內(nèi)外紅外隱身材料研究進展[J].激光與紅外,2015,45(11):1285-1291.
[5]張改平,屈紹波.紅外隱身技術(shù)研究新進展[J].航天電子對抗,2012,28(5):41-43.
[6]秦娜.裝甲車輛在紅外隱身措施下的仿真評估[D].南京:南京理工大學(xué),2016.
[7]唐建楊.基于Hilbert_Huang變換的隱身目標(biāo)涂層損傷檢測研究[D].成都:電子科技大學(xué),2016.
[8]劉斌,張翔,周建江,廖文和.一種可實現(xiàn)雷達隱身的微小衛(wèi)星外形設(shè)計(英文)[J].宇航學(xué)報,2011,32(12):2628-2634.
[9]張志剛,劉麗梅,等.Matlab與數(shù)學(xué)實驗[M].北京:中國鐵道出版社,2004.
Radar Detection Probability of Software Platform Design Based on MATLAB GUI Patent Technology
Song Wengui1,2Shen Anwei1
(1.College of Aeronautics and Astronautics Engineering,Air Force Engineering University,Xi'an Shanxi 710043;2.Shijiazhuang Haishan Industrial Development Corporation,Shijiazhuang Hebei 020208)
The radar stealth performance of aircraft is an important part of aircraft stealth performance,is the research focus of understanding at this stage,in today's world,all countries are committed to the research of this area,the results of the study will have a decisive impact on the outcome of the air combat.Calculation of aircraft stealth performance,using ideas from simple to complex,first of all the basic principle of radar detection were analysis and calculation,then the RCS unchanged when the detection probability and variable detection probability are calculated respectively.For the problem of the quantitative relationship between the detection probability between the detection probability and the detection distance,we introduce the intermediate output signal to noise ratio,thus the indirect function relation is established,and the problem is simplified.
stealth performance of aircraft radar;swerling model;SNR;MATLAB;software programmer
TN957.51
A
1003-5168(2016)07-0065-04
2016-7-1
宋文瑰(1983-),女,碩士研究生,工程師,研究方向:機載設(shè)備修理。