JP2013236127A

JP2013236127A - 無線基地局

Info

Publication number: JP2013236127A
Application number: JP2012105456A
Authority: JP
Inventors: Kazuhiko Fukawa; 和彦府川; Satoshi Suyama; 聡須山; Hiroshi Suzuki; 博鈴木; Yusuke Owatari; 裕介大渡; Nobuhiko Miki; 信彦三木; Yosuke Sano; 洋介佐野; Yukihiko Okumura; 幸彦奥村
Original assignee: NTT Docomo Inc; Tokyo Institute of Technology NUC
Current assignee: NTT Docomo Inc; Tokyo Institute of Technology NUC
Priority date: 2012-05-02
Filing date: 2012-05-02
Publication date: 2013-11-21
Anticipated expiration: 2032-05-02
Also published as: JP6029848B2

Abstract

【課題】同一チャネル干渉を抑圧し、ユーザ端末で受信される信号の受信品質の劣化を抑えられる、伝送路情報のフィードバック不要のプリコーディング方式を実行する無線基地局を提供する。
【解決手段】無線基地局は、上りリンクの受信信号のベクトルの積を平均化して上りリンクの自己相関行列を求め、その自己相関行列から下りリンクの自己相関行列を推定する。無線基地局は、推定された下りリンクの自己相関行列のうち、所望ユーザ端末の下りリンクの自己相関行列を固有値分解して固有ベクトルを求め、下りリンク自己相関行列生成部で推定された下りリンクの自己相関行列のうち、全ユーザ端末または他ユーザ端末の下りリンクの自己相関行列に基づいて、他ユーザ端末への送信信号に起因する同一チャネル干渉を抑圧するための干渉抑圧行列を生成し、固有ベクトルと干渉抑圧行列に基づいて、プリコーディング行列を生成するプリコーディング行列変換部を備える。
【選択図】図７

Description

本発明は、無線基地局に関するものである。

携帯電話システム等の無線通信において、周波数帯域を広げずに伝送速度を高める技術として、複数の送信アンテナと複数の受信アンテナを用いて空間多重伝送を行うMIMO(Multiple Input Multiple Output)伝送が注目を集めている。このMIMO伝送の更なる特性改善を図る送信技術として、送信側でチャネルインパルス応答等の伝送路情報が既知の場合、伝送路に応じて送信信号に線形処理を施すプリコーディング技術が知られている。プリコーディングは、送信信号の指向性を高めてメインビーム（メインローブ）を送信先に向けることができる。

上りリンクと下りリンクで異なる周波数帯域を使用するFDD(Frequency Division Duplex)の場合、送信側で下りリンクの伝送路情報が直接得られないので、プリコーディングを行うためには、受信側から送信側へ下りリンクの伝送路情報をフィードバックする必要がある。

以下では、複数のユーザ端末が通信チャネルを共有するマルチユーザMIMOについて、伝送路情報のフィードバックを用いるMIMO送信機の説明を行う。送信側を基地局、受信側を複数のユーザ端末とし、基地局からユーザ端末への下りリンクにおいてプリコーディングを用いた通信を行い、ユーザ端末から基地局への上りリンクにおいて下りリンク伝送路情報のフィードバックを行うものとする。

図１に、従来のMIMOシングルキャリア伝送の無線送受信機（無線基地局）で、ユーザ端末からフィードバックされた下りリンクの伝送路情報を用いてプリコーディングを行う構成を示す。説明を簡単にするため、下りリンクのフェージングはフラットフェージングと想定する。無線送受信機の送信アンテナ（アダプティブ・アレー・アンテナ）１０の数はN_T、ユーザ端末の数はKと想定し、第kユーザ端末(1≦k≦K)の受信アンテナの数はN_rkと想定する。以下、無線送受信機の受信アンテナの数は、送信アンテナ１０の数N_Tと等しいと想定する。無線送受信機に送信アンテナ１０と受信アンテナを兼ねる複数の送受信アンテナを設けてもよい。

まず、時刻iにおいて、第kユーザ端末へのN_rk次元の変調信号ベクトル

(1≦k≦K)が入力端子S_k(1≦k≦K)を介して乗算部１２に入力される。変調信号ベクトル

から

は、ベースバンド信号である変調信号を要素とするベクトルである。以下ではベースバンド信号は全て、同相（in-phase）成分を実部、直交（quadrature-phase）成分を虚部とする複素表示で表わすことにする。乗算部１２は、これらの変調信号にプリコーディング行列生成部１４が生成するプリコーディング行列を乗算して、送信信号を生成する。時刻iは、変調のシンボル周期分、次の時刻i＋１と離れた離散時刻である。

送信信号は下記の式（１）で表される。

ここで、

はN_T次元の送信信号ベクトルであり、

は、プリコーディング行列生成部１４が生成するプリコーディング行列である。

このプリコーディング行列はN_T行N_r列であり、N_rは受信側であるすべてのＫ個のユーザ端末の受信アンテナの総数であり、次式で表される。

また、式（１）の

は、前記の変調信号ベクトル

の集合から得られるN_r次元変調信号ベクトルであって、次式で表される。

ここで上付文字のTは、転置（transposition）を示す。

プリコーディング行列生成部１４には、ユーザ端末からフィードバックされた下りリンクの伝送路情報が入力される。下りリンクの伝送路情報は、N_r行N_T列のチャネル行列

であり、このチャネル行列の各要素は下りリンクのインパルス応答である。プリコーディング行列生成部１４は、このチャネル行列を用いて、プリコーディング行列

を生成し、乗算部１２へプリコーディング行列を入力する。

線形プリコーディングには、ゼロフォーシング(Zero Forcing、ZF)規範プリコーディングと、ＭＭＳＥ（最小平均二乗誤差）規範プリコーディングがある。

ZFプリコーディングの場合、プリコーディング行列は次式で得られる。

ここで、上付文字のＨは、複素共役転置（complex conjugate transposition）すなわちエルミート転置（Hermitian transpose）を示す。γは送信電力を一定にするための規格化定数である。

この場合のプリコーディング行列

は、チャネル行列

の一般逆行列（generalized inverse）に相当するので、ZFは、チャネルの逆特性を変調信号に乗算することにより、ユーザ端末で他のユーザ端末への変調信号が漏れ込まないようにすることができる。即ち、ユーザ端末で、同一周波数での他のユーザ端末への信号に起因する干渉をゼロにできる。

この無線送受信機は、送信アンテナ１０（１０_１〜１０_NT）の数と同じN_T個の直交変調器１６（１６_１〜１６_NT）を有する。乗算部１２が生成するN_T次元の送信信号ベクトル

の第１要素は直交変調器１６_１へ供給され、送信信号ベクトルの第N_T要素は直交変調器１６_NTへ供給される。つまり、N_T次元の送信信号ベクトルの要素は、N_T個の直交変調器１６にそれぞれ入力される。

各直交変調器１６は、発振器１８が生成するRF周波数の搬送波を、端子Minを介して入力される送信信号の要素で変調し、その直交変調器１６に対応する送信アンテナ１０によって無線送信する。

図２は、各直交変調器１６の構成を示す。各直交変調器１６は、D/A変換器２０_I、D/A変換器２０_Q、乗算器２２_I、乗算器２２_Q、移相器２４、加算器２６およびアンプ２８を備える。端子Minから送信信号の要素が入力される。この送信信号の要素はディジタル信号であり、同相成分がD/A変換器２０_Iでアナログ信号に変換され、直交成分がD/A変換器２０_Qでアナログ信号に変換される。端子CinからRF周波数の搬送波が入力される。乗算器２２_Iでは同相成分のアナログ信号と搬送波が乗算される。乗算器２２_Qでは直交成分のアナログ信号と、移相器２４で位相が90度回転した搬送波が乗算される。加算器２６は、これらの乗算結果を足し合わせ、変調された搬送波を生成する。変調された搬送波はアンプ２８で増幅され、端子Moutから出力される。

図１に示す無線送受信機におけるプリコーディング行列生成部１４は、ユーザ端末からフィードバックされたチャネル行列に基づいて、プリコーディング行列を生成する。しかしながら、チャネル行列のような情報量が大きい伝送路情報をユーザ端末からフィードバックすると、上りリンクをそのために使用し、上りリンクの伝送効率を劣化させてしまう。そこで、このフィードバックを不要とし、上りリンクの受信信号から伝送路情報を抽出する方法が検討されている。

図３は、この方法に適するプリコーディング行列生成部１４Ａの構成を示す。プリコーディング行列生成部１４Ａは、到来方向（direction-of-arrival、DOA）推定部３０とプリコーディング行列推定部３２を備える。

プリコーディング行列生成部１４Ａにおいては、まず、チャネル行列の代わりに、複数のユーザ端末から受信される上りリンクの受信信号のベクトルが入力される。DOA推定部３０は、上りリンクの受信信号のベクトルを用いて、各ユーザ端末の上りリンクのパスの到来角を推定する。この推定にはMUSIC(MUltiple SIgnal Classification)などの到来方向推定アルゴリズムが用いられる（非特許文献１参照）。

プリコーディング行列推定部３２は、DOA推定部３０で推定された各ユーザ端末の上りリンクの到来角の推定値に基づいて、プリコーディング行列

を生成する。プリコーディング行列推定部３２は、各ユーザ端末について、上りリンクのパスの角度と下りリンクのパスの角度とが一致するという仮定の下で、プリコーディング行列を下記の式（２）に従って計算する。

ここで、

ただし、式（２）は、各ユーザ端末の受信アンテナの数N_rk = 1、受信側であるすべてのＫ個のユーザ端末の受信アンテナの総数N_r = Kという想定の下で使用される。ここで、

は、式（２）で使用されるK行N_T列の行列である。

は、無線送受信機の送信アンテナ１０のN_T次元アレー・レスポンス・ベクトルであり、

は第kユーザ端末の到来角推定値であり、γ’は送信電力を一定にするための規格化定数である。

図４は、図１の無線送受信機におけるプリコーディング行列生成部１４の代わりに使用されうる従来の他のプリコーディング行列生成部１４Ｂの構成を示す。プリコーディング行列生成部１４Ｂは、自己相関行列推定部３４、自己相関行列変換部３６、固有値分解部３８、および固有ベクトル選択部３９を備える。

プリコーディング行列生成部１４Ｂにおいては、まず、複数のユーザ端末から受信される上りリンクの受信信号のベクトルが入力される。離散時刻iにおける第kユーザ端末からの受信信号は、N_T次元ベクトル

で表わすことができる。自己相関行列推定部３４は、受信信号のベクトルの積、

を時間平均し、上りリンクの自己相関行列（auto-correlation matrix）を求める。

自己相関行列変換部３６は、この上りリンクの自己相関行列から下りリンクの自己相関行列を推定する。この推定法は、各ユーザ端末について上りリンクのパスの角度と下りリンクのパスの角度とが一致するという仮定の下で、上りリンクと下りリンクのアレー・レスポンス・ベクトルの違いを考慮して、上りリンクの自己相関行列を線形変換することで、下りリンクの自己相関行列を推定する（非特許文献２参照）。

固有値分解部３８は下りリンクの自己相関行列を固有値分解（eigenvalue decomposition）し、N_T次元の固有ベクトルを生成する。固有ベクトル選択部３９は、固有値の大きい固有ベクトル（eigenvector）を選択し、プリコーディング行列の列ベクトルに設定し、得られたプリコーディング行列に最後に規格化定数を乗算して最終的にプリコーディング行列

を求める（非特許文献３参照）。

S. Haykin, "Adaptive Filter Theory (3rd Edition)", Prentice-Hall, U.S.A., December 27, 1995. Y-C Liang and F. P. S. Chin, "Downlink channel covariance matrix (DCCM) estimation and its applications in wireless DS-CDMA systems", IEEE Journal on Selected Areas in Communications, vol. 19, no. 2, pp. 222-232, February 2001. L. Dong, "Open-loop beamforming for frequency-division duplex mobile wireless access", IEEE Transactions on Vehicular Technology, Vol. 56, No. 4, pp. 1845-1849, July 2007.

図４に示すプリコーディング行列生成部１４Ｂによりプリコーディング行列を求めると、ユーザ端末の受信信号電力を高めることができる。しかし、他のユーザ端末への信号に起因する干渉を抑圧する効果は得られず、同一周波数での干渉（同一チャネル干渉）のためユーザ端末で受信される信号の受信品質が劣化する。

そこで、本発明は、同一チャネル干渉を抑圧し、ユーザ端末で受信される信号の受信品質の劣化を抑えられる、伝送路情報のフィードバック不要のプリコーディング方式を実行する無線基地局を提供する。

本発明に係る無線基地局は、上りリンクの受信信号のベクトルの積を平均化して上りリンクの自己相関行列を求め、前記上りリンクの自己相関行列から下りリンクの自己相関行列を推定する下りリンク自己相関行列生成部と、前記下りリンク自己相関行列生成部で推定された前記下りリンクの自己相関行列のうち、所望ユーザ端末の下りリンクの自己相関行列を固有値分解して固有ベクトルを求めるとともに、前記下りリンク自己相関行列生成部で推定された前記下りリンクの自己相関行列のうち、全ユーザ端末または他ユーザ端末の下りリンクの自己相関行列に基づいて、他ユーザ端末への送信信号に起因する同一チャネル干渉を抑圧するための干渉抑圧行列を生成し、前記固有ベクトルと前記干渉抑圧行列に基づいて、プリコーディング行列を生成するプリコーディング行列変換部と、変調信号ベクトルに前記プリコーディング行列を乗算して、下りリンクの送信信号ベクトルを生成する送信信号生成部とを備える。

本発明によれば、下りリンク自己相関行列生成部で推定された下りリンクの自己相関行列のうち、全ユーザ端末または他ユーザ端末の下りリンクの自己相関行列に基づいて、他ユーザ端末への送信信号に起因する同一チャネル干渉を抑圧するための干渉抑圧行列が生成される。この干渉抑圧行列と、所望ユーザ端末の下りリンクの自己相関行列を固有値分解して得られた固有ベクトルに基づいて、プリコーディング行列を生成することによって、他ユーザ端末への送信信号に起因する同一チャネル干渉を抑圧し、所望ユーザ端末で受信される信号の受信品質の劣化を抑えることができる。

前記下りリンク自己相関行列生成部は、前記上りリンクの自己相関行列に線形変換を施すことにより、前記下りリンクの自己相関行列を推定する下りリンク自己相関行列変換部を備えてもよい。

このように、下りリンクの自己相関行列を推定する過程において、上りリンクの自己相関行列に線形変換を施すことにより、精度良く下りリンクの自己相関行列を推定することができる。

前記下りリンク自己相関行列生成部は、前記上りリンクの自己相関行列から、各ユーザ端末について上りリンクのパスの到来角分布を推定する到来角分布推定部と、
前記到来角分布から前記下りリンクの自己相関行列を推定する自己相関行列演算部とを備えてもよい。

このように到来角分布を推定し、到来角分布から下りリンクの自己相関行列を推定することにより、マルチパス伝搬環境でパスに角度広がりがある場合でも、精度良く下りリンクの自己相関行列を推定することができる。

前記下りリンク自己相関行列生成部は、前記上りリンクの自己相関行列から、各ユーザ端末について上りリンクのパスの到来角分布を推定するとともに、各ユーザ端末について前記到来角分布を所定の種類の確率分布で近似させて、各ユーザ端末について前記確率分布の平均値と分散を推定する到来方向推定部と、前記平均値と分散から前記下りリンクの自己相関行列を推定する自己相関行列演算部とを備えてもよい。

このように到来角分布を推定し、到来角分布から下りリンクの自己相関行列を推定することにより、マルチパス伝搬環境でパスに角度広がりがある場合でも、精度良く下りリンクの自己相関行列を推定することができる。また、各ユーザ端末について前記到来角分布を所定の種類の確率分布で近似させて、各ユーザ端末について前記確率分布の平均値と分散を推定することにより、上りリンクの受信信号に含まれる雑音の影響で誤差が大きくなる場合でも、誤差を軽減することができる。

前記プリコーディング行列変換部は、前記下りリンク自己相関行列生成部で推定された前記下りリンクの自己相関行列のうち、所望ユーザ端末の下りリンクの自己相関行列を固有値分解して固有ベクトルを求める固有値分解部と、前記固有ベクトルのうち、固有値が大きい方の１以上の固有ベクトルを選択し、選択された固有ベクトルの位相を回転させる固有ベクトル選択位相回転部と、前記下りリンク自己相関行列生成部で推定された前記下りリンクの自己相関行列のうち、全ユーザ端末または他ユーザ端末の下りリンクの自己相関行列の逆行列である前記干渉抑圧行列を生成する逆行列演算部と、前記固有ベクトル選択位相回転部で得られた位相が回転させられた固有ベクトルに、前記逆行列演算部で得られた前記干渉抑圧行列を乗算し、前記プリコーディング行列を生成する乗算部とを備えてもよい。

このように、所望ユーザ端末に関する固有ベクトルに、全ユーザ端末または他ユーザ端末の下りリンクの自己相関行列の逆行列（干渉抑圧行列）を乗算し、プリコーディング行列を生成することにより、他ユーザ端末への送信信号に起因する同一チャネル干渉を抑圧することができる。

前記固有ベクトル選択位相回転部は、前記固有値分解部で求められた前記固有ベクトルのうち、前記所望ユーザ端末の受信アンテナの数と等しい数の固有ベクトルを選択してもよい。

このように、所望ユーザ端末の受信アンテナの数と等しい数の固有ベクトルを固有ベクトル選択位相回転部が選択することにより、乗算部で生成されるプリコーディング行列は、ユーザ端末あたり複数の送信信号を送信するのに適合させられる。

前記プリコーディング行列変換部は、前記下りリンク自己相関行列生成部で推定された前記下りリンクの自己相関行列のうち、他ユーザ端末の下りリンクの自己相関行列を固有値分解して、固有値分解結果から固有値が小さい方の複数の第１の固有ベクトルを選択する第１の固有値分解部と、前記下りリンク自己相関行列生成部で推定された前記下りリンクの自己相関行列のうち、所望ユーザ端末の下りリンクの自己相関行列を、前記第１の固有ベクトルを用いて変換し、変換された下りリンクの自己相関行列を固有値分解し、固有値分解結果から固有値が大きい方の１以上の第２の固有ベクトルを選択する第２の固有値分解部と、前記第１の固有値分解部で選択された前記第１の固有ベクトルを含む行列である前記干渉抑圧行列と前記第２の固有値分解部で選択された前記第２の固有ベクトルを含む行列である所望信号増強行列とを乗算し、前記プリコーディング行列を生成する乗算部とを備えてもよい。

ここで、干渉抑圧行列は、他ユーザ端末に対応し、固有値が小さい方の複数の第１の固有ベクトルを含むため、他ユーザ端末への送信信号に起因する同一チャネル干渉を抑圧する働きを有する。他方、所望信号増強行列は、所望ユーザ端末に対応し、固有値が大きい１以上の第２の固有ベクトルを含むため、所望ユーザ端末にとっての所望信号を強める働きを有する。干渉抑圧行列と所望信号増強行列を乗算し、プリコーディング行列を生成することにより、他ユーザ端末への送信信号に起因する同一チャネル干渉を抑圧することができる。

前記第２の固有値分解部は、前記下りリンク自己相関行列生成部で推定された前記下りリンクの自己相関行列のうち、所望ユーザ端末の下りリンクの自己相関行列を、前記第１の固有ベクトルを用いて変換し、変換された下りリンクの自己相関行列を固有値分解し、固有値分解結果から前記所望ユーザ端末の受信アンテナの数と等しい数の、固有値が大きい方の前記第２の固有ベクトルを選択してもよい。

このように、所望ユーザ端末の受信アンテナの数と等しい数の第２の固有ベクトルを第２の固有値分解部が選択することにより、乗算部で生成されるプリコーディング行列は、ユーザ端末あたり複数の送信信号を送信するのに適合させられる。

従来のプリコーディングを用いたMIMOシングルキャリア無線送受信機（無線基地局）の構成を示すブロック図である。図１の無線送受信機における直交変調器の構成を示すブロック図である。図１の無線送受信機におけるプリコーディング行列生成部の代わりに使用されうる従来のプリコーディング行列生成部の構成を示すブロック図である。図１の無線送受信機におけるプリコーディング行列生成部の代わりに使用されうる従来の他のプリコーディング行列生成部の構成を示すブロック図である。本発明の実施の形態に係るプリコーディングを用いたMIMOシングルキャリア無線送受信機（無線基地局）の構成を示すブロック図である。図５の無線送受信機におけるプリコーディング行列生成部の構成を示すブロック図である。図６のプリコーディング行列生成部で使用される自己相関行列生成部の一例を示すブロック図である。図６のプリコーディング行列生成部で使用される自己相関行列生成部の他の例を示すブロック図である。図６のプリコーディング行列生成部で使用される自己相関行列生成部の他の例を示すブロック図である。図６のプリコーディング行列生成部で使用されるプリコーディング行列変換部の一例を示すブロック図である。図６のプリコーディング行列生成部で使用されるプリコーディング行列変換部の他の例を示すブロック図である。本発明の実施の形態に係る無線通信システムを示す概略図である。

以下、本発明に係る様々な実施の形態について、マルチユーザMIMOかつシングルキャリア伝送を例に説明する。

図５は、本発明の実施の形態に係るプリコーディングを用いたMIMOシングルキャリア無線送受信機（無線基地局）の構成を示すブロック図である。図５は、無線基地局のうち無線送信に関する部分のみを示す。この無線基地局は、送信アンテナ１０、乗算部（送信信号生成部）１２、プリコーディング行列生成部１４０、直交変調器１６および発振器１８を備える。この無線基地局は、図１の無線基地局と類似する構成を有するが、プリコーディング行列生成部１４の代わりにプリコーディング行列生成部１４０を有する。送信アンテナ１０、乗算部１２、直交変調器１６および発振器１８は、図１に示す背景技術の無線基地局のそれらと同じでよい。無線基地局は、送信アンテナ１０の数N_Tと等しい数の受信アンテナ１１（図１２参照）を有する。無線基地局に送信アンテナ１０と受信アンテナを兼ねる複数の送受信アンテナを設けてもよい。

プリコーディング行列生成部１４０および乗算部１２は、無線基地局の図示しないＣＰＵ（central processing unit）がコンピュータプログラムを実行し、そのコンピュータプログラムに従って機能することによって実現される機能ブロックである。後述するプリコーディング行列生成部１４０の構成要素も同様である。但し、ＣＰＵが実行する各機能は、ＣＰＵの代わりに、ハードウェアで実行してもよいし、例えばFPGA（Field Programmable Gate Array）、DSP（Digital Signal Processor）等のプログラマブルロジックデバイスで実行してもよい。

図６は、本発明の実施の形態に係るプリコーディング行列生成部の構成を示すブロック図である。プリコーディング行列生成部１４０は、下りリンク自己相関行列生成部４０およびプリコーディング行列変換部６０を備える。

プリコーディング行列生成部１４０においては、まず、複数のユーザ端末から受信される上りリンクの受信信号のベクトルが、下りリンク自己相関行列生成部４０へ入力される。下りリンク自己相関行列生成部４０は、上りリンクの受信信号のベクトルの積を平均化して上りリンクの自己相関行列を求め、上りリンクの自己相関行列から下りリンクの自己相関行列を推定する。下りリンク自己相関行列生成部４０は、推定された下りリンクの自己相関行列をプリコーディング行列変換部６０へ供給する。

プリコーディング行列変換部６０は、下りリンク自己相関行列生成部４０で推定された下りリンクの自己相関行列のうち、所望ユーザ端末の下りリンクの自己相関行列を固有値分解して固有ベクトルを求めるとともに、下りリンク自己相関行列生成部４０で推定された下りリンクの自己相関行列のうち、全ユーザ端末または他ユーザ端末の下りリンクの自己相関行列に基づいて、他ユーザ端末への送信信号に起因する同一チャネル干渉を抑圧するための干渉抑圧行列を生成し、固有ベクトルと干渉抑圧行列に基づいて、プリコーディング行列を生成する。プリコーディング行列の生成においては、プリコーディング行列変換部６０は送信電力が一定になるよう規格化を行う。下りリンクの自己相関行列から固有ベクトル及び干渉抑圧行列を生成することで、後述するように、同一周波数での他のユーザ端末への信号に起因する干渉を抑圧するようプリコーディング行列を求めることができる。

図７は、図６の下りリンク自己相関行列生成部４０の構成の一例を示す。図７の下りリンク自己相関行列生成部４０は、上りリンク自己相関行列推定部４２および下りリンク自己相関行列変換部４４を備える。

図７の下りリンク自己相関行列生成部４０においては、まず、複数のユーザ端末から受信される上りリンクの受信信号のベクトルが、上りリンク自己相関行列推定部４２へ入力される。上りリンク自己相関行列推定部４２は上りリンクの受信信号のベクトルの積を平均化して上りリンクの自己相関行列を求める。この上りリンク自己相関行列推定部４２は、図４を参照して背景技術に関連して説明した自己相関行列推定部３４と同じでよく、精度良く上りリンクの自己相関行列を推定できる。

下りリンク自己相関行列変換部４４は、上りリンク自己相関行列推定部４２で推定された上りリンクの自己相関行列に線形変換を施すことで、下りリンクの自己相関行列を推定する。この下りリンク自己相関行列変換部４４は、図４を参照して背景技術に関連して説明した自己相関行列変換部３６と同じでよい。このように、下りリンクの自己相関行列を推定する過程において、上りリンクの自己相関行列に線形変換を施すことにより、精度良く下りリンクの自己相関行列を推定することができる。

図８は、図６の下りリンク自己相関行列生成部４０の構成の他の一例を示す。図８の下りリンク自己相関行列生成部４０は、上りリンク自己相関行列推定部４２、到来角分布推定部４８および下りリンク自己相関行列演算部５０を備える。

図８の下りリンク自己相関行列生成部４０においては、まず、複数のユーザ端末から受信される上りリンクの受信信号のベクトルが、上りリンク自己相関行列推定部４２へ入力される。上りリンク自己相関行列推定部４２は上りリンクの受信信号のベクトルの積を平均化して上りリンクの自己相関行列を求める。この上りリンク自己相関行列推定部４２は、図４を参照して背景技術に関連して説明した自己相関行列推定部３４と同じでよく、精度良く上りリンクの自己相関行列を推定できる。

到来角分布推定部４８は、上りリンク自己相関行列推定部４２で推定された上りリンクの自己相関行列から、上りリンクの複数のパス（複数のユーザ端末に相当する）の到来角分布を推定する。下りリンク自己相関行列演算部５０は、到来角分布推定部４８で推定された到来角分布から下りリンクの自己相関行列を推定する。

到来角分布推定部４８において、到来角分布はMVDR (Minimum Variance Distortionless Response)に基づき、推定される。以下では数式を用いてMVDRを説明する。第kユーザ端末からの到来角分布をp_k(φ)とし、第kユーザ端末の上りリンクの自己相関行列を

とする。

MVDRに従って、到来角分布推定部４８は、p_k(φ)を下記の式（３）に従って計算すなわち推定する。

ここで、

は、無線送受信機（無線基地局）の受信アンテナ１１（図１２参照）のN_T次元アレー・レスポンス・ベクトルであり、Cは規格化定数であって、次の式を満足するように求められている。

次に、下りリンク自己相関行列演算部５０が実行する下りリンクの自己相関行列の推定について、数式を用いて説明する。第kユーザ端末の下りリンクの自己相関行列を

とする。各ユーザ端末について上りリンクのパスの角度と下りリンクのパスの角度とが一致すると仮定すると、
自己相関行列

は下記の式（４）に従って計算することができる。

なお、

は、無線送受信機の送信アンテナ１０のN_T次元アレー・レスポンス・ベクトルである。

ここで、P_kは第kユーザ端末の平均受信電力であり、送信電力制御が行われていれば既知の値である。式（３）の到来角分布p_k(φ)を式（４）に代入することにより、自己相関行列

を求めることができる。

このように、図８に示す下りリンク自己相関行列生成部４０は、到来角分布を求め、下りリンクの自己相関行列を推定するので、マルチパス伝搬環境でパスに角度広がりがある場合でも、精度良く下りリンクの自己相関行列を推定することができる。

図３に示す背景技術のプリコーディング行列生成部１４Ａが使用する式（２）は、マルチパス伝搬環境でパスの角度に広がりが無いもしくは非常に小さいという想定を前提とする。したがって、式（２）は、マルチパス伝搬環境でパスの角度に広がりがあったとしても、到来角の平均値しか用いていないことになる。そのため、このプリコーディング行列生成部１４Ａで得られたプリコーディング行列を用いて無線送信を行うと、ユーザ端末で受信された信号の受信品質が劣化する。図８に示す下りリンク自己相関行列生成部４０は、この問題を解決する。

図９は、図６の下りリンク自己相関行列生成部４０の構成の他の一例を示す。図９の下りリンク自己相関行列生成部４０は、上りリンク自己相関行列推定部４２、到来方向（direction-of-arrival、DOA）推定部５４および下りリンク自己相関行列演算部５６を備える。

図９の下りリンク自己相関行列生成部４０においては、まず、複数のユーザ端末から受信される上りリンクの受信信号のベクトルが、上りリンク自己相関行列推定部４２へ入力される。上りリンク自己相関行列推定部４２は上りリンクの受信信号のベクトルの積を平均化して上りリンクの自己相関行列を求める。この上りリンク自己相関行列推定部４２は、図４を参照して背景技術に関連して説明した自己相関行列推定部３４と同じでよく、精度良く上りリンクの自己相関行列を推定できる。

DOA推定部５４は、上りリンク自己相関行列推定部４２で推定された上りリンクの自己相関行列から、各ユーザ端末について上りリンクのパスの到来角の平均値と分散を推定し、これらを下りリンク自己相関行列演算部５６に供給する。具体的には、DOA推定部５４は、まず図８の到来角分布推定部４８と同様に、式（３）に従って第kユーザ端末の到来角分布p_k(φ)を求める。さらに、DOA推定部５４は、到来角分布p_k(φ)を所定の種類の確率分布（例えば正規分布）で近似して、その確立分布の平均値φ_kと分散σ_k ²を求める。

下りリンク自己相関行列演算部５６は、この到来角の平均値と分散から、下りリンクの自己相関行列

を推定すなわち生成する。具体的には、式（４）のp_k(φ)を平均値φ_k、分散σ_k ²の正規分布に置換える。即ち、下りリンクの自己相関行列は下記の式（５）に従って計算される。

図９に示す下りリンク自己相関行列生成部４０によれば、到来角分布p_k(φ)に、上りリンクの受信信号に含まれる雑音の影響で誤差が大きくなる場合、正規分布とマッチングをとり近似することで、この誤差を軽減できる。したがって、マルチパス伝搬環境でパスに角度広がりがある場合でも、さらに精度良く下りリンクの自己相関行列を推定することができる。図９に示す下りリンク自己相関行列生成部４０のDOA推定部５４は、到来角分布を正規分布で近似するが、ラプラス分布等の他の確率密度関数でマッチングをとり近似してもよい。

図３に示す背景技術のプリコーディング行列生成部１４Ａが使用する式（２）は、マルチパス伝搬環境でパスの角度に広がりが無いもしくは非常に小さいという想定を前提とする。したがって、式（２）は、マルチパス伝搬環境でパスの角度に広がりがあったとしても、到来角の平均値しか用いていないことになる。そのため、このプリコーディング行列生成部１４Ａで得られたプリコーディング行列を用いて無線送信を行うと、ユーザ端末で受信された信号の受信品質が劣化する。図９に示す下りリンク自己相関行列生成部４０は、この問題を解決する。

図１０は、図６のプリコーディング行列変換部６０の構成の一例を示す。図１０のプリコーディング行列変換部６０は、固有値分解部６２、固有ベクトル選択位相回転部６４、逆行列演算部６６および乗算正規化部６８を備える。

図１０のプリコーディング行列変換部６０においては、図７から図９のいずれかの下りリンク自己相関行列生成部４０によって推定された下りリンクの自己相関行列

が入力される。

固有値分解部６２は、下りリンク自己相関行列生成部４０によって推定された下りリンクの自己相関行列のうち、所望ユーザ端末（第kユーザ端末）の下りリンクの自己相関行列を固有値分解して固有ベクトルを求める。固有ベクトル選択位相回転部６４は、固有値が大きい方のいくつかの固有ベクトル（第kユーザ端末の受信アンテナの数と等しい数の固有ベクトル）を選択し、それらの固有ベクトルの位相を回転させて、乗算正規化部６８に供給する。

一方、逆行列演算部６６は、下りリンク自己相関行列生成部４０によって推定された下りリンクの自己相関行列のうち、全ユーザ端末または他ユーザ端末の下りリンクの自己相関行列の逆行列（干渉抑圧行列）を生成し、これを乗算正規化部６８に供給する。乗算正規化部（乗算部）６８は、固有ベクトル選択位相回転部６４で得られた位相が回転させられた固有ベクトルと逆行列演算部６６の逆行列を乗算することで、プリコーディング行列を生成する。乗算正規化部６８は、さらに送信電力が一定になるようその乗算結果に規格化定数を乗算してプリコーディング行列を生成してもよい。

図１０のプリコーディング行列変換部６０で求められるプリコーディング行列

を具体的に説明する。このプリコーディング行列は、複数のユーザ端末のためのプリコーディング行列の集合として、下記の式（６）で表される。

ここで、

は第kユーザ端末への送信信号のためのプリコーディング行列であって、N_T行N_rk列の行列である。

は、下記の式（７）に従って計算される。

ここで、gは送信電力を一定にするための規格化定数であり乗算正規化部６８により与えられる。

は逆行列演算部６６で生成される逆行列であり、

は固有ベクトル選択位相回転部６４で得られた位相が回転させられた固有ベクトルである。

ηは正の定数であり、

は、N_T行N_T列の単位行列である。但し、この定数と単位行列の積の加算は、本発明にとって必ずしも不可欠ではない（式（８）でも同じである）。

Lは、第kユーザ端末の受信アンテナの数である（L = N_rk）。添字k1は全ユーザ端末に関する総和を計算するに付加されており、付加されたパラメータの意味を変えるものではない（式（８）でも同じである）。また、θ_kq(1≦q≦L) は位相回転量であって、固有ベクトル選択位相回転部６４により与えられる。jは虚数単位である。

は、下りリンクの自己相関行列

の正規化固有ベクトルの内、q番目に固有値が大きい正規化固有ベクトルである。

第kユーザ端末への送信信号のためのプリコーディング行列

は、式（７）の代わりに、下記の式（８）に従って計算してもよい。

は逆行列演算部６６で生成される逆行列である。

式（７）は全ユーザ端末に対応しており、式（８）は、k1≠kから理解されるように、第kユーザ端末を除くすべてのユーザ端末（他ユーザ端末、所望ユーザ端末以外のすべての端末）に対応する。

式（７）または式（８）において、正規化固有ベクトル

に、送信電力に相当する正の重み付け係数を乗算することもできる。

さらに、式（７）または式（８）において、正規化固有ベクトル

をパスの到来角の方向を向いたアレー・レスポンス・ベクトル

に置換えることもできる。

式（７）および式（８）で逆行列が乗算されている。これは、他ユーザ端末への送信信号から所望ユーザ端末への送信信号への漏れ込みを抑える、即ち同一チャネル干渉を抑える効果を達成する。図４を参照して説明した背景技術のプリコーディング行列生成部１４Ｂでは、このような効果は達成できない。

このように、図１０に示すプリコーディング行列変換部６０は、同一チャネル干渉を抑える効果を達成することができる。また、式（７）と式（８）で、L = N_rk個（第kユーザ端末の受信アンテナの数）の固有ベクトル

を選択していることから明らかなように、ユーザ端末あたり複数の送信信号を送信することができる。図３に示す背景技術のプリコーディング行列生成部１４Ａが使用する式（２）は、各ユーザ端末の受信アンテナの数N_rkが１であるという想定の下で使用されるため、あるユーザ端末の受信アンテナの数N_rkが２以上の場合に、そのユーザ端末に向けて複数の送信信号を送ることができない。図１０に示すプリコーディング行列変換部６０はこの問題を解決する。

図１１は、図６のプリコーディング行列変換部６０の構成の他の一例を示す。図１１のプリコーディング行列変換部６０は、第１の固有値分解部７０、第２の固有値分解部７２および乗算正規化部７４を備える。

図１１のプリコーディング行列変換部６０においては、図７から図９のいずれかの下りリンク自己相関行列生成部４０によって推定された下りリンクの自己相関行列

が入力される。

第１の固有値分解部７０は、下りリンク自己相関行列生成部４０によって推定された下りリンクの自己相関行列のうち、他ユーザ端末の下りリンクの自己相関行列を固有値分解して、固有値分解結果から固有値が小さい方のいくつかの固有ベクトル（第１の固有ベクトル）を選択する。

第２の固有値分解部７２は、下りリンク自己相関行列生成部４０によって推定された下りリンクの自己相関行列のうち、所望ユーザ端末の下りリンクの自己相関行列を、第１の固有値分解部７０で選択された複数の第１の固有ベクトルを用いて変換し、変換された下りリンクの自己相関行列を固有値分解し、固有値分解結果から固有値が大きい方のいくつかの固有ベクトル（第２の固有ベクトル）を選択する。選択される第２の固有ベクトルの数は、所望ユーザ端末の受信アンテナの数と等しい。

乗算正規化部（乗算部）７４は、第１の固有値分解部７０で選択された固有値が小さい方のいくつかの固有ベクトル（第１の固有ベクトル）を含む第１の行列（干渉抑圧行列）と、第２の固有値分解部７２で選択された固有値が大きい方のいくつかの固有ベクトル（第２の固有ベクトル）を含む第２の行列（所望信号増強行列）とを乗算することで、プリコーディング行列を生成する。乗算正規化部７４は、さらに送信電力が一定になるようその乗算結果を規格化してプリコーディング行列を生成してもよい。

図１１のプリコーディング行列変換部６０で求められるプリコーディング行列

を具体的に説明する。このプリコーディング行列も、複数のユーザ端末のためのプリコーディング行列の集合として、上記の式（６）で表される。
ここで、

は、下記の式（９）に従って計算される。

式（９）において、

は、N_T行Q_k列の行列（第１の行列）である(Q_k≧N_rk)。この第１の行列は、下記の式（１０）で定めるN_T行N_T列の行列

の正規化固有ベクトルの内、固有値の小さいQ_k個の正規化固有ベクトル（第１の固有ベクトル）を列ベクトルとして持つ。式（１０）において、*は複素共役（complex conjugation）を表わす。添字k1は複数ユーザ端末に関する総和を計算するに付加されており、付加されたパラメータの意味を変えるものではない。

このように、第１の行列（干渉抑圧行列）

は、第１の固有値分解部７０により生成される行列である。

また、式（９）において、

は、Q_k行N_rk列の行列（第２の行列）である。この第２の行列は、下記の式（１１）で定めるQ_k行Q_k列の行列

の正規化固有ベクトルの内、固有値が大きい方のN_rk個の正規化固有ベクトル（第２の固有ベクトル）を列ベクトルとして持つ。

このように、第２の行列（所望信号増強行列）

は、第２の固有値分解部７２により生成される行列である。

式（９）において、

は、N_rk行N_rk列の対角行列であり、対角成分は各送信信号の送信電力の平方根であり、規格化定数に相当する。これは、乗算正規化部７４により、規格化のために与えられる。但し、これは本発明にとって必ずしも不可欠ではない。

図１１に示すプリコーディング行列変換部６０が使用する式（９）において、第１の行列（干渉抑圧行列）

は、第kユーザ端末にとっての干渉信号を抑圧するための行列である。その理由は次の通りである。この第１の行列の基礎となる行列

は、式（１０）のk1≠kから理解されるように、第kユーザ端末を除くすべてのユーザ端末（他ユーザ端末）に対応する。そして、第１の行列

は、基礎となる行列の正規化固有ベクトルの内、固有値が小さい正規化固有ベクトル（第１の固有ベクトル）を列ベクトルとして有するため、信号を弱めるように働く。したがって、第１の行列

は、他ユーザ端末への送信信号に起因する同一チャネル干渉を大幅に抑圧することができる。

図１１に示すプリコーディング行列変換部６０が使用する式（９）において、第２の行列（所望信号増強行列）

は、第kユーザ端末にとっての所望信号を強める成分である。その理由は次の通りである。この第２の行列の基礎となる行列

は、第kユーザ端末（所望ユーザ端末）に対応する。そして、第２の行列

は、基礎となる行列の正規化固有ベクトルの内、固有値が大きい正規化固有ベクトル（第２の固有ベクトル）を列ベクトルとして有するため、信号を強めるように働く。したがって、第２の行列

は、第kユーザ端末の受信電力を高めることができる。

以上の通り、図１１に示すプリコーディング行列変換部６０によれば、同一チャネル干渉を大幅に抑圧することができるとともに、第kユーザ端末の受信電力を高めることができる。

また、所望ユーザ端末の受信アンテナの数N_rkと等しい数の第２の固有ベクトルを第２の固有値分解部７２が選択することにより、乗算正規化部７４で生成されるプリコーディング行列は、ユーザ端末あたり複数の送信信号を送信するのに適合させられる。式（９）で求められる各ユーザ端末のためのプリコーディング行列の列数はN_rk個（第kユーザ端末の受信アンテナの数）であるので、ユーザ端末あたり複数の送信信号を送信することができる。図３に示す背景技術のプリコーディング行列生成部１４Ａが使用する式（２）は、各ユーザ端末の受信アンテナの数N_rkが１であるという想定の下で使用されるため、あるユーザ端末の受信アンテナの数N_rkが２以上の場合に、そのユーザ端末に向けて複数の送信信号を送ることができない。図１０に示すプリコーディング行列変換部６０はこの問題を解決する。

以上、本発明を実施するための最良の形態についてマルチユーザMIMOのシングルキャリア伝送を例に説明したが、OFDM(Orthogonal Frequency Division Multiplexing)伝送の場合でも本発明を実施しうる。OFDMの場合、各サブキャリアは狭い帯域幅を有するので、サブキャリア内の実質的なフェージングはフラットフェージングとなる。したがって、各サブキャリアについて上記の操作を行えばよい。

図１２は、本発明の実施の形態に係る無線通信システムを示す。無線通信システムは、複数の無線基地局１を備え、各無線基地局１はその無線基地局のセル２にあるユーザ端末（移動局）５と通信する。各無線基地局１は、複数の送信アンテナ（アダプティブ・アレー・アンテナ）１０を備え、各ユーザ端末５は1以上の受信アンテナ５Ａを備える。したがって、各無線基地局１はMIMO送信が可能である。また、各無線基地局１は、複数の受信アンテナ１１を備える。

ユーザ端末５は例えば携帯電話（UMTS（Universal Mobile Telecommunications System）LTE（Long Term Evolution）でのUE（User Equipment）である。各無線基地局１は、UMTS LTEでのeNB（evoluted Node B）であってもよい。あるいは、各無線基地局１は、無線LAN（Local Area Network）またはWiMAX（Worldwide Interoperability for Microwave Access）でのアクセスポイントであってもよい。

以上、本発明をその好適な実施の形態を参照しながら詳細に図示して説明したが、請求の範囲に記載された本発明の区域内で、形式および細部に関する様々な変更が可能であることは当業者であれば理解できることだろう。かかる変更、代替、修正も本発明の範囲に含まれるものであると出願人は意図している。

１無線基地局、２セル、５ユーザ端末、５Ａ受信アンテナ、１０送信アンテナ、１１受信アンテナ、１２乗算部（送信信号生成部）、１６直交変調器、１８発振器、３０ DOA推定部、１４０プリコーディング行列生成部、４０下りリンク自己相関行列生成部、６０プリコーディング行列変換部、４２上りリンク自己相関行列推定部、４４下りリンク自己相関行列変換部、４２自己相関行列推定部、４８到来角分布推定部、５０下りリンク自己相関行列演算部、５４到来方向（direction-of-arrival、DOA）推定部、５６下りリンク自己相関行列演算部、４２自己相関行列推定部、６２固有値分解部、６４固有ベクトル選択位相回転部、６６逆行列演算部、６８乗算正規化部（乗算部）、７０第１の固有値分解部、７２第２の固有値分解部、７４乗算正規化部（乗算部）

Claims

上りリンクの受信信号のベクトルの積を平均化して上りリンクの自己相関行列を求め、前記上りリンクの自己相関行列から下りリンクの自己相関行列を推定する下りリンク自己相関行列生成部と、
前記下りリンク自己相関行列生成部で推定された前記下りリンクの自己相関行列のうち、所望ユーザ端末の下りリンクの自己相関行列を固有値分解して固有ベクトルを求めるとともに、前記下りリンク自己相関行列生成部で推定された前記下りリンクの自己相関行列のうち、全ユーザ端末または他ユーザ端末の下りリンクの自己相関行列に基づいて、他ユーザ端末への送信信号に起因する同一チャネル干渉を抑圧するための干渉抑圧行列を生成し、前記固有ベクトルと前記干渉抑圧行列に基づいて、プリコーディング行列を生成するプリコーディング行列変換部と、
変調信号ベクトルに前記プリコーディング行列を乗算して、下りリンクの送信信号ベクトルを生成する送信信号生成部と
を備える無線基地局。
前記下りリンク自己相関行列生成部は、前記上りリンクの自己相関行列に線形変換を施すことにより、前記下りリンクの自己相関行列を推定する下りリンク自己相関行列変換部を備えることを特徴とする請求項１に記載の無線基地局。
前記下りリンク自己相関行列生成部は、前記上りリンクの自己相関行列から、各ユーザ端末について上りリンクのパスの到来角分布を推定する到来角分布推定部と、
前記到来角分布から前記下りリンクの自己相関行列を推定する自己相関行列演算部とを備えることを特徴とする請求項１に記載の無線基地局。
前記下りリンク自己相関行列生成部は、前記上りリンクの自己相関行列から、各ユーザ端末について上りリンクのパスの到来角分布を推定するとともに、各ユーザ端末について前記到来角分布を所定の種類の確率分布で近似させて、各ユーザ端末について前記確率分布の平均値と分散を推定する到来方向推定部と、
前記平均値と分散から前記下りリンクの自己相関行列を推定する自己相関行列演算部とを備えることを特徴とする請求項１に記載の無線基地局。
前記プリコーディング行列変換部は、前記下りリンク自己相関行列生成部で推定された前記下りリンクの自己相関行列のうち、所望ユーザ端末の下りリンクの自己相関行列を固有値分解して固有ベクトルを求める固有値分解部と、
前記固有ベクトルのうち、固有値が大きい方の１以上の固有ベクトルを選択し、選択された固有ベクトルの位相を回転させる固有ベクトル選択位相回転部と、
前記下りリンク自己相関行列生成部で推定された前記下りリンクの自己相関行列のうち、全ユーザ端末または他ユーザ端末の下りリンクの自己相関行列の逆行列である前記干渉抑圧行列を生成する逆行列演算部と、
前記固有ベクトル選択位相回転部で得られた位相が回転させられた固有ベクトルに、前記逆行列演算部で得られた前記干渉抑圧行列を乗算し、前記プリコーディング行列を生成する乗算部とを備えることを特徴とする請求項１から請求項４のいずれか１項に記載の無線基地局。
前記固有ベクトル選択位相回転部は、前記固有値分解部で求められた前記固有ベクトルのうち、前記所望ユーザ端末の受信アンテナの数と等しい数の固有ベクトルを選択することを特徴とする請求項５に記載の無線基地局。
前記プリコーディング行列変換部は、前記下りリンク自己相関行列生成部で推定された前記下りリンクの自己相関行列のうち、他ユーザ端末の下りリンクの自己相関行列を固有値分解して、固有値分解結果から固有値が小さい方の複数の第１の固有ベクトルを選択する第１の固有値分解部と、
前記下りリンク自己相関行列生成部で推定された前記下りリンクの自己相関行列のうち、所望ユーザ端末の下りリンクの自己相関行列を、前記第１の固有ベクトルを用いて変換し、変換された下りリンクの自己相関行列を固有値分解し、固有値分解結果から固有値が大きい方の１以上の第２の固有ベクトルを選択する第２の固有値分解部と、
前記第１の固有値分解部で選択された前記第１の固有ベクトルを含む行列である前記干渉抑圧行列と前記第２の固有値分解部で選択された前記第２の固有ベクトルを含む行列である所望信号増強行列とを乗算し、前記プリコーディング行列を生成する乗算部とを備えることを特徴とする請求項１から請求項４のいずれか１項に記載の無線基地局。
前記第２の固有値分解部は、前記下りリンク自己相関行列生成部で推定された前記下りリンクの自己相関行列のうち、所望ユーザ端末の下りリンクの自己相関行列を、前記第１の固有ベクトルを用いて変換し、変換された下りリンクの自己相関行列を固有値分解し、固有値分解結果から前記所望ユーザ端末の受信アンテナの数と等しい数の、固有値が大きい方の前記第２の固有ベクトルを選択することを特徴とする請求項７に記載の無線基地局。