JP2003152674A - 周波数同期方法及びこれを用いたｏｆｄｍ受信装置 - Google Patents
周波数同期方法及びこれを用いたｏｆｄｍ受信装置Info
- Publication number
- JP2003152674A JP2003152674A JP2001345466A JP2001345466A JP2003152674A JP 2003152674 A JP2003152674 A JP 2003152674A JP 2001345466 A JP2001345466 A JP 2001345466A JP 2001345466 A JP2001345466 A JP 2001345466A JP 2003152674 A JP2003152674 A JP 2003152674A
- Authority
- JP
- Japan
- Prior art keywords
- frequency offset
- frequency
- signal
- unit
- correlation value
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Landscapes
- Synchronisation In Digital Transmission Systems (AREA)
Abstract
(57)【要約】
【課題】 広い周波数同期引き込み範囲を実現し、周波
数オフセットの推定誤差を小さくできる周波数同期方法
とそれを用いた受信装置を提供する。 【解決手段】 送信側では固定信号系列を複数繰り返し
たプリアンブルを送信し、受信側では周波数同期確立の
ため前記送信側で送信されたプリアンブルの受信信号に
基づいて周波数オフセットを補正するＯＦＤＭ伝送シス
テムに適用され、受信信号中の同期用固定信号系列か
ら、同期引き込み周波数範囲および周波数オフセット推
定精度の異なる誤差信号を複数生成し、これらの大小比
較および正負の判定により周波数オフセット推定に用い
る誤差信号を選択する。
数オフセットの推定誤差を小さくできる周波数同期方法
とそれを用いた受信装置を提供する。 【解決手段】 送信側では固定信号系列を複数繰り返し
たプリアンブルを送信し、受信側では周波数同期確立の
ため前記送信側で送信されたプリアンブルの受信信号に
基づいて周波数オフセットを補正するＯＦＤＭ伝送シス
テムに適用され、受信信号中の同期用固定信号系列か
ら、同期引き込み周波数範囲および周波数オフセット推
定精度の異なる誤差信号を複数生成し、これらの大小比
較および正負の判定により周波数オフセット推定に用い
る誤差信号を選択する。
Description
【０００１】
【発明の属する技術分野】本発明は、周波数同期方法及
びこれを用いたＯＦＤＭ受信装置に関し、特に直交周波
数分割多重（ＯＦＤＭ）方式における周波数同期方法及
び受信装置に利用して有効な技術に関するものである。
びこれを用いたＯＦＤＭ受信装置に関し、特に直交周波
数分割多重（ＯＦＤＭ）方式における周波数同期方法及
び受信装置に利用して有効な技術に関するものである。
【０００２】
【従来の技術】ＯＦＤＭ伝送方式は、互いに直交する多
数の搬送波を用いるディジタル変調方式であり、他のサ
ービスに妨害を与えにくく、妨害を受けにくい、周波数
利用効率が比較的良いなどの特徴があるが、周波数オフ
セットによる信号歪が大きく、高精度の周波数同期が必
要となる。従来のＯＦＤＭ用周波数同期手法は、例えば
文献"A Fast Synchronization Scheme of OFDM Signals
for High-Rate Wireless LAN" T.Onizawa et.al, IEIC
E TRANS. COMMUN. VOL.E82-B. NO.2等に示されている。
数の搬送波を用いるディジタル変調方式であり、他のサ
ービスに妨害を与えにくく、妨害を受けにくい、周波数
利用効率が比較的良いなどの特徴があるが、周波数オフ
セットによる信号歪が大きく、高精度の周波数同期が必
要となる。従来のＯＦＤＭ用周波数同期手法は、例えば
文献"A Fast Synchronization Scheme of OFDM Signals
for High-Rate Wireless LAN" T.Onizawa et.al, IEIC
E TRANS. COMMUN. VOL.E82-B. NO.2等に示されている。
【０００３】以下、従来の周波数同期手法を説明する。
この周波数同期手法を用いて周波数同期を行うための送
信パケットフォーマットを図５に示す。送信パケット１
０はプリアンブル１０１とデータ１０２からなり、プリ
アンブル１０１を用いて周波数オフセットを推定し、周
波数補正することにより、周波数同期を行う。プリアン
ブル１０１の構成を図２に示す。プリアンブル１０１
は、１６サンプル長の同一の固定信号系列２０，２１，
２２，２３，２４，２５，２６，２７，２８および２９
からなる。
この周波数同期手法を用いて周波数同期を行うための送
信パケットフォーマットを図５に示す。送信パケット１
０はプリアンブル１０１とデータ１０２からなり、プリ
アンブル１０１を用いて周波数オフセットを推定し、周
波数補正することにより、周波数同期を行う。プリアン
ブル１０１の構成を図２に示す。プリアンブル１０１
は、１６サンプル長の同一の固定信号系列２０，２１，
２２，２３，２４，２５，２６，２７，２８および２９
からなる。
【０００４】次に、従来の周波数同期方法を、図３を用
いて説明する。図３において、ｓ（ｎ）は受信されたＯ
ＦＤＭ信号であり、Ｒｅ（ｓ（ｎ））およびＩｍ（ｓ
（ｎ））はそれぞれｓ（ｎ）の実数部および虚数部を表
す。また、ｏｕｔ（ｎ）は周波数補正後の信号を表す。
ここで、ｎはサンプリング周期をＴとする受信側での離
散時刻を表す。誤差信号生成部３０において、プリアン
ブル１０１の受信信号ｓ（ｎ）から、周波数オフセット
量に応じた誤差信号Ｅ（ｎ）を計算する。次に、誤差信
号変換部３１にて、誤差信号Ｅ（ｎ）から周波数オフセ
ットの推定値Δｆ ESTを求める。正弦波発生部３２で
は、ΔｆESTの周波数で振動する複素数の正弦波ｃ
（ｎ）を出力する。このｃ（ｎ）と、入力信号ｓ（ｎ）
を、複素乗算部３３にて乗算することで、周波数オフセ
ット補正されたｏｕｔ（ｎ）を得る。
いて説明する。図３において、ｓ（ｎ）は受信されたＯ
ＦＤＭ信号であり、Ｒｅ（ｓ（ｎ））およびＩｍ（ｓ
（ｎ））はそれぞれｓ（ｎ）の実数部および虚数部を表
す。また、ｏｕｔ（ｎ）は周波数補正後の信号を表す。
ここで、ｎはサンプリング周期をＴとする受信側での離
散時刻を表す。誤差信号生成部３０において、プリアン
ブル１０１の受信信号ｓ（ｎ）から、周波数オフセット
量に応じた誤差信号Ｅ（ｎ）を計算する。次に、誤差信
号変換部３１にて、誤差信号Ｅ（ｎ）から周波数オフセ
ットの推定値Δｆ ESTを求める。正弦波発生部３２で
は、ΔｆESTの周波数で振動する複素数の正弦波ｃ
（ｎ）を出力する。このｃ（ｎ）と、入力信号ｓ（ｎ）
を、複素乗算部３３にて乗算することで、周波数オフセ
ット補正されたｏｕｔ（ｎ）を得る。
【０００５】誤差信号生成部３０の構成を図４に示す。
誤差信号生成部３０では、遅延部３０１で、入力信号を
６４サンプル遅延させ、遅延信号ｓ（ｎ−６４）を得
る。このｓ（ｎ−６４）に対し、複素共役部３０２にて
複素共役ｓ（ｎ−６４）*を求め、複素乗算部３０３に
て、入力信号ｓ（ｎ）と複素共役部３０２から出力され
たｓ（ｎ−６４）*の積Ｓ（ｎ）を求める。受信したプ
リアンブルは６４サンプル周期の繰り返し信号であるの
で、周波数オフセットがない状態では、ｓ（ｎ）とｓ
（ｎ−６４）はノイズ成分を除いて同一となり、Ｓ
（ｎ）は実数となる。しかし、周波数オフセットΔｆが
ある状態では、ｓ（ｎ−６４）は式（１）のようにな
り、ｓ（ｎ）とｓ（ｎ−６４）の間で位相差が生じる。
誤差信号生成部３０では、遅延部３０１で、入力信号を
６４サンプル遅延させ、遅延信号ｓ（ｎ−６４）を得
る。このｓ（ｎ−６４）に対し、複素共役部３０２にて
複素共役ｓ（ｎ−６４）*を求め、複素乗算部３０３に
て、入力信号ｓ（ｎ）と複素共役部３０２から出力され
たｓ（ｎ−６４）*の積Ｓ（ｎ）を求める。受信したプ
リアンブルは６４サンプル周期の繰り返し信号であるの
で、周波数オフセットがない状態では、ｓ（ｎ）とｓ
（ｎ−６４）はノイズ成分を除いて同一となり、Ｓ
（ｎ）は実数となる。しかし、周波数オフセットΔｆが
ある状態では、ｓ（ｎ−６４）は式（１）のようにな
り、ｓ（ｎ）とｓ（ｎ−６４）の間で位相差が生じる。
【０００６】
ｓ（ｎ−６４）＝ｓ（ｎ）ｅｘｐ（ｊ２π・６４ＴΔｆ） ・・・・・（１）
このため、Ｓ（ｎ）は式（２）のように、位相成分を持つ。
Ｓ（ｎ）＝| Ｓ（ｎ）|2ｅｘｐ（ｊ２π・６４ＴΔｆ） ・・・・・・（２）
【０００７】積算部３０４は、上記方法にて求められた
Ｓ（ｎ）を、６４サンプルにわたって積算し、相関値Σ
Ｓ（ｎ）を求める。また除算部３０５は、積算部３０４
の出力ΣＳ（ｎ）について、式（３）により誤差信号Ｅ
（ｎ）を得る。 Ｅ（ｎ９＝Ｉｍ（ΣＳ（ｎ））／Ｒｅ（ΣＳ（ｎ）） ・・・・・・・（３） 除算部３０５から出力される誤差信号Ｅ（ｎ）は、式
（４）、式（５）のように対応付ける事ができる。 Ｅ（ｎ）＝ｔａｎ（ΔθEST) ・・・・・・・・（４） ΔθEST ＝２π・６４ＴΔｆEST ・・・・・・・・（５）
Ｓ（ｎ）を、６４サンプルにわたって積算し、相関値Σ
Ｓ（ｎ）を求める。また除算部３０５は、積算部３０４
の出力ΣＳ（ｎ）について、式（３）により誤差信号Ｅ
（ｎ）を得る。 Ｅ（ｎ９＝Ｉｍ（ΣＳ（ｎ））／Ｒｅ（ΣＳ（ｎ）） ・・・・・・・（３） 除算部３０５から出力される誤差信号Ｅ（ｎ）は、式
（４）、式（５）のように対応付ける事ができる。 Ｅ（ｎ）＝ｔａｎ（ΔθEST) ・・・・・・・・（４） ΔθEST ＝２π・６４ＴΔｆEST ・・・・・・・・（５）
【０００８】ΔθESTは、受信したプリアンブル１０１
から推定されたｓ（ｎ）とｓ（ｎ−６４）の位相ずれで
あり、ΔｆESTは周波数オフセット推定値である。した
がって、誤差信号Ｅ（ｎ）から、式（６）の演算を行う
ことで、周波数オフセット推定値ΔｆESTが得られる。
図３の誤差信号変換部３１は、式（６）の演算を行うも
のである。 ΔｆEST＝ｔａｎ-1（Ｅ（ｎ））／（２π・６４Ｔ） ・・・・・・・（６） 以上が従来の周波数同期手法の説明である。
から推定されたｓ（ｎ）とｓ（ｎ−６４）の位相ずれで
あり、ΔｆESTは周波数オフセット推定値である。した
がって、誤差信号Ｅ（ｎ）から、式（６）の演算を行う
ことで、周波数オフセット推定値ΔｆESTが得られる。
図３の誤差信号変換部３１は、式（６）の演算を行うも
のである。 ΔｆEST＝ｔａｎ-1（Ｅ（ｎ））／（２π・６４Ｔ） ・・・・・・・（６） 以上が従来の周波数同期手法の説明である。
【０００９】
【発明が解決しようとする課題】従来の周波数同期手法
は、前述したように繰り返し送信される既知の信号列を
用いて周波数同期を行うものであり、誤差信号は式
（４）、式（５）のように、周波数オフセット値の正接
として求められる。正接関数の周期性から、誤差信号の
位相成分で表されるΔθESTがとりうる値の範囲は、−
π／２≦ΔθEST≦π／２に限定される。よって、式
（５）より、推定可能な周波数オフセット値の範囲は、
式（７）で表される。 −１／（２５６・Ｔ）≦Δｆ≦１／（２５６・Ｔ） ・・・・・・・・（７）
は、前述したように繰り返し送信される既知の信号列を
用いて周波数同期を行うものであり、誤差信号は式
（４）、式（５）のように、周波数オフセット値の正接
として求められる。正接関数の周期性から、誤差信号の
位相成分で表されるΔθESTがとりうる値の範囲は、−
π／２≦ΔθEST≦π／２に限定される。よって、式
（５）より、推定可能な周波数オフセット値の範囲は、
式（７）で表される。 −１／（２５６・Ｔ）≦Δｆ≦１／（２５６・Ｔ） ・・・・・・・・（７）
【００１０】すなわち、同期引き込み周波数範囲は±１
／（２５６・Ｔ）に限定され、これ以上の周波数オフセ
ットが存在する場合は周波数同期が不可能となるという
第１の課題がある。例えば、サンプリング周波数を１／
Ｔ＝２０ＭＨｚとした場合、同期引き込み周波数範囲は
±７８．５ＫＨｚとなる。
／（２５６・Ｔ）に限定され、これ以上の周波数オフセ
ットが存在する場合は周波数同期が不可能となるという
第１の課題がある。例えば、サンプリング周波数を１／
Ｔ＝２０ＭＨｚとした場合、同期引き込み周波数範囲は
±７８．５ＫＨｚとなる。
【００１１】さらに、周波数オフセットが同期引き込み
周波数範囲の限界近くの値を持つ場合、ΔθESTの大き
さはπ／２に近い値となる。このとき、除算部３０５の
入力ΣＳ（ｎ）のうち、実数成分Ｒｅ（ΣＳ（ｎ））の
絶対値の大きさが非常に小さいものとなる。このため、
伝送路ノイズによって、実数成分の符号が容易に反転
し、得られたΔｆESTの符号も反転してしまうことか
ら、周波数オフセットの推定誤差が大きくなるという第
２の課題がある。
周波数範囲の限界近くの値を持つ場合、ΔθESTの大き
さはπ／２に近い値となる。このとき、除算部３０５の
入力ΣＳ（ｎ）のうち、実数成分Ｒｅ（ΣＳ（ｎ））の
絶対値の大きさが非常に小さいものとなる。このため、
伝送路ノイズによって、実数成分の符号が容易に反転
し、得られたΔｆESTの符号も反転してしまうことか
ら、周波数オフセットの推定誤差が大きくなるという第
２の課題がある。
【００１２】本発明は、上記のような課題を解決し、広
い周波数同期引き込み範囲を実現し、周波数オフセット
の推定誤差を小さく抑えることのできる周波数同期方法
とそれを用いた受信装置を提供するとを目的とする。こ
の発明の前記ならびにそのほかの目的と新規な特徴は、
本明細書の記述および添付図面から明らかになるであろ
う。
い周波数同期引き込み範囲を実現し、周波数オフセット
の推定誤差を小さく抑えることのできる周波数同期方法
とそれを用いた受信装置を提供するとを目的とする。こ
の発明の前記ならびにそのほかの目的と新規な特徴は、
本明細書の記述および添付図面から明らかになるであろ
う。
【００１３】
【課題を解決するための手段】本願において開示される
発明のうち代表的なものの概要を簡単に説明すれば、下
記の通りである。すなわち、遅延量が異なる複数の相関
値計算部と、それぞれの相関値計算部から得られた相関
値より求める誤差信号と、設定された判定閾値の大小比
較を行い、周波数同期に使用する誤差信号を選択する。
前記複数の相関値の属する複素平面上の象限を検出する
ことで、周波数オフセット推定に用いる相関値を選択
し、選択された相関値と該相関値の属する複素平面上の
象限から周波数オフセット推定値を求める。上述のそれ
ぞれの誤差信号生成部において、遅延量が等しい複数の
相関値を計算し、これらを平均化して誤差信号を求め
る。上述のそれぞれの誤差信号生成部における除算部に
おいて、除算前に実数成分の入力の絶対値を求め、この
絶対値をもって虚数成分の除算を行う。
発明のうち代表的なものの概要を簡単に説明すれば、下
記の通りである。すなわち、遅延量が異なる複数の相関
値計算部と、それぞれの相関値計算部から得られた相関
値より求める誤差信号と、設定された判定閾値の大小比
較を行い、周波数同期に使用する誤差信号を選択する。
前記複数の相関値の属する複素平面上の象限を検出する
ことで、周波数オフセット推定に用いる相関値を選択
し、選択された相関値と該相関値の属する複素平面上の
象限から周波数オフセット推定値を求める。上述のそれ
ぞれの誤差信号生成部において、遅延量が等しい複数の
相関値を計算し、これらを平均化して誤差信号を求め
る。上述のそれぞれの誤差信号生成部における除算部に
おいて、除算前に実数成分の入力の絶対値を求め、この
絶対値をもって虚数成分の除算を行う。
【００１４】
【発明の実施の形態】以下に添付図面を参照して、本発
明に係る周波数同期方法およびこれを用いる受信装置の
実施の形態を詳細に説明する。
明に係る周波数同期方法およびこれを用いる受信装置の
実施の形態を詳細に説明する。
【００１５】（実施例１）本発明に係る周波数オフセッ
ト推定部の構成例を図１に示す。この周波数オフセット
推定部は、受信信号ｓ（ｎ）が入力され、これからそれ
ぞれ６４サンプル、３２サンプル、１６サンプル遅れて
受信する信号ｓ（ｎ−６４）、ｓ（ｎ−３２）、ｓ（ｎ
−１６）とｓ（ｎ）との相関値を求める相関値計算部５
１，５２，５３と、相関値計算部５１，５２，５３から
出力される相関値ΣＳ（ｎ）６４，ΣＳ（ｎ）３２，Σ
Ｓ（ｎ）１６を入力して、これらの相関値に対する誤差
信号Ｅ（ｎ）６４，Ｅ（ｎ）３２，Ｅ（ｎ）１６を求
め、周波数オフセット推定誤差が最小となるように、Ｅ
（ｎ）６４，Ｅ（ｎ）３２，Ｅ（ｎ）１６の中から推定
に用いる誤差信号Ｅ（ｎ）を選択して、さらにこのＥ
（ｎ）からΔθESTを求め、ΔθESTから周波数オフセッ
ト推定値ΔfESTを計算して出力する判定・演算部５４か
ら構成される。
ト推定部の構成例を図１に示す。この周波数オフセット
推定部は、受信信号ｓ（ｎ）が入力され、これからそれ
ぞれ６４サンプル、３２サンプル、１６サンプル遅れて
受信する信号ｓ（ｎ−６４）、ｓ（ｎ−３２）、ｓ（ｎ
−１６）とｓ（ｎ）との相関値を求める相関値計算部５
１，５２，５３と、相関値計算部５１，５２，５３から
出力される相関値ΣＳ（ｎ）６４，ΣＳ（ｎ）３２，Σ
Ｓ（ｎ）１６を入力して、これらの相関値に対する誤差
信号Ｅ（ｎ）６４，Ｅ（ｎ）３２，Ｅ（ｎ）１６を求
め、周波数オフセット推定誤差が最小となるように、Ｅ
（ｎ）６４，Ｅ（ｎ）３２，Ｅ（ｎ）１６の中から推定
に用いる誤差信号Ｅ（ｎ）を選択して、さらにこのＥ
（ｎ）からΔθESTを求め、ΔθESTから周波数オフセッ
ト推定値ΔfESTを計算して出力する判定・演算部５４か
ら構成される。
【００１６】図１において、相関値計算部５１，５２，
５３で、遅延量がそれぞれ６４サンプル、３２サンプ
ル、１６サンプルとなる遅延信号ｓ（ｎ−６４）、ｓ
（ｎ−３２）、ｓ（ｎ−１６）と、受信信号ｓ（ｎ）と
の相関値ΣＳ（ｎ）６４，ΣＳ（ｎ）３２，ΣＳ（ｎ）
１６を計算する。これらの計算結果を判定・演算部５４
に入力し、判定・演算部５４では、ΣＳ（ｎ）６４，Σ
Ｓ（ｎ）３２，ΣＳ（ｎ）１６の３種類の相関値に対応
する誤差信号Ｅ（ｎ）６４，Ｅ（ｎ）３２，Ｅ（ｎ）１
６を求め、Ｅ（ｎ）６４，Ｅ（ｎ）３２，Ｅ（ｎ）１６
の中から最終的な周波数オフセット推定に用いる誤差信
号を選択し、選択された誤差信号から位相ずれΔθEST
を計算し、選択された誤差信号に応じて、式（８）を用
いて周波数オフセット推定値ΔfESTを出力する。ここで
Ｎは、選択された相関値に対応する遅延サンプル数の値
であり、本実施例では１６，３２，６４のいずれかの値
を取る。
５３で、遅延量がそれぞれ６４サンプル、３２サンプ
ル、１６サンプルとなる遅延信号ｓ（ｎ−６４）、ｓ
（ｎ−３２）、ｓ（ｎ−１６）と、受信信号ｓ（ｎ）と
の相関値ΣＳ（ｎ）６４，ΣＳ（ｎ）３２，ΣＳ（ｎ）
１６を計算する。これらの計算結果を判定・演算部５４
に入力し、判定・演算部５４では、ΣＳ（ｎ）６４，Σ
Ｓ（ｎ）３２，ΣＳ（ｎ）１６の３種類の相関値に対応
する誤差信号Ｅ（ｎ）６４，Ｅ（ｎ）３２，Ｅ（ｎ）１
６を求め、Ｅ（ｎ）６４，Ｅ（ｎ）３２，Ｅ（ｎ）１６
の中から最終的な周波数オフセット推定に用いる誤差信
号を選択し、選択された誤差信号から位相ずれΔθEST
を計算し、選択された誤差信号に応じて、式（８）を用
いて周波数オフセット推定値ΔfESTを出力する。ここで
Ｎは、選択された相関値に対応する遅延サンプル数の値
であり、本実施例では１６，３２，６４のいずれかの値
を取る。
【００１７】
ΔｆEST = ΔθEST／２πＴＮ ・・・・・・・・（８）
相関値計算部５１，５２，５３の構成を図６に示す。相
関値計算部５１は、受信信号ｓ（ｎ）を６４サンプル遅
延させてｓ（ｎ−６４）を出力する遅延部５１１と、ｓ
（ｎ−６４）の複素共役ｓ（ｎ−６４）*を得る複素共
役部５１２と、ｓ（ｎ）とｓ（ｎ−６４）*の複素乗算
を行って積Ｓ（ｎ）６４を得る複素乗算部５１３と、Ｓ
（ｎ）６４を６４サンプルにわたって積算し、相関値Σ
Ｓ（ｎ）６４を出力する積算部５１４から構成される。
関値計算部５１は、受信信号ｓ（ｎ）を６４サンプル遅
延させてｓ（ｎ−６４）を出力する遅延部５１１と、ｓ
（ｎ−６４）の複素共役ｓ（ｎ−６４）*を得る複素共
役部５１２と、ｓ（ｎ）とｓ（ｎ−６４）*の複素乗算
を行って積Ｓ（ｎ）６４を得る複素乗算部５１３と、Ｓ
（ｎ）６４を６４サンプルにわたって積算し、相関値Σ
Ｓ（ｎ）６４を出力する積算部５１４から構成される。
【００１８】また、相関値計算部５２は、受信信号ｓ
（ｎ）を３２サンプル遅延させてｓ（ｎ−３２）を出力
する遅延部５２１と、ｓ（ｎ−３２）の複素共役ｓ（ｎ
−３２）*を得る複素共役部５２２と、ｓ（ｎ）とｓ
（ｎ−３２）*の複素乗算を行って積Ｓ（ｎ）３２を得
る複素乗算部５２３と、Ｓ（ｎ）３２を３２サンプルに
わたって積算し、相関値ΣＳ（ｎ）３２を出力する積算
部５２４から構成される。
（ｎ）を３２サンプル遅延させてｓ（ｎ−３２）を出力
する遅延部５２１と、ｓ（ｎ−３２）の複素共役ｓ（ｎ
−３２）*を得る複素共役部５２２と、ｓ（ｎ）とｓ
（ｎ−３２）*の複素乗算を行って積Ｓ（ｎ）３２を得
る複素乗算部５２３と、Ｓ（ｎ）３２を３２サンプルに
わたって積算し、相関値ΣＳ（ｎ）３２を出力する積算
部５２４から構成される。
【００１９】また、相関値計算部５３は、受信信号ｓ
（ｎ）を１６サンプル遅延させてｓ（ｎ−１６）を出力
する遅延部５３１と、ｓ（ｎ−１６）の複素共役ｓ（ｎ
−１６）*を得る複素共役部５３２と、ｓ（ｎ）とｓ
（ｎ−１６）*の複素乗算を行って積Ｓ（ｎ）１６を得
る複素乗算部５３３と、Ｓ（ｎ）１６を１６サンプルに
わたって積算し、相関値ΣＳ（ｎ）１６を出力する積算
部５３４から構成される。
（ｎ）を１６サンプル遅延させてｓ（ｎ−１６）を出力
する遅延部５３１と、ｓ（ｎ−１６）の複素共役ｓ（ｎ
−１６）*を得る複素共役部５３２と、ｓ（ｎ）とｓ
（ｎ−１６）*の複素乗算を行って積Ｓ（ｎ）１６を得
る複素乗算部５３３と、Ｓ（ｎ）１６を１６サンプルに
わたって積算し、相関値ΣＳ（ｎ）１６を出力する積算
部５３４から構成される。
【００２０】これらの相関値計算部の動作を、相関値計
算部５２を例にとって説明する。上記遅延部５２１で
は、入力信号ｓ（ｎ）に対して３２サンプルの遅延を施
し、ｓ（ｎ−３２）を得る。上記複素共役部５２２でｓ
（ｎ−３２）の複素共役ｓ（ｎ−３２）*を求めて、上
記複素乗算部５２３にてｓ（ｎ）とｓ（ｎ−３２）*の
乗算を行い、積Ｓ（ｎ）３２を得る。上記積算部５２４
では、このＳ（ｎ）３２を３２サンプルにわたって積算
し、相関値ΣＳ（ｎ）３２を求める。上記相関値計算部
５１および５３においても、遅延量と積算サンプル数が
それぞれ６４サンプル、１６サンプルとなる以外は同一
の動作を行う。
算部５２を例にとって説明する。上記遅延部５２１で
は、入力信号ｓ（ｎ）に対して３２サンプルの遅延を施
し、ｓ（ｎ−３２）を得る。上記複素共役部５２２でｓ
（ｎ−３２）の複素共役ｓ（ｎ−３２）*を求めて、上
記複素乗算部５２３にてｓ（ｎ）とｓ（ｎ−３２）*の
乗算を行い、積Ｓ（ｎ）３２を得る。上記積算部５２４
では、このＳ（ｎ）３２を３２サンプルにわたって積算
し、相関値ΣＳ（ｎ）３２を求める。上記相関値計算部
５１および５３においても、遅延量と積算サンプル数が
それぞれ６４サンプル、１６サンプルとなる以外は同一
の動作を行う。
【００２１】次に、図１中の判定・演算部５４の構成に
ついて、図７を用いて説明する。判定・演算部５４は、
入力した３種類の相関値、ΣＳ（ｎ）６４、ΣＳ（ｎ）
３２、ΣＳ（ｎ）１６のそれぞれについて、位相検出の
ために、虚数部を実数部で除算して、それぞれに対する
誤差信号Ｅ（ｎ）６４，Ｅ（ｎ）３２，Ｅ（ｎ）１６を
出力する除算部５４１、５４２、５４３と、除算部５４
１、５４２、５４３から出力された誤差信号から、周波
数オフセット推定に用いる誤差信号を選択し、これをＥ
（ｎ）として出力する誤差信号選択部５４４と、Ｅ
（ｎ）にアークタンジェント演算を行い、Ｅ（ｎ）から
ΔθESTを求めて出力するアークタンジェント演算部５
４５と、ΔθESTとＮから、式（８）を用いて周波数オ
フセット推定値ΔfESTを計算する演算部５４６から構成
される。なお、Ｎは誤差信号選択部５４４により、演算
部５４６へ供給される。
ついて、図７を用いて説明する。判定・演算部５４は、
入力した３種類の相関値、ΣＳ（ｎ）６４、ΣＳ（ｎ）
３２、ΣＳ（ｎ）１６のそれぞれについて、位相検出の
ために、虚数部を実数部で除算して、それぞれに対する
誤差信号Ｅ（ｎ）６４，Ｅ（ｎ）３２，Ｅ（ｎ）１６を
出力する除算部５４１、５４２、５４３と、除算部５４
１、５４２、５４３から出力された誤差信号から、周波
数オフセット推定に用いる誤差信号を選択し、これをＥ
（ｎ）として出力する誤差信号選択部５４４と、Ｅ
（ｎ）にアークタンジェント演算を行い、Ｅ（ｎ）から
ΔθESTを求めて出力するアークタンジェント演算部５
４５と、ΔθESTとＮから、式（８）を用いて周波数オ
フセット推定値ΔfESTを計算する演算部５４６から構成
される。なお、Ｎは誤差信号選択部５４４により、演算
部５４６へ供給される。
【００２２】相関値の位相成分を求めるために、相関値
の虚数部を実数部で除算する。得られた３種類の相関値
ΣＳ（ｎ）６４、ΣＳ（ｎ）３２、ΣＳ（ｎ）１６につ
いてこの除算をそれぞれ除算部５４１，５４２，５４３
を用いて行い、それぞれの遅延量に対する誤差信号Ｅ
（ｎ）６４，Ｅ（ｎ）３２，Ｅ（ｎ）１６を出力する。
誤差信号選択部５４４では、３種類の誤差信号Ｅ（ｎ）
６４，Ｅ（ｎ）３２，Ｅ（ｎ）１６を入力し、最も周波
数オフセット推定精度が良好となるものを選択して、Ｅ
（ｎ）としてアークタンジェント演算部５４５に出力す
る。
の虚数部を実数部で除算する。得られた３種類の相関値
ΣＳ（ｎ）６４、ΣＳ（ｎ）３２、ΣＳ（ｎ）１６につ
いてこの除算をそれぞれ除算部５４１，５４２，５４３
を用いて行い、それぞれの遅延量に対する誤差信号Ｅ
（ｎ）６４，Ｅ（ｎ）３２，Ｅ（ｎ）１６を出力する。
誤差信号選択部５４４では、３種類の誤差信号Ｅ（ｎ）
６４，Ｅ（ｎ）３２，Ｅ（ｎ）１６を入力し、最も周波
数オフセット推定精度が良好となるものを選択して、Ｅ
（ｎ）としてアークタンジェント演算部５４５に出力す
る。
【００２３】また、ΔθESTから周波数オフセット推定
値ΔｆESTを求めるため、選択された誤差信号に対する
遅延サンプル数Ｎを、演算部５４６に出力する。このよ
うにして得られたＥ（ｎ）をアークタンジェント演算部
５４５に入力し、アークタンジェント演算部５４５で
は、誤差信号選択部５４４から出力されたＥ（ｎ）に対
して式（９）の演算を行い、ΔθESTを求め、これを演
算部５４６に出力する。 ΔθEST = ｔａｎ-1（Ｅ（ｎ）） ・・・・・・・・（９）
値ΔｆESTを求めるため、選択された誤差信号に対する
遅延サンプル数Ｎを、演算部５４６に出力する。このよ
うにして得られたＥ（ｎ）をアークタンジェント演算部
５４５に入力し、アークタンジェント演算部５４５で
は、誤差信号選択部５４４から出力されたＥ（ｎ）に対
して式（９）の演算を行い、ΔθESTを求め、これを演
算部５４６に出力する。 ΔθEST = ｔａｎ-1（Ｅ（ｎ）） ・・・・・・・・（９）
【００２４】演算部５４６では、アークタンジェント演
算部５４５より入力したΔθESTと、誤差信号選択部５
４４より入力した遅延サンプル数Ｎを用いて、式（８）
の演算を行い、ΔfESTを求める。
算部５４５より入力したΔθESTと、誤差信号選択部５
４４より入力した遅延サンプル数Ｎを用いて、式（８）
の演算を行い、ΔfESTを求める。
【００２５】次に、誤差信号選択部５４４の、誤差信号
選択アルゴリズムについて説明する。Ｅ（ｎ）は式
（３）および式（４）で表されることから、ΔθESTが
−π／２≦ΔθEST≦＋π／２の範囲にある場合、Ｅ
（ｎ）の値は周波数オフセットが大きくなると共に大き
くなる。また、式（４）、式（５）、式（６）から、遅
延量を少なくすると、引き込み周波数範囲は拡大する。
一方、遅延量を多くすると、ノイズ成分が抑圧され、引
き込み周波数範囲内の推定精度が良好になる。これらの
性質を利用して、Ｅ（ｎ）１６とＥ（ｎ）３２の値の大
きさを観測し、ΔθESTを求めるために使用する誤差信
号Ｅ（ｎ）を選択する。
選択アルゴリズムについて説明する。Ｅ（ｎ）は式
（３）および式（４）で表されることから、ΔθESTが
−π／２≦ΔθEST≦＋π／２の範囲にある場合、Ｅ
（ｎ）の値は周波数オフセットが大きくなると共に大き
くなる。また、式（４）、式（５）、式（６）から、遅
延量を少なくすると、引き込み周波数範囲は拡大する。
一方、遅延量を多くすると、ノイズ成分が抑圧され、引
き込み周波数範囲内の推定精度が良好になる。これらの
性質を利用して、Ｅ（ｎ）１６とＥ（ｎ）３２の値の大
きさを観測し、ΔθESTを求めるために使用する誤差信
号Ｅ（ｎ）を選択する。
【００２６】誤差信号選択アルゴリズムを図８に示す。
まず、判定閾値を設定し、処理ブロック８１にてＥ
（ｎ）１６の絶対値との大小比較を行う。Ｅ（ｎ）１６
の絶対値が判定閾値より大きい場合、周波数オフセット
はＥ（ｎ）３２で推定した場合の引き込み周波数範囲内
にないと判断し、最も引き込み周波数範囲の広いＥ
（ｎ）１６を用いて推定を行う（処理ブロック８５）。
また、判定閾値の方がＥ（ｎ）１６と比べて大きいか等
しい場合、次に処理ブロック８２にてＥ（ｎ）３２の絶
対値と判定閾値の大小比較を行う。
まず、判定閾値を設定し、処理ブロック８１にてＥ
（ｎ）１６の絶対値との大小比較を行う。Ｅ（ｎ）１６
の絶対値が判定閾値より大きい場合、周波数オフセット
はＥ（ｎ）３２で推定した場合の引き込み周波数範囲内
にないと判断し、最も引き込み周波数範囲の広いＥ
（ｎ）１６を用いて推定を行う（処理ブロック８５）。
また、判定閾値の方がＥ（ｎ）１６と比べて大きいか等
しい場合、次に処理ブロック８２にてＥ（ｎ）３２の絶
対値と判定閾値の大小比較を行う。
【００２７】ここで、Ｅ（ｎ）３２の絶対値の方が判定
閾値よりも大きい場合、周波数オフセットはＥ（ｎ）６
４で推定した場合の引き込み周波数範囲内にないと判断
し、二番目に引き込み範囲の広いＥ（ｎ）３２を用いて
推定を行う（処理ブロック８４）。また、処理ブロック
８２において判定閾値の方がＥ（ｎ）３２より大きいか
等しい場合、周波数オフセットはＥ（ｎ）６４で推定し
た場合の引き込み周波数範囲内にあると判断し、最も推
定精度が良好となるＥ（ｎ）６４を用いて推定を行う
（処理ブロック８３）。
閾値よりも大きい場合、周波数オフセットはＥ（ｎ）６
４で推定した場合の引き込み周波数範囲内にないと判断
し、二番目に引き込み範囲の広いＥ（ｎ）３２を用いて
推定を行う（処理ブロック８４）。また、処理ブロック
８２において判定閾値の方がＥ（ｎ）３２より大きいか
等しい場合、周波数オフセットはＥ（ｎ）６４で推定し
た場合の引き込み周波数範囲内にあると判断し、最も推
定精度が良好となるＥ（ｎ）６４を用いて推定を行う
（処理ブロック８３）。
【００２８】（実施例２）本発明の周波数同期方法に係
る相関値計算部の別の実施の形態について述べる。相関
値計算部５２は、図９に示すような構成とすることもで
きる。図９に示す相関値計算部では、前述図６の遅延部
５２１、複素共役部５２２、複素乗算部５２３、積算部
５２４に加えて、受信信号に対して６４サンプルの遅延
を与えて、ｓ（ｎ−６４）を出力する遅延部５２５、ｓ
（ｎ−６４）に対して３２サンプルの遅延を与え、ｓ
（ｎ−９６）を出力する遅延部５２６、ｓ（ｎ−９６）
の複素共役を取り、ｓ（ｎ−９６）*を出力する複素共
役部５２７、ｓ（ｎ−６４）と、ｓ（ｎ−９６）*を複
素乗算して、積Ｓ（ｎ）３２''を出力する複素乗算部５
２８、Ｓ（ｎ）３２''を３２サンプルにわたって積算し
て、相関値ΣＳ（ｎ）３２''を出力する積算部５２９、
および、積算部５２４の出力である相関値ΣＳ（ｎ）３
２' と、積算部５２９の出力である相関値ΣＳ（ｎ）３
２''を平均化して、推定に用いる相関値ΣＳ（ｎ）３２
を出力する平均化部５２０で構成される。
る相関値計算部の別の実施の形態について述べる。相関
値計算部５２は、図９に示すような構成とすることもで
きる。図９に示す相関値計算部では、前述図６の遅延部
５２１、複素共役部５２２、複素乗算部５２３、積算部
５２４に加えて、受信信号に対して６４サンプルの遅延
を与えて、ｓ（ｎ−６４）を出力する遅延部５２５、ｓ
（ｎ−６４）に対して３２サンプルの遅延を与え、ｓ
（ｎ−９６）を出力する遅延部５２６、ｓ（ｎ−９６）
の複素共役を取り、ｓ（ｎ−９６）*を出力する複素共
役部５２７、ｓ（ｎ−６４）と、ｓ（ｎ−９６）*を複
素乗算して、積Ｓ（ｎ）３２''を出力する複素乗算部５
２８、Ｓ（ｎ）３２''を３２サンプルにわたって積算し
て、相関値ΣＳ（ｎ）３２''を出力する積算部５２９、
および、積算部５２４の出力である相関値ΣＳ（ｎ）３
２' と、積算部５２９の出力である相関値ΣＳ（ｎ）３
２''を平均化して、推定に用いる相関値ΣＳ（ｎ）３２
を出力する平均化部５２０で構成される。
【００２９】図９に示す相関値計算部では、図６に示す
相関値計算部５２の動作を、受信信号ｓ（ｎ）に対して
だけではなく、受信信号を６４サンプル遅延させて得ら
れた信号ｓ（ｎ−６４）についても行い、平均化部５２
０にて、得られた二つの相関値を加算し、定数２で除算
することで相関値が平均化され、伝送路ノイズの影響を
軽減し、推定精度を向上できる。ただし、相関値計算部
５２の出力ΣＳ（ｎ）３２に対しては、図１中の判定・
演算部５４にて、虚数部を実数部で除算するため、平均
化部５２０では、ΣＳ（ｎ）３２' とΣＳ（ｎ）３２''
の加算を行うだけで、平均化することと等価になる。こ
の処理の説明図を図２１に示す。
相関値計算部５２の動作を、受信信号ｓ（ｎ）に対して
だけではなく、受信信号を６４サンプル遅延させて得ら
れた信号ｓ（ｎ−６４）についても行い、平均化部５２
０にて、得られた二つの相関値を加算し、定数２で除算
することで相関値が平均化され、伝送路ノイズの影響を
軽減し、推定精度を向上できる。ただし、相関値計算部
５２の出力ΣＳ（ｎ）３２に対しては、図１中の判定・
演算部５４にて、虚数部を実数部で除算するため、平均
化部５２０では、ΣＳ（ｎ）３２' とΣＳ（ｎ）３２''
の加算を行うだけで、平均化することと等価になる。こ
の処理の説明図を図２１に示す。
【００３０】図２中の固定信号系列２０と２１からなる
３２サンプルの信号系列Ｓ１（ｎ）と、図２中の固定信
号系列２２と２３からなる３２サンプルの信号系列Ｓ２
（ｎ）との相関値ΣＳ（ｎ）３２' を処理ブロック７１
にて求め、さらに、図２中の固定信号系列２４と２５か
らなる３２サンプル信号系列Ｓ３（ｎ）と、図２中の固
定信号系列２６と２７からなる３２サンプルの信号系列
Ｓ４（ｎ）との相関値ΣＳ（ｎ）３２''を処理ブロック
７２にて求める。処理ブロック７３にて、ΣＳ（ｎ）３
２' とΣＳ（ｎ）３２''を平均化（加算のみでも良い）
して、周波数オフセット推定に用いる相関値ΣＳ（ｎ）
３２を得る。以上が本発明に係る第二の実施例である。
３２サンプルの信号系列Ｓ１（ｎ）と、図２中の固定信
号系列２２と２３からなる３２サンプルの信号系列Ｓ２
（ｎ）との相関値ΣＳ（ｎ）３２' を処理ブロック７１
にて求め、さらに、図２中の固定信号系列２４と２５か
らなる３２サンプル信号系列Ｓ３（ｎ）と、図２中の固
定信号系列２６と２７からなる３２サンプルの信号系列
Ｓ４（ｎ）との相関値ΣＳ（ｎ）３２''を処理ブロック
７２にて求める。処理ブロック７３にて、ΣＳ（ｎ）３
２' とΣＳ（ｎ）３２''を平均化（加算のみでも良い）
して、周波数オフセット推定に用いる相関値ΣＳ（ｎ）
３２を得る。以上が本発明に係る第二の実施例である。
【００３１】（実施例３）本発明の周波数同期方法に係
る、判定・演算部の別の実施の形態について、図１０を
用いて説明する。図１０において、アークタンジェント
演算部５４５および演算部５４６の動作は図７に示すも
のと同一であり、同一の符号を付し、説明を省略する。
る、判定・演算部の別の実施の形態について、図１０を
用いて説明する。図１０において、アークタンジェント
演算部５４５および演算部５４６の動作は図７に示すも
のと同一であり、同一の符号を付し、説明を省略する。
【００３２】図１０では、判定・演算部５４は上記アー
クタンジェント演算部５４５と、演算部５４６と、周波
数オフセット推定に用いる相関値を選択する相関値選択
部５４７と、相関値選択部５４７で選択された相関値の
虚数部を実数部で除算する除算部５４８にて構成され
る。
クタンジェント演算部５４５と、演算部５４６と、周波
数オフセット推定に用いる相関値を選択する相関値選択
部５４７と、相関値選択部５４７で選択された相関値の
虚数部を実数部で除算する除算部５４８にて構成され
る。
【００３３】さらに、相関値選択部５４７の構成を図１
１に示す。相関値選択部５４７は、入力したΣＳ（ｎ）
３２の実数部Ｒｅ（ΣＳ（ｎ））３２とΣＳ（ｎ）１６
の実数部Ｒｅ（ΣＳ（ｎ））１６にそれぞれ判定閾値を
乗算して、積Ｒｅ（ΣＳ（ｎ））３２ＴＨとＲｅ（ΣＳ
（ｎ））１６ＴＨをそれぞれ出力する乗算部８０１、８
０２と、Ｉｍ（ΣＳ（ｎ））３２Ｉｍ（ΣＳ（ｎ））１
６と、Ｒｅ（ΣＳ（ｎ））３２ＴＨとＲｅ（ΣＳ
（ｎ））１６ＴＨを入力して、ΣＳ（ｎ）６４、ΣＳ
（ｎ）３２、ΣＳ（ｎ）１６のうちのどれを周波数オフ
セット推定に用いるΣＳ（ｎ）とするかとＮを決定する
判定部８０３から構成される。
１に示す。相関値選択部５４７は、入力したΣＳ（ｎ）
３２の実数部Ｒｅ（ΣＳ（ｎ））３２とΣＳ（ｎ）１６
の実数部Ｒｅ（ΣＳ（ｎ））１６にそれぞれ判定閾値を
乗算して、積Ｒｅ（ΣＳ（ｎ））３２ＴＨとＲｅ（ΣＳ
（ｎ））１６ＴＨをそれぞれ出力する乗算部８０１、８
０２と、Ｉｍ（ΣＳ（ｎ））３２Ｉｍ（ΣＳ（ｎ））１
６と、Ｒｅ（ΣＳ（ｎ））３２ＴＨとＲｅ（ΣＳ
（ｎ））１６ＴＨを入力して、ΣＳ（ｎ）６４、ΣＳ
（ｎ）３２、ΣＳ（ｎ）１６のうちのどれを周波数オフ
セット推定に用いるΣＳ（ｎ）とするかとＮを決定する
判定部８０３から構成される。
【００３４】次に、判定部８０３における判定アルゴリ
ズムを、図１２を用いて説明する。まず、処理ブロック
８０４にて、Ｉｍ（ΣＳ（ｎ））１６の絶対値とＲｅ
（ΣＳ（ｎ））１６ＴＨの絶対値の大小比較を行う。Ｉ
ｍ（ΣＳ（ｎ））１６の絶対値の方がＲｅ（ΣＳ
（ｎ））１６ＴＨの絶対値よりも大きい場合、周波数オ
フセットはΣＳ（ｎ）３２で推定した場合の引き込み周
波数範囲内にないと判断し、最も引き込み周波数範囲の
広いΣＳ（ｎ）１６を用いて推定を行う（処理ブロック
８０８) 。
ズムを、図１２を用いて説明する。まず、処理ブロック
８０４にて、Ｉｍ（ΣＳ（ｎ））１６の絶対値とＲｅ
（ΣＳ（ｎ））１６ＴＨの絶対値の大小比較を行う。Ｉ
ｍ（ΣＳ（ｎ））１６の絶対値の方がＲｅ（ΣＳ
（ｎ））１６ＴＨの絶対値よりも大きい場合、周波数オ
フセットはΣＳ（ｎ）３２で推定した場合の引き込み周
波数範囲内にないと判断し、最も引き込み周波数範囲の
広いΣＳ（ｎ）１６を用いて推定を行う（処理ブロック
８０８) 。
【００３５】Ｒｅ（ΣＳ（ｎ））１６ＴＨの絶対値の方
が、Ｉｍ（ΣＳ（ｎ））１６の絶対値よりも大きいか等
しい場合は、次に、処理ブロック８０５にて、Ｉｍ（Σ
Ｓ（ｎ））３２の絶対値とＲｅ（ΣＳ（ｎ））３２ＴＨ
の絶対値の大小比較を行う。Ｉｍ（ΣＳ（ｎ））３２の
絶対値の方が、Ｒｅ（ΣＳ（ｎ））３２ＴＨの絶対値よ
りも大きい場合、周波数オフセットはΣＳ（ｎ）６４で
推定した場合の引き込み周波数範囲内にないと判断し、
二番目に引き込み周波数範囲の広いΣＳ（ｎ）３２を用
いて推定を行う（処理ブロック８０７）。
が、Ｉｍ（ΣＳ（ｎ））１６の絶対値よりも大きいか等
しい場合は、次に、処理ブロック８０５にて、Ｉｍ（Σ
Ｓ（ｎ））３２の絶対値とＲｅ（ΣＳ（ｎ））３２ＴＨ
の絶対値の大小比較を行う。Ｉｍ（ΣＳ（ｎ））３２の
絶対値の方が、Ｒｅ（ΣＳ（ｎ））３２ＴＨの絶対値よ
りも大きい場合、周波数オフセットはΣＳ（ｎ）６４で
推定した場合の引き込み周波数範囲内にないと判断し、
二番目に引き込み周波数範囲の広いΣＳ（ｎ）３２を用
いて推定を行う（処理ブロック８０７）。
【００３６】Ｒｅ（ΣＳ（ｎ））３２ＴＨの絶対値の方
がＩｍ（ΣＳ（ｎ））３２の絶対値よりも大きいか等し
い場合は、周波数オフセットはΣＳ（ｎ）６４で推定し
た場合の引き込み周波数範囲内にあると判断し、最も推
定精度が良好となるΣＳ（ｎ）６４を用いて推定を行う
（処理ブロック８０６）。この処理は、相関値の実数部
に判定閾値を乗じて該相関値の虚数部との大小比較を行
うもので、相関値の虚数部を該相関値の実数部で除算し
て、判定閾値と大小比較することと等価である。以上が
本発明に係る第三の実施例である。
がＩｍ（ΣＳ（ｎ））３２の絶対値よりも大きいか等し
い場合は、周波数オフセットはΣＳ（ｎ）６４で推定し
た場合の引き込み周波数範囲内にあると判断し、最も推
定精度が良好となるΣＳ（ｎ）６４を用いて推定を行う
（処理ブロック８０６）。この処理は、相関値の実数部
に判定閾値を乗じて該相関値の虚数部との大小比較を行
うもので、相関値の虚数部を該相関値の実数部で除算し
て、判定閾値と大小比較することと等価である。以上が
本発明に係る第三の実施例である。
【００３７】（実施例４）本発明の周波数同期方法に係
る判定・演算部の別の実施の形態について、図１３を用
いて説明する。図１３の構成において、判定・演算部５
４は、周波数オフセット推定に用いる相関値を判定する
相関値選択部５４９と、選択された相関値の虚数部Ｉｍ
（ΣＳ（ｎ））を実数部Ｒｅ（ΣＳ（ｎ））で除算して
誤差信号Ｅ（ｎ）を求める除算部５４８と、得られた誤
差信号Ｅ（ｎ）にアークタンジェント演算を施してΔθ
ESTを求めるアークタンジェント演算部５４５と、Δθ
ESTと以下で説明する象限情報とＮとから周波数オフセ
ットΔｆESTを計算する位相周波数変換部５４０から構
成される。ここで、除算部５４８の動作およびアークタ
ンジェント演算部５４５の動作は図１０に示すものと同
じであり、同一の符号を付し、説明を省略する。
る判定・演算部の別の実施の形態について、図１３を用
いて説明する。図１３の構成において、判定・演算部５
４は、周波数オフセット推定に用いる相関値を判定する
相関値選択部５４９と、選択された相関値の虚数部Ｉｍ
（ΣＳ（ｎ））を実数部Ｒｅ（ΣＳ（ｎ））で除算して
誤差信号Ｅ（ｎ）を求める除算部５４８と、得られた誤
差信号Ｅ（ｎ）にアークタンジェント演算を施してΔθ
ESTを求めるアークタンジェント演算部５４５と、Δθ
ESTと以下で説明する象限情報とＮとから周波数オフセ
ットΔｆESTを計算する位相周波数変換部５４０から構
成される。ここで、除算部５４８の動作およびアークタ
ンジェント演算部５４５の動作は図１０に示すものと同
じであり、同一の符号を付し、説明を省略する。
【００３８】以下、図１３に示す判定・演算部５４の動
作について説明する。まず、相関値計算部５１、５２、
５３でそれぞれ得られた相関値ΣＳ（ｎ９６４、ΣＳ
（ｎ）３２、ΣＳ（ｎ）１６は、相関値選択部５４９に
入力される。相関値選択部５４９では、入力した３種類
の相関値の中から、推定に用いる相関値を選択して、除
算部５４８へ出力し、さらに選択された相関値が属する
複素平面上の象限を決定し、象限情報として位相周波数
変換部５４０へ出力する。さらに、選択された相関値が
ΣＳ（ｎ）６４、ΣＳ（ｎ）３２、ΣＳ（ｎ）１６のい
ずれかに応じて、Ｎとして６４、３２、１６のいずれか
を出力する。
作について説明する。まず、相関値計算部５１、５２、
５３でそれぞれ得られた相関値ΣＳ（ｎ９６４、ΣＳ
（ｎ）３２、ΣＳ（ｎ）１６は、相関値選択部５４９に
入力される。相関値選択部５４９では、入力した３種類
の相関値の中から、推定に用いる相関値を選択して、除
算部５４８へ出力し、さらに選択された相関値が属する
複素平面上の象限を決定し、象限情報として位相周波数
変換部５４０へ出力する。さらに、選択された相関値が
ΣＳ（ｎ）６４、ΣＳ（ｎ）３２、ΣＳ（ｎ）１６のい
ずれかに応じて、Ｎとして６４、３２、１６のいずれか
を出力する。
【００３９】位相周波数変換部５４０では、アークタン
ジェント演算部５４５から得られたΔθESTと、相関値
選択部５４９から得られた遅延量Ｎと、相関値選択部５
４９より得られた象限情報を用いて、ΔθESTからΔｆ
ESTを計算し、出力する。
ジェント演算部５４５から得られたΔθESTと、相関値
選択部５４９から得られた遅延量Ｎと、相関値選択部５
４９より得られた象限情報を用いて、ΔθESTからΔｆ
ESTを計算し、出力する。
【００４０】相関値選択部５４９における、推定に使用
する相関値の選択アルゴリズムを図１４を用いて説明す
る。まず、処理ブロック１１１にて、Ｒｅ（ΣＳ
（ｎ））１６の正負を判定する。これが負の場合、ΣＳ
（ｎ）１６を用いて推定する（処理ブロック１１５）。
正または０の場合は処理ブロック１１２にて、Ｒｅ（Σ
Ｓ（ｎ））３２の正負を判定する。Ｒｅ（ΣＳ（ｎ））
３２が負の場合、ΣＳ（ｎ）３２を用いて推定する（処
理ブロック１１４）。Ｒｅ（ΣＳ（ｎ））３２が正また
は０の場合はΣＳ（ｎ）６４を用いて推定する（処理ブ
ロック１１３）。
する相関値の選択アルゴリズムを図１４を用いて説明す
る。まず、処理ブロック１１１にて、Ｒｅ（ΣＳ
（ｎ））１６の正負を判定する。これが負の場合、ΣＳ
（ｎ）１６を用いて推定する（処理ブロック１１５）。
正または０の場合は処理ブロック１１２にて、Ｒｅ（Σ
Ｓ（ｎ））３２の正負を判定する。Ｒｅ（ΣＳ（ｎ））
３２が負の場合、ΣＳ（ｎ）３２を用いて推定する（処
理ブロック１１４）。Ｒｅ（ΣＳ（ｎ））３２が正また
は０の場合はΣＳ（ｎ）６４を用いて推定する（処理ブ
ロック１１３）。
【００４１】また、相関値選択部５４９における象限情
報決定方法を図１５、図１６、図１７を用いて説明す
る。ここで、図１５はΣＳ（ｎ）６４を用いて推定する
場合の方法、図１６はΣＳ（ｎ）３２を用いて推定する
場合の方法、図１７はΣＳ（ｎ）１６を用いて推定する
場合の方法である。
報決定方法を図１５、図１６、図１７を用いて説明す
る。ここで、図１５はΣＳ（ｎ）６４を用いて推定する
場合の方法、図１６はΣＳ（ｎ）３２を用いて推定する
場合の方法、図１７はΣＳ（ｎ）１６を用いて推定する
場合の方法である。
【００４２】例えば、ΣＳ（ｎ）６４を用いて推定する
場合、Ｒｅ（ΣＳ（ｎ））６４とＩｍ（ΣＳ（ｎ））６
４の正負から象限情報を得る。Ｉｍ（ΣＳ（ｎ））６４
が正または０の場合、ΣＳ（ｎ）６４から得られるΔθ
EST(以下、Δθ６４と記す）の範囲は０≦Δθ６４≦π
であり、このとき、ΣＳ（ｎ）６４の属する複素平面上
の象限は第１象限または第２象限となる。
場合、Ｒｅ（ΣＳ（ｎ））６４とＩｍ（ΣＳ（ｎ））６
４の正負から象限情報を得る。Ｉｍ（ΣＳ（ｎ））６４
が正または０の場合、ΣＳ（ｎ）６４から得られるΔθ
EST(以下、Δθ６４と記す）の範囲は０≦Δθ６４≦π
であり、このとき、ΣＳ（ｎ）６４の属する複素平面上
の象限は第１象限または第２象限となる。
【００４３】また、Ｉｍ（ΣＳ（ｎ））６４が負の場
合、Δθ６４の範囲は−π＜Δθ６４＜０であり、ΣＳ
（ｎ）６４の属する複素平面上の象限は第３象限または
第４象限となる。また、Ｒｅ（ΣＳ（ｎ））６４が正ま
たは０の場合、Δθ６４の範囲は−π／２≦Δθ６４≦
π／２であり、このとき、ΣＳ（ｎ）６４の属する複素
平面上の象限は第１象限または第４象限となる。
合、Δθ６４の範囲は−π＜Δθ６４＜０であり、ΣＳ
（ｎ）６４の属する複素平面上の象限は第３象限または
第４象限となる。また、Ｒｅ（ΣＳ（ｎ））６４が正ま
たは０の場合、Δθ６４の範囲は−π／２≦Δθ６４≦
π／２であり、このとき、ΣＳ（ｎ）６４の属する複素
平面上の象限は第１象限または第４象限となる。
【００４４】一方、Ｒe （ΣＳ（ｎ））６４が負の場
合、Δθ６４の範囲は−π＜Δθ６４＜−π／２または
π／２＜Δθ６４＜πであり、このとき、ΣＳ（ｎ）６
４の属する複素平面上の象限は第２象限または第３象限
となる。以上のことから、図１５に示すように、Ｒｅ
（ΣＳ（ｎ））６４が正または０の場合は、Ｉｍ（ΣＳ
（ｎ））６４が正または０ならば第１象限、Ｉｍ（ΣＳ
（ｎ））６４が負ならば第４象限、Ｒｅ（ΣＳ（ｎ））
６４が負の場合は、Ｉｍ（ΣＳ（ｎ））６４が正または
０ならば第２象限、Ｉｍ（ΣＳ（ｎ））６４が負ならば
第３象限と判定することができる。
合、Δθ６４の範囲は−π＜Δθ６４＜−π／２または
π／２＜Δθ６４＜πであり、このとき、ΣＳ（ｎ）６
４の属する複素平面上の象限は第２象限または第３象限
となる。以上のことから、図１５に示すように、Ｒｅ
（ΣＳ（ｎ））６４が正または０の場合は、Ｉｍ（ΣＳ
（ｎ））６４が正または０ならば第１象限、Ｉｍ（ΣＳ
（ｎ））６４が負ならば第４象限、Ｒｅ（ΣＳ（ｎ））
６４が負の場合は、Ｉｍ（ΣＳ（ｎ））６４が正または
０ならば第２象限、Ｉｍ（ΣＳ（ｎ））６４が負ならば
第３象限と判定することができる。
【００４５】また、ΣＳ（ｎ）３２を用いて推定する場
合、図１６に示すように、Ｒｅ（ΣＳ（ｎ））３２とＩ
ｍ（ΣＳ（ｎ））３２の正負から象限情報を得る。すな
わち、Ｒｅ（ΣＳ（ｎ））３２が正または０の場合、Ｉ
ｍ（ΣＳ（ｎ））３２が正または０であれば第１象限、
Ｉｍ（ΣＳ（ｎ））３２が負であれば第４象限と判定す
る。また、Ｒｅ（ΣＳ（ｎ））３２が負の場合、Ｉｍ
（ΣＳ（ｎ））３２が正または０であれば第２象限、Ｉ
ｍ（ΣＳ（ｎ））３２が負であれば第３象限と判定す
る。
合、図１６に示すように、Ｒｅ（ΣＳ（ｎ））３２とＩ
ｍ（ΣＳ（ｎ））３２の正負から象限情報を得る。すな
わち、Ｒｅ（ΣＳ（ｎ））３２が正または０の場合、Ｉ
ｍ（ΣＳ（ｎ））３２が正または０であれば第１象限、
Ｉｍ（ΣＳ（ｎ））３２が負であれば第４象限と判定す
る。また、Ｒｅ（ΣＳ（ｎ））３２が負の場合、Ｉｍ
（ΣＳ（ｎ））３２が正または０であれば第２象限、Ｉ
ｍ（ΣＳ（ｎ））３２が負であれば第３象限と判定す
る。
【００４６】さらに、ΣＳ（ｎ）１６を用いて推定する
場合、図１７に示すように、Ｒｅ（ΣＳ（ｎ））１６と
Ｉｍ（ΣＳ（ｎ））１６の正負から象限情報を得る。す
なわち、Ｒｅ（ΣＳ（ｎ））１６が正または０の場合、
Ｉｍ（ΣＳ（ｎ））１６が正または０であれば第１象
限、Ｉｍ（ΣＳ（ｎ））１６が負であれば第４象限と判
定する。また、Ｒｅ（ΣＳ（ｎ））１６が負の場合、Ｉ
ｍ（ΣＳ（ｎ））１６が正または０であれば第２象限、
Ｉｍ（ΣＳ（ｎ））１６が負であれば第３象限と判定す
る。
場合、図１７に示すように、Ｒｅ（ΣＳ（ｎ））１６と
Ｉｍ（ΣＳ（ｎ））１６の正負から象限情報を得る。す
なわち、Ｒｅ（ΣＳ（ｎ））１６が正または０の場合、
Ｉｍ（ΣＳ（ｎ））１６が正または０であれば第１象
限、Ｉｍ（ΣＳ（ｎ））１６が負であれば第４象限と判
定する。また、Ｒｅ（ΣＳ（ｎ））１６が負の場合、Ｉ
ｍ（ΣＳ（ｎ））１６が正または０であれば第２象限、
Ｉｍ（ΣＳ（ｎ））１６が負であれば第３象限と判定す
る。
【００４７】次に、位相周波数変換部５４０の動作を説
明する。まず、相関値選択部５４９より得られた象限情
報から、選択された相関値の属する複素平面上の象限を
決定する。選択された相関値が複素平面上の第４象限ま
たは第１象限にある場合、すなわちΔθESTが−π／２
≦ΔθEST≦π／２の範囲にある場合は、これまでに述
べた方法、すなわち数８の演算によりΔｆESTが求めら
れる。ΔθESTが上記範囲内にないとき、すなわち、選
択された相関値が複素平面上の第２象限または第３象限
にある場合は、数８の演算ののち、オフセット値を加え
てΔｆESTとする。以上が本発明に係る第四の実施例で
ある。
明する。まず、相関値選択部５４９より得られた象限情
報から、選択された相関値の属する複素平面上の象限を
決定する。選択された相関値が複素平面上の第４象限ま
たは第１象限にある場合、すなわちΔθESTが−π／２
≦ΔθEST≦π／２の範囲にある場合は、これまでに述
べた方法、すなわち数８の演算によりΔｆESTが求めら
れる。ΔθESTが上記範囲内にないとき、すなわち、選
択された相関値が複素平面上の第２象限または第３象限
にある場合は、数８の演算ののち、オフセット値を加え
てΔｆESTとする。以上が本発明に係る第四の実施例で
ある。
【００４８】（実施例５）本発明の周波数同期方法に係
る除算部の別の実施の形態について、図１８を用いて説
明する。図１８において、除算部５４８は、Ｒｅ（ΣＳ
（ｎ））の絶対値 |Ｒｅ（ΣＳ（ｎ））|を求める絶対
値算出部９０１と、Ｉｍ（ΣＳ（ｎ））を、絶対値算出
部の出力 |Ｒｅ（ΣＳ（ｎ））|で除算し、誤差信号Ｅ
（ｎ）を求める除算部９０２から構成される。
る除算部の別の実施の形態について、図１８を用いて説
明する。図１８において、除算部５４８は、Ｒｅ（ΣＳ
（ｎ））の絶対値 |Ｒｅ（ΣＳ（ｎ））|を求める絶対
値算出部９０１と、Ｉｍ（ΣＳ（ｎ））を、絶対値算出
部の出力 |Ｒｅ（ΣＳ（ｎ））|で除算し、誤差信号Ｅ
（ｎ）を求める除算部９０２から構成される。
【００４９】前記実施例１乃至３の方法でΔθESTを求
めた場合、理想的な状態では、ΔθE STの取りうる範囲
が−π／２≦ΔθEST≦π／２に限られることから、相
関値ΣＳ（ｎ）の実数部は伝送路ノイズ成分を除けば常
に正または０となる。しかし、実際の受信信号にはノイ
ズ成分が加わるため、相関値の実数部の値が小さい場
合、負の値になることがある。
めた場合、理想的な状態では、ΔθE STの取りうる範囲
が−π／２≦ΔθEST≦π／２に限られることから、相
関値ΣＳ（ｎ）の実数部は伝送路ノイズ成分を除けば常
に正または０となる。しかし、実際の受信信号にはノイ
ズ成分が加わるため、相関値の実数部の値が小さい場
合、負の値になることがある。
【００５０】このことを、除算部５４８を例にとって、
図１９を用いて説明する。例えば、相関値ΣＳ（ｎ）が
図１９中の点９１であるとき、ノイズの影響で、求めら
れた相関値が実際には点９２となり、実数部Ｒｅ（ΣＳ
（ｎ））が負の値をとる場合がある。しかし、ΔθEST
の取りうる範囲は−π／２≦ΔθEST≦π／２であるた
め、Ｉｍ（ΣＳ（ｎ９）をＲｅ（ΣＳ（ｎ））で除算し
た結果、周波数オフセット推定に用いられる相関値は点
９３となる。このとき、周波数オフセット推定値は、真
の周波数オフセットと比較して符号自身が反転したもの
となり、推定誤差は非常に大きくなる。
図１９を用いて説明する。例えば、相関値ΣＳ（ｎ）が
図１９中の点９１であるとき、ノイズの影響で、求めら
れた相関値が実際には点９２となり、実数部Ｒｅ（ΣＳ
（ｎ））が負の値をとる場合がある。しかし、ΔθEST
の取りうる範囲は−π／２≦ΔθEST≦π／２であるた
め、Ｉｍ（ΣＳ（ｎ９）をＲｅ（ΣＳ（ｎ））で除算し
た結果、周波数オフセット推定に用いられる相関値は点
９３となる。このとき、周波数オフセット推定値は、真
の周波数オフセットと比較して符号自身が反転したもの
となり、推定誤差は非常に大きくなる。
【００５１】この影響を低減するため、除算部５４８
を、本実施の形態である図１８のような構成にし、絶対
値 |Ｒｅ（ΣＳ(n))|でＩｍ（ΣＳ（ｎ））を除算し
て、誤差信号Ｅ（ｎ）を得ることにより、Ｒｅ（ΣＳ
（ｎ））の符号反転によるＥ（ｎ）の符号反転、すなわ
ちΔｆESTの符号反転を防止することができ、推定精度
の劣化を抑えることができる。この構成を、他の除算部
５４１、５４２、５４３にも適用することができる。以
上が本発明に係る第五の実施例である。
を、本実施の形態である図１８のような構成にし、絶対
値 |Ｒｅ（ΣＳ(n))|でＩｍ（ΣＳ（ｎ））を除算し
て、誤差信号Ｅ（ｎ）を得ることにより、Ｒｅ（ΣＳ
（ｎ））の符号反転によるＥ（ｎ）の符号反転、すなわ
ちΔｆESTの符号反転を防止することができ、推定精度
の劣化を抑えることができる。この構成を、他の除算部
５４１、５４２、５４３にも適用することができる。以
上が本発明に係る第五の実施例である。
【００５２】（実施例６）本発明に係る周波数同期方法
を用いたＯＦＤＭ受信装置の構成について説明する。図
２０はＯＦＤＭ信号受信装置のブロック図であり、アン
テナ６０から入力された受信信号の周波数変換を行うＲ
Ｆ部６１と、アナログ信号を標本量子化してディジタル
信号ｓ（ｎ）に変換するＡ／Ｄ変換器６２と、受信パケ
ットの開始タイミングを検出するタイミング同期部６３
と、ｓ（ｎ）から周波数オフセットを推定して推定値Δ
ｆESTを出力する周波数オフセット推定部６４と、Δｆ
ESTの周波数で振動する正弦波を受信信号に乗算して周
波数オフセットを補正する周波数オフセット補正部６５
と、ＦＦＴ演算を行うＦＦＴ演算部６６と、誤り訂正を
行うビタビデコーダ６７とを有する。
を用いたＯＦＤＭ受信装置の構成について説明する。図
２０はＯＦＤＭ信号受信装置のブロック図であり、アン
テナ６０から入力された受信信号の周波数変換を行うＲ
Ｆ部６１と、アナログ信号を標本量子化してディジタル
信号ｓ（ｎ）に変換するＡ／Ｄ変換器６２と、受信パケ
ットの開始タイミングを検出するタイミング同期部６３
と、ｓ（ｎ）から周波数オフセットを推定して推定値Δ
ｆESTを出力する周波数オフセット推定部６４と、Δｆ
ESTの周波数で振動する正弦波を受信信号に乗算して周
波数オフセットを補正する周波数オフセット補正部６５
と、ＦＦＴ演算を行うＦＦＴ演算部６６と、誤り訂正を
行うビタビデコーダ６７とを有する。
【００５３】周波数同期部６８は、周波数オフセット推
定部６４と、周波数オフセット補正部６５から構成され
る。本実施例では、周波数オフセット推定部６４に、上
記１乃至５のいずれかに記載した周波数オフセット推定
方法を適用する。
定部６４と、周波数オフセット補正部６５から構成され
る。本実施例では、周波数オフセット推定部６４に、上
記１乃至５のいずれかに記載した周波数オフセット推定
方法を適用する。
【００５４】図２０において、ＲＦ部６１に受信ＯＦＤ
Ｍ信号が入力される。ＲＦ部６１では、入力される受信
ＯＦＤＭ信号を搬送波にほぼ近いローカル信号によりア
ナログ複素ベースバンド信号に変換する。Ａ／Ｄ変換器
６２は、ＲＦ部６１にて得られたアナログ複素ベースバ
ンド信号を標本量子化して、ディジタル受信信号ｓ
（ｎ）を出力する。
Ｍ信号が入力される。ＲＦ部６１では、入力される受信
ＯＦＤＭ信号を搬送波にほぼ近いローカル信号によりア
ナログ複素ベースバンド信号に変換する。Ａ／Ｄ変換器
６２は、ＲＦ部６１にて得られたアナログ複素ベースバ
ンド信号を標本量子化して、ディジタル受信信号ｓ
（ｎ）を出力する。
【００５５】ｓ（ｎ）はタイミング同期部63および周波
数オフセット推定部６４に入力する。タイミング同期部
６３では、受信信号におけるプリアンブルの開始位置を
検出し、周波数オフセット推定部６４およびＦＦＴ演算
部６６の動作開始タイミングを決定する。周波数オフセ
ット推定部64では、信号ｓ（ｎ）を入力し、タイミング
検出部６３にて示されたタイミングで周波数オフセット
を推定し、推定値Δｆ ESTを出力する。
数オフセット推定部６４に入力する。タイミング同期部
６３では、受信信号におけるプリアンブルの開始位置を
検出し、周波数オフセット推定部６４およびＦＦＴ演算
部６６の動作開始タイミングを決定する。周波数オフセ
ット推定部64では、信号ｓ（ｎ）を入力し、タイミング
検出部６３にて示されたタイミングで周波数オフセット
を推定し、推定値Δｆ ESTを出力する。
【００５６】さらに、周波数オフセット補正部６５で
は、周波数オフセット推定部６４で得られたΔｆESTを
入力し、これに対応する正弦波を発生して、受信信号ｓ
（ｎ）と複素乗算し、周波数オフセットを補正する。周
波数オフセット補正された信号に対して、ＦＦＴ演算部
６６でＦＦＴ演算を施し、ビタビデコーダ６７で誤り訂
正を行う。
は、周波数オフセット推定部６４で得られたΔｆESTを
入力し、これに対応する正弦波を発生して、受信信号ｓ
（ｎ）と複素乗算し、周波数オフセットを補正する。周
波数オフセット補正された信号に対して、ＦＦＴ演算部
６６でＦＦＴ演算を施し、ビタビデコーダ６７で誤り訂
正を行う。
【００５７】図２０のＯＦＤＭ受信装置は、本発明に係
る周波数同期方法を用いることで、従来より広い周波数
引き込み範囲の周波数同期、または、従来より小さい周
波数オフセットの推定誤差の周波数同期を実現でき、良
好なＯＦＤＭの復調信号を得ることができる。
る周波数同期方法を用いることで、従来より広い周波数
引き込み範囲の周波数同期、または、従来より小さい周
波数オフセットの推定誤差の周波数同期を実現でき、良
好なＯＦＤＭの復調信号を得ることができる。
【００５８】なお、上記１乃至５の実施例では遅延量を
６４サンプル、３２サンプル、１６サンプルの３種類と
しているが、これに限るものではない。すなわち、一般
的に、プリアンブル１０１の構成が、Ｌサンプル長の同
一の固定信号系列のＲ回の繰り返しの場合、遅延量は
Ｌ，２Ｌ…，（Ｒ／２）Ｌとすることができる。以上が
本発明に係る第六の実施例である。
６４サンプル、３２サンプル、１６サンプルの３種類と
しているが、これに限るものではない。すなわち、一般
的に、プリアンブル１０１の構成が、Ｌサンプル長の同
一の固定信号系列のＲ回の繰り返しの場合、遅延量は
Ｌ，２Ｌ…，（Ｒ／２）Ｌとすることができる。以上が
本発明に係る第六の実施例である。
【００５９】
【発明の効果】本発明によれば、大まかな周波数オフセ
ットの大きさを推定し、周波数オフセットが十分小さい
と判定されれば、より同期引き込み周波数範囲が狭い
（推定精度の良い）誤差信号を用いて周波数オフセット
を推定する。これにより、周波数オフセットが小さい場
合には高精度の周波数同期を確立でき、また、周波数オ
フセットが大きい場合の周波数同期、即ち広い周波数引
き込み範囲を実現できる。
ットの大きさを推定し、周波数オフセットが十分小さい
と判定されれば、より同期引き込み周波数範囲が狭い
（推定精度の良い）誤差信号を用いて周波数オフセット
を推定する。これにより、周波数オフセットが小さい場
合には高精度の周波数同期を確立でき、また、周波数オ
フセットが大きい場合の周波数同期、即ち広い周波数引
き込み範囲を実現できる。
【００６０】さらに本発明によれば、相関値を複数求め
てその平均値を用いて周波数オフセット推定の誤差信号
を計算することにより、周波数オフセット推定における
ノイズ成分の影響を抑制することができ、周波数オフセ
ット推定精度が向上し、高精度の周波数同期が確立でき
るようになる。
てその平均値を用いて周波数オフセット推定の誤差信号
を計算することにより、周波数オフセット推定における
ノイズ成分の影響を抑制することができ、周波数オフセ
ット推定精度が向上し、高精度の周波数同期が確立でき
るようになる。
【００６１】さらに本発明によれば、相関値の実数部の
絶対値をもって除算を行うことにより、伝送路ノイズに
よる誤差信号の符号の反転を防ぐことができることか
ら、周波数オフセットが同期引き込み周波数範囲の限界
近くの値を持つ場合でも、推定精度の劣化を抑制するこ
とが可能となる。
絶対値をもって除算を行うことにより、伝送路ノイズに
よる誤差信号の符号の反転を防ぐことができることか
ら、周波数オフセットが同期引き込み周波数範囲の限界
近くの値を持つ場合でも、推定精度の劣化を抑制するこ
とが可能となる。
【００６２】さらに本発明によれば、相関値からΔθ
ESTの属する象限も検出することにより、周波数オフセ
ット推定引き込み範囲を拡大することが可能となる。
ESTの属する象限も検出することにより、周波数オフセ
ット推定引き込み範囲を拡大することが可能となる。
【図１】本発明に係る周波数オフセット推定部のブロッ
ク図である。
ク図である。
【図２】本発明に係るプリアンブルの構成図である。
【図３】従来の周波数同期回路の構成を示すブロック図
である。
である。
【図４】従来の周波数同期回路における誤差信号生成部
のブロック図である。
のブロック図である。
【図５】ＯＦＤＭシステムにおいて伝送されるパケット
の構成図である。
の構成図である。
【図６】本発明に係る周波数オフセット推定部中の相関
値計算部のブロック図である。
値計算部のブロック図である。
【図７】本発明に係る周波数オフセット推定部のブロッ
ク図である。
ク図である。
【図８】本発明に係る誤差信号選択アルゴリズムを示す
フローチャート図である。
フローチャート図である。
【図９】本発明に係る周波数オフセット推定部中の相関
値計算部のブロック図である。
値計算部のブロック図である。
【図１０】本発明に係る周波数オフセット推定部中の判
定・演算部のブロック図である。
定・演算部のブロック図である。
【図１１】本発明に係る周波数オフセット推定部中の相
関値選択部のブロック図である。
関値選択部のブロック図である。
【図１２】本発明に係る相関値選択アルゴリズムを示す
フローチャート図である。
フローチャート図である。
【図１３】本発明に係る周波数オフセット推定部中の判
定・演算部のブロック図である。
定・演算部のブロック図である。
【図１４】本発明に係る相関値選択アルゴリズムを示す
フローチャート図である。
フローチャート図である。
【図１５】本発明に係る相関値の属する象限決定方法を
示す説明図である。
示す説明図である。
【図１６】本発明に係る相関値の属する象限決定方法を
示す説明図である。
示す説明図である。
【図１７】本発明に係る相関値の属する象限決定方法を
示す説明図である。
示す説明図である。
【図１８】本発明に係る周波数オフセット推定回路中の
除算部のブロック図である。
除算部のブロック図である。
【図１９】相関値計算における伝送路ノイズの影響を示
す説明図である。
す説明図である。
【図２０】本発明に係るＯＦＤＭ信号受信装置の一実施
例のブロック図である。
例のブロック図である。
【図２１】本発明に係る周波数オフセット推定部中の相
関値計算の説明図である。
関値計算の説明図である。
５１、５２、５３……相関値計算部、５４……判定・演
算部、２０、２１、２２、２３、２４、２５、２６、２
７、２８、２９……プリアンブルにおける固定信号系
列、３０……誤差信号生成部、３１……誤差信号変換
部、３２……正弦波発生部、３３……複素乗算部、３０
１……遅延部、３０２……複素共役部、３０３……複素
乗算部、３０４……積算部、３０５……除算部、１０…
…ＯＦＤＭシステムにおいて伝送されるパケット、１０
１……ＯＦＤＭシステムにおいて伝送されるパケットの
プリアンブル部、１０２……ＯＦＤＭシステムにおいて
伝送されるパケットのデータ部、５１１、５２１、５３
１……遅延部、５１２、５２２、５３２……複素共役
部、５１３、５２３、５３３……複素乗算部、５１４、
５２４、５３４……積算部、５４１、５４２、５４３…
…除算部、５４４……誤差信号選択部、５４５……アー
クタンジェント演算部、５４６……演算部、８１……Ｅ
（ｎ）１６に対する判定分岐処理、８２……Ｅ（ｎ）３
２に対する判定分岐処理、８３、８４、８５……誤差信
号判定結果、５２５、５２６……遅延部、５２７……複
素共役部、５２８……複素乗算部、５２９……積算部、
５２０……平均化部、５４７……相関値選択部、５４８
……除算部、８０１、８０２……乗算部、８０３……判
定部、８０４…… ΣＳ（ｎ）１６に対する判定分岐処
理、８０５……ΣＳ（ｎ）３２に対する判定分岐処理、
８０６、８０７、８０８……相関値選択結果、５４９…
…相関値選択部、５４０……位相周波数変換部、１１１
……ΣＳ（ｎ）１６に対する判定分岐処理、１１２……
ΣＳ（ｎ）３２に対する判定分岐処理、１１３、１１
４、１１５……相関値選択結果、９０１……絶対値算出
部、９０２……除算部、９１……理想状態における相関
値を表す点、９２……ノイズの影響を受けた状態で求め
られた相関値を表す点、９３……ノイズの影響を受けて
実数部が負数になった場合に推定に用いられる相関値を
表す点、６０……アンテナ、６１……ＲＦ部、６２……
Ａ／Ｄ変換器、６３……タイミング同期部、６４……周
波数オフセット推定部、６５……周波数オフセット補正
部、６６……ＦＦＴ演算部、６７……ビタビデコーダ、
６８……周波数同期部、７１、７２……相関値計算処理
ブロック、７３……平均化処理ブロック。
算部、２０、２１、２２、２３、２４、２５、２６、２
７、２８、２９……プリアンブルにおける固定信号系
列、３０……誤差信号生成部、３１……誤差信号変換
部、３２……正弦波発生部、３３……複素乗算部、３０
１……遅延部、３０２……複素共役部、３０３……複素
乗算部、３０４……積算部、３０５……除算部、１０…
…ＯＦＤＭシステムにおいて伝送されるパケット、１０
１……ＯＦＤＭシステムにおいて伝送されるパケットの
プリアンブル部、１０２……ＯＦＤＭシステムにおいて
伝送されるパケットのデータ部、５１１、５２１、５３
１……遅延部、５１２、５２２、５３２……複素共役
部、５１３、５２３、５３３……複素乗算部、５１４、
５２４、５３４……積算部、５４１、５４２、５４３…
…除算部、５４４……誤差信号選択部、５４５……アー
クタンジェント演算部、５４６……演算部、８１……Ｅ
（ｎ）１６に対する判定分岐処理、８２……Ｅ（ｎ）３
２に対する判定分岐処理、８３、８４、８５……誤差信
号判定結果、５２５、５２６……遅延部、５２７……複
素共役部、５２８……複素乗算部、５２９……積算部、
５２０……平均化部、５４７……相関値選択部、５４８
……除算部、８０１、８０２……乗算部、８０３……判
定部、８０４…… ΣＳ（ｎ）１６に対する判定分岐処
理、８０５……ΣＳ（ｎ）３２に対する判定分岐処理、
８０６、８０７、８０８……相関値選択結果、５４９…
…相関値選択部、５４０……位相周波数変換部、１１１
……ΣＳ（ｎ）１６に対する判定分岐処理、１１２……
ΣＳ（ｎ）３２に対する判定分岐処理、１１３、１１
４、１１５……相関値選択結果、９０１……絶対値算出
部、９０２……除算部、９１……理想状態における相関
値を表す点、９２……ノイズの影響を受けた状態で求め
られた相関値を表す点、９３……ノイズの影響を受けて
実数部が負数になった場合に推定に用いられる相関値を
表す点、６０……アンテナ、６１……ＲＦ部、６２……
Ａ／Ｄ変換器、６３……タイミング同期部、６４……周
波数オフセット推定部、６５……周波数オフセット補正
部、６６……ＦＦＴ演算部、６７……ビタビデコーダ、
６８……周波数同期部、７１、７２……相関値計算処理
ブロック、７３……平均化処理ブロック。
─────────────────────────────────────────────────────
フロントページの続き
(72)発明者 大久保 隆志
東京都青梅市新町六丁目16番地の３ 株式
会社日立製作所デバイス開発センタ内
Ｆターム(参考） 5K022 DD01 DD18 DD33 DD42 DD43
5K047 AA03 CC01 HH15 HH53 MM13
Claims (5)
- 【請求項１】 送信側では固定信号系列を複数回繰り返
したプリアンブルを送信し、受信側では周波数同期確立
のため前記送信側で送信されたプリアンブルの受信信号
に基づいて周波数オフセットを補正するＯＦＤＭ伝送シ
ステムに適用され、 当該システムの受信側で受信信号とその遅延信号の相関
値を求め、得られた相関値の位相成分を誤差信号とし、
この誤差信号から該周波数オフセットを推定する周波数
同期方法において、 該プリアンブルの受信信号とその遅延信号の間の遅延量
の異なる複数の相関値を求め、これら複数の相関値それ
ぞれに基づいて複数の誤差信号を計算し、該複数の相関
値もしくは該複数の誤差信号に基づいて周波数オフセッ
ト推定に用いる誤差信号を選択し、周波数オフセットを
推定することを特徴とする周波数同期方法。 - 【請求項２】 請求項１において、 前記複数の相関値のうち少なくとも一つの相関値につい
て、複素平面上のいずれの象限にあるかを検出すること
で、前記複数の相関値の中から周波数オフセット推定値
を求めるための相関値を選択し、この選択された相関値
と、該検出された象限の値に基づいて周波数オフセット
推定値を求めることを特徴とする周波数同期方法。 - 【請求項３】 請求項１又は２において、 前記誤差信号のうち少なくとも一つは、前記相関値の虚
数部をその実数部の絶対値で除算した値であることを特
徴とする周波数同期方法。 - 【請求項４】 請求項１乃至３のいずれかにおいて、 誤差信号のうち少なくとも一つは、受信信号とその遅延
信号の間の遅延量が等しい複数の相関値を計算し、これ
らの平均値に基づいて計算されることを特徴とする周波
数同期方法。 - 【請求項５】 請求項１乃至４の周波数同期方法を用い
て周波数オフセットを補正する回路を備えてなることを
特徴とするＯＦＤＭ受信装置。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2001345466A JP3945623B2 (ja) | 2001-11-12 | 2001-11-12 | 周波数同期方法及びこれを用いたｏｆｄｍ受信装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2001345466A JP3945623B2 (ja) | 2001-11-12 | 2001-11-12 | 周波数同期方法及びこれを用いたｏｆｄｍ受信装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2003152674A true JP2003152674A (ja) | 2003-05-23 |
| JP3945623B2 JP3945623B2 (ja) | 2007-07-18 |
Family
ID=19158788
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2001345466A Expired - Fee Related JP3945623B2 (ja) | 2001-11-12 | 2001-11-12 | 周波数同期方法及びこれを用いたｏｆｄｍ受信装置 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP3945623B2 (ja) |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| KR100601939B1 (ko) * | 2004-01-16 | 2006-07-14 | 삼성전자주식회사 | Ｏｆｄｍ 시스템에서의 초기 주파수 동기 방법 및 장치 |
| JP2007528643A (ja) * | 2004-01-28 | 2007-10-11 | クゥアルコム・インコーポレイテッド | 無線通信ネットワークにおける周波数取得のための装置及び方法 |
| JPWO2005101711A1 (ja) * | 2004-04-14 | 2008-03-06 | 松下電器産業株式会社 | 受信装置 |
| JP2011160034A (ja) * | 2010-01-29 | 2011-08-18 | Fujitsu Semiconductor Ltd | 信号処理装置，信号処理方法とそれを有する受信装置 |
| US8433005B2 (en) | 2004-01-28 | 2013-04-30 | Qualcomm Incorporated | Frame synchronization and initial symbol timing acquisition system and method |
| US8724447B2 (en) | 2004-01-28 | 2014-05-13 | Qualcomm Incorporated | Timing estimation in an OFDM receiver |
-
2001
- 2001-11-12 JP JP2001345466A patent/JP3945623B2/ja not_active Expired - Fee Related
Cited By (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| KR100601939B1 (ko) * | 2004-01-16 | 2006-07-14 | 삼성전자주식회사 | Ｏｆｄｍ 시스템에서의 초기 주파수 동기 방법 및 장치 |
| JP2007528643A (ja) * | 2004-01-28 | 2007-10-11 | クゥアルコム・インコーポレイテッド | 無線通信ネットワークにおける周波数取得のための装置及び方法 |
| JP2010220238A (ja) * | 2004-01-28 | 2010-09-30 | Qualcomm Inc | 無線通信ネットワークにおける周波数取得のための装置及び方法 |
| US8433005B2 (en) | 2004-01-28 | 2013-04-30 | Qualcomm Incorporated | Frame synchronization and initial symbol timing acquisition system and method |
| US8724447B2 (en) | 2004-01-28 | 2014-05-13 | Qualcomm Incorporated | Timing estimation in an OFDM receiver |
| JPWO2005101711A1 (ja) * | 2004-04-14 | 2008-03-06 | 松下電器産業株式会社 | 受信装置 |
| JP4583374B2 (ja) * | 2004-04-14 | 2010-11-17 | パナソニック株式会社 | 受信装置 |
| JP2011160034A (ja) * | 2010-01-29 | 2011-08-18 | Fujitsu Semiconductor Ltd | 信号処理装置，信号処理方法とそれを有する受信装置 |
Also Published As
| Publication number | Publication date |
|---|---|
| JP3945623B2 (ja) | 2007-07-18 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP3214159B2 (ja) | キャリア検出器 | |
| JP4125715B2 (ja) | Ｏｆｄｍシステムでの初期周波数の同期方法及び装置 | |
| US6628606B1 (en) | Coarse frequency offset estimator in orthogonal frequency division multiplexing receiver and method thereof | |
| US9712316B2 (en) | Reception apparatus, phase error estimation method, and phase error correction method | |
| JP4263119B2 (ja) | Ｏｆｄｍシステムでの初期周波数の同期方法および装置 | |
| US20060146962A1 (en) | Method and device for frame detection and synchronizer | |
| WO2000065756A1 (fr) | Recepteur de communication par paquets ofdm | |
| JP4279027B2 (ja) | Ｏｆｄｍ復調方法及び半導体集積回路 | |
| US20040202234A1 (en) | Low-complexity and fast frequency offset estimation for OFDM signals | |
| JP2005102169A (ja) | Ｏｆｄｍ受信装置およびｏｆｄｍ受信方法 | |
| US8503594B2 (en) | Phase tracking in communications systems | |
| JP3905541B2 (ja) | 遅延プロファイル推定装置及び相関器 | |
| JP3546804B2 (ja) | Ｏｆｄｍパケット通信用受信装置 | |
| JP3945623B2 (ja) | 周波数同期方法及びこれを用いたｏｆｄｍ受信装置 | |
| US7583770B2 (en) | Multiplex signal error correction method and device | |
| KR100458630B1 (ko) | 위상편이 및 주파수편이 결합추정방법 및 그 장치 | |
| JP2000059332A (ja) | 遅延プロファイル検出回路及び検出方法 | |
| JP2818148B2 (ja) | Ｏｆｄｍ復調装置 | |
| JPH09246917A (ja) | 周波数誤差推定装置 | |
| KR20020086161A (ko) | 직교 주파수 분할 다중 전송 시스템의 훈련 심볼 결정방법 및 주파수 옵셋 추정과 동기를 위한 장치 및 방법 | |
| JP2003188846A (ja) | 通信装置および周波数偏差推定方法 | |
| JP2000358011A (ja) | 受信装置 | |
| JP6037992B2 (ja) | 受信装置及び方法 | |
| JP3086144B2 (ja) | バースト復調器 | |
| JP2001313629A (ja) | Ｏｆｄｍ通信装置 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20041109 |
|
| A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20060703 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20060927 |
|
| A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20061121 |
|
| TRDD | Decision of grant or rejection written | ||
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20070404 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20070404 |
|
| R150 | Certificate of patent or registration of utility model |
Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100420 Year of fee payment: 3 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110420 Year of fee payment: 4 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110420 Year of fee payment: 4 |
|
| S111 | Request for change of ownership or part of ownership |
Free format text: JAPANESE INTERMEDIATE CODE: R313111 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110420 Year of fee payment: 4 |
|
| R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120420 Year of fee payment: 5 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120420 Year of fee payment: 5 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130420 Year of fee payment: 6 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20140420 Year of fee payment: 7 |
|
| S531 | Written request for registration of change of domicile |
Free format text: JAPANESE INTERMEDIATE CODE: R313531 |
|
| R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |
|
| LAPS | Cancellation because of no payment of annual fees |