JP2003101449A - アダプティブアレーアンテナのウエイト計算方法、ウエイト計算装置、ウエイト計算システム及びウエイト計算プログラムを記録した記録媒体 - Google Patents

Abstract

(57)【要約】 （修正有） 【課題】 ウエイト算出のための計算量を格段に削減し
たアダプティブアレーアンテナのウエイト計算方法・装
置・記録媒体を提供する。 【解決手段】 ウエイト計算器３１０は、Ｎ_ele個の参
照信号対雑音電力比ＤＮ_e及びＮ_ele個の複素インパルス
応答ｈ_e(k)を入力し、まずＮ_ele×Ｎ_ele次の行列ＮとＨ
を計算し、正方行列Ｎのｐ行ｑ列の要素をＮ（ｐ，ｑ）
と表すと、Ｎの対角成分は式Ｎ（ｅ，ｅ）＝１／ＤＮ_e
を用いて算出し、その他の成分は０とする。また、Ｈの
要素は、次式により算出する。但しＨ_p ^*はＨ_pの複素共
役を表す。 更に、Ｒ_xx＝Ｈ＋ＮなるＮ_ele×Ｎ_ele次の行列Ｒ_xxを計
算し、Ｎ_ele次の列ベクトルＲ_xdを計算し、ｅ行目の要
素をＲ_xd（e）と表すと、式Ｒ_xd＝Ｈ_e ^*（０）を用いて
算出する。そして、Ｒ_xxとＲ_xdを用いて目的のウエイト
Ｗ_e（ｅ＝１，２，…，Ｎ_ele）を計算する。

Description

【発明の詳細な説明】

【０００１】

【発明の属する技術分野】本発明は、アダプティブアレ
ーアンテナのウエイト計算方法、ウエイト計算装置、ウ
エイト計算システム及びウエイト計算プログラムを記録
した記録媒体に関する。

【０００２】

【従来の技術】デジタル無線通信システムにおいて、周
波数選択性フェージングを軽減する方式として、アレー
アンテナの指向性を制御し遅延波を抑圧する、アダプテ
ィブアレーアンテナ方式がある。アダプティブアレーア
ンテナ方式は、複数のアンテナに適切な重み付けを行う
ことによって非希望波を抑圧し、希望波の信号だけを取
り出すようにアレーの出力を最適に制御する技術であ
る。

【０００３】アダプティブアレーアンテナを制御する方
式として、“抑圧パイロット信号を用いたビームフォー
ミング法”電子情報通信学会論文誌Vol. J81-B-1 No.11
pp.661（1998-11発行）に示されるような抑圧パイロッ
ト信号を用いたビームフォーミング方式（以下、方式Ａ
と呼ぶ）がある。この方式では、伝送信号の中に含まれ
ているパイロット信号を用いて伝搬路のインパルス応答
を測定し、これを基にしてアレーアンテナのウエイトを
計算し、制御している。具体例として、上記方式Ａの受
信方式について詳細を述べる。

【０００４】図１４は、上記方式Ａの受信器のアダプテ
ィブアレーアンテナの制御に関する部分の構成を示す図
である。図１４において、１０１はアレーアンテナ、１
０２はＲＦ，ＩＦ，ダウンコンバータ，Ａ／Ｄ変換器，
ＬＰＦ等からなるダウンコンバート部、１０３はデジタ
ルベースバンド信号、１０４は複素インパルス応答測定
器、１０５は指向性制御器、１０６は遅延波の抑圧され
た信号、１１０はウエイト計算器、１１１は複素インパ
ルス応答、１１２はウエイトである。

【０００５】４本のアンテナ素子からなるアレーアンテ
ナ１０１で受信した信号は、ＲＦ，ＩＦ，ダウンコンバ
ータ，Ａ／Ｄ変換器，ＬＰＦ等からなるダウンコンバー
ト部１０２を経てデジタルベースバンド信号１０３に変
換される。受信信号には、パイロット信号が含まれてい
るため、複素インパルス応答測定器１０４にて逆拡散処
理及び同期加算処理を行うことによって、伝搬路の複素
インパルス応答１１１が作成される。これを用いてウエ
イト計算器１１０によりウエイト１１２が計算される。
ウエイト計算方法についての詳細は後述する。指向性制
御器１０５は、ウエイト１１２を用いてアンテナの指向
性を制御し、遅延波の抑圧された信号１０６を復調部へ
と出力する。

【０００６】図１５は、上記ウエイト計算器１１０の構
成を示す図である。図１５において、ウエイト計算器１
１０は、擬似送信信号発生器２０１、擬似受信信号発生
器２０３、及び最小二乗法計算器２０５から構成され
る。２０２は擬似送信信号、２０４は擬似受信信号発生
器２０３により発生した擬似受信信号である。

【０００７】図１６は、上記受信器の信号処理のタイム
チャートを示す図である。一般に無線伝搬路の状態は時
間的に変化するため、ウエイトもこれに応じて更新する
必要がある。上記方式Ａでは、伝送信号をフレーム化
し、１フレーム分のパイロット信号を用いてインパルス
応答ｈ_e(k)を測定し、これを用いてウエイトＷ_eを算出
している。このような作業を毎フレーム行いウエイトを
更新している。

【０００８】以上の構成において、図１５に示すよう
に、入力は複素インパルス応答ｈ_e(k)１１１であり、出
力はウエイトＷ_e１１２である。ここで添え字“e”はア
ンテナ素子を表す。方式Ａでは、アレーアンテナが４本
のアンテナ素子から構成されているため、１フレーム分
の信号処理を行うために、ｅ＝１，２，３，４ なる４
個のアンテナ素子に対応する４個の複素インパルス応答
ｈ_e(k)を入力し、４個のウエイトＷ_eを算出、出力す
る。また、“ｋ”は（遅延）時間を表す整数である。ｋ
＝１は１シンボル時間を表す。

【０００９】擬似送信信号発生器２０１では、ランダム
なＱＰＳＫ（Quadrature Phase Shift Keying）信号を
擬似送信信号ｄ(k)２０２として生成する。擬似受信信
号発生器２０３では、複素インパルス応答ｈ_e(k)をタッ
プ係数にもつフィルタに擬似送信信号ｄ(k)を通すこと
によって、アンテナ素子“ｅ”に対応する擬似受信信号
Ｘ_e(k)２０４を生成する。また、このとき雑音ｎ_e(k)も
加える。雑音の電力は、擬似受信信号Ｘ_e(k)２０４のＳ
Ｎ比と図１４のデジタルベースバンド信号１０３にて示
されるベースバンド信号のＳＮ比が同じになるようにあ
らかじめ調整しておく。これらの処理を全てのアンテナ
素子ｅ＝１，２，３，４ について行う。擬似受信信号
発生器における信号処理を数式を用いて表現すると、次
式（１）となる。 Ｘ_e(k)＝ｄ(k)※ｈ_e(k) ＋ｎ_e(k) …（１） 但し“※”はｋに関する畳み込み演算を表す。

【００１０】最小二乗法計算器２０５では、擬似送信信
号ｄ(k)と擬似受信信号Ｘ_e(k)（ｅ＝１，２，３，４）
を入力し、ＲＬＳ（Recursive Least Square：再帰最小
二乗）法を用いてウエイトＷ_e（ｅ＝１，２，３，４）
１１２を計算、出力する。また、ここでＲＬＳ法のルー
プ回数は５０回としている。よって最小二乗法計算器２
０５は、１フレーム当たりｋ＝１，２，…，５０ なる
５０組のｄ(k)とＸ_e(k)（ｅ＝１，２，３，４）とを入
力し、これらの情報を基に統計処理を行い、最適なＷ_e
を計算、出力していると言うこともできる。上記方式Ａ
ではＲＬＳ法が採用されているが、最小二乗法を実現す
る他のアルゴリズム、例えばＳＭＩ法（逆行列演算法）
等を採用することも可能である。

【００１１】

【発明が解決しようとする課題】しかしながら、このよ
うな従来のアダプティブアレーアンテナのウエイト計算
方法にあっては以下のような問題点があった。図１５に
示す方式Ａのウエイト計算器１１０では、統計処理を行
ってウエイトを算出している。よって、ウエイトを高精
度に算出するためには統計に用いるサンプル数を増大
（例えば、方式Ａでは５０サンプル）させる必要があ
り、計算量が大きくなるという課題がある。また、ＲＬ
Ｓ法よりも比較的計算量の少ないＳＭＩ法を用いても同
様である。

【００１２】本発明は、このような課題に鑑みてなされ
たものであって、ウエイト算出のための計算量を格段に
削減することのできるアダプティブアレーアンテナのウ
エイト計算方法、ウエイト計算装置、ウエイト計算シス
テム及びウエイト計算プログラムを記録した記録媒体を
提供することを目的としている。

【００１３】

【課題を解決するための手段】本発明は、アダプティブ
アレーアンテナを用いて受信した受信信号から伝搬路の
インパルス応答を測定し、該伝搬路のインパルス応答か
ら前記アンテナのウエイトを計算するウエイト計算装置
及び方法において、あらかじめ設定された数Ｎ_eleに対
して、Ｎ_ele個の参照信号対雑音電力比ＤＮ_e（ｅ＝１，
２，…，Ｎ_ele）を入力し、Ｎ_ele個の複素インパルス応
答ｈ_e(k)（ｅ＝１，２，…，Ｎ_ele）を入力し、ｐ行ｑ
列の要素Ｎ（ｐ，ｑ）が次式で表されるＮ_ele次の正方
行列Ｎを計算する手段と、

【００１４】

【数３４】 但しＤＮ_p ^-1はＤＮ_pの逆数を表す ｐ行ｑ列の要素Ｈ（ｐ，ｑ）が次式で表されるＮ_ele次
の正方行列Ｈを計算する手段と、

【００１５】

【数３５】 但しｈ_p ^*(k)はｈ_p(k)の複素共役を表す 式Ｒ_xx＝Ｈ＋Ｎで表されるＮ_ele次の正方行列Ｒ_xxを計
算する手段と、ｐ行の要素Ｒ_xd(p)が次式で表されるＮ
_ele次の列ベクトルＲ_xdを計算する手段と、 Ｒ_xd(p)＝ｈ_p ^*(0) 但しｈ_p ^*(0)はｈ_p(0)の複素共役を表す 連立一次方程式Ｒ_xxＷ＝Ｒ_xdを満足するＮ_ele次の列ベ
クトルＷを計算する手段と、前記列ベクトルＷのｐ行目
の要素Ｗ(p)を、ウエイトＷ_pとし、Ｎ_ele個のウエイト
Ｗ_p（ｐ＝１，２，…，Ｎ_ele）を出力する出力手段とを
備えることを特徴としている。

【００１６】本発明は、アダプティブアレーアンテナを
用いて受信した受信信号から伝搬路のインパルス応答を
測定し、該伝搬路のインパルス応答から前記アンテナの
ウエイトを計算するウエイト計算装置及び方法におい
て、あらかじめ設定された数Ｎ_eleに対して、Ｎ_ele個の
参照信号対雑音電力比の逆数ＤＮ_e ^-1（ｅ＝１，２，
…，Ｎ_ele）を入力し、Ｎ_ele個の複素インパルス応答ｈ
_e(k)（ｅ＝１，２，…，Ｎ_ele）を入力し、ｐ行ｑ列の
要素Ｎ（ｐ，ｑ）が次式で表されるＮ_ele次の正方行列
Ｎを計算する手段と、

【００１７】

【数３６】 但しＤＮ_p ^-1はＤＮ_pの逆数を表す ｐ行ｑ列の要素Ｈ（ｐ，ｑ）が次式で表されるＮ_ele次
の正方行列Ｈを計算する手段と、

【００１８】

【数３７】 但しｈ_p ^*(k)はｈ_p(k)の複素共役を表す 式Ｒ_xx＝Ｈ＋Ｎで表されるＮ_ele次の正方行列Ｒ_xxを計
算する手段と、ｐ行の要素Ｒ_xd(p)が次式で表されるＮ
_ele次の列ベクトルＲ_xdを計算する手段と、 Ｒ_xd(p)＝ｈ_p ^*(0) 但しｈ_p ^*(0)はｈ_p(0)の複素共役を表す 連立一次方程式Ｒ_xxＷ＝Ｒ_xdを満足するＮ_ele次の列ベ
クトルＷを計算する手段と、前記列ベクトルＷのｐ行目
の要素Ｗ(p)を、ウエイトＷ_pとし、Ｎ_ele個のウエイト
Ｗ_p（ｐ＝１，２，…，Ｎ_ele）を出力する出力手段とを
備えることを特徴としている。

【００１９】本発明は、アダプティブアレーアンテナを
用いて受信した受信信号から伝搬路のインパルス応答を
測定し、該伝搬路のインパルス応答から前記アンテナの
ウエイトを計算するウエイト計算装置及び方法におい
て、あらかじめ設定された数Ｎ_eleに対して、Ｎ_ele個の
参照信号対雑音電力比ＤＮ_e（ｅ＝１，２，…，Ｎ_ele）
を入力し、Ｎ_ele個の複素インパルス応答ｈ_e(k)（ｅ＝
１，２，…，Ｎ_ele）を入力し、ｐ行ｑ列の要素Ｎ
（ｐ，ｑ）が次式で表されるＮ_ele次の正方行列Ｎを計
算する手段と、

【００２０】

【数３８】 但しＤＮ_p ^-1はＤＮ_pの逆数を表す ｐ行ｑ列の要素Ｈ（ｐ，ｑ）が次式で表されるＮ_ele次
の正方行列Ｈを計算する手段と、

【００２１】

【数３９】 但しｈ_p ^*(k)はｈ_p(k)の複素共役を表す 式Ｒ_xx＝Ｈ＋Ｎで表されるＮ_ele次の正方行列Ｒ_xxを計
算する手段と、ｐ行の要素Ｒ_xd(p)が次式で表されるＮ
_ele次の列ベクトルＲ_xdを計算する手段と、 Ｒ_xd(p)＝ｈ_p ^*(0) 但しｈ_p ^*(0)はｈ_p(0)の複素共役を表す 式Ｗ＝Ｒ_xx ^-1・Ｒ_xd（但しＲ_xx ^-1はＲ_xxの逆行列を表
す）で表されるＮ_ele次の列ベクトルＷを計算する手段
と、前記列ベクトルＷのｐ行目の要素Ｗ(p)を、ウエイ
トＷ_pとし、Ｎ_ele個のウエイトＷ_p（ｐ＝１，２，…，
Ｎ_ele）を出力する出力手段とを備えることを特徴とし
ている。

【００２２】本発明は、アダプティブアレーアンテナを
用いて受信した受信信号から伝搬路のインパルス応答を
測定し、該伝搬路のインパルス応答から前記アンテナの
ウエイトを計算するウエイト計算装置及び方法におい
て、あらかじめ設定された数Ｎ_eleに対して、Ｎ_ele個の
参照信号対雑音電力比の逆数ＤＮ_e ^-1（ｅ＝１，２，
…，Ｎ_ele）を入力し、Ｎ_ele個の複素インパルス応答ｈ
_e(k)（ｅ＝１，２，…，Ｎ_ele）を入力し、ｐ行ｑ列の
要素Ｎ（ｐ，ｑ）が次式で表されるＮ_ele次の正方行列
Ｎを計算する手段と、

【００２３】

【数４０】 但しＤＮ_p ^-1はＤＮ_pの逆数を表す ｐ行ｑ列の要素Ｈ（ｐ，ｑ）が次式で表されるＮ_ele次
の正方行列Ｈを計算する手段と、

【００２４】

【数４１】 但しｈ_p ^*(k)はｈ_p(k)の複素共役を表す 式Ｒ_xx＝Ｈ＋Ｎで表されるＮ_ele次の正方行列Ｒ_xxを計
算する手段と、ｐ行の要素Ｒ_xd(p)が次式で表されるＮ
_ele次の列ベクトルＲ_xdを計算する手段と、 Ｒ_xd(p)＝ｈ_p ^*(0) 但しｈ_p ^*(0)はｈ_p(0)の複素共役を表す 式Ｗ＝Ｒ_xx ^-1・Ｒ_xd（但しＲ_xx ^-1はＲ_xxの逆行列を表
す）で表されるＮ_ele次の列ベクトルＷを計算する手段
と、前記列ベクトルＷのｐ行目の要素Ｗ(p)を、ウエイ
トＷ_pとし、Ｎ_ele個のウエイトＷ_p（ｐ＝１，２，…，
Ｎ_ele）を出力する出力手段とを備えることを特徴とし
ている。

【００２５】本発明は、アダプティブアレーアンテナを
用いて受信した受信信号から伝搬路のインパルス応答を
測定し、該伝搬路のインパルス応答から前記アンテナの
ウエイトを計算するウエイト計算装置及び方法におい
て、あらかじめ設定された数Ｎ_ele、閾値ＴＨに対し
て、Ｎ_ele個の参照信号対雑音電力比ＤＮ_e（ｅ＝１，
２，…，Ｎ_ele）を入力し、Ｎ_ele個の複素インパルス応
答ｈ_e(k)（ｅ＝１，２，…，Ｎ_ele）を入力し、ｐ行ｑ
列の要素Ｎ（ｐ，ｑ）が次式で表されるＮ_ele次の正方
行列Ｎを計算する手段と、

【００２６】

【数４２】 但しＤＮ_p ^-1はＤＮ_pの逆数を表す ｐ行ｑ列の要素Ｈ（ｐ，ｑ）が次式で表されるＮ_ele次
の正方行列Ｈを計算する手段と、

【００２７】

【数４３】 但しｈ_p ^*(k)はｈ_p(k)の複素共役を表す 式Ｒ_xx＝Ｈ＋Ｎで表されるＮ_ele次の正方行列Ｒ_xxを計
算する手段と、ｐ行の要素Ｒ_xd(p)が次式で表されるＮ
_ele次の列ベクトルＲ_xdを計算する手段と、 Ｒ_xd(p)＝ｈ_p ^*(0) 但しｈ_p ^*(0)はｈ_p(0)の複素共役を表す 式Ｗ１＝Ｒ_xx ^-1・Ｒ_xd（但しＲ_xx ^-1はＲ_xxの逆行列を表
す）で表されるＮ_ele次の列ベクトルＷ１を計算する手
段と、ｅ＝１，２，…，Ｎ_eleの中から最大の｜Ｈ_e(0)
｜（但し｜Ｈ_e(0)｜はＨ_e(0)の絶対値を表す）を持つｅ
＝ｅ_maxを探す手段と、ｐ行の要素Ｗ０(p)が次式で表さ
れるＮ_ele次の列ベクトルＷ０を計算する手段と、

【００２８】

【数４４】 Ｒ_xxの行列式Δを計算する手段と、行列式Δの絶対値｜
Δ｜を前記閾値ＴＨと比較する手段と、｜Δ｜≧ＴＨの
場合、Ｗ１のｐ行目の要素Ｗ１(p)を、ウエイトＷ_pと
し、｜Δ｜＜ＴＨの場合、Ｗ０のｐ行目の要素Ｗ０(p)
を、ウエイトＷ_pとし、Ｎ_ele個のウエイトＷ_p（ｐ＝
１，２，…，Ｎ_ele）を出力する出力手段とを備えるこ
とを特徴としている。

【００２９】本発明は、アダプティブアレーアンテナを
用いて受信した受信信号から伝搬路のインパルス応答を
測定し、該伝搬路のインパルス応答から前記アンテナの
ウエイトを計算するウエイト計算装置及び方法におい
て、あらかじめ設定された数Ｎ_ele、閾値ＴＨに対し
て、Ｎ_ele個の参照信号対雑音電力比の逆数ＤＮ_e ^-1（ｅ
＝１，２，…，Ｎ_ele）を入力し、Ｎ_ele個の複素インパ
ルス応答ｈ_e(k)（ｅ＝１，２，…，Ｎ_ele）を入力し、
ｐ行ｑ列の要素Ｎ（ｐ，ｑ）が次式で表されるＮ_ele次
の正方行列Ｎを計算する手段と、

【００３０】

【数４５】 但しＤＮ_p ^-1はＤＮ_pの逆数を表す ｐ行ｑ列の要素Ｈ（ｐ，ｑ）が次式で表されるＮ_ele次
の正方行列Ｈを計算する手段と、

【００３１】

【数４６】 但しｈ_p ^*(k)はｈ_p(k)の複素共役を表す 式Ｒ_xx＝Ｈ＋Ｎで表されるＮ_ele次の正方行列Ｒ_xxを計
算する手段と、ｐ行の要素Ｒ_xd(p)が次式で表されるＮ
_ele次の列ベクトルＲ_xdを計算する手段と、 Ｒ_xd(p)＝ｈ_p ^*(0) 但しｈ_p ^*(0)はｈ_p(0)の複素共役を表す 式Ｗ１＝Ｒ_xx ^-1・Ｒ_xd（但しＲ_xx ^-1はＲ_xxの逆行列を表
す）で表されるＮ_ele次の列ベクトルＷ１を計算する手
段と、ｅ＝１，２，…，Ｎ_eleの中から最大の｜Ｈ_e(0)
｜（但し｜Ｈ_e(0)｜はＨ_e(0)の絶対値を表す）を持つｅ
＝ｅ_maxを探す手段と、ｐ行の要素Ｗ０(p)が次式で表さ
れるＮ_ele次の列ベクトルＷ０を計算する手段と、

【００３２】

【数４７】 Ｒ_xxの行列式Δを計算する手段と、行列式Δの絶対値｜
Δ｜を前記閾値ＴＨと比較する手段と、｜Δ｜≧ＴＨの
場合、Ｗ１のｐ行目の要素Ｗ１(p)を、ウエイトＷ_pと
し、｜Δ｜＜ＴＨの場合、Ｗ０のｐ行目の要素Ｗ０(p)
を、ウエイトＷ_pとし、Ｎ_ele個のウエイトＷ_p（ｐ＝
１，２，…，Ｎ_ele）を出力する出力手段とを備えるこ
とを特徴としている。

【００３３】他の観点において本発明は、アダプティブ
アレーアンテナを用いて受信した受信信号から伝搬路の
インパルス応答を測定するインパルス応答測定装置と、
伝搬路のインパルス応答から前記アンテナのウエイトを
計算するウエイト計算装置と、計算されたウエイトを用
いて前記アンテナの指向性を制御し、遅延波の抑圧され
た信号を復調部へと出力する指向性制御装置とを備える
ウエイト計算システムにおいて、前記ウエイト計算装置
は、請求項７乃至１２のいずれか一項に記載のウエイト
計算装置を用いることを特徴としている。

【００３４】さらに、本発明は、アダプティブアレーア
ンテナを用いて受信した受信信号から伝搬路のインパル
ス応答を測定し、該伝搬路のインパルス応答から前記ア
ンテナのウエイトを計算するウエイト計算方法におい
て、あらかじめ設定された数Ｎ_eleに対して、Ｎ_ele個の
参照信号対雑音電力比ＤＮ_e（ｅ＝１，２，…，Ｎ_ele）
を入力し、Ｎ_ele個の複素インパルス応答ｈ_e(k)（ｅ＝
１，２，…，Ｎ_ele）を入力し、ｐ行ｑ列の要素Ｎ
（ｐ，ｑ）が次式で表されるＮ_ele次の正方行列Ｎを計
算し、

【００３５】

【数４８】 但しＤＮ_p ^-1はＤＮ_pの逆数を表す ｐ行ｑ列の要素Ｈ（ｐ，ｑ）が次式で表されるＮ_ele次
の正方行列Ｈを計算し、

【００３６】

【数４９】 但しｈ_p ^*(k)はｈ_p(k)の複素共役を表す 式Ｒ_xx＝Ｈ＋Ｎで表されるＮ_ele次の正方行列Ｒ_xxを計
算し、ｐ行の要素Ｒ_xd(p)が次式で表されるＮ_ele次の列
ベクトルＲ_xdを計算し、 Ｒ_xd(p)＝ｈ_p ^*(0) 但しｈ_p ^*(0)はｈ_p(0)の複素共役を表す 連立一次方程式Ｒ_xx・Ｗ＝Ｒ_xdを満足するＮ_ele次の列
ベクトルＷを計算し、前記列ベクトルＷのｐ行目の要素
Ｗ(p)を、ウエイトＷ_pとし、Ｎ_ele個のウエイトＷ_p（ｐ
＝１，２，…，Ｎ_ele）を出力するウエイト計算方法を
コンピュータに実行させるためのプログラムを記録した
ことを特徴とするコンピュータ読み取り可能な記録媒体
である。

【００３７】また、本発明は、アダプティブアレーアン
テナを用いて受信した受信信号から伝搬路のインパルス
応答を測定し、該伝搬路のインパルス応答から前記アン
テナのウエイトを計算するウエイト計算方法において、
あらかじめ設定された数Ｎ_ele、閾値ＴＨに対して、Ｎ
_ele個の参照信号対雑音電力比ＤＮ_e（ｅ＝１，２，…，
Ｎ_ele）を入力し、Ｎ_ele個の複素インパルス応答ｈ_e(k)
（ｅ＝１，２，…，Ｎ_ele）を入力し、ｐ行ｑ列の要素
Ｎ（ｐ，ｑ）が次式で表されるＮ_ele次の正方行列Ｎを
計算し、

【００３８】

【数５０】 但しＤＮ_p ^-1はＤＮ_pの逆数を表す ｐ行ｑ列の要素Ｈ（ｐ，ｑ）が次式で表されるＮ_ele次
の正方行列Ｈを計算し、

【００３９】

【数５１】 但しｈ_p ^*(k)はｈ_p(k)の複素共役を表す 式Ｒ_xx＝Ｈ＋Ｎで表されるＮ_ele次の正方行列Ｒ_xxを計
算し、ｐ行の要素Ｒ_xd(p)が次式で表されるＮ_ele次の列
ベクトルＲ_xdを計算し、 Ｒ_xd(p)＝ｈ_p ^*(0) 但しｈ_p ^*(0)はｈ_p(0)の複素共役を表す 式Ｗ１＝Ｒ_xx ^-1・Ｒ_xd（但しＲ_xx ^-1はＲ_xxの逆行列を表
す）で表されるＮ_ele次の列ベクトルＷ１を計算し、ｅ
＝１，２，…，Ｎ_eleの中から最大の｜Ｈ_e(0)｜（但し
｜Ｈ_e(0)｜はＨ_e(0)の絶対値を表す）を持つｅ＝ｅ_max
を探し、ｐ行の要素Ｗ０(p)が次式で表されるＮ_ele次の
列ベクトルＷ０を計算し、

【００４０】

【数５２】 Ｒ_xxの行列式Δを計算し、行列式Δの絶対値｜Δ｜を前
記閾値ＴＨと比較し、｜Δ｜≧ＴＨの場合、Ｗ１のｐ行
目の要素Ｗ１(p)を、ウエイトＷ_pとし、｜Δ｜＜ＴＨの
場合、Ｗ０のｐ行目の要素Ｗ０(p)を、ウエイトＷ_pと
し、Ｎ_ele個のウエイトＷ_p（ｐ＝１，２，…，Ｎ_ele）
を出力するウエイト計算方法をコンピュータに実行させ
るためのプログラムを記録したことを特徴とするコンピ
ュータ読み取り可能な記録媒体である。

【００４１】

【発明の実施の形態】以下、添付図面を参照しながら本
発明の好適なアダプティブアレーアンテナのウエイト計
算方法、ウエイト計算装置及びウエイト計算システムの
実施の形態について詳細に説明する。

【００４２】第１の実施の形態 図１は、本発明の第１の実施の形態の受信器のアダプテ
ィブアレーアンテナの制御に関する部分の構成を示す図
である。

【００４３】図１において、３０１はアレーアンテナ、
３０２はＲＦ，ＩＦ，ダウンコンバータ，Ａ／Ｄ変換
器，ＬＰＦ等からなるダウンコンバート部、３０３はデ
ジタルベースバンド信号、３０４は複素インパルス応答
測定器（CIRE:Complex ImpulseResponse Estimator）
（インパルス応答測定装置）、３０５は指向性制御器
（指向性制御装置）、３０６は遅延波の抑圧された信
号、３１０はウエイト計算器（ウエイト計算装置）、３
１１は複素インパルス応答ｈ_e(k)、３１２はウエイトＷ
_e、３１３は参照信号対雑音電力比ＤＮ_eである。

【００４４】アレーアンテナ３０１は、Ｎ_ele＝４本の
アンテナ素子から構成される。これらのアンテナで受信
された信号は、ダウンコンバート部３０２においてシン
ボルレートの４倍のサンプリングレート４００ＭＳａｍ
ｐｌｅｓ／ｓｅｃのデジタル複素ベースバンド信号３０
３に変換される。ダウンコンバート部３０２からは参照
信号対雑音電力比ＤＮ_e（ｅ＝１，２，…，Ｎ_ele）３１
３が出力される。ダウンコンバート部３０２の構成につ
いては図３により後述する。受信信号には、パイロット
信号が含まれているため、複素インパルス応答測定器３
０４にて逆拡散処理及び同期加算処理を行うことによっ
て、伝搬路の複素インパルス応答３１１が作成される。
これを用いてウエイト計算器３１０によりウエイトＷ_e
３１２が計算される。本発明は、ウエイト計算方法に特
徴があり、ウエイト計算の詳細については図５により後
述する。指向性制御器３０５は、ウエイトＷ_e３１２を
用いてアンテナの指向性を制御し、遅延波の抑圧された
信号３０６を復調部へと出力する。

【００４５】図２は、本発明にかかるウエイト計算方法
をハードウェアで行う場合のウエイト計算器３１０の構
成を示す図である。図２において、ウエイト計算器３１
０は、複素インパルス応答ｈ_e(k)からＲ_xdベクトルを計
算するＲ_xdベクトル計算器４０１、複素インパルス応答
ｈ_e(k)からＨ行列を計算するＨ行列計算器４０２、参照
信号対雑音電力比ＤＮ_eからＮ行列を計算するＮ行列計
算器４０３、Ｈ行列とＮ行列からＲ_xx行列を計算するＲ
_xx行列計算器４０４、及びＲ_xx行列とＲ_xdベクトルから
ウエイトＷ_eを計算する連立一次方程式計算器４０５か
ら構成される。

【００４６】図３は、上記ダウンコンバート部３０２の
構成を示す図である。図３において、ダウンコンバート
部３０２は、ＲＦ，ＩＦ，ダウンコンバータ，Ａ／Ｄ変
換器，ＬＰＦ等を備え、アンテナからの信号５０２を入
力し、増幅、周波数変換、Ａ／Ｄ変換、フィルタリング
等を行い、デジタル複素ベースバンド信号５０３を出力
する。ダウンコンバート部３０２の中には、ＡＧＣ（Au
tomatic Gain Control）部５０１がＮ_ele個あり、アン
テナ素子ｅに対応する部分から参照信号対雑音電力比Ｄ
Ｎ_eを出力する。

【００４７】図４は、本実施の形態で用いる伝送信号を
送信する送信器の構成を示す図である。本実施の形態で
用いる伝送信号は、データチャネルと、それに比べて電
力が十分に抑圧されたパイロットチャネルから構成され
る。データチャネル及びパイロットチャネルは、同じキ
ャリア上にあり、同じ帯域幅を持つ。パイロットチャネ
ルは、受信器において複素インパルス応答を測定し、ア
レーアンテナの指向性を制御するために用いる。

【００４８】図４において、送信器６００は、同相チャ
ンネル（I-ch）用ＰＮ（Pseudo Noise）系列発生器６０
１、直交チャンネル（Q-ch）用ＰＮ系列発生器６０２、
送信データ（１００Ｍ ｂｐｓ）６０８を元にベースバ
ンド信号であるデータチャネルを生成するＱＰＳＫ（Qu
adrature Phase Shift Keying）変調器６０３、ＰＮ系
列発生器６０１，６０２により生成したＭ（Maximum Le
ngth Shift Resister）系列からパイロットチャネルを
生成するＱＰＳＫ変調器６０４、ＱＰＳＫ変調器６０４
で生成した信号を減衰する減衰器６０５、データチャネ
ル及び抑圧パイロットチャネルを加算する加算器６０
６、及びアップコンバート部６０７を備えて構成され
る。ＰＮ系列として８段のＭ系列を用いているため、Ｐ
Ｎ系列の周期は２５５シンボルとなる。また、データチ
ャネル、抑圧パイロットチャネル共に信号レートは１０
０Ｍ Ｓｙｍｂｏｌｓ／ｓｅｃである。アップコンバー
ト部６０７は、帯域制限、周波数変換、増幅等の処理を
行った後、伝送信号として送信する。

【００４９】以下、アダプティブアレーアンテナのウエ
イト計算方法について説明する。本発明は、ウエイト計
算方法に特徴があり、ウエイト計算はコンピュータによ
るソフトウェアで行っても良いし、ウエイト計算の手順
を全てハードウェアで行っても良い。ハードウェアで行
う場合のウエイト計算システムは、前記図１乃至図４で
示される。

【００５０】まず、本実施の形態のウエイト計算器につ
いて述べる。図５は、本実施の形態のウエイト計算器３
１０の信号処理手順を示すフローチャートであり、複素
インパルス応答からウエイトを求める手順を示す。図
中、Ｓはフローの各ステップを示す。本実施の形態で
は、アレーアンテナ３０１のアンテナ素子の数Ｎ
_eleは、Ｎ_ele＝４である。

【００５１】入力信号は複素インパルス応答ｈ_e(k)（ｋ
は整数、ｅ＝１，２，…，Ｎ_ele）３１１と参照信号対
雑音電力比ＤＮ_e（ｅ＝１，２，…，Ｎ_ele）３１３であ
り、出力信号はウエイトＷ_e（ｅ＝１，２，…，Ｎ_ele）
３１２である（ステップＳ１００）。参照信号対雑音電
力比ＤＮ_eとは、従来技術で述べた方式Ａ（図１５参
照）において擬似送信信号“Ｄ(k)”２０２の電力と擬
似受信信号発生器２０３にて付加される雑音Ｎ_e(k)の電
力の比に相当するものである。ウエイト計算器３１０内
部で行う信号処理がスタートすると、まず、ステップＳ
１０１でｐ行ｑ列の要素Ｎ（ｐ，ｑ）が次式で表される
Ｎ_ele次の正方行列Ｎを計算する。

【００５２】

【数５３】 但しＤＮ_p ^-1はＤＮ_pの逆数を表す 次いで、ステップＳ１０２でｐ行ｑ列の要素Ｈ（ｐ，
ｑ）が次式で表されるＮ _ele次の正方行列Ｈを計算す
る。

【００５３】

【数５４】 但しｈ_p ^*(k)はｈ_p(k)の複素共役を表す すなわち、まずＮ_ele×Ｎ_ele次の行列ＮとＨを計算す
る。行列Ｎのｐ行ｑ列の要素をＮ（ｐ，ｑ）と表すと、
Ｎの対角成分は式Ｎ（ｅ，ｅ）＝１／ＤＮ_eを用いて算
出し、その他の成分は０とする。また、Ｈの要素は、
〔数５４〕に示すように算出する。但し、「ｈ_p ^*」はｈ
_pの複素共役を表す。ステップＳ１０３では、これらの
行列を用いて、次式（２）で示されるＮ_ele×Ｎ_ele次の
正方行列Ｒ_xxを計算する。 Ｒ_xx＝Ｈ＋Ｎ …（２） 次いで、ステップＳ１０４でｐ行の要素Ｒ_xd(p)が次式
で表されるＮ_ele次の列ベクトルＲ_xdを計算する。ここ
で、ｅ行目の要素をＲ_xd（ｅ）と表すと、式Ｒ_xd（ｅ）
＝ｈ_e ^*(0)を用いて算出する。

【００５４】

【数５５】 但しｈ_p ^*(0)はｈ_p(0)の複素共役を表す

【００５５】次いで、ステップＳ１０５でＲ_xxとＲ_xdを
用いて目的のウエイトＷ_e（ｅ＝１，２，…，Ｎ_ele）を
計算して本フローを終了する。具体的には、ウエイトＷ
_e（ｅ＝１，２，…，Ｎ_ele）をｅ行目の要素にもつＮ
_ele次の列ベクトルをＷとすると、連立一次方程式Ｒ_xx
・Ｗ＝Ｒ_xdを満足するＷを算出する。この方法として、
Gauss-Jordanの消去法等を用いることができる。

【００５６】本実施の形態のウエイト計算器３１０によ
れば、従来技術の方式Ａのウエイト計算器１１０に比べ
て計算手順が簡略化され、計算量を減少させることがで
きる。以上のウエイト計算の手順を全てハードウェアで
行う場合のウエイト計算器３１０のブロック図が図２で
あり、図５の信号処理手順に対応している。

【００５７】図２示すように、Ｒ_xdベクトル計算器４０
１では、複素インパルス応答ｈ_e(k)からＲ_xdベクトルを
計算する。Ｈ行列計算器４０２では、複素インパルス応
答ｈ _e(k)からＨ行列を計算する。Ｎ行列計算器４０３で
は、参照信号対雑音電力比ＤＮ_eからＮ行列を計算す
る。Ｒ_xx行列計算器４０４では、Ｈ行列とＮ行列からＲ
_{x x}行列を計算する。連立一次方程式計算器４０５では、
Ｒ_xx行列とＲ_xdベクトルからウエイトＷ_eを計算する。

【００５８】次に、本実施の形態で用いる伝送信号につ
いて述べる。伝送信号は、データチャネルと、それに比
べて電力が十分に抑圧されたパイロットチャネルから構
成される。これらは同じキャリア上にあり、同じ帯域幅
を持つ。パイロットチャネルは受信機において複素イン
パルス応答を測定し、アレーアンテナの指向性を制御す
るために用いる。

【００５９】前記図４に示す送信器６００のブロック図
を用いて伝送信号の発生方法を説明する。データチャネ
ル６０９は、２００Ｍｂｐｓの送信データ６０８を元
に、ＱＰＳＫ変調器６０３により生成する。一方、パイ
ロットチャネルはＰＮ系列発生器６０１，６０２により
生成したＭ系列を元に、ＱＰＳＫ変調器６０４により生
成したものを減衰器６０５で減衰して生成する。データ
チャネル、パイロットチャネル共に信号レートは１００
Ｍ Ｓｙｍｂｏｌｓ／ｓｅｃであり、これらを加算器６
０６で加算し、アップコンバート部６０７により帯域制
限、周波数変換、増幅等の処理を行った後、伝送信号と
して送信する。ベースバンドにおける帯域制限前の伝送
信号をＳ_B(t)、データチャネルの信号をＳ_Da(t)、パイ
ロットチャネルの信号をＳ_Pi(t)とすると、これらは次
式で示される。

【００６０】

【数５６】 ここで、ｄ_Ik、ｄ_Qkはデータチャネルのｋ番目のシンボ
ルの同相及び直交成分（±１）、Ｔ_sは１シンボル長
（＝１０ｎｓｅｃ）、δ(t)はディラックのデルタ関
数、ＭはＭ系列の周期（＝２５５）、βはデータチャネ
ルとパイロットチャネルの電力比（＝０．０４）、
ｍ_Ik、ｍ_Qkはパイロットチャネルの同相及び直交成分
（±１）でありそれぞれ別のＭ系列のｋ番目の値であ
る。

【００６１】次に、上記ウエイト計算器３１０以外の部
分の動作について説明する。前記図１において、アレー
アンテナ３０１はＮ_ele＝４本のアンテナ素子から構成
され、これらのアンテナで受信された信号は、ダウンコ
ンバート部３０２においてシンボルレートの４倍のサン
プリングレート４００ＭＳａｍｐｌｅｓ／ｓｅｃのデジ
タル複素ベースバンド信号３０３に変換される。

【００６２】ダウンコンバート部３０２からは参照信号
対雑音電力比ＤＮ_e（ｅ＝１，２，…，Ｎ_ele）３１３も
出力する。この値の決定方法については後述する。複素
インパルス応答測定器３０４では、入力したベースバン
ド信号３０３をパイロットチャネルに対応するマッチド
フィルタに通す。これによってまだデータチャネルが十
分抑圧されていない瞬時複素インパルス応答が、Ｍ系列
の周期の２５５シンボル毎に次々と得られる。これらの
連続する１６個について、同期加算を行い、データチャ
ネルの十分抑圧された複素インパルス応答を生成する。
次に生成した応答から最大電力を持つパスを探す。すな
わち、生成した応答をｈ_e（τ）とすると、

【００６３】

【数５７】 を計算し、σ（τ）を最大にするτ＝τ_maxを求める。
さらに、最大電力を持つパスを原点に移動し、シンボル
時間間隔で値を間引く。ここまでの結果をｈ_e(k)とする
と、ｈ_e(k)は次式（３）で示される。 ｈ_e(k)＝ｈ_e（τ_max ＋ ｋＴ_s） …（３） さらに、有効パス数をＬ＝１２とし、電力の大きなパス
から順番にＬ個分だけを選択し、後は値を０とする。す
なわち、

【００６４】

【数５８】 を計算し、σ(k)の大きなものから順番にＬ個分選択し
たものを、σ（ｋ₁），σ（ｋ₂），… ，σ（ｋ_L）とす
ると、ｋ₁，ｋ₂，…，ｋ_L以外のｋについて ｈ₁(k)＝ｈ₂(k)＝…＝ｈ_Nele(k)＝０ …（４） 式（４）に示す処理を行い、これを複素インパルス応答
として出力する。指向性制御器３０５では、入力したベ
ースバンド信号３０３にウエイトＷ_e３１２による重み
付けを行い、加算して出力することによって、アレーア
ンテナ３０１の指向性を制御する。つまり、入力ベース
バンド信号をＸ_e（ｅ＝１，２，…，Ｎ_ele）、出力をｙ
と表すと、次式で示される演算を行う。

【００６５】

【数５９】

【００６６】ダウンコンバート部３０２の動作を図３を
参照して説明する。ダウンコンバート部３０２は、アン
テナからの信号を入力し、増幅、周波数変換、Ａ／Ｄ変
換、フィルタリング等を行い、デジタル複素ベースバン
ド信号３０３を出力するものである。ダウンコンバート
部３０２の中には、ＡＧＣ部５０１がＮ_ele個あり、ア
ンテナ素子ｅに対応する部分から参照信号対雑音電力比
ＤＮ_eを出力する。ＤＮ_eは、以下の方法によって決定す
る。例としてｅ＝１について採り上げる。その他のｅに
ついても同様に決定できる。

【００６７】基本的な考え方は、アンテナ入力５０２に
遅延波がなくフェージングも受けていない伝送信号を入
力し、複素ベースバンド出力５０３におけるデータチャ
ネルの電力と雑音電力の比をＳＮＲ_e、複素インパルス
応答測定器３０４の出力ｈ_e(0)の値をＡ_e、とした時、
次式（５）を満足するようにＤＮ_eを出力するようにす
ればよい。

【００６８】 ＤＮ_e・|Ａ_e｜²＝ＳＮＲ_e …（５） 具体的には、送信機のＲＦ出力の信号を減衰器を通して
有線でアンテナ入力５０２に入力し、複素ベースバンド
出力５０３におけるデータチャネルの電力と雑音電力の
比をＳＮＲ１_e＝１００（２０ｄＢ）になるように減衰
器を調整し、この時のＡＧＣ部５０１の利得Ｇ１_e、複
素インパルス応答測定器３０４の出力ｈ_e(0)の値Ａ１_e
を測定する。これらの値を用いて、ＡＧＣ部５０１の利
得がＧ_eのときの出力すべきＤＮ_e（Ｇ_e）の値を、次式
（６）により決定する。 ＤＮ_e（Ｇ_e）＝（ＳＮＲ１_e・Ｇ１_e）／（｜Ａ１_e｜²・Ｇ_e） …（６）

【００６９】このようにしておけば、前記式（５）は満
足される。なぜならば、複素ベースバンド出力５０３に
おけるＡＧＣ部５０１の利得を除いた雑音電力をＮ１_e
とすると、上述した減衰器の設定において複素ベースバ
ンド出力５０３のパイロットチャネルの電力は、Ｎ１_e
・Ｇ１_e・ＳＮＲ１_e・βであり、複素インパルス応答測
定器３０４の利得は｜Ａ１_e｜²／（Ｎ１_e・Ｇ１_e・ＳＮ
Ｒ１_e・β）である。

【００７０】よって、一般的に利得がＧ_e、複素ベース
バンド出力５０２３におけるデータチャネル電力対雑音
電力比がＳＮＲ_eの時、複素インパルス応答測定器３０
４の出力ｈ_e(0)の値をＡ_eとすると、|Ａ_e｜²＝（Ｎ１_e
・Ｇ_e・ＳＮＲ_e・β）・| Ａ１_e｜²／（Ｎ１_e・Ｇ１_e・
ＳＮＲ１_e・β）であり、これを上記式（６）の両辺に
掛ければ、前記式（５）となるためである。

【００７１】以上述べたように、本実施の形態のウエイ
ト計算システムは、４本のアンテナ素子からなるアレー
アンテナ３０１、ＲＦ，ＩＦ，ダウンコンバータ，Ａ／
Ｄ変換器，ＬＰＦ等からなるダウンコンバート部３０
２、逆拡散処理及び同期加算処理を行うことによって、
伝搬路の複素インパルス応答を測定する複素インパルス
応答測定器３０４、伝搬路の複素インパルス応答ｈ
_e（ｋ）からアレーアンテナ３０１のウエイトＷ_eを計算
するウエイト計算器３１０と、ウエイトＷ_eを用いてア
ンテナの指向性を制御し、遅延波の抑圧された信号３０
６を復調部へと出力する指向性制御器３０５を備え、ウ
エイト計算器３１０は、あらかじめ設定された数Ｎ_ele
に対して、Ｎ_ele個の参照信号対雑音電力比ＤＮ_e（ｅ＝
１，２，…，Ｎ_{e le}）及びＮ_ele個の複素インパルス応答
ｈ_e(k)（ｅ＝１，２，…，Ｎ_ele）を入力し、まずＮ_ele
×Ｎ_ele次の行列ＮとＨを計算し、正方行列Ｎのｐ行ｑ
列の要素をＮ（ｐ，ｑ）と表すと、Ｎの対角成分は式Ｎ
（ｅ，ｅ）＝１／ＤＮ_eを用いて算出し、その他の成分
は０とする。また、Ｈの要素は、次式により算出する。
但し、Ｈ_p ^*はＨ_pの複素共役を表す。

【００７２】

【数６０】

【００７３】これらの行列を用いて、Ｒ_xx＝Ｈ＋Ｎなる
Ｎ_ele×Ｎ_ele次の行列Ｒ_xxを計算し、Ｎ_ele次の列ベク
トルＲ_xdを計算し、ｅ行目の要素をＲ_xd（e）と表す
と、式Ｒ _xd＝Ｈ_e ^*（０）を用いて算出する。そして、Ｒ
_xxとＲ_xdを用いて目的のウエイトＷ_e（ｅ＝１，２，
…，Ｎ_ele）を計算し、ウエイトＷ_e（ｅ＝１，２，…，
Ｎ_ele）をｅ行目の要素にもつＮ_ele次の列ベクトルを、
連立一次方程式Ｒ_xx・Ｗ＝Ｒ_{x d}により算出する。このた
め、従来技術の方式Ａのように統計処理を行ってウエイ
トを算出するウエイト計算方式では、ウエイトを高精度
に算出するために統計に用いるサンプル数を増大（例え
ば、方式Ａでは５０サンプル）させる必要があり、計算
量が大きくなる欠点があったが、本ウエイト計算では統
計処理行わずに高精度にウエイトを算出できるため、ウ
エイト算出のための計算量を格段に削減することができ
る。

【００７４】以下、この効果について詳細に説明する。
図６及び図７は、本実施の形態のウエイト計算器と従来
の方式Ａのウエイト計算器との計算量を示す図であり、
図６は、アンテナエレメント数Ｎ_eleが４個の場合、図
７は、アンテナエレメント数Ｎ_eleが７個の場合を示
す。

【００７５】横軸は、上記方式Ａにおいて計算に用いる
サンプル数である。本発明では統計的な計算を実際には
行わないためサンプル数の概念が無い。縦軸は、計算量
であり、複素数の加減乗除算、実数の加減乗除算すべて
を計算量１としている。パラメータＬは、入力する複素
インパルス応答のうち、有効なパスの数である。上記方
式Ａのウエイト計算アルゴリズムとしてはＲＬＳ方式、
ＳＭＩ方式があるが、これら両方の場合について計算量
を示した。

【００７６】図６及び図７に示すように、本発明によっ
て、計算量が大幅に削減されていることが分かる。一般
にウエイト計算アルゴリズムとしては、ＲＬＳ方式がよ
く用いられていること、ウエイト計算に用いるサンプル
数は、アンテナ素子の数の２倍以上、すなわちＮ_eleが
４、７の時、それぞれ８、１４以上が採用されているこ
とから、このような場合においては、図６及び図７から
本発明によって計算量が約９０％以上削減されているこ
とが分かる。

【００７７】第２の実施の形態 第１の実施の形態では、Ｎ_ele次の列ベクトルＷを計算
する手段として、連立一次方程式Ｒ_xx・Ｗ＝Ｒ_xdを満足
するＮ_ele次の列ベクトルＷを計算していた。本実施の
形態は、この部分を更に具体化した信号処理手順を提供
するものである。

【００７８】図８は、本発明にかかるウエイト計算をハ
ードウェアで行う場合のウエイト計算器の構成を示す図
である。本実施の形態のウエイト計算器は、前記図１の
ウエイト計算器３１０に置き換えたウエイト計算システ
ムである。本実施の形態の説明にあたり図２と同一構成
部分には同一符号を付している。

【００７９】図８において、ウエイト計算器７１０（ウ
エイト計算装置）は、複素インパルス応答ｈ_e(k)からＲ
_xdベクトルを計算するＲ_xdベクトル計算器４０１、複素
インパルス応答ｈ_e(k)からＨ行列を計算するＨ行列計算
器４０２、参照信号対雑音電力比ＤＮ_eからＮ行列を計
算するＮ行列計算器４０３、Ｈ行列とＮ行列からＲ_xx行
列を計算するＲ_xx行列計算器４０４、Ｒ_xx行列の逆行列
Ｒ_xx ^-1を計算する逆行列計算器７０１、及びＲ_xx ^-1行列
とＲ_xdベクトルからウエイトＷ_eを計算する行列乗算器
７０２から構成される。

【００８０】すなわち、本ウエイト計算器７１０は、前
記図２のＲ_xx行列とＲ_xdベクトルからウエイトＷ_eを計
算する連立一次方程式計算器４０５に換えて、Ｒ_xx行列
の逆行列Ｒ_xx ^-1を計算する逆行列計算器７０１及びＲ_xx
^-1行列とＲ_xdベクトルからウエイトＷ_eを計算する行列
乗算器７０２が設けられている。

【００８１】以下、本実施の形態のウエイト計算方法に
ついて説明する。本発明は、ウエイト計算方法に特徴が
あり、ウエイト計算はコンピュータによるソフトウェア
で行っても良いし、ウエイト計算の手順を全てハードウ
ェアで行っても良い。ハードウェアで行う場合のウエイ
ト計算システムは、前記図８で、ソフトウェアで行う場
合は図９で示される。

【００８２】図９は、本実施の形態のウエイト計算器７
１０の信号処理手順を示すフローチャートであり、複素
インパルス応答からウエイトを求める手順を示す。図５
に示すフローと同一処理手順には、同一ステップ番号を
付して重複部分の説明を省略する。図９のステップＳ１
０４でｐ行の要素Ｒ_xd(p)が次式で表されるＮ_ele次の列
ベクトルＲ_xdを計算すると、ステップＳ２０１に進む。
ステップＳ２０１では、次式に従ってＲ_xx行列の逆行列
Ｒ_xx ^-1を計算し、Ｒ_xx ^-1行列とＲ_xdベクトルからウエイ
トＷ_eを計算して本フローを終了する。

【００８３】

【数６１】 但しＲ_xx ^-1はＲ_xxの逆行列を表す

【００８４】次に、図９に示される本実施の形態のウエ
イト計算器７１０の信号処理手順を用いて算出したウエ
イトが、従来の方式Ａを用いて算出したウエイトと同等
のものであることを以下に示す。

【００８５】方式Ａの最小二乗法計算器２０５（図１
５）では、擬似送信信号ｄ(k)と擬似受信信号Ｘ_e(k)を
入力し、ＲＬＳ法を用いてウエイトＷ_eを計算してい
る。また、ＲＬＳ法の代わりにＳＭＩ法を用いることも
可能である。これらは結局、最小二乗法を計算するアル
ゴリズムである。最小二乗法は、Ｎ_ele×Ｎ_ele次のＲ_xx
行列

【００８６】

【数６２】 を、次式（７）によって計算し、 Ｒ_xx（ｐ，ｑ）＝Ｅ[Ｘ_p ^*(k)・Ｘq(k)] …（７） Ｎ_ele次の列ベクトルＲ_xd

【数６３】 を、次式（８）によって計算し、 Ｒ_xd(e)＝Ｅ[Ｘ_e ^*(k)・ｄ(k)] …（８）

【数６４】 によって計算し、これらを用いてウエイトＷ_eをｅ行目
の要素にもつＮ_ele次の列ベクトルＷを、次式（９）に
よって求めるものである。 Ｗ＝Ｒ_xx ^-1・Ｒ_xd …（９） 但し、Ｅ[・] はアンサンブル平均を表す。

【００８７】ＲＬＳ法やＳＭＩ法ではアンサンブル平均
を計算する代わりに、有限個のサンプルの平均、すなわ
ち方式Ａでは５０シンボルの時間平均を計算することに
よって最小二乗法を計算している。ところで、前記式
（１）には畳み込み演算が含まれるため、式（１）を上
記式（７）、式（８）に代入しても計算式が複雑になり
すぎてこれ以上計算手順は簡略化されない。

【００８８】そこで本実施の形態では、前記図１５に示
す擬似送信信号ｄ(k)２０２と擬似受信信号発生器２０
３にて付加される雑音ｎ_e(k)（ｅ＝１，２，…，
Ｎ_ele）の自己相関関数と相互相関関数の性質をここで
さらに導入し、計算手順を簡略化する。具体的には、擬
似送信信号ランダムな系列であるから、電力をＤとする
と、次式（１０）を満足し、

【００８９】

【数６５】 また擬似送信信号と雑音の関係、異なるアンテナ素子に
対応する雑音同士の関係は互いに無相関であるから、雑
音電力をＮ_eとすると、次式（１１）を満足する。

【００９０】

【数６６】 よって、前記式（１）を式（７）、式（８）に代入した
ものにさらにこれらの関係を導入すると、式（７）、式
（８）は、式（１２）、式（１３）のような形に簡略化
できる。

【００９１】

【数６７】 これらの簡略化の結果が本発明の計算手順となってい
る。上記式（１２）を計算する手順が図９のフローのス
テップＳ１０１、Ｓ１０２、Ｓ１０３であり、式（１
３）を計算する手順がステップＳ１０４となっている。

【００９２】本実施の形態の計算手順を用いたときの結
果と、従来技術の方式Ａの結果の違いは、方式Ａが５０
サンプルを用いた最小二乗法であるのに対して、本実施
の形態は無限個のサンプルを用いた最小二乗法となって
いる点である。これは、式（１０）、式（１１）の等式
において無限個のサンプルを考慮しているためである。
このように本発明では計算量が減少できると言う利点と
共に、性能的にもより良いものとなっている。

【００９３】以上述べたように、本実施の形態のウエイ
ト計算方法及び装置は、ウエイト計算器７１０が、連立
一次方程式Ｒ_xx・Ｗ＝Ｒ_xdに代えて、式Ｗ＝Ｒ_xx ^-1・Ｒ
_xd（但しＲ_xx ^-1はＲ_xdの逆行列を表す）で表されるＮ
_ele次の列ベクトルＷを計算するので、第１の実施の形
態と同様に、従来の方式Ａのウエイト計算器１１０に比
べて計算手順が簡略化され、計算量を大幅に減少させる
ことができる。特に、本実施の形態では計算量が減少で
きると言う利点と共に、性能的にもより良いものとなっ
ている。

【００９４】第３の実施の形態 第１及び第２の実施の形態では、行列式の絶対値にかか
わらず一律にウエイト計算を行っていた。本実施の形態
は、行列式Δの絶対値が閾値ＴＨ以下の場合は仮ウエイ
ト計算を行うようにしたものである。

【００９５】図１０は、本発明にかかるウエイト計算を
ハードウェアで行う場合のウエイト計算器の構成を示す
図である。本実施の形態のウエイト計算器は、前記図１
のウエイト計算器３１０に置き換えたウエイト計算シス
テムである。本実施の形態の説明にあたり図８と同一構
成部分には同一符号を付している。

【００９６】図１０において、ウエイト計算器８１０
（ウエイト計算装置）は、複素インパルス応答ｈ_e(k)か
らＲ_xdベクトルを計算するＲ_xdベクトル計算器４０１、
複素インパルス応答ｈ_e(k)からＨ行列を計算するＨ行列
計算器４０２、参照信号対雑音電力比ＤＮ_eからＮ行列
を計算するＮ行列計算器４０３、Ｈ行列とＮ行列からＲ
_{x x}行列を計算するＲ_xx行列計算器４０４、Ｒ_xx行列の逆
行列Ｒ_xx ^-1を計算する逆行列計算器７０１、Ｒ_xx ^-1行列
とＲ_xdベクトルからウエイトＷ_eを計算する行列乗算器
７０２、Ｒ_xx行列の行列式を計算する行列式計算器８０
１、行列式の絶対値が閾値以上のとき行列乗算器７０２
からの信号を出力し、閾値より小さいとき仮ウエイト計
算器８０３からの信号を出力する出力選択器８０２、及
び行列式の絶対値が閾値以下の時のウエイトを計算する
仮ウエイト計算器８０３から構成される。このように、
本ウエイト計算器８１０は、図８のウエイト計算器７１
０に行列式計算器８０１、出力選択器８０２及び仮ウエ
イト計算器８０３が追加された構成となっている。

【００９７】以下、本実施の形態のウエイト計算方法に
ついて説明する。上述したように、ウエイト計算はコン
ピュータによるソフトウェアで行っても良いし、ウエイ
ト計算の手順を全てハードウェアで行っても良い。ハー
ドウェアで行う場合のウエイト計算システムは、前記図
１０で、ソフトウェアで行う場合は図１１及び図１２で
示される。

【００９８】図１１及び図１２は、本実施の形態のウエ
イト計算器８１０の信号処理手順を示すフローチャート
であり、複素インパルス応答からウエイトを求める手順
を示す。図９に示すフローと同一処理手順には、同一ス
テップ番号を付して重複部分の説明を省略する。図１１
のステップＳ１０３でＮ_ele×Ｎ_ele次の正方行列Ｒ_xxを
計算すると、ステップＳ３０１に進む。

【００９９】ステップＳ３０１では、行列Ｒ_xxの行列式
Δを計算し、図１２のステップＳ３０２で行列式Δの絶
対値｜Δ｜があらかじめ設定された閾値ＴＨ以上か否か
を判別する。行列Ｒ_xxの行列式Δの絶対値｜Δ｜があら
かじめ設定された閾値ＴＨ以上のときは、ステップＳ１
０４以下に進んで前記図９の手順と同じ処理を行う。

【０１００】上記ステップＳ３０２で行列Ｒ_xxの行列式
Δの絶対値｜Δ｜があらかじめ設定された閾値ＴＨより
小さいときは、ステップＳ３０３で最大の｜Ｈ_e(0)｜を
持つｅ＝ｅ_maxを探す。次いで、ステップＳ３０４で目
的のウエイトＷ_e（ｅ＝１，２，…，Ｎ_ele）を計算して
本フローを終了する。具体的には、ｅ＝ｅ_max なるＷ_e
についてはＷ_e＝Ｈ_e ^*(0)／（｜Ｈ_e(0)｜² + Ｄ
Ｎ_e ^-1）、ｅ≠ｅ_max なるＷ_eについてはＷ_e＝０と決定
する。以上のウエイト計算の手順を全てハードウェアで
実現した場合のウエイト計算器８１０のブロック図が前
記図１０である。

【０１０１】本実施の形態のウエイト計算器８１０によ
れば、従来技術の方式Ａのウエイト計算器１１０に比べ
て計算手順が簡略化され、計算量を減少させることがで
きる。また、本実施の形態では、行列Ｒ_xxの行列式Δの
絶対値｜Δ｜と閾値ＴＨとを比較し、｜Δ｜が閾値ＴＨ
より小さいときは第２の実施の形態のウエイトＷ_e計算
を行わず、最大の｜Ｈ_e(0)｜を持つｅ_maxから直ちにウ
エイトＷ_eを決定しているので、この部分の計算量を減
少させることができる。

【０１０２】第４の実施の形態 上記各実施の形態に対して図５、図９及び図１１のステ
ップＳ１００の手順の代わりとなる手順を別例として示
す。

【０１０３】図１３は、図５、図９及び図１１の信号処
理手順の手順のステップＳ１００の代わりとなる手順を
示すステップＳ４００である。図１３において、ステッ
プＳ４００では、参照信号対雑音電力比ＤＮ_e（ｅ＝
１，２，…，Ｎ_ele）ではなく、その逆数ＤＮ_e ^-1（ｅ＝
１，２，…，Ｎ_ele）を用意するデータとしている。こ
れによって、図５、図９及び図１１の各ステップＳ１０
１において行っているＤＮ_e ^-1を算出する演算を省略す
ることができる。

【０１０４】なお、上記各実施の形態に係るウエイト計
算方法を、上述したようなビームフォーミング方式の送
受信機器に適用することもできるが、勿論これには限定
されず、無線通信を行う装置であればスペクトル拡散方
式を含む全ての通信装置に適用可能である。また、本ウ
エイト計算システムが、無線端末や送受信器のデータ通
信機能として組み込まれたものでもよい。

【０１０５】また、本実施の形態ではウエイト計算方
法，ウエイト計算システムという名称を用いたが、これ
は説明の便宜上であり、ビームフォーミング方式、アダ
プティブアレー制御方式等であってもよいことは勿論で
ある。また、上記各実施の形態において、アンテナ素子
の数及び種類、あるいは参照信号対雑音電力比ＤＮ_e、
複素インパルス応答ｈ_e（ｋ）等を適当な値に設定して
示したが、これらの設定値は式の導出過程で明らかなよ
うに任意に設定することができ、さまざまな規格のアダ
プティブアレーアンテナにおいて適用可能である。

【０１０６】さらに、上記ウエイト計算システムを構成
する各回路部、例えばダウンコンバート部、複素インパ
ルス応答測定器やウエイト計算器を構成する行列計算
器、乗算器の種類、数及び接続方法などは前述した実施
の形態に限られない。上記ウエイト計算の手順は、ハー
ドウェア又はソフトウェアのいずれにおいて実現しても
よい。

【０１０７】以上説明したウエイト計算システムでは、
この処理装置を機能させるためのプログラムでも実現さ
れる。このプログラムはコンピュータで読み取り可能な
記録媒体に格納されている。本発明では、この記録媒体
として、メインメモリそのものがプログラムメディアで
あってもよいし、また外部記憶装置としてプログラム読
み取り装置が設けられ、そこに記録媒体を挿入すること
で読み取り可能なプログラムメディアであってもよい。
いずれの場合においても、格納されているプログラムは
ＣＰＵがアクセスして実行させる構成であってもよい
し、あるいはいずれの場合もプログラムを読み出し、読
み出されたプログラムは、図示されていないプログラム
記憶エリアにダウンロードされて、そのプログラムが実
行される方式であってもよい。このダウンロード用のプ
ログラムは予め本体装置に格納されているものとする。

【０１０８】ここで、上記プログラムメディアは、本体
と分離可能に構成される記録媒体であり、例えばＰＣカ
ード（ＳＲＡＭカード）のほか、磁気テープやカセット
テープ等のテープ系、フロッピー（登録商標）ディスク
やハードディスク等の磁気ディスクやＣＤ−ＲＯＭ／Ｍ
Ｏ／ＭＤ／ＤＶＤ等の光ディスクのディスク系、ＩＣカ
ード（メモリカードを含む）／光カード等のカード系、
あるいはマスクＲＯＭ、ＥＰＲＯＭ、ＥＥＰＲＯＭ、フ
ラッシュＲＯＭ等による半導体メモリを含めた固定的に
プログラムを担持する媒体であってもよい。

【０１０９】また、本実施の形態においては、インター
ネットを含む通信ネットワークを接続可能なシステム構
成であることから、通信ネットワークからプログラムを
ダウンロードするように流動的にプログラムを担持する
媒体であってもよい。なお、このように通信ネットワー
クからプログラムをダウンロードする場合には、そのダ
ウンロード用のプログラムは予め本体装置に格納してお
くか、あるいは別な記録媒体からインストールされるも
のであってもよい。また、記録媒体に格納されている内
容としてはプログラムに限定されず、データであっても
よい。

【０１１０】

【発明の効果】以上、詳述したように、本発明によれ
ば、ウエイト算出のための計算量を従来の方式に比べ約
９０％以上に削減することのできるアダプティブアレー
アンテナのウエイト計算方法、ウエイト計算装置、ウエ
イト計算システムを実現することができる。

【図面の簡単な説明】

【図１】本発明の第１の実施の形態の受信器のアダプテ
ィブアレーアンテナの制御に関する部分の構成を示す図
である。

【図２】本実施の形態のウエイト計算方法をハードウェ
アで行う場合のウエイト計算器の構成を示す図である。

【図３】本実施の形態のウエイト計算システムのダウン
コンバート部の構成を示す図である。

【図４】本実施の形態で用いる伝送信号を送信する送信
器の構成を示す図である。

【図５】本実施の形態のウエイト計算器の信号処理手順
を示すフローチャートであり、

【図６】本実施の形態のウエイト計算器と従来の方式Ａ
のウエイト計算器との計算量を示す図である。

【図７】本実施の形態のウエイト計算器と従来の方式Ａ
のウエイト計算器との計算量を示す図である。

【図８】本発明の第２の実施の形態のウエイト計算方法
をハードウェアで行う場合のウエイト計算器の構成を示
す図である。

【図９】本実施の形態のウエイト計算器の信号処理手順
を示すフローチャートである。

【図１０】本発明の第３の実施の形態のウエイト計算方
法をハードウェアで行う場合のウエイト計算器の構成を
示す図である。

【図１１】本実施の形態のウエイト計算器の信号処理手
順を示すフローチャートである。

【図１２】本実施の形態のウエイト計算器の信号処理手
順を示すフローチャートである。

【図１３】図５、図９及び図１１の信号処理手順の手順
のステップＳ１００の代わりとなる手順を示すステップ
を示す図である。

【図１４】従来技術の方式Ａの受信器のアダプティブア
レーアンテナの制御に関する部分の構成を示す図であ
る。

【図１５】従来技術の方式Ａの受信器のウエイト計算器
の構成を示す図である。

【図１６】従来技術の方式Ａの受信器の信号処理のタイ
ムチャートを示す図である。

【符号の説明】

３０１ アレーアンテナ ３０２ ダウンコンバート部 ３０３ デジタル複素ベースバンド信号 ３０４ 複素インパルス応答測定器（インパルス応答測
定装置） ３０５ 指向性制御器（指向性制御装置） ３０６ 指向性制御器出力デジタル複素ベースバンド信
号 ３１０，７１０，８１０ ウエイト計算器（ウエイト計
算装置） ３１１ 複素インパルス応答 ３１２ ウエイト ３１３ 参照信号対雑音電力比 ４０１ Ｒ_xdベクトル計算器 ４０２ Ｈ行列計算器 ４０３ Ｎ行列計算器 ４０４ Ｒ_xx行列計算器 ４０５ 連立一次方程式計算器 ５０１ ＡＧＣ部でアンテナ素子１に対応するもの ５０２ アンテナ入力でアンテナ素子１に対応するもの ５０３ デジタル複素ベースバンド信号でアンテナ素子
１に対応するもの ６０１，６０２ ＰＮ系列発生器 ６０３ ＱＰＳＫ変調器 ６０４ ＱＰＳＫ変調器 ６０５ 減衰器 ６０６ 加算器 ６０７ アップコンバート部 ６０８ 送信データ ６０９ データチャネル ６１０ パイロットチャネル ７０１ 逆行列計算器 ７０２ 行列乗算器 ８０１ 行列式計算器 ８０２ 出力選択器 ８０３ 仮ウエイト計算器

Claims (15)

【特許請求の範囲】
1. 【請求項１】 アダプティブアレーアンテナを用いて受
   信した受信信号から伝搬路のインパルス応答を測定し、
   該伝搬路のインパルス応答から前記アンテナのウエイト
   を計算するウエイト計算方法において、 あらかじめ設定された数Ｎ_eleに対して、 Ｎ_ele個の参照信号対雑音電力比ＤＮ_e（ｅ＝１，２，
   …，Ｎ_ele）を入力し、 Ｎ_ele個の複素インパルス応答ｈ_e(k)（ｅ＝１，２，
   …，Ｎ_ele）を入力し、 ｐ行ｑ列の要素Ｎ（ｐ，ｑ）が次式で表されるＮ_ele次
   の正方行列Ｎを計算し、 【数１】 但しＤＮ_p ^-1はＤＮ_pの逆数を表す ｐ行ｑ列の要素Ｈ（ｐ，ｑ）が次式で表されるＮ_ele次
   の正方行列Ｈを計算し、 【数２】 但しｈ_p ^*(k)はｈ_p(k)の複素共役を表す 式Ｒ_xx＝Ｈ＋Ｎで表されるＮ_ele次の正方行列Ｒ_xxを計
   算し、 ｐ行の要素Ｒ_xd(p)が次式で表されるＮ_ele次の列ベクト
   ルＲ_xdを計算し、 Ｒ_xd(p)＝ｈ_p ^*(0) 但しｈ_p ^*(0)はｈ_p(0)の複素共役を表す 連立一次方程式Ｒ_xx・Ｗ＝Ｒ_xdを満足するＮ_ele次の列
   ベクトルＷを計算し、 前記列ベクトルＷのｐ行目の要素Ｗ(p)を、ウエイトＷ_p
   とし、 Ｎ_ele個のウエイトＷ_p（ｐ＝１，２，…，Ｎ_ele）を出
   力することを特徴とするウエイト計算方法。
2. 【請求項２】 アダプティブアレーアンテナを用いて受
   信した受信信号から伝搬路のインパルス応答を測定し、
   該伝搬路のインパルス応答から前記アンテナのウエイト
   を計算するウエイト計算方法において、 あらかじめ設定された数Ｎ_eleに対して、 Ｎ_ele個の参照信号対雑音電力比の逆数ＤＮ_e ^-1（ｅ＝
   １，２，…，Ｎ_ele）を入力し、 Ｎ_ele個の複素インパルス応答ｈ_e(k)（ｅ＝１，２，
   …，Ｎ_ele）を入力し、 ｐ行ｑ列の要素Ｎ（ｐ，ｑ）が次式で表されるＮ_ele次
   の正方行列Ｎを計算し、 【数３】 但しＤＮ_p ^-1はＤＮ_pの逆数を表す ｐ行ｑ列の要素Ｈ（ｐ，ｑ）が次式で表されるＮ_ele次
   の正方行列Ｈを計算し、 【数４】 但しｈ_p ^*(k)はｈ_p(k)の複素共役を表す 式Ｒ_xx＝Ｈ＋Ｎで表されるＮ_ele次の正方行列Ｒ_xxを計
   算し、ｐ行の要素Ｒ_xd(p)が次式で表されるＮ_ele次の列
   ベクトルＲ_xdを計算し、 Ｒ_xd(p)＝ｈ_p ^*(0) 但しｈ_p ^*(0)はｈ_p(0)の複素共役を表す 連立一次方程式Ｒ_xx・Ｗ＝Ｒ_xdを満足するＮ_ele次の列
   ベクトルＷを計算し、 前記列ベクトルＷのｐ行目の要素Ｗ(p)を、ウエイトＷ_p
   とし、 Ｎ_ele個のウエイトＷ_p（ｐ＝１，２，…，Ｎ_ele）を出
   力することを特徴とするウエイト計算方法。
3. 【請求項３】 アダプティブアレーアンテナを用いて受
   信した受信信号から伝搬路のインパルス応答を測定し、
   該伝搬路のインパルス応答から前記アンテナのウエイト
   を計算するウエイト計算方法において、 あらかじめ設定された数Ｎ_eleに対して、 Ｎ_ele個の参照信号対雑音電力比ＤＮ_e（ｅ＝１，２，
   …，Ｎ_ele）を入力し、 Ｎ_ele個の複素インパルス応答ｈ_e(k)（ｅ＝１，２，
   …，Ｎ_ele）を入力し、 ｐ行ｑ列の要素Ｎ（ｐ，ｑ）が次式で表されるＮ_ele次
   の正方行列Ｎを計算し、 【数５】 但しＤＮ_p ^-1はＤＮ_pの逆数を表す ｐ行ｑ列の要素Ｈ（ｐ，ｑ）が次式で表されるＮ_ele次
   の正方行列Ｈを計算し、 【数６】 但しｈ_p ^*(k)はｈ_p(k)の複素共役を表す 式Ｒ_xx＝Ｈ＋Ｎで表されるＮ_ele次の正方行列Ｒ_xxを計
   算し、 ｐ行の要素Ｒ_xd(p)が次式で表されるＮ_ele次の列ベクト
   ルＲ_xdを計算し、 Ｒ_xd(p)＝ｈ_p ^*(0) 但しｈ_p ^*(0)はｈ_p(0)の複素共役を表す 式Ｗ＝Ｒ_xx ^-1・Ｒ_xd（但しＲ_xx ^-1はＲ_xxの逆行列を表
   す）で表されるＮ_ele次の列ベクトルＷを計算し、 前記列ベクトルＷのｐ行目の要素Ｗ(p)を、ウエイトＷ_p
   とし、 Ｎ_ele個のウエイトＷ_p（ｐ＝１，２，…，Ｎ_ele）を出
   力することを特徴とするウエイト計算方法。
4. 【請求項４】 アダプティブアレーアンテナを用いて受
   信した受信信号から伝搬路のインパルス応答を測定し、
   該伝搬路のインパルス応答から前記アンテナのウエイト
   を計算するウエイト計算方法において、 あらかじめ設定された数Ｎ_eleに対して、 Ｎ_ele個の参照信号対雑音電力比の逆数ＤＮ_e ^-1（ｅ＝
   １，２，…，Ｎ_ele）を入力し、 Ｎ_ele個の複素インパルス応答ｈ_e(k)（ｅ＝１，２，
   …，Ｎ_ele）を入力し、 ｐ行ｑ列の要素Ｎ（ｐ，ｑ）が次式で表されるＮ_ele次
   の正方行列Ｎを計算し、 【数７】 但しＤＮ_p ^-1はＤＮ_pの逆数を表す ｐ行ｑ列の要素Ｈ（ｐ，ｑ）が次式で表されるＮ_ele次
   の正方行列Ｈを計算し、 【数８】 但しｈ_p ^*(k)はｈ_p(k)の複素共役を表す 式Ｒ_xx＝Ｈ＋Ｎで表されるＮ_ele次の正方行列Ｒ_xxを計
   算し、 ｐ行の要素Ｒ_xd(p)が次式で表されるＮ_ele次の列ベクト
   ルＲ_xdを計算し、 Ｒ_xd(p)＝ｈ_p ^*(0) 但しｈ_p ^*(0)はｈ_p(0)の複素共役を表す 式Ｗ＝Ｒ_xx ^-1・Ｒ_xd（但しＲ_xx ^-1はＲ_xxの逆行列を表
   す）で表されるＮ_ele次の列ベクトルＷを計算し、 前記列ベクトルＷのｐ行目の要素Ｗ(p)を、ウエイトＷ_p
   とし、 Ｎ_ele個のウエイトＷ_p（ｐ＝１，２，…，Ｎ_ele）を出
   力することを特徴とするウエイト計算方法。
5. 【請求項５】 アダプティブアレーアンテナを用いて受
   信した受信信号から伝搬路のインパルス応答を測定し、
   該伝搬路のインパルス応答から前記アンテナのウエイト
   を計算するウエイト計算方法において、 あらかじめ設定された数Ｎ_ele、閾値ＴＨに対して、 Ｎ_ele個の参照信号対雑音電力比ＤＮ_e（ｅ＝１，２，
   …，Ｎ_ele）を入力し、 Ｎ_ele個の複素インパルス応答ｈ_e(k)（ｅ＝１，２，
   …，Ｎ_ele）を入力し、 ｐ行ｑ列の要素Ｎ（ｐ，ｑ）が次式で表されるＮ_ele次
   の正方行列Ｎを計算し、 【数９】 但しＤＮ_p ^-1はＤＮ_pの逆数を表す ｐ行ｑ列の要素Ｈ（ｐ，ｑ）が次式で表されるＮ_ele次
   の正方行列Ｈを計算し、 【数１０】 但しｈ_p ^*(k)はｈ_p(k)の複素共役を表す 式Ｒ_xx＝Ｈ＋Ｎで表されるＮ_ele次の正方行列Ｒ_xxを計
   算し、 ｐ行の要素Ｒ_xd(p)が次式で表されるＮ_ele次の列ベクト
   ルＲ_xdを計算し、 Ｒ_xd(p)＝ｈ_p ^*(0) 但しｈ_p ^*(0)はｈ_p(0)の複素共役を表す 式Ｗ１＝Ｒ_xx ^-1・Ｒ_xd（但しＲ_xx ^-1はＲ_xxの逆行列を表
   す）で表されるＮ_ele次の列ベクトルＷ１を計算し、 ｅ＝１，２，…，Ｎ_eleの中から最大の｜Ｈ_e(0)｜（但
   し｜Ｈ_e(0)｜はＨ_e(0)の絶対値を表す）を持つｅ＝ｅ
   _maxを探し、 ｐ行の要素Ｗ０(p)が次式で表されるＮ_ele次の列ベクト
   ルＷ０を計算し、 【数１１】 Ｒ_xxの行列式Δを計算し、 行列式Δの絶対値｜Δ｜を前記閾値ＴＨと比較し、 ｜Δ｜≧ＴＨの場合、Ｗ１のｐ行目の要素Ｗ１(p)を、
   ウエイトＷ_pとし、 ｜Δ｜＜ＴＨの場合、Ｗ０のｐ行目の要素Ｗ０(p)を、
   ウエイトＷ_pとし、 Ｎ_ele個のウエイトＷ_p（ｐ＝１，２，…，Ｎ_ele）を出
   力することを特徴とするウエイト計算方法。
6. 【請求項６】 アダプティブアレーアンテナを用いて受
   信した受信信号から伝搬路のインパルス応答を測定し、
   該伝搬路のインパルス応答から前記アンテナのウエイト
   を計算するウエイト計算方法において、 あらかじめ設定された数Ｎ_ele、閾値ＴＨに対して、 Ｎ_ele個の参照信号対雑音電力比の逆数ＤＮ_e ^-1（ｅ＝
   １，２，…，Ｎ_ele）を入力し、 Ｎ_ele個の複素インパルス応答ｈ_e(k)（ｅ＝１，２，
   …，Ｎ_ele）を入力し、 ｐ行ｑ列の要素Ｎ（ｐ，ｑ）が次式で表されるＮ_ele次
   の正方行列Ｎを計算し、 【数１２】 但しＤＮ_p ^-1はＤＮ_pの逆数を表す ｐ行ｑ列の要素Ｈ（ｐ，ｑ）が次式で表されるＮ_ele次
   の正方行列Ｈを計算し、 【数１３】 但しｈ_p ^*(k)はｈ_p(k)の複素共役を表す 式Ｒ_xx＝Ｈ＋Ｎで表されるＮ_ele次の正方行列Ｒ_xxを計
   算し、 ｐ行の要素Ｒ_xd(p)が次式で表されるＮ_ele次の列ベクト
   ルＲ_xdを計算し、 Ｒ_xd(p)＝ｈ_p ^*(0) 但しｈ_p ^*(0)はｈ_p(0)の複素共役を表す 式Ｗ１＝Ｒ_xx ^-1・Ｒ_xd（但しＲ_xx ^-1はＲ_xxの逆行列を表
   す）で表されるＮ_ele次の列ベクトルＷ１を計算し、 ｅ＝１，２，…，Ｎ_eleの中から最大の｜Ｈ_e(0)｜（但
   し｜Ｈ_e(0)｜はＨ_e(0)の絶対値を表す）を持つｅ＝ｅ
   _maxを探し、 ｐ行の要素Ｗ０(p)が次式で表されるＮ_ele次の列ベクト
   ルＷ０を計算し、 【数１４】 Ｒ_xxの行列式Δを計算し、 行列式Δの絶対値｜Δ｜を前記閾値ＴＨと比較し、 ｜Δ｜≧ＴＨの場合、Ｗ１のｐ行目の要素Ｗ１(p)を、
   ウエイトＷ_pとし、 ｜Δ｜＜ＴＨの場合、Ｗ０のｐ行目の要素Ｗ０(p)を、
   ウエイトＷ_pとし、 Ｎ_ele個のウエイトＷ_p（ｐ＝１，２，…，Ｎ_ele）を出
   力することを特徴とするウエイト計算方法。
7. 【請求項７】 アダプティブアレーアンテナを用いて受
   信した受信信号から伝搬路のインパルス応答を測定し、
   該伝搬路のインパルス応答から前記アンテナのウエイト
   を計算するウエイト計算装置において、 あらかじめ設定された数Ｎ_eleに対して、 Ｎ_ele個の参照信号対雑音電力比ＤＮ_e（ｅ＝１，２，
   …，Ｎ_ele）を入力し、 Ｎ_ele個の複素インパルス応答ｈ_e(k)（ｅ＝１，２，
   …，Ｎ_ele）を入力し、 ｐ行ｑ列の要素Ｎ（ｐ，ｑ）が次式で表されるＮ_ele次
   の正方行列Ｎを計算する手段と、 【数１５】 但しＤＮ_p ^-1はＤＮ_pの逆数を表す ｐ行ｑ列の要素Ｈ（ｐ，ｑ）が次式で表されるＮ_ele次
   の正方行列Ｈを計算する手段と、 【数１６】 但しｈ_p ^*(k)はｈ_p(k)の複素共役を表す 式Ｒ_xx＝Ｈ＋Ｎで表されるＮ_ele次の正方行列Ｒ_xxを計
   算する手段と、 ｐ行の要素Ｒ_xd(p)が次式で表されるＮ_ele次の列ベクト
   ルＲ_xdを計算する手段と、 Ｒ_xd(p)＝ｈ_p ^*(0) 但しｈ_p ^*(0)はｈ_p(0)の複素共役を表す 連立一次方程式Ｒ_xxＷ＝Ｒ_xdを満足するＮ_ele次の列ベ
   クトルＷを計算する手段と、 前記列ベクトルＷのｐ行目の要素Ｗ(p)を、ウエイトＷ_p
   とし、 Ｎ_ele個のウエイトＷ_p（ｐ＝１，２，…，Ｎ_ele）を出
   力する出力手段とを備えることを特徴とするウエイト計
   算装置。
8. 【請求項８】 アダプティブアレーアンテナを用いて受
   信した受信信号から伝搬路のインパルス応答を測定し、
   該伝搬路のインパルス応答から前記アンテナのウエイト
   を計算するウエイト計算装置において、 あらかじめ設定された数Ｎ_eleに対して、 Ｎ_ele個の参照信号対雑音電力比の逆数ＤＮ_e ^-1（ｅ＝
   １，２，…，Ｎ_ele）を入力し、 Ｎ_ele個の複素インパルス応答ｈ_e(k)（ｅ＝１，２，
   …，Ｎ_ele）を入力し、 ｐ行ｑ列の要素Ｎ（ｐ，ｑ）が次式で表されるＮ_ele次
   の正方行列Ｎを計算する手段と、 【数１７】 但しＤＮ_p ^-1はＤＮ_pの逆数を表す ｐ行ｑ列の要素Ｈ（ｐ，ｑ）が次式で表されるＮ_ele次
   の正方行列Ｈを計算する手段と、 【数１８】 但しｈ_p ^*(k)はｈ_p(k)の複素共役を表す 式Ｒ_xx＝Ｈ＋Ｎで表されるＮ_ele次の正方行列Ｒ_xxを計
   算する手段と、 ｐ行の要素Ｒ_xd(p)が次式で表されるＮ_ele次の列ベクト
   ルＲ_xdを計算する手段と、 Ｒ_xd(p)＝ｈ_p ^*(0) 但しｈ_p ^*(0)はｈ_p(0)の複素共役を表す 連立一次方程式Ｒ_xxＷ＝Ｒ_xdを満足するＮ_ele次の列ベ
   クトルＷを計算する手段と、 前記列ベクトルＷのｐ行目の要素Ｗ(p)を、ウエイトＷ_p
   とし、 Ｎ_ele個のウエイトＷ_p（ｐ＝１，２，…，Ｎ_ele）を出
   力する出力手段とを備えることを特徴とするウエイト計
   算装置。
9. 【請求項９】 アダプティブアレーアンテナを用いて受
   信した受信信号から伝搬路のインパルス応答を測定し、
   該伝搬路のインパルス応答から前記アンテナのウエイト
   を計算するウエイト計算装置において、 あらかじめ設定された数Ｎ_eleに対して、 Ｎ_ele個の参照信号対雑音電力比ＤＮ_e（ｅ＝１，２，
   …，Ｎ_ele）を入力し、 Ｎ_ele個の複素インパルス応答ｈ_e(k)（ｅ＝１，２，
   …，Ｎ_ele）を入力し、 ｐ行ｑ列の要素Ｎ（ｐ，ｑ）が次式で表されるＮ_ele次
   の正方行列Ｎを計算する手段と、 【数１９】 但しＤＮ_p ^-1はＤＮ_pの逆数を表す ｐ行ｑ列の要素Ｈ（ｐ，ｑ）が次式で表されるＮ_ele次
   の正方行列Ｈを計算する手段と、 【数２０】 但しｈ_p ^*(k)はｈ_p(k)の複素共役を表す 式Ｒ_xx＝Ｈ＋Ｎで表されるＮ_ele次の正方行列Ｒ_xxを計
   算する手段と、 ｐ行の要素Ｒ_xd(p)が次式で表されるＮ_ele次の列ベクト
   ルＲ_xdを計算する手段と、 Ｒ_xd(p)＝ｈ_p ^*(0) 但しｈ_p ^*(0)はｈ_p(0)の複素共役を表す 式Ｗ＝Ｒ_xx ^-1・Ｒ_xd（但しＲ_xx ^-1はＲ_xxの逆行列を表
   す）で表されるＮ_ele次の列ベクトルＷを計算する手段
   と、 前記列ベクトルＷのｐ行目の要素Ｗ(p)を、ウエイトＷ_p
   とし、 Ｎ_ele個のウエイトＷ_p（ｐ＝１，２，…，Ｎ_ele）を出
   力する出力手段とを備えることを特徴とするウエイト計
   算装置。
10. 【請求項１０】 アダプティブアレーアンテナを用いて
    受信した受信信号から伝搬路のインパルス応答を測定
    し、該伝搬路のインパルス応答から前記アンテナのウエ
    イトを計算するウエイト計算装置において、 あらかじめ設定された数Ｎ_eleに対して、 Ｎ_ele個の参照信号対雑音電力比の逆数ＤＮ_e ^-1（ｅ＝
    １，２，…，Ｎ_ele）を入力し、 Ｎ_ele個の複素インパルス応答ｈ_e(k)（ｅ＝１，２，
    …，Ｎ_ele）を入力し、 ｐ行ｑ列の要素Ｎ（ｐ，ｑ）が次式で表されるＮ_ele次
    の正方行列Ｎを計算する手段と、 【数２１】 但しＤＮ_p ^-1はＤＮ_pの逆数を表す ｐ行ｑ列の要素Ｈ（ｐ，ｑ）が次式で表されるＮ_ele次
    の正方行列Ｈを計算する手段と、 【数２２】 但しｈ_p ^*(k)はｈ_p(k)の複素共役を表す 式Ｒ_xx＝Ｈ＋Ｎで表されるＮ_ele次の正方行列Ｒ_xxを計
    算する手段と、 ｐ行の要素Ｒ_xd(p)が次式で表されるＮ_ele次の列ベクト
    ルＲ_xdを計算する手段と、 Ｒ_xd(p)＝ｈ_p ^*(0) 但しｈ_p ^*(0)はｈ_p(0)の複素共役を表す 式Ｗ＝Ｒ_xx ^-1・Ｒ_xd（但しＲ_xx ^-1はＲ_xxの逆行列を表
    す）で表されるＮ_ele次の列ベクトルＷを計算する手段
    と、 前記列ベクトルＷのｐ行目の要素Ｗ(p)を、ウエイトＷ_p
    とし、 Ｎ_ele個のウエイトＷ_p（ｐ＝１，２，…，Ｎ_ele）を出
    力する出力手段とを備えることを特徴とするウエイト計
    算装置。
11. 【請求項１１】 アダプティブアレーアンテナを用いて
    受信した受信信号から伝搬路のインパルス応答を測定
    し、該伝搬路のインパルス応答から前記アンテナのウエ
    イトを計算するウエイト計算装置において、 あらかじめ設定された数Ｎ_ele、閾値ＴＨに対して、 Ｎ_ele個の参照信号対雑音電力比ＤＮ_e（ｅ＝１，２，
    …，Ｎ_ele）を入力し、 Ｎ_ele個の複素インパルス応答ｈ_e(k)（ｅ＝１，２，
    …，Ｎ_ele）を入力し、 ｐ行ｑ列の要素Ｎ（ｐ，ｑ）が次式で表されるＮ_ele次
    の正方行列Ｎを計算する手段と、 【数２３】 但しＤＮ_p ^-1はＤＮ_pの逆数を表す ｐ行ｑ列の要素Ｈ（ｐ，ｑ）が次式で表されるＮ_ele次
    の正方行列Ｈを計算する手段と、 【数２４】 但しｈ_p ^*(k)はｈ_p(k)の複素共役を表す 式Ｒ_xx＝Ｈ＋Ｎで表されるＮ_ele次の正方行列Ｒ_xxを計
    算する手段と、ｐ行の要素Ｒ_xd(p)が次式で表されるＮ
    _ele次の列ベクトルＲ_xdを計算する手段と、 Ｒ_xd(p)＝ｈ_p ^*(0) 但しｈ_p ^*(0)はｈ_p(0)の複素共役を表す 式Ｗ１＝Ｒ_xx ^-1・Ｒ_xd但しＲ_xx ^-1はＲ_xxの逆行列を表
    す）で表されるＮ_ele次の列ベクトルＷ１を計算する手
    段と、 ｅ＝１，２，…，Ｎ_eleの中から最大の｜Ｈ_e(0)｜（但
    し｜Ｈ_e(0)｜はＨ_e(0)の絶対値を表す）を持つｅ＝ｅ
    _maxを探す手段と、 ｐ行の要素Ｗ０(p)が次式で表されるＮ_ele次の列ベクト
    ルＷ０を計算する手段と、 【数２５】 Ｒ_xxの行列式Δを計算する手段と、 行列式Δの絶対値｜Δ｜を前記閾値ＴＨと比較する手段
    と、 ｜Δ｜≧ＴＨの場合、Ｗ１のｐ行目の要素Ｗ１(p)を、
    ウエイトＷ_pとし、 ｜Δ｜＜ＴＨの場合、Ｗ０のｐ行目の要素Ｗ０(p)を、
    ウエイトＷpとし、 Ｎ_ele個のウエイトＷ_p（ｐ＝１，２，…，Ｎ_ele）を出
    力する出力手段とを備えることを特徴とするウエイト計
    算装置。
12. 【請求項１２】 アダプティブアレーアンテナを用いて
    受信した受信信号から伝搬路のインパルス応答を測定
    し、該伝搬路のインパルス応答から前記アンテナのウエ
    イトを計算するウエイト計算装置において、 あらかじめ設定された数Ｎ_ele、閾値ＴＨに対して、 Ｎ_ele個の参照信号対雑音電力比の逆数ＤＮ_e ^-1（ｅ＝
    １，２，…，Ｎ_ele）を入力し、 Ｎ_ele個の複素インパルス応答ｈ_e(k)（ｅ＝１，２，
    …，Ｎ_ele）を入力し、 ｐ行ｑ列の要素Ｎ（ｐ，ｑ）が次式で表されるＮ_ele次
    の正方行列Ｎを計算する手段と、 【数２６】 但しＤＮ_p ^-1はＤＮ_pの逆数を表す ｐ行ｑ列の要素Ｈ（ｐ，ｑ）が次式で表されるＮ_ele次
    の正方行列Ｈを計算する手段と、 【数２７】 但しｈ_p ^*(k)はｈ_p(k)の複素共役を表す 式Ｒ_xx＝Ｈ＋Ｎで表されるＮ_ele次の正方行列Ｒ_xxを計
    算する手段と、p行の要素Ｒ_xd(p)が次式で表されるＮ
    _ele次の列ベクトルＲ_xdを計算する手段と、 Ｒ_xd(p)＝ｈ_p ^*(0) 但しｈ_p ^*(0)はｈ_p(0)の複素共役を表す 式Ｗ１＝Ｒ_xx ^-1・Ｒ_xd（但しＲ_xx ^-1はＲ_xxの逆行列を表
    す）で表されるＮ_ele次の列ベクトルＷ１を計算する手
    段と、 ｅ＝１，２，…，Ｎ_eleの中から最大の｜Ｈ_e(0)｜（但
    し｜Ｈ_e(0)｜はＨ_e(0)の絶対値を表す）を持つｅ＝ｅ
    _maxを探す手段と、 ｐ行の要素Ｗ０(p)が次式で表されるＮ_ele次の列ベクト
    ルＷ０を計算する手段と、 【数２８】 Ｒ_xxの行列式Δを計算する手段と、 行列式Δの絶対値｜Δ｜を前記閾値ＴＨと比較する手段
    と、 ｜Δ｜≧ＴＨの場合、Ｗ１のｐ行目の要素Ｗ１(p)を、
    ウエイトＷ_pとし、 ｜Δ｜＜ＴＨの場合、Ｗ０のｐ行目の要素Ｗ０(p)を、
    ウエイトＷ_pとし、 Ｎ_ele個のウエイトＷ_p（ｐ＝１，２，…，Ｎ_ele）を出
    力する出力手段とを備えることを特徴とするウエイト計
    算装置。
13. 【請求項１３】 アダプティブアレーアンテナを用いて
    受信した受信信号から伝搬路のインパルス応答を測定す
    るインパルス応答測定装置と、 伝搬路のインパルス応答から前記アンテナのウエイトを
    計算するウエイト計算装置と、 計算されたウエイトを用いて前記アンテナの指向性を制
    御し、遅延波の抑圧された信号を復調部へと出力する指
    向性制御装置とを備えるウエイト計算システムにおい
    て、 前記ウエイト計算装置は、請求項７乃至１２のいずれか
    一項に記載のウエイト計算装置を用いることを特徴とす
    るウエイト計算システム。
14. 【請求項１４】 アダプティブアレーアンテナを用いて
    受信した受信信号から伝搬路のインパルス応答を測定
    し、該伝搬路のインパルス応答から前記アンテナのウエ
    イトを計算するウエイト計算方法において、 あらかじめ設定された数Ｎ_eleに対して、 Ｎ_ele個の参照信号対雑音電力比ＤＮ_e（ｅ＝１，２，
    …，Ｎ_ele）を入力し、 Ｎ_ele個の複素インパルス応答ｈ_e(k)（ｅ＝１，２，
    …，Ｎ_ele）を入力し、 ｐ行ｑ列の要素Ｎ（ｐ，ｑ）が次式で表されるＮ_ele次
    の正方行列Ｎを計算し、 【数２９】 但しＤＮ_p ^-1はＤＮ_pの逆数を表す ｐ行ｑ列の要素Ｈ（ｐ，ｑ）が次式で表されるＮ_ele次
    の正方行列Ｈを計算し、 【数３０】 但しｈ_p ^*(k)はｈ_p(k)の複素共役を表す 式Ｒ_xx＝Ｈ＋Ｎで表されるＮ_ele次の正方行列Ｒ_xxを計
    算し、 ｐ行の要素Ｒ_xd(p)が次式で表されるＮ_ele次の列ベクト
    ルＲ_xdを計算し、 Ｒ_xd(p)＝ｈ_p ^*(0) 但しｈ_p ^*(0)はｈ_p(0)の複素共役を表す 連立一次方程式Ｒ_xx・Ｗ＝Ｒ_xdを満足するＮ_ele次の列
    ベクトルＷを計算し、 前記列ベクトルＷのｐ行目の要素Ｗ(p)を、ウエイトＷ_p
    とし、 Ｎ_ele個のウエイトＷ_p（ｐ＝１，２，…，Ｎ_ele）を出
    力するウエイト計算方法をコンピュータに実行させるた
    めのプログラムを記録したことを特徴とするコンピュー
    タ読み取り可能な記録媒体。
15. 【請求項１５】 アダプティブアレーアンテナを用いて
    受信した受信信号から伝搬路のインパルス応答を測定
    し、該伝搬路のインパルス応答から前記アンテナのウエ
    イトを計算するウエイト計算方法において、 あらかじめ設定された数Ｎ_ele、閾値ＴＨに対して、 Ｎ_ele個の参照信号対雑音電力比ＤＮ_e（ｅ＝１，２，
    …，Ｎ_ele）を入力し、 Ｎ_ele個の複素インパルス応答ｈ_e(k)（ｅ＝１，２，
    …，Ｎ_ele）を入力し、 ｐ行ｑ列の要素Ｎ（ｐ，ｑ）が次式で表されるＮ_ele次
    の正方行列Ｎを計算し、 【数３１】 但しＤＮ_p ^-1はＤＮ_pの逆数を表す ｐ行ｑ列の要素Ｈ（ｐ，ｑ）が次式で表されるＮ_ele次
    の正方行列Ｈを計算し、 【数３２】 但しｈ_p ^*(k)はｈ_p(k)の複素共役を表す 式Ｒ_xx＝Ｈ＋Ｎで表されるＮ_ele次の正方行列Ｒ_xxを計
    算し、 ｐ行の要素Ｒ_xd(p)が次式で表されるＮ_ele次の列ベクト
    ルＲ_xdを計算し、 Ｒ_xd(p)＝ｈ_p ^*(0) 但しｈ_p ^*(0)はｈ_p(0)の複素共役を表す 式Ｗ１＝Ｒ_xx ^-1・Ｒ_xd（但しＲ_xx ^-1はＲ_xxの逆行列を表
    す）で表されるＮ_ele次の列ベクトルＷ１を計算し、 ｅ＝１，２，…，Ｎ_eleの中から最大の｜Ｈ_e(0)｜（但
    し｜Ｈ_e(0)｜はＨ_e(0)の絶対値を表す）を持つｅ＝ｅ
    _maxを探し、 ｐ行の要素Ｗ０(p)が次式で表されるＮ_ele次の列ベクト
    ルＷ０を計算し、 【数３３】 Ｒ_xxの行列式Δを計算し、 行列式Δの絶対値｜Δ｜を前記閾値ＴＨと比較し、 ｜Δ｜≧ＴＨの場合、Ｗ１のｐ行目の要素Ｗ１(p)を、
    ウエイトＷ_pとし、 ｜Δ｜＜ＴＨの場合、Ｗ０のｐ行目の要素Ｗ０(p)を、
    ウエイトＷ_pとし、 Ｎ_ele個のウエイトＷ_p（ｐ＝１，２，…，Ｎ_ele）を出
    力するウエイト計算方法をコンピュータに実行させるた
    めのプログラムを記録したことを特徴とするコンピュー
    タ読み取り可能な記録媒体。