楊婧，李銀伢，戚國慶，盛安冬

（南京理工大學(xué)，南京210014）

純方位目標(biāo)運(yùn)動分析可觀測性研究*

楊婧，李銀伢，戚國慶，盛安冬

（南京理工大學(xué)，南京210014）

純方位目標(biāo)運(yùn)動分析（BOTMA）僅利用方位信息實現(xiàn)對目標(biāo)狀態(tài)參數(shù)的估計，是一種有效的無源被動定位跟蹤方法。純方位系統(tǒng)中的可觀測性條件、目標(biāo)跟蹤與估計策略、觀測器最優(yōu)機(jī)動軌跡構(gòu)成了BOTMA的核心研究內(nèi)容?？捎^測性問題是BOTMA首先必須解決的關(guān)鍵問題，可觀測性條件是后續(xù)目標(biāo)定位與跟蹤的前提和基礎(chǔ)。介紹了純方位系統(tǒng)中可觀測性的基本概念，從幾何方法、線性與非線性方法、數(shù)值分析方法等角度，對BOTMA的可觀測性研究成果進(jìn)行了系統(tǒng)的總結(jié)和評述，最后對這一領(lǐng)域的研究提出了新的展望。

可觀測性，BOTMA，被動跟蹤

0 引言

純方位目標(biāo)定位與跟蹤是最經(jīng)典的被動式無源定位跟蹤技術(shù)，具有隱蔽性好、作用距離遠(yuǎn)、抗干擾能力強(qiáng)等優(yōu)點(diǎn)，在無線電探測、被動聲納、紅外等領(lǐng)域有著廣泛的應(yīng)用背景［1-4］。純方位目標(biāo)運(yùn)動分析（Bearings-Only Target Motion Analysis，BOTMA）指的是，利用傳感器被動地測量目標(biāo)的方位信息，實時估計出目標(biāo)的位置和運(yùn)動（速度、航向）參數(shù)，達(dá)到對目標(biāo)進(jìn)行被動定位與跟蹤的目的。BOTMA旨在掌握運(yùn)動目標(biāo)參數(shù)（狀態(tài)）隨時間的變化規(guī)律，建立和預(yù)測目標(biāo)運(yùn)動狀態(tài)的歷程和趨勢，為后續(xù)武器系統(tǒng)的進(jìn)一步動作，諸如攻擊要素的獲取奠定基礎(chǔ)。BOTMA由于其固有的非線性和弱可觀測性，使得該領(lǐng)域的研究具有極大的難度和挑戰(zhàn)性。BOTMA的研究內(nèi)容主要包括三個方面：目標(biāo)可觀測性條件，觀測器最優(yōu)機(jī)動軌跡以及目標(biāo)定位跟蹤算法［5］。其中，可觀測性研究是后續(xù)研究的前提和基礎(chǔ)。在純方位系統(tǒng)中，BOTMA的可觀測性問題本質(zhì)上是非線性系統(tǒng)的可觀測性問題。與線性系統(tǒng)不同，非線性系統(tǒng)具有局部可觀測性和弱可觀測性等特點(diǎn)，其研究過程涉及李代數(shù)等繁復(fù)的數(shù)學(xué)內(nèi)容［6］，使得該類系統(tǒng)可觀測性的分析研究更為困難。

圖1 純方位系統(tǒng)及其可觀測性研究趨勢

圖1給出了SCI和EI兩大數(shù)據(jù)庫中涉及純方位系統(tǒng)及其可觀測性研究的檢索結(jié)果。其中，“條件一”標(biāo)識曲線代表純方位系統(tǒng)研究的檢索結(jié)果，檢索關(guān)鍵詞為“bearings-only”；“條件二”標(biāo)識曲線代表純方位系統(tǒng)可觀測性研究的檢索結(jié)果，檢索關(guān)鍵詞為“bearings-only&observability”。從圖1可大致看出，純方位系統(tǒng)及其可觀測性問題的研究起源于上世紀(jì)六七十年代［7-8］，在經(jīng)歷了相對平緩的發(fā)展過程后，近十年來逐漸成為研究熱點(diǎn)。而分析現(xiàn)有文獻(xiàn)可知，可觀測性問題的研究工作在八九十年代取得了較大進(jìn)展，許多經(jīng)典成果為后續(xù)的研究奠定了基礎(chǔ)［9-11］。

1 可觀測性概念

可觀測性概念是卡爾曼在20世紀(jì)60年代首先提出的，表征系統(tǒng)狀態(tài)能否由其輸出完全反映的問題［12］。

在純方位系統(tǒng)中，觀測器測量的是目標(biāo)相對于觀測器的偏轉(zhuǎn)角度，因此，在笛卡爾坐標(biāo)系下，系統(tǒng)的量測方程是非線性的。以圖2所示二維情形為例，設(shè)目標(biāo)和觀測器的坐標(biāo)分別為和o），目標(biāo)和觀測器的狀態(tài)向量xT（t）和xo（t）分別為

設(shè)β（t）為t時刻觀測器測得的目標(biāo)方位角，則系統(tǒng)的狀態(tài)方程與量測方程分別為

其中，u（t）為控制量，表示觀測器相對于目標(biāo)的機(jī)動；rx和ry分別為目標(biāo)相對觀測器在x軸、y軸方向上的位置坐標(biāo)分量；w（t），v（t）分別為狀態(tài)噪聲和量測噪聲。

圖2 目標(biāo)與觀測器運(yùn)動態(tài)勢

記t時刻目標(biāo)狀態(tài)的信息集合為

相應(yīng)地，t=t0初始時刻的目標(biāo)參數(shù)集合為

2 可觀測性研究

可觀測性分析是BOTMA中最核心也是最基礎(chǔ)的研究內(nèi)容。一方面，具有良好可觀測性的純方位系統(tǒng)，可以大幅度提高觀測器對目標(biāo)狀態(tài)參數(shù)估計的可靠性和準(zhǔn)確性。

2.1 幾何分析方法

用幾何方法進(jìn)行目標(biāo)定位通?？傻玫脚c實際運(yùn)動參數(shù)相關(guān)聯(lián)的結(jié)論，比較直觀易懂。測向交叉定位法是一種發(fā)展較早，應(yīng)用廣泛的無源定位方法。文獻(xiàn)［13］對二維和三維目標(biāo)的可觀測性進(jìn)行了直觀的分析。在二維情況下，觀測器每次測量只得到方位角信息；在三維情形下，一次測量可得目標(biāo)的方位角和高低角。其結(jié)論可歸納為，目標(biāo)可觀測要求觀測器與目標(biāo)有相對運(yùn)動。文獻(xiàn)［14］介紹了雙站測向交叉定位法，基本原理與單站情形相同，同時給出了單站對勻速目標(biāo)航向的計算方法。文獻(xiàn)［15］提出線面相交定位法，比傳統(tǒng)方法具有更高的定位精度。文獻(xiàn)［16］用初等幾何方法，得到了用方位信息計算勻速直線運(yùn)動目標(biāo)距離的公式。為保證可觀測性，觀測器采用折線運(yùn)動方式。文獻(xiàn)［17］也在觀測器折線運(yùn)動條件下分析了純方位系統(tǒng)的可觀測性，實現(xiàn)了對勻速目標(biāo)的距離和速度的估計。

2.2 線性系統(tǒng)分析方法

在笛卡爾坐標(biāo)系下，BOTMA中的測量方程是非線性的，給后續(xù)分析帶來諸多不便。對測量方程進(jìn)行偽線性化處理可以將測量方程轉(zhuǎn)化為線性形式，然后運(yùn)用線性系統(tǒng)理論分析可觀測性。不考慮量測誤差，對式（3）中量測方程進(jìn)行代數(shù)變換可得

用矩陣向量形式可表示為

上式即為偽觀測方程，C（t）稱為偽觀測矩陣，其維數(shù)由狀態(tài)向量決定，y（t）是偽觀測值。相應(yīng)地，對式（3）中的狀態(tài)方程作線性化處理，可得偽線性化后的純方位系統(tǒng)動態(tài)方程為

式中，x（t）?Rn，y（t）?Rm，A（t）?Rn×n，B（t）?Rn×p，C（t）?Rm×n。因此，由式（3）描述的系統(tǒng)可觀測性問題就轉(zhuǎn)化為式（6）所示線性系統(tǒng)的可觀測性問題。

下面介紹相關(guān)的分析方法。

判據(jù)1（Gram矩陣判據(jù)）［18］考慮式（6）描述的線性時變連續(xù)系統(tǒng)，系統(tǒng)在時刻t0完全能觀測的充分必要條件是存在一個有限時刻t0?Tt，t1＞t0，使如下定義的Gram矩陣：

為非奇異，其中Φ（t，t0）為A（t）狀態(tài)轉(zhuǎn)移矩陣。

應(yīng)當(dāng)注意，上述判據(jù)中，Gram矩陣只包含C（t）和Φ（t），從控制理論的角度看，可觀測性與輸入無關(guān)［19］。但在BOTMA問題中，輸入u（t）代表觀測器相對于目標(biāo)機(jī)動，u（t）的變化會影響測得的方位角信息β（t），而觀測矩陣C（t）=［cosβ（t）-sinβ（t）0 0…］，所以輸入u（t）會影響B(tài)OTMA可觀測性能。

可觀測性的Gram矩陣判據(jù)需要先計算狀態(tài)轉(zhuǎn)移矩陣，再求Gram矩陣，計算困難且計算量較大，故有以下簡化的判據(jù)［19-20］。

判據(jù)2線性時變系統(tǒng)（6）在時刻t0完全能觀測的充分條件是存在一個有限時刻t1?Tt，t1＞t0，使

其中，M（t）=C（t）Φ（t，t0），Φ（t，t0）為A（t）的狀態(tài)轉(zhuǎn)移矩陣。

與判據(jù)2等價的判據(jù)3［19-20］如下：

判據(jù)3設(shè)線性時變系統(tǒng)（6）中，矩陣A（t），C（t）都為（n-1）階連續(xù)可導(dǎo)的，則系統(tǒng)狀態(tài)在時刻完全能觀測的充分條件是存在一個有限時刻t1?Tt，t1＞t0，使

其中

由判據(jù)1可推導(dǎo)如下等價的判據(jù)4和判據(jù)5［19-20］：

判據(jù)4系統(tǒng)（6）在時刻t0完全能觀測的充要條件是，存在一個有限時刻t1?Tt，t1＞t0，使矩陣C（t1）Φ（t1，t0）的N個列線性無關(guān)。

判據(jù)5系統(tǒng)（6）在時刻t0完全能觀測的充要條件是，對任意的非零矢量ξ，存在t?［t0，tf］，使得C（t）Φ（t，t0）ξ≠0。

Aidala［9］、Hammel［10］等利用上述線性系統(tǒng)的可觀測性判據(jù)對BOTMA可觀測性進(jìn)行分析，其方法和結(jié)論被廣泛引用。文獻(xiàn)［9］對二維情況下的勻速直線運(yùn)動目標(biāo)，將量測方程偽線性化，運(yùn)用判據(jù)2推導(dǎo)出了目標(biāo)可觀測條件：

其中，ao（t）是觀測器加速度矢量，ao（t）和v（t0）分別是目標(biāo)相對于觀測器的初始位置和速度矢量，α（t）為任意標(biāo)量時間函數(shù)。文獻(xiàn)［10］將結(jié)論推廣到了三維目標(biāo)，得到了形式類似的可觀測條件，不同的是推導(dǎo)過程中運(yùn)用了判據(jù)3，簡化了步驟。文獻(xiàn)［11］將研究對象從勻速運(yùn)動目標(biāo)擴(kuò)展到N階動態(tài)運(yùn)動目標(biāo)，用代數(shù)方法得到了可觀測的充分必要條件，推導(dǎo)過程極為簡單，避免了求解微分方程，同時指出文獻(xiàn)［9］和文獻(xiàn)［10］的可觀測判據(jù)只是必要非充分條件。而文獻(xiàn)［21］根據(jù)判據(jù)4也得出了N階動態(tài)運(yùn)動目標(biāo)不可觀測的充要條件，其形式是觀測器相對于目標(biāo)的機(jī)動向量，即系統(tǒng)輸入u關(guān)于目標(biāo)初始狀態(tài)的表達(dá)式，在實際中代表觀測器的無效機(jī)動。筆者還進(jìn)一步證明了文獻(xiàn)［9］和文獻(xiàn)［10］的結(jié)論是目標(biāo)為一階的特殊情況。同樣基于線性系統(tǒng)方法，文獻(xiàn)［22］運(yùn)用判據(jù)5，經(jīng)過簡單的代數(shù)推導(dǎo)即可得文獻(xiàn)［9-11］的結(jié)論，說明了以往線性系統(tǒng)方法在解決BOTMA可觀測性問題時所呈現(xiàn)出的繁瑣過程并不是線性系統(tǒng)方法本身的原因，而是與可觀測性判據(jù)的選用有關(guān)，其結(jié)論使得線性系統(tǒng)方法得到了廣泛應(yīng)用［23-29］。文獻(xiàn)［30］認(rèn)為這種方法不局限于對N階動態(tài)運(yùn)動目標(biāo)的分析，還可處理更廣泛一類目標(biāo)的情況，比如對于一階延遲目標(biāo)加速模型，當(dāng)觀測器的運(yùn)動也滿足一階延遲模型時，要使得目標(biāo)可觀測，只須使觀測器的機(jī)動時間常數(shù)不同于目標(biāo)，觀測器不必作比目標(biāo)更高階次的機(jī)動運(yùn)動，這與N階動態(tài)運(yùn)動目標(biāo)的可觀測條件是截然不同的。該結(jié)論同樣適用于隨機(jī)環(huán)境下，目標(biāo)加速度為Singer模型及其改進(jìn)形式的情況。此外，文獻(xiàn)［3，31］利用向量變換和線性方程組解的性質(zhì)推導(dǎo)出N階動態(tài)運(yùn)動目標(biāo)的不可觀測條件，也簡化了推導(dǎo)過程，為線性系統(tǒng)框架內(nèi)的可觀測性研究提供了新的思路。

基于線性系統(tǒng)分析方法的BOTMA可觀測性分析研究，可以得到觀測器的無效機(jī)動表達(dá)式，即可觀測條件的解析形式，相對幾何分析方法而言，研究方法在數(shù)學(xué)上更為嚴(yán)謹(jǐn)。然而，該類方法盡管在可觀測性條件上給出了解析表達(dá)式，卻一般是以目標(biāo)和觀測器的相對運(yùn)動參數(shù)（或其導(dǎo)數(shù)）而不是以便于直接獲取的方位量測來表示的。在實際工程應(yīng)用中，常希望建立直接利用方位量測數(shù)據(jù)，無需求導(dǎo)和求解目標(biāo)相對位置與運(yùn)動參數(shù)的BOTMA可觀測性判據(jù)，以簡化操作，提高工程實用性。

2.3 非線性系統(tǒng)分析方法

純方位系統(tǒng)本質(zhì)上是非線性系統(tǒng)，利用非線性系統(tǒng)分析方法對BOTMA可觀測問題可以進(jìn)行更廣泛、更深入的研究。與線性系統(tǒng)的Gram矩陣判據(jù)類似，Lee和Dunn等人提出如下可觀測性判據(jù)［32］：

判據(jù)6對于如下非線性系統(tǒng)：

如果對于凸集S?Rn上的所有x0，都有

是正定的，則系統(tǒng)在S上是完全可觀測的。其中

將上述判據(jù)與線性系統(tǒng)Gram矩陣判據(jù)比較可知，在非線性系統(tǒng)中用相應(yīng)的Jacobin矩陣代替了線性系統(tǒng)的狀態(tài)轉(zhuǎn)移矩陣和測量矩陣。

將非線性系統(tǒng)弱可觀測和局部可觀測的概念相結(jié)合，還可得出與線性系統(tǒng)類似的局部弱可觀測性的秩判據(jù)如下［2，33］：

判據(jù)7設(shè)非線性系統(tǒng)為

其中，x?Xn?Rn，y?Rm?？捎^測性矩陣定義為

對x0?Xn，如果rank（Q）=n，則稱系統(tǒng)在x0處是局部弱可觀測的；如果對任意x?Xn，都有rank（Q）= n，則稱系統(tǒng)是局部弱可觀測的。其中，N是如下定義的非線性算子：

上述非線性系統(tǒng)分析方法在BOTMA的可觀測性研究中得到廣泛應(yīng)用［34-37］。特別是文獻(xiàn)［38］從實際應(yīng)用出發(fā)，證明了非線性離散系統(tǒng)可觀測性判據(jù)的“一致性”要求在系統(tǒng)初始狀態(tài)和輸入屬于緊集時并非必須滿足。只要系統(tǒng)狀態(tài)兩兩可區(qū)分，且對于每個狀態(tài)變量和輸入，存在一定長度的觀測序列使得可觀測性判據(jù)的秩條件成立，那么整個系統(tǒng)就是可觀測的。這種將可觀測性判據(jù)“化整為零”的思路是符合實際需求的。

非線性系統(tǒng)方法在一定程度上可以看作是基于Gram矩陣的線性系統(tǒng)方法的推廣，因此，適用于線性和非線性目標(biāo)模型，使可觀測性問題的研究擴(kuò)展到了更廣闊的領(lǐng)域。在數(shù)學(xué)原理上，用矩陣的正定或滿秩條件作為可觀測性判據(jù)。不足的是，當(dāng)處理高階目標(biāo)時會產(chǎn)生較大的計算量。

2.4 數(shù)值分析方法

前面介紹的可觀測性分析方法一般適用于判斷系統(tǒng)全體狀態(tài)變量在某一時刻是否可觀測，是一種總體的定性分析。對于時不變系統(tǒng)，若系統(tǒng)是可觀測的，并不能得知系統(tǒng)各狀態(tài)的可觀測程度；而對于時變系統(tǒng)，則需要分析從0時刻到t時刻的可觀測性矩陣的秩或正定性，計算量很大。基于線性代數(shù)理論的數(shù)值方法可以對可觀測性進(jìn)行更為深入細(xì)致的定量分析。

文獻(xiàn)［39］先計算系統(tǒng)可觀測性矩陣Q，然后將正定對稱陣QTQ單位化后的特征值λi定義為其所對應(yīng)的狀態(tài)向量或狀態(tài)向量線性組合的可觀測度。當(dāng)0≤λi≤1，且λi越大，其特征向量所對應(yīng)的狀態(tài)向量或狀態(tài)向量線性組合的可觀測度越好；λi=0，則所對應(yīng)的狀態(tài)向量或狀態(tài)向量組合不可觀測。文獻(xiàn)［40-41］對可觀測性矩陣進(jìn)行奇異值分解，并用奇異值大小判斷可觀測度。當(dāng)可觀測性矩陣是實矩陣時，有QH=QT，這與文獻(xiàn)［39］的方法實質(zhì)上是相同的。文獻(xiàn)［42］根據(jù)線性代數(shù)理論定義了系統(tǒng)可觀測矩陣的偽逆矩陣，將狀態(tài)向量解的表達(dá)式分為固定部分和自由部分，自由部分決定系統(tǒng)的可觀測性，給出了計算不可觀測程度的公式。

文獻(xiàn)［43］用系統(tǒng)可觀測矩陣的條件數(shù)定義可觀測度，條件數(shù)越大，則可觀測度越差；若條件數(shù)是無窮大，則系統(tǒng)不可觀測；條件數(shù)等于1時，可觀測度最好。特別對于時變系統(tǒng)，文獻(xiàn)［43］通過計算其局部可觀測矩陣的條件數(shù)定量地分析系統(tǒng)在某些時刻的可觀測性能，從而可對系統(tǒng)在任意時段內(nèi)的可觀測性進(jìn)行定量分析。文獻(xiàn)［44］基于條件數(shù)給出一種非線性系統(tǒng)可觀測性度量方法，并指出系統(tǒng)可觀測度依賴于狀態(tài)模型和觀測模型，與觀測數(shù)據(jù)本身無關(guān)，可作為比較觀測模型優(yōu)劣的一個性能準(zhǔn)則。

可觀測矩陣的奇異值（或特征值）分析方法不僅可以定性判斷整個系統(tǒng)是否可觀測，還可定量分析系統(tǒng)各狀態(tài)的可觀測度；而基于條件數(shù)的方法適用于研究時變系統(tǒng)在某一時間段內(nèi)的可觀測性，這些都與實際工程中系統(tǒng)可觀測性及后續(xù)研究的需求相符。

2.5 基于分段勻速直線運(yùn)動模型的分析方法

在工程實際中，復(fù)雜的系統(tǒng)模型不利于問題的分析和求解，對模型進(jìn)行簡化操作是一種有效的解決方法。文獻(xiàn)［45］認(rèn)為，一些線性時變系統(tǒng)用分段定常系統(tǒng)（Piece-Wise Constant System，PWCS）模型近似時，可保留系統(tǒng)特性，引入了簡化的可觀測性矩陣分析PWCS系統(tǒng)的可觀測性。文獻(xiàn)［46］改進(jìn)了文獻(xiàn)［45］的方法，減少了計算量。以上述研究為依據(jù)，文獻(xiàn)［47］用分段勻速直線運(yùn)動模型描述目標(biāo)和觀測器的運(yùn)動，分析了多種情形下純方位系統(tǒng)的可觀測性問題。其主要結(jié)論有：（1）如果目標(biāo)做勻速直線運(yùn)動，觀測器運(yùn)動軌跡由多段勻速直線運(yùn)動構(gòu)成，則目標(biāo)是可觀測的，當(dāng)且僅當(dāng)觀測器每一段軌跡測得的目標(biāo)方位變化率不全為零。（2）如果目標(biāo)和觀測器的運(yùn)動軌跡都由多段勻速直線運(yùn)動構(gòu)成，那么目標(biāo)是可觀測的，當(dāng)且僅當(dāng)系統(tǒng)滿足下列條件：①目標(biāo)和觀測器的所有機(jī)動時刻都不同且目標(biāo)的機(jī)動時刻已知；②觀測器軌跡至少包括2段勻速直線運(yùn)動；③不存在方位變化率為零的分段。（3）如果目標(biāo)和觀測器的所有機(jī)動時刻都不同，目標(biāo)的機(jī)動時刻未知且觀測器軌跡沒有方位變化率為零的分段，那么目標(biāo)是可觀測的。

文獻(xiàn)［47］認(rèn)為如果目標(biāo)運(yùn)動軌跡由多段勻速直線運(yùn)動構(gòu)成，觀測器做勻速直線運(yùn)動，那么目標(biāo)是不可觀測的。文獻(xiàn)［48-50］的研究則表明如果觀測器的速度矢量vo和目標(biāo)軌跡中相鄰的兩段速度矢量vs，1，vs，2滿足voT（vs，1-vs，2）≠0，那么此時目標(biāo)可觀測。文獻(xiàn)［51-52］進(jìn)一步研究了目標(biāo)軌跡包含勻速轉(zhuǎn)彎（Constant Turn，CT）運(yùn)動，而觀測器仍做勻速直線運(yùn)動的情形，結(jié)合文獻(xiàn)［48］有如下結(jié)論：①只要觀測器不是靜止的，則任何做CT運(yùn)動的目標(biāo)都是可觀測的；②如果觀測器靜止，則存在測量方位角與真正目標(biāo)同位相似的假目標(biāo)；③BOTMA可以區(qū)分CT運(yùn)動目標(biāo)和CV目標(biāo)，即②中假目標(biāo)不會是CV目標(biāo)；④如果目標(biāo)運(yùn)動軌跡由多段CV和CT運(yùn)動組成，那么在至少有一段為CT運(yùn)動的條件下，目標(biāo)的整個軌跡是可觀測的。

2.6 有噪條件下的可觀測性分析方法

BOTMA的可觀測性研究一般只考慮無噪聲影響的確定性系統(tǒng)，文獻(xiàn)［4］則將實際觀測量視為受隨機(jī)噪聲干擾的隨機(jī)變量，認(rèn)為純方位系統(tǒng)實質(zhì)上是狀態(tài)方程連續(xù)、量測方程離散的非線性隨機(jī)系統(tǒng)，其觀測性能更適宜用“好”和“差”來描述?？紤]到估計量的Fisher信息矩陣（Fisher Information Matrix，F(xiàn)IM）是衡量系統(tǒng)估計量好壞的性能指標(biāo)，因此，可作為判定純方位系統(tǒng)可觀測性的依據(jù)。當(dāng)FIM不可逆時，系統(tǒng)是不可觀測的；若FIM可逆，但具有大條件數(shù)，此時系統(tǒng)雖是可觀測的，但觀測性“較差”，反之若FIM可逆且不存在大條件數(shù)，此時系統(tǒng)是可觀測的，且具有“好”的觀測性。文獻(xiàn)［53-54］還引入高斯條件下狀態(tài)不確定超橢球體積，作為系統(tǒng)可觀測度大小的度量。文獻(xiàn)［55-56］對FIM的可逆性與非線性系統(tǒng)在無噪條件下的可觀測性之間的關(guān)系進(jìn)行了討論，結(jié)果表明：如果在高斯獨(dú)立噪聲條件下，非線性系統(tǒng)的FIM在x0處可逆，那么相應(yīng)的無噪系統(tǒng)在x0處是局部可觀測的。應(yīng)當(dāng)注意該結(jié)論的逆命題不成立，即系統(tǒng)在某點(diǎn)處可觀測并不能保證FIM在該點(diǎn)的可逆性。而對于線性系統(tǒng)，局部可觀測性與全局可觀測性是等價的。

FIM可以看作對系統(tǒng)整體可觀測性的定量分析，利用數(shù)值方法則可對每個狀態(tài)變量分別進(jìn)行定量的可觀測性分析，2.4節(jié)已對此詳細(xì)總結(jié)，這里不再贅述。需要說明的是，文獻(xiàn)［43，57］將系統(tǒng)可觀測矩陣的條件數(shù)定義為可觀測度，看似與文獻(xiàn)［4，53］的可觀測度定義不同，實際上二者原理相通［58］。FIM與可觀測矩陣相比，前者包含了測量噪聲的協(xié)方差，適合于分析有噪條件下的可觀測性能，而后者適用于理想的無噪條件。

2.7 可觀測性分析其他方法

近年來，對純方位系統(tǒng)的可觀測性問題，出現(xiàn)了一些新的研究思路，進(jìn)一步拓展了BOTMA可觀測性分析方法，豐富了BOTMA可觀測性理論。文獻(xiàn)［59］討論一般目標(biāo)的可觀測性問題。這里“一般”的含義是指目標(biāo)的運(yùn)動假設(shè)突破了以往勻速直線運(yùn)動等形式的限制，擴(kuò)展為更一般各種復(fù)雜運(yùn)動。指出在二維純方位系統(tǒng)中，系統(tǒng)可觀測的必要條件是觀測器作比目標(biāo)更為復(fù)雜的運(yùn)動，并且獨(dú)立方位數(shù)目等于待定參數(shù)數(shù)目。這一結(jié)論不僅有助于對純方位系統(tǒng)本質(zhì)的認(rèn)識，對后續(xù)跟蹤算法和觀測器機(jī)動軌跡研究等都具有重要意義。連續(xù)和離散系統(tǒng)是純方位系統(tǒng)的兩種不同的表現(xiàn)形式，對于線性連續(xù)定常系統(tǒng)離散化后是否依然可觀測，文獻(xiàn)［18］給出了判斷依據(jù)。但對線性時變系統(tǒng)和非線性系統(tǒng)而言，對應(yīng)的可觀測判據(jù)類比或轉(zhuǎn)換比較困難，不利于實際應(yīng)用。文獻(xiàn)［60］證明了離散純方位系統(tǒng)僅依賴方位測量的可觀測判據(jù)，文獻(xiàn)［61］則在連續(xù)系統(tǒng)下推導(dǎo)了該判據(jù)，表達(dá)形式與離散系統(tǒng)判據(jù)類似，但作者研究僅限于一階目標(biāo)，結(jié)論有待進(jìn)一步推廣。

對于海上運(yùn)動目標(biāo)，傳統(tǒng)分析方法［62-63］一般忽略了聲音在海水中的傳播時間，認(rèn)為聲速為無限值。文獻(xiàn)［64］對此提出了質(zhì)疑，考慮到聲音的傳播需要一定時間，所以方位的測量必定與聲速有關(guān)。由于海水中的聲速受諸多因素影響，筆者將其設(shè)定為一個假想的平均速度。當(dāng)目標(biāo)做勻速直線運(yùn)動且平均聲速為有限常數(shù)時，證明了目標(biāo)在觀測器不機(jī)動時也是可觀測的，且在相同的觀測器機(jī)動條件下，目標(biāo)可觀測性要優(yōu)于認(rèn)為聲速為無限值的情形。文章的研究結(jié)果表明，聲速有限的假設(shè)包含了目標(biāo)的狀態(tài)信息，對于可觀測性分析有重要價值。

基于圖論的分析方法［65-66］常用于電力系統(tǒng)的狀態(tài)可觀測性分析，具有計算量小，速度快等優(yōu)點(diǎn)。對于多站純方位系統(tǒng)目標(biāo)運(yùn)動分析，其系統(tǒng)拓?fù)浣Y(jié)構(gòu)的復(fù)雜性遠(yuǎn)小于電力系統(tǒng)。文獻(xiàn)［67-68］研究機(jī)器人僅依靠方位量測的協(xié)同定位問題，定義了相對位置測量圖（Relative Position Measurement Graph，RPMG），結(jié)合非線性系統(tǒng)可觀測判據(jù)建立了RPMG性質(zhì)與可觀測性的關(guān)系，從而得出可觀測條件。該類方法在處理基于傳感網(wǎng)絡(luò)的BOTMA可觀測性分析上，提供了新的研究思路和技術(shù)途徑。

3 結(jié)束語

BOTMA的可觀測性問題本質(zhì)上屬于非線性系統(tǒng)解的存在唯一性問題，但早期的研究者多采用線性系統(tǒng)的分析方法，對量測方程進(jìn)行偽線性化處理。非線性系統(tǒng)理論的發(fā)展使可觀測性問題的研究有了新的理論依據(jù)，研究對象也不再局限于N階動態(tài)運(yùn)動目標(biāo)。盡管后來又有其他研究方法的出現(xiàn)，但線性和非線性系統(tǒng)方法仍然得到了最廣泛的應(yīng)用。20世紀(jì)八九十年代是可觀測性研究發(fā)展的重要時期，出現(xiàn)了大量經(jīng)典的研究成果，為后續(xù)的研究奠定了基礎(chǔ)。而后來的學(xué)者們主要從兩方面拓展豐富了前人的研究：一是尋找更為有效的方法，簡化分析過程，用最簡單的方法推導(dǎo)出已有的結(jié)論并將其融會貫通［69］；二是考慮實際應(yīng)用需求，對可觀測性的判據(jù)進(jìn)行改進(jìn)或是調(diào)整觀測條件［70-73］。

綜上所述，BOTMA可觀測性的研究重點(diǎn)表現(xiàn)在以下兩個方面：一是建立可觀測性判據(jù)，判斷所給系統(tǒng)是否可觀測或研究其可觀測程度；二是利用可觀測性判據(jù)，研究可觀測條件，并由此提出對觀測器的數(shù)目、運(yùn)動方式等方面的要求，以期改善系統(tǒng)可觀測性能。前者力求判據(jù)的簡單易用，后者則追求可觀測條件的直觀易懂。

從現(xiàn)階段可觀測性研究情況來看，可從以下幾方面對BOTMA可觀測性問題進(jìn)行更深入的研究：

（1）連續(xù)系統(tǒng)和離散系統(tǒng)的可觀測性研究方法不同，所得到的判據(jù)形式也不同，研究它們之間是否有類比或相互轉(zhuǎn)換的可能，以使可觀測理論更為完善，為實際應(yīng)用提供便利。

（2）一些可觀測條件是目標(biāo)和觀測器相對運(yùn)動狀態(tài)的表達(dá)式，而目標(biāo)與觀測器的相對運(yùn)動通常是未知的，需要更為直觀、實用的判據(jù)形式。諸如，能否建立直接基于方位量測數(shù)據(jù)，無須求導(dǎo)和求解觀測器與目標(biāo)的相對運(yùn)動參數(shù)的可觀測性判據(jù)?

（3）在有噪條件下，基于FIM的可觀測性判定方法僅適用于低階目標(biāo)，對于高階目標(biāo)會產(chǎn)生很大的計算量，而對于更一般性的運(yùn)動目標(biāo)，相應(yīng)的可觀測性分析理論和方法仍有待深入研究。

（4）純方位系統(tǒng)的可觀測性和觀測器與目標(biāo)的相對運(yùn)動形式密切相關(guān)。由于目標(biāo)運(yùn)動的未知性和不合作性，觀測器本身的合理有效機(jī)動是保證純方位系統(tǒng)可觀測的首要條件。從工程應(yīng)用角度出發(fā)，考慮觀測器最大機(jī)動速度、機(jī)動次數(shù)、傳感器探測盲區(qū)、戰(zhàn)術(shù)意圖、武器效能等實際因素制約下的觀測器機(jī)動模式，將后續(xù)目標(biāo)跟蹤估計算法與觀測器工程化機(jī)動軌跡研究緊密結(jié)合起來，從定性分析過渡至定量分析，尋求更具實用性和可操作性的可觀測性分析方法，將是未來該領(lǐng)域研究的重要課題。

（5）可觀測性研究的進(jìn)展總是與研究分析方法的革新密切相關(guān)，因此，進(jìn)行理論創(chuàng)新，尋求更為有效的研究途徑將是未來該領(lǐng)域研究的重要方面。

［1］劉忠，周豐，石章松，等.純方位目標(biāo)運(yùn)動分析［M］.北京：國防工業(yè)出版社，2009.

［2］潘泉，戴冠中.被動式跟蹤可觀測性分析的非線性系統(tǒng)方法［J］.信息與控制，1997，26（3）：168-173.

［3］Xu Z，Sheng A，Guo Z.The Analysis of the Unobservable Criteria on the Single Platform Bearings-only Target Tracking System［C］//IEEE International Conference on Networking，Sensing and Control，2008.

［4］權(quán)宏偉，彭冬亮，薛安克.純方位角目標(biāo)運(yùn)動分析的可觀測性研究［J］.火力與指揮控制，2009，34（10）：43-46.

［5］李洪瑞.水下目標(biāo)運(yùn)動分析關(guān)鍵技術(shù)研究［D］.南京：南京理工大學(xué)，2009.

［6］丁玉琴.幾類非線性系統(tǒng)觀測器設(shè)計研究［D］.濟(jì)南：山東大學(xué)，2007.

［7］Kolb R C，Hollister F H．Bearing-only Target Motion Estimation［C］//Proceedings of 1st Asilomar Conference on Circuits and Systems，1967：935-946．

［8］Murphy D J．Noise Bearings-only Target Motion Analysis［D］.Ph.D.Dissertation of Northwestern University，1970.

［9］Nardone S C，Aidala V J.Observability Criteria for Bearingsonly Target Motion Analysis［J］.IEEE Transactions on Aerospace and Electronic Systems，1981，17（2）：162-166.

［10］Hammel S E，Aidala V J.Observability Requirements for Three Dimensional Tracking Via Angle Measurements［J］. IEEE Transactions on Aerospace and Electronic Systems，1985，21（2）：200-206.

［11］Fogel E，Gavish M.Nth-order Dynamics Target Observability from Angle Measurements［J］.IEEE Transactions on Aerospace and Electronic Systems，1988，24（3）：305-308.

［12］Kalman R E.Contributions to the Theory of Optimal Control［J］.Bol.Soc.Mat.Mexicana，1960，5（2）：102-119.

［13］呂文亭，黃亮，王澍初，等.無人機(jī)單站無源定位中的可觀測性分析［J］.艦船電子工程，2013，32（11）：44-45.

［14］孫仲康，郭福成，馮道旺，等.單站無源定位跟蹤技術(shù)［M］.北京：國防工業(yè)出版社，2008.

［15］楊柏勝.多站被動傳感器定位技術(shù)研究［D］.西安：西安電子科技大學(xué)，2006.

［16］Gad A S，Mojica F，F(xiàn)arooq M.Geometric Approach to Target Tracking Motion Analysis in Bearing-only Tracking［C］//International Society for Optics and Photonics，2002. AeroSense 2002：1-12.

［17］Rao S K.Correspondence-Comments on Discrete-time Observability and Estimability Analysis for Bearings-only Target Motion Analysis［J］.IEEE Transactions on AerospaceandElectronicSystems，1998，34（4）：1361-1367.

［18］李先允.現(xiàn)代控制理論基礎(chǔ)［M］.北京：機(jī)械工業(yè)出版社，2007.

［19］Chen C T，王紀(jì)文.線性系統(tǒng)理論與設(shè)計［M］.北京：科學(xué)出版社，1988.

［20］鄧新蒲.單站無源定位可觀測性評述［J］.中國工程科學(xué)，2007，9（11）：54-62.

［21］Payne A N.Observability Conditions for Angles-only Tracking［C］//Twenty-Second Asilomar Conference on Signals，Systems and Computers.IEEE，1988.

［22］Becker K.Simple Linear Theory Approach to TMA Observability［J］.IEEE transactions on Aerospace and Electronic Systems，1993，29（2）：575-578.

［23］許兆鵬，韓樹平.多基陣純方位非機(jī)動目標(biāo)跟蹤可觀測性研究［J］.傳感器與微系統(tǒng)，2012，30（12）：57-59.

［24］Zhu W，Xu Z，Li B，et al.Research on the Observability of Bearings only Target Tracking Based on Multiple Sonar Sensors［C］//2012 Second International Conference on Instrumentation，Measurement，Computer，Communication and Control（IMCCC）.IEEE，2012：631-634.

［25］Zhu W，Xu Z，Li B，et al.Research on the Observability of Acceleration Based on Multiple Sonar Sensors［C］//2012 International Symposium on Information Science and Engineering（ISISE）.IEEE，2012：3-6.

［26］黃耀光，李建新，高博.基于空頻域信息的單站無源定位可觀測分析［J］.指揮控制與仿真，2012，34（5）：46-49.

［27］馮道旺，李宗華，周一宇，等.一種單站無源定位方法及其可觀測性分析［J］.國防科技大學(xué)學(xué)報，2004，26（1）：68-71.

［28］熊茂濤，張勇，吳欽章.基于純方位量測的三維運(yùn)動目標(biāo)可觀測性研究［J］.現(xiàn)代防御技術(shù)，2009，37（4）：30-34.

［29］Jauffret C，Pillon D.Observability in Passive Target Motion Analysis［J］.IEEE Transactions on Aerospace and Electronic Systems，1996，32（4）：1290-1300.

［30］Song T L.Observability of Target Tracking with Bearings-only Measurements［J］.IEEE Transactions on AerospaceandElectronicSystems，1996，32（4）：1468-1472.

［31］許志剛.TMA中不完全量測估計與觀測平臺機(jī)動模式研究［D］.南京：南京理工大學(xué)，2009.

［32］Lee T S，Dunn K P，Chang C B.On Observability and Unbiased Estimation of Nonlinear Systems［M］.System Modeling and Optimization.Springer Berlin Heidelberg，1982：258-266.

［33］Hermann R，Krener A J.Nonlinear Controllability and Observability［J］.IEEE Transactions on Automatic Control，1977，22（5）：728-740.

［34］喬梁.外輻射源角度測量定位算法可觀測性分析［J］.青島科技大學(xué)學(xué)報：自然科學(xué)版，2011，32（4）：426-428.

［35］曲毅，劉忠，彭鵬菲.基于垂直線陣無源定位的可觀測性研究［J］.艦船科學(xué)技術(shù)，2009，31（6）：91-94.

［36］喬梁.信源定位的可觀測性及跟蹤技術(shù)研究［D］.哈爾濱：哈爾濱工程大學(xué)，2010.

［37］Yi Y Z，Zhong K S.Observability Analysis of Single Passive Observer［C］//Proceedings of the IEEE 1995 National Aerospace andElectronics Conference.IEEE，1995.

［38］Hanba S.On the“Uniform”O(jiān)bservability of Discrete-time Nonlinear Systems［J］.IEEE Transactions on Automatic Control，2009，54（8）：1925-1928.

［39］馮紹軍，袁信.觀測度及其在卡爾曼濾波器設(shè)計中的應(yīng)用［J］.中國慣性技術(shù)學(xué)報，1999，7（2）：18-21.

［40］Yang B，Yao Y，He F.Passive Ranging Based on Observability Analysisand Receding Horizon Filter［J］.Tsinghua Science&Technology，2009，14（s2）：32-37.

［41］Chang X，Cui P，Cui H.Observability Analysis of Autonomous Navigation system of Interplanetary Spacecraft［C］//Chinese Control and Decision Conference，2009. IEEE，2009：766-770.

［42］Dong J L，Mo B.The Method of System Observability Analysis Using Pseudo-inverse of System Observability Matrix［C］//2013 32nd Chinese Control Conference（CCC）. IEEE，2013：55-59.

［43］劉準(zhǔn)，陳哲.條件數(shù)在系統(tǒng)可觀測性分析中的應(yīng)用研究［J］.系統(tǒng)仿真學(xué)報，2004，16（7）：1552-1555.

［44］常曉華，崔平遠(yuǎn)，王曉明，等.基于條件數(shù)的能觀性度量方法及在自主導(dǎo)航系統(tǒng)中的應(yīng)用［J］.宇航學(xué)報，2010，31（5）：1331-1337.

［45］Goshen-meskin D，Bar-itzhack I.Observability Analysis of Piece-wise Constant Systems.I.Theory［J］.IEEE Transactions on Aerospace and Electronic Systems，1992，28（4）：1056-1067.

［46］Kong X，Guo M，Dong J.An improved PWCS Approach on Observabilityanalysis of Linear Time-varying System［C］//Proceedings of the 21st Annual International Conference on Chinese Control and Decision Conference.IEEE Press，2009：816-820.

［47］Le Cadre J P，Jauffret C.Discrete-time Observability and Estimability Analysis for Bearings-only Target Motion Analysis［J］.IEEE Transactions on Aerospace and Electronic Systems，1997，33（1）：178-201.

［48］Jauffret C，Pillon D，Pignol A C.Bearings-only Maneuvering Target Motion Analysis from a Nonmaneuvering Platform［J］.IEEE Transactions on Aerospace and Electronic Systems，2010，46（4）：1934-1949.

［49］Jauffret C，Pillon D，Pignol A C.Bearings-only TMA without Observer Maneuver［C］//2008 11th International Conference on Information Fusion.IEEE，2008：1-8.

［50］Jauffret C，Pillon D，Pignol A C.Leg-by-leg Bearings-only TMA without Observer Maneuver［C］//Proceedings of the International Conference on Information Fusion.Seattle，USA：2009.

［51］Clavard J，Pillon D，Pignol A C，et al.Bearings-only Target Motion Analysis of a Source in a Circular Constant Speed Motion from A Non-maneuvering Platform［C］//Proceedings of the 14th International Conference on Information Fusion（FUSION）.2011：1-8.

［52］Clavard J，Pillon D，Pignol A C，et al.Target Motion Analysis of a Source in a Constant Turn from a Nonmaneuvering Observer［J］.IEEE Transactions on Aerospace and Electronic Systems，2013，49（3）：1760-1780.

［53］吳京，王博，李駿，等.空間目標(biāo)天基光學(xué)跟蹤可觀測性分析［J］.湖南大學(xué)學(xué)報：自然科學(xué)版，2009，36（11）：32-37.

［54］Fang X，Yan W，Wang Y，et al.Observability Analysis of Underwater Localization Based on Fisher Information Matrix［C］//2013 32nd Chinese Control Conference（CCC）. IEEE，2013：4866-4869.

［55］Jauffret C.Observability and Fisher Information Matrix in NonlinearRegression［J］.IEEETransactionson Aerospace and Electronic Systems，2007，43（2）：756-759.

［56］Pignol A C，Jauffret C，Pillon D.Properties of Range-only Target Motion Analysis［C］//2013 16th International Conference on Information Fusion（FUSION）.IEEE，2013：1693-1698.

［57］夏佩倫.純方位目標(biāo)跟蹤系統(tǒng)的可觀測程度［J］.火力與指揮控制，1992，17（2）：23-29.

［58］孫仲康，周一宇，何黎星，等.單多基地有源無源定位技術(shù)［M］.北京：國防工業(yè)出版社，1996.

［59］夏佩倫.一般僅方位目標(biāo)跟蹤系統(tǒng)可觀測的必要條件［J］.電光與控制，2011，18（6）：1-4.

［60］李洪瑞，盛安冬.離散純方位系統(tǒng)的可觀測性判據(jù)［J］.控制理論與應(yīng)用，2009，26（5）：570-572.

［61］李洪瑞，盛安冬.連續(xù)純方位系統(tǒng)的可觀測性分析［J］.兵工學(xué)報，2009，30（11）：1446-1450.

［62］董志榮.純方位系統(tǒng)定位與跟蹤的本載體最優(yōu)軌線方程及其最優(yōu)軌線［J］.指揮控制與仿真，2007，29（1）：7-15.

［63］石章松，劉忠.單站純方位目標(biāo)跟蹤系統(tǒng)可觀測性分析［J］火力與指揮控制，2007，32（2）：26-29.

［64］李長文，趙建昕.聲速有限條件下的純方位方法［J］.指揮控制與仿真，2012，34（5）：94-98.

［65］Parlangeli G N G.Controllability and Observability of Grid Graphs Via Reduction and Symmetries［J］.arXiv Preprint arXiv，2012（129）：1203-1204.

［66］Doostmohammadian M，Khan U A.Topology Design in Net-worked Estimation：A Generic Approach［C］//American Control Conference（ACC），2013.

［67］Sharma R，Beard R W，Taylor C N，et al.Graph-based Observability Analysis of Bearing-only Cooperative Localization［J］.IEEE Transactions on Robotics，2012，28（2）：522-529.

［68］Sharma R，Quebe S，Beard R W，et al.Bearing-only Cooperative Localization［J］.Journal of Intelligent&Robotic Systems，2013，72（3-4）：429-440.

［69］Li H，Chen L，Sheng A.Two-step Modeling and Observability for Bearings-only Tracking［C］//International Conference on Mechatronics and Automation，2009.

［70］劉忠，鄧聚龍.多傳感器系統(tǒng)純方位定位與可觀測性分析［J］.火力與指揮控制，2004，29（5）：79-83.

［71］劉忠.單靜止站純方位系統(tǒng)的可觀測性與目標(biāo)參數(shù)估計算法［J］.海軍工程大學(xué)學(xué)報，2004，13（5）：20-23.

［72］劉忠，鄧聚龍.等速運(yùn)動觀測站純方位系統(tǒng)的可觀測性［J］.火力與指揮控制，2004，29（6）：51-54.

［73］Ferdowsi M H.Observability Conditions for Target States with Bearing-only Measurements in Three-dimensional Case［C］//Computer Aided Control System Design，2006 IEEE International Conference on Control Applications，2006 IEEE International Symposium on Intelligent Control. IEEE，2006.

Overview of Observability of Bearings-only Target Motion Analysis

YANG Jing，LI Yin-ya，QI Guo-qing，SHENG An-dong
（Nanjing University of Science and Technology，Nanjing210014，China）

Bearings-Only Target Motion Analysis（BOTMA）is an efficient method for passive localization and tracking due to the superiority of which uses only bearing measurements to estimate parameters of target motion.The researches of bearings-only systems concentrate on three aspects：observability conditions，target tracking and estimation strategies and optimal observer maneuver trajectories.The observability problem is the key issue that is to be handled firstly because the corresponding observability conditions are the basis and precondition of follow-up researches of target localization and tracking.First，the basic concept of observability in bearings-only system is introduced. Then from the research methodology perspective，the existing BOTMA observability research results，including geometric methods，linear and nonlinear system based methods and numerical methods，etc.，are summarized and reviewed.Finally，several new opinions on future researches in this area are presented.

observability，BOTMA，passive tracking

TP271

A

1002-0640（2015）12-0001-08

2014-11-05

2015-01-15

國家自然科學(xué)基金（61273076，61104186）；江蘇省自然科學(xué)基金資助項目（BK2012801）

楊婧（1990-），女，山西運(yùn)城人，碩士研究生。研究方向：水下純方位系統(tǒng)的可觀測性。