(57)【要約】【課題】アンテナで受信した電界信号の位相を用いずに、デコンボリューション法とプロニー法を用いて電波源情報を直接求めることができる観測装置、観測方法、及びプログラムを提供すること。【解決手段】観測装置１０は、所定の方位角内でビームが走査されるアンテナからの受信電力信号の平方根信号及び走査されるビームの方位角の信号に基づき、平方根信号を方位角に関してフーリエ変換し、フーリエ変換した第１の方位角周波数信号を、アンテナのアンテナ電力パターンの平方根信号又はアンテナ電界パターン信号を方位角に関してフーリエ変換した第２の方位角周波数信号で除算し、除算した信号を、引数に実数部と虚数部とを有する指数関数でプロニー法によりフィッティングする信号処理部１５を具備する。【選択図】図１
2020060543.tif 000002
【特許請求の範囲】
【請求項１】
所定の方位角内でビームが走査されるアンテナからの受信電力信号の平方根信号及び前記走査されるビームの方位角の信号に基づき、前記平方根信号を方位角に関してフーリエ変換し、前記フーリエ変換した第１の方位角周波数信号を、前記アンテナのアンテナ電力パターンの平方根信号又はアンテナ電界パターン信号を前記方位角に関してフーリエ変換した第２の方位角周波数信号で除算し、前記除算した信号を、引数に実数部と虚数部とを有する指数関数でプロニー法によりフィッティングする信号処理部
を具備する観測装置。
【請求項２】
前記信号処理部は、前記フィッティングにより得られた電波源候補のうち、前記実数部の絶対値が所定値以下の電波源候補を抽出する
請求項１に記載の観測装置。
【請求項３】
前記信号処理部は、前記抽出した電波源候補から実際の電波源を特定するための情報を、前記指数関数から取得する
請求項２に記載の観測装置。
【請求項４】
前記抽出した電波源候補から実際の電波源を特定するための情報は、前記電波源候補の相対振幅の値、実数部の値及び方位角の値である
請求項３に記載の観測装置。
【請求項５】
前記抽出された電波源候補及び前記実際の電波源を特定するための情報を表示する表示部
を更に具備する、請求項３又は４に記載の観測装置。
【請求項６】
前記信号処理部は、プロニー法によりフィッティングする前に、前記除算した信号にローパスフィルタ処理を施す
請求項１から５のうちいずれか１項に記載の観測装置。
【請求項７】
前記所定の方位角内でビームが走査されるアンテナとして、ビームアンテナ及び当該ビームアンテナを回転するアンテナ回転装置で構成されたアンテナシステム、電子的アンテナビーム走査装置を有するアンテナシステム及びビームアンテナ及び当該ビームアンテナを移動するアンテナ移動装置で構成されたアンテナシステムのうち、いずれかのアンテナシステム
を更に具備する、請求項１から６のうちいずれか１項に記載の観測装置。
【請求項８】
前記アンテナに送信電力を供給する送信部と、
送受信の切り替えを行うスイッチ部と
を更に具備する、請求項１から７のうちいずれか１項に記載の観測装置。
【請求項９】
所定の方位角内でビームが走査されるアンテナからの受信電力信号の平方根信号及び前記走査されるビームの方位角の信号に基づき、前記平方根信号を方位角に関してフーリエ変換するステップと、
前記フーリエ変換した第１の方位角周波数信号を、前記アンテナのアンテナ電力パターンの平方根信号又はアンテナ電界パターン信号を方位角に関してフーリエ変換した第２の方位角周波数信号で除算するステップと、
前記除算した信号を、引数に実数部と虚数部とを有する指数関数でプロニー法によりフィッティングするステップと、
を備える、観測方法。
【請求項１０】
所定の方位角内でビームが走査されるアンテナからの受信電力信号の平方根信号及び前記走査されるビームの方位角の信号に基づき、前記受信電力信号の平方根信号を方位角に関してフーリエ変換するステップと、
前記フーリエ変換した第１の方位角周波数信号を、前記アンテナのアンテナ電力パターンの平方根信号又はアンテナ電界パターン信号を方位角に関してフーリエ変換した第２の方位角周波数信号で除算するステップと、
前記除算した信号を、引数に実数部と虚数部とを有する指数関数でプロニー法によりフィッティングするステップと
をコンピュータに実行させるプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、例えばレーダによる観測に使われる観測装置、観測方法及びプログラムに関する。
【背景技術】
【０００２】
例えばレーダにおいて観測対象物を観測する場合は、アンテナの回転等によってアンテナビームを回転させながら、アンテナで受信した電力信号の強さを観測する方法が用いられている。
【０００３】
この際、回転範囲内で受信信号のピーク点をサーチすることで、電波源情報(位置（方位角）、強度)を求めることが一般的に実施されており、直接電波源情報を抽出することができないという問題がある。本発明者は、デコンボリューション法とプロニー法を用いて電波源情報を直接求める手法を提唱した（特許文献１及び非特許文献１参照）。この手法は理論的に電波源情報を求める方法であり、数値シミュレーションで電波源情報が得られることを確認している。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】特開２０００−１１４８５０号公報
【非特許文献】
【０００５】
【非特許文献１】ＳＥＺＡＩ，Ｔ．"Ｉｍｐｒｏｖｅｍｅｎｔ  ｏｆ  Ｄｉｒｅｃｔｉｏｎ  Ｒｅｓｏｌｕｔｉｏｎ  ｏｆ  Ｒｏｔａｔｉｏｎａｌ  Ｒａｄａｒ  ｂｙ  Ｐｒｏｎｙ  Ａｎａｌｙｓｉｓ  ０ｆ  Ｄｅｃｏｎｖｏｌｖｅｄ    Ａｎｔｅｎｎａ  Ｏｕｔｐｕｔ"  ．  電子情報通信学会技術研究報告  （ＩＥＩＣＥ  ＴＥＣＨＮＩＣＡＬ  ＲＥＰＯＲＴ）［ＩＳＳＮ：０９１３−５６８５］．  Ｖｏｌｕｍｅ  １１７，Ｎｕｍｂｅｒ  １８２．  ２０１７，０８，１７．  Ｐａｇｅｓ  ３７−４２＊Ｓｅｃｔｉｏｎｓ  ２−４＊
【発明の概要】
【発明が解決しようとする課題】
【０００６】
本発明者が提唱した手法では、アンテナで受信した電界信号とアンテナの電界パターン信号を使用する。電界信号及び電界パターン信号は複素数信号であるため、電界信号及び電界パターン信号の振幅と位相の取得が必要となる。従来の観測装置においては、振幅値の取得は一般的であるが、位相については取得していない場合がほとんどである。  
アンテナパターンの測定では、汎用設備を使用して測定することが多いため、位相の取得は困難ではない。  
これに対して、例えばレーダや通信のようにアンテナに到来する電波の周波数が正確にわかっている場合、アンテナで受信した電界信号の位相を取得するためには、観測装置の受信部を改修し、位相取得機能を追加するか、新しい受信部を準備することが必要になる。例えば電波の探索のようにアンテナに到来する電波の周波数が正確にわかっていない場合、アンテナで受信した電界信号の位相を取得するためには、装置用アンテナに隣接してリファレンスアンテナを設置し、リファレンスアンテナで受信した電波を基準にして位相を取得することになる。このように、位相取得機能を追加することにより受信部、あるいは装置全体が複雑になり、寸法、重量、コストが増大するという課題があった。
【０００７】
このような事情に鑑み、本発明の目的は、アンテナで受信した電界信号の位相を用いずに、デコンボリューション法とプロニー法を用いて電波源情報を直接求めることができる観測装置、観測方法、及びプログラムを提供することにある。
【課題を解決するための手段】
【０００８】
上記問題を解決するため、本発明の一形態に係る観測装置は、所定の方位角内でビームが走査されるアンテナからの受信電力信号の平方根信号及び前記走査されるビームの方位角の信号に基づき、前記平方根信号を方位角に関してフーリエ変換し、前記フーリエ変換した第１の方位角周波数信号を、前記アンテナのアンテナ電力パターンの平方根信号又はアンテナ電界パターン信号を方位角に関してフーリエ変換した第２の方位角周波数信号で除算し、前記除算した信号を、引数に実数部と虚数部とを有する指数関数でプロニー法によりフィッティングする信号処理部を具備する。
【０００９】
デコンボリューション法処理の過程で得られる電波源分布方位角周波数関数は、電波源の分布関数を方位角に関してフーリエ変換したものである。電波源の分布は複数の離散電波源の和であるとモデル化することが可能である。離散した電波源はそれぞれを点波源と考えることができるため、それをフーリエ変換したものは方位角周波数上では振幅が一定の指数関数で表現される。従って、電波源分布方位角周波数関数は振幅が一定の指数関数の和として表現される。  
それ故、特許文献１で開示した手法、すなわち電波源分布方位角周波数関数をプロニー法により指数関数でフィッティングすることにより、電波源分布に関する情報を直接取得することが可能である。また、非特許文献１で開示した手法により、電波源分布方位角周波数関数から指数関数成分を抽出する際に、引数の実数部の値により、指数関数の識別を行い、雑音に起因する指数関数を除外することが可能である。  
電波源分布関数は方位角の関数で、振幅は複素数である。しかし、電波源観測の際に位相情報が必要となることはほとんどない。そのため、アンテナからの受信電力信号の平方根信号及びアンテナのアンテナ電力パターンの平方根信号又はアンテナ電界パターン信号を使用して特許文献１及び非特許文献１の手法で処理をすることにより、電波源分布関数の振幅を実数とみなした場合の結果を取得しても、実用上、支障となることはない。
【００１０】
前記信号処理部は、前記フィッティングにより得られた電波源候補のうち、前記実数部の絶対値が所定値以下の電波源候補を抽出することが好ましい。
【００１１】
前記信号処理部は、前記抽出した電波源候補から実際の電波源を特定するための情報を前記指数関数から取得することが好ましい。
【００１２】
前記抽出した電波源候補から実際の電波源を特定するための情報は、前記電波源候補の相対振幅の値、実数部の値及び方位角の値であることがより好ましい。
【００１３】
本発明に係る観測装置は、前記抽出された電波源候補及び前記実際の電波源を特定するための情報を表示する表示部を更に具備することが好ましい。
【００１４】
前記信号処理部は、プロニー法によりフィッティングする前に、前記除算した信号にローパスフィルタ処理を施すことが好ましい。
【００１５】
本発明に係る観測装置は、前記所定の方位角内でビームが走査されるアンテナとして、ビームアンテナ及び当該ビームアンテナを回転するアンテナ回転装置で構成されたアンテナシステム、電子的アンテナビーム走査装置を有するアンテナシステム及びビームアンテナ及び当該ビームアンテナを移動するアンテナ移動装置で構成されたアンテナシステムのうち、いずれかのアンテナシステムを更に具備することが好ましい。
【００１６】
本発明に係る観測装置は、前記アンテナに送信電力を供給する送信部と、送受信の切り替えを行うスイッチ部とを更に具備することが好ましい。
【００１７】
本発明の一形態に係る観測方法は、所定の方位角内でビームが走査されるアンテナからの受信電力信号の平方根信号及び前記走査されるビームの方位角の信号に基づき前記受信電力信号の平方根信号を方位角に関してフーリエ変換するステップと、前記フーリエ変換した第１の方位角周波数信号を、前記アンテナのアンテナ電力パターンの平方根信号又はアンテナ電界パターン信号を方位角に関してフーリエ変換した第２の方位角周波数信号で除算するステップと、前記除算した信号を、引数に実数部と虚数部とを有する指数関数でプロニー法によりフィッティングするステップとを備える。
【００１８】
本発明の一形態に係るプログラムは、所定の方位角内でビームが走査されるアンテナからの受信電力信号の平方根信号及び前記走査されるビームの方位角の信号に基づき、前記受信電力信号の平方根信号を方位角に関してフーリエ変換するステップと、前記フーリエ変換した第１の方位角周波数信号を、前記アンテナのアンテナ電力パターンの平方根信号又はアンテナ電界パターン信号を方位角に関してフーリエ変換した第２の方位角周波数信号で除算するステップと、前記除算した信号を、引数に実数部と虚数部とを有する指数関数でプロニー法によりフィッティングするステップとをコンピュータに実行させる。
【発明の効果】
【００１９】
本発明により、アンテナで受信した電界信号の位相を用いずに、デコンボリューション法とプロニー法を用いて電波源情報を直接求めることができる。
【図面の簡単な説明】
【００２０】
【図１】本発明の一実施形態に係る観測装置を示すブロック図である。
【図２】図１に示した観測装置の信号処理部の構成を示すブロック図である。
【図３】本発明の一実施形態に係る電波源情報のシミュレーション結果の一例を示すグラフである。
【図４】非特許文献１に示された手法に係る電波源情報のシミュレーション結果の一例を示すグラフである。
【発明を実施するための形態】
【００２１】
以下、図面を参照しながら、本発明の実施形態を説明する。  
図１は、本発明に一実施形態に係る観測装置を示すブロック図である。この実施形態は、本発明に係る観測装置をレーダ装置に適用したものである。  
図１に示すように、観測装置１０は、アンテナシステム１１と、送信部１２と、受信部１３と、スイッチ部１４と、信号処理部１５と、表示部１６とを具備する。  
この観測装置１０は、レーダ装置であり、対象物に向けて電波を送信し、その反射波を観測することにより、対象物の電波源情報を測る装置である。
【００２２】
アンテナシステム１１は、所定の方位内でアンテナビームを走査するものであり、ビームアンテナ１１１及び当該ビームアンテナ１１１を回転するアンテナ回転装置１１２で構成される。アンテナ回転装置１１２は、ビームアンテナ１１１の分解能を向上させる方向にビームアンテナ１１１を回転させる。アンテナシステム１１は、アンテナ回転装置１１２の回転に応じて走査されるビームの方位角の信号を出力する。ビームアンテナとしては、ダイポール系のアンテナ、ホーンアンテナ、パラボラアンテナ、アレイアンテナなど如何なるアンテナでも用いることができる。  
送信部１２は、ビームアンテナ１１１に送信電力を供給する。  
受信部１３は、ビームアンテナ１１１の受信電力信号をその平方根信号に変換する。  
スイッチ部１４は、送受信の切り替えを行う。  
信号処理部１５は、スイッチ部１４を介して受信部１３から受信電力信号の平方根信号及びアンテナシステム１１から方位角の信号を入力し、所定の信号処理を施し、電波源候補、その電波源候補の相対強度、後述する指数関数引数の実数部の値及び位置（方位角）に関する情報などを出力する。  
表示部１６は、信号処理部１５からの出力に基づき電波源候補の相対強度、指数関数引数の実数部の値及び位置（方位角）に関する情報の相対強度及び方位角などに関する情報を表示する。
【００２３】
図２は、上記の信号処理部１５の構成を示すブロック図である。  
図２に示すように、信号処理部１５は、フーリエ変換部１５１と、アンテナパターン信号保持部１５２と、除算部１５３と、ローパスフィルタ部１５４と、指数関数抽出部１５５と、指数関数選択部１５６と、電波源情報取得部１５７とを有する。
【００２４】
フーリエ変換部１５１は、入力した受信電力信号の平方根信号ｅ（θ）を方位角に関してフーリエ変換した方位角周波数信号Ｅ（ω）を求める。ここで、θは方位角、ωは方位角周波数である。
アンテナパターン信号保持部１５２は、ビームアンテナ１１１のレーダ装置で使用するアンテナ電力パターンの平方根信号ｇ（θ）を方位角に関してフーリエ変換した方位角周波数信号Ｇ（ω）を保持する。  
除算部１５３は、方位角周波数信号Ｅ（ω）を、方位角周波数信号Ｇ（ω）で除算する。除算値は、電波源分布を方位角に関してフーリエ変換した方位角周波数信号Ｆ（ω）である。すなわち、除算部１５３は、  
Ｆ（ω）＝Ｅ（ω）／Ｇ（ω）
を求める。
【００２５】
ローパスフィルタ部１５４は、除算部１５３からの出力信号に対してローパスフィルタ処理を行う。アンテナ電力パターンの平方根信号をフーリエ変換した信号には、物理的に方位角周波数上の帯域が存在するため、その帯域を超えた成分は持ち得ないが、除算部１５３で信号処理を行って取得した除算信号（電波源分布方位角周波数信号）には、上記帯域を超えた方位角周波数成分が含まれる。これは実際の処理と理想的な処理との差により発生する。この帯域を超えた信号成分を、そのまま用いて次段の処理を行うと、信号を劣化させることになる。従って、ここでローパスフィルタ処理を行い、上記帯域を超えた信号成分を除去し、信号の劣化を防止している。
【００２６】
指数関数抽出部１５５は、除算部１５３で除算し、ローパスフィルタ処理を行った信号を、引数に実数部αｉと虚数部ｊθｉとを有する指数関数（ｅｘｐ（αｉ−ｊθｉ））でプロニー法によりフィッティングする。  
指数関数選択部１５６は、予めαｉの大きさに閾値を設け、閾値以下のαｉを有する指数関数（電波源候補）を選択する。  
電波源情報取得部１５７は、指数関数選択部１５６で選択した指数関数から電波源情報を取得する。電波源情報は、各電波源候補の相対強度、位置（方位角）の値である。
【００２７】
次に、このように構成した観測装置１０の動作を説明する。  
アンテナ回転装置１１２によりビームアンテナ１１１を回転させながら、送信部１２からビームアンテナ１１１に送信電力を供給し、ビームアンテナ１１１から電波を送信する。
【００２８】
ビームアンテナ１１１から送信された電波が散乱体（対象物）で反射されて戻ってくると、ビームアンテナ１１１は受信電波を出力する。スイッチ部１４を介した受信電波は受信部１３で受信電力信号の平方根信号に変換され、次いで信号処理部１５に入力される。
【００２９】
信号処理部１５では、以下の信号処理が実行される。  
・受信電力信号の平方根信号ｅ（θ）を方位角に関してフーリエ変換した方位角周波数信号Ｅ（ω）を求める。  
・この方位角周波数信号Ｅ（ω）を、ビームアンテナ１１１のレーダ装置におけるアンテナ電力パターンの平方根信号ｇ（θ）を方位角に関してフーリエ変換した方位角周波数信号Ｇ（ω）で除算する。  
・除算信号に対してローパスフィルタ処理を行う。  
・ローパスフィルタ処理を行った信号を、引数に実数部αｉと虚数部ｊθｉとを有する指数関数でプロニー法によりフィッティングする。  
・予めαｉの大きさに閾値を設け、閾値以下のαｉを有する指数関数（電波源候補）を選択する。  
・選択した指数関数から電波源情報を取得する。  
各電波源候補の相対強度、位置（方位角）の値が表示部１６により表示される。表示部１６は表形式などでこれらの情報を表示する。
【００３０】
次に、本発明に係る効果を確認するために行った数値シミュレーション結果を説明する。  
本発明に係るアンテナとしてアンテナ長が波長の５０．４２倍の一様分布の開口面アンテナ（アンテナビーム幅：１度）を用い、このアンテナを−４５度から＋４５度まで回転させながら電波を送信し、アンテナから等距離上の−０．５度及び＋０．５度の方向に位置する電波的に同一の点散乱体で散乱され、このアンテナに戻ってくる電波を信号対雑音比が２０ｄＢの環境で取得した際に得られる電波情報源を検討した。
【００３１】
ローパスフィルタ処理は方位角周波数における絶対値で例えば０．８（度−１）内の信号を通過させる。  
ローパスフィルタ処理を行ったＦ（ω）を、引数に実数部αｉと虚数部ｊθｉとを有する指数関数でプロニー法によりフィッティングし、実数部αｉの閾値が絶対値で０．０１以下の指数関数を選択し、その指数関数の電波源情報、及び指数関数引数の実数部αｉ(電波源強度の上位5情報)を表１に示す。
【００３２】
【表１】
2020060543.tif 000003
また、シミュレーションで得られた電波情報源を図示した結果を図３に示す。  
これにより、本発明により得られる電波源情報より、実際の電波源に近い２つの方位にほぼ同じ強度の電波源が存在することが分かる。また、これらの方向から離れた方向に多数の電波源の存在を示す情報も示されるが、強度が約１３ｄＢ以上小さいため、雑音によるものと判断することができる。
【００３３】
比較例として、非特許文献１示された手法により、上記と同様の数値シミュレーションを行った結果を表２と図４に示す。
【００３４】
【表２】
2020060543.tif 000004
表１と表２を比較すると、本発明で得られる電波源の強度は非特許文献１で得られる強度とほぼ同じである。本発明で得られる電波源の方位の精度は非特許文献１で得られる電波源の方位の精度よりも悪い。図３と図４を比較すると、本発明では電波源以外の出力の数は、非特許文献１よりも多くなる。このように本発明は非特許文献１と比較すると、電波源情報を取得する性能は低下するが、方位誤差は最悪でも０．０３５度と高精度であり、電波源以外による出力は雑音によるものと容易に判断できる。よって、本発明においても、電波源情報は精度よく取得できる。
【００３５】
本発明は、上記の実施形態には限定されず、様々な変形が可能であり、それらも本発明の技術的思想の範囲に属する。  
上記の実施形態では、アンテナパターン信号保持部１５２がビームアンテナ１１１のレーダ装置で使用するアンテナ電力パターンの平方根信号ｇ（θ）を方位角に関してフーリエ変換した方位角周波数信号Ｇ（ω）を保持し、除算部１５３が方位角周波数信号Ｅ（ω）を、方位角周波数信号Ｇ（ω）で除算し、方位角周波数信号Ｆ（ω）を求めていたが、本発明では、アンテナ電力パターンの平方根信号ｇ（θ）に代えてアンテナのアンテナ電界パターン信号を用いてもよい。
【００３６】
上記の実施形態においては、アンテナビームの指向方向を、観測電波源方向に移動させるための手段として、機械的なアンテナ回転装置を用いたものを示したが、アンテナビームの指向方向を移動させる手段としては、電子的なアンテナビーム走査手段や、航空機、人工衛星等のような移動プラットフォームに搭載してアンテナビームを移動させる手段等を用いることができる。
【００３７】
また、上記の実施形態は、観測装置がアンテナシステムや表示部などを有するものであったが、観測装置が外部のアンテナシステムから受信電力信号や角度信号を入力する構成であっても構わない。また、観測装置が外部の表示装置に情報を出力し、表示装置に情報を表示させる構成であってもよい。また、観測装置は情報を表示することなく、指数関数から電波源情報に基づき所定の判断処理を実行してもよい。
【００３８】
また、上記の実施形態では、アンテナパターン信号保持部がビームアンテナのアンテナ電力パターンの平方根信号ｇ（θ）を方位角に関してフーリエ変換した方位角周波数信号Ｇ（ω）を保持するものであったが、ビームアンテナのアンテナ電力パターンの平方根信号ｇ（θ）が動的に変動するような場合にはその都度入力し、方位角周波数信号Ｇ（ω）を算出するように構成してもよい。
【００３９】
上記の実施形態で例示した実数部αｉの閾値は適応的に可変するものであってもよい。
【符号の説明】
【００４０】
１０  観測装置
１１  アンテナシステム
１２  送信部
１３  受信部
１４  スイッチ部
１５  信号処理部
１６  表示部
１５１  フーリエ変換部
１５２  アンテナパターン信号保持部
１５３  除算部
１５４  ローパスフィルタ部
１５５  指数関数抽出部
１５６  指数関数選択部
１５７  電波源情報取得部
【図１】
2020060543.tif 000005
【図２】
2020060543.tif 000006
【図３】
2020060543.tif 000007
【図４】
2020060543.tif 000008