JP2016161269A - Target state quantity estimation device - Google Patents
Target state quantity estimation device Download PDFInfo
- Publication number
- JP2016161269A JP2016161269A JP2015043108A JP2015043108A JP2016161269A JP 2016161269 A JP2016161269 A JP 2016161269A JP 2015043108 A JP2015043108 A JP 2015043108A JP 2015043108 A JP2015043108 A JP 2015043108A JP 2016161269 A JP2016161269 A JP 2016161269A
- Authority
- JP
- Japan
- Prior art keywords
- observation
- noise matrix
- target
- state quantity
- matrix
- 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.)
- Granted
Links
- 239000011159 matrix material Substances 0.000 claims abstract description 95
- 238000004364 calculation method Methods 0.000 claims abstract description 69
- 238000000034 method Methods 0.000 claims abstract description 47
- 238000013213 extrapolation Methods 0.000 claims abstract description 28
- 230000001133 acceleration Effects 0.000 claims description 4
- 230000002123 temporal effect Effects 0.000 abstract 5
- 230000001052 transient effect Effects 0.000 description 6
- 230000006698 induction Effects 0.000 description 3
- 238000013178 mathematical model Methods 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 230000007423 decrease Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
Images
Landscapes
- Radar Systems Or Details Thereof (AREA)
- Aiming, Guidance, Guns With A Light Source, Armor, Camouflage, And Targets (AREA)
Abstract
Description
この発明は、飛しょう体を目標に会合させるにあたり、目標の位置・速度・加速度等の状態量を推定するための目標状態量推定装置に関するものである。 The present invention relates to a target state quantity estimating device for estimating a state quantity such as a target position, velocity, and acceleration when a flying object is associated with a target.
従来、飛しょう体等に搭載されているシーカによって目標の状態量を推定する場合、カルマンフィルタが広く用いられてきた。シーカから得られる目標の観測結果には観測誤差が含まれており、それを取り除く手段として、目標の運動モデルに基づいて状態量を推定するカルマンフィルタは非常に有効である。カルマンフィルタには、線形カルマンフィルタ(例えば、非特許文献1のP.163〜P.184参照)や、非線形カルマンフィルタ(例えば、非特許文献1のP.373〜P.388参照)等の種類がある。また、従来、目標の状態量を推定する装置が開示されている(例えば、特許文献1参照)。 Conventionally, a Kalman filter has been widely used to estimate a target state quantity with a seeker mounted on a flying object or the like. The observation result of the target obtained from the seeker includes an observation error, and a Kalman filter that estimates the state quantity based on the target motion model is very effective as a means for removing it. The Kalman filter includes types such as a linear Kalman filter (for example, see P.163 to P.184 of Non-Patent Document 1) and a nonlinear Kalman filter (for example, see P.373 to P.388 of Non-Patent Document 1). Conventionally, an apparatus for estimating a target state quantity has been disclosed (see, for example, Patent Document 1).
カルマンフィルタの特性は、目標の運動モデルに加えて、シーカによる観測誤差の大きさを表す観測ノイズ行列と、モデル誤差の大きさを表すプロセスノイズ行列の値によって決定される。このうち観測ノイズ行列については、シーカのS/N比や観測誤差モデルなどによって一意に決定される。それに対してプロセスノイズ行列は、目標の運動が予測できないことから一意に決定することができず、一般的に設計パラメータとなる。 The characteristics of the Kalman filter are determined by the observation noise matrix that represents the magnitude of the observation error due to the seeker and the value of the process noise matrix that represents the magnitude of the model error, in addition to the target motion model. Of these, the observation noise matrix is uniquely determined by the Seeker S / N ratio and the observation error model. In contrast, the process noise matrix cannot be uniquely determined because the target motion cannot be predicted, and is generally a design parameter.
カルマンフィルタに求められる応答特性は、シーカの観測値に基づいて誘導を開始してから会合に至るまでの間に変化する。誘導の開始直後においては、フィルタの推定誤差による誘導弾の余分な機動が会合の妨げとなるため、定常特性が重視される。一方で会合直前においては、フィルタの時間遅れが誘導時間の不足に直結するため、過渡特性が重視される。
定常特性と過渡特性はトレードオフの関係にある。プロセスノイズ行列を設定するにあたっては、その点を考慮する必要がある。
The response characteristic required for the Kalman filter changes between the start of guidance and the meeting based on the observed value of the seeker. Immediately after the start of the guidance, the excessive movement of the guided bullet due to the estimation error of the filter hinders the association, so that steady characteristics are emphasized. On the other hand, immediately before the meeting, the time delay of the filter is directly linked to the shortage of the induction time, so the transient characteristics are emphasized.
Steady state characteristics and transient characteristics are in a trade-off relationship. This must be taken into account when setting the process noise matrix.
レーダー反射断面積(RCS)が大きく、観測値のS/N比が高い目標の場合は、観測誤差が小さいため、プロセスノイズ行列の値がある一種類に固定されていても、定常特性と過渡特性に関する性能指標を両方満たすことができる。
しかし、ステルス機のようにRCSが小さく、観測値のS/N比が低い目標の場合は、観測誤差が大きいために、定常特性と過渡特性は大きく劣化する。この場合、プロセスノイズ行列の値がある一種類に固定されている状態では、定常特性と過渡特性に関する性能指標を同時に満たすことが難しいという課題があった。
For targets with a large radar reflection cross-section (RCS) and high S / N ratio, the observation error is small, so even if the process noise matrix value is fixed to a certain type, steady-state characteristics and transients Both performance indicators regarding characteristics can be met.
However, in the case of a target having a small RCS and a low S / N ratio of the observed value, such as a stealth machine, since the observation error is large, the steady-state characteristic and the transient characteristic are greatly deteriorated. In this case, in the state where the value of the process noise matrix is fixed to a certain type, there is a problem that it is difficult to satisfy the performance index related to the steady characteristic and the transient characteristic at the same time.
この発明は係る課題を解決するためになされたものであり、シーカの観測値のS/N比が低い状況にあっても、目標の運動を高精度に推定する目標運動量推定装置を提供することを目標とする。 The present invention has been made to solve such a problem, and provides a target momentum estimation device that estimates a target motion with high accuracy even in a situation where the S / N ratio of the observed value of a seeker is low. To the goal.
この発明に係る目標状態量推定装置は、目標の状態量をカルマンフィルタにより推定する目標状態量推定装置であって、飛しょう体に搭載され、目標までの距離と方向を観測した観測量と、前記観測量のS/N比を出力するシーカと、前記観測量と前記S/N比に基づき、前記観測量の誤差の大きさを表す観測ノイズ行列を算出し、前記観測ノイズ行列を観測更新計算装置に出力する観測ノイズ行列計算装置と、前記観測量と、前記観測ノイズ行列と、時間外挿計算装置が出力する目標の状態量(−)と誤差共分散行列(−)に基づき観測更新計算を行い、観測更新計算後の目標の状態量(+)と誤差共分散行列(+)を出力する観測更新計算装置と、更新計算後の前記目標の状態量(+)と誤差共分散行列(+)と、前記目標の運動モデルの誤差の大きさであるプロセスノイズ行列に基づき、時間外挿計算により目標の状態量(−)と誤差共分散行列(−)を算出し、出力する時間外挿計算装置と、前記時間外挿計算装置が算出した前記目標の状態量(−)に基づき、前記プロセスノイズ行列を算出し、算出した前記プロセスノイズ行列を前記時間外挿計算装置に出力するプロセスノイズ行列計算装置とから構成される。 A target state quantity estimation device according to the present invention is a target state quantity estimation apparatus that estimates a target state quantity by a Kalman filter, and is mounted on a flying object, and an observed quantity obtained by observing the distance and direction to the target, Based on the seeker that outputs the S / N ratio of the observation amount, the observation noise matrix that represents the magnitude of the observation amount error based on the observation amount and the S / N ratio, and the observation update calculation of the observation noise matrix Observation update matrix calculation device output to device, observation update calculation based on observation amount, observation noise matrix, target state quantity (-) and error covariance matrix (-) output by time extrapolation calculation device And an observation update calculation device that outputs the target state quantity (+) and error covariance matrix (+) after the observation update calculation, and the target state quantity (+) and error covariance matrix (after the update calculation) +) And the target exercise model A time extrapolation calculation device that calculates and outputs a target state quantity (−) and an error covariance matrix (−) by time extrapolation based on a process noise matrix that is the magnitude of the error, and the time extrapolation calculation The process noise matrix calculation device is configured to calculate the process noise matrix based on the target state quantity (−) calculated by the device and to output the calculated process noise matrix to the time extrapolation calculation device.
この発明に係る目標状態量推定装置によれば、シーカの観測値のS/N比が低い状況にあっても目標の運動を高精度に推定することができ、目標を高精度に追尾することができる。 According to the target state quantity estimation device according to the present invention, the target motion can be estimated with high accuracy even in a situation where the S / N ratio of the observed value of the seeker is low, and the target can be tracked with high accuracy. Can do.
実施の形態1.
図1は、本実施の形態に係る目標状態量推定装置100の構成を示す図である。
目標状態量推定装置100は、飛しょう体に搭載され、目標の距離・方向等を観測するためのシーカ1と、観測更新計算装置2と、時間外挿計算装置3と、観測ノイズ行列計算装置4と、プロセスノイズ行列計算装置5から構成される。
以下、図1を参照しながら、実施の形態1における目標状態量推定装置の動作を詳しく説明する。
Embodiment 1 FIG.
FIG. 1 is a diagram showing a configuration of target state
A target state
Hereinafter, the operation of the target state quantity estimation apparatus in the first embodiment will be described in detail with reference to FIG.
シーカ1は、目標が反射する電波を受信、あるいは目標が放射する赤外線を受信すること等により、目標の運動を観測する。
シーカ1は、観測された目標の観測量6を観測更新計算装置2へ出力する。
The seeker 1 observes the movement of the target by receiving radio waves reflected by the target or receiving infrared rays radiated by the target.
The seeker 1 outputs the observed target observation amount 6 to the observation update calculation device 2.
観測更新計算装置2は、シーカ1から得られる観測量6と、時間外挿計算装置3から得られる「状態量(−)・誤差共分散行列(−)」8を用い、観測ノイズ行列計算装置4から得られる観測ノイズ行列9の値に基づいて観測更新計算を行う。
The observation update calculation device 2 uses the observation amount 6 obtained from the seeker 1 and the “state quantity (−) / error covariance matrix (−)” 8 obtained from the time
観測更新計算の処理を下記の式(1)に示す。
シーカ1はその装置の特性上、観測量6を離散的に出力するため、観測更新計算装置2は、シーカ1から観測量6が入力されるたびに計算されるものとする。
The observation update calculation process is shown in the following equation (1).
Since the seeker 1 outputs the observation amount 6 discretely due to the characteristics of the apparatus, the observation update calculation device 2 is calculated every time the observation amount 6 is input from the seeker 1.
観測方程式:h
観測量:z
観測ノイズ行列:R
状態量(観測更新計算前):x−
状態量(観測更新計算後):x+
誤差共分散行列(観測更新計算前):P−
誤差共分散行列(観測更新計算後):P+
とする。
Observed quantity: z
Observation noise matrix: R
State quantity (before observation update calculation): x −
State quantity (after observation update calculation): x +
Error covariance matrix (before observation update calculation): P −
Error covariance matrix (after observation update calculation): P +
And
観測方程式hは、シーカの特性をあらかじめ適切に同定した数式モデルから構成される。観測ノイズ行列は、シーカから得られる観測値のS/N比7や観測誤差モデルなどをもとに、適切に設定されるものとする。誤差共分散行列の初期値は、あらかじめ想定される誤差共分散に応じて適切に設定されるものとする。
The observation equation h is composed of a mathematical model in which the characteristics of the seeker are appropriately identified in advance. The observation noise matrix is appropriately set based on the S /
以上の計算で得られた状態量(観測更新計算後)と誤差共分散行列(観測更新計算後)を「状態量(+)・誤差共分散行列(+)」10として時間外挿計算装置3へ出力する。
The time
時間外挿計算装置3は、観測更新計算装置2からの入力「状態量(+)・誤差共分散行列(+)」10と、プロセスノイズ行列計算装置5からの入力「プロセスノイズ行列」11を用いて、時間外挿計算を行う。
The time
下記の式(2)に時間外挿計算の処理を示す。 The time extrapolation calculation processing is shown in the following equation (2).
運動方程式:f
プロセスノイズ行列:Q
現在時間:tk
とする。
Process noise matrix: Q
Current time: tk
And
運動方程式fは、対象をあらかじめ適切に同定した数式モデルから構成される。 The equation of motion f is composed of a mathematical model in which an object is appropriately identified in advance.
時間外挿計算の結果として得られる状態量(時間外挿計算後)・誤差共分散行列(時間外挿計算後)を「状態量(−)・誤差共分散行列(−)」8として、観測更新計算装置2とプロセスノイズ行列計算装置5へ出力する。 State quantity (after time extrapolation calculation) and error covariance matrix (after time extrapolation calculation) obtained as a result of time extrapolation calculation are observed as "state quantity (-) and error covariance matrix (-)" 8 Output to the update calculation device 2 and the process noise matrix calculation device 5.
本実施の形態に係るプロセスノイズ行列計算装置5は、時間外挿計算装置3から得られた「状態量(−)・誤差共分散行列(−)」8を用いて、プロセスノイズ行列11の値を計算する。
プロセスノイズ行列11の値は、下記の式(3)に示す処理で計算される。
The process noise matrix calculation device 5 according to the present embodiment uses the “state quantity (−) / error covariance matrix (−)” 8 obtained from the time
The value of the
プロセスノイズ行列:Q
プロセスノイズ行列初期値:Q0
ゲイン行列関数:G
とする。
Process noise matrix initial value: Q0
Gain matrix function: G
And
Gは、状態量に基づいてプロセスノイズ行列の各行のゲインを算出する関数である。
各時刻において、状態量(−)を基にゲイン行列を算出し、プロセスノイズ行列初期値に乗じる。このように、状態量を反映したプロセスノイズ行列Qを算出することで、その場に適したフィルタ特性を得ることが可能になる。
G is a function for calculating the gain of each row of the process noise matrix based on the state quantity.
At each time, a gain matrix is calculated based on the state quantity (−) and multiplied by the process noise matrix initial value. Thus, by calculating the process noise matrix Q reflecting the state quantity, it is possible to obtain filter characteristics suitable for the situation.
一例として、状態量から計算される「予想会合時間」を用いて、プロセスノイズ行列の算出を行う場合を説明する。
予想会合時間は飛しょう体が目標に会合するまでにかかると予想される時間を表し、誘導開始時点で最も大きく、誘導中は単調減少し、会合時に0となる。
そのため、ゲイン行列Gを予想会合時間に基づいて算出することで、(1)誘導開始直後では定常特性を重視した特性、(2)会合直前では過渡特性を重視した特性というように、その場で求められる性能指標を満足したフィルタ特性を得ることができる。
As an example, a case where the process noise matrix is calculated using “expected meeting time” calculated from the state quantity will be described.
The expected meeting time represents the time that the flying object is expected to take to meet the target, and is the largest at the start of induction, decreases monotonically during induction, and becomes zero at the meeting.
Therefore, by calculating the gain matrix G based on the expected meeting time, (1) characteristics that emphasize steady characteristics immediately after the start of guidance, and (2) characteristics that emphasize transient characteristics immediately before the meeting, on the spot. A filter characteristic satisfying the required performance index can be obtained.
なお、プロセスノイズ行列初期値Q0は、目標の最大旋回加速度などをもとにあらかじめ設定される初期値である。
この計算は時間外挿計算装置3から新たな値が送信される毎に行われ、式(3)により計算されたプロセスノイズ行列Qは、プロセスノイズ行列11として時間外挿計算装置3に送信される。
The process noise matrix initial value Q0 is an initial value set in advance based on a target maximum turning acceleration or the like.
This calculation is performed every time a new value is transmitted from the time
以上の数式1〜数式3の処理を繰り返し計算し、「状態量(−)・誤差共分散行列(−)」8を毎サンプリング計算する。サンプリング毎に、「状態量(−)・誤差共分散行列(−)」8を目標状態量推定装置の最終的な出力とする。 The processing of the above Equations 1 to 3 is repeatedly calculated, and the “state quantity (−) / error covariance matrix (−)” 8 is calculated every sampling. For each sampling, “state quantity (−) / error covariance matrix (−)” 8 is set as the final output of the target state quantity estimation apparatus.
このように、本実施の形態に係る目標状態量推定装置によれば、RCSが小さく、観測値のS/N比が低い目標の運動を推定する状況にあっても、フィルタ特性をその場に応じて変化させることができる。
その結果、例えば、誘導開始直後では定常特性を重視した特性、会合直前では過渡特性を重視した特性というように、その場で求められる性能指標を満足したフィルタ特性を得ることができる。これにより、目標の位置や速度、旋回加速度等の状態量を高精度に推定することができ、ひいてはその推定された状態量を用いて、目標を高精度に追尾することができる。
As described above, according to the target state quantity estimating device according to the present embodiment, even in a situation where the motion of a target with a small RCS and a low S / N ratio of an observed value is estimated, the filter characteristics can be placed on the spot. It can be changed accordingly.
As a result, it is possible to obtain a filter characteristic that satisfies the performance index required on the spot, such as a characteristic that emphasizes the steady characteristic immediately after the start of guidance and a characteristic that emphasizes the transient characteristic immediately before the meeting. As a result, state quantities such as the target position, speed, and turning acceleration can be estimated with high accuracy, and the target can be tracked with high accuracy by using the estimated state quantities.
なお、この目標状態量推定装置を用いて、目標追尾装置を構成することが可能となる。 In addition, it becomes possible to comprise a target tracking apparatus using this target state quantity estimation apparatus.
1 シーカ、2 観測更新計算装置、3 時間外挿計算装置、4 観測ノイズ行列計算装置、5 プロセスノイズ行列計算装置、6 観測量、7 S/N比、8 状態量(−)・誤差共分散行列(−)、9 観測ノイズ行列、10 状態量(+)・誤差共分散行列(+)、11 プロセスノイズ行列、100 目標状態量推定装置 1 Seeker, 2 Observation update calculator, 3 Time extrapolation calculator, 4 Observation noise matrix calculator, 5 Process noise matrix calculator, 6 Observed quantity, 7 S / N ratio, 8 State quantity (-) and error covariance Matrix (−), 9 observation noise matrix, 10 state quantity (+) / error covariance matrix (+), 11 process noise matrix, 100 target state quantity estimation device
Claims (5)
飛しょう体に搭載され、目標までの距離と方向を観測した観測量と、前記観測量のS/N比を出力するシーカと、
前記観測量と前記S/N比に基づき、前記観測量の誤差の大きさを表す観測ノイズ行列を算出し、前記観測ノイズ行列を観測更新計算装置に出力する観測ノイズ行列計算装置と、
前記観測量と、前記観測ノイズ行列と、時間外挿計算装置が出力する目標の状態量(−)と誤差共分散行列(−)に基づき観測更新計算を行い、観測更新計算後の目標の状態量(+)と誤差共分散行列(+)を出力する観測更新計算装置と、
更新計算後の前記目標の状態量(+)と誤差共分散行列(+)と、前記目標の運動モデルの誤差の大きさであるプロセスノイズ行列に基づき、時間外挿計算により目標の状態量(−)と誤差共分散行列(−)を算出し、出力する時間外挿計算装置と、
前記時間外挿計算装置が算出した前記目標の状態量(−)に基づき、前記プロセスノイズ行列を算出し、算出した前記プロセスノイズ行列を前記時間外挿計算装置に出力するプロセスノイズ行列計算装置と、
から構成されることを特徴とする目標状態量推定装置。 A target state quantity estimation device for estimating a target state quantity by a Kalman filter,
An observable that is mounted on a flying object and that observes the distance and direction to the target, and a seeker that outputs the S / N ratio of the observable,
An observation noise matrix calculation device that calculates an observation noise matrix that represents a magnitude of an error of the observation amount based on the observation amount and the S / N ratio, and outputs the observation noise matrix to an observation update calculation device;
The observation update calculation is performed based on the observation amount, the observation noise matrix, the target state quantity (−) and the error covariance matrix (−) output by the time extrapolation calculation device, and the target state after the observation update calculation is performed. An observation update calculation device that outputs a quantity (+) and an error covariance matrix (+);
Based on the state quantity (+) and error covariance matrix (+) of the target after the update calculation, and the process noise matrix that is the magnitude of the error of the target motion model, the target state quantity ( -) And a time extrapolation calculation device that calculates and outputs an error covariance matrix (-),
A process noise matrix calculation device that calculates the process noise matrix based on the target state quantity (−) calculated by the time extrapolation calculation device, and outputs the calculated process noise matrix to the time extrapolation calculation device; ,
A target state quantity estimation device comprising:
ここで
プロセスノイズ行列:Q
プロセスノイズ行列初期値:Q0
ゲイン行列関数:G
である。 2. The process noise matrix is a matrix represented by the following equation, and the process noise matrix initial value Q0 is an initial value set in advance based on the target maximum turning acceleration. Target state quantity estimation device.
Where process noise matrix: Q
Process noise matrix initial value: Q0
Gain matrix function: G
It is.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2015043108A JP6248967B2 (en) | 2015-03-05 | 2015-03-05 | Target state quantity estimation device |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2015043108A JP6248967B2 (en) | 2015-03-05 | 2015-03-05 | Target state quantity estimation device |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2016161269A true JP2016161269A (en) | 2016-09-05 |
JP6248967B2 JP6248967B2 (en) | 2017-12-20 |
Family
ID=56846715
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2015043108A Active JP6248967B2 (en) | 2015-03-05 | 2015-03-05 | Target state quantity estimation device |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP6248967B2 (en) |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH08338698A (en) * | 1995-06-14 | 1996-12-24 | Toshiba Corp | Radio wave induction control device |
JPH11218398A (en) * | 1998-02-02 | 1999-08-10 | Mitsubishi Electric Corp | Guidance controller for missile |
US20050035900A1 (en) * | 2003-08-14 | 2005-02-17 | Zaugg David A. | System and method for tracking beam-aspect targets with combined Kalman and particle filters |
US7180443B1 (en) * | 2005-03-16 | 2007-02-20 | Lockheed Martin Corporation | Reduced state estimator for systems with physically bounded parameters |
JP2008076096A (en) * | 2006-09-19 | 2008-04-03 | Mitsubishi Electric Corp | Tracking method and its system |
JP2011002266A (en) * | 2009-06-17 | 2011-01-06 | Nec Corp | Target tracking processor and method for correcting error covariance matrix used therefor |
JP2011220625A (en) * | 2010-04-12 | 2011-11-04 | Mitsubishi Electric Corp | Missile guiding device |
JP2013195104A (en) * | 2012-03-16 | 2013-09-30 | Mitsubishi Electric Corp | Target state quantity estimation device and flying-object guidance control apparatus |
JP2014092840A (en) * | 2012-11-01 | 2014-05-19 | Nec Corp | Missile guidance system, missile guidance signal generation circuit, missile guidance method, and missile guidance program |
JP2015197403A (en) * | 2014-04-03 | 2015-11-09 | 三菱電機株式会社 | target tracking device |
-
2015
- 2015-03-05 JP JP2015043108A patent/JP6248967B2/en active Active
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH08338698A (en) * | 1995-06-14 | 1996-12-24 | Toshiba Corp | Radio wave induction control device |
JPH11218398A (en) * | 1998-02-02 | 1999-08-10 | Mitsubishi Electric Corp | Guidance controller for missile |
US20050035900A1 (en) * | 2003-08-14 | 2005-02-17 | Zaugg David A. | System and method for tracking beam-aspect targets with combined Kalman and particle filters |
US7180443B1 (en) * | 2005-03-16 | 2007-02-20 | Lockheed Martin Corporation | Reduced state estimator for systems with physically bounded parameters |
JP2008076096A (en) * | 2006-09-19 | 2008-04-03 | Mitsubishi Electric Corp | Tracking method and its system |
JP2011002266A (en) * | 2009-06-17 | 2011-01-06 | Nec Corp | Target tracking processor and method for correcting error covariance matrix used therefor |
JP2011220625A (en) * | 2010-04-12 | 2011-11-04 | Mitsubishi Electric Corp | Missile guiding device |
JP2013195104A (en) * | 2012-03-16 | 2013-09-30 | Mitsubishi Electric Corp | Target state quantity estimation device and flying-object guidance control apparatus |
JP2014092840A (en) * | 2012-11-01 | 2014-05-19 | Nec Corp | Missile guidance system, missile guidance signal generation circuit, missile guidance method, and missile guidance program |
JP2015197403A (en) * | 2014-04-03 | 2015-11-09 | 三菱電機株式会社 | target tracking device |
Also Published As
Publication number | Publication date |
---|---|
JP6248967B2 (en) | 2017-12-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105549049B (en) | A kind of adaptive Kalman filter algorithm applied to GPS navigation | |
KR101998236B1 (en) | Apparatus for compensating impact angle error in real-time | |
JP2015197403A (en) | target tracking device | |
RU2013101745A (en) | METHOD FOR RADAR DETECTION OF BALLISTIC OBJECTIVE MANEUVER AT PASSIVE TRAJECTORY SECTION | |
JP5577807B2 (en) | Flying object guidance device | |
RU2006101467A (en) | METHOD FOR OPERATION OF INFORMATION-COMPUTING SYSTEM OF ROCKET AND DEVICE FOR ITS IMPLEMENTATION | |
JP6248967B2 (en) | Target state quantity estimation device | |
US20050242242A1 (en) | Methods and systems for guiding an object to a target using an improved guidance law | |
Chaitanya et al. | Convergence issues of taylor series method in determining unknown target location using hyperbolic multilateration | |
JP5760772B2 (en) | Flying object guidance control device | |
Garapati Vaishnavi et al. | Underwater bearings-only tracking using particle filter | |
Kumar et al. | Target tracking using adaptive Kalman Filter | |
CN111043914B (en) | Method and system for acquiring bullet sight line angular rate applied to strapdown seeker | |
Hexner et al. | Stochastic approach to optimal guidance with uncertain intercept time | |
RU2014144061A (en) | METHOD FOR CONTROLLING AN ANTENNA INERTIAL DRIVE, PROVIDING SUSTAINABLE SUPPORT OF INTENSE MANEUVERING AND HIGH-SPEED AIRCRAFT FROM AN EXTENDED YESISTAN YEAR-RESISTANCE | |
Kumar et al. | Adaptive extended kalman filter for ballistic missile tracking | |
JP2013195104A (en) | Target state quantity estimation device and flying-object guidance control apparatus | |
Nakagawa et al. | Disturbance-compensated intercept guidance using line-of-sight rate information | |
Jahan et al. | Cannon fired ball with relative velocity | |
RU2632476C2 (en) | Method for detecting maneuver of ballistic object by sampling products of distance and radial speed and device for its implementation | |
RU2592747C1 (en) | Method of firing at moving target | |
Garg et al. | Error estimation in ballistic missile trajectory using Kalman Filter | |
Kim et al. | JPDAS multi-target tracking algorithm for cluster bombs tracking | |
Kim et al. | Modified receding horizon guidance law with information on small accurate time-to-go | |
Fu et al. | Simulation Analysis on Optimal Guidance Law of Maneuvering Target Based on Multiple Model Adaptive Estimator |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20170120 |
|
TRDD | Decision of grant or rejection written | ||
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20171024 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20171026 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20171106 |
|
R151 | Written notification of patent or utility model registration |
Ref document number: 6248967 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R151 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |