(57)【特許請求の範囲】
【請求項１】
連続波からなる送信信号を、放射自在な送信アンテナ、第１チャンネルから第Ｋチャンネルまでの複数のアンテナ素子からなる受信アンテナ、前記複数のアンテナ素子で受信される受信信号と前記送信信号をミキシングして前記複数のアンテナ素子に対応した複数チャンネル分のビート信号を得るミキサ、前記ミキサで得られたビート信号を所定のサンプリング周波数でサンプリングして前記複数のアンテナ素子に対応した複数チャンネル分の受信データを得る、Ａ／Ｄ変換器、前記サンプリングされた前記複数チャンネル分の受信データを、各チャンネルについて時間方向に複数の短時間データに切り出す、短時間データ切出し部、前記各チャンネルについて、前記複数の短時間データの周波数スペクトルを算出する、周波数スペクトル算出部、前記周波数スペクトルから、干渉波の干渉成分周波数を検出する、干渉周波数検出部、検出された各チャンネルの干渉成分周波数における周波数スペクトル成分に基づいて、位相補正テーブルを算出する位相補正テーブル算出部、前記位相補正テーブルに基づいて前記受信データ中の干渉方向成分を抑圧する干渉成分除去手段、該干渉方向成分が抑圧された前記受信データに基づいて物標の距離、相対速度などを検出する、物標検出部、を有する、電子走査式レーダ装置において、
前記位相補正テーブル算出部は、
前記干渉周波数検出部で得られる前記干渉成分周波数における周波数スペクトル成分で構成される複素数ベクトルデータを振幅項と位相項に分解演算処理するデータ分解部、
前記分解された複素数ベクトルデータの前記位相項を方位０度に揃える補正テーブルを演算生成する補正テーブル生成部を有し、
前記干渉成分除去手段は、
時刻ｔにおける前記第１〜第Ｋチャンネルの前記受信データxc[t]に、前記補正テーブルを作用させ、各チャンネルの前記干渉成分の位相を方位０度に揃えると共に、射影行列(I−ａ（０）．ａ（０）Ｔ)を掛けて（Ｉは単位行列(サイズＫ)、a(0)は０度方向のモードベクトル（サイズK））、該位相が方位０度に揃った各チャンネルの前記受信データxc[t]から干渉成分を除去し、更に前記干渉方向成分が除去された各チャンネルの受信データの位相を元に戻すために、前記補正テーブルの複素共役を作用させることで行う、干渉方向成分の抑圧処理を、
式（１９）を用いて、
（１）          
前記受信データxc[t]に対して、補正テーブルの共役転置Hosei[t]Hを乗じる演算を行い、
（２）          
（１）の結果に対して、１／Ｋ・補正テーブルHosei[t]を乗じる演算を行い、
（３）          
次いで、受信データxc[t]から（２）の結果を減じる演算を行って、
前記干渉方向成分の抑圧処理と等価な演算を行う干渉抑圧部、
【数１３】
0005478010.tif 000016
を有し、
更に、前記干渉方向成分が抑圧された前記受信データをマージするバッファ部を有し、
前記復元されたデータに基づいて、前記物標の距離、相対速度などを検出することを特徴とする電子走査式レーダ装置。
【請求項２】
前記ミキサと前記複数のアンテナ素子間に設けられ、前記複数のアンテナ素子を選択的に前記ミキサに接続する切替え器が設けられたことを特徴とする請求項１記載の電子走査式レーダ装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、電子走査式レーダ装置に係わり、特に、車載用のFM-CW方式またはＣＷ方式の電子走査式レーダ装置において、受信信号に含まれる干渉信号を抑圧することのできる電子走査式レーダ装置に関する。
【背景技術】
【０００２】
図１は、FM-CWレーダ方式における送受信信号と、ミキシング処理の原理を示すタイムチャート、図２は、対向車線を有する道路環境の一例を示す平面図、図３は、従来のレーダ装置において、他車からの干渉信号が受信された場合の、自車での信号処理状態示す図、図４は、全チャンネル同時受信の場合と、時分割（切替え）受信の場合の、各チャンネルのサンプリング値と、その際の想定される干渉成分信号の状態（折り返し前信号と折り返し後信号）を示す図である。
【０００３】
自動車の衝突事故防止や車間制御のために、先行する車両などの前方物標に対する距離・速度・方位を計測する車載レーダが開発されている。
【０００４】
前方物標に対する距離と相対速度を計測する手法としては、信号処理回路構成が簡易であるなどの理由からFM-CWレーダ方式が採用される。FM-CW方式では、図１（ａ）に示すように、送信アンテナより直線的に周波数が変化する信号Ｓ１を送信する。これが物標に反射してきた信号Ｓ２を受信し、受信信号Ｓ２と送信信号Ｓ１のミキシングを、図１（ｂ）に示すように、行う。これにより、送受信信号の周波数差（ビート周波数ｆｂ）を成分とするビート信号Ｓ３が生成される。このビート周波数は物標から往復伝播遅延時間Δｔに比例しており、ここから距離を換算することができる。
【０００５】
方位を計測する手法としては、短時間で全方位の走査処理が可能ものとして電子走査方式がある。電子走査方式では、対象からの反射波をある規則により配置された複数のアンテナ素子（アレーアンテナ）で受信する。この受信データのチャンネル間には、各アンテナに対する物標の方位、各アンテナの配置位置及び受信信号周波数によって決定される時間差が生じている。この時間差（または位相差）から物標の方位検出ができる。たとえばそのような手法としてディジタルビームフォーミング（DBF）が知られている。DBFでは受信データをＡＤ変換器でディジタル化した後、各チャネルとベクトルデータ（モードベクトル）との相関をとることで方位検出をおこなうことができる（例えば、非特許文献１参照）。
【０００６】
以上のように、電子走査方式では、複数のアンテナ素子での同時受信データが必要となる。しかしながら、アンテナ素子ごとにＡＤ変換器を用意する構成では、装置が複雑・高価になるため、例えば図５に示すように、各アンテナ素子６とＡＤ変換器１３の間に設けられた切替え器７により、時分割受信する構成のものが提案されている。（たとえば、特許文献１参照）。
【０００７】
このような時分割受信方式では、切替えによる遅延時間τ[k]が各チャンネルに生じる。ここでkはチャンネル番号をあらわす。もし、この切替えによる遅延時間τ[k]がビート周波数ｆｂの周期1/ｆｂに対して無視できるほど小さければ(τ[k]<<1／ｆｂ)、すべてのチャンネルを同時に受信したとみなして処理可能である。しかしながら、コスト面の都合などの理由から、比較的低速な切替え器（VCO2の駆動周波数が比較的低い形態）を採用することがある。この場合、遅延時間が無視できなくなる。各チャンネルにおける受信信号の位相の誤差が大きくなると、対象物（物標）方位の検出が不正確となるため、（１）式によって表される位相Δφ［κ］を各チャンネルで補正することが好ましい。
【０００８】
【数１】
0005478010.tif 000002
【０００９】
この位相補正処理により、時分割受信の場合でも正確な方位検出が可能となる。ところで、車載レーダ装置搭載車両が多数往来するような、例えば図２に示すような道路環境では、対向車線を走行する車両に搭載されたレーダからの電波Ｒｘ２が自車に搭載したレーダへ混入することで、自車において放射した電波Ｔｘの物標での反射波Ｒｘ１と対向車からの電波Ｒｘ２との干渉が生じる。特に他車送信アンテナからの直接波は電力レベルが比較的大きく、計測精度に及ぼす影響が大きい。
【００１０】
このような状況では、受信信号中に含まれる干渉成分を抑圧することが有効である。たとえば、特定方位からの成分を抑圧するフィルタを利用して干渉成分を抑圧する方法が提案されている。（例えば、非特許文献２参照）
【００１１】
しかしながら、前述の時分割切替え型のレーダでは、他車からの干渉信号Ｒｘ２成分の方位が適切に求められない場合がある。以下この理由を説明する。
【００１２】
図３（ａ）に示すように、他車からの干渉信号Ｒｘ２成分の変調方式がFM-CW方式やCW方式の場合、ミキシング後の信号Ｒｘ２の干渉成分は、図３（ｂ）に示すように、周波数が変動する連続信号となる。さらにＡＤ変換器によるサンプリングを実施すると、サンプリング周波数Fの1/2倍（いわゆるナイキスト周波数）以上のビート周波数は、図３（ｃ）に示すように、折り返す成分となって現れる。
【００１３】
このように周波数が時間変動しているものに対しては、複数のアンテナ素子を切替え器により切替えるチャンネル切替え方式における、各チャンネルおける受信信号の位相補正量の決定および位相補正処理が困難となる。もしFM−CW変調の傾きが平行に近い条件であれば、干渉成分の周波数変動が比較的緩やかとなるが、その周波数は折り返している可能性があり、位相補量を一意に決定することが困難である。図４に、サンプリング後のデータが同一であっても折返し前周波数によって位相のズレ量が異なる例を示す。図４（ｂ）の切替え受信の場合（切替え遅延時間τ）、点線で示すラインが折り返し前の信号ＢＳであり、実線で示すラインが折り返し後の信号ＡＳとなる。各信号ＢＳ、ＡＳの周期が異なるために、サンプリングの値に対する必要な位相補正量が、各信号ＢＳ、ＡＳで大きく異なることが分かる。
【００１４】
そこで、サンプリングされたビート信号を、各アンテナ素子について時間方向に複数の短時間データに切り出し、その短時間データの周波数スペクトルから、干渉波の干渉成分周波数を検出した後、干渉波の折り返し前周波数の複数の候補に基づいてを位相補正を行ない、更に、補正後の信号のデジタルビームフォーミング結果を用いて、干渉成分の一番尤もらしい到来方位を推定し、推定された干渉成分の到来方位から、当該干渉成分を抑圧するフィルタを作用させ、干渉成分を抑圧する技術が提案されている（例えば、特許文献２参照）
【００１５】
【非特許文献１】菊間信良著「アレーアンテナによる適応信号処理」，科学技術出版，1998 年）
【非特許文献２】論文：Adaptive Mainbeam Jamming Suppressionfor Multi-Function Radars, T.J. Nohara 他著」）
【特許文献１】特開平１１−２３１０４０号公報
【特許文献２】特開２００７−２３２３８３号公報
【発明の開示】
【発明が解決しようとする課題】
【００１６】
しかし、この場合、干渉成分の到来方位を検出する場合には、干渉成分の折り返し前周波数の候補の数だけ、位相補正およびDBF(デジタルビームフォーミング)による方位推定を繰り返す必要がある。特に切替受信でカット周波数を高くとる必要がある(折り返してくる周波数帯が広くなる)場合、折り返し前周波数の候補の数も多くなり、演算回数が増大し、リアルタイム処理が必要なレーダ装置において短時間での処理を可能とするためには、高性能、従って複雑かつ高価な演算手段（演算プログラム及び演算素子）が必要となる不都合がある。
【００１７】
そこで、本発明は、前述の問題点を解決するため、FM−CW方式又はＣＷ方式の電子走査型レーダにおいて、簡易な構成で干渉抑圧が可能な、電子走査式レーダ装置を提供することを目的とする。
【課題を解決するための手段】
【００１８】
請求項１の発明は、連続波からなる送信信号（Ｔｘ）を、放射自在な送信アンテナ（５）、第１チャンネルから第Ｋチャンネルまでの複数のアンテナ素子（６）からなる受信アンテナ、前記複数のアンテナ素子で受信される受信信号（Ｒｘ）と前記送信信号をミキシングして前記複数のアンテナ素子（６）に対応した複数チャンネル分（例えば、Ｋチャンネル分）のビート信号（Ｓ３）を得るミキサ（１０）、前記ミキサで得られたビート信号（Ｓ３）を所定のサンプリング周波数（ｆＳ）でサンプリングして前記複数のアンテナ素子に対応した複数チャンネル分の受信データ（ＤＴ１）を得る、Ａ／Ｄ変換器（１３）、前記サンプリングされた前記複数チャンネル分の受信データ（ＤＴ１）を、各チャンネルについて時間方向に複数の短時間データ（ＳＤ）に切り出す、短時間データ切出し部（１９）、前記各チャンネルについて、前記複数の短時間データ（ＳＤ）の周波数スペクトルを算出する、周波数スペクトル算出部（２０）、前記周波数スペクトルから、干渉波の干渉成分周波数を検出する、干渉周波数検出部（１８）、検出された各チャンネルの干渉成分周波数における周波数スペクトル成分に基づいて、位相補正テーブルを算出する位相補正テーブル算出部（２１）、前記位相補正テーブルに基づいて前記受信データ中の干渉方向成分を抑圧する干渉成分除去手段（２６）、該干渉方向成分が抑圧された前記受信データ（ＤＴ１）に基づいて物標の距離、相対速度などを検出する、物標検出部（１７）、を有する、電子走査式レーダ装置（１）において、
前記位相補正テーブル算出部は、
前記干渉周波数検出部で得られる前記干渉成分周波数における周波数スペクトル成分で構成される複素数ベクトルデータを振幅項と位相項に分解演算処理するデータ分解部、
前記分解された複素数ベクトルデータの前記位相項を方位０度に揃える補正テーブルを演算生成する補正テーブル生成部を有し、
前記干渉成分除去手段（２６）は、
時刻ｔにおける前記第１〜第Ｋチャンネルの前記受信データxc[t]に、前記補正テーブルを作用させ、各チャンネルの前記干渉成分の位相を方位０度に揃えると共に、射影行列(I−ａ（０）．ａ（０）Ｔ)を掛けて（Ｉは単位行列(サイズＫ)、a(0)は０度方向のモードベクトル（サイズK））、該位相が方位０度に揃った各チャンネルの前記受信データxc[t]から干渉成分を除去し、更に前記干渉方向成分が除去された各チャンネルの受信データの位相を元に戻すために、前記補正テーブルの複素共役を作用させることで行う、干渉方向成分の抑圧処理を、
式（１９）を用いて、
（１）          
受信データxc[t]に対して、補正テーブルの共役転置Hosei[t]Hを乗じる演算を行い、
（２）          
（１）の結果に対して、１／Ｋ・補正テーブルHosei[t]を乗じる演算を行い、
（３）          
次いで、受信データxc[t]から（２）の結果を減じる演算行って、
前記干渉方向成分の抑圧処理と等価な演算を行う干渉抑圧部、
【数１３】
0005478010.tif 000003
を有し、
更に、前記干渉方向成分が抑圧された前記受信データをマージするバッファ部（２７）を有し、
前記復元されたデータに基づいて、前記物標の距離、相対速度などを検出することを特徴として構成される。
【００１９】
請求項２の発明は、前記ミキサ（１０）と前記複数のアンテナ素子（６）間に設けられ、前記複数のアンテナ素子（６）を選択的に前記ミキサに接続する切替え器（７）が設けられて構成される。
【発明の効果】
【００２１】
請求項１の発明によれば、対向車などからのＦＭ−ＣＷ、ＣＷレーダ波などを干渉波として受信したような場合などにおいて、受信データの干渉成分の位相を方位０度に揃えた後、当該受信データから干渉成分のみを除去するようにしたので、従来のように、干渉成分の折り返し前周波数の候補の数だけのデジタルビームフォーミング処理を行うといった、演算量の大きな処理を行う必要が無くなり、簡易な構成で干渉抑圧が可能な、電子走査式レーダ装置を提供することができる。
【００２２】
請求項２の発明によれば、切替え遅延時間が問題となる複数のアンテナ素子（６）を切替え器（７）で選択的にミキサ（１０）に接続する、チャンネル切替え形の電子走査式レーダ装置（１）に適用することが可能となる。
【００２３】
また、請求項１の発明によれば、複素ベクトルの位相項を方位０度に揃えるように補正テーブルを演算生成すればよいので、補正テーブルの算出のためのコストが少なくてすむ。
【００２４】
なお、括弧内の番号等は、図面における対応する要素を示す便宜的なものであり、従って、本記述は図面上の記載に限定拘束されるものではない。
【発明を実施するための最良の形態】
【００２５】
以下、図面に基づき、本発明の実施例を説明する。
【００２６】
図５は、本発明による電子走査式レーダ装置の１実施例を示すブロック図、図６は、短時間データ切り出し処理の内容を示す模式図、図７は、周波数スペクトル算出処理の内容を示す模式図、図８は、各時刻における干渉信号の瞬時ビート周波数を示す模式図、図９は、本発明の基本思想を示す概念図、図１０は、本発明による干渉信号除去の流れを示すフローチャート（一例）である。
【００２７】
図５は、本発明の一実施形態である電子走査式レーダ装置１を示す構成図である。このレーダ装置１は、連続波（ＣＷ）に周波数変調（ＦＭ）を掛けた送信信号Ｔｘを用いるＦＭ−ＣＷレーダ装置である。このレーダ装置１は、自動車に搭載されるいわゆる車載用レーダ装置であり、前方を走行する車輌（物標）までの距離やその相対速度などを検知するものである。このレーダ装置１の検知結果は、車輌走行の制御情報等に利用される。送信電波にはマイクロ波が用いられている。
【００２８】
このレーダ装置１では、切換え器７を利用することにより、ＲＦアンプ９やミキサ１０などのアナログデバイスを全体で一組だけ備えた構成となっている。レーダ装置１は、送受信部４を有しており、送受信部４は、中心周波数がｆ０（たとえば７６ＧＨｚ）の発振器２と、アンプ３と、送信アンテナ５とを備えている。発振器２は、図示しない変調用の直流電源から出力される制御電圧によって、周波数ｆ０の搬送波に対して周波数変調幅ΔＦの三角波変調を掛けた信号、すなわち周波数ｆ０±ΔＦ／２の被変調波（送信信号Ｔｘ）を出力する。被変調波はアンプ３で増幅され、送信アンテナ５から電磁波として放射される。なお、送信信号Ｔｘの一部は受信検波用のローカル信号としてミキサ１０に出力される。
【００２９】
送受信部４に設けられた受信用アレーアンテナ８は、第１チャネル（♯１）から第Ｋチャネル（♯Ｋ）までの各チャネルに対応するＫ個のアレーアンテナ素子６を備えている。切替え器７は、Ｋ個の入力端子と１個の出力端子とを有し、各入力端子にはアレーアンテナ８の各アレーアンテナ素子６が１個づつ接続されている。出力端子は入力端子のいずれか一つと接続されるものであり、切換信号（クロック信号）により、その接続は周期的に切替えられる。接続切替えは、回路上で電気的に行われる。
【００３０】
受信信号Ｒｘは切替器７で周期１／ｆｓｗで時分割多重化される。ここで、切替えの順番はランダムに行うものとする。時分割多重化された信号は、ＲＦアンプ９で増幅され、ミキサ１０により分配された送信信号Ｔｘとミキシングされる。このミキシングにより受信信号Ｒｘはダウンコンバートされ、図１（ｂ）に示すように、送信信号Ｔｘと受信信号Ｒｘとの差信号であるビート信号Ｓ３が生成される。受信信号ＲＸ及び送信信号ＴＸに基づいてビート信号Ｓ３を得る処理の詳細は、例えば特開平１１−１３３１４２号公報などで述べられている公知技術なので、本明細書ではその詳細な説明は省略する。
【００３１】
ところで、三角波変調ＦＭ−ＣＷ方式では、相対速度が零のときのビート周波数をｆｒ、相対速度に基づくドップラ周波数をｆｄ、周波数が増加する区間（アップ区間）のビート周波数をｆｂ１、周波数が減少する区間（ダウン区間）のビート周波数をｆｂ２とすると、
ｆｂ１＝ｆｒ−ｆｄ …（２）
ｆｂ２＝ｆｒ＋ｆｄ …（３）
が成り立つ。
【００３２】
従って、変調サイクルのアップ区間とダウン区間のビート周波数ｆｂ１およびｆｂ２を別々に測定すれば、次式（４）及び（５）からｆｒおよびｆｄを求めることができる。
ｆｒ＝（ｆｂ１＋ｆｂ２）／２ …（４）
ｆｄ＝（ｆｂ２−ｆｂ１）／２ …（５）
ｆｒおよびｆｄが求まれば、目標物の距離Ｒと速度Ｖを次の（６）（７）式により求めることができる。
Ｒ＝（Ｃ／（４・ΔＦ・ｆｍ））・ｆｒ …（６）
Ｖ＝（Ｃ／（２・ｆ０））・ｆｄ …（７）
ここに、Ｃは光の速度、ｆｍはＦＭ変調周波数である。
【００３３】
生成されたビート信号Ｓ３は、アンプ１１、ローパスフィルタ１２を経由して、A／D変換器１３にてサンプリング周波数ｆｓでN個のデータとしてサンプリング量子化される。そして式（８）のようなK（チャンネル）×N個の受信データＤＴ１としてバッファ部１４へ蓄積され、物標検出部１７に出力される。
【００３４】
【数２】
0005478010.tif 000004
【００３５】
物標検出部１７は、図５に示すように、干渉抑圧部３０，ビート周波数検出部３１，位相補正部３２及び方位検出部３３を有しており、干渉抑圧部３０は、短時間データ切り出し部１９，周波数スペクトル算出部２０，干渉周波数検出部１８，位相補正テーブル算出部２１、干渉成分除去部２６及びバッファ部２７を有している。
【００３６】
短時間データ切り出し部１９では、図６に示すように、各アレーアンテナ素子６に対応する各チャンネルについて、それぞれ時間方向にN個蓄積された受信データＲＤを、次式のようなM個づつの短いデータＳＤに切り出す。なお、データＳＤを切り出す際の、各データＳＤのずらし数は、式（９）の場合は、１であるが、ずらし数は、１以上の適宜な数に可変とすることが出来る。
【００３７】
【数３】
0005478010.tif 000005
【００３８】
次に、周波数スペクトル算出部２０では、短時間に切り出したデータに対し、図７及び式（１０）に示すように、離散フーリエ変換を行い周波数領域へのデータに変換して、周波数スペクトルを演算算出する。
【００３９】
【数４】
0005478010.tif 000006
【００４０】
干渉周波数検出部１８では、Kチャンネル分の離散フーリエ変換後電力の平均を求め、図８に示すように、さらにその周波数方向のピークを検出しそのピークの平均電力レベルが最大となる周波数をもって、各時刻tでの干渉成分の瞬時ビート周波数（干渉成分周波数）として求める。これを
【００４１】
【数５】
0005478010.tif 000007
と表記する。
【００４２】
【数６】
0005478010.tif 000008
【００４３】
ここで、求められた干渉成分の瞬時ビート周波数ｆBA[t]（干渉成分周波数）における1からK chの複素周波数スペクトルは、式（１３）で表される。
【００４４】
【数７】
0005478010.tif 000009
【００４５】
図９に示すように、周波数スペクトル算出部２０で算出されたビート周波数スペクトルから、干渉周波数検出部１８で得られる干渉成分周波数fBAの成分を取り出したもの（式（１３））は、ほぼ干渉成分だけで構成される。このことから、本発明では、それらの信号から、従来のように、折り返し前周波数の複数の候補を生成して、それぞれの候補に対応した位相補正量を求め、そのかなかで一番尤もらしい干渉成分の方位を求める大規模な演算処理を行わなくても、各チャンネルの受信信号中の干渉成分のみを補正（位相成分を揃える）してからその主成分（直流成分）を除去する処理を行うことで、簡単に干渉成分を除去することが出来るであろうとの知見に基づいて処理を行う。
【００４６】
即ち、本発明による処理を模式的に示すと、図１０に示すように、干渉成分が混入した受信信号は、ステップＳ１に示すように、各受信チャンネルＣｈ１，Ｃｈ２，Ｃｈ３……について、干渉成分ＫＣと本来の物標からの反射波からなる所望成分ＳＣがその位相及び大きさが異なる形で混在している。なお、図１０において、各チャンネルにおける干渉成分ＫＣ及び所望成分ＳＣは、その矢印の長さが信号の大きさを表し、矢印の向きが、位相を示している。なお、図１０の表示は、本発明の概念を分かり易く表示するための模式図であり、各信号の大きさや方位は、必ずしも実際の信号状態を反映するものではない。
【００４７】
次に、ステップＳ２で、各チャンネルにおける受信信号の干渉成分ＫＣの位相をそろえ、次いでステップＳ３で当該干渉成分ＫＣを除去する処理を行うことで、簡単に受信信号中の干渉成分ＫＣを除去することができる。その後、ステップＳ４で、各チャンネルの位相を基に戻すと、干渉成分ＫＣが除去された受信信号を得ることが出来る。
【００４８】
以下、干渉抑圧処理の詳細について記述する。まず、式（１３）で示される干渉成分周波数における周波数スペクトル成分で構成される複素数ベクトルデータＹ［］［ｔ］（ｆＢＡ（ｔ））を、振幅項（ａｉ）と位相項（ｅｊθｉ）に分解演算処理する。
【００４９】
【数８】
0005478010.tif 000010
このとき、位相補正テーブル算出部２１は、振幅項（ａｉ）と位相項（ｅｊθｉ）に分解された受信データの位相項を１にする補正テーブルを演算生成する。即ち、時刻ｔにおける補正テーブルHosei[t]、即ち、図１０におけるステップＳ２において、各チャンネルにおける受信信号の干渉成分（信号）の位相をそろえるための式（１５）は以下のようなものとなる。
【００５０】
【数９】
0005478010.tif 000011
この補正テーブルは、式（１４）で示した受信データの位相項（ｅｊθｉ）を１にする、即ち、受信データにおける干渉信号の位相を揃える働きをする（図１０のステップＳ２に相当）。即ち、各チャンネルの受信データの位相を、各チャンネルについて、干渉成分（信号）の位相を基準に、それら干渉成分（信号）の位相が揃うように、各チャンネルの受信データを演算処理する。このように、従来は折り返し周波数の候補の数だけ、位相補正およびDBFを繰り返して干渉方位を検出後に補正テーブルを算出する必要があったが、提案方式は、干渉方位を求める必要が無くなり、１回の処理で補正テーブルを求めることができる。
【００５１】
こうして、各受信信号における干渉信号の位相をそろえるための補正テーブルHosei[t]が求められたところで、時刻tにおける１〜Ｋチャンネルの受信データをXc[t]とすると、式（１６）のようになる。
【００５２】
【数１０】
0005478010.tif 000012
【００５３】
干渉成分除去部２６では、以下の式（１７）によって干渉抑圧を実施する。
【００５４】
【数１１】
0005478010.tif 000013
Ｉは単位行列(サイズＫ)、a(0)は０度方向のモードベクトル（サイズK）である（式（１８））。式（１７）で、＊は複素共役を表す。
【００５５】
【数１２】
0005478010.tif 000014
干渉成分除去部２６で、式（１６）に示す受信データXc[t]に補正テーブルを作用させることで、Xc[t]中の干渉成分の位相が揃う(位相項が、例えば全て１、従って方位０度（位相が揃えば、必ずしも方位０度に限定されるわけではない）になる、図１０のステップＳ２、)こととなる。次に受信信号中の直流成分(方位０度に相当)の信号を除去する射影行列(I−ａ（０）．ａ（０）Ｔ)を掛けると、受信信号中の干渉成分が除去される（図１０のステップＳ３）こととなる。最後に図１０のステップ２で各チャンネルの受信データの位相補正した分を、ステップＳ４で元に戻す（図１０のステップＳ４）。これにより、受信データから干渉成分のみが除去されることとなる。
【００５６】
式（１７）は以下の式（１９）と等価であるので、実際はこちらを計算する。その際、演算コスト削減のために、後ろから順番に計算する。式（１９）で、Ｈは共役転置を表す。
【００５７】
【数１３】
0005478010.tif 000015
バッファ部２７では、干渉信号成分が抑圧されたデータｘC[t]を元のデータ数分蓄積して、後段のビート周波数検出部３１へ送る。この状態で、後段のビート周波数検出部３１へは、図５の送受信部５のバッファ部１４に蓄積されたビート信号から、干渉波成分が除去（抑圧）された形の信号が適切に出力される。
【００５８】
物標検出部１７の干渉抑圧部３０で干渉成分が抑圧されたビート信号は、ビート周波数検出部３１、位相補正部３２及び方位検出部３３で公知の処理が施され、自車と先行車両などの物標との距離、相対速度、方位などが演算され、更に、図５に示す、物標追従処理部３５において、時間的な追跡処理を行って前方の車両を検出するなどの演算処理を行う。なお、物標追従処理部３５における詳しい処理内容については、特開２００３−２７０３４１号公報などにその詳細が述べられている公知技術なので、本明細書でははその説明を省略する。また、ビート周波数検出部３１、位相補正部３２及び方位検出部３３での処理は、非特許文献１等に詳細に述べられており、公知の手法となっているので、ここでは、その詳細な説明は省略する。
【００５９】
なお、本実施例では、これらの処理部とその動作内容をマイクロプロセッサやディジタルシグナルプロセッサー等で動作する信号処理ソフトウエアとしての実現を想定して説明を行うが、FPGAやLSI等の半導体デバイス上の集積回路としての実現も可能である。
【００６０】
なお、上述の実施例は、切替え器７で複数のアレーアンテナ素子６を切り替えて単一のＡ／Ｄ変換器１３でビート信号Ｓ３を量子化する構成としたが、本発明は、切替器７を用いずに、各アレーアンテナ素子についてそれぞれＡ／Ｄ変換器１３を接続して、各チャンネルの同時受信をする構成であっても、受信信号に混入した干渉成分を除去することにおいて、同様に適用することが出来る。
【００６１】
このように、対向車などからのＦＭ−ＣＷ、ＣＷレーダ波などを干渉波として受信したような場合などにおいて、干渉成分周波数の時間変動があっても、当該干渉波の短時間データを切り出して処理する構成としたため、切り出した時間範囲では周波数がほとんど変わらず、周波数スペクトル算出部２０と干渉周波数検出部１８によって、その時間区間での干渉成分周波数を検出することができる。
【産業上の利用可能性】
【００６２】
本発明は、車載用のFM-CW方式またはＣＷ方式の電子走査式レーダ装置に利用することが出来る。
【図面の簡単な説明】
【００６３】
【図１】図１は、FM-CWレーダ方式における送受信信号と、ミキシング処理の原理を示すタイムチャート。
【図２】図２は、対向車線を有する道路環境の一例を示す平面図。
【図３】図３は、従来のレーダ装置において、他車からの干渉信号か受信された場合の、自車での信号処理状態示す図。
【図４】図４は、全チャンネル同時受信の場合と、時分割（切替え）受信の場合の、各チャンネルのサンプリング値と、その際の想定される干渉成分信号の状態（折り返し前信号と折り返し後信号）を示す図。
【図５】図５は、本発明による電子走査式レーダ装置の１実施例を示すブロック図。
【図６】図６は、短時間データ切り出し処理の内容を示す模式図。
【図７】図７は、周波数スペクトル算出処理の内容を示す模式図。
【図８】図８は、各時刻における干渉信号の瞬時ビート周波数を示す模式図。
【図９】図９は、本発明の基本思想を示す概念図。
【図１０】図１０は、本発明による干渉信号除去の流れを示すフローチャート（一例）。
【符号の説明】
【００６４】
１……電子走査式レーダ装置
５……送信アンテナ
６……アンテナ素子
７……切替えスイッチ
１０……ミキサ
１３……Ａ／Ｄ変換器
１７……物標検出部
１８……干渉周波数検出部
１９……短時間データ切出し部
２０……周波数スペクトル算出部
２１……位相補正テーブル算出部
２６……干渉成分除去手段（干渉成分除去部）
２７……バッファ部
Ｓ３、ＲＤ……ビート信号
ＲＸ……受信信号
ＴＸ……送信信号
【図１】
0005478010.tif 000017
【図２】
0005478010.tif 000018
【図３】
0005478010.tif 000019
【図４】
0005478010.tif 000020
【図５】
0005478010.tif 000021
【図６】
0005478010.tif 000022
【図７】
0005478010.tif 000023
【図８】
0005478010.tif 000024
【図９】
0005478010.tif 000025
【図１０】
0005478010.tif 000026