JP3640891B2 - 到来波方向推定装置及び到来波方向推定方法 - Google Patents
到来波方向推定装置及び到来波方向推定方法 Download PDFInfo
- Publication number
- JP3640891B2 JP3640891B2 JP2001019010A JP2001019010A JP3640891B2 JP 3640891 B2 JP3640891 B2 JP 3640891B2 JP 2001019010 A JP2001019010 A JP 2001019010A JP 2001019010 A JP2001019010 A JP 2001019010A JP 3640891 B2 JP3640891 B2 JP 3640891B2
- Authority
- JP
- Japan
- Prior art keywords
- correction amount
- arrival
- antenna elements
- deviation
- reception
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
Images
Landscapes
- Variable-Direction Aerials And Aerial Arrays (AREA)
- Mobile Radio Communication Systems (AREA)
- Radio Transmission System (AREA)
Description
【０００１】
【発明の属する技術分野】
本発明は、移動通信や構内無線やレーダ装置などに利用される到来波方向推定装置及び到来波方向推定方法に関する。
【０００２】
【従来の技術】
近年、送受共用のアダプティブアンテナ装置に対する需要が高まっており、例えばＴＤＭＡ／ＴＤＤの通信システムにおいては、送受の周波数が一致することから、送受で同一のパターンを形成するアダプティブアンテナ装置を利用することが有効である。
【０００３】
図１２は従来の送受共用のアダプティブアンテナ装置の構成例を示したブロック図である。アダプティブアンテナを形成する複数のアンテナ素子３１ａ、…、３１ｎがあり、個々のアンテナ素子には信号を受信する受信アナログ部３２ａ、…、３２ｎと、信号を送信する送信アナログ部３３ａ、…、３３ｎが切り替えで接続されている。
【０００４】
受信アナログ部３２ａ、…、３２ｎの受信信号は受信パターン形成回路３４に入力されて最終的な受信信号になる。送信信号は送信パターン形成回路３５に入力され、個々の送信アナログ部３３ａ、…、３３ｎに分配される。
【０００５】
このような送受信共用のアダプティブアンテナ装置は、個々のアンテナ素子３１ａ、…、３１ｎに対する受信信号に基づいて、個々のアンテナ素子３１ａ、…、３１ｎに対応する複素振幅の重み付けを受信パターン形成回路３４で決定し、この重みを共通のアンテナ素子に対応する受信信号に付与し、或いは前記重み付けを送信パターン形成回路３５で送信信号に付与することで、所望のアンテナパターンを形成している。従って原理的には、受信と送信で同一のパターンを形成することができることになる。
【０００６】
しかしながら、受信と送信で同一のパターンを形成することは現実的には困難である。その理由は、多くの場合、アダプティブアンテナ内のアンテナ素子毎の送信アナログ部３３、受信アナログ部３２の通過特性（通過振幅特性、通過位相特性）に偏差（ばらつき）があることに起因している。
【０００７】
このような偏差があると、図１３に示すように、受信のパターン５０と送信パターン６０が異なり、所望方向への利得が低くなったり、干渉ユーザーへの利得が高くなったりしてしまう。
【０００８】
或いは、受信アダプティブアンテナ装置において、方向拘束付のアダプティブアルゴリズムを用いる場合がある。この場合にも、受信アナログ部の通過特性に偏差（ばらつき）があると、最適な解を得ることができなくなる。
【０００９】
更に、受信アダプティブアンテナ装置はレーダ装置の一部として、電波の到来方向を推定することに使用されるが、受信アナログ部の通過特性に偏差（ばらつき）があると、正確に電波の到来方向を推定することができなくなる。
【００１０】
上記のような偏差を取り除く方法として、一般的には、キャリブレーションが行われる。キャリブレーションとは、アダプティブアンテナ装置の製作が終了した時点で、すべての受信アナログ部と送信アナログ部の偏差を測定し、補正量 （偏差量に同じ）を予め算出しておく。そして、実際にアダプティブアンテナ装置を運用するときには、補正量を考慮して各素子に設定する振幅位相（重み）を算出すればよい。
【００１１】
しかしながら、最初に求めた補正量が運用中に変化してしまう場合がある。例えば、昼夜の外気温度変化の違いによるアナログ部の部品（増幅器、周波数変換機、フィルタなど）の特性変化、そして、それらを接続するケーブルの長さの伸縮が考えられる。このような変化が生じると、受信アダプティブアンテナ装置はその性能が大幅に劣化して到来方向を正確に推定できなくなり、更に、これにより加入者容量が大幅に劣化して、最悪的には、通信の切断にまで及ぶ場合がある。つまり、キャリブレーションは一度では不十分であり、運用中にリアルタイムにキャリブレーションすることが要求される。
【００１２】
このような問題を解決する方法として、図１４は特開平２０００−９１８３３の発明例を示したもので、アンテナ素子４１ａ、４１ｂ、４１ｃで受信した受信信号はビームフォーマー４２を通ってビームスペースデータに変換される。ビームフォーマーの出力（ビームスペースデータ）は、レベル調整器４４ｂ，４４ｃにより補正される。補正されたデータは、逆変換回路によって再び、アンテナ素子で受信した受信信号（これは補正されている）に変換される。
【００１３】
位相補正量演算部４３は各、ビームフォーマーの出力間の相対的な位相差を検出する。ここでは、到来波がひとつであり、アンテナ素子４１ａ、４１ｂ、…、４１ｃの配置位置が決まっていれば、ビームフォーマーの出力間の位相差は到来角度によって決まってくる。そこで、位相補正量演算部４３では、検出した位相差と本来得られる値との差を補正量として算出し、その補正量によりアンテナ素子４１ｂ、４１ｃの受信レベルをレベル調整器４４ｂ、４４ｃにより調整して、到来波の推定方向の精度を高めることができる。しかも、この従来例では、運用中に上記の補正をリアルタイムに行うことができるメリットがある。
【００１４】
しかしながら、到来波の数は、ひとつに限定され、また、位相のみの補正になっているので、振幅の変動があった場合には、偏差を保証することができない。特に、アダプティブアンテナ装置では、複数の到来波を受信している状態でパターン形成を行うものである。到来してくる電波の数をひとつにしなければいけないので、通常の通信中にキャリブレーションを行うことができない。
【００１５】
以上説明したように、従来のアダプティブアンテナ装置に対して発明されたキャリブレーション方法は、リアルタイムに行えるものであるが、到来波の数がひとつに限定され、位相しか補正できないという問題点があった。
【００１６】
また、以上のアダプティブアンテナ装置に関する従来の問題点は、受信アダプティブ装置と同様の構成を一部に持つアンテナレーダ装置における到来方向の推定精度を向上させる場合においても、まったく同様な問題点となる。
【００１７】
【発明が解決しようとする課題】
以上説明したように従来の到来波方向推定装置では、回路の特性偏差を取り除く方法として、一般的にキャリブレーションが行われるが、最初に求めた補正量が運用中に変化してしまう場合があるため、リアルタイムにキャリブレーションを行なわなければならない。しかし、その場合、到来波はひとつに限定され、しかも位相のみの補正しかできず、しかも複数の到来波がある運用中にキャリブレーションを行うことができないため、実用上不十分なものであった。
【００１８】
本発明は、上述の如き従来の課題を解決するためになされたもので、運用中の受信データを用いてリアルタイムに高精度なキャリブレーションを行うことができ、到来波方向推定精度を向上させることができる到来波方向推定装置及び到来波方向推定方法を提供することを目的とする。
【００１９】
【課題を解決するための手段】
上記目的を達成するために、本発明の特徴は、到来波方向推定装置において、複数のアンテナ素子と、前記複数のアンテナ素子に対応して設けられ、対応するアンテナ素子より伝達される受信信号を前記アンテナ素子対応で処理する複数の受信部と、前記複数の受信部から出力される受信信号から前記複数の受信部から出力される信号間の相対的な偏差を推定する補正量算出手段と、前記複数の受信部から出力される受信信号を前記補正量算出手段により得られる偏差量を用いて補正し、これら補正された受信信号を用いて前記複数のアンテナ素子に到来する少なくともひとつの電波の到来方向を推定する到来方向推定手段と、を具備し、前記補正量算出手段は、ひとつ或いは複数の電波の到来方向と、複数の受信部から出力される信号間の振幅位相偏差と、複数の受信部から出力される信号間の相互結合係数の３種類の未知変数を求める処理を行い、ステップ１として振幅位相偏差と相互結合係数を固定して到来方向を推定し、ステップ２として到来方向と相互結合係数を固定して振幅位相偏差を推定し、ステップ３として到来方向と振幅位相偏差を固定して相互結合係数を推定し、前記ステップ１からステップ３の処理を推定結果が収束するまで繰り返し実行することにより、前記偏差量を求めることにある。
【００２０】
本発明の他の特徴は、到来波方向推定装置において、複数のアンテナ素子と、前記複数のアンテナ素子に対応して設けられ、対応するアンテナ素子より伝達される受信信号を前記アンテナ素子対応で処理する複数の受信部と、前記複数の受信部から出力される受信信号から前記複数の受信部から出力される信号間の相対的な偏差を推定する補正量算出手段と、前記複数の受信部から出力される受信信号を前記補正量算出手段により得られる偏差量を用いて補正し、これら補正された受信信号を用いて前記複数のアンテナ素子に到来する少なくともひとつの電波の到来方向を推定する到来方向推定手段と、を具備し、前記補正量算出手段は、ステップ１として振幅位相偏差と相互結合係数を固定して到来方向を推定し、ステップ２として到来方向と相互結合係数を固定して振幅位相偏差を推定し、ステップ３として振幅位相偏差と相互結合係数を固定して到来方向を推定し、ステップ４として到来方向と振幅位相偏差を固定して相互結合係数を推定し、前記ステップ１からステップ４の処理を推定結果が収束するまで繰り返し実行することにより、前記偏差量を求めることにある。
【００２１】
本発明の他の特徴は、到来波方向推定装置において、複数のアンテナ素子と、前記複数のアンテナ素子に対応して設けられ、対応するアンテナ素子より伝達される受信信号を前記アンテナ素子対応で処理する複数の受信部と、前記複数の受信部から出力される受信信号から前記複数の受信部から出力される信号間の相対的な偏差を推定する補正量算出手段と、前記複数の受信部から出力される受信信号を前記補正量算出手段により得られる偏差量を用いて補正し、これら補正された受信信号を用いて前記複数のアンテナ素子に到来する少なくともひとつの電波の到来方向を推定する到来方向推定手段と、を具備し、前記到来波方向推定装置は、偏差量算出のためのひとつ或いは複数の電波の送信装置を有しており、前記電波の送信装置の電波発信位置は、前記到来波推定装置に対して相対的な方向が既知であり、前記補正量算出手段は、前記送信装置から発せられた既知の方向から到来する電波の受信信号を用いて、ステップ１として到来方向と相互結合係数を固定して振幅位相偏差を推定し、ステップ２として到来方向と振幅位相偏差を固定して相互結合係数を推定し、前記ステップ１とステップ２の処理を推定結果が収束するまで繰り返し実行することにより、前記偏差量を求めることにある。
【００２２】
本発明の他の特徴は、到来波方向推定装置において、複数のアンテナ素子と、前記複数のアンテナ素子に対応して設けられ、対応するアンテナ素子より伝達される受信信号を前記アンテナ素子対応で処理する複数の受信部と、前記複数の受信部から出力される受信信号から前記複数の受信部から出力される信号間の相対的な偏差を推定する補正量算出手段と、前記複数の受信部から出力される受信信号を前記補正量算出手段により得られる偏差量を用いて補正し、これら補正された受信信号を用いて前記複数のアンテナ素子に到来する少なくともひとつの電波の到来方向を推定する到来方向推定手段と、を具備し、前記補正量算出手段は補正量記憶装置を有しており、偏差量を算出する毎に得られた偏差量を補正量記憶装置に記憶して記憶内容を更新し、前記補正量記億装置は、前記補正量算出装置によって算出された偏差量と、そのときの外気気温の両方を記憶することにある。
【００２３】
本発明の他の特徴は、到来波方向推定装置において、複数のアンテナ素子と、前記複数のアンテナ素子に対応して設けられ、対応するアンテナ素子より伝達される受信信号を前記アンテナ素子対応で処理する複数の受信部と、前記複数の受信部から出力される受信信号から前記複数の受信部から出力される信号間の相対的な偏差を推定する補正量算出手段と、前記複数の受信部から出力される受信信号を前記補正量算出手段により得られる偏差量を用いて補正し、これら補正された受信信号を用いて前記複数のアンテナ素子に到来する少なくともひとつの電波の到来方向を推定する到来方向推定手段と、を具備し、前記補正量算出手段は補正量記憶装置を有しており、偏差量を算出する毎に得られた偏差量を補正量記憶装置に記憶して記憶内容を更新し、前記補正量記憶装置は、前記補正量算出装置によって算出された偏差量と、そのときの受信アナログ部の周りの気温の両方を記憶することにある。
【００２４】
本発明の他の特徴は、到来波方向推定装置において、複数のアンテナ素子と、前記複数のアンテナ素子に対応して設けられ、対応するアンテナ素子より伝達される受信信号を前記アンテナ素子対応で処理する複数の受信部と、前記複数の受信部から出力される受信信号から前記複数の受信部から出力される信号間の相対的な偏差を推定する補正量算出手段と、前記複数の受信部から出力される受信信号を前記補正量算出手段により得られる偏差量を用いて補正し、これら補正された受信信号を用いて前記複数のアンテナ素子に到来する少なくともひとつの電波の到来方向を推定する到来方向推定手段と、を具備し、前記補正量算出手段は補正量記憶装置を有しており、偏差量を算出する毎に得られた偏差量を補正量記憶装置に記憶して記憶内容を更新し、前記補正量記憶装置は、前記補正量算出手段によって算出された偏差量と、そのときの時刻、日時の両方を記憶することにある。
【００３４】
【発明の実施の形態】
以下、本発明の実施形態を図面に基づいて説明する。図１は、本発明の到来波方向推定装置の第１の実施形態に係る構成を示したブロック図である。到来波方向推定装置は、複数のアンテナ素子１ａ、…、１ｎと、複数のアンテナ素子１ａ、…、１ｎに対応して設けられ、対応するアンテナ素子より伝達される受信信号を処理する複数の受信アナログ部２ａ、…、２ｎと、受信アナログ部２ａ、…、２ｎの出力信号をデジタル変換するＡ／Ｄ変換器３ａ、…、３ｎと、複数の受信アナログ部２ａ、…、２ｎの相対的な偏差を推定する補正量算出手段４と、補正量算出手段４の出力する補正量を用いて複数のアンテナ素子１ａ、…、１ｎに到来するひとつ或いは複数の電波の到来方向を推定する到来方向推定手段５とから構成される。
【００３５】
ここで、受信アナログ部２（受信アナログ部２ａ、…、２ｎの総称）は、所望信号を増幅する低雑音増幅器２１と、所望の周波数帯域のみを抽出するフィルタ２２と、周波数帯をＲＦ（電波周波数）からＢＢ（ベースバンド）又はＩＦ（中間周波数）へ変換するための周波数変換器２３とから構成されている。この受信アナログ部２は、受信信号を抽出し、周波数変換し、増幅する機能を有するものであれば、他のいかなる構成を有するものに置き換えて構わない。
【００３６】
次に本実施形態の動作について説明する。補正量計算手段４は、ひとつ或いは複数の到来してくる電波の受信信号を用い、ＭＵＳＩＣ（Multiple Signal Classification）アルゴリズムで推定される複数の電波の到来方向情報を利用して補正量を算出する第１の実施例の方法を採用している。
【００３７】
ここで、複数のアンテナ素子１ａ、１ｎはどのようなものであっても良い。例えば、ダイポールアンテナを用いれば、水平面内で概略無指向性が得られるため、ひとつの到来波方向推定装置で３６０度をカバーできる。例えば、ホーンアンテナを用いれば各アンテナは指向性を有しており、限定した領域のみを到来波方向推定手段５でカバーすることになる。こうすると、領域が限定されるものの、後方からの到来波の影響を受けずに電波の方向を推定することができる。ここで、アンテナの放射パターンは予め測定しておき、到来方向推定手段５と補正量算出手段４において使用するものとする。
【００３８】
次に受信アナログ部２であるが、ここでは、受信信号が低雑音増幅器２１で増幅され、フィルタ２２により所望の周波数帯域のみが抽出され、周波数変換器２３により周波数がベースバンド（ＢＢ）または中間周波（ＩＦ）へ変換される。ここで、各受信アナログ部２ａ、…、２ｎは、偏差（ばらつき）を有しており、ブランチ毎に異なった通過特性を有することとなる。また、受信アナログ部２ａ、…、２ｎ間には、シールドが完全でない場合に相互結合が生じることがある。つまり、ある受信アナログ部で受信した信号が、他の受信アナログ部へ漏れこむ現象が発生する。
【００３９】
次に各受信信号は、各Ａ／Ｄ変換器３ａ、…、３ｎによってデジタル信号に変換され、到来方向推定手段５において処理され、到来波の方向が推定される。
【００４０】
一方、補正量算出手段４は、各Ａ／Ｄ変換器３ａ、…、３ｎによってデジタル信号に変換された受信信号を用いて補正量を算出する。この補正量が正確に求まったならば、どのような偏差があったとしても（ケーブルが切断され振幅が０となった場合を除く）、それを補正できるものであり、精度良く到来波方向を推定できるようになる。また、到来波方向推定手段５は、どのような方法によって方向を推定するものでも構わない。ＭＵＳＩＣアルゴリズムを用いても、ＥＳＰＲＩＴ（Estimation of Signal Parameters via Rotational In rariance Techniques）アルゴリズムを用いても、モノパルス側角であっても、ビームフォーマー法でも構わない。
【００４１】
図２は受信アナログ部２の偏差を回路モデルで表した図である。複数のアンテナ１ａ、…、１ｎで受信した信号は、Γｉｉ（ｉ＝１〜ｎ）（ｎは素子数）倍される。ここでΓｉｉは通過特性をあらわす複素数であり、振幅と位相の偏差を意味している。つまり、全てのΓｉｉがｌｅｘｐ（ｊ０゜）ならば、通過特性はすべての受信アナログ部２ａ、…、２ｎで同一であることを意味している。一方、Γｉｉが受信アナログ部２ａ、…、２ｎによって異なる場合には、通過特性を補正する必要がある。
【００４２】
次に、Γｉｉ倍された受信信号は、Ｃｉｊ（ｉ＝１〜ｎ，ｊ＝１〜ｎ）倍されて、他の受信アナログ部２に漏れこむ。Ｃｉｊとは受信アナログ部ｉ番目とｊ番目の相互結合係数である。ここで、Ｃｉｉは常に１である。また、Ｃｉｊ（ｉ≠ｊ）＝０ならば、結合は発生していない状態を意味している。このように、受信アナログ部２ａ、…、２ｎでは、振幅位相偏差と相互結合の二つの補正をしなければならない未知変数がある。
【００４３】
尚、振幅位相偏差と相互結合の２つは受信アナログ部２ａ、…、２ｎのどの部分に発生しているかは問題としない。つまり、受信アナログ部２の全体の偏差を求めることが重要であるからである。また、振幅位相偏差と相互結合の回路モデルの順番も問題ではない。なぜならば、順番を逆にしたとしても、受信アナログ部２の全体の偏差としては等価なものが得られるからである。
【００４４】
次に補正量算出手段４であるが、ひとつ或いは複数の到来してくる電波の受信信号を用いるものである。つまり、任意の到来電波を使用することができるので、アダプティブアンテナ装置を運用中にリアルタイムにキャリブレーションを行うことができ、また、キャリブレーション用の特別な送信源を必要としない。そして、ＭＵＳＩＣアルゴリズムで推定される複数の電波の到来方向情報を利用して補正量を算出する。
【００４５】
一般的に、キャリブレーションが正確に行われていたならば、ＭＵＳＩＣアルゴリズムによって推定される到来方向においては、図３の実線に示すように、鋭いピークを有するものとなる。なお、ＭＵＳＩＣスペクトラムは、アンテナの放射パターンと、受信信号と、補正量（キャリブレーションデータ）から算出されるものである。しかしながら、偏差がある場合には、図３の破線に示すように推定方向のピークが鈍り、推定精度が劣化してしまう。
【００４６】
ここでの補正量算出手段４は、ＭＵＳＩＣによって得られるＭＵＳＩＣスペクトラムのピークを鋭くするように、補正量（図２におけるΓｉｉとＣｉｊ）を算出するものである。算出する方法は、いかなる方法によっても良い。例えば、考えうる全ての組み合わせの補正量からもっともスペクトラムが鋭くなるものを選んでもよい。または、ジェネティックアルゴリズム（遺伝的手法）によっても良いし、ニューロを用いた方法でも良い。以下に、Γを求めるための計算方法の一例を示す。なお、この計算方法は文献（B.Friedlander,et al.,"Direction finding in the presence of mutual coupling,"IEEE Trans.Antenna and Propagation,vol.39,no.3,pp.273-284,March 1991.)を参考にしている。
【００４７】
複数のアンテナ（ここではｎ素子）１ａ、…、１ｎで受信した受信ベクトルは以下の式で与えられる。
【００４８】
Ｘ（ｔ）＝ＡＳ（ｔ）＋Ｎ（ｔ）
ここでＸ（ｔ）は受信ベクトル、Ａは到来方向におけるアレイマニフォルド、Ｓ（ｔ）は信号ベクトル、Ｎ（ｔ）はノイズベクトルである。
【００４９】
そして、偏差のある受信アナログ部２を通過した後の受信ベクトルＸ“（ｔ）は以下のようになる。
【００５０】
Ｘ“（ｔ）＝ＣΓＡＳ（ｔ）＋Ｎ（ｔ）
Γ＝ｄｉａｇ｛Γ１１、Γ２２、…、Γｎｎ｝
Ｃは、要素Ｃｉｊがｉ番目とｊ番目の受信アナログ部２の相互結合係数を意味する行列である。
【００５１】
この例では、振幅位相Γｉｉを算出する方法を記述する。本計算では、以下の評価関数Ｊｃを最小化する。この評価関数は、到来方向におけるＭＵＳ１Ｃスペクトラムの合計値の逆数を表している。つまり、ＭＵＳＩＣスペクトラムの値の最大化を目的とするが、計算の都合上、逆数を用いて最小化を目指す。
【００５２】
【数１】
ここで、Ｎは到来波の数、θｍはｍ番目の波の到来方向、ａ（θｍ）はθｍ方
【外１】
前）の相関行列の雑音固有ベクトルから構成される行列である。
【００５３】
また、δ＝［Γ１１，Γ２２，…，Γｍｍ］Ｔ，Ｑ１（ｍ）＝ｄｉａｇ｛ａ （θｍ）｝である。但し、Ｔは転置。
【００５４】
このとき、Ｊｃを最小化するΓｉｉは以下の式で与えられる。
【００５５】
【数２】
ＭＵＳＩＣのスペクトラムを鋭くするようにする補正量は、実際の補正量を表すものとなる。そこで、この補正量を到来方向推定手段５に反映することで、最適な到来方向を推定することができる。即ち、到来方向推定手段５は補正量算出手段４で算出された補正量を受け取ると、その逆数をＡ／Ｄ変換器３ａ、…、３ｎの出力に掛けて、アンテナの通過特性であるΓｉｉや他の受信アナログ部２との相互結合係数Ｃｉｊにより生じる受信誤差を解消して、入力受信レベルを対応するアンテナ素子の本来の受信レベルと同じにする。
【００５６】
以上説明したように、本実施形態の到来波方向推定装置では、複数の電波の受信信号のみを用いて、ＭＵＳＩＣアルゴリズムで得られる到来方向情報を利用し、補正量を算出（つまりキャリブレーション）できる。ここで、複数の電波の受信信号は、キャリブレーション用として特別に用意する必要はなく、実使用中の受信信号を用いるだけでよい。
【００５７】
また、振幅位相、そして相互結合までもが、リアルタイムにキャリブレーションできるようになる。更に、新たなアナログ部品を必要とすることもなく、装置を安価に構成することができる。
【００５８】
次に補正量算出手段４における第２の実施例の算出方法について図４のフローチャートを参照して説明する。始めに、振幅位相Γｉｉと相互結合係数Ｃｉｊに初期値を与える（ステップ４０１）。初期値としては、前回のキャリブレーションデータが存在するならばその値を、まったく未知な場合には、偏差や相互結合がない条件で与えればよい。次に、振幅位相偏差と相互結合係数を固定して到来方向をＭＵＳＩＣアルゴリズムで推定する（ステップ４０２）。ここでは、振幅位相Γｉｉと相互結合係数Ｃｉｊは固定である。
【００５９】
次に到来方向と相互結合係数を固定して振幅位相偏差Γｉｉを推定する（ステップ４０３）。ここでは、到来方向と相互結合係数は固定である。Γｉｉは、ＭＵＳＩＣスペクトラムのピークが到来方向においてより鋭くなる（大きな値になる）ように設定される。
【００６０】
更に、到来方向と振幅位相偏差を固定して相互結合係数Ｃｉｊを推定する（ステップ４０４）。ここでは、到来方向と振幅位相は固定である。Ｃｉｊは、ＭＵＳＩＣスペクトラムのピークが到来方向において鋭くなる（大きな値になる）ように設定される。このステップ４０２からステップ４０４の処理を繰り返し計算して、推定結果が収束して未知変数の値が変化しなくなったら計算を終了し（ステップ４０５）、補正量を到来方向推定手段５に出力する（ステップ４０６）。尚、ステップ４０３と４０４の処理における計算方法は、いかなる手法を用いても良い。
【００６１】
この実施例は、効率的に補正量を算出する方法について記述している。未知パラメータとしては、電波の到来方向、振幅位相偏差、相互結合係数の３種類である。これらを同時に推定することは未知パラメータの数が膨大になり、リアルタイムにキャリブレーションを行う場合に問題になる。そこで、これら３種類の未知変数を同時に算出するのではなく、順番に算出し、結果が収束するまで繰り返し求めていくものである。
【００６２】
このアルゴリズムを用いた場合の補正量算出結果の一例を示す。アンテナは６素子の円形アレイで、誤差として図５に示す振幅位相を設定した。このときの上記アルゴリズムを用い、ステップ４０２からステップ４０４の処理を繰り返し行った結果を図５に示す。尚、振幅位相の初期値としては、誤差のない値を用いている。また、振幅位相を求める計算は、第１の実施例で記述した方法で行った。図５より、初期値として誤差のない条件を用いているにも拘らず、補正量が計算できていることがわかる。値としては、１００％真の値にはなっていないが、十分キャリブレーションデータとして使用できる値になっている。
【００６３】
以上のように、本実施例によれば、未知変数を繰り返し計算で求めるため、効率的に計算することができるようになり、リアルタイムでのキャリブレーションを行うに際して高価な高速の計算機を用いなくともよく、装置を安価に構成することに大いに貢献することとなる。
【００６４】
次に補正量算出手段の第３の実施例の算出方法について図６のフローチャートを参照して説明する。始めに、振幅位相Γｉｉと相互結合係数Ｃｉｊに初期値を与える（ステップ６０１）。初期値としては、前回のキャリブレーションデータが存在するならばその値を、まったく未知な場合には、偏差や相互結合がない条件で与えればよい。次に、振幅位相偏差と相互結合係数を固定して到来方向をＭＵＳＩＣアルゴリズムで推定する（ステップ６０２）。ここでは、振幅位相Γｉｉと相互結合係数Ｃｉｊは固定である。次に到来方向と相互結合係数を固定して振幅位相偏差Γｉｉを推定する（ステップ６０３）。ここでは、到来方向と相互結合係数は固定である。Γｉｉは、ＭＵＳＩＣスペクトラムのピークが到来方向においてより鋭くなる（大きな値になる）ように設定される。次に、振幅位相偏差と相互結合係数を固定して到来方向をＭＵＳＩＣアルゴリズムで推定する（ステップ６０４）。ここでは、振幅位相Γｉｉと相互結合係数Ｃｉｊは固定である。
【００６５】
更に、到来方向と振幅位相偏差を固定して相互結合係数Ｃｉｊを推定する（ステップ６０５）。ここでは、到来方向と振幅位相は固定である。Ｃｉｊは、ＭＵＳＩＣスペクトラムのピークが到来方向において鋭くなる（大きな値になる）ように設定される。このステップ６０２からステップ６０５の処理を繰り返し計算して、推定結果が収束して未知変数の値が変化しなくなったら計算を終了し（ステップ６０６）、補正量を到来方向推定手段５に出力する（ステップ６０７）。
【００６６】
本実施例では、図４で示した第２の実施例に対して、計算順序を変え、特に到来波推定方向の算出をステップ６０２とステップ６０４の処理で２回行っている。
【００６７】
本実施例においては、図６に示したように、到来方向の推定→振幅位相算出→到来方向の推定→相互結合係数算出→…としている。このようにすると、到来方向推定のステップがひとつ増えて計算量が増加するように見えるが、繰り返し回数が減少する効果が期待できる。以下にその理由を説明する。
【００６８】
図４では、ステップ４０３で補正量が算出されるが、補正後のデータを用いて到来方向を推定し直すと、方向が変化している場合がある。つまり、図４に示したようなステップを踏むと、ステップ４０２で推定した到来方向と、ステップ４０３で補正したデータを用いて改めて推定した到来方向は異なっている場合がある。ここで、このままステップ４０４で補正を行うが、ステップ４０２の到来方向のデータと、ステップ４０３の補正データは整合が取れていないために、ステップ４０４で求まる補正量には誤差が含まれている場合がある。
【００６９】
そこで、本実施例では、振幅位相の補正量を算出した後に、ステップ６０４で到来方向を推定する。そして、相互結合係数の算出を行う。その結果、相互結合係数の算出においては、誤差を含まない形で求めることができるので、繰り返し回数を少なくして補正量を求めることができるようになる。
【００７０】
以上説明したように、本実施例によれば、補正量を算出した後に毎回到来方向を推定し直すので、誤差のない補正量算出が実現でき、その結果、少ない繰り返し回数で算出結果が収束し、リアルタイムでのキャリブレーションをより早く行うことができる。
【００７１】
また、振幅位相補正量が相互結合補正量に比ベて明らかに大きい場合には、ステップ６０２、６０３、６０２、６０３、６０２、６０４、６０２、６０３、６０２、６０３、６０２、６０４、…、と、振幅位相の補正量を求める回数を増やすことで、更に繰り返し回数を減少させる効果も期待できる。
【００７２】
図７は本発明の到来波方向推定装置の第２の実施形態の構成を示した図である。本例の到来波方向推定装置１００は、図７に示すように、補正量算出のためのひとつ或いは複数の電波の送信装置９ａ、９ｂを有している。また、前記電波の送信装置９ａ、９ｂは、前記到来波推定装置１００に対しての相対的な方向が既知となっている。尚、本例の到来波方向推定装置１００は図１に示した第１の実施形態と同様であるため、以降、図１を借用して説明する。
【００７３】
図８は図７に示した到来波方向推定装置の補正量算出手段４の補正量算出方法の実施例を示したフローチャートである。本例の補正量算出手段４は、図７に示すように既知の方向から到来する送信装置９ａ、９ｂから送信された電波の受信信号を用いて、振幅位相Γｉｉと相互結合係数Ｃｉｊを繰り返し計算して補正量を算出する。
【００７４】
始めに、振幅位相Γｉｉと相互結合係数Ｃｉｊに初期値を与える（ステップ８０１）。初期値としては、前回のキャリブレーションデータが存在するならばその値を、まったく未知な場合には、偏差や相互結合がない条件で与えればよい。
【００７５】
次に到来方向（これは既知である）と相互結合係数を固定して振幅位相偏差Γｉｉを推定する（ステップ８０２）。ここでは、到来方向と相互結合係数は固定である。ΓｉｉはＭＵＳＩＣスペクトラムのピークが到来方向においてより鋭くなる（大きな値になる）ように設定される。
【００７６】
更に、到来方向と振幅位相偏差を固定して相互結合係数Ｃｉｊを推定する（ステップ８０３）。ここでは、到来方向と振幅位相は固定である。Ｃｉｊは、ＭＵＳＩＣスペクトラムのピークが到来方向において鋭くなる（大きな値になる）ように設定される。このステップ８０２からステップ８０３の処理を繰り返し計算して、推定結果が収束して未知変数の値が変化しなくなったら計算を終了し（ステップ８０４）、補正量を到来方向推定手段５に出力する（ステップ８０５）。
【００７７】
本実施例では、効率的に補正量を算出する方法について記述している。本例では、既知の方向の波源からの電波の受信信号を用いているため、未知パラメータとしては、振幅位相偏差、相互結合係数の２種類だけである。つまり、未知パラメータであった到来方向が既知の条件となり、この条件で補正量を推定する。従って、方向推定の計算がまつたく不要になるので処理量も減少する。また、方向が既知であるので、最終的に得られる補正量の精度の向上が期待できる。
【００７８】
図９に方向が既知とした場合の補正量の算出結果を示す。シミュレーション条件は図５と同じである。図５に示した未知の方向からの受信信号を用いた場合に比べて、精度良く補正量を算出できていることが分かる。
【００７９】
尚、本例は、ひとつのみの送信装置を用いても、補正量は算出できる。しかし、複数の送信装置を用いることで、送信装置の位置ずれ、送信装置と到来波推定装置の間の伝搬環境の変化、送信装置の故障、などに対する対応が可能となり、補正量を安定且つ更に精度良く算出することができる。
【００８０】
以上説明したように、本実施形態によれば、補正量を算出するための電波の送信装置を用いて補正量を算出することで、処理量を大幅に少なくでき、また、精度の高い補正量が得られるキャリブレーションをリアルタイムに行うことができる。
【００８１】
図１０は本発明の到来波方向推定装置の第３の実施形態の構成を示したブロック図である。但し、図１に示した第１の実施形態と同様の部分には同一符号を用い、且つ適宜その説明を省略する。到来波方向推定装置は、複数のアンテナ素子１ａ、…、１ｎと、複数のアンテナ素子１ａ、…、１ｎに対応して設けられ、対応するアンテナ素子より伝達される受信信号を処理する複数の受信アナログ部２ａ、…、２ｎと、受信アナログ部２ａ、…、２ｎの出力信号をデジタル変換するＡ／Ｄ変換器３ａ、…、３ｎと、複数の受信アナログ部２ａ、…、２ｎの相対的な偏差を推定する補正量算出手段４と、補正量算出手段４の出力する補正量を用いて複数のアンテナ素子１ａ、…、１ｎに到来するひとつ或いは複数の電波の到来方向を推定する到来方向推定手段５及び補正量が算出される毎に、新たな補正量を記憶して記憶内容を更新する補正量記憶装置６を有している。
【００８２】
次に本実施形態の動作について説明する。本例の到来波方向推定動作、特に補正量算出手段４による補正量の算出動作は第１の実施形態のそれと同様である。異なる動作は、補正量算出手段４により算出された補正量を補正量記憶装置６によって記憶しているところにある。
【００８３】
本実施形態でも、補正量はリアルタイムに補正することが可能であるが、短時間の変動が小さい場合には、ある一定の時間間隔で補正すれば十分である。この場合には、補正量記憶装置６に記憶された補正量を用いて到来波方向の推定を間欠的に行えばよく、装置及び本装置を含むシステムの負荷を減らすことができる。また、補正量を計算する時の初期値として、補正量記憶装置６内の補正量を用いれば、補正量算出アルゴリズムの収束時間を短くすることができる。
【００８４】
以下、上記した補正量記憶装置６内の補正量を用いて、到来波の推定を行う動作の実施例について説明する。
【００８５】
本例の補正量算出手段４は、補正量記億装置６に記憶された過去の補正量を利用して、次の補正量を予測する。
【００８６】
ここで、補正量をリアルタイムに算出することは、精度を向上させる上では重要である。しかし、補正量の算出回数を減らし、その空いた時間には、補正量を算出するための計算リソースを到来方向推定アルゴリズムに利用し、方向推定に係る計算時間を短くすることはリアルタイム処理にとっては非常に重要である。
【００８７】
そこで、本例では、補正量記憶装置６に記憶された過去の補正量から現時点での補正量を予測し、その補正量を用いて電波の到来方向を推定する。
【００８８】
補正量を予測する方法の具体例を以下に示す。第１の例として補正量を線形補間で予測する。図１１は振幅位相の振幅を線形補間で予測する方法を示した例である。横軸は補正回数、縦軸は補正振幅量を表している。過去の２回の振幅補正量から現在の補正量を一次の線形補間で予測する。ここで、過去の３回およびそれ以上のデータを用いて２次の近似から予測することもできる。
【００８９】
第２の例として、補正量算出手段４によって算出された補正量と、そのときの外気気温の両方を補正量記憶装置６に記憶することで予測することもできる。受信アナログ部２の偏差が外気気温に相関がある場合がある。つまり、外気気温が分かっていれば、受信アナログ部２の偏差が予測できる場合がある。このような場合には、外気気温をモニターすることで、補正量を予測することもできる。または、受信アナログ部２が外気気温とは関わりがなく、受信アナログ部２を含んでいる装置内の気温と関わりがある場合がある。つまり、装置内の機器の発生する熱によって外気と装置内の温度差が異なる場合がある。この場合には、装置内の気温と補正量を補正量記憶装置６に記憶することが有効である。
【００９０】
第３の例として、補正量算出手段４によって算出された補正量と、そのときの時刻の両方を補正量記憶装置６に記憶する。例えば、構内無線や、移動体通信の基地局では、使用時間と使用頻度の間には相関がある。つまり、工場内の構内無線では、昼間に使用される頻度が高く、屋外の基地局では、社会人の仕事が終わった時間帯に使用される頻度が高い。使用する頻度が高いと、機器内の発生する熱量も多くなり、装置内の温度が上昇する。つまり、時刻と気温の間にも相関がある。そして、先に説明したように、気温と補正量の間にも相関がある。つまり、時刻と気温の間に相関がある。つまり、時刻と補正量を補正量記憶装置６に記憶しておくことは有効なことである。また、時刻によって外気の気温に変化がある。また、季節によって外気の差がある。そこで、時刻のほかにも、日時も記憶することは有効なことである。
【００９１】
また、記憶された離散的な時刻や日時と補正量の値から、線形補間によって記憶されていない時刻や日時における補正量を予測することも有効である。
【００９２】
同様に、記憶された離散的な気温と補正量の値から、線形補間によって記憶されていない気温における補正量を予測することも有効である。
【００９３】
更に、上記した気温、装置内部の温度、時刻や日時により、補正量算出時の初期値の設定をより適切なものにすることもできる。
【００９４】
尚、上記の各実施形態では、到来方向推定装置での説明であったが、同様に、アダプティブアンテナ装置の受信系の補正としてもまったく同様に用いることができる。
【００９５】
また、上記の各実施形態で得られる補正量と受信で用いたアダプティブアンテナ装置の重みから、送信に用いるべき重みを逆算することも勿論可能である。このときには、別に用意した送信用のキャリブレーションデータを必要とする。特にＴＤＭＡ／ＴＤＤシステムでは、送受で同一のパターンが要求されるので、有効である。以上の説明では、相互結合係数の補正量を算出していたが、受信アナログ部のシールドが精度良く構成されている場合には、相互結合係数を無視することができる。このような場合には、振幅位相の補正量のみを算出すればよい。また、以上の説明に用いた複数の実施形態を組み合わせて到来方向推定装置を構成することも勿論、可能である。
【００９６】
更に、本発明は上記実施形態に限定されることなく、その要旨を逸脱しない範囲において、具体的な構成、機能、作用、効果において、他の種々の形態によっても実施することができる。
【００９７】
【発明の効果】
以上詳細に説明したように、本発明によれば、運用中の受信データを用いてリアルタイムに高精度なキャリブレーションを行うことができ、到来波方向推定精度を向上させることができる。また、キャリプレーション用の送信装置を用いることで、さらに、精度の良い補正を行うことができるようになる。また、過去の補正量を記憶しておくことで、現時点での補正量を線形補間などで簡易に予測することができる。
【図面の簡単な説明】
【図１】本発明の到来波方向推定装置の第１の実施形態に係る構成を示したブロック図である。
【図２】図１に示した受信アナログ部の偏差を回路モデルで表した等価回路図である。
【図３】図１に示した補正量算出手段で用いるＭＵＳＩＣスペクトラムの特性を示した特性図である。
【図４】図１に示した補正量算出手段における補正量算出方法の第２の実施例を示したフローチャートである。
【図５】図１に示した補正量算出手段により算出される振幅位相Γｉｉの設定値と計算値を示した表図である。
【図６】図１に示した補正量算出手段における補正量算出方法の第３の実施例を示したフローチャートである。
【図７】本発明の到来波方向推定装置の第２の実施形態に係る構成を示したブロック図である。
【図８】図７に示した補正量算出手段における補正量算出方法の実施例を示したフローチャートである。
【図９】図７に示した補正量算出手段により算出される振幅位相Γｉｉの設定値と計算値を示した表図である。
【図１０】本発明の到来波方向推定装置の第３の実施形態に係る構成を示したブロック図である。
【図１１】図１０の補正量算出手段により求められる振幅量の線形補間方法を説明する図である。
【図１２】従来のアダプティブアンテナ装置の構成例を示したブロック図である。
【図１３】図１２の装置における送信パターンと受信パターン例を示した図である。
【図１４】従来の到来波方向推定装置の構成例を示したブロック図である。
【符号の説明】
１ａ、１ｎ アンテナ素子
２ａ、２ｎ 受信アナログ部
３ａ、３ｎ Ａ／Ｄ変換器
４ 補正量算出手段
５ 到来方向推定手段
６ 補正量記憶装置
９ａ、９ｎ 送信装置
２１ 低雑音増幅器
２２ フィルタ
２３ 周波数変換器
１００ 到来波方向推定装置
【発明の属する技術分野】
本発明は、移動通信や構内無線やレーダ装置などに利用される到来波方向推定装置及び到来波方向推定方法に関する。
【０００２】
【従来の技術】
近年、送受共用のアダプティブアンテナ装置に対する需要が高まっており、例えばＴＤＭＡ／ＴＤＤの通信システムにおいては、送受の周波数が一致することから、送受で同一のパターンを形成するアダプティブアンテナ装置を利用することが有効である。
【０００３】
図１２は従来の送受共用のアダプティブアンテナ装置の構成例を示したブロック図である。アダプティブアンテナを形成する複数のアンテナ素子３１ａ、…、３１ｎがあり、個々のアンテナ素子には信号を受信する受信アナログ部３２ａ、…、３２ｎと、信号を送信する送信アナログ部３３ａ、…、３３ｎが切り替えで接続されている。
【０００４】
受信アナログ部３２ａ、…、３２ｎの受信信号は受信パターン形成回路３４に入力されて最終的な受信信号になる。送信信号は送信パターン形成回路３５に入力され、個々の送信アナログ部３３ａ、…、３３ｎに分配される。
【０００５】
このような送受信共用のアダプティブアンテナ装置は、個々のアンテナ素子３１ａ、…、３１ｎに対する受信信号に基づいて、個々のアンテナ素子３１ａ、…、３１ｎに対応する複素振幅の重み付けを受信パターン形成回路３４で決定し、この重みを共通のアンテナ素子に対応する受信信号に付与し、或いは前記重み付けを送信パターン形成回路３５で送信信号に付与することで、所望のアンテナパターンを形成している。従って原理的には、受信と送信で同一のパターンを形成することができることになる。
【０００６】
しかしながら、受信と送信で同一のパターンを形成することは現実的には困難である。その理由は、多くの場合、アダプティブアンテナ内のアンテナ素子毎の送信アナログ部３３、受信アナログ部３２の通過特性（通過振幅特性、通過位相特性）に偏差（ばらつき）があることに起因している。
【０００７】
このような偏差があると、図１３に示すように、受信のパターン５０と送信パターン６０が異なり、所望方向への利得が低くなったり、干渉ユーザーへの利得が高くなったりしてしまう。
【０００８】
或いは、受信アダプティブアンテナ装置において、方向拘束付のアダプティブアルゴリズムを用いる場合がある。この場合にも、受信アナログ部の通過特性に偏差（ばらつき）があると、最適な解を得ることができなくなる。
【０００９】
更に、受信アダプティブアンテナ装置はレーダ装置の一部として、電波の到来方向を推定することに使用されるが、受信アナログ部の通過特性に偏差（ばらつき）があると、正確に電波の到来方向を推定することができなくなる。
【００１０】
上記のような偏差を取り除く方法として、一般的には、キャリブレーションが行われる。キャリブレーションとは、アダプティブアンテナ装置の製作が終了した時点で、すべての受信アナログ部と送信アナログ部の偏差を測定し、補正量 （偏差量に同じ）を予め算出しておく。そして、実際にアダプティブアンテナ装置を運用するときには、補正量を考慮して各素子に設定する振幅位相（重み）を算出すればよい。
【００１１】
しかしながら、最初に求めた補正量が運用中に変化してしまう場合がある。例えば、昼夜の外気温度変化の違いによるアナログ部の部品（増幅器、周波数変換機、フィルタなど）の特性変化、そして、それらを接続するケーブルの長さの伸縮が考えられる。このような変化が生じると、受信アダプティブアンテナ装置はその性能が大幅に劣化して到来方向を正確に推定できなくなり、更に、これにより加入者容量が大幅に劣化して、最悪的には、通信の切断にまで及ぶ場合がある。つまり、キャリブレーションは一度では不十分であり、運用中にリアルタイムにキャリブレーションすることが要求される。
【００１２】
このような問題を解決する方法として、図１４は特開平２０００−９１８３３の発明例を示したもので、アンテナ素子４１ａ、４１ｂ、４１ｃで受信した受信信号はビームフォーマー４２を通ってビームスペースデータに変換される。ビームフォーマーの出力（ビームスペースデータ）は、レベル調整器４４ｂ，４４ｃにより補正される。補正されたデータは、逆変換回路によって再び、アンテナ素子で受信した受信信号（これは補正されている）に変換される。
【００１３】
位相補正量演算部４３は各、ビームフォーマーの出力間の相対的な位相差を検出する。ここでは、到来波がひとつであり、アンテナ素子４１ａ、４１ｂ、…、４１ｃの配置位置が決まっていれば、ビームフォーマーの出力間の位相差は到来角度によって決まってくる。そこで、位相補正量演算部４３では、検出した位相差と本来得られる値との差を補正量として算出し、その補正量によりアンテナ素子４１ｂ、４１ｃの受信レベルをレベル調整器４４ｂ、４４ｃにより調整して、到来波の推定方向の精度を高めることができる。しかも、この従来例では、運用中に上記の補正をリアルタイムに行うことができるメリットがある。
【００１４】
しかしながら、到来波の数は、ひとつに限定され、また、位相のみの補正になっているので、振幅の変動があった場合には、偏差を保証することができない。特に、アダプティブアンテナ装置では、複数の到来波を受信している状態でパターン形成を行うものである。到来してくる電波の数をひとつにしなければいけないので、通常の通信中にキャリブレーションを行うことができない。
【００１５】
以上説明したように、従来のアダプティブアンテナ装置に対して発明されたキャリブレーション方法は、リアルタイムに行えるものであるが、到来波の数がひとつに限定され、位相しか補正できないという問題点があった。
【００１６】
また、以上のアダプティブアンテナ装置に関する従来の問題点は、受信アダプティブ装置と同様の構成を一部に持つアンテナレーダ装置における到来方向の推定精度を向上させる場合においても、まったく同様な問題点となる。
【００１７】
【発明が解決しようとする課題】
以上説明したように従来の到来波方向推定装置では、回路の特性偏差を取り除く方法として、一般的にキャリブレーションが行われるが、最初に求めた補正量が運用中に変化してしまう場合があるため、リアルタイムにキャリブレーションを行なわなければならない。しかし、その場合、到来波はひとつに限定され、しかも位相のみの補正しかできず、しかも複数の到来波がある運用中にキャリブレーションを行うことができないため、実用上不十分なものであった。
【００１８】
本発明は、上述の如き従来の課題を解決するためになされたもので、運用中の受信データを用いてリアルタイムに高精度なキャリブレーションを行うことができ、到来波方向推定精度を向上させることができる到来波方向推定装置及び到来波方向推定方法を提供することを目的とする。
【００１９】
【課題を解決するための手段】
上記目的を達成するために、本発明の特徴は、到来波方向推定装置において、複数のアンテナ素子と、前記複数のアンテナ素子に対応して設けられ、対応するアンテナ素子より伝達される受信信号を前記アンテナ素子対応で処理する複数の受信部と、前記複数の受信部から出力される受信信号から前記複数の受信部から出力される信号間の相対的な偏差を推定する補正量算出手段と、前記複数の受信部から出力される受信信号を前記補正量算出手段により得られる偏差量を用いて補正し、これら補正された受信信号を用いて前記複数のアンテナ素子に到来する少なくともひとつの電波の到来方向を推定する到来方向推定手段と、を具備し、前記補正量算出手段は、ひとつ或いは複数の電波の到来方向と、複数の受信部から出力される信号間の振幅位相偏差と、複数の受信部から出力される信号間の相互結合係数の３種類の未知変数を求める処理を行い、ステップ１として振幅位相偏差と相互結合係数を固定して到来方向を推定し、ステップ２として到来方向と相互結合係数を固定して振幅位相偏差を推定し、ステップ３として到来方向と振幅位相偏差を固定して相互結合係数を推定し、前記ステップ１からステップ３の処理を推定結果が収束するまで繰り返し実行することにより、前記偏差量を求めることにある。
【００２０】
本発明の他の特徴は、到来波方向推定装置において、複数のアンテナ素子と、前記複数のアンテナ素子に対応して設けられ、対応するアンテナ素子より伝達される受信信号を前記アンテナ素子対応で処理する複数の受信部と、前記複数の受信部から出力される受信信号から前記複数の受信部から出力される信号間の相対的な偏差を推定する補正量算出手段と、前記複数の受信部から出力される受信信号を前記補正量算出手段により得られる偏差量を用いて補正し、これら補正された受信信号を用いて前記複数のアンテナ素子に到来する少なくともひとつの電波の到来方向を推定する到来方向推定手段と、を具備し、前記補正量算出手段は、ステップ１として振幅位相偏差と相互結合係数を固定して到来方向を推定し、ステップ２として到来方向と相互結合係数を固定して振幅位相偏差を推定し、ステップ３として振幅位相偏差と相互結合係数を固定して到来方向を推定し、ステップ４として到来方向と振幅位相偏差を固定して相互結合係数を推定し、前記ステップ１からステップ４の処理を推定結果が収束するまで繰り返し実行することにより、前記偏差量を求めることにある。
【００２１】
本発明の他の特徴は、到来波方向推定装置において、複数のアンテナ素子と、前記複数のアンテナ素子に対応して設けられ、対応するアンテナ素子より伝達される受信信号を前記アンテナ素子対応で処理する複数の受信部と、前記複数の受信部から出力される受信信号から前記複数の受信部から出力される信号間の相対的な偏差を推定する補正量算出手段と、前記複数の受信部から出力される受信信号を前記補正量算出手段により得られる偏差量を用いて補正し、これら補正された受信信号を用いて前記複数のアンテナ素子に到来する少なくともひとつの電波の到来方向を推定する到来方向推定手段と、を具備し、前記到来波方向推定装置は、偏差量算出のためのひとつ或いは複数の電波の送信装置を有しており、前記電波の送信装置の電波発信位置は、前記到来波推定装置に対して相対的な方向が既知であり、前記補正量算出手段は、前記送信装置から発せられた既知の方向から到来する電波の受信信号を用いて、ステップ１として到来方向と相互結合係数を固定して振幅位相偏差を推定し、ステップ２として到来方向と振幅位相偏差を固定して相互結合係数を推定し、前記ステップ１とステップ２の処理を推定結果が収束するまで繰り返し実行することにより、前記偏差量を求めることにある。
【００２２】
本発明の他の特徴は、到来波方向推定装置において、複数のアンテナ素子と、前記複数のアンテナ素子に対応して設けられ、対応するアンテナ素子より伝達される受信信号を前記アンテナ素子対応で処理する複数の受信部と、前記複数の受信部から出力される受信信号から前記複数の受信部から出力される信号間の相対的な偏差を推定する補正量算出手段と、前記複数の受信部から出力される受信信号を前記補正量算出手段により得られる偏差量を用いて補正し、これら補正された受信信号を用いて前記複数のアンテナ素子に到来する少なくともひとつの電波の到来方向を推定する到来方向推定手段と、を具備し、前記補正量算出手段は補正量記憶装置を有しており、偏差量を算出する毎に得られた偏差量を補正量記憶装置に記憶して記憶内容を更新し、前記補正量記億装置は、前記補正量算出装置によって算出された偏差量と、そのときの外気気温の両方を記憶することにある。
【００２３】
本発明の他の特徴は、到来波方向推定装置において、複数のアンテナ素子と、前記複数のアンテナ素子に対応して設けられ、対応するアンテナ素子より伝達される受信信号を前記アンテナ素子対応で処理する複数の受信部と、前記複数の受信部から出力される受信信号から前記複数の受信部から出力される信号間の相対的な偏差を推定する補正量算出手段と、前記複数の受信部から出力される受信信号を前記補正量算出手段により得られる偏差量を用いて補正し、これら補正された受信信号を用いて前記複数のアンテナ素子に到来する少なくともひとつの電波の到来方向を推定する到来方向推定手段と、を具備し、前記補正量算出手段は補正量記憶装置を有しており、偏差量を算出する毎に得られた偏差量を補正量記憶装置に記憶して記憶内容を更新し、前記補正量記憶装置は、前記補正量算出装置によって算出された偏差量と、そのときの受信アナログ部の周りの気温の両方を記憶することにある。
【００２４】
本発明の他の特徴は、到来波方向推定装置において、複数のアンテナ素子と、前記複数のアンテナ素子に対応して設けられ、対応するアンテナ素子より伝達される受信信号を前記アンテナ素子対応で処理する複数の受信部と、前記複数の受信部から出力される受信信号から前記複数の受信部から出力される信号間の相対的な偏差を推定する補正量算出手段と、前記複数の受信部から出力される受信信号を前記補正量算出手段により得られる偏差量を用いて補正し、これら補正された受信信号を用いて前記複数のアンテナ素子に到来する少なくともひとつの電波の到来方向を推定する到来方向推定手段と、を具備し、前記補正量算出手段は補正量記憶装置を有しており、偏差量を算出する毎に得られた偏差量を補正量記憶装置に記憶して記憶内容を更新し、前記補正量記憶装置は、前記補正量算出手段によって算出された偏差量と、そのときの時刻、日時の両方を記憶することにある。
【００３４】
【発明の実施の形態】
以下、本発明の実施形態を図面に基づいて説明する。図１は、本発明の到来波方向推定装置の第１の実施形態に係る構成を示したブロック図である。到来波方向推定装置は、複数のアンテナ素子１ａ、…、１ｎと、複数のアンテナ素子１ａ、…、１ｎに対応して設けられ、対応するアンテナ素子より伝達される受信信号を処理する複数の受信アナログ部２ａ、…、２ｎと、受信アナログ部２ａ、…、２ｎの出力信号をデジタル変換するＡ／Ｄ変換器３ａ、…、３ｎと、複数の受信アナログ部２ａ、…、２ｎの相対的な偏差を推定する補正量算出手段４と、補正量算出手段４の出力する補正量を用いて複数のアンテナ素子１ａ、…、１ｎに到来するひとつ或いは複数の電波の到来方向を推定する到来方向推定手段５とから構成される。
【００３５】
ここで、受信アナログ部２（受信アナログ部２ａ、…、２ｎの総称）は、所望信号を増幅する低雑音増幅器２１と、所望の周波数帯域のみを抽出するフィルタ２２と、周波数帯をＲＦ（電波周波数）からＢＢ（ベースバンド）又はＩＦ（中間周波数）へ変換するための周波数変換器２３とから構成されている。この受信アナログ部２は、受信信号を抽出し、周波数変換し、増幅する機能を有するものであれば、他のいかなる構成を有するものに置き換えて構わない。
【００３６】
次に本実施形態の動作について説明する。補正量計算手段４は、ひとつ或いは複数の到来してくる電波の受信信号を用い、ＭＵＳＩＣ（Multiple Signal Classification）アルゴリズムで推定される複数の電波の到来方向情報を利用して補正量を算出する第１の実施例の方法を採用している。
【００３７】
ここで、複数のアンテナ素子１ａ、１ｎはどのようなものであっても良い。例えば、ダイポールアンテナを用いれば、水平面内で概略無指向性が得られるため、ひとつの到来波方向推定装置で３６０度をカバーできる。例えば、ホーンアンテナを用いれば各アンテナは指向性を有しており、限定した領域のみを到来波方向推定手段５でカバーすることになる。こうすると、領域が限定されるものの、後方からの到来波の影響を受けずに電波の方向を推定することができる。ここで、アンテナの放射パターンは予め測定しておき、到来方向推定手段５と補正量算出手段４において使用するものとする。
【００３８】
次に受信アナログ部２であるが、ここでは、受信信号が低雑音増幅器２１で増幅され、フィルタ２２により所望の周波数帯域のみが抽出され、周波数変換器２３により周波数がベースバンド（ＢＢ）または中間周波（ＩＦ）へ変換される。ここで、各受信アナログ部２ａ、…、２ｎは、偏差（ばらつき）を有しており、ブランチ毎に異なった通過特性を有することとなる。また、受信アナログ部２ａ、…、２ｎ間には、シールドが完全でない場合に相互結合が生じることがある。つまり、ある受信アナログ部で受信した信号が、他の受信アナログ部へ漏れこむ現象が発生する。
【００３９】
次に各受信信号は、各Ａ／Ｄ変換器３ａ、…、３ｎによってデジタル信号に変換され、到来方向推定手段５において処理され、到来波の方向が推定される。
【００４０】
一方、補正量算出手段４は、各Ａ／Ｄ変換器３ａ、…、３ｎによってデジタル信号に変換された受信信号を用いて補正量を算出する。この補正量が正確に求まったならば、どのような偏差があったとしても（ケーブルが切断され振幅が０となった場合を除く）、それを補正できるものであり、精度良く到来波方向を推定できるようになる。また、到来波方向推定手段５は、どのような方法によって方向を推定するものでも構わない。ＭＵＳＩＣアルゴリズムを用いても、ＥＳＰＲＩＴ（Estimation of Signal Parameters via Rotational In rariance Techniques）アルゴリズムを用いても、モノパルス側角であっても、ビームフォーマー法でも構わない。
【００４１】
図２は受信アナログ部２の偏差を回路モデルで表した図である。複数のアンテナ１ａ、…、１ｎで受信した信号は、Γｉｉ（ｉ＝１〜ｎ）（ｎは素子数）倍される。ここでΓｉｉは通過特性をあらわす複素数であり、振幅と位相の偏差を意味している。つまり、全てのΓｉｉがｌｅｘｐ（ｊ０゜）ならば、通過特性はすべての受信アナログ部２ａ、…、２ｎで同一であることを意味している。一方、Γｉｉが受信アナログ部２ａ、…、２ｎによって異なる場合には、通過特性を補正する必要がある。
【００４２】
次に、Γｉｉ倍された受信信号は、Ｃｉｊ（ｉ＝１〜ｎ，ｊ＝１〜ｎ）倍されて、他の受信アナログ部２に漏れこむ。Ｃｉｊとは受信アナログ部ｉ番目とｊ番目の相互結合係数である。ここで、Ｃｉｉは常に１である。また、Ｃｉｊ（ｉ≠ｊ）＝０ならば、結合は発生していない状態を意味している。このように、受信アナログ部２ａ、…、２ｎでは、振幅位相偏差と相互結合の二つの補正をしなければならない未知変数がある。
【００４３】
尚、振幅位相偏差と相互結合の２つは受信アナログ部２ａ、…、２ｎのどの部分に発生しているかは問題としない。つまり、受信アナログ部２の全体の偏差を求めることが重要であるからである。また、振幅位相偏差と相互結合の回路モデルの順番も問題ではない。なぜならば、順番を逆にしたとしても、受信アナログ部２の全体の偏差としては等価なものが得られるからである。
【００４４】
次に補正量算出手段４であるが、ひとつ或いは複数の到来してくる電波の受信信号を用いるものである。つまり、任意の到来電波を使用することができるので、アダプティブアンテナ装置を運用中にリアルタイムにキャリブレーションを行うことができ、また、キャリブレーション用の特別な送信源を必要としない。そして、ＭＵＳＩＣアルゴリズムで推定される複数の電波の到来方向情報を利用して補正量を算出する。
【００４５】
一般的に、キャリブレーションが正確に行われていたならば、ＭＵＳＩＣアルゴリズムによって推定される到来方向においては、図３の実線に示すように、鋭いピークを有するものとなる。なお、ＭＵＳＩＣスペクトラムは、アンテナの放射パターンと、受信信号と、補正量（キャリブレーションデータ）から算出されるものである。しかしながら、偏差がある場合には、図３の破線に示すように推定方向のピークが鈍り、推定精度が劣化してしまう。
【００４６】
ここでの補正量算出手段４は、ＭＵＳＩＣによって得られるＭＵＳＩＣスペクトラムのピークを鋭くするように、補正量（図２におけるΓｉｉとＣｉｊ）を算出するものである。算出する方法は、いかなる方法によっても良い。例えば、考えうる全ての組み合わせの補正量からもっともスペクトラムが鋭くなるものを選んでもよい。または、ジェネティックアルゴリズム（遺伝的手法）によっても良いし、ニューロを用いた方法でも良い。以下に、Γを求めるための計算方法の一例を示す。なお、この計算方法は文献（B.Friedlander,et al.,"Direction finding in the presence of mutual coupling,"IEEE Trans.Antenna and Propagation,vol.39,no.3,pp.273-284,March 1991.)を参考にしている。
【００４７】
複数のアンテナ（ここではｎ素子）１ａ、…、１ｎで受信した受信ベクトルは以下の式で与えられる。
【００４８】
Ｘ（ｔ）＝ＡＳ（ｔ）＋Ｎ（ｔ）
ここでＸ（ｔ）は受信ベクトル、Ａは到来方向におけるアレイマニフォルド、Ｓ（ｔ）は信号ベクトル、Ｎ（ｔ）はノイズベクトルである。
【００４９】
そして、偏差のある受信アナログ部２を通過した後の受信ベクトルＸ“（ｔ）は以下のようになる。
【００５０】
Ｘ“（ｔ）＝ＣΓＡＳ（ｔ）＋Ｎ（ｔ）
Γ＝ｄｉａｇ｛Γ１１、Γ２２、…、Γｎｎ｝
Ｃは、要素Ｃｉｊがｉ番目とｊ番目の受信アナログ部２の相互結合係数を意味する行列である。
【００５１】
この例では、振幅位相Γｉｉを算出する方法を記述する。本計算では、以下の評価関数Ｊｃを最小化する。この評価関数は、到来方向におけるＭＵＳ１Ｃスペクトラムの合計値の逆数を表している。つまり、ＭＵＳＩＣスペクトラムの値の最大化を目的とするが、計算の都合上、逆数を用いて最小化を目指す。
【００５２】
【数１】
【外１】
【００５３】
また、δ＝［Γ１１，Γ２２，…，Γｍｍ］Ｔ，Ｑ１（ｍ）＝ｄｉａｇ｛ａ （θｍ）｝である。但し、Ｔは転置。
【００５４】
このとき、Ｊｃを最小化するΓｉｉは以下の式で与えられる。
【００５５】
【数２】
【００５６】
以上説明したように、本実施形態の到来波方向推定装置では、複数の電波の受信信号のみを用いて、ＭＵＳＩＣアルゴリズムで得られる到来方向情報を利用し、補正量を算出（つまりキャリブレーション）できる。ここで、複数の電波の受信信号は、キャリブレーション用として特別に用意する必要はなく、実使用中の受信信号を用いるだけでよい。
【００５７】
また、振幅位相、そして相互結合までもが、リアルタイムにキャリブレーションできるようになる。更に、新たなアナログ部品を必要とすることもなく、装置を安価に構成することができる。
【００５８】
次に補正量算出手段４における第２の実施例の算出方法について図４のフローチャートを参照して説明する。始めに、振幅位相Γｉｉと相互結合係数Ｃｉｊに初期値を与える（ステップ４０１）。初期値としては、前回のキャリブレーションデータが存在するならばその値を、まったく未知な場合には、偏差や相互結合がない条件で与えればよい。次に、振幅位相偏差と相互結合係数を固定して到来方向をＭＵＳＩＣアルゴリズムで推定する（ステップ４０２）。ここでは、振幅位相Γｉｉと相互結合係数Ｃｉｊは固定である。
【００５９】
次に到来方向と相互結合係数を固定して振幅位相偏差Γｉｉを推定する（ステップ４０３）。ここでは、到来方向と相互結合係数は固定である。Γｉｉは、ＭＵＳＩＣスペクトラムのピークが到来方向においてより鋭くなる（大きな値になる）ように設定される。
【００６０】
更に、到来方向と振幅位相偏差を固定して相互結合係数Ｃｉｊを推定する（ステップ４０４）。ここでは、到来方向と振幅位相は固定である。Ｃｉｊは、ＭＵＳＩＣスペクトラムのピークが到来方向において鋭くなる（大きな値になる）ように設定される。このステップ４０２からステップ４０４の処理を繰り返し計算して、推定結果が収束して未知変数の値が変化しなくなったら計算を終了し（ステップ４０５）、補正量を到来方向推定手段５に出力する（ステップ４０６）。尚、ステップ４０３と４０４の処理における計算方法は、いかなる手法を用いても良い。
【００６１】
この実施例は、効率的に補正量を算出する方法について記述している。未知パラメータとしては、電波の到来方向、振幅位相偏差、相互結合係数の３種類である。これらを同時に推定することは未知パラメータの数が膨大になり、リアルタイムにキャリブレーションを行う場合に問題になる。そこで、これら３種類の未知変数を同時に算出するのではなく、順番に算出し、結果が収束するまで繰り返し求めていくものである。
【００６２】
このアルゴリズムを用いた場合の補正量算出結果の一例を示す。アンテナは６素子の円形アレイで、誤差として図５に示す振幅位相を設定した。このときの上記アルゴリズムを用い、ステップ４０２からステップ４０４の処理を繰り返し行った結果を図５に示す。尚、振幅位相の初期値としては、誤差のない値を用いている。また、振幅位相を求める計算は、第１の実施例で記述した方法で行った。図５より、初期値として誤差のない条件を用いているにも拘らず、補正量が計算できていることがわかる。値としては、１００％真の値にはなっていないが、十分キャリブレーションデータとして使用できる値になっている。
【００６３】
以上のように、本実施例によれば、未知変数を繰り返し計算で求めるため、効率的に計算することができるようになり、リアルタイムでのキャリブレーションを行うに際して高価な高速の計算機を用いなくともよく、装置を安価に構成することに大いに貢献することとなる。
【００６４】
次に補正量算出手段の第３の実施例の算出方法について図６のフローチャートを参照して説明する。始めに、振幅位相Γｉｉと相互結合係数Ｃｉｊに初期値を与える（ステップ６０１）。初期値としては、前回のキャリブレーションデータが存在するならばその値を、まったく未知な場合には、偏差や相互結合がない条件で与えればよい。次に、振幅位相偏差と相互結合係数を固定して到来方向をＭＵＳＩＣアルゴリズムで推定する（ステップ６０２）。ここでは、振幅位相Γｉｉと相互結合係数Ｃｉｊは固定である。次に到来方向と相互結合係数を固定して振幅位相偏差Γｉｉを推定する（ステップ６０３）。ここでは、到来方向と相互結合係数は固定である。Γｉｉは、ＭＵＳＩＣスペクトラムのピークが到来方向においてより鋭くなる（大きな値になる）ように設定される。次に、振幅位相偏差と相互結合係数を固定して到来方向をＭＵＳＩＣアルゴリズムで推定する（ステップ６０４）。ここでは、振幅位相Γｉｉと相互結合係数Ｃｉｊは固定である。
【００６５】
更に、到来方向と振幅位相偏差を固定して相互結合係数Ｃｉｊを推定する（ステップ６０５）。ここでは、到来方向と振幅位相は固定である。Ｃｉｊは、ＭＵＳＩＣスペクトラムのピークが到来方向において鋭くなる（大きな値になる）ように設定される。このステップ６０２からステップ６０５の処理を繰り返し計算して、推定結果が収束して未知変数の値が変化しなくなったら計算を終了し（ステップ６０６）、補正量を到来方向推定手段５に出力する（ステップ６０７）。
【００６６】
本実施例では、図４で示した第２の実施例に対して、計算順序を変え、特に到来波推定方向の算出をステップ６０２とステップ６０４の処理で２回行っている。
【００６７】
本実施例においては、図６に示したように、到来方向の推定→振幅位相算出→到来方向の推定→相互結合係数算出→…としている。このようにすると、到来方向推定のステップがひとつ増えて計算量が増加するように見えるが、繰り返し回数が減少する効果が期待できる。以下にその理由を説明する。
【００６８】
図４では、ステップ４０３で補正量が算出されるが、補正後のデータを用いて到来方向を推定し直すと、方向が変化している場合がある。つまり、図４に示したようなステップを踏むと、ステップ４０２で推定した到来方向と、ステップ４０３で補正したデータを用いて改めて推定した到来方向は異なっている場合がある。ここで、このままステップ４０４で補正を行うが、ステップ４０２の到来方向のデータと、ステップ４０３の補正データは整合が取れていないために、ステップ４０４で求まる補正量には誤差が含まれている場合がある。
【００６９】
そこで、本実施例では、振幅位相の補正量を算出した後に、ステップ６０４で到来方向を推定する。そして、相互結合係数の算出を行う。その結果、相互結合係数の算出においては、誤差を含まない形で求めることができるので、繰り返し回数を少なくして補正量を求めることができるようになる。
【００７０】
以上説明したように、本実施例によれば、補正量を算出した後に毎回到来方向を推定し直すので、誤差のない補正量算出が実現でき、その結果、少ない繰り返し回数で算出結果が収束し、リアルタイムでのキャリブレーションをより早く行うことができる。
【００７１】
また、振幅位相補正量が相互結合補正量に比ベて明らかに大きい場合には、ステップ６０２、６０３、６０２、６０３、６０２、６０４、６０２、６０３、６０２、６０３、６０２、６０４、…、と、振幅位相の補正量を求める回数を増やすことで、更に繰り返し回数を減少させる効果も期待できる。
【００７２】
図７は本発明の到来波方向推定装置の第２の実施形態の構成を示した図である。本例の到来波方向推定装置１００は、図７に示すように、補正量算出のためのひとつ或いは複数の電波の送信装置９ａ、９ｂを有している。また、前記電波の送信装置９ａ、９ｂは、前記到来波推定装置１００に対しての相対的な方向が既知となっている。尚、本例の到来波方向推定装置１００は図１に示した第１の実施形態と同様であるため、以降、図１を借用して説明する。
【００７３】
図８は図７に示した到来波方向推定装置の補正量算出手段４の補正量算出方法の実施例を示したフローチャートである。本例の補正量算出手段４は、図７に示すように既知の方向から到来する送信装置９ａ、９ｂから送信された電波の受信信号を用いて、振幅位相Γｉｉと相互結合係数Ｃｉｊを繰り返し計算して補正量を算出する。
【００７４】
始めに、振幅位相Γｉｉと相互結合係数Ｃｉｊに初期値を与える（ステップ８０１）。初期値としては、前回のキャリブレーションデータが存在するならばその値を、まったく未知な場合には、偏差や相互結合がない条件で与えればよい。
【００７５】
次に到来方向（これは既知である）と相互結合係数を固定して振幅位相偏差Γｉｉを推定する（ステップ８０２）。ここでは、到来方向と相互結合係数は固定である。ΓｉｉはＭＵＳＩＣスペクトラムのピークが到来方向においてより鋭くなる（大きな値になる）ように設定される。
【００７６】
更に、到来方向と振幅位相偏差を固定して相互結合係数Ｃｉｊを推定する（ステップ８０３）。ここでは、到来方向と振幅位相は固定である。Ｃｉｊは、ＭＵＳＩＣスペクトラムのピークが到来方向において鋭くなる（大きな値になる）ように設定される。このステップ８０２からステップ８０３の処理を繰り返し計算して、推定結果が収束して未知変数の値が変化しなくなったら計算を終了し（ステップ８０４）、補正量を到来方向推定手段５に出力する（ステップ８０５）。
【００７７】
本実施例では、効率的に補正量を算出する方法について記述している。本例では、既知の方向の波源からの電波の受信信号を用いているため、未知パラメータとしては、振幅位相偏差、相互結合係数の２種類だけである。つまり、未知パラメータであった到来方向が既知の条件となり、この条件で補正量を推定する。従って、方向推定の計算がまつたく不要になるので処理量も減少する。また、方向が既知であるので、最終的に得られる補正量の精度の向上が期待できる。
【００７８】
図９に方向が既知とした場合の補正量の算出結果を示す。シミュレーション条件は図５と同じである。図５に示した未知の方向からの受信信号を用いた場合に比べて、精度良く補正量を算出できていることが分かる。
【００７９】
尚、本例は、ひとつのみの送信装置を用いても、補正量は算出できる。しかし、複数の送信装置を用いることで、送信装置の位置ずれ、送信装置と到来波推定装置の間の伝搬環境の変化、送信装置の故障、などに対する対応が可能となり、補正量を安定且つ更に精度良く算出することができる。
【００８０】
以上説明したように、本実施形態によれば、補正量を算出するための電波の送信装置を用いて補正量を算出することで、処理量を大幅に少なくでき、また、精度の高い補正量が得られるキャリブレーションをリアルタイムに行うことができる。
【００８１】
図１０は本発明の到来波方向推定装置の第３の実施形態の構成を示したブロック図である。但し、図１に示した第１の実施形態と同様の部分には同一符号を用い、且つ適宜その説明を省略する。到来波方向推定装置は、複数のアンテナ素子１ａ、…、１ｎと、複数のアンテナ素子１ａ、…、１ｎに対応して設けられ、対応するアンテナ素子より伝達される受信信号を処理する複数の受信アナログ部２ａ、…、２ｎと、受信アナログ部２ａ、…、２ｎの出力信号をデジタル変換するＡ／Ｄ変換器３ａ、…、３ｎと、複数の受信アナログ部２ａ、…、２ｎの相対的な偏差を推定する補正量算出手段４と、補正量算出手段４の出力する補正量を用いて複数のアンテナ素子１ａ、…、１ｎに到来するひとつ或いは複数の電波の到来方向を推定する到来方向推定手段５及び補正量が算出される毎に、新たな補正量を記憶して記憶内容を更新する補正量記憶装置６を有している。
【００８２】
次に本実施形態の動作について説明する。本例の到来波方向推定動作、特に補正量算出手段４による補正量の算出動作は第１の実施形態のそれと同様である。異なる動作は、補正量算出手段４により算出された補正量を補正量記憶装置６によって記憶しているところにある。
【００８３】
本実施形態でも、補正量はリアルタイムに補正することが可能であるが、短時間の変動が小さい場合には、ある一定の時間間隔で補正すれば十分である。この場合には、補正量記憶装置６に記憶された補正量を用いて到来波方向の推定を間欠的に行えばよく、装置及び本装置を含むシステムの負荷を減らすことができる。また、補正量を計算する時の初期値として、補正量記憶装置６内の補正量を用いれば、補正量算出アルゴリズムの収束時間を短くすることができる。
【００８４】
以下、上記した補正量記憶装置６内の補正量を用いて、到来波の推定を行う動作の実施例について説明する。
【００８５】
本例の補正量算出手段４は、補正量記億装置６に記憶された過去の補正量を利用して、次の補正量を予測する。
【００８６】
ここで、補正量をリアルタイムに算出することは、精度を向上させる上では重要である。しかし、補正量の算出回数を減らし、その空いた時間には、補正量を算出するための計算リソースを到来方向推定アルゴリズムに利用し、方向推定に係る計算時間を短くすることはリアルタイム処理にとっては非常に重要である。
【００８７】
そこで、本例では、補正量記憶装置６に記憶された過去の補正量から現時点での補正量を予測し、その補正量を用いて電波の到来方向を推定する。
【００８８】
補正量を予測する方法の具体例を以下に示す。第１の例として補正量を線形補間で予測する。図１１は振幅位相の振幅を線形補間で予測する方法を示した例である。横軸は補正回数、縦軸は補正振幅量を表している。過去の２回の振幅補正量から現在の補正量を一次の線形補間で予測する。ここで、過去の３回およびそれ以上のデータを用いて２次の近似から予測することもできる。
【００８９】
第２の例として、補正量算出手段４によって算出された補正量と、そのときの外気気温の両方を補正量記憶装置６に記憶することで予測することもできる。受信アナログ部２の偏差が外気気温に相関がある場合がある。つまり、外気気温が分かっていれば、受信アナログ部２の偏差が予測できる場合がある。このような場合には、外気気温をモニターすることで、補正量を予測することもできる。または、受信アナログ部２が外気気温とは関わりがなく、受信アナログ部２を含んでいる装置内の気温と関わりがある場合がある。つまり、装置内の機器の発生する熱によって外気と装置内の温度差が異なる場合がある。この場合には、装置内の気温と補正量を補正量記憶装置６に記憶することが有効である。
【００９０】
第３の例として、補正量算出手段４によって算出された補正量と、そのときの時刻の両方を補正量記憶装置６に記憶する。例えば、構内無線や、移動体通信の基地局では、使用時間と使用頻度の間には相関がある。つまり、工場内の構内無線では、昼間に使用される頻度が高く、屋外の基地局では、社会人の仕事が終わった時間帯に使用される頻度が高い。使用する頻度が高いと、機器内の発生する熱量も多くなり、装置内の温度が上昇する。つまり、時刻と気温の間にも相関がある。そして、先に説明したように、気温と補正量の間にも相関がある。つまり、時刻と気温の間に相関がある。つまり、時刻と補正量を補正量記憶装置６に記憶しておくことは有効なことである。また、時刻によって外気の気温に変化がある。また、季節によって外気の差がある。そこで、時刻のほかにも、日時も記憶することは有効なことである。
【００９１】
また、記憶された離散的な時刻や日時と補正量の値から、線形補間によって記憶されていない時刻や日時における補正量を予測することも有効である。
【００９２】
同様に、記憶された離散的な気温と補正量の値から、線形補間によって記憶されていない気温における補正量を予測することも有効である。
【００９３】
更に、上記した気温、装置内部の温度、時刻や日時により、補正量算出時の初期値の設定をより適切なものにすることもできる。
【００９４】
尚、上記の各実施形態では、到来方向推定装置での説明であったが、同様に、アダプティブアンテナ装置の受信系の補正としてもまったく同様に用いることができる。
【００９５】
また、上記の各実施形態で得られる補正量と受信で用いたアダプティブアンテナ装置の重みから、送信に用いるべき重みを逆算することも勿論可能である。このときには、別に用意した送信用のキャリブレーションデータを必要とする。特にＴＤＭＡ／ＴＤＤシステムでは、送受で同一のパターンが要求されるので、有効である。以上の説明では、相互結合係数の補正量を算出していたが、受信アナログ部のシールドが精度良く構成されている場合には、相互結合係数を無視することができる。このような場合には、振幅位相の補正量のみを算出すればよい。また、以上の説明に用いた複数の実施形態を組み合わせて到来方向推定装置を構成することも勿論、可能である。
【００９６】
更に、本発明は上記実施形態に限定されることなく、その要旨を逸脱しない範囲において、具体的な構成、機能、作用、効果において、他の種々の形態によっても実施することができる。
【００９７】
【発明の効果】
以上詳細に説明したように、本発明によれば、運用中の受信データを用いてリアルタイムに高精度なキャリブレーションを行うことができ、到来波方向推定精度を向上させることができる。また、キャリプレーション用の送信装置を用いることで、さらに、精度の良い補正を行うことができるようになる。また、過去の補正量を記憶しておくことで、現時点での補正量を線形補間などで簡易に予測することができる。
【図面の簡単な説明】
【図１】本発明の到来波方向推定装置の第１の実施形態に係る構成を示したブロック図である。
【図２】図１に示した受信アナログ部の偏差を回路モデルで表した等価回路図である。
【図３】図１に示した補正量算出手段で用いるＭＵＳＩＣスペクトラムの特性を示した特性図である。
【図４】図１に示した補正量算出手段における補正量算出方法の第２の実施例を示したフローチャートである。
【図５】図１に示した補正量算出手段により算出される振幅位相Γｉｉの設定値と計算値を示した表図である。
【図６】図１に示した補正量算出手段における補正量算出方法の第３の実施例を示したフローチャートである。
【図７】本発明の到来波方向推定装置の第２の実施形態に係る構成を示したブロック図である。
【図８】図７に示した補正量算出手段における補正量算出方法の実施例を示したフローチャートである。
【図９】図７に示した補正量算出手段により算出される振幅位相Γｉｉの設定値と計算値を示した表図である。
【図１０】本発明の到来波方向推定装置の第３の実施形態に係る構成を示したブロック図である。
【図１１】図１０の補正量算出手段により求められる振幅量の線形補間方法を説明する図である。
【図１２】従来のアダプティブアンテナ装置の構成例を示したブロック図である。
【図１３】図１２の装置における送信パターンと受信パターン例を示した図である。
【図１４】従来の到来波方向推定装置の構成例を示したブロック図である。
【符号の説明】
１ａ、１ｎ アンテナ素子
２ａ、２ｎ 受信アナログ部
３ａ、３ｎ Ａ／Ｄ変換器
４ 補正量算出手段
５ 到来方向推定手段
６ 補正量記憶装置
９ａ、９ｎ 送信装置
２１ 低雑音増幅器
２２ フィルタ
２３ 周波数変換器
１００ 到来波方向推定装置
Claims (6)
- 複数のアンテナ素子と、
前記複数のアンテナ素子に対応して設けられ、対応するアンテナ素子より伝達される受信信号を前記アンテナ素子対応で処理する複数の受信部と、
前記複数の受信部から出力される受信信号から前記複数の受信部から出力される信号間の相対的な偏差を推定する補正量算出手段と、
前記複数の受信部から出力される受信信号を前記補正量算出手段により得られる偏差量を用いて補正し、これら補正された受信信号を用いて前記複数のアンテナ素子に到来する少なくともひとつの電波の到来方向を推定する到来方向推定手段と、を具備し、
前記補正量算出手段は、ひとつ或いは複数の電波の到来方向と、複数の受信部から出力される信号間の振幅位相偏差と、複数の受信部から出力される信号間の相互結合係数の３種類の未知変数を求める処理を行い、ステップ１として振幅位相偏差と相互結合係数を固定して到来方向を推定し、ステップ２として到来方向と相互結合係数を固定して振幅位相偏差を推定し、ステップ３として到来方向と振幅位相偏差を固定して相互結合係数を推定し、前記ステップ１からステップ３の処理を推定結果が収束するまで繰り返し実行することにより、前記偏差量を求めることを特徴とする到来波方向推定装置。 - 複数のアンテナ素子と、
前記複数のアンテナ素子に対応して設けられ、対応するアンテナ素子より伝達される受信信号を前記アンテナ素子対応で処理する複数の受信部と、
前記複数の受信部から出力される受信信号から前記複数の受信部から出力される信号間の相対的な偏差を推定する補正量算出手段と、
前記複数の受信部から出力される受信信号を前記補正量算出手段により得られる偏差量を用いて補正し、これら補正された受信信号を用いて前記複数のアンテナ素子に到来する少なくともひとつの電波の到来方向を推定する到来方向推定手段と、を具備し、
前記補正量算出手段は、ステップ１として振幅位相偏差と相互結合係数を固定して到来方向を推定し、ステップ２として到来方向と相互結合係数を固定して振幅位相偏差を推定し、ステップ３として振幅位相偏差と相互結合係数を固定して到来方向を推定し、ステップ４として到来方向と振幅位相偏差を固定して相互結合係数を推定し、前記ステップ１からステップ４の処理を推定結果が収束するまで繰り返し実行することにより、前記偏差量を求めることを特徴とする到来波方向推定装置。 - 複数のアンテナ素子と、
前記複数のアンテナ素子に対応して設けられ、対応するアンテナ素子より伝達される受信信号を前記アンテナ素子対応で処理する複数の受信部と、
前記複数の受信部から出力される受信信号から前記複数の受信部から出力される信号間の相対的な偏差を推定する補正量算出手段と、
前記複数の受信部から出力される受信信号を前記補正量算出手段により得られる偏差量を用いて補正し、これら補正された受信信号を用いて前記複数のアンテナ素子に到来する少なくともひとつの電波の到来方向を推定する到来方向推定手段と、を具備し、
前記到来波方向推定装置は、偏差量算出のためのひとつ或いは複数の電波の送信装置を有しており、前記電波の送信装置の電波発信位置は、前記到来波推定装置に対して相対的な方向が既知であり、
前記補正量算出手段は、前記送信装置から発せられた既知の方向から到来する電波の受信信号を用いて、ステップ１として到来方向と相互結合係数を固定して振幅位相偏差を推定し、ステップ２として到来方向と振幅位相偏差を固定して相互結合係数を推定し、前記ステップ１とステップ２の処理を推定結果が収束するまで繰り返し実行することにより、前記偏差量を求めることを特徴とする到来波方向推定装置。 - 複数のアンテナ素子と、
前記複数のアンテナ素子に対応して設けられ、対応するアンテナ素子より伝達される受信信号を前記アンテナ素子対応で処理する複数の受信部と、
前記複数の受信部から出力される受信信号から前記複数の受信部から出力される信号 間の相対的な偏差を推定する補正量算出手段と、
前記複数の受信部から出力される受信信号を前記補正量算出手段により得られる偏差量を用いて補正し、これら補正された受信信号を用いて前記複数のアンテナ素子に到来する少なくともひとつの電波の到来方向を推定する到来方向推定手段と、を具備し、
前記補正量算出手段は補正量記憶装置を有しており、偏差量を算出する毎に得られた偏差量を補正量記憶装置に記憶して記憶内容を更新し、
前記補正量記億装置は、前記補正量算出装置によって算出された偏差量と、そのときの外気気温の両方を記憶することを特徴とする到来波方向推定装置。 - 複数のアンテナ素子と、
前記複数のアンテナ素子に対応して設けられ、対応するアンテナ素子より伝達される受信信号を前記アンテナ素子対応で処理する複数の受信部と、
前記複数の受信部から出力される受信信号から前記複数の受信部から出力される信号間の相対的な偏差を推定する補正量算出手段と、
前記複数の受信部から出力される受信信号を前記補正量算出手段により得られる偏差量を用いて補正し、これら補正された受信信号を用いて前記複数のアンテナ素子に到来する少なくともひとつの電波の到来方向を推定する到来方向推定手段と、を具備し、
前記補正量算出手段は補正量記憶装置を有しており、偏差量を算出する毎に得られた偏差量を補正量記憶装置に記憶して記憶内容を更新し、
前記補正量記憶装置は、前記補正量算出装置によって算出された偏差量と、そのときの受信アナログ部の周りの気温の両方を記憶することを特徴とする到来波方向推定装置。 - 複数のアンテナ素子と、
前記複数のアンテナ素子に対応して設けられ、対応するアンテナ素子より伝達される受信信号を前記アンテナ素子対応で処理する複数の受信部と、
前記複数の受信部から出力される受信信号から前記複数の受信部から出力される信号間の相対的な偏差を推定する補正量算出手段と、
前記複数の受信部から出力される受信信号を前記補正量算出手段により得られる偏差量を用いて補正し、これら補正された受信信号を用いて前記複数のアンテナ素子に到来する少なくともひとつの電波の到来方向を推定する到来方向推定手段と、を具備し、
前記補正量算出手段は補正量記憶装置を有しており、偏差量を算出する毎に得られた偏差量を補正量記憶装置に記憶して記憶内容を更新し、
前記補正量記憶装置は、前記補正量算出手段によって算出された偏差量と、そのときの時刻、日時の両方を記憶することを特徴とする到来波方向推定装置。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2001019010A JP3640891B2 (ja) | 2001-01-26 | 2001-01-26 | 到来波方向推定装置及び到来波方向推定方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2001019010A JP3640891B2 (ja) | 2001-01-26 | 2001-01-26 | 到来波方向推定装置及び到来波方向推定方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2002221561A JP2002221561A (ja) | 2002-08-09 |
| JP3640891B2 true JP3640891B2 (ja) | 2005-04-20 |
Family
ID=18884957
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2001019010A Expired - Fee Related JP3640891B2 (ja) | 2001-01-26 | 2001-01-26 | 到来波方向推定装置及び到来波方向推定方法 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP3640891B2 (ja) |
Families Citing this family (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP4954112B2 (ja) * | 2008-02-18 | 2012-06-13 | 三菱電機株式会社 | アレー測角装置 |
| JP4794613B2 (ja) | 2008-10-15 | 2011-10-19 | 三菱電機株式会社 | 信号波到来角度測定装置 |
| RU2681942C1 (ru) * | 2018-03-28 | 2019-03-14 | Акционерное общество "Особое конструкторское бюро Московского энергетического института" | Способ фазовой пеленгации и фазовый пеленгатор |
| JP7248454B2 (ja) * | 2019-02-28 | 2023-03-29 | 株式会社デンソーテン | レーダ装置および補正値算出方法 |
| KR102243445B1 (ko) * | 2020-09-24 | 2021-04-22 | 한화시스템 주식회사 | 안테나 장치, 이를 구비하는 능동 위상배열 레이더, 및 이의 오차 교정방법 |
-
2001
- 2001-01-26 JP JP2001019010A patent/JP3640891B2/ja not_active Expired - Fee Related
Also Published As
| Publication number | Publication date |
|---|---|
| JP2002221561A (ja) | 2002-08-09 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| KR101019521B1 (ko) | 어레이 안테나 전송링크의 조정장치 및 방법 | |
| EP2549812B1 (en) | RF Fingerprinting for Location Estimate | |
| US6496140B1 (en) | Method for calibrating a smart-antenna array radio transceiver unit and calibrating system | |
| JP5367843B2 (ja) | 無線信号処理装置及び無線装置 | |
| JP2005197772A (ja) | アダプティブアレイアンテナ装置 | |
| JPH0653727A (ja) | アレイアンテナを有する無線周波送信機をビーム形成する方法 | |
| US8818291B2 (en) | Method and device for antenna calibration | |
| US20130044024A2 (en) | Method for iterative estimation of global parameters | |
| KR102001394B1 (ko) | 로그-영역 안테나 어레이 보간에 기반한 수신신호의 도래각 추정 방법과 이를 위한 장치 | |
| US20190173592A1 (en) | Beamforming calibration system and method | |
| KR20130049803A (ko) | 어레이 안테나의 교정 장치 및 교정 방법 | |
| JP2015012568A (ja) | 指向性制御装置、および指向性制御方法 | |
| JP4015180B2 (ja) | 到来方向推定装置、到来方向推定方法及び障害物推定装置 | |
| JP6292566B2 (ja) | 推定方法およびそれを利用した推定装置 | |
| JP3640891B2 (ja) | 到来波方向推定装置及び到来波方向推定方法 | |
| CN110429994A (zh) | 基于偶数根天线的均匀圆阵幅相误差自校正装置及方法 | |
| CN103248412B (zh) | 一种卫星导航多波束形成的方法 | |
| JP4014941B2 (ja) | 到来方向推定装置及び到来方向推定方法 | |
| JP6628596B2 (ja) | アレーアンテナ装置及び無線装置 | |
| JP2004147079A (ja) | 無線システム | |
| JP2002368664A (ja) | 通信装置およびその送信アレーアンテナ校正方法 | |
| JP2007504729A (ja) | 無線通信システムにおける信号の伝送方法、ならびに相応する送信局及び受信局 | |
| WO2005041443A1 (ja) | 無線通信装置 | |
| US10826596B2 (en) | Placement of calibration earth stations for ground based beam forming in a satellite system | |
| CN111181617A (zh) | 一种发射波束的形成方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20040618 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20040720 |
|
| A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20040917 |
|
| TRDD | Decision of grant or rejection written | ||
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20050104 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20050119 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20080128 Year of fee payment: 3 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090128 Year of fee payment: 4 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100128 Year of fee payment: 5 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110128 Year of fee payment: 6 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120128 Year of fee payment: 7 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130128 Year of fee payment: 8 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130128 Year of fee payment: 8 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20140128 Year of fee payment: 9 |
|
| LAPS | Cancellation because of no payment of annual fees |