(57)【要約】【課題】周波数誤差の推定の精度と範囲を向上させる。【解決手段】周波数誤差検出装置１５０は、既知シンボルに対する受信データの誤差であるシンボル誤差の複数のラグの自己相関の総和である自己相関和を、複数のラグのうちの最大ラグを小から大に変化させて求め、その自己相関和から、搬送波に対する周波数誤差の推定値を求める。累積加算装置１５４は、周波数誤差の推定値の累積加算を行い、その累積加算により得られる累積加算値を、周波数誤差の最終的な推定値として出力する。本発明は、例えば、搬送波をディジタル変調した変調信号を受信する受信装置等に適用できる。【選択図】図９
2008167230.tif 000002
【特許請求の範囲】
【請求項１】
搬送波をディジタル変調した変調信号をベースバンドの信号に変換することにより得られる、前記搬送波と同相のI成分と、前記搬送波と直交するQ成分とからなる受信データを処理するデータ処理装置において、
既知のシンボルである既知シンボルに対する前記受信データの誤差であるシンボル誤差の複数のラグの自己相関の総和である自己相関和を、前記複数のラグのうちの最大のラグである最大ラグを小から大に変化させて求め、その自己相関和から、前記搬送波に対する周波数誤差の推定値を求める周波数誤差検出手段と、
前記周波数誤差の推定値の累積加算を行い、その累積加算により得られる累積加算値を、前記周波数誤差の最終的な推定値として出力する累積加算手段と
を備えるデータ処理装置。
【請求項２】
前記周波数誤差の最終的な推定値に基づいて、周波数誤差を補正する周波数補正手段をさらに備える
請求項１に記載のデータ処理装置。
【請求項３】
前記周波数誤差検出手段は、
既知シンボルの前記受信データと、前記既知シンボルの複素共役とを乗算することにより、前記既知シンボルに対する前記受信データのシンボル誤差を求める誤差演算手段と、
既知シンボルの系列である既知シンボル系列について、前記シンボル誤差の複数のラグの自己相関の総和である自己相関和を求める相関演算手段と、
複数の既知シンボル系列それぞれについて、前記複数のラグの自己相関和から、前記周波数誤差の第１の推定値を求める周波数誤差演算手段と、
複数の既知シンボル系列それぞれについて求められた周波数誤差の第１の推定値を平均化することにより、前記周波数誤差の第２の推定値を求める推定値平均化手段と
を有し、
前記累積加算手段は、前記周波数誤差の第２の推定値の累積加算を行い、その累積加算により得られる累積加算値を、前記周波数誤差の最終的な推定値として出力する
請求項１に記載のデータ処理装置。
【請求項４】
前記相関演算手段は、
前記シンボル誤差の所定のラグの自己相関を求める複数の相関部と、
前記複数の相関部で求められた自己相関を加算することにより、既知シンボル系列についての、前記複数のラグの自己相関和を求める演算部と
を有し、
前記複数の相関部のうちの、前記演算部での加算の対象とする自己相関を求める相関部を指定する制御手段をさらに備える
請求項３に記載のデータ処理装置。
【請求項５】
前記周波数誤差検出手段は、
既知シンボルの前記受信データと、前記既知シンボルの複素共役とを乗算することにより、前記既知シンボルに対する前記受信データのシンボル誤差を求める誤差演算手段と、
既知シンボルの系列である既知シンボル系列について、前記シンボル誤差の複数のラグの自己相関の総和である自己相関和を求める相関演算手段と、
複数の既知シンボル系列それぞれについて求められた、前記複数のラグの自己相関和を平均化し、平均自己相関和を求める自己相関和平均化手段と、
前記平均自己相関和から、前記周波数誤差の推定値を求める周波数誤差演算手段と
を有する
請求項１に記載のデータ処理装置。
【請求項６】
前記相関演算手段は、
前記シンボル誤差の所定のラグの自己相関を求める複数の相関部と、
前記複数の相関部で求められた自己相関を加算することにより、既知シンボル系列についての、前記複数のラグの自己相関和を求める演算部と
を有し、
前記複数の相関部のうちの、前記演算部での加算の対象とする自己相関を求める相関部を指定する制御手段をさらに備える
請求項５に記載のデータ処理装置。
【請求項７】
搬送波をディジタル変調した変調信号をベースバンドの信号に変換することにより得られる、前記搬送波と同相のI成分と、前記搬送波と直交するQ成分とからなる受信データを処理するデータ処理方法において、
既知のシンボルである既知シンボルに対する前記受信データの誤差であるシンボル誤差の複数のラグの自己相関の総和である自己相関和を、前記複数のラグのうちの最大のラグである最大ラグを小から大に変化させて求め、その自己相関和から、前記搬送波に対する周波数誤差の推定値を求め、
前記周波数誤差の推定値の累積加算を行い、その累積加算により得られる累積加算値を、前記周波数誤差の最終的な推定値として出力する
ステップを含むデータ処理方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、データ処理装置及びデータ処理方法に関し、特に、例えば、周波数誤差の推定の精度と範囲を向上させることができるようにするデータ処理装置及びデータ処理方法に関する。
【背景技術】
【０００２】
例えば、BS(Broadcasting Satellite)ディジタル放送を行う送信装置では、搬送波が、送信対象のデータで、ディジタル変調され、すなわち、例えば、BPSK(Binary PSK(Phase Shift Keying))，QPSK(Quadrature PSK)、又は8PSK等でディジタル変調され、その結果得られる変調信号が送信される。
【０００３】
また、BSディジタル放送では、フレームと呼ばれる単位で、データが送信される。
【０００４】
フレームは、その先頭から、ユニークなデータW1のシンボル、TMCC(Transmission Multiplexing Configuration Control)のシンボル、ユニークなデータW2又はW3のシンボル、主信号のシンボルとバースト信号のシンボルとのセットが順次配置されて構成される。
【０００５】
データW1,W2、及びW3のシンボルは、ユニークで既知のシンボル（既知シンボル）であり、BSディジタル放送を受信する受信装置において、同期（フレーム同期）を確立するのに用いられる。
【０００６】
TMCCは、変調信号における主信号の変調方式や符号化方式等を含む制御情報である。
【０００７】
主信号は、画像データをMPEG(Moving Picture Experts Group)符号化した符号化データ等の、送信すべき本来の情報である。
【０００８】
バースト信号のシンボルは、フレーム周期のPRBS(Pseudo Random Bit/Binary Sequence)系列であり、既知シンボルである。バースト信号のシンボルは、受信装置において、低C/N(Carrier to Noise Ratio)であっても、同期（搬送波同期）を確立することができるように、フレームに、間欠的に配置される。
【０００９】
すなわち、受信装置では、変調信号に、搬送波に相当する信号（以下、適宜、搬送波相当信号という）を乗算することにより、変調信号を、搬送波と同相のI成分と、搬送波と直交するQ成分とからなるベースバンドの復調信号に復調するが、受信装置で用いられる搬送波相当信号と、変調信号を送信してくる送信装置で用いられる搬送波との間には、一般に、誤差があり、その誤差に起因して、受信装置で得られる復調信号のシンボルは、I成分を表すI軸と、Q成分を表すQ軸とで規定されるIQ平面において回転する。
【００１０】
受信装置では、以上のような復調信号のシンボルの回転を補償するために、搬送波同期が確立される。
【００１１】
搬送波同期を確立する場合において、シンボルが、IQ平面上の信号点どうしの距離が比較的近いQPSK変調、又は8PSK変調をされているときには、C/Nが低いと、搬送波同期を確立することが困難になることがある。
【００１２】
そこで、変調信号には、既知シンボルで、しかも、IQ平面上での信号点どうしの距離が遠いBPSK変調がされたバースト信号のシンボルが、間欠的に配置され、このバースト信号のシンボルにより、C/Nが低くても、受信装置において、搬送波同期を確立することができるようになっている。
【００１３】
ここで、ディジタル変調を含む変調では、搬送波を、送信の対象のデータにしたがって変調することにより、変調信号が得られるのであるが、本明細書では、説明の便宜上、搬送波を、送信の対象のデータにしたがって変調することを、送信の対象のデータを変調するともいう。
【００１４】
上述したように、受信装置において、変調信号の復調に用いられる搬送波相当信号と、変調信号を送信してくる送信装置で用いられる搬送波との間には、一般に、誤差があり、この誤差、すなわち、受信装置の搬送波相当信号の、送信装置の搬送波に対する誤差には、周波数誤差と位相誤差とがある。
【００１５】
受信装置では、既知シンボルを用いて周波数誤差が推定（検出）され、その推定値に基づいて、周波数誤差が補正される。
【００１６】
周波数誤差を推定する推定アルゴリズムとしては、L&R(M. Luise and R. Reggiannini)アルゴリズム（L&R検出器）がある（例えば、非特許文献１を参照）。
【００１７】
L&Rアルゴリズムでは、送信装置から送信されてくる変調信号をベースバンドの信号に変換することにより得られる、I成分とQ成分とからなるアナログの復調信号r(t)が、式（１）で表される。
【００１８】
【数１】
2008167230.tif 000003
・・・（１）
【００１９】
ここで、式（１）において、tは、時刻を表し、r(t)は、復調信号を表す。また、eは、ネイピア数(Napier's constant)を表し、jは、虚数単位(imaginary unit)を表す。さらに、πは、円周率を、vは、周波数誤差を、θは、位相誤差を、それぞれ表す。また、kは、サンプリングの時刻（サンプル点）を表し、c(k)は、kサンプル目の、送信対象のデータ（シンボル）を表す。さらに、g(t)は、送信装置でかけられる波形整形フィルタ（の係数）を、τは、時間のジッタを、w(t)は、ノイズを表す。
【００２０】
受信装置では、例えば、復調信号r(t)が、受信フィルタ（マッチフィルタ）g(-t)でフィルタリングされることにより波形整形され、さらに、シンボルの周期に等しい時間Tごとにサンプリング（A/D(Analog/Digital)変換）され、ディジタル信号の受信データy(k)とされる。この受信データy(k)は、式（２）で表される。
【００２１】
【数２】
2008167230.tif 000004
・・・（２）
【００２２】
ここで、式（２）において、n(k)は、ノイズを表す。
【００２３】
いま、式（２）で表される受信データy(k)を、既知シンボルの系列である既知シンボル系列の受信データであるとし、１の既知シンボル系列が、L0個のシンボルから構成されることとする。なお、L0を、既知シンボル系列長という。
【００２４】
また、１の既知シンボル系列の、先頭からk+1番目(k=0,1,・・・,L0)の既知シンボルは、c(k)で表され、その既知シンボルc(k)に対応する受信データが、y(k)で表されることとする。
【００２５】
さらに、既知シンボルc(k)の複素共役を、c*(k)と表すこととすると、c(k)c*(k)は1であるから、受信データy(k)と、既知シンボルc(k)の複素共役c*(k)とを乗算することにより、既知シンボルc(k)に対する受信データy(k)の誤差（以下、適宜、シンボル誤差という）z(k)が、式（３）のように求められる。
【００２６】
【数３】
2008167230.tif 000005
・・・（３）
【００２７】
ここで、式（３）において、n'(k)は、n(k)c*(k)である。
【００２８】
また、ラグがmの、シンボル誤差z(k)の自己相関を、R(m)と表すこととすると、自己相関R(m)は、式（４）で表される。
【００２９】
【数４】
2008167230.tif 000006
・・・（４）
【００３０】
ここで、式（４）において、z*(m)は、シンボル誤差z(m)の複素共役であり、n''(m)は、式（３）のn'(k)が関係するノイズ成分を表す。また、ラグmは、1,2,・・・,Nの値をとり、Nは、ラグの最大値（最大ラグ）で、既知シンボル系列長L0より小さい整数である。
【００３１】
１の既知シンボル系列について、複数のラグの自己相関R(m)、すなわち、ラグmが1,2,・・・,Nの自己相関R(1),R(2),・・・,R(N)の総和（以下、適宜、自己相関和ともいう）ΣR(m)は、式（５）で近似することができる。
【００３２】
【数５】
2008167230.tif 000007
・・・（５）
【００３３】
ここで、式（５）では、ノイズ成分n''(m)の平均を、0に近似している。
【００３４】
式（５）の、mを1からNに変化させてのej2πmvTのサメーションであるΣej2πmvTは、式x0+x1+・・・+xN=(xN+1-1)/(x-1)を利用して、式（６）で表すことができる。
【００３５】
【数６】
2008167230.tif 000008
・・・（６）
【００３６】
式（５）に、式（６）を代入することにより、式（５）の自己相関和ΣR(m)は、式（７）で表される。
【００３７】
【数７】
2008167230.tif 000009
・・・（７）
【００３８】
したがって、実軸を横軸とするとともに、虚軸を縦軸とする複素平面(complex plane)において、自己相関和ΣR(m)を表す点と原点とを結ぶ線分が、実軸となす角度（偏角 (argument) ）を、arg[ΣR(m)]と表すこととすると、周波数誤差vは、式（７）の両辺のarg[]をとることにより、式（８）にしたがって求めることができる。
【００３９】
【数８】
2008167230.tif 000010
・・・（８）
【００４０】
なお、図１に示すような不確定性をなくすため、式（７）右辺のsin(πNvT)/sin(πvT)は、正である必要があり、そのためには、式（９）が成立する必要がある。
【００４１】
【数９】
2008167230.tif 000011
・・・（９）
【００４２】
式（９）から、式（８）によって求めることができる周波数誤差vの範囲は、式（１０）で表される範囲となる。
【００４３】
【数１０】
2008167230.tif 000012
・・・（１０）
【００４４】
図２は、以上のようなL&Rアルゴリズムにより周波数誤差を推定する、BSディジタル放送の、従来の受信装置の一例の構成を示している。
【００４５】
受信装置では、変調信号が復調されることにより、式（１）で表されるベースバンドの復調信号r(t)に変換され、その復調信号r(t)が、マッチフィルタ１に供給される。
【００４６】
マッチフィルタ１では、復調信号r(t)がフィルタリングされることにより波形整形され、A/D変換部２に出力される。A/D変換部２は、マッチフィルタ１の出力をサンプリングし、式（２）で表される受信データy(k)を出力する。
【００４７】
A/D変換部２が出力する受信データy(k)は、周波数誤差検出装置３に供給される。
【００４８】
周波数誤差検出装置３は、誤差演算部１１、相関演算部１２、及び周波数誤差演算部１３から構成され、L&Rアルゴリズムにより周波数誤差を推定する。
【００４９】
すなわち、誤差演算部１１は、演算部２１から構成され、演算部２１は、A/D変換部２が出力する受信データy(k)のうちの、既知シンボルc(k)の受信データy(k)と、既知シンボルc(k)の複素共役c*(k)とを乗算することにより、式（３）で表されるシンボル誤差z(k)を求め、相関演算部１２に供給する。
【００５０】
相関演算部１２は、N個の相関部３１1，３１2，・・・，３１Nと、演算部３２とから構成される。
【００５１】
相関部３１mは、式（４）にしたがい、既知シンボル系列について、演算部２１から供給されるシンボル誤差z(k)の、ラグがmの自己相関R(m)を求め、演算部３２に供給する。
【００５２】
演算部３２は、N個の相関部３１1ないし３１Nそれぞれから供給される自己相関R(1)ないしR(N)の総和を求め、その結果得られる、式（５）の自己相関和ΣR(m)を、周波数誤差演算部１３に供給する。
【００５３】
周波数誤差演算部１３は、演算部３２からの自己相関和ΣR(m)を用い、式（８）にしたがって、周波数誤差の推定値vを求めて出力する。
【００５４】
受信装置では、以上のようにして得られる周波数誤差の推定値vに基づき、周波数誤差を補正するAFC(Automatic Frequency Control)が行われる。
【００５５】
なお、遅延検波をして、大まかな周波数偏差を推定し、この周波数偏差に基づいて、受信信号の周波数偏差を補償した後、その補償後の受信信号に対し、再度、遅延検波を行うことにより、精度の高い周波数偏差の推定を行うAFC装置がある（例えば、特許文献１を参照）。
【００５６】
【特許文献１】特許第3230567号
【非特許文献１】Umberto Mengali, Aldo N. D'Andrea, "Synchronization Techniques for Digital Receivers (Applications of Communications Theory) ", Plenum Pub Corp (1997/11)
【発明の開示】
【発明が解決しようとする課題】
【００５７】
上述したL&Rアルゴリズムは、フィードフォワード型であり、また、周波数誤差が大きくなるに連れ、受信フィルタであるマッチフィルタ１（図２）において、受信信号r(t)が減衰し、歪みが生じる。その結果、周波数誤差の推定の精度が劣化することがある。
【００５８】
また、L&Rアルゴリズムでは、周波数誤差vを推定する範囲が、式（１０）で表される範囲に制限される。
【００５９】
すなわち、L&Rアルゴリズムでは、最大ラグNが大である場合には、既知シンボル系列において、遠く離れている既知シンボルに対するシンボル誤差z(k)どうしの自己相関を用いて、周波数誤差の推定値vが求められるため、ノイズやマルチパスが存在する悪い環境でも、周波数誤差を高精度で推定することができる。しかしながら、最大ラグNが大である場合には、式（１０）により、周波数誤差の推定の範囲が狭くなる。
【００６０】
一方、最大ラグNを小とすることにより、周波数誤差の推定の範囲を広くすることができるが、その反面、推定の精度が劣化する。
【００６１】
本発明は、このような状況に鑑みてなされたものであり、周波数誤差の推定の精度と範囲を向上させることができるようにするものである。
【課題を解決するための手段】
【００６２】
本発明の一側面のデータ処理装置は、搬送波をディジタル変調した変調信号をベースバンドの信号に変換することにより得られる、前記搬送波と同相のI成分と、前記搬送波と直交するQ成分とからなる受信データを処理するデータ処理装置であり、既知のシンボルである既知シンボルに対する前記受信データの誤差であるシンボル誤差の複数のラグの自己相関の総和である自己相関和を、前記複数のラグのうちの最大のラグである最大ラグを小から大に変化させて求め、その自己相関和から、前記搬送波に対する周波数誤差の推定値を求める周波数誤差検出手段と、前記周波数誤差の推定値の累積加算を行い、その累積加算により得られる累積加算値を、前記周波数誤差の最終的な推定値として出力する累積加算手段とを備える。
【００６３】
本発明の一側面のデータ処理方法は、搬送波をディジタル変調した変調信号をベースバンドの信号に変換することにより得られる、前記搬送波と同相のI成分と、前記搬送波と直交するQ成分とからなる受信データを処理するデータ処理方法であり、既知のシンボルである既知シンボルに対する前記受信データの誤差であるシンボル誤差の複数のラグの自己相関の総和である自己相関和を、前記複数のラグのうちの最大のラグである最大ラグを小から大に変化させて求め、その自己相関和から、前記搬送波に対する周波数誤差の推定値を求め、前記周波数誤差の推定値の累積加算を行い、その累積加算により得られる累積加算値を、前記周波数誤差の最終的な推定値として出力するステップを含む。
【００６４】
以上のような本発明の一側面においては、既知のシンボルである既知シンボルに対する前記受信データの誤差であるシンボル誤差の複数のラグの自己相関の総和である自己相関和が、前記複数のラグのうちの最大のラグである最大ラグを小から大に変化させて求められ、その自己相関和から、前記搬送波に対する周波数誤差の推定値が求められる。そして、前記周波数誤差の推定値の累積加算が行われ、その累積加算により得られる累積加算値が、前記周波数誤差の最終的な推定値として出力される。
【発明の効果】
【００６５】
本発明の一側面によれば、周波数誤差を推定することができ、特に、周波数誤差の推定の精度と範囲を向上させることができる。
【発明を実施するための最良の形態】
【００６６】
以下に本発明の実施の形態を説明するが、本発明の構成要件と、明細書又は図面に記載の実施の形態との対応関係を例示すると、次のようになる。この記載は、本発明をサポートする実施の形態が、明細書又は図面に記載されていることを確認するためのものである。従って、明細書又は図面中には記載されているが、本発明の構成要件に対応する実施の形態として、ここには記載されていない実施の形態があったとしても、そのことは、その実施の形態が、その構成要件に対応するものではないことを意味するものではない。逆に、実施の形態が構成要件に対応するものとしてここに記載されていたとしても、そのことは、その実施の形態が、その構成要件以外の構成要件には対応しないものであることを意味するものでもない。
【００６７】
本発明の一側面のデータ処理装置は、
搬送波をディジタル変調した変調信号をベースバンドの信号に変換することにより得られる、前記搬送波と同相のI成分と、前記搬送波と直交するQ成分とからなる受信データを処理するデータ処理装置（例えば、図３の受信装置）であり、
既知のシンボルである既知シンボルに対する前記受信データの誤差であるシンボル誤差の複数のラグの自己相関の総和である自己相関和を、前記複数のラグのうちの最大のラグである最大ラグを小から大に変化させて求め、その自己相関和から、前記搬送波に対する周波数誤差の推定値を求める周波数誤差検出手段（例えば、図９の周波数誤差検出装置１５０や、図１５の周波数誤差検出装置２００）と、
前記周波数誤差の推定値の累積加算を行い、その累積加算により得られる累積加算値を、前記周波数誤差の最終的な推定値として出力する累積加算手段（例えば、図９や図１５の累積加算装置１５４）と
を備える。
【００６８】
一側面のデータ処理装置には、
前記周波数誤差の最終的な推定値に基づいて、周波数誤差を補正する周波数補正手段（例えば、図９や図１５の周波数補正部１０３）をさらに設けることができる。
【００６９】
前記周波数誤差検出手段には、
既知シンボルの前記受信データと、前記既知シンボルの複素共役とを乗算することにより、前記既知シンボルに対する前記受信データのシンボル誤差を求める誤差演算手段（例えば、図１３の誤差演算部１８１）と、
既知シンボルの系列である既知シンボル系列について、前記シンボル誤差の複数のラグの自己相関の総和である自己相関和を求める相関演算手段（例えば、図１３の相関演算部１８２）と、
複数の既知シンボル系列それぞれについて、前記複数のラグの自己相関和から、前記周波数誤差の第１の推定値を求める周波数誤差演算手段（例えば、図１３の周波数誤差演算部１８３）と、
複数の既知シンボル系列それぞれについて求められた周波数誤差の第１の推定値を平均化することにより、前記周波数誤差の第２の推定値を求める推定値平均化手段（例えば、図１３の平均部１８４）と
を設け、
前記累積加算手段は、前記周波数誤差の第２の推定値の累積加算を行い、その累積加算により得られる累積加算値を、前記周波数誤差の最終的な推定値として出力することができる。
【００７０】
前記相関演算手段には、
前記シンボル誤差の所定のラグの自己相関を求める複数の相関部（例えば、図１３の相関部１８５1ないし１８５N）と、
前記複数の相関部で求められた自己相関を加算することにより、既知シンボル系列についての、前記複数のラグの自己相関和を求める演算部（例えば、図１３の演算部１８６）と
を設け、
一側面のデータ処理装置には、
前記複数の相関部のうちの、前記演算部での加算の対象とする自己相関を求める相関部を指定する制御手段（例えば、図９の制御装置１５８）をさらに設けることができる。
【００７１】
前記周波数誤差検出手段には、
既知シンボルの前記受信データと、前記既知シンボルの複素共役とを乗算することにより、前記既知シンボルに対する前記受信データのシンボル誤差を求める誤差演算手段（例えば、図１４の誤差演算部１９１）と、
既知シンボルの系列である既知シンボル系列について、前記シンボル誤差の複数のラグの自己相関の総和である自己相関和を求める相関演算手段（例えば、図１４の相関演算部１９２）と、
複数の既知シンボル系列それぞれについて求められた、前記複数のラグの自己相関和を平均化し、平均自己相関和を求める自己相関和平均化手段（例えば、図１４の平均部１９３）と、
前記平均自己相関和から、前記周波数誤差の推定値を求める周波数誤差演算手段（例えば、図１４の周波数誤差演算部１９４）と
を設けることができる。
【００７２】
前記相関演算手段には、
前記シンボル誤差の所定のラグの自己相関を求める複数の相関部（例えば、図１４の相関部１９５1ないし１９５N）と、
前記複数の相関部で求められた自己相関を加算することにより、既知シンボル系列についての、前記複数のラグの自己相関和を求める演算部（例えば、図１４の演算部１９６）と
を設け、
一側面のデータ処理装置には、
前記複数の相関部のうちの、前記演算部での加算の対象とする自己相関を求める相関部を指定する制御手段（例えば、図９の制御装置１５８）をさらに設けることができる。
【００７３】
本発明の一側面のデータ処理装置は、
搬送波をディジタル変調した変調信号をベースバンドの信号に変換することにより得られる、前記搬送波と同相のI成分と、前記搬送波と直交するQ成分とからなる受信データを処理するデータ処理方法であり、
既知のシンボルである既知シンボルに対する前記受信データの誤差であるシンボル誤差の複数のラグの自己相関の総和である自己相関和を、前記複数のラグのうちの最大のラグである最大ラグを小から大に変化させて求め、その自己相関和から、前記搬送波に対する周波数誤差の推定値を求め（例えば、図１２のステップＳ４３）、
前記周波数誤差の推定値の累積加算を行い、その累積加算により得られる累積加算値を、前記周波数誤差の最終的な推定値として出力する（例えば、図１２のステップＳ４４）
ステップを含む。
【００７４】
以下、図面を参照して、本発明の実施の形態について説明する。
【００７５】
図３は、本発明を適用した受信装置の一実施の形態の構成例を示している。
【００７６】
図３において、受信装置は、アンテナ１０１，A/D変換部１０２、周波数補正部１０３、クロック同期部１０４、受信フィルタ１０５、等化部１０６、位相同期部１０７、位相補正部１０８、誤り訂正部１０９、及び周波数同期部１１０等から構成される。
【００７７】
アンテナ１０１は、例えば、BSディジタル放送の放送波を受信し、IF(Intermediate Frequency)信号となった変調信号、すなわち、BSディジタル放送を行う放送局において、搬送波をディジタル変調した変調信号のIF信号を、A/D変換部１０２に供給する。
【００７８】
A/D変換部１０２は、アンテナ１０１からの変調信号をA/D変換し、周波数補正部１０３に供給する。
【００７９】
周波数補正部１０３は、A/D変換部１０２からの変調信号に、BSディジタル放送の送信装置で用いられる搬送波に相当する搬送波相当信号を乗算することにより、変調信号を、I成分とQ成分とからなるベースバンドの復調信号に変換し、クロック同期部１０４に供給する。ここで、周波数補正部１０３がクロック同期部１０４に供給する復調信号が、前述の式（１）の復調信号r(t)のディジタル値（サンプル値）に相当する。
【００８０】
また、周波数補正部１０３には、周波数同期部１１０から、搬送波相当信号の、BSディジタル放送の送信装置で用いられる搬送波に対する周波数誤差の推定値vが供給される。周波数補正部１０３は、周波数同期部１１０からの周波数誤差の推定値vに基づき、搬送波相当信号の、送信装置で用いられる搬送波に対する周波数誤差を補正する。
【００８１】
クロック同期部１０４は、周波数補正部１０３からの復調信号に基づき、フレーム同期等をとり、受信フィルタ１０５に供給する。
【００８２】
受信フィルタ１０５は、前述の図２のマッチフィルタ１に相当し、クロック同期部１０４からの復調信号をフィルタリングすることにより、復調信号を波形整形し、その結果得られる受信データを、等化部１０６及び周波数同期部１１０に供給する。
【００８３】
ここで、受信フィルタ１０５が等化部１０６及び周波数同期部１１０に供給する受信データが、前述の式（２）の受信データy(k)に相当する。
【００８４】
等化部１０６は、受信フィルタ１０５からの受信データを対象に等化処理を行い、位相同期部１０７、及び位相補正部１０８に供給する。
【００８５】
位相同期部１０７は、受信データとしてのシンボル系列を構成するシンボルのIQ平面上の信号点の、そのシンボルの本来の信号点に対する位相誤差を推定し、その結果得られる位相誤差の推定値を、位相補正部１０８に供給する。
【００８６】
位相補正部１０８は、位相同期部１０７からの位相誤差の推定値に基づき、等化部１０６からの受信データとしてのシンボル系列を構成するシンボルの位相誤差を補正し、誤り訂正部１０９に供給する。
【００８７】
誤り訂正部１０９は、位相補正部１０８からの受信データに対して、ビタビ復号、及びリードソロモン復号等の誤り訂正処理(FEC(Forward Error Correction))を施し、その結果得られるFECデータを出力する。このFECデータは、例えば、図示せぬMPEGのデコーダに供給されてデコードされる。
【００８８】
一方、周波数同期部１１０には、上述したように、受信フィルタ１０５から、前述の式（２）で表される受信データy(k)が供給される。
【００８９】
周波数同期部１１０は、周波数誤差を検出する周波数誤差検出機能を有し、受信フィルタ１０５からの受信データy(k)を用いて、搬送波相当信号の、BSディジタル放送の送信装置で用いられる搬送波に対する周波数誤差を推定し、その結果得られる推定値vを、周波数補正部１０３に供給する。
【００９０】
上述したように、周波数補正部１０３では、周波数同期部１１０からの周波数誤差の推定値vに基づき、その推定値vだけ、周波数誤差の補正が行われる。
【００９１】
なお、図３において、周波数補正部１０３、クロック同期部１０４、受信フィルタ１０５、及び周波数同期部１１０が、AFC装置を構成している。
【００９２】
次に、図４は、図３の周波数同期部１１０の第１の構成例を示している。
【００９３】
なお、図中、図２の周波数誤差検出装置３と対応する部分については、同一の符号を付してあり、以下では、その説明は適宜省略する。すなわち、図４の周波数同期部１１０は、誤差演算部１１、相関演算部１２、及び周波数誤差演算部１３を有する点で、図２の周波数誤差検出装置３と共通し、平均部１２１及び制御部１２２をさらに有する点で、図２の周波数誤差検出装置３と相違する。
【００９４】
平均部１２１は、周波数誤差演算部１３から供給される周波数誤差の第１の推定値を平均化することにより、周波数誤差の第２の推定値を求め、周波数補正部１０３（図３）に供給する。
【００９５】
制御部１２２は、相関演算部１２に対して、最大ラグD(≦N)を表す最大ラグ情報を供給することにより、相関演算部１２が有するN個の相関部３１1ないし３１Nのうちの、演算部３２での加算の対象となる自己相関、すなわち、式（５）左辺の総和（サメーション）の対象となる自己相関R(m)を求める相関部３１1ないし３１Dを指定する。
【００９６】
また、制御部１２２は、平均部１２１に対して、平均化の対象とする第１の推定値の個数である平均数Pを表す平均数情報を供給することにより、その平均数P、すなわち、平均化の対象とする第１の推定値の個数を指定する。
【００９７】
次に、図５を参照して、図４の周波数同期部１１０の処理について説明する。
【００９８】
図５は、式（２）で表される受信データy(k)の時系列を示している。
【００９９】
受信データy(k)には、図５に示すように、複数の既知シンボル系列が含まれる。いま、受信データy(k)におけるi番目の既知シンボル系列を、第i既知シンボル系列ということとするとともに、第i既知シンボル系列（の受信データ）の、式（３）で表されるシンボル誤差z(k)の、ラグがmの自己相関を、Ri(m)と表すこととすると、第i既知シンボル系列について求められる自己相関Ri(m)は、式（１１）で表される。
【０１００】
【数１１】
2008167230.tif 000013
・・・（１１）
【０１０１】
ここで、式（１１）において、Dは、制御部１２２（図４）が最大ラグ情報によって指定する最大ラグを表し、既知シンボル系列長L0未満の整数（L0-1以下の整数）である。なお、図４において、相関演算部１２を構成する相関部３１1ないし３１Nの数Nは、例えば、取り得る値の最大値であるL0-1であるとする。
【０１０２】
また、式（１１）において、Pは、制御部１２２が平均数情報によって指定する平均数を表す。
【０１０３】
図４において、相関演算部１２では、制御部１２２からの最大ラグ情報によって指定される最大ラグDだけの個数の相関部３１1ないし３１Dにおいて、それぞれ、第i既知シンボル系列についての、ラグmが1ないしDの自己相関Ri(1)ないしRi(D)が求められ、演算部３２において、その自己相関Ri(1)ないしRi(D)の総和である自己相関和ΣRi(m)(=Ri(1)+Ri(2)+・・・+Ri(D))が求められる。
【０１０４】
相関演算部１２は、第i既知シンボル系列についての自己相関和ΣRi(m)を求めると、その自己相関和ΣRi(m)を、周波数誤差演算部１３に供給する。
【０１０５】
周波数誤差演算部１３は、相関演算部１２から供給される、第i既知シンボル系列についての自己相関和ΣRi(m)から、式（８）にしたがい、周波数誤差の推定値を求め、平均部１２１に供給する。
【０１０６】
ここで、周波数誤差演算部１３は、第i既知シンボル系列についての自己相関和ΣRi(m)から、式（８）にしたがって求める周波数誤差の推定値を、第１の推定値といい、viと表す。
【０１０７】
平均部１２１は、周波数誤差演算部１３から供給される周波数誤差の第１の推定値viを一時記憶し、制御部１２２からの平均数情報が表す平均数Pだけの第１の推定値v1,v2,・・・,vPを記憶すると、そのP個の第１の推定値v1,v2,・・・,vP、つまり、P個の既知シンボル系列それぞれについて求められた周波数誤差の第１の推定値v1,v2,・・・,vPを、式v=(v1+v2+・・・+vP)/Pにしたがって平均化することにより、周波数補正部１０３（図３）での補正に用いられる、周波数誤差の推定値vを求める。
【０１０８】
ここで、平均部１２１において、P個の既知シンボル系列それぞれについて求められた周波数誤差の第１の推定値v1ないしvPを平均化することにより求められる周波数誤差の推定値vを、第２の推定値という。
【０１０９】
第２の推定値vは、結局、式（１２）にしたがって求められることになる。
【０１１０】
【数１２】
2008167230.tif 000014
・・・（１２）
【０１１１】
次に、図６のフローチャートを参照して、図４の周波数同期部１１０の処理について、さらに説明する。
【０１１２】
制御部１２２は、最大ラグ情報を相関演算部１２に供給するとともに、平均数情報を平均部１２１に供給する。
【０１１３】
ここで、最大ラグ情報が表す最大ラグD、及び平均数情報が表す平均数Pは、例えば、ユーザの操作や、C/N、誤り訂正部１０９（図３）で得られるFECデータの誤り率等に応じて決定することができる。すなわち、例えば、C/NやFECデータの誤り率が悪い場合には、最大ラグD、及び平均数Pを大とすることができ、逆に、C/NやFECデータの誤り率が良い場合には、最大ラグD、及び平均数Pを小とすることができる。
【０１１４】
ステップＳ１１において、誤差演算部１１は、受信フィルタ１０５（図３）から供給される式（２）の、既知シンボルc(k)の受信データy(k)と、既知シンボルc(k)の複素共役c*(k)とを乗算することにより、式（３）のシンボル誤差z(k)を求め、処理は、ステップＳ１２に進む。
【０１１５】
ステップＳ１２では、相関演算部１２は、誤差演算部１１からのシンボル誤差z(k)を用い、制御部１２２からの最大ラグ情報によって指定される最大ラグDだけの個数の相関部３１1ないし３１Dにおいて、それぞれ、既知シンボル系列についての、ラグmが1ないしDの自己相関Ri(1)ないしRi(D)を求める。さらに、ステップＳ１２では、相関演算部１２は、演算部３２において、相関部３１1ないし３１Dでそれぞれ求められた自己相関Ri(1)ないしRi(D)の総和である自己相関和ΣRi(m)を求め、周波数誤差演算部１３に供給する。
【０１１６】
なお、相関演算部１２は、P個の既知シンボル系列である第１既知シンボル系列ないし第P既知シンボル系列それぞれについての自己相関和ΣR1(m)ないしΣRP(m)を求め、周波数誤差演算部１３に供給する。
【０１１７】
ステップＳ１２の後、処理は、ステップＳ１３に進み、周波数誤差演算部１３は、相関演算部１２から供給される、P個の第１既知シンボル系列ないし第P既知シンボル系列それぞれについて求められた自己相関和ΣR1(m)ないしΣRP(m)から、式（８）にしたがい、P個の周波数誤差の第１の推定値v1ないしvPを求め、平均部１２１に供給して、処理は、ステップＳ１４に進む。
【０１１８】
ステップＳ１４では、平均部１２１は、周波数誤差演算部１３から供給されるP個の周波数誤差の第１の推定値v1ないしvPを、式v=(v1+v2+・・・+vP)/Pにしたがって平均化することにより、周波数誤差の第２の推定値vを求めて、周波数補正部１０３に供給する。
【０１１９】
以上のように、複数であるP個の既知シンボル系列それぞれについて求められた周波数誤差の第１の推定値v1ないしvPを平均化することにより、周波数誤差の第２の推定値vを求めるので、既知シンボル系列長L0が長い既知シンボル系列を用いて求めた周波数誤差の推定値と等価な第２の推定値vを求めることができ、したがって、既知シンボル系列長L0が制限されていても、周波数誤差の推定の精度を向上させることができる。また、周波数誤差の推定の結果のばらつきを抑制することができる。
【０１２０】
次に、図７は、図３の周波数同期部１１０の第２の構成例を示している。
【０１２１】
なお、図中、図２の周波数誤差検出装置３と対応する部分については、同一の符号を付してあり、以下では、その説明は適宜省略する。すなわち、図７の周波数同期部１１０は、誤差演算部１１、及び相関演算部１２を有する点で、図２の周波数誤差検出装置３と共通し、平均部１３１、及び制御部１３３をさらに有する点、並びに、周波数誤差演算部１３に代えて周波数誤差演算部１３２が設けられている点で、図２の周波数誤差検出装置３と相違する。
【０１２２】
平均部１３１は、周波数誤差演算部１３から供給される自己相関和を、制御部１３３から供給される平均数情報が表す平均数Pによって指定される個数だけ平均化することにより、平均自己相関和を求めて、周波数誤差演算部１３２に供給する。
【０１２３】
周波数誤差演算部１３２は、自己相関和に代えて、平均部１３１から供給される平均自己相関和を用いて、周波数誤差の推定値vを求め、周波数補正部１０３（図３）に供給する。
【０１２４】
制御部１３３は、図４の制御部１２２と同様に、相関演算部１２に対して、最大ラグD(≦N)を表す最大ラグ情報を供給することにより、相関演算部１２が有するN個の相関部３１1ないし３１Nのうちの、演算部３２での加算の対象となる自己相関、すなわち、式（５）左辺の総和（サメーション）の対象となる自己相関R(m)を求める相関部３１1ないし３１Dを指定する。
【０１２５】
また、制御部１３３は、平均部１３１に対して、平均化の対象とする自己相関和の個数である平均数Pを表す平均数情報を供給することにより、平均数P、すなわち、平均化の対象とする自己相関和の個数を指定する。
【０１２６】
以上のように構成される周波数同期部１１０では、相関演算部１２が、図４で説明した場合と同様に、制御部１３３からの最大ラグ情報によって指定される最大ラグDだけの個数の相関部３１1ないし３１Dにおいて、上述の式（１１）にしたがい、それぞれ、第i既知シンボル系列についての、ラグmが1ないしDの自己相関Ri(1)ないしRi(D)が求められ、演算部３２において、その自己相関Ri(1)ないしRi(D)の総和である自己相関和ΣRi(m)が求められる。
【０１２７】
相関演算部１２は、第i既知シンボル系列についての自己相関和ΣRi(m)を求めると、その自己相関和ΣRi(m)を、平均部１３１に供給する。
【０１２８】
平均部１３１は、相関演算部１２から供給される自己相関和ΣRi(m)を一時記憶し、制御部１３３からの平均数情報が表す平均数Pだけの自己相関和ΣR1(m),ΣR2(m),・・・,ΣRP(m)を記憶すると、そのP個の自己相関和ΣR1(m),ΣR2(m),・・・,ΣRP(m)、つまり、P個の既知シンボル系列それぞれについて求められた自己相関和ΣR1(m),ΣR2(m),・・・,ΣRP(m)を、式R=(ΣR1(m),ΣR2(m),・・・,ΣRP(m))/Pにしたがって平均化することにより、平均自己相関和Rを求めて、周波数誤差演算部１３２に供給する。
【０１２９】
周波数誤差演算部１３２は、１つの既知シンボル系列についての自己相関和ΣR(m)に代えて、平均部１３１から供給される平均自己相関和Rを用い、式（８）にしたがって、周波数誤差の推定値vを求める。
【０１３０】
したがって、図７の周波数同期部１１０では、周波数誤差の推定値vが、式（１３）にしたがって求められる。
【０１３１】
【数１３】
2008167230.tif 000015
・・・（１３）
【０１３２】
次に、図８のフローチャートを参照して、図７の周波数同期部１１０の処理について説明する。
【０１３３】
制御部１３３は、最大ラグ情報を相関演算部１２に供給するとともに、平均数情報を平均部１３１に供給する。
【０１３４】
ここで、最大ラグ情報が表す最大ラグD、及び平均数情報が表す平均数Pは、図６の場合と同様に、ユーザの操作や、C/N、誤り訂正部１０９（図３）で得られるFECデータの誤り率等に応じて決定することができる。
【０１３５】
ステップＳ２１において、誤差演算部１１は、受信フィルタ１０５（図３）から供給される式（２）の、既知シンボルc(k)の受信データy(k)と、既知シンボルc(k)の複素共役c*(k)とを乗算することにより、式（３）のシンボル誤差z(k)を求め、処理は、ステップＳ２２に進む。
【０１３６】
ステップＳ２２では、相関演算部１２は、誤差演算部１１からのシンボル誤差z(k)を用い、制御部１３３からの最大ラグ情報によって指定される最大ラグDだけの個数の相関部３１1ないし３１Dにおいて、それぞれ、既知シンボル系列についての、ラグmが1ないしDの自己相関Ri(1)ないしRi(D)を求める。さらに、ステップＳ２２では、相関演算部１２は、演算部３２において、相関部３１1ないし３１Dでそれぞれ求められた自己相関Ri(1)ないしRi(D)の総和である自己相関和ΣRi(m)を求め、平均部１３１に供給する。
【０１３７】
なお、相関演算部１２は、P個の既知シンボル系列である第１既知シンボル系列ないし第P既知シンボル系列それぞれについての自己相関和ΣR1(m)ないしΣRP(m)を求め、平均部１３１に供給する。
【０１３８】
ステップＳ２２の後、処理は、ステップＳ２３に進み、平均部１３１は、相関演算部１２から供給される、P個の第１既知シンボル系列ないし第P既知シンボル系列それぞれについて求められた、P個の自己相関和ΣR1(m)ないしΣRP(m)を、式R=(ΣR1(m),ΣR2(m),・・・,ΣRP(m))/Pにしたがって平均化することにより、平均自己相関和Rを求め、周波数誤差演算部１３２に供給して、処理は、ステップＳ２４に進む。
【０１３９】
ステップＳ２４では、周波数誤差演算部１３２は、平均部１３１からの平均自己相関和Rから、周波数誤差の推定値vを求め、すなわち、自己相関和ΣR(m)に代えて、平均自己相関和Rを用いて、式（８）を計算することにより、周波数誤差の推定値vを求め、周波数補正部１０３に供給する。
【０１４０】
以上のように、複数であるP個の既知シンボル系列それぞれについて求められた自己相関和ΣR1(m)ないしΣRP(m)を平均化することにより、平均自己相関和Rを求め、さらに、その平均自己相関和Rから、周波数誤差の推定値vを求めるので、図４の場合と同様に、既知シンボル系列長L0が長い既知シンボル系列を用いて求めた周波数誤差の推定値と等価な周波数誤差の推定値vを求めることができ、したがって、既知シンボル系列長L0が制限されていても、周波数誤差の推定の精度を向上させることができる。また、周波数誤差の推定の結果のばらつきを抑制することができる。
【０１４１】
次に、図９は、図３の周波数同期部１１０の第３の構成例を示している。
【０１４２】
図９では、周波数同期部１１０は、周波数誤差検出装置１５０、累積加算装置１５４、及び、制御装置１５８から構成されている。
【０１４３】
周波数誤差検出装置１５０は、スイッチ１５１、Q個の検出部１５２1ないし１５２Q、及びスイッチ１５３から構成され、式（３）のシンボル誤差z(k)の、複数のラグの自己相関R(m)の総和である自己相関和ΣR(m)を、複数のラグのうちの最大ラグDを小から大に変化させて求め、その自己相関和ΣR(m)から、周波数誤差の推定値vを求めて、累積加算装置１５４に供給する。
【０１４４】
すなわち、スイッチ１５１は、制御装置１５８の制御にしたがって、Q個の検出部１５２1ないし１５２Qのうちの、１の検出部１５２q(q=1,2,・・・,Q)を選択し、その検出部１５２qに、受信フィルタ１０５から供給される受信データy(k)を供給する。
【０１４５】
検出部１５２qは、スイッチ１５１から供給される、既知シンボル系列の受信データy(k)から、例えば、L&Rアルゴリズムにより、周波数誤差の推定値v(q)を求めて出力する。
【０１４６】
スイッチ１５３は、制御装置１５８の制御にしたがって、Q個の検出部１５２1ないし１５２Qのうちの、スイッチ１５１が選択するのと同一の１の検出部１５２qを選択し、その検出部１５２qが出力する周波数誤差の推定値v(q)を、累積加算装置１５４に供給する。
【０１４７】
累積加算装置１５４は、演算部１５５、セレクタ１５６、及び記憶部１５７から構成され、周波数誤差検出装置１５０から供給される周波数誤差の推定値v(q)の累積加算を行い、その累積加算により得られる累積加算値を、周波数誤差の最終的な推定値vとして出力する。
【０１４８】
すなわち、演算部１５５は、記憶部１５７の記憶値vmemに、周波数誤差検出装置１５０から供給される周波数誤差の推定値v(q)を加算し、その結果得られる加算値vmem+v(q)をセレクタ１５６に供給する。
【０１４９】
セレクタ１５６には、演算部１５５から加算値vmem+v(q)が供給される他、記憶部１５７から記憶値vmemが供給されるようになっている。セレクタ１５６は、制御装置１５８の制御にしたがい、演算部１５５からの加算値vmem+v(q)、又は、記憶部１５７からの記憶値vmemを選択し、周波数誤差の最終的な推定値vとして、記憶部１５７と、周波数補正部１０３とに供給する。
【０１５０】
記憶部１５７は、初期値を0として、セレクタ１５６から供給される周波数誤差の最終的な推定値vを、記憶値vmemとして、上書きの形で記憶する。
【０１５１】
制御装置１５８は、周波数誤差検出装置１５０を構成するスイッチ１５１及び１５３を制御する。また、制御装置１５８は、周波数誤差検出装置１５０を構成する検出部１５２1ないし１５２Qを制御する。さらに、制御装置１５８は、累積加算装置１５４を構成するセレクタ１５６を制御する。
【０１５２】
以上のように構成される周波数同期部１１０では、制御装置１５８は、周波数誤差検出装置１５０の検出部１５２1から検出部１５２Qを順次選択するように、スイッチ１５１及び１５３を制御し、これにより、周波数誤差検出装置１５０から累積加算装置１５４には、検出部１５２1ないし検出部１５２Qで順次得られる周波数誤差の推定値v(1),v(2),・・・,v(Q)が、順次供給される。
【０１５３】
また、制御装置１５８は、周波数誤差検出装置１５０の検出部１５２1から検出部１５２Qを順次選択するように、スイッチ１５１及び１５３を制御している間は、演算部１５５の加算値vmem+v(q)を選択するように、セレクタ１５６を制御し、検出部１５２Qで求められた周波数誤差の推定値v(Q)と記憶部１５７の記憶値vmemとを加算した加算値vmem+v(q)を選択した後は、記憶部１５７の記憶値vmemを選択するように、セレクタ１５６を制御する。
【０１５４】
したがって、制御装置１５８は、まず、検出部１５２1を選択するように、スイッチ１５１及び１５３を制御する。この場合、検出部１５２1において、周波数誤差の推定値v(1)が求められ、演算部１５５に供給される。
【０１５５】
また、演算部１５５には、記憶部１５７の記憶値vmemが供給されるが、制御装置１５８は、検出部１５２1を選択するように、スイッチ１５１及び１５３を制御するとき、記憶部１５７の記憶値vmemを、初期値である0に初期化する。したがって、いまの場合、演算部１５５では、記憶部１５７の記憶値vmemである0に、検出部１５２1で求められた周波数誤差の推定値v(1)が加算され、その加算値vmem+v(1)=v(1)が、周波数誤差の最終的な推定値vとして、セレクタ１５６を介して、記憶部１５７と周波数補正部１０３に供給される。
【０１５６】
そして、記憶部１５７では、その周波数誤差の最終的な推定値v=v(1)が、記憶値vmemとして記憶される。
【０１５７】
その後、制御装置１５８は、検出部１５２2を選択するように、スイッチ１５１及び１５３を制御する。この場合、検出部１５２2において、周波数誤差の推定値v(2)が求められ、演算部１５５に供給される。
【０１５８】
また、演算部１５５には、記憶部１５７の記憶値vmem=v(1)が供給される。演算部１５５では、記憶部１５７の記憶値vmemに、検出部１５２2で求められた周波数誤差の推定値v(2)が加算され、その加算値vmem+v(2)=v(1)+v(2)、すなわち、周波数誤差の推定値v(1)及びv(2)の累積加算値が、周波数誤差の最終的な推定値vとして、セレクタ１５６を介して、記憶部１５７と周波数補正部１０３に供給される。
【０１５９】
そして、記憶部１５７では、その周波数誤差の最終的な推定値v=v(1)+v(2)が、記憶値vmemとして記憶される。
【０１６０】
その後、制御装置１５８は、検出部１５２3を選択するように、スイッチ１５１及び１５３を制御する。この場合、検出部１５２3において、周波数誤差の推定値v(3)が求められ、演算部１５５に供給される。
【０１６１】
また、演算部１５５には、記憶部１５７の記憶値vmem=v(1)+v(2)が供給される。演算部１５５では、記憶部１５７の記憶値vmemに、検出部１５２3で求められた周波数誤差の推定値v(3)が加算され、その加算値vmem+v(3)=v(1)+v(2)+v(3)、すなわち、周波数誤差の推定値v(1),v(2),v(3)の累積加算値が、周波数誤差の最終的な推定値vとして、セレクタ１５６を介して、記憶部１５７と周波数補正部１０３に供給される。
【０１６２】
以下、同様にして、累積加算装置１５４では、周波数誤差の推定値v(1),v(2),・・・,v(Q)が、順次、累積加算される。
【０１６３】
そして、周波数誤差の推定値v(1),v(2),・・・,v(Q)のすべての累積加算をした累積加算値v(1)+v(2)+・・・+v(Q)が、周波数誤差の最終的な推定値vとして、セレクタ１５６を介して、記憶部１５７と周波数補正部１０３に供給され、記憶部１５７において、その周波数誤差の最終的な推定値v=v(1)+v(2)+・・・+v(Q)が、記憶値vmemとして記憶されると、制御装置１５８は、記憶部１５７の記憶値vmemを選択するように、セレクタ１５６を制御する。
【０１６４】
セレクタ１５６は、制御装置１５８の制御にしたがい、記憶部１５７の記憶値vmem=v(1)+v(2)+・・・+v(Q)を選択し、周波数誤差の最終的な推定値vとして、記憶部１５７と周波数補正部１０３に供給する。
【０１６５】
記憶部１５７は、セレクタ１５６からの周波数誤差の最終的な推定値v=v(1)+v(2)+・・・+v(Q)を、記憶値vmemとして記憶し、したがって、以降は、記憶部１５７の記憶値vmem=v(1)+v(2)+・・・+v(Q)が、周波数誤差の最終的な推定値vとして、周波数補正部１０３に供給される。
【０１６６】
図１０は、図９の検出部１５２qの第１の構成例を示している。
【０１６７】
検出部１５２qは、誤差演算部１７１、相関演算部１７２、及び周波数誤差演算部１７３から構成され、上述したように、スイッチ１５１（図９）から供給される、既知シンボル系列の受信データy(k)から、L&Rアルゴリズムにより、周波数誤差の推定値v(q)を求めて出力する。
【０１６８】
すなわち、誤差演算部１７１は、図２の誤差演算部１１と同様に、スイッチ１５１からの既知シンボルc(k)の受信データy(k)と、既知シンボルc(k)の複素共役c*(k)とを乗算することにより、式（３）で表されるシンボル誤差z(k)を求め、相関演算部１７２に供給する。
【０１６９】
相関演算部１７２は、N個の相関部１７４1，１７４2，・・・，１７４Nと、演算部１７５とから構成され、図２の相関演算部１２と同様に、誤差演算部１７１からのシンボル誤差z(k)を用い、式（５）の自己相関和ΣR(m)を求めて、周波数誤差演算部１７３に供給する。
【０１７０】
すなわち、相関部１７４mは、式（４）にしたがい、既知シンボル系列について、誤差演算部１７１から供給されるシンボル誤差z(k)の、ラグがmの自己相関R(m)を求め、演算部１７５に供給する。
【０１７１】
なお、図９の制御装置１５８は、相関演算部１７２に対して、最大ラグDq(≦N)を表す最大ラグ情報を供給するようになっており、これにより、N個の相関部１７４1ないし１７４Nのうちの、演算部１７５での加算の対象となる自己相関、すなわち、式（５）左辺の総和（サメーション）の対象となる自己相関R(m)を求める相関部１７４1ないし１７４Dqを指定する。
【０１７２】
相関演算部１７２では、N個の相関部１７４1ないし１７４Nのうちの、制御装置１５８からの最大ラグ情報が表す最大ラグDqの数の相関部１７４1ないし１７４Dqにおいて、それぞれ、自己相関R(1),R(2),・・・,R(Dq)を求め、演算部１７５に供給する。
【０１７３】
演算部１７５は、N個の相関部１７４1ないし１７４Dqそれぞれから供給される自己相関R(1)ないしR(Dq)の総和である、式（５）の自己相関和ΣR(m)を求め、周波数誤差演算部１７３に供給する。
【０１７４】
周波数誤差演算部１７３は、演算部１７５からの自己相関和ΣR(m)を用い、式（８）にしたがって、周波数誤差の推定値v(q)を求めて出力する。
【０１７５】
次に、図１１を参照して、図９の周波数同期部１１０の処理について説明する。
【０１７６】
なお、ここでは、説明を簡単にするために、周波数誤差検出装置１５０が、２つの検出部１５２1及び１５２2から構成されることとする。
【０１７７】
図１１は、周波数誤差検出装置１５０を構成する検出部１５２1と１５２2のSカーブ、すなわち、周波数誤差の推定特性を示している。
【０１７８】
図１１において、LRは、理想的な周波数誤差の推定特性を、L1は、検出部１５２1の周波数誤差の推定特性を、L2は、検出部１５２2の周波数誤差の推定特性を、それぞれ表している。
【０１７９】
ここで、図１１において、横軸は、周波数誤差の真値を表し、縦軸は、周波数誤差の真値に対して推定される、周波数誤差の推定値を示している。
【０１８０】
制御装置１５８は、上述したように、検出部１５２q（図１０）を構成する相関演算部１７２に対して、最大ラグDq(≦N)を表す最大ラグ情報を供給するが、最大ラグD1,D2,・・・,DQは、式D1<D2<・・・<DQを満たすようになっている。
【０１８１】
ここで、L&Rアルゴリズムでは、前述したように、最大ラグDqが大である場合には、既知シンボル系列において、遠く離れている既知シンボルに対するシンボル誤差z(k)どうしの自己相関を用いて、周波数誤差の推定値vが求められるため、ノイズやマルチパスが存在する悪い環境でも、周波数誤差を高精度で推定することができる。しかしながら、最大ラグDqが大である場合には、前述の式（１０）において最大ラグを表すNが大となるから、周波数誤差の推定の範囲が狭くなる。
【０１８２】
一方、最大ラグDqを小とすることにより、周波数誤差の推定の範囲を広くすることができるが、その反面、推定の精度が劣化する。
【０１８３】
したがって、周波数誤差検出装置１５０が、２つの検出部１５２1と１５２2で構成される場合において、式D1<D2の関係がある最大ラグD1を表す最大ラグ情報が検出部１５２1に供給されるとともに、最大ラグD2を表す最大ラグ情報が検出部１５２2に供給されると、検出部１５２1では、図１１に示すように、推定の範囲が広いが、推定の精度が、理想的な推定特性LRに比較して劣化している推定特性L1にしたがった周波数誤差の推定値v(1)が求められる。
【０１８４】
一方、検出部１５２2では、図１１に示すように、推定の範囲が狭いが、推定の精度が、理想的な推定特性LRに近い推定特性L2にしたがった周波数誤差の推定値v(2)が求められる。
【０１８５】
上述したように、周波数誤差検出装置１５０では、式（３）のシンボル誤差z(k)の、複数のラグの自己相関R(m)の総和である自己相関和ΣR(m)を、複数のラグのうちの最大ラグDを小から大に変化させて求め、その自己相関和ΣR(m)から、周波数誤差の推定値vを求めて、累積加算装置１５４に供給する。そして、累積加算装置１５４では、周波数誤差検出装置１５０から供給される周波数誤差の推定値v(q)の累積加算が行われ、その累積加算により得られる累積加算値が、周波数誤差の最終的な推定値vとして出力される。
【０１８６】
したがって、周波数誤差検出装置１５０では、まず、検出部１５２1において、最大ラグDとして、小さい値D1を用いて、自己相関和ΣR(m)が求められ、その自己相関和ΣR(m)から、推定の精度は悪いが、推定の範囲が広い推定特性L1にしたがった周波数誤差の推定値v(1)が求められる。
【０１８７】
そして、周波数補正部１０３では、周波数誤差の推定値v(1)を、周波数誤差の最終的な推定値として、その推定値に基づいて、周波数誤差の補正が行われる。
【０１８８】
上述したように、周波数誤差の推定値v(1)は、精度が悪いので、周波数誤差の推定値v(1)に基づく周波数誤差の補正を行うだけでは、図１１に示すように、ある程度の周波数誤差△vR1が残留する。
【０１８９】
その後、周波数誤差検出装置１５０では、検出部１５２2において、最大ラグDとして、大きな値D2を用いて、自己相関和ΣR(m)が求められ、その自己相関和ΣR(m)から、推定の範囲は狭いが、推定の精度が高い推定特性L2にしたがった周波数誤差の推定値v(2)が求められる。
【０１９０】
すなわち、検出部１５２2において求められる周波数誤差の推定値v(2)は、推定の範囲は狭いが、推定の精度が高いので、図１１に示すように、周波数誤差の推定値v(1)に基づく周波数誤差の補正を行った後に残留する周波数誤差△vR1を、精度良く表す値となる。
【０１９１】
そして、累積加算装置１５４では、先に周波数誤差の補正に用いられた周波数誤差の最終的な推定値としての周波数誤差の推定値v(1)に、検出部１５２2において求められた周波数誤差の推定値v(2)が累積加算され、周波数補正部１０３では、その累積加算値v(1)+v(2)を、周波数誤差の最終的な推定値として、その推定値に基づいて、周波数誤差の補正が行われる。
【０１９２】
以上のように、図９の周波数同期部１１０では、周波数誤差検出装置１５０において、シンボル誤差z(k)の、複数のラグの自己相関R(m)の総和である自己相関和ΣR(m)を、複数のラグのうちの最大ラグDを小から大に変化させて求め、その自己相関和ΣR(m)から、周波数誤差の推定値vを求める。さらに、累積加算装置１５４において、周波数誤差検出装置１５０で求められた周波数誤差の推定値v(q)の累積加算を行い、その累積加算により得られる累積加算値を、周波数誤差の最終的な推定値vとして出力する。したがって、周波数誤差の推定の精度と範囲を向上させることができる。
【０１９３】
その結果、例えば、受信フィルタ１０５（図９）において、受信データy(k)の歪みが生じている場合や、受信環境が劣悪な場合のように、大きな周波数誤差が生じる場合であっても、その周波数誤差を、正確に補正することができる。
【０１９４】
次に、図１２のフローチャートを参照して、図９の周波数同期部１１０の処理について、さらに説明する。
【０１９５】
ステップＳ４１において、制御装置１５８は、周波数誤差検出装置１５０に設定するパラメータとして、例えば、最大ラグD1,D2,・・・,DQを、式D1<D2<・・・<DQを満たすように決定し、その最大ラグDqを表す最大ラグ情報を、検出部１５２qに供給して、処理は、ステップＳ４２に進む。
【０１９６】
ステップＳ４２では、制御装置１５８は、周波数誤差の推定値v(q)を求める検出部１５２qを指定するための変数qを1に初期化して、処理は、ステップＳ４３に進む。
【０１９７】
ステップＳ４３では、制御装置１５８は、検出部１５２qを選択するように、スイッチ１５１及び１５３を制御し、これにより、受信データy(k)が、受信フィルタ１０５から、スイッチ１５１を介して、検出部１５２qに供給される。
【０１９８】
検出部１５２qは、受信フィルタ１０５からスイッチ１５１を介して供給される、既知シンボル系列の受信データy(k)を用い、シンボル誤差z(k)の、最大ラグをDqとする複数のラグの自己相関R(m)の総和である自己相関和ΣR(m)(=R(1)+R(2)+・・・+R(Dq))を求め、その自己相関和ΣR(m)から、式（８）にしたがい、周波数誤差の推定値v(q)を求める。そして、検出部１５２qは、周波数誤差の推定値v(q)を、累積加算装置１５４に供給して、処理は、ステップＳ４３からステップＳ４４に進む。
【０１９９】
ステップＳ４４では、累積加算装置１５４において、演算部１５５が、記憶部１５７の記憶値vmemと、検出部１５２qからの周波数誤差の推定値v(q)とを加算し、その結果得られる加算値vmem+v(q)を、セレクタ１５６に供給して、処理は、ステップＳ４５に進む。
【０２００】
ステップＳ４５では、制御装置１５８が、変数qが、検出部１５２1ないし１５２Qの総数Qより大きいかどうかを判定する。
【０２０１】
ステップＳ４５において、変数qが、総数Qより大きくないと判定された場合、処理は、ステップＳ４６に進み、制御装置１５８は、演算部１５５で求められた加算値vmem+v(q)を選択するように、セレクタ１５６を制御する。
【０２０２】
セレクタ１５６は、制御装置１５８の制御にしたがい、演算部１５５で求められた加算値vmem+v(q)を選択し、周波数誤差の最終的な推定値vとして、記憶部１５７と周波数補正部１０３に供給する。記憶部１５７は、セレクタ１５６からの周波数誤差の最終的な推定値vを、記憶値vmemとして新たに記憶する。
【０２０３】
その後、処理は、ステップＳ４６からステップＳ４７に進み、制御装置１５８は、変数qを1だけインクリメントして、処理は、ステップＳ４３に戻り、次の既知シンボル系列の受信データy(k)を対象として、以下、同様の処理が繰り返される。
【０２０４】
また、ステップＳ４５において、変数qが、総数Qより大きいと判定された場合、処理は、ステップＳ４８に進み、制御装置１５８は、記憶部１５７の記憶値vmemを選択するように、セレクタ１５６を制御する。
【０２０５】
セレクタ１５６は、制御装置１５８の制御にしたがい、記憶部１５７の記憶値vmemを選択し、周波数誤差の最終的な推定値vとして、記憶部１５７と周波数補正部１０３に供給する。記憶部１５７は、セレクタ１５６からの周波数誤差の最終的な推定値vを、記憶値vmemとして新たに記憶する。
【０２０６】
その後、処理は、ステップＳ４８に戻り、以下、同様の処理が繰り返される。
【０２０７】
次に、図１３は、図９の検出部１５２qの第２の構成例を示している。
【０２０８】
図１３では、検出部１５２qは、誤差演算部１８１、相関演算部１８２、周波数誤差演算部１８３、及び平均部１８４から構成され、図４の場合と同様にして、周波数誤差の推定値v(q)を求めて出力する。
【０２０９】
すなわち、図１３において、誤差演算部１８１、相関演算部１８２、周波数誤差演算部１８３、及び平均部１８４は、それぞれ、図４の誤差演算部１１、相関演算部１２、周波数誤差演算部１３、及び平均部１２１と同様の処理を行う。
【０２１０】
具体的には、誤差演算部１８１は、受信フィルタ１０５（図９）から供給される式（２）の、既知シンボルc(k)の受信データy(k)と、既知シンボルc(k)の複素共役c*(k)とを乗算することにより、式（３）のシンボル誤差z(k)を求め、相関演算部１８２に供給する。
【０２１１】
相関演算部１８２は、N個の相関部１８５1，１８５2，・・・，１８５Nと、演算部１８６とから構成され、図２の相関演算部１２と同様に、誤差演算部１８１からのシンボル誤差z(k)を用い、式（１１）の自己相関和ΣRi(m)を求めて、周波数誤差演算部１８３に供給する。
【０２１２】
すなわち、相関部１８５mは、式（１１）にしたがい、既知シンボル系列について、誤差演算部１８１から供給されるシンボル誤差z(k)の、ラグがmの自己相関Ri(m)を求め、演算部１８６に供給する。
【０２１３】
なお、図９の制御装置１５８は、相関演算部１８２に対して、最大ラグDq(≦N)を表す最大ラグ情報を供給するようになっており、これにより、N個の相関部１８５1ないし１８５Nのうちの、演算部１８６での加算の対象となる自己相関、すなわち、式（１２）右辺の総和（サメーション）の対象となる自己相関Ri(m)を求める相関部１８５1ないし１８５Dqを指定する。
【０２１４】
相関演算部１８２では、N個の相関部１８５1ないし１８５Nのうちの、制御装置１５８からの最大ラグ情報が表す最大ラグDqの数の相関部１８５1ないし１８５Dqにおいて、それぞれ、ラグmが１,2,・・・,Dqの自己相関Ri(1),Ri(2),・・・,Ri(Dq)を求め、演算部１８６に供給する。
【０２１５】
演算部１８６は、N個の相関部１８５1ないし１８５Dqそれぞれから供給される自己相関Ri(1)ないしRi(Dq)の総和である、式（１２）の自己相関和ΣRi(m)を求め、周波数誤差演算部１８３に供給する。
【０２１６】
周波数誤差演算部１８３は、相関演算部１８２の演算部１８６から供給される自己相関和ΣRi(m)から、式（８）にしたがい、周波数誤差の第１の推定値viを求め、平均部１８４に供給する。
【０２１７】
そして、相関演算部１８２において、P個の既知シンボル系列である第１既知シンボル系列ないし第P既知シンボル系列それぞれについての自己相関和ΣR1(m)ないしΣRP(m)が求められ、さらに、周波数誤差演算部１８３において、P個の第１既知シンボル系列ないし第P既知シンボル系列それぞれについて求められた自己相関和ΣR1(m)ないしΣRP(m)から、式（８）にしたがい、P個の周波数誤差の第１の推定値v1ないしvPが求められると、平均部１８４は、そのP個の周波数誤差の第１の推定値v1ないしvPを、式v(q)=(v1+v2+・・・+vP)/Pにしたがって平均化することにより、周波数誤差の第２の推定値v(q)を求めて、累積加算装置１５４（図９）に供給する。
【０２１８】
なお、図９の制御装置１５８は、平均部１８４に対して、平均数Pを表す平均数情報を供給するようになっており、これにより、平均化の対象とする第１の推定値の個数である平均数Pを指定する。
【０２１９】
以上のように、図９の周波数同期部１１０において、シンボル誤差の、複数のラグの自己相関の総和である自己相関和を、複数のラグのうちの最大ラグを小から大に変化させて求め、その自己相関和から、周波数誤差の推定値v(q)を求めて、その推定値v(q)の累積加算を行い、その累積加算により得られる累積加算値を、周波数誤差の最終的な推定値として出力する場合において、複数であるP個の既知シンボル系列それぞれについて求められた周波数誤差の第１の推定値v1ないしvPを平均化して、累積加算の対象となる周波数誤差の第２の推定値v(q)を求めることにより、周波数誤差の推定の精度を、より向上させることができる。
【０２２０】
次に、図１４は、図９の検出部１５２qの第３の構成例を示している。
【０２２１】
図１４では、検出部１５２qは、誤差演算部１９１、相関演算部１９２、平均部１９３、及び周波数誤差演算部１９４から構成され、図７の場合と同様にして、周波数誤差の推定値v(q)を求めて出力する。
【０２２２】
すなわち、図１４において、誤差演算部１９１、相関演算部１９２、平均部１９３、及び周波数誤差演算部１９４は、それぞれ、図７の誤差演算部１１、相関演算部１２、平均部１３１、及び周波数誤差演算部１３２と同様の処理を行う。
【０２２３】
具体的には、誤差演算部１９１は、受信フィルタ１０５（図９）から供給される式（２）の、既知シンボルc(k)の受信データy(k)と、既知シンボルc(k)の複素共役c*(k)とを乗算することにより、式（３）のシンボル誤差z(k)を求め、相関演算部１９２に供給する。
【０２２４】
相関演算部１９２は、N個の相関部１９５1，１９５2，・・・，１９５Nと、演算部１９６とから構成され、図２の相関演算部１２と同様に、誤差演算部１９１からのシンボル誤差z(k)を用い、式（１１）の自己相関和ΣRi(m)を求めて、平均部１９３に供給する。
【０２２５】
すなわち、相関部１９５mは、式（１１）にしたがい、既知シンボル系列について、誤差演算部１９１から供給されるシンボル誤差z(k)の、ラグがmの自己相関Ri(m)を求め、演算部１９６に供給する。
【０２２６】
なお、図９の制御装置１５８は、相関演算部１９２に対して、最大ラグDq(≦N)を表す最大ラグ情報を供給するようになっており、これにより、N個の相関部１９５1ないし１９５Nのうちの、演算部１９６での加算の対象となる自己相関、すなわち、式（１３）右辺の総和（サメーション）の対象となる自己相関Ri(m)を求める相関部１９５1ないし１９５Dqを指定する。
【０２２７】
相関演算部１９２では、N個の相関部１９５1ないし１９５Nのうちの、制御装置１５８からの最大ラグ情報が表す最大ラグDqの数の相関部１９５1ないし１９５Dqにおいて、それぞれ、ラグmが１,2,・・・,Dqの自己相関Ri(1),Ri(2),・・・,Ri(Dq)を求め、演算部１９６に供給する。
【０２２８】
演算部１９６は、N個の相関部１９５1ないし１９５Dqそれぞれから供給される自己相関Ri(1)ないしRi(Dq)の総和である、式（１３）の自己相関和ΣRi(m)を求める。
【０２２９】
相関演算部１９２は、以上のようにして、第i既知シンボル系列についての自己相関和ΣRi(m)を求めると、その自己相関和ΣRi(m)を、平均部１９３に供給する。
【０２３０】
平均部１９３には、相関演算部１９２から自己相関和ΣRi(m)が供給される他、制御装置１５８（図９）から、平均数Pを表す平均数情報が供給される。
【０２３１】
平均部１９３は、相関演算部１９２から、平均数情報が表す平均数Pだけの自己相関、すなわち、P個の第１既知シンボル系列ないし第P既知シンボル系列それぞれについて求められた自己相関和ΣR1(m),ΣR2(m),・・・,ΣRP(m)が供給されると、そのP個の自己相関和ΣR1(m),ΣR2(m),・・・,ΣRP(m)を、式R=(ΣR1(m)+ΣR2(m)+・・・+ΣRP(m))/Pにしたがって平均化することにより、平均自己相関和Rを求めて、周波数誤差演算部１９４に供給する。
【０２３２】
周波数誤差演算部１９４は、１つの既知シンボル系列についての自己相関和ΣR(m)に代えて、平均部１９３から供給される平均自己相関和Rを用い、式（８）にしたがって、周波数誤差の推定値v(q)を求め、すなわち、式（１３）で表される推定値v(q)を求め、累積加算装置１５４（図９）に供給する。
【０２３３】
以上のように、図９の周波数同期部１１０において、シンボル誤差の、複数のラグの自己相関の総和である自己相関和を、複数のラグのうちの最大ラグを小から大に変化させて求め、その自己相関和から、周波数誤差の推定値v(q)を求めて、その推定値v(q)の累積加算を行い、その累積加算により得られる累積加算値を、周波数誤差の最終的な推定値として出力する場合において、複数であるP個の既知シンボル系列それぞれについて求められた自己相関和ΣR1(m)ないしΣRP(m)を平均化することにより、平均自己相関和Rを求め、さらに、その平均自己相関和Rから、周波数誤差の推定値v(q)を求めることにより、周波数誤差の推定の精度を、より向上させることができる。
【０２３４】
次に、図１５は、図３の周波数同期部１１０の第４の構成例を示している。
【０２３５】
なお、図中、図９の場合と対応する部分については、同一の符号を付してあり、以下では、その説明は、適宜省略する。
【０２３６】
すなわち、図１５の周波数同期部１１０は、累積加算装置１５４を有する点で、図９の場合と共通するが、周波数誤差検出装置１５０と制御装置１５８に代えて、それぞれ、周波数誤差検出装置２００と制御装置２０２が設けられている点で、図９の場合と相違する。
【０２３７】
周波数誤差検出装置２００は、１つの検出部２０１から構成される。検出部２０１は、例えば、図９の検出部１５２qと同様に、図１０や、図１３、図１４に示したように構成され、制御装置２０２の制御にしたがって、受信フィルタ１０５から供給される受信データy(k)から、周波数誤差の推定値v(q)を求め、累積加算装置１５４に供給する。
【０２３８】
制御装置２０２は、図９の制御装置１５８と同様に、累積加算装置１５４を制御する他、最大ラグDqと平均数Pを適宜変更し、その最大ラグDqを表す最大ラグ情報と、平均数Pを表す平均数情報を、周波数誤差検出装置２００の検出部２０１に供給することにより、検出部２０１を制御する。
【０２３９】
次に、図１６のフローチャートを参照して、図１５の周波数同期部１１０の処理について説明する。
【０２４０】
なお、ここでは、検出部２０１が、例えば、図１０に示したように構成されることとする。
【０２４１】
ステップＳ６１において、制御装置２０２は、変数qを1に初期化して、処理は、ステップＳ６２に進む。
【０２４２】
ステップＳ６２において、制御装置２０２は、周波数誤差検出装置２００の検出部２０１に設定するパラメータとして、例えば、最大ラグDqを、式D1<D2<・・・<DQを満たすように決定し、その最大ラグDqを表す最大ラグ情報を、検出部２０１に供給する。
【０２４３】
ここで、検出部２０１が、例えば、図１３や図１４に示したように構成される場合には、制御装置２０２は、周波数誤差検出装置２００の検出部２０１に設定するパラメータとして、最大ラグ情報の他、平均数Pを表す平均数情報を、検出部２０１に供給する。
【０２４４】
その後、処理は、ステップＳ６２からステップＳ６３に進み、周波数誤差検出装置２００の検出部２０１は、受信フィルタ１０５から供給される、既知シンボル系列の受信データy(k)を用い、シンボル誤差z(k)を求め、そのシンボル誤差z(k)の、最大ラグmを制御装置２０２からの最大ラグ情報が表す最大ラグDqとする、Dq個の自己相関R(1),R(2),・・・,R(Dq)を求める。さらに、ステップＳ６３では、検出部２０１は、Dq個の自己相関R(1)ないしR(Dq)の総和である自己相関和ΣR(m)を求め、その自己相関和ΣR(m)から、式（８）にしたがい、周波数誤差の推定値v(q)を求める。そして、検出部２０１は、周波数誤差の推定値v(q)を、累積加算装置１５４に供給して、処理は、ステップＳ６３からステップＳ６４に進む。
【０２４５】
ステップＳ６４では、累積加算装置１５４において、演算部１５５が、記憶部１５７の記憶値vmemと、検出部２０１からの周波数誤差の推定値v(q)とを加算し、その結果得られる加算値vmem+v(q)を、セレクタ１５６に供給して、処理は、ステップＳ６５に進む。
【０２４６】
ステップＳ６５では、制御装置２０２が、変数qが、所定の値Qより大きいかどうかを判定する。ここで、所定の値Qとしては、図９の周波数同期部１１０の周波数誤差検出装置１５０を構成する検出部１５２1ないし１５２Qの総数Qに相当する値が、あらかじめ設定される。
【０２４７】
ステップＳ６５において、変数qが、所定の値Qより大きくないと判定された場合、処理は、ステップＳ６６に進み、制御装置２０２は、演算部１５５で求められた加算値vmem+v(q)を選択するように、セレクタ１５６を制御する。
【０２４８】
セレクタ１５６は、制御装置２０２の制御にしたがい、演算部１５５で求められた加算値vmem+v(q)を選択し、周波数誤差の最終的な推定値vとして、記憶部１５７と周波数補正部１０３に供給する。記憶部１５７は、セレクタ１５６からの周波数誤差の最終的な推定値vを、記憶値vmemとして新たに記憶する。
【０２４９】
その後、処理は、ステップＳ６６からステップＳ６７に進み、制御装置２０２は、変数qを1だけインクリメントして、処理は、ステップＳ６３に戻り、次の既知シンボル系列の受信データy(k)を対象として、以下、同様の処理が繰り返される。
【０２５０】
また、ステップＳ６５において、変数qが、所定の値Qより大きいと判定された場合、処理は、ステップＳ６８に進み、制御装置２０２は、記憶部１５７の記憶値vmemを選択するように、セレクタ１５６を制御する。
【０２５１】
セレクタ１５６は、制御装置２０２の制御にしたがい、記憶部１５７の記憶値vmemを選択し、周波数誤差の最終的な推定値vとして、記憶部１５７と周波数補正部１０３に供給する。記憶部１５７は、セレクタ１５６からの周波数誤差の最終的な推定値vを、記憶値vmemとして新たに記憶する。
【０２５２】
その後、処理は、ステップＳ６８に戻り、以下、同様の処理が繰り返される。
【０２５３】
以上のように、検出部２０１に、最大ラグ情報を、適宜供給して、その最大ラグ情報が表す値Dqを最大ラグとする自己相関の自己相関和を求めさせる場合には、周波数誤差検出装置２００を、１つの検出部２０１で構成することができ、装置の規模、及び消費電力を削減することができる。
【０２５４】
なお、制御部１２２（図４）は、CPU(Central Processing Unit)やメモリ等を有するマイクロコンピュータで構成し、制御部１２２が行う処理は、マイクロコンピュータにプログラムを実行させることで行わせることができる。制御部１３３（図７）、制御装置１５８（図９）、及び制御装置２０２（図１５）についても同様である。
【０２５５】
マイクロコンピュータに実行させるプログラムは、そのマイクロコンピュータを構成するメモリにあらかじめインストールしておくことができる。
【０２５６】
また、マイクロコンピュータに実行させるプログラムは、フレキシブルディスク、CD-ROM(Compact Disc Read Only Memory)，MO(Magneto Optical)ディスク，DVD(Digital Versatile Disc)、磁気ディスク、半導体メモリなどのリムーバブル記録媒体に記録し、いわゆるパッケージソフトウエアとして提供することができる。
【０２５７】
さらに、プログラムは、インターネットその他のネットワークを経由してダウンロードすることができる。
【０２５８】
以上、本発明を、BS放送を受信する受信装置に適用した場合について説明したが、本発明は、その他、例えば、搬送波をディジタル変調した変調信号を受信する受信装置等に適用することができる。
【０２５９】
なお、本発明の実施の形態は、上述した実施の形態に限定されるものではなく、本発明の要旨を逸脱しない範囲において種々の変更が可能である。
【図面の簡単な説明】
【０２６０】
【図１】周波数誤差を求める際の不確定性を示す図である。
【図２】従来の受信装置の一例の構成を示すブロック図である。
【図３】本発明を適用した受信装置の一実施の形態の構成例を示すブロック図である。
【図４】周波数同期部１１０の第１の構成例を示すブロック図である。
【図５】受信データの時系列を示す図である。
【図６】周波数同期部１１０の処理を説明するフローチャートである。
【図７】周波数同期部１１０の第２の構成例を示すブロック図である。
【図８】周波数同期部１１０の処理を説明するフローチャートである。
【図９】周波数同期部１１０の第３の構成例を示すブロック図である。
【図１０】検出部１５２qの第１の構成例を示すブロック図である。
【図１１】周波数誤差の推定特性を示す図である。
【図１２】周波数同期部１１０の処理を説明するフローチャートである。
【図１３】検出部１５２qの第２の構成例を示すブロック図である。
【図１４】検出部１５２qの第３の構成例を示すブロック図である。
【図１５】周波数同期部１１０の第４の構成例を示すブロック図である。
【図１６】周波数同期部１１０の処理を説明するフローチャートである。
【符号の説明】
【０２６１】
１１  誤差演算部，  １２  相関演算部，  １３  周波数誤差演算部，  ３１1ないし３１N  相関部，  ３２  演算部，  １０１  アンテナ，  １０２  A/D変換部，  １０３  周波数補正部，  １０４  クロック同期部，  １０５  受信フィルタ，  １０６  等化部，  １０７  位相同期部，  １０８  位相補正部，  １０９  誤り訂正部，  １１０  周波数同期部，  １２１  平均部，  １２２  制御部，  １３１  平均部，  １３２  周波数誤差演算部，  １３３  制御部，  １５０  周波数誤差検出装置，  １５１  スイッチ，  １５２1ないし１５２Q  検出部，  １５３  スイッチ，  １５４  累積加算装置，  １５５  演算部，  １５６  セレクタ，  １５７  記憶部，  １５８  制御装置，  １７１  誤差演算部，  １７２  相関演算部，  １７３  周波数誤差演算部，  １７４1ないし１７４N  相関部，  １７５  演算部，  １８１  誤差演算部，  １８２  相関演算部，  １８３  周波数誤差演算部，  １８４  平均部，  １８５1ないし１８５N  相関部，  １８６  演算部，  １９１  誤差演算部，  １９２  相関演算部，  １９３  平均部，  １９４  周波数誤差演算部，  １９５1ないし１９５N  相関部，  １９６  演算部，  ２００  周波数誤差検出装置，  ２０１  検出部，  ２０２  制御装置
【図１】
2008167230.tif 000016
【図２】
2008167230.tif 000017
【図３】
2008167230.tif 000018
【図４】
2008167230.tif 000019
【図５】
2008167230.tif 000020
【図６】
2008167230.tif 000021
【図７】
2008167230.tif 000022
【図８】
2008167230.tif 000023
【図９】
2008167230.tif 000024
【図１０】
2008167230.tif 000025
【図１１】
2008167230.tif 000026
【図１２】
2008167230.tif 000027
【図１３】
2008167230.tif 000028
【図１４】
2008167230.tif 000029
【図１５】
2008167230.tif 000030
【図１６】
2008167230.tif 000031