多傳感器空間配準(zhǔn)的病態(tài)性及其穩(wěn)健估計(jì)

潘江懷， 喬慧

(江蘇自動(dòng)化研究所， 江蘇 連云港 222006)

多傳感器空間配準(zhǔn)的病態(tài)性及其穩(wěn)健估計(jì)

潘江懷， 喬慧

(江蘇自動(dòng)化研究所， 江蘇 連云港 222006)

為提高多傳感器空間配準(zhǔn)的精度和穩(wěn)定性，以最小二乘空間配準(zhǔn)算法為基礎(chǔ)，分析多傳感器空間配準(zhǔn)病態(tài)性產(chǎn)生的原因，指出系統(tǒng)誤差估計(jì)Fisher信息矩陣的病態(tài)性是影響空間配準(zhǔn)質(zhì)量的主要因素。運(yùn)用條件指標(biāo)對(duì)病態(tài)性程度進(jìn)行評(píng)價(jià)，并根據(jù)信息矩陣的奇異值呈階梯型分布特點(diǎn)，提出一種混合的奇異值修正穩(wěn)健估計(jì)算法。該算法通過(guò)只對(duì)較小的條件指標(biāo)進(jìn)行修正，實(shí)現(xiàn)對(duì)系統(tǒng)誤差觀測(cè)方程中不確定成分的有效抑制和確定成分的保留，能較好解決多傳感器空間配準(zhǔn)系統(tǒng)誤差實(shí)時(shí)估計(jì)問(wèn)題；通過(guò)兩種不同場(chǎng)景下的系統(tǒng)誤差估計(jì)結(jié)果證明，該算法的使用可顯著地提高系統(tǒng)誤差參數(shù)估計(jì)的精度和穩(wěn)定性。

信息處理技術(shù)； 空間配準(zhǔn)； 病態(tài)性； 穩(wěn)健估計(jì)； 奇異值分解； 多傳感器

0 引言

在多傳感器融合跟蹤系統(tǒng)中，傳感器探測(cè)的系統(tǒng)誤差是影響融合精度的最主要因素[1]，系統(tǒng)誤差的存在將直接對(duì)數(shù)據(jù)融合中的跟蹤門(mén)計(jì)算、目標(biāo)航跡狀態(tài)估計(jì)和協(xié)方差計(jì)算、航跡相關(guān)波門(mén)及準(zhǔn)則產(chǎn)生不利影響[2]。若直接將各傳感器的探測(cè)數(shù)據(jù)進(jìn)行融合，則將由于系統(tǒng)誤差的存在而使跟蹤結(jié)果惡化，導(dǎo)致目標(biāo)一批變多批或航跡呈鋸齒、航跡斷航，甚至不如單傳感器的跟蹤效果[3]?？臻g配準(zhǔn)是指消除探測(cè)源/探測(cè)系統(tǒng)所存在的系統(tǒng)誤差或?qū)⑵淇刂圃谌菰S范圍內(nèi)的處理過(guò)程。在多平臺(tái)多傳感器系統(tǒng)中，一般是指將分布在不同地理位置的傳感器目標(biāo)探測(cè)信息轉(zhuǎn)換到一個(gè)統(tǒng)一坐標(biāo)系或融合中心過(guò)程中，使由于傳感器探測(cè)、平臺(tái)定位不準(zhǔn)以及坐標(biāo)轉(zhuǎn)換等原因產(chǎn)生的系統(tǒng)誤差(配準(zhǔn)誤差)最小化[4]，以實(shí)現(xiàn)同一目標(biāo)位置信息在不同傳感器中的迭合和正確交接，確保后續(xù)數(shù)據(jù)融合正確進(jìn)行。因此，在多傳感器融合跟蹤系統(tǒng)中，必須對(duì)傳感器進(jìn)行空間配準(zhǔn)[5]。

在多傳感器系統(tǒng)誤差估計(jì)過(guò)程中，影響偏差可觀測(cè)性的因素有很多，偏差的種類、偏差之間相互耦合的程度、目標(biāo)與傳感器之間的幾何關(guān)系等多種因素都會(huì)造成偏差估計(jì)可觀測(cè)性的變化，從而導(dǎo)致偏差估計(jì)可觀測(cè)度難以進(jìn)行定量的描述[6]。偏差估計(jì)精度除了受到估計(jì)算法的影響，主要是由系統(tǒng)可觀測(cè)度決定，因此如何在未知可觀測(cè)度條件下定量地確定系統(tǒng)的配準(zhǔn)效果就成為一個(gè)難題[7]?，F(xiàn)有的空間配準(zhǔn)算法如歐洲和北美洲廣泛使用SAGE系統(tǒng)監(jiān)視系統(tǒng)中的實(shí)時(shí)質(zhì)量控制(RTQC)法[4]、考慮傳感器測(cè)量噪聲的最小二乘(LS)配準(zhǔn)法[5]及廣義LS(GLS)配準(zhǔn)法[6]、精確極大似然(EML)配準(zhǔn)法[7]、極大似然(MLR)配準(zhǔn)法[8]及文獻(xiàn)[9-14]提出的空間配準(zhǔn)算法，在可觀測(cè)度較高、目標(biāo)狀態(tài)觀測(cè)完整、傳感器布站合理且累計(jì)了一定的觀測(cè)數(shù)據(jù)條件下，通常都可以獲得比較好的配準(zhǔn)結(jié)果。當(dāng)可觀測(cè)度較低時(shí)，偏差的估計(jì)值將發(fā)生巨大的波動(dòng)[15-16]。這是因?yàn)橛^測(cè)數(shù)據(jù)中緩變的觀測(cè)信息被劇烈變化的隨機(jī)噪聲湮沒(méi)，觀測(cè)數(shù)據(jù)的變化主要由于隨機(jī)誤差引起，導(dǎo)致了估計(jì)結(jié)果的波動(dòng)，若直接根據(jù)這種估計(jì)值去修正目標(biāo)航跡，則后續(xù)的融合結(jié)果將必然會(huì)受到質(zhì)疑[17]。

本文基于LS空間配準(zhǔn)算法，分析了LS系統(tǒng)誤差估計(jì)的病態(tài)性產(chǎn)生原因，并根據(jù)系統(tǒng)誤差估計(jì)的Fisher信息矩陣奇異值分布為階梯型分布特點(diǎn)，提出一種新的奇異值修正(CSVD)算法，對(duì)奇異特征值進(jìn)行修正。該算法可兼顧系統(tǒng)誤差參數(shù)估計(jì)的分辨率與方差，實(shí)現(xiàn)了對(duì)系統(tǒng)誤差觀測(cè)方程中不確定成分的有效抑制和確定成分的有效保留，能有效解決多傳感器空間配準(zhǔn)系統(tǒng)誤差實(shí)時(shí)估計(jì)問(wèn)題，并從兩種不同的仿真場(chǎng)景進(jìn)行了驗(yàn)證比較。結(jié)果表明：當(dāng)信息矩陣的條件數(shù)明顯分梯形分布時(shí)，該方法可顯著提高系統(tǒng)誤差估計(jì)結(jié)果的精度和穩(wěn)定性。

1 多傳感器系統(tǒng)誤差的LS估計(jì)

對(duì)于不同平臺(tái)的n部三坐標(biāo)雷達(dá)而言，需要估計(jì)的系統(tǒng)誤差參數(shù)有3n個(gè)之多，(Δri,Δβi,Δεi)分別表示為第i個(gè)傳感器的距離系統(tǒng)誤差、方位角系統(tǒng)誤差、俯仰角系統(tǒng)誤差，i=1，…，n. 為簡(jiǎn)化推導(dǎo)，下面以兩部傳感器為例對(duì)空間配準(zhǔn)的算法進(jìn)行分析。如圖1所示，假設(shè)兩傳感器的大地坐標(biāo)為(L1,B1,H1)和(L2,B2,H2)，兩傳感器在地心地固(ECEF)坐標(biāo)系中的位置為(xp1,yp1,zp1)、(xp2,yp2,zp2)，k時(shí)刻兩傳感器對(duì)目標(biāo)距離、方位角和俯仰角的探測(cè)為(r1(k),β1(k),ε1(k))、(r2(k),β2(k),ε2(k))，轉(zhuǎn)換到兩傳感器地理坐標(biāo)系下的位置為(x1(k),y1(k),z1(k))、(x2(k),y2(k),z2(k)). 以下省略時(shí)刻k.

圖1 有系統(tǒng)誤差下的多平臺(tái)傳感器探測(cè)Fig.1 Multi-sensor detection with bias

若沒(méi)有任何誤差，兩傳感器同一時(shí)刻對(duì)同一目標(biāo)的探測(cè)在ECEF坐標(biāo)系中應(yīng)該迭合，即

(1)

式中：TR1和TR2分別為兩傳感器地理坐標(biāo)系到ECEF坐標(biāo)系的旋轉(zhuǎn)矩陣；(xi,t,yi,t,zi,t)為目標(biāo)在傳感器i(i=1，2)地理坐標(biāo)系下的真值：

(2)

(3)

兩傳感器對(duì)同一目標(biāo)的探測(cè)在ECEF坐標(biāo)系中的坐標(biāo)之差如下：

(4)

簡(jiǎn)記為

Z=Ab，

(5)

式中：A=[TR1Jf1,-TR2Jf2]；b=[Δr1,Δβ1,Δε1,Δr2,Δβ2,Δε2]T.

當(dāng)有k(k>2)個(gè)量測(cè)值時(shí)，(5)式的LS解為

=(ATA)-1ATZ.

(6)

定義系統(tǒng)誤差參數(shù)b的Fisher信息矩陣為

J=ATA.

(7)

2 LS系統(tǒng)誤差估計(jì)的病態(tài)性

LS估計(jì)系統(tǒng)誤差的條件是估計(jì)參數(shù)b的信息矩陣J可逆，即A列滿秩。當(dāng)A的列線性相關(guān)或弱相關(guān)時(shí)，信息矩陣是病態(tài)的。

2.1 常見(jiàn)的病態(tài)場(chǎng)景

1) 傳感器平臺(tái)相距較近。當(dāng)兩傳感器相距較近時(shí)，兩傳感器對(duì)同一目標(biāo)的探測(cè)距離、方位角、俯仰角近似相等，即有(r1(k),β1(k),ε1(k))≈(r2(k),β2(k),ε2(k))，而 (7) 式中矩陣A的系數(shù)完全由傳感器對(duì)目標(biāo)探測(cè)的距離、方位角、俯仰角所決定，將導(dǎo)致矩陣A的列向量具有相關(guān)性，從而使信息矩陣J是病態(tài)的。

2) 觀測(cè)數(shù)據(jù)密集情況。當(dāng)目標(biāo)相對(duì)于傳感器平臺(tái)靜止不動(dòng)或者活動(dòng)范圍很小時(shí)，傳感器在觀測(cè)時(shí)間內(nèi)對(duì)目標(biāo)的距離、方位角、俯仰角探測(cè)近似相等，即(r1(1),β1(1),ε1(1))≈…≈(r1(k),β1(k),ε1(k))，此時(shí)也會(huì)導(dǎo)致矩陣A的列向量具有相關(guān)性，從而使信息矩陣J是病態(tài)的。

2.2 病態(tài)性的評(píng)估

一般地，通過(guò)條件數(shù)對(duì)方程病態(tài)性進(jìn)行評(píng)估。對(duì)于(5)式，小擾動(dòng)引起的估計(jì)相對(duì)系統(tǒng)誤差有如下關(guān)系：

(8)

當(dāng)A非奇異時(shí)，A條件數(shù)定義如下：

(9)

式中：λl為矩陣A的第l個(gè)特征值。

可見(jiàn)，系數(shù)矩陣A的最大最小特征值越分散，cond{A}越大，解向量b的相對(duì)誤差越大。條件數(shù)是一個(gè)相對(duì)數(shù)，其值的大小表示了最小特征根相對(duì)于最大特征根“小”的程度，即條件數(shù)從某種程度上度量了病態(tài)性的嚴(yán)重程度。但條件數(shù)不能判定觀測(cè)矩陣A的數(shù)據(jù)列之間存在幾個(gè)相關(guān)關(guān)系，也不能確切斷定出每個(gè)相關(guān)性到底存在于哪些數(shù)據(jù)列之間，這些信息對(duì)于全面把握病態(tài)性產(chǎn)生的機(jī)理和實(shí)質(zhì)，以及判定哪個(gè)系統(tǒng)誤差參數(shù)的可觀測(cè)性高低具有重要的參考價(jià)值。

定義條件指標(biāo)：

(10)

稱ηj法矩陣的第j個(gè)條件指標(biāo)，λ1為最大的特征值，顯然ηj≥1.

條件指標(biāo)直接反映了具體的系統(tǒng)誤差參數(shù)項(xiàng)估計(jì)受測(cè)量噪聲干擾的影響程度。若出現(xiàn)高條件指標(biāo)，則矩陣A的該列元有微小擾動(dòng)時(shí)，該系統(tǒng)誤差參數(shù)項(xiàng)的估計(jì)結(jié)果將引起相當(dāng)大的變化，從而觀測(cè)矩陣A的數(shù)據(jù)列之間存在相關(guān)關(guān)系。大量模擬研究表明[18-19]，若病態(tài)性很弱，則條件指標(biāo)小于100；若病態(tài)性較強(qiáng)，則條件指標(biāo)在100～1 000之間；若病態(tài)性很嚴(yán)重，則條件指標(biāo)在1 000以上。

3 多傳感器系統(tǒng)誤差CSVD穩(wěn)健估計(jì)算法

由系統(tǒng)誤差參數(shù)b的Fisher信息矩陣定義可知，J是一個(gè)實(shí)對(duì)稱矩陣，根據(jù)矩陣奇異值分解定理，對(duì)于任意實(shí)對(duì)稱矩陣Jm×m，存在正交矩陣Vm×m使得：

J=VSVT，

(11)

J的Mooer-Penrose廣義逆為

J+=VS-1VT，

(12)

系統(tǒng)誤差參數(shù)估計(jì)的方差為

(13)

向量形式：

(14)

TSVD法就是將較小的奇異值舍棄，即Σ=diag(λ1,λ2,…,λr)舍棄較小的，使其為0，以減少參數(shù)解的方差、增加解算的穩(wěn)定性。若令D=diag(λ1,λ2,…,λr1) (r1

(15)

(15)式向量形式表示為

(16)

在具體采用TSVD法時(shí)，主要需要考慮按什么原則去截?cái)唷?/p>

(17)

參數(shù)估計(jì)的向量形式為

(18)

參數(shù)估計(jì)的方差為

(19)

從(18)式和(19)式可以看出，適當(dāng)選擇α可以減少估計(jì)的方差、增加估計(jì)的穩(wěn)定性。但上述修正方法存在如下缺陷：1)當(dāng)奇異值分布呈均勻下降型時(shí)，MSVD法對(duì)不可靠的成分進(jìn)行了抑制，它不但修正了較小的奇異值，也修正了相對(duì)較大的奇異值，從而使模型中的確定成分發(fā)生了畸變，降低了分辨率；2)TSVD法將小于方差門(mén)限值的奇異值全部截掉，將引起分辨率的明顯變壞，因?yàn)楸Ａ舻暮徒氐舻某煞植顒e不明顯。

在進(jìn)行奇異值分解之后，相對(duì)大的奇異值及相應(yīng)的左右特征向量表示模型參數(shù)中比較肯定的和可靠的部分，很小的奇異值和相應(yīng)的特征向量則表示不可靠的部分。TSVD法在處理病態(tài)問(wèn)題的奇異值時(shí)，刪除了小的奇異值及對(duì)應(yīng)的特征向量，實(shí)際上刪掉了模型參數(shù)中不可靠的部分，以減少解的方差。但是這種做法同時(shí)也可能嚴(yán)重?fù)p害了解估計(jì)的分辨率。

對(duì)于多傳感器空間配準(zhǔn)問(wèn)題，同時(shí)對(duì)傳感器的測(cè)距和測(cè)向系統(tǒng)誤差進(jìn)行估計(jì)時(shí)，其觀測(cè)矩陣的特征值的大小順序呈階梯狀態(tài)的分布，如(20)式所示：

η1≤η2≤η3<η4≤η5≤η6.

(20)

理想的做法是修正奇異值較小的，而對(duì)較大的奇異值不進(jìn)行修正或給予不同于較小的奇異值的修正。本文綜合了TSVD法和MSVD法的優(yōu)點(diǎn)，提出了一種新的CSVD方案，即對(duì)大的奇異值盡量趨近于用TSVD法進(jìn)行修正，對(duì)小的奇異值趨近于用MSVD法進(jìn)行修正。具體做法為：對(duì)超過(guò)門(mén)限值t的奇異值保持不變，對(duì)小于門(mén)限值t的奇異值進(jìn)行修正，修正公式為

(21)

相應(yīng)的修改后的奇異值為

Σδ=diag(δ1,δ2,…,δr).

(22)

相應(yīng)的配準(zhǔn)偏差估計(jì)結(jié)果為

(23)

1)配準(zhǔn)估計(jì)結(jié)果的分辨率和方差。偏差估計(jì)結(jié)果(23)式的分辨率為

(24)

當(dāng)λp≥t時(shí)，(23)式解估計(jì)的方差為

(25)

由方差公式(24)式可以看出，與LS法估計(jì)的方差相比，新的解估計(jì)的方差減小了。本文的方法兼具了MSVD法和TSVD法的優(yōu)點(diǎn)，通過(guò)降低分辨率減小了解的方差，解估計(jì)的質(zhì)量得到了提高。

2)最小奇異值門(mén)限t的確定。判斷病態(tài)觀測(cè)方程式的Fisher信息矩陣ATA的條件數(shù)大于103，可以得到

(26)

設(shè)奇異值λq(q=1,2,…,r)滿足

(27)

則選擇最小奇異值門(mén)限t為

t=λq.

(28)

4 仿真驗(yàn)證

4.1 可觀測(cè)性較好的情況

目標(biāo)1的初始位置[121.05°，26.6°，2 km]，運(yùn)動(dòng)速度[121.2°，26.5°,0 km]；平臺(tái)1初始位置[121.2°，26.5°，0 km]，運(yùn)動(dòng)速度[0 m/s，0 m/s，0 m/s]；平臺(tái)2初始位置[121.4°，26.5°，0 km]，運(yùn)動(dòng)速度[0 m/s，0 m/s，0 m/s]；傳感器1探測(cè)系統(tǒng)誤差：Δr1=-500 m,Δβ1=0.9°,Δε1=0.5°，噪聲σr1=100 m,σβ1=0.3°,σε1=0.3°；傳感器2探測(cè)系統(tǒng)誤差：Δr2=400 m,Δβ2=0.9° ,Δε2=-0.5°，噪聲σr2=100 m,σβ2=0.3° ,σε2=0.3°.

仿真場(chǎng)景如圖2所示，仿真結(jié)果如圖3～圖9所示。

圖2 仿真場(chǎng)景示意圖Fig.2 Schematic diagram of simulation scene

圖3 條件指標(biāo)變化情況Fig.3 Changes of condition indexes

圖4 CSVD修正之后條件指標(biāo)變化情況Fig.4 Changes of modified condition indexes

圖5 兩傳感器距離系統(tǒng)誤差估計(jì)結(jié)果Fig.5 Estimated range deviations of two sensors

圖6 兩傳感器方位角系統(tǒng)誤差估計(jì)結(jié)果Fig.6 Estimated azimuth deviations of two sensors

圖7 兩傳感器俯仰角系統(tǒng)誤差估計(jì)結(jié)果Fig.7 Estimated elevation deviations of two sensors

圖8 CSVD實(shí)時(shí)配準(zhǔn)前后兩傳感器探測(cè)結(jié)果Fig.8 Detection results of two sensors before and after real time CSVD registration

圖9 實(shí)時(shí)配準(zhǔn)前后兩傳感器誤差比較Fig.9 Comparison of detection errors of two sensors before and after real time CSVD registration

4.2 可觀測(cè)性較差的情況

目標(biāo)1初始位置[120.95°，26.6°，2 km]，運(yùn)動(dòng)速度[0 m/s，100 m/s，0 m/s]；平臺(tái)1初始位置[121.2°，26.5°，0 km]，運(yùn)動(dòng)速度[0 m/s，0 m/s，0 m/s]；平臺(tái)2初始位置[121.4°，26.5°，0 km]，運(yùn)動(dòng)速度[0 m/s，0 m/s，0 m/s]；傳感器1探測(cè)系統(tǒng)誤差：Δr1=-500 m,Δβ1=0.9°,Δε1=0.5°，噪聲σr1=100 m,σβ1=0.3°,σε1=0.3°；傳感器2探測(cè)系統(tǒng)誤差：Δr2=400 m, Δβ2=0.9°,Δε2=-0.5°，噪聲σr2=100 m,σβ2=0.3°,σε2=0.3°.

仿真場(chǎng)景如圖10所示，仿真結(jié)果如圖11～圖16所示。

圖10 仿真場(chǎng)景示意圖Fig.10 Schematic diagram of simulation scene

圖11 條件指標(biāo)變化情況Fig.11 Change of condition indexes

圖12 CSVD修正之后條件指標(biāo)變化情況Fig.12 Changes of modified condition indexes

圖13 兩傳感器距離系統(tǒng)誤差估計(jì)結(jié)果Fig.13 Estimated range deviations of two sensors

圖14 兩傳感器方位角系統(tǒng)誤差估計(jì)結(jié)果Fig.14 Estimated azimuth deviations of two sensors

圖15 兩傳感器俯仰角系統(tǒng)誤差估計(jì)結(jié)果Fig.15 Estimated elevation deviations of two sensors

圖16 CSVD實(shí)時(shí)配準(zhǔn)前后兩傳感器探測(cè)結(jié)果Fig.16 Detection results of two sensors before and after real time CSVD registration

圖17 實(shí)時(shí)配準(zhǔn)前后兩傳感器誤差比較Fig.17 Comparison of detection errors of two sensors before and after real time CSVD registration

對(duì)比圖3和圖17中兩種場(chǎng)景的條件指標(biāo)可以看出：1)第1種場(chǎng)景的條件指標(biāo)明顯比第2種場(chǎng)景的要低；2)兩個(gè)場(chǎng)景之間一開(kāi)始均存在3個(gè)變量之間的相關(guān)性，開(kāi)始階段均只能估計(jì)出相對(duì)的系統(tǒng)誤差；3)隨著數(shù)據(jù)的積累，第1種場(chǎng)景的可觀測(cè)性明顯變好，當(dāng)數(shù)據(jù)量超過(guò)100時(shí)所有的待估計(jì)參數(shù)之間是完全獨(dú)立的，此時(shí)使用經(jīng)典的算法均可得到有效的結(jié)果。

5 結(jié)論

在多傳感器空間配準(zhǔn)中，信息矩陣的病態(tài)性是影響系統(tǒng)誤差參數(shù)估計(jì)質(zhì)量的主要因素之一。本文針對(duì)LS法系統(tǒng)誤差估計(jì)過(guò)程中由于信息矩陣的病態(tài)性導(dǎo)致估計(jì)結(jié)果無(wú)效的情況，分析LS法系統(tǒng)誤差估計(jì)的病態(tài)性產(chǎn)生原因，并根據(jù)Fisher信息矩陣奇異值的分布為階梯型分布特點(diǎn)，提出一種新的CSVD方案對(duì)奇異特征值進(jìn)行修正，該算法兼顧參數(shù)估計(jì)的分辨率與方差，實(shí)現(xiàn)了對(duì)觀測(cè)方程中不確定成分的有效抑制和確定成分的有效保留，有效解決了多傳感器空間配準(zhǔn)系統(tǒng)誤差實(shí)時(shí)估計(jì)問(wèn)題，并從兩種不同的仿真場(chǎng)景進(jìn)行了驗(yàn)證比較，結(jié)果表明：當(dāng)信息矩陣的條件數(shù)明顯分梯形分布時(shí)，該方法可顯著提高系統(tǒng)誤差估計(jì)結(jié)果的精度和穩(wěn)定性。

References)

[1] Bo Y, Chen Z, Yin M, et al. Improved different dimen-sional sensors combined space registration algorithm[J]. Mathematical Problems in Engineering, 2015(1):1-9.

[2] 程志剛. 多雷達(dá)組網(wǎng)系統(tǒng)空間配準(zhǔn)方法研究[J]. 艦船科學(xué)技術(shù), 2015, 37(12):131-134.

CHENG Zhi-gang. Research on radar space registration method for multi-radar network system[J]. Ship Science and Technology, 2015, 37(12):131-134.(in Chinese)

[3] 胡勤振, 蘇洪濤, 劉子威. 配準(zhǔn)誤差下的多基地雷達(dá)目標(biāo)檢測(cè)算法[J]. 電子與信息學(xué)報(bào), 2017,39(1):87-94.

HU Qin-zhen, SU Hong-tao, LIU Zi-wei. Target detection algorithm for multistatic radar with registration errors[J]. Journal of Electronics & Information Technology, 2017,39(1):87-94.(in Chinese)

[4] 史江林, 周學(xué)平, 劉長(zhǎng)海,等. 基于EM算法的聯(lián)合數(shù)據(jù)關(guān)聯(lián)與空間配準(zhǔn)[J]. 現(xiàn)代雷達(dá), 2016, 38(1):51-55.

SHI Jiang-lin, ZHOU Xue-ping, LIU Chang-hai, et al. Joint data association and spatial registration based on EM[J]. Modern Radar, 2016, 38(1):51-55.(in Chinese)

[5] 方峰, 蔡遠(yuǎn)利. 機(jī)動(dòng)目標(biāo)多傳感器組網(wǎng)空間配準(zhǔn)方法[J]. 固體火箭技術(shù), 2016, 39(4):574-579.

FANG Feng, CAI Yuan-li. Multi-sensor space registration for maneuvering target tracking [J]. Journal of Solid Rocket Technology, 2016, 39(4):574-579. (in Chinese)

[6] Okello N N, Chall S L.A joint radar registration and track-to-track fusion for distributed trackers[J]. IEEE Transactions on Aerospace and Electronic Systems, 2004, 40(3):808-823.

[7] Herman S M, Poore A B.Nonlinear least-squares estimation for radar and navigation biases[J]. Proceedings of SPIE, Signal and Data Processing of Small Targets, 2006, 6236:1701-1717.

[8] Okello N N, Ristic B. Maximum likelihood registration for multiple dissimilar radars[J]. IEEE Transactions on Aerospace and Electronic Systems, 2003, 39(3):1235-1240.

[9] 潘江懷, 李洪梅, 何佳洲.三維傳感器組網(wǎng)偏差估計(jì)方法[J].火力與指揮控制, 2008, 4(4):104-106.

PAN Jiang-huai, LI Hong-mei, HE Jia-zhou.Bias estimation in 3D sensor networks[J]. Fire Control and Command Control, 2008, 4(4):104-106.(in Chinese)

[10] 宋強(qiáng), 崔亞奇, 何友.反饋式多目標(biāo)多傳感器系統(tǒng)誤差融合估計(jì)技術(shù)[J]. 宇航學(xué)報(bào), 2011, 32(1):115-122.

SONG Qiang, CUI Ya-qi, HE You.Multi-target multi-sensor systematic error fusion estimation technique with feedback[J]. Journal of Astronautics, 2011, 32(1):115-122 .(in Chinese)

[11] 王成飛, 王航宇, 石章松, 等. 基于高精度導(dǎo)航設(shè)備的海上多傳感器配準(zhǔn)算法[J].控制理論與應(yīng)用, 2011, 28(4):497-503.

WANC Cheng-fei, WANG Hang-yu, SHI Zhang-song, et al. Maritime sensor registration algorithm based on high-precision navigation equipment[J]. Control Theory and Applications, 2011, 28(4):497-503.(in Chinese)

[12] 崔亞奇, 熊偉, 何友.基于MLR的機(jī)動(dòng)平臺(tái)傳感器誤差配準(zhǔn)算法[J].航空學(xué)報(bào), 2012, 33(1):118-128.

CUI Ya-qi, XIONG Wei, HE You.Mobile platform sensor registration algorithm based on MLR[J]. Acta Aeronautica et Astronautica Sinica, 2012, 33(1):118-128.(in Chinese)

[13] 連峰, 韓崇昭, 彭一峰, 等.基于廣義似然比的自適應(yīng)在線配準(zhǔn)算法[J]. 控制與決策, 2009, 24(1):23-28.

LIAN Feng, HAN Chong-zhao, PENG Yi-feng, et al.Adaptive on-line registration algorithm based on GLR[J]. Control and Decision, 2009, 24(1):23-28.(in Chinese)

[14] 吳衛(wèi)華, 江晶. WGS-84坐標(biāo)系下多空基無(wú)源傳感器最大似然配準(zhǔn)[J]. 系統(tǒng)工程與電子技術(shù), 2015, 37(2): 304-309.

WU Wei-hua, JIANG Jing. Maximun likelihood registration for passive sensors of multiple airborne platforms in WGS-84[J]. Systems Engineering and Electronics, 2015, 37(2):304-309.(in Chinese)

[15] Huang D, Leung H.A pseudo-measurement approach to simultaneous registration and track fusion[J]. IEEE Transactions on Aerospace and Electronic System, 2012, 48(3):2315-2331.

[16] Steafano F, Alfonso F, Fulvio G, et al.Least squares estimation and Cramer-Rao type lower bounds for relative sensor registration process[J]. IEEE Transactions on Signal Processing, 2011, 59(3):1075-1087.

[17] Boer Y，Driessen J N. Multi-target particle filter track before detect application[J]. IEEE Proceedings Radar, Sonar and Navigation, 2004, 151(6):351-357.

[18] Ong H T.Sensor registration using airlanes maximum likelihood solution[J]. Proceedingd of SPIE, Signal and Data Processing of Small Targets, 2003,5204:390-401.

[19] Watson G A, Rice T R.Sensor alignment and compensation for composite tracking[J]. Proceedings of SPIE, Signal and Data Processing of Small Targets, 2002, 4728: 354-367.

Ill-conditionAnalysisandRobustEstimationforMulti-sensorSpatialRegistration

PAN Jiang-huai, QIAO Hui
(Jiangsu Automation Research Institute, Lianyungang 222006，Jiangsu， China)

In order to improve the accuracy and stability of multi-sensor spatial registration, the causes of the ill condition of the multi sensor spatial registration are analyzed based on LS estimation for multi-sensor spatial registration. The analysis result shows that the ill-condition of the Fisher information matrix is the key factor which influences the quality of system error parameter estimation. The conditional index of information matrix is used to measure the ill-condition of spatial registration. A new algorithm of compound modified singular value decomposition (CSVD) for robust estimation is proposed according to the characteristics of trapezia distribution of the information matrix singular values. The algorithm can effectively suppress uncertainty components and retain deterministic components in the observation equation by modifing only minor conditional index, and can effectively solve the multi sensor real time estimation of spatial registration. Two different scenarios of the estimatied results show that CSVD method can be used to significantly improve the estimation accuracy and stability of system error parameters.

information processing technology; spatial registration; ill-condition; robust estimation; singular value decomposition; multi-sensor

2016-03-02

國(guó)家“973”計(jì)劃項(xiàng)目(613101)

潘江懷(1982—)，男，高級(jí)工程師。E-mail：panjianghuai@163.com

TN953+.6

A

1000-1093(2017)10-1965-09

10.3969/j.issn.1000-1093.2017.10.012