(57)【要約】【課題】波源の帯域幅が狭い場合にも波源相互相関値を低減することができる受信波の相関行列算出装置を提供する。【解決手段】受信波の相関行列算出装置は、アレーアンテナ２４ａにより電波を受信して相関行列を算出する。受信波の相関行列算出装置は、受信波の方向行列の位相を回転し、位相回転された方向行列を含む相関行列を周波数平均して平均化された相関行列を算出する。【選択図】図１
2010266294.tif 000002
【特許請求の範囲】
【請求項１】
複数の受信素子により電波又は音波を受信して相関行列を算出する受信波の相関行列算出装置において、
受信波の方向行列の位相を回転する位相回転手段と、
位相回転された方向行列を含む相関行列を平均して平均化された相関行列を算出する平均手段とを備える
ことを特徴とする受信波の相関行列算出装置。
【請求項２】
請求項１に記載の受信波の相関行列算出装置において、
受信波の波源のおおよその方向を推定する初期方向推定手段と、
方向推定された受信波の方向行列を算出する方向行列算出手段とを備える
ことを特徴とする受信波の相関行列算出装置。
【請求項３】
請求項１に記載の受信波の相関行列算出装置において、
受信波の波源のおおよその位置を推定する初期位置推定手段と、
位置推定された受信波の方向行列を算出する方向行列算出手段とを備える
ことを特徴とする受信波の相関行列算出装置。
【請求項４】
複数の受信素子により電波を受信して相関行列を算出する受信波の相関行列算出方法において、
受信波の方向行列の位相を回転する位相回転と、
位相回転された方向行列を含む相関行列を平均して平均化された相関行列を算出する平均処理とを行う
ことを特徴とする受信波の相関行列算出方法。
【請求項５】
請求項１又は２に記載の受信波の相関行列算出装置が算出した相関行列に基づいて受信波の波源の方向を推定する波源方向推定手段を備える
ことを特徴とする波源方向推定装置。
【請求項６】
請求項１又は３に記載の受信波の相関行列算出装置が算出した相関行列に基づいて受信波の波源の位置を推定する波源位置推定手段を備える
ことを特徴とする波源位置推定装置。
【請求項７】
請求項２に記載の受信波の相関行列算出装置が算出した相関行列に基づいて前記波源の方向を推定する波源方向推定手段を備える波源方向推定装置であって、
前記波源方向推定手段は、前記初期方向推定手段による方向推定から前記相関行列の算出までの処理を複数回繰り返させ、算出された相関行列に基づいて前記波源の方向を推定する
ことを特徴とする波源方向推定装置。
【請求項８】
請求項３に記載の受信波の相関行列算出装置が算出した相関行列に基づいて前記波源の位置を推定する波源位置推定手段を備える波源位置推定装置であって、
前記波源位置推定手段は、前記初期位置推定手段による位置推定から前記相関行列の算出までの処理を複数回繰り返させ、算出された相関行列に基づいて前記波源の位置を推定する
ことを特徴とする波源位置推定装置。
【発明の詳細な説明】
【技術分野】
【０００１】
この発明は、受信波の相関行列算出装置及び受信波の相関行列算出方法、並びに波源方向推定装置及び波源位置推定装置に関する。
【背景技術】
【０００２】
近年、地上デジタル放送機器、携帯電話、電子キー等の様々な無線機器の使用が日常化してきており、希望の相手からの電波を受信するために伝播特性測定の必要性が高まっている。また、反射、回折、散乱により電波の伝播経路が複数となるパルチパスではそれぞれの経路の電波が互いに干渉して受信レベルが激しく変動する。
【０００３】
例えば、車両においては、無線通信により車両に固有のＩＤコードを発信可能な電子キーを車両のキーとして使用する電子キーシステムが搭載されている。この電子キーシステムでは、電子キーから発信されたＩＤコードを車両が受信すると、車両が電子キーのＩＤコードと自身に予め登録されているＩＤコードとを照合するＩＤ照合を行い、これらＩＤコードが一致してＩＤ照合が成立すれば、車載機器の制御、例えば車両のドアロックの施解錠やエンジンの始動が許可又は実行される。この電子キーシステムでは、車両室内において、マルチパスによりコヒーレントな複数波が同時に入射し、波源と受信アレーアンテナとの距離が短く近傍波源として取り扱う必要がある。また、電子キーの位置が車内と車外とのいずれにあるかで車載機器の制御を切り分けている場合には、車両からのＩＤコードの返信要求を車内と車外とで異なるタイミングで発信して対応している。そこで、車両が受信した電子キーからの電波の波源の位置を推定できれば、ＩＤコードの返信要求を異なるタイミングで発信する必要がなくなるため、波源位置推定が求められている。
【０００４】
波源位置推定を行う方法としては、ＭＵＳＩＣ（ＭＵｌｔｉｐｌｅ  ＳＩｇｎａｌ  Ｃｌａｓｓｉｆｉｃａｔｉｏｎ）法やＥＳＰＲＩＴ（Ｅｓｔｉｍａｔｉｏｎ  ｏｆ  Ｓｉｇｎａｌ  Ｐａｒａｍｅｔｅｒｓ  ｖｉｓａ  Ｒｏｔａｔｉｏｎａｌ  Ｉｎｖａｒｉａｎｃｅ  Ｔｅｃｈｎｉｑｕｅｓ）法等の固有値展開に基づく方法がある。しかしながら、車両室内のようなコヒーレント波の推定を行うためには、相関行列のランクを回復するために空間平均法、周波数平均法等の前処理を必要とする（例えば、非特許文献１参照）。
【０００５】
上記のＭＵＳＩＣ法では、コヒーレント波には元々対応していないため、コヒーレント波の位相関係は受信位置で異なるので、受信点を平行移動させて相関値を求める空間平均法を用いて波源相互相関値を低減していた。ところが、球面波を受信した場合には等間隔アレーであるにも関わらず、素子ごとの位相差が異なる仮想的な不等間隔アレーとして処理することとなるため、等間隔アレーを前提としている空間平均法をそのまま適用することができない。そこで、補間ＭＵＳＩＣ法という方法が提案されているが、平均化処理のためにアンテナ素子を多数用意する必要があった。
【先行技術文献】
【非特許文献】
【０００６】
【非特許文献１】菊間信良著「アレーアンテナによる適応信号処理」科学技術出版
【発明の概要】
【発明が解決しようとする課題】
【０００７】
一方、波源相互相関値を低減する方法として周波数平均法を用いることが考えられる。この周波数平均法においては、周波数範囲が広い場合に周波数ごとの方向ベクトルの差が大きくなるため、推定精度が劣化するという問題点があった。この問題を解決する方法として、ＣＳＳ（Ｃｏｈｅｒｅｎｔ  Ｓｉｇｎａｌ−Ｓｕｂｓｐａｃｅ）法が提案されている。ＣＳＳ法はアレー形状に拠らないため、近傍波源に対しても適用可能であり、空間平均法のように平均化処理のためにアレー全体のアンテナ素子数を増やす必要がなく、ＣＳＳ法は空間平均法と比べ、アンテナ素子数が削減できることとなり、特に平面アレーにおいては素子数削減効果が大きいことが特徴である。しかしながら、実際に装置に適用するにあたっては電波法によって広い帯域の電波を発信することは困難であり、狭い帯域幅の電波をやりとりすることとなる。発明者らは、ＣＳＳ法について検討を行い、波源の帯域幅が狭い場合には波源相互相関値が十分に低減されず、ＣＳＳ法による前処理では推定精度が劣化することを確認した。また、波源相互相関値を低減させて、信号を読み出せば、通信品質を向上させることができる。そこで、波源の帯域幅が狭い場合にも波源相互相関値を低減することができる相関行列算出装置及び相関行列算出方法が求められていた。また、算出された相関行列に基づいて波源の方向や位置を推定する波源方向又は波源位置推定装置が求められていた。
【０００８】
この発明は、こうした実情に鑑みてなされたものであり、その目的は、波源の帯域幅が狭い場合にも波源相互相関値を低減することができる受信波の相関行列算出装置及び受信波の相関行列算出方法、並びに波源方向推定装置及び波源位置推定装置を提供することにある。
【課題を解決するための手段】
【０００９】
以下、上記目的を達成するための手段及びその作用効果について説明する。
請求項１に記載の発明は、複数の受信素子により電波又は音波を受信して相関行列を算出する受信波の相関行列算出装置において、受信波の方向行列の位相を回転する位相回転手段と、位相回転された方向行列を含む相関行列を平均して平均化された相関行列を算出する平均手段とを備えることをその要旨としている。
【００１０】
同構成によれば、複数の受信素子により電波又は音波を受信して相関行列を算出するにあたって、位相回転手段により受信波の方向行列の位相を回転し、平均手段により周波数平均して平均化された相関行列を算出する。従来のＣＳＳ（Ｃｏｈｅｒｅｎｔ  Ｓｉｇｎａｌ−Ｓｕｂｓｐａｃｅ）法による相関行列の算出では、波源の帯域幅が狭い場合には伝播による位相回転が小さく、波源相互相関値が十分に低減されなかった。しかし、本装置では、位相回転手段により方向行列に位相回転が生じるため、波源相互相関値を低減することができる。
【００１１】
請求項２に記載の発明は、請求項１に記載の受信波の相関行列算出装置において、受信波の波源のおおよその方向を推定する初期方向推定手段と、方向推定された受信波における方向行列を算出する方向行列算出手段とを備えることをその要旨としている。
【００１２】
同構成によれば、複数の受信素子により電波又は音波を受信して相関行列を算出するにあたって、まず初期方向推定手段により従来から行われている、例えばＭＵＳＩＣ（ＭＵｌｔｉｐｌｅ  ＳＩｇｎａｌ  Ｃｌａｓｓｉｆｉｃａｔｉｏｎ）法等により受信波の波源の方向を大まかに推定し、方向行列算出手段により方向推定された受信波の方向行列を算出する。そして、位相回転手段によりこの方向行列の位相を回転し、平均手段により平均して平均化された相関行列を算出する。従来のＣＳＳ法による相関行列の算出では、波源の帯域幅が狭い場合には伝播による位相回転が小さく、波源相互相関値が十分に低減されなかった。しかし、本装置では、位相回転手段により方向行列に位相回転が生じるため、波源相互相関値を低減することができる。
【００１３】
請求項３に記載の発明は、請求項１に記載の受信波の相関行列算出装置において、受信波の波源のおおよその位置を推定する初期位置推定手段と、位置推定された受信波の方向行列を算出する方向行列算出手段とを備えることを要旨としている。
【００１４】
同構成によれば、複数の受信素子により電波又は音波を受信して相関行列を算出するにあたって、まず初期方向推定手段により従来から行われている、例えばＭＵＳＩＣ法等により受信波の波源の位置を大まかに推定し、方向行列算出手段により位置推定された受信波の方向行列を算出する。そして、位相回転手段によりこの方向行列の位相を回転し、平均手段により平均して平均化された相関行列を算出する。従来のＣＳＳ法による相関行列の算出では、波源の帯域幅が狭い場合には伝播による位相回転が小さく、波源相互相関値が十分に低減されなかった。しかし、本装置では、位相回転手段により方向行列に位相回転が生じるため、波源相互相関値を低減することができる。
【００１５】
請求項４に記載の発明は、複数の受信素子により電波を受信して相関行列を算出する受信波の相関行列算出方法において、受信波の方向行列の位相を回転する位相回転と、位相回転された方向行列を含む相関行列を平均して平均化された相関行列を算出する平均処理とを行うことをその要旨としている。
【００１６】
同方法によれば、複数の受信素子により電波を受信して相関行列を算出するにあたって、位相回転により受信波の方向行列の位相を回転し、平均処理により平均して相関行列を算出する。従来のＣＳＳ法による相関行列の算出では、波源の帯域幅が狭い場合には伝播による位相回転が小さく、波源相互相関値が十分に低減されなかった。しかし、本方法では、位相回転により方向行列に位相の回転が生じるため、波源相互相関値を低減することができる。
【００１７】
請求項５に記載の発明は、請求項１又は２に記載の受信波の相関行列算出装置が算出した相関行列に基づいて受信波の波源の方向を推定する波源方向推定手段を備えることを特徴とする波源方向推定装置であることをその要旨としている。
【００１８】
同構成によれば、波源相互相関値を低減した相関行列に基づいて、波源方向推定手段によりＭＵＳＩＣ法等により波源方向を推定する。よって、波源相互相関値が低減された相関行列を用いるため、高精度な波源の方向推定が可能となる。
【００１９】
請求項６に記載の発明は、請求項１又は３に記載の受信波の相関行列算出装置が算出した相関行列に基づいて受信波の波源の位置を推定する波源位置推定手段を備えることを特徴とする波源位置推定装置であることをその要旨としている。
【００２０】
同構成によれば、波源相互相関値を低減した相関行列に基づいて、波源位置推定手段によりＭＵＳＩＣ法等により波源位置を推定する。よって、波源相互相関値が低減された相関行列を用いるため、高精度な波源の位置推定が可能となる。
【００２１】
請求項７に記載の発明は、請求項２に記載の受信波の相関行列算出装置が算出した相関行列に基づいて前記波源の方向を推定する波源方向推定手段を備える波源方向推定装置であって、前記波源方向推定手段は、前記初期方向推定手段による方向推定から前記相関行列の算出までの処理を複数回繰り返させ、算出された相関行列に基づいて前記波源の方向を推定することをその要旨としている。
【００２２】
同構成によれば、波源相互相関値を低減した相関行列に基づいて、波源位置推定手段によりＭＵＳＩＣ法等により波源位置を推定する。そして、この推定された波源方向を初期方向として相関行列の算出処理を繰り返し行うことによって、繰り返し回数に伴って精度がさらに向上されて高精度な波源の方向推定が可能となる。
【００２３】
請求項８に記載の発明は、請求項３に記載の受信波の相関行列算出装置が算出した相関行列に基づいて前記波源の位置を推定する波源位置推定手段を備える波源位置推定装置であって、前記波源位置推定手段は、前記初期位置推定手段による位置推定から前記相関行列の算出までの処理を複数回繰り返させ、算出された相関行列に基づいて前記波源の位置を推定することをその要旨としている。
【００２４】
同構成によれば、波源相互相関値を低減した相関行列に基づいて、波源位置推定手段によりＭＵＳＩＣ法等により波源位置を推定する。そして、この推定された波源位置を初期位置として相関行列の算出処理を繰り返し行うことによって、繰り返し回数に伴って精度がさらに向上されて高精度な波源の位置推定が可能となる。
【発明の効果】
【００２５】
本発明によれば、波源の帯域幅が狭い場合にも波源相互相関値を低減することができる受信波の相関行列算出装置及び受信波の相関行列算出方法、並びに波源方向推定装置及び波源位置推定装置を提供することができる。
【図面の簡単な説明】
【００２６】
【図１】電子キーシステムの構成を示すブロック図。
【図２】車両の通信エリアを示す上面図。
【図３】（ａ）アレーアンテナと波源との位置関係を示す図、（ｂ）波源面（ｚ＝ｚｌ）を示す図、（ｃ）アレー面（ｚ＝０）を示す図。
【発明を実施するための形態】
【００２７】
以下、本発明にかかる相関行列算出装置及び相関行列算出方法を電子キーシステムに適用した一実施形態について図１〜図３を参照して説明する。
図１に示されるように、車両２には、例えば運転者が実際に車両キーを操作しなくてもドアロックの施解錠やエンジンの始動及び停止等の車両動作を行うことが可能な電子キーシステムが搭載されている。電子キーシステムは、キー固有のＩＤコードを無線通信で発信可能な電子キー１が車両キーとして使用されている。電子キーシステムは、車両２からＩＤコード返信要求としてリクエスト信号Ｓｒｑを発信させ、このリクエスト信号Ｓｒｑを電子キー１が受信すると、それに応答する形で電子キー１が自身のＩＤコードを含ませたＩＤコード信号Ｓｉｄを狭域無線通信により車両２に返信し、電子キー１のＩＤコードが車両２のＩＤコードと一致すると、ドアロックの施解錠やエンジンの始動及び停止が許可又は実行されるシステムである。
【００２８】
電子キーシステムを以下に説明すると、車両２には、電子キー１との間で狭域無線通信を行う際にＩＤ照合を行う照合ＥＣＵ（Ｅｌｅｃｔｒｏｎｉｃ  Ｃｏｎｔｒｏｌ  Ｕｎｉｔ）２１と、車両２の状態を管理するメインボディＥＣＵ３１とが設けられている。照合ＥＣＵ２１には、車両２の室内天井に埋設されて（図２参照）車外及び車内にＵＨＦ（Ｕｌｔｒａ  Ｈｉｇｈ  Ｆｒｅｑｕｅｎｃｙ）帯の信号を発信及び受信可能な通信機２４が接続されている。通信機２４には、複数のアンテナ素子からなるアレーアンテナ２４ａが搭載されている。照合ＥＣＵ２１には、例えばドアロック施解錠等を管理するメインボディＥＣＵ３１が車内ＬＡＮ（Ｌｏｃａｌ  Ａｒｅａ  Ｎｅｔｗｏｒｋ）３０を介して接続されている。
【００２９】
一方、電子キー１には、車両２との間で電子キーシステムに準じた無線通信を行う際のコントロールユニットとして通信制御部１１が設けられている。通信制御部１１は、固有のキーコードとしてＩＤコードが記憶されたメモリ１１ａを備えている。通信制御部１１には、外部で発信されたＵＨＦ帯の信号を受信可能な受信部１２と、通信制御部１１の指令に従いＵＨＦ帯の信号を発信可能な発信部１３とが接続されている。
【００３０】
照合ＥＣＵ２１は、通信機２４からＵＨＦ帯のリクエスト信号Ｓｒｑを間欠的に発信させ、車両２を含む周辺に通信エリアＡを形成する（図２参照）。電子キー１がこの通信エリアＡに入り込んでリクエスト信号Ｓｒｑを受信部１２で受信すると、電子キー１はリクエスト信号Ｓｒｑに応答する形で、自身のメモリ１１ａに登録されたＩＤコードを含ませたＵＨＦ帯のＩＤコード信号Ｓｉｄを発信部１３から返信する。照合ＥＣＵ２１は、通信機２４でＩＤコード信号Ｓｉｄを受信すると、ＩＤコード信号Ｓｉｄを発信した電子キー１の位置を位置推定部２１ｂにより推定する。位置推定部２１ｂが受信波の相関行列算出装置であるとともに波源位置推定装置に相当する。照合ＥＣＵ２１は、自身のメモリ２１ａに登録されたＩＤコードと電子キー１のＩＤコードとを照らし合わせてＩＤ照合を行う。照合ＥＣＵ２１は、電子キー１が車外（Ａ−Ａｉ）に位置すると推定して照合が成立する（車外照合成立）と、メモリ２１ａに車外照合フラグを一定時間立てて、メインボディＥＣＵ３１を介してドアロック装置３８によるドアロックの施解錠を許可又は実行する。
【００３１】
また、電子キーシステムには、エンジン始動停止操作の際に実際の車両キー操作を必要とせずに単なるスイッチ操作のみでエンジン始動停止操作を行うことが可能な機能としてワンプッシュエンジンスタートシステムがある。このワンプッシュエンジンスタートシステムを以下に説明すると、車両２には、照合ＥＣＵ２１のＩＤ照合成立結果を基に、エンジンの点火制御及び燃料噴射制御を行うエンジンＥＣＵ３２が設けられている。エンジンＥＣＵ３２は、車内ＬＡＮ３０を通じて照合ＥＣＵ２１等の各種ＥＣＵに接続されている。車両２の運転席には、車両２の電源状態（電源ポジション）を切り換える際に操作されるエンジンスイッチ３５が設けられている。
【００３２】
照合ＥＣＵ２１は、車外照合が成立してドアロックが解錠された後、ドアが開けられて運転者が乗車したことを例えばカーテシスイッチ３７で認識すると、通信機２４からリクエスト信号Ｓｒｑを発信して通信エリアＡを形成する。照合ＥＣＵ２１は、電子キー１がこの通信エリアＡに入り込んで返信してきたＩＤコード信号Ｓｉｄを通信機２４で受信すると、ＩＤコード信号Ｓｉｄを発信した電子キー１の位置を位置推定部２１ｂにより推定する。照合ＥＣＵ２１は、自身に登録されたＩＤコードと電子キー１のＩＤコードとを照らし合わせてＩＤ照合を行う。照合ＥＣＵ２１は、電子キー１が車内（Ａｉ）に位置すると推定して照合が成立する（車内照合成立）と、メモリ２１ａに車内照合フラグを立てる。このとき、メインボディＥＣＵ３１は、エンジンスイッチ３５が押圧操作されると、車両２の電源状態をオフ状態からＡＣＣ（Ａｃｃｅｓｓｏｒｙ）オン状態やＩＧ（Ｉｇｎｉｔｉｏｎ）オン状態に切り換え可能とする。
【００３３】
メインボディＥＣＵ３１は、エンジンが停止している際にブレーキペダルが踏み込み操作された状態でエンジンスイッチ３５が押圧操作されたことを検出すると、照合ＥＣＵ２１に車内照合成立の有無を確認する。メインボディＥＣＵ３１は、照合ＥＣＵ２１から車内照合が成立の旨の通知を受ければ車内照合が成立済みであることを認識して、停止状態のエンジンを始動すべくエンジンＥＣＵ３２に起動信号を出力する。起動信号を受け付けたエンジンＥＣＵ３２は、車内照合結果を確認し、点火制御及び燃料噴射制御を開始してエンジンを始動する。一方、メインボディＥＣＵ３１は、エンジンが稼働している際にエンジンスイッチ３５が押圧操作されたことを検出すると、車両２が停止していることを条件にエンジンを停止状態にする。
【００３４】
次に、照合ＥＣＵ２１の位置推定部２１ｂによる位置推定方法について図３を参照して説明する。本発明では、ＭＵＳＩＣ（ＭＵｌｔｉｐｌｅ  ＳＩｇｎａｌ  Ｃｌａｓｓｉｆｉｃａｔｉｏｎ）法により位置推定を行い、前処理としてＣＳＳ（Ｃｏｈｅｒｅｎｔ  Ｓｉｇｎａｌ−Ｓｕｂｓｐａｃｅ）法を改良したＶＳＳＰ（Ｖｉｒｔｕａｌ  ＳＳＰ：Ｓｐａｔｉａｌ  Ｓｍｏｏｔｈｉｎｇ  Ｐｒｏｃｅｓｓｉｎｇ）法を行う。よって、ＭＵＳＩＣ法、ＣＳＳ法、ＶＳＳＰ法の順番で説明する。なお、位置推定部２１ｂが、初期位置推定手段、方向行列算出手段、位相回転手段、平均手段、及び波源位置推定手段として機能する。
【００３５】
＜平面アレーにおける球面波モードベクトルＭＵＳＩＣ法＞
図３（ａ）、図３（ｂ）、図３（ｃ）は（ｘ，ｙ，ｚ）＝（ｘｌ，ｙｌ，ｚｌ）＝βｌに位置する第ｌ波の波源とアレーアンテナとの位置関係を示した図である。Ｌ波の到来波がｚ＝０のｘｙ平面上のＮｘ×Ｎｙ素子長方形アレーアンテナに入射した場合、Ｇ個の周波数ｆｇ（ｇ＝１，２，…，Ｇ）における各アンテナの受信信号は式（１）のように表すことができる。
【００３６】
【数１】
2010266294.tif 000003
ここで、Ｆｌ（ｆｇ）は周波数ｆｇにおけるｌ番目の波源の振幅である。式（２）はｌ番目の波源と（ｎｘ，ｎｙ）番のアレー素子との距離である。式（３）はｌ番目の波源の周波数ｆｇにおける伝搬による位相遅れである。
【００３７】
式（１）を用いて入力ベクトルＸ（ｆｇ）を次のように定義する。
【００３８】
【数２】
2010266294.tif 000004
このとき、周波数ｆｇにおける振幅ベクトルＦ（ｆｇ）、位相回転行列Φ（ｆｇ）、平面アレーの方向行列Ａ（ｆｇ）をそれぞれ次のように定義する。
【００３９】
【数３】
2010266294.tif 000005
ただし、式（１０）はｌ番目の波源の周波数ｆｇにおける（ｎｘ，ｎｙ）番のアレー素子の基準位置に対する位相遅れであり、φ（ｆｇ，βｌ）は周波数ｆｇにおけるｌ番目の波源から基準位置までの伝搬による位相遅れである。なお、Tは転置を表し、diag［…］は対角行列を表す。この時、平面アレーの入力ベクトルは次のように表すことができる。
【００４０】
【数４】
2010266294.tif 000006
ここで、Ｗ（ｆｇ）はＮ（＝Ｎｘ×Ｎｙ）次元の雑音ベクトルである。更に、式（１３）より、入力ベクトルＸ（ｆｇ）の相関行列Ｒｘｘ（ｆｇ）は次式のように表される。
【００４１】
【数５】
2010266294.tif 000007
ここで、Ｅ［・］は期待値を表す。各アンテナの内部雑音は無相関で、スナップショット数が無限である場合、内部雑音電力をσ２とすると、Ｗ（ｆｇ）＝σ２Ｉ（Ｉ：単位行列）となり、ＭＵＳＩＣ法による推定結果に影響を与えない。実際にはスナップショット数は有限であり、Ｗ（ｆｇ）は単位行列の実数倍とはならず、推定結果に誤差が生じる。
【００４２】
また、全波源がコヒーレントであるとすると、各波源の位相差は周波数によらず一定と考えられるため次式のようになる。
【００４３】
【数６】
2010266294.tif 000008
ここで、ｆｃは中心周波数を表す。よって、全波源がコヒーレントである場合、相関行列Ｒｘｘ（ｆｇ）は次式のように表すことができる。
【００４４】
【数７】
2010266294.tif 000009
＜ＣＳＳ法＞
式（１８）で導かれたＲｘｘ（ｆｇ）は、波源がコヒーレントである場合にはrank（Ｓ（ｆｃ））＝１であるため、ＭＵＳＩＣ法を適用することができない。よって、平均化処理によりランク回復を図る。コヒーレントな広帯域信号を扱う上でよく知られる方法であるコヒーレント信号部分空間（ＣＳＳ）法により、平均化処理を行う。
【００４５】
まず、次のようにＮ次のベクトルＹ（ｆｇ）を定義する。
【００４６】
【数８】
2010266294.tif 000010
ここで、Ｔ（ｆｇ）は焦点行列と呼ばれ、Ｎ×Ｎ次の変換行列であり、次の関係式を満たすよう選ばれる。
【００４７】
【数９】
2010266294.tif 000011
よって、Ａ（ｆｇ），Ａ（ｆｃ）を得るために、まず初期推定を行いおおよその波源位置を知る必要がある。初期推定の方法としては、周波数平均法を前処理としたＭＵＳＩＣ法等が挙げられる。ただし、システムによってはおおよその波源位置が既知である場合も考えられる。このような場合には、必ずしも初期推定を行う必要は無い。また、後述する変換誤差の影響を考慮する必要が無い場合には、適当な波源位置を想定しＡ（ｆｇ），Ａ（ｆｃ）を決定することで、初期推定を行わないようにしてもよい。
【００４８】
このようにして求められたＴ（ｆｇ），Ｙ（ｆｇ）を用い、次式で与えられる相関行列Ｒｙｙバーを求める。
【００４９】
【数１０】
2010266294.tif 000012
よって、Ｇ≧ＬとするとΦ（ｆｇ）（ｇ＝１，…，Ｇ）の独立性が高い場合には式（２３）のランクがＬとなり、ＭＵＳＩＣ法が適用可能となる。しかし、周波数帯域幅が十分でない場合、Φ（ｆｇ）は周波数によらずほぼ同一となり、各波源の相互相関値が十分に低下しない。この場合、雑音行列Ｗバーによる推定誤差が大きくなる。
【００５０】
＜ＶＳＳＰ法＞
ＶＳＳＰ法は、基本的には周波数帯域幅が０であっても波源相互相関値を低下させることが可能である。よって、ここでは受信波の中心周波数ｆｃの受信信号のみを用いることとする。また、ＶＳＳＰ法では、ＣＳＳ法における変換行列Ｔ（ｆｇ）が満たすべき式（２０）を、式（２５）を満たす変換行列Ｔｎ（ｎ＝１，２，…，Ｎ）に変更する。
【００５１】
【数１１】
2010266294.tif 000013
ここで、Ｄｎは方向行列Ａ（ｆｃ）の各列の位相回転を表す対角行列である。この時、変換された受信信号Ｙｎ、相関行列Ｒｙｙバーは次式の通りとなる。
【００５２】
【数１２】
2010266294.tif 000014
以上より、Ｄｎを適切に選ぶことで各波源の相互相関値は低下し、推定精度が向上する。例えば、Ｄｎを次式のように選べば各波源の相互相関値はすべて０となる。なお、δは任意の実定数である。
【００５３】
【数１３】
2010266294.tif 000015
＜初期推定誤差と繰り返し推定＞
通常、初期推定値には誤差が含まれる。この誤差を含む初期推定位置β'ｌより得られる方向行列をＡ'（ｆｃ）とし、これらを用いて算出されるＶＳＳＰ法における変換行列をＴｎ'とする。このＴｎ'を用い真の波源位置を表す方向行列Ａ（ｆc）を変換すると次式のように表される。
【００５４】
【数１４】
2010266294.tif 000016
式（３２）では、変換誤差行列Ｐｎが発生する。一方、ＣＳＳ法においてはβ'ｌより得られる方向行列Ａ'（ｆｃ）、Ａ'（ｆｇ）を用いて算出される変換行列Ｔ'（ｆｇ）を用い、真の波源位置を表す方向行列Ａ（ｆｇ）を変換すると、次式のようになる。
【００５５】
【数１５】
2010266294.tif 000017
しかしながら、やはり同様に変換誤差行列が発生する。よって、Ｔ'ｎを用いてＶＳＳＰ法を適用した相関行列Ｒｙｙ’バーは、熱雑音が無いと仮定すると次式のように表すことができる。
【００５６】
【数１６】
2010266294.tif 000018
Ｐバー（式（３５））は変換誤差のために発生する項であり、熱雑音と同様に推定誤差の原因となる。よって、各波源の相互相関値が高い場合にはＰバーに起因する推定誤差が大きくなる。条件によっては、Ｒｙｙ’バー（式（３４））を用いて推定し得られる位置β''ｌは初期推定位置β'ｌよりも真の波源位置から遠くなる可能性がある。ＣＳＳ法では各波源の相互相関値が十分に低減できない可能性があるが、ＶＳＳＰ法では各波源の相互相関値を０とすることが可能であるため、β''ｌがβ'ｌより真の波源位置から遠くなる可能性を低減できる。
【００５７】
また、β''ｌがβ'ｌより真の波源位置に近くなった場合、β''ｌを初期推定値として再びＶＳＳＰ法を適用すると、変換誤差行列Ｐｎ，Ｐバーは小さくなり、推定結果β'''ｌは更に真の波源位置に近づくと考えられる。つまり、ＶＳＳＰ法を繰り返し適用することにより、さらなる推定精度を向上できる。
【００５８】
さて、狭帯域かつコヒーレントな近傍波源を、平面アレーによりＭＵＳＩＣ法で位置推定する場合に有効な前処理の手法として、ＶＳＳＰ法を説明した。初期推定誤差が無いと仮定した場合の位置推定精度は、比帯域が狭くＣＳＳ法では推定精度が劣化する状況においても、ＶＳＳＰ法では良好な推定精度が得られる。次に、初期推定誤差が存在する場合において、比帯域が狭い場合にはＣＳＳ法に対してＶＳＳＰ法の推定精度が良好である。
【００５９】
以上、説明した実施形態によれば、以下の作用効果を奏することができる。
（１）アレーアンテナ２４ａを備えた通信機２４により電波を受信して、位置推定部２１ｂにより相関行列を算出するにあたって、まず従来から行われている、周波数平均法を前処理として用いたＭＵＳＩＣ法により受信波の波源の方向を大まかに推定し、受信波の方向行列を算出する。そして、この方向行列の位相を回転し、平均処理して平均化された相関行列を算出する。従来のＣＳＳ法による相関行列の算出では、波源の帯域幅が狭い場合には伝播による位相回転が小さく、波源相互相関値が十分に低減されなかった。しかし、本装置では、位相回転させて方向行列に位相回転が生じるため、波源相互相関値を低減することができる。
【００６０】
（２）受信波の単一の周波数における相関行列を用いることにより、Φ（ｆg）による位相回転を考慮せずに、波源を無相関化できる。
（３）方向行列の位相を回転させることにより、波源相互相関値を低減した相関行列に基づいて波源位置を推定する。よって、高精度な波源の位置推定が可能となる。
【００６１】
（４）推定された波源の位置を初期推定位置として相関行列の算出処理を繰り返し行うことによって、繰り返し回数に伴って精度がさらに向上されて高精度な波源の位置推定が可能となる。
【００６２】
なお、上記実施形態は、これを適宜変更した以下の形態にて実施することができる。
・上記実施形態において、式（２０），（２１）を満たす変換行列Ｔ（ｆｎ）は次式により得ることができる。
【００６３】
【数１７】
2010266294.tif 000019
ここで、Ｕ（ｆｎ），Ｖ（ｆｎ）は、式（３７）を特異値分解した場合の左特異ベクトル、右特異ベクトルである。
【００６４】
・上記実施形態において、変換行列Ｔ（ｆｇ）は式（２１）を満たさなくてもよい。
・上記実施形態において、式（２０）を満たし、式（２１）を満たさない変換行列Ｔ（ｆｇ）は、例えば次式により得ることができる。なお、＋は一般逆行列を表す。
【数１８】
2010266294.tif 000020
【００６５】
・上記実施形態において、Ｄｎは式（２９）が無相関となる行列の例であり、その他の方法により算出してもよい。また、Ｄｎは、波源が無相関とならなくても、波源相互相関値が低減されるようなものであればよい。
【００６６】
・上記実施形態では、波源位置推定をＭＵＳＩＣ法により行ったが、ＭＵＳＩＣ法に限らず、ＥＳＰＲＩＴ法等の他の波源位置推定方法により行ってもよい。
・上記実施形態では、初期位置推定を、周波数平均法を前処理としたＭＵＳＩＣ法により行ったが、ビームフォーマ法等の他の波源位置推定方法により行ってもよい。
【００６７】
・上記実施形態では、受信波の中心周波数の相関行列を用いるようにしたが、受信波の中心周波数に限らず、受信波の他の周波数の相関行列を用いてもよい。
・上記実施形態では、推定された波源方向を初期方向として相関行列の算出処理を繰り返し行うようにしたが、所望の推定精度に合わせて適宜変更可能である。また、１回で所望の推定精度を得ることができれば、繰り返さなくてもよい。
【００６８】
・上記実施形態では、波源位置推定装置を電子キーシステムの電子キー１の位置推定に適用したが、電子キーシステムに限らず、種々の通信機器に適用してもよい。
・上記実施形態では、波源位置推定装置に初期推定手段及び方向行列算出手段を備えたが、初期推定手段及び方向行列算出手段を省略してもよい。
【００６９】
・また、波源位置推定手段を備えた波源位置推定装置に限らず、波源方向推定手段を備えた波源方向推定装置として種々の通信機器に適用してもよい。
・また、波源方向推定手段を備えた波源方向推定装置は、初期方向推定手段及び方向行列算出手段を設けてもよいが、初期方向推定手段及び方向行列算出手段を省略してもよい。
【００７０】
・また、波源位置推定装置及び波源方向推定装置に限らず、相関行列算出装置として種々の通信機器に適用してもよい。この場合、波源相互相関値が低減された信号を読み出せば、通信品質を向上させることができる。
【００７１】
・波源位置推定装置、波源方向推定装置、及び相関行列算出装置等の装置に限らず、波源位置推定方法、波源方向推定方法、及び相関行列算出方法を、アレーアンテナを備えた種々の通信機器の演算に適用してもよい。
【００７２】
・上記構成において、図１に示されるように、電子キー１に施錠スイッチ１９ａ及び解錠スイッチ１９ｂを設け、この施錠スイッチ１９ａ及び解錠スイッチ１９ｂを操作することで、車両２へ施錠信号Ｓｌ又は解錠信号Ｓｕｌを発信して、車両２のドアロックの解錠又は施錠を可能としてもよい。
【００７３】
・電子キーシステムで使用する各通信の周波数はどの周波数の電波を使用してもよい。
・上記実施形態において、波源は電波を発信する装置であったが、これに限らず、例えばスピーカーのような音波を発信する装置であっても良い。
【００７４】
・上記実施形態において、受信素子は電波を受信する装置であったが、これに限らず、例えばマイクロホンのような音波を受信する装置であっても良い。
・上記実施形態では、電波に適用したが、音波に適用してもよい。
【００７５】
次に、前記実施形態から把握できる技術的思想をその効果と共に記載する。
（イ）請求項４に記載の受信波の相関行列算出方法において、受信波の波源のおおよその方向を推定する初期方向推定と、方向推定された受信波の方向行列を算出する方向行列算出とを行うことを特徴とする受信波の相関行列算出方法。
【００７６】
同方法によれば、まず初期方向推定を従来から行われている、例えば周波数平均法を前処理に用いたＭＵＳＩＣ法等により受信波の波源の方向を大まかに推定し、方向行列算出により方向推定された受信波の各周波数における方向行列を算出する。そして、前記方向行列を用いて位相回転により受信波の方向行列の位相を回転することにより、回転の際に発生する誤差を低減することができる。
【００７７】
（ロ）請求項４に記載の受信波の相関行列算出方法において、受信波の波源のおおよその位置を推定する初期位置推定と、位置推定された受信波の方向行列を算出する方向行列算出とを行うことを特徴とする受信波の相関行列算出方法。
【００７８】
同方法によれば、まず初期方向推定を従来から行われている、例えば周波数平均法を前処理に用いたＭＵＳＩＣ法等により受信波の波源の位置を大まかに推定し、方向行列算出により位置推定された受信波の各周波数における方向行列を算出する。そして、前記方向行列を用いて位相回転により受信波の方向行列の位相を回転することにより、回転の際に発生する誤差を低減することができる。
【００７９】
（ハ）請求項４又は付記（イ）に記載の受信波の相関行列算出方法が算出した相関行列に基づいて受信波の波源の方向を推定する波源方向推定を行うことを特徴とする波源方向推定方法。
【００８０】
同方法によれば、波源相互相関値を低減した相関行列に基づいて、ＭＵＳＩＣ法等により波源方向を推定する。よって、波源相互相関値が低減された相関行列を用いるため、高精度な波源の方向推定が可能となる。
【００８１】
（ニ）請求項４又は付記（ロ）に記載の受信波の相関行列算出方法が算出した相関行列に基づいて受信波の波源の位置を推定する波源位置推定を行うことを特徴とする波源位置推定方法。
【００８２】
同方法によれば、波源相互相関値を低減した相関行列に基づいて、ＭＵＳＩＣ法等により波源位置を推定する。よって、波源相互相関値が低減された相関行列を用いるため、高精度な波源の位置推定が可能となる。
【００８３】
（ホ）付記（イ）に記載の受信波の相関行列算出方法が算出した相関行列に基づいて前記波源の方向を推定する波源方向推定を行う波源方向推定方法であって、前記波源方向推定は、前記初期方向推定手段による方向推定から前記相関行列の算出までの処理を複数回繰り返させ、算出された相関行列に基づいて前記波源の方向を推定することを特徴とする波源方向推定方法。
【００８４】
同方法によれば、波源相互相関値を低減した相関行列に基づいて、ＭＵＳＩＣ法等により波源方向を推定する。そして、この推定された波源方向を初期方向として相関行列の算出処理を繰り返し行うことによって、繰り返し回数に伴って精度がさらに向上されて高精度な波源の方向推定が可能となる。
【００８５】
（ヘ）付記（ロ）に記載の受信波の相関行列算出方法が算出した相関行列に基づいて前記波源の位置を推定する波源位置推定を行う波源位置推定方法であって、前記波源位置推定は、前記初期位置推定手段による位置推定から前記相関行列の算出までの処理を複数回繰り返させ、算出された相関行列に基づいて前記波源の位置を推定することをその要旨としている。
【００８６】
同方法によれば、波源相互相関値を低減した相関行列に基づいて、ＭＵＳＩＣ法等により波源位置を推定する。そして、この推定された波源位置を初期位置として相関行列の算出処理を繰り返し行うことによって、繰り返し回数に伴って精度がさらに向上されて高精度な波源の位置推定が可能となる。
【符号の説明】
【００８７】
１…電子キー、２…車両、１１…通信制御部、１１ａ…メモリ、１２…受信部、１３…発信部、２１…照合ＥＣＵ、２１ａ…メモリ、２１ｂ…位置推定部、２４…通信機、２４ａ…アレーアンテナ、３５…エンジンスイッチ、３８…ドアロック装置、Ａ…通信エリア、Ａｉ…車内エリア、Ｓｉｄ…ＩＤコード信号、Ｓｒｑ…リクエスト信号。
【図１】
2010266294.tif 000021
【図２】
2010266294.tif 000022
【図３】
2010266294.tif 000023