(57)【要約】  追跡の測定値を極座標系から直交座標系に変換することによって発生するレーダ測定値の測定誤差共分散の変化を基に、α−β追跡フィルタのα及びβフィルタ利得を求めるための追跡指数を計算する追跡指数部、及び追跡指数変化率に対するα及びβフィルタ利得の変化率の偏微分係数と追跡指数変化量とを基に、α及びβフィルタ利得の変化量を計算し、α及びβフィルタ利得を更新するフィルタ利得更新部を具備するα−β追跡フィルタを有するレーダ追跡装置である。
2010505111.tif 000002
【特許請求の範囲】
【請求項１】
α−β追跡フィルタを有するレーダ追跡装置において、
前記α−β追跡フィルタは、
追跡の測定値を極座標系から直交座標系に変換することによって発生するレーダ測定値の測定誤差共分散の変化を基に、前記α−β追跡フィルタのα及びβフィルタ利得を求めるための追跡指数を計算する追跡指数部と、
前記追跡指数変化率に対する前記α及びβフィルタ利得の変化率の偏微分係数と前記追跡指数変化量とを基に、前記α及びβフィルタ利得の変化量を計算し、前記α及びβフィルタ利得を更新するフィルタ利得更新部と、
を具備することを特徴とするレーダ追跡装置。
【請求項２】
前記レーダ測定値の測定誤差共分散の変化は、レーダ測定誤差の標準偏差の関数を基に示されることを特徴とする請求項１に記載のレーダ追跡装置。
【請求項３】
前記追跡指数は、
レーダ測定誤差の標準偏差に対する前記追跡の運動モデルに係る工程誤差の標準偏差の比例値と、離散化間隔とに比例することを特徴とする請求項１に記載のレーダ追跡装置。
【請求項４】
2010505111.tif 000091
【請求項５】
α−β追跡フィルタを利用するレーダ追跡方法において、
追跡の測定値を極座標系から直交座標系に変換することによって発生するレーダ測定値の測定誤差共分散の変化を基に、前記α−β追跡フィルタのα及びβフィルタ利得を求めるための追跡指数を計算する追跡指数の計算段階と、
前記追跡指数変化率に対する前記α及びβフィルタ利得の変化率の偏微分係数と前記追跡指数変化量とを基に、前記α及びβフィルタ利得の変化量を計算し、前記α及びβフィルタ利得を更新するフィルタ利得の更新段階と、
を含むことを特徴とするレーダ追跡方法。
【請求項６】
前記レーダ測定値の測定誤差共分散の変化は、レーダ測定誤差の標準偏差の関数を基に示されることを特徴とする請求項５に記載のレーダ追跡方法。
【請求項７】
請求項５または請求項６に記載の方法をコンピュータで実行させるためのプログラムを記録したコンピュータで読み取り可能な記録媒体。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、レーダ追跡装置及びその方法に係り、さらに詳細には、測定値の誤差特性変化を考慮してα−β利得を更新するレーダ追跡装置及びその方法に関する。
【０００２】
すなわち本発明は、直交座標系に変換されたレーダ測定値から測定誤差共分散を計算し、これを利用して追跡指数の変化量、α−β利得変化量計算及びα−β利得を更新するレーダ追跡装置及びその方法に関する。
【背景技術】
【０００３】
レーダを利用した追跡装置は、軍事用目的だけではなく、港湾及び大型船舶などで主に使われてきた。しかし、最近には、車両の衝突防止分野にまでその利用分野が拡大されている。
【０００４】
レーダ追跡装置には、α−βフィルタ、カルマン（Kalman）フィルタなどが使われており、特に、数十〜数百の標的を同時に追跡しなければならないＴＷＳ（Track−While−Scan）追跡装置では、α−βフィルタが広範囲に使われている。
【０００５】
α−βフィルタは、他の追跡方法に比べて演算量が非常に少ないために、リアルタイムで多数の標的を追跡せねばならないＴＷＳ追跡装置に適するためである。しかし、直交座標系に変換されたレーダ測定値は、標的の距離及び方位によって測定誤差の特性が変わるので、これを考慮したα−β利得の更新が必要である。
【発明の概要】
【発明が解決しようとする課題】
【０００６】
本発明がなそうとする技術的課題は、レーダ測定値を利用して物標に対する追跡を行う場合、レーダ測定値の誤差特性変化を考慮して追跡フィルタの利得を更新することによって、さらに向上した追跡性能を提供するところにある。
【課題を解決するための手段】
【０００７】
前記技術的課題をなすための本発明によるレーダ追跡装置の一実施例は、追跡の測定値を極座標系から直交座標系に変換することによって発生するレーダ測定値の測定誤差共分散の変化を基に、α−β追跡フィルタのα及びβフィルタ利得を求めるための追跡指数を計算する追跡指数部、及び追跡指数変化率に対するα及びβフィルタ利得の変化率の偏微分係数と追跡指数変化量とを基にα及びβフィルタ利得の変化量を計算し、α及びβフィルタ利得を更新するフィルタ利得更新部を具備するα−β追跡フィルタを有する。
【０００８】
前記技術的課題をなすための本発明によるレーダ追跡装置の一実施例は、追跡の測定値を極座標系から直交座標系に変換することによって発生するレーダ測定値の測定誤差共分散の変化を基に、α−β追跡フィルタのα及びβフィルタ利得を求めるための追跡指数を計算する追跡指数の計算段階、及び追跡指数変化率に対するα及びβフィルタ利得の変化率の偏微分係数と追跡指数変化量とを基にα及びβフィルタ利得の変化量を計算し、α及びβフィルタ利得を更新するフィルタ利得の更新段階を有する。
【発明の効果】
【０００９】
本発明は、レーダ測定値を利用した追跡装置において、測定誤差の特性変化によって追跡フィルタの利得を更新することによって、向上した追跡性能を提供する。
【００１０】
このように本発明は、軍事用レーダ；港湾などの船舶の航路誘導及び統制のためのレーダ；航空機、船舶及び車両などの衝突防止のためのレーダのような多様な目的のレーダに使われうる。
【図面の簡単な説明】
【００１１】
【図１】本発明の望ましい一実施例によるレーダ追跡装置を示す図面である。
【図２】本発明の望ましい一実施例によるレーダ追跡装置で使用する追跡フィルタ部を示す図面である。
【図３】本発明の望ましい一実施例によるレーダ追跡方法のフローチャートである。
【図４】本発明の他の望ましい一実施例によるレーダ追跡方法のフローチャートである。
【図５】本発明の望ましい一実施例によるレーダ測定誤差特性の変化による追跡指数の変化を示すグラフである。
【図６】本発明の望ましい一実施例による更新されたα利得パラメータを示すグラフである。
【図７】本発明の望ましい一実施例による更新されたβ利得パラメータを示すグラフである。
【図８】本発明の望ましい一実施例による更新されたα−β利得を利用した追跡フィルタの位置推定誤差を示すグラフである。
【図９】本発明の望ましい一実施例による更新されたα−β利得を利用した追跡フィルタの速度推定誤差を示すグラフである。
【発明を実施するための形態】
【００１２】
以下、添付した図面を使用しつつ、本発明によるレーダ追跡装置及びその方法について詳細に説明する。
【００１３】
図１は、本発明の望ましい一実施例によるレーダ追跡装置である。図１を参照すれば、本発明によるレーダ追跡装置は、座標変換部１０１、追跡−物標連繋部１０２、追跡管理部１０３、α−β追跡フィルタ１０４及びゲーティング部１０５を備える。
【００１４】
座標変換部１０１は、極座標系で観測された物標の位置を、直交座標系での位置に変換する。レーダ観測値は極座標系で得られるので、追跡フィルタで使用する座標系と同じ座標系（直交座標系）への変換が必要である。
【００１５】
追跡−物標連繋部１０２は、探知された物標の位置と予測された標的の位置とを利用し、距離誤差の自乗が最小になる標的と物標とを連繋する。物標は、レーダスキャン映像から抽出された「object」である。一般的に物標という用語は、軍隊で使用する表現であり、軍隊で物標は、レーダによって探知された未確認物体を意味する。
【００１６】
標的は、物標のうちで追跡しなければならない価値のある関心物標をいう。主に、任務遂行中である味方の船舶、飛行機や操業中である船舶または威嚇要因になると推定される物標が標的である。
【００１７】
標的を追跡して標的の現在位置及び航速方向が分かる。これを介して味方の円滑な任務遂行のための誘導（guidance）、操業中である船舶の場合には、軍事境界線侵犯、操業区域離脱または離脱可能性などに対して判断できる。
【００１８】
すなわち、船舶などの誘導を行うことができ、多数の船舶が密集した区域では、衝突回避などのための自動衝突予防補助装置（ＡＲＰＡ：Automatic Radar Plotting Aids）機能に追跡情報を使用する。
【００１９】
追跡管理部１０３は、新たに探知された物標に対する追跡の生成を行い、物標と関連しない追跡に対して追跡の維持及び削除を管理する。
【００２０】
α−β追跡フィルタ１０４は、物標と関連した追跡に対して、物標を測定値として追跡の更新及び予測を行う。α−β追跡フィルタ１０４は、既追跡物標の予測位置提供のために追跡予測を行う。
【００２１】
ゲーティング部１０５は、予測された標的の位置とレーダによって観測された物標との連繋のための境界を設定する。
【００２２】
図２は、本発明の望ましい一実施例によるレーダ追跡装置で使用する追跡フィルタ部である。図２を参照すれば、本発明のレーダ追跡装置で使用するα−β追跡フィルタ２００は、追跡指数部２１０及びフィルタ利得更新部２２０を備える。
【００２３】
追跡指数部２１０は、追跡の測定値を極座標系から直交座標系に変換することによって発生するレーダ測定値の測定誤差共分散の変化を基に、α−β追跡フィルタのα及びβフィルタ利得を求めるための追跡指数を計算する。
【００２４】
追跡指数部２１０でのレーダ測定値の測定誤差共分散の変化は、レーダ測定誤差の標準偏差の関数を基に示す。これについて詳細に説明すれば、次の通りである。
【００２５】
α−β追跡フィルタは、レーダによって観測された物標の位置情報から測定誤差共分散を計算する。
【００２６】
まず、最適α−β追跡フィルタの運動及び観測モデルは、次の通りである。離散時間で位置と速度とを状態変数とする標的の運動モデルは、次の通り示すことができる。
【００２７】
【数１】
2010505111.tif 000003
【００２８】
2010505111.tif 000004
【００２９】
Ｔ秒の間隔で観測される標的の位置に係る測定方程式は、次の通り示すことができる。
【００３０】
【数２】
2010505111.tif 000005
【００３１】
2010505111.tif 000006
【００３２】
α−β追跡フィルタの予測式（式３）と改善式（式４）は、次の通り示すことができる。
【００３３】
【数３】
2010505111.tif 000007
【００３４】
【数４】
2010505111.tif 000008
【００３５】
2010505111.tif 000009
【００３６】
【数５】
2010505111.tif 000010
【００３７】
最適のα−β利得は、追跡指数から次の通り求めることができる。
【００３８】
【数６】
2010505111.tif 000011
【００３９】
【数７】
2010505111.tif 000012
【００４０】
従って、α−β追跡フィルタの利得行列は、次の通り示すことができる。
【００４１】
【数８】
2010505111.tif 000013
【００４２】
また、変換された測定値を利用する直交座標系で設計された分離型追跡フィルタの利得行列は、距離と方位とを分けて次の通り示すことができる。
【００４３】
【数９】
2010505111.tif 000014
【００４４】
ここで、ＤとＤ2は、（式１０）及び（式１１）の通りである。
【００４５】
【数１０】
2010505111.tif 000015
【００４６】
【数１１】
2010505111.tif 000016
【００４７】
2010505111.tif 000017
【００４８】
ＭＵＣＭ（Modified Unbiased Converted Measurements）を利用して直交座標系に変換されたレーダ測定値に対する測定誤差の共分散行列は、次の通り示すことができる。
【００４９】
【数１２】
2010505111.tif 000018
【００５０】
ここで、
【００５１】
【数１３】
2010505111.tif 000019
【００５２】
【数１４】
2010505111.tif 000020
【００５３】
【数１５】
2010505111.tif 000021
【００５４】
【数１６】
2010505111.tif 000022
【００５５】
【数１７】
2010505111.tif 000023
【００５６】
2010505111.tif 000024
【００５７】
（式１２）と（式１３）とから分かるように、ｘ，ｙ軸に対する測定誤差共分散は、標的の距離によって変わるので、あらかじめ計算された追跡指数及びα−β利得が不適切になって追跡性能が低下しうる。従って、これを考慮して、追跡フィルタの利得を更新せねばならない。
【００５８】
これまで、ＭＵＣＭを利用して直交座標系に変換されたレーダ測定値に対する測定誤差の共分散行列を示した。
【００５９】
しかし、変換方法をＭＵＣＭのみに制限する必要はなく、ＭＵＣＭだけではなく、古典的な変換（Classical Conversion）、Debiased Conversion、不偏変換（Unbiased Conversion）を含むいなかる方法で変換しても、直交座標系に変換した測定値の測定誤差共分散は、変わるようになる。
【００６０】
次に、追跡指数の変化量を計算する。近似的α−β利得更新のために、追跡指数の変化量を計算する。
【００６１】
まず、ｘ，ｙ軸に対する測定誤差共分散の変化を考慮し、（式５）を次の通り示すことができる。
【００６２】
【数１８】
2010505111.tif 000025
【００６３】
2010505111.tif 000026
【００６４】
フィルタ利得更新部２２０は、追跡指数変化率に対するα及びβフィルタ利得の変化率の偏微分係数と追跡指数変化量とを基に、α及びβフィルタ利得の変化量を計算し、α及びβフィルタ利得を更新する。
【００６５】
追跡指数での偏微分係数を利用し、δα，δβを計算する。追跡指数の変化率によるα−β利得の変化率は、次の通り示すことができる。
【００６６】
【数１９】
2010505111.tif 000027
【００６７】
【数２０】
2010505111.tif 000028
【００６８】
追跡指数に対するα−β利得の偏微分方程式は、（式６）と（式７）とから、次の通り誘導できる。
【００６９】
【数２１】
2010505111.tif 000029
【００７０】
【数２２】
2010505111.tif 000030
【００７１】
ｋ番目スキャンとｋ＋１番目スキャンとの間の追跡指数と、α−β利得の変化量は、次の通り示すことができる。
【００７２】
【数２３】
2010505111.tif 000031
【００７３】
【数２４】
2010505111.tif 000032
【００７４】
【数２５】
2010505111.tif 000033
【００７５】
（式１９）ないし（式２５）から、α−β利得の変化量は、次の通り与えられる。
【００７６】
【数２６】
2010505111.tif 000034
【００７７】
【数２７】
2010505111.tif 000035
【００７８】
従って、ｋ＋１番目の更新周期で、公称点に対する追跡指数の変化量と、α−β利得の変化量とを利用した更新されたα−β利得は、次の通り示すことができる。
【００７９】
【数２８】
2010505111.tif 000036
【００８０】
【数２９】
2010505111.tif 000037
【００８１】
【数３０】
2010505111.tif 000038
【００８２】
また、α−β追跡フィルタの利得行列は、次の通り与えられる。
【００８３】
【数３１】
2010505111.tif 000039
【００８４】
（式１８）ないし（式３１）まで追跡フィルタの利得を計算する過程は、追跡指数に対するα−β利得の偏微分が定数であると仮定すれば、非常に簡略に整理でき、（式２６）と（式２７）は、次の通り示すことができる。
【００８５】
【数３２】
2010505111.tif 000040
【００８６】
【数３３】
2010505111.tif 000041
【００８７】
2010505111.tif 000042
【００８８】
【数３４】
2010505111.tif 000043
【００８９】
【数３５】
2010505111.tif 000044
【００９０】
2010505111.tif 000045
【００９１】
計算されたα−β利得の変化量を利用し、α−β利得を更新する。計算されたα−β利得の変化量を利用し、次の通りα−β利得を更新する。また、更新されたα−β利得を利用し、追跡を更新する。
【００９２】
図３は、本発明の望ましい一実施例によるレーダ追跡方法のフローチャートである。図３を参照すれば、α−β追跡フィルタを利用するレーダ追跡方法は、追跡指数の計算段階（Ｓ３０１）及びフィルタ利得の更新段階（Ｓ３０２）を含む。
【００９３】
追跡指数の計算段階（Ｓ３０１）は、追跡の測定値を極座標系から直交座標系に変換することによって発生するレーダ測定値の測定誤差共分散の変化を基に、α−β追跡フィルタのα及びβフィルタ利得を求めるための追跡指数を計算する。
【００９４】
追跡指数の計算段階（Ｓ３０１）で、レーダ測定値の測定誤差共分散の変化は、レーダ測定誤差の標準偏差の関数を基に示すことができる。
【００９５】
フィルタ利得の更新段階（Ｓ３０２）は、追跡指数変化率に対するα及びβフィルタ利得の変化率の偏微分係数と追跡指数変化量とを基に、α及びβフィルタ利得の変化量を計算し、α及びβフィルタ利得を更新する。
【００９６】
図４は、本発明の他の望ましい一実施例によるレーダ追跡方法で、測定値誤差特性の変化を考慮した追跡フィルタの利得更新方法に係る説明である。
【００９７】
α−β追跡フィルタは、レーダによって観測された物標の位置情報から測定誤差共分散を計算する。
【００９８】
まず、最適α−β追跡フィルタの運動及び観測モデルは、次の通りである。離散時間で、位置と速度とを状態変数とする標的の運動モデルは、次の通り示すことができる。
【００９９】
【数３６】
2010505111.tif 000046
【０１００】
2010505111.tif 000047
【０１０１】
Ｔ秒の間隔で観測される標的の位置に係る測定方程式は、次の通り示すことができる。
【０１０２】
【数３７】
2010505111.tif 000048
【０１０３】
2010505111.tif 000049
【０１０４】
α−β追跡フィルタの予測式（式３）と改善式（式４）は、次の通り示すことができる。
【０１０５】
【数３８】
2010505111.tif 000050
【０１０６】
【数３９】
2010505111.tif 000051
【０１０７】
2010505111.tif 000052
【０１０８】
【数４０】
2010505111.tif 000053
【０１０９】
最適のα−β利得は、追跡指数から次の通り求めることができる。
【０１１０】
【数４１】
2010505111.tif 000054
【０１１１】
【数４２】
2010505111.tif 000055
【０１１２】
従って、α−β追跡フィルタの利得行列は、次の通り示すことができる。
【０１１３】
【数４３】
2010505111.tif 000056
【０１１４】
また、変換された測定値を利用する直交座標系で設計された分離型追跡フィルタの利得行列は、距離と方位とを分けて次の通り示すことができる。
【０１１５】
【数４４】
2010505111.tif 000057
【０１１６】
ここで、ＤとＤ2は、（式１０）及び（式１１）の通りである。
【０１１７】
【数４５】
2010505111.tif 000058
【０１１８】
【数４６】
2010505111.tif 000059
【０１１９】
である。
【０１２０】
2010505111.tif 000060
【０１２１】
ＭＵＣＭを利用し、直交座標系に変換されたレーダ測定値に対する測定誤差の共分散行列は、次の通り示すことができる。
【０１２２】
【数４７】
2010505111.tif 000061
【０１２３】
ここで、
【０１２４】
【数４８】
2010505111.tif 000062
【０１２５】
【数４９】
2010505111.tif 000063
【０１２６】
【数５０】
2010505111.tif 000064
【０１２７】
【数５１】
2010505111.tif 000065
【０１２８】
【数５２】
2010505111.tif 000066
【０１２９】
2010505111.tif 000067
【０１３０】
（式１２）と（式１３）とから分かるように、ｘ，ｙ軸に対する測定誤差共分散は、標的の距離によって変わるので、あらかじめ計算された追跡指数及びα−β利得が不適切になって追跡性能が低下しうる。従って、これを考慮して、追跡フィルタの利得を更新せねばならない。
【０１３１】
これまで、ＭＵＣＭを利用して直交座標系に変換されたレーダ測定値に対する測定誤差の共分散行列を示した。しかし、変換方法をＭＵＣＭだけに制限する必要はなく、ＭＵＣＭだけではなく、古典的な変換、Debiased Conversion）、不偏変換を含むいかなる方法で変換しても、直交座標系に変換した測定値の測定誤差共分散は変わるようになる。
【０１３２】
次に、追跡指数の変化量を計算する（Ｓ４０２）。近似的α−β利得更新のために、追跡指数の変化量を計算する。
【０１３３】
まず、ｘ，ｙ軸に対する測定誤差共分散の変化を考慮し、（式５）を次の通り示すことができる。
【０１３４】
【数５３】
2010505111.tif 000068
【０１３５】
2010505111.tif 000069
【０１３６】
追跡指数での偏微分係数を利用し、δα，δβを計算する（Ｓ４０３）。追跡指数の変化率によるα−β利得の変化率は、次の通り示すことができる。
【０１３７】
【数５４】
2010505111.tif 000070
【０１３８】
【数５５】
2010505111.tif 000071
【０１３９】
追跡指数に対するα−β利得の偏微分方程式は、（式６）と（式７）とから、次の通り誘導できる。
【０１４０】
【数５６】
2010505111.tif 000072
【０１４１】
【数５７】
2010505111.tif 000073
【０１４２】
ｋ番目スキャンとｋ＋１番目スキャンとの間の追跡指数と、α−β利得の変化量は、次の通り示すことができる。
【０１４３】
【数５８】
2010505111.tif 000074
【０１４４】
【数５９】
2010505111.tif 000075
【０１４５】
【数６０】
2010505111.tif 000076
【０１４６】
（式１９）ないし（式２５）から、α−β利得の変化量は、次の通り与えられる。
【０１４７】
【数６１】
2010505111.tif 000077
【０１４８】
【数６２】
2010505111.tif 000078
【０１４９】
従って、ｋ＋１番目の更新周期で、公称点に対する追跡指数の変化量と、α−β利得の変化量とを利用した更新されたα−β利得は、次の通り示すことができる。
【０１５０】
【数６３】
2010505111.tif 000079
【０１５１】
【数６４】
2010505111.tif 000080
【０１５２】
【数６５】
2010505111.tif 000081
【０１５３】
また、α−β追跡フィルタの利得行列は、次の通り与えられる。
【０１５４】
【数６６】
2010505111.tif 000082
【０１５５】
（式１８）ないし（式３１）まで追跡フィルタの利得を計算する過程は、追跡指数に対するα−β利得の偏微分が定数であると仮定すれば、非常に簡略に整理でき、（式２６）と（式２７）は、次の通り示すことができる。
【０１５６】
【数６７】
2010505111.tif 000083
【０１５７】
【数６８】
2010505111.tif 000084
【０１５８】
2010505111.tif 000085
【０１５９】
【数６９】
2010505111.tif 000086
【０１６０】
【数７０】
2010505111.tif 000087
【０１６１】
2010505111.tif 000088
【０１６２】
計算されたα−β利得の変化量を利用し、α−β利得を更新する（Ｓ４０４）。計算されたα−β利得の変化量を利用し、次の通りα−β利得を更新する。また、更新されたα−β利得を利用して追跡を更新する。
【０１６３】
図５は、本発明の望ましい一実施例によるレーダ測定誤差特性の変化による追跡指数の変化を示す。すなわち、図５は、観測誤差共分散の変化に対する追跡指数の変化量を示したものである。
【０１６４】
標的の運動モデルは、等速運動モデルを利用し、工程誤差の標準偏差は、０．９８ｍ／ｓ2に設定した。レーダ観測値は、距離と方位とに対してそれぞれ５０ｍと１°との標準偏差を有する平均が０である白色ノイズと仮定して生成した。また、更新周期は５秒、すなわち、アンテナの回転率を１２ｒｐｍ（revolution per minute）に設定した。
【０１６５】
図６は、本発明の望ましい一実施例による更新されたα利得パラメータを示す。図６を参照すれば、既存方法による公称α利得パラメータと、本発明の一実施例による更新されたα利得パラメータ及び最適α利得パラメータとを示している。
【０１６６】
図６は、最適追跡指数を利用して計算したα利得と、近似法を利用して計算された追跡指数を利用したα利得とを示したものであり、後者は、本発明によるα−β追跡フィルタに使われる。
【０１６７】
円は、公称点での最適α利得を意味する。標的の運動モデルは等速運動モデルを利用し、工程誤差の標準偏差は、０．９８に設定した。レーダ観測値は、距離と方位とに対してそれぞれ５０ｍと１°との標準偏差を有する平均が０である白色ノイズと仮定して生成した。また、更新周期は５秒、すなわち、アンテナの回転率を１２ｒｐｍに設定した。
【０１６８】
2010505111.tif 000089
【０１６９】
図７は、本発明の望ましい一実施例による更新されたβ利得パラメータを示す。図７を参照すれば、既存方法による公称β利得パラメータと、本発明の一実施例による更新されたβ利得パラメータ及び最適β利得パラメータとを示している。
【０１７０】
図７は、最適追跡指数を利用して計算したβ利得と、近似法を利用して計算された追跡指数を利用したβ利得とを示したものであり、後者は、本発明によるα−β追跡フィルタに使われる。
【０１７１】
円は、公称点での最適β利得を意味する。標的の運動モデルは等速運動モデルを利用し、工程誤差の標準偏差は、０．９８に設定した。レーダ観測値は、距離と方位とに対してそれぞれ５０ｍと１°との標準偏差を有する平均が０である白色ノイズと仮定して生成した。また、更新周期は５秒、すなわち、アンテナの回転率を１２ｒｐｍに設定した。
【０１７２】
2010505111.tif 000090
【０１７３】
図８は、本発明の望ましい一実施例による更新されたα−β利得を利用した追跡フィルタの位置推定誤差を示す。図８を参照し、既存方法と本発明とによる更新されたα−β利得を利用した追跡フィルタの位置推定誤差に対するＲＭＳ値の比較する。
【０１７４】
図８は、位置と速度との推定値に対するＲＭＳ（Root Mean Square）誤差を示したものであり、実線は、一般的なα−β追跡フィルタの位置推定誤差に対するＲＭＳ値であり、点線は、提案した方式によるα−β追跡フィルタの位置推定誤差に対するＲＭＳ値である。
【０１７５】
本発明によるα−β追跡フィルタが、一般的なα−β追跡フィルタに比べて大きく改善された位置追跡性能を提供することが分かる。
【０１７６】
点線は、本発明によるα−β利得更新を利用したα−β追跡フィルタを利用した標的の位置推定誤差をＲＭＳ値で示したものであり、実線は、既存方法である公称点でのα−β利得を利用して追跡した場合の標的の位置推定誤差ＲＭＳ値で示したものである。
【０１７７】
図９は、本発明の望ましい一実施例による更新されたα−β利得を利用した追跡フィルタの速度推定誤差を示す。図９を参照し、既存方法と本発明とによる更新されたα−β利得を利用した追跡フィルタの速度推定誤差に対するＲＭＳ値を比較する。
【０１７８】
図９は、位置と速度との推定値に対するＲＭＳ誤差を示したものであり、実線は、一般的なα−β追跡フィルタの位置推定誤差に対するＲＭＳ値であり、点線は、提案した方式によるα−β追跡フィルタの位置推定誤差に対するＲＭＳ値である。
【０１７９】
点線は、本発明によるα−β利得更新を利用したα−β追跡フィルタを利用した標的の位置推定誤差をＲＭＳ値で示したものであり、実線は、既存方法である公称点でのα−β利得を利用して追跡した場合の標的の位置推定誤差ＲＭＳ値で示したものである。
【０１８０】
本発明によるレーダ観測誤差特性変化を考慮したα−β追跡フィルタが、一般的なα−β追跡フィルタに比べて大きく改善された位置追跡性能を提供することが分かる。
【０１８１】
本発明はまた、コンピュータで読み取り可能な記録媒体にコンピュータで読み取り可能なコードでもって具現することが可能である。コンピュータで読み取り可能な記録媒体は、コンピュータシステムによって読み取り可能なデータが保存されるあらゆる種類の記録装置を含む。
【０１８２】
コンピュータで読み取り可能な記録媒体の例としては、ＲＯＭ（Read-Only Memory）、ＲＡＭ（Random-Access Memory）、ＣＤ−ＲＯＭ、磁気テープ、フロッピーディスク、光データ保存装置などがあり、またキャリアウェーブ（例えば、インターネットを介した伝送）の形態で具現されるものも含む。また、コンピュータで読み取り可能な記録媒体は、ネットワークに連結されたコンピュータシステムに分散され、分散方式でコンピュータで読み取り可能なコードが保存されて実行されうる。
【０１８３】
以上、本発明についてその望ましい実施例を中心に説明した。本発明が属する技術分野で当業者は、本発明が本発明の本質的な特性から外れない範囲で変形された形態で具現しうることを理解することができるであろう。従って、開示された実施例は、限定的な観点ではなく、説明的な観点から考慮されねばならない。本発明の範囲は、前述の説明ではなくして特許請求の範囲に示されており、それと同等な範囲内にあるあらゆる差異点は、本発明に含まれるものと解釈されるものである。
【産業上の利用可能性】
【０１８４】
本発明は、レーダ測定値を利用して物標に対する追跡を行う場合、レーダ測定値の誤差特性変化を考慮して追跡フィルタの利得を更新することによって、さらに向上した追跡性能を提供する。
【図１】
2010505111.tif 000092
【図２】
2010505111.tif 000093
【図３】
2010505111.tif 000094
【図４】
2010505111.tif 000095
【図５】
2010505111.tif 000096
【図６】
2010505111.tif 000097
【図７】
2010505111.tif 000098
【図８】
2010505111.tif 000099
【図９】
2010505111.tif 000100
【国際調査報告】
2010505111.tif 000101
2010505111.tif 000102