(57)【特許請求の範囲】
【請求項１】
電波によって物体を検知するための物体検知装置であって、
前記物体に向けて、送信信号となる電波を照射する、送信部と、
複数の受信アンテナによって、前記物体で反射された電波を受信信号として受信し、更に、複数の前記受信アンテナそれぞれで受信された前記受信信号毎に、当該受信信号を用いて、中間周波数信号を生成する、受信部と、
前記中間周波数信号それぞれを合成して得られるビームパターンによる虚像の発生が抑制されるように、サンプリング時刻を決定し、決定したサンプリング時刻で前記中間周波数信号に対してサンプリングを行うことによって、物体の位置検出用の中間周波数信号を生成し、
生成した前記位置検出用の中間周波数信号を用いて、前記物体を検知する、演算部と、
を備え、
前記送信部は、前記物体に向けて周波数が時間に応じて変化する電波を、前記送信信号として照射し、
前記受信部は、前記受信アンテナの各々で受信した前記受信信号に前記送信信号をミキシングして前記中間周波数信号を生成し、
前記受信アンテナが、複数備えられ、
隣接する前記受信アンテナの間隔が、前記送信信号の周波数の最小値と最大値とで定められる上限値以下に設定されている、
ことを特徴とする物体検知装置。
【請求項２】
請求項１に記載の物体検知装置であって、
前記受信部が、基準受信アンテナと測定受信アンテナとを備え、
前記演算部は、
前記測定受信アンテナで受信された受信信号から生成された中間周波数信号と、前記基準受信アンテナで受信された受信信号から生成された中間周波数信号とに対して、サンプリングを実行し、
得られたサンプリング値について、前記測定受信アンテナで受信された受信信号から生成された中間周波数信号を、前記基準受信アンテナで受信された受信信号から生成された中間周波数信号によって、規格化して、前記位置検出用の中間周波数信号として、規格化中間周波数信号を生成し、
更に前記規格化中間周波数信号から、前記物体の位置分布を表す評価関数を算出し、算出した前記評価関数を用いて、前記物体の位置及び形状を検知する、
ことを特徴とする物体検知装置。
【請求項３】
請求項２に記載の物体検知装置であって、
前記演算部は、前記測定受信アンテナで受信された受信信号から生成された中間周波数信号と、前記基準受信アンテナで受信された受信信号から生成された中間周波数信号とに対し、検知される前記物体の位置を制限する処理を実行し、
前記処理が実行された前記中間周波数信号を用いて、前記サンプリングを実行する、
ことを特徴とする物体検知装置。
【請求項４】
請求項２または３に記載の物体検知装置であって、
前記演算部は、前記規格化中間周波数信号から相関行列を算出し、算出した前記相関行列から、前記評価関数を算出する、
ことを特徴とする物体検知装置。
【請求項５】
請求項４に記載の物体検知装置であって、
前記測定受信アンテナが、複数備えられ、且つ、前記基準受信アンテナの位置を基準にして、少なくとも二つの方向それぞれに沿って配置されており、
前記演算部は、前記方向毎に、当該方向に設置された前記測定受信アンテナから、前記規格化中間周波数信号を生成し、更に、前記方向毎に生成された前記規格化中間周波数信号の直積から、前記相関行列を算出し、そして、算出した前記相関行列から、前記物体の多次元の位置分布を表す前記評価関数を算出する、
ことを特徴とする物体検知装置。
【請求項６】
請求項４または５に記載の物体検知装置であって、
前記演算部は、前記規格化中間周波数信号の範囲を複数選択し、選択した前記複数の範囲を定めた複数の前記規格化中間周波数信号それぞれから、前記相関行列を算出し、更に前記相関行列の平均から、前記物体の位置分布を表す前記評価関数を算出する、
ことを特徴とする物体検知装置。
【請求項７】
電波によって物体を検知するための方法であって、
（ａ）前記物体に向けて、送信信号となる電波を照射する、ステップと、
（ｂ）複数の受信アンテナによって、前記物体で反射された電波を受信信号として受信し、更に、複数の前記受信アンテナそれぞれで受信された前記受信信号毎に、当該受信信号を用いて、中間周波数信号を生成する、ステップと、
（ｃ）前記中間周波数信号それぞれを合成して得られるビームパターンによる虚像の発生が抑制されるように、サンプリング時刻を決定し、決定したサンプリング時刻で前記中間周波数信号に対してサンプリングを行うことによって、物体の位置検出用の中間周波数信号を生成する、ステップと、
（ｄ）生成した前記位置検出用の中間周波数信号を用いて、前記物体を検知する、ステップと、を有し、
前記（ａ）のステップにおいて、前記物体に向けて周波数が時間に応じて変化する電波を、前記送信信号として照射し、
前記（ｂ）のステップにおいて、前記受信アンテナの各々で受信した前記受信信号に前記送信信号をミキシングして前記中間周波数信号を生成し、
前記受信アンテナが、複数備えられ、
前記（ｂ）のステップにおいて、隣接する前記受信アンテナの間隔が、前記送信信号の周波数の最小値と最大値とで定められる上限値以下に設定されている、
ことを特徴とする物体検知方法。
【請求項８】
請求項７に記載の物体検知方法であって、
前記（ｂ）のステップにおいて、前記受信アンテナとして、基準受信アンテナと測定受信アンテナとを用い、
前記（ｃ）のステップにおいて、
前記測定受信アンテナで受信された受信信号から生成された中間周波数信号と、前記基準受信アンテナで受信された受信信号から生成された中間周波数信号とに対して、サンプリングを実行し、
得られたサンプリング値について、前記測定受信アンテナで受信された受信信号から生成された中間周波数信号を、前記基準受信アンテナで受信された受信信号から生成された中間周波数信号によって、規格化して、前記位置検出用の中間周波数信号として、規格化中間周波数信号を生成し、
前記（ｄ）のステップにおいて、前記規格化中間周波数信号から、前記物体の位置分布を表す評価関数を算出し、算出した前記評価関数を用いて、前記物体の位置及び形状を検知する、
ことを特徴とする物体検知方法。
【請求項９】
請求項８に記載の物体検知方法であって、
前記（ｃ）のステップにおいて、
前記測定受信アンテナで受信された受信信号から生成された中間周波数信号と、前記基準受信アンテナで受信された受信信号から生成された中間周波数信号とに対し、検知される前記物体の位置を制限する処理を実行し、
前記処理が実行された前記中間周波数信号を用いて、前記サンプリングを実行する、
ことを特徴とする物体検知方法。
【請求項１０】
請求項８または９に記載の物体検知方法であって、
前記（ｄ）のステップにおいて、前記規格化中間周波数信号から相関行列を算出し、算出した前記相関行列から、前記評価関数を算出する、
ことを特徴とする物体検知方法。
【請求項１１】
請求項１０に記載の物体検知方法であって、
前記測定受信アンテナが、複数備えられ、且つ、前記基準受信アンテナの位置を基準にして、少なくとも二つの方向それぞれに沿って配置されており、
前記（ｃ）のステップにおいて、前記方向毎に、当該方向に設置された前記測定受信アンテナから、前記規格化中間周波数信号を生成し、
前記（ｄ）のステップにおいて、前記方向毎に生成された前記規格化中間周波数信号の直積から、前記相関行列を算出し、そして、算出した前記相関行列から、前記物体の多次元の位置分布を表す前記評価関数を算出する、
ことを特徴とする物体検知方法。
【請求項１２】
請求項１０または１１に記載の物体検知方法であって、
前記（ｄ）のステップにおいて、前記規格化中間周波数信号の範囲を複数選択し、選択した前記複数の範囲を定めた複数の前記規格化中間周波数信号それぞれから、前記相関行列を算出し、更に前記相関行列の平均から、前記物体の位置分布を表す前記評価関数を算出する、
ことを特徴とする物体検知方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、電波を検知対象物に照射し、対象物からの反射ないし放射された電波を検知する事で検知対象物の存在を認識又は識別するための、物体検知装置、及び物体検知方法に関し、更には、これらを実現するためのプログラムに関する。
【０００２】
電波（マイクロ波、ミリ波、テラヘルツ波など）は、光に比べて、物体を透過する能力の点で優れている。電波の透過能力を活用し、衣服下、鞄内の物品を画像化して検査を行う装置が実用化されている。また、電波の透過能力を活用して、衛星又は航空機から雲を透過して地表を画像化するリモートセンシング技術も実用化されている。【０００３】
電波を用いたイメージング装置（物体検知装置）としては、いくつかの方式が提案されている。一つは、図１２で概念図を示したアレイアンテナ方式である。図１２は、一般的なアレイアンテナ方式の概念を示す概念図である。図１２に示すように、アレイアンテナ方式において、測定装置は、送信機２１１と、受信機２０１とを備えている。【０００４】
また、図１２に示すように、送信機２１１は、送信アンテナ２１２を備えている。また、受信機２０１は、受信アンテナ２０２１、２０２２、・・・、及び２０２Ｎを備えている（Ｎは受信アンテナの数）。【０００５】
送信機２１１は、送信アンテナ２１２からＲＦ信号（電波）２１３を検知対象物２０４１、２０４２、・・・、２０４Ｋ（Ｋは対象物の数）に向けて照射する。ＲＦ信号（電波）２１３は検知対象物２０４１、２０４２、・・・、２０４Ｋにおいて反射され、反射波２０３１、２０３２、・・・、２０３Ｋがそれぞれ発生する。発生した反射波２０３１、２０３２、・・・、２０３Ｋは、受信アンテナ２０２１、２０２２、・・・、２０２Ｎにおいて受信される。【０００６】
受信機２０１は、受信した反射波２０３１、２０３２、・・・、２０３Ｋに基づいて、検知対象物２０４１、２０４２、・・・、２０４Ｋから反射されている電波強度を算出する。その電波強度の分布を画像化する事で、検知対象物２０４１、２０４２、・・・、２０４Ｋの像を得る事ができる。【０００７】
アレイアンテナ方式では、図１３で示すように、受信機２０１における受信アンテナ２０２１、２０２２、・・・、２０２Ｎそれぞれ毎に、移相器２０６１、２０６２、・・・、２０６Ｎがそれぞれ備えられている。図１３は、図１２に示した受信機の構成をより具体的に示す図である。図１３に示すように、移相器２０６１、２０６２、・・・、２０６Ｎは、受信アンテナ２０２１、２０２２、・・・、２０２Ｎで受信した到来波２０８１、２０８２、・・・、２０８Ｎに対し、位相回転Φ１、Φ２、・・・、ΦＮをそれぞれ加える。位相回転Φ１、Φ２、・・・、ΦＮが加えられた到来波２０８１、２０８２、・・・、２０８Ｎは、加算器２０７で加算される。【０００８】
移相器２０６１、２０６２、・・・、２０６Ｎと、加算器２０７とは、アナログ回路で実装される事もあれば、ソフトウェアによるデジタル処理で実装される事もある。アレイアンテナ方式では、移相器２０６１、２０６２、・・・、２０６Ｎによる位相回転Φ１、Φ２、・・・、ΦＮの設定により、アレイアンテナの指向性が制御される。受信アンテナ２０２の指向性をｇ（θ）、受信アンテナ２０２ｎで受信した到来波２０８ｎ（ｎ＝１，２，・・・，Ｎ）の振幅と位相をそれぞれａｎおよびφｎとした場合、アレイアンテナの指向性Ｅ（θ）は、以下の式（１）のように計算される。【０００９】
【数１】
0006939981.tif 000002
【００１０】
なお、式（１）において、アレイアンテナの指向性Ｅ（θ）から受信アンテナ２０２の指向性ｇ（θ）を除去した指向性成分ＡＦ（θ）は、アレイファクターと呼ばれる。アレイファクターＡＦ（θ）が、アレイアンテナを形成した事による指向性の効果を表す。受信アンテナ２０２ｎ（ｎ＝１，２，・・・，Ｎ）で受信した信号はｇ（θ）ａｎｅｘｐ（ｊφｎ）であり、移相器２０６ｎの位相回転Φｎを受けた信号ｇ（θ）ａｎｅｘｐ（ｊφｎ）ｅｘｐ（ｊΦｎ）が、ｎ＝１，２，・・・，Ｎに渡って加算器２０７で加算されて得られた信号が式（１）の指向性Ｅ（θ）として得られる。【００１１】
受信アンテナ２０２ｎ（ｎ＝１，２，・・・，Ｎ）で入射角θの到来波２０８ｎを受信する際、到来波２０８ｎの位相φｎは、−２π・ｎ・ｄ・ｓｉｎθ／λで与えられる（ｎ＝１，２，・・・，Ｎ）。ここで、ｄは受信アンテナ２０２ｎ（ｎ＝１，２，・・・，Ｎ）の間隔であり、λは到来波２０８１、２０８２、・・・、２０８Ｎの波長である。【００１２】
式（１）において、振幅ａｎがｎによらず一定とした場合、移相器２０６ｎの位相回転Φｎ（ｎ＝１，２，・・・，Ｎ）を到来波２０８ｎの位相φｎと等しくなるように設定すると、アレイファクターＡＦ（θ）は角度θの方向において最大となる。この事がすなわち、移相器２０６ｎの位相回転Φｎによるアレイアンテナの指向性の制御法を示している。【００１３】
アレイアンテナ方式による電波イメージング装置の例は、特許文献１〜３において開示されている。【００１４】
特許文献１及び２に開示されたアレイアンテナ方式では、移相器２０６１、２０６２、・・・、２０６Ｎは、受信機２０１に内蔵され、且つ、受信アンテナ２０２１、２０２２、・・・、２０２Ｎに接続されている。これら移相器２０６１、２０６２、・・・、２０６Ｎにより、受信アンテナ２０２１、２０２２、・・・、２０２Ｎで形成される受信アレイアンテナの指向性が制御される。【００１５】
つまり、ビーム状に形成された受信アレイアンテナ（２０２１、２０２２、・・・、２０２Ｎ）の指向性を変化させ、検知対象物２０４１、２０４２、・・・、２０４Ｋのそれぞれに受信アレイアンテナ（２０２１、２０２２、・・・、２０２Ｎ）の指向性ビームを向ける事で、検知対象物２０４１、２０４２、・・・、２０４Ｋから反射されている電波強度が算出される。【００１６】
特許文献３に開示されたアレイアンテナ方式では、受信アレイアンテナ（２０２１、２０２２、・・・、２０２Ｎ）の周波数依存性が利用されて、受信アレイアンテナ（２０２１、２０２２、・・・、２０２Ｎ）の指向性が制御されている。特許文献３に開示されたアレイアンテナ方式でも、検知対象物２０４１、２０４２、・・・、２０４Ｋのそれぞれに受信アレイアンテナ（２０２１、２０２２、・・・、２０２Ｎ）の指向性ビームを向ける事で、検知対象物２０４１、２０４２、・・・、２０４Ｋから反射されている電波強度が算出される。この点で、特許文献３に開示されたアレイアンテナ方式も、特許文献１又は２に開示されたアレイアンテナ方式と共通である。【００１７】
ところで、アレイアンテナ方式では、対象物２０４の虚像を抑制するため、受信アンテナ２０２１、２０２２、・・・、２０２Ｎの各アンテナの間隔は、受信機２０１において受信される反射波２０３１、２０３２、・・・、２０３Ｋの波長λの半分以下で無ければならないという制約がある。【００１８】
アレイアンテナ方式における受信アンテナ２０２１、２０２２、・・・、２０２Ｎの間隔の制約を緩和する技術が特許文献４において開示されている。特許文献４では、ＲＦ周波数が異なる複数のパルス信号が対象物に向けて送信され、対象物から反射された異なるＲＦ周波数の信号が複数の受信アンテナで受信される。それらの受信アンテナで受信された異なるＲＦ周波数のパルス信号は、位相が揃えられ、その後、（コヒーレントに）合成される。このような合成処理により、ＲＦ周波数に応じた位置に仮想的なアンテナを形成する事ができ、結果、実際の受信アンテナの間隔を波長λの半分以上に拡張する事が可能となる。【００１９】
電波を用いたイメージング装置の他の方式として、図１４で概念図を示した合成開口レーダー（Ｓｙｎｔｈｅｔｉｃ  ａｐｅｒｔｕｒｅ  ｒａｄａｒ，  ＳＡＲ）方式がある。図１４は、一般的な合成開口レーダー方式の概念を示す概念図である。図１４に示すように、合成開口レーダー方式において、測定装置は、送信機３１１と、受信機３０１とを備えている。送信機３１１は、送信アンテナ３１２を備えている。また、受信機３０１は、受信アンテナ３０２を備えている。【００２０】
送信機３１１は、送信アンテナ３１２から、ＲＦ信号（電波）３１３を、検知対象物３０４１、３０４２、・・・、３０４Ｋ（Ｋは対象物の数）に向けて照射する。ＲＦ信号（電波）３１３は、検知対象物３０４１、３０４２、・・・、３０４Ｋにおいて反射され、反射波３０３１、３０３２、・・・、３０３Ｋがそれぞれ発生する。この場合、受信機３０１は、受信機３０１１、３０１２、・・・、３０１Ｎの位置に移動しながら、受信アンテナ３０２１、３０２２、・・・、３０２Ｎ（の位置）において反射波３０３１、３０３２、・・・、３０３Ｋを受信する。【００２１】
更に、この場合、受信アンテナ３０２１、３０２２、・・・、３０２Ｎは、図１２で示したアレイアンテナ方式における受信アンテナ２０２１、２０２２、・・・、２０２Ｎと同じく、Ｎ本のアンテナによる受信アレイアンテナ（仮想アレイアンテナ）を形成する。【００２２】
従って、図１２で示したアレイアンテナ方式と同じく、図１４で示した合成開口レーダー方式においても、受信機３０１は、受信した反射波３０３１、３０３２、・・・、３０３Ｋに基づいて、検知対象物３０４１、３０４２、・・・、３０４Ｋから反射されている電波強度を算出する。また、合成開口レーダー方式においても、その電波強度の分布を画像化する事で、検知対象物３０４１、３０４２、・・・、３０４Ｋの像を得る事ができる。【００２３】
また、合成開口レーダー方式による電波イメージング装置の例は、特許文献５〜７において開示されている。
【先行技術文献】
【特許文献】
【００２４】
【特許文献１】特表２０１３−５２８７８８号公報
【特許文献２】特開２０１５−０１４６１１号公報
【特許文献３】特許第５０８０７９５号
【特許文献４】特開２０１５−２３０２１６号公報
【特許文献５】特許第４６５３９１０号公報
【特許文献６】特表２０１１−５１３７２１号公報
【特許文献７】特開２０１５−０３６６８２号公報
【発明の概要】
【発明が解決しようとする課題】
【００２５】
ここで、上記特許文献１〜３に開示されたアレイアンテナ方式（図１２参照）の問題点について説明する。まず、アレイアンテナ方式の問題点は、必要となる受信アンテナ２０２１、２０２２、・・・、２０２Ｎの数が増加すると、それに付随して受信機の数Ｎも増加して、結果、装置のコスト、サイズ、及び重量が大きくなる点である。【００２６】
上記の問題点について以下に具体的に説明する。まず、アレイアンテナ方式の場合、受信アンテナ２０２１、２０２２、・・・、２０２Ｎの各アンテナの間隔は、受信機２０１において受信される反射波２０３１、２０３２、・・・、２０３Ｋの波長λの半分以下で無ければならない、という条件がある。この条件を満たせない場合、生成した画像において、検知対象物２０４１、２０４２、・・・、２０４Ｋが存在しない位置に虚像が発生するという問題が生じる。反射波２０３１、２０３２、・・・、２０３Ｋがミリ波である場合、波長は数ｍｍ程度であるので、アンテナの間隔は極めて小さい。【００２７】
次に、アレイアンテナ方式の場合、画像の分解能は、受信アレイアンテナ（２０２１、２０２２、・・・、２０２Ｎ）の指向性ビーム幅Δθで決まる。受信アレイアンテナ（２０２１、２０２２、・・・、２０２Ｎ）の指向性ビームの幅Δθは、Δθ〜λ／Ｄにて与えられる。ここで、Ｄは受信アレイアンテナ（２０２１、２０２２、・・・、２０２Ｎ）の開口サイズであり、両端に存在する受信アンテナ２０２１と受信アンテナ２０２Ｎとの間の距離に対応する。従って、衣服の下にある物品、鞄の中にある物品の画像化において実用的な分解能を得るには、受信アレイアンテナ（２０２１、２０２２、・・・、２０２Ｎ）の開口サイズＤを数十ｃｍから数ｍ程度取る必要がある。【００２８】
このように、アレイアンテナ方式には２つの条件がある。即ち、受信アンテナ２０２１、２０２２、・・・、２０２Ｎの各アンテナの間隔は波長λの半分以下（数ｍｍ以下）で無ければならないという事と、両端に存在する受信アンテナ２０２１と受信アンテナ２０２Ｎとの間の距離は少なくとも数十ｃｍ程度必要であるという事とである。このため、アレイアンテナ方式を採用すると、一列当たりに必要なアンテナの数Ｎは、数百個程度となる。【００２９】
また、実際の電波を用いたイメージング装置は、２次元画像を形成する必要がある。このため、図１５に示すように縦方向と横方向とに、それぞれＮ個ずつの受信アンテナ２０２を敷き詰める必要があり、全体で必要なアンテナの数はＮ２個となる。従って、全体で必要なアンテナの数と、それに付随する受信機の数とは、数万個程度となる。図１５は、一般的なアレイアンテナ方式のイメージング装置におけるアンテナのレイアウトの例を示す構成図である。【００３０】
このように、従来からのアレイアンテナ方式のイメージング装置では、大量のアンテナと受信機とが必要となるため、コストは非常に高いものになる。また、数十ｃｍ〜数ｍ四方の領域にアンテナが敷き詰められるので、装置のサイズ及び重量は非常に大きなものとなる。【００３１】
一方、上記の特許文献４に開示されたアレイアンテナ方式は、ＲＦ周波数の異なる複数のパルス信号を送信し、ＲＦ周波数に応じた位置に仮想的なアンテナを形成する事で、実際の受信アンテナの間隔を波長λの半分以上に拡張しながら虚像の抑制を実現している。【００３２】
ところが、特許文献４に開示されたアレイアンテナ方式には、高分解能を実現するためにＲＦ周波数の範囲を広く設定しなければならないという問題点がある。具体的には、特許文献４に開示されたアレイアンテナ方式では、ＲＦ周波数ｆ０、ｆ１、・・・、ｆＮのパルス信号を送信し、且つ間隔ｄを空けて設置された２本の受信アンテナで受信を行った場合、ｄ・ｆｎ/ｆ０（ｎ＝０、１、・・・、Ｎ）の位置に仮想アンテナが設置される。この時の実効的な開口サイズＤは、ｄ・（ｆＮ−ｆ０）/ｆ０となる。実効的な開口サイズＤから、角度分解能Δθは以下の式（２）で与えられる。【００３３】
【数２】
0006939981.tif 000003
【００３４】
上記式（２）から分かるように、良好な角度分解能Δθを得るには、ＲＦ周波数の設定範囲ｆＮ−ｆ０を大きく取る必要がある。しかしながら、電波法の制約上、ＲＦ周波数の設定範囲ｆＮ−ｆ０を大きく取る事は困難であり、ｆ０/(ｆＮ−ｆ０)の値は、通常１０〜１００程度となる。ここで、特許文献４の方式における受信アンテナの間隔ｄと、特許文献１〜３に開示されたアレイアンテナ方式の開口サイズＤとが同じあるとする。そうすると、特許文献４に開示されたアレイアンテナ方式における角度分解能は、特許文献１〜３に開示されたアレイアンテナ方式の角度分解能に比べ、ｆ０/(ｆＮ−ｆ０)倍、即ち、１０〜１００倍程度に劣化する事となる。【００３５】
更に、特許文献４に開示されたアレイアンテナ方式では、パルス信号が用いられているが、距離分解能を得るためには、パルス信号のＲＦ帯域幅は、数百ＭＨｚから数ＧＨｚに設定されている必要がある。加えて、特許文献４において採用されているパルスレーダ方式では、ベースバンド信号の帯域幅はＲＦ帯域幅と同一である。そのため、ベースバンド信号のサンプリング周波数を数百ＭＨｚから数ＧＨｚにまで高速化しなければならず、この事は演算量の増加につながる。また、一般にベースバンド信号の帯域幅に比例した雑音がベースバンド信号に混入するが、パルスレーダ方式では、ベースバンド信号の帯域幅が広い分、雑音による信号品質の劣化も問題となる。【００３６】
続いて、上記特許文献５〜７に開示された合成開口レーダー方式の問題点について以下に説明する。合成開口レーダー方式の問題点は、装置を機械的に動かす必要があり、そのために走査時間の短縮が難しいという事である。この事は、物品又は人をイメージング装置で検査する場合に、単位時間あたりに検査できる対象物の数が限られるという問題につながる。また、特に、特許文献７に開示されたイメージング装置では、受信機を動かすための機械的な機構が必要となるので、よりいっそう装置のサイズ及び重量が増大してしまう。【００３７】
以上のように、一般的に、電波を用いたイメージング装置では、装置のコスト、サイズ、重量が非常に大きなものになるという問題が発生する。一方、この問題を解決しようとすると、要求される分解能及び信号品質が得られず、精度が悪いという問題が発生してしまう。このため、イメージング装置を実際に使用できる用途及び機会は、限定されてしまうという問題が発生する。また、対象物を検査する速度が限られたものになるという問題も発生する。【００３８】
本発明の目的の一例は、上記問題を解消し、電波による物体のイメージングにおいて、画質の低下を抑制しつつ、装置の低コスト化及び小型化を図り得る、物体検知装置、物体検知方法、及びプログラムを提供することにある。
【課題を解決するための手段】
【００３９】
上記目的を達成するため、本発明の一側面における物体検知装置は、電波によって物体を検知するための物体検知装置であって、
前記物体に向けて、送信信号となる電波を照射する、送信部と、
複数の受信アンテナによって、前記物体で反射された電波を受信信号として受信し、更に、複数の前記受信アンテナそれぞれで受信された前記受信信号毎に、当該受信信号を用いて、中間周波数信号を生成する、受信部と、
前記中間周波数信号それぞれを合成して得られるビームパターンによる虚像の発生が抑制されるように、サンプリング時刻を決定し、決定したサンプリング時刻で前記中間周波数信号に対してサンプリングを行うことによって、物体の位置検出用の中間周波数信号を生成し、
生成した前記位置検出用の中間周波数信号を用いて、前記物体を検知する、演算部と、
を備えている、
ことを特徴とする。【００４０】
また、上記目的を達成するため、本発明の一側面における物体検知方法は、電波によって物体を検知するための方法であって、
（ａ）前記物体に向けて、送信信号となる電波を照射する、ステップと、
（ｂ）複数の受信アンテナによって、前記物体で反射された電波を受信信号として受信し、更に、複数の前記受信アンテナそれぞれで受信された前記受信信号毎に、当該受信信号を用いて、中間周波数信号を生成する、ステップと、
（ｃ）前記中間周波数信号それぞれを合成して得られるビームパターンによる虚像の発生が抑制されるように、サンプリング時刻を決定し、決定したサンプリング時刻で前記中間周波数信号に対してサンプリングを行うことによって、物体の位置検出用の中間周波数信号を生成する、ステップと、
（ｄ）生成した前記位置検出用の中間周波数信号を用いて、前記物体を検知する、ステップと、を有している、
ことを特徴とする。【００４１】
更に、上記目的を達成するため、本発明の一側面におけるプログラムは、
物体に向けて、送信信号となる電波を照射する、送信部と、複数の受信アンテナによって、前記物体で反射された電波を受信信号として受信し、更に、複数の前記受信アンテナそれぞれで受信された前記受信信号毎に、当該受信信号を用いて、中間周波数信号を生成する、受信部と、コンピュータと、を備えた物体検知装置において、
前記コンピュータに、
（ａ）前記中間周波数信号それぞれを合成して得られるビームパターンによる虚像の発生が抑制されるように、サンプリング時刻を決定し、決定したサンプリング時刻で前記中間周波数信号に対してサンプリングを行うことによって、物体の位置検出用の中間周波数信号を生成する、ステップと、
（ｂ）生成した前記位置検出用の中間周波数信号を用いて、前記物体を検知する、ステップと、を実行させる命令を含む、
ことを特徴とする。
【発明の効果】
【００４２】
以上のように、本発明によれば、電波による物体のイメージングにおいて、画質の低下を抑制しつつ、装置の低コスト化及び小型化を図ることができる。
【図面の簡単な説明】
【００４３】
【図１】図１は、本発明の実施の形態１における物体検知装置の概略構成を示す図である。
【図２】図２は、本発明の実施の形態１における物体検知装置の送信部及び受信部の構成を具体的に示す図である。
【図３】図３は、本発明の実施の形態１における物体検知装置の送信部及び受信部の他の例の構成を具体的に示す図である。
【図４】図４は、本発明の実施の形態１における物体検知装置の送信部及び受信部の他の例の構成を具体的に示す図である。
【図５】図５は、本発明の実施の形態１における物体検知装置が照射する送信信号の一例を示す図である
【図６】図６は、本発明の実施の形態１において検知の対象となる物体と検知を行うための送信アンテナ及び受信アンテナとの位置関係の一例を示す図である。
【図７】図７は、本発明の実施の形態１において検知の対象となる物体と検知を行うための送信アンテナ及び受信アンテナとの位置関係の他の例を示す図である。
【図８】図８は、本発明の実施の形態１における物体検知装置の動作を示すフロー図である。
【図９】図９は、本発明の実施の形態１で用いられる測距受信アンテナの配置と所望の時間サンプリング点の取り方を示す図である。
【図１０】図１０は、本発明の実施の形態２において検知の対象となる物体と検知を行うための送信アンテナ及び受信アンテナとの位置関係の一例を示す図である。
【図１１】図１１は、本発明の実施の形態１及び２における物体検知装置を実現するコンピュータの一例を示すブロック図である。
【図１２】図１２は、一般的なアレイアンテナ方式の概念を示す概念図である。
【図１３】図１３は、図１２に示した受信機の構成をより具体的に示す図である。
【図１４】図１４は、一般的な合成開口レーダー方式の概念を示す概念図である。
【図１５】図１５は、一般的なアレイアンテナ方式のイメージング装置におけるアンテナのレイアウトの例を示す構成図である。
【発明を実施するための形態】
【００４４】
以下、本発明による送信装置および送信方法の好適な実施形態について添付図を参照して説明する。なお、以降に示す各図面において、同一または相当部分の部位については、同一符号を付して示すこととし、その説明は繰り返さないことにする。
【００４５】
（実施の形態）
以下、本発明の実施の形態１における、物体検知装置、物体検知方法、及びプログラムについて、図１〜図９を参照しながら説明する。本実施の形態１は、小型のレーダー装置を用いながら、検知対象となる物体の位置だけでなく、この物体の幅といった形状に関する情報も検知できる例について示す。
【００４６】
［装置構成］
最初に、図１を用いて、本実施の形態１における物体検知装置の構成について説明する。図１は、本発明の実施の形態１における物体検知装置の概略構成を示す図である。
【００４７】
図１に示す、本実施の形態１における物体検知装置１０００は、電波によって物体を検知するための装置である。図１に示すように、物体検知装置１０００は、送信部１１０１と、受信部１１０２と、演算装置１２１１とを備えている。
【００４８】
また、本実施の形態では、送信部１１０１と受信部１１０２とで送受信装置１００１が構成されており、物体検知装置１０００は、主に、送受信装置１００１と演算装置１２１１とで構成されている。また、演算装置１２１１は、物体検知装置１０００において、演算部として機能する。
【００４９】
送信部１１０１は、対象物配置面１００５上に存在している物体であって、検知対象となる物体（以下、「対象物」と表記する）１００３に向けて、送信信号となる電波を照射する。
【００５０】
受信部１１０２は、複数の受信アンテナ（後述の図２参照）によって、対象物１００３で反射された電波を受信信号として受信する。また、受信部１１０２は、更に、受信アンテナそれぞれで受信された受信信号毎に、各受信信号を用いて、中間周波数信号（以下「ＩＦ（Intermediate Frequency）信号」と表記する。）を生成する。
【００５１】
具体的には、図１に示すように、送信部１１０１は、受信部１１０２に向けて、端子１２０８を経由して送信信号を出力する。受信部１１０２は、対象物１００３で反射された受信信号と、端子１２０８を経由して出力されてきた送信信号とを、ミキシングして、ＩＦ信号を生成する。また、送信部１１０１は、生成したＩＦ信号を演算装置１２１１に出力する。
【００５２】
演算装置（演算部）１２１１は、ＩＦ信号それぞれを合成して得られるビームパターンによる虚像の発生が抑制されるように、サンプリング時刻を決定する。また、演算装置１２１１は、決定したサンプリング時刻でＩＦ信号に対してサンプリングを行うことによって、対象物１００３の位置検出用のＩＦ信号（以下「位置検出用ＩＦ信号」と表記する。）を生成する。そして、演算部１２１１は、位置検出用ＩＦ信号を用いて、対象物１００３を検知する。
【００５３】
このように本実施の形態１では、ＩＦ信号を、適切なサンプリング時刻でサンプリングすることによって、得られたサンプリングデータを、仮想的なアンテナで受信された受信信号からのＩＦ信号として用いることができる。つまり、本実施の形態１では、実際には存在しない仮想的なアンテナを、一般的なアレイアンテナ方式における受信アンテナと置き換えることができ、実際よりも多数の受信アンテナを備えていた場合と同様の高分解能を実現できる。このため、本実施の形態１によれば、電波による物体のイメージングにおいて、画質の低下を抑制しつつ、装置の低コスト化及び小型化を図ることができる。
【００５４】
また、図１においては、送信部１１０１及び受信部１１０２は、それぞれ一つのみが図示されているが、送信部１１０１及び受信部１１０２は、実際には複数備えられていても良い。送信部１１０１及び受信部１１０２が複数備えられている場合は、複数の受信部１１０２それぞれは、いずれかの送信部１１０１に対応している。
【００５５】
続いて、図１に加えて、図２から図４までを用いて、本実施の形態１における物体検知装置１０００の構成について更に具体的に説明する。図２は、本発明の実施の形態１における物体検知装置の送信部及び受信部の構成を具体的に示す図である。図３は、本発明の実施の形態１における物体検知装置の送信部及び受信部の他の例の構成を具体的に示す図である。図４は、本発明の実施の形態１における物体検知装置の送信部及び受信部の他の例の構成を具体的に示す図である。
【００５６】
図２に示すように、本実施の形態１では、送受信装置１００１において、送信部１１０１は、発振器１２０１と、送信アンテナ１２０２とを備えている。また、受信部１１０２は、受信アンテナ１２０３と、ミキサ１２０４と、インターフェイス回路１２０５とを備えている。更に、図１でも示したように、送信部１１０１と受信部１１０２とは、端子１２０８を介して接続されている。
【００５７】
送信部１１０１において、発振器１２０１は、ＲＦ信号（電波）を生成する。発振器１２０１で生成されたＲＦ信号は、送信アンテナ１２０２から送信信号として送信され、対象物１００３に照射される。対象物１００３で反射された電波は、受信部１１０２において、受信アンテナ１２０３によって受信される。
【００５８】
ミキサ１２０４は、発振器１２０１から端子１２０８を経由して入力されてきたＲＦ信号と受信アンテナ１２０３で受信された電波（受信信号）とを、ミキシングする事で、ＩＦ信号を生成する。ミキサ１２０４で生成されたＩＦ信号は、インターフェイス回路１２０５を経由して、演算装置１２１１へと送信される。インターフェイス回路１２０５は、アナログ信号であるＩＦ信号を、演算装置１２１１で扱えるデジタル信号に変換する機能を持ち、得られたデジタル信号を演算装置１２１１へと出力する。
【００５９】
また、図２に示した例では、一つの送信部１１０１には、一つの送信アンテナ１２０２が備えられているが、本実施の形態は、この態様に限定されない。本実施の形態では、例えば、図３に示すように、一つの送信部１１０１において、複数の送信アンテナ１２０２が備えられていても良い。
【００６０】
具体的には、図３の例では、送信部１１０１は、一つの発振器１２０１と、複数の送信アンテナ１２０２とを備えている。また、送信部１１０１は、送信アンテナ１２０２毎に設けられた可変移相器１２０６も備え、各送信アンテナ１２０２は、可変移相器１２０６を介して、発振器１２０１に接続されている。各可変移相器１２０６は、発振器１２０１から送信アンテナ１２０２の各々に供給される送信信号の位相を制御する事で、送信アンテナ１２０２の指向性の制御を行なっている。
【００６１】
本実施の形態1では、例えば、図４に示すように、送信アンテナ１２０２の代わりに、送受信共用アンテナ１２１０が備えられていてもよい。具体的には、図４の例では、送信部１１０１において、送受信共用アンテナ１２１０とサーキュレータ１２０９が備えられている。サーキュレータ１２０９は、発振器１２０１から供給される送信信号を送受信共用アンテナ１２１０に出力する。送受信共用アンテナ１２１０は、サーキュレータ１２０９から出力された送信信号を対象物１００３に照射する。送受信共用アンテナ１２１０で受信した対象物１００３からの反射信号は、サーキュレータ１２０９を介してミキサ１２０４へと出力される。
【００６２】
また、本実施の形態１では、受信部１１０２は、基準受信アンテナと測定受信アンテナとを備えることができる。この場合、演算装置１２１１は、測定受信アンテナで受信された受信信号から生成された中間周波数信号と、前記基準受信アンテナで受信された受信信号から生成された中間周波数信号とに対して、サンプリングを実行する。
【００６３】
そして、演算装置１２１１は、得られたサンプリング値について、測定受信アンテナで受信された受信信号から生成されたＩＦ信号を、基準受信アンテナで受信された受信信号から生成されたＩＦ信号によって、規格化して、位置検出用ＩＦ信号を生成する。
【００６４】
また、演算装置１２１１は、測定受信アンテナで受信された受信信号から生成されたＩＦ信号と、基準受信アンテナで受信された受信信号から生成されたＩＦ信号とに対し、対象物１００３の位置（物体検知装置１０００までの距離）を制限する処理を実行することもできる。この場合は、演算装置１２１１は、処理が実行されたＩＦ信号を用いて、サンプリングを実行する。
【００６５】
更に、本実施の形態では、演算装置１２１１は、位置検出用ＩＦ信号から相関行列を算出し、算出した相関行列から、対象物１００３の位置分布を表す評価関数を算出する。そして、演算装置１２１１は、算出した評価関数を用いて、対象物１００３の位置及び形状を検知する。
【００６６】
ここで、図５を用いて、本実施の形態において物体に照射される送信信号について説明する。図５は、本発明の実施の形態１における物体検知装置が照射する送信信号の一例を示す図である。
【００６７】
本実施の形態１において、発振器１２０１で生成されるＲＦ信号は、図５に示すように、周期Ｔchirpで、ＲＦ周波数がｆｍｉｎからｆｍｉｎ＋ＢＷ（＝ｆｍａｘ）に時間的に変化する、ＦＭＣＷ信号であることが好ましい。なお、ｆｍｉｎは、ＲＦ周波数の最小値であり、ＢＷはＲＦ信号の帯域幅である。また、ＲＦ周波数の時間変化は連続的と離散的のいずれでも良い。
【００６８】
また、ここで、図６及び図７を用いて、対象物１００３と、送信アンテナ１２０２及び受信アンテナ１２０３との位置関係について説明する。図６は、本発明の実施の形態１において検知の対象となる物体と検知を行うための送信アンテナ及び受信アンテナとの位置関係の一例を示す図である。図７は、本発明の実施の形態１において検知の対象となる物体と検知を行うための送信アンテナ及び受信アンテナとの位置関係の他の例を示す図である。
【００６９】
図６の例では、受信アンテナ１２０３は、基準受信アンテナ１２０３０と測定受信アンテナ１２０３ｍ（ｍ＝１、２、・・・、Ｍ）とで構成されている。
【００７０】
図６で示すように、送信アンテナ１２０２と基準受信アンテナ１２０３０を焦点とする楕円状の対象物配置面１００５上にＫ個の対象物１００３ｋ（ｋ＝１、２、・・・、Ｋ）が存在する状況を考える。また、送信アンテナ１２０２から対象物１００３ｋを経由して基準受信アンテナ１２０３０に至るまでの経路長を２Ｒ０とする。また、送信アンテナ１２０２と受信アンテナ１２０３ｍ（ｍ＝０、１、２、・・・、Ｍ）との間隔をｄｍとする。対象物１００３ｋ（ｋ＝１、２、・・・、Ｋ）からの反射波は、各受信アンテナ１２０３ｍ（ｍ＝１、２、・・・、Ｍ）において振幅skで受信されるものとする。
【００７１】
また、図６の例では、送信アンテナ１２０２と基準受信アンテナ１２０３０とは、個別に設置されている。一方、図７に示す例では、送信アンテナ１２０２と基準受信アンテナ１２０３０との代わりに、送受信共用アンテナ１２１０が用いられている。送受信共用アンテナ１２１０は、送信アンテナ１２０２及び基準受信アンテナ１２０３０両方の機能を有し、両者としての役割を果たす。また、図７の例では、対象物１００３ｋ（ｋ＝１、２、・・・、Ｋ）は、送受信共用アンテナ１２１０を中心とする半径Ｒ０の円上に配置されているとする。
【００７２】
［装置動作］
次に、本実施の形態１における物体検知装置１０００の動作について図８及び図９を用いて説明する。図８は、本発明の実施の形態１における物体検知装置の動作を示すフロー図である。以下の説明においては、適宜図１〜図７を参酌する。また、本実施の形態１では、物体検知装置１０００を動作させることによって、物体検知方法が実施される。よって、本実施の形態１における物体検知方法の説明は、以下の物体検知装置１０００の動作説明に代える。
【００７３】
図８に示すように、最初に、送受信装置１００１において、送信部１１０１から、対象物１００３に向けて、送信信号となる電波を照射する（ステップＡ１）。また、送信部１１０１は、送信信号となる電波の照射と同時に、端子１２０８を介して、送信信号を受信部１１０２に出力する。
【００７４】
次に、送受信装置１００１において、受信部１１０２の受信アンテナ１２０３が、対象物１００３から反射された電波を、受信信号として受信する（ステップＡ２）。
【００７５】
次に、受信部１１０２は、各受信アンテナ１２０３で受信した受信信号それぞれに、送信部１１０１で生成された送信信号をそれぞれミキシングして、ＩＦ信号を生成する（ステップＡ３）。
【００７６】
次に、演算装置１２１１が、ステップＡ３で各受信アンテナ１２０３の受信信号から得たＩＦ信号それぞれから、対象物１００３の物体検知装置１０００までの距離を特定の値に制約したＩＦ信号を生成する（ステップＡ４）。具体的には、演算装置１２１１は、測定受信アンテナで受信された受信信号から生成されたＩＦ信号と、基準受信アンテナで受信された受信信号から生成されたＩＦ信号とに対し、対象物１００３の検知範囲を制限する処理を実行する。
【００７７】
次に、演算装置１２１１は、測定受信アンテナ１２０３ｍ（ｍ＝１、２、・・・、Ｍ）のＩＦ信号と、基準受信アンテナ１２０３０のＩＦ信号とを用いて、位置検出用ＩＦ信号を生成する（ステップＡ５）。
【００７８】
具体的には、演算装置１２１１は、ＩＦ信号それぞれを合成して得られるビームパターンによる虚像の発生が抑制されるように、サンプリング時刻を決定する。そして、演算装置１２１１は、決定したサンプリング時刻の値が得られるように、補間処理を行って、測定受信アンテナ１２０３ｍ（ｍ＝１、２、・・・、Ｍ）のＩＦ信号の値と、基準受信アンテナ１２０３０のＩＦ信号の値とを算出する。その後、演算装置１２１１は、算出した値を用いて、測定受信アンテナＩＦ信号を基準受信アンテナＩＦ信号によって規格化して、規格化ＩＦ信号を生成する。この生成された規格化ＩＦ信号が、位置検出用ＩＦ信号となる。
【００７９】
次に、演算装置１２１１は、規格化ＩＦ信号から評価関数を算出し、さらに評価関数から対象物１００３の位置（距離および角度）を算出する（ステップＡ６）。
【００８０】
次に、演算装置１２１１は、予め定められている距離全てについて処理を実施しているかどうかを判定する（ステップＡ７）。ステップＡ７の判定の結果、Ｙｅｓである場合は、演算装置１２１１は動作を終了する。一方、ステップＡ７の判定の結果、Ｎｏである場合は、演算装置１２１１は、再度ステップＡ４を実行する。
【００８１】
続いて、図８に示したステップＡ１〜Ａ７のうちステップＡ３〜Ａ７について、より詳細に説明する。
【００８２】
［ステップＡ３］
まず、各受信アンテナ１２０３で受信した受信信号に、送信部１１０１で生成された送信信号をそれぞれミキシングしてＩＦ信号を生成する、ステップＡ３の詳細について説明する。
【００８３】
図６又は図７に示した配置において、基準受信アンテナ１２０３０又は送受信共用アンテナ１２１０の受信信号から算出される基準ＩＦ信号（ＩＦ０（ｔ））は、以下の式（３）で与えられる。
【００８４】
【数３】
0006939981.tif 000004
【００８５】
上記式（３）において、αは、図５で示したチャープ信号のチャープ率を示しており、α＝ＢＷ/Ｔchirpで与えられる。cは光速、tは時刻を示している。
【００８６】
また、測定受信アンテナ１２０３ｍ（ｍ＝１、２、・・・、Ｍ）の受信信号から算出される測定ＩＦ信号（ＩＦｍ（ｔ））は、以下の式（４）で与えられる。以下の式（４）において、Δｄｍ = ｄｍ−ｄ０である。
【００８７】
【数４】
0006939981.tif 000005
【００８８】
［ステップＡ４］
次に、演算装置１２１１が、ステップＡ３で各受信アンテナ１２０３の受信信号から得たＩＦ信号に基づいて、対象物１００３の物体検知装置１０００までの距離を特定の値に制約したＩＦ信号を生成する、ステップＡ４の詳細について説明する。
【００８９】
図６又は図７に示した例では、送信アンテナ１２０２から対象物１００３を経由して基準受信アンテナ１２０３０に至るまでの経路長が２Ｒ０となる楕円状（又は円状）の対象物配置面１００５上に、対象物１００３が配置されている。実際の状況では、距離Ｒ０で指定される対象物配置面１００５以外の位置に対象物が存在する事もある。そこでステップＡ４では、演算装置１２１１は、距離Ｒ０で指定される対象物配置面１００５上の位置に存在する対象物１００３からの信号のみを抽出する。
【００９０】
例えば、対象物配置面１００５において、距離Ｒ０の位置に、対象物１００３が存在しているとする。この場合、対象物１００３により生成される式（３）の基準ＩＦ信号（ＩＦ０（ｔ））は、２αＲ０/cの周波数を有することになる。従って、基準受信アンテナ１２０３０又は送受信共用アンテナ１２１０で得た受信ＩＦ信号の時間波形をフーリエ変換し、２αＲ０/cの周波数の信号成分のみを抽出すれば、対象物１００３の位置を所望の位置（距離範囲）に制約したＩＦ信号が得られる。
【００９１】
上記式（４）の測定ＩＦ信号（ＩＦｍ（ｔ））は、「２α（Ｒ０−Δｄｍ）/c」から「２α（Ｒ０＋Δｄｍ）/c」までの周波数を有する。従って、測定受信アンテナ１２０３ｍ（ｍ＝１、２、・・・、Ｍ）で得た測定ＩＦ信号の時間波形をフーリエ変換し、２α（Ｒ０−Δｄｍ）/cから２α（Ｒ０＋Δｄｍ）/cまでの周波数の信号成分のみを抽出すれば、距離の値がＲ０に制約されたＩＦ信号が得られる。
【００９２】
［ステップＡ５］
次に、測定受信アンテナＩＦ信号を、基準受信アンテナＩＦ信号によって規格化して、規格化ＩＦ信号を生成する、ステップＡ５の詳細について説明する。
【００９３】
図９は、本実施の形態１で用いられる測定受信アンテナ１２０３ｍ（ｍ＝１、２、・・・、Ｍ）の配置と所望の時間サンプリング点の取り方を示す図である。図９の下側の図は、従来方式の受信アンテナ２０２の配置も示している。
【００９４】
背景技術の欄で述べたとおり、従来方式のアンテナアレイでは、虚像の発生を防ぐために受信アンテナ２０２ｎ（ｎ＝１、２、・・・、Ｎ）の間隔を、受信波の波長λの半分以下としなければならない。また、高分解能の実現のために、アレイ開口サイズ（受信アンテナ２０２１から２０２Ｎまでの長さ）を大きく取る必要があり、結果として多数の受信アンテナ２０２が必要となっていた。
【００９５】
本実施の形態１では、測定ＩＦ信号（ＩＦｍ（ｔ））を基準ＩＦ信号（ＩＦ０（ｔ））で除算して得られる規格化ＩＦ信号ｒｍ（ｔ）の時間波形のサンプリングデータを仮想的なアンテナ列１３０１ｍ（ｍ＝１、２、・・・、Ｍ）とする。そして仮想的なアンテナ列１３０１ｍ（ｍ＝１、２、・・・、Ｍ）を、従来方式の受信アンテナ２０２ｎ（ｎ＝１、２、・・・、Ｎ）の代用とする事で、実際の受信アンテナ１２０３ｍ（ｍ＝１、２、・・・、Ｍ）の数Ｍを小さく収めながら、高分解能を実現する。
【００９６】
本実施の形態１では、虚像の発生を防ぐために、時間波形のサンプリング時刻および実際の受信アンテナ１２０３ｍ（ｍ＝１、２、・・・、Ｍ）の位置を適正に設定する必要がある。以下では、虚像の発生を防ぐ時間波形のサンプリング時刻および実際の受信アンテナ１２０３ｍ（ｍ＝１、２、・・・、Ｍ）の位置の設定例を示す。
【００９７】
ステップＡ５では、測定ＩＦ信号（ＩＦｍ（ｔ））を基準ＩＦ信号（ＩＦ０（ｔ））で除算する事で、式（５）で与えられる規格化ＩＦ信号ｒｍ（ｔ）が生成される。
【００９８】
【数５】
0006939981.tif 000006
【００９９】
受信アンテナ１２０３ｍ（ｍ＝１、２、・・・、Ｍ）で得られる規格化ＩＦ信号ｒｍ（ｔ）を、サンプリング時刻ｔ（ｍ,１），ｔ（ｍ,2）,・・・，ｔ（ｍ,p(m)）において取得する。ｐ（ｍ）は、受信アンテナ１２０３ｍで得られる規格化ＩＦ信号のサンプリング点数である。全てのアンテナのサンプリング点数の総和ｐ（１）＋ｐ（２）＋・・・＋ｐ（Ｍ）をＮとする。
【０１００】
規格化ＩＦ信号ベクトルｒを以下のように定義する。
r = [r1,1,・・・, r1,p(1),・・・,rM,1,・・・, rM,p(M)]T,
ここで、rm（ｔ（ｍ,n））= rm,nとしている。rは、Ｎ次ベクトルである。規格化ＩＦ信号ベクトルrのn番目の成分rnは、以下の式（６）及び（７）のように表される。
【０１０１】
【数６】
0006939981.tif 000007
【０１０２】
【数７】
0006939981.tif 000008
【０１０３】
上記式（７）における時刻tnは、n番目のサンプリング時刻である。アンテナ位置Δdmは、nの値によって変化する。
【０１０４】
次に、虚像の発生条件について説明する。角度θ0の方向に対象物１００３が存在し、且つ、以下の式（８）の条件が満たされる時、対象物１００３が本来は存在しない角度θaの方向に、対象物の虚像が発生する。
【０１０５】
【数８】
0006939981.tif 000009
【０１０６】
次に、前記の虚像を排除するため、以下の式（９）の条件が設定される。
【０１０７】
【数９】
0006939981.tif 000010
【０１０８】
上記式（９）において、Δtmは、受信アンテナ１２０３ｍ（ｍ＝１、２、・・・、Ｍ）で得られる規格化ＩＦ信号ｒｍ（ｔ）のサンプリング周期である。角度θ0およびθaが−πから＋πの範囲で任意の値を取っても上記式（９）が満たされるサンプリング周期Δtmの条件は、以下の式（１０）で与えられる。
【０１０９】
【数１０】
0006939981.tif 000011
【０１１０】
上記式（１０）が示すように、本実施の形態１では、上限を受信アンテナ１２０３ｍ（ｍ＝１、２、・・・、Ｍ）の位置Δdmに応じて定めると、この上限に収まるようにサンプリング周期Δtmを設定する事が望ましい。
【０１１１】
また、隣接する受信アンテナ１２０３ｍおよび１２０３ｍ＋１において、上記式（７）の位相φn(θ)の範囲が連続となるように、以下の式（１１）の条件が設定される。
【０１１２】
【数１１】
0006939981.tif 000012
【０１１３】
本実施の形態１では、式（１１）の条件が満たされるように、アンテナ位置Δdm(ｍ＝１、２、・・・、Ｍ−１)を定める事が望ましい態様である。
【０１１４】
なお、隣接する受信アンテナ１２０３ｍと１２０３ｍ＋１との間隔を、Δd’m(＝Δdm+1−Δdm)とした場合、上記式(１１)は、以下の式(１２)と等価である。
【０１１５】
【数１２】
0006939981.tif 000013
【０１１６】
すなわち、隣接する受信アンテナ１２０３ｍと１２０３ｍ＋１との間隔は、最小ＲＦ周波数fminと最大ＲＦ周波数fmaxとで定まる上限値以下に収める事が、本実施の形態１では好ましい態様である。
【０１１７】
ステップＡ５では、演算装置１２１１は、まず測定ＩＦ信号（ＩＦｍ（ｔ））と基準ＩＦ信号（ＩＦ０（ｔ））とのそれぞれに対し、補間処理を行い、上記式（１０）を満たすサンプリング周期Δtmのサンプリング時刻における値を算出する。次に、演算装置１２１１は、補間処理により測定ＩＦ信号（ＩＦｍ（ｔ））と基準ＩＦ信号（ＩＦ０（ｔ））とのサンプリング時刻を揃えた後に、上記式(５)で示すように測定ＩＦ信号（ＩＦｍ（ｔ））を基準ＩＦ信号（ＩＦ０（ｔ））で除算して規格化ＩＦ信号ｒｍ（ｔ）(ｍ＝１、２、・・・、Ｍ)を算出する。
【０１１８】
［ステップＡ６］
次に、規格化ＩＦ信号から評価関数を算出するステップＡ６の詳細について説明する。
【０１１９】
ステップＡ６では、演算装置１２１１、規格化ＩＦ信号の範囲を複数選択し、選択した複数の範囲を定めた複数の規格化ＩＦ信号それぞれから、相関行列を算出し、そして、相関行列の平均から、対象物１００３の位置分布を表す評価関数を算出する。この点について以下に説明する。
【０１２０】
まず、ステップＡ６では、演算装置１２１１は、規格化ＩＦ信号ベクトルｒ=[r1,r2,・・・,rn,・・・,rN]からサブアレイベクトルを構成する。
【０１２１】
具体的には、v番目のサブアレイは、v番目からv+W−1番目までのサブアレイの受信信号、即ち、rv= [rv,rv+1,・・・,rv+W-1]Tで構成される。Wは、各サブアレイを構成するサンプリング点数に相当する。よって、演算装置１２１１は、v番目のサブアレイから計算される相関行列Rcol(v)を以下の式（１３）のように計算する。下記の式(１３)における添字Hは、ベクトルの複素共役転置を表す。
【０１２２】
【数１３】
0006939981.tif 000014
【０１２３】
全てのサブアレイの相関行列Rcol(v) (v=1,2,・・・,V)の平均をRallとする。サブアレイの数Vは対象物の数K以上とする。
【０１２４】
上記の方法では、異なるサブアレイの受信信号間では相関が弱まるという性質を利用する事で、反射間の相関に起因する問題が回避される。更に、ステップＡ６では、演算装置１２１１は、以下の式(１４)で与えられる方向ベクトルa(θ)と、規格化ＩＦ信号から算出した相関行列Rallとを用いて、式(１５)から式(１７)で与えられる評価関数のいずれかを計算する。
【０１２５】
【数１４】
0006939981.tif 000015
【０１２６】
【数１５】
0006939981.tif 000016
【０１２７】
【数１６】
0006939981.tif 000017
【０１２８】
【数１７】
0006939981.tif 000018
【０１２９】
ここでベクトルｅｋ（ｋ＝Ｋ＋１，・・・，Ｗ）は、相関行列Ｒallの固有ベクトルの内、その固有値がノイズ電力に等しいものである。式(１４)内の位相φ(θ)は、上述の式(７)で与えられるものである。
【０１３０】
式（１５）から式（１７）で与えられる評価関数は、対象物１００３ｋ（ｋ＝１,2,・・・,Ｋ）が存在する角度θkにおいてピークを生じる。従って、演算装置１２１１は、式（１５）から式（１７）の評価関数のピークが得られる角度θｋから、対象物１００３ｋの位置を検知する。さらに、演算装置１２１１は、評価関数の値から、対象物１００３ｋ（ｋ＝１,2,・・・,Ｋ）における反射強度ｓｋの分布を見積もる事ができる。
【０１３１】
上記のステップＡ６までの処理により、対象物１００３ｋ（ｋ＝１,2,・・・,Ｋ）の位置を表す指標として距離Ｒ０と角度θｋとが定まる。
【０１３２】
［ステップＡ７］
次に、ステップＡ７について説明する。ステップＡ１〜ステップＡ６までのステップにおいて、距離を特定の値Ｒ０に限定した上で、対象物１００３ｋの存在角度θｋが、算出されている。そこで、ステップＡ７では、演算装置１２１１は、予め定められている距離全てについて処理を実施しているかどうかを判定する。
【０１３３】
判定の結果、距離全てについて処理を実施していない場合は、演算装置１２１１は、ステップＡ４に戻り、限定する距離を別の値に変更した上で、対象物１００３ｋの存在角度θｋを算出する。このように、ステップＡ４からステップＡ７までを繰り返す事で、演算装置１２１１は、予め定められている距離全てについて、対象物１００３ｋの存在角度θｋを算出する。
【０１３４】
一方、判定の結果、距離全てについて処理を実施している場合は、演算装置１２１１は、ステップＡ４からステップＡ６の処理を完了する。また、物体検知装置１０００は、予め定められている全ての距離と角度とで対象物１００３ｋの分布をスキャンできた時点で、その動作を終了する。
【０１３５】
［プログラム］
本実施の形態１におけるプログラムは、コンピュータに、図７に示すステップＡ１〜Ａ７を実行させるプログラムであれば良い。このプログラムを物体検知装置に備えられているコンピュータにインストールし、実行することによって、コンピュータのプロセッサが演算装置１２１１として機能する。この結果、本実施の形態１における物体検知装置１０００及び物体検知方法を実現することができる。
【０１３６】
［装置性能］
次に、本実施の形態における物体検知装置１０００の性能、詳しくは、対象物１００３ｋの反射強度ｓｋの分解能について説明する。
【０１３７】
上述の式（１５）から式（１７）の評価関数から得られる反射強度ｓｋの分解能は、本実施形態１における図９で示した受信アンテナ配置を持つアンテナアレイで形成されるビーム幅で定まる。例えば、式（１５）および（１６）に対応するビームステアリング型の方式であれば、反射強度ｓｋの分解能は、アンテナアレイのビーム幅と同程度となる。また式（１７）に対応するヌルステアリング型の方式であれば、反射強度ｓｋの分解能は、アンテナアレイのビーム幅の数分の一程度となる。
【０１３８】
本実施形態１においては、図９で示した受信アンテナ配置を持つアンテナアレイで形成されるビーム幅が算出される。
【０１３９】
また、アンテナアレイの受信信号が、上述の式(６)で与えられる時、ビームを角度θ０方向に向けた時のアンテナアレイ指向性(アレイファクターＡＦ（θ）)は、以下の式（１８）で与えられる。
【０１４０】
【数１８】
0006939981.tif 000019
【０１４１】
ここで、ｕ＝２πΔｑ（ｓｉｎθ−ｓｉｎθ０）である。なお、Δｑ＝αΔｔｍΔｄｍ/ｃ （ｍ＝１，２、・・・、Ｍ）である。Δｑはアンテナ番号ｍによらず一定になるように、アンテナ位置Δｄｍに応じて各受信アンテナに対応するサンプリング周期Δｔｍを設定する。
【０１４２】
上記式（１８）で与えられるアレイファクター｜ＡＦ（θ）｜は角度θ０において最大値を持ち、かつそのピーク幅すなわちビーム幅ΔθはおよそΔθ〜１/ＮΔｑ（ｒａｄ）で与えられる。
【０１４３】
即ち、ビーム幅Δθは、全てのアンテナのサンプリング点数の総和Ｎと、サンプリング周期とアンテナ位置とで決まるパラメータΔｑに依存する。なお、−πから＋πの範囲で虚像を発生させない条件である式（１０）の条件は、Δｑ？≦１/２と等価である。−πから＋πの範囲で虚像を発生させない条件下でビーム幅（角度分解能）を最良にするのはΔｑ＝１/２であり、その場合のビーム幅（角度分解能）はΔθ〜２/Ｎ（ｒａｄ）となる。
【０１４４】
ビーム幅Δθ、即ち、分解能は、アンテナ数Ｍに直接依存しない。即ち、本実施の形態１では、原理的には、一本の基準受信アンテナ１２０３０と、一本の測定受信アンテナ１２０３1と、の２本の受信アンテナのみで構成が可能である。この事は、多数のアンテナを用いる従来方式に比べて、大幅にアンテナ数を削減できる事を意味する。更に、本実施の形態１では、サンプリング点数の総和Ｎを増やすだけで、アンテナ数を増やさなくても任意に分解能の向上が可能である。
【０１４５】
また、ビーム幅Δθ、即ち、分解能は、ＲＦ信号の帯域幅ＢＷに直接依存しない。即ち、本実施の形態は、原理的には、ＲＦ信号の帯域幅を任意に狭く設定しても動作が可能である。この事は、電波法で使用できるＲＦ信号の帯域幅が制約されている状況においても、本実施の形態１を実現可能である事を意味する。更に、本実施の形態では、サンプリング点数の総和Ｎを増やすだけで、ＲＦ信号の帯域幅ＢＷを拡大しなくても、任意に分解能の向上が可能である。
【０１４６】
本実施の形態１では、ＦＭＣＷ信号の使用が想定されている。この時、ＩＦ信号の周波数は、４ＢＷ・Ｒ０/ｃＴchirpで与えられ、実用的な状況において、ＩＦ周波数は数十ｋＨｚから数百ｋＨｚ程度となる。ベースバンド信号の帯域幅が、数百ＭＨｚから数ＧＨｚと広帯域となるパルス方式に比べ、本実施の形態では、サンプリングレートと演算量とが削減され、且つ、雑音による信号品質劣化の問題が回避される。
【０１４７】
［実施の形態１における効果］
以上のように、本実施の形態における物体検知装置１０００及び物体検知方法によれば、一般的なアレイアンテナ方式と比較した場合、ＩＦ信号の時間波形のサンプリングデータを仮想的なアンテナとして用いることができる。このため、本実施の形態１によれば、仮想的なアンテナを、一般的なアレイアンテナ方式における受信アンテナと置き換えることができるので、受信アンテナの数を削減しながら、高分解能を実現できる。
【０１４８】
また、ＲＦ周波数の異なる複数のパルス信号を送信し、ＲＦ周波数に応じた位置に仮想的なアンテナを形成する従来方式と比較した場合を検討する。従来方式では、高分解能を実現するためにＲＦ周波数の範囲を広く設定する必要があるのに対し、本実施の形態１では、ＩＦ信号の時間波形のサンプリングデータを仮想的なアンテナとする事で、電波法の制約でＲＦ信号の帯域幅が取れない状況でも高分解能の動作を実現できる。更に、パルス信号を用いる従来方式と比較して、本実施の形態１では、ＩＦ信号が狭帯域であるため、サンプリングレートと演算量は少量で良く、また雑音による信号品質劣化の問題も回避できる。
【０１４９】
また、従来からの合成開口レーダー方式には、受信機を機械的に動かす必要があり、この結果、物体の検知及び検査のための時間が長くなるという問題があった。一方、本実施の形態１では、受信機の位置ではなく、受信周波数を電子的に走査すればよいので、合成開口レーダー方式に比べて、物体の検知及び検査のための時間を短縮化できる。
【０１５０】
即ち、本実施の形態１における物体検知装置及び物体検知方法においては、一般的なアレイアンテナ方式よりも、必要なアンテナ及びそれに付随する受信機の数を削減する事ができるので、装置のコスト、サイズ、重量を削減できるという効果を奏する。更に、本実施の形態１は、電波法の制約でＲＦ信号の帯域幅が取れない状況でも、高分解能の動作を実現できるという効果を奏する。
【０１５１】
更に、本実施の形態１における物体検知装置及び物体検知方法は、ＩＦ信号のサンプリングレートと演算量とを少量に抑え、雑音による信号品質劣化の問題を回避できるという効果も奏する。また、本実施の形態１における物体検知装置及び物体検知方法は、一般的な合成開口レーダー方式と異なり、装置を機械的に動かす必要がないため、物体検知及び検査の時間を短縮できるという効果も奏する。
【０１５２】
（実施の形態２）
次に、本発明の実施の形態２における、物体検知装置、物体検知方法、及びプログラムについて、図１０を用いて説明する。
【０１５３】
上述した実施の形態１においては、距離Ｒ０と角度θとの２つの変数により、対象物１００３の２次元的な位置が測定されている。これに対し、本実施の形態２では、対象物１００３の３次元的な位置が測定される。以下に具体的に説明する。
【０１５４】
［装置構成］
本実施の形態２における物体検知装置は、実施の形態１において図１から図５に示した物体検知装置と同様の構成を備えている。従って、本実施の形態２における物体検知装置の構成についての説明は、省略する。
但し、本実施の形態２においては、演算装置１２１１での処理において、実施の形態１と異なる部分がある。よって、以下に、実施の形態１との相違点について説明する。
【０１５５】
まず、図１０を用いて、本実施の形態２における、対象物１００３と、送受信共用アンテナ１２１０及び受信アンテナ１２０３との位置関係について説明する。図１０は、本発明の実施の形態２において検知の対象となる物体と検知を行うための送信アンテナ及び受信アンテナとの位置関係の一例を示す図である。
【０１５６】
図１０に示すように、本実施の形態２では、送受信共用アンテナ１２１０は、原点に設置される。また、ｘ軸上に受信アンテナ１２０３ｘ１、１２０３ｘ２、・・・、１２０３ｘＭｘがそれぞれ設置され、ｙ軸上に受信アンテナ１２０３ｙ１、１２０３ｙ２、・・・、１２０３ｙＭｙがそれぞれ設置される。
【０１５７】
このように、実施の形態１では、ｘ軸上のみに受信アンテナ１２０３が配置されていたのに対し、本実施の形態２では、ｘ軸上に加えて、ｙ軸上にも受信アンテナ１２０３が配置されている。本実施の形態２は、受信アンテナの配置の点で、実施の形態１と異なっている。
【０１５８】
また、図１０の例において、送受信共用アンテナ１２１０とｘ軸上に配置された受信アンテナ１２０３ｘｍ（ｍ＝０、１、２、・・・、Ｍｘ）との間隔をΔｄｘｍとする。更に、送受信共用アンテナ１２１０とｙ軸上に配置された受信アンテナ１２０３ｙｍ（ｍ＝０、１、２、・・・、Ｍｙ）との間隔をΔｄｙｍとする。
【０１５９】
そして、図１０の例では、対象物１００３ｋ（ｋ＝１、２、・・・、Ｋ）は、送受信共用アンテナ１２１０を中心とする半径Ｒ０の円上に配置されている状況を考える。また、対象物１００３ｋ（ｋ＝１、２、・・・、Ｋ）からの反射波は、各受信アンテナ１２０３において振幅skで受信されるものとする。また、送受信共用アンテナ１２１０と対象物１００３ｋとを結ぶベクトルとを想定する。この想定したベクトルとｘ軸との成す角をθｘｋとし、同ベクトルとｙ軸との成す角をθｙｋとする。
【０１６０】
［装置動作］
次に、本実施の形態２における物体検知装置の動作について説明する。また、本実施の形態２においても、物体検知装置を動作させることによって、物体検知方法が実施される。よって、本実施の形態２における物体検知方法の説明も、以下の物体検知装置の動作説明に代える。
【０１６１】
また、本実施の形態２においても、実施の形態１において図８に示したステップＡ１〜Ａ７が実行される。但し、本実施の形態２におけるステップＡ３〜Ａ７の詳細は、実施の形態１と異なっている。よって、以下に、本実施の形態におけるステップＡ３〜Ａ７の詳細について説明する。
【０１６２】
［ステップＡ３］
まず、各受信アンテナ１２０３で受信した受信信号に、送信部１１０１で生成された送信信号をそれぞれミキシングしてＩＦ信号を生成する、ステップＡ３の詳細について説明する。
【０１６３】
図１０に示した配置において、送受信共用アンテナ１２１０の受信信号から算出される基準ＩＦ信号（ＩＦ０（ｔ））は、実施の形態１で示した上述の式（３）で与えられる。
【０１６４】
ｘ軸に配置された測定受信アンテナ１２０３ｘｍ（ｍ＝１、２、・・・、Ｍｘ）の受信信号から算出される測定ＩＦ信号（ＩＦｘｍ（ｔ））は、以下の式（１９）で与えられる。
【０１６５】
【数１９】
0006939981.tif 000020
【０１６６】
またｙ軸に配置された測定受信アンテナ１２０３ｙｍ（ｍ＝１、２、・・・、Ｍｙ）の受信信号から算出される測定ＩＦ信号（ＩＦｙｍ（ｔ））は、以下の式（２０）で与えられる。
【０１６７】
【数２０】
0006939981.tif 000021
【０１６８】
［ステップＡ４］
次に、演算装置１２１１が、ステップＡ３で各受信アンテナ１２０３の受信信号から得たＩＦ信号に基づいて、対象物１００３の物体検知装置までの距離を特定の値に制約したＩＦ信号を生成する、ステップＡ４の詳細について説明する。
【０１６９】
図１０に示した例では、送受信共用アンテナ１２１０からの距離がＲ０となる球面状の対象物配置面１００５上に、対象物１００３が配置されている。実際の状況では、距離Ｒ０で指定される対象物配置面１００５上以外の位置に対象物が存在する事もある。そこでステップＡ４では、演算装置１２１１は、距離Ｒ０で指定される対象物配置面１００５上の位置に存在する対象物１００３からの信号のみを抽出する。
【０１７０】
例えば、対象物配置面１００５において、距離Ｒ０の位置に、対象物１００３が存在しているとする。この場合、対象物１００３により生成される式（３）の基準ＩＦ信号（ＩＦ０（ｔ））は、２αＲ０/cの周波数を有することになる。従って、基準受信アンテナ１２０３０又は送受信共用アンテナ１２１０で得た受信ＩＦ信号の時間波形をフーリエ変換し、２αＲ０/cの周波数の信号成分のみを抽出すれば、対象物１００３の位置を所望の位置（距離範囲）に制約したＩＦ信号が得られる。
【０１７１】
ｘ軸上に配置した受信アンテナ１２０３ｘｍ（ｍ＝１、２、・・・、Ｍｘ）で得られる式（１９）の測定ＩＦ信号（ＩＦｘｍ（ｔ））は、２α（Ｒ０−Δｄｘｍ）/cから２α（Ｒ０＋Δｄｘｍ）/cまでの周波数を有する。従って、測定受信アンテナ１２０３ｍ（ｍ＝１、２、・・・、Ｍ）で得た測定ＩＦ信号の時間波形をフーリエ変換し、２α（Ｒ０−Δｄｘｍ）/cから２α（Ｒ０＋Δｄｘｍ）/cまでの周波数の信号成分のみを抽出すれば、距離の値がＲ０に制約されたＩＦ信号が得られる。
【０１７２】
ｙ軸上に配置した受信アンテナ１２０３ｙｍ（ｍ＝１、２、・・・、Ｍｙ）で得られる式（２０）の測定ＩＦ信号（ＩＦｙｍ（ｔ））も同様である。フーリエ変換した後に２α（Ｒ０−Δｄｙｍ）/cから２α（Ｒ０＋Δｄｙｍ）/cまでの周波数の信号成分のみを抽出すれば、所望の距離範囲に制約したＩＦ信号が得られる。
【０１７３】
［ステップＡ５］
次に、測定受信アンテナＩＦ信号を、基準受信アンテナＩＦ信号によって規格化して、規格化ＩＦ信号を生成する、ステップＡ５の詳細について説明する。
【０１７４】
本実施の形態２では、ｘ軸上に配置した測定受信アンテナ１２０３ｘｍ（ｍ＝１、２、・・・、Ｍｘ）とｙ軸上に配置した測定受信アンテナ１２０３ｙｍ（ｍ＝１、２、・・・、Ｍｙ）とのそれぞれにおいて、実施の形態１と同じく、図９で示した内容に沿って受信アンテナ１２０３の位置と時間サンプリング点とが設定される。
【０１７５】
更に、本実施の形態２では、ｘ軸上に配置された測定受信アンテナ１２０３ｘｍで得られる測定ＩＦ信号（ＩＦｘｍ（ｔ））を、基準ＩＦ信号（ＩＦ０（ｔ））で除算して、規格化ＩＦ信号ｒｘｍ（ｔ）が生成される。規格化ＩＦ信号ｒｘｍ（ｔ）は、以下の式（２１）で与えられる。
【０１７６】
【数２１】
0006939981.tif 000022
【０１７７】
同様に、本実施の形態２では、ｙ軸上に配置された測定受信アンテナ１２０３ｙｍで得られる測定ＩＦ信号（ＩＦｙｍ（ｔ））を、基準ＩＦ信号（ＩＦ０（ｔ））で除算して、規格化ＩＦ信号ｒｙｍ（ｔ）が生成される。規格化ＩＦ信号ｒｙｍ（ｔ）は、以下の式（２２）で与えられる。
【０１７８】
【数２２】
0006939981.tif 000023
【０１７９】
x軸上に配置された受信アンテナ１２０３xm（ｍ＝１、２、・・・、Mx）で得られる規格化ＩＦ信号ｒxm（ｔ）を、サンプリング時刻ｔ（xm,１）, ｔ（xm,2）,・・・、ｔ（xm,px(m)）において取得する。px(m)は、受信アンテナ１２０３xmで得られる規格化ＩＦ信号のサンプリング点数である。全てのアンテナのサンプリング点数の総和px(1)+px(2)+・・・+px(Mx)をNxとする。
【０１８０】
同様に、y軸上に配置された受信アンテナ１２０３ym（ｍ＝１、２、・・・、My）で得られる規格化ＩＦ信号ｒym（ｔ）を、サンプリング時刻ｔ(ym,1), ｔ(ym,2),・・・、ｔ（ym,py(m)）において取得する。py(m)は、受信アンテナ１２０３ymで得られる規格化ＩＦ信号のサンプリング点数である。全てのアンテナのサンプリング点数の総和py(1)+py(2)+・・・+py(My)をNyとする。
【０１８１】
x軸上に配置された受信アンテナ１２０３xm（ｍ＝１、２、・・・、Mx）に関する規格化ＩＦ信号ベクトルrxは、以下の式（２３）のように定義される。式（２３）において、rxm（ｔ（ｍ,n））= rxm,nとしている。rxはＮx次ベクトルである。
【０１８２】
【数２３】
0006939981.tif 000024
【０１８３】
同様に、y軸上に配置された受信アンテナ１２０３ym（ｍ＝１、２、・・・、My）に関する規格化ＩＦ信号ベクトルryは、以下の式（２４）のように定義される。式（２４）においても、rym（ｔ（ｍ,n））= rym,nとしている。ryはＮy次ベクトルである。
【０１８４】
【数２４】
0006939981.tif 000025
【０１８５】
なお、本実施の形態では、x軸上に配置された受信アンテナ１２０３xm（ｍ＝１、２、・・・、Mx）で得られるサンプリング点数の総和Ｎｘと、ｙ軸上に配置された受信アンテナ１２０３ym（ｍ＝１、２、・・・、My）で得られるサンプリング点数の総和Ｎyとは、等しい事が好ましい以下では、Ｎ＝Ｎｘ＝Ｎｙとする。
【０１８６】
受信アンテナ１２０３xm及び受信アンテナ１２０３ymで得られる、規格化ＩＦ信号ｒxm（ｔ）及びｒym（ｔ）のサンプリング周期を、それぞれΔtxmおよびΔtymとする。虚像の発生を抑制するため、実施の形態１と同様に、サンプリング周期に対して、以下の式（２５）に示す条件が設定される。
【０１８７】
【数２５】
0006939981.tif 000026
【０１８８】
また、隣接する受信アンテナで位相の範囲が連続となるように、実施の形態１と同じく、受信アンテナの位置に以下の式（２６）及び（２７）に示す条件が設定される。
【０１８９】
【数２６】
0006939981.tif 000027
【０１９０】
【数２７】
0006939981.tif 000028
【０１９１】
ステップＡ５では、演算装置１２１１は、まずｘ軸上に配置された測定受信アンテナ１２０３xmで得られた、測定ＩＦ信号（ＩＦxm（ｔ））と基準ＩＦ信号（ＩＦ０（ｔ））とのそれぞれに対して、補間処理を行う。そして、演算装置１２１１は、上述の式（２５）を満たすサンプリング周期Δtxmのサンプリング時刻における値を算出する。次に、演算装置１２１１は、補間処理により測定ＩＦ信号（ＩＦxm（ｔ））と基準ＩＦ信号（ＩＦ０（ｔ））とのサンプリング時刻を揃えた後に、上述の式(２１)で示すように測定ＩＦ信号（ＩＦxm（ｔ））を基準ＩＦ信号（ＩＦ０（ｔ））で除算して規格化ＩＦ信号ｒxm（ｔ）(ｍ＝１、２、・・・、Ｍx)を算出する。
【０１９２】
同様に、演算装置１２１１は、y軸上に配置された測定受信アンテナ１２０３ymで得た測定ＩＦ信号（ＩＦym（ｔ））と基準ＩＦ信号（ＩＦ０（ｔ））とのそれぞれに対しても、補間処理を行う。そして、演算装置１２１１は、上述の式（２５）を満たすサンプリング周期Δtymのサンプリング時刻における値を算出する。次に、演算装置１２１１は、補間処理により測定ＩＦ信号（ＩＦym（ｔ））と基準ＩＦ信号（ＩＦ０（ｔ））とのサンプリング時刻を揃えた後に、上述の式(２２)で示すように測定ＩＦ信号（ＩＦym（ｔ））を基準ＩＦ信号（ＩＦ０（ｔ））で除算して規格化ＩＦ信号ｒym（ｔ）(ｍ＝１、２、・・・、Ｍy)を算出する。
【０１９３】
［ステップＡ６］
次に、規格化ＩＦ信号から評価関数を算出するステップＡ６の詳細について説明する。
【０１９４】
まず、ステップＡ６では、演算装置１２１１は、x軸上に配置した受信アンテナ１２０３xm及びy軸上に配置した受信アンテナ１２０３ymで得られる、規格化ＩＦ信号ベクトルrx=[rx1,rx2,・・・, rxn,・・・,rxN]T及びry=[ry1,ry2,・・・, ryn,・・・,ryN]Tから、それぞれサブアレイベクトルを構成する。具体的には、v番目のサブアレイは、v番目からv+W−1番目までのサブアレイの受信信号、即ちrxv = [rxv,rx(v+1),・・・,rx(v+W-1)]Tとryv= [ryv,ry(v+1),・・・,ry(v+W-1)]Tとで構成される。Wは、各サブアレイを構成するサンプリング点数に相当する。
【０１９５】
次に、演算装置１２１１は、サブアレイベクトルrxvとryvの直積を取り、以下の式（２８）に示す規格化IF信号rxyvを生成する。
【０１９６】
【数２８】
0006939981.tif 000029
【０１９７】
規格化IF信号rxyvは、rxvとryvとの要素の全ての組み合わせ毎の積を要素とした（N2×1）次ベクトルとなる。
【０１９８】
次に、演算装置１２１１は、v番目のサブアレイから計算される相関行列Rcol(v)を以下の式（２９）のように計算する。
【０１９９】
【数２９】
0006939981.tif 000030
【０２００】
全てのサブアレイの相関行列Rcol(v) (v=1,2,・・・,V)の平均をRallとする。サブアレイの数Vは対象物の数K以上とする。
【０２０１】
上記の方法では、異なるサブアレイの受信信号間では相関が弱まるという性質を利用する事で、反射間の相関に起因する問題が回避される。
【０２０２】
次に、演算装置１２１１は、以下の式(３０)で与えられる方向ベクトルa(θx,θy)を定義する。
【０２０３】
【数３０】
0006939981.tif 000031
【０２０４】
また、上記式（３０）における、時刻txn及びtynは、それぞれ規格化ＩＦ信号ベクトルrxおよびryのn番目の成分に対応するサンプリング時刻である。またアンテナ位置Δdmx及びΔdmyは、nの値によって変化する。
【０２０５】
次に、式(３１)から式(３３)で与えられる評価関数のいずれかを計算する。
【０２０６】
【数３１】
0006939981.tif 000032
【０２０７】
【数３２】
0006939981.tif 000033
【０２０８】
【数３３】
0006939981.tif 000034
【０２０９】
上記式（３１）から式（３３）で与えられる評価関数は、対象物１００３ｋ（ｋ＝１,2,・・・,Ｋ）が存在する角度(θxk,θyk）においてピークを生じる。従って式（３１）から式（３３）の評価関数のピークが得られる角度(θxk,θyk）から、対象物１００３ｋの位置を検知できる。更に、評価関数の値から、対象物１００３ｋ（ｋ＝１,2,・・・,Ｋ）における反射強度ｓｋの分布を見積もる事ができる。
【０２１０】
ステップＡ１〜ステップＡ６までの処理で、対象物１００３ｋの距離Ｒ０と角度(θxk,θyk）とが測定されるので、それらのデータから対象物１００３ｋの３次元的な位置を算出できる。
【０２１１】
［ステップＡ７］
次に、ステップＡ７について説明する。ステップＡ１〜ステップＡ６までのステップにおいて、距離を特定の値Ｒ０に限定した上で、対象物１００３ｋの存在角度(θxk,θyk）が、算出されている。そこで、ステップＡ７では、演算装置１２１１は、予め定められている距離全てについて処理を実施しているかどうかを判定する。
【０２１２】
判定の結果、距離全てについて処理を実施していない場合は、演算装置１２１１は、ステップＡ４に戻り、限定する距離を別の値に変更した上で、対象物１００３ｋの存在角度(θxk,θyk）を算出する。このように、ステップＡ４からステップＡ７までを繰り返す事で、演算装置１２１１は、予め定められている距離全てについて、対象物１００３ｋの存在角度(θxk,θyk）を算出する。
【０２１３】
一方、判定の結果、距離全てについて処理を実施している場合は、演算装置１２１１は、ステップＡ４からステップＡ６の処理を完了する。また、物体検知装置１０００は、予め定められている全ての距離と角度とで対象物１００３ｋの分布をスキャンできた時点で、その動作を終了する。
【０２１４】
［プログラム］
本実施の形態２におけるプログラムも、コンピュータに、図７に示すステップＡ１〜Ａ７を実行させるプログラムであれば良い。このプログラムを物体検知装置に備えられているコンピュータにインストールし、実行することによって、コンピュータのプロセッサが演算装置１２１１として機能する。この結果、本実施の形態２における物体検知装置１０００及び物体検知方法を実現することができる。
【０２１５】
［装置性能及び実施の形態２における効果］
本実施の形態２においても、実施の形態１において述べた性能及び効果を得ることができる。また、本実施の形態２によれば、対象物１００３の３次元空間における位置を測定でき、実施の形態１よりも多くの場面に適用可能となる。
【０２１６】
（物理構成）
ここで、実施の形態１及び２におけるプログラムを実行することによって、物体検知装置を実現するコンピュータ（演算装置）について図１１を用いて説明する。図１１は、本発明の実施の形態１及び２における物体検知装置を実現するコンピュータの一例を示すブロック図である。
【０２１７】
図１１に示すように、コンピュータ１１０は、ＣＰＵ（Central Processing Unit）１１１と、メインメモリ１１２と、記憶装置１１３と、入力インターフェイス１１４と、表示コントローラ１１５と、データリーダ／ライタ１１６と、通信インターフェイス１１７とを備える。これらの各部は、バス１２１を介して、互いにデータ通信可能に接続される。なお、コンピュータ１１０は、ＣＰＵ１１１に加えて、又はＣＰＵ１１１に代えて、ＧＰＵ（Graphics Processing Unit）、又はＦＰＧＡ（Field-Programmable Gate Array）を備えていても良い。
【０２１８】
ＣＰＵ１１１は、記憶装置１１３に格納された、本実施の形態におけるプログラム（コード）をメインメモリ１１２に展開し、これらを所定順序で実行することにより、各種の演算を実施する。メインメモリ１１２は、典型的には、ＤＲＡＭ（Dynamic Random Access Memory）等の揮発性の記憶装置である。また、本実施の形態におけるプログラムは、コンピュータ読み取り可能な記録媒体１２０に格納された状態で提供される。なお、本実施の形態におけるプログラムは、通信インターフェイス１１７を介して接続されたインターネット上で流通するものであっても良い。
【０２１９】
また、記憶装置１１３の具体例としては、ハードディスクドライブの他、フラッシュメモリ等の半導体記憶装置が挙げられる。入力インターフェイス１１４は、ＣＰＵ１１１と、キーボード及びマウスといった入力機器１１８との間のデータ伝送を仲介する。表示コントローラ１１５は、ディスプレイ装置１１９と接続され、ディスプレイ装置１１９での表示を制御する。なお、コンピュータ１１０は、ＣＰＵ１１１に加えて、又はＣＰＵ１１１に代えて、ＧＰＵ（Graphics Processing Unit）、又はＦＰＧＡ（Field-Programmable Gate Array）を備えていても良い。
【０２２０】
データリーダ／ライタ１１６は、ＣＰＵ１１１と記録媒体１２０との間のデータ伝送を仲介し、記録媒体１２０からのプログラムの読み出し、及びコンピュータ１１０における処理結果の記録媒体１２０への書き込みを実行する。通信インターフェイス１１７は、ＣＰＵ１１１と、他のコンピュータとの間のデータ伝送を仲介する。
【０２２１】
また、記録媒体１２０の具体例としては、ＣＦ（Compact Flash（登録商標））及びＳＤ（Secure Digital）等の汎用的な半導体記憶デバイス、フレキシブルディスク（Flexible Disk）等の磁気記録媒体、又はＣＤ−ＲＯＭ（Compact DiskRead Only Memory）などの光学記録媒体が挙げられる。
【０２２２】
なお、本実施の形態における物体検知装置は、プログラムがインストールされたコンピュータではなく、各部に対応したハードウェアを用いることによっても実現可能である。更に、物体検知装置は、一部がプログラムで実現され、残りの部分がハードウェアで実現されていてもよい。
【０２２３】
（実施の形態１及び２における効果の要約）
以下において、本実施の形態の効果を要約する。
【０２２４】
特許文献１から特許文献３に記載された一般的なアレイアンテナ方式と、本実施の形態とを比較した場合、前者のアレイアンテナ方式は多数のアンテナを必要とする。その一方で、本実施の形態では、ＩＦ信号の時間波形のサンプリングデータを仮想的なアンテナとして実際の受信アンテナと置き換える事で、受信アンテナの数を削減しながら高分解能を実現できる。
【０２２５】
ＲＦ周波数の異なる複数のパルス信号を送信し、ＲＦ周波数に応じた位置に仮想的なアンテナを形成する特許文献４に記載の従来方式では、高分解能を実現するためにＲＦ周波数の範囲を広く設定しなければならない。しかしながら、この従来方式には、電波法で使用できるＲＦ周波数の範囲が制約されているため実現が困難という問題点があった。
【０２２６】
一方で、本実施の形態では、ＲＦ周波数に応じた位置に仮想的なアンテナを形成する代わりに、ＩＦ信号の時間波形のサンプリングデータを仮想的なアンテナとする事で、電波法の制約でＲＦ信号の帯域幅が取れない状況でも高分解能の動作を実現できる。
【０２２７】
更に、パルス方式に基づく特許文献４に記載の従来方式では、ベースバンド信号の帯域幅が数百ＭＨｚから数ＧＨｚと広帯域であるため、高速なサンプリングレートと多大な演算量が必要となる事、及び雑音による信号品質劣化の問題があった。一方で、本実施の形態では、ＦＭＣＷ方式に基づくため、ＩＦ信号の帯域幅は、数十ｋＨｚから数百ｋＨｚと狭く、サンプリングレートと演算量とは少量でよく、また雑音による信号品質劣化の問題も回避できる。
【０２２８】
合成開口レーダー方式と本実施の形態とを比較した場合、合成開口レーダー方式は受信機３０１を機械的に動かす必要があり（図１４参照）、このため、物体の検知及び検査のための時間が長くなるという問題があった。一方、本実施の形態では、受信機の位置ではなくＩＦ信号の時間波形のサンプリングデータを電子的に走査すればよいので、合成開口レーダー方式に比べて、物体の検知及び検査のための時間を短縮できる。更に、本実施の形態では、アンテナを機械的に動かす装置が不要であるので、装置を小型化でき、電子走査で機構が摩滅する事はないので、機械走査に比べて装置寿命とメンテナンスコストとを改善できるという利点がある。
【０２２９】
即ち、本実施の形態における物体検知装置及び物体検知方法においては、一般的なアレイアンテナ方式よりも、必要なアンテナの数、及びそれに付随する受信機の数を削減する事ができるので、装置のコスト、サイズ、重量を削減できるという効果を奏する。更に、電波法の制約でＲＦ信号の帯域幅が取れない状況でも高分解能の動作を実現できるという効果を奏する。
【０２３０】
更に、ＩＦ信号のサンプリングレートと演算量とを少量に抑え、雑音による信号品質劣化の問題も回避できるという効果を奏する。また、本実施の形態における物体検知装置及び物体検知方法においては、一般的な合成開口レーダー方式と異なり、装置を機械的に動かす必要がないため、物体検知及び検査の時間を短縮できるという効果を奏する。
【０２３１】
上述した実施の形態の一部又は全部は、以下に記載する（付記１）〜（付記２４）によって表現することができるが、以下の記載に限定されるものではない。
【０２３２】
（付記１）
電波によって物体を検知するための物体検知装置であって、
前記物体に向けて、送信信号となる電波を照射する、送信部と、
複数の受信アンテナによって、前記物体で反射された電波を受信信号として受信し、更に、複数の前記受信アンテナそれぞれで受信された前記受信信号毎に、当該受信信号を用いて、中間周波数信号を生成する、受信部と、
前記中間周波数信号それぞれを合成して得られるビームパターンによる虚像の発生が抑制されるように、サンプリング時刻を決定し、決定したサンプリング時刻で前記中間周波数信号に対してサンプリングを行うことによって、物体の位置検出用の中間周波数信号を生成し、
生成した前記位置検出用の中間周波数信号を用いて、前記物体を検知する、演算部と、
を備えている、
ことを特徴とする物体検知装置。
【０２３３】
（付記２）
付記１に記載の物体検知装置であって、
前記受信部が、基準受信アンテナと測定受信アンテナとを備え、
前記演算部は、
前記測定受信アンテナで受信された受信信号から生成された中間周波数信号と、前記基準受信アンテナで受信された受信信号から生成された中間周波数信号とに対して、サンプリングを実行し、
得られたサンプリング値について、前記測定受信アンテナで受信された受信信号から生成された中間周波数信号を、前記基準受信アンテナで受信された受信信号から生成された中間周波数信号によって、規格化して、前記位置検出用の中間周波数信号として、規格化中間周波数信号を生成し、
更に前記規格化中間周波数信号から、前記物体の位置分布を表す評価関数を算出し、算出した前記評価関数を用いて、前記物体の位置及び形状を検知する、
ことを特徴とする物体検知装置。
【０２３４】
（付記３）
付記２に記載の物体検知装置であって、
前記演算部は、前記測定受信アンテナで受信された受信信号から生成された中間周波数信号と、前記基準受信アンテナで受信された受信信号から生成された中間周波数信号とに対し、検知される前記物体の位置を制限する処理を実行し、
前記処理が実行された前記中間周波数信号を用いて、前記サンプリングを実行する、
ことを特徴とする物体検知装置。
【０２３５】
（付記４）
付記２または３に記載の物体検知装置であって、
前記演算部は、前記規格化中間周波数信号から相関行列を算出し、算出した前記相関行列から、前記評価関数を算出する、
ことを特徴とする物体検知装置。
【０２３６】
（付記５）
付記４に記載の物体検知装置であって、
前記測定受信アンテナが、複数備えられ、且つ、前記基準受信アンテナの位置を基準にして、少なくとも二つの方向それぞれに沿って配置されており、
前記演算部は、前記方向毎に、当該方向に設置された前記測定受信アンテナから、前記規格化中間周波数信号を生成し、更に、前記方向毎に生成された前記規格化中間周波数信号の直積から、前記相関行列を算出し、そして、算出した前記相関行列から、前記物体の多次元の位置分布を表す前記評価関数を算出する、
ことを特徴とする物体検知装置。
【０２３７】
（付記６）
付記４または５に記載の物体検知装置であって、
前記演算部は、前記規格化中間周波数信号の範囲を複数選択し、選択した前記複数の範囲を定めた複数の前記規格化中間周波数信号それぞれから、前記相関行列を算出し、更に前記相関行列の平均から、前記物体の位置分布を表す前記評価関数を算出する、
ことを特徴とする物体検知装置。
【０２３８】
（付記７）
付記１〜６のいずれかに記載の物体検知装置であって、
前記送信部は、前記物体に向けて周波数が時間に応じて変化する電波を、前記送信信号として照射し、
前記受信部は、前記受信アンテナの各々で受信した前記受信信号に前記送信信号をミキシングして前記中間周波数信号を生成する、
ことを特徴とする物体検知装置。
【０２３９】
（付記８）
付記７に記載の物体検知装置であって、
前記受信アンテナが、複数備えられ、
隣接する前記受信アンテナの間隔が、前記送信信号の周波数の最小値と最大値とで定められる上限値以下に設定されている、
ことを特徴とする物体検知装置。
【０２４０】
（付記９）
電波によって物体を検知するための方法であって、
（ａ）前記物体に向けて、送信信号となる電波を照射する、ステップと、
（ｂ）複数の受信アンテナによって、前記物体で反射された電波を受信信号として受信し、更に、複数の前記受信アンテナそれぞれで受信された前記受信信号毎に、当該受信信号を用いて、中間周波数信号を生成する、ステップと、
（ｃ）前記中間周波数信号それぞれを合成して得られるビームパターンによる虚像の発生が抑制されるように、サンプリング時刻を決定し、決定したサンプリング時刻で前記中間周波数信号に対してサンプリングを行うことによって、物体の位置検出用の中間周波数信号を生成する、ステップと、
（ｄ）生成した前記位置検出用の中間周波数信号を用いて、前記物体を検知する、ステップと、を有している、
ことを特徴とする物体検知方法。
【０２４１】
（付記１０）
付記９に記載の物体検知方法であって、
前記（ｂ）のステップにおいて、前記受信アンテナとして、基準受信アンテナと測定受信アンテナとを用い、
前記（ｃ）のステップにおいて、
前記測定受信アンテナで受信された受信信号から生成された中間周波数信号と、前記基準受信アンテナで受信された受信信号から生成された中間周波数信号とに対して、サンプリングを実行し、
得られたサンプリング値について、前記測定受信アンテナで受信された受信信号から生成された中間周波数信号を、前記基準受信アンテナで受信された受信信号から生成された中間周波数信号によって、規格化して、前記位置検出用の中間周波数信号として、規格化中間周波数信号を生成し、
前記（ｄ）のステップにおいて、前記規格化中間周波数信号から、前記物体の位置分布を表す評価関数を算出し、算出した前記評価関数を用いて、前記物体の位置及び形状を検知する、
ことを特徴とする物体検知方法。
【０２４２】
（付記１１）
付記１０に記載の物体検知方法であって、
前記（ｃ）のステップにおいて、
前記測定受信アンテナで受信された受信信号から生成された中間周波数信号と、前記基準受信アンテナで受信された受信信号から生成された中間周波数信号とに対し、検知される前記物体の位置を制限する処理を実行し、
前記処理が実行された前記中間周波数信号を用いて、前記サンプリングを実行する、
ことを特徴とする物体検知方法。
【０２４３】
（付記１２）
付記１０または１１に記載の物体検知方法であって、
前記（ｄ）のステップにおいて、前記規格化中間周波数信号から相関行列を算出し、算出した前記相関行列から、前記評価関数を算出する、
ことを特徴とする物体検知方法。
【０２４４】
（付記１３）
付記１２に記載の物体検知方法であって、
前記測定受信アンテナが、複数備えられ、且つ、前記基準受信アンテナの位置を基準にして、少なくとも二つの方向それぞれに沿って配置されており、
前記（ｃ）のステップにおいて、前記方向毎に、当該方向に設置された前記測定受信アンテナから、前記規格化中間周波数信号を生成し、
前記（ｄ）のステップにおいて、前記方向毎に生成された前記規格化中間周波数信号の直積から、前記相関行列を算出し、そして、算出した前記相関行列から、前記物体の多次元の位置分布を表す前記評価関数を算出する、
ことを特徴とする物体検知方法。
【０２４５】
（付記１４）
付記１２または１３に記載の物体検知方法であって、
前記（ｄ）のステップにおいて、前記規格化中間周波数信号の範囲を複数選択し、選択した前記複数の範囲を定めた複数の前記規格化中間周波数信号それぞれから、前記相関行列を算出し、更に前記相関行列の平均から、前記物体の位置分布を表す前記評価関数を算出する、
ことを特徴とする物体検知方法。
【０２４６】
（付記１５）
付記９〜１４のいずれかに記載の物体検知方法であって、
前記（ａ）のステップにおいて、前記物体に向けて周波数が時間に応じて変化する電波を、前記送信信号として照射し、
前記（ｂ）のステップにおいて、前記受信アンテナの各々で受信した前記受信信号に前記送信信号をミキシングして前記中間周波数信号を生成する、
ことを特徴とする物体検知方法。
【０２４７】
（付記１６）
付記１５に記載の物体検知方法であって、
前記受信アンテナが、複数備えられ、
前記（ｂ）のステップにおいて、隣接する前記受信アンテナの間隔が、前記送信信号の周波数の最小値と最大値とで定められる上限値以下に設定されている、
ことを特徴とする物体検知方法。
【０２４８】
（付記１７）
物体に向けて、送信信号となる電波を照射する、送信部と、複数の受信アンテナによって、前記物体で反射された電波を受信信号として受信し、更に、複数の前記受信アンテナそれぞれで受信された前記受信信号毎に、当該受信信号を用いて、中間周波数信号を生成する、受信部と、コンピュータと、を備えた物体検知装置において、
前記コンピュータに、
（ａ）前記中間周波数信号それぞれを合成して得られるビームパターンによる虚像の発生が抑制されるように、サンプリング時刻を決定し、決定したサンプリング時刻で前記中間周波数信号に対してサンプリングを行うことによって、物体の位置検出用の中間周波数信号を生成する、ステップと、
（ｂ）生成した前記位置検出用の中間周波数信号を用いて、前記物体を検知する、ステップと、を実行させる命令を含む、
ことを特徴とするプログラム。
【０２４９】
（付記１８）
付記１７に記載のプログラムであって、
前記受信部が、前記受信アンテナとして、基準受信アンテナと測定受信アンテナとを備えており、
前記（ａ）のステップにおいて、
前記測定受信アンテナで受信された受信信号から生成された中間周波数信号と、前記基準受信アンテナで受信された受信信号から生成された中間周波数信号とに対して、サンプリングを実行し、
得られたサンプリング値について、前記測定受信アンテナで受信された受信信号から生成された中間周波数信号を、前記基準受信アンテナで受信された受信信号から生成された中間周波数信号によって、規格化して、前記位置検出用の中間周波数信号として、規格化中間周波数信号を生成し、
前記（ｂ）のステップにおいて、前記規格化中間周波数信号から、前記物体の位置分布を表す評価関数を算出し、算出した前記評価関数を用いて、前記物体の位置及び形状を検知する、
ことを特徴とするプログラム。
【０２５０】
（付記１９）
付記１８に記載のプログラムであって、
前記（ａ）のステップにおいて、
前記測定受信アンテナで受信された受信信号から生成された中間周波数信号と、前記基準受信アンテナで受信された受信信号から生成された中間周波数信号とに対し、検知される前記物体の位置を制限する処理を実行し、
前記処理が実行された前記中間周波数信号を用いて、前記サンプリングを実行する、
ことを特徴とするプログラム。
【０２５１】
（付記２０）
付記１８または１９に記載のプログラムであって、
前記（ｂ）のステップにおいて、前記規格化中間周波数信号から相関行列を算出し、算出した前記相関行列から、前記評価関数を算出する、
ことを特徴とするプログラム。
【０２５２】
（付記２１）
付記２０に記載のプログラムであって、
前記測定受信アンテナが、複数備えられ、且つ、前記基準受信アンテナの位置を基準にして、少なくとも二つの方向それぞれに沿って配置されており、
前記（ａ）のステップにおいて、前記方向毎に、当該方向に設置された前記測定受信アンテナから、前記規格化中間周波数信号を生成し、
前記（ｂ）のステップにおいて、前記方向毎に生成された前記規格化中間周波数信号の直積から、前記相関行列を算出し、そして、算出した前記相関行列から、前記物体の多次元の位置分布を表す前記評価関数を算出する、
ことを特徴とするプログラム。
【０２５３】
（付記２２）
付記２０または２１に記載のプログラムであって、
前記（ｂ）のステップにおいて、前記規格化中間周波数信号の範囲を複数選択し、選択した前記複数の範囲を定めた複数の前記規格化中間周波数信号それぞれから、前記相関行列を算出し、更に前記相関行列の平均から、前記物体の位置分布を表す前記評価関数を算出する、
ことを特徴とするプログラム。
【０２５４】
（付記２３）
付記１７〜２２のいずれかに記載のプログラムであって、
前記送信部が、前記物体に向けて周波数が時間に応じて変化する電波を、前記送信信号として照射し、
前記受信部が、前記受信アンテナの各々で受信した前記受信信号に前記送信信号をミキシングして前記中間周波数信号を生成する、
ことを特徴とするプログラム。
【０２５５】
（付記２４）
付記２３に記載のプログラムであって、
前記受信アンテナが、複数備えられ、
隣接する前記受信アンテナの間隔が、前記送信信号の周波数の最小値と最大値とで定められる上限値以下に設定されている、
ことを特徴とするプログラム。
【０２５６】
以上、本発明の好適な実施形態の構成を説明した。しかし、前述の各特許文献等に開示されている内容は、本発明に引用をもって繰り込むことも可能とする。本発明の全開示（特許請求の範囲を含む）の枠内において、さらにその基本的技術思想に基づいて、実施の形態の変更・調整が可能である。
【０２５７】
また、本発明の特許請求の範囲の枠内において種々の開示要素の多様な組み合わせあるいは選択も可能である。すなわち、本発明は、特許請求の範囲を含む全開示、技術的思想にしたがって、当業者であればなし得ることが可能な各種変形、修正を含むことは勿論である。
【０２５８】
更に、本願発明は上記実施の形態に限定されるものではない。本願発明の構成や詳細には、本願発明のスコープ内で当業者が理解し得る様々な変更をすることができる。
【産業上の利用可能性】
【０２５９】
以上のように、本発明によれば、電波による物体のイメージングにおいて、画質の低下を抑制しつつ、装置の低コスト化及び小型化を図ることができる。本発明は、レーダー装置、衣服の下及びカバンの中の物体を検査するイメージング装置等に有用である。
【符号の説明】
【０２６０】
１１０  コンピュータ
１１１  ＣＰＵ
１１２  メインメモリ
１１３  記憶装置
１１４  入力インターフェイス
１１５  表示コントローラ
１１６  データリーダ／ライタ
１１７  通信インターフェイス
１１８  入力機器
１１９  ディスプレイ装置
１２０  記録媒体
１２１  バス
【図１】
0006939981.tif 000035
【図２】
0006939981.tif 000036
【図３】
0006939981.tif 000037
【図４】
0006939981.tif 000038
【図５】
0006939981.tif 000039
【図６】
0006939981.tif 000040
【図７】
0006939981.tif 000041
【図８】
0006939981.tif 000042
【図９】
0006939981.tif 000043
【図１０】
0006939981.tif 000044
【図１１】
0006939981.tif 000045
【図１２】
0006939981.tif 000046
【図１３】
0006939981.tif 000047
【図１４】
0006939981.tif 000048
【図１５】
0006939981.tif 000049