JPH0453393B2 - - Google Patents
Info
- Publication number
- JPH0453393B2 JPH0453393B2 JP60167007A JP16700785A JPH0453393B2 JP H0453393 B2 JPH0453393 B2 JP H0453393B2 JP 60167007 A JP60167007 A JP 60167007A JP 16700785 A JP16700785 A JP 16700785A JP H0453393 B2 JPH0453393 B2 JP H0453393B2
- Authority
- JP
- Japan
- Prior art keywords
- target
- value
- motion
- time interval
- observation
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Lifetime
Links
Landscapes
- Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
- Radar Systems Or Details Thereof (AREA)
Description
〔産業上の利用分野〕
この発明は、例えば飛行目標等の追尾目標をレ
ーダ等の観測手段により追尾するため、実際の観
測値等から次の予測状態値等を得るための追尾フ
イルタに関するものである。 〔従来の技術〕 第４図は、レーダにより飛行目標を追尾する従
来の追尾フイルタの構成図であり、図において１
は図示しないレーダにより観測され、数値化され
た観測値の入力線、２は同じくレーダにより与え
られる観測誤差分散値の入力線、３は極座標系か
ら直交座標系への座標変換器、４は観測誤差分散
値と予測誤差分散値を入力して推定ゲインを出力
する推定ゲイン計算器、５は予測状態値と観測値
との残差と、予測状態値と、推定ゲインを入力し
て推定状態値を出力する状態推定器、６は推定ゲ
インと予測誤差分散値を入力して推定誤差分散値
を出力する推定誤差分散計算器、７は推定状態値
を入力して予測状態値を出力する状態予測器、８
は推定誤差分散値を入力して予測誤差分散値を出
力する予測誤差分散計算器、９は予測状態値出力
線、１０は予測誤差分散値出力線、１１は１時刻
分の時間遅延を意味する遅延要素１２は直交座標
系から極座標系への座標変換器である。 次に、飛行目標の運動と観測過程をモデル化
し、従来の追尾フイルタの動作を説明する。 〓K＋１＝〓〓K＋〓K (1) ここで、添字Ｋは時刻Ｋにおける値であることを
示し、太文字はその値がベクトルまたは行列であ
ることを示す。式中、〓Kは目標の状態値であり、
直交座標系における３方向（ｘ，ｙ，ｚ）の位置
と速度により構成される６次ベクトルである。〓
は６×６の遷移行列、〓はシステム雑音として扱
われる加速度雑音の影響であり、６次ベツトルと
なる。これらは、それぞれ次に示すように表され
る。 〓K＝Ｘ方向の位置 Ｘ方向の速度 ｙ方向の位置 ｙ方向の速度 ｚ方向の位置 ｚ方向の速度K ，〓＝1T0000 010000 001T00 000100 00001T 000001 〓K＝加速度雑音がｘ方向の位置に及す影響 加速度雑音がｘ方向の速度に及す影響 加速度雑音がｙ方向の位置に及す影響 加速度雑音がｙ方向の速度に及す影響 加速度雑音がｚ方向の位置に及す影響 加速度雑音がｚ方向の速度に及す影響K ここで、Ｔは観測時間間隔（以下サンプリング
時間と記す）であり、この従来の追尾フイルタで
は常に一定である。 上記システム雑音〓Kの統計的性質は次のとお
りである。
ーダ等の観測手段により追尾するため、実際の観
測値等から次の予測状態値等を得るための追尾フ
イルタに関するものである。 〔従来の技術〕 第４図は、レーダにより飛行目標を追尾する従
来の追尾フイルタの構成図であり、図において１
は図示しないレーダにより観測され、数値化され
た観測値の入力線、２は同じくレーダにより与え
られる観測誤差分散値の入力線、３は極座標系か
ら直交座標系への座標変換器、４は観測誤差分散
値と予測誤差分散値を入力して推定ゲインを出力
する推定ゲイン計算器、５は予測状態値と観測値
との残差と、予測状態値と、推定ゲインを入力し
て推定状態値を出力する状態推定器、６は推定ゲ
インと予測誤差分散値を入力して推定誤差分散値
を出力する推定誤差分散計算器、７は推定状態値
を入力して予測状態値を出力する状態予測器、８
は推定誤差分散値を入力して予測誤差分散値を出
力する予測誤差分散計算器、９は予測状態値出力
線、１０は予測誤差分散値出力線、１１は１時刻
分の時間遅延を意味する遅延要素１２は直交座標
系から極座標系への座標変換器である。 次に、飛行目標の運動と観測過程をモデル化
し、従来の追尾フイルタの動作を説明する。 〓K＋１＝〓〓K＋〓K (1) ここで、添字Ｋは時刻Ｋにおける値であることを
示し、太文字はその値がベクトルまたは行列であ
ることを示す。式中、〓Kは目標の状態値であり、
直交座標系における３方向（ｘ，ｙ，ｚ）の位置
と速度により構成される６次ベクトルである。〓
は６×６の遷移行列、〓はシステム雑音として扱
われる加速度雑音の影響であり、６次ベツトルと
なる。これらは、それぞれ次に示すように表され
る。 〓K＝Ｘ方向の位置 Ｘ方向の速度 ｙ方向の位置 ｙ方向の速度 ｚ方向の位置 ｚ方向の速度K ，〓＝1T0000 010000 001T00 000100 00001T 000001 〓K＝加速度雑音がｘ方向の位置に及す影響 加速度雑音がｘ方向の速度に及す影響 加速度雑音がｙ方向の位置に及す影響 加速度雑音がｙ方向の速度に及す影響 加速度雑音がｚ方向の位置に及す影響 加速度雑音がｚ方向の速度に及す影響K ここで、Ｔは観測時間間隔（以下サンプリング
時間と記す）であり、この従来の追尾フイルタで
は常に一定である。 上記システム雑音〓Kの統計的性質は次のとお
りである。
しかしながら、上記のような従来の追尾フイル
タでは、サンプリング時間Ｔ、及び、予測におけ
る目標の運動特性の仮定値が一定であつたため、
例えば飛行目標が旋回運動を始めた場合には、変
化の激しい運動特性を仮定して速いサンプリン
グ・レートで観測しなければならないのにもかか
わらず、それができないために目標を見失つた
り、反対に、飛行目標が等速直線運動をしている
場合には、変化が微小な運動特性を仮定すれば粗
いサンプリング・レートでも追尾が可能であるの
にもかかわらず、それができないたに必要以上の
観測をしてしまうなどという問題があつた。 この発明は、かかる問題点を解決するためにな
されたもので、等速直線運動や旋回運動等の追尾
目標の運動状態に応じて観測時間間隔及び運動特
性の仮定値を変えることができるようにして、観
測回数の削減が図れるとともに、追尾目標を見失
うことのない追尾フイルタを得ることを目的とす
る。 〔問題点を解決するための手段〕 この発明に係る追尾フイルタは、予測状態値と
観測値との残差と、観測誤差分散値と、予測誤差
分散値とを入力して追尾目標の運動状態を検出す
る検出手段と、観測時間間隔指示信号により観測
時間間隔を切り換える第１切換手段と、上記検出
手段の検出結果にもとづき上記第１切換手段に観
測時間間隔指示信号を出力する観測時間間隔指示
手段と、上記検出結果にもとづき運動特性の仮定
値を切り換える第２切換手段とを備えたものであ
る。 〔作用〕 この発明においては、追尾目標の運動状態を検
出し、それにしたがつてサンプリング時間と、目
標の運動特性の仮定値とを適応的に切り換えるこ
とにより、目標が等速直線運動等を行つている場
合には観測頻度を節約し、目標が旋回運動等を行
つている場合には目標の運動変化に即応してこを
見失わないように追尾を行う。 〔実施例〕 以下、この発明の一実施例を図について説明す
る。 第１図はこの発明の一実施例を示す構成図であ
り、１〜１２は上記従来装置と全く同一のもので
ある。１３は旋回検出器であり、予測状態値と実
際の観測値との残差、観測誤差分散値、及び予測
誤差分散値を入力して目標の旋回状態を検出する
検出手段である。１４は観測時間間隔指示手段と
してのサンプリング指示器である。１５は状態予
測器７の入力切換スイツチ、１６は予測誤差分散
計算器８の入力切換スイツチ、１７はサンプリン
グ時間指示出力線であり、これらは、上記旋回検
出器１３の検出結果にしたがつて上記サンプリン
グ指示器１４より出力されるサンプリング指示信
号によりサンプリング時間を適応的に変化させる
第１切換手段である。１８は目標運動特性仮定値
の切換スイツチであり、イとロが目標の運動特性
の仮定値であつて、それぞれ等速直線運動、旋回
運動に対応している。これらは、旋回検出器１３
の検出結果により目標の運動特性の仮定値を切り
換える第２切換手段である。 次に、上記実施例の動作を第２図及び第３図を
参照しながら説明する。 まず、基本サンプリング時間Ｔというものを考
え、実際のサンプリング時間はこの整数倍である
ものとする。したがつて、時間Ｔ毎に、観測は行
われるか行われないかであり、上記実施例の追尾
フイルタも基本サンプリング時間Ｔ毎にデータを
更新する。 旋回検出器１３は、レーダにより観測が行われ
たときのみ動作する。そして、このとき、飛行目
標の予測状態値と実際の観測値の残差〓Kと、観
測誤差分散値〓Kと、予測誤差分散値〓Kを入力
し、目標の旋回状態を検出する。この検出には、
x2検定を行う。観測雑音〓Kが３方向で独立な正
規白色雑音なら、残差〓Kも３方向で独立な正規
白色雑音となり、その３方向の正規化された二乗
和ZKは、自由度３のx2分布にしたがう。 〓K＝Ｅ〓〓K〓T〓＝〓K＋〓〓K〓T (9) ZK＝〓T K〓-1〓K (10) ここで、次の仮定H0，H1を設ける。 H0：目標は等速直線運動を行つている。 H1：目標は等速直線運動を行つていない。 この仮定H1に対し、H0の検定を行うための棄
却域Ｃを次のように定める。一般に、自由度φの
x2分布の上側100α％の点はx2（φ，α）と書か
れ、数表として与えられている。これを利用すれ
ば、危険度ｄで検出する際の棄却域Ｃは Z(o)＝x2（３，ｄ） （11） Ｃ ＝〓ZK｜ZK＞Z(o)〓 （12） となる。すなわち、第10式で得られたZKの値が、
指定した危険度ｄで決まる値Z(o)よりも大きな場
合、H1「目標は旋回運動を行つている」と判断
し、出力をhK＝１とする。そうでない場合には、
H0「目標は旋回運動を行つていない」と判断して
出力をhK＝０とする。また、観測が行われず、こ
の旋回検出器１３が動作しない場合には、この旋
回検出信号をhK＝０となつているようにする。 サンプリング指示器１４は、この旋回検出信号
hKを入力し、これにより次の２つのモードを選択
してサンプリング指示信号SKを出力する。 oAモード（等速直線運動対応モード）……サ
ンプリング時間NTで目標を観測し、微小な運動
変化を仮定してフイルタリングを行う。 oBモード（旋回運動対応モード）……サンプ
リングリング時間Ｔで目標を観測し、激しい運動
変化を仮定してフイルタリングを行う。 ここで、Ｎは２以上の整数である。この選択
は、過去N′T間で、旋回が検出されたかどうかに
よる。つまり、過去N′T間で、旋回が検出されて
いれば上記のＢモードを選び、検出されていなけ
ればＡモードを選ぶ。ここで、N′は１以上の整
数であり、 １≦N′≦Ｎ 〓 である。そして、選んだＡモード、Ｂモードにし
たがつて、サンプリング指示信号SKを、時間Ｔ
毎に出力する。この信号の意味は、SK＝０「サン
プリングなし」、SK＝１「サンプリングあり」と
いう指示である。 上記、旋回検出信号hKとサンプリング指示信号
SKの関係を第２図に示す。 また、上記のＡモード、Ｂモードにより、目標
の運動特性の仮定値として、フイルタリングにお
ける加速度雑音スペクトラル密度ｇを変数gKとし
て切り換える。これにより〓も変数〓Kとなる。 前記サンプリング指示信号SKは、状態予測器
７の入力切換スイツチ１５、予測誤差分散計算器
８の入力切換スイツチ１６、サンプリング時間指
示出力線１７に伝えられて、それぞれ以下のよう
な働きをする。 状態予測器７の入力切換スイツチ１５は、サン
プリング指示信号SKにより、次回（時刻ｋ＋１）
の状態予測器７への入力を切り換える。SK＝１
のとき、スイツチ１５は状態推定器５側に接続
し、入力として推定状態値〓K＋１を選ぶ。 SK＝０のときには、遅延要素１１側に接続し、
予測状態値〓K＋１を選ぶ。 予測誤差分散計算器８の入力切換スイツチ１６
も、サンプリング指示信号SKにより、次回（時
刻ｋ＋１）の予測誤差分散計算器８への入力を切
り換える。SK＝１のときスイツチ１６は推定誤
差分散計算器６側に接続し、入力として推定誤差
分散値〓＋１を選ぶ。SK＝０のときには、遅延
要素１１側に接続し、予測誤差分散値〓K＋１を
選ぶ。 さらに、このサンプリング指示信号SKは、サ
ンプリング時間指示出力線１７からレーダに送ら
れ、次回（時刻ｋ＋１）にその目標を観測するか
どうかを指示する。 一方、目標の運動特性仮定値の切換スイツチ１
８は、旋回検出信号hKにより、次回（時刻ｋ＋
１）の予測誤差分散計算器８における加速度雑音
のスペクトラル密度〓を、外部から切り換える。
この値は、目標の運動特性の仮定値であり、この
値を大きくとれば、「目標は激しい運動変化を起
す」という仮定の下に予測が行われ、この値を小
さくとれば、「目標は微小な運動変化しか起さな
い」という仮定の下で予測が行われる。そこで、
この値を２種類用意しておいて、それらを、旋回
の検出・不検出にしたがつて、サンプリング時間
とともに使い分ける。その２種類の値とは、次の
２つであり、第１図中イ、ロのように、外部にあ
らかじめ設定されている。 〓(o)：目標は微小な運動変化しか起さないとい
う仮定の下の加速度雑音スペクトラル密度 〓(1)：目標は激しい運動変化を起すという仮定
の下の加速度雑音スペクトラル密度 つまり、目標の運動特性仮定値の切替スイツチ
１８は、前記Ｂモードのとき、〓K＋１＝〓(1)と
し、前記Ａモードのとき〓K＋１＝〓(o)とする。 ここで、従来のフイルタと異なり、加速度雑音
のスペクトラル密度を〓Kとして、変数として扱
つていることは重要である。これにより、〓も変
数〓Kとなり、予測誤差分散計算器８は、外部か
らの指示で〓K、すなわち〓Kを切り換えて、第８
式の代わりに次の計算を行うことになる。 〓K＋１＝〓〓K〓T＋〓K （8A） ここで〓Kは
タでは、サンプリング時間Ｔ、及び、予測におけ
る目標の運動特性の仮定値が一定であつたため、
例えば飛行目標が旋回運動を始めた場合には、変
化の激しい運動特性を仮定して速いサンプリン
グ・レートで観測しなければならないのにもかか
わらず、それができないために目標を見失つた
り、反対に、飛行目標が等速直線運動をしている
場合には、変化が微小な運動特性を仮定すれば粗
いサンプリング・レートでも追尾が可能であるの
にもかかわらず、それができないたに必要以上の
観測をしてしまうなどという問題があつた。 この発明は、かかる問題点を解決するためにな
されたもので、等速直線運動や旋回運動等の追尾
目標の運動状態に応じて観測時間間隔及び運動特
性の仮定値を変えることができるようにして、観
測回数の削減が図れるとともに、追尾目標を見失
うことのない追尾フイルタを得ることを目的とす
る。 〔問題点を解決するための手段〕 この発明に係る追尾フイルタは、予測状態値と
観測値との残差と、観測誤差分散値と、予測誤差
分散値とを入力して追尾目標の運動状態を検出す
る検出手段と、観測時間間隔指示信号により観測
時間間隔を切り換える第１切換手段と、上記検出
手段の検出結果にもとづき上記第１切換手段に観
測時間間隔指示信号を出力する観測時間間隔指示
手段と、上記検出結果にもとづき運動特性の仮定
値を切り換える第２切換手段とを備えたものであ
る。 〔作用〕 この発明においては、追尾目標の運動状態を検
出し、それにしたがつてサンプリング時間と、目
標の運動特性の仮定値とを適応的に切り換えるこ
とにより、目標が等速直線運動等を行つている場
合には観測頻度を節約し、目標が旋回運動等を行
つている場合には目標の運動変化に即応してこを
見失わないように追尾を行う。 〔実施例〕 以下、この発明の一実施例を図について説明す
る。 第１図はこの発明の一実施例を示す構成図であ
り、１〜１２は上記従来装置と全く同一のもので
ある。１３は旋回検出器であり、予測状態値と実
際の観測値との残差、観測誤差分散値、及び予測
誤差分散値を入力して目標の旋回状態を検出する
検出手段である。１４は観測時間間隔指示手段と
してのサンプリング指示器である。１５は状態予
測器７の入力切換スイツチ、１６は予測誤差分散
計算器８の入力切換スイツチ、１７はサンプリン
グ時間指示出力線であり、これらは、上記旋回検
出器１３の検出結果にしたがつて上記サンプリン
グ指示器１４より出力されるサンプリング指示信
号によりサンプリング時間を適応的に変化させる
第１切換手段である。１８は目標運動特性仮定値
の切換スイツチであり、イとロが目標の運動特性
の仮定値であつて、それぞれ等速直線運動、旋回
運動に対応している。これらは、旋回検出器１３
の検出結果により目標の運動特性の仮定値を切り
換える第２切換手段である。 次に、上記実施例の動作を第２図及び第３図を
参照しながら説明する。 まず、基本サンプリング時間Ｔというものを考
え、実際のサンプリング時間はこの整数倍である
ものとする。したがつて、時間Ｔ毎に、観測は行
われるか行われないかであり、上記実施例の追尾
フイルタも基本サンプリング時間Ｔ毎にデータを
更新する。 旋回検出器１３は、レーダにより観測が行われ
たときのみ動作する。そして、このとき、飛行目
標の予測状態値と実際の観測値の残差〓Kと、観
測誤差分散値〓Kと、予測誤差分散値〓Kを入力
し、目標の旋回状態を検出する。この検出には、
x2検定を行う。観測雑音〓Kが３方向で独立な正
規白色雑音なら、残差〓Kも３方向で独立な正規
白色雑音となり、その３方向の正規化された二乗
和ZKは、自由度３のx2分布にしたがう。 〓K＝Ｅ〓〓K〓T〓＝〓K＋〓〓K〓T (9) ZK＝〓T K〓-1〓K (10) ここで、次の仮定H0，H1を設ける。 H0：目標は等速直線運動を行つている。 H1：目標は等速直線運動を行つていない。 この仮定H1に対し、H0の検定を行うための棄
却域Ｃを次のように定める。一般に、自由度φの
x2分布の上側100α％の点はx2（φ，α）と書か
れ、数表として与えられている。これを利用すれ
ば、危険度ｄで検出する際の棄却域Ｃは Z(o)＝x2（３，ｄ） （11） Ｃ ＝〓ZK｜ZK＞Z(o)〓 （12） となる。すなわち、第10式で得られたZKの値が、
指定した危険度ｄで決まる値Z(o)よりも大きな場
合、H1「目標は旋回運動を行つている」と判断
し、出力をhK＝１とする。そうでない場合には、
H0「目標は旋回運動を行つていない」と判断して
出力をhK＝０とする。また、観測が行われず、こ
の旋回検出器１３が動作しない場合には、この旋
回検出信号をhK＝０となつているようにする。 サンプリング指示器１４は、この旋回検出信号
hKを入力し、これにより次の２つのモードを選択
してサンプリング指示信号SKを出力する。 oAモード（等速直線運動対応モード）……サ
ンプリング時間NTで目標を観測し、微小な運動
変化を仮定してフイルタリングを行う。 oBモード（旋回運動対応モード）……サンプ
リングリング時間Ｔで目標を観測し、激しい運動
変化を仮定してフイルタリングを行う。 ここで、Ｎは２以上の整数である。この選択
は、過去N′T間で、旋回が検出されたかどうかに
よる。つまり、過去N′T間で、旋回が検出されて
いれば上記のＢモードを選び、検出されていなけ
ればＡモードを選ぶ。ここで、N′は１以上の整
数であり、 １≦N′≦Ｎ 〓 である。そして、選んだＡモード、Ｂモードにし
たがつて、サンプリング指示信号SKを、時間Ｔ
毎に出力する。この信号の意味は、SK＝０「サン
プリングなし」、SK＝１「サンプリングあり」と
いう指示である。 上記、旋回検出信号hKとサンプリング指示信号
SKの関係を第２図に示す。 また、上記のＡモード、Ｂモードにより、目標
の運動特性の仮定値として、フイルタリングにお
ける加速度雑音スペクトラル密度ｇを変数gKとし
て切り換える。これにより〓も変数〓Kとなる。 前記サンプリング指示信号SKは、状態予測器
７の入力切換スイツチ１５、予測誤差分散計算器
８の入力切換スイツチ１６、サンプリング時間指
示出力線１７に伝えられて、それぞれ以下のよう
な働きをする。 状態予測器７の入力切換スイツチ１５は、サン
プリング指示信号SKにより、次回（時刻ｋ＋１）
の状態予測器７への入力を切り換える。SK＝１
のとき、スイツチ１５は状態推定器５側に接続
し、入力として推定状態値〓K＋１を選ぶ。 SK＝０のときには、遅延要素１１側に接続し、
予測状態値〓K＋１を選ぶ。 予測誤差分散計算器８の入力切換スイツチ１６
も、サンプリング指示信号SKにより、次回（時
刻ｋ＋１）の予測誤差分散計算器８への入力を切
り換える。SK＝１のときスイツチ１６は推定誤
差分散計算器６側に接続し、入力として推定誤差
分散値〓＋１を選ぶ。SK＝０のときには、遅延
要素１１側に接続し、予測誤差分散値〓K＋１を
選ぶ。 さらに、このサンプリング指示信号SKは、サ
ンプリング時間指示出力線１７からレーダに送ら
れ、次回（時刻ｋ＋１）にその目標を観測するか
どうかを指示する。 一方、目標の運動特性仮定値の切換スイツチ１
８は、旋回検出信号hKにより、次回（時刻ｋ＋
１）の予測誤差分散計算器８における加速度雑音
のスペクトラル密度〓を、外部から切り換える。
この値は、目標の運動特性の仮定値であり、この
値を大きくとれば、「目標は激しい運動変化を起
す」という仮定の下に予測が行われ、この値を小
さくとれば、「目標は微小な運動変化しか起さな
い」という仮定の下で予測が行われる。そこで、
この値を２種類用意しておいて、それらを、旋回
の検出・不検出にしたがつて、サンプリング時間
とともに使い分ける。その２種類の値とは、次の
２つであり、第１図中イ、ロのように、外部にあ
らかじめ設定されている。 〓(o)：目標は微小な運動変化しか起さないとい
う仮定の下の加速度雑音スペクトラル密度 〓(1)：目標は激しい運動変化を起すという仮定
の下の加速度雑音スペクトラル密度 つまり、目標の運動特性仮定値の切替スイツチ
１８は、前記Ｂモードのとき、〓K＋１＝〓(1)と
し、前記Ａモードのとき〓K＋１＝〓(o)とする。 ここで、従来のフイルタと異なり、加速度雑音
のスペクトラル密度を〓Kとして、変数として扱
つていることは重要である。これにより、〓も変
数〓Kとなり、予測誤差分散計算器８は、外部か
らの指示で〓K、すなわち〓Kを切り換えて、第８
式の代わりに次の計算を行うことになる。 〓K＋１＝〓〓K〓T＋〓K （8A） ここで〓Kは
以上のように、この発明による追尾フイルタ
は、予測状態値と観測値との残差と、観測誤差分
散値と、予測誤差分散値とを入力して追尾目標の
運動状態を検出する検出手段と、この検出結果に
もとづき観測時間間隔を切換える第１切替手段
と、上記検出結果にもとづき運動特性の仮定値を
切り換える第２切換手段とを備えたことにより、
追尾目標の運動状態に応じて観測時間間隔と運動
特性の仮定値とを変えることができるので、観測
回数の削減が図れるとともに、追尾目標を見失う
ことのない追尾フイルタが得られるという効果が
ある。
は、予測状態値と観測値との残差と、観測誤差分
散値と、予測誤差分散値とを入力して追尾目標の
運動状態を検出する検出手段と、この検出結果に
もとづき観測時間間隔を切換える第１切替手段
と、上記検出結果にもとづき運動特性の仮定値を
切り換える第２切換手段とを備えたことにより、
追尾目標の運動状態に応じて観測時間間隔と運動
特性の仮定値とを変えることができるので、観測
回数の削減が図れるとともに、追尾目標を見失う
ことのない追尾フイルタが得られるという効果が
ある。
第１図はこの発明による追尾フイルタの一実施
例を示す構成図、第２図は上記実施例の旋回検出
信号とサンプリング指示信号のタイミング図、第
３図はその動作説明のためのフロー・チヤート、
第４図は従来の追尾フイルタを示す構成図であ
る。 図において、１は観測値入力線、２は観測誤差
分散値入力線、３は極座標系から直交座標系への
座標交換器、４は推定ゲイン計算器、５は状態推
定器、６は推定誤差分散計算器、７は状態予測
器、８は予測誤差分散計算器、９は予測状態出力
線、１０は予測誤差分散出力線、１１は遅延要
素、１２は直交座標系から極座標系への座標変換
器、１３は旋回検出器（検出手段）、１４はサン
プリング指示器、１５は状態予測器７の入力切換
スイツチ、１６は予測誤差分散計算器８の入力切
換スイツチ、１７はサンプリング時間指示出力
線、上記１４〜１７は第１切換手段、１８は目標
の運動特性仮定値の切換スイツチ（第２切換手
段）、イ、ロは運動特性仮定値である。なお、図
中同一符号は同一または相当部分を示す。
例を示す構成図、第２図は上記実施例の旋回検出
信号とサンプリング指示信号のタイミング図、第
３図はその動作説明のためのフロー・チヤート、
第４図は従来の追尾フイルタを示す構成図であ
る。 図において、１は観測値入力線、２は観測誤差
分散値入力線、３は極座標系から直交座標系への
座標交換器、４は推定ゲイン計算器、５は状態推
定器、６は推定誤差分散計算器、７は状態予測
器、８は予測誤差分散計算器、９は予測状態出力
線、１０は予測誤差分散出力線、１１は遅延要
素、１２は直交座標系から極座標系への座標変換
器、１３は旋回検出器（検出手段）、１４はサン
プリング指示器、１５は状態予測器７の入力切換
スイツチ、１６は予測誤差分散計算器８の入力切
換スイツチ、１７はサンプリング時間指示出力
線、上記１４〜１７は第１切換手段、１８は目標
の運動特性仮定値の切換スイツチ（第２切換手
段）、イ、ロは運動特性仮定値である。なお、図
中同一符号は同一または相当部分を示す。
Claims (1)
- 【特許請求の範囲】 １ 追尾目標を指定時間間隔で観測する観測手段
より上記追尾目標の観測値と観測誤差分散値とを
入力し、これらの値と上記追尾目標の予測におけ
る運動特性の仮定値とにもとづき害追尾目標の予
測状態値と予測誤差分散値とを出力するカルマン
型追尾フイルタにおいて、上記予測状態値と観測
値との残差と上記観測誤差分散値と上記予測誤差
分散値とを入力して上記追尾目標の運動状態を検
出する検出手段と、上記検出手段の検出結果にも
とづき観測時間間隔指示信号を出力する観測時間
間隔指示手段と、上記観測時間間隔指示手段から
の観測時間間隔指示信号により観測時間間隔を切
り換える第１切換手段と、上記検出結果にもとづ
き上記運動特性の仮定値を切り換える第２切換手
段とを備えたことを特徴とする追尾フイルタ。 ２ 追尾目標を飛行目標とし、観測手段をレーダ
とし、検出手段は上記飛行目標が等速直線運動か
旋回運動かの運動状態を検出し、第１切換手段
は、上記運動状態にもとづく観測時間間隔指示手
段からの指示信号により、飛行目標が等速直線運
動のときは長い観測時間間隔に、旋回運動のとき
は短い観測時間間隔に切り換え、第２切換手段
は、上記運動状態を示す信号にもとづき、飛行目
標が等速直線運動のときは運動特性の小なる仮定
値に、旋回運動のときは運動特性の大なる仮定値
に切替えることを特徴とする特許請求の範囲第１
項記載の追尾フイルタ。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP60167007A JPS6227677A (ja) | 1985-07-29 | 1985-07-29 | 追尾フイルタ |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP60167007A JPS6227677A (ja) | 1985-07-29 | 1985-07-29 | 追尾フイルタ |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPS6227677A JPS6227677A (ja) | 1987-02-05 |
| JPH0453393B2 true JPH0453393B2 (ja) | 1992-08-26 |
Family
ID=15841648
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP60167007A Granted JPS6227677A (ja) | 1985-07-29 | 1985-07-29 | 追尾フイルタ |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JPS6227677A (ja) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2500347B2 (ja) * | 1992-08-21 | 1996-05-29 | 防衛庁技術研究本部長 | ３個の１軸磁力計を用いた移動物体の位置局限装置 |
| JP4977492B2 (ja) * | 2007-02-26 | 2012-07-18 | 日本電気株式会社 | 目標追尾装置及び目標追尾方法 |
Family Cites Families (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS58135978A (ja) * | 1982-02-09 | 1983-08-12 | Mitsubishi Electric Corp | 追尾フイルタ |
| JPS58191520A (ja) * | 1982-05-04 | 1983-11-08 | Mitsubishi Electric Corp | デイジタルフイルタ |
| JPS59222777A (ja) * | 1983-06-02 | 1984-12-14 | Mitsubishi Electric Corp | 追尾装置 |
-
1985
- 1985-07-29 JP JP60167007A patent/JPS6227677A/ja active Granted
Also Published As
| Publication number | Publication date |
|---|---|
| JPS6227677A (ja) | 1987-02-05 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Jilkov et al. | Design and comparison of mode-set adaptive IMM algorithms for maneuvering target tracking | |
| Blair | Fixed-gain two-stage estimators for tracking maneuvering targets | |
| JP4928171B2 (ja) | レーダ装置及びレーダ画像表示方法 | |
| JPH0453393B2 (ja) | ||
| SE503641C2 (sv) | Fasstyrt radarsystem för spårföljning | |
| JP3850950B2 (ja) | ドップラソナー | |
| JPH0453394B2 (ja) | ||
| JP3674952B2 (ja) | 射撃統制装置 | |
| Bar-Shalom et al. | Precision target tracking for small extended objects | |
| GB1535338A (en) | Radar systems | |
| JP2003149328A (ja) | レーダの目標相関装置 | |
| Nomoto et al. | An adaptive tracking filter controlling sampling frequency | |
| JPH03245081A (ja) | 主センサ決定システム | |
| JPS59220666A (ja) | 目標位置予測装置 | |
| JP3432737B2 (ja) | 追尾装置 | |
| Efe et al. | A tracking algorithm for both highly maneuvering and nonmaneuvering targets | |
| JPS6227678A (ja) | 追尾フイルタ | |
| JPS61191977A (ja) | 追尾フイルタ | |
| JPS61184476A (ja) | 追尾フイルタ | |
| JPH06217974A (ja) | 超音波診断装置 | |
| JPS624668B2 (ja) | ||
| JPH028668B2 (ja) | ||
| JPS63196816A (ja) | 船舶用運動予測方式 | |
| RU2103713C1 (ru) | Способ управления и система для его осуществления | |
| Karan et al. | Bearings-only tracking with sea trial sonar data from multiple asynchronous sonobuoys |