"受信信号の搬送波周波数と局部発振器の出力信号の周波数との差を推定する周波数オフセット推定方法であって、所定のサンプリング周波数で予めサンプリングされた受信信号を離散フーリエ変換し、複数個の周波数成分を持つ周波数スペクトルを出力するステップと、前記周波数スペクトルの平均電力を算出するステップと、前記平均電力、または前記平均電力を定数倍した電力に、所定の値を加算してしきい値を算出するステップと、前記しきい値に基づいて、前記周波数スペクトルの各周波数成分の電力を1ビット量子化するステップと、前記各周波数成分の周波数と前記1ビット量子化された各周波数成分の電力とを掛け合わせ、この掛け合わせた積の和を求め、この積の和を前記1ビット量子化された前記周波数スペクトルの各周波数成分の電力の和で除算して重心周波数を算出するステップとを含む周波数オフセット推定方法。受信信号の搬送波周波数と局部発振器の出力信号の周波数との差を推定する周波数オフセット推定方法であって、所定のサンプリング周波数で予めサンプリングされた受信信号を離散フーリエ変換し、複数個の周波数成分を持つ周波数スペクトルを出力するステップと、前記周波数スペクトルの複数個の周波数成分をその電力の大きさで順に並べ替えるステップと、前記並べ替えられた周波数成分のうちの特定の順位の電力値または特定の順位の電力値を定数倍した電力にあらかじめ定められた値を加算してしきい値を算出するステップと、前記しきい値に基づいて、前記周波数スペクトルの各周波数成分の電力を1ビット量子化するステップと、前記各周波数成分の周波数と前記1ビット量子化された各周波数成分の電力とを掛け合わせ、この掛け合わせた積の和を求め、この積の和を前記1ビット量子化された前記周波数スペクトルの各周波数成分の電力の和で除算して重心周波数を算出するステップとを含む周波数オフセット推定方法。受信信号の搬送波周波数と局部発振器の出力信号の周波数との差を推定する周波数オフセット推定方法であって、所定のサンプリング周波数で予めサンプリングされた受信信号を離散フーリエ変換し、複数個の周波数成分を持つ周波数スペクトルを出力するステップと、前記周波数スペクトルの複数個の周波数成分をその電力の大きさで順に並べ替えるステップと、前記並べ替えられた複数個の周波数成分の電力の算術平均によりしきい値を算出するステップと、前記しきい値に基づいて、前記周波数スペクトルの各周波数成分の電力を1ビット量子化するステップと、前記各周波数成分の周波数と前記1ビット量子化された各周波数成分の電力とを掛け合わせ、この掛け合わせた積の和を求め、この積の和を前記1ビット量子化された前記周波数スペクトルの各周波数成分の電力の和で除算して重心周波数を算出するステップとを含む周波数オフセット推定方法。前記しきい値を算出するステップは、NS=2mであり、また、k=0,1,2,...,m-1であるとき、NS個の周波数成分を電力の大きさで降順に並べ替えたときの先頭から2k番目の電力値と2k1番目との電力値の平均値を前記しきい値とする請求項3に記載の周波数オフセット推定方法。受信信号の搬送波周波数と局部発振器の出力信号の周波数との差を推定する周波数オフセット推定方法であって、所定のサンプリング周波数で予めサンプリングされた受信信号を離散フーリエ変換し、複数個の周波数成分を持つ周波数スペクトルを出力するステップと、予め設定されたしきい値に基づいて、前記周波数スペクトルの各周波数成分の電力を1ビット量子化するステップと、前記各周波数成分の周波数と前記1ビット量子化された各周波数成分の電力とを掛け合わせ、この掛け合わせた積の和を求め、この積の和を前記1ビット量子化された前記周波数スペクトルの各周波数成分の電力の和で除算して重心周波数を算出するステップとを含む周波数オフセット推定方法。受信信号の搬送波周波数と局部発振器の出力信号の周波数との差を推定する周波数オフセット推定装置であって、所定のサンプリング周波数で予めサンプリングされた受信信号を離散フーリエ変換し、複数個の周波数成分を持つ周波数スペクトルを出力する離散フーリエ変換部と、前記離散フーリエ変換部により変換された周波数スペクトルの平均電力を算出する平均電力算出部と、前記平均電力算出部により算出された前記平均電力、または前記平均電力を定数倍した電力に、所定の値を加算してしきい値を算出するしきい値電力算出部と、前記しきい値電力算出部により算出された前記しきい値に基づいて、前記周波数スペクトルの各周波数成分の電力を1ビット量子化する1ビット量子化部と、前記離散フーリエ変換部により変換された周波数スペクトルの各周波数成分の周波数と前記1ビット量子化部により1ビット量子化された各周波数成分の電力とを掛け合わせ、この掛け合わせた積の和を求め、この積の和を前記1ビット量子化された前記周波数スペクトルの各周波数成分の電力の和で除算して重心周波数を算出する重心算出部とを備える周波数オフセット推定装置。受信信号の搬送波周波数と局部発振器の出力信号の周波数との差を推定する周波数オフセット推定装置であって、所定のサンプリング周波数で予めサンプリングされた受信信号を離散フーリエ変換し、複数個の周波数成分を持つ周波数スペクトルを出力する離散フーリエ変換部と、前記離散フーリエ変換部により変換された周波数スペクトルの複数個の周波数成分をその電力の大きさで順に並べ替える電力値整列部と、前記電力値整列部により並べ替えられた前記周波数成分のうちの特定の順位の電力値または特定の順位の電力値を定数倍した電力にあらかじめ定められた値を加算してしきい値を算出するしきい値電力算出部と、前記しきい値電力算出部により算出された前記しきい値に基づいて、前記周波数スペクトルの各周波数成分の電力を1ビット量子化する1ビット量子化部と、前記離散フーリエ変換部により変換された周波数スペクトルの各周波数成分の周波数と前記1ビット量子化部により1ビット量子化された各周波数成分の電力とを掛け合わせ、この掛け合わせた積の和を求め、この積の和を前記1ビット量子化された前記周波数スペクトルの各周波数成分の電力の和で除算して重心周波数を算出する重心算出部とを備える周波数オフセット推定装置。受信信号の搬送波周波数と局部発振器の出力信号の周波数との差を推定する周波数オフセット推定装置であって、所定のサンプリング周波数で予めサンプリングされた受信信号を離散フーリエ変換し、複数個の周波数成分を持つ周波数スペクトルを出力する離散フーリエ変換部と、前記離散フーリエ変換部により変換された周波数スペクトルの複数個の周波数成分をその電力の大きさで順に並べ替える電力値整列部と、前記電力値整列部により並べ替えられた複数個の周波数成分の電力の算術平均によりしきい値を算出するしきい値電力算出部と、前記しきい値電力算出部により算出された前記しきい値に基づいて、前記周波数スペクトルの各周波数成分の電力を1ビット量子化する1ビット量子化部と、前記離散フーリエ変換部により変換された周波数スペクトルの各周波数成分の周波数と前記1ビット量子化部により1ビット量子化された各周波数成分の電力とを掛け合わせ、この掛け合わせた積の和を求め、この積の和を前記1ビット量子化された前記周波数スペクトルの各周波数成分の電力の和で除算して重心周波数を算出する重心算出部とを備える周波数オフセット推定装置。前記しきい値電力算出部は、NS=2mであり、また、k=0,1,2,...,m-1であるとき、NS個の周波数成分を電力の大きさで降順に並べ替えたときの先頭から2k番目の電力値と2k1番目との電力値の平均値を前記しきい値とする請求項8に記載の周波数オフセット推定装置。受信信号の搬送波周波数と局部発振器の出力信号の周波数との差を推定する周波数オフセット推定装置であって、所定のサンプリング周波数で予めサンプリングされた受信信号を離散フーリエ変換し、複数個の周波数成分を持つ周波数スペクトルを出力する離散フーリエ変換部と、予め設定されたしきい値に基づいて、前記周波数スペクトルの各周波数成分の電力を1ビット量子化する1ビット量子化部と、前記離散フーリエ変換部により変換された周波数スペクトルの各周波数成分の周波数と前記1ビット量子化部により1ビット量子化された各周波数成分の電力とを掛け合わせ、この掛け合わせた積の和を求め、この積の和を前記1ビット量子化された前記周波数スペクトルの各周波数成分の電力の和で除算して重心周波数を算出する重心算出部とを備える周波数オフセット推定装置。しかしながら、数式1により周波数オフセットを算出する従来の周波数オフセット推定方法を、ハードウェアのデジタルシグナルプロセッサDSPに実装する場合、除算演算の回路規模が大きくなるという問題がある。一般に、周波数オフセット推定方法をDSPに実装する際に、除算器に浮動小数点演算を用いた場合には、固定小数点演算に比べて著しく回路規模が増大し、計算速度も遅くなるため、固定小数点演算を行う必要がある。固定小数点演算は、割る数が2のべき乗であれば、ビットシフトのみの簡単な演算で処理できるが、割る数が2のべき乗以外の場合には、複雑な演算となる。特に、割る数のダイナミックレンジが大きい場合、オーバーフローや、アンダーフローを生じさせないことと、演算精度を高く取ることとを両立させるような小数点の位置を決めることが困難であり、余裕を持った桁取りを行うと、回路規模が大きくなる。数式1の分母は、入力信号の周波数成分の電力の和であり、通信システムにより定まる入力信号電力のダイナミックレンジに対応して、様々な実数の値を取るため、回路規模が大きくなるという問題がある。本発明は、このような事情を考慮してなされたものであり、その目的は、周波数スペクトルの形状を利用して周波数オフセットを算出する際に、除算演算の割る数を2のべき乗、または限定された範囲の整数値にすることができるため、ハードウェア化した場合の回路規模を小さくすることができる周波数オフセット推定方法および周波数オフセット推定装置を提供することにある。上述した課題を解決するために、本発明は、受信信号の搬送波周波数と局部発振器の出力信号の周波数との差を推定する周波数オフセット推定方法であって、所定のサンプリング周波数で予めサンプリングされた受信信号を離散フーリエ変換し、複数個の周波数成分を持つ周波数スペクトルを出力するステップと、前記周波数スペクトルの平均電力を算出するステップと、前記平均電力、または前記平均電力を定数倍した電力に、所定の値を加算してしきい値を算出するステップと、前記しきい値に基づいて、前記周波数スペクトルの各周波数成分の電力を1ビット量子化するステップと、前記各周波数成分の周波数と前記1ビット量子化された各周波数成分の電力とを掛け合わせ、この掛け合わせた積の和を求め、この積の和を前記1ビット量子化された前記周波数スペクトルの各周波数成分の電力の和で除算して重心周波数を算出するステップとを含む周波数オフセット推定方法である。本発明は、受信信号の搬送波周波数と局部発振器の出力信号の周波数との差を推定する周波数オフセット推定方法であって、所定のサンプリング周波数で予めサンプリングされた受信信号を離散フーリエ変換し、複数個の周波数成分を持つ周波数スペクトルを出力するステップと、前記周波数スペクトルの複数個の周波数成分をその電力の大きさで順に並べ替えるステップと、前記並べ替えられた周波数成分のうちの特定の順位の電力値または特定の順位の電力値を定数倍した電力にあらかじめ定められた値を加算してしきい値を算出するステップと、前記しきい値に基づいて、前記周波数スペクトルの各周波数成分の電力を1ビット量子化するステップと、前記各周波数成分の周波数と前記1ビット量子化された各周波数成分の電力とを掛け合わせ、この掛け合わせた積の和を求め、この積の和を前記1ビット量子化された前記周波数スペクトルの各周波数成分の電力の和で除算して重心周波数を算出するステップとを含む周波数オフセット推定方法である。本発明は、受信信号の搬送波周波数と局部発振器の出力信号の周波数との差を推定する周波数オフセット推定方法であって、所定のサンプリング周波数で予めサンプリングされた受信信号を離散フーリエ変換し、複数個の周波数成分を持つ周波数スペクトルを出力するステップと、前記周波数スペクトルの複数個の周波数成分をその電力の大きさで順に並べ替えるステップと、前記並べ替えられた複数個の周波数成分の電力の算術平均によりしきい値を算出するステップと、前記しきい値に基づいて、前記周波数スペクトルの各周波数成分の電力を1ビット量子化するステップと、前記各周波数成分の周波数と前記1ビット量子化された各周波数成分の電力とを掛け合わせ、この掛け合わせた積の和を求め、この積の和を前記1ビット量子化された前記周波数スペクトルの各周波数成分の電力の和で除算して重心周波数を算出するステップとを含む周波数オフセット推定方法である。本発明において、前記しきい値を算出するステップは、NS=2mであり、また、k=0,1,2,...,m-1であるとき、NS個の周波数成分を電力の大きさで降順に並べ替えたときの先頭から2k番目の電力値と2k1番目との電力値の平均値を前記しきい値とするようにしても良い。本発明は、受信信号の搬送波周波数と局部発振器の出力信号の周波数との差を推定する周波数オフセット推定方法であって、所定のサンプリング周波数で予めサンプリングされた受信信号を離散フーリエ変換し、複数個の周波数成分を持つ周波数スペクトルを出力するステップと、予め設定されたしきい値に基づいて、前記周波数スペクトルの各周波数成分の電力を1ビット量子化するステップと、前記各周波数成分の周波数と前記1ビット量子化された各周波数成分の電力とを掛け合わせ、この掛け合わせた積の和を求め、この積の和を前記1ビット量子化された前記周波数スペクトルの各周波数成分の電力の和で除算して重心周波数を算出するステップとを含む周波数オフセット推定方法である。また、上述した課題を解決するために、本発明は、受信信号の搬送波周波数と局部発振器の出力信号の周波数との差を推定する周波数オフセット推定装置であって、所定のサンプリング周波数で予めサンプリングされた受信信号を離散フーリエ変換し、複数個の周波数成分を持つ周波数スペクトルを出力する離散フーリエ変換部と、前記離散フーリエ変換部により変換された周波数スペクトルの平均電力を算出する平均電力算出部と、前記平均電力算出部により算出された前記平均電力、または前記平均電力を定数倍した電力に、所定の値を加算してしきい値を算出するしきい値電力算出部と、前記しきい値電力算出部により算出された前記しきい値に基づいて、前記周波数スペクトルの各周波数成分の電力を1ビット量子化する1ビット量子化部と、前記離散フーリエ変換部により変換された周波数スペクトルの各周波数成分の周波数と前記1ビット量子化部により1ビット量子化された各周波数成分の電力とを掛け合わせ、この掛け合わせた積の和を求め、この積の和を前記1ビット量子化された前記周波数スペクトルの各周波数成分の電力の和で除算して重心周波数を算出する重心算出部とを備える周波数オフセット推定装置である。本発明は、受信信号の搬送波周波数と局部発振器の出力信号の周波数との差を推定する周波数オフセット推定装置であって、所定のサンプリング周波数で予めサンプリングされた受信信号を離散フーリエ変換し、複数個の周波数成分を持つ周波数スペクトルを出力する離散フーリエ変換部と、前記離散フーリエ変換部により変換された周波数スペクトルの複数個の周波数成分をその電力の大きさで順に並べ替える電力値整列部と、前記電力値整列部により並べ替えられた前記周波数成分のうちの特定の順位の電力値または特定の順位の電力値を定数倍した電力にあらかじめ定められた値を加算してしきい値を算出するしきい値電力算出部と、前記しきい値電力算出部により算出された前記しきい値に基づいて、前記周波数スペクトルの各周波数成分の電力を1ビット量子化する1ビット量子化部と、前記離散フーリエ変換部により変換された周波数スペクトルの各周波数成分の周波数と前記1ビット量子化部により1ビット量子化された各周波数成分の電力とを掛け合わせ、この掛け合わせた積の和を求め、この積の和を前記1ビット量子化された前記周波数スペクトルの各周波数成分の電力の和で除算して重心周波数を算出する重心算出部とを備える周波数オフセット推定装置である。本発明は、受信信号の搬送波周波数と局部発振器の出力信号の周波数との差を推定する周波数オフセット推定装置であって、所定のサンプリング周波数で予めサンプリングされた受信信号を離散フーリエ変換し、複数個の周波数成分を持つ周波数スペクトルを出力する離散フーリエ変換部と、前記離散フーリエ変換部により変換された周波数スペクトルの複数個の周波数成分をその電力の大きさで順に並べ替える電力値整列部と、前記電力値整列部により並べ替えられた複数個の周波数成分の電力の算術平均によりしきい値を算出するしきい値電力算出部と、前記しきい値電力算出部により算出された前記しきい値に基づいて、前記周波数スペクトルの各周波数成分の電力を1ビット量子化する1ビット量子化部と、前記離散フーリエ変換部により変換された周波数スペクトルの各周波数成分の周波数と前記1ビット量子化部により1ビット量子化された各周波数成分の電力とを掛け合わせ、この掛け合わせた積の和を求め、この積の和を前記1ビット量子化された前記周波数スペクトルの各周波数成分の電力の和で除算して重心周波数を算出する重心算出部とを備える周波数オフセット推定装置である。本発明において、前記しきい値電力算出部は、NS=2mであり、また、k=0,1,2,...,m-1であるとき、NS個の周波数成分を電力の大きさで降順に並べ替えたときの先頭から2k番目の電力値と2k1番目との電力値の平均値を前記しきい値とするようにしても良い。本発明は、受信信号の搬送波周波数と局部発振器の出力信号の周波数との差を推定する周波数オフセット推定装置であって、所定のサンプリング周波数で予めサンプリングされた受信信号を離散フーリエ変換し、複数個の周波数成分を持つ周波数スペクトルを出力する離散フーリエ変換部と、予め設定されたしきい値に基づいて、前記周波数スペクトルの各周波数成分の電力を1ビット量子化する1ビット量子化部と、前記離散フーリエ変換部により変換された周波数スペクトルの各周波数成分の周波数と前記1ビット量子化部により1ビット量子化された各周波数成分の電力とを掛け合わせ、この掛け合わせた積の和を求め、この積の和を前記1ビット量子化された前記周波数スペクトルの各周波数成分の電力の和で除算して重心周波数を算出する重心算出部とを備える周波数オフセット推定装置である。この発明によれば、周波数スペクトルの形状を利用して周波数オフセットを算出する回路をハードウェア化した場合の回路規模を小さくすることができる。以下、本発明の一実施形態を、図面を参照して説明する。なお、以下の各図において、同一のまたは対応する構成には同一の符号を用いている。本発明の各実施形態では、周波数スペクトルの形状を利用して周波数オフセットを算出する際に、予め閾値電力を決めて、閾値電力と各周波数成分の信号電力を比較することで、各周波数成分の信号電力を1ビット量子化し、量子化された信号電力から重心周波数を求めることを特徴とする。これにより、重心周波数は、周波数オフセットに等しくなるので、周波数オフセットの推定が可能になるという効果を奏する。また、重心を演算除算する際の分母が1からNSまでの整数値になるので回路規模を小さくすることが可能になるという効果を奏する。また、閾値電力を求める際に、中央値medianを用いることで、除算の分母を2のべき乗とすることができるので、大幅に回路規模を削減することが可能になる。以下、これらの特徴についてさらに詳細に説明する。A.第1実施形態図1は、本発明の第1実施形態による周波数オフセット推定装置の構成例を示すブロック図である。図1において、周波数オフセット推定装置4は、直並列変換回路5、離散フーリエ変換回路DFT6、平均電力算出回路9、しきい値電力算出回路10、1ビット量子化器11、重心算出回路7により構成される。また、周波数オフセット推定装置4の周辺回路として、第1のA/D変換器1、第2のA/D変換器2、合成回路3が周波数オフセット推定装置4に接続されている。周波数オフセット推定装置4に入力された信号は、1つの入力ポートとNS個の出力ポートとを持つ直並列変換回路5により直並列変換された後に、離散フーリエ変換回路6により所定のサンプリング周波数で予めサンプリングされた受信信号が離散フーリエ変換され、NS個の周波数成分を持つ周波数スペクトルが出力される。離散フーリエ変換のサイズNSは、一般に、2のべき乗の整数である。離散フーリエ変換DFTを実行するためのアルゴリズムとしては、高速フーリエ変換FFTなどを用いることができる。離散フーリエ変換回路6の出力信号は、2分岐されて、一方は1ビット量子化器11に入力され、他方は平均電力算出回路9に入力される。この平均電力算出回路9では、離散フーリエ変換回路6により出力される周波数スペクトルの平均電力を算出する。しきい値電力算出回路10では、この平均電力を元にして、しきい値電力を算出する。1ビット量子化器11は、このしきい値電力を用いて離散フーリエ変換回路6の出力信号を1ビット量子化する。1ビット量子化器11はコンパレータとして動作する。重心算出回路7は、各周波数成分の周波数と1ビット量子化された各周波数成分の電力とに基づいて重心周波数本第1実施形態では周波数オフセットに等しいを算出する。図2は、第1実施形態の動作を示し、周波数オフセットがないときの周波数スペクトルを示す概念図である。図3は、本第1実施形態の動作を示し、周波数オフセットがあるときの周波数スペクトルを示す概念図である。図2、図3は、共に離散フーリエ変換回路6の出力信号と1ビット量子化器11の出力信号であり、受信信号の周波数スペクトルと1ビット量子化した周波数スペクトルを表している。離散フーリエ変換回路6の出力信号である周波数スペクトルのうち、周波数fiの周波数成分の信号電力をPi、平均電力算出回路9で算出される平均電力をPAVE、しきい値電力算出回路10で算出されるしきい値電力をPTH、1ビット量子化器11の出力信号のうち、周波数fiの周波数成分の信号電力をPQiとすると、PAVE、PTH、及びPQiは、次式2、3、4となる。この場合の重心周波数fCGは、次式5で与えられる。但し、FSは、サンプリング周波数である。iは、1からNSまでの整数値を取る変数であり、数式6により周波数fiは、-1/21/NSFSから1/2FSまでの値を取る。周波数スペクトルが中央の周波数に関して対称であれば、数式5は、ほぼ0となる。今、周波数fiの周波数成分の信号電力は、Piであり、周波数オフセットΔfが加えられた場合の周波数スペクトルにおいては、fiΔfの周波数成分の信号電力がPiとなる。また、この関係は、1ビット量子化後も保持される。従って、周波数オフセットΔfがある場合の重心周波数fCGは次式7となる。すなわち、重心周波数fCGは、周波数オフセットΔfに等しくなる。本第1実施形態では、平均電力を算出する数式2の分母NSは、2のべき乗であるため、回路規模を小さくすることができる。しきい値電力を算出する数式3は、平均電力または平均電力を定数倍した電力に予め定められた値を加算する演算であり、除算が不要であるため、回路規模は、小さくできる。また、重心周波数fCGを算出する数式5の分母は、周波数スペクトルの形状や、入力信号電力の大小にかかわらず、1からNSまでの範囲の整数値に限定されるため、回路規模を小さくすることができる。図4ないし図6は、本第1実施形態の動作をシミュレーションした結果を示す概念図である。図4は、シンボルレート31.5GHzのQPSK四位相偏移変調変調信号に5GHzの周波数オフセットを与えた時の周波数スペクトルである。図5は、平均電力をしきい値として1ビット量子化した時の周波数スペクトルであり、周波数オフセットに対応して周波数スペクトルが中心から右にシフトしていることが分かる。図6では、周波数オフセットの設定値を横軸にとり、本第1実施形態による周波数オフセット推定結果を縦軸にとって、シミュレーション結果をプロットしている。1ビット量子化による演算量削減を行っても、周波数オフセットを精度良く推定できていることが分かる。なお、上述した第1実施形態において、第1のA/D変換器1、及び第2のA/D変換器2に入力される受信信号の変調方式は、OFDMであっても、シングルキャリア変調であってもよい。B.第2実施形態次に、本発明の第2実施形態について説明する。図7は、本発明の第2実施形態による周波数オフセット推定装置の構成例を示すブロック図である。図7において、参照符号12は電力値整列回路である。電力値整列回路12では、出力される周波数スペクトルの周波数成分毎の電力を大きさの昇順、または降順に並べ替える。しきい値電力算出回路10では、並べ替えられた電力値のうち、予め定めた順番の電力値を用いてしきい値を算出する。例えば、周波数スペクトルの周波数成分の個数をNS、1以上NS以下の任意の整数をkとし、周波数成分を電力の大きさで昇順または降順に並べたときに、先頭からk番目の値またはk番目の値を定数倍した電力にあらかじめ定められた値を加算して得られる値をしきい値電力に設定することができる。すなわち、しきい値電力は、受信信号電力と相関を持った値に設定すれば良く、第1実施形態では平均電力、本第2実施形態では周波数スペクトルの周波数成分の値を元に算出している。特に、NS個の周波数成分を電力の大きさで降順に並べ替えたときの先頭からNS/2番目の電力値とNS/21番目との電力値の平均値、すなわち中央値medianを1ビット量子化器11が用いるしきい値電力PTHとした場合、数式4の右辺の1と0は等しい個数となり、NS/2個ずつとなる。従って、数式5の分母は、NS/2となり、この値は2のべき乗であるため、回路規模を大幅に削減することができる。同様に、NS個の周波数成分を降順に並べ替えたときの先頭からNS/4番目の電力値とNS/41番目との電力値の平均値を1ビット量子化器11が用いるしきい値電力PTHとした場合、数式4の右辺の1はNS/4個、0は3NS/4個となる。従って、数式5の分母は、NS/4となり、この値も2のべき乗であるため、回路規模を大幅に削減することができる。すなわち、NS=2mであり、また、k=0,1,2,...,m-1であるとき、周波数成分を電力の大きさの降順に並べ替えたときの先頭から2k番目の電力値と2k1番目との電力値の平均値をしきい値とすれば、数式5の分母は、2のべき乗となり、回路規模を大幅に削減することができる。kの値の選択は、変調方式や、伝送路の帯域制限素子によって定まる周波数スペクトルの形状に合わせて選択すればよい。図8、図9は、本第2実施形態の動作をシミュレーションした結果を示す概念図である。図8には、シンボルレート31.5GHzのQPSK変調信号に5GHzの周波数オフセットを与え、中央値をしきい値として1ビット量子化したときの周波数スペクトルを示し、周波数オフセットに対応して周波数スペクトルが中心から右にシフトしていることが分かる。図9では、周波数オフセットの設定値を横軸にとり、本第2実施形態による周波数オフセット推定結果を縦軸にとって、シミュレーション結果をプロットしており、周波数オフセットを精度良く推定できていることが分かる。また、第1、及び第2の実施形態の説明において、周波数オフセット推定には、離散フーリエ変換後の周波数スペクトルの周波数成分をすべて用いて重心を算出することとしていた。しかし、要求される推定精度と回路規模とのトレードオフにより、周波数成分を数個おきに間引いて重心を算出することも可能である。また、周波数オフセット推定装置4,4aに入力される信号レベルが常に一定であるように、自動利得制御回路などの他の回路で制御されていれば、しきい値電力を予め設定しておき、運用中は変更せずに動作させることが可能である。この場合、数式3の定数aを0とし、瞬時的な受信信号電力とは無関係な定数bのみでしきい値電力を決めれば良い。定数bは、受信信号電力が変動する場合と同様に、周波数オフセット推定装置4,4aに入力される電力またはこの電力を定数倍した電力にあらかじめ定められた値を加算して得られる値や、周波数スペクトルの周波数成分を電力の大きさで昇順または降順に並べたときに、先頭からk番目の値またはk番目の値を定数倍した電力にあらかじめ定められた値を加算して得られる値から決めることができる。なお、この場合、図1に示した平均電力算出回路9およびしきい値電力算出回路10、ならびに、図7に示した電力値整列回路12およびしきい値電力算出回路10は省略することができる。すなわち、1ビット量子化器11は、予め設定されたしきい値電力を用いて離散フーリエ変換回路6の出力信号を1ビット量子化する。以上説明したように、本発明の各実施形態によれば、周波数スペクトルの形状を利用する周波数オフセット推定方法において除算演算の割る数を2のべき乗、または限定された範囲の整数値にすることができるため、ハードウェア化した場合の回路規模を小さくできる。以上、本発明の実施形態について図面を参照して詳述したが、具体的な構成はこれら実施形態に限られるものではなく、本発明の要旨を逸脱しない範囲の設計等構成の付加、省略、置換、およびその他の変更も含まれる。本発明は前述した説明によって限定されることはなく、添付の請求の範囲によってのみ限定される。"