"複数のアンテナを用いて無線通信を行う無線装置において、受信した受信信号点から、取り得る送信信号点に伝搬路状態の影響を反映させたレプリカを減算する減算部と、2次関数を2のべき乗で区切った線形補間式を用いて、前記減算部によって算出される値から、前記受信信号点と前記レプリカとのメトリックを計算する計算部と、を有することを特徴とする無線装置。前記計算部は、前記減算部によって算出される値の2進数における桁数を検出する検出部を有し、前記線形補間式を用いて、前記減算部によって算出される値と、前記検出部で検出される桁数とから、前記受信信号点と前記レプリカとのメトリックを計算することを特徴とする請求項1記載の無線装置。前記計算部は、前記減算部によって算出される値から、前記検出部によって検出される桁数から1を引いた数のビットを抽出する抽出部をさらに有し、前記線形補間式を用いて、前記減算部によって算出される値、前記検出部で検出される桁数、および前記抽出部から抽出されるビットから、前記受信信号点と前記レプリカとのメトリックを計算することを特徴とする請求項2記載の無線装置。前記線形補間式は、前記減算部によって算出される値をx、前記xの桁数をnとして、2n13|x|-2nで示され、前記計算部は、前記xの|x|を算出する絶対値算出部と、前記|x|をビットシフトして2|x|を算出する第1のシフト部と、前記|x|と前記2|x|とを加算して3|x|を算出する加算部と、前記|x|から前記nを検出する検出部と、前記nから2nを出力する第2のシフト部と、前記3|x|から前記2nを減算して3|x|-2nを算出する減算部と、前記3|x|-2nをビットシフトして2n13|x|-2nを算出する第3のシフト部と、を有することを特徴とする請求項1記載の無線装置。前記線形補間式は、前記減算部によって算出される値をx、前記xの桁数をnとして、2n1{|x|2|x|-2n1}で示され、前記計算部は、前記xの|x|を算出する絶対値算出部と、前記|x|から前記nを検出する検出部と、前記|x|から下位n-1ビットを抽出して|x|-2n1を出力する抽出部と、前記|x|-2n1をビットシフトして2|x|-2n1を出力する第1のシフト部と、前記|x|と前記2|x|-2n1とを加算して|x|2|x|-2n1を算出する加算部と、前記|x|2|x|-2n1をビットシフトして2n1{|x|2|x|-2n1}を出力する第2のシフト部と、を有することを特徴とする請求項1記載の無線装置。複数のアンテナを用いて無線通信を行う無線装置のメトリック計算方法において、受信した受信信号点から、取り得る送信信号点に伝搬路状態の影響を反映させたレプリカを減算し、2次関数を2のべき乗で区切った線形補間式を用いて、前記受信信号点から前記レプリカを減算した値から、前記受信信号点と前記レプリカとのメトリックを計算する、ことを特徴とするメトリック計算方法。しかし、2乗ユークリッド距離の計算は、2乗計算を行うために乗算器を要し、多数のレプリカ検索のため大量の乗算処理を要する。そのため、回路規模が増大し、消費電力が大きくなるという問題点があった。本件はこのような点に鑑みてなされたものであり、回路規模の増大を抑制し、消費電力を低減することができる無線装置およびメトリック計算方法を提供することを目的とする。上記課題を解決するために、複数のアンテナを用いて無線通信を行う無線装置が提供される。この無線装置は、受信した受信信号点から、取り得る送信信号点に伝搬路状態の影響を反映させたレプリカを減算する減算部と、2次関数を2のべき乗で区切った線形補間式を用いて、前記減算部によって算出される値から、前記受信信号点と前記レプリカとのメトリックを計算する計算部と、を有する。開示の装置および方法によれば、回路規模の増大を抑制し、消費電力を低減することができる。以下、実施の形態を、図面を参照して詳細に説明する。[第1の実施の形態]図1は、第1の実施の形態に係る無線装置のブロック図である。図1に示すように、無線装置1は、減算部1aおよび計算部1bを有している。図1には、無線装置1にデータを送信する送信装置2も示している。無線装置1および送信装置2は、それぞれ複数のアンテナを有し、MIMOにより無線通信を行う。無線装置1の減算部1aは、送信装置2から受信した受信信号点から、取り得る送信信号点に伝搬路状態の影響を反映させたレプリカを減算する。受信信号点およびレプリカは、IQ座標により、複素数で示される。従って、受信した受信信号点からレプリカを減算するということは、レプリカに対する受信信号点のベクトルを算出しているということになる。減算部1aからは、次の式1に示す形式の値が得られる。なお、jは、虚数単位を示す。ajb...1計算部1bは、2次関数を2のべき乗で区切った線形補間式を用いて、減算部1aによって算出される、受信信号点からレプリカを減算した値から、受信信号点とレプリカとのメトリックを計算する。ここで、受信信号点とレプリカとのメトリックを2乗ユークリッド距離で求めると、次の式2で示される。a2b2...2一方、計算部1bは、上記したように2次関数を2のべき乗で区切った線形補間式を用いてメトリックを計算する。例えば、2次関数fx=x2を2のべき乗で区切った線形補間式は、次の式3によって示される。なお、nは、xを2進数で示したときの桁数を示している。fx=2n13|x|-2n...3従って、計算部1bの計算するメトリックは、次の式4で示される。2k13|a|-2k2l13|b|-2l...4kは、aを2進数で示したときの桁数を示している。lは、bを2進数で示したときの桁数を示している。なお、式4の第1項は、式2の第1項の近似計算を行っている。式4の第2項は、式2の第2項の近似計算を行っている。例えば、減算部1aから出力される式1のa,bは、16ビットで出力されるとする。a,bのそれぞれの値は、000101...101、0000011...001であるとする。前記の2進数の左側をMSBMostSignificantBit、右側をLSBLeastSignificantBitとする。なお、以下で2進数を示した場合、左側をMSB、右側をLSBとする。この場合、aの桁数は、MSB16ビット目から14ビット目までが0であるので、13桁となる。bの桁数は、MSBから12ビット目までが0であるので、11桁となる。従って、式4のk,lは、k=13およびl=11となり、計算部1bの計算するメトリックは、次の式5に示すようになる。2123|a|-2132103|b|-211...5なお、aは、上記した000101...101であり、bは0000011...001である。すなわち、計算部1bは、2次関数を2のべき乗で区切った線形補間式を用いて、2乗ユークリッド距離を近似計算する。2次関数を2のべき乗で区切った線形補間式について説明する。図2は、2次関数を2のべき乗で区切った線形補間式を説明する図である。図2に示す点線のグラフW1は、fx=x2のグラフを示している。実線のグラフW2~W4は、グラフW1を2のべき乗2n1≦x<2nで区切った線形補間式のグラフを示している。例えば、グラフW2は、グラフW1を28と29で区切った区間の線形補間式を示している。グラフW3は、グラフW1を29と210で区切った区間の線形補間式を示している。グラフW4は、グラフW1を210と211で区切った区間の線形補間式を示している。なお、図2では、0~21,21~22,...,27~28および211~212,...で区切った線形補間式も存在するが、その図示を省略している。式3の桁数nは、2次関数を2のべき乗2n1≦x<2nで区切った線形補間式を選択していると言える。例えば、n=9の場合、線形補間式は、グラフW2に示すようになる。n=10の場合、線形補間式は、グラフW3に示すようになる。n=11の場合、線形補間式は、グラフW4に示すようになる。例えば、520の2乗を、式3の線形補間式で求める近似するとする。520は、2進数で示すと、‘1000001000’である。桁数nは、10桁である。従って、520を求める線形補間式は、グラフW3n=10となり、次の式6で示される。fx=293|x|-210...6xに520を代入すると、式6の値は、274432となる。なお、520の2乗は、270400である。2次関数を2のべき乗で区切った線形補間式は、式3で示される。式3に示す3|x|の乗算は、2|x||x|に分解でき、2|x|は、|x|を1ビットシフトすることにより演算できる。また、2n1と3|x|-2nの乗算は、3|x|-2nをn-1ビットシフトすることにより演算できる。すなわち、式3の線形補間式は、乗算器をビットシフトに置き換えることによってメトリックを計算できる。つまり、無線装置1は、乗算器を用いないで2乗ユークリッド距離の精度のよい近似値を算出でき、回路規模の増大を抑制し、消費電力を低減することができる。このように、無線装置1の減算部1aは、受信した受信信号点からレプリカを減算する。そして、計算部1bは、2次関数を2のべき乗で区切った線形補間式を用いて、減算部1aによって算出される値から、受信信号点とレプリカとのメトリックを計算するようにした。これにより、無線装置1は、回路規模の増大を抑制し、消費電力を低減することができる。また、無線装置1は、2次関数を2のべき乗で区切った線形補間式を用いてメトリックを算出するので、精度のよいメトリックを算出できる。[第2の実施の形態]次に、第2の実施の形態を、図面を参照して詳細に説明する。図3は、第2の実施の形態に係る移動通信システムを示した図である。図3に示すように、移動通信システムは、無線装置10および送信装置20を有している。無線装置10および送信装置20は、それぞれ複数のアンテナを有し、MIMOにより無線通信を行う。無線装置10は、例えば、携帯電話機などの無線端末である。送信装置20は、例えば、基地局である。無線装置10および送信装置20は、例えば、LTEによる無線通信を行う。図4は、送信装置のブロック図である。図4に示すように、送信装置20は、データ分割部21、シンボル変調部22a,22b、パイロット多重部23a,23b、無線部24a,24b、およびアンテナ25a,25bを有している。データ分割部21には、無線装置10に無線送信するデータが入力される。データ分割部21は、入力されたデータを分割し、シンボル変調部22a,22bに出力する。例えば、データ分割部21には、データD1,D2,D3,D4,...が入力されるとする。データ分割部21は、データD1,D3,...をシンボル変調部22aに出力し、データD2,D4,...をシンボル変調部22bに出力する。シンボル変調部22a,22bは、データ分割部21から出力されるデータをシンボル変調し、パイロット多重部23a,23bに出力する。例えば、シンボル変調部22a,22bは、データ分割部21から出力されるデータを、QPSKQuadraturePhaseShiftKeyingやM値のQAMQuadratureAmplitudeModulationによってシンボル変調する。パイロット多重部23a,23bは、シンボル変調部22a,22bから出力されるシンボルに、パイロット信号パイロットシンボルを多重する。無線装置10は、データシンボルに多重されるパイロット信号によって、電波伝搬路のチャネル状態を推定することができる。無線部24a,24bは、パイロット多重部23a,23bから出力されるシンボルを無線周波数の無線信号に変換し、アンテナ25a,25bに出力する。なお、図4では、送信装置20は、2本のアンテナ25a,25bで無線信号を送信する例を示しているが、アンテナは、3本以上であってもよい。例えば、アンテナがn本の場合、シンボル変調部、パイロット多重部、および無線部は、n本のアンテナに対応して、n個ずつ存在する。データ分割部21は、入力されるデータをn個に分割して、n個のシンボル変調部に出力する。図5は、無線装置のブロック図である。図5に示すように、無線装置10は、アンテナ11a,11b、無線部12a,12b、シンボル分離部13a,13b、チャネル推定部14、およびMLD処理部15を有している。無線部12a,12bは、アンテナ11a,11bで受信された無線信号の周波数をベースバンド周波数に変換し、ベースバンド信号を抽出する。シンボル分離部13a,13bは、無線部12a,12bから出力されるベースバンド信号から、データのシンボルと、パイロット信号のシンボルとに分離する。シンボル分離部13a,13bは、分離したデータのシンボルをMLD処理部15に出力し、分離したパイロット信号のシンボルをチャネル推定部14に出力する。シンボル分離部13a,13bから出力されるデータのシンボル受信信号点は、伝達関数行列Hの要素をHijとして、次の式7によって表される。iは、無線装置10のアンテナ11a,11bに対応する。アンテナ11aをi=1、アンテナ11bをi=2とする。jは、アンテナ25a,25bに対応する。アンテナ25aをj=1、アンテナ25bをj=2とする。Ntxは、送信側のアンテナの総数を示す。送信装置20は、2本のアンテナ25a,25bを備えているので、Ntxは2となる。sjは、送信装置20のアンテナjから送信された信号点を示す。niは、アンテナiで受信した信号に含まれる外部雑音や熱雑音などのノイズを示す。チャネル推定部14は、シンボル分離部13a,13bから出力されるパイロット信号を用いて、アンテナ11a,11bごとに、送信装置20のアンテナ25a,25bごとにおける伝搬路のチャネル値を算出する。すなわち、チャネル推定部14は、伝達関数行列Hの各成分を算出する。例えば、チャネル推定部14は、アンテナ25aとアンテナ11aとにおける伝搬路のチャネル値と、アンテナ25bとアンテナ11aとにおける伝搬路のチャネル値とを算出する。また、チャネル推定部14は、アンテナ25aとアンテナ11bとにおける伝搬路のチャネル値と、アンテナ25bとアンテナ11bとにおける伝搬路のチャネル値とを算出する。MLD処理部15には、チャネル推定部14で算出されたチャネル値と、シンボル分離部13a,13bから出力されるデータのシンボルとが入力される。MLD処理部15は、MLDによって、チャネル推定部14で算出されたチャネル値と、シンボル分離部13a,13bから出力されるデータのシンボルとから、送信装置20が送信したデータを復調する。なお、図5では、無線装置10は、2本のアンテナ11a,11bで無線信号を受信する例を示しているが、アンテナは、3本以上であってもよい。例えば、アンテナがn本の場合、無線部およびシンボル分離部は、n本のアンテナに対応して、n個ずつ存在する。以下では、説明を簡単にするために、送信装置20が2本のアンテナ25a,25bを備え、無線装置10が2本のアンテナ11a,11bを備えた2×2MIMOの例について説明するが、もちろん、n×mMIMOについても適用できる。図6は、MLD処理部のブロック図である。図6に示すように、MLD処理部15は、レプリカ生成部15a、減算部15b、メトリック計算部15c、および最小メトリック候補選択部15dを有している。レプリカ生成部15aには、チャネル推定部14で算出されたチャネル値が入力される。レプリカ生成部15aは、チャネル推定部14で算出されたチャネル値を考慮したレプリカを生成する。例えば、レプリカ生成部15aは、送信装置20から送信され得る複数の送信信号点に伝達関数行列Hを乗算したレプリカを生成する。すなわち、レプリカ生成部15aは、伝搬路の状態に応じて、送信装置20から送信され得る複数の送信信号点の振幅を変更または回転してレプリカを生成する。具体的には、送信装置20のアンテナ25a,25bから送信され得る信号点をsjkとする。jは、上記したように、送信装置20のアンテナ25a,25bに対応する。kは、アンテナjに適用されている変調方式によって決まる値である。例えば、アンテナjにQPSKが適用されている場合、kは、k=1,2,...,4の値をとり、アンテナjから送信され得る信号点sjkは、IQ軸上のそれぞれにおいて次の値をとる。また、例えば、アンテナjに16QAMが適用されている場合、kは、k=1,2,...,16の値をとり、アンテナjから送信され得る信号点sjkは、IQ軸上のそれぞれにおいて次の値をとる。また、アンテナjに64QAMが適用されている場合、kは、k=1,2,...,64の値をとり、アンテナjから送信され得る信号点sjkは、IQ軸上のそれぞれにおいて次の値をとる。レプリカ生成部15aは、sjkに伝達関数行列Hを乗算してレプリカを生成する。例えば、アンテナjに64QAMの変調方式が適用されているとすると、レプリカ生成部15aは、k=1,2,...,64の64個の信号点sjkに対して、Hijsjkを計算する。Hijは、伝達関数行列Hの要素を示す。iは、上記したように、無線装置10のアンテナ11a,11bに対応する。jは、上記したように、送信装置20のアンテナ25a,25bに対応する。ここで、上記のa~cに示す各変調方式における信号点sjkを、それぞれ、212,1012,4212倍することを考える。すると、各変調方式のIQ軸上のそれぞれの値実数部および虚数部の値は、次の値をとる。QPSK:±1...d16QAM:±1,±3...e64QAM:±1,±3,±5,±7...f従って、レプリカ生成処理は、Hijsjkの乗算を含むが、成分に分解すれば、乗算回数を低減することができる。例えば、64QAMにおいて、Hijsjkの実数部は、次の式8によって算出することができるReHijsjk=ReHijResjk-ImHijImsjk...8式8の右辺第1項は、次の式9に示すように変形することができる。式9の最右辺の4212sjkは、上記のfの値をとる。従って、レプリカ生成部15aは、±A,±3A,±5A,±7Aを計算することにより、式8の右辺第1項を算出できる。3A,5A,7Aは、それぞれ、3A=2AA,5A=4AA,7A=8A-Aと分解できる。2A,4A,8Aの乗算は、ビットシフトで実現できる。また、-A,-3A,-5A,-7Aは、A,3A,5A,Aを符号反転することにより算出できる。従って、レプリカ生成部15aは、例えば、上記のfに示す定数をメモリなどの記憶装置に記憶し、Aの計算を1回行えば、後は、ビットシフト、加算処理、および符号反転により、式8の右辺第1項を計算できる。つまり、レプリカ生成部15aは、1回の乗算で式8の右辺第1項を計算できる。なお、レプリカ生成部15aは、式8の右辺第2項も上記と同様にして計算することができる。また、レプリカ生成部15aは、次の式10によって、64QAMにおけるHijsjkの虚数部を算出できる。ImHijsjk=ReHijImsjkImHijResjk...10レプリカ生成部15aは、上記の実数部の場合と同様にして式10の値を算出できる。また、上記では、64QAMのレプリカ生成について説明したが、その他の変調方式の場合も同様にしてレプリカを生成することができる。また、上記は、送信装置20から送信され得る信号点が均等な確率で発生する場合に、sjkの電力が1となるような振幅で送信されているとして説明した。sjkの電力が1でない場合、レプリカ生成部15aは、対応する振幅で規格化することにより、上記と同様にしてレプリカを生成することができる。減算部15bは、シンボル分離部13a,13bから出力される受信信号点から、レプリカ生成部15aで生成されたレプリカを減算する。受信信号点からのレプリカの減算は、次の式11で示される。式11は、複素数であり、次の式12の形式で示される。xijyi...12xi,yiは、s1k1,s2k2のバリエーションについて計算される。例えば、s1k1,s2k2がともに64QAMで変調されているとすると、s1k1,s2k2のそれぞれのk1,k2は、64通りあるので、s1k,s2kの組み合わせは、64×64通りある。なお、この計算は、例えば、QRM-MLDやシンボルランキングに基づく方法などによって、検索対象となるレプリカの組み合わせを絞り込み、計算数を低減することができる。すなわち、減算部15bは、無線装置10のアンテナiにおけるxi,yiを算出してメトリック計算部15cに出力する。メトリック計算部15cは、2次関数を2のべき乗で区切った線形補間式を用いて、減算部15bによって算出される値xi,yiから、シンボル受信信号点とレプリカとのメトリックを計算する。ここで、受信信号点とレプリカとのメトリックを2乗ユークリッド距離で求めると、次の式13で示される。式13に示すように、2乗ユークリッド距離では乗算を要する。すなわち、2ユークリッド距離を求めるには、例えば、xiとxiとを乗算する乗算器を要する。一方、メトリック計算部15cは、上記したように、2次関数を2のべき乗で区切った線形補間式を用いて、減算部15bによって算出される値xi,yiから、受信信号点とレプリカとのメトリックを計算する。従って、以下で詳細に説明するが、メトリック計算部15cは、乗算器を用いることなく、受信信号点とレプリカとのメトリックを計算することができる。最小メトリック候補選択部15dは、メトリック計算部15cの計算するメトリックが最小となる送信信号点s1k1,s2k2の組み合わせを選択する。図7は、メトリック計算部のブロック図である。図7に示すように、メトリック計算部15cは、演算部30,40および加算部50を有している。演算部30には、減算部15bから出力されるxiが入力され、演算部40には、減算部15bから出力されるyiが入力される。以下では、xi,yiを単にx,yと示すことがある。演算部30は、2次関数を2のべき乗で区切った線形補間式によって、x2の近似値を算出する。演算部40は、2次関数を2のべき乗で区切った線形補間式によって、y2の近似値を算出する。2次関数を2のべき乗で区切った線形補間式は、例えば、式3で示される。加算部50は、演算部30,40から出力されるx2,y2の近似値を加算する。つまり、加算部50は、式13に示す2乗ユークリッド距離の近似値を出力する。演算部30は、絶対値算出部31、シフト部32,35,37、加算部33、桁検出部34、および減算部36を有している。演算部30は、前記の各部によって、式3の演算を行う。絶対値算出部31には、減算部15bから出力されるxが入力される。絶対値算出部31は、入力されるxの絶対値|x|を算出する。シフト部32は、絶対値算出部31から出力される|x|を、2倍するように1ビットシフトする。加算部33には、絶対値算出部31から出力される|x|と、シフト部32から出力される2|x|とが入力される。加算部33は、入力される|x|と2|x|とを加算する。すなわち、加算部33は、式3の3|x|を算出する。桁検出部34は、絶対値算出部31から出力される|x|の桁数を検出する。桁検出部34は、例えば、|x|の桁数を、|x|のビット幅のビットフラグで出力する。例えば、桁検出部34は、|x|の桁数に対応するビットを立てた、|x|のビット幅のビットフラグを出力する。具体的には、絶対値算出部31から出力される|x|を、16ビットのデータとする。桁検出部34は、16ビットデータ|x|の桁数を示す16ビット幅のビットフラグを出力する。この桁検出部34から出力される16ビット幅のビットフラグの1ビット目は、16ビットデータ|x|の1桁目に対応し、桁検出部34から出力される16ビット幅のビットフラグの2ビット目は、16ビットデータ|x|の2桁目に対応し、以下、同様にして、桁検出部34から出力される16ビット幅のビットフラグの16ビット目は、16ビットデータ|x|の16桁目に対応するとする。桁検出部34に入力される|x|を、例えば、000101...10116ビットとする。この場合、桁検出部34からは、13ビット目に1が立ち、その他のビットは0の16ビット幅のビットフラグが出力される。シフト部35は、桁検出部34から出力される、|x|の桁数を示すビットフラグを2倍するように1ビットシフトする。桁検出部34からは、各ビットが|x|の桁数に対応したビットフラグで出力される。従って、ビットフラグの値を2倍するということは、例えば、|x|の桁数がnであるとすると、シフト部35からは、2nの値が出力されることになる。例えば、上記例に従い、桁検出部34からは、13ビット目に1が立ったビットフラグが出力されたとする。この値を2のべき乗で示すと、212となる。これを2倍するようにビットシフトすると、213となる。すなわち、|x|の桁数をnとすると、桁検出部34からは、2n1の値が出力される。シフト部35は、これを2倍し、2nの値を出力する。つまり、シフト部35は、式3の右辺の括弧内の2nを算出する。減算部36は、加算部33から出力される3|x|から、シフト部35から出力される2nを減算する。すなわち、減算部36は、式3の‘3|x|-2n’を算出する。シフト部37には、桁検出部34から出力される|x|の桁数と、減算部36から出力される‘3|x|-2n’とが入力される。シフト部37は、減算部36から出力される‘3|x|-2n’が、2の‘桁数-1’乗倍となるようにビットシフトする。すなわち、|x|の桁数をnとすると、シフト部37は、‘3|x|-2n’を2n1倍する。つまり、シフト部37は、式3の右辺を算出する。以上から、演算部30からは、次の式14に示す値が出力される。2n13|x|-2n...14なお、nは、|x|の桁数を示す。演算部40は、演算部30と同様の部位を有している。従って、演算部40からは、次の式15に示す値が出力される。2m13|y|-2m...15なお、mは、|y|の桁数を示す。加算部50は、演算部30から出力される値と、演算部40から出力される値とを加算する。次の式16は、加算部50から出力される値を示している。2n13|x|-2n2m13|y|-2m...16式16の第1項は、式13の|x|2に対応し、式16の第2項は、式13の|y|2に対応する。このようにして、メトリック計算部15cからは、受信信号点とレプリカとのメトリックが出力される。図8は、桁検出部のブロック図である。図8に示すように、桁検出部34は、OR部61a,61b,...,61l,61m,61nおよびXOR部62a,62b,...,62l,62m,62nを有している。OR部61aには、|x|のMSBと、それより1ビット桁の小さいビットとが入力される。OR部61aは、|x|のMSBと、それより1ビット桁の小さいビットとのOR演算を行う。XOR部62aには、|x|のMSBと、OR部61aのOR演算結果とが入力される。XOR部62aは、|x|のMSBと、OR部61aのOR演算結果とのXOR演算を行う。OR部61bには、OR部61aのOR演算結果と、|x|のMSBより2ビット桁の小さいビットとが入力される。OR部61bは、OR部61aのOR演算結果と、|x|のMSBより2ビット桁の小さいビットとのOR演算を行う。XOR部62bには、OR部61aのOR演算結果と、OR部61bのOR演算結果とが入力される。XOR部62bは、OR部61aのOR演算結果と、OR部61bのOR演算結果とのXOR演算を行う。以下同様にして、残りのOR部およびXOR部もOR演算およびXOR演算を行う。なお、OR部61nには、|x|のLSBが入力される。|x|のMSBからあるビットで初めて1が立つと、そのビットに対応するOR部は、1を出力する。そして、そのビットに対応するXORも1を出力する。1が立ったビットより下位のビットに対応するOR部は、上位のビットに対応するOR部から1が入力されるので、入力される|x|のビットが0であっても1であっても1を出力する。1が立ったビットより下位のビットに対応するXOR部は、対応するOR部から1が入力され、また、上位のビットに対応するOR部から1が入力されるので0を出力する。すなわち、桁検出部34は、|x|のMSBからあるビットで初めて1が立つと1を出力し、それより下位のビットでは0を出力する。なお、桁検出部34は、|x|のMSBが1である場合、1をそのまま出力し、MSBより下位のビットでは0を出力する。図7に示すメトリック計算部15cは、絶対値算出部31、シフト部32,35,37、加算部33、および減算部36を有し、乗算器を含んでない。また、桁検出部34は、図8に示すように、乗算器を含むことなく形成することができる。すなわち、メトリック計算部15cは、乗算部を備えることなく、受信信号点とレプリカとのメトリックを計算することができる。よって、メトリック計算部15cは、回路規模の増大を抑制し、消費電力を低減することができる。図9は、線形補間式、2乗ユークリッド距離、およびマンハッタン距離でメトリックを算出した場合に対する誤り率のシミュレーション結果を示した図である。図9の横軸はSNRSignalNoiseRatioを示し、縦軸は無線装置10の受信したデータの誤り率を示している。図9に示す実線の波形W11は、2乗ユークリッド距離でメトリックを算出した場合の、データの誤り率のシミュレーション結果を示している。1点鎖線の波形W12は、2次関数を2のべき乗で区切った線形補間式を用いてメトリックを算出した場合の、データの誤り率のシミュレーション結果を示している。2点鎖線の波形W13は、マンハッタン距離を用いてメトリックを算出した場合の、データの誤り率のシミュレーション結果を示している。波形W13に示すように、マンハッタン距離でメトリックを算出した場合は、2乗ユークリッド距離でメトリックを算出した場合より誤り率が劣化している。一方、波形W12に示すように、2次関数を2のべき乗で区切った線形補間式を用いてメトリックを算出した場合は、マンハッタン距離より誤り率がよく、2乗ユークリッド距離と同程度の誤り率を得ている。図10は、無線装置の動作を示したフローチャートである。[ステップS1]レプリカ生成部15aは、チャネル推定部14で算出されたチャネル値と、送信装置20から送信され得る複数の送信信号点とに基づいてレプリカを生成する。[ステップS2]減算部15bは、受信された受信信号点から、レプリカ生成部15aで生成されたレプリカを減算する。[ステップS3]メトリック計算部15cは、2次関数を2のべき乗で区切った線形補間式を用いて、減算部15bによって算出される値から、受信信号点とレプリカとのメトリックを計算する。[ステップS4]最小メトリック候補選択部15dは、メトリック計算部15cの計算するメトリックが最小となる送信信号点の組み合わせを選択する。このように、無線装置10の減算部15bは、受信した受信信号点からレプリカを減算する。そして、メトリック計算部15cは、2次関数を2のべき乗で区切った線形補間式を用いて、減算部15bによって算出された値から、受信信号点とレプリカとのメトリックを計算するようにした。これにより、無線装置10は、回路規模の増大を抑制し、消費電力を低減することができる。また、メトリック計算部15cは、2次関数を2のべき乗で区切った線形補間式を用いてメトリックを算出するので、2乗ユークリッド距離と同程度の誤り率を得ることができる。[第3の実施の形態]次に、第3の実施の形態を、図面を参照して詳細に説明する。第3の実施の形態では、メトリック計算部の別の例について説明する。なお、第3の実施の形態に係る無線装置のブロックは、図5に示したブロックと同様である。また、MLD処理部15は、図6に示したブロックと同様である。式3に示す2次関数を2のべき乗で区切った線形補間式は、次の式17に示すように変形することができる。fx=2n13|x|-2n=2n1|x|2|x|-2n=2n1{|x|2|x|-2n1}...17式17の第2項‘|x|-2n1’は、2n1≦x<2nであることから、|x|を2進数で示した場合の下位n-1ビットを抽出したものである。例えば、|x|は、16ビットのデータで、00...0010010100であるとする。前記データのMSBから9ビット目までは0とする。|x|の桁数は、8桁n=8である。この場合、‘|x|-2n1’の値は、下位7ビットを抽出した10100となる7,6ビットは0。第3の実施の形態のメトリック計算部15cは、式17を計算する。図11は、第3の実施の形態に係るメトリック計算部のブロック図である。図11に示すように、メトリック計算部15cは、演算部70,80および加算部90を有している。演算部70には、減算部15bから出力されるxが入力され、演算部80には、減算部15bから出力されるyが入力される。演算部70は、2次関数を2のべき乗で区切った線形補間式によって、x2の近似値を算出する。演算部70は、2次関数を2のべき乗で区切った線形補間式によって、y2の近似値を算出する。2次関数を2のべき乗で区切った線形補間式は、式17で示される。加算部90は、演算部70,80から出力されるx2,y2の近似値を加算する。つまり、加算部50は、式13に示す2乗ユークリッド距離の近似値を出力する。演算部70は、絶対値算出部71、桁検出部72、ビット抽出部73、シフト部74,76、および加算部75を有している。絶対値算出部71には、減算部15bから出力されるxが入力される。絶対値算出部71は、入力されるxの絶対値|x|を算出する。桁検出部72は、絶対値算出部71から出力される|x|の桁数を検出する。桁検出部72は、例えば、|x|の桁数を、|x|のビット幅のビットフラグで出力する。例えば、桁検出部72は、|x|の桁数に対応するビットを立てた、|x|のビット幅のビットフラグを出力する。桁検出部72は、図8に示したブロックと同様のブロックを有している。ビット抽出部73は、|x|の下位n-1ビットを抽出する。nは、|x|の桁数を示す。例えば、|x|の桁数がnであるとすると、桁検出部72からは、2n1の値が出力されることになる。ビット抽出部73は、例えば、|x|と2n1とのXORを演算することにより、|x|の下位n-1ビットを抽出することができる。なお、|x|の下位n-1ビットを抽出した値は、|x|-2n1であり、ビット抽出部73は、式17の|x|-2n1を算出している。シフト部74は、ビット抽出部73から出力される‘|x|-2n1’の値が2倍となるように1ビットシフトする。すなわち、シフト部74は、式17の2|x|-2n1を算出する。加算部75は、絶対値算出部71から出力される|x|と、シフト部74から出力される2|x|-2n1とを加算する。すなわち、加算部75は、式17の|x|2|x|-2n1を算出する。シフト部76は、加算部75から出力される値‘|x|2|x|-2n1’が2n1倍となるようにn-1ビットシフトする。すなわち、シフト部76は、式17の‘2n1{|x|2|x|-2n1}’を算出する。演算部80は、演算部70と同様の部位を有している。従って、演算部80からは、次の式18に示す値が出力される。2m1{|y|2|y|-2m1}...18なお、mは、|y|の桁数を示す。加算部90は、演算部70から出力される値と、演算部80から出力される値とを加算する。次の式19は、加算部90の出力の値を示している。2n1{|x|2|x|-2n1}2m1{|y|2|y|-2m1}...19式19の第1項は、式13の|x|2に対応し、式19の第2項は、式13の|y|2に対応する。このようにして、メトリック計算部15cからは、受信信号点とレプリカとのメトリックが出力される。図11に示すメトリック計算部15cは、図7に示すメトリック計算部15cに対し、減算部36が不要である。従って、図11に示すメトリック計算部15cは、図7に示すメトリック計算部15cよりさらに回路規模の増大を抑制し、消費電力を低減することができる。このように、無線装置10のビット抽出部73は、減算部15bによって算出される値から、桁検出部72によって検出される桁数から1を引いた数のビットを抽出する。そして、メトリック計算部15cは、線形補間式を用いて、減算部15bによって算出される値、桁検出部72で検出される桁数、およびビット抽出部73から抽出されるビットから、受信信号点とレプリカとのメトリックを計算するようにした。これにより、無線装置10は、回路規模の増大を抑制し、消費電力を低減することができる。以下、無線装置10のハードウェア構成例について説明する。図12は、無線装置のハードウェア構成例を示した図である。図12に示すように、無線装置10は、RFRadioFrequency部101a,101b、ADCAnalogtoDigitalConverter102a,102b、およびLSILargeScaleIntegration103を有している。RF部101a,101bは、例えば、1チップの半導体装置で形成されている。ADC102a,102bは、RF部101a,101bから出力されるアナログ信号をデジタル信号に変換する。ADC102a,102bは、例えば、1チップの半導体装置で形成されている。図5に示した無線部12a,12bは、例えば、図12に示すRF部101a,101bおよびADC102a,102bによって、その機能が実現される。LSI103は、BBベースバンド処理部103aおよびCPUCentralProcessingUnit103bを有している。図1に示した減算部1aおよび計算部1bは、例えば、図12に示すBB処理部103aによって、その機能が実現される。また、図5に示したシンボル分離部13a,13b、チャネル推定部14、およびMLD処理部15は、例えば、図12に示すBB処理部103aによって、その機能が実現される。CPU103bは、通信制御や上位レイヤの通信処理を行う。"