(57)【要約】【課題】ＦＭ−ＣＷレーダ装置において、レーダの受信信号に混入したインパルス状の干渉波をビート信号から除去する。【解決手段】ビート信号を連続する３以上の被調整・注目・被調整スイープデータ列とし、被調整スイープデータ列の振幅と位相を前記注目スイープデータ列の振幅と位相に等しくなるようにそれぞれ調整して調整済スイープデータ列を得、調整済・注目・調整済スイープデータ列の対応するデータの中央値を出力用スイープデータ列とする。【選択図】      図３
2010025901.tif 000002
【特許請求の範囲】
【請求項１】
所定の変調周期で時間の経過とともに線形的に周波数変調した連続波の送信信号を送信し、空間的に隔てた物標にて反射された信号を含む受信信号を受信して、送信信号と受信信号の周波数差によって物標までの距離を求められることを利用したＦＭ−ＣＷレーダ装置において、
送信信号と受信信号の周波数差を周波数成分とするビート信号を規定の標本化クロックに同期してディジタル信号に変換するＡＤ変換部と、
前記ＡＤ変換部によって標本化されたビート信号を、前記変調周期毎のスイープデータ列として、連続するＮのスイープデータ列（Ｎは、３以上の整数）を格納するスイープバッファと、
前記スイープバッファに格納された連続するＮのスイープデータ列の１つを注目スイープデータ列とし、該注目スイープデータ列以外の複数Ｎ−１のスイープデータ列を被調整スイープデータ列とし、該Ｎ−１の被調整スイープデータ列の振幅と位相を前記注目スイープデータ列の振幅と位相に等しくなるようにそれぞれ調整してＮ−１の調整済スイープデータ列を得る振幅・位相調整部と、
前記振幅・位相調整部から出力された前記注目スイープデータ列及び前記Ｎ−１の調整済スイープデータ列の対応するデータの実部と虚部をそれぞれ小さい順に並べ替え、実部の中央値を実部とし虚部の中央値を虚部とする複素数を１スイープデータ列に含まれるすべてのサンプリング時刻のデータについて出力するセレクタ回路とを備えて、
前記注目スイープデータ列に含まれる妨害信号を低減することを特徴とする、ＦＭ−ＣＷレーダ装置。
【請求項２】
前記振幅・位相調整部は、
前記スイープバッファに格納された連続する注目スイープデータ列及びＮ−１の被調整スイープデータ列をそれぞれフーリエ変換するフーリエ変換部と、
前記フーリエ変換部による連続する注目スイープデータ列及びＮ−１の被調整スイープデータ列のフーリエ変換結果を、前記スイープバッファの注目スイープデータ列及びＮ−１の被調整スイープデータ列と対応付けできるように、格納するフーリエバッファと、
前記フーリエバッファに格納された前記注目スイープデータ列のフーリエ変換結果の複素絶対値を計算する複素絶対値演算部と、
前記複素絶対値演算部での複素絶対値演算によって計算された注目スイープデータ列のフーリエ変換結果の複素絶対値を一時的に格納する複素絶対値バッファと、
前記複素絶対値バッファに格納された複素絶対値を探索してその複素絶対値が最大になる振幅最大周波数を前記フーリエ変換部の周波数分解能よりも細かい周波数分解能で算出する最大値検出部と、
前記フーリエバッファに格納された連続する注目スイープデータ列及びＮ−１の被調整スイープデータ列のフーリエ変換結果に対して前記振幅最大周波数に対応する各関数値を計算する複素関数補間処理部と、
前記複素関数補間処理部で計算された各関数値を用いて、前記Ｎ−１の被調整スイープデータ列の振幅と位相が、前記注目スイープデータ列の振幅と位相に一致するように補正するビート信号補正部と、を有することを特徴とする、請求項１に記載のＦＭ−ＣＷレーダ装置。
【請求項３】
前記ビート信号補正部は、
前記複素関数補間処理部で計算された関数値の内の前記Ｎ−１の被調整スイープデータ列にそれぞれ対応する関数値の逆数をそれぞれに計算するＮ−１の被調整スイープ用複素逆数回路と、
前記Ｎ−１の被調整スイープ用複素逆数回路によってそれぞれに計算された関数値の逆数と、前記注目スイープデータ列に対応する計算された関数値とをそれぞれ乗算することによって前記Ｎ−１の被調整スイープデータ列にそれぞれ対応する複素補正値を計算し、該複素補正値をそれぞれ対応する前記被調整スイープデータ列に乗算して、前記Ｎ−１の被調整スイープデータ列の振幅と位相を前記注目スイープデータ列の振幅と位相に一致させる複素乗算回路と、を有することを特徴とする、請求項２に記載のＦＭ−ＣＷレーダ装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、レーダの受信信号に混入する干渉波を除去するようにしたＦＭ−ＣＷレーダ装置に関する。
【背景技術】
【０００２】
パルスレーダ装置では、送信信号を物標で反射させて戻ってきた信号に比べ, 例えば他レーダの送信信号を直接受信した信号は非常に大きな振幅をもっているから、レーダの受信信号に混入した干渉波を除去するために、振幅比較によって干渉を除去する技術が用いられている。
【０００３】
このパルスレーダでは、連続したパルス送信タイミングに対応する受信データ列 (スイープデータ) を保持しておき、同一レンジの受信信号の振幅を連続スイープ間で比較する。その中で振幅が突出していれば干渉であると判断されるため、よく用いられる手法としては、連続スイープ間で振幅比較し、前回のスイープにより得られた受信波と今回のスイープにより得られた受信波との相関を求め、相関の低い反射波を干渉波と見なして除去するスイープ相関処理を行うことによって干渉を除去することが行われている（特許文献１）。
【特許文献１】特開平８−３６０５５号公報
【発明の開示】
【発明が解決しようとする課題】
【０００４】
しかし、ＦＭ−ＣＷレーダ装置は、パルスレーダ装置とは異なり、アンテナの指向方向に存在する物標で反射した信号が途切れることなく続いている。当然、レーダ装置の近傍に反射体が存在する状態では、受信期間の全体にわたって大振幅の信号が受信されるため、干渉が入ってきたとしても振幅だけではそれを判断できない。
【０００５】
そこで、本発明は、レーダの受信信号に混入したインパルス状の干渉波を、ビート信号から除去するようにしたＦＭ−ＣＷレーダ装置を提供することを目的とする。
【課題を解決するための手段】
【０００６】
請求項１に記載のＦＭ−ＣＷレーダ装置は、所定の変調周期で時間の経過とともに線形的に周波数変調した連続波の送信信号を送信し、空間的に隔てた物標にて反射された信号を含む受信信号を受信して、送信信号と受信信号の周波数差によって物標までの距離を求められることを利用したＦＭ−ＣＷレーダ装置において、
送信信号と受信信号の周波数差を周波数成分とするビート信号を規定の標本化クロックに同期してディジタル信号に変換するＡＤ変換部と、
前記ＡＤ変換部によって標本化されたビート信号を、前記変調周期毎のスイープデータ列として、連続するＮのスイープデータ列（Ｎは、３以上の整数）を格納するスイープバッファと、
前記スイープバッファに格納された連続するＮのスイープデータ列の１つを注目スイープデータ列とし、該注目スイープデータ列以外の複数Ｎ−１のスイープデータ列を被調整スイープデータ列とし、該Ｎ−１の被調整スイープデータ列の振幅と位相を前記注目スイープデータ列の振幅と位相に等しくなるようにそれぞれ調整してＮ−１の調整済スイープデータ列を得る振幅・位相調整部と、
前記振幅・位相調整部から出力された前記注目スイープデータ列及び前記Ｎ−１の調整済スイープデータ列の対応するデータの実部と虚部をそれぞれ小さい順に並べ替え、実部の中央値を実部とし虚部の中央値を虚部とする複素数を１スイープデータ列に含まれるすべてのサンプリング時刻のデータについて出力するセレクタ回路とを備えて、
前記注目スイープデータ列に含まれる妨害信号を低減することを特徴とする。
【０００７】
請求項２に記載のＦＭ−ＣＷレーダ装置は、請求項１に記載のＦＭ−ＣＷレーダ装置において、
前記振幅・位相調整部は、
前記スイープバッファに格納された連続する注目スイープデータ列及びＮ−１の被調整スイープデータ列をそれぞれフーリエ変換するフーリエ変換部と、
前記フーリエ変換部による連続する注目スイープデータ列及びＮ−１の被調整スイープデータ列のフーリエ変換結果を、前記スイープバッファの注目スイープデータ列及びＮ−１の被調整スイープデータ列と対応付けできるように、格納するフーリエバッファと、
前記フーリエバッファに格納された前記注目スイープデータ列のフーリエ変換結果の複素絶対値を計算する複素絶対値演算部と、
前記複素絶対値演算部での複素絶対値演算によって計算された注目スイープデータ列のフーリエ変換結果の複素絶対値を一時的に格納する複素絶対値バッファと、
前記複素絶対値バッファに格納された複素絶対値を探索してその複素絶対値が最大になる振幅最大周波数を前記フーリエ変換部の周波数分解能よりも細かい周波数分解能で算出する最大値検出部と、
前記フーリエバッファに格納された連続する注目スイープデータ列及びＮ−１の被調整スイープデータ列のフーリエ変換結果に対して前記振幅最大周波数に対応する各関数値を計算する複素関数補間処理部と、
前記複素関数補間処理部で計算された各関数値を用いて、前記Ｎ−１の被調整スイープデータ列の振幅と位相が、前記注目スイープデータ列の振幅と位相に一致するように補正するビート信号補正部と、を有することを特徴とする。
【０００８】
請求項３に記載のＦＭ−ＣＷレーダ装置は、請求項２に記載のＦＭ−ＣＷレーダ装置において、
前記ビート信号補正部は、
前記複素関数補間処理部で計算された関数値の内の前記Ｎ−１の被調整スイープデータ列にそれぞれ対応する関数値の逆数をそれぞれに計算するＮ−１の被調整スイープ用複素逆数回路と、
前記Ｎ−１の被調整スイープ用複素逆数回路によってそれぞれに計算された関数値の逆数と、前記注目スイープデータ列に対応する計算された関数値とをそれぞれ乗算することによって前記Ｎ−１の被調整スイープデータ列にそれぞれ対応する複素補正値を計算し、該複素補正値をそれぞれ対応する前記被調整スイープデータ列に乗算して、前記Ｎ−１の被調整スイープデータ列の振幅と位相を前記注目スイープデータ列の振幅と位相に一致させる複素乗算回路と、を有することを特徴とする。
【発明の効果】
【０００９】
本発明のＦＭ−ＣＷレーダ装置によれば、ビート信号を連続する３以上の被調整・注目・被調整スイープデータ列とし、複数の被調整スイープデータ列の振幅と位相を前記注目スイープデータ列の振幅と位相に等しくなるようにそれぞれ調整して調整済スイープデータ列を得、調整済・注目・調整済スイープデータ列の対応するデータの中央値を出力用スイープデータ列とするから、レーダの受信信号に混入した妨害信号を著しく低減することができる。
【発明を実施するための最良の形態】
【００１０】
本発明は、ＦＭ−ＣＷレーダ装置の受信信号にインパルス状の干渉波が混入した場合に、ビート信号からそのインパルス状干渉波分を除去しようとするものである。図１は、本発明のＦＭ−ＣＷレーダ装置の構成例を示す図である。
【００１１】
図１において、信号発生器２０１は所定の変調周期で時間の経過とともに線形的に周波数変調した連続波が線形的に変化する周波数変調信号（送信周波数信号）を生成するもので、電圧制御による発振器（ＶＣＯ）またはシンセサイザによって実現される。周波数変換器２０３は、局部発振器２０２で発生される無線周波数信号を用いて、信号発生器２０１の周波数変調信号を基本周波数帯から無線周波数帯へ周波数変換する。周波数変換器２０３で周波数変換された周波数変調信号は、電力増幅器２０４で送信電力レベルまで増幅されて、送信アンテナ２０５から送信信号として放射される。
【００１２】
受信アンテナ２０６は、空間的に隔てた物標にて反射された受信信号を受信する。この受信信号には、物標による反射信号の外に、例えばインパルス状の干渉波が混入している。送信アンテナと受信アンテナは同一アンテナでもよい。受信アンテナで受信された受信信号は、レーダの探知性能を劣化させないような低雑音増幅器２０７で増幅された後、周波数変換器２０８で周波数変換器２０３からの周波数変調信号（送信周波数信号）とを混合し、受信信号周波数と送信周波数の差の周波数をもつビート信号を取り出す。Ａ／Ｄ変換器２０８は、ビート信号を規定の標本化クロックに同期してディジタル信号（標本化されたビート信号）に変換し、直交検波回路２１０でそのビート信号をＩ信号とＱ信号に分離する。
【００１３】
妨害信号除去回路１００は本発明のＦＭ−ＣＷレーダ装置の特徴部分であり、ビート信号からインパルス状干渉波分を除去するための回路である。妨害信号除去回路１００で干渉波成分が著しく低減されたビート信号は、フーリエ変換部２１１で周波数（距離）軸における信号に変換され、更に対数変換部２１２で振幅の対数を計算し（いわゆるＬＯＧ検波）、レーダ表示装置２１３にレーダー映像を表示するために供給される。この図１において、一重矢尻の矢印は実数信号であることを示し、二重矢尻の矢印は複素信号であることを示している。
【００１４】
本発明で妨害信号除去回路１００は図１に示すようなＦＭ−ＣＷレーダの内部に組み込まれて使用され、特に、ビート信号を周波数軸(距離に対応する軸)における信号に変換するためのフーリエ変換の直前に配置される。なお、図１の例では、直交検波回路２１０はＡ／Ｄ変換の直後に配置されている。これは、論理回路またはソフトウェアによるヒルベルトフィルタによって直交検波することを想定している。直交検波は必ずしもこの位置である必要はなく、アナログによって直交検波する例も可能である。本発明は、ビート信号を周波数軸の信号に変換する前に、重畳している干渉を取り除き、その後は、通常のＦＭ−ＣＷの信号処理が実行され、レーダ表示装置にレーダ映像が表示されるようになっている。
【００１５】
図２は、本発明の基本構成を示す図である。本発明においては、標本化されたビート信号を、変調周期毎のスイープデータ列として、連続するＮのスイープデータ列（Ｎは、３以上の整数）をスイープバッファ１０に格納する。このスイープバッファ１０に格納された連続するＮのスイープデータ列の１つを注目スイープデータ列とし、その注目スイープデータ列以外の複数Ｎ−１のスイープデータ列を被調整スイープデータ列とし、振幅・位相調整部２０でこれらＮ−１の被調整スイープデータ列の振幅と位相を注目スイープデータ列の振幅と位相に等しくなるようにそれぞれ調整してＮ−１の調整済スイープデータ列を得る。そして、セレクタ回路３０で、振幅・位相調整部２０から出力された注目スイープデータ列及びＮ−１の調整済スイープデータ列の対応するデータの実部と虚部をそれぞれ小さい順に並べ替え、実部の中央値を実部とし虚部の中央値を虚部とする複素数を１スイープデータ列に含まれるすべてのサンプリング時刻のデータについて出力する。これにより、注目スイープデータ列に含まれる妨害信号を著しく低減する。
【００１６】
図３は、本発明の概念図である。図２、図３及び実施例では、例示のために、スイープデータ列のＮを３とし、中央のスイープデータ列を注目スイープとして掃引ｎと表記し、中央の１つ前のスイープデータ列を第１の被調整スイープとして掃引ｎー１と表記し、中央の１つ後のスイープデータ列を第２の被調整スイープとして掃引ｎ＋１と表記する。なお、Ｎは３以上であれば良く、また注目スイープはＮのスイープデータ列のうちのいずれか１つのスイープデータ列であればよい。
【００１７】
本発明では、連続するスイープのデータを比較することによって干渉を除去する。ただし、パルスレーダのような振幅比較による方法では不十分である。ＦＭ−ＣＷレーダは受信信号に対して周波数解析をすることによって物標までの距離を算出する方式であるので、位相情報ができるだけ壊れないようにする必要がある。
【００１８】
まず、ビート信号を互いに９０度の位相差をもつI信号とQ信号に分離しておく。これら分離された信号はＩ信号を実部としＱ信号を虚部として１つの複素信号として取り扱う。干渉除去に関しては、連続したスイープ間でＩ信号同士とＱ信号同士で比較をして干渉を判断すればよい。アンテナにはビーム幅があるため、連続するスイープの信号振幅は大きな変化はないが異なる値である。一方、物標が運動している場合、スイープ間で物標までの距離が変化することにより、ビート信号の位相が変化する。そのため、信号を比較する前にスイープ間で振幅と位相が等しくなるように信号補正する必要がある。
【００１９】
これを図３を用いて説明すると、この図は、注目スイープであるスイープｎに含まれる干渉を除去する処理について描いている。保持されている３本のスイープには異なる箇所にインパルス状の干渉が混入している。このインパルス状の干渉は、例えば他のパルスレーダから放射され、反射信号に重畳された干渉波によってビート信号に混入したものが想定される。
【００２０】
保持された３本のスイープデータはそれぞれ振幅と位相が異なるが、振幅・位相補正処理を用いてこれらをすべて同一振幅、同位相の信号に変換する。振幅・位相補正処理によって３本のスイープデータは同一の信号となるはずであり、その中に他とは異なる値が混入していた場合にそれが干渉ということである。連続する３スイープのビート信号を用いて干渉除去を実行する。ビート信号には物標までの距離に対応する周波数成分の信号に干渉が重ねあわされている。即ち、ビート信号は、一般的には、図３のように干渉を含み、しかも、物標の方位とビーム幅に関連してスイープごとに振幅がわずかに異なり、物標の相対速度のためドップラ効果が発生しビート信号の位相がスイープごとで異なる。補正された３本のスイープデータを正負の符号つきの値のまま小さい順に並べたとき、干渉のような異常値は、必ず最も大きい値か最も小さい値になっているはずである。よって、３本のスイープデータの中央値（２番目の値）を出力すれば干渉が除去される。
【００２１】
さて、振幅・位相補正処理の実現についてはフーリエ変換を用いればよい。例えば、角周波数ω0のビート信号を時刻の関数として、数式１のように表されるとする。ただし、ｉは虚数単位(すなわち、ｉ2＝−１)、Ａｏは信号振幅、φは時刻ゼロにおける信号位相 (初期位相) である。
【００２２】
【数１】
2010025901.tif 000003
【００２３】
この信号を時刻ゼロからＴにわたってフーリエ変換すると、数式２のようになる。
【００２４】
【数２】
2010025901.tif 000004
【００２５】
数式２のＦ(r)（ω）は、よく知られているようにはシンク関数と呼ばれる関数形状をしている。この関数は複素数であり、その絶対値はω＝ω０のとき、すなわち、ビート信号の周波数に対応する場所で最大値となる。その絶対値の最大値はＡ０Ｔとなり、信号振幅に比例する。また、その絶対値最大の場所における位相はφ、すなわち、ビート信号の初期位相と同じ値である。したがって、保持したスイープすべてについてフーリエ変換をした結果に対し、その絶対値が最大となる場所を求め、その場所における絶対値を調べれば振幅の補正値が計算でき, 位相を調べれば, 位相の補正値を計算できる。なお、インパルス成分はフーリエ変換されることにより全ての周波数成分に同一振幅で分散されるから、インパルス状の干渉はフーリエ変換の後には考慮しなくても良い。
【００２６】
このようなフーリエ変換による周波数解析では、ハニング窓のような窓関数を信号に乗じた後でフーリエ変換をすることが多い。窓関数なしでフーリエ変換をした場合、周波数軸上のサイドローブが大きいため、複数の周波数成分が存在する場合に他成分のサイドローブによって必要な情報が影響を受けてしまうからである。当然、ＦＭ−ＣＷレーダのビート信号には複数の周波数成分が入っていることが多いから、振幅・位相補正のためのフーリエ変換にも窓関数（例えば、数式３）をかけておいたほうがよいと思われる。
【００２７】
【数３】
2010025901.tif 000005
【００２８】
例えば、数式３で表されるハニング窓を、ビート信号ｆ（ｔ）に掛けた場合のフーリエ変換を計算してみると、数式４となる。
【００２９】
【数４】
2010025901.tif 000006
【００３０】
この場合も、フーリエ変換の絶対値はω＝ω０のときに最大となる。その最大値はＡ０Ｔ／２、そのときのフーリエ変換Ｆ(h)（ω）の位相はφである。よって、ハニング窓関数をかけたときも窓なしと同様の関係があるため、フーリエ変換を比較することによって振幅・位相補正処理が可能になる。
【００３１】
例として、連続する３スイープのビート信号を用いて干渉除去する場合を説明する。干渉除去に用いる３スイープのビート信号を、数式５で表す。
【００３２】
【数５】
2010025901.tif 000007
【００３３】
これら３スイープのビート信号は同一物標からの信号を想定しているため、ビート周波数ω0はこれらのスイープで共通した値である。振幅は、アンテナの放射パターンにより隣接するスイープ間で急激に変化しないが異なる値となる。初期位相も、物標が運動していると仮定すれば、ドップラ周波数のため掃引ごとに異なる値となる。そうなると、掃引ごとに物標までの距離が変化するのでビート周波数も微妙に変化するはずであるが、その変化は無視できるくらい小さな値である。例えば、搬送波の周波数を２４ＧＨｚ、１掃引における周波数変調時間を２ミリ秒、掃引間隔を２ミリ秒、周波数偏移幅を２０ＭＨｚとする。そのとき、速さ２０ｍ／ｓで接近する物標を見ているとする。この場合、各スイープの間に物標は４０ｍｍずつ近づいてくる。搬送波の波長が１２．５ｍｍであるので初期位相は掃引ごとに、（−４０×２）／１２．５×２π＝−４０．１９［ｒａｄ］ずつ回転する。一方、ビート周波数の掃引ごとの変化量は２．６７Ｈｚとなる。この送信条件では、ビート周波数は数ｋＨｚであり、しかも、レーダ信号処理における周波数分解能が５００Ｈｚなので、ビート周波数はほぼ一定であると考えてもよい。
【００３４】
ところで、連続した３掃引のビート信号ｆn-1（ｔ）、ｆｎ（ｔ）、ｆｎ＋１（ｔ）をフーリエ変換した結果をＦｎ−１（ω）、Ｆｎ（ω）、Ｆｎ＋１（ω）とする。上で説明したように、フーリエ変換するときにハニング窓をかけているとすると、数式６が成り立つ。
【００３５】
【数６】
2010025901.tif 000008
【００３６】
これら数式６の値を使って、数式７で表される補正した値を計算する。これらの値は、ｇｎ−１（ｔ）＝ｇｎ＋１（ｔ）＝ｆｎ（ｔ）、となるはずである。
【００３７】
【数７】
2010025901.tif 000009
【００３８】
現実的なビート信号の場合、信号に含まれる雑音や重畳した干渉信号の影響のために、厳密にはｇｎ−１（ｔ）＝ｇｎ＋１（ｔ）＝ｆｎ（ｔ）なる関係は成り立たない。しかし、干渉のように大きな外乱がない状態では、ｇｎ−１（ｔ）とｆｎ（ｔ）とｇｎ＋１（ｔ）はほぼ等しいはずである。干渉がある場合には、干渉が存在する箇所だけ、ビート信号の実部または虚部または両方が大きく異なる値となっている。つまり、ｇｎ−１（ｔ）とｆｎ（ｔ）とｇｎ＋１（ｔ）の実部と虚部を、それぞれ小さい順に並び替えたとき、最も大きい値と最も小さい値は干渉の影響を受けている可能性がある。しかし、ビート信号の同一位置において、干渉が連続掃引にわたって現れないという条件では、小さい順に並べ替えた２番目の値は干渉の影響を受けていない。すなわち、実部を小さい順に並べ替えた２番目と虚部を小さい順に並べ替えた２番目を実部および虚部とする複素数は、ビート信号から干渉を取り除いた値に他ならない。
【００３９】
以上のような干渉除去を行うための振幅・位相調整部２０の構成例が図４に、また、セレクタ回路３０の構成例が図５に示されている。
【００４０】
図４においても二重矢尻の矢印は複素信号（Ｉ信号とＱ信号で表現される）であり、一重矢尻の矢印は実数信号である。振幅・位相調整部２０はビート信号のフーリエ変換を利用して掃引ｎ−１と掃引ｎ＋１の振幅と位相を掃引ｎと一致するように調整する回路である。
【００４１】
この図４の構成では、スイープバッファ１０に保存された３スイープのビート信号それぞれをフーリエ変換部［Ｓｉｚｅ−４Ｎ  ＦＥＴ］２１−１〜２１−３に入力してフーリエ変換を計算する。ただし、このフーリエ変換部２１−１〜２１−３は、絶対値が最大になる周波数を高精度で計算するため、ビート信号の後ろに信号長の３倍の長さにゼロを追加して得られる4倍長の系列をフーリエ変換するようになっている。このようにゼロを追加して長さを拡張することで、周波数軸における計算分解能を細かくすることができ、絶対値が最大となる周波数を特定しやすくなる。
【００４２】
計算分解能を細かくするにはフーリエ変換する系列長の拡張倍率を４倍とは言わず、さらに長くするとよいと思われるが、８倍、１６倍のように倍率を上げてもそれに見合った効果は現れない。後に説明する [最大値検出部] において、２次補間処理していることによって、それほど拡張倍率を上げないうちに十分な効果が得られるのである。発明者が実施した評価によると、拡張倍率２倍では誤差が大きく、４倍で十分に収束することがわかっている。それ以上倍率を上げても、フーリエ変換の演算量の増大に見合う効果が得られないことがわかっている。よって、効果的な拡張倍率は４倍または８倍である。
【００４３】
フーリエ変換結果は、一時格納場所であるフーリエバッファ２２に供給され、掃引ｎ−１、掃引ｎ、掃引ｎ＋１のフーリエ変換結果がそれぞれバッファ２２−１、２２−２、２２−３に格納される。この実施例で保持されているスイープデータ列（掃引ｎ−１、掃引ｎ、掃引ｎ＋１）のうち、スイープｎが干渉除去をする注目スイープである。
【００４４】
注目スイープｎのフーリエ変換結果に含まれるすべての周波数に対して、絶対値計算部［ａｂｓ（）］２３によって複素絶対値を計算する。この絶対値計算部２３は複素数の絶対値、すなわち実部と虚部の自乗和の平方根を計算する演算部である。絶対値計算によって得られた実数系列は、複素絶対値バッファ２４で一時的に保持され、最大値検出部［ｐｅａｋ−ｐｏｓ（）］２５に入力される。最大値検出部２５で、フーリエ変換の絶対値が最大となる周波数ｘが計算され、出力される。この周波数ｘは、フーリエ変換で得られる離散的な周波数ではなく、補間処理によって得られた高精度な値である。最大値検出部２５が２次補間によってフーリエ変換の絶対値が最大になる周波数を離散周波数の分解能よりも細かい小数点以下の精度で周波数ｘを決定する。
【００４５】
最大値検出部２５によって得られた周波数ｘは、複素関数補間処理部［ｆ（ｘ）］２６−１〜２６−３に入力され、保存されている３スイープそれぞれに対して、周波数ｘにおけるフーリエ変換の値(複素数)を計算する。この複素関数補間処理部［ｆ（ｘ）］２６−１〜２６−２では、フーリエバッファ２２から与えられたデータ系列に対して補間演算によって、小数点以下の精度で指定された周波数ｘに対する高精度の複素量、即ち関数値、を演算する。
【００４６】
複素関数補間処理部２６−１の関数値ｚ−＝Ａ−ｅｘｐｉθ−は複素逆数回路２７−１で逆数に変換され、複素乗算回路２８−１で、複素関数補間処理部２６−２の関数値ｚ＝Ａ−ｅｘｐｉθと乗算され、第１被調整スイープデータ列である掃引ｎ−１に対応する複素補正値ｚ／ｚ−＝（Ａ／Ａ−）ｅｘｐｉ（θ−θ−）を得る。同様に、複素関数補間処理部２６−３の関数値ｚ＋＝Ａ＋ｅｘｐｉθ＋は複素逆数回路２７−３で逆数に変換され、複素乗算回路２８−３で、複素関数補間処理部２６−２の関数値ｚ＝Ａｅｘｐｉθと乗算され、第２被調整スイープデータ列である掃引ｎ＋１に対応する複素補正値ｚ／ｚ＋＝（Ａ／Ａ＋）ｅｘｐｉ（θ−θ＋）を得る。
【００４７】
複素補正値ｚ／ｚ−＝（Ａ／Ａ−）ｅｘｐｉ（θ−θ−）は、複素乗算回路２９−１で、第１被調整スイープデータ列である掃引ｎ−１に乗算されて、振幅・位相補正された第１調整済スイープデータ列である掃引ｎ−１（補正値）を得る。同様に、複素補正値ｚ／ｚ＋＝（Ａ／Ａ＋）ｅｘｐｉ（θ−θ＋）は、複素乗算回路２９−３で、第２被調整スイープデータ列である掃引ｎ＋１に乗算されて、振幅・位相補正された第２調整済スイープデータ列である掃引ｎ＋１（補正値）を得る。この部分の演算は、数式７に記述したとおりの演算になっている。これら複素逆数回路２７−１，２７−３及び複素乗算回路２８−１〜２９−３が、ビート信号補正部を構成する。
【００４８】
図５のセレクタ回路３０は、入力された複素信号の実数Ｉを小さい順に並べるソート回路３１と、入力された複素信号の虚数Ｑを小さい順に並べるソート回路３２とを備えている。
【００４９】
振幅・位相補正された３つのスイープ掃引ｎ−１（補正値）、掃引ｎ、掃引ｎ＋１（補正値）がセレクタ３０に入力され、３スイープの実部と虚部をそれぞれ比較し、中央値を出力する。したがって、セレクタ３０から出力されるスイープデータは、注目スイープｎから干渉を取り除いた値になっている。
【００５０】
図６は、振幅・位相調整部２０からフーリエ変換回路を削減した構成例を示しており、図５とはフーリエ変換回路が異なるのみで他の構成は同じである。図５で説明した例では、スイープバッファ１０に保持されたすべてのビート信号のフーリエ変換を計算するためにフーリエ変換回路が３つ実装されている。しかし、図６では、ビート信号の入力に同期して、最近のビート信号をフーリエ変換回路に入力し、その結果を保持するようにして、振幅・位相調整部２０のフーリエ変換回路を１つ実装する構成としている。
【００５１】
本発明の振幅・位相調整部２０に含まれる最大値検出部２３は、フーリエ変換結果Ｆ（ω）の絶対値が最大になる周波数ωを計算する。しかし、実際の計算においては、周波数ωは整数値しかとることができず、フーリエ変換Ｆ（ω）はその離散的な波数ωに対する値しか計算されていない。しかし、フーリエ変換の絶対値を最大にする角周波数とはビート信号に含まれる信号の周波数であるため、ビート信号の位相を正確に合わせる目的上、角周波数の特定は小数点以下の精度で計算することが望ましい。
【００５２】
そのために行う以下の計算は、離散的周波数において、フーリエ変換の絶対値が最大となる付近でフーリエ変換の絶対値を２次補間することによるアルゴリズムを用いている。そのため、最大値探索に用いるフーリエ変換は単純にビート系列ｆ（ｔ）（ただし、ｔ＝１，２，・・・Ｎ−１)をフーリエ変換した値を使うよりも、例えばｆ（Ｎ）〜ｆ（４Ｎ−１）＝０のようにもとの系列の長さの３倍の長さのゼロを後ろに追加してデータ数を４倍にする４倍ゼロ拡張した系列を用いたほうが、２次補間の性能がよい。そのように４倍ゼロ拡張されたビート系列を用いた場合、フーリエ変換に現れる周波数はω＝０，１，２，・・・４Ｎ−１のような整数値をとる。しかしながら、ビート信号の周波数は整数値であるとは限らないから、小数部を含む精度で周波数を計算する必要がある。以下に、高精度で周波数を特定するための手順を、その処理概念を示す図７を参照して説明する。
【００５３】
手順1．フーリエ変換の絶対値｜Ｆ（ω）｜を最大にする整数周波数ω＝ωmを探索する。
【００５４】
手順２．探索された整数周波数について、ωm＝０ならばωm←１、またωm＝４Ｎ−１ならばωm←４Ｎ−２のように整数周波数を補正する。なお、ωm＝０、４Ｎ−１にならないように、例えばビート信号に適切な特性の帯域通過フィルタを通過させるなどの前処理を施すことが良い。
【００５５】
手順３．次の数式８の代入式によって係数α２，α１を求める。この係数はω＝ωmの近傍において、｜Ｆ（ω）｜＝α２（ω−ωm）2＋α１（ω−ωm）＋｜Ｆ（ωm）｜のような２次補間をするための係数である。
【００５６】
【数８】
2010025901.tif 000010
【００５７】
手順４．二次補間された式から、｜Ｆ（ω）｜を最大にする周波数ω０は、次の代入式［ω０←ωｍ−α１／（２・α２）］によって計算される。なお、既に説明したように、フーリエバッファに保存されるフーリエ変換結果は、４倍ゼロ拡張された系列に対してフーリエ変換した結果である。
【００５８】
周波数の単位について補足すると、ここでの周波数は長さ４Ｎの系列をフーリエ変換したときの周波数分解能を単位とする周波数で表現されている。ビート信号入力における標本化周波数ｆｓとすれば、周波数分解能はｆｓ／４Ｎとなる。よって、ここで記述した周波数ωをＭＫＳＡ単位(すなわちＨｚ)に変換した値をｆｂとするなら、ｆｂ＝ｆｓ・ω／４Ｎ、によって計算することができる。
【００５９】
また、振幅・位相調整部２０に含まれる複素関数補間処理部［ｆ（ｘ）］２６−１〜２６−３は、小数点以下の精度で指定された周波数ｘについて、複素数で表現されるフーリエ変換の結果を計算する処理であり、この処理は、ビート信号の振幅・位相補正のために必要である。この処理も最大値検出部２５と同様、対象とする周波数の近傍に存在する離散的周波数３点を用いて２次補間することにより、正確なフーリエ変換の値Ｆ（ω）を算出する。以下に、その手順を説明する。
【００６０】
手順１．対応する周波数ω(小数部含む)に最も近い整数周波数ωmを得る。具体的には、ωm←ｆｌｏｏｒ（ω＋０．５）とすればよい。ここで、ｆｌｏｏｒ（ω）はω以下である最大の整数を得る演算である。
【００６１】
手順２．整数化された周波数ωmについて、ωm＝０ならばωm←１、またωm＝４Ｎ−１ならばωm←４Ｎ−２のように整数周波数を補正する。最大値検出部の手順２と同様である。
【００６２】
手順３．次の数式９の代入式によって係数α２，α１を求める。この係数はω＝ωmの近傍において、Ｆ（ω）＝α２（ω−ωm）2＋α１（ω−ωm）＋Ｆ（ωm）のような２次補間をするための係数である。なお、最大値検出部とは異なり、ここで求める係数α２，α１は複素数である。
【００６３】
【数９】
2010025901.tif 000011
【００６４】
手順４．手順３に記述したように、求めた係数を用いて、以下の代入式、Ｆ（ω）←α２（ω−ωm）2＋α１（ω−ωm）＋Ｆ（ωm）、によって任意周波数に対するフーリエ変換Ｆ（ω）を計算する。なお、周波数の単位については、最大値検出部における補足事項と同じ。
【００６５】
本発明は、レーダ物標が単一の場合に効果を奏することはもちろんであるが、レーダ物標が複数の場合においても効果を得ることができる。以下、本発明によって得ることができる効果を、計算機シミュレーションによって示す。
【００６６】
このシミュレーションにおいて、表１に示す３つの物標を想定した。これらの物標は、振幅、距離(ビート周波数)および相対速度(ドップラ周波数)が異なるものとしている。
【００６７】
表１
振幅       ビート周波数      ドップラ周波数
[x2p rad/SWEEP]    [rad/SWEEP]
Target 1          1.0                10.8             1.0
Target 2         0.95                45.7             0.7
Target 3         0.03               141.4           −0.3
【００６８】
表１にしたがって作成したビート信号に図８（ａ）（ｂ）（ｃ）のような干渉を重畳させた。図８のグラフは、横軸が時刻(サンプリング周期を単位とする)、縦軸がビート信号の波動関数(電圧)である。この図８（ａ）（ｂ）（ｃ）には３スイープのビート信号(Iチャンネルのみ)を描いているが、そのうち、上から２番目の図８（ｂ）のスイープが干渉除去の対象となる注目スイープである。各スイープには２〜３サンプル幅のインパルス状干渉が重畳している。
【００６９】
本発明によって注目スイープを干渉除去すると、図８（ｄ）の結果を得ることができる。この図８（ｄ）のグラフは、図８（ｂ）の干渉を含むグラフから干渉が除去されており、除去された干渉以外の場所ではほとんど重なり合っている。
【００７０】
図９のグラフは、最終的なレーダ出力における本発明の効果を示すグラフである。ＦＭ−ＣＷレーダの出力はビート信号をフーリエ変換した結果である。このグラフは、横軸をビート周波数(距離に対応する)、縦軸を信号強度(何らかの正規化をした値、単位はデシベル)である。この図は、「本発明の適用なし」と「本発明の適用あり」を重ねて表示している。本発明を適用しない場合、インパルス性の干渉によって雑音レベルが非常に高くなっている。この現象は、デルタ関数のフーリエ変換が１になることを考えると理解できるはずである。それに対し、本発明を適用した場合、雑音レベルには約２０ｄＢの改善が見られている。正確に言うと、干渉によって雑音レベルが２０ｄＢ劣化していたということである。特に、ｔａｒｇｅｔ３が干渉除去をしなければ、雑音に埋もれて見えなくなっているのに対し、干渉除去をすることによって十分に見えるようになっている。同一物標からの反射信号の電力はよく知られているように距離の4乗に反比例するから、２０ｄＢの雑音劣化は、レーダの最大探知距離を３１．６％まで劣化させてしまう。本発明は、そのような著しい性能劣化の原因となる干渉を効果的に取り除くことができる。
【図面の簡単な説明】
【００７１】
【図１】本発明のＦＭ−ＣＷレーダ装置の構成例を示す図
【図２】本発明の干渉除去の基本構成を示す図
【図３】本発明の干渉除去の概念を示す図
【図４】振幅・位相調整部２０の構成例を示す図
【図５】セレクタ回路３０の構成例を示す図
【図６】振幅・位相調整部２０のフーリエ変換回路を削減した構成例を示す図
【図７】最大値検出部の処理概念を示す図
【図８】本発明の干渉除去効果を説明するためのビート信号を示す図
【図９】本発明の効果を説明するためのシミュレーション結果を示す図
【符号の説明】
【００７２】
１００：妨害信号除去回路、１０：スイープバッファ、２０：振幅・位相調整部、
３０：セレクタ回路、２１，２１−１〜２１−３：フーリエ変換部（フーリエ変換回路）、
２２，２２−１〜２２−３：フーリエバッファ、２３：複素絶対値演算部（ａｂｓ（））、
２４：複素絶対値バッファ、２５：最大値検出部（ｐｅａｋ−ｐｏｓ（））、
２６−１〜２６−３：複素関数補間処理部（ｆ（ｘ））、
２１−１，２７−３：複素逆数回路、２８−１〜２９−１：複素乗算回路、
３１、３２：ソート回路
【図１】
2010025901.tif 000012
【図２】
2010025901.tif 000013
【図３】
2010025901.tif 000014
【図４】
2010025901.tif 000015
【図５】
2010025901.tif 000016
【図６】
2010025901.tif 000017
【図７】
2010025901.tif 000018
【図８】
2010025901.tif 000019
【図９】
2010025901.tif 000020