李 源， 張京國(guó)， 梁 謙

（中國(guó)空空導(dǎo)彈研究院，河南 洛陽(yáng)471009）

0 引言

近炸引信虛擬樣機(jī)是在計(jì)算機(jī)上建立的表示引信系統(tǒng)（包括引戰(zhàn)配合特性）及其使用環(huán)境的完整、詳細(xì)的模型，它能夠逼真地反映引信在測(cè)試、試驗(yàn)和作戰(zhàn)使用中的性能，并且它能夠用來(lái)評(píng)價(jià)引信的設(shè)計(jì)方案和評(píng)估產(chǎn)品的性能［1］。由于虛擬樣機(jī)的以上優(yōu)點(diǎn)，隨著計(jì)算機(jī)技術(shù)的發(fā)展，近炸引信虛擬樣機(jī)設(shè)計(jì)已經(jīng)成為引信設(shè)計(jì)中的一個(gè)重要階段。

作為激光引信的一個(gè)重要組成部分，激光引信虛擬樣機(jī)必然要對(duì)光學(xué)系統(tǒng)進(jìn)行建模。傳統(tǒng)的光學(xué)仿真往往采用商業(yè)的光學(xué)軟件來(lái)實(shí)現(xiàn)，但對(duì)于激光引信而言，其采用的光學(xué)系統(tǒng)形式不同于傳統(tǒng)光學(xué)系統(tǒng)，特別是用來(lái)進(jìn)行光電轉(zhuǎn)換的探測(cè)器往往具有不規(guī)則面形，而光學(xué)軟件只能針對(duì)光學(xué)系統(tǒng)的某個(gè)功能進(jìn)行仿真以及優(yōu)化，對(duì)目標(biāo)特性以及不規(guī)則面形探測(cè)器都不能進(jìn)行仿真?；诖朔N原因，本文采用統(tǒng)計(jì)建模與光學(xué)追跡的方法，利用Matlab軟件編程對(duì)此光學(xué)系統(tǒng)進(jìn)行虛擬樣機(jī)建模。

1 發(fā)射光學(xué)系統(tǒng)的數(shù)學(xué)建模

1.1 發(fā)射光學(xué)系統(tǒng)參數(shù)

發(fā)射光學(xué)系統(tǒng)子午面視場(chǎng)角α；

發(fā)射光學(xué)系統(tǒng)弧矢面視場(chǎng)角β；

激光源：量子阱陣列GaAs半導(dǎo)體激光器。

圖1所示為光學(xué)發(fā)射機(jī)的組成示意圖。

圖1 光學(xué)發(fā)射機(jī)的組成示意圖

1.2 發(fā)射光束的數(shù)學(xué)統(tǒng)計(jì)建模

按照發(fā)射光學(xué)系統(tǒng)要求，對(duì)激光器發(fā)出的激光束進(jìn)行整形設(shè)計(jì)［2］，如圖2所示。

圖2 發(fā)射光學(xué)系統(tǒng)

激光器發(fā)出的激光經(jīng)過(guò)光學(xué)系統(tǒng)整形后，在子午方向光束能量P（α）是高斯分布；在弧矢方向由于激光器發(fā)光結(jié)是陣列分布，且又經(jīng)光學(xué)系統(tǒng)擴(kuò)束，所以光束能量P（β）在此方向即不是高斯分布，又不是理想的均勻分布。這就給光學(xué)發(fā)射機(jī)的激光光束建模帶來(lái)了難度。所以，在本文中借助試驗(yàn)參數(shù)，運(yùn)用統(tǒng)計(jì)建模法將多個(gè)測(cè)試結(jié)果進(jìn)行融合，最終得到發(fā)射光學(xué)系統(tǒng)的能量空間分布。

圖3是光學(xué)發(fā)射機(jī)在子午面以及弧矢面的能量分布示意圖。

圖3 光學(xué)發(fā)射機(jī)在子午面和弧矢面的能量分布示意圖

上述情況僅考慮了光束能量的空間分布，而在本文中還要考慮激光束功率隨時(shí)間的變化。本文所研究的激光引信采用的是脈沖工作體制，直接向半導(dǎo)體激光器注入電流，改變發(fā)射波形。激光光束功率P（t）隨時(shí)間的變化如圖4所示。

圖4 光學(xué)發(fā)射機(jī)光束能量隨時(shí)間分布示意圖

假設(shè)由光學(xué)系統(tǒng)光面發(fā)射到外部空間的功率為PS，將此功率在子午面與弧矢面內(nèi)進(jìn)行離散化，離散化后的單位發(fā)射光束與目標(biāo)反射面的截面稱為單位反射面元，如圖5、圖6所示。

由上圖可知，單位面元面積 ATAR＝

圖5 光學(xué)發(fā)射機(jī)在子午面的能量均分示意圖

圖6 光學(xué)發(fā)射機(jī)在弧矢面的能量均分示意圖

（π／4）R2θαθβ，其中R 為發(fā)射機(jī)到目標(biāo)間的距離；θα、θβ分別為單位發(fā)射光束在子午面以及弧矢面的角寬度。照射到單位面元上的功率PTAR為

式中：A（α，β，t）為光束能量PS隨子午面視場(chǎng)角、弧矢面視場(chǎng)角以及時(shí)間變化的系數(shù)；τ 為發(fā)射機(jī)到目標(biāo)間的大氣透射比。

由目標(biāo)單位反射面元反射到半球空間內(nèi)的回波功率PSP為

式中：ρ為目標(biāo)的反射比；θ為探測(cè)方向與朗伯面表面法線的夾角。離反射源距離為R 的空間點(diǎn)上的輻照度HR 由下式表示：

2 接收光學(xué)系統(tǒng)的數(shù)學(xué)建模

2.1 接收光學(xué)系統(tǒng)參數(shù)

接收光學(xué)系統(tǒng)子午面視場(chǎng)角α；

接收光學(xué)系統(tǒng)弧矢面視場(chǎng)角β；

接收口徑L×H。

2.2 接收光學(xué)系統(tǒng)仿真及其數(shù)學(xué)建模

本文所研究的是以二次曲柱面反射鏡為主體形式的光學(xué)系統(tǒng)，如圖7所示。

圖7 接收光學(xué)系統(tǒng)光線追跡圖

入射光線經(jīng)過(guò)接收光學(xué)系統(tǒng)聚焦成像在反射鏡的焦平面上，形成與反射鏡口徑等寬度的一段線像，圖8是成像光斑的能量分布圖。

為了減小探測(cè)器的動(dòng)態(tài)范圍，本光學(xué)系統(tǒng)采用的是異形探測(cè)器。由圖7可知，入射光線經(jīng)過(guò)光學(xué)系統(tǒng)在探測(cè)器上的成像位置與其在接收口徑上的入射位置有關(guān)，由于所需建模的接收光學(xué)系統(tǒng)采用的是不規(guī)則面形的探測(cè)器，這就使得光學(xué)設(shè)計(jì)軟件不能進(jìn)行精確的成像光斑與探測(cè)器相對(duì)位置關(guān)系的仿真?；谝陨显?，此接收光學(xué)系統(tǒng)采用Matlab軟件進(jìn)行建模。

圖8 接收光學(xué)系統(tǒng)成像光斑能量分布圖

將引信接收窗口在兩正交方向離散化，均分為多個(gè)單位接收面元，如圖9所示。

圖9 引信窗口接收面元與目標(biāo)反射面元的關(guān)系示意圖

引信窗口單位接收面元AW上接收到的回波功率PW由下式表示：

將公式（1）、（2）、（3）、（4）、（5）整理后得到

由圖7可以看出接收光學(xué)系統(tǒng)的入射光線是一個(gè)在某一時(shí)刻并攜帶了光束能量的空間矢量→S（PW，t），它是一個(gè)多維變量，將此矢量在接收光學(xué)系統(tǒng)的子午面、弧矢面投影，得到傳統(tǒng)定義的子午面入射角α、弧矢面入射角β。在進(jìn)行光學(xué)系統(tǒng)的建模時(shí)，可以分別在這兩個(gè)投影面上進(jìn)行光線追跡，得到成像光斑位置。

2.3 接收光學(xué)系統(tǒng)成像數(shù)學(xué)建模

假設(shè)在接收窗口點(diǎn)A 處，且子午面入射角為3°，弧矢面入射角為45°時(shí)，經(jīng)過(guò)仿真可以精確的判斷入射光線經(jīng)過(guò)光學(xué)系統(tǒng)后是否成像在探測(cè)器的光敏面上，如圖10所示。

如果成像光斑落在探測(cè)器的光敏面上，根據(jù)PIN 硅光電二極管的光電轉(zhuǎn)換特性［3］，可以得到

探測(cè)器的輸出電流信號(hào)及其大小

圖10 反射光束經(jīng)過(guò)光學(xué)系統(tǒng)后成像在探測(cè)器上的位置示意圖

式中：A 是探測(cè)器的響應(yīng)度，單位為A／W；A0是接收光學(xué)系統(tǒng)總效率。由于現(xiàn)階段引信的激光回波能量一般為幾個(gè)μw ～幾十個(gè)μw，完全在探測(cè)器的線性響應(yīng)范圍內(nèi)，所以就可以將某一時(shí)刻所有有效的成像光斑所產(chǎn)生的信號(hào)線性迭加，得

由上式可以得到在某一時(shí)刻t探測(cè)器的輸出電流信號(hào)It，將此信號(hào)在時(shí)間軸上展開(kāi)計(jì)算得到I（t），它就是探測(cè)器與一個(gè)發(fā)射脈沖相對(duì)應(yīng)的輸出電流隨時(shí)間變化的曲線。

3 結(jié)論

激光引信的光學(xué)系統(tǒng)形式以及所要探測(cè)的目標(biāo)特性都較為復(fù)雜，本文采用的這種統(tǒng)計(jì)建模與Matlab軟件編程相結(jié)合的光學(xué)系統(tǒng)虛擬樣機(jī)建模方法合理有效，從最終的仿真結(jié)果來(lái)看，比較符合物理樣機(jī)的現(xiàn)實(shí)情況，真實(shí)的反應(yīng)了光學(xué)系統(tǒng)的特性，滿足激光引信虛擬樣機(jī)的要求。

［1］ 張京國(guó)，周宗海，劉建新.近炸引信虛擬樣機(jī)技術(shù)研究［J］.制導(dǎo)與引信，2006，（3）.

［2］ 趙為黨，楊李茗，等.激光光束整形的設(shè)計(jì)和研究［J］.紅外與激光工程，1999，28（1）.

［3］ 安毓英，曾曉東.光電探測(cè)原理［M］.西安：西安電子科技大學(xué)出版社，2004.