JPH11500219A - 射撃統制装置 - Google Patents

Abstract

(57)【要約】 本発明は射撃統制装置、特にターゲット行路の周期的成分を予測する方法に関する。効率よく高速で予測装置をセッテングすることは、はじめに周期を予測し、次に最新の測定に基づいて周期的成分の振幅を予測することにより実現される。

Description

【発明の詳細な説明】 射撃統制装置 本発明は、少なくともターゲットセンサと、ターゲットセンサに結合した計算 手段とを具え、前記ターゲットセンサを、飛行機のようなターゲットを検出し、 測定したターゲット位置のような測定したターゲット位置に関連する信号Ｐ（ｔ ）を発生するように構成し、前記計算手段を、予測される周波数決定値Ｗ、及び ターゲットの予測される位置、速度又は加速度のようなターゲットの状態と関連 し前記信号Ｐ（ｔ）から取り出したターゲット信号Ｓ（ｔ）の周期成分の予測さ れる振幅Ａを計算するように構成した射撃統制装置に関するものである。 本発明は、飛行機のようなターゲットの測定した位置から取り出され予測され る位置、速度又は加速度のようなターゲット状態と関連するターゲット信号Ｓ（ ｔ）の周期成分の周波数決定値Ｗ及び振幅Ａを予測する方法に関するものである 。 この場合、周波数決定値Ｗはターゲット信号の周波数又は周期のいずれかを含 んでいる。ターゲット信号Ｓ（ｔ）は予測され又は測定したターゲット位置、タ ーゲット速度又はターゲットの加速度の線形な結合を含むことができる。 この形式の射撃統制装置及び方法は既知である。射撃統制装置に連結されてい る銃装置についてのリード角の計算の観点に関して、ターゲット行路に関係する 値を予測することがしばしば必要である。これらの値は予測したターゲット位置 だけでなくターゲット速度及びターゲット加速度を含んでいる。ターゲットは、 例えば補足しがたい進路をとる飛行機又はミサイルとすることができる。このよ うな場合、実際的なことから明らかなように、飛行機はウィービング行路を進行 する。ウィービング行路の周期又は周波数及び振幅が既知であれば、直後のター ゲット位置又はターゲット速度のようなターゲットの状態を正確に予測すること ができる。一方、実際においては、進路の周波数又は振幅のいずれも知られてい ない。これにもかかわらず、予測した周波数及び振幅に基づいて追従することが できる。 射撃統制装置にとっては高速での動作が必須である。従って、短い時間期間中 に限られた情報量でターゲットの進路を正確に予測することが必要となる。この 分野で既知のカルマンフィルタリング技術はこの課題を解決するのに十分ではな い。この理由は、カルマンフィルタリングは設定するのに時間がかかってしまう からである。ターゲットがウィービングパターンで進路を追従し、そのパターン の周波数が連続して変化する場合、カルマンフィルタは応答が遅すぎターゲット の進路の周波数及び振幅の両方について正確で新しい予測を直ちに発生しえない 。これは、カルマンフィルタが周波数及び振幅の先に行った予測に応じて周波数 及び振幅の両方の新しい予測を常時行うという事実から発生する。一方、振幅か ら独立して周波数を連続して高精度に予測することが考えられる。予測した周波 数に基づき、次に正確な振幅予測を行うことができる。この方法により、カルマ ンフィルタを適用する場合よりも直ちに正確な予測が得られる。 本発明の射撃統制装置は上述した課題に対する解決策を提供するものであり、 従って、本発明の射撃統制装置は、計算手段を、前記予測される周波数決定値Ｗ を少なくともほぼ連続して計算できると共に、前記予測した周波数決定値Ｗを連 続して用いて、少なくともほぼ連続して予測される振幅Ａを計算できるように構 成したことを特徴とする。 本発明の方法は、予測される周波数決定値Ｗを少なくともほぼ連続して計算し 、前記予測される周波数決定値Ｗを連続して用い、前記予測される振幅Ａを少な くともほぼ連続して計算することを特徴とする。 ターゲット進路の予測アルゴリズムにおいて、予測されるターゲット加速度を 用いることによりより高い精度を得ることができる。 従って、本発明の射撃統制装置の別の有益な実施例は、ターゲット信号Ｓ（ｔ ）がターゲットの予測される加速度を含むことを特徴とする。 従って、本発明の射撃統制装置のさらに別の有益な実施例は、ターゲット信号 Ｓ（ｔ）が予測される加速度信号を含むことを特徴とする。 ターゲット信号Ｓ（ｔ）が周波数ｗ₀の周期を有しているか否かを高速で効率 よく判別する方法は、この周波数ｗ₀について高いゲインを有しそれ以外の周波 数に対して低いゲインを有するバントパスフィルタによりターゲット信号をフィ ルタリングすることである。このフィルタは、計算時間をほとんど必要としない ために、好ましくは周波数ｗ₀において低い減衰率及び共振周波数を有する２次 フィルタで構成する。急峻なロールオフが得られる高次のフィルタで構成するこ ともできる。信号Ｓ（ｔ）の周波数は周波数ｗ₀に接近しているので、フィルタ の出力信号はより大きな絶対値を有することになる。従って、信号Ｓ（ｔ）の絶 対値と関連するフィルタ出力信号の絶対値を決定することにより、ｗ₀付近の周 波数について信号Ｓ（ｔ）中に周波数ｗ₀が含まれるか否か決定することができ る。適当なフィルタ特性を選択することにより、感度を増大させ又は低下させる ことができる。一方、広い周波数帯域についてターゲット信号がどの周波数を有 しているかチェックしようとする場合、単一のフィルタでは不十分である。この 理由は、ｗ₀付近にない周波数について信号Ｓ（ｔ）にｗ₀が含まれている否かチ ェックできないためである。この課題は、異なる周波数に対して感度を有する複 数のフィルタに同時に信号Ｓ（ｔ）を供給することにより解決できる。Ｓ（ｔ） の周期は各フィルタの出力信号を比較することにより直ちに決定することができ る。 従って、本発明の射撃統制装置の別の有益な実施例は、計算手段がＮ個のフィ ルタＦ_i、ｉ＝１，..．，Ｎ，を具え、各フィルタＦ_iがフィルタＦ_iに属する周 波数決定値Ｗ_iの信号に主として応答し、前記ターゲット信号Ｓ（ｔ）が各フィ ルタＦ_iの入力部に存在し、前記計算手段を、Ｎ個のフィルタＦ_iからの出力信号 に基づいて前記予測される周波数決定値Ｗを計算するように構成したことを特徴 とする。 従って、本発明の方法の別の有益な実施例は、ターゲット信号をＮ個のフィル タＦ_i、ｉ＝１，..．，Ｎ，に供給し、各フィルタＦ_iがフィルタＦ_iに属する周 波数決定値Ｗ_iの信号に主として応答し、Ｎ個のフィルタＦ_iからの出力信号に基 づいて前記予測する周波数決定値Ｗを決定することを特徴とする。 各フィルタの出力信号を比較できるようにするためには、各フィルタの応答周 波数、すなわち共振周波数における各フィルタのゲインを等しくすることが有益 である。 従って、本発明の射撃統制装置の別の有益な実施例は、フィルタＦ_iに属する 周波数決定値Ｗ_iを有するシヌソイダル入力信号についての各フィルタＦ_iが、 少なくともほぼ同一のゲインを有することを特徴とする。 本発明の方法の別の有益な実施例は、フィルタＦ_iに属する周波数決定値Ｗ_iを 有するシヌソイダル入力信号について各フィルタＦ_iが、少なくともほぼ同一の 信号ゲインを有することを特徴とする。 各フィルタの順次の共振周波数を決定するに際し、２個の順次の共振周波数間 で対数スケールで一定の間隔が得られるように順次の共振周波数間で一定の比率 を選択することが好ましい。このように構成することにより、比較的高い周波数 において２個の順次の共振周波数の２個の順次の信号形状間の識別が顕著になる 。 従って、本発明の射撃統制装置の別の有益な実施例は、周波数決定値Ｗ_iが対 数スケールで少なくともほぼ等間隔の周波数ω_iを含むことを特徴とする。 従って、本発明の方法の別の有益な実施例は、周波数決定値Ｗ_iが対数スケー ルで少なくともほぼ等間隔の周波数ω_iを含むことを特徴とする。 信号Ｓ（ｔ）の周波数を決定するためには、フィルタの各出力信号の振幅を決 定し、次に出力振幅が最大になるフィルタを決定することが有益である。出力振 幅信号は、出力信号の絶対値を求めその後ローパスフィルタを介して絶対値信号 を発生させることにより決定する。信号Ｓ（ｔ）のの予測される周波数はこのフ ィルタの共振周波数の式で表される。この方法により、予測される周波数は高速 で効率よく決定することができる。 従って、本発明の射撃統制装置の別の有益な実施例は、計算手段を、フィルタ Ｆ_iの出力信号の振幅Ａ_iを決定するように構成すると共に、前記周波数決定値Ｗ_i 及び振幅値Ａ_iを用いて予測される周波数決定値Ｗを決定するように構成したこ とを特徴とする。 従って、本発明の方法の別の有益な実施例は、フィルタＦ_iの出力信号の振幅 Ａ_iを決定し、前記周波数決定値Ｗ_i及び振幅値Ａ_iを用いて予測される周波数決 定値Ｗを決定することを特徴とする。 必要な場合、各フィルタの振幅間で補間を行うことができ、Ｓ（ｔ）の予測さ れる周波数を補間が全体として最大を示す周波数に等しいものとする。この結果 、より精度の高い予測が得られ、より少ない数のフィルタを用いることができ、 計算時間を節約することができる。 従って、本発明の射撃統制装置の別の有益な実施例は、計算手段が、前記周波 数決定値Ｗ_iを横座標として用い前記振幅Ａ_iを周波数決定値Ｗ_iに属する縦座標 として用いて縦座標Ａ_i間で補間を行うように構成され、前記計算手段が、前記 補間が全体として最大値を示す周波数決定値から予測する周波数決定値を取り出 すように構成したことを特徴とする。 従って、本発明の方法の別の有益な実施例は、前記周波数決定値Ｗ_iを横座標 として用い前記振幅Ａ_iを周波数決定値Ｗ_iに属する縦座標として用いて縦座標Ａ_i 間で補間を行い、補間が全体として最大値を示す周波数決定値から予測される 周波数決定値を取り出すことを特徴とする。 Ｓ（ｔ）の周波数の良好な予測は既知であり、Ｓ（ｔ）の振幅の良好な予測は この予測される周波数を用いることにより得ることができる。そのように行う明 確な方法は、最小２乗予測器によりこの周波数の位相が互いに９０°ずれた２個 のシヌソイダル信号の振幅を信号Ｓ（ｔ）に適合させることによりこれら２個の シヌソイダル信号を信号をＳ（ｔ）に適合させることである。この方法は、シヌ ソイダル信号の値は連続して計算しなければならず、計算時間が長くなる欠点が ある。予測される周波数の矩形波形を有する信号の振幅を適合させることはより 有益である。この理由は、矩形波信号の値の計算は、その計算時間が極めて僅か であるからである。最小２乗予測器により矩形波の信号をターゲット信号Ｓ（ｔ ）に適合させるためには、矩形波の正端縁部分すなわち矩形波の半分の期間中及 び矩形波の負の端縁部分すなわち矩形波の半分の期間中にターゲット信号Ｓ（ｔ ）について時間積分処理を行う必要がある。このためには、サンプルしたターゲ ット信号Ｓ（ｔ）を決定し、サンプルした矩形波をサンプルしたターゲット信号 に適合させることが有益である。さらに、ターゲット信号を位置信号のような以 前にサンプルされた信号から取り出すことも可能である。上述した積分は加算す ることにも変更される。この場合、予測されることが可能な最も低い周波数の全 期間中にターゲット信号を周期的なバッファメモリに増分的に記憶することが有 益である。すなわち、バッファメモリが完全に満たされた場合、場合メモリの最 初のメモリ素子に戻る。このようにして、各サンプリング時刻について、信号Ｓ （ｔ）のサンプリング時刻について加算処理が行われる。多数のサンプリング時 刻 により形成されるある期間にわたる加算は、その期間の終了時刻のバッファメモ リの数値からその期間の開始時刻におけるバッファメモリの数値を減算すること により直接得ることができる。バッファメモリの内容は、サンプリング時刻につ いて加算されたターゲット信号値が必要となる各計算処理に用いることができる 。 従って、本発明の射撃統制装置の別の有益な実施例は、サンプリングしたター ゲット信号Ｓ（ｔ）を決定する手段を設け、周期的なバッファメモリを設け、前 記ターゲット信号Ｓ（ｔ）のサンプリング値を前記バッファメモリの順次のメモ リ素子に増分的に記憶する手段を設けたことを特徴とする。 本発明の方法の有益な実施例は、ターゲット信号Ｓ（ｔ）がサンプリングした 信号を含み、前記ターゲット信号Ｓ（ｔ）のサンプリング値をバッファメモリの 順次のメモリ素子に増分的に記憶することを特徴とする。 増分的なターゲット信号値をバッファメモリに記憶すれば、これらの値を用い て予測される振幅Ａを決定するのに有益である。この構成を実施するためには、 単に数個の計算工程を必要とするにすぎない。 従って、本発明の射撃統制装置の別の有益な実施例は、計算手段を周期的なバ ッファメモリに記憶した値に基づいて予測される振幅Ａを決定するように構成し たことを特徴とする。 本発明の方法の別の有益な実施例は、予測される振幅を、周期的なバッファメ モリに記憶した値及び予測される周波数決定値に基づいて決定することを特徴と する。 サンプリングされた信号Ｓ（ｔ）を決定するに際し、例えばクロック信号によ り一定のサンプリング時間を設けることが有益であり、一方、信号Ｓ（ｔ）を一 定のサンプリング時間でサンプルされた位置信号から取り出すこともできる。こ れにより、サンプルされた信号と共に行われる計算処理が簡単になる。 従って、本発明の射撃統制装置の別の有益な実施例は、計算手段に、前記周期 的なバッファメモリに記憶した値に時間バリディティを割り当てる手段を設け、 順次の値の時間バリディティ間の差をほぼ一定にしたことを特徴とする。 従って、本発明の方法の別の実施例は、周期的なバッファメモリに記憶した値 が増分的なターゲット信号ｇ（ｔ_i）を形成し、これらの値が少なくともほぼ等 間隔の離散的な時刻ｔ_iに時間バリディティを有することを特徴とする。 サンプリング時刻間の増分的なターゲット信号値を決定する手段を設けること により射撃統制装置を一層改良することができる。これにより、サンプルしたタ ーゲット信号の加算を決定することができる。 従って、本発明の射撃統制装置の有益な実施例は、ｔが時間を表すものとした 場合に、前記計算手段を、前記周期的なバッファメモリに記憶されている値ｇ（ ｔ_i）に基づいて増分ターゲット信号値ｇ（ｔ）を決定するように構成し、前記 時間時刻ｔ_iの離散点間に位置することを特徴とする。 従って、本発明方法の有益な実施例は、ｔが時刻を表し、この時刻が離散的な 時刻ｔ_iとした場合に、前記周期的なバッファメモリに記憶されている値ｇ（ｔ_i ）に基づいて増分ターゲット信号値ｇ（ｔ）を決定することを特徴とする。 増分ターゲット信号値ｇを決定する方法は、２個のサンプリング時刻間で線形 補間することであり、この方法は計算容量をほとんど必要とせず十分な精度が得 られる。 本発明の射撃統制装置の別の有益な実施例は、計算手段を、前記周期的なバッ ファメモリに記憶されている増分ターゲット信号値ｇ（ｔ_i）間の線形補間によ り前記ターゲット信号値ｇ（ｔ）を決定するように構成したことを特徴とする。 本発明の方法の別の有益な実施例は、周期的なバッファメモリに記憶されてい る増分ターゲット信号値ｇ（ｔ_i）間の線形補間により前記ターゲット信号値ｇ （ｔ）を決定することを特徴とする。 バッファメモリに記憶した最新のサンプルに関連していかなる時刻の増分ター ゲット信号をも決定することができるので、予測される周期に対応する波長を有 する矩形波をサンプルされたターゲット信号に高速でいかなる瞬時においても適 合させることができる。このようにする場合、少なくともほぼ１／４の周期だけ 位相がずれた２個の矩形波を毎回適合させることができる。この理由は位相情報 もこのようにして発生するからである。バッファメモリに記憶されている基準サ ンプルは完全な１周期だけ戻った時刻の矩形波を適合する基準として作用する。 勿論、数周期後退することも可能であるが、これにはより多くのバッファ容量及 び計算時間が必要である。正の端縁側の半周期及び負の端縁側の半周期について のサンプルされたターゲット信号の半周期毎の加算は、各矩形波についてこの半 周期の終了時の全加算からスタート時の全加算を減算することにより決定するこ とができ、その場合各矩形波について２個の加算を互いに減算することができる 。この結果、各矩形波について各矩形波の予測される振幅が得られ、２個を予測 される振幅Ａに結合することができる。加算はバッファメモリにおいてすでに得 られるので、この処理は迅速に行うことができる。 従って、本発明の射撃統制装置の別の有益な実施例は、計算手段を、予測され る周期Ｔを決定すると共に、次に ａ₁＝（ｘ₁−２ｘ₃＋ｘ₅）／Ｔ ａ₂＝（ｘ₁−２ｘ₂＋２ｘ₄−ｘ₅）／Ｔ を決定するように構成する。 ここで、 ｘ₁は少なくともほぼｋ₁ｇ（ｔ_r）に等しい、 ｘ₂は少なくともほぼｋ₁ｇ（ｔ_r−Ｔ／４）に等しい、 ｘ₃は少なくともほぼｋ₁ｇ（ｔ_r−Ｔ／２）に等しい、 ｘ₄は少なくともほぼｋ₁ｇ（ｔ_r−３Ｔ／４）に等しい、 ｘ₅は少なくともほぼｋ₁ｇ（ｔ_r−Ｔ）に等しい、 また、ｔ_rは周期的バッファメモリに導入された時間パリディティを表わし、 ｋ₁は数１のようなスケーリングファクタを表わす。 従って、本発明の方法の別の有益な実施例は、予測される周期Ｔを決定すると 共に、次に ａ₁＝（ｘ₁−２ｘ₃＋ｘ₅）／Ｔ ａ₂＝（ｘ₁−２ｘ₂＋２ｘ₄−ｘ₅）／Ｔ を決定するように構成する。 ここで、 ｘ₁は少なくともほぼｋ₁ｇ（ｔ_r）に等しい、 ｘ₂は少なくともほぼｋ₁ｇ（ｔ_r−Ｔ／４）に等しい、 ｘ₃は少なくともほぼｋ₁ｇ（ｔ_r−Ｔ／２）に等しい、 ｘ₄は少なくともほぼｋ₁ｇ（ｔ_r−３Ｔ／４）に等しい、 ｘ₅は少なくともほぼｋ₁ｇ（ｔ_r−Ｔ）に等しい、 また、ｔ_rは周期的バッファメモリに導入された時間パリディティを表わし、 ｋ₁は数１のようなスケーリングファクタを表わす。 ここで、最新に導入されたサンプルを基準として用いることが有益である。こ の理由は、このようにすれば、できるだけ最新の振幅予測が得られるからである 。 本発明の射撃統制装置の別の有益な実施例は、ｔ_rが周期的なバッファメモリ に最も新しく導入された値の時間バリディティを表すことを特徴とする。 本発明の方法の別の有益な実施例は、ｔ_rが周期的なバッファメモリに最も新 しく導入された値の時間バリディティを表すことを特徴とする。 予測される振幅Ａを計算する際、１／４周期だけ位相のずれた２個の矩形波を 直交関数とすると、各矩形波についてピタゴラスの定理を予測される振幅に適用 できることを利用する。所望の場合、スケーリングファクタを適用することもで きる。 従って、本発明の射撃統制装置の有益な実施例は、予測される振幅Ａをｋ₂× √（ｋ₃ａ₁＋ｋ₃ａ₂ ²）に等しくすることを特徴とする。ここで、ｋ₂及びｋ₃は 数値１のようなスケーリングファクタを表わす。 矩形波からスタートする場合、スケーリングファクタを正しく選択することに より良好な近似としてシヌソイダル信号をターゲット信号に適合させることがで きる。ターゲット信号Ｓ（ｔ）が実際に純粋のシヌソイダルの形態を有すれば、 この方法により予測される振幅は正確なものとなる。 予測される周期Ｔを決定すると共に、次に ａ₁＝（ｘ₁−２ｘ₃＋ｘ₅）／Ｔ ａ₂＝（ｘ₁−２ｘ₂＋２ｘ₄−ｘ₅）／Ｔ を決定するように構成する。 ここで、 ｘ₁は少なくともほぼｋ₁ｇ（ｔ_r）に等しい、 ｘ₂は少なくともほぼｋ₁ｇ（ｔ_r−Ｔ／４）に等しい、 ｘ₃は少なくともほぼｋ₁ｇ（ｔ_r−Ｔ／２）に等しい、 ｘ₄は少なくともほぼｋ₁ｇ（ｔ_r−３Ｔ／４）に等しい、 ｘ₅は少なくともほぼｋ₁ｇ（ｔ_r−Ｔ）に等しい、 また、ｔ_rは周期的バッファメモリに導入された時間パリディティを表わし、 ｋ₁は数１のようなスケーリングファクタを表わす。 π／２に少なくともほぼ等しくしたことを特徴とする。 本発明の方法の別の有益な実施例は、予測される振幅Ａをｋ₂×√（ｋ₃ａ₁＋ ｋ₃ａ₂ ²）に等しくすることを特徴とする。ここで、ｋ₂及びｋ₃は数値１のよう なスケーリングファクタを表わす。 矩形波からスタートする場合、スケーリングファクタを正しく選択することに より良好な近似としてシヌソイダル信号をターゲット信号に適合させることがで きる。ターゲット信号Ｓ（ｔ）が実際に純粋のシヌソイダルの形態を有すれば、 この方法により予測される振幅は正確なものとなる。 予測される周期Ｔを決定すると共に、次に ａ₁＝（ｘ₁−２ｘ₃＋ｘ₅）／Ｔ ａ₂＝（ｘ₁−２ｘ₂＋２ｘ₄−ｘ₅）／Ｔ を決定するように構成する。 ここで、 ｘ₁は少なくともほぼｋ₁ｇ（ｔ_r）に等しい、 ｘ₂は少なくともほぼｋ₁ｇ（ｔ_r−Ｔ／４）に等しい、 ｘ₃は少なくともほぼｋ₁ｇ（ｔ_r−Ｔ／２）に等しい、 ｘ₄は少なくともほぼｋ₁ｇ（ｔ_r−３Ｔ／４）に等しい、 ｘ₅は少なくともほぼｋ₁ｇ（ｔ_r−Ｔ）に等しい、 また、ｔ_rは周期的バッファメモリに導入された時間パリディティを表わし、 ｋ₁は数１のようなスケーリングファクタを表わす。 以下図面を参照して本発明の射撃統制装置及びその方法を詳細に説明する。 図１は本発明の射撃統制装置及びこれに連結されている砲器装置の実施例を示 す。 図２はターゲット信号の周波数及び振幅が予測される射撃統制装置の一部の実 施例を示す。 図３は周波数が予測される射撃統制装置の一部の実施例を示す。 図４はサンプルされたターゲット信号及びこれらの信号に適合される矩形波を 表すグラフである。 本発明による「射撃統制装置」の実施例を図１に示す。射撃統制装置はターゲ ットセンサ１を具え、本例では追跡レーダのように図示したが、光学センサのよ うな他のセンサを用いることもできる。追跡レーダは、ウィービング行路（Weav ing path）３を追従するターゲット２を追跡するように取り付ける。ターゲット センサはターゲット位置信号４を計算「ユニットに供給する。ターゲットセサン はターゲットの速度信号及び／又はターゲット加速度信号センサ計算ユニット５ に供給することができる。計算ユニット５は方位角及び砲目線信号６を銃装置７ に供給し、銃装置７はターゲットを破壊するミサイル装置又は他のいかなる装置 で構成することができる。計算ユニット５はウィービング行路３の周期成分の周 波数を予測できるように構成する。 図２は本発明によるターゲット射撃統制装置に適用するための計算ユニット５ の実施例を示す。本例のターゲット信号４はサンプリング周期ΔＴでサンプルす る。一方、供給されるターゲット信号がすでにサンプルされたものも用いること ができる。サンプリングされたターゲット信号は、サンプリングされたターゲッ ト信号に基づいてターゲット加速度信号９を予測する評価器８に供給する。次に ターゲット加速度信号を周波数評価器１０に供給し、最新の時刻のサンプリング 点において有効なターゲット加速度信号の周波数を時間的にほぼ連続して又は少 なくとも各サンプリング周期に亘たるターゲットのコースを考慮する。循環予測 器又はＦＩＲ（有限インパルス応答）フィルタを用いることができる。ＦＩＲフ ィルタを用いる場合、過去のある時間周期を超える信号は予測結果に影響を与え ない。周波数の予測において、予測された振幅の使用は不要である。 周波数予測器は、予測した周波数及び／又は周期１１を発生し、これを振幅予 測器１２に供給する。この振幅予測器１２は、予測したターゲット加速度信号９ と予測した周期１１に基いて予測される振幅及び位相１３を決定する。次に、予 測した周期１１と振幅及び位相１３とをターゲット経路予測器１４に供給する。 この予測器１４は、予測した振幅、周波数及び位相に対応する振幅、周波数及び 位相を有する例えばシヌソイダルな経路の変化の仮定に基いて将来のターゲット 位置を予測する。一方、別の周期的な経路変化も用いることができる。 予測したターゲットのコース、並びに銃装置の位置、気象条件及び弾道データ のような射撃統制に重要な他のデータに基いて銃装置用の方位角及び砲目線値６ を発生させる。 図３は周波数予測器１０の可能な実施例を表わす。周波数予測器１０は並列動 作する複数のフィルタ１５を具え、各フィルタはそれ自身の固定された応答周波 数を有する。さらに、フィルタは２次フィルタを具えるが、高次のフィルタも採 用することができる。各フィルタについて高い選択性が選ばれる場合、多数の個 別のフィルタが必要になる。低い選択度の場合、数個の個別のフィルタが必要で あるが、これは評価器の精度を犠牲にすることになる。合理的な精度が得られる と共にフィルタの数が多過ぎない有益な折衷案は、１４〜１６個の２次フィルタ により得られ、この場合最も低い共振周波数は約０．１Hzであり最も高い共振周 波数は約１Hzである。絶対値捕獲器１６において、フィルタ１５からの出力信号 の絶対値を決定する。この絶対値信号は、リップルを平滑にするフィルタ１７に 供給する。入力信号９がある周波数を有する場合、この特有の周波数付近の共振 周波数を有するフィルタは大きな絶対値出力信号を有することになる。濾波され た絶対値信号は補間器１８に供給する。この補間器１８は、各サンプリング点毎 にフィルタの出力信号間において補間し、この補間において周波数軸が横軸とし て作用し振幅軸は縦軸として作用する。補間された信号１９をブロック２０に供 給し、このブロックにおいて補間信号１９の最大値及びこの最大値が存在する周 波数を決定する。この場合、周波数は予測した周波数となる。フィルタ１５の共 振周波数は対数スケールで等間隔で位置することが好ましい。 別の可能な実施例において、フィルタ１５をシヌソイダル信号の相関器で構成 する。この場合、絶対値捕獲器１６及びフィルタ１７は削除することができる。 固定周波数ω_iのシヌソイダル信号を用いて各フィルタ１５毎に相関がとられる ので、これらサンプリング点におけるシヌソイダル信号の値は予めアレイに記憶 することができ、これにより計算時間が節約される。フィルタ１５は相関値を補 間器１８に供給する。 別の有益な実施例において、フィルタ１５は循環形最小２乗予測器で構成し、 この予測器において各フィルタ毎にシヌソイダル信号の振幅及び位相が各フィル タについての固定した周波数によりターゲット信号に適合される。この循環形最 小２乗評価器は、留数ｒｅｓ_iの２乗の時間積分を行なうミニマイザ(minimizer) で構成する。ここで、パラメータｂ₁，ｂ₂及びｂ₃について、 res_i(t)＝ｓ(t)−b₁sin(ω_it)−b₂sin(ω_it)−ｂ₃ とする。ここで、ω_iは個々のフィルタについて固定された周波数である。ミニ マイザにおいて、指数フォーゲット関数(exponential forget function)を用い る。各フィルタについて、指数フォーゲット関数は、各フィルタの固定周波数の 周期である逆数にほぼ等しい時定数を有している。留数ｒｅｓ_i，ｉ＝１，..., Ｎ及び周波数ω_i，ｉ＝１_...,Ｎにより、例えば以下の式で表わされる品質ファ クタを決定することができる。 この品質ファクタは補間器１８に直接供給する。ここで、ｒｅｓ_i ²はフィルタ１ ５に属し、ある経過した時間期間に亘って平均化した予測留数である。これは、 好ましくは１／ω_i程度とする。ターゲット信号の周波数が約ω_iの場合、各フィ ルタは最高の品質ファクタを発生する。この場合、絶対値捕獲器１６及びフィル タ１７も削除することができる。 図４は同一の時間軸を有する４個のグラフにおいて矩形波及びオフセットをい かにして信号ｓ（ｔ）に適合させるかを示す。信号ｓ（ｔ）は周期ΔＴでサンプ ルした。一番上のグラフはサンプルされた信号ｓ（ｔ）を示し、時間ｔ_kのサン プリング点間の信号Ｓ（ｔ）は電流のホールド回路に亘って一定に維持する。３ 個の下側のグラフは矩形波信号ｑ₁（ｔ）、ｑ₂（ｔ）及びオフセットｑ₃（ｔ） をそれぞれ示す。これらの矩形波は時刻ｔ_kのサンプリング点から開始するよう に適合され、時間的に戻った矩形波の周期Ｔに亘って最新のサンプルに等しくす る。４個のグラフの細い縦の線により示すように、周期Ｔはサンプリング時間Δ Ｔの倍数に必ずしも等しくはない。矩形波の振幅ａ₁及びａ₂並びにオフセット値 ａ₃は最小２乗の基準により適合している。 一般的に、パラメータａ₁，ａ₂及びａ₃を有するパラメータ化された信号ｑ（ ｔ）は、最小２乗基準をａ₁，ａ₂及びａ₃に最小化することにより信号Ｓ（ｔ） に適合することができる。 これを本発明による射撃統制装置に用いる。本発明の実施例において、ｑ（ｔ ）は以下の３個のパラメータ化した信号の和で構成する。 ｑ（ｔ）＝ａ₁ｑ₁（ｔ）＋ａ₂ｑ₂（ｔ）＋ａ₃ｑ₃（ｔ） 信号ｑ₁（ｔ）及びｑ₂（ｔ）は、図４に示すように、１に等しい振幅において 周期的な矩形波を有すると共に先に予測した周期Ｔに等しい周期を有する。信号 ｑ₁（ｔ）及びｑ₂（ｔ）は１／４周期だけ位相がずれている。信号ｑ₃（ｔ）は 時間に対して一定であり１に等しい。最小２乗基準の最小化により決定したパラ メータａ₁及びｑ₂はＳ（ｔ）の振幅の目安であり、パラメータａ₃は基準に比例 するＳ（ｔ）のオフセットの目安となる。本発明の実施例においては、信号Ｓ（ ｔ）のサンプリングの時間有効性に等しいある時刻ｔ_rから、信号ｑ（ｔ）が１ 個の完全な周期に亘って信号Ｓ（ｔ）に適合する。所望の場合、１個の完全な周 期Ｔよりも長い後退が実際的に可能である。より長い時間周期の戻りを考慮すれ ば、振幅予測器のより強いダンペニング効果（dampening effect）が得られるこ とになる。 振幅予測法の別の利点は、ダンペニング効果が予測される周期Ｔに依存するこ とである。この理由は、より短い周期Ｔの場合後退はより長い周期の場合よりも 時間的に長くならないからである。周期Ｔが短いことは、ターゲット行路に急激 な変化があることを意味する。急激に変化する場合、振幅予測器を急激に更新す ることが望ましく、これは上述した方法により実現される。 最小２乗基準を証明する。 上式は以下の式より最小にされる。 従って、パラメータａ₁，ａ₂及びａ₃は上記積分を解くことにより得られる。 ここで、信号ｑ_m（ｔ），ｎ＝１，...,３は値１又は−１だけをとる事実を利用 することが一般的である。 従って、信号を利用することが可能になる。 さらに、ｔ₀は固定された開始時刻であり、ｔ_kは終了時刻である。ｓ（ｔ_i） がサンプリングされ零次ホールド回路により維持されるので、トラペゾイド式を 用いることが可能であるが、高次積分アルブリズムを用いることもできる。ｚ（ ｔ_k）は、ｋ＝０，１，２，..．，についてメモリ素子のアレイに記憶すること ができ、循環して計算することができる。この理由は、ｚ(ｔ_k+1)＝ｚ（ｔ_k）＋ ｓ（ｔ_k）となるからである。 このようにして、ｓ（ｔ_k）がメモリ素子のアレイに増分的に記憶される。ｔ が離散的なサンプリング時刻（ｔ_k）と(ｔ_k+1)との間に位置する場合のｇ（ｔ） の値は、これらサンプリング時刻間の線形補間により計算することができる。ま た、ｇ（ｔ）を最も近い時刻の値ｇ（ｔ_k）に等しいものとすることができ、こ の場合精度が若干低下する。 式（１），（２）及び（３）の積分は、ｑ_n（ｔ）＝１の部分とｑ_n（ｔ）＝− １の部分の副積分に分割することができる。これらの副積分は好ましく補間によ り計算した信号ｇ（ｔ）を用いることにより決定することができ、ｔについて副 積分毎の境界が挿入される。従って、以下の式が得られる。 ａ₁＝(１/T)(ｇ(t_r)−２ｇ(t_r−T/2)＋ｇ(t_r−T)) ａ₂＝(１/T)(ｇ(t_r)−２ｇ(t_r−T/4)＋２ｇ(t_r−3T/4)−g(t_r−T) ａ₃＝(１/T)(ｇ(t_r)−ｇ(t_r−T)) これより、時刻ｔ_rに対して１周期Ｔよりも僅かに長い期間戻った時刻におけ る値ｚ（ｔ_k）をメモリ素子アレイに記憶する必要がある。そのために、周期的 なバッファメモリを用いることができる。バッファメモリが完全に満たされ場合 、最初のバッファ素子に戻り、１周期Ｔ以上前の時刻に位置する以前の値ｚ（ｔ_k ）の上に書き込む。バッファメモリの数を規定する際、振幅を予測できるだけ の最長の周期Ｔを考慮する必要がある。 上述した方法は、僅かな加算と乗算が必要になるだけであり、処理速度に対し て有利になる利点がある。シヌソイダル信号を適合させることが望まれる場合、 振幅予測についてサンプリング時刻におけるサイン値の計算が必要となり、この 場合上述した方法よりも長い時間がかかってしまう。

───────────────────────────────────────────────────── フロントページの続き (81)指定国 ＥＰ(ＡＴ，ＢＥ，ＣＨ，ＤＥ， ＤＫ，ＥＳ，ＦＲ，ＧＢ，ＧＲ，ＩＥ，ＩＴ，ＬＵ，Ｍ Ｃ，ＮＬ，ＰＴ，ＳＥ)，ＡＭ，ＡＵ，ＢＢ，ＢＧ，Ｂ Ｒ，ＢＹ，ＣＡ，ＣＮ，ＣＺ，ＥＥ，ＧＥ，ＨＵ，ＩＳ ，ＪＰ，ＫＥ，ＫＧ，ＫＰ，ＫＲ，ＫＺ，ＬＫ，ＬＲ， ＬＴ，ＬＶ，ＭＤ，ＭＧ，ＭＮ，ＭＷ，ＭＸ，ＮＯ，Ｎ Ｚ，ＰＬ，ＲＯ，ＲＵ，ＳＤ，ＳＧ，ＳＩ，ＳＫ，ＴＪ ，ＴＭ，ＴＲ，ＴＴ，ＵＡ，ＵＧ，ＵＳ，ＵＺ，ＶＮ

Claims (1)

1. 【特許請求の範囲】 １．少なくともターゲットセンサと、ターゲットセンサに結合した計算手段とを 具え、前記ターゲットセンサを、飛行機のようなターゲットを検出し、測定した ターゲット位置のような測定したターゲット位置に関連する信号Ｐ（ｔ）を発生 するように構成し、前記計算手段を、予測される周波数決定値Ｗ、及びターゲッ トの予測される位置、速度又は加速度のようなターゲットの状態と関連し前記信 号Ｐ（ｔ）から取り出したターゲット信号Ｓ（ｔ）の周期成分の予測される振幅 Ａを計算するように構成した射撃統制装置において、前記計算手段を、前記予測 される周波数決定値Ｗを少なくともほぼ連続して計算できると共に、前記予測し た周波数決定値Ｗを連続して用いて、少なくともほぼ連続して予測される振幅Ａ を計算できるように構成したことを特徴とする射撃統制装置。 ２．請求項１に記載の射撃統制装置において、前記ターゲット信号Ｓ（ｔ）が予 測したターゲットの加速度を含むことを特徴とする射撃統制装置。 ３．請求項１又は２に記載の射撃統制装置において、前記計算手段がＮ個のフィ ルタＦ_i、ｉ＝１，..．，Ｎ，を具え、各フィルタＦ_iがフィルタＦ_iに属する周 波数決定値Ｗ_iの信号に主として応答し、前記ターゲット信号Ｓ（ｔ）が各フィ ルタＦ_iの入力部に存在し、前記計算手段を、Ｎ個のフィルタＦ_iからの出力信号 に基づいて前記予測した周波数決定値Ｗを計算するように構成したことを特徴と する射撃統制装置。 ４．請求項３に記載の射撃統制装置において、フィルタＦ_iに属する周波数決定 値Ｗ_iを有するシヌソイダル入力信号についての各フィルタＦ_iが、少なくともほ ぼ同一のゲインを有することを特徴とする射撃統制装置。 ５．請求項３又は４に記載の射撃統制装置において、前記周波数決定値Ｗ_iが、 少なくとも対数スケールでほぼ等間隔の複数の周波数ω_iを含むことを特徴とす る射撃統制装置。 ６．請求項３から５までのいずれか１項に記載の射撃統制装置において、前記計 算手段を、フィルタＦ_iの出力信号の振幅Ａ_iを決定するように構成すると共に、 前記周波数決定値Ｗ_i及び振幅値Ａ_iを用いて予測される周波数決定値Ｗ を決定するように構成したことを特徴とする射撃統制装置。 ７．請求項６に記載の射撃統制装置において、前記計算手段が、前記周波数決定 値Ｗ_iを横座標として用い前記振幅Ａ_iを周波数決定値Ｗ_iに属する縦座標として 用いて縦座標Ａ_i間で補間を行うように構成され、前記計算手段が、前記補間が 全体として最大値を示す周波数決定値から予測する周波数決定値を取り出すよう に構成したことを特徴とする射撃統制装置。 ８．請求項１から７までのいずれ１項に記載の射撃統制装置において、サンプリ ングしたターゲット信号Ｓ（ｔ）を決定する手段を設け、周期的なバッファメモ リを設け、前記ターゲット信号Ｓ（ｔ）のサンプリング値を前記バッファメモリ の順次のメモリ素子に増分的に記憶する手段を設けたことを特徴とする射撃統制 装置。 ９．請求項８に記載の射撃統制装置において、前記計算手段を、前記周期的なバ ッファメモリに記憶した値に基づいて予測される振幅Ａを決定するように構成し たことを特徴とする射撃統制装置。 10．請求項９に記載の射撃統制装置において、前記計算手段に、前記周期的なバ ッファメモリに記憶した値に時間バリディティを割り当てる手段を設け、順次の 値の時間バリディティ間の差をほぼ一定にしたことを特徴とする射撃統制装置。 11．請求項１０に記載の射撃統制装置において、ｔが時間を表すものとした場合 に、前記計算手段を、前記周期的なバッファメモリに記憶されている値ｇ（ｔ_i ）に基づいて増分ターゲット信号値ｇ（ｔ）を決定するように構成し、前記時間 時刻ｔ_iの離散点間に位置することを特徴とする射撃統制装置。 12．請求項１１に記載の射撃統制装置において、前記計算手段を、前記周期的な バッファメモリに記憶されている増分ターゲット信号値ｇ（ｔ_i）間の線形補間 により前記ターゲット信号値ｇ（ｔ）を決定するように構成したことを特徴とす る射撃統制装置。 13．請求項１１又は１２に記載の射撃統制装置において、前記計算手段を、予測 した周期Ｔを決定すると共に、 ｘ₁をｋ₁ｇ（ｔ_r）に少なくともほぼ等くし、 ｘ₂をｋ₁ｇ（ｔ_r−Ｔ／４）に少なくともほぼ等くし、 ｘ₃をｋ₁ｇ（ｔ_r−Ｔ／２）に少なくともほぼ等くし、 ｘ₄をｋ₁ｇ（ｔ_r−３Ｔ／４）に少なくともほぼ等くし、 ｘ₅をｋ₁ｇ（ｔ_r−Ｔ）に少なくともほぼ等くし、 ｔ_rが前記周期的バッファメモリに導入された時間パリディティを表わし、ｋ_Lが 数値１のようなスケーリングファクタを表わすものとした場合に、 ａ₁＝（ｘ₁−２ｘ₃＋ｘ₅）／Ｔ ａ₂＝（ｘ₁−２ｘ₂＋２ｘ₄−ｘ₅）／Ｔ の値を決定するように構成したことを特徴とする射撃統制装置。 14．請求項１３に記載の射撃統制装置において、ｔ_r前記周期的なバッファメモ リに記憶した最も新しい値の時間バリディティを表すことを特徴とする射撃統制 装置。 15．請求項１３又は１４に記載の射撃統制装置において、前記計算手段を、ｋ₂ 及びｋ₃を数値１のようなスケーリングファクタを表すものとした場合に、前記 予測した振幅Ａをを式 ｋ₂×√（ｋ₃ａ₁ ²＋ｋ₃ａ₂ ²） で表すように構成したことを特徴とする射撃統制装置。 ２に少なくともほぼ等しくしたことを特徴とする射撃統制装置。 17．飛行機のようなターゲットの測定した位置から取り出され予測される位置、 速度又は加速度のようなターゲット状態と関連するターゲット信号Ｓ（ｔ）の周 期成分の周波数決定値Ｗ及び振幅Ａを予測するに当たり、予測される周波数決定 値Ｗを少なくともほぼ連続して計算し、前記予測される周波数決定値Ｗを連続し て用い、前記予測される振幅Ａを少なくともほぼ連続して計算することを特徴と する方法。 18．請求項１７に記載の方法において、前記ターゲット信号Ｓ（ｔ）が予測した 加速度信号を含むことを特徴とする方法。 19．請求項１８に記載の方法において、前記ターゲット信号をＮ個のフィルタＦ_i 、ｉ＝１，..．，Ｎ，に供給し、各フィルタＦ_iがフィルタＦ_iに属する周波数 決定値Ｗ_iの信号に主として応答し、Ｎ個のフィルタＦ_iからの出力信号に基づい て前記予測される周波数決定値Ｗを決定することを特徴とする方法。 20．請求項１９に記載の方法において、フィルタＦ_iに属する周波数決定値Ｗ_iを 有するシヌソイダル入力信号についての各フィルタＦ_iが、少なくともほぼ同一 の信号ゲインを有することを特徴とする方法。 21．請求項１９又は２０に記載の方法において、前記周波数決定値Ｗ_iが、少な くとも対数スケールでほぼ等間隔の複数の周波数ω_iを含むことを特徴とする方 法。 22．請求項１９から２１までのいずれか１項に記載の方法において、フィルタＦ_i の出力信号の振幅Ａ_iを決定し、周波数決定値Ｗ_i及び振幅値Ａ_iを用いて予測さ れる周波数決定値Ｗを決定することを特徴とする方法。 23．請求項２２に記載の方法において、前記周波数決定値を横座標として用い前 記振幅Ａ_iを前記周波数決定値Ｗ_iに属する縦座標として用いて縦座標Ａ_i間で補 間を行い、前記補間が全体として最大値を示す周波数決定値から予測される周波 数決定値Ｗを取り出すことを特徴とする方法。 24．請求項１７から２３までのいずれ１項に記載の方法において、前記ターゲッ ト信号Ｓ（ｔ）がサンプリングした信号を含み、前記ターゲット信号Ｓ（ｔ）の サンプリング値をバッファメモリの順次のメモリ素子に増分的に記憶することを 特徴とする方法。 25．請求項２４に記載の方法において、前記周期的なバッファメモリに記憶した 値に基づいて予測した振幅Ａを決定することを特徴とする方法。 26．請求項２５に記載の方法において、前記周期的なバッファメモリに記憶した 値が増分的なターゲット信号ｇ（ｔ_i）を形成し、これらの値が少なくともほぼ 等間隔の離散的な時刻ｔ_iに時間バリディティを有することを特徴とする方法。 27．請求項２６に記載の方法において、ｔが時刻を表し、この時刻が離散的な時 刻ｔ₁とした場合に、前記周期的なバッファメモリに記憶されている値ｇ（ｔ_i ）に基づいて増分ターゲット信号値ｇ（ｔ）を決定するすることを特徴とする 方法。 28．請求項２７に記載の方法において、前記周期的なバッファメモリに記憶され ている増分ターゲット信号値ｇ（ｔ_i）間の線形補間により前記ターゲット信号 値ｇ（ｔ）を決定することを特徴とする方法。 29．請求項２７又は２８に記載の方法において、前記計算手段を、予測した周期 Ｔを決定すると共に、 ｘ₁をｋ₁ｇ（ｔ_r）に少なくともほぼ等くし、 ｘ₂をｋ₁ｇ（ｔ_r−Ｔ／４）に少なくともほぼ等くし、 ｘ₃をｋ₁ｇ（ｔ_r−Ｔ／２）に少なくともほぼ等くし、 ｘ₄をｋ₁ｇ（ｔ_r−３Ｔ／４）に少なくともほぼ等くし、 ｘ₅をｋ₁ｇ（ｔ_r−Ｔ）に少なくともほぼ等くし、 ｔ_rが前記周期的バッファメモリに導入された時間パリディティを表わし、ｋ_L が数値１のようなスケーリングファクタを表わすものとした場合に、 ａ₁＝（ｘ₁−２ｘ₃＋ｘ₅）／Ｔ ａ₂＝（ｘ₁−２ｘ₂＋２ｘ₄−ｘ₅）／Ｔ の値を決定するように構成したことを特徴とする射撃統制装置。 30．請求項２９に記載の方法において、ｔ_rが前記周期的なバッファメモリに記 憶した最も新しい値の時間バリディティを表すことを特徴とする方法。 31．請求項２９又は３０に記載の方法において、ｋ₂及びｋ₃を数値１のようなス ケーリングファクタを表すものとした場合に、前記予測した振幅Ａをを式ｋ₂× √（ｋ₃ａ₁ ²＋ｋ₃ａ₂ ²）で表すことを特徴とする方法。 くともほぼ等しくしたことを特徴とする方法。