"自局宛の送信信号を送信する他局の送信アンテナと自局の受信アンテナとの間の自チャネルと、他局宛の送信信号を送信する他局の送信アンテナと前記受信アンテナとの間の干渉チャネルと、を推定する推定部と、前記推定部による前記自チャネルの推定値と前記干渉チャネルの推定値とを行方向に結合したチャネル行列を生成する生成部と、前記生成部によって生成されたチャネル行列を、直交行列と、前記自チャネルの成分が前記干渉チャネルの成分よりゼロ成分の多い行に配置された三角行列と、に分解する分解部と、前記受信アンテナによって受信された受信信号を、前記分解部によって得られた直交行列に基づいて直交化する直交化部と、前記分解部によって得られた三角行列の行ごとの成分をゼロ成分のより多い行から順に取得し、前記行ごとの成分を取得するごとに、取得した行ごとの成分と前記直交化部によって直交化された受信信号とに基づいて前記自局宛の送信信号の値のメトリックを累積的に算出する算出部と、前記算出部によって算出されたメトリックに基づいて前記自局宛の送信信号を復号する復号部と、を備えることを特徴とする受信装置。前記算出部は、前記行ごとの成分を、前記自局宛の送信信号を送信する他局の送信アンテナの数だけ取得することを特徴とする請求項1に記載の受信装置。前記生成部は、前記自チャネルの推定値を前記干渉チャネルの推定値より右側に配置したチャネル行列を生成し、前記分解部は、前記チャネル行列を直交行列および上三角行列に分解し、前記算出部は、前記行ごとの成分を前記上三角行列の下段から順に取得することを特徴とする請求項1または2に記載の受信装置。前記推定部による各推定値に基づく雑音電力を推定する雑音推定部と、前記生成部によって生成されたチャネル行列の列数に等しい行数および列数の単位行列に前記雑音推定部によって推定された雑音電力の偏差を乗算した行列を、前記チャネル行列の下側に追加して前記チャネル行列を拡張する拡張部と、を備え、前記分解部は、前記拡張部によって拡張されたチャネル行列を分解することを特徴とする請求項1~3のいずれか一つに記載の受信装置。前記推定部は、自局宛の送信信号を送信する他局の複数の送信アンテナと前記受信アンテナとの間の各自チャネルを推定し、他局宛の送信信号を送信する他局の複数の送信アンテナと前記受信アンテナとの間の各干渉チャネルを推定し、前記生成部は、前記自チャネルの各推定値を前記自チャネルの各推定値のベクトルノルムに基づいて並べ、前記干渉チャネルの各推定値を前記干渉チャネルの各推定値のベクトルノルムに基づいて並べたチャネル行列を生成することを特徴とする請求項1~4のいずれか一つに記載の受信装置。前記推定部による各推定値に基づく雑音電力を推定する雑音推定部と、前記生成部によって生成されたチャネル行列と同じ列数および前記自局宛の送信信号の数に等しい行数を有し、対角成分の値が前記雑音推定部によって推定された雑音電力の偏差であり非対角成分の値がゼロである行列を、前記チャネル行列の下側に追加して前記チャネル行列を拡張する拡張部と、を備え、前記分解部は、前記拡張部によって拡張されたチャネル行列を分解することを特徴とする請求項1~3のいずれか一つに記載の受信装置。前記推定部による各推定値に基づく信号雑音比を推定する信号雑音比推定部を備え、前記拡張部は、前記信号雑音比推定部によって推定された信号雑音比と閾値との比較結果に基づいて前記チャネル行列の拡張の有無を切り替え、前記分解部は、前記拡張部によって拡張されたチャネル行列、または前記拡張部によって拡張されなかったチャネル行列を分解することを特徴とする請求項4または6に記載の受信装置。自局宛の送信信号を送信する他局の送信アンテナと自局の受信アンテナとの間の自チャネルと、他局宛の送信信号を送信する他局の送信アンテナと前記受信アンテナとの間の干渉チャネルと、を推定し、前記自チャネルの推定値と、前記干渉チャネルの推定値と、を行方向に結合したチャネル行列を生成し、生成した前記チャネル行列を、直交行列と、前記自チャネルの成分が前記干渉チャネルの成分よりゼロ成分の多い行に配置された三角行列と、に分解し、前記受信アンテナによって受信された受信信号を、前記直交行列に基づいて直交化し、前記三角行列の行ごとの成分をゼロ成分のより多い行から順に取得し、前記行ごとの成分を取得するごとに、取得した行ごとの成分と直交化した前記受信信号とに基づいて前記自局宛の送信信号の値のメトリックを累積的に算出し、算出した前記メトリックに基づいて前記自局宛の送信信号を復号する、ことを特徴とする受信方法。しかしながら、上述した従来技術では、他局宛の送信信号の変調方式が不明であるシステムにおいては、他局宛の送信信号の値の候補変調多値数も不明であり、MLD法による信号分離が利用できない。このため、MLD法を利用することによって受信特性を向上させることができないという問題がある。本発明は、上述した従来技術による問題点を解消するため、受信特性を向上させることができる受信装置および受信方法を提供することを目的とする。上述した課題を解決し、目的を達成するため、本発明の一側面によれば、自局宛の送信信号を送信する他局の送信アンテナと自局の受信アンテナとの間の自チャネルを推定し、他局宛の送信信号を送信する他局の送信アンテナと前記受信アンテナとの間の干渉チャネルを推定し、前記自チャネルの推定値と、前記干渉チャネルの推定値と、を行方向に結合したチャネル行列を生成し、生成した前記チャネル行列を、直交行列と、前記自チャネルの成分が前記干渉チャネルの成分よりゼロ成分の多い行に配置された三角行列と、に分解し、前記受信アンテナによって受信された受信信号を、前記直交行列に基づいて直交化し、前記三角行列の行ごとの成分をゼロ成分のより多い行から順に取得し、前記行ごとの成分を取得するごとに、取得した行ごとの成分と直交化した前記受信信号とに基づいて前記自局宛の送信信号の値のメトリックを累積的に算出し、算出した前記メトリックに基づいて前記自局宛の送信信号を復号する、受信装置および受信方法が提案される。本発明の一側面によれば、受信特性を向上させることができるという効果を奏する。以下に添付図面を参照して、本発明にかかる受信装置および受信方法の実施の形態を詳細に説明する。実施の形態1図1は、実施の形態1にかかる通信システムの一例を示す図である。図1に示すように、実施の形態1にかかる通信システム100は、送信機110と、ユーザ1~KKは1以上の自然数に対応するK個の受信機121~12Kと、を含んでいる。通信システム100は、複数の送信アンテナを備える送信機110から、複数の受信アンテナを備える受信機121~12Kへ複数のデータストリームを同時に無線送信するMU-MIMOの通信システムである。たとえば、送信機110は、移動体通信システムにおける基地局に適用することができる。また、受信機121~12Kは、移動体通信システムにおける移動局に適用することができる。受信機121~12Kの各アンテナは、自局宛の送信信号だけでなく、他局宛の送信信号も受信する。受信機121~12Kの少なくともいずれかは、受信した各送信信号のうちの自局宛の送信信号をMLD法により分離して復号する。送信機の構成図2は、送信機の構成例を示す図である。図2に示す送信機200は、図1に示した送信機110の構成例である。図2に示すように、送信機200は、スケジューリング部210と、ユーザ1~Kに対応するK個のユーザデータ生成部221~22Kと、K個の誤り訂正符号化部231~23Kと、K個の変調部241~24Kと、ストリームマッピング部250と、Nt個Ntは2以上の自然数の送信部261~26Ntと、Nt本の送信アンテナ271~27Ntと、を備えている。したがって、送信機200が多重可能な最大のデータストリーム数はNtである。また、送信機200が同時にスケジュール可能なユーザ数は最大でNtである。<スケジューリング部>スケジューリング部210は、受信側のユーザ1~Kとの通信についてスケジューリングを行う。たとえば、スケジューリング部210は、ユーザ1~Kのチャネル品質に応じて、スケジュールするユーザを決定し、スケジュールされるユーザごとのデータストリーム数、符号化率および変調方式などを決定する。スケジューリング部210は、ユーザ1~Kのスケジューリング結果をそれぞれユーザデータ生成部221~22Kへ出力する。<ユーザデータ生成部>ユーザデータ生成部221~22Kは、スケジューリング部210から出力されたスケジューリング結果に基づいて、それぞれユーザ1~Kへ送信するユーザデータを生成する。ユーザデータ生成部221~22Kは、生成したユーザデータをそれぞれ誤り訂正符号化部231~23Kへ出力する。<誤り訂正符号化部>誤り訂正符号化部231~23Kは、それぞれユーザデータ生成部221~22Kから出力されたユーザデータの誤り訂正符号化を行う。誤り訂正符号化部231~23Kは、たとえば、それぞれユーザ1~Kのチャネル品質に応じた符号化率によって誤り訂正符号化を行ってもよい。誤り訂正符号化部231~23Kは、誤り訂正符号化を行ったユーザデータをそれぞれ変調部241~24Kへ出力する。<変調部>変調部241~24Kは、それぞれ誤り訂正符号化部231~23Kから出力されたユーザデータに基づく変調を行う。変調部241~24Kは、たとえば、それぞれユーザ1~Kのチャネル品質に応じた変調方式によって変調を行ってもよい。変調部241~24Kは、変調を行ったユーザデータをストリームマッピング部250へ出力する。<ストリームマッピング部>ストリームマッピング部250は、変調部241~24Kから出力されたユーザ1~Kのユーザデータを、送信アンテナ271~27Ntにより送信される各データストリームにマッピング割り当てする。ストリームマッピング部250は、マッピング結果に基づいて、ユーザ1~Kのユーザデータをそれぞれ送信部261~26Ntのうちの対応する送信部へ出力する。<送信部および送信アンテナ>送信部261~26Ntのそれぞれは、ストリームマッピング部250から出力されたユーザデータを無線周波数にアップコンバートする。送信部261~26Ntは、アップコンバートしたユーザデータをそれぞれ送信アンテナ271~27Ntへ出力する。これにより、送信アンテナ271~27Ntからユーザデータが無線信号により同時送信される。送信アンテナ271~27Ntからの送信信号をそれぞれ送信信号x1~xNtとする。受信機の構成図3は、受信機の構成例を示す図である。図3に示す受信機300は、図1に示した受信機121~12Kのそれぞれの構成例である。ここでは、一例としてユーザ1の受信機121の構成例としての受信機300について。図3に示す例では、図2に示した送信機200においてユーザ1にスケジュールされたデータストリームの数をL≦Nrとする。すなわち、送信機200の送信アンテナ271~27Ntのうちの1~L番目の送信アンテナにユーザ1へのユーザデータが割り当てられたとする。受信機300は、たとえば、送信アンテナ271~27Ntの数Ntと、送信アンテナ271~27Ntのうちの受信機300を宛先とする送信信号を送信する送信信号の数Lと、をあらかじめ送信機200から取得する。送信機200は、たとえば、ストリームマッピング部250によるマッピングの結果としてNtおよびLを受信機300へ通知する。また、受信機300は、送信機200が受信機300へ送信する送信信号の変調方式を示す情報をあらかじめ送信機200から取得する。送信機200は、たとえば、ストリームマッピング部250によるマッピングの結果として、受信機300へ送信する送信信号の変調方式を受信機300へ通知する。図3に示すように、受信機300は、Nr本Nrは1以上の自然数の受信アンテナ311~31Nrと、Nr個の受信部321~32Nrと、自チャネル推定部331と、干渉チャネル推定部332と、チャネル行列生成部333と、雑音電力推定部334と、チャネル行列拡張部335と、QR分解部336と、ユニタリ変換部337と、MLD処理部338と、誤り訂正復号部339と、を備えている。受信部321~32Nr、自チャネル推定部331、干渉チャネル推定部332、チャネル行列生成部333、雑音電力推定部334、チャネル行列拡張部335、QR分解部336、ユニタリ変換部337、MLD処理部338および誤り訂正復号部339は、たとえばFPGAFieldProgrammableGateArrayなどの電子回路によって実現することができる。<受信アンテナ>受信アンテナ311~31Nrのそれぞれは、送信機200の送信アンテナ271~27Ntから無線信号により送信された送信信号x1~xNtを受信する。すなわち、受信アンテナ311~31Nrのそれぞれが受信する送信信号には、受信機300を宛先とする送信信号x1~xLだけでなく、受信機300を宛先としない送信信号xL1~xNtも含まれる。受信アンテナ311~31Nrは、受信した信号受信信号をそれぞれ受信部321~32Nrへ出力する。<受信部>受信部321~32Nrは、それぞれ受信アンテナ311~31Nrから出力された受信信号をベースバンド信号に変換する。受信部321~32Nrのそれぞれは、ベースバンド信号に変換した受信信号を、自チャネル推定部331、干渉チャネル推定部332およびユニタリ変換部337へ出力する。受信部321~32Nrから出力される受信信号をそれぞれ受信信号y1~yNrとする。<自チャネル推定部>自チャネル推定部331は、受信部321~32Nrから出力される受信信号y1~yNrに含まれるパイロット信号などに基づいて自チャネルを推定する。具体的には、自チャネル推定部331は、受信アンテナ311~31Nrのうちのm1≦m≦Nr番目の受信アンテナと、送信アンテナ271~27Ntのうちのn1≦n≦L番目の送信アンテナと、の間の各チャネルを推定する。チャネルの推定は、たとえば、チャネルの振幅および位相の変動量の推定である。自チャネル推定部331は、自チャネルの推定値を自チャネル推定値ベクトルhnとしてチャネル行列生成部333および雑音電力推定部334へ出力する。<干渉チャネル推定部>干渉チャネル推定部332は、受信部321~32Nrから出力される受信信号y1~yNrに含まれるパイロット信号などに基づいて自チャネルへの干渉チャネルを推定する。具体的には、干渉チャネル推定部332は、受信アンテナ311~31Nrのうちのm番目の受信アンテナと、送信アンテナ271~27Ntのうちのユーザ1以外へ信号を送信している可能性のある送信アンテナと、の間のチャネルを推定する。ユーザ1以外へ信号を送信している可能性のある送信アンテナは、たとえば、ユーザ1へ信号を送信していない送信アンテナである。すなわち、ユーザ1以外へ信号を送信している可能性のある送信アンテナは、送信アンテナ271~27NtのうちのkL1≦k≦Nt番目の送信アンテナである。干渉チャネル推定部332は、チャネルの推定値を、干渉チャネル推定値ベクトルhkとしてチャネル行列生成部333および雑音電力推定部334へ出力する。<チャネル行列生成部>チャネル行列生成部333は、自チャネル推定部331から出力された自チャネル推定値ベクトルhnと、干渉チャネル推定部332から出力された干渉チャネル推定値ベクトルhkと、を行方向行の各成分の並び方向に結合したチャネル行列Hを生成する。たとえば、チャネル行列生成部333は、自チャネル推定値ベクトルhnが干渉チャネル推定値ベクトルhkより右側に配置されるようにチャネル行列Hを生成する。チャネル行列生成部333によって生成されるチャネル行列Hは、たとえば下記1式によって表される。また、送信信号ベクトルをxx1~xNt、雑音をnn1~nNrとすると、受信機300の受信信号yy1~yNrはy=Hxnとなる。したがって、受信信号yはたとえば下記2式によって表される。上述のように、チャネル行列Hにおいては自チャネル推定値ベクトルhnが干渉チャネル推定値ベクトルhkより右側に配置されている。このため、上記2式に示すように、送信信号ベクトルxにおいては、自局宛の送信信号x1,x2,...,xLが他局宛の送信信号xL1,xL2,...,xNtより下段になる。これにより、後述するMLD処理部338において、自局宛の送信信号x1,x2,...,xLが他局宛の送信信号xL1,xL2,...,xNtより先に処理される。チャネル行列生成部333は、生成したチャネル行列Hをチャネル行列拡張部335へ出力する。このように、チャネル行列生成部333は、チャネル行列Hに対応する送信信号ベクトルxの所定方向ここでは下段から上段へ向かう方向の順序において、自局宛の送信信号が他局宛の送信信号より先になるようにチャネル行列Hを生成する。<雑音電力推定部>雑音電力推定部334は、自チャネル推定部331から出力された自チャネル推定値ベクトルhnと、干渉チャネル推定部332から出力された干渉チャネル推定値ベクトルhkと、に基づいて、受信アンテナあたりの雑音電力σ2を推定する雑音推定部である。雑音電力推定部334は、推定した雑音電力σ2をチャネル行列拡張部335へ出力する。<チャネル行列拡張部>チャネル行列拡張部335は、チャネル行列生成部333から出力されたチャネル行列Hを、たとえば下記3式に示すように拡張して拡張チャネル行列Heとする。INtはチャネル行列Hの列数Ntに等しい行数および列数Nt×Ntの単位行列である。σINtは、INtに雑音電力σ2の偏差を乗算した行列である。チャネル行列拡張部335は、チャネル行列生成部333から出力されたチャネル行列Hの下側にσINtを追加することによって拡張を行う。したがって、拡張チャネル行列HeはNrNt×Ntの行列である。チャネル行列拡張部335は、チャネル行列Hを拡張した拡張チャネル行列HeをQR分解部336へ出力する。たとえば、Nt=Nr=4、L=2の場合は、チャネル行列Hは下記4式のようになり、拡張チャネル行列Heは下記5式のようになる。<QR分解部>QR分解部336は、チャネル行列拡張部335から出力された拡張チャネル行列HeをQR分解する。たとえば、QR分解部336は、下記6式に示すように、拡張チャネル行列HeをNrNt×NrNtのユニタリ行列Qe直交行列とNrNt×Ntの上三角行列Reとに分解する。上三角行列Reは、自チャネルの成分が干渉チャネルの成分よりゼロ成分の多い行下段に配置された三角行列となる。QR分解部336は、QR分解により得られたユニタリ行列Qeの左上Nr×Nrを取り出した部分行列Qのエルミート共役QHをユニタリ変換部337へ出力する。また、QR分解部336は、QR分解により得られた上三角行列Reの上Nr行を取り出した部分行列RをMLD処理部338へ出力する。これにより、ユニタリ変換部337において受信信号yの直交化が可能になる。<ユニタリ変換部>ユニタリ変換部337は、受信部321~32Nrから出力された受信信号y=y1~yNrの直交化を行う直交化部である。具体的には、ユニタリ変換部337は、QR分解部336から出力されたエルミート共役QHを受信信号yに乗算z=QHyし、直交化した受信信号zz1~zNrを得る。ユニタリ変換部337は、直交化した受信信号zをMLD処理部338へ出力する。MLD処理部338は、ユニタリ変換部337から出力された受信信号zに対するMLD法による処理を、QR分解部336から出力された部分行列Rに基づいて行う。受信信号zは、たとえば下記7式に示すように表すことができる。下記7式、受信信号z、部分行列R、エルミート共役QHに基づいて送信信号xを求めることができる。<MLD処理部>MLD処理部338は、ユニタリ変換部337から出力された受信信号zと、QR分解部336から出力された部分行列Rと、に基づいて、自局宛の各送信信号の値の各候補のメトリックを累積的に算出する算出部である。具体的には、MLD処理部338は、部分行列R三角行列の行ごとの成分を、ゼロ成分のより多い行ここでは下段から順に、自局宛の送信信号を送信する他局の送信アンテナの数Lだけ取得する。そして、MLD処理部338は、行ごとの成分を取得するごとに、取得した行ごとの成分と受信信号zとに基づいて自局宛の送信信号の値のメトリックを累積的に算出する。また、MLD処理部338は、自局宛の各送信信号のメトリックの算出が終わると、算出した累積メトリックに基づいて、自局宛の各送信信号が示す値の候補のLLRLogLikelihoodRatio:対数尤度比を誤り訂正復号部339へ出力する。したがって、MLD処理部338は、他局宛の各送信信号のメトリックについては算出しなくてもよい。<誤り訂正復号部>誤り訂正復号部339は、MLD処理部338から出力されたLLRに基づく軟判定により、送信信号xが示す値の誤り訂正復号を行う。誤り訂正復号部339は、誤り訂正復号により得られたユーザデータを出力する。MLD処理部の具体例図4は、MLD処理部の具体例を示す図である。ここでは一例として、MLD処理部338が、QR分解部336によるQR分解の結果に基づいてMアルゴリズムによるシンボル候補の絞り込みを行うQRM-MLDアルゴリズムを用いる場合について説明する。Mアルゴリズムにおいては、受信アンテナ311~31Nrに対応する各ステージにおいて、N個のシンボル候補の中からMN≦M個のシンボル候補を選択することにより順次シンボル候補の絞り込みを行う。MLD処理部338は、ステージ処理制御部410と、Nr個のステージ処理部421~42Nrと、LLR算出部430と、を備えている。MLD処理部338には、ユニタリ変換部337から出力された受信信号zと、QR分解部336から出力された部分行列Rと、が入力される。<ステージ処理制御部>ステージ処理制御部410は、自局宛の送信信号x1~xLまでのステージの処理が実行され、他局宛の送信信号xL1~xNtまでのステージの処理が実行されないようにステージ処理部421~42Nrの動作を制御する。具体的には、ステージ処理制御部410は、第1ステージ~第Lステージに対応するステージ処理部421~42Lを動作するように設定する。また、ステージ処理制御部410は、第L1ステージ~第Nrステージに対応するステージ処理部42L1~42Nrは動作しないように設定する。また、ステージ処理制御部410は、第1ステージ~第L-1ステージに対応するステージ処理部421~42L-1の出力選択部403の出力先をシンボル候補選択部404に設定する。また、ステージ処理制御部410は、第Lステージに対応するステージ処理部42Lの出力選択部403の出力先をLLR算出部430に設定する。<ステージ処理部>ステージ処理部421~42Nrは、それぞれMLD法における第1ステージ~第Nrステージの処理を行う。具体的には、ステージ処理部421~42Nrのそれぞれは、二乗ユークリッド距離算出部401と、累積メトリック算出部402と、出力選択部403と、シンボル候補選択部404と、を備えている。ただし、ステージ処理部421は累積メトリック算出部402を備えていなくてもよい。また、ステージ処理部42Nrはシンボル候補選択部404を備えていなくてもよい。<第1ステージ>まず、第1ステージに対応するステージ処理部421の二乗ユークリッド距離算出部401、出力選択部403およびシンボル候補選択部404について説明する。ステージ処理部421は、受信信号zおよび部分行列Rに基づいて、送信信号xLが示す値の候補シンボル候補を選択する。具体的には、二乗ユークリッド距離算出部401は、送信信号xLが示す値の各候補CLii=1,2,...,mLについて、受信信号zNrとの二乗ユークリッド距離を算出する。mLは、送信信号xLの変調多値数たとえば位相点の数である。たとえば変調方式がQPSKQuadraturePhaseShiftKeyingである場合はmL=4となる。たとえば、二乗ユークリッド距離算出部401は、下記8式によって二乗ユークリッド距離d1cLiを算出する。二乗ユークリッド距離算出部401は、算出した二乗ユークリッド距離d1cLiを、第1ステージまでの累積メトリックとして出力選択部403へ出力する。出力選択部403は、ステージ処理制御部410からの設定により、二乗ユークリッド距離算出部401から出力された二乗ユークリッド距離d1cLiをシンボル候補選択部404へ出力する。シンボル候補選択部404は、送信信号xLが示す値の各候補CLiの中から、出力選択部403から出力された二乗ユークリッド距離d1cLiが小さい順にS1個の候補cLkを選択するk=1,2,...,S1。S1は、たとえばあらかじめ設定された固定値である。シンボル候補選択部404は、選択した候補cLkをステージ処理部422へ出力する。また、シンボル候補選択部404は、出力選択部403から出力された二乗ユークリッド距離d1cLiのうちの、選択した候補cLkに対応する二乗ユークリッド距離を第1ステージまでの累積メトリックd1cLkとしてステージ処理部422へ出力する。<第2ステージ>つぎに、第2ステージに対応するステージ処理部422の二乗ユークリッド距離算出部401、累積メトリック算出部402、出力選択部403およびシンボル候補選択部404について説明する。ステージ処理部422は、ステージ処理部421から出力された送信信号xLの候補cLkおよび累積メトリックd1cLkに基づいて、送信信号xL,xL1が示す値の組み合わせの候補を選択する。具体的には、二乗ユークリッド距離算出部401は、送信信号xL,xL1が示す値の各候補CL1ii=1,2,...,mL1について、受信信号zNr1との二乗ユークリッド距離を算出する。ただし、送信信号xLの値については、ステージ処理部421から出力された候補cLkの値とする。たとえば、二乗ユークリッド距離算出部401は、下記9式によって二乗ユークリッド距離d2cLk,cL1iを算出する。二乗ユークリッド距離算出部401は、算出した二乗ユークリッド距離d2cLk,cL1iを累積メトリック算出部402へ出力する。累積メトリック算出部402は、第2ステージまでのS1mL1個の累積メトリックE2を算出する。たとえば、累積メトリック算出部402は、下記10式のように累積メトリックE2を算出する。すなわち、累積メトリック算出部402は、ステージ処理部421からの第1ステージまでの累積メトリックd1cLと、二乗ユークリッド距離算出部401からの第2ステージの二乗ユークリッド距離d2cLk,cL1iと、を加算する。累積メトリック算出部402は、算出した累積メトリックE2を出力選択部403へ出力する。出力選択部403は、ステージ処理制御部410からの設定により、累積メトリック算出部402から出力された累積メトリックE2をシンボル候補選択部404へ出力する。シンボル候補選択部404は、送信信号xL,xL1が示す値の組み合わせの各候補CL1iの中から、出力選択部403から出力された累積メトリックE2が小さい順にS2個の候補CL1kを選択するk=1,2,...,S2。ただし、送信信号xLの値の候補については、ステージ処理部421から出力された候補cLkとする。S2は、たとえばあらかじめ設定された固定値である。シンボル候補選択部404は、選択した候補CL1kをステージ処理部423へ出力する。また、シンボル候補選択部404は、出力選択部403から出力された累積メトリックE2のうちの、選択した候補CL1kに対応する累積メトリックE2を第2ステージまでの累積メトリックd2cL1kとしてステージ処理部423へ出力する。<第3ステージ~第L-1ステージ>つぎに、第3ステージ~第L-1ステージに対応するステージ処理部423~42L-1について説明する。ステージ処理部423~42L-1は、ステージ処理部422と同様の処理を行う。たとえば、ステージ処理部42L-1は、前段のステージ処理部42L-2から出力された送信信号の候補および累積メトリックに基づいて、それぞれ送信信号xL~x2の値の組み合わせの候補C2を選択する。ステージ処理部42L-1のシンボル候補選択部404は、選択した候補C2kをステージ処理部42Lへ出力する。また、ステージ処理部42L-1のシンボル候補選択部404は、選択した候補C2kに対応する累積メトリックEL1を第L-1ステージまでの累積メトリックdL1c2kとしてステージ処理部42Lへ出力する。<第Lステージ>つぎに、第Lステージに対応するステージ処理部42Lの二乗ユークリッド距離算出部401、累積メトリック算出部402および出力選択部403について説明する。ステージ処理部42Lは、ステージ処理部42L-1から出力された送信信号x2の候補CLiおよび累積メトリックdL1c2kに基づいて、送信信号xL~x1の値の組み合わせの候補C1を選択する。具体的には、二乗ユークリッド距離算出部401は、送信信号xL~x1の値の組み合わせの各候補CL1ii=1,2,...,mL1について、受信信号zNrL1との二乗ユークリッド距離を算出する。ただし、送信信号xL~x2の値については、ステージ処理部42L-1から出力された候補C2kの値とする。たとえば、二乗ユークリッド距離算出部401は、下記11式によって二乗ユークリッド距離dLcLk,cL1k,...c2k,c1iを算出する。二乗ユークリッド距離算出部401は、算出した二乗ユークリッド距離dLcLk,cL1k,...c2k,c1iを累積メトリック算出部402へ出力する。累積メトリック算出部402は、第Lステージまでの累積メトリックELを算出する。たとえば、累積メトリック算出部402は、下記12のようにS1mL1個の累積メトリックELを算出する。すなわち、累積メトリック算出部402は、ステージ処理部421からの累積メトリックdL1c2kと、二乗ユークリッド距離算出部401からの二乗ユークリッド距離dLcLk,cL1k,...c2k,c1iと、を加算する。累積メトリック算出部402は、算出した累積メトリックELを出力選択部403へ出力する。出力選択部403は、ステージ処理制御部410からの設定により、累積メトリック算出部402から出力された累積メトリックELをLLR算出部430へ出力する。<第L1ステージ~第Nrステージ>第L1ステージ~第Nrステージに対応するステージ処理部42L1~42Nrは、ステージ処理制御部410からの設定により、図4に示す例では動作しない。<LLRの算出>LLR算出部430は、ステージ処理部42Lから出力された累積メトリックELに基づいて、送信信号xが示す値のLLRビットLLRを算出する。LLRは、たとえば、信号が0となる信頼度情報尤度と信号が1となる信頼度情報尤度との比の対数値である。LLR算出部430は、算出したLLRを誤り訂正復号部339へ出力する。図4では、QRM-MLDアルゴリズムを例に説明したが、シンボル候補の選択アルゴリズムはQRM-MLDアルゴリズムに限られない。たとえば、シンボル候補の選択アルゴリズムとしてASESSAdaptiveSElectionofSurvivingSymbolreplicacandidates:適応生き残りシンボルレプリカ候補選択法アルゴリズムなどを用いてもよい。なお、QRM-MLDアルゴリズムについてはたとえば非特許文献K.J.KimandJ.Yue,JointchannelestimationanddatadetectionalgorithmsforMIMO-OFDMsystemsinProc.Thirty-SixthAsilomarConferenceonSignals,SystemsandComputers,pp.1857-1861,Nov.2002に記載されている。また、ASESSアルゴリズムについてはたとえば非特許文献K.Higuchi,H.Kawai,N.MaedaandM.Sawahashi,AdaptiveSelectionofSurvivingSymbolReplicaCandidatesBasedonMaximumReliabilityinQRM-MLDforOFCDMMIMOMultiplexingProc.ofIEEEGlobecom2004,pp.2480-2486,Nov.2004に記載されている。このように、実施の形態1にかかる受信機300によれば、自局宛の送信信号がMLD法で先に処理されるようにチャネル行列を生成することで、MLD法において他局宛の送信信号を処理しなくても、自局宛の送信信号の累積メトリックを算出することができる。これにより、他局宛の送信信号が示す値の候補が不明であってもMLD法を実施できるため、他局の変調方式が不明なシステムにおいても受信特性を向上させることができる。たとえば、受信機300は、自チャネルの推定値を干渉チャネルの推定値より右側に配置したチャネル行列を生成し、チャネル行列をユニタリ行列および上三角行列に分解し、送信信号ベクトルの下段側から順にメトリックを算出する。これにより、MLD法において、自局宛の送信信号を他局宛の送信信号より先に処理させることができる。これに対して、受信機300は、自チャネルの推定値を干渉チャネルの推定値より左側に配置したチャネル行列を生成してもよい。この場合は、受信機300は、チャネル行列をユニタリ行列および下三角行列に分解し、送信信号ベクトルxの上段側から順にメトリックを算出する。これにより、MLD法において、自局宛の送信信号を他局宛の送信信号より先に処理させることができる。また、たとえば上記特許文献4のように、同時にスケジュールされる全てのユーザの変調方式を同じにすることによりMLD法を使用する構成に比べて、変調方式の制限がない。このため、端末の性能や通信環境に応じた柔軟なスケジューリングが可能になるため、スループットを向上させることができる。また、推定した雑音電力に基づいてチャネル行列を拡張することで、チャネル行列に含まれる雑音成分を低減し、後段の処理における雑音強調を抑えることができる。このため、自局宛の送信信号の復号の精度を向上させ、受信特性を向上させることができる。実施の形態2実施の形態2にかかる受信機300のチャネル行列生成部333は、自チャネル推定値ベクトルhnに含まれる推定値間の順序をベクトルノルムに基づいて並べたチャネル行列Hを生成する。また、チャネル行列生成部333は、干渉チャネル推定値ベクトルhkに含まれる推定値間の順序をベクトルノルムに基づいて並べたチャネル行列Hを生成する。具体的には、自チャネル推定値ベクトルhnに含まれる各推定値を、チャネル行列Hに対して、下記13式に示すベクトルノルムPnの大きい順に右になるように配置する。これにより、自チャネル推定値ベクトルhnに含まれる各推定値は、SNRが大きいほどチャネル行列Hの右側に配置される。また、チャネル行列生成部333は、干渉チャネル推定値ベクトルhkに含まれる各推定値を、チャネル行列Hに対して、下記14式に示すベクトルノルムPkの大きい順に右になるように配置する。これにより、干渉チャネル推定値ベクトルhkに含まれる各推定値は、SNRが小さいほどチャネル行列Hの右側に配置される。これにより、チャネル行列Hに対応する送信信号ベクトルxにおいては、SNRが大きいチャネルの送信信号ほど下段に配置される。これにより、MLD処理部338において、自局宛の送信信号x1,x2,...,xLのうちのSNRが大きいチャネルの送信信号が先に処理される。このため、推定精度の高い送信信号から先に値の候補を絞り込んでいくことができるため、値の候補の絞り込みの精度が向上し、受信特性を向上させることができる。このように、実施の形態2にかかる受信機300によれば、自チャネルの各推定値および干渉チャネルの各推定値をそれぞれベクトルノルムに基づいて並べた拡張行列を生成することができる。これにより、MLD法においてSNRSignalNoiseRatio:信号雑音比が大きいチャネルの送信信号が先に処理され、推定精度の高い送信信号から先に値の候補を絞り込んでいくことができる。このため、値の候補の絞り込みの精度が向上し、受信特性を向上させることができる。実施の形態3実施の形態3にかかるチャネル行列拡張部335は、非対角成分が0で対角成分がσのL×Ntの行列をチャネル行列Hの下側に追加することによってチャネル行列Hを拡張する。具体的には、チャネル行列拡張部335は、下記15式のようにチャネル行列Hを拡張することによって拡張チャネル行列H’eを得る。ただし、ΣL×Ntは下記16式を満たす。すなわち、拡張チャネル行列H’eはNrL×Ntの行列となる。たとえば、Nt=Nr=4、L=2の場合は、拡張チャネル行列H’eは下記17式のようになる。チャネル行列拡張部335は、チャネル行列Hを拡張した拡張チャネル行列H’eをQR分解部336へ出力する。QR分解部336は、チャネル行列拡張部335から出力された拡張チャネル行列H’eをQR分解する。たとえば、QR分解部336は、下記18式に示すように、拡張チャネル行列H’eをNrL×NrLのユニタリ行列Q’eとNrL×Ntの上三角行列R’eとに分解する。QR分解部336は、QR分解により得られたユニタリ行列Q’eの左上Nr×Nrを取り出した部分行列Qのエルミート共役QHをユニタリ変換部337へ出力する。また、QR分解部336は、QR分解により得られた上三角行列R’eの上Nr行を取り出した部分行列RをMLD処理部338へ出力する。この場合は、受信信号zは、たとえば下記19式に示すように表すことができる。このように、実施の形態3にかかる受信機300によれば、チャネル行列と同じ列数および自局宛の送信信号の数に等しい行数を有し、対角成分の値が雑音電力の偏差であり、非対角成分の値がゼロである行列によってチャネル行列を拡張することができる。これにより、少ない行数の追加によりチャネル行列を拡張し、後段の計算量を低減することができる。このため、後段の処理における雑音強調を抑えるとともに、受信機300の消費電力を低減することができる。実施の形態4図5は、実施の形態4にかかる受信機の構成例を示す図である。図5において、図3に示した部分と同様の部分については同一の符号を付して説明を省略する。図3に示すように、受信機300は、図3に示した構成に加えて、SNR推定部501と、拡張制御部502と、を備えている。チャネル行列拡張部335は、拡張制御部502の制御にしたがって、チャネル行列Hの拡張の有無を切り替える。具体的には、チャネル行列拡張部335は、チャネル行列Hの拡張を行う場合は拡張により得られた拡張チャネル行列Heを出力する。また、チャネル行列拡張部335は、チャネル行列Hの拡張を行わない場合はチャネル行列Hを出力する。QR分解部336は、チャネル行列拡張部335から出力される拡張チャネル行列Heまたはチャネル行列HのQR分解を行う。自チャネル推定部331は、自チャネル推定値ベクトルhnをSNR推定部501へも出力する。干渉チャネル推定部332は、干渉チャネル推定値ベクトルhkをSNR推定部501へも出力する。雑音電力推定部334は、推定した雑音電力σ2をSNR推定部501へも出力する。SNR推定部501は、受信機300の受信信号のSNRを推定する信号雑音比推定部である。具体的には、SNR推定部501は、自チャネル推定部331からの自チャネル推定値ベクトルhnと、干渉チャネル推定部332からの干渉チャネル推定値ベクトルhkと、雑音電力推定部334からの雑音電力σ2と、に基づいてSNRを推定する。SNR推定部501は、推定したSNRを拡張制御部502へ出力する。拡張制御部502は、SNR推定部501から出力されたSNRと所定の閾値とを比較する。そして、拡張制御部502は、SNRが閾値以下である場合はチャネル行列Hの拡張を行うようにチャネル行列拡張部335を制御する。また、拡張制御部502は、SNRが閾値より大きい場合はチャネル行列Hの拡張を行わないようにチャネル行列拡張部335を制御する。図6は、SNRに対するBERの特性の一例を示す図である。図6において、横軸は受信機300の受信信号のSNRを示し、縦軸は受信機300の受信信号のBERBitErrorRate:ビット誤り率を示している。特性601は、MMSE法を用いると仮定した場合のSNRに対するBERの特性を参考として示している。特性602は、チャネル行列拡張部335においてチャネル行列Hの拡張を行う場合におけるSNRに対するBERの特性を示している。特性603は、チャネル行列拡張部335においてチャネル行列Hの拡張を行わない場合におけるSNRに対するBERの特性を示している。また、特性601~603のそれぞれは、受信信号がQPSKたとえばEnhancedVehicularAmodelであり、最大ドップラー周波数fdが70[Hz]である場合の各特性を示している。SNRが約12[dB]より低い場合では、チャネル行列Hの拡張なしの特性603はMMSE法の特性601よりBERが劣化している。一方、SNRが約12[dB]より高い場合では、チャネル行列Hの拡張なしの特性603でもMMSE法の特性601よりも優れたBER受信特性を得られる。したがって、たとえば図6に示した例において、MMSE法より良好な特性が得られれば十分である場合は、所定の閾値を12[dB]としてもよい。これにより、SNRが閾値以下である場合は、MMSE法より優れた受信特性を得つつ、チャネル行列Hの拡張を行わないことで処理量を低減することができる。このため、受信機300の消費電力を低減することができる。このように、実施の形態4にかかる受信機300によれば、SNRと閾値との比較結果に基づいてチャネル行列の拡張の有無を切り替えることで、信号の受信状態がよい場合はチャネル行列Hの拡張を行わないようにすることができる。これにより、処理量を低減することができる。実施の形態5図7は、実施の形態5にかかる受信機の構成例を示す図である。図7において、図3と同様の構成については同一の符号を付して説明を省略する。図7に示すように、実施の形態5にかかる受信機300は、図3に示した雑音電力推定部334およびチャネル行列拡張部335を省いた構成としてもよい。この場合は、チャネル行列生成部333は、生成したチャネル行列HをQR分解部336へ出力する。QR分解部336は、チャネル行列生成部333から出力されたチャネル行列HをQR分解する。この場合も、実施の形態1にかかる受信機300と同様に、自局宛の送信信号がMLD法で先に処理されるようにチャネル行列を生成することで、MLD法において他局宛の送信信号を処理しなくても、自局宛の送信信号の累積メトリックを算出することができる。これにより、他局宛の送信信号が示す値の候補が不明であってもMLD法を実施できるため、他局の変調方式が不明なシステムにおいても受信特性を向上させることができる。実施の形態6図8は、実施の形態6にかかる通信システムの一例を示す図である。図8において、図1と同様の構成については同一の符号を付して説明を省略する。図8に示すように、実施の形態6にかかる通信システム100は、複数の送信機111,112を含んでいてもよい。たとえば、送信機111,112はそれぞれ一本の送信アンテナを有しており、送信機111,112はそれぞれ受信機121,122へ信号を送信してもよい。この場合は、受信機121の受信アンテナは、受信機121を宛先とする送信機111からの送信信号だけでなく、受信機122を宛先とする送信機112からの送信信号も受信する。このように、通信システム100には、たとえば、受信機121を宛先とする送信信号を送信する第1送信アンテナと、受信機121とは異なる受信機たとえば受信機122を宛先とする送信信号を送信する第2送信アンテナと、が含まれている。第1送信アンテナおよび第2送信アンテナは、同一の送信機たとえば図1の送信機110に備えられていてもよいし、別々の送信機たとえば図8の送信機111,112に備えられていてもよい。また、第1送信アンテナおよび第2送信アンテナのそれぞれは、1つでもよいし、複数でもよい。受信機121の受信アンテナは、1つでもよいし、複数でもよい。また、受信機121の受信アンテナは、第1送信アンテナからの送信信号と、第2送信アンテナからの送信信号と、を受信する。受信機121は、受信した各信号のうちの第1送信アンテナからの送信信号を復号する。以上説明したように、受信装置および受信方法によれば、受信特性を向上させることができる。上述した各実施の形態に関し、さらに以下の付記を開示する。付記1自局宛の送信信号を送信する他局の送信アンテナと自局の受信アンテナとの間の自チャネルと、他局宛の送信信号を送信する他局の送信アンテナと前記受信アンテナとの間の干渉チャネルと、を推定する推定部と、前記推定部による前記自チャネルの推定値と前記干渉チャネルの推定値とを行方向に結合したチャネル行列を生成する生成部と、前記生成部によって生成されたチャネル行列を、直交行列と、前記自チャネルの成分が前記干渉チャネルの成分よりゼロ成分の多い行に配置された三角行列と、に分解する分解部と、前記受信アンテナによって受信された受信信号を、前記分解部によって得られた直交行列に基づいて直交化する直交化部と、前記分解部によって得られた三角行列の行ごとの成分をゼロ成分のより多い行から順に取得し、前記行ごとの成分を取得するごとに、取得した行ごとの成分と前記直交化部によって直交化された受信信号とに基づいて前記自局宛の送信信号の値のメトリックを累積的に算出する算出部と、前記算出部によって算出されたメトリックに基づいて前記自局宛の送信信号を復号する復号部と、を備えることを特徴とする受信装置。付記2前記算出部は、前記行ごとの成分を、前記自局宛の送信信号を送信する他局の送信アンテナの数だけ取得することを特徴とする付記1に記載の受信装置。付記3前記生成部は、前記自チャネルの推定値を前記干渉チャネルの推定値より右側に配置したチャネル行列を生成し、前記分解部は、前記チャネル行列を直交行列および上三角行列に分解し、前記算出部は、前記行ごとの成分を前記上三角行列の下段から順に取得することを特徴とする付記1または2に記載の受信装置。付記4前記生成部は、前記自チャネルの推定値を前記干渉チャネルの推定値より左側に配置したチャネル行列を生成し、前記分解部は、前記チャネル行列を直交行列および下三角行列に分解し、前記算出部は、前記行ごとの成分を前記下三角行列の上段から順に取得することを特徴とする付記1または2に記載の受信装置。付記5前記推定部による各推定値に基づく雑音電力を推定する雑音推定部と、前記生成部によって生成されたチャネル行列の列数に等しい行数および列数の単位行列に前記雑音推定部によって推定された雑音電力の偏差を乗算した行列を、前記チャネル行列の下側に追加して前記チャネル行列を拡張する拡張部と、を備え、前記分解部は、前記拡張部によって拡張されたチャネル行列を分解することを特徴とする付記1~4のいずれか一つに記載の受信装置。付記6前記推定部は、自局宛の送信信号を送信する他局の複数の送信アンテナと前記受信アンテナとの間の各自チャネルを推定し、他局宛の送信信号を送信する他局の複数の送信アンテナと前記受信アンテナとの間の各干渉チャネルを推定し、前記生成部は、前記自チャネルの各推定値を前記自チャネルの各推定値のベクトルノルムに基づいて並べ、前記干渉チャネルの各推定値を前記干渉チャネルの各推定値のベクトルノルムに基づいて並べたチャネル行列を生成することを特徴とする付記1~5のいずれか一つに記載の受信装置。付記7前記推定部による各推定値に基づく雑音電力を推定する雑音推定部と、前記生成部によって生成されたチャネル行列と同じ列数および前記自局宛の送信信号の数に等しい行数を有し、対角成分の値が前記雑音推定部によって推定された雑音電力の偏差であり非対角成分の値がゼロである行列を、前記チャネル行列の下側に追加して前記チャネル行列を拡張する拡張部と、を備え、前記分解部は、前記拡張部によって拡張されたチャネル行列を分解することを特徴とする付記1~4のいずれか一つに記載の受信装置。付記8前記推定部による各推定値に基づく信号雑音比を推定する信号雑音比推定部を備え、前記拡張部は、前記信号雑音比推定部によって推定された信号雑音比と閾値との比較結果に基づいて前記チャネル行列の拡張の有無を切り替え、前記分解部は、前記拡張部によって拡張されたチャネル行列、または前記拡張部によって拡張されなかったチャネル行列を分解することを特徴とする付記5または7に記載の受信装置。付記9前記算出部は、前記他局宛の送信信号が示す値のメトリックについては算出しないことを特徴とする付記1~8のいずれか一つに記載の受信装置。付記10自局宛の送信信号を送信する他局の送信アンテナと自局の受信アンテナとの間の自チャネルと、他局宛の送信信号を送信する他局の送信アンテナと前記受信アンテナとの間の干渉チャネルと、を推定し、前記自チャネルの推定値と、前記干渉チャネルの推定値と、を行方向に結合したチャネル行列を生成し、生成した前記チャネル行列を、直交行列と、前記自チャネルの成分が前記干渉チャネルの成分よりゼロ成分の多い行に配置された三角行列と、に分解し、前記受信アンテナによって受信された受信信号を、前記直交行列に基づいて直交化し、前記三角行列の行ごとの成分をゼロ成分のより多い行から順に取得し、前記行ごとの成分を取得するごとに、取得した行ごとの成分と直交化した前記受信信号とに基づいて前記自局宛の送信信号の値のメトリックを累積的に算出し、算出した前記メトリックに基づいて前記自局宛の送信信号を復号する、ことを特徴とする受信方法。"