"未知のシンボル速度で変調された受信信号から直交検波されたIQ信号を周波数領域の信号に変換し、周波数領域の信号の中でピーク電力が得られる周波数値を特定するとともに、上記周波数値からシンボル速度を粗推定するシンボル粗推定器と、上記IQ信号から各シンボル内のナイキスト位相を検出して、各シンボル内のナイキスト位相から周波数偏差を算出し、上記周波数偏差を上記シンボル粗推定器により粗推定されたシンボル速度に加算するシンボル精推定器とを備えたシンボル推定回路。シンボル粗推定器は、IQ信号に重畳されているノイズを除去する帯域制限フィルタを有し、上記帯域制限フィルタによるノイズ除去後のIQ信号を周波数領域の信号に変換することを特徴とする請求項1記載のシンボル推定回路。シンボル粗推定器は、帯域制限フィルタによるノイズ除去後のIQ信号に対する非線形処理を実施して、非線形処理後のIQ信号を周波数領域の信号に変換し、周波数領域の信号の中でピーク電力が得られる周波数値を特定する第1のシンボル情報検出器と、上記帯域制限フィルタによるノイズ除去後のIQ信号のゼロクロス点の検出処理を実施し、その検出結果を周波数領域の信号に変換し、周波数領域の信号の中でピーク電力が得られる周波数値を特定する第2のシンボル情報検出器と、上記帯域制限フィルタによるノイズ除去後のIQ信号の周波数変動の検出処理を実施し、その検出結果に対する非線形処理を実施して周波数領域の信号に変換し、周波数領域の信号の中でピーク電力が得られる周波数値を特定する第3のシンボル情報検出器と、上記第1、第2及び第3のシンボル情報検出器により特定された周波数値の中から、最も高いピーク電力が得られる周波数値を選択する周波数値選択部と、上記周波数値選択部により選択された周波数値からシンボル速度を算出するシンボル速度算出部とから構成されていることを特徴とする請求項2記載のシンボル推定回路。シンボル精推定器は、帯域制限フィルタによるノイズ除去後のIQ信号から各シンボル内のナイキスト位相を検出するBTR検出回路と、上記BTR検出回路により検出された各シンボル内のナイキスト位相から周波数偏差を算出して、上記周波数偏差をシンボル速度算出部により算出されたシンボル速度に加算し、その加算結果をシンボル速度の精推定結果として出力するシンボル偏差検出回路とから構成されていることを特徴とする請求項3記載のシンボル推定回路。未知のシンボル速度で変調された受信信号から直交検波されたIQ信号を周波数領域の信号に変換し、周波数領域の信号の中でピーク電力が得られる周波数値を特定するとともに、上記周波数値からシンボル速度を粗推定するシンボル粗推定器と、上記IQ信号から各シンボル内のナイキスト位相を検出して、各シンボル内のナイキスト位相から周波数偏差を算出し、上記周波数偏差を上記シンボル粗推定器により粗推定されたシンボル速度に加算するシンボル精推定器と、上記シンボル精推定器により周波数偏差が加算されたシンボル速度を用いて、上記受信信号の復調処理を行う復調器とを備えた復調回路。従来の変調型式識別回路は以上のように構成されているので、受信信号の中心周波数や、シンボル速度などの諸元が事前に分からなければ、受信信号の変調型式を推定することができず、受信信号を復調することができないなどの課題があった。この発明は上記のような課題を解決するためになされたもので、未知のシンボル速度で変調された受信信号を復調するに際して、そのシンボル速度を推定することができるシンボル推定回路を得ることを目的とする。また、この発明は、シンボル推定回路により推定されたシンボル速度を用いて、変調型式が未知の受信信号を復調することができる復調回路を得ることを目的とする。この発明に係るシンボル推定回路は、未知のシンボル速度で変調された受信信号から直交検波されたIQ信号を周波数領域の信号に変換し、周波数領域の信号の中でピーク電力が得られる周波数値を特定するとともに、その周波数値からシンボル速度を粗推定するシンボル粗推定器と、そのIQ信号から各シンボル内のナイキスト位相を検出して、各シンボル内のナイキスト位相から周波数偏差を算出し、その周波数偏差をシンボル粗推定器により粗推定されたシンボル速度に加算するシンボル精推定器とを設けるようにしたものである。この発明によれば、未知のシンボル速度で変調された受信信号から直交検波されたIQ信号を周波数領域の信号に変換し、周波数領域の信号の中でピーク電力が得られる周波数値を特定するとともに、その周波数値からシンボル速度を粗推定するシンボル粗推定器と、そのIQ信号から各シンボル内のナイキスト位相を検出して、各シンボル内のナイキスト位相から周波数偏差を算出し、その周波数偏差をシンボル粗推定器により粗推定されたシンボル速度に加算するシンボル精推定器とを設けるように構成したので、未知のシンボル速度で変調された受信信号を復調するに際して、そのシンボル速度を推定することができる効果がある。実施の形態1.図1はこの発明の実施の形態1によるシンボル推定回路を実装している復調回路を示す構成図である。図1において、シンボル推定回路1は未知のシンボル速度で変調された受信信号を復調するに際して、そのシンボル速度を推定する装置である。シンボル粗推定器2は未知のシンボル速度で変調された受信信号から直交検波されたIQ信号を入力し、そのIQ信号に重畳されているノイズガウス雑音を除去する帯域制限フィルタ2aを実装している。シンボル粗推定器2は帯域制限フィルタ2aによるノイズ除去後のIQ信号を周波数領域の信号に変換し、周波数領域の信号の中でピーク電力が得られる周波数値を特定するとともに、その周波数値からシンボル速度を粗推定する処理を実施する。帯域制限フィルタ2aは例えば99法や3dB法を実施することで占有帯域幅を推定する図示せぬ推定回路等から出力された占有帯域幅推定値Fwまたは、既知の情報として、スペアナ等で測定された占有帯域幅値を入力し、その占有帯域幅推定値Fwのα倍例えば、約1.3倍=1.3×Fwのカットオフ周波数でIQ信号を通過させるものである。これにより、シンボル帯域の推定に必要なIQ信号がフィルタリングで肩落ちしない程度にガウス雑音の除去が行われる。シンボル精推定器3は帯域制限フィルタ2aによるノイズ除去後のIQ信号から各シンボル内のナイキスト位相を検出して、各シンボル内のナイキスト位相から周波数偏差を算出し、その周波数偏差をシンボル粗推定器2により粗推定されたシンボル速度に加算する処理を実施する。復調器4はシンボル推定回路1のシンボル精推定器3により周波数偏差が加算されたシンボル速度を用いて、上記受信信号の復調処理を行う。図2はこの発明の実施の形態1によるシンボル推定回路1のシンボル粗推定器2を示す構成図である。図2において、シンボル情報検出器10はIQ信号の非線形処理とFFT高速フーリエ変換処理を実施することで、周波数領域の信号の中でピーク電力MAXPAVE1ωが得られる周波数値ωPT1を特定する処理を実施する。なお、シンボル情報検出器10は第1のシンボル情報検出器を構成している。シンボル情報検出器20はIQ信号のゼロクロス点の検出処理とFFT処理を実施することで、周波数領域の信号の中でピーク電力MAXPAVE2ωが得られる周波数値ωPT2を特定する処理を実施する。なお、シンボル情報検出器20は第2のシンボル情報検出器を構成している。シンボル情報検出器30はIQ信号の周波数変動の検出処理、非線形処理及びFFT処理を実施することで、周波数領域の信号の中でピーク電力MAXPAVE3ωが得られる周波数値ωPT3を特定する処理を実施する。なお、シンボル情報検出器30は第3のシンボル情報検出器を構成している。シンボル情報検出器10のLPF部11は直交検波されたIQ信号SRnに対してフィルタリング処理を実施する。このフィルタリング処理は、例えば、ルートナイキストフィルタ例えば、カットオフ周波数fc=3dB:占有帯域幅推定値Fw×1.3を用いるものであり、帯域制限フィルタ2aのフィルタリング処理に相当する。非線形部12はLPF部11によるフィルタリング処理後のIQ信号S’Rnに対する非線形処理2乗処理を実施し、非線形処理後のIQ信号fnを出力する。FFT部13は非線形部12による非線形処理後のIQ信号fnに対してハニング窓関数係数Wnを乗算したのち、Nポイント例えば、1024ポイントの複素FFT処理を実施する。FFT時間領域平均部14はFFT部13による複数個のFFT結果Fωに対して、複数回の時間加算処理を行うことで、0~1023の周波数値ωの範囲で、周波数領域信号の電力PAVE1ωを算出する。ピーク電力検出部15はFFT時間領域平均部14により算出された周波数領域信号の電力PAVE1ωの中で、最も高いピーク電力MAXPAVE1ωを検出して、そのピーク電力MAXPAVE1ωが得られる周波数値ωPAVE1をシンボル周波数として出力する処理を実施する。ラグランジェ補間部16はピーク電力検出部15から出力されたシンボル周波数ωPAVE1に対する補間処理例えば、ラグランジェ2次補間処理を実施し、補間処理後のシンボル周波数ωPT1及びピーク電力MAXPAVE1ωを周波数値選択部40に出力する。シンボル情報検出器20の周波数偏差検出部21は直交検波されたIQ信号SRnから、例えば逓倍法などを実施することで、IQ信号SRnの中心周波数Δfcの推定を行う。ここでは、周波数偏差検出部21が中心周波数Δfcを推定しているが、既知情報として中心周波数Δfcが事前に与えられるようにしてもよい。周波数回転部22は周波数偏差検出部21により推定された中心周波数Δfc分だけ、IQ信号SRnの周波数を回転させる処理を実施する。LPF部23は周波数回転部22による周波数回転処理後のIQ信号SROTnに対してフィルタリング処理を実施する。このフィルタリング処理は、例えば、ルートナイキストフィルタ例えば、カットオフ周波数fc=3dB:占有帯域幅推定値Fw×1.3を用いるものであり、帯域制限フィルタ2aのフィルタリング処理に相当する。ゼロクロス検出部24はLPF部23によるフィルタリング処理後のIQ信号S’ROTnと1シンボル前のIQ信号S’ROTn-1との間のゼロクロス点の検出処理を実施して、パルス波のIQ信号fnを生成する。FFT部25はゼロクロス検出部24により生成されたパルス波のIQ信号fnに対してハニング窓関数係数Wnを乗算したのち、Nポイント例えば、1024ポイントの複素FFT処理を実施する。FFT時間領域平均部26はFFT部25による複数個のFFT結果Fωに対して、複数回の時間加算処理を行うことで、0~1023の周波数値ωの範囲で、周波数領域信号の電力PAVE2ωを算出する。ピーク電力検出部27はFFT時間領域平均部26により算出された周波数領域信号の電力PAVE2ωの中で、最も高いピーク電力MAXPAVE2ωを検出して、そのピーク電力MAXPAVE2ωが得られる周波数値ωPAVE2をシンボル周波数として出力する処理を実施する。ラグランジェ補間部28はピーク電力検出部27から出力されたシンボル周波数ωPAVE2に対する補間処理例えば、ラグランジェ2次補間処理を実施し、補間処理後のシンボル周波数ωPT2及びピーク電力MAXPAVE2ωを周波数値選択部40に出力する。シンボル情報検出器30のLPF部31は直交検波されたIQ信号SRnに対してフィルタリング処理を実施する。このフィルタリング処理は、例えば、ルートナイキストフィルタ例えば、カットオフ周波数fc=3dB:占有帯域幅推定値Fw×1.3を用いるものであり、帯域制限フィルタ2aのフィルタリング処理に相当する。周波数ディスクリミネータ部32はLPF部31によるフィルタリング処理後のIQ信号S’Rnの位相情報θnを微分演算することで、IQ信号S’Rnの瞬時周波数変動を検出する。非線形部33は周波数ディスクリミネータ部32の検出結果に対する非線形処理2乗処理を実施し、非線形処理後の信号fnを出力する。FFT部34は非線形部33による非線形処理後の信号fnに対してハニング窓関数係数Wnを乗算したのち、Nポイント例えば、1024ポイントの複素FFT処理を実施する。FFT時間領域平均部35はFFT部34による複数個のFFT結果Fωに対して、複数回の時間加算処理を行うことで、0~1023の周波数値ωの範囲で、周波数領域信号の電力PAVE3ωを算出する。ピーク電力検出部36はFFT時間領域平均部35により算出された周波数領域信号の電力PAVE3ωの中で、最も高いピーク電力MAXPAVE3ωを検出して、そのピーク電力MAXPAVE3ωが得られる周波数値ωPAVE3をシンボル周波数として出力する処理を実施する。ラグランジェ補間部37はピーク電力検出部36から出力されたシンボル周波数ωPAVE3に対する補間処理例えば、ラグランジェ2次補間処理を実施し、補間処理後のシンボル周波数ωPT3及びピーク電力MAXPAVE3ωを周波数値選択部40に出力する。周波数値選択部40はシンボル情報検出器10,20,30から出力された補間処理後のシンボル周波数ωPT1,ωPT2,ωPT3の中から、最も高いピーク電力が得られるシンボル周波数周波数値ωPTを選択する処理を実施する。シンボル速度算出部50は周波数値選択部40により選択されたシンボル周波数ωPTからシンボル速度fRを算出する処理を実施する。図3はこの発明の実施の形態1によるシンボル推定回路1のシンボル精推定器3を示す構成図である。図3において、非線形処理によるBTRBitTimingRecovery検出回路61は周波数値選択部40により選択されたシンボル周波数ωPTが、シンボル情報検出器10から出力された補間処理後のシンボル周波数ωPT1又はシンボル情報検出器30から出力された補間処理後のシンボル周波数ωPT3である場合、帯域制限フィルタ2aによるノイズ除去後のIQ信号S’Rnに対して非線形処理を実施することで、各シンボル内のナイキスト位相θNnを検出する処理を実施する。ゼロクロスによるBTR検出回路62は周波数値選択部40により選択されたシンボル周波数ωPTが、シンボル情報検出器20から出力された補間処理後のシンボル周波数ωPT2である場合、帯域制限フィルタ2aによるノイズ除去後のIQ信号S’Rnのゼロクロス点の検出処理を実施することで、各シンボル内のナイキスト位相θNnを検出する処理を実施する。シンボル偏差検出回路63はBTR検出回路61又はBTR検出回路62により検出された各シンボル内のナイキスト位相θNnから周波数偏差ΔfRを算出して、その周波数偏差ΔfRをシンボル粗推定器2のシンボル速度算出部50により算出されたシンボル速度fRに加算し、その加算結果fRΔfRをシンボル速度の精推定結果として復調器4に出力する処理を実施する。次に動作について説明する。シンボル粗推定器2は、未知のシンボル速度で変調された受信信号から直交検波されたIQ信号SRnを入力すると、帯域制限フィルタ2aを用いて、そのIQ信号SRnに重畳されているノイズガウス雑音を除去する。シンボル粗推定器2は、帯域制限フィルタ2aによるノイズ除去後のIQ信号を周波数領域の信号に変換し、周波数領域の信号の中でピーク電力が得られるシンボル周波数ωPT周波数値を特定するとともに、そのシンボル周波数ωPTからシンボル速度fRを粗推定する。シンボル粗推定器2では、FFTポイント精度例えば、1024ポイントのFFT処理を実施する場合、1000ppm程度の偏差精度による粗いシンボル精度の推定処理を行う。以下、シンボル粗推定器2の処理内容を具体的に説明する。シンボル情報検出器10は、未知のシンボル速度で変調された受信信号から直交検波されたIQ信号SRnを入力すると、そのIQ信号SRnの非線形処理とFFT処理を実施することで、周波数領域の信号の中でピーク電力MAXPAVE1ωが得られる周波数値ωPT1を特定する。即ち、シンボル情報検出器10のLPF部11は、例えば、ルートナイキストフィルタ例えば、カットオフ周波数fc=3dB:占有帯域幅推定値Fw×1.3を用いて、直交検波されたIQ信号SRnに対するフィルタリング処理を実施し、フィルタリング処理後のIQ信号S’Rnを非線形部12に出力する。S’Rn=S’RInjS’RQn1ただし、jは複素数である。非線形部12は、LPF部11がIQ信号SRnに対するフィルタリング処理を実施すると、下記の式2に示すように、フィルタリング処理後のIQ信号S’Rnに対する非線形処理2乗処理を実施して、非線形処理後のIQ信号fnをFFT部13に出力する。fn=S’Rn×S’Rn=S’RIn×S’RInS’RQn×S’RQn2ただし、は共役を示す記号である。FFT部13は、非線形部12から非線形処理後のIQ信号fnを受けると、下記の式3に示すように、その非線形処理後のIQ信号fnに対してハニング窓関数係数Wnを乗算したのち、Nポイント例えば、1024ポイントの複素FFT処理を実施することで、複数個のFFT結果FωをFFT時間領域平均部14に出力する。FFT時間領域平均部14は、FFT部13による複数個=m個のFFT結果Fωに対して、複数回の時間加算処理を行うことで、0~1023の周波数値ωの範囲で、周波数領域信号の電力PAVE1ωを算出する。即ち、FFT時間領域平均部14は、Fω導出時のIQ信号fnのn値を、例えば256ずつスライドさせながら、下記の式4に示すように、4m-1回の加算処理を行う。ピーク電力検出部15は、FFT時間領域平均部14により算出された周波数領域信号の電力PAVE1ωの中で、最も高いピーク電力MAXPAVE1ωを検出する。ピーク電力検出部15は、ピーク電力MAXPAVE1ωを検出すると、そのピーク電力MAXPAVE1ωが得られる周波数値ωPAVE1をシンボル周波数としてラグランジェ補間部16に出力する。ただし、ピーク電力を検出する際、周波数値の範囲ωは、実際にピーク電力が検出されることが期待される範囲に限定して探索するようにしてもよい。ラグランジェ補間部16は、ピーク電力検出部15からシンボル周波数ωPAVE1を受けると、そのシンボル周波数ωPAVE1に対する補間処理を実施して、補間処理後のシンボル周波数ωPT1及びピーク電力MAXPAVE1ωを周波数値選択部40に出力する。例えば、下記の式5に示すようなラグランジェ2次補間処理を実施する。図4はシンボル情報検出器10の検出結果である補間処理後のシンボル周波数ωPT1と、ピーク電力MAXPAVE1ωとを示す説明図である。図4の例では、FFTポイント番号が256の周波数がシンボル周波数ωPT1であり、約-12dBの電力が得られている。シンボル情報検出器20は、未知のシンボル速度で変調された受信信号から直交検波されたIQ信号SRnを入力すると、そのIQ信号SRnのゼロクロス点の検出処理とFFT処理を実施することで、周波数領域の信号の中でピーク電力MAXPAVE2ωが得られる周波数値ωPT2を特定する。即ち、シンボル情報検出器20の周波数偏差検出部21は、直交検波されたIQ信号SRnから、例えば逓倍法などを実施することで、そのIQ信号SRnの中心周波数Δfcの推定を行う。ここでは、周波数偏差検出部21が中心周波数Δfcを推定しているが、既知情報として中心周波数Δfcが事前に与えられるようにしてもよい。周波数回転部22は、周波数偏差検出部21が中心周波数Δfcを推定すると、下記の式6に示すように、その中心周波数Δfc分だけ、IQ信号SRnの周波数を回転させる。SROTn=SRn×exp-2πΔfc/Fs×n6ただし、FsはデジタルIQ信号をサンプリングしたA/DAnalogtoDigitalConvertorサンプリング周波数、nは時間サンプル番号n=0,,1023である。LPF部23は、例えば、ルートナイキストフィルタ例えば、カットオフ周波数fc=3dB:占有帯域幅推定値Fw×1.3を用いて、周波数回転部22による周波数回転処理後のIQ信号SROTnに対するフィルタリング処理を実施し、フィルタリング処理後のIQ信号S’ROTnをゼロクロス検出部24に出力する。S’ROTn=S’ROTInjS’ROTQn7ただし、jは複素数である。ゼロクロス検出部24は、LPF部23からフィルタリング処理後のIQ信号S’ROTnを受けると、そのIQ信号S’ROTnと1シンボル前のIQ信号S’ROTn-1との間のゼロクロス点の検出処理を実施して、パルス波のIQ信号fnを生成する。fn=doutInj×doutQn8ただし、jは複素数である。即ち、ゼロクロス検出部24は、Ich信号及びQch信号について下記の計算を行うことにより、パルス波のIQ信号fnを生成する。ただし、下記の計算において、InZeroDet1はゼロクロス有、InZeroDet0はゼロクロス無を示している。FFT部25は、ゼロクロス検出部24がパルス波のIQ信号fnを生成すると、下記の式9に示すように、そのパルス波のIQ信号fnに対してハニング窓関数係数Wnを乗算したのち、Nポイント例えば、1024ポイントの複素FFT処理を実施することで、複数個のFFT結果FωをFFT時間領域平均部26に出力する。FFT時間領域平均部26は、FFT部25による複数個=m個のFFT結果Fωに対して、複数回の時間加算処理を行うことで、0~1023の周波数値ωの範囲で、周波数領域信号の電力PAVE2ωを算出する。即ち、FFT時間領域平均部26は、Fω導出時のIQ信号fnのn値を、例えば256ずつスライドさせながら、下記の式10に示すように、4m-1回の加算処理を行う。ピーク電力検出部27は、FFT時間領域平均部26により算出された周波数領域信号の電力PAVE2ωの中で、最も高いピーク電力MAXPAVE2ωを検出する。ピーク電力検出部27は、ピーク電力MAXPAVE2ωを検出すると、そのピーク電力MAXPAVE2ωが得られる周波数値ωPAVE2をシンボル周波数としてラグランジェ補間部28に出力する。ただし、ピーク電力を検出する際、周波数値の範囲ωは、実際にピーク電力が検出されることが期待される範囲に限定して探索するようにしてもよい。ラグランジェ補間部28は、ピーク電力検出部27からシンボル周波数ωPAVE2を受けると、そのシンボル周波数ωPAVE2に対する補間処理を実施して、補間処理後のシンボル周波数ωPT2及びピーク電力MAXPAVE2ωを周波数値選択部40に出力する。例えば、下記の式11に示すようなラグランジェ2次補間処理を実施する。シンボル情報検出器30は、未知のシンボル速度で変調された受信信号から直交検波されたIQ信号SRnを入力すると、そのIQ信号SRnの周波数変動の検出処理、非線形処理及びFFT処理を実施することで、周波数領域の信号の中でピーク電力MAXPAVE3ωが得られる周波数値ωPT3を特定する。即ち、シンボル情報検出器30のLPF部31は、例えば、ルートナイキストフィルタ例えば、カットオフ周波数fc=3dB:占有帯域幅推定値Fw×1.3を用いて、直交検波されたIQ信号SRnに対するフィルタリング処理を実施し、フィルタリング処理後のIQ信号S’Rnを周波数ディスクリミネータ部32に出力する。S’Rn=S’RInjS’RQn12ただし、jは複素数である。周波数ディスクリミネータ部32は、LPF部31からフィルタリング処理後のIQ信号S’Rnを受けると、下記の式13に示すように、そのIQ信号S’Rnの位相情報θnを微分演算することで、IQ信号S’Rnの瞬時周波数変動を検出する。非線形部33は、周波数ディスクリミネータ部32の検出結果fnに対する非線形処理2乗処理を実施して、非線形処理後の信号fnをFFT部34に出力する。FFT部34は、非線形部33から非線形処理後の信号fnを受けると、下記の式14に示すように、その非線形処理後の信号fnに対してハニング窓関数係数Wnを乗算したのち、Nポイント例えば、1024ポイントの複素FFT処理を実施することで、複数個のFFT結果FωをFFT時間領域平均部35に出力する。FFT時間領域平均部35は、FFT部34による複数個=m個のFFT結果Fωに対して、複数回の時間加算処理を行うことで、0~1023の周波数値ωの範囲で、周波数領域信号の電力PAVE3ωを算出する。即ち、FFT時間領域平均部35は、Fω導出時の信号fnのn値を、例えば256ずつスライドさせながら、下記の式15に示すように、4m-1回の加算処理を行う。ピーク電力検出部36は、FFT時間領域平均部35により算出された周波数領域信号の電力PAVE3ωの中で、最も高いピーク電力MAXPAVE3ωを検出する。ピーク電力検出部36は、ピーク電力MAXPAVE3ωを検出すると、そのピーク電力MAXPAVE3ωが得られる周波数値ωPAVE3をシンボル周波数としてラグランジェ補間部37に出力する。ただし、ピーク電力を検出する際、周波数値の範囲ωは、実際にピーク電力が検出されることが期待される範囲に限定して探索するようにしてもよい。ラグランジェ補間部37は、ピーク電力検出部36からシンボル周波数ωPAVE3を受けると、そのシンボル周波数ωPAVE3に対する補間処理を実施して、補間処理後のシンボル周波数ωPT3及びピーク電力MAXPAVE3ωを周波数値選択部40に出力する。例えば、下記の式16に示すようなラグランジェ2次補間処理を実施する。周波数値選択部40は、シンボル情報検出器10,20,30から補間処理後のシンボル周波数ωPT1,ωPT2,ωPT3と、ピーク電力MAXPAVE1ω,MAXPAVE2ω,MAXPAVE3ωを受けると、ピーク電力MAXPAVE1ω,MAXPAVE2ω,MAXPAVE3ωを比較して、補間処理後のシンボル周波数ωPT1,ωPT2,ωPT3の中から、最も高いピーク電力が得られるシンボル周波数ωPTを選択する。例えば、MAXPAVE2ω>MAXPAVE3ω>MAXPAVE1ωであれば、補間処理後のシンボル周波数ωPT2を選択し、MAXPAVE3ω>MAXPAVE1ω>MAXPAVE2ωであれば、補間処理後のシンボル周波数ωPT3を選択する。シンボル速度算出部50は、周波数値選択部40がシンボル周波数ωPTを選択すると、下記の式17に示すように、そのシンボル周波数ωPTからシンボル速度fRを算出し、そのシンボル速度fRを粗推定結果としてシンボル精推定器3に出力する。また、周波数値選択部40により選択されたシンボル周波数ωPTをシンボル精推定器3に出力する。fR=ωPT/N×Fs17ただし、NはFFTのポイント数、FsはA/Dのサンプリング周波数である。シンボル精推定器3は、シンボル粗推定器2からシンボル速度fR及びシンボル周波数ωPTを受けると、帯域制限フィルタ2aによるノイズ除去後のIQ信号S’Rnから各シンボル内のナイキスト位相θNnを検出して、各シンボル内のナイキスト位相θNnから周波数偏差ΔfRを算出し、その周波数偏差ΔfRをシンボル粗推定器2により粗推定されたシンボル速度fRに加算する。即ち、シンボル精推定器3の非線形処理によるBTR検出回路61は、シンボル粗推定器2から出力されたシンボル周波数ωPTが、シンボル情報検出器10から出力された補間処理後のシンボル周波数ωPT1である場合、あるいは、シンボル情報検出器30から出力された補間処理後のシンボル周波数ωPT3である場合、帯域制限フィルタ2aによるノイズ除去後のIQ信号S’Rnに対して非線形処理を実施することで、各シンボル内のナイキスト位相θNnを検出する。非線形処理によるBTR検出回路61がナイキスト位相θNnを検出する処理自体は公知の技術であるため詳細な説明を省略する。ゼロクロスによるBTR検出回路62は、シンボル粗推定器2から出力されたシンボル周波数ωPTが、シンボル情報検出器20から出力された補間処理後のシンボル周波数ωPT2である場合、帯域制限フィルタ2aによるノイズ除去後のIQ信号S’Rnのゼロクロス点の検出処理を実施することで、各シンボル内のナイキスト位相θNnを検出する。ゼロクロスによるBTR検出回路62がナイキスト位相θNnを検出する処理自体は公知の技術であるため詳細な説明を省略する。ここで、図5はBTR検出回路61,62により検出される各シンボル内のナイキスト位相θNn及び周波数偏差ΔfRであるシンボル速度偏差の一例を示す説明図である。シンボル偏差検出回路63は、BTR検出回路61又はBTR検出回路62が各シンボル内のナイキスト位相θNnを検出すると、下記の式18に示すように、そのナイキスト位相θNnから周波数偏差ΔfRを算出する。ΔfR=θN-θN/360×Fs18ただし、FsはA/Dのサンプリング周波数である。図5の例では、nAシンボル目のナイキスト位相θNnAが140、nBシンボル目のナイキスト位相θNnBが230であるため、周波数偏差ΔfRは、下記のようになる。ΔfR=θN-θN/360×Fs=230140360×FsHzシンボル偏差検出回路63は、ナイキスト位相θNnから周波数偏差ΔfRを算出すると、その周波数偏差ΔfRをシンボル粗推定器2のシンボル速度算出部50により算出されたシンボル速度fRに加算し、その加算結果fRΔfRをシンボル速度の精推定結果として復調器4に出力する。これにより、シンボル粗推定器2の粗推定結果では、FFT1024ptの精度約1000ppmであるが、シンボル精推定器3の精推定結果では、約50ppm以下の偏差に改善される。復調器4は、シンボル推定回路1のシンボル精推定器3の精推定結果であるシンボル速度fRΔfRを用いて、上記受信信号の復調処理を行う。復調器4の復調処理自体は公知の技術であるため詳細な説明を省略する。以上で明らかなように、この実施の形態1によれば、未知のシンボル速度で変調された受信信号から直交検波されたIQ信号SRnに重畳されているノイズを除去する帯域制限フィルタ2aを実装し、帯域制限フィルタ2aによるノイズ除去後のIQ信号S’Rnを周波数領域の信号に変換し、周波数領域の信号の中でピーク電力が得られるシンボル周波数ωPTを特定するとともに、そのシンボル周波数ωPTからシンボル速度fRを粗推定するシンボル粗推定器2と、帯域制限フィルタ2aによるノイズ除去後のIQ信号S’Rnから各シンボル内のナイキスト位相θNnを検出して、各シンボル内のナイキスト位相θNnから周波数偏差ΔfRを算出し、その周波数偏差ΔfRをシンボル粗推定器2により粗推定されたシンボル速度fRに加算するシンボル精推定器3とを設けるように構成したので、未知のシンボル速度で変調された受信信号を復調するに際して、高精度にシンボル速度を推定することができる効果を奏する。なお、本願発明はその発明の範囲内において、実施の形態の任意の構成要素の変形、もしくは実施の形態の任意の構成要素の省略が可能である。"