旋轉(zhuǎn)折射光楔對控制模型在機(jī)載光電探測設(shè)備中的應(yīng)用分析

時(shí)光煜，李明月，孫輝輝

(中航工業(yè)洛陽電光設(shè)備研究所，河南洛陽 471009)

旋轉(zhuǎn)折射光楔對控制模型在機(jī)載光電探測設(shè)備中的應(yīng)用分析

時(shí)光煜，李明月，孫輝輝

(中航工業(yè)洛陽電光設(shè)備研究所，河南洛陽 471009)

在現(xiàn)代機(jī)載光電探測設(shè)備中，旋轉(zhuǎn)折射光楔對起著擴(kuò)大探測視場和搜索范圍的重要作用。首先基于幾何光學(xué)和光的標(biāo)量衍射理論對旋轉(zhuǎn)折射光楔對的折光特性進(jìn)行了數(shù)學(xué)分析;然后結(jié)合其在機(jī)載光電探測設(shè)備中的功能要求，研究了一種旋轉(zhuǎn)光楔對控制模型，并基于該控制模型設(shè)計(jì)了旋轉(zhuǎn)折射光楔對的控制算法;最后，通過仿真實(shí)驗(yàn)證明了該控制算法的正確性和可靠性。

機(jī)械光電探測設(shè)備;旋轉(zhuǎn)折射光楔對;控制模型;控制算法

0 引言

機(jī)載光電探測設(shè)備主要用于對空、對地探測，將飛機(jī)航線前方或左右景物形成圖像供平顯或多功能顯示器顯示，以增強(qiáng)飛行員對空地態(tài)勢的感知能力。為了擴(kuò)大觀察范圍和實(shí)現(xiàn)對目標(biāo)的快速搜索定位，要求機(jī)載光電探測設(shè)備具備方位掃描功能和在“上”、“下”、“左”、“右”、“中”5個(gè)探測方位之間實(shí)現(xiàn)快速切換的功能。旋轉(zhuǎn)折射光楔對(下文簡稱“光楔對”)是較常使用的光學(xué)元件，通過控制光楔對的旋轉(zhuǎn)，在光路中能起到偏折光線的作用。由于其具有控制靈活、機(jī)械運(yùn)動平穩(wěn)、振動噪聲小、掃描效率高，結(jié)構(gòu)緊湊等優(yōu)點(diǎn)，光楔對已被廣泛應(yīng)用于性能先進(jìn)的機(jī)載光電產(chǎn)品中［1］。

1 光楔理論概述

光楔就是頂角很小的折射棱鏡，在光路中能起偏折光線的作用，如圖1所示。

圖1 光楔對光線的偏折作用Fig．1The refractive property of light for dual-wedge prism system

光楔的工作面是兩個(gè)折射面，兩折射面的交線稱為折射棱，兩折射面間的二面角稱為折射棱鏡折射角，用α表示。垂直于折射棱的平面稱為主截面。由于沿折射棱方向的分量對光線的折射沒有貢獻(xiàn)，光線只在平行于主截面的方向發(fā)生折射［2］，為簡化討論，本文只討論光線在主截面內(nèi)的折射情況。

在圖1中，光線AB沿主截面方向入射到折射棱鏡上，經(jīng)兩折射面的折射，出射光線DE與入射光線AB的夾角δ稱為偏向角。設(shè)棱鏡折射率為n，光線在兩折射面上的入射角和折射角分別為I1、I1'和I2、I2'，可得入射角I1和偏向角δ之間的關(guān)系式［3－4］為

當(dāng)I1為有限小時(shí)，由于α通常很小，可將式(1)簡化為

兩個(gè)相同的光楔一起使用，稱為光楔對。將其繞同一軸相對轉(zhuǎn)動，即一個(gè)光楔逆時(shí)針旋轉(zhuǎn)φ角，另一個(gè)光楔同時(shí)順時(shí)針旋轉(zhuǎn)φ角，折射光束就會在一定的范圍內(nèi)變動，且有如下關(guān)系式［5－6］成立將式(3)整理后得:

式(4)描述了光楔對產(chǎn)生的總偏向角δ與轉(zhuǎn)角φ之間的函數(shù)關(guān)系。

2 光楔對控制模型

在機(jī)載光電探測設(shè)備中，根據(jù)光楔對能對入射光線產(chǎn)生偏折作用的原理，通過控制光楔對的相對旋轉(zhuǎn)，可以實(shí)現(xiàn)對于光楔導(dǎo)向范圍內(nèi)任意方向的探測功能。

在旋轉(zhuǎn)折射光楔對控制模型中，輸入量為機(jī)載光電設(shè)備的探測方位矢量(u，ξ)，其中u，ξ分別表示方位、俯仰角度;輸出量是光楔對的旋轉(zhuǎn)角度(φ0，φ1)，其中φ0，φ1表示光楔對的旋轉(zhuǎn)角度。

2．1 坐標(biāo)系的建立

建立如圖2所示的“下左前”坐標(biāo)系，即X軸方向垂直向下，Y軸方向水平向左，Z軸方向水平向前。設(shè)圖2中光楔對的主截面位于ZOX平面中，折射棱鏡折射角為α，折射棱為n，光楔對當(dāng)前位置為初始零位。光線L2位于機(jī)載光電探測設(shè)備的外部，自右向左入射光楔對，出射光線L0平行于OZ軸向。

在圖2所示的坐標(biāo)系中，光線矢量可以有兩種表示方法:1)使用光線的方位角和俯仰角(u，ξ)來表示，優(yōu)點(diǎn)是物理意義清晰明確，便于終端用戶接受、操作; 2)使用光線與OZ軸向的夾角及其在XOY平面投影與OX軸夾角(δ，β)來表示，優(yōu)點(diǎn)是便于建立數(shù)學(xué)模型，簡化數(shù)值運(yùn)算。其中，δ為非負(fù)數(shù)，β從Z軸箭頭方向看進(jìn)去逆時(shí)針為正。

圖2 旋轉(zhuǎn)折射光楔對控制模型中坐標(biāo)系的建立Fig．2The coordinate system of control model of dual-wedge prism

2．2 光楔對的旋轉(zhuǎn)角度計(jì)算

由于在建模過程中對于矢量光線的表示方法不同(控制模型輸入?yún)?shù)采用第1)種表示方法，而光楔對轉(zhuǎn)角的計(jì)算中使用第2)種表示方法)，因此，需要進(jìn)行光線的矢量變換，以統(tǒng)一表示形式，方便數(shù)值運(yùn)算。由圖2可知有如下函數(shù)關(guān)系式成立。

將式(5)整理后得

式(6)、式(7)描述了光線矢量從(u，ξ)到(δ，β)的轉(zhuǎn)換關(guān)系。

結(jié)合第1節(jié)的光楔理論知識，對于給定的輸入?yún)⒘?δ，β)，可將待求解的光楔對的轉(zhuǎn)動角度(φ0，φ1)重寫為

其中:φ可由式(4)求得，不妨取正值，表示光楔對的相對轉(zhuǎn)角部分;β可由式(7)得到，表示光楔對的相同轉(zhuǎn)角部分。將φ、β代入式(8)可得:式(9)、式(10)描述了從光楔對控制模型中的輸入(u，ξ)到輸出(φ0，φ1)之間的函數(shù)映射關(guān)系。

2．3 約束方程

光楔對的光束導(dǎo)向范圍是一個(gè)關(guān)于旋轉(zhuǎn)軸對稱的圓錐立體角，即折射光線只能位于這個(gè)圓錐立體角之內(nèi)，可用最大導(dǎo)向角δmax表示光楔對對于光線的導(dǎo)向能力。δmax的大小取決于光楔的折射角α和材料的折射率n，可表示為

3 光楔對控制算法

為了能夠?qū)︼w機(jī)的前方進(jìn)行大范圍搜索，擴(kuò)大探測范圍，要求機(jī)載光電探測設(shè)備具備掃描功能和對特定方位的快速觀察功能，將詳細(xì)討論實(shí)現(xiàn)上述功能的光楔對控制算法。

3．1 方位掃描控制算法

方位掃描控制算法首先根據(jù)當(dāng)前機(jī)載光電設(shè)備的探測俯仰角度ξ和光楔對的導(dǎo)向范圍，確定掃描的水平軌跡;然后將掃描軌跡離散化，取出N個(gè)等步長的位置點(diǎn)，并計(jì)算出每個(gè)點(diǎn)處的方位矢量(u，ξ);緊接著根據(jù)式(6)、式(7)分別對這N個(gè)點(diǎn)進(jìn)行矢量變換，即實(shí)現(xiàn)從(u，ξ)到(δ，β)的變換;接下來將(δ，β)代入到式(9)、式(10)中，求解出光楔對的轉(zhuǎn)動角度(φ0，φ1);最后發(fā)送指令給光楔對控制電機(jī)，控制光楔對旋轉(zhuǎn)到指定角度(φ0，φ1)。光楔對的方位掃描控制算法流程圖參見圖3。

圖3 方位掃描控制算法流程圖Fig．3The flow chart of azimuth scanning algorithm

3．2 快速觀察控制算法

快速觀察算法就是在如圖4所示的“上”、“下”、“左”、“右”、“中”5個(gè)探測方位之間實(shí)現(xiàn)快速切換，滿足光楔對的旋轉(zhuǎn)角度盡量小的要求。另外，還希望兩個(gè)光楔的旋轉(zhuǎn)方向能夠盡量相反，以改善光楔對旋轉(zhuǎn)過程中視頻畫面的視覺效果。本文采用人工智能理論中狀態(tài)空間搜索的方法，從初始狀態(tài)節(jié)點(diǎn)出發(fā)，根據(jù)上述的規(guī)則約束，在狀態(tài)空間中搜索出一條解路徑，從而得出光楔對的最佳轉(zhuǎn)動角度［7－8］。

圖4 快速觀察示意圖Fig．4The sketch of quick watching

設(shè)光楔對的當(dāng)前位置為P0(φ0，φ1)，目標(biāo)位置為Pd(φd0，φd1)。光楔i(i=0，1)的當(dāng)前位置φi可以重寫為

其中:ri為整數(shù);θi∈［0，2π)。

為了滿足光楔轉(zhuǎn)動角度盡量小的要求，rdi取值空間為{ri－1，ri，ri+1}，相應(yīng)φdi取值空間記為(φ1di，φ2di，φ3di)。這樣，目標(biāo)位置的狀態(tài)空間便有9種不同的排列組合，分別記為Pijd(φid0，φjd1)(i，j=1，2，3)。

如果進(jìn)行全狀態(tài)搜索，至少需要進(jìn)行9次運(yùn)算，運(yùn)算量大，相應(yīng)的運(yùn)算時(shí)間較長。本文采用搜索深度為2級的有代價(jià)寬度優(yōu)先搜索方法［9－10］，如圖5所示，僅需6次運(yùn)算就可以找到最優(yōu)解。

圖5 狀態(tài)空間搜索樹Fig．5A state space searching tree

在圖5中順次將每個(gè)節(jié)點(diǎn)編號，其中節(jié)點(diǎn)1表示光楔對初始角度(φ0，φ1)。第一級求解光楔0的轉(zhuǎn)動角度，定義式(9)為其代價(jià)函數(shù);第二級求解光楔1的轉(zhuǎn)動角度，定義式(10)為其代價(jià)函數(shù)。建立Open表，用于記錄搜索過程中尚未考查的節(jié)點(diǎn)和新生成的節(jié)點(diǎn)，Open表的節(jié)點(diǎn)域可定義如下，示意圖見圖6。

圖6 Open表節(jié)點(diǎn)域定義示意Fig．6Node threshold of OPEN table

其中:j表示當(dāng)前節(jié)點(diǎn)編號;i表示父節(jié)點(diǎn)編號;gj表示節(jié)點(diǎn)j的代價(jià)函數(shù)值;Sj表示節(jié)點(diǎn)狀態(tài)，取值空間為(0，1)，1表示該節(jié)點(diǎn)被選用，0表示未被選用。

狀態(tài)空間搜索算法如下:

1)初始化Open表，把初始節(jié)點(diǎn)(編號1)放入Open表中，且令g1=φ0，S1=1;

2)將Open表的第一個(gè)節(jié)點(diǎn)移出，并依據(jù)圖5所示的搜索樹擴(kuò)展子節(jié)點(diǎn)2、3、4，依次放入Open表中;

3)依次將節(jié)點(diǎn)2、3、4從Open表中移出，計(jì)算函數(shù)值gj(j=2，3，4)，選取|g1－gj|最小的節(jié)點(diǎn)將其狀態(tài)標(biāo)識Sj置1，設(shè)該節(jié)點(diǎn)為m，其余節(jié)點(diǎn)的Sj清零;

4)將節(jié)點(diǎn)m擴(kuò)展子節(jié)點(diǎn)，并將新擴(kuò)展的子節(jié)點(diǎn)依次放入Open表中;

5)依次從Open表中移出節(jié)點(diǎn)，計(jì)算函數(shù)值gj，選取|gm－gj|值最小的兩個(gè)進(jìn)行比較，如果差別小于某特定值(例如小于10°)，則選取滿足(gj－φ1) (gm·g1)＜0的節(jié)點(diǎn)，否則選取|gm－gj|值最小的節(jié)點(diǎn)，設(shè)選中的節(jié)點(diǎn)為n，令Sn=1，其余節(jié)點(diǎn)的Sj清零;

6)將狀態(tài)標(biāo)識Sj=1(j=1，2，…，13)的節(jié)點(diǎn)從搜索樹中順次找出，即成功找到了解路徑。

4 數(shù)值仿真及其分析

仿真實(shí)驗(yàn)中選取光楔材料的折射率n=3．42，光楔頂角α=5．2°。由約束方程(13)可得最大導(dǎo)向角δmax=25．1680°。

首先，對光楔對的方位掃描控制算法進(jìn)行Matlab數(shù)值仿真，分別選取俯仰角ξ為±11．46°(約合±0．2 rad)，方位角在導(dǎo)向范圍內(nèi)往返掃描，掃描步長len設(shè)為2°(約合0．0349 rad)，掃描起始方位角位置為u= 0、方向?yàn)樨?fù)。仿真結(jié)果參見圖7、圖8，圖中星號表示光楔0的轉(zhuǎn)角φ0，菱形表示光楔1的轉(zhuǎn)角φ1，實(shí)線表示方位角度u。實(shí)驗(yàn)結(jié)果證明光楔運(yùn)動平滑連續(xù)，控制效果良好。

圖7 光楔對掃描仿真圖(俯仰角為11．46°)Fig．7The simulation diagram of dual-wedge prism system with pitch angle as 11．46°

圖8 光楔對掃描仿真圖(俯仰角為－11．46°)Fig．8The simulation diagram of dual-wedge prism system with pitch angle as－11．46°

然后，將光楔對快速觀察控制算法進(jìn)行仿真。在仿真實(shí)驗(yàn)中，設(shè)“上”、“下”、“左”、“右”、“中”5個(gè)快速觀察位置的坐標(biāo)(u，ξ)分別為(0，－18°)、(0，18°)、(21°，0)、(－21°，0)、(0，0)，快速觀察位置的切換次序?yàn)椤爸小稀隆蟆摇小?，使用光楔對的快速觀察控制算法，在狀態(tài)空間中搜索出相應(yīng)的解路徑，仿真結(jié)果參見表1。

表1 快速觀察控制算法得到的解路徑Table 1The result of the quick watching algorithm

由表1所示的仿真結(jié)果可知，在各種快速觀察位置的切換過程中，控制效果符合預(yù)期的要求。

5 結(jié)論

文章首先基于幾何光學(xué)和光的標(biāo)量衍射理論對旋轉(zhuǎn)折射光楔對的折光特性進(jìn)行了數(shù)學(xué)分析，然后結(jié)合其在機(jī)載光電探測設(shè)備中的功能要求，研究了一種旋轉(zhuǎn)光楔對控制模型，并基于該控制模型設(shè)計(jì)了旋轉(zhuǎn)折射光楔對的控制算法。最后，通過仿真實(shí)驗(yàn)證明了該控制算法是正確、有效的。

［1］郁道銀，談恒英．工程光學(xué)［M］．北京:機(jī)械工業(yè)出版社，2006．

［2］賀磊，袁家虎，李展，等．雙光楔高精度角度發(fā)生器設(shè)計(jì)［J］．光電工程，2002，29(6):46-49．

［3］玻恩M，沃耳夫E．光學(xué)原理［M］．北京:科學(xué)出版社，1978．

［4］徐斌．雙光楔實(shí)現(xiàn)外腔半導(dǎo)體激光器波長精密調(diào)諧研究［J］．激光與紅外，2008，38(5):444-446．

［5］韋中超，熊言威，莫瑋，等．旋轉(zhuǎn)雙光楔折射特性與二維掃描軌跡的分析［J］．應(yīng)用光學(xué)，2009，30(6):939-943．

［6］楊輝．光測設(shè)備引導(dǎo)數(shù)據(jù)的誤差分析［J］．電光與控制，2010，17(7):83-86．

［7］王萬森．人工智能原理及其應(yīng)用［M］．北京:電子工業(yè)出版社，2007．

［8］冷勁夫．單源路徑搜索算法的實(shí)現(xiàn)及比較［J］．地理空間信息，2010(3):154-156．

［9］馬蓉平．搜索推理技術(shù)與寬度優(yōu)先搜索的實(shí)現(xiàn)［J］．赤峰學(xué)院學(xué)報(bào):自然科學(xué)版，2009，25(3):42-43．

［10］陸釗．基于狀態(tài)空間表示的邏輯問題及其寬度優(yōu)先搜索算法實(shí)現(xiàn)［J］．中國科教創(chuàng)新導(dǎo)刊，2007(467): 82-83．

Application of a Rotary Dual-Wedge Prism System Control Model in Airborne Photoelectronic Equipment

SHI Guangyu，LI Mingyue，SUN Huihui
(Luoyang Institute of Electro-Optical Equipment，AVIC，Luoyang 471009，China)

In airborne photoelectronic detection equipment，rotary dual-wedge prism system plays an important role for enlarging the detection scope and searching range．First，the refractive property of the dual-wedge prism system was analyzed according to the theories of geometrical optics and scalar diffraction．Then，a control model was established for the dual-wedge prism system considering the functional requirements of it in airborne photoelectronic equipment，and the control algorithms were given based on the model．The correctness and reliability of the control model and control algorithms were proved through simulation．

airborne photoelectric detection equipment;rotary dual-wedge prism;control model;control algorithm

V271．4

A

1671－637X(2011)12－0052－04

2010－12－03

2011－03－02

時(shí)光煜(1984—)，男，河南輝縣人，碩士生，研究方向?yàn)榭刂评碚撆c控制工程。