張振鐸，張 彬

（中國科學院 長春光學精密儀器與物理研究所，吉林 長春 130033）

引 言

靶場光測設備是對彈道導彈、飛行器、衛(wèi)星等武器設備進行精確觀測的有效工具，它通過角度測量和交誤差分析實施對空間目標精確定位以完成外彈道的觀測。測量中的各種誤差源直接決定定位精度［1，2］。但誤差并不是越小越好還要受到研制成本等多種因素的制約，同時各種誤差因素對光測系統(tǒng)而言，其重要性和影響方式有所不同，對各種誤差源進行精確分析對靶場光測設備有重要的實際意義。一方面在設計的初期研制過程中可以對設備測量不確定度有影響的各種因素進行嚴格的過程控制，同時對優(yōu)化材料選用、設備制造等多種環(huán)節(jié)進行控制，以更低的研制成本更好地滿足用戶的使用要求［3－5］；另一方面通過對光測設備進行合理的布站以提高光測設備測量精度［6，7］。光測設備的精度模型和光測設備的布站是兩個相互聯(lián)系的問題，對光測設備的布站優(yōu)化必須以精確的誤差模型為基礎，同時對光測設備的誤差分析也需要給出在一定布站方式下對特定目標軌跡的精度分析。

文中以光電經(jīng)緯儀為例，對光電經(jīng)緯儀建立了精確的Verilog－A模型，并在此基礎上使用 Monte Carlo法對光電經(jīng)緯儀的各種誤差源的影響進行了詳細的統(tǒng)計學分析和研究，并針對彈道導彈不同的布站方式進行數(shù)值分布，給出能適應不同靶場地形的通用的布站優(yōu)化方法。

1 光電經(jīng)緯儀的誤差模型

文中提出的誤差分析方法具有一定的通用性，使用通用的模擬系統(tǒng)建模語言Verilog－A以對光電經(jīng)緯儀建立數(shù)學模型，并在數(shù)學模型中對測量精度有影響的各種誤差源進行建模。Verilog－A提供了層次化的模擬系統(tǒng)模型構(gòu)架，可以一定的數(shù)學表達式在抽象的層次上對系統(tǒng)建模，以方便系統(tǒng)的分析設計和驗證。Verilog－A支持直流、交流、瞬態(tài)、Monte Carlo法等多種分析方法。提供給用戶一定的設計參數(shù)，用戶使用Verilog－A模型可以對設備的性能進行仿真分析進一步指導光測設備的應用。對各種光測設備建立通用的Verilog－A模型具有非常重要的實際意義。光電經(jīng)緯儀的Verilog－A模型分為基本模型和誤差模型，如圖1所示。

圖1 評估光電經(jīng)緯儀測量誤差的VerilogA模型框圖Fig.1 The verilog A model for theodolite measure error evaluation

其中基本模型實現(xiàn)光電經(jīng)緯儀的基本功能，即實現(xiàn)空間目標3維坐標（x，y，z）到俯仰角E和方位角A的轉(zhuǎn)換：

誤差模型使用坐標變換的方法，光電跟蹤經(jīng)過了基座坐標系、水平軸坐標系、照準部坐標系、相面坐標系，分別引入豎軸誤差、橫軸誤差、照準軸誤差和傳感器誤差，同時編碼器也會引入編碼器誤差，在坐標系X1（x，y，z）經(jīng)過旋轉(zhuǎn)變換轉(zhuǎn)到坐標系X2（x′，y′，z′），旋轉(zhuǎn)θ角的旋轉(zhuǎn)變換矩陣為M（θ），則同一點在X1中坐標與X2中坐標的關系為：

由于旋轉(zhuǎn)變換矩陣的逆矩陣為：

得到：

光電跟蹤依次經(jīng)過了基座坐標系、水平軸坐標系、照準部坐標系、相面坐標系等旋轉(zhuǎn)變換。各誤差源所引入的旋轉(zhuǎn)變換矩陣定義為Mi（i＝1，2…n），則在相面坐標系中單位光軸方向矢量（1，0，0）T在基座坐標系中的坐標可以表達為：考慮到方位角和俯仰角旋轉(zhuǎn)方向的定義，各坐標系間旋轉(zhuǎn)矩陣如表1所示

其中rot x，rot z，rot y分別為繞x，z，y軸的旋轉(zhuǎn)變換矩陣，在旋轉(zhuǎn)角很小時可以表示為：

表1 光電經(jīng)緯儀誤差模型中旋轉(zhuǎn)變換Tab.1 The rotation transformation used for theodolite intercross measure model

坐標變換法是一種對光測設備誤差源進行表達的通用的方法，使用此方法可以方便地對光電經(jīng)緯儀建立包含各種誤差源的準確的數(shù)學模型，并進一步建立光電經(jīng)緯儀雙站交匯測量模型如圖2所示，輸出dx，dz，dy分別為x，z，y三個方向上的交匯測量誤差。其中交匯測量模塊在異側(cè)布站時使用K公式，在同側(cè)布站時使用L公式，在Verilog－A模型中需要對除法運算分母做限幅處理以避免在平衡點處出現(xiàn)收斂問題。

圖2 光電經(jīng)緯儀交匯測量模塊Fig.2 Theodolite intercross measure model

2 基于Verilog－A模型的誤差分析法

針對上述光電經(jīng)緯儀的數(shù)學模型進行統(tǒng)計分析，基本模型參數(shù)如表2所示。使用兩種基本的誤差分析方法，即最壞情況法分析和Monte Carlo法分析。計算了各種單項誤差對系統(tǒng)性能的影響?，F(xiàn)使用的分析方法具有一定的通用性，可用于各種光測儀器的誤差計算中。

2.1 最壞情況法分析法

在最壞情況下假定所有的誤差源都取最大誤差，由于在實際情況下所有的誤差源不是同時取得最大誤差數(shù)值，因此最壞情況分析將得到較為悲觀的預測數(shù)值，但這種分析方法可以快速估計各單項誤差對系統(tǒng)誤差影響程度和影響方式，同時給出理論上最大誤差。使用此模型仿真了各單項誤差都取最大數(shù)值的最壞情況，誤差主要影響方位角A，而對俯仰角E影響較小。其中編碼器誤差直接加入到方位角和俯仰角誤差中，而傳感器誤差、照準差、橫軸差對系統(tǒng)誤差的影響相對較為復雜，可以通過掃描相關誤差參數(shù)計算出各誤差參數(shù)對系統(tǒng)性能的影響。

表2 光電經(jīng)緯儀模型參數(shù)選取Tab.2 Parameter selection for intercross measure model

使用該模型仿真了因俯仰角變化而引起的傳感器誤差、軸系照準差和水平軸誤差A分量的變化曲線，如圖3所示。從圖3中可以看出，隨著俯仰角的增加，三種誤差源快速增加，文獻表明照準差和水平軸誤差的δE分別與sec（E）和tan（E）成正比，而水平軸誤差與tan（E）與仿真結(jié)果相一致。同時可以看出在俯仰角較大時照準差對系統(tǒng)影響約為橫軸誤差的2倍。其中圖3（a）表示了不同誤差像素數(shù)目對系統(tǒng)誤差的影響，傳感器誤差將隨著誤差像素數(shù)目的提高而增加，誤差像素主要來自量化誤差、細分誤差和拖尾誤差。

圖3 俯仰角變化對傳感器誤差、軸系照準差和橫軸誤差的影響Fig.3 The influence of collimate error，vertical axes error，horizontal axes error

圖4表示了豎軸誤差隨方位角和俯仰角變化曲線，其中左圖表示豎軸方位誤差隨方位角和俯仰角的變化，可以看出豎軸差與方位角呈現(xiàn)正弦變化規(guī)律，且在豎軸傾斜方向上取得最大值，而隨俯仰角E以tan（E）的規(guī)律快速增加；右圖豎軸俯仰誤差與俯仰角E無關且隨方位角呈現(xiàn)正弦規(guī)律變化。此仿真結(jié)果與文獻［1］論述相一致。圖5表示總誤差源隨著方位角和俯仰角變化曲線，其中總方位誤差隨方位角A以正弦形式變化，其原因主要來自豎軸誤差的影響，總方位誤差隨俯仰角E快速增加，主要受到橫軸差、照準差以及傳感器誤差的影響。

2.2 Monte Carlo法分析法

在實際情況下各個單項是以一定概率分布的形式出現(xiàn)的，誤差合成不是簡單的線性疊加關系，傳統(tǒng)的方法是基于統(tǒng)計學的t分布和χ2分布不確定度分析法，這種方法在處理測量誤差傳遞時是基于線性化近似模型，同時假設各種誤差源間是相互獨立的。由于光測設備本身是一個復雜的非線性系統(tǒng)，同時誤差源間也不是完全獨立的，因此決定了這種方法具有一定的局限性。Monte Carlo法是一種通用的誤差分析工具，它將設備的各項誤差源表達為一定概率分布函數(shù)的形式，以相同概率分布產(chǎn)生隨機數(shù)進行仿真計算，從而得到各種誤差對系統(tǒng)的影響。系統(tǒng)的輸出呈現(xiàn)出一定統(tǒng)計分布，通過使用MatLab統(tǒng)計學工具箱擬合的方法，可以得到系統(tǒng)輸出的準確的統(tǒng)計分布函數(shù)，并以一定的數(shù)字特征如期望和方差的形式來表達。Monte Carlo算法的準確度主要取決于采樣點的數(shù)目，可以通過合理的選用采樣點的數(shù)目以達到所需的計算精度。

圖4 豎軸誤差隨方位角和俯仰角的變化曲線Fig.4 The relationship of vertical axes error with azimuth and pitch angle

圖5 最壞情況下經(jīng)緯儀總誤差隨方位角和俯仰角變化曲線Fig.5 The total error with azimuth and pitch angle in the worst case

圖6 Monte Carlo誤差分析的一般步驟Fig.6 The general step for Monte Carlo error analysis

圖7 Monte Carlo分析結(jié)果的正態(tài)分布擬合Fig.7 The normal distribution fitting result

圖6表示了對系統(tǒng)進行Monte Carlo法分析的一般步驟。在Monte Carlo法分析中假設各誤差源符合正態(tài)分布，精度范圍為±δ，其中δ為正態(tài)分布標準差，各誤差精度范圍如表1所示。圖7為方位角135°，俯仰角為45°時使用Monte Carlo法分析得到統(tǒng)計直方圖以及使用MatLab統(tǒng)計工具箱擬合的結(jié)果，擬合得到正態(tài)分布δ＝8.53。使用此方法進一步分析了各種誤差源對系統(tǒng)性能的影響如表3所示，是在考慮各種單項誤差的情況下計算得到的系統(tǒng)的方位角誤差和俯仰角誤差。從分析結(jié)果看，照準差和傳感器誤差對總誤差的影響最大。

3 基于Verilog－A模型的光電經(jīng)緯儀布站優(yōu)化

表3 A＝135°，E＝45°情況下 Monte Carlo法分析結(jié)果Tab.3 The Monte Carlo analytical result when A＝135°and E＝45°

基于Monte Carlo法的精度分析可以進一步應用于光電經(jīng)緯儀的布站設計中，可以針對特定的地形特點和彈道軌跡使用Monte Carlo法仿真光電經(jīng)緯儀的交匯測量過程，彈道軌跡公式如下：

現(xiàn)選擇彈道發(fā)射后2s到7s進行理論彈道軌跡優(yōu)化，彈道發(fā)射方向為X方向，彈道方程參數(shù)為v0＝4.5km／s，θ＝0.13，采用異側(cè)布站方式，兩站點坐標為（10km，20km，0）和（10km，－20km，0），使用K公式進行交匯測量，在只考慮經(jīng)緯儀的測量誤差的情況下，測量得到的x，z，y三個方向的誤差曲線隨時間變化曲線如圖8所示。對其在仿真時間內(nèi)的平均誤差繪制統(tǒng)計直方圖，并使用MatLab統(tǒng)計工具箱進行正態(tài)分布擬合得到如圖9的計算結(jié)果，x，z，y三個方向的正態(tài)分布方差分別為1.168m、2.756m和0.515m。計算表明在垂直于飛行方向的z方向上誤差最大，而在高度y方向下誤差最小。

圖8 仿真x，z，y方向測量誤差隨時間變化曲線Fig.8 The time variable measure error in the direction of x，z，y

圖9 x，z，y方向平均誤差直方圖Fig.9 The average measure error histogram in the direction of x，z，y

在異側(cè)布站情形下，設置站點z坐標分別為＋20km和－20km，使用Monte Carlo法仿真得到站點x坐標變化對交匯測量誤差的影響，如圖10所示，圖中4條曲線分別對應于x，z方向站點定位誤差d分別為0m，1m，2m和3m時系統(tǒng)三個方向上的測量誤差，從計算結(jié)果可知站點定位誤差在垂直于飛行方向（z方向）上影響較大，且在站點定位誤差小于1m時，經(jīng)緯儀本身誤差起主要作用；而當站點定位誤并差大于1m時，站點定位誤差將成為主要的誤差源。同時注意到當站點x坐標為20km時三個方向交匯測量誤差均呈現(xiàn)最小值，為最優(yōu)布站位置。仿真結(jié)果表明在飛行段x方向中部位置布站可以獲得最小測量誤差。

圖10 異側(cè)布站時不同站點位置精度下站點x方向距離與測量誤差關系Fig.10 The relationship between measure error xdistance for different side station position optimization

設置站點x坐標位＋20km，－20km時，使用Monte Carlo法仿真了站點z坐標變化對交匯測量誤差的影響，如圖11所示，其中z坐標變化范圍設置為1～20km，從仿真結(jié)果可以看出：

（1）x，z方向的誤差隨z方向站點距離減小而減小。

（2）而y方向誤差隨z方向站點距離呈現(xiàn)拋物線形態(tài)，且受布站位置誤差的影響較大。在不存在站點定位誤差時，站點z坐標為2.25km處，y方向誤差最??；在布站定位誤差為1m時，站點z坐標為5m處，y方向誤差最?。徽军c定位誤差為3m時，站點z坐標為5m處，y方向誤差最小。定義y方向誤差最小的點為最優(yōu)布站點。

（3）最優(yōu)布站點z方向距離定義為z＊，當站點z向距離z＜z＊時，y方向誤差隨z減少快速增加；當站點z向距離z＞z＊時，x方向和z方向誤差隨z增加快速增加。

圖11 異側(cè)布站時不同站點位置精度下站點z方向距離與測量誤差關系Fig.11 The relationship between measure error z distance for different side station position optimization

文獻［1］表明對定點測量時，當站點交匯角在90°左右交匯測量的精度最高，使用Monte Carlo法計算表明針對特定的彈道軌跡，在考慮經(jīng)緯儀的各項誤差源以及布站位置誤差時，交匯測量精度將呈現(xiàn)更為復雜的關系，Monte Carlo法可以為光電經(jīng)緯儀的布站優(yōu)化提供指導。

使用Monte Carlo法進一步繪制了站點位置（x，z）變化時異側(cè)對稱布站情形下三維誤差分布曲面如圖12所示，使用此方法可以快速計算得到最優(yōu)布站位置為（5km，20km，0）附近。如果給定一定的靶場地形信息，使用此方法可以針對特定靶場地形給出最優(yōu)布站方式。

圖12 站點位置精度1m時，x，z，y方向測量誤差分布Fig.12 The x，z，y measure error distribution with 1mthe station position precision

4 結(jié) 論

對光電經(jīng)緯儀建立了包含各種誤差源的準確的數(shù)學模型，使用此模型進行最壞情況分析和Monte Carlo法分析，分析了各種誤差源影響系統(tǒng)性能的程度和方式。進一步使用此模型對雙站異側(cè)布站情形進行了優(yōu)化，分析表明針對沿x方向發(fā)射的彈道軌跡，飛行段x方向中部位置布站可以獲得最小測量誤差；布站z方向存在最優(yōu)位置z＊，此處y方向交匯測量誤差最小，當布站z向距離z＜z＊時，y方向誤差隨z減少快速增加；當布站z向距離z＞z＊時，x方向和z方向誤差隨z增加快速增加。提出的基于Monte Carlo法的布站優(yōu)化可以進一步推廣到特定靶場地形的情況，對經(jīng)緯儀總體設計及布站方式的選擇具有一定的理論指導意義。

［1］趙學顏，李迎春.靶場光學測量［D］.北京：裝備指揮技術學院，2001：35－38.

［2］王家騏.光學儀器總體設計［D］.長春：長春光機所研究生部，2003：27－32.

［3］杜俊峰.光電經(jīng)緯儀測量精度指標的確定［J］.應用光學，2006，5（6）：506－509.

［4］張 寧，沈湘衡.應用跟蹤誤差等效模型評價光電經(jīng)緯儀跟蹤性能［J］.光學精密工程，2003，18（3）：677－684.

［5］李 慧，沈湘衡.光電經(jīng)緯儀數(shù)字化模型的建立方法及其應用［J］.光電子技術，2007，27（2）：97－100.

［6］侯宏錄，周德云.光電經(jīng)緯儀異面交會測量及組網(wǎng)布站優(yōu)化設計［J］.光子學報，2008，37（5）：1023－1028.

［7］郭麗華，朱元昌.基于遺傳算法的光電經(jīng)緯儀布站優(yōu)化設計［J］.儀器儀表學報，2010，5（4）：741－746.