"複数の送信アンテナを備え、少なくとも一つのアンテナを有する複数の端末装置と無線通信を行なう基地局装置であって、前記各端末装置との間の伝搬路情報に関連付けられた第一の制御情報に基づいて、前記伝搬路情報を示す第一の伝搬路行列を取得する取得部と、前記第一の伝搬路行列および前記第一の伝搬路行列に変換行列が乗算された第二の伝搬路行列に基づいて、前記各端末装置宛ての送信データに対して、それぞれプリコーディングを行なうプリコーディング部と、を備え、前記プリコーディング後の複数の送信データを同一無線リソースに空間多重して送信することを特徴とする基地局装置。modulo演算を含む非線形信号処理が可能な第一の端末装置に対する送信データと、線形信号処理のみが可能な第二の端末装置に対する送信データとを、同一無線リソースに空間多重して送信することを特徴とする請求項1記載の基地局装置。前記変換行列は、行列式が1または-1であり、すべての構成要素がガウス整数であることを特徴とする請求項1または請求項2記載の基地局装置。前記プリコーディング部は、前記変換行列、前記第一の伝搬路行列および前記第二の伝搬路行列に基づいて、線形フィルタ行列および摂動ベクトルを算出し、前記送信データに対して、前記算出した摂動ベクトルを加算し、前記線形フィルタを乗算し、前記変換行列は、少なくとも一つの単位行ベクトルを含み、前記摂動ベクトルは少なくとも一つの0を含むことを特徴とする請求項1から請求項3のいずれかに記載の基地局装置。前記プリコーディング部は、前記変換行列、前記第一の伝搬路行列および前記第二の伝搬路行列に基づいて、線形フィルタ行列および摂動ベクトルを算出し、前記送信データに対して、前記算出した摂動ベクトルを加算し、前記線形フィルタ行列を乗算し、前記摂動ベクトルと前記変換行列の逆行列との積で算出されるベクトルは、少なくとも一つの0を含むことを特徴とする請求項1から請求項3のいずれかに記載の基地局装置。前記プリコーディング部は、前記変換行列、前記第一の伝搬路行列および前記第二の伝搬路行列に基づいて、線形フィルタ行列および摂動ベクトルを算出し、前記摂動ベクトルと前記変換行列の逆行列との積で算出される第一の摂動ベクトルおよび第二の摂動ベクトルを加算することで得られるベクトルが、少なくとも一つの0を含み、前記送信データに対して前記第一の摂動ベクトルおよび前記第二の摂動ベクトルを加算し、前記線形フィルタ行列を乗算することを特徴とする請求項1から請求項3のいずれかに記載の基地局装置。前記第一の伝搬路行列および前記第二の伝搬路行列に基づいて、前記各端末装置における干渉電力を測定し、前記測定した干渉電力値に基づいて、前記各端末装置を第一の端末装置または第二の端末装置のいずれか一方に分類し、前記分類結果に関連付けた第二の制御情報を、前記各端末装置に通知することを特徴とする請求項1記載の基地局装置。前記第二の制御情報は、前記各端末装置との間で共有するMCSModulationandCodingSchemeセットに関連付けられた情報であることを特徴とする請求項7記載の基地局装置。請求項1から請求項8のいずれかに記載の基地局装置から無線信号を受信し、前記空間多重された送信データから自局宛てのデータを検出することを特徴とする移動局装置。複数の送信アンテナを備え、少なくとも一つのアンテナを有する複数の端末装置と無線通信を行なう基地局装置の制御プログラムであって、前記各端末装置との間の伝搬路情報に関連付けられた第一の制御情報に基づいて、前記伝搬路情報を示す第一の伝搬路行列を取得する処理と、前記第一の伝搬路行列および前記第一の伝搬路行列に変換行列が乗算された第二の伝搬路行列に基づいて、前記各端末装置宛ての送信データに対して、それぞれプリコーディングを行なう処理と、前記プリコーディング後の複数の送信データを同一無線リソースに空間多重して送信する処理と、の一連の処理を、コンピュータに実行させることを特徴とする制御プログラム。複数の送信アンテナを備えた基地局装置に実装されることにより、前記基地局装置に複数の機能を発揮させる集積回路であって、少なくとも一つのアンテナを有する複数の端末装置と無線通信を行なう機能と、前記各端末装置との間の伝搬路情報に関連付けられた第一の制御情報に基づいて、前記伝搬路情報を示す第一の伝搬路行列を取得する機能と、前記第一の伝搬路行列および前記第一の伝搬路行列に変換行列が乗算された第二の伝搬路行列に基づいて、前記各端末装置宛ての送信データに対して、それぞれプリコーディングを行なう機能と、前記プリコーディング後の複数の送信データを同一無線リソースに空間多重して送信する機能と、の一連の機能を、前記基地局装置に発揮させることを特徴とする集積回路。LRATHP技術は空間多重される全端末装置がmodulo演算を行なえることを前提とした技術であり、THPのように第2の端末装置宛の送信信号にmodulo演算を行なわないだけでは、二つの端末装置を混在して空間多重させることが出来ない。LRATHPにおいて第1の端末装置と第2の端末装置とを高効率に多重させられる技術については開示されていないのが実情である。本発明は、上述した課題に鑑み、modulo演算が行なえる第1の端末装置とmodulo演算が行なえない第2の端末装置とを、LRATHPに基づく非線形MUMIMOにおいても空間多重させることにより、周波数利用効率の改善に寄与できる基地局装置、移動局装置、制御プログラムおよび集積回路を提供することを目的とする。1上記の目的を達成するために、本発明は、以下のような手段を講じた。すなわち、本発明の基地局装置は、複数の送信アンテナを備え、少なくとも一つのアンテナを有する複数の端末装置と無線通信を行なう基地局装置であって、前記各端末装置との間の伝搬路情報に関連付けられた第一の制御情報に基づいて、前記伝搬路情報を示す第一の伝搬路行列を取得する取得部と、前記第一の伝搬路行列および前記第一の伝搬路行列に変換行列が乗算された第二の伝搬路行列に基づいて、前記各端末装置宛ての送信データに対して、それぞれプリコーディングを行なうプリコーディング部と、を備え、前記プリコーディング後の複数の送信データを同一無線リソースに空間多重して送信することを特徴としている。このように、基地局装置は、第一の伝搬路行列および第一の伝搬路行列に変換行列が乗算された第二の伝搬路行列に基づいて、各端末装置宛ての送信データに対して、それぞれプリコーディングを行なうので、LRA-THPに基づく非線形MU-MIMO空間多重伝送においても、modulo演算が行なえる端末装置とmodulo演算が行なえない端末装置とを空間多重することが可能となり、周波数利用効率の改善に寄与できる。2また、本発明の基地局装置において、modulo演算を含む非線形信号処理が可能な第一の端末装置に対する送信データと、線形信号処理のみが可能な第二の端末装置に対する送信データとを、同一無線リソースに空間多重して送信することを特徴としている。このように、基地局装置は、modulo演算を含む非線形信号処理が可能な第一の端末装置に対する送信データと、線形信号処理のみが可能な第二の端末装置に対する送信データとを、同一無線リソースに空間多重して送信するので、周波数利用効率の改善に寄与できる。3また、本発明の基地局装置において、前記変換行列は、行列式が1または-1であり、すべての構成要素がガウス整数であることを特徴としている。このように、変換行列が、行列式が1または-1であり、すべての構成要素がガウス整数であるので、基地局装置は、第一の伝搬路行列を、より直交性の高い行列に変換することができる。4また、本発明の基地局装置において、前記プリコーディング部は、前記変換行列、前記第一の伝搬路行列および前記第二の伝搬路行列に基づいて、線形フィルタ行列および摂動ベクトルを算出し、前記送信データに対して、前記算出した摂動ベクトルを加算し、前記線形フィルタを乗算し、前記変換行列は、少なくとも一つの単位行ベクトルを含み、前記摂動ベクトルは少なくとも一つの0を含むことを特徴としている。このように、変換行列が、少なくとも一つの単位行ベクトルを含み、摂動ベクトルは少なくとも一つの0を含むので、ある端末装置においてmodulo演算を必要としない伝送を実現することができる。5また、本発明の基地局装置において、前記プリコーディング部は、前記変換行列、前記第一の伝搬路行列および前記第二の伝搬路行列に基づいて、線形フィルタ行列および摂動ベクトルを算出し、前記送信データに対して、前記算出した摂動ベクトルを加算し、前記線形フィルタ行列を乗算し、前記摂動ベクトルと前記変換行列の逆行列との積で算出されるベクトルは、少なくとも一つの0を含むことを特徴としている。このように、摂動ベクトルと変換行列の逆行列との積で算出されるベクトルが、少なくとも一つの0を含むので、ある端末装置においてmodulo演算を必要としない伝送を実現することができる。6また、本発明の基地局装置において、前記プリコーディング部は、前記変換行列、前記第一の伝搬路行列および前記第二の伝搬路行列に基づいて、線形フィルタ行列および摂動ベクトルを算出し、前記摂動ベクトルと前記変換行列の逆行列との積で算出される第一の摂動ベクトルおよび第二の摂動ベクトルを加算することで得られるベクトルが、少なくとも一つの0を含み、前記送信データに対して前記第一の摂動ベクトルおよび前記第二の摂動ベクトルを加算し、前記線形フィルタ行列を乗算することを特徴としている。このように、摂動ベクトルと変換行列の逆行列との積で算出される第一の摂動ベクトルおよび第二の摂動ベクトルを加算することで得られるベクトルが、少なくとも一つの0を含むので、ある端末装置においてmodulo演算を必要としない伝送を実現することができる。7また、本発明の基地局装置において、前記第一の伝搬路行列および前記第二の伝搬路行列に基づいて、前記各端末装置における干渉電力を測定し、前記測定した干渉電力値に基づいて、前記各端末装置を第一の端末装置または第二の端末装置のいずれか一方に分類し、前記分類結果に関連付けた第二の制御情報を、前記各端末装置に通知することを特徴としている。このように、基地局装置は、測定した干渉電力値に基づいて、各端末装置を第一の端末装置または第二の端末装置のいずれか一方に分類するので、modulo演算を行なう端末装置と行なわない端末装置とを空間多重させることができる。8また、本発明の基地局装置において、前記第二の制御情報は、前記各端末装置との間で共有するMCSModulationandCodingSchemeセットに関連付けられた情報であることを特徴としている。このように、第二の制御情報が、各端末装置との間で共有するMCSModulationandCodingSchemeセットに関連付けられた情報であるので、基地局装置は端末装置に対して、新たに制御信号により通知せずに、第二の制御情報を通知することができる。9また、本発明の移動局装置は、上記1から8のいずれかに記載の基地局装置から無線信号を受信し、前記空間多重された送信データから自局宛てのデータを検出することを特徴としている。このように、上記1から8のいずれかに記載の基地局装置から無線信号を受信するので、基地局装置は、LRA-THPに基づく非線形MU-MIMO空間多重伝送においても、modulo演算が行なえる端末装置とmodulo演算が行なえない端末装置とを空間多重することが可能となり、周波数利用効率の改善に寄与できる。10また、本発明の制御プログラムは、複数の送信アンテナを備え、少なくとも一つのアンテナを有する複数の端末装置と無線通信を行なう基地局装置の制御プログラムであって、前記各端末装置との間の伝搬路情報に関連付けられた第一の制御情報に基づいて、前記伝搬路情報を示す第一の伝搬路行列を取得する処理と、前記第一の伝搬路行列および前記第一の伝搬路行列に変換行列が乗算された第二の伝搬路行列に基づいて、前記各端末装置宛ての送信データに対して、それぞれプリコーディングを行なう処理と、前記プリコーディング後の複数の送信データを同一無線リソースに空間多重して送信する処理と、の一連の処理を、コンピュータに実行させることを特徴としている。このように、基地局装置は、第一の伝搬路行列および第一の伝搬路行列に変換行列が乗算された第二の伝搬路行列に基づいて、各端末装置宛ての送信データに対して、それぞれプリコーディングを行なうので、LRA-THPに基づく非線形MU-MIMO空間多重伝送においても、modulo演算が行なえる端末装置とmodulo演算が行なえない端末装置とを空間多重することが可能となり、周波数利用効率の改善に寄与できる。11また、本発明の集積回路は、複数の送信アンテナを備えた基地局装置に実装されることにより、前記基地局装置に複数の機能を発揮させる集積回路であって、少なくとも一つのアンテナを有する複数の端末装置と無線通信を行なう機能と、前記各端末装置との間の伝搬路情報に関連付けられた第一の制御情報に基づいて、前記伝搬路情報を示す第一の伝搬路行列を取得する機能と、前記第一の伝搬路行列および前記第一の伝搬路行列に変換行列が乗算された第二の伝搬路行列に基づいて、前記各端末装置宛ての送信データに対して、それぞれプリコーディングを行なう機能と、前記プリコーディング後の複数の送信データを同一無線リソースに空間多重して送信する機能と、の一連の機能を、前記基地局装置に発揮させることを特徴としている。このように、集積回路が、第一の伝搬路行列および第一の伝搬路行列に変換行列が乗算された第二の伝搬路行列に基づいて、各端末装置宛ての送信データに対して、それぞれプリコーディングを行なうので、基地局装置は、LRA-THPに基づく非線形MU-MIMO空間多重伝送においても、modulo演算が行なえる端末装置とmodulo演算が行なえない端末装置とを空間多重することが可能となり、周波数利用効率の改善に寄与できる。本発明によれば、LRATHPに基づく非線形MUMIMO空間多重伝送においても、modulo演算が行なえる端末装置とmodulo演算が行なえない端末装置とを空間多重することが可能となるため、周波数利用効率の改善に寄与できる。以下、図面を参照して本発明の無線通信システムを適用した場合における実施形態について説明する。なお、本実施形態において説明した事項は、発明を理解するための一態様であり、実施形態に限定して発明の内容が解釈されるものではない。[1.基本技術]本発明の実施形態の説明をする前に、背景技術であるLRATHPを用いる非線形MUMIMO空間多重伝送以降、LRATHPMUMIMO、もしくは単にLRATHPと記述するについて説明する。Nt本の送信アンテナを有する基地局装置送信装置とも呼ぶに対して、1本の受信アンテナを有する端末装置移動局装置、もしくは受信装置とも呼ぶがU個接続している通信を対象とし、Nt=Uであるものとする。なお、全ての端末装置はmodulo演算が行なえる第1の端末装置であるものとする。はじめに基地局装置の送信アンテナ109nn=1,,Ntと端末装置300u間u=1,,Uの複素チャネル利得をhu,nとしたとき、伝搬路行列Hをと定義する。LRATHP方式では、基地局装置は、伝搬路行列Hを予め知っている必要がある。ここで、伝搬路行列に関連付けられた情報は、各端末装置から理想的に通知されており、基地局装置は、伝搬路行列Hを理想的に知っているものとしている。[1.1基地局装置]図16は、基地局装置100の構成を示すブロック図である。図16を用いて基地局装置100における信号処理について説明する。各端末装置宛の送信データはチャネル符号化部101において、チャネル符号化が行なわれたのち、データ変調部103において、QPSK、16QAM等にデータ変調される。チャネル符号化における符号化率や、データ変調方式は、常に一定のものを用いても良いし、各端末装置の受信品質に応じて適応的に変更するようにしても良い。適応的に変更する場合、基地局装置100と端末装置は、予め受信品質毎に設定される符号化率とデータ変調方式がテーブル化されているModulationandCodingSchemeMCSセットを共有しておき、端末装置から通知される受信品質に関連付けられた制御情報に基づき、符号化率等を決定すれば良い。以下の説明では、符号化率およびデータ変調方式は適切に決定されているものとして説明する。データ変調部103からの出力は参照信号多重部105に入力され、各端末装置において伝搬路推定を行なうための既知参照信号系列が参照信号多重部105において多重される。各端末装置宛の参照信号については、受信した端末装置において分離可能なように、それぞれが直交するように多重されるものとする。以下の説明では、参照信号は任意の無線リソースに理想的に配置されたものとし、端末装置では上記既知参照信号系列により、理想的に伝搬路推定が行なわれるものとする。参照信号多重部105の出力は、Precoding部107に入力される。図17は、Precoding部107の詳細な構成について説明するためのブロック図である。Precoding部107には各端末装置宛の送信シンボル{du;u=1,,U}を要素とする送信シンボルベクトルd=[d1,,dU]Tと伝搬路行列Hが入力される。はじめに、入力された伝搬路行列Hが格子基底縮小部201に入力され、格子基底縮小技術が適用される。本発明で対象とする格子基底縮小技術は、行列変換技術の一つであり、入力された行列Hに対して、unimodular行列Tを乗算し、より直交性の高い行列Gに変換する技術である。なお、unimodular行列とは、構成要素が全てガウス整数であり、行列式の絶対値が1となる行列を指す。格子基底縮小が行なわれた伝搬路行列をGとすると、G=THと表現される。従来のLRATHPに適したunimodulor行列の算出方法についてはLenstraらが提案したLLLアルゴリズムがある。以下では、任意のアルゴリズムにより、unimodular行列が算出されたものとして説明を行なう。格子基底縮小部201からは行列Tが出力される。伝搬路情報Hと格子基底縮小部201から出力されたTは線形フィルタ生成部203に入力され、線形フィルタWが生成される。線形フィルタWは、格子基底縮小された伝搬路行列Gに対して、GWが単位下三角行列となるような線形フィルタである。Gはユニタリ行列Qと下三角行列Lを用いて、G=LQのようにLQ分解することができるから、GWを単位下三角行列とする線形フィルタWはW=QH{diagL}-1となる。ここで、diag{A}は行列Aの対角成分のみを抽出した対角行列、A-1は行列Aの逆行列を表す。なお、GWを単位下三角行列とする線形フィルタは残留干渉を0とするZF規範に基づく線形フィルタである。ZF規範ではなく、送信信号と受信信号の平均二乗誤差を最小とするMMSE規範に基づく線形フィルタを用いても良い。この場合、H’=[HαI]で定義される拡張伝搬路行列を伝搬路行列と見なしてZF規範の線形フィルタを計算すれば良い。なお、Iは単位行列である。またαは残留干渉を制御する係数であり、所望の伝送品質に応じて任意に設定されるが、例えば端末装置で観測される雑音電力の標準偏差σに設定すれば良い。送信シンボルベクトルdにWを乗算することにより、送信信号ベクトルs=[s1,...,sNt]Tはs=Wdと計算できる。仮に送信信号ベクトルsを基地局装置100から送信したとする。端末装置300uで受信される受信信号を{ru;u=1,...,U}としたとき、各端末装置の受信信号を要素とする受信信号ベクトルr=[r1,...,rU]Tは式21で与えられる。ここで、η=[η1,...,ηU]Tは各端末装置で加わる雑音を表すが、以下では、簡単のため、雑音項は無視して説明する。式21は式22のように展開できる。ここで、ai,jは単位下三角行列L{diagL}-1の第i行j列成分を表す。式22より、T-1を無視して考えた場合、第1の端末装置は自身の送信シンボルのみを受信できるが、第2の端末装置の受信信号には、第1の端末装置の送信シンボルが干渉を与えていることが分かる。つまり、端末装置300uの受信信号には第1~u1端末装置の送信シンボルが干渉として含まれている。そこで、Precoding部107に入力された送信シンボルベクトルdは事前干渉抑圧部205に入力されることで、この各端末装置で観測される干渉成分が、予め減算されることになる。例えば、第1の端末装置への送信信号として、d2でなくx2=d2-a2,1d1を送信するものとする。式22のd2の部分にx2を代入すれば、第2の端末装置はd1からの干渉を受けずに、自身の送信シンボルd2のみを受信することが可能であることが分かる。以下、同様に端末装置300u宛の送信信号から、予め第1~u1端末装置宛の送信シンボルを減算してから送信を行なうことで、全端末装置が干渉を受けずに通信を行なうことが可能となる。端末装置300u宛の事前干渉抑圧が為された信号のことを送信符号xuと呼ぶ事とする。xuは式31で与えられることになる。このとき、送信符号ベクトルx=[x1,,xNt]Tはx=GW-1dと表現できる。xを式22のdに代入することで受信信号は式23で与えられる。しかし、受信信号にはT-1が残るため、さらに送信シンボルベクトルdをd^=[d^1,,d^Nt]T=Tdに置き換える。以上の信号処理により全端末装置は干渉を受けずに通信を行なうことができる。しかし、unimodular行列Tは必ずしも直交行列ではないから、Tが乗算された送信シンボルベクトルd^の大きさは元々の送信シンボルベクトルdより大きくなってしまう場合がある。また、式31で示される干渉抑圧により、送信符号ベクトルxの大きさが、送信シンボルベクトルd^より大きくなる場合もある。そこで、LRATHPでは、各端末装置宛の信号から干渉を取り除く際に、modulo演算と呼ばれる信号処理を行なう。modulo演算ModMxは、ある入力xに対して、その出力が-Mより大きく、かつM以下に収まるようにするものである。入力xが複素数であった場合、その出力の実部と虚部がそれぞれ-Mより大きく、かつM以下に収まる。ここでMをmodulo幅と呼ぶこととする。実際に、送信符号x2=d^2-a2,1d^1にmodulo演算を適用した場合、その出力は式32で与えられる。ここで、floorxは実数xを超えない最大の整数を返す関数であり、床関数とも呼ばれる。また、RecおよびImcはそれぞれ、複素数cの実数および虚数を返す関数である。また、zt,2のことを摂動項と呼ぶこととする。modulo演算を行なう場合、送信符号xuは式33で与えられることになる。以上説明してきたmodulo演算をx1=d^1から順番にxUまでの全ての送信符号に施す。送信符号ベクトルxは式4で与えられる。ここで、zt=[zt,1,...,zt,Nt]Tであり、以降ztのことを摂動ベクトルと呼ぶ。数式上は行列演算で表現されているが、実際には、式31の干渉抑圧と式32のmodulo演算を各送信符号に対して逐次的に行なうことで送信符号が生成される。以上説明してきたように、事前干渉抑圧部205には、送信シンボルベクトルd、伝搬路情報H、unimodular行列Tおよび線形フィルタWが入力され、事前干渉抑圧部205出力として送信符号ベクトルxが出力されることとなる。事前干渉抑圧部205出力xは線形フィルタ乗算部207に入力される。線形フィルタ乗算部207では、入力されたxに対して、線形フィルタW=QH{diagL}-1および電力正規化係数βが乗算され、送信信号ベクトルs=βWxが算出される。ここで、であり、Psは摂動ベクトルが付与された送信シンボルベクトルTd2Mztの共分散行列を表す。trAは行列Aに対するトレース演算を表す。Precoding部107からは線形フィルタ乗算部207出力sが最終的に出力される。その後、基地局装置100では、Precoding部107の出力sが対応する各送信アンテナ109の無線送信部111に入力される。無線送信部111において、ベースバンド帯の送信信号が無線周波数RF帯の送信信号に変換される。無線送信部111の出力信号は、各送信アンテナ109よりそれぞれ送信される。[1.2端末装置]図18は、端末装置300の構成を示すブロック図である。図18を用いて端末装置300における信号処理について説明する。端末装置300では、アンテナ301で受信された信号が無線受信部303に入力され、ベースバンド帯の信号に変換される。ベースバンド帯に変換された信号は、参照信号分離部305に入力される。参照信号分離部305では、受信信号はデータ系列と既知参照信号系列とに分離され、データ系列は伝搬路補償部307に入力され、既知参照信号系列は伝搬路推定部309に入力される。伝搬路推定部309では、入力された既知参照信号系列を用いて伝搬路推定が行なわれる。各端末装置300宛の既知参照信号系列はそれぞれ直交するように基地局装置100より送信されているから、端末装置300では、伝搬路情報を推定することができる。端末装置300uにおいて推定される伝搬路情報は基地局装置100のPrecoding部107にて用いられる伝搬路行列Hの第u行成分に該当する。推定された伝搬路情報はそれぞれ伝搬路補償部307およびフィードバック情報生成部311に入力される。フィードバック情報生成部311では、各端末装置300がフィードバックする伝搬路情報形式に応じて、基地局装置100にフィードバックする情報を生成する。本発明においては、伝搬路情報形式については何かに限定されるものではない。例えば、推定された伝搬路情報について、有限ビット数にて量子化を行ない、その量子化情報をフィードバックする方法が考えられる。フィードバック情報生成部311で生成された情報は、無線送信部313に入力され、基地局装置100に向けて通知される。基地局装置100では制御情報取得部115により、端末装置300より送信された情報を取得し、その情報がCSI取得部117に入力される。そして通知されてきた情報と、その伝搬路情報形式に基づいて、伝搬路情報Hが算出され、HはPrecoding部107に入力されることになる。一方、受信データ系列については伝搬路補償部307に入力され、伝搬路補償が行なわれる。ベースバンド帯に変換された受信信号は式6で与えられる。ここで、z^t=T-1ztであり,z^tの構成要素は全てガウス整数となる。式6より、各端末装置300の受信信号に対して、他の端末装置300宛の送信シンボルは干渉を与えていないことが分かる。端末装置300uの伝搬路補償部307では、入力された受信信号から電力正規化項βを除算する。そのため、端末装置300はβの値を把握する必要があるが、βは基地局装置100より送信される既知参照信号より推定することが可能である。例えば、基地局装置100が伝搬路情報推定のために送信する既知参照信号の一部に対して、受信データ系列と同様のPrecoding処理を行なうことで、βを推定できる。以下では、βは理想的に推定できるものとして説明していく。電力正規化項が除算されたのち、基地局装置100のPrecoding部107で用いられたmodulo演算と同じmodulo幅を有するmodulo演算が適用される。端末装置300uのModulo出力は式7で与えられる。ここで、zr,uは端末装置300uの受信機で適用されるmodulo演算により生成される摂動項であり、z^t,uzr,u=0となる。よって、最終的に伝搬路補償部307から出力されるのはduということになる。その後、伝搬路補償部307の出力はデータ復調部315およびチャネル復号部317に入力され、データ復調、チャネル復号がそれぞれ適用されたのち、各端末装置300宛の送信データが検出される。以上が、LRATHPMUMIMOの送受信処理の説明となる。空間多重される端末が、全てmodulo演算が行なえる第1の端末装置300であった場合、上述してきた空間多重方法により線形MUMIMOやTHPMUMIMOと比較して、優れた伝送特性を実現できる。しかし、実際の無線通信システムにおいては、第1の端末装置300と、modulo演算が行なえない第2の端末装置300とが混在する場合が考えられる。ところで、LRAを行なわないTHPMUMIMOは上記説明の中で、unimodular行列Tを単位行列Iに置き換えることで実現できる。この場合、式32で行なわれるmodulo演算の適用の有無については、送信電力の強調を許容する限りにおいては、端末装置300毎に変えることが可能である。ここで、端末装置300iがmodulo演算を行なえない第2の端末装置300であるものとする。この場合、端末装置300i宛の送信信号に対してmodulo演算を行なわない、すなわちzt,i=0とすることで、端末装置300iが第2の端末装置300であったとしてもMUMIMO伝送を行なうことが可能である。なお、全端末装置300に対してmodulo演算を行なわない、すなわちzt=0U×10m×nは全ての構成要素が0であるm行n列の零行列とすれば、THPMUMIMOは線形MUMIMOと等価となる。しかし、LRATHPでは、第2の端末装置300宛の送信信号に対してmodulo演算を行なわないだけでは、第2の端末装置300を空間多重することが出来ない。それは、式6において現れるz^t=T-1ztのためである。Tは伝搬路行列に応じて決定される行列であるから、THPのようにzt,i=0としても、必ずしもz^t,i=0とはならない。よって、従来のLRATHPで第1の端末装置300と第2の端末装置300とを空間多重するためには、基地局装置100においてすべての端末装置300宛の送信信号に対するmodulo演算を行なわない必要がある。しかし、それは線形MUMIMOと等価であり、周波数利用効率の改善は望めないことを意味している。そこで本実施形態では、LRATHPにおいて、第1の端末装置300と第2の端末装置300とを空間多重する方法を開示する。[2.第1の実施形態]第1の実施形態においては、LRATHPにおいて、modulo演算を行なう第1の端末装置300と、modulo演算を行なわない第2の端末装置300とを空間多重させる場合を対象とする。Nt本の送信アンテナ109を有する基地局装置送信装置とも呼ぶに対して、1本のアンテナ301を有する端末装置300移動局装置もしくは受信装置とも呼ぶがU個接続している通信を対象とし、端末装置300のうち、端末装置300iが第2の端末装置300であるものとする。残りの端末装置300は全て第1の端末装置300である。なお、各端末装置300が有するアンテナ301数や、接続されている端末装置300に含まれる第1の端末装置300と第2の端末装置300の割合についてはこれに限ったものではなく、複数アンテナ301を有する端末装置300が存在する場合や、複数の第2の端末装置300が存在する場合にも、本発明は適用可能である。また、以下の説明では、基地局装置は接続している端末装置300が第1の端末装置300か第2の端末装置300かを把握しているものとしている。図1は、本発明の第1の実施形態に係るシステム全体の信号処理についてその概要を説明するフローチャートである。はじめに、基地局装置100は、送信データを生成するとともにステップS101、接続している第1の端末装置300および第2の端末装置300に対して、既知の参照信号系列を送信するステップS102。次いで、各端末装置300は受信された既知参照信号系列に基づき、伝搬路推定を行なうステップS103。各端末装置300は推定された伝搬路情報を基地局装置100に通知するが、通知方法については、特定の方法に限定されるものではなく、基地局装置100が伝搬路情報を把握できるのであればどのような方法でも構わないステップS104。例えば、推定された伝搬路情報を有限ビット長に量子化した情報を通知すれば良い。また、基地局装置100と端末装置300間で予め複数の線形フィルタが記載されたコードブックを共有しておき、伝搬路情報と最も相関が高い線形フィルタのインデックスを通知しても良い。基地局装置100は、各端末装置300より伝搬路情報を取得し、その情報に基づき、各端末装置300宛の送信データに対して、後で述べる送信符号化プリコーディングを施すステップS105。プリコーディング後のデータを同一無線リソースに空間多重して、接続されている複数の端末装置300宛に送信するステップS106。複数端末装置300宛のデータが空間多重された信号を受信した端末装置300では、それぞれ前述の伝搬路情報等に基づき、自局の所望データを検出するがステップS108、第1の端末装置300では、modulo演算を含む非線形演算処理を行なう必要があるのに対してステップS107、第2の端末装置300ではmodulo演算は必要としない。[2.1基地局装置100]図2は、本発明の第1の実施形態に係る基地局装置100の構成を示すブロック図である。図2に示すように、基地局装置100は、チャネル符号化部101と、データ変調部103と、参照信号多重部105と、Precoding部プリコーディング部401と、無線送信部111と、無線受信部113と、制御情報取得部115と、CSI取得部取得部117と、送信アンテナ109とを含んで構成されている。各端末装置300宛の送信データ系列はチャネル符号化部101において、チャネル符号化が行なわれたのち、データ変調部103において、QPSK、16QAM等にデータ変調される。データ変調部103からの出力は参照信号多重部105に入力され、各端末装置300において伝搬路推定を行なうための既知参照信号系列が参照信号多重部105において多重される。各端末装置300宛の参照信号については、受信した端末装置300において分離可能なように、それぞれが直交するように多重されるものとする。以下の説明では、参照信号は任意の無線リソースに理想的に配置されたものとし、端末装置300では上記既知参照信号系列により、理想的に伝搬路推定が行なわれるものとする。参照信号多重部105の出力は、Precoding部401に入力される。なお、各端末装置300より通知される伝搬路情報についても、CSI取得部117よりPrecoding部401に入力されることとなる。図3は、本発明の第1の実施形態に係るPrecoding部401の装置構成を示すブロック図である。図3に示すように、Precoding部401は、線形フィルタ生成部203と、格子基底縮小部501と、事前干渉抑圧部と、線形フィルタ乗算部207とを含んで構成されている。Precoding部401で行なわれる信号処理は、図17に示したPrecoding部401で行なわれる信号処理とほぼ同じであるが、格子基底縮小部501および事前干渉抑圧部503に行なわれる信号処理が異なっている。はじめに、格子基底縮小部501で行なわれる信号処理について説明する。格子基底縮小部501では入力された伝搬路行列Hに基づいて、Hの直交性を高めるunimodular行列Tが算出される。格子基底縮小技術に適したTの算出アルゴリズムとして、LLLアルゴリズムが良く知られている。図4は、LLLアルゴリズムのアルゴリズムを示す表である。算出されたTに基づいて、事前干渉抑圧部503および線形フィルタ生成部203における信号処理が行なわれる。なお、unimodular行列Tの構成要素はガウス整数であることが望ましいが、伝送特性の劣化を許容できる場合、実数に属する数字が構成要素となっても構わない。しかし、従来のLLLアルゴリズムにおいて算出されるunimodular行列Tでは、空間多重される全端末装置300がmodulo演算を行なえる第1の端末装置300である必要がある。そこで、本実施形態においては、従来のLLLアルゴリズムを改良することで、第2の端末装置300の空間多重を可能とする。LRATHPにて送信された信号を、端末装置300iがmodulo演算を行なわずに復調するためには、式6においてz^t,i=0である必要がある。z^t,iはベクトルT-1ztの第i成分であるから、行列T-1の第i行とztとの積で表現される。ここで、第k成分が1であり、他の構成要素が全て0であるような要素数Uの行ベクトルを第k単位行ベクトルeU,kと定義する。ここで、行列の持つ性質により、U行U列の行列Tの第k行がeU,iであるとき、行列Tの逆行列T-1は第i行にeU,kを有する行列となる。このとき、zt,k=0とすれば、端末装置300iではmodulo演算を必要としない伝送とできる。つまり、端末装置300iでmodulo演算を必要としない伝送を実現するためには、行列Tのいずれかの行に第i単位行ベクトルeU,iが存在すれば良いことになる。LLLアルゴリズムは入力された行列Hに対して、unimodular行列Tを算出しG=HTという直交性の高い行列を算出する。このアルゴリズムをLRATHPに用いる場合、入力される値は伝搬路行列Hのエルミート転置行列HHであり、このとき算出されるunimodular行列をTHと考えると、行列GH=HHTHが算出され、直交性の高い行列GはG=THとして出力される。なお、入力される値は伝搬路行列Hの転置行列HT、もしくは逆行列H-1のいずれかでも良い。実際のLLLアルゴリズムで出力されるのはTHである一方で、本実施形態の目的はTのいずれかの行にeU,iを存在させることにある。よって、LLLアルゴリズムで算出される行列TH図4におけるTのいずれかの列ベクトルがeU,iHであれば良い。LLLアルゴリズムは、図4から分かるように、4行目から10行目までの1準直交化処理、11行目の2条件比較、12行目から16行目までの3入れ替え作業の3つの処理の繰り返しで構成されている。ここで、unimodular行列Tが直接算出されているのは、1の準直交化処理の部分である。unimodular行列の初期値は単位行列Iであるから、第i列を対象とした準直交化処理を行なわなければ、最終的に算出されるTにはeU,iHが存在することになる。なお、2の条件比較および3の入れ替え作業により、第i列が違う列例えば第j列と交換される場合がある。この場合、新しく第j列がeU,iHとなる。つまり、繰り返し処理の中で、元々第i列成分であった列ベクトルがどの列に移動したのかを常に把握しておく必要がある。以下の説明では、アルゴリズムが終了した時点でeU,iHは第v列に存在するものとする。図5は、本発明の第1の実施形態に係るLLLアルゴリズムのアルゴリズムを示す表である。入力は伝搬路行列Hに加えて、第2の端末装置300の端末装置番号を示すiであり、出力は、unimodular行列Tと最終的にeU,iHが存在する列番号を示すvである。なお、第2の端末装置300数が複数で有る場合は準直交化を行なわない列が複数となることになる。その場合、3の入れ替え処理時に覚えておくべき列番号も複数となる。本実施形態では、格子基底縮小部501からは、図5で表されるアルゴリズムにより、unimodular行列Tとvという二つの情報が出力される。図3に戻り、線形フィルタ生成部203にはunimodular行列Tと伝搬路行列Hが入力され、線形フィルタWが生成される。線形フィルタWは、格子基底縮小された伝搬路行列Gに対して、GWが単位下三角行列となるような線形フィルタである。Gはユニタリ行列Qと下三角行列Lを用いて、G=LQのようにLQ分解することができるから、GWを単位下三角行列とする線形フィルタWはW=QH{diagL}-1となる。生成された線形フィルタは事前干渉抑圧部503および線形フィルタ乗算部207に入力される。事前干渉抑圧部503には、送信信号ベクトルd、伝搬路情報H、unimodular行列T、eU,iHが存在する列番号を示すv、および線形フィルタWが入力され、事前干渉抑圧部503出力としてxが出力されることとなる。xの算出方法は基本技術で述べたように、式4を実現するものとなるが、このとき、d^vに対してはmodulo演算を行なわない。つまり、式4において、zt=[zt,1,...,zt,v-1,0,zt,v1,...,zt,Nt]Tとなることになる。具体的には、d^vに対する干渉抑圧は式31に基づき、それ以外の送信シンボルに対しては式33に基づいて行なわれることになる。第2の端末装置300が複数である場合、vに該当する数字も複数となるが、この場合も、vに該当する全ての摂動ベクトルを0とすれば良い。事前干渉抑圧部503より出力されるxはその後、線形フィルタ乗算部207に入力され、線形フィルタWおよび電力正規化係数βが乗算され、送信信号ベクトルsが出力され、これがPrecoding部401出力となる。図6は、本発明の第1の実施形態に係るPrecoding部401における信号処理を示すフローチャートである。まず、Precoding部401は、送信シンボルベクトルdおよび伝搬路情報Hを取得するステップS201。次に、格子基底縮小部501では入力された伝搬路行列Hに基づいて、Hの直交性を高めるunimodular行列Tが算出されるステップS202。伝搬路情報Hと格子基底縮小部501から出力されたTは線形フィルタ生成部203に入力され、線形フィルタWが生成されるステップS203。続いて、Precoding部401に入力された送信シンボルベクトルdは事前干渉抑圧部503に入力されることで、各端末装置300で観測される干渉成分が、予め減算されるステップS204。この時、第1の端末装置300宛てには式33に基づいて行なわれ、第2の端末装置300宛てには式31に基づいて行なわれる。次に、線形フィルタ乗算部207では、入力されたxに対して、線形フィルタWが乗算され、送信信号ベクトルsが算出されるステップS205。図2に戻り、基地局装置100では、Precoding部401の出力sが対応する各送信アンテナ109の無線送信部111に入力される。無線送信部111において、ベースバンド帯の送信信号が無線周波数RF帯の送信信号に変換される。無線送信部111の出力信号は、各送信アンテナ109よりそれぞれ送信される。[2.2端末装置300]図7は、本発明の第1の実施形態に係る第1の端末装置300の構成を示すブロック図である。図8は、本発明の第1の実施形態に係る第2の端末装置300の構成を示すブロック図である。図7および図8に示すように、端末装置300はアンテナ301と、無線受信部303と、参照信号分離部305と、伝搬路推定部309と、フィードバック情報生成部311と、無線送信部313と、伝搬路補償部601と、データ復調部315と、チャネル復号部317とを含んで構成されている。本実施形態においては、端末装置300iが第2の端末装置300に含まれ、残りの端末装置300については全て第1の端末装置300に含まれるものとする。第1の端末装置300と第2の端末装置300の装置構成は同じであるが、それぞれの伝搬路補償部601において行なわれる信号処理が異なる。各端末装置300においては、アンテナ301で受信された信号が、無線受信部303に入力され、無線受信部303において、ベースバンド帯の信号に変換される。ベースバンド帯に変換された信号は、参照信号分離部305に入力される。参照信号分離部305では、受信信号はデータ系列と既知参照信号系列とに分離され、データ系列は伝搬路補償部601に入力され、既知参照信号系列は伝搬路推定部309に入力される。伝搬路推定部309では、入力された既知参照信号系列を用いて伝搬路推定が行なわれる。各端末装置300宛の既知参照信号系列はそれぞれ直交するように基地局装置100より送信されているから、端末装置300では、伝搬路情報を推定することができる。端末装置300uにおいて推定される伝搬路情報は基地局装置100のPrecoding部401にて用いられる伝搬路行列Hの第u行成分に該当する。推定された伝搬路情報はそれぞれ伝搬路補償部601およびフィードバック情報生成部311に入力される。フィードバック情報生成部311では、各端末装置300がフィードバックする伝搬路情報形式に応じて、基地局装置100にフィードバックする情報を生成する。本発明においては、伝搬路情報形式については何かに限定されるものではない。例えば、推定された伝搬路情報を有限ビット長に量子化した情報を通知すれば良い。また、基地局装置100と端末装置300間で予め複数の線形フィルタが記載されたコードブックを共有しておき、伝搬路情報と最も相関が高い線形フィルタのインデックスを通知しても良い。フィードバック情報生成部311で生成された情報は、無線送信部313に入力され、基地局装置100に向けて通知される。基地局装置100ではCSI取得部117において、通知されてきた情報と、その伝搬路情報形式に基づいて、伝搬路情報Hを取得し、取得された伝搬路情報がPrecoding部401に入力されることになる。一方、受信データ系列については伝搬路補償部601に入力され、伝搬路補償が行なわれる。ベースバンド帯に変換された受信信号は式6で与えられるものと同じ形式である。端末装置300uの伝搬路補償部601では、入力された受信信号から電力正規化項βを除算する。そのため、端末装置300はβの値を把握する必要があるが、βは基地局装置100より送信される既知参照信号より推定することが可能である。例えば、基地局装置100が伝搬路情報推定のために送信する既知参照信号の一部に対して、受信データ系列と同様のPrecoding処理を行なうことで、βを推定できる。以下では、βは理想的に推定できるものとして説明していく。なお、端末装置300から基地局装置100への伝搬路情報のフィードバック精度が低い場合や、伝搬路の時間および周波数の選択性が激しい場合、端末装置300-uの受信信号において、希望信号duの振幅および位相が変動する。この振幅および位相の変動は電力正規化項βと同様に推定可能であるから、βを除算する際に、同時に振幅および位相の変動分の補償を行えば良い。電力正規化項が除算されたのち、基地局装置100のPrecoding部401で用いられたmodulo演算と同じmodulo幅を有するmodulo演算が適用される。しかし、本実施形態の基地局装置100のPrecoding部401において算出されるunimodulor行列Tと、modulo演算の適用方法により、z^tのうち、端末装置300iに該当する成分、つまり、z^t,iについては0となる。そのため、第2の端末装置300に含まれる端末装置300iの伝搬路補償部701においてはmodulo演算を必要とせず、電力正規化係数βの除算のみが行なわれる。一方で、第1の端末装置300に含まれる端末装置300については、基本技術と同様に、電力正規化係数βの除算後にmodulo演算が行なわれる事となる。その後、伝搬路補償部601、701の出力はデータ復調部315およびチャネル復号部317に入力され、データ復調、チャネル復号がそれぞれ適用されたのち、各端末装置300の送信データが検出される。以上の処理により、LRATHPMUMIMOにおいても、第1の端末装置300と第2の端末装置300とを高効率に空間多重させることが可能となる。ところで、本実施形態によるMUMIMOでは、第1の端末装置300と第2の端末装置300とを空間多重させているが、このことは、modulo演算を行なう端末装置300と行なわない端末装置300とを空間多重できることを意味している。modulo演算はIUIの増加に伴う送信電力の増大を抑圧するため、伝送特性の改善に寄与できる。しかし、IUIが小さい場合には、modulo演算に起因するmodulo損失やprecoding損失により伝送特性がmodulo演算を行なわない場合より劣化する場合があり、特にQPSK変調時に劣化しやすい。そのために、空間多重されている端末装置300が全て第1の端末装置300であった場合でも、IUIの状況に応じて、modulo演算の適用の有無を適応的に切り替える伝送方法が取られる場合がある。このときに、観測されるIUIが少ない端末を第2の端末装置300と見なすことにより、modulo演算を行なわないQPSK変調信号と、modulo演算を行なうQPSK変調信号とを本実施形態の方法により空間多重させることが可能となる。modulo演算の切り替えは、実際のプリコーディング時に測定できる干渉成分例えば式31の第2項の電力を測定し、ある閾値を超えたらmodulo演算を行ない、超えなければmodulo演算を行なわないように制御すれば良い。閾値については、与えられた符号化率等に基づき、予めBER測定などを行ない決定すれば良い。なお、このとき端末装置300は受信された信号に対して、modulo演算の有無を知る必要がある。基地局装置100が端末装置300に対して、新たに制御信号により通知しても良いが、基地局装置100と端末装置300とが共有するMCSに関連付けて通知する方法が考えられる。例えば、QPSK変調であれば、常にmodulo演算を行なうように制御しても良いし、modulo演算を行なうQPSKと行なわないQPSKをそれぞれMCSの項目に含めるようなシステムとしても良い。本実施形態においては、伝送方式もしくはアクセス方式については制限を与えていない。例えば、LTEの下りリンク伝送に採用されている直交周波数分割多重アクセスOFDMA方式に適用することが可能である。この場合は、サブキャリア毎に本実施形態を適用すれば良く、また複数サブキャリアを一纏めとしたリソースブロック毎に本実施形態を適用しても良い。同様に、シングルキャリアベースのアクセス方式例えばシングルキャリア周波数分割多重アクセスSCFDMA方式などに適用することも可能であり、周波数成分毎に適用しても良いし、送信電力の強調を回避するために、全周波数帯域に渡って同一のプリコーディングを行なうようにしても良い。第1の実施形態によれば、第1の端末装置300と、第2の端末装置300とが混在している無線通信システムにおいて、LRATHPによる高効率な空間多重が可能となるため、システム全体の周波数利用効率の改善に寄与できる。[3.第2の実施形態]第1の実施形態においては、格子基底縮小部501で算出されるunimodular行列Tと、事前干渉抑圧部における信号処理を従来技術と異なるものとすることで、第1の端末装置300と第2の端末装置300とをLRATHPにより空間多重させる方法を明らかにした。しかし、第1の実施形態の方法では、LLLアルゴリズムの収束速度が、従来の方法と比較して、著しく劣化してしまう場合がある。そこで、第2の実施形態では、格子基底縮小部に入力される伝搬路行列に予めオーダリングを施すことにより、第1の端末装置300と第2の端末装置300とをLRATHPにより空間多重させながら、LLLアルゴリズムの収束性の低下を回避する方法を対象とする。第1の実施形態と同様に、Nt本の送信アンテナ109を有する基地局装置100送信装置とも呼ぶに対して、1本のアンテナ301を有する端末装置300受信装置とも呼ぶがU個接続している通信を対象とし、端末装置300のうち、端末装置300iが第2の端末装置300であるものとする。なお、各端末装置300が有するアンテナ301数や、端末装置300に含まれる第1の端末装置300と第2の端末装置300の割合については、これに限ったものではなく、複数アンテナ301を有する端末装置300が存在したり、複数の第2の端末装置300が存在する場合にも、本発明は適用可能である。[3.1基地局装置100]第2の実施形態に係る基地局装置100構成は図2とほぼ同じであり、異なるのはPrecoding部における信号処理のみであるため、以下では、第2の実施形態に係るPrecoding部について説明する。図9は、本発明の第2の実施形態に係るPrecoding部801の構成を示すブロック図である。図9に示すように、Precoding部801は、線形フィルタ生成部203と、格子基底縮小部803と、事前干渉抑圧部805と、線形フィルタ乗算部207に加えて、オーダリング部807とを含んで構成されている。Precoding部801では、入力される伝搬路行列Hが初めにオーダリング部807に入力されてオーダリング処理が施される。ここでいうオーダリングとは、行列の列同士もしくは行同士を入れ替えることを意味する。オーダリング済みの伝搬路行列をHpとすると、Hpは順列行列Πと伝搬路行列Hを用いて、Hp=ΠHと表現できる。本実施形態においては、第2の端末装置300に該当する伝搬路情報が、オーダリング済み伝搬路行列Hpの一番上になるように制御する。例えば、本実施形態が対象としているように、端末装置300iが第2の端末装置300であった場合、順列行列Πとして、eU,iが第1行に配置されるような行列を用いれば良い。Πの他の行については、自由に入れ替えて良いから、順列行列としては全部でU-1通りの選択肢が考えられることとなる。いずれを用いても本実施形態の目的は達せられるが、一般にLLLアルゴリズムでは、入力行列の列ベクトルの大きさは、左から昇順に並ぶ方がアルゴリズムの収束特性を向上させる面から望ましい。そのような順列行列を求める方法として、例えば、初めに伝搬路行列Hに対して、列ベクトルの大きさが左から昇順になるようなソート付きQR分解を行なうことで、ある順列行列Πを生成する。その後、Πの各行成分のうちeU,iを有する行を一番上に配置しなおせば良い。オーダリング部807出力として、オーダリング済みの伝搬路行列Hpと選択されたオーダリング行列Πがそれぞれ出力される。オーダリング部807出力Hpは次いで、格子基底縮小部803に入力され、格子基底縮小処理が施されることになる。図10は、本発明の第2の実施形態に係る格子基底縮小部803において用いられるアルゴリズムを示す表である。図10と第1の実施形態で対象とした図5の違いは、第2の端末装置300に係る伝搬路行列成分本実施形態では、第2の端末装置300は端末装置300iであるから、入力行列Hの第i列成分が該当するに対する信号処理である。第1の実施形態では、第i列成分に対して、準直交化処理こそ行なわないものの、2条件比較、および3入れ替え作業については適用している。この処理により、端末装置300iを考慮した格子基底縮小が行なえるため、たとえ第2の端末装置300が混在したとしても、良好な伝送特性を実現できる。しかし、LLLアルゴリズムは行列の全成分に対して、1~3の信号処理を行なうことを前提としているため、図5の方法では、特に第2の端末装置300数が複数である場合に、アルゴリズムの収束に膨大な時間を要してしまう。そこで、第2の実施形態が用いるアルゴリズムは、図10に示されているように、入力された行列に対して、所定の列番号vより信号処理を開始し、vより小さい列番号の列ベクトルについては、一切の信号処理を行なわないものとしている。これは、事前のオーダリングにより、第2の端末装置300に関わる伝搬路行列成分が第1行アルゴリズムに入力されるのはHHであるから、アルゴリズム中では第1列を指すに置換されているためである。この場合、アルゴリズムの収束性は図4で示されている従来のLLLアルゴリズムとほぼ同等とすることができる。本実施形態の場合、第2の端末装置300数は1個であるから、vとして3が入力されることになる。つまり、vには従来のLLLアルゴリズムのkの初期値である2に第2の端末装置300数が加算されたものとなる。格子基底縮小部803からは最終的にunimodular行列Tが出力される。次いで、線形フィルタ生成部203において、オーダリング部807出力Hpと格子基底縮小部803出力Tに基づいて、線形フィルタWが生成される。生成される線形フィルタWは第1の実施形態と同様に、格子基底縮小された伝搬路行列G=THp=TΠHを単位下三角行列に変換するものであるから、Gに対するLQ分解により得られる下三角行列Lとユニタリ行列QからW=QH{diagL}-1と算出できる。線形フィルタ生成部203出力Wは、事前干渉抑圧部805に入力される。事前干渉抑圧部805には、他にオーダリング部807出力Πと格子基底縮小部803出力Tとオーダリング済みの伝搬路行列Hpおよび、送信シンボルベクトルdが入力され、事前干渉抑圧が行なわれる。実際の信号処理は第1の実施形態とほぼ同様であり、式4のような送信信号ベクトルxを作り出すものであるが、若干処理が異なる。まず、第1の実施形態では、dにTを乗算したTdを新たな送信シンボルベクトルd^として干渉抑圧を行なっていたが、第2の実施形態においては、d^=TΠdを新たな送信シンボルベクトルとして干渉抑圧を行なう。そして、第1の実施形態においては、modulo演算の有無を、unimodular行列が有する単位行ベクトルの位置に応じて決定していたのに対して、第2の実施形態においては、d^のうち、第1成分から、第2の端末装置300数だけの成分までについてはmodulo演算を行なわないように制御する。例えば、本実施形態では第2の端末装置300数は1であるから、d^のうち、d^1に対してのみ、modulo演算を行なわない、つまり式31に基づいて干渉抑圧を行なうことになる。第2の端末装置300数がmであれば、d^1~d^mに対してmodulo演算を行なわないということである。それ以外の送信シンボルについては式33に基づいて干渉抑圧が行なわれる。以上の処理により、最終的に事前干渉抑圧部805より出力される送信符号ベクトルxは式8で与えられる。ただし、zt,1=0となる。事前干渉抑圧部805出力xは線形フィルタ乗算部207に入力され、線形フィルタWが乗算される。図11は、本発明の第2の実施形態に係るPrecoding部801の信号処理を示すフローチャートである。まず、Precoding部801は、送信シンボルベクトルdおよび伝搬路情報Hを取得するステップS301。次に、入力される伝搬路行列Hが初めにオーダリング部807に入力されてオーダリング処理が施されるステップS302。次に、格子基底縮小部803では入力された伝搬路行列Hpと図10記載のアルゴリズムに基づいて、Hpの直交性を高めるunimodular行列Tが算出されるステップS303。伝搬路情報Hpと格子基底縮小部501から出力されたTは線形フィルタ生成部203に入力され、線形フィルタWが生成されるステップS304。続いて、Precoding部801に入力された送信シンボルベクトルdは事前干渉抑圧部805に入力されることで、各端末装置300で観測される干渉成分が、予め減算されるステップS305。この時、第1の端末装置300宛てには式33に基づいて行なわれ、第2の端末装置300宛てには式31に基づいて行なわれる。次に、線形フィルタ乗算部207では、入力されたxに対して、線形フィルタWが乗算され、送信信号ベクトルsが算出されるステップS306。以降の信号処理については、Precoding部801出力に対する基地局装置100の信号処理を含めて図2と同じであるため、説明は省略する。[3.2端末装置300]第2の実施形態に係る端末装置300構成は、第1の端末装置300および第2の端末装置300ともに、それぞれ図7および図8と同じであり、行なわれる信号処理も同じであるから説明は省略する。第2の実施形態では、格子基底縮小部803において行なわれる格子基底縮小アルゴリズムを変更することで、格子基底縮小アルゴリズムの収束速度を低下させない方法を対象とした。本実施形態の方法によれば、伝搬路の直交性は若干低下するものの、特に第2の端末装置300数が多い場合に、格子基底縮小アルゴリズムを第1の実施形態の方法と比較して、高速に収束させることが可能である。[4.第3の実施形態]これまで説明してきた方法はいずれも、格子基底縮小部803で算出されるunimodular行列Tを従来LRATHPで用いられてきたTとは異なるものを使うことにより、第1の端末装置300と第2の端末装置300とを空間多重することを可能としてきた。しかし、Tの生成に制限を加えることは、伝搬路の直交性や、格子基底縮小アルゴリズムの収束速度を著しく低下させてしまう場合がある。ところで、LRATHPは伝搬路に対する準直交化処理と、THPによる事前干渉抑圧処理の二つの信号処理を組み合わせたものと捉える事ができ、これまでの方法は、準直交化処理の方を主に変更することで、空間多重を実現させてきたと言える。第3の実施形態では、THPによる事前干渉抑圧処理を変更することにより、空間多重を実現させる方法を対象とする。第1の実施形態と同様に、Nt本の送信アンテナ109を有する基地局装置100送信装置とも呼ぶに対して、1本のアンテナ301を有する端末装置300受信装置とも呼ぶがU個接続している通信を対象とし、端末装置300のうち、端末装置300iが第2の端末装置300であるものとする。なお、各端末装置300が有するアンテナ301数や、端末装置300に含まれる第1の端末装置300と第2の端末装置300の割合については、これに限ったものではなく、複数アンテナ301を有する端末装置300が存在したり、複数の第2の端末装置300が存在する場合にも、本発明は適用可能である。[4.1基地局装置100]第3の実施形態に係る基地局装置100構成は図2とほぼ同じであり、異なるのはPrecoding部801における信号処理のみであるため、以下では、第3の実施形態に係るPrecoding部について説明する。図12は、本発明の第3の実施形態に係るPrecoding部901の構成を示すブロック図である。装置構成は図9とほぼ同じであるが、第2の実施形態とはオーダリング部903の位置が異なり、また摂動ベクトル制御部905が追加される構成となる。はじめに、入力された伝搬路行列Hが格子基底縮小部907に入力され、格子基底縮小処理が施されることになる。第1および第2の実施形態においては、格子基底縮小部907で算出されるunimoular行列Tに制限を与えていたが、本実施形態においては、算出されるTはunimodular行列でありさえすればなんでも良い。例えば、図4のLLLアルゴリズムをそのまま適用し、unimoular行列Tを算出すれば良い。算出された行列Tはその後線形フィルタ生成部909と事前干渉抑圧部911に入力される。次いで、線形フィルタ生成部909には、格子基底縮小部907出力Tと伝搬路行列Hおよびオーダリング部903出力Πが入力される。オーダリング部903における信号処理については後述するが、初めはΠとして単位行列Iが入力されることになる。線形フィルタ生成部909では、入力された情報に基づき見掛け上の伝搬路行列G=ΠTHを生成する。第2の実施形態とは、unimodular行列Tと順列行列Πを乗算する順番が異なる。生成される線形フィルタWは伝搬路行列Gに対して、GWを単位下三角行列とするものであるから、第2の実施形態と同様に、Gに対するLQ分解により得られる下三角行列Lとユニタリ行列QからW=QH{diagL}-1と算出できる。線形フィルタ生成部909出力Wは、事前干渉抑圧部911に入力される。事前干渉抑圧部911には、他に伝搬路行列Hとオーダリング部903出力Πと格子基底縮小部907出力Tおよび、送信シンボルベクトルdが入力され、事前干渉抑圧が行なわれる。ここで行なわれる事前干渉抑圧は、伝搬路行列がHからGに、見掛け上の送信シンボルベクトルがd^=ΠTdに変更される以外は、基本技術で述べた事前干渉抑圧部911における信号処理と同じであり、式9のような送信符号ベクトルxを生成するものである。ここで、式9で表現される送信信号に線形フィルタWを乗算したものを基地局装置100より伝送したとすると、受信信号は式10で与えられることになる。ここで、第2の端末装置300である端末装置300iにおいて、modulo演算を行なわずとも信号復調をするためには、受信信号を構成する各成分のうち、列ベクトルT-1Π-1ztの第i列成分が0で有りさえすればよい。ここで、摂動ベクトルztは式31および式32から分かるように、伝搬路行列に応じて決まる残留IUIの大きさ式22中のai,ji>jが該当によって変化するから、伝搬路行列が変化すれば、摂動ベクトルの値も変わる。そこで、事前干渉抑圧部911出力xは一度、摂動ベクトル制御部905に入力され、第2の端末装置300で観測されるmodulo演算の摂動項T-1Π-1ztが0であるかどうかを測定し、0であれば、摂動ベクトル制御部905入力xをそのまま出力として、線形フィルタ乗算部に出力し、摂動項T-1Π-1ztが0でなければ、入力されたxを破棄し、オーダリング部903に対して、新たな順列行列を生成するよう制御信号を出力する。オーダリング部903における具体的な順列行列の探査方法としては、空間多重数がUである場合、選択可能な順列行列はU個だけ存在するから、それらすべてに関してT-1Π-1ztを算出し、第2の端末装置300に関わる成分が0となる順列行列を探せば良い。以上の処理により摂動ベクトル制御部905において適切な送信符号ベクトルxが得られた後、送信符号ベクトルxを線形フィルタ乗算部に入力する。線形フィルタ乗算部では、電力正規化などが行なわれ、Precoding部901出力sが生成される。図13は、本発明の第3の実施形態に係るPrecoding部901の信号処理について示すフローチャートである。まず、Precoding部901は、送信シンボルベクトルdおよび伝搬路情報Hを取得するステップS401。次に、格子基底縮小部907では入力された伝搬路行列Hと図4記載のアルゴリズムに基づいて、Hの直交性を高めるunimodular行列Tが算出されるステップS402。次に、オーダリング部903により、オーダリング行列Πとして単位行列Iが入力されるステップS403。伝搬路情報Hと格子基底縮小部501から出力されたTとオーダリング行列Πは線形フィルタ生成部203に入力され、線形フィルタWが生成されるステップS404。続いて、Precoding部901に入力された送信シンボルベクトルdは事前干渉抑圧部911に入力されることで、各端末装置300で観測される干渉成分が、予め減算されるステップS405。この時、常に式33に基づいて行なわれる。摂動ベクトル制御部905は、第2の端末装置300で観測される摂動ベクトルが0であるかどうかを判断しステップS406、0でなければステップS406:No、オーダリング部903に対して、オーダリング行列Πを更新するよう制御信号を出力しステップS407、ステップS404に戻る。0であればステップS406:Yes、線形フィルタ乗算部207では、入力されたxに対して、線形フィルタWが乗算され、送信信号ベクトルsが算出されるステップS408。ただし、全ての順列行列を探査したとしても、T-1Π-1ztの第2の端末装置300に関わる成分が0となる順列行列が抽出されるとは限らない。この場合、全ての順列行列の中で、該当要素を最も小さくできる順列行列を選択すれば良いが、その場合、第2の端末装置300の伝送特性は残留する摂動ベクトルの大きさに比例して劣化する。しかし、変調方式がQPSK等の位相変調であった場合、付与されている摂動ベクトルが送信データと同じ象限に存在すれば例えば、QPSK変調時に第2の端末装置300宛の送信信号が、2-1/21jである場合に、摂動ベクトルが正の整数p,qを用いてpjqと表現できる状態、データ復調には影響を与えないから、そのような摂動ベクトルとなる順列行列を選択しても良い。同様の事は、16QAMを用いる場合でも、信号候補点の位置によって摂動ベクトルの付与が信号復調に影響しない場合があることを利用できる。また、IUIが小さい場合にQPSK変調はmodulo演算を行なうと伝送特性は劣化してしまうから、空間多重されている端末装置300が全て第1の端末装置300であったとしても、QPSK変調を行なう端末装置300が混在している場合に、QPSK変調を行なう端末装置300宛の信号に付与される摂動ベクトルの存在事象が送信データと一致するような順列行列を選択し、端末装置300側ではmodulo演算を行なわないように制御しても良い。また、適切な順列行列が発見されない場合、各端末装置300宛の送信信号に対する干渉抑圧を行なわない方法も考えられる。つまり、IUIを完全抑圧に行なうためには、式31で表現される干渉抑圧を行なう必要があるが、一部の干渉を残留させれば、その後式32で付与される摂動ベクトルも変わることとなる。そこで、順列行列だけではなく、事前干渉抑圧部911における干渉抑圧を変化させることで、目的を達するように制御しても良い。その場合、残留IUIを許容する端末装置300の伝送特性は劣化してしまうことになるから、残留IUIを許容させる端末装置300は、元々受信品質に優れた端末装置300を選択することが望ましい。Precoding部901以外の信号処理については、全て第2の実施形態と同様であるから説明は省略する。[4.2端末装置300]第3の実施形態に係る端末装置300構成は、第1の端末装置300および第2の端末装置300ともに、それぞれ図7および図8と同様であり、行なわれる信号処理も第1および第2の実施形態で説明したものと同じであるから説明は省略する。第3の実施形態では、格子基底縮小部907において行なわれる格子基底縮小アルゴリズムを変化させず、事前干渉抑圧部911における信号処理を変更することで、第1の端末装置300と第2の端末装置300とを空間多重させる方法を対象とした。第1や第2の実施形態と異なり、事前干渉抑圧の方法を変更するため、空間多重された端末装置300間で伝送品質が一定とならないという問題はあるものの、格子基底縮小は完全に行なわれるため、伝搬路の直交性を維持しやすく、また収束アルゴリズムの収束速度に影響を与えずに空間多重を行なうことが可能となる。[5.第4の実施形態]これまで説明してきた方法はいずれも、格子基底縮小部907で算出されるunimodular行列Tと伝搬路行列Hに基づきTHPを行なうLRATHPを前提とした方法である。ところで、THPを含む非線形プリコーディングは全てVectorperturbationVPと呼ばれる技術に帰着することが知られている。VP技術の演算量を削減したものが、THPやLRATHPであると言える。第4の実施形態では、VPの考え方を応用して、格子基底縮小を用いる非線形MUMIMOで第1の端末装置300と第2の端末装置300とを空間多重させる方法を対象とする。第1の実施形態と同様に、Nt本の送信アンテナ109を有する基地局装置100送信装置とも呼ぶに対して、1本のアンテナ301を有する端末装置300受信装置とも呼ぶがU個接続している通信を対象とし、端末装置300のうち、端末装置300iが第2の端末装置300であるものとする。なお、各端末装置300が有するアンテナ301数や、端末装置300に含まれる第1の端末装置300と第2の端末装置300の割合については、これに限ったものではなく、複数アンテナ301を有する端末装置300が存在したり、複数の第2の端末装置300が存在する場合にも、本発明は適用可能である。[5.1基地局装置100]第4の実施形態に係る基地局装置100構成は図2とほぼ同じであり、異なるのはPrecoding部における信号処理のみであるため、以下では、第4の実施形態に係るPrecoding部について説明する。Precoding部の説明の前に、VP技術について簡単に説明する。式4で与えられるLRATHPの送信信号において、unimodular行列Tを単位行列Iに置き換えたものが、THPMUMIMOの送信信号ベクトルとなり、その後線形フィルタが乗算された送信信号sは式111で与えられる。ここで、摂動ベクトルztを零ベクトルに置きかえることで、ZF規範に基づく線形MUMIMOの送信信号に一致する。つまり、THPMUMIMOは線形MUMIMOの送信信号に摂動ベクトルを付与することにより送信電力を軽減したものと言える。ここで、THPでは式32に基づいて、各端末装置300で独立に摂動ベクトルを求めているから各端末装置300の所要送信電力を一定にできる。VPは空間多重され送信される送信信号自体の送信電力を最小とする摂動ベクトルを求める方法であり、その送信信号は式112のように表現できる。しかし、候補摂動ベクトルは無数に存在し、空間多重数が増加するにつれて、候補数は指数関数的に増加していくため、VPをそのまま用いるのは非現実的である。ところで、LRATHPの線形フィルタ乗算後の送信信号は式113で与えられることになる。基本的には式111と同様であるが、摂動ベクトルにT-1が乗算された形となる。この場合、ztはTHPにより送信符号毎に独立に求められたものであるが、T-1は与えられた伝搬路行列を直交化させるためのものであり、伝搬路行列全体を考慮して算出された行列であるから、z^t=T-1ztは送信信号全体の送信電力をztよりも大幅に抑圧できる摂動ベクトルとなる。よって、この摂動ベクトルを基にVPの摂動ベクトルを探査することで、演算量の大幅な削減が行なえることが知られている。本実施形態では、この考えに基づき、第1の端末装置300と第2の端末装置300とを空間多重させるプリコーディングを行なう。図14は、本発明の第4の実施形態に係るPrecoding部1001の構成を示すブロック図である。図14に示すように、Precoding部1001は、線形フィルタ生成部1003と、格子基底縮小部1005と、摂動ベクトル制御部1007とを含んで構成されている。はじめに、格子基底縮小部1005に伝搬路行列Hが入力され、格子基底縮小が行なわれる。実際の信号処理は基本技術で述べた方法と同じであり、例えば、図4に示したLLLアルゴリズムに基づいてunimodular行列Tを求めれば良い。格子基底縮小部1005より出力された行列Tは線形フィルタ生成部1003に入力される。線形フィルタ生成部1003には格子基底縮小部1005出力Tと伝搬路行列Hが入力されて、線形フィルタが生成される。線形フィルタ生成部1003において生成される線形フィルタWは格子基底縮小された伝搬路行列G=THを対角行列にするフィルタであり、例えば、Gの逆行列G-1を用いれば良い。また、MMSE規範に基づく線形フィルタを用いても良い。生成された線形フィルタWは摂動ベクトル制御部1007に入力されることになる。摂動ベクトル制御部1007には、送信シンボルベクトルdと線形フィルタWおよびunimodular行列Tが入力され、摂動ベクトルが算出される。本実施形態では、摂動ベクトルとして、第2の端末装置300に属する端末装置300iにおいて、modulo演算を必要としないような摂動ベクトルを探査する。LRATHPにおいて、THPによる事前干渉抑圧を行なわなかった場合の送信信号もやはり式113のように表現することができる。摂動ベクトル制御部1007では、はじめに式113で表現される信号を算出する。このとき、摂動ベクトルz^t=T-1ztを第1の摂動ベクトルと呼ぶこととする。しかし、第2の端末装置300に属する端末装置300i宛の送信符号xiには、z^t,iが付与される事となり、このままでは端末装置300iではmodulo演算を行なわないと信号復調することが出来ない。そこで、摂動ベクトル制御部1007では、z^tに加えて、新たな摂動ベクトルzTを付与する。zTを第2の摂動ベクトルと呼ぶ。つまり、最終的な送信信号を式114のようにすることを意味する。このとき、新たに付与する摂動ベクトルzTの第i成分zT,iについては、zT,iz^t,i=0を満たすように制御する。このように摂動ベクトルを新たに付与することにより、端末装置300iではmodulo演算を用いずとも信号復調することができる。zT,i以外の摂動ベクトルについては、特に制限は存在しないが、適切な摂動ベクトルを選択しないと、伝送特性は大幅に劣化してしまう。伝送特性を劣化させないためには、zT,i以外の摂動ベクトルについては、送信信号電力を最小にするように選択する必要がある。つまり、求めるべき摂動ベクトルzTは式12を満たすものとなる。以上の信号処理により、摂動ベクトルzTを算出したのち、電力正規化係数βを算出する。摂動ベクトル制御部1007からは電力正規化係数βが乗算された送信信号ベクトルsが出力され、Precoding部1001出力として出力される。このとき、sは式13で与えられる。図15は、本発明の第4の実施形態に係るPrecoding部1001の信号処理について示すフローチャートである。まず、Precoding部1001は、送信シンボルベクトルdおよび伝搬路情報Hを取得するステップS501。次に、格子基底縮小部1005では入力された伝搬路行列Hと図4記載のアルゴリズムに基づいて、Hの直交性を高めるunimodular行列Tが算出されるステップS502。伝搬路情報Hと格子基底縮小部1005から出力されたTは線形フィルタ生成部1003に入力され、線形フィルタWが生成されるステップS503。次に、摂動ベクトル制御部1007は、第1の摂動ベクトルを生成するステップS504。続いて、摂動ベクトル制御部1007は、伝搬路情報Hと格子基底縮小部1005から出力されたTと線形フィルタWに基づいて第2の摂動ベクトルを生成するステップS505。摂動ベクトル制御部1007は、dに第1および第2の摂動ベクトルを加算し、送信符号ベクトルxを生成するステップS506。更に、摂動ベクトル制御部1007は、送信符号ベクトルxに線形フィルタWを乗算し、送信信号ベクトルsを生成するステップS507。以降の基地局装置100の信号処理については、図2と同じであるから説明は省略する。[5.2端末装置300]第4の実施形態に係る端末装置300構成は、第1の端末装置300および第2の端末装置300ともに、それぞれ図7および図8と同様であり、行なわれる信号処理も第1から第3の実施形態と同じであるから説明は省略する。第4の実施形態では、格子基底縮小に基づく非線形MUMIMO伝送において、摂動ベクトルを適切に制御することにより、第1の端末装置300と第2の端末装置300とを空間多重させる方法を対象とした。第1から第3の実施形態で対象としたように、摂動ベクトルを端末装置300毎に算出しないため、演算量は大幅に増加してしまうものの、第2の端末装置300を空間多重させることにより発生する伝送特性の劣化を最小限に抑えることが可能となる。<全実施形態共通>[変形例]以上、この発明の実施形態について図面を参照して詳述してきたが、具体的な構成はこの実施形態に限られるものではなく、この発明の要旨を逸脱しない範囲の設計等も特許請求の範囲に含まれる。本発明に関わる端末装置300および基地局装置100で動作するプログラムは、本発明に関わる上記実施形態の機能を実現するように、CPU等を制御するプログラムコンピュータを機能させるプログラムである。そして、これら装置で取り扱われる情報は、その処理時に一時的にRAMに蓄積され、その後、各種ROMやHDDに格納され、必要に応じてCPUによって読み出し、修正書き込みが行なわれる。プログラムを格納する記録媒体としては、半導体媒体例えば、ROM、不揮発性メモリカード等、光記録媒体例えば、DVD、MO、MD、CD、BD等、磁気記録媒体例えば、磁気テープ、フレキシブルディスク等等のいずれであってもよい。また、ロードしたプログラムを実行することにより、上述した実施形態の機能が実現されるだけでなく、そのプログラムの指示に基づき、オペレーティングシステムあるいは他のアプリケーションプログラム等と共同して処理することにより、本発明の機能が実現される場合もある。また市場に流通させる場合には、可搬型の記録媒体にプログラムを格納して流通させたり、インターネット等のネットワークを介して接続されたサーバコンピュータに転送したりすることができる。この場合、サーバコンピュータの記憶装置も本発明に含まれる。また、上述した実施形態における端末装置300および基地局装置100の一部、または全部を典型的には集積回路であるLSIとして実現してもよい。端末装置300および基地局装置100の各機能ブロックは個別にプロセッサ化してもよいし、一部、または全部を集積してプロセッサ化してもよい。また、集積回路化の手法はLSIに限らず専用回路、または汎用プロセッサで実現しても良い。また、半導体技術の進歩によりLSIに代替する集積回路化の技術が出現した場合、当該技術による集積回路を用いることも可能である。以上、この発明の実施形態を、図面を参照して詳述してきたが、具体的な構成はこの実施形態に限られるものではなく、この発明の要旨を逸脱しない範囲の設計等も特許請求の範囲に含まれる。本発明は無線通信装置に利用可能である。"