球直坐標變換非對稱誤差分析

李慧敏 馬 可 范學禎

(西安電子工程研究所 西安 710100)

0 引言

對于常見的雷達、光電等主動探測設(shè)備,信號收發(fā)一般放置在一起,稱作單站測量,單個接收站可獨自獲取目標的距離、角度等信息。對于收發(fā)分置的探測設(shè)備來說,目標參數(shù)的獲取方式只要不依賴于多個接收站之間的配合,也可稱作單站測量設(shè)備。

對于單站測量設(shè)備來說,其原始測量值在球坐標系下獲取,各測量值的誤差一般是互相獨立的,而濾波、外推等算法則在直角坐標系下進行,由于球坐標與直角坐標之間的非線性映射關(guān)系[1-2],導致測量誤差互相耦合,空間分布產(chǎn)生畸變。

濾波的前提是隨機變量的期望值等于目標真實位置[3-5],但若期望偏離真實值,則濾波值必然偏離目標真實位置。對于比較常見的直線運動軌跡,如果測量誤差在直線兩側(cè)的分布不對稱,將導致目標運動軌跡濾波值產(chǎn)生單側(cè)偏移,產(chǎn)生系統(tǒng)偏差。這種系統(tǒng)偏差對于一些需要外推遭遇點的應(yīng)用場景來說[1,6],誤差將進一步放大,因此必須引起重視。

本文從球坐標系與直角坐標系的坐標變換出發(fā),對直角坐標系下誤差分布的非對稱性進行了分析計算。

1 球坐標系與直角坐標系的坐標變換

以兩坐標雷達為例,球坐標系到直角坐標系的坐標轉(zhuǎn)換公式見式(1)所示。

(1)

其中R、θ為球坐標系下的原始測量值,分別代表目標距離和角度,x、y為直角坐標系下的坐標值。坐標逆變換公式為

(2)

兩種坐標系之間的非線性映射關(guān)系如圖1所示。

圖1 直角坐標系與球坐標系之間的非線性映射

我們描述物體的運動軌跡是“直線”的,一般均默認其處于直角坐標系。直線軌跡轉(zhuǎn)換到球坐標系下后,其距離-角度曲線一般都不再是直線。

直線軌跡可分為兩種:截距為0和截距不為0。

1)截距為0時,x、y的關(guān)系式為y=cx,根據(jù)坐標逆變換公式可推導得球坐標系下的關(guān)系式為θ=tan-1c,可知此時角度θ為恒定值,距離-角度關(guān)系體現(xiàn)為直線。

2)截距不等于0時,x、y的關(guān)系式為y=cx+d,且d≠0,將式(1)帶入可得Rsinθ=cRcosθ+d,兩邊同時對R求導數(shù),整理后可得

(3)

(4)

圖2 直角坐標系與球坐標關(guān)系示意圖

2 球坐標系下的誤差分布

(3)

在球坐標系下,該函數(shù)關(guān)于點(R0,θ0)呈對稱分布,如圖3所示。

圖3中的圓環(huán)為概率分布等高線。在(R,θ)二維平面上,任意一條通過點(R0,θ0)的直線θ=a·R+b,均可將p(R,θ)等分,即滿足

(6)

圖3 二維聯(lián)合概率密度函數(shù)

根據(jù)第2節(jié)的結(jié)論:截距d=0時,θ=tan-1c。此時在(R,θ)坐標系下表現(xiàn)為一條直線,如圖4所示。

圖4 概率密度函數(shù)關(guān)于θ=tan-1c對稱分布

因此θ=tan-1c將概率密度函數(shù)等分為兩部分,誤差在直線兩側(cè)的分布概率是相等的。

根據(jù)第2節(jié)的結(jié)論:截距d≠0時,在(R,θ)坐標系下,R和θ的關(guān)系表現(xiàn)為一條經(jīng)過點(R0,θ0)的單調(diào)曲線,假定其函數(shù)表達式記為θ=g(R)。若c、d均大于0,曲線形狀如下圖所示。

圖5 概率密度函數(shù)關(guān)于θ=g(R)不對稱分布

由于θ=g(R)是單調(diào)的,因此必然可以找到一條直線,在(R0,θ0)處與θ=g(R)相切。由此可得出以下結(jié)論:

1)對于曲線θ=g(R)上的任意一點(Ri,θi),其切線將(Ri,θi)處的分布概率等分為兩部分。

2)對于曲線θ=g(R)上的任意一點(Ri,θi),其分布概率在曲線θ=g(R)兩側(cè)的分布概率不相等。

3)隨著R的增大,曲線θ=g(R)在點(Ri,θi)處的斜率越接近于0,從而和其切線越來越靠近,因此曲線兩側(cè)的分布概率逐漸趨近相等。

由此可知,在直角坐標系下進行濾波時,目標在真實軌跡兩側(cè)的分布概率是不相等的,因此導致期望偏向真實軌跡的一側(cè),引起系統(tǒng)誤差。當目標距離較近,且航路截徑和目標距離相當時,其影響將比較明顯;目標距離越遠,該誤差逐漸減小,其影響可忽略。

3 仿真計算

假定在直角坐標系下有一條直線軌跡,其坐標值x,y滿足關(guān)系y=0.1x+500,x的取值范圍為[100,500]。設(shè)定測距誤差σR=5m,測角誤差σθ=1°。

通過蒙特卡洛仿真,模擬某一特定測量點的誤差散布。例如,直角坐標系下,某一點處的散布誤差“點云”如圖6所示。

圖6 直角坐標系下的直線航跡與誤差分布

其中實線為目標的直線軌跡,“點云”為軌跡上某點處的測量值散布誤差。對應(yīng)的球坐標系下的散布誤差為圖7所示。

圖7 球坐標系下直線航跡與誤差分布

統(tǒng)計點云處于軌跡上方(參考縱軸正向)的點的個數(shù),可以反映出誤差的單邊分布概率。對于圖6中的情況,軌跡上方的單邊分布概率約為49.6%,沿軌跡上各點依次通過蒙特卡洛仿真手段,分析單邊分布概率,仿真結(jié)果如圖8所示。

圖8 單邊分布概率的最小二乘擬合

圖8中,細實線為通過蒙特卡洛仿真得出的單邊分布概率。對仿真概率進行最小二乘擬合,近似得出實際的單邊分布概率,如圖8中粗實線所示。

可見,對于直線軌跡上的各點,實際的單邊分布概率始終小于0.5,說明在設(shè)定的仿真條件下,軌跡上方的單邊分布概率始終小于軌跡下方的單邊分布概率,由于概率分布不對稱,屬于“偏心”分布,將導致期望值偏向軌道下方一側(cè)。根據(jù)圖8的仿真結(jié)果,可得出軌跡上每一點處目標的位置期望值,從而得到預(yù)測的目標軌跡。實際軌跡與預(yù)測軌跡對比如圖9所示。

圖9 直角坐標系下的實際軌跡與預(yù)測軌跡

預(yù)測軌跡和實際軌跡的差別較小,二者幾乎重疊在一起,將圖9左圖黑色框線內(nèi)部放大,如右圖所示。可見由于誤差的非對稱分布,預(yù)測軌跡始終處于實際軌跡下方。預(yù)測軌跡偏離實際軌跡的差值變化趨勢如圖10所示。

圖10 直角坐標系下預(yù)測軌跡與實際軌跡的偏差

目標距離越近,二者的角度系統(tǒng)偏差越大,最大接近0.03°,目標距離越遠,偏差越小。對于遠程單站測量設(shè)備來說,誤差的非對稱分布導致的角度系統(tǒng)偏差一般很小,因此在試驗數(shù)據(jù)中一般無法覺察到它的存在。對于近距探測設(shè)備來說,該誤差隨著目標距離的減小而迅速增大,需要予以考慮。

4 結(jié)束語

球坐標系下獲取的原始數(shù)據(jù),通過坐標變換到直角坐標系下進行處理時,誤差分布不再對稱,導致濾波的期望值為有偏估計,引起濾波系統(tǒng)誤差。以直線航路為例進行了仿真,結(jié)果表明該誤差隨距離的減小而增大,因此對遠程探測設(shè)備影響較小,對近距探測設(shè)備影響較大。

對于曲線航路來說,可分段等效為直線航路,因此其誤差分布在局部也是非對稱的,同樣,其非對稱性在近距更為明顯。