孫 濤，謝曉方，唐 江，張龍杰

（1.海軍航空工程學(xué)院 兵器科學(xué)與技術(shù)系，山東 煙臺 264001；2.海軍航空工程學(xué)院 研究生管理大隊，山東 煙臺 264001）

空艦導(dǎo)彈純方位攻擊是指僅根據(jù)目標方位信息，采取導(dǎo)彈自控飛行時間裝定最小值，向目標初始方位發(fā)射導(dǎo)彈的攻擊方式。這種攻擊方式可以更好地適應(yīng)復(fù)雜的電磁環(huán)境［1］和惡劣的信息條件。由于采用了被動探測設(shè)備，這種攻擊方式也更有利于提高發(fā)射平臺的隱蔽性。目前，美國的最新一代戰(zhàn)機上均安裝了被動探測設(shè)備［2－3］，主動雷達僅在必要時短時間開機；俄羅斯和以色列等國在這方面也有很強的技術(shù)實力。被動探測手段的重要性日益突出。

被動探測設(shè)備通常只具有測向能力，如何在這種條件下科學(xué)處理目標方位數(shù)據(jù)，利用空艦導(dǎo)彈有效進行目標搜索攻擊，已成為航空武器系統(tǒng)與運用工程領(lǐng)域的熱點問題。

在院校教學(xué)中，仿真實驗有利于靈活配置實驗內(nèi)容，使理論教學(xué)和實踐教學(xué)相融合［4－5］。為了使學(xué)生更好地理解純方位目標指示數(shù)據(jù)處理的一系列模型和方法，有必要構(gòu)建綜合仿真實驗系統(tǒng)，為準確、直觀地研究純方位目標指示條件下數(shù)據(jù)處理方法提供條件。

1 實驗系統(tǒng)功能需求

仿真實驗系統(tǒng)主要支持研究性實驗，服務(wù)于探究性教學(xué)和實驗比較法教學(xué)［6－8］，具備以下主要功能：

（1）支持主要的純方位目標指示數(shù)據(jù)處理方法仿真，包括幾何法、概率法、角速率法等，提供參數(shù)設(shè)置接口，實驗用戶可以自行修改目標、載機等參數(shù)；

（2）針對各類方法仿真結(jié)果，解算空艦導(dǎo)彈火控參數(shù)，支持目標捕捉概率的仿真計算；

（3）支持仿真過程記錄和結(jié)果比對，滿足實驗教學(xué)與理論教學(xué)相結(jié)合的需求。

（4）提供模型開發(fā)接口，既能夠為研究生等高級實驗用戶提供自定義仿真計算模型的實驗條件，又支持本科生研究型實驗，培養(yǎng)學(xué)生的創(chuàng)新意識和科研品質(zhì)［9］；

2 主要仿真模型分析

2.1 幾何法處理模型

幾何法的基本原理與測向定位法（bearing－only，BO）［10－11］相同，通過在2個以上不同位置對目標方向的多次被動探測，解得目標的位置信息，從而為空艦導(dǎo)彈提供目標指示。設(shè)機載被動探測設(shè)備的最大探測距離為Rmax，目標方向指示的角誤差為Δα，則目標指示為圖1所示的虛線扇形區(qū)域。

圖1 基于被動測向的目標定位

在一定時間內(nèi)對2個以上目標指示扇區(qū)的相交運算，可以為導(dǎo)彈計算出更加準確的目標區(qū)。隨著探測次數(shù)的增加，目標區(qū)域不斷縮小。通常情況下，在不同位置進行的被動探測很難做到同時進行，由此會在目標區(qū)計算中引入誤差。先后在P1和P2兩個不同位置進行被動探測的時間間隔為Δt，目標艦艇最大可能的機動速度為Vm，則誤差補償方法是將搜索區(qū)沿P2探測的正反方向分別延伸Δt·Vm，如圖2所示。

圖2 測向時間間隔補償

設(shè)導(dǎo)彈到達目標區(qū)之前的自控飛行時間為tZK，則期間艦艇的最大機動距離為tZK·Vm，補償方法是按tZK·Vm將搜索區(qū)向外擴展［12］，如圖3所示。

圖3 導(dǎo)彈自控飛行時間補償

2.2 概率法處理模型

載機在不同位置P1，P2，…，Pi，…，Pn，通過無源測向設(shè)備進行方位探測，分別得到目標相對于不同觀測點的方位角φ1，φ2，…，φi，…，φn。通過以上信息估計目標位置。設(shè)目標T的真實位置為η＝（xT，yT），載機的測向點的位置坐標分別為（xp1，yp1），（xp2，yp2），…，（xpn，ypn），對于任意測向點pi（1≤i≤n），觀測值為

δi是第i次測向誤差。

設(shè)載機各次測向噪聲不相關(guān)，為零均值，標準差為σi的隨機噪聲。表示方位角真值，則有），對于距離r的目標，φi的概率密度函數(shù)為：

其中，Rmax為無源探測設(shè)備的最大探測距離，Rmin為無源探測設(shè)備的最小探測距離，則目標在探測區(qū)域的概率密度函數(shù)表示為：

若對同一目標進行N次測量，各次測量相互獨立，則聯(lián)合概率密度函數(shù)為：

2.3 方位角速率法處理模型

本方法適用于純方位目標的快速距離估計和定位誤差估計，適用條件是在數(shù)據(jù)采集及處理期間載機保持平飛，并且被動傳感器能夠被目標連續(xù)照射。

假設(shè)載機與目標實際距離為r，載機飛行速度vP，載機方位指示分辨率εmax，將εmax視為載機對目標方位探測的量化單位。令Δn為單位時間Δt內(nèi)載機的目標方位角示值變化量對應(yīng)的量化單位個數(shù)?？紤]到載機可能無法始終沿以目標點為中心的圓弧切向飛行，需要進行飛行方向補償：

其中，θ為目標相對于載機飛行方向的方位角。該式表明，在已知載機飛行速度和載機方位指示分辨率的條件下，目標相對于載機的距離與方位示值變化率成反比關(guān)系。

對于空艦導(dǎo)彈，目標運動速度相對于載機速度很小，利用對稱原理，視載機不動，而目標艦艇相對于載機以速度vP運動（見圖4）。

圖4 載機－目標相對運動與目標方位測量

圖中θ變化間隔為最小量化單位εmax，下標n為測量周期序號。

估計結(jié)果與方位示值變化周期的正確性和精度有直接關(guān)系。然而，載機上能夠測得的只有各次測量之間的平均角速度，不能反映被測目標當(dāng)前的實際角速度。相鄰兩測量時刻，目標方位角速度ωn，ωn＋1存在如下關(guān)系：

平均角速率比測量點當(dāng)前實際值小，導(dǎo)致目標距離預(yù)測結(jié)果偏遠。因此利用平均角速度作為目標距離估計的模型。

載機平飛條件下，如果把角速度看作目標方位角的函數(shù)，則有：

又因為，rn＋1＝a／cosθn＋1，所以，第n＋1次測量時，目標距離估計為：

2.4 目標捕獲仿真

在目標捕獲概率方面，主要需考慮導(dǎo)彈進入搜索區(qū)的方向、自控終點的位置、所用的彈量，依照目標區(qū)內(nèi)部目標概率密度分布情況進行積分求解（見圖5）。而在仿真結(jié)果方面，則主要考慮對目標區(qū)的搜索覆蓋率和用彈效費比。

3 仿真程序總體流程

綜合仿真系統(tǒng)總體流程如圖6所示。其中，幾何法、概率法、方位角速率法計算模型作為開放資源預(yù)先存儲，實驗用戶可直接調(diào)用，也可自行開發(fā)設(shè)計；初始仿真參數(shù)和導(dǎo)彈攻擊參數(shù)可由用戶設(shè)置。

圖5 目標區(qū)搜索

圖6 載機純方位無源定位算法流程圖

4 部分仿真實驗效果

4.1 幾何法

設(shè)置被動探測角誤差10°，載機運動速度為272 m／s，仿真結(jié)果如圖7所示。圖中飛機圖標表示載機平臺進行各次被動測向的位置，多邊形表示經(jīng)過多次探測之后，仿真計算得到的目標區(qū)。利用幾何法得到的目標區(qū)可能是任意形狀的凸多邊行，也可提取如圖所示的矩形目標區(qū)，同時可給出載機完成所有方位探測所需要的機動距離和飛行時間。

圖7 提取的目標搜索區(qū)

4.2 概率法

選取載機初次探測位置為原點，φ＝0°，Rmax＝400 km，Rmin＝10km，δ＝3°，則單次測量目標概率密度分布如圖8所示。

圖8 單次測量目標概率分布

設(shè)載機對目標方位進行3次探測，將坐標系原點取在第二次探測位置，載機沿x軸正向水平飛行。假設(shè)目標位于距原點r＝100km，φ＝90°處，探測間隔距離為L＝10km。概率法仿真結(jié)果如圖9所示。

圖9 目標聯(lián)合離散概率密度分布

4.3 方位角速率法

不同情況下目標估計距離r與方位角示值變化率仿真及對比效果如圖10所示。

圖10 方位指示數(shù)據(jù)變化周期與目標距離的關(guān)系

令vP＝300m／s，＝10－10rad2／s2（平均角速度方差），εmax＝π／512，＝4×10－4rad2，＝10 m2／s2，＝3.765×10－9?？捎嬎憔嚯x估計方差隨θn＋1，rn＋1變化情況，仿真結(jié)果如圖11所示。

圖11 距離估計方差（無源雷達）

4.4 目標捕獲

導(dǎo)彈自控終點誤差3.5km，搜索扇面張角80°，搜索距離20km，采用2枚導(dǎo)彈對目標區(qū)進行搜索攻擊，仿真效果如圖12所示。

圖12 目標區(qū)搜索仿真結(jié)果

5 結(jié)束語

綜合仿真實驗系統(tǒng)提供了各種處理方法的標準仿真模型，統(tǒng)一了仿真實驗環(huán)境和結(jié)果表示形式，使得各種方法的處理效果具有可比性。實驗用戶不必將主要精力放在仿真環(huán)境構(gòu)建、模型驗證、仿真程序編制等方面，而是通過不同條件下的仿真參數(shù)設(shè)置獲得直觀的仿真結(jié)果并進行對比分析，從而獲得更好的教學(xué)效果。

（

）

［1］田錦昌，袁健全，孔玉，等.空艦導(dǎo)彈制導(dǎo)雷達適應(yīng)復(fù)雜電磁環(huán)境方法分析［J］.航天電子對抗，2010，26（6）：12－13.

［2］芮國勝，康健.機載被動定位系統(tǒng)仿真研究［J］.計算機仿真，2003，20（9）：27－28.

［3］楊建華，陳思興.國外機載快速無源定位技術(shù)最新發(fā)展［C］／／雷達跟蹤技術(shù)研討會論文集.2001：38－40.

［4］凌霖，蘇胤杰，晉春.軟件仿真在通信原理實驗教學(xué)中的應(yīng)用［J］.實驗技術(shù)與管理，2011，28（4）：83－85.

［5］劉伯紅，閻英，方義秋.高校計算機實踐教學(xué)質(zhì)量保障體系改革探索與實踐［J］.實驗室研究與探索，2012，31（12）：121－123.

［6］朱昌平，施鈴泉，王斌，等.講授好“卓越計劃”高頻電路理論與實踐課程的探索［J］.實驗技術(shù)與管理，2012，29（10）：142－145.

［7］馮毅平，張光新，榮岡.基于知識系統(tǒng)化重構(gòu)的研究型實驗教學(xué)［J］.實驗技術(shù)與管理，2012，29（10）：154－156.

［8］徐凱，徐果微.基于“實驗比較法”的電氣信息類研究生培養(yǎng)教學(xué)探究［J］.實驗室研究與探索，2012，31（9）：135－138.

［9］張曄.基礎(chǔ)實驗中科研能力的培養(yǎng)［J］.實驗室研究與探索，2012，31（12）：91－92.

［10］Taek L Song，Tae Yoon Um.Practical guidance for homing missiles with Bearing－Only measurements［J］.IEEE Transactions OD Aerospace and E1ectronic Systems.1996，32（1）：434－443.

［11］Hammel S E，Aidala V J.Observability requirements for three－dimensional tracking via angle measurements［J］.IEEE Transactions on Aerospace and Electronic Systems.1985，21（2）：200－207.

［12］孫濤，謝曉方，梁捷，等.純方位目指條件下空艦導(dǎo)彈目標搜索策略研究［J］.彈箭與制導(dǎo)學(xué)報，2011，31（5）：72－74.