JP5867968B2 - 無線受信装置及び無線受信方法並びに無線システム - Google Patents

無線受信装置及び無線受信方法並びに無線システム Download PDF

Info

Publication number
JP5867968B2
JP5867968B2 JP2013526893A JP2013526893A JP5867968B2 JP 5867968 B2 JP5867968 B2 JP 5867968B2 JP 2013526893 A JP2013526893 A JP 2013526893A JP 2013526893 A JP2013526893 A JP 2013526893A JP 5867968 B2 JP5867968 B2 JP 5867968B2
Authority
JP
Japan
Prior art keywords
covariance matrix
eigenvector
signal
matrix
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.)
Active
Application number
JP2013526893A
Other languages
English (en)
Other versions
JPWO2013018716A1 (ja
Inventor
次夫 丸
次夫 丸
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
NEC Corp
Original Assignee
NEC Corp
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by NEC Corp filed Critical NEC Corp
Priority to JP2013526893A priority Critical patent/JP5867968B2/ja
Publication of JPWO2013018716A1 publication Critical patent/JPWO2013018716A1/ja
Application granted granted Critical
Publication of JP5867968B2 publication Critical patent/JP5867968B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04B—TRANSMISSION
    • H04B7/00—Radio transmission systems, i.e. using radiation field
    • H04B7/02—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas
    • H04B7/04—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas
    • H04B7/08—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas at the receiving station
    • H04B7/0837—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas at the receiving station using pre-detection combining
    • H04B7/0842—Weighted combining
    • H04B7/0848—Joint weighting
    • H04B7/0854—Joint weighting using error minimizing algorithms, e.g. minimum mean squared error [MMSE], "cross-correlation" or matrix inversion

Landscapes

  • Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Radio Transmission System (AREA)

Description

本発明は固定局間のマイクロ波無線通信システムのマイクロ波無線受信に関し、特にパイロット信号やリファレンス信号或いはトレーニング信号と呼ばれる既知の系列信号を用いない手法であるブラインド処理によるスペースダイバーシチに関する。
マイクロ波無線通信システムでは、フェージングの防止策としてスペースダイバーシチ方式（Space Diversity：以下、ＳＤ方式と呼ぶ。）を採用することが有効である。
ここで、ＳＤ方式とは、適切な距離をおいて配置された複数のアンテナへの受信波が同時に劣化する可能性が低いことを利用して、フェージングによる電波伝搬特性の劣化を補償する技術である。ＳＤ方式に関して図１０及び図１１を参照して具体的に説明する。
図１０はＳＤ方式を採用する場合の具体的な実装例である。図１０を参照すると、本具体例は、固定局５１００と固定局５２００の二つの固定局を有している。
そして、固定局５１００から送信された信号は固定局５２００で二基のアンテナにより受信され、固定局５２００での受信処理によってダイバーシチ処理が施される。
図１１を参照すると固定局５２００での受信処理部が示されている。
受信処理は、図１１に示す様に、各アンテナ（アンテナ５２１１、アンテナ５２１２）に対応して設けられた受信データ処理があり、処理された各ブランチの信号を合成することでダイバーシチ処理が完了する。ここでブランチとは各アンテナに対応して設けられた受信データ処理用の回路や処理のことを指す。図１１の例でいえば、ブランチ１は、アンテナ５２１１、ミキサ５２２１及びＡ／Ｄ変換器５２３１により実現される。また、ブランチ２は、アンテナ５２１２、ミキサ５２２２及びＡ／Ｄ変換器５２３２により実現される。また、ダイバーシチ合成部５２４０により各ブランチの信号が合成される。なお、図１１に示すのは本具体例の説明において特に重要な部分のみであり、フィルタ等の他の構成要素や直交復調用の９０度位相成分に関しては図示を省略する。
このような、ＳＤ方式に関する技術としては、例えば特許文献１に記載の技術等が挙げられる。
特開平０７−０３８５８３６号公報
上述のように、ＳＤ方式はフェージングの防止策として一般的に用いられている。近年はトラフィックの増加に伴いマイクロ波無線通信システムの大容量化が進み、高次の多値変調信号を高速に伝送するようになってきており、一般的なＳＤ方式では対応できない問題が生じてきた。この問題について説明する。
マイクロ波帯では、例えば、図１２に示す様に、多重伝搬路（干渉波の発生要因）による周波数選択性フェージングが発生する。高速で更に高次の多値変調信号を使用すると、この周波数選択性フェージングにより誤り率特性が著しく劣化する。
この誤り率特性の著しい劣化を防止するために従来から周波数選択性フェージング対策として適応等化処理が行われてきた。一般に適応アルゴリズムではタップ係数の収束特性を得る為にトレーニング信号を伝送する。特に時分割多元接続(ＴＤＭＡ：Time Division Multiple Access)の移動通信などのバースト伝送では、トレーニング信号の伝送が必須である。
一方、固定マイクロ波無線通信システムなどでは、既存システムの信号構成上、トレーニング信号等を入れられない場合がある。この場合は、パイロット信号やリファレンス信号或いはトレーニング信号と呼ばれる既知の系列信号を用いない手法であるブラインド処理が有効である。
また、例えパイロット信号やリファレンス信号或いはトレーニング信号を入れられる場合であったとしても、これらの信号は、既知の系列信号で実際の情報の伝達に何ら寄与しない為、伝送効率を上げるためにもブラインド処理を行うことが好ましい。
具体的に説明すると、移動通信等では、１０％程度の無線リソースをパイロット信号やリファレンス信号に割り当てている。これをそのまま１０２４ＱＡＭで伝送している固定マイクロ波無線通信システムに当てはめたとすると、５１２ＱＡＭ相当の伝送速度になってしまう。すなわち、折角送信電力を倍増して高速伝送を実現したのに、情報伝達に対する実効速度がそれに伴わないことになるのである。
このような前提のもとで、ブラインド処理で実現したＳＤ方式について更に検討する。ブラインド処理で実現したＳＤ方式としては、同相合成（In-Phase：以下「ＩＰ合成」と呼ぶ。）と呼ばれる手法がある。
ＩＰ合成法は、ブランチの片方の信号に低周波の位相変調をかけて、合成後のＡＭ信号成分により位相の同相状態を検出し、その検出結果を用いて位相器へのフィードバック制御を行うというものである。
ＩＰ合成法は、低い周波数の振幅成分を検出に用いている。この点、高次の多値変調が可能なＱＡＭはその変調方式に起因する振幅成分の変動を内包している。その為、多値のＱＡＭへのＩＰ合成法の適用は問題がある。また後述するが、ＩＰ合成法は、周波数選択性フェージングにより誤り率特性が著しく劣化するという問題点もある。つまり、ＩＰ合成法は如何なる場合であっても有効な手法であるとはいえない。
他方、周波数選択性フェージングを考慮し、且つ、ブラインド処理で実現したＳＤ方式として最小振幅偏差合成（Minimum Dispersion Combiner：以下、「ＭＩＤ合成」と呼ぶ。）と呼ばれる方法も存在する。
ここで、ＩＰ合成法と異なる部分で、ＭＩＤ合成法の特徴を示す干渉波相殺合成動作について、説明のし易さから極端な例を用いて説明する。
図１３は、ＭＩＤ合成法の特徴を示す干渉波相殺合成動作を二つのアンテナで受信した例で示した図である。異なる電波伝搬条件の二つの例（図１３中の、例１及び例２）で説明している。また図１３では受信波を回転フェザーのベクトル表記で示している。
図１３、ブランチ１（Ｂｒａｎｃｈ １）は、例えば図１２における受信局である基地局５００２の上側のアンテナに対する受信信号を示しているものとする。また、同様に図１３、ブランチ２（Ｂｒａｎｃｈ ２）は、図１２における受信局である基地局５００２の下側のアンテナに対する受信信号を示しているということにする。
各ブランチにおける受信波は図１２に示されている様に、直接波（図１２中では、「ＬＯＳ ｒａｙ」と表記する、図１３中では実線にて表記する。）に遅延反射波（図１２中、「Ｒｅｆｌｅｃｔｅｄ ｒａｙ」と表記する図１３中では破線にて表記する。）が加わったものである。
干渉波相殺合成動作は、干渉波であるそれぞれのブランチの遅延反射波の位相を逆位相になるようにブランチ内にある位相器を調整して合成を行うものである。従って、図１３上側例１の右側に示されている様に、干渉波成分であるブランチ１の遅延反射波成分とブランチ２の遅延反射波成分が相殺され、残るブランチ１の直接波成分とブランチ２の直接波成分が合成されダイバーシチ処理が完結する。
図１３上側例１のようにこの処理が理想的に行われたならば、遅延反射波成分が完全に無くなるので、周波数選択性フェージングによる干渉が結果的に解消される。しかしながらこの方式は以下の問題を包含している。
マイクロ波無線通信システムの電波伝搬環境は、アンテナの設置状況、反射対象である大地、あるいは海峡を跨ぐ場合には海面反射の状況によって色々と変わる。上述した図１３の上側の例１は幸運な電波伝搬環境の場合である。つまり、干渉波成分である遅延反射波成分を相殺する様にブランチの位相器を調整した結果各ブランチの直接波が同じ方向を向くようになり、合成後の希望波電力が増大してダイバーシチ利得が得られた幸運な例である。一方、図１３の下側の例２は厳しい条件の場合である。干渉波成分である遅延反射波成分を相殺するようにブランチの位相器を調整した結果各ブランチの直接波が反対方向に近い状態になってしまっている。その結果、合成後の希望波電力が減少して、周波数選択性は解消されてもダイバーシチ利得が逆に減ってしまった例である。
以上は説明のし易さから極端な形で、ＭＩＤ合成法を説明したが、ここまで極端でないとしても電波伝搬環境によってダイバーシチ効果が大きく異なるという点でＭＩＤ合成法には問題があることは明らかである。
そこで、本発明は、実際の情報の伝達に何ら寄与しないパイロット信号やリファレンス信号或いはトレーニング信号など伝送効率を低下させる既知の系列を用いることのないブラインド処理を実現し、且つ、高速で高次の多値変調を用いる大容量マイクロ波無線通信に厳しい遅延干渉波が存在する周波数選択性フェージング環境下であったとしても、無線伝搬環境による影響の少ない常に良好なダイバーシチ効果を発揮することが可能な無線受信装置及び無線受信方法並びに無線システムを提供することを目的とする。
本発明の第１の観点によれば、複数のアンテナで受信した受信信号毎に対応するブランチを設け、当該各ブランチの受信信号を合成することによりスペースダイバーシチを実現する無線受信装置において、前記各ブランチの受信信号を重み付け係数に基づくウェイト処理によりダイバーシチ合成する際に、前記ウェイト処理対象とする受信信号の受信信号ベクトルから得られた共分散行列が単位行列となる様な修正行列を生成し、当該修正行列によって前記受信信号を変換し、当該変換後の受信信号の受信信号ベクトルから得られた共分散行列の最大固有値に対応した固有ベクトルを算出し、当該固有ベクトルを前記重み付け係数として前記ダイバーシチ合成を行うことを特徴とする無線受信装置が提供される。
本発明の第２の観点によれば、複数のアンテナで受信した受信信号毎に対応するブランチを設け、当該各ブランチの受信信号を合成することによりスペースダイバーシチを実現する無線受信装置において、前記各ブランチの受信信号から干渉波信号を検出する干渉波信号検出部と、前記干渉波信号から得られた共分散行列が単位行列となる様な修正行列を生成し、前記受信信号の受信信号ベクトルから得られた共分散行列を前記修正行列により修正し、当該修正後の共分散行列の最大固有値に対応した固有ベクトルを用いて重み付け係数を生成し、当該該重み付け係数に基づくウェイト処理によりダイバーシチ合成を行うブラインド処理部と、を備えることを特徴とする無線受信装置が提供される。
本発明の第３の観点によれば、信号を無線送信する無線送信装置と、前記送信された信号を複数のアンテナにより受信する無線受信装置と、を含んだ無線通信システムにおいて、
前記無線送信装置と前記無線受信装置間の無線通信はブラインド処理を実現した通信方式に準拠しており、前記無線受信装置は上記本発明の第１の観点又は第２の観点により提供される記載の無線受信装置である、ことを特徴とする無線通信システムが提供される。
本発明の第４の観点によれば、複数のアンテナで受信した受信信号毎に対応するブランチを設け、当該各ブランチの受信信号を合成することによりスペースダイバーシチを実現する無線受信装置が行う無線受信方法において、前記各ブランチの受信信号を重み付け係数に基づくウェイト処理によりダイバーシチ合成する際に、前記ウェイト処理対象とする受信信号の受信信号ベクトルから得られた共分散行列が単位行列となる様な修正行列を生成し、当該修正行列によって前記受信信号を変換し、当該変換後の受信信号の受信信号ベクトルから得られた共分散行列の最大固有値に対応した固有ベクトルを算出し、当該固有ベクトルを前記重み付け係数として前記ダイバーシチ合成を行うことを特徴とする無線受信方法が提供される。
本発明の第５の観点によれば、複数のアンテナで受信した受信信号毎に対応するブランチを設け、当該各ブランチの受信信号を合成することによりスペースダイバーシチを実現する無線受信装置が行う無線受信方法において、前記各ブランチの受信信号から干渉波信号を検出し、前記干渉波信号から得られた共分散行列が単位行列となる様な修正行列を生成し、前記受信信号の受信信号ベクトルから得られた共分散行列を前記修正行列により修正し、当該修正後の共分散行列の最大固有値に対応した固有ベクトルを用いて重み付け係数を生成し、当該該重み付け係数に基づくウェイト処理によりダイバーシチ合成を行う、ことを特徴とする無線受信方法が提供される。
本発明によれば、受信信号ベクトルから得られた共分散行列の最大固有値に対応した固有ベクトルを算出して、該固有ベクトルを重み付け係数として用いることから、伝送効率を低下させる既知の系列を用いることのないブラインド処理を実現し、且つ、高速で高次の多値変調を用いる大容量マイクロ波無線通信に厳しい遅延干渉波が存在する周波数選択性フェージング環境下であったとしても、無線伝搬環境による影響の少ない常に良好なダイバーシチ効果を発揮することが可能となる。
固有ベクトル空間で表したウェイトベクトルについて表す図である。 最大固有値の固有ベクトル 本発明の実施形態１における固有ベクトル合成方式のシミュレーション結果を表す図である。 白色化により特性劣化を緩和した固有ベクトルによるダイバーシチ合成を行う実施形態である、本発明の実施形態２の基本的構成を表す図である。 本発明の実施形態２における白色化による固有ベクトル合成方式のシミュレーション結果を表す図である。 本発明の干渉波白色化により劣化を緩和した固有ベクトル合成を行う実施形態である、本発明の実施形態３の基本的構成を表す図である。 本発明の実施形態２及び３における白色化固有ベクトル合成方式のシミュレーション結果を表す図（１／２）である。 本発明の実施形態２及び３における白色化固有ベクトル合成方式のシミュレーション結果を表す図（２／２）である。 図７及び８で用いた伝搬環境のベクトル表記とその時のタップ係数表す図（１／２）である。 図７及び８で用いた伝搬環境のベクトル表記とその時のタップ係数表す図（２／２）である。 ＳＤ方式を用いたマイクロ波無線通信システムの具体例の構成を表す図である。 ＳＤ方式を用いたマイクロ波無線通信システムの具体例における受信処理部の構成を表す図である。 干渉波の発生要因である多重伝搬路について説明するための図である。 ＭＩＤ合成法について説明するための図である。
１１１、１１２、５２１１、５２１２ アンテナ
１２１、１２２、５２２１、５２２２ ミキサ
１３１、１３２、５２３１、５２３２ Ａ／Ｄ変換器
１４１、１４２ タップ処理部
１５０、５２５０ 局部発振器
１６０ ブラインド処理部
１７０ ダイバーシチ合成部
２１１、２２２ 白色化・タップ処理部
３１０ 干渉波信号検出部
３２０ 白色化ブラインド処理部
１０００、２０００、３０００ 受信装置
５００１ 固定局（受信）
５００２ 固定局（送信）
次に、本発明の実施形態について図面を参照して詳細に説明する。
最初に、本発明の実施形態の概略について述べる。
上述で説明した様に、一般的なブラインド処理でＩＰ（In-Phase）合成と呼ばれる方式がある。このＩＰ合成方式は片方の信号に低周波の位相変調をかけ、合成後のＡＭ信号により位相のフィードバック制御を行う必要があり、低周波の振幅成分を用いている。そのため定振幅ではない多値変調のＱＡＭにはこのままでは使えない。すなわち、ＩＰ合成方式はＰＳＫ（phase shift keying）の様にＡＭ変動の無い定振幅変調に適した方法であって、高密度変調には適さない。そこで、本発明の実施形態であるマイクロ波無線受信装置は、別のアプローチを採用して高密度変調でも使える様にした点に特徴がある。
更に、一般的なＭＩＤ合成法の特徴を示す干渉波相殺合成は、確かに周波数選択性は解消されたがダイバーシチ利得が逆に減ってしまった。また、電波伝搬環境によってダイバーシチ効果が大きく異なるという点でＭＩＤ合成法には問題がある。本発明の実施形態であるマイクロ波無線受信装置は、これらの問題に対処した点に特徴がある。
以下の説明では、先ず本発明の各実施形態の原理的な側面を定量的に説明する為に、定量解析から始めてブラインド処理の概念で構築する。そして、本実施形態における定量解析を示す。最後に、具体的例を示し、その具体例に対して本実施形態の定量的解析結果を用いたシミュレーションにより、各実施形態の効果を定量的に示す。
［実施形態１］
まず、本発明の基本的構成を示す実施形態である実施形態１について説明する。実施形態１は、最大固有値の固有ベクトルによるダイバーシチ合成方式に関する実施形態である。
最初に最大固有値の固有ベクトルによるダイバーシチ合成方式の原理について説明する。
以下の説明では、送信信号をｓ（ｔ）、ウェイトベクトルをＷ、チャネルベクトルをｈ、雑音ベクトルをｎとする。この場合、信号対雑音電力比（ＳＮＲ：signal-to-noise ratio）γ（以下、適宜「ＳＮＲγ」と呼ぶ）は、
・・・（式１）
と表せる。ここで、
∵
とした。
ＲｓとＲｎとはそれぞれ共分散行列を表す。具体的には、Ｒｓは送信信号の共分散行列を表す。また、Ｒｎは雑音の共分散行列を表す。また、集合平均に対して、統計的独立性から、
とした。
Ｉは単位行列である。送信信号の共分散行列であるＲｓは階数（以下、適宜「ランク」と呼ぶ）１である。一方、雑音の共分散行列であるＲｎはフルランクである。
ＳＮＲγは（式１）より分母及び分子のそれぞれにウェイトベクトルｗが含まれている。そのため、ＳＮＲγはウェイトベクトルｗの大きさには関係しない。そこで、ＳＮＲの最大化問題を次の拘束条件の下で考える。
・・・（式２）
これより、（式１）の分母は
となり一定となる。従って、ＳＮＲγの最大化は、分子のみの最大化に帰着する。
ここで、本実施形態ではブラインド型を用いることを想定している。そのため、受信側でチャネルベクトルｈを求めることが出来ない。そこで評価関数の要素として下記の数式により表される入力信号ベクトルｘを用いることにする。
この入力信号ベクトルｘにウェイトベクトルｗを掛けた
に対する評価関数として、
・・・（式３）
を採用することにする。
上記の（式３）における、Ｒは入力信号の共分散行列であり、
・・・（式４）
即ち、評価関数Ｊの最大化はｗＨ・Ｒｓ・ｗの最大化となり、（式１）の分子の最大化となり、ＳＮＲγの最大化と同じになる。
［最適化問題］
ｗＨ・ｗ＝１の制約条件のもとで、評価関数Ｊを最大化する解ｗを調べる。
ラグランジュの未定乗数法（Lagrange-Multiplier Method）を用いて上述の（式３）を変形すると、
ｗＨに関して偏微分を取ることによって極解が求まるから、
・・・（式５）
λＬはラグランジュの未定乗数で、解ｗは次の条件を満たすことになる。
・・・（式６）
これは固有ベクトルと固有値の定義そのものに他ならない。よって、解ｗは入力信号の共分散行列Ｒの固有ベクトルであり、ラグランジュの未定乗数は固有値となる。
一方、エルミート（Hermitian）行列である共分散行列Ｒは、その正規化された固有ベクトルβｊと、その固有値λｊによってスペクトル分解すると、
・・・（式７）
となり固有ベクトルはＮ個存在する。このＮ個の中から最大固有値λ１に対応する固有ベクトルβ１を求めればそれがｗ＝β１となる。
今、任意のウェイトベクトルｗに対して、（式３）の評価関数は、（式７）を用いて、
・・・（式８）
ここで、Ｒはエルミート行列であるから、その固有値λｊは全て必ず実数で非負である。
そこで、
・・・（式９）
を満足するウェイトベクトルｗに対してスカラー値（直交系である｛βｎ｝への射影係数）
・・・（式１０）
とすると、
で（式８）は、（式９）を条件として、以下の様になる。
・・・（式１１）
この式１１は超楕円を示す。即ちｗを固有ベクトルにより直交展開した
で、その直交軸βｎに対する射影係数が
となって、直交系｛βｎ｝での幾何学的関係を示している。
従って、例えば一つの固有ベクトルが形作る軸に対する端点を見ると、極限の値であるから（式１１）における他の軸成分は全て零となり、全ての値ｃ２をこの軸に集めた形となり、その固有ベクトルをβ１とすると、
、
に対し、
となり、
に対するβ１軸上の端点となる。即ち、
を満足するウェイトベクトルは固有ベクトルで形成された軸に対してその軸の固有値の平方根に逆比例した値で与えられる超楕円上に存在することになる。
以上の内容を表現のし易さからウェイトベクトルｗの要素は複素数であるが実数のように扱って描くことにする。
直感的理解の説明のし易さから二次元の実数の要素の場合で示した簡単な例を図１に示す。
図１より、
上の範囲は
で、軸が固有ベクトルで与えられる楕円となる。
で
は、
となってβ１方向での最適な値を示し、
は
となってβ２方向での値を示す。
即ち固有値の平方根に反比例する長さで軸が固有ベクトルで与えられる楕円上の
がｃ２値を持つ等高線となる。制約条件ｗＨ・ｗ＝１は図１中の円筒となる。最大固有値を与える固有ベクトルを求める問題は、等高線ｃ２の値を変えることによって出来る超放物曲線とこの円筒が交わる曲線上で、
が最大となるｗ＝β１を求める問題と同じことになる。
次に、最大固有値に対応する固有ベクトルの算出手段の一例について説明する。
本例は、冪乗法による最大固有値に対応する固有ベクトル算出法である。
（式７）式より、入力信号の共分散行列Ｒは
で表される。固有ベクトルβｊは正規化されている。そのｋ乗をとると
、且つ、
であるから、
となる。ここで、全てのｊに対して
ならば、
で、
ｋを十分に大きくすると
第二項は消えて、
となる。即ちＲｋはランク１の行列に限りなく近づく。
一方、βと直交関係にない任意のベクトルをｖ０とする。この時新しいベクトルｖ１を
とする。
もし、Ｒのランクが１であるなら、
であるから、
となり、ｖ１は固有ベクトルβ１に比例する。
今、ｖ０とＲｋ・ｖ０を以下の様に正規化された直交基底ベクトルβによって展開する。
、
全てのｊに対して
λ１＞λｊ
１≠ｊならば、
でｋを十分に大きくすると第二項は消えて、
となり、最大固有値に対応する固有ベクトルが得られる。この時、ｖ０がβ１と直交しているとα１が零となってしまうので、βと直交関係にない任意のベクトルをｖ０とする必要がある。
また、上式よりλ１／λｊの比率で解が収束することが分かる。
以上の原理に基づいて、次のアルゴリズムによる処理手段が得られる。
［処理手段］
離散時刻（ｉ−１）の時点で、入力信号の共分散行列Ｒ（ｉ−１）の推定値が既に求められているとする。時刻ｉにおける初期値として入力信号ベクトル
を用いて共分散行列
の推定を
雑音が存在しない場合の共分散行列は（式１）下よりランク１の
で
に対応して固有ベクトルがｈに相当するから、例として、
が推定出来たとして、
を用いると、正規化した
は、
として与えられる。
上式において、λは忘却係数である、１≧λ＞０である。
なお、一般に、離散時刻ｉ−１の推定が無いときは以下の様に設定する。
ｋ＝１
その後、正の小さな実数εに対して
が成立している間、以下のステップを繰り返す。
Ｓｔｅｐ１・・・
Ｓｔｅｐ２・・・
ｋ＝ｋ＋１Ｓｔｅｐ１へ
条件が成立すると、次式を実行して繰り返しを終了する。
収束条件
を満たした時、
は、時刻ｉにおける共分散行列
の最大固有ベクトルとなる。
このアルゴリズムは、最大固有値に対する他の固有値（主に第２最大固有値）の比率で収束速度が決まる。雑音が無いときには、ランク１となるので、収束が最も速い。
以上説明した様に、複数のブランチの受信信号を重み付け係数に基づいてウェイト処理により合成する際、受信信号ベクトルから得られた共分散行列Ｒの最大固有値λに対応した固有ベクトルβを共分散行列Ｒから算出し、当該算出した固有ベクトルβを前記重み付け係数として用いた時、ダイバーシチ合成後の信号対雑音電力比（ＳＮＲ）γが最大となり良好なダイバーシチ利得が得られる。
しかも、上述の方法はトレーニングシーケンスなどの伝送効率を低下させる信号を使わないブラインド処理であることに注意を要する。すなわち、上述の方法では、トレーニング信号などの伝送効率を低下させる信号を使うことなく、良好なダイバーシチ利得を得ることが出来る。
更に、上述の方法は共分散行列の最大固有値に対応した固有ベクトルβを共分散行列Ｒから算出する際、共分散行列Ｒの冪乗演算によって最大固有値に対する固有ベクトルを算出し、得られた固有ベクトルを重み付け係数として用いたことを特徴としている。よって、固有ベクトル算出の為に、固有値分解等の複雑な方法ではなく、簡単な積和演算で実現出来、回路規模を小さく出来ることが分かる。
上述した動作を実行するための受信装置の構成例を図２に示す。
図２を参照すると本実施形態の受信装置１０００は、アンテナ１１１、アンテナ１１２、ミキサ１２１、ミキサ１２２、Ａ／Ｄ変換器１３１、Ａ／Ｄ変換器１３２、タップ処理部１４１、タップ処理部１４２、局部発振器１５０、ブラインド処理部１６０及びダイバーシチ合成部１７０を有する。
図２において、上側の系統（アンテナ１１１、ミキサ１２１、Ａ／Ｄ変換器１３１及びタップ処理部１４１）がブランチ１である。一方、図２において、下側の系統（アンテナ１１２、ミキサ１２２、Ａ／Ｄ変換器１３２及びタップ処理部１４２）がブランチ２である。なお、ブランチ１及びブランチ２を構成するこれら各部を具体的にどのような回路や部材により実現するかは本実施形態の要旨ではない。各部は任意の方式及び任意の形状により実現することが可能である。また、図２に示すのは本実施形態の説明において特に重要な部分のみであり、フィルタ等の他の構成要素や直交復調用の９０度位相成分に関しては図示を省略する。
続いて、受信装置１０００が有する各部の動作について説明する。まず、アンテナ１１１及び１１２にて信号が受信される。そして、アンテナ１１１及び１１２でそれぞれ受信した信号は、波数変換用のミキサ１２１、ミキサ１２２に入力される。
周波数変換用のミキサ１２１、ミキサ１２２は、局部発振器（local oscillator）１５０で生成した信号を用いることでアンテナ１１１及び１１２から入力された信号をそれぞれダウンコンバートする。なお、図２では１つの局部発振器１５０を両ブランチで共有する例で書かれているが、ブランチ毎に別々の局部発振器を用いるようにしてもよい。
そして、周波数変換用のミキサ１２１、ミキサ１２２でダウンコンバートされた信号は、アナログ−デジタル変換器であるＡ／Ｄ変換器１３１、１３２でそれぞれデジタル信号へと変換される。
Ａ／Ｄ変換器１３１、１３２で変換された信号は、重み付け係数に基づいてウェイトの乗算を行うタップ処理部１４１、１４２へそれぞれ入力される。
ここで、本実施形態のブラインド処理部１６０による重み付け計数の算出について具体的に説明する。
上述した、共分散行列の演算の式［数４２］から［数４８］の処理を行うのがブラインド１６０処理部である。式［数４２］のＸに相当するのが、本実施形態の場合は、
となる。これはブランチ１からブラインド１６０処理部へ入力される受信信号ｘ１とブランチ２からブラインド１６０処理部へ入力される受信信号ｘ２の両受信信号から信号ベクトルＸを作ったことに相当する。また、処理の結果である式［数４８］は、この具体例の場合は、
となる。これは、上述の説明から受信信号ベクトルＸの共分散行列Ｒの最大固有に対応する固有ベクトルである。ブラインド処理部１６０の入力は［数５０］式のｘ１、ｘ２である。また、ブラインド処理部１６０の出力は、［数５１］式のｗ１、ｗ２である。
そして、タップ処理部１４１はブラインド処理部１６０から送られてきたｗ１を基に
なる演算を行う。
タップ処理部１４１の出力の同相成分がｙ１Ｉ、直交成分がｙ１Ｑである。
また、同様にタップ処理部１４２はブラインド処理部１６０から送られてきたｗ２を基に
なる演算を行う。タップ処理部１４２の出力の同相成分がｙ２Ｉ、直交成分がｙ２Ｑである。
この各タップ処理部によって乗算されたタップ出力は、ダイバーシチ合成部１７０へ送られる。そして、ダイバーシチ合成部１７０においてブランチ１と２の信号が合成され、ダイバーシチ処理が完結する。即ち上述の［数７］式に相当する操作が行われることとなる。
具体的に記載すると、
となる。
上述した本実施形態では、固有ベクトルβを重み付け係数として用いているので、ダイバーシチ合成後の信号対雑音電力比（ＳＮＲ）γが最大となり良好なダイバーシチ利得が得られるという効果を奏する。
なお、ブラインド処理部１６０に入力する同相成分ｙ１Ｉ直交成分ｙ１Ｑや同相成分ｙ２Ｉ直交成分ｙ２Ｑはモニタ用である。
上述した本実施形態の効果をシミュレーションにより確認した結果が図３に示されている。図３では、縦軸はシンボル誤り率、横軸はアンテナ当たりのＳＮＲで、チャネルは直接波のみの伝搬環境、多値変調として６４ＱＡＭを使用している。
図３中の「ｅ：本発明の固有ベクトル合成」が本実施形態を用いた場合のシミュレーション結果に相当する。
なお、図３には本実施形態との比較の為、一般的な技術で説明したＩＰ合成法、等化機能を持たせたもの、ダイバーシチ機能が無いもの等の特性も載せている。
具体的に検討すると、それぞれの中で本実施形態を用いた場合の固有ベクトル合成が一番良い特性を示している。
一方、本実施形態に肉薄する「ｄ：ＩＰ合成法」は、背景技術として説明した様に多値のＱＡＭの適用には問題がある。そこでこのシミュレーションでは位相が分かっているものと仮定して行った理想状態の特性であり、実現性には問題がある。
しかし、本実施形態の固有ベクトルは完全ブラインド処理で行われている。比較対象「ａ：ダイバーシチ無し」は、ダイバーシチ無しでＣＨ１のみの場合の特性であり、この特性との比較から本実施形態の固有ベクトル合成ｅは良好なダイバーシチ利得を有していることがわかる。
他方、比較対象「ｂ：ダイバーシチ無しでＬＭＳ等化機能有り」はダイバーシチ無しであって、ＬＭＳ等化機能がある場合である。このシミュレーション環境は直接波のみで大地反射波等の干渉波が存在しない。従って等化機能を必要としない。同じ比較対象のａと比べて特性が劣化しているのは、等化機能を司るＬＭＳステップサイズμによる適応処理で最適点からのμによる揺らぎの為である。同様のことは、「ｃ：ダイバーシチ有りで、等化機能有り」でも言える。ｃは、ダイバーシチ有りで、等化機能を有する比較対象である。ｃは、ダイバーシチ利得を得ているもののμの揺らぎの為に特性が劣化している。
このように、上述したａ〜ｄの何れの方法も本実施形態の固有ベクトル合成ｅより特性が劣っている。すなわち、本実施形態は一般的な技術と比べて、有利な効果を奏することが分かる。
［実施形態２］
上述した実施形態１の特性は、直接波のみで大地反射波等の干渉波が存在しない伝搬環境での特性である。上述したように実施形態１の固有ベクトル合成のダイバーシチは一般的な技術で説明したＩＰ合成法や等化機能を持たせたものよりも特性が良かった。しかし大地反射波等の遅延干渉波が存在する伝播環境では急激に特性が劣化する。そこで、実施形態２として、白色化処理による特性劣化を緩和する実施形態を説明する。すなわち、実施形態２は、白色化処理により特性劣化を緩和して最大固有値の固有ベクトルにより合成した方式に関するものである。
まず、白色化処理の原理について説明する。
今、各ブランチの受信信号ベクトルをＸｉとする。その共分散行列をＲｉｉとすると、
そこで、
なる修正行列Ａｉにより、例えば
とすると、
となる。即ち、共分散行列Ｒｉｉを単位行列にするように修正行列Ａｉは作用する。
そこで、受信信号ベクトルＸｉをこの修正行列Ａｉで変換すると、ベクトル
の要素間は無相関となる。例えばブランチ１の修正行列によって変換された受信信号をｘ１とし、ブランチ２の修正行列によって変換された受信信号をｘ２とする。受信信号ベクトルＸは、
となる。これは上述のブランチ１の受信信号ｘ１とブランチ２の受信信号ｘ２から受信信号ベクトルＸを作ったことと同じである。これを基にウェイトを計算すると、上述より処理の結果である例えば式［数０４８］を用いて
が得られる。ダイバーシチ合成後の信号ｙは
となる。
上述した動作を実行するための構成例を図４に示す。
図４を参照すると、本実施形態の受信装置２０００は、アンテナ１１１、アンテナ１１２、ミキサ１２１、ミキサ１２２、Ａ／Ｄ変換器１３１、Ａ／Ｄ変換器１３２、白色化・タップ処理部２１１、白色化・タップ処理部２１２、局部発振器１５０、ブラインド処理部１６０及びダイバーシチ合成部１７０を有する。
図２に示される実施形態１の構成例と共通する符号が振られている部分での処理は上述と同じなので説明を省略する。
一方、本実施形態特有の箇所である白色化・タップ処理部２１１と白色化・タップ処理部２１２については、以下説明する。
ブランチ１の白色化・タップ処理部２１１を代表として説明する。この点、ブランチ２の白色化・タップ処理部２１２における処理は白色化・タップ処理部２１１と同様であるので説明を省略する。
ブランチ１の白色化・タップ処理部２１１は、式［数５５］を用いて共分散行列Ｒ１１を計算する。実際の処理は例えば［数４２］を用いれば良い。そして［数５６］の関係から修正行列Ａ１を算出する。更に、［数５７］を用いて修正行列Ａ１により変換された受信信号ｘ１を算出し、ブラインド処理部１６０へ送る。ブランチ２も同様である。
ブラインド処理部１６０は、式［数５８］により受信信号ベクトルＸを生成して、図２に示される実施形態１の構成例の場合と同様に、最大固有地に対応した固有ベクトルを算出し［数６１］として重み付け係数ｗ１及びｗ２を白色化・タップ処理部２１１と白色化・タップ処理部２１２へ送る。
タップ処理部２１１はブラインド処理部１６０から送られてきたｗ１を基に
なる演算を行う。タップ処理部２１２はブラインド処理部１６０から送られてきたｗ２を基に
なる演算を行う。
この乗算されたタップ出力は、ダイバーシチ合成部１７０へ送られブランチ１と２の信号が合成されダイバーシチ処理が完結する。即ち上述の［数６２］式に相当する操作が行われているのである。具体的に書くと、
となる。
本実施形態では、固有ベクトルβを重み付け係数として用いているので、ダイバーシチ合成後の信号対雑音電力比（ＳＮＲ）γが最大となり良好なダイバーシチ利得が得られる。
更に、本実施形態では、白色化・タップ処理部２１１と白色化・タップ処理部２１２による白色化処理があるので、大地反射波等の遅延干渉波が存在する伝播環境での急激な特性劣化を緩和する効果が期待できる。
上述した本実施形態の効果をシミュレーションにより確認した結果が図５に示されている。図５で、縦軸はシンボル誤り率、横軸は遅延干渉波の遅延時間である。遅延時間が大きくなる程干渉波の影響が厳しい。
図５中の「ｆ：本発明の固有ベクトル合成（白色化有り）」が本実施形態を用いた場合のシミュレーション結果に相当する。
なお、図５には本実施形態との比較の為、一般的な技術で説明したＩＰ合成法、等化機能を持たせたもの、ダイバーシチ機能が無いもの等の特性も載せている。
具体的に検討すると、一般的な技術として紹介した比較対象「ｄ：ＩＰ合成法」の場合、干渉波の遅延が少し増えただけで特性が極端に劣化していることが分かる。
また上述の「ｅ：固有ベクトル合成（白色化処理無し）」の場合も特性が極端に劣化している。これに対して本実施形態の「ｆ：本発明の固有ベクトル合成（白色化有り）」の場合、等化機能の有る比較対象
「ｂ：ダイバーシチ無しでＬＭＳ等化機能有り」や「ｃ：ダイバーシチ有りで、等化機能有り」には及ばないが、白色化処理により特性劣化を緩和していることが分かる。尚この例はかなり厳しい干渉波を想定した場合であり、後述する様に別の伝搬環境では良好なダイバーシチ効果を得ている。
［実施形態３］
次に、厳しい干渉波が存在する伝搬環境で急激に特性が劣化することについての対策として、実施形態２とは別の白色化処理による特性劣化を緩和する実施形態である実施形態３について説明する。本実施形態では受信信号から干渉波信号を検出する手段を用いる。
すなわち、実施形態３は、白色化処理により特性劣化を緩和して最大固有値の固有ベクトルにより合成した方式（干渉波信号検出手段を用いた場合）に関するものである。
まず、本実施形態での白色化処理の原理について説明する。
上述の式［数１］でＳＮＲγは
となっているが、遅延干渉波が存在する伝播環境では以下のＳＩＮＲ（Signal−to−Interference plus Noise power Ratio）γとして考える必要がある。
ここで、ＲＩは干渉波信号ベクトルＸＩの共分散行列であり以下の式で与えられる。
そこで、
なる修正行列ＡＩにより、
とすると、［数６７］の分母は、
となる。
即ち、干渉波信号ベクトルＸＩの共分散行列ＲＩを単位行列になる様に修正行列ＡＩは作用する。干渉波信号ベクトルＸＩをこの修正行列ＡＩで変換すると、ベクトル
の要素間は無相関となる。尚、干渉波信号ベクトルＸＩは雑音成分も含んでいる。
従って、この修正行列ＡＩで変換した後の［数６７］は、
ここで、
とすると、
となり、
ＳＩＮＲγを最大化するｗは、修正行列ＡＩを用いて共分散行列Ｒを修正した修正共分散行列
の最大固有値に対応した固有ベクトルとなる。例えば、修正行列ＡＩによって変換された受信信号ベクトルＸは、元の受信信号ベクトルをＵとすると、
となる。ここでｕ１はブランチ１の受信信号、ｕ２はブランチ２の受信信号である。従って、
の関係がある。ウェイトは［数７４］の最大固有値に対応した固有ベクトルで
とすると、ダイバーシチ合成後の信号ｙは
となる。
上述した動作を実行するための動作を行っている構成例を図６に示す。
図６を参照すると、本実施形態の受信装置３０００は、アンテナ１１１、アンテナ１１２、ミキサ１２１、ミキサ１２２、Ａ／Ｄ変換器１３１、Ａ／Ｄ変換器１３２、局部発振器１５０、干渉波信号検出部３１０及び白色化ブラインド処理部３２０を有する。
図２に示される実施形態１の構成例と共通する符号が振られている部分での処理は上述と同じなので説明を省略する。
一方、本実施形態特有の箇所である干渉波信号検出部３１０と白色化ブラインド処理部３２０については、以下説明する。
干渉波検出部３１０は各ブランチの受信信号より干渉波を検出する。
干渉波検出部３１０によって検出された干渉波は信号ベクトルの形に生成して白色化ブラインド処理部３２０へ送られる。この送られる干渉波信号ベクトルが［数６８］のＸＩに相当する。
白色化ブラインド処理部３２０では入力された干渉波信号ベクトルＸＩを基にして［数６８］の演算を行う。実際には［数４２］の様な演算を用いても良い。これによって干渉波信号ベクトルから共分散行列ＲＩを得る。
更に、白色化ブラインド処理部３２０では［数６９］の関係より修正行列ＡＩを得る。この修正行列ＡＩは［数７０］の関係より共分散行列ＲＩを単位行列にもっていく様に作用することが分かる。白色化ブラインド処理部３２０には、ブランチ１の受信信号ｕ１をＡ／Ｄ変換器１３１から、ブランチ２の受信信号ｕ２をＡ／Ｄ変換器１３２から得ている。
そして、この受信信号ｕ１及び受信信号ｕ２を受信信号ベクトルＵとして統合し、例えば［数４２］の様な演算を用いて共分散行列Ｒを生成する。上の修正行列ＡＩを用いてこの共分散行列Ｒを［数７４］により修正する。これを修正共分散行列と呼ぶことにする。
加えて、白色化ブラインド処理部３２０は、修正共分散行列の最大固有値に対応した固有ベクトルを上述と同じ手順により得ると同時に、修正行列ＡＩを用いて受信信号ベクトルＵを［数７５］により変換し、変換された受信信号ベクトルＸを得る。
上記固有ベクトルは式［数７７］により重み付け係数として使われるウェイトベクトルＷとなる。変換された受信信号ベクトルＸとウェイトベクトルＷにより式［数７８］の演算が行われ、ダイバーシチ合成後の信号ｙを得て出力してダイバーシチ処理が完結する。
以上の結果、［数７３］のＳＩＮＲを最大化する修正共分散行列の最大固有ベクトルβを重み付け係数として用いているので、ダイバーシチ合成後の信号対干渉波雑音電力比（ＳＩＮＲ）γが最大となり、良好なダイバーシチ利得が期待出来る。即ち干渉波信号検出部３１０と白色化ブラインド処理部３２０によって、［数６７］の分母を白色化してその分子に反映させているので、大地反射波等の遅延干渉波が存在する劣悪な伝播環境でも急激な特性劣化を起こすことがない。
上述した本実施形態の効果をシミュレーションにより確認した結果が図７及び８に示されている。
図７及び８中の「ｇ：本発明の固有ベクトル合成（干渉波信号検出手段を用いた白色化）」が実施形態３を用いた場合のシミュレーション結果に相当する。
加えて、図７及び８では同時に一つ前の構成である実施形態２を用いた場合のシミュレーション結果を「ｆ：本発明の固有ベクトル合成（白色化有り）」として示す。
共に同じ二波モデルによる遅延干渉波の伝搬環境でのシミュレーション結果である。
なお、実施形態１の白色化を用いない固有ベクトル合成では６４ＱＡＭのコンスタレーションが干渉波に埋もれてしまって特性にならないので図７及び８では図示を省略している。即ち、白色化を用いた実施形態２及び実施形態３は、白色化処理により特性劣化を緩和していることが分かる。
同シミュレーションで使われている二波モデルによる遅延干渉波の様子を詳しく見ることにする。
図９ａは図７のシミュレーションで用いた伝搬環境をベクトル表記したものである。同図より直接波成分はブランチ１と２とも同じ方向を示している。一方、直接波成分とは逆に干渉波成分はブランチ１と２で逆相である。但し、ブランチ１の干渉波成分の大きさはブランチ２の０．２５倍となっている。
従って干渉波を相殺しようとした場合、第二のタップ要素ｗ２の大きさを第一のタップ要素ｗ１の０．２５とすることになる。しかし、このようにすると、ブランチ２の直接波の大きさが減少してダイバーシチ利得が減少する。
図９ａの下の表は、その時の比を、シミュレーション値に基づいてＳＮＲ毎に観測した値である。同表より、ＳＮＲが高いとき、干渉波成分を相殺するように０．２５となっているが、ＳＮＲが低くなるにつれて０．２５より高い値となって、ダイバーシチ利得を得る方向へと作用している。
即ち、白色化処理はただ干渉波成分を相殺する様に作用するだけでなく、ダイバーシチ利得を考慮する様に作用していることが分かる。
図９ｂは図８のシミュレーションで用いた伝搬環境をベクトル表記したものである。図９ｂより干渉成分を相殺するには第一のタップ要素の位相と第二のタップ要素の位相の差を大きくすることになる。
シミュレーション結果はＳＮＲが高いうちは位相差をつけて逆相となるが、ＳＮＲが低い時はその差が小さくなっていた。しかし、それでも実施形態３の場合、各ブランチの干渉波を相殺する様に作用する。そのため、この例で使われた伝搬環境の場合、実施形態２より特性が劣っているが、より厳しい伝搬環境では逆転して実施形態３の方が特性が良くなる。
［算出の簡略化］
以上では、直接波のみで大地反射波等の干渉波が存在しない伝搬環境での実施形態１の効果や、特性が急激に劣化する大地反射波等の遅延干渉波が存在する伝播環境で白色化処理により特性劣化を緩和する実施形態２や３の効果について説明してきた。ここからは、各構成で用いられてきた最大固有値の固有ベクトルを簡単に算出する方法を説明する。以下で述べる算出方法は、上述した各実施形態で利用することが可能である。
［二つのアンテナによりスペースダイバーシチを実現する際に用いる代数的演算手法］
二つのアンテナによりスペースダイバーシチを実現する際に用いる共分散行列の最大固有値に対応した固有ベクトルを、共分散行列から算出する代数的演算手法について説明する。
共分散行列Ｒはエルミート行列だから
とする。この固有値λは、
によって計算出来、
を得る。この内大きい方の固有値λを基に固有ベクトルβを
但し、ｘは定数であるものとする。
として計算出来る。
何とならば、固有ベクトルβの定義により
で、
となるからである。
［二つのアンテナによりスペースダイバーシチを実現する際に用いる更に簡略化した手法］
式［数７９］の共分散行列は更に簡略化出来、受信データが定常なら、
になる。従って、式［数８０］の固有値は、
となる。大きい方の固有値を選択するから、当然、
従って、式［数８１］の固有ベクトルβは、
但し、ｘは定数であるものとする。
として計算出来る。
即ち、最大固有値に対応する固有ベクトルを、最大固有値を求めることなく、共分散行列より直接得ることが出来、非常に簡単な処理で実現することが可能となる。
以上の説明では受信装置のアンテナが２つの場合を想定していた。次に、受信装置のアンテナが三つの場合で三次元の行列に対する最大固有値の固有ベクトルを簡単に算出する方法を説明する。
［三つのアンテナによりスペースダイバーシチを実現する際に用いる代数的演算手法］
三つのアンテナによりスペースダイバーシチを実現する際に用いる共分散行列の最大固有値に対応した固有ベクトルを、共分散行列から算出する代数的演算手法について説明する。
共分散行列Ｒはエルミート行列だから
とする。この固有値λは、
によって計算出来る。上式を展開すると、
ここで、
とすると、
そこで、
とすると、上式は、
整理していくと、
ここで、式［数９１］より、
更に、
であるから、
一方、
従って、［数９８］は、
となる。そこで、
と置くと、式［数９５］の固有値の条件は、
を得る。ここで、
とすると、ｆ（ｙ）は原点を通る三次特性となる。極値は微分して、
より
の点に有り、その時の値は、
従って、この範囲内
にｆ（ｙ）があれば、ｙの取り得る三個の実数の値が存在し、それが［数１０１］、［数９３］の実数による変換を経て固有値となる。
エルミート行列の固有値は必ず実数となるから、根ｙは実数である必要がある。この範囲外になるとｙの取る根は虚数を含むことになるから、エルミート行列の固有値となりえない。［数１０２］の関係より、ｆ（ｙ）の取る値の絶対値は１以下でなければこの条件は満たされない。従って、端点として
より、左側端点はｙ＝−１、右側端点はｙ＝１となる。
の範囲にあれば全ての固有値が実数となりエルミート行列の固有値となる。これを
なるｚで表す。すると、
従って、
を得る。これより、
ここで、ｎは任意の整数。更に、
とすると、
従って、
となる。固有値λは、
とすると、式［数１０１］、［数９３］より、
であるから、結局三つの固有値λ１、λ２、λ３は、
となる。
この内大きい方の固有値λを基に固有ベクトルβを
但し、ｘは定数であるものとする。
として計算出来る。
何とならば、固有ベクトルβの定義により
で、
となるからである。ここで、第二要素は受信データの定常性から、
とした。更に第三要素は、式［数９０］より、
なる関係を用いた。
更に、上記固有ベクトルが零になる時も考慮しなければならない。その時は、
の第１要素より、
となる。 ・・・（式ａ）
更に第２要素から、この関係を用いて、
・・・（式ｂ）
更に第３要素から、
・・・（式ｃ）
となる。（これは（式ｂ）と同じ。また固有値の条件から下記の様に（式ｂ）、（式ｃ）を得ることが出来る。）
即ち（式ａ）を用いて、式［数８９］を式［数１２０］の基に以下の様に書き直す。
そこで、以下の式を固有ベクトルβとして用いる。
（
の場合）
但し、ｘは定数であるものとする。
何とならば、固有ベクトルβの定義により
で、
となるからである。ここで、第一要素は上記の条件から、
なる関係を用いた。更に、受信データの定常性から、
である。
更に、第２要素は、（式ｂ）の
を用いた。
更に、第３要素は、［数１２７］を用いると、
となり、ｙに対する係数は（式ｂ）より
ｘに対する係数は（式ｂ）より、
従って、
を用いた。
このように３次の共分散行列に対しても各構成で用いられてきた最大固有値の固有ベクトルを代数的演算によって簡単に算出することが出来る。
以上説明した本発明の実施形態は、以下に示すような多くの効果を奏する。
第１の効果は実際の情報の伝達に何ら寄与しないパイロット信号やリファレンス信号或いはトレーニング信号など伝送効率を低下させる既知の系列を用いることのないブラインド処理で、高速で高次の多値変調を用いる大容量マイクロ波無線通信に厳しい遅延干渉波が存在する周波数選択性フェージング環境下でも、良好なダイバーシチ利得を得ることが出来ることである。
その理由は、白色化処理を用いるからである。
第２の効果はマイクロ波無線受信装置に用いられている固有値ベクトルを得る為に、固有値分解等の複雑な方法ではなく、簡単な積和演算で実現出来、複雑度を削減出来ることである。また、それに伴い回路規模を小さく出来ることである。
その理由は、本発明の実施形態では代数的演算手法を用いるからである。
なお、上述した実施形態は、本発明の好適な実施形態ではあるが、上記実施形態のみに本発明の範囲を限定するものではなく、本発明の要旨を逸脱しない範囲において種々の変更を施した形態での実施が可能である。
また、本発明の実施形態であるマイクロ波無線受信装置は、ハードウェアにより実現することもできるが、コンピュータをそのマイクロ波無線受信装置として機能させるためのプログラムをコンピュータがコンピュータ読み取り可能な記録媒体から読み込んで実行することによっても実現することができる。
また、本発明の実施形態によるマイクロ波無線受信方法は、ハードウェアにより実現することもできるが、コンピュータにその方法を実行させるためのプログラムをコンピュータがコンピュータ読み取り可能な記録媒体から読み込んで実行することによっても実現することができる。
この点、マイクロ波無線受信装置を、装置の全部又は一部として動作させ、あるいは、上述の処理を実行させるためのプログラムは、様々なタイプの非一時的なコンピュータ可読媒体(non-transitory computer readable medium)を用いて格納され、コンピュータに供給することができる。非一時的なコンピュータ可読媒体は、様々なタイプの実体のある記録媒体(tangible storage medium)を含む。非一時的なコンピュータ可読媒体の例は、磁気記録媒体（例えば、フレキシブルディスク、磁気テープ、ハードディスクドライブ）、光磁気記録媒体（例えば、光磁気ディスク）、ＣＤ−ＲＯＭ(Read Only Memory)、ＣＤ−Ｒ、ＣＤ−Ｒ／Ｗ、半導体メモリ（例えば、マスクＲＯＭ、ＰＲＯＭ(Programmable ROM)、ＥＰＲＯＭ(Erasable PROM)、フラッシュＲＯＭ、ＲＡＭ(random access memory)）を含む。また、プログラムは、様々なタイプの一時的なコンピュータ可読媒体(transitory computer readable medium)によってコンピュータに供給されてもよい。一時的なコンピュータ可読媒体の例は、電気信号、光信号、及び電磁波を含む。一時的なコンピュータ可読媒体は、電線及び光ファイバ等の有線通信路、又は無線通信路を介して、プログラムをコンピュータに供給できる。
上記の実施形態の一部又は全部は、以下の付記のようにも記載されうるが、以下には限られない。
（付記１） 複数のアンテナで受信した受信信号毎に対応するブランチを設け、当該各ブランチの受信信号を合成することによりスペースダイバーシチを実現する無線受信装置において、
前記各ブランチの受信信号を重み付け係数に基づくウェイト処理によりダイバーシチ合成する際に、
前記ウェイト処理対象とする受信信号の受信信号ベクトルから得られた共分散行列が単位行列となる様な修正行列を生成し、当該修正行列によって前記受信信号を変換し、当該変換後の受信信号の受信信号ベクトルから得られた共分散行列の最大固有値に対応した固有ベクトルを算出し、当該固有ベクトルを前記重み付け係数として前記ダイバーシチ合成を行うことを特徴とする無線受信装置。
（付記２） 複数のアンテナで受信した受信信号毎に対応するブランチを設け、当該各ブランチの受信信号を合成することによりスペースダイバーシチを実現する無線受信装置において、
前記各ブランチの受信信号から干渉波信号を検出する干渉波信号検出部と、
前記干渉波信号から得られた共分散行列が単位行列となる様な修正行列を生成し、前記受信信号の受信信号ベクトルから得られた共分散行列を前記修正行列により修正し、当該修正後の共分散行列の最大固有値に対応した固有ベクトルを用いて重み付け係数を生成し、当該該重み付け係数に基づくウェイト処理によりダイバーシチ合成を行うブラインド処理部と、
を備えることを特徴とする無線受信装置。
（付記３） 付記１又は２に記載の無線受信装置において、
前記共分散行列の最大固有値に対応した固有ベクトルの算出は、当該共分散行列の冪乗演算によって最大固有値に対する固有ベクトルを算出することにより行われることを特徴とする無線受信装置。
（付記４） 付記１又は２に記載の無線受信装置において、
前記共分散行列の最大固有値に対応した固有ベクトルの算出は、当該共分散行列から代数的演算によって最大固有値を算出することにより行われることを特徴とする無線受信装置。
（付記５） 付記４に記載の無線受信装置において、
前記アンテナの数は２つであり、
前記共分散行列の最大固有値に対応した固有ベクトルを該共分散行列から算出する際、
共分散行列Ｒを
とした時、固有値λを
の如く計算し、この内大きい方の固有値λを基に固有ベクトルβを
但し、ｘは定数であるものとする、
の如く計算することを特徴とする無線受信装置。
（付記６） 付記４に記載の無線受信装置において、
前記アンテナの数は３つであり、
前記共分散行列の最大固有値に対応した固有ベクトルを該共分散行列から算出する際、
共分散行列Ｒを
とした時、
なる変数の元に、あるいは直接的に、
固有値λを
の如く計算し、この内一番大きい固有値λを基に固有ベクトルβを
（
の場合）
但し、ｘは定数であるものとする、
或いは、
（
の場合）
但し、ｘ、ｙは定数であるものとする、
の如く計算することを特徴とする無線受信装置。
（付記７） 付記４に記載の無線受信装置において、
前記アンテナの数は２つであり、
前記共分散行列の最大固有値に対応した固有ベクトルを該共分散行列から算出する際、
共分散行列Ｒを
とした時、固有ベクトルβを
但し、ｘは定数であるものとする、
の如く計算することを特徴とする無線受信装置。
（付記８） 信号を無線送信する無線送信装置と、前記送信された信号を複数のアンテナにより受信する無線受信装置と、を含んだ無線通信システムにおいて、
前記無線送信装置と前記無線受信装置間の無線通信はブラインド処理を実現した通信方式に準拠しており、
前記無線受信装置は付記１乃至７の何れか１に記載の無線受信装置である、
ことを特徴とする無線通信システム。
（付記９） 複数のアンテナで受信した受信信号毎に対応するブランチを設け、当該各ブランチの受信信号を合成することによりスペースダイバーシチを実現する無線受信装置が行う無線受信方法において、
前記各ブランチの受信信号を重み付け係数に基づくウェイト処理によりダイバーシチ合成する際に、
前記ウェイト処理対象とする受信信号の受信信号ベクトルから得られた共分散行列が単位行列となる様な修正行列を生成し、当該修正行列によって前記受信信号を変換し、当該変換後の受信信号の受信信号ベクトルから得られた共分散行列の最大固有値に対応した固有ベクトルを算出し、当該固有ベクトルを前記重み付け係数として前記ダイバーシチ合成を行うことを特徴とする無線受信方法。
（付記１０） 複数のアンテナで受信した受信信号毎に対応するブランチを設け、当該各ブランチの受信信号を合成することによりスペースダイバーシチを実現する無線受信装置が行う無線受信方法において、
前記各ブランチの受信信号から干渉波信号を検出する干渉波信号検出し、
前記干渉波信号から得られた共分散行列が単位行列となる様な修正行列を生成し、前記受信信号の受信信号ベクトルから得られた共分散行列を前記修正行列により修正し、当該修正後の共分散行列の最大固有値に対応した固有ベクトルを用いて重み付け係数を生成し、当該該重み付け係数に基づくウェイト処理によりダイバーシチ合成を行う、ことを特徴とする無線受信方法。
（付記１１） 複数のアンテナで受信した受信信号毎に対応するブランチを設け、当該各ブランチの受信信号を合成することによりスペースダイバーシチを実現する無線受信装置が行う無線受信方法において、
前記各ブランチの受信信号から干渉波信号を検出する干渉波信号検出ステップと、
前記干渉波信号から得られた共分散行列が単位行列となる様な修正行列を生成し、前記受信信号の受信信号ベクトルから得られた共分散行列を前記修正行列により修正し、当該修正後の共分散行列の最大固有値に対応した固有ベクトルを用いて重み付け係数を生成し、当該該重み付け係数に基づくウェイト処理によりダイバーシチ合成を行うブラインド処理ステップと、
を備えることを特徴とする無線受信方法。
（付記１２） 付記１０又は１１に記載の無線受信方法において、
前記共分散行列の最大固有値に対応した固有ベクトルの算出は、当該共分散行列の冪乗演算によって最大固有値に対する固有ベクトルを算出することにより行われることを特徴とする無線受信方法。
（付記１３） 付記１０又は１１に記載の無線受信方法において、
前記共分散行列の最大固有値に対応した固有ベクトルの算出は、当該共分散行列から代数的演算によって最大固有値を算出することにより行われることを特徴とする無線受信方法。
（付記１４） 付記１３に記載の無線受信方法において、
前記アンテナの数は２つであり、
前記共分散行列の最大固有値に対応した固有ベクトルを該共分散行列から算出する際、
共分散行列Ｒを
とした時、固有値λを
の如く計算し、この内大きい方の固有値λを基に固有ベクトルβを
但し、ｘは定数であるものとする、
の如く計算することを特徴とする無線受信方法。
（付記１５） 付記１３に記載の無線受信方法において、
前記アンテナの数は３つであり、
前記共分散行列の最大固有値に対応した固有ベクトルを該共分散行列から算出する際、
共分散行列Ｒを
とした時、
なる変数の元に、あるいは直接的に、
固有値λを
の如く計算し、この内一番大きい固有値λを基に固有ベクトルβを
（
の場合）
但し、ｘは定数であるものとする、
或いは、
（
の場合）
但し、ｘ、ｙは定数であるものとする、
の如く計算することを特徴とする無線受信方法。
（付記１６） 付記１３に記載の無線受信方法において、
前記アンテナの数は２つであり、
前記共分散行列の最大固有値に対応した固有ベクトルを該共分散行列から算出する際、
共分散行列Ｒを
とした時、固有ベクトルβを
但し、ｘは定数であるものとする、
の如く計算することを特徴とする無線受信方法。
本願は、日本の特願２０１１−１６６７５０（２０１１年７月２９日に出願）に基づいたものであり、又、特願２０１１−１６６７５０に基づくパリ条約の優先権を主張するものである。特願２０１１−１６６７５０の開示内容は、特願２０１１−１６６７５０を参照することにより本明細書に援用される。
本発明の代表的な実施の形態が詳細に述べられたが、様々な変更(changes)、置き換え(substitutions)及び選択(alternatives)が請求項で定義された発明の精神と範囲から逸脱することなくなされることが理解されるべきである。また、仮にクレームが出願手続きにおいて補正されたとしても、クレームされた発明の均等の範囲は維持されるものと発明者は意図する。

Claims (10)

  1. 複数のアンテナで受信した受信信号毎に対応するブランチを設け、当該各ブランチの受信信号を合成することによりスペースダイバーシチを実現する無線受信装置において、
    前記各ブランチの受信信号を重み付け係数に基づくウェイト処理によりダイバーシチ合成する際に、
    前記ウェイト処理対象とする受信信号の受信信号ベクトルから得られた共分散行列が単位行列となる様な修正行列を生成し、当該修正行列によって前記受信信号を変換し、当該変換後の受信信号の受信信号ベクトルから得られた共分散行列の最大固有値に対応した固有ベクトルを算出し、当該固有ベクトルを前記重み付け係数として前記ダイバーシチ合成を行うことを特徴とする無線受信装置。
  2. 複数のアンテナで受信した受信信号毎に対応するブランチを設け、当該各ブランチの受信信号を合成することによりスペースダイバーシチを実現する無線受信装置において、
    前記各ブランチの受信信号から干渉波信号を検出する干渉波信号検出部と、
    前記干渉波信号から得られた共分散行列が単位行列となる様な修正行列を生成し、前記受信信号の受信信号ベクトルから得られた共分散行列を前記修正行列により修正し、当該修正後の共分散行列の最大固有値に対応した固有ベクトルを用いて重み付け係数を生成し、当該該重み付け係数に基づくウェイト処理によりダイバーシチ合成を行うブラインド処理部と、
    を備えることを特徴とする無線受信装置。
  3. 請求項１又は２に記載の無線受信装置において、
    前記共分散行列の最大固有値に対応した固有ベクトルの算出は、当該共分散行列の冪乗演算によって最大固有値に対する固有ベクトルを算出することにより行われることを特徴とする無線受信装置。
  4. 請求項１又は２に記載の無線受信装置において、
    前記共分散行列の最大固有値に対応した固有ベクトルの算出は、当該共分散行列から代数的演算によって最大固有値を算出することにより行われることを特徴とする無線受信装置。
  5. 請求項４に記載の無線受信装置において、
    前記アンテナの数は２つであり、
    前記共分散行列の最大固有値に対応した固有ベクトルを該共分散行列から算出する際、
    共分散行列Ｒを
    とした時、固有値λを
    の如く計算し、この内大きい方の固有値λを基に固有ベクトルβを
    但し、ｘは定数であるものとする、
    の如く計算することを特徴とする無線受信装置。
  6. 請求項４に記載の無線受信装置において、
    前記アンテナの数は３つであり、
    前記共分散行列の最大固有値に対応した固有ベクトルを該共分散行列から算出する際、
    共分散行列Ｒを
    とした時、
    なる変数の元に、あるいは直接的に、
    固有値λを
    の如く計算し、この内一番大きい固有値λを基に固有ベクトルβを
    （
    の場合）
    但し、ｘは定数であるものとする、
    或いは、
    （
    の場合）
    但し、ｘ、ｙは定数であるものとする、
    の如く計算することを特徴とする無線受信装置。
  7. 請求項４に記載の無線受信装置において、
    前記アンテナの数は２つであり、
    前記共分散行列の最大固有値に対応した固有ベクトルを該共分散行列から算出する際、
    共分散行列Ｒを
    とした時、固有ベクトルβを
    但し、ｘは定数であるものとする、
    の如く計算することを特徴とする無線受信装置。
  8. 信号を無線送信する無線送信装置と、前記送信された信号を複数のアンテナにより受信する無線受信装置と、を含んだ無線通信システムにおいて、
    前記無線送信装置と前記無線受信装置間の無線通信はブラインド処理を実現した通信方式に準拠しており、
    前記無線受信装置は請求項１乃至７の何れか１項に記載の無線受信装置である、
    ことを特徴とする無線通信システム。
  9. 複数のアンテナで受信した受信信号毎に対応するブランチを設け、当該各ブランチの受信信号を合成することによりスペースダイバーシチを実現する無線受信装置が行う無線受信方法において、
    前記各ブランチの受信信号を重み付け係数に基づくウェイト処理によりダイバーシチ合成する際に、
    前記ウェイト処理対象とする受信信号の受信信号ベクトルから得られた共分散行列が単位行列となる様な修正行列を生成し、当該修正行列によって前記受信信号を変換し、当該変換後の受信信号の受信信号ベクトルから得られた共分散行列の最大固有値に対応した固有ベクトルを算出し、当該固有ベクトルを前記重み付け係数として前記ダイバーシチ合成を行うことを特徴とする無線受信方法。
  10. 複数のアンテナで受信した受信信号毎に対応するブランチを設け、当該各ブランチの受信信号を合成することによりスペースダイバーシチを実現する無線受信装置が行う無線受信方法において、
    前記各ブランチの受信信号から干渉波信号を検出し、
    前記干渉波信号から得られた共分散行列が単位行列となる様な修正行列を生成し、前記受信信号の受信信号ベクトルから得られた共分散行列を前記修正行列により修正し、当該修正後の共分散行列の最大固有値に対応した固有ベクトルを用いて重み付け係数を生成し、当該該重み付け係数に基づくウェイト処理によりダイバーシチ合成を行う、ことを特徴とする無線受信方法。
JP2013526893A 2011-07-29 2012-07-27 無線受信装置及び無線受信方法並びに無線システム Active JP5867968B2 (ja)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2013526893A JP5867968B2 (ja) 2011-07-29 2012-07-27 無線受信装置及び無線受信方法並びに無線システム

Applications Claiming Priority (4)

Application Number Priority Date Filing Date Title
JP2011166750 2011-07-29
JP2011166750 2011-07-29
JP2013526893A JP5867968B2 (ja) 2011-07-29 2012-07-27 無線受信装置及び無線受信方法並びに無線システム
PCT/JP2012/069198 WO2013018716A1 (ja) 2011-07-29 2012-07-27 無線受信装置及び無線受信方法並びに無線システム

Publications (2)

Publication Number Publication Date
JPWO2013018716A1 JPWO2013018716A1 (ja) 2015-03-05
JP5867968B2 true JP5867968B2 (ja) 2016-02-24

Family

ID=47629239

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2013526893A Active JP5867968B2 (ja) 2011-07-29 2012-07-27 無線受信装置及び無線受信方法並びに無線システム

Country Status (2)

Country Link
JP (1) JP5867968B2 (ja)
WO (1) WO2013018716A1 (ja)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP7486305B2 (ja) * 2019-11-22 2024-05-17 日本無線株式会社 無線受信装置
JP7289600B2 (ja) * 2019-11-27 2023-06-12 日本無線株式会社 無線受信装置
CN112615652B (zh) * 2020-12-17 2022-06-10 吉首大学 一种适用于高阶mimo的基于特征值检验的半盲频谱感知方法
CN116827367A (zh) * 2023-07-06 2023-09-29 厦门纵行信息科技有限公司 一种强干扰源下的干扰消除接收机及接收方法

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP4403916B2 (ja) * 2004-08-10 2010-01-27 株式会社豊田中央研究所 重み係数演算装置
JP2006270847A (ja) * 2005-03-25 2006-10-05 Mitsubishi Electric Corp アンテナ装置
JP2008312188A (ja) * 2007-05-14 2008-12-25 Mitsubishi Electric Corp アダプティブアンテナ
JP5183376B2 (ja) * 2008-09-04 2013-04-17 三菱電機株式会社 無線受信装置および復調方法

Also Published As

Publication number Publication date
JPWO2013018716A1 (ja) 2015-03-05
WO2013018716A1 (ja) 2013-02-07

Similar Documents

Publication Publication Date Title
US8638260B2 (en) Transmitter beamforming steering matrix processing and storage
US7830988B2 (en) MIMO decoder and MIMO decoding method
WO2010025203A1 (en) System and method for wireless communications
JP5867968B2 (ja) 無線受信装置及び無線受信方法並びに無線システム
CN110870216B (zh) 一种波束成形方法及设备
Park et al. Regularized interference alignment based on weighted sum-MSE criterion for MIMO interference channels
JP2017204855A (ja) 受信機回路
US20220209830A1 (en) Wireless communication system, wireless communication method, transmitting station device and receiving station device
CN114830607B (zh) 用于大型多维无线系统的无线x2x接入方法和接收器
WO2016148262A1 (ja) 通信装置、方法及びシステムと端末とプログラム
JP6166641B2 (ja) 重み係数算出装置、アダプティブアレー受信装置、及び、重み係数算出方法
CN110138413A (zh) 一种空时域自适应宽线性降秩波束形成方法
US9853836B2 (en) Apparatus and method for signal detection in a wireless communication system
JP5340196B2 (ja) 送信ウェイト決定方法
WO2019125240A1 (en) Steering vector weighting
KR101556482B1 (ko) 다중안테나 간섭채널에서 slnr 최대화를 통한 간섭정렬 방법 및 장치
KR20180080768A (ko) 복수의 무선 신호들을 컴바이닝하는 방법 및 장치
KR102425095B1 (ko) 주파수 분할-다중 입력 다중 출력 방식을 지원하는 통신 시스템에서 빔 포밍 동작을 수행하는 장치 및 방법
JP6254962B2 (ja) 基地局装置、無線通信方法及び無線通信システム
JP5812010B2 (ja) 無線通信装置
CN114301567B (zh) 一种基于人工噪声的通信方法及设备
US11870517B2 (en) Wireless communication system, wireless communication method, transmitting station device and receiving station device
JP4901643B2 (ja) アダプティブアレーアンテナ装置およびプログラム
Anit et al. Hybrid Analog and Digital Beamforming for SU-MIMO Frameworks
WO2017002323A1 (ja) 無線通信システムと基地局装置とプリコーダ決定方法と記録媒体

Legal Events

Date Code Title Description
A02 Decision of refusal

Free format text: JAPANESE INTERMEDIATE CODE: A02

Effective date: 20141113

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20160104

R150 Certificate of patent or registration of utility model

Ref document number: 5867968

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150