JP2007028636A - サンプリング周波数誤差推定装置、サンプリング周波数誤差推定方法、直交周波数分割多重システムの受信機及びその位相歪み補正方法 - Google Patents

サンプリング周波数誤差推定装置、サンプリング周波数誤差推定方法、直交周波数分割多重システムの受信機及びその位相歪み補正方法 Download PDF

Info

Publication number
JP2007028636A
JP2007028636A JP2006197325A JP2006197325A JP2007028636A JP 2007028636 A JP2007028636 A JP 2007028636A JP 2006197325 A JP2006197325 A JP 2006197325A JP 2006197325 A JP2006197325 A JP 2006197325A JP 2007028636 A JP2007028636 A JP 2007028636A
Authority
JP
Japan
Prior art keywords
error
phase change
phase
pilot subcarrier
sampling
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
Application number
JP2006197325A
Other languages
English (en)
Other versions
JP4320026B2 (ja
Inventor
Yun-Young Kim
倫 永 金
Jae-Ho Roh
在 皓 盧
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Samsung Electronics Co Ltd
Original Assignee
Samsung Electronics Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Samsung Electronics Co Ltd filed Critical Samsung Electronics Co Ltd
Publication of JP2007028636A publication Critical patent/JP2007028636A/ja
Application granted granted Critical
Publication of JP4320026B2 publication Critical patent/JP4320026B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L27/00Modulated-carrier systems
    • H04L27/26Systems using multi-frequency codes
    • H04L27/2601Multicarrier modulation systems
    • H04L27/2647Arrangements specific to the receiver only
    • H04L27/2655Synchronisation arrangements
    • H04L27/2668Details of algorithms
    • H04L27/2673Details of algorithms characterised by synchronisation parameters
    • H04L27/2675Pilot or known symbols
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L25/00Baseband systems
    • H04L25/02Details ; arrangements for supplying electrical power along data transmission lines
    • H04L25/0202Channel estimation
    • H04L25/0224Channel estimation using sounding signals
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L27/00Modulated-carrier systems
    • H04L27/01Equalisers
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L27/00Modulated-carrier systems
    • H04L27/26Systems using multi-frequency codes
    • H04L27/2601Multicarrier modulation systems
    • H04L27/2647Arrangements specific to the receiver only
    • H04L27/2655Synchronisation arrangements
    • H04L27/2657Carrier synchronisation
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L27/00Modulated-carrier systems
    • H04L27/26Systems using multi-frequency codes
    • H04L27/2601Multicarrier modulation systems
    • H04L27/2647Arrangements specific to the receiver only
    • H04L27/2655Synchronisation arrangements
    • H04L27/2668Details of algorithms
    • H04L27/2673Details of algorithms characterised by synchronisation parameters
    • H04L27/2676Blind, i.e. without using known symbols
    • H04L27/2679Decision-aided

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Power Engineering (AREA)
  • Synchronisation In Digital Transmission Systems (AREA)

Abstract

【課題】パイロット副搬送波の平均位相変化率ではなく、パイロット副搬送波ごとに位相変化率を計算してサンプリング周波数の誤差を正確に推定できるOFDMシステムに適用されるサンプリング周波数誤差推定装置及びその方法を提供することにある。
【解決手段】本発明によるサンプリング周波数誤差推定装置は、サンプリングされ周波数領域に変換されたサンプル受信信号に含まれたパイロット副搬送波の位相変化を計算する位相計算部610と、前記位相変化によって前記パイロット副搬送波のサンプルがシフトされる位相変化到達時間を計算する判断部620と、前記パイロット副搬送波を利用してサンプリング時間の誤差を計算する時間誤差計算部と、前記位相変化到達時間及び前記サンプリング時間誤差に基づいてサンプリング周波数の誤差を計算する周波数誤差計算部と、を含む。
【選択図】図4

Description

本発明は、OFDMシステムでのサンプリング周波数誤差推定装置、サンプリング周波数誤差推定方法、直交周波数分割多重システムの受信機及びその位相歪み補正方法に関し、より詳細には受信信号に含まれたパイロット副搬送波ごとに位相変化を推定してサンプリング周波数の誤差を推定するOFDMシステムでのサンプリング周波数誤差推定装置、サンプリング周波数誤差推定方法、直交周波数分割多重システムの受信機及びその位相歪み補正方法に関する。
直交周波数分割多重化方式であるOFDM方式は、直列形態で入力されるデータ列を所定のブロック単位の並列データに変換した後、並列化されたシンボルを互いに直交である相違する搬送波周波数に多重化(Multiplexing)し、広帯域伝送を多数の狭帯域並列伝送に変える方式である。こうしたOFDM方式は、無線通信環境下で多重経路フェーディング(fading)に強く、高速データ伝送が可能である。
OFDMシステムにおいて、サンプリング周波数の誤差追跡(tracking)はパイロット副搬送波(pilot subcarrier)がよく利用されている。このパイロット副搬送波を利用する方式では、受信機が同期化を実行するためにデータを使用できるように、パイロットシンボルデータ(これは送信機と受信機に認識されているデータである)が送られる。このパイロット副搬送波は2つの集合に分けられる。つまり、負の副搬送波に位置したパイロットと正の副搬送波に位置したパイロットである。サンプリング周波数の誤差は、誤差とパイロット副搬送波のインデックスによって発生する位相回転(phase rotation)との間の線形的関係に対する情報を使うことで推定される。
図1A及び図1Bは、OFDMシステムにおけるサンプリング周波数誤差の影響を説明するための図面である。図1Aは、0.5の初期時間の誤差のある30ppm(parts per million)の場合における、位相及び振幅歪みを示す図面である。そして、図1Bは、0.5の初期時間の誤差と1−サンプルシンボル時間の誤差のある30ppmの場合における、位相及び振幅歪みを示す図面である。なお、図1A及び図1Bにおいて、横軸は位相の大きさを表わし、縦軸は振幅歪みの大きさを表わす。
図2A及び図2Bは、従来のOFDMシステムにおけるサンプリング周波数の誤差推定を説明するための図面である。
図2Aは、パイロット副搬送波間の位相差を利用してサンプリング周波数の誤差を推定する場合を説明するための図面である。同図によると、パイロット搬送波間の位相差は次の数式によって求めることができる。
Figure 2007028636
ここで、Φm,pは、m番目のOFDMシンボルでのp番目のパイロット副搬送波を、Φm,p+1は、p+1番目のパイロット副搬送波の位相を示す。そして、Δkは、パイロット副搬送波間の間隔、Δtは、サンプリング時間誤差、Nは、FFTウィンドウサイズ(Fast Fourier Transform Window Size)を示す。
[数1]をサンプリング時間誤差に関して示すと次の通りである。
Figure 2007028636
[数2]に示すように、サンプリング時間の誤差を計算して位相の傾きを推定することで、サンプリング周波数の誤差を推定することができる。
図2Bは、時間に応じたパイロット副搬送波間の位相差の変化推定を利用してサンプリング周波数の誤差を推定する場合を説明するための図面である。
同図によると、パイロット搬送波間の位相差の変化は次の数式によって求めることができる。
Figure 2007028636
ここで、εp,p+1は、l番目とl+1番目のOFDMシンボルでp番目とp+1番目のパイロット副搬送波間の位相差を利用して計算された位相の傾きの変化量を示す。ここで、Δkは、パイロット副搬送波間の間隔、Δtはサンプリング時間の誤差を示す。
このとき、l+1番目のOFDMシンボルの代わりにDだけ遅延されたシンボルとl番目のシンボルでp番目、p+1番目のパイロットとの間の位相差を利用して位相の傾きの変化量を推定することもできる。
図3は、従来のサンプリング周波数の誤差推定方法によるシンボルインデックスと推定した位相との関係を示す図面である。ここで、「S」は1−サンプルがシフトされる地点で、サンプル一つを追加もしくは削除してサンプリング周波数の誤差によって発生した歪みを補償する必要がある地点である。
図3に示すように、「L」領域では時間(OFDM SYMBOL INDEX)が増加するにつれサンプリング時間の誤差の位相(ESTIMATED PHASE [RADIAN])が線形的に増加している。しかし、「L」領域以後では時間が経過するにつれ位相が−πもしくはπ以上を持つパイロット副搬送波が存在し(図3で示される「E」を参照)、パイロット副搬送波間の位相の線形性が成り立たない。従って、パイロット副搬送波のインデックスによる位相の平均変化量を利用する場合、誤ったサンプリング周波数の誤差を推定することになる。
そして、「S」は多数のパイロット副搬送波によって計算される平均位相傾きが1つのサンプルによって異なっており、雑音の影響によりサンプルを追加もしくは削除しなければならないときの点を示している。そのため、標準偏差が非常に大きくなっている。つまり、パイロット副搬送波の平均位相傾きを利用する場合、1−サンプルを追加もしくは削除する点を正確に検出することができない。
韓国特許出願公開第2004−014571号明細書 特開平11−331120号公報 米国特許出願公開第2003−0210645号明細書
本発明の目的は、パイロット副搬送波の平均位相変化率ではなく、パイロット副搬送波ごとに位相変化率を計算してサンプリング周波数の誤差を正確に推定できるOFDMシステムに適用されるサンプリング周波数誤差推定装置及びその方法を提供することにある。
前記目的を達成するための本発明に係るサンプリング周波数誤差推定装置は、サンプリングされ周波数領域に変換されたサンプル受信信号に含まれるパイロット副搬送波ごとの位相変化を計算する位相計算部と、前記位相変化によって前記パイロット副搬送波のサンプルがシフトする位相変化到達時間を計算する判断部と、前記パイロット副搬送波に基づいてサンプリング時間の誤差を計算する時間誤差計算部と、前記位相変化到達時間及び前記サンプリング時間誤差に基づいてサンプリング周波数の誤差を計算する周波数誤差計算部と、を含むことを特徴としている。
ここで、前記周波数領域に変換された前記サンプル受信信号のうち、前記パイロット副搬送波を保存する保存部を更に含むのが好ましい。
また、前記判断部は、前記パイロット副搬送波の位相がπ以上もしくは−π以下となるシンボルインデックスを検出して前記位相変化到達時間を計算するのが好ましい。
一方、前記目的を達成するための本発明に係る直交周波数分割多重システムの受信機は、予め設定されたサンプリング周波数に基づいて受信信号をサンプリングするアナログ/デジタル(AD)変換部と、前記サンプリングされた受信信号を周波数領域に変換する高速フーリエ変換(FET)部と、前記周波数領域に変換された受信信号に含まれるパイロット副搬送波ごとの位相変化を計算する位相計算部と、前記位相変化によって前記パイロット副搬送波のサンプルがシフトする位相変化到達時間を計算する判断部と、前記パイロット副搬送波に基づいてサンプリング時間の誤差を計算する時間誤差計算部と、前記位相変化到達時間及び前記サンプリング時間の誤差に基づいてサンプリング周波数の誤差を計算する周波数誤差計算部と、を含むサンプリング周波数誤差推定装置を備えることを特徴としている。
ここで、前記判断部は、前記パイロット副搬送波の位相がπ以上もしくは−π以下となるシンボルインデックスを検出して前記位相変化到達時間を計算するのが好ましい。
また、計算した前記サンプリング時間の誤差に基づいて前記サンプリングされた受信信号の位相歪みを補償する等化部を更に含むのが好ましい。
そして、前記サンプリング周波数の誤差に基づいて1−サンプルを削除もしくは追加する削除/追加部を更に含むのが好ましい。
一方、前記目的を達成するための本発明に係るサンプリング周波数誤差推定方法は、位相計算部により、サンプリングされ周波数領域に変換された受信信号に含まれたパイロット副搬送波ごとの位相変化を計算する段階と、周波数誤差計算部により、前記位相変化によってそれぞれの前記パイロット副搬送波のサンプルがシフトする位相変化到達時間に基づいてサンプリング周波数の誤差を計算する段階と、を含むことを特徴としている。
ここで、前記サンプリング周波数の誤差を計算する段階は、前記位相計算部により、サンプリングされ周波数領域に変換された前記受信信号に含まれたパイロット副搬送波ごとの位相変化を計算する段階と、判断部により、前記位相変化によって前記パイロット副搬送波のサンプルがシフトされる位相変化到達時間を計算する段階と、時間誤差計算部により、前記パイロット副搬送波を利用してサンプリング時間の誤差を計算する段階と、周波数誤差計算部により、前記位相変化到達時間及び前記サンプリング時間誤差に基づいて前記サンプリング周波数の誤差を計算する段階と、を含むのが好ましい。
また、保存部により、サンプリングされ周波数領域に変換された前記受信信号のうち、前記パイロット副搬送波を保存する段階を更に含むのが好ましい。
そして、前記位相変化到達時間を計算する段階は、前記判断部により、前記パイロット副搬送波の位相がπ以上もしくは−π以下となるシンボルインデックスを検出して前記位相変化到達時間を計算するのが好ましい。
一方、前記目的を達成するための本発明による直交周波数分割多重システムの受信機における位相歪み補正方法は、アナログ/デジタル(AD)変換部により、予め設定されたサンプリング周波数に基づいて受信信号をサンプリングする段階と、高速フーリエ変換(FET)部により、サンプリングされた受信信号を周波数領域に変換する段階と、位相計算部により、サンプリングされ周波数領域に変換された前記受信信号に含まれたパイロット副搬送波ごとの位相変化を計算する段階と、周波数誤差計算部により、前記位相変化によって前記パイロット副搬送波のサンプルがシフトする位相変化到達時間に基づいてサンプリング周波数の誤差を計算する段階と、を含むことを特徴としている。
ここで、前記サンプリング周波数の誤差を計算する段階は、前記位相計算部により、サンプリングされ周波数領域に変換された前記受信信号に含まれたパイロット副搬送波ごとの位相変化を計算する段階と、判断部により、前記位相変化によって前記パイロット副搬送波のサンプルがシフトする位相変化到達時間を計算する段階と、時間誤差計算部により、前記パイロット副搬送波を利用してサンプリング時間の誤差を計算する段階と、前記周波数誤差計算部により、前記位相変化到達時間及び前記サンプリング時間誤差に基づいて前記サンプリング周波数の誤差を計算する段階と、を含むのが好ましい。
また、保存部により、サンプリングされ周波数領域に変換された前記受信信号の前記パイロット副搬送波を保存する段階を更に含むのが好ましい。
また、前記位相変化到達時間は、前記判断部により、前記サンプリングされ周波数領域に変換された前記受信信号の前記パイロット副搬送波の位相がπ以上もしくは−π以下となるシンボルインデックスを検出して計算されるのが好ましい。
また、等化部により、前記サンプリング周波数の誤差に基づいて前記サンプリングされた受信信号の位相歪みを補正する段階を更に含む。
また、削除/追加部により、前記サンプリング周波数の誤差に基づいて1−サンプルを削除もしくは追加する段階を更に含むのが好ましい。
本発明によると、パイロット副搬送波ごとの位相変化を計算することで、正確なサンプリング周波数の誤差が推定できるだけではなく、サンプリング周波数の誤差による位相歪みを正確に補償することができる。
以下、添付の図面に基づいて本発明の好適な実施形態を詳述する。
なお、各図面の構成要素に参照符号を付加することにおいて、同一構成要素に対し、たとえ他の図面上に表示されるとしても、できる限り同一符号を付加することとした。
また、本発明を説明することにおいて、関連する公知構成或は機能に対する詳細な説明が本発明の要旨を不明瞭にすると判断される場合は、そのような詳細な説明を省略しすることとした。
OFDMシステムにおいて、サンプリング周波数誤差による回転は2つの方式によって補償(Correction)することができる。一つは、受信機ADC(Analog−to−digital Converter)のサンプリング周波数を変えて修正する方式の同期化サンプリング方式であり、もう一つは、副搬送波をまた回転させることにより、DFT(Discrete Fourier Transform)処理以後に回転現象を補償する方式の非同期化サンプリング方式である。
以下では、一定のサンプリング周波数に基づいて受信信号をサンプリングした後、DFT処理以後に回転現象を補償する非同期化サンプリング方式に適用されるサンプリング周波数の誤差を推定するサンプリング周波数誤差推定装置について説明する。
図4は、本発明の一実施形態に係るサンプリング周波数誤差推定装置が適用されるOFDMシステムのブロック図である。
同図に示すように、本発明のサンプリング周波数誤差推定装置が適用されるOFDMシステムは、AD変換部100、削除/追加部(rob/stuff unit)200、FFT部(Fast Fourier Transform Unit)300、等化部(Equalizer)400、保存部500、及び周波数誤差推定部(Frequency offset estimator)600を含む。
AD変換部100は、受信される信号を発信器(Oscillator)(図示せず)から生成された一定のサンプリング周波数に応じてサンプリングする。
FFT部300は、サンプリングした前記受信信号を周波数領域に変換する。
保存部500は、周波数領域に変換した受信信号のうちパイロット副搬送波(pilot subcarrier)を保存する。OFDM受信機は送信機および受信機が認識しているデータであるパイロットシンボル(pilot symbol)が伝送され、同期化を行い、サンプリング周波数の誤差推定時に用いられるパイロット副搬送波が別途保存される。
周波数誤差推定部600は、位相計算部610、判断部620、及び誤差計算部630を含み、保存されたパイロット副搬送波ごとの位相変化に基づいてサンプリング周波数誤差及びサンプリング時間誤差を推定(計算)する。
位相計算部610は、周波数領域に変換された受信信号に含まれた複数のパイロット副搬送波の位相変化をそれぞれ計算する。ここで、パイロット副搬送波は、サンプリング周波数の誤差に応じて互いに異なる位相変化率を有するため、パイロット副搬送波ごとに位相変化を計算する。
判断部620は、それぞれのパイロット副搬送波で急激な位相変化を発生する時間を計算する。パイロット副搬送波の位相がπ以上もしくは−π以下となる時間、つまり1−サンプルがシフトする時間である位相変化到達時間を計算する。サンプリング周波数の誤差は、パイロット副搬送波のインデックスによって発生する位相回転(phase rotation)と線形的関係にある。従って、一つのシンボルにおいて外側に位置する最も上のパイロット副搬送波の位相がπ以上もしくは−π以下となる時間に到達する。
誤差計算部630は、図示しない時間誤差計算部により、パイロット副搬送波のインデックスに基づいてサンプリング時間誤差を計算し、図示しない周波数誤差計算部により、計算されたサンプリング時間誤差及び推定した位相変化到達時間に基づいてサンプリング周波数誤差を計算する。
等化部400は、周波数誤差推定部600で推定されたサンプリング時間誤差を利用して位相歪みを補償する。
削除/追加部200は、周波数誤差推定部600で推定されたサンプリング周波数の誤差に基づいてコピーされたサンプルを追加もしくは信号から一つのサンプルを削除する。サンプリング周波数の誤差によってサンプリング瞬間の変化が次第にサンプリング周期より大きくなるため、1−サンプルシフトが発生する位置を検出しサンプルを追加もしくは削除してサンプリング周波数の誤差による歪みを補償する。
図5は、本発明の一実施形態に係るOFDMシステムに適用されるサンプリング周波数の誤差推定方法を説明するための図面である。
同図に示すように、まず受信信号のうち送信機および受信機が認識しているデータであるパイロット副搬送波の位相変化を計算する(S910)。受信信号は、一定のサンプリング周波数によってサンプリングされた後、FFT部300によって周波数領域に変換される。周波数領域に変換された受信信号のうち、パイロット副搬送波は、パイロット副搬送波のインデックスごとに位相変化が計算される。
続いて、判断部620は、それぞれのパイロット副搬送波ごとに位相がπ以上もしくは−π以下となるかを判断する(S920)。パイロット副搬送波は、サンプリング周波数の誤差に応じて互いに異なる位相変化率を有し、サンプリング周波数の誤差とパイロット副搬送波のインデックスによって発生する位相回転(phase rotation)とは線形的関係にある。
従って、位相がπ以上もしくは−π以下となる位置はパイロット副搬送波のインデックスごとに相違し、一つのシンボルにおいて外側に位置する最も上のパイロット副搬送波の位相がπ以上もしくは−π以下となる時間に到達する。
そして、S920段階での判断結果、パイロット副搬送波の位相がπ以上もしくは−π以下となる場合、位相がπ以上もしくは−π以下となる位相変化到達時間を計算する(S930)。位相変化到達時間は、位相がπ以上もしくは−π以下となる位置のシンボルインデックス、サンプリング時間区間、及びシンボル当たりのサンプル数を利用して計算する。このとき、位相がπ以上もしくは−π以下となる時間は、1−サンプルがシフトされる時間である。1−サンプルがシフトする位相変化到達時間を正確に検出することで、検出した位相変化到達時間に基づいてサンプルを一つ追加もしくは削除し、サンプリング周波数の誤差による歪みを補償することができる。
一方、S920段階での判断結果、パイロット副搬送波の位相がπ以上もしくは−π以下ではない場合、シンボルインデックスを増加させながらパイロット副搬送波の位相変化を計算する。
次に、サンプリング時間の誤差を計算する(S940)。サンプリング時間の誤差は、パイロット副搬送波のインデックスとFFTウィンドウサイズに基づいて求めることができる。
次に、S930段階で計算された位相変化到達時間と、S940段階で計算されたサンプリング時間誤差とに基づいてサンプリング周波数の誤差を計算する(S950。)サンプリング周波数の誤差は、サンプリング時間誤差を位相変化到達時間で割ること(割算)で求めることができる。
本発明では、このように、相違する位相変化率を持つパイロット副搬送波ごとに1−サンプルがシフトされる時間とサンプリング時間の誤差を計算し、サンプリング周波数の誤差を計算する。従って、従来のようにパイロット副搬送波の平均位相変化を推定する場合と異なり、独立のパイロット副搬送波の位相変化を計算するので、1−サンプルを追加もしくは削除する時点を正確に推定することができる。
このとき、多数のパイロット副搬送波のサンプリング周波数の誤差を計算し、サンプリング周波数の誤差による位相変化を補償することができるが、最初に1−サンプルがシフトされる最も上のパイロット副搬送波のサンプリング周波数の誤差のみ利用して位相変化を補償することができる。
図6は、図5のサンプリング周波数の誤差推定方法を詳説するための図面である。同図において、横軸はOFDMシンボルインデックスを示し、縦軸はパイロット副搬送波のインデックスごとに推定される位相を示す。
図6に示すように、OFDM信号の受信時における各パイロット副搬送波に応じて発生する回転角の量は、パイロット副搬送波のインデックスとOFDMシンボルインデックスに応じて異なる。この回転角は、外側のパイロット副搬送波部分において最も大きく、連続するOFDMシンボルで増加する。つまり、サンプリング誤差は一般的に非常に小さいが、シンボルインデックスが大きくなるにつれ回転角がだんだん大きくなり、正確な復調が不可能になることがある。
従って、#1乃至#6のパイロット副搬送波(図6において「PILOT ♯1・・・PILOT ♯6」)のうち、外側のパイロット副搬送波である#1のパイロット副搬送波の位相がπ以上もしくは−π以下となる位相変化到達時間に先に到達する。ここで、Aは、#1のパイロット副搬送波が位相変化到達時間に到達する場合のOFDMシンボルインデックスを示す。
#1のパイロット副搬送波が位相変化到達時間に到達する時のOFDMシンボルインデックスは「182」である。従って、S930段階で説明したように、パイロット副搬送波の位相がπ以上もしくは−π以下となる位置の「182」、シンボル当たりのサンプル数、サンプリング時間区間を利用してパイロット副搬送波の位相変化到達時間を計算すると次の通りである。
Figure 2007028636
ここで、Δt1はm#1のパイロット副搬送波の位相変化到達時間であり、「182」は、位相変化到達時間でのOFDMシンボルインデックスを示す。そして、Nは、シンボル当たりのサンプル数であり、Tは、サンプリング時間区間を示す。
#1のパイロット副搬送波の位相変化到達時間を計算する方法と同様の方法で、#2乃至#4のパイロット副搬送波の位相変化到達時間を次のように示すことができる。
Figure 2007028636
ここで、Δt、Δt、Δtは、それぞれ#2乃至#4のパイロット副搬送波の位相変化到達時間を示す。
そして、それぞれのパイロット副搬送波のサンプリング時間誤差は、S940段階で述べたように、パイロット副搬送波のインデックスとFFTウィンドウサイズを利用して計算すれば次の通りである。
Figure 2007028636
ここで、δ1は、#1のパイロット副搬送波のサンプリング時間誤差を示し、Nは、FFTウィンドウサイズを示す。そして、「55」は、#1のパイロット副搬送波のインデックスを示す。
#1のパイロット副搬送波のサンプリング時間誤差を計算する方法と同様の方法で、#2乃至#4のパイロット副搬送波のサンプリング時間誤差を次のように示すことができる。
Figure 2007028636
ここで、δ、δ、δは、それぞれ#2乃至#4のパイロット副搬送波のサンプリング時間の誤差を示し、「45」、「35」、「25」は、それぞれ#2乃至#4のパイロット副搬送波のインデックスを示す。
そして、サンプリング周波数の誤差は、S950段階で説明したように、パイロット副搬送波の位相変化到達時間とサンプリング時間誤差を利用して計算される。従って、サンプリング周波数の誤差は、次の数式により求めることができる。
Figure 2007028636
ここで、ΔfSは、サンプリング周波数の誤差、Δδは、サンプリング時間誤差を示し、Δtは、パイロット副搬送波の位相変化到達時間を示す。
数8に示すように、それぞれのパイロット副搬送波に基づいてサンプリング周波数の誤差を個別に計算し、一つのサンプルを追加もしくは削除してサンプリング周波数の誤差による位相変化を補償することができる。或いは、多数のパイロット副搬送波のうち、最初に位相変化到達時間に至る最も上のパイロット副搬送波のサンプリング周波数の誤差のみに基づいて位相変化を補償することができる。
図6において、シンボル当たりのサンプル数が165である場合、[数4]、[数6]、[数8]を利用して計算される最も上のパイロット副搬送波である#1のパイロット副搬送波のサンプリング周波数の誤差は−40ppm(parts per million)である。従って、#1のパイロット副搬送波で250サンプル当たり一つのサンプルを追加することにより、サンプリング周波数の誤差による位相歪みを補償することができる。
図7A及び図7Bは、本発明に係るサンプリング周波数の誤差推定方法及び従来のサンプリング周波数誤差推定方法の性能をシミュレーションした結果を示すグラフである。図7Aは、シグナル対ノイズ比(Signal to Noise Ratio;SNR[dB])の変化による誤差平均(Means Square Error;MSE)を示す図面であり、図7Bは、シグナル対ノイズ比の変化によるパケット誤差率(Packet Error Rate;PER)を示す図面である。
誤差が0ppm、−10ppm、−20ppm、−40ppmの場合、図7Aのような誤差平均を求めることができる。図7Aに示すように、本発明のサンプリング周波数の誤差推定方法を適用することによって、本発明では(propsed)、誤差平均が従来(conventional)に比べて減少させることができる。
また、誤差が−10ppm、−20ppm、−40ppmの場合、図7Bのようなパケット誤差率を求めることができる。図7Bに示すように、本発明のサンプリング周波数の誤差推定方法を適用することによって、本発明では(propsed)、パケット誤差率も従来(conventional)に比べて減少させることができる。
なお、このときの超広帯域無線(Ultra Wide Band;UWB)OFDMシステム下でのサンプリングレートは、528MHz、OFDMシンボルサイズは165に設定し、データ率は53.3Mbpsに設定する。そして、ペイロードサイズは1024byteに設定し、多重経路チャネルモデル(multipath Channel Model)はCM1に設定している。
以上、図面に基づいて本発明の好適な実施形態を図示および説明してきたが、本発明の保護範囲は、前述の実施形態に限定されるものではなく、特許請求の範囲に記載された発明とその均等物にまで及ぶものである。
OFDMシステムにおいてサンプリング周波数誤差の影響を説明するための図面である。 OFDMシステムにおいてサンプリング周波数誤差の影響を説明するための図面である。 従来のOFDMシステムにおいてサンプリング周波数の誤差推定を説明するための図面である。 従来のOFDMシステムにおいてサンプリング周波数の誤差推定を説明するための図面である。 従来のサンプリング周波数の誤差推定方法によるシンボルインデックスと推定した位相との関係を示す図面である。 本発明の一実施形態に係るサンプリング周波数の誤差推定装置が適用されるOFDMシステムのブロック図である。 本発明の一実施形態に係るOFDMシステムに適用されるサンプリング周波数の誤差推定方法を説明するための図面である。 図5のサンプリング周波数の誤差推定方法を詳説するための図面である。 本発明に係るサンプリング周波数の誤差推定方法及び従来のサンプリング周波数の誤差推定方法の性能をシミュレーションした結果を示すグラフである。 本発明に係るサンプリング周波数の誤差推定方法及び従来のサンプリング周波数の誤差推定方法の性能をシミュレーションした結果を示すグラフである。
符号の説明
100 AD変換部
200 削除/追加部
300 FFT部
400 等化部
500 保存部
600 推定部
610 位相計算部
620 判断部
630 誤差計算部

Claims (17)

  1. サンプリングされ周波数領域に変換されたサンプル受信信号に含まれるパイロット副搬送波ごとの位相変化を計算する位相計算部と、
    前記位相変化によって前記パイロット副搬送波のサンプルがシフトする位相変化到達時間を計算する判断部と、
    前記パイロット副搬送波に基づいてサンプリング時間の誤差を計算する時間誤差計算部と、
    前記位相変化到達時間及び前記サンプリング時間誤差に基づいてサンプリング周波数の誤差を計算する周波数誤差計算部と、
    を含むことを特徴とするサンプリング周波数誤差推定装置。
  2. 前記周波数領域に変換された前記サンプル受信信号のうち、前記パイロット副搬送波を保存する保存部を更に含むことを特徴とする請求項1に記載のサンプリング周波数誤差推定装置。
  3. 前記判断部は、前記パイロット副搬送波の位相がπ以上もしくは−π以下となるシンボルインデックスを検出して前記位相変化到達時間を計算することを特徴とする請求項1に記載のサンプリング周波数誤差推定装置。
  4. 予め設定されたサンプリング周波数に基づいて受信信号をサンプリングするアナログ/デジタル(AD)変換部と、
    前記サンプリングされた受信信号を周波数領域に変換する高速フーリエ変換(FET)部と、
    前記周波数領域に変換された受信信号に含まれるパイロット副搬送波ごとの位相変化を計算する位相計算部と、
    前記位相変化によって前記パイロット副搬送波のサンプルがシフトする位相変化到達時間を計算する判断部と、
    前記パイロット副搬送波に基づいてサンプリング時間の誤差を計算する時間誤差計算部と、
    前記位相変化到達時間及び前記サンプリング時間の誤差に基づいてサンプリング周波数の誤差を計算する周波数誤差計算部と、
    を含むサンプリング周波数誤差推定装置を備えることを特徴とする直交周波数分割多重システムの受信機。
  5. 前記判断部は、前記パイロット副搬送波の位相がπ以上もしくは−π以下となるシンボルインデックスを検出して前記位相変化到達時間を計算することを特徴とする請求項4に記載の直交周波数分割多重システムの受信機。
  6. 計算した前記サンプリング時間の誤差に基づいて前記サンプリングされた受信信号の位相歪みを補償する等化部を更に含むことを特徴とする請求項4に記載の直交周波数分割多重システムの受信機。
  7. 前記サンプリング周波数の誤差に基づいて1−サンプルを削除もしくは追加する削除/追加部を更に含むことを特徴とする請求項4に記載の直交周波数分割多重システムの受信機。
  8. 位相計算部により、サンプリングされ周波数領域に変換された受信信号に含まれたパイロット副搬送波ごとの位相変化を計算する段階と、
    周波数誤差計算部により、前記位相変化によってそれぞれの前記パイロット副搬送波のサンプルがシフトする位相変化到達時間に基づいてサンプリング周波数の誤差を計算する段階と、
    を含むことを特徴とするサンプリング周波数誤差推定方法。
  9. 前記サンプリング周波数の誤差を計算する段階は、
    前記位相計算部により、サンプリングされ周波数領域に変換された前記受信信号に含まれたパイロット副搬送波ごとの位相変化を計算する段階と、
    判断部により、前記位相変化によって前記パイロット副搬送波のサンプルがシフトされる位相変化到達時間を計算する段階と、
    時間誤差計算部により、前記パイロット副搬送波を利用してサンプリング時間の誤差を計算する段階と、
    周波数誤差計算部により、前記位相変化到達時間及び前記サンプリング時間誤差に基づいて前記サンプリング周波数の誤差を計算する段階と、
    を含むことを特徴とする請求項8に記載のサンプリング周波数誤差推定方法。
  10. 保存部により、サンプリングされ周波数領域に変換された前記受信信号のうち前記パイロット副搬送波を保存する段階を更に含むことを特徴とする請求項8に記載のサンプリング周波数誤差推定方法。
  11. 前記位相変化到達時間を計算する段階は、前記判断部により、前記パイロット副搬送波の位相がπ以上もしくは−π以下となるシンボルインデックスを検出して前記位相変化到達時間を計算することを特徴とする請求項8に記載のサンプリング周波数誤差推定方法。
  12. アナログ/デジタル(AD)変換部により、予め設定されたサンプリング周波数に基づいて受信信号をサンプリングする段階と、
    高速フーリエ変換(FET)部により、サンプリングされた受信信号を周波数領域に変換する段階と、
    位相計算部により、サンプリングされ周波数領域に変換された前記受信信号に含まれたパイロット副搬送波ごとの位相変化を計算する段階と、
    周波数誤差計算部により、前記位相変化によって前記パイロット副搬送波のサンプルがシフトする位相変化到達時間に基づいてサンプリング周波数の誤差を計算する段階と、
    を含むことを特徴とする直交周波数分割多重システムの受信機での位相歪み補正方法。
  13. 前記サンプリング周波数の誤差を計算する段階は、
    前記位相計算部により、サンプリングされ周波数領域に変換された前記受信信号に含まれたパイロット副搬送波ごとの位相変化を計算する段階と、
    判断部により、前記位相変化によって前記パイロット副搬送波のサンプルがシフトする位相変化到達時間を計算する段階と、
    時間誤差計算部により、前記パイロット副搬送波を利用してサンプリング時間の誤差を計算する段階と、
    前記周波数誤差計算部により、前記位相変化到達時間及び前記サンプリング時間誤差に基づいて前記サンプリング周波数の誤差を計算する段階と、
    を含むことを特徴とする請求項12に記載の直交周波数分割多重システムの受信機での位相歪み補正方法。
  14. 保存部により、サンプリングされ周波数領域に変換された前記受信信号の前記パイロット副搬送波を保存する段階を更に含むことを特徴とする請求項12に記載の直交周波数分割多重システムの受信機での位相歪み補正方法。
  15. 前記位相変化到達時間は、前記判断部により、前記サンプリングされ周波数領域に変換された前記受信信号の前記パイロット副搬送波の位相がπ以上もしくは−π以下となるシンボルインデックスを検出して計算されることを特徴とする請求項12に記載の直交周波数分割多重システムの受信機での位相歪み補正方法。
  16. 等化部により、前記サンプリング周波数の誤差に基づいて前記サンプリングされた受信信号の位相歪みを補正する段階を更に含むことを特徴とする請求項12に記載の直交周波数分割多重システムの受信機での位相歪み補正方法。
  17. 削除/追加部により、前記サンプリング周波数の誤差に基づいて1−サンプルを削除もしくは追加する段階を更に含むことを特徴とする請求項12に記載の直交周波数分割多重システムの受信機での位相歪み補正方法。
JP2006197325A 2005-07-19 2006-07-19 サンプリング周波数誤差推定装置、サンプリング周波数誤差推定方法、直交周波数分割多重システムの受信機及びその位相歪み補正方法 Active JP4320026B2 (ja)

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
KR1020050065526A KR100719112B1 (ko) 2005-07-19 2005-07-19 Ofdm 시스템에 적용되는 샘플링 주파수 오차 추정장치및 그 방법

Publications (2)

Publication Number Publication Date
JP2007028636A true JP2007028636A (ja) 2007-02-01
JP4320026B2 JP4320026B2 (ja) 2009-08-26

Family

ID=37054798

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2006197325A Active JP4320026B2 (ja) 2005-07-19 2006-07-19 サンプリング周波数誤差推定装置、サンプリング周波数誤差推定方法、直交周波数分割多重システムの受信機及びその位相歪み補正方法

Country Status (5)

Country Link
US (1) US7702024B2 (ja)
EP (1) EP1746793A2 (ja)
JP (1) JP4320026B2 (ja)
KR (1) KR100719112B1 (ja)
CN (1) CN100591064C (ja)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2015089125A (ja) * 2013-10-31 2015-05-07 富士通株式会社 マルチキャリア変調信号の定時同期化装置及び方法
JP2017517939A (ja) * 2014-04-24 2017-06-29 富士通株式会社 チャネル等化追跡装置、方法及び受信機

Families Citing this family (23)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8363757B1 (en) * 1999-10-12 2013-01-29 Qualcomm Incorporated Method and apparatus for eliminating the effects of frequency offsets in a digital communication system
CN101136881B (zh) * 2006-08-31 2012-10-10 世意法(北京)半导体研发有限责任公司 一种ofdma上行链路接收机和估计载波频偏的方法
US20080117953A1 (en) * 2006-11-20 2008-05-22 Yu-Min Chuang Frequency offset estimation circuit of receiving device of wireless transmission system
KR100874919B1 (ko) * 2007-03-12 2008-12-19 삼성전자주식회사 샘플링 오프셋 추정 장치 및 샘플링 오프셋 추정 방법
CN101277284B (zh) * 2007-03-29 2012-12-26 深圳赛意法微电子有限公司 在drm中的采样时钟频率偏移量估计的方法及集成芯片
US8045628B2 (en) * 2007-10-18 2011-10-25 Nokia Corporation Digital video broadcast service discovery
CN101551790B (zh) * 2008-04-03 2014-11-05 中兴通讯股份有限公司 应用于通信领域的快速傅立叶变换实现方法及装置
CN101267423B (zh) * 2008-05-09 2011-07-20 哈尔滨工业大学 在超宽带系统中利用分数傅立叶域估计时钟频偏的方法
CN101789807B (zh) * 2008-07-18 2014-04-16 俊茂微电子(上海)有限公司 频偏估计方法、装置以及通信设备
KR101018140B1 (ko) * 2009-02-02 2011-02-25 성균관대학교산학협력단 Ofdm 시스템의 샘플링 주파수 동기 장치
WO2010099808A1 (en) 2009-03-06 2010-09-10 Telefonaktiebolaget Lm Ericsson (Publ) A method of and equipment for compensating carrier frequency offset in an orthogonal frequency division multiplexing wireless radio transmission system
WO2011051029A1 (en) * 2009-10-29 2011-05-05 Sony Corporation Receiving apparatus for high frequency imaging system
US8660013B2 (en) 2010-04-12 2014-02-25 Qualcomm Incorporated Detecting delimiters for low-overhead communication in a network
CN102790737B (zh) * 2011-05-17 2017-11-28 中兴通讯股份有限公司 一种系统的同步方法及装置
WO2013028553A1 (en) * 2011-08-19 2013-02-28 Lumedyne Technologies Incorporated Time domain switched analog-to digital converter apparatus and methods
CN102333370B (zh) * 2011-08-31 2015-11-18 华南理工大学 Ofdm系统中基于抽除/补偿机制的抽样时钟同步方法
CN103095638B (zh) * 2012-12-19 2016-06-08 西安电子科技大学 一种多径衰落信道下ofdm系统的采样频率偏移盲估算方法
CN103905361A (zh) * 2012-12-25 2014-07-02 普天信息技术研究院有限公司 一种ofdm系统中的采样频率同步方法
CN103281282B (zh) * 2013-06-13 2016-04-06 湖南国科微电子股份有限公司 一种dvb-t2系统中频偏估计的方法
US9160599B2 (en) * 2013-08-12 2015-10-13 Samsung Electronics Co., Ltd Method and apparatus for channel smoothing and estimation in OFDM system
CN105610755B (zh) * 2015-12-21 2020-05-26 航天恒星科技有限公司 一种突发信号的频偏估计方法及装置
CA3070694C (en) * 2018-01-22 2022-05-24 Radius Co., Ltd. Receiver method and transmission method
CN110031675A (zh) * 2019-04-19 2019-07-19 南京大学 一种数采系统实际采样频率的测量方法

Family Cites Families (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH11205275A (ja) 1998-07-09 1999-07-30 Toshiba Corp Ofdm受信装置
KR100313860B1 (ko) * 1998-07-30 2001-12-12 윤종용 Ofdm 전송 방식에서 미세 주파수 복원장치 및 그 방법
JP2000151545A (ja) 1998-11-06 2000-05-30 Sony Corp 復調装置および方法、並びに提供媒体
JP3851017B2 (ja) 1999-04-06 2006-11-29 日本放送協会 Ofdm復調装置
KR100488802B1 (ko) * 2002-12-09 2005-05-12 한국전자통신연구원 직교 주파수 분할 다중화 무선 통신 시스템에서의 반송파주파수 오차와 샘플링 주파수 오차 추적 방법 및 그 장치
TWI252656B (en) * 2003-03-21 2006-04-01 Realtek Semiconductor Corp Sampling clock compensation device of multi-carrier system and method thereof
US7333548B2 (en) * 2003-11-10 2008-02-19 Mediatek Inc. Phase and frequency drift compensation in Orthogonal Frequency Division Multiplexing systems
KR100983273B1 (ko) * 2003-12-22 2010-09-24 엘지전자 주식회사 Dmb 수신기의 주파수 복구 장치 및 방법
CN1719818B (zh) 2005-07-08 2010-06-23 广州海格通信集团股份有限公司 一种正交频分复用系统的采样频率偏差跟踪信号处理装置及方法

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2015089125A (ja) * 2013-10-31 2015-05-07 富士通株式会社 マルチキャリア変調信号の定時同期化装置及び方法
JP2017517939A (ja) * 2014-04-24 2017-06-29 富士通株式会社 チャネル等化追跡装置、方法及び受信機
US9979572B2 (en) 2014-04-24 2018-05-22 Fujitsu Limited Channel equalization and tracking apparatus and method and receiver

Also Published As

Publication number Publication date
JP4320026B2 (ja) 2009-08-26
KR20070010721A (ko) 2007-01-24
US7702024B2 (en) 2010-04-20
CN1901527A (zh) 2007-01-24
KR100719112B1 (ko) 2007-05-17
US20070019763A1 (en) 2007-01-25
EP1746793A2 (en) 2007-01-24
CN100591064C (zh) 2010-02-17

Similar Documents

Publication Publication Date Title
JP4320026B2 (ja) サンプリング周波数誤差推定装置、サンプリング周波数誤差推定方法、直交周波数分割多重システムの受信機及びその位相歪み補正方法
KR100770924B1 (ko) 무선 통신 시스템에서 주파수 오차 보상 장치 및 방법
US8345734B2 (en) Time error estimation for data symbols
JP4159030B2 (ja) Ofdmを用いる無線ネットワーク用のタイミング同期方法
KR100699490B1 (ko) 샘플링 주파수 오프셋 추정방법 및 이 방법이 적용되는ofdm 시스템
US7110387B1 (en) System and method for compensating timing error using pilot symbol in OFDM/CDMA communication system
US6798738B1 (en) FFT window position recovery apparatus and method for OFDM system receiver
JP4263119B2 (ja) Ofdmシステムでの初期周波数の同期方法および装置
KR100719111B1 (ko) Ofdm 시스템에 적용되는 위상잡음 보상장치 및 그 방법
JP4125715B2 (ja) Ofdmシステムでの初期周波数の同期方法及び装置
US20040208267A1 (en) Frequency synchronization apparatus and method for OFDM systems
EP1964346B1 (en) Method and a system for estimating a symbol time error in a broadband transmission system
US7529179B1 (en) Joint maximum likelihood estimation of integer carrier frequency offset and channel in OFDM systems
US7792202B2 (en) Apparatus and method for estimating timing offset of OFDM symbol, and method of recovering symbol timing of OFDM symbol
US8665976B2 (en) Blind symbol synchronization scheme for OFDM system
JP3555843B2 (ja) 復調装置
KR20060067034A (ko) 직교 주파수 분할 다중화 통신시스템에서 주파수 및 시간옵셋 추정 방법과 그를 이용한 장치
US8023576B2 (en) Apparatus and method to estimate sampling offset
JP2007036549A (ja) Ofdm受信装置
JP2001069113A (ja) マルチキャリア同期検波回路
KR20160116994A (ko) 샘플링 클락 옵셋을 보상하는 방법 및 그 장치
JP6419016B2 (ja) Ofdm受信装置及びチップ
JP2005323077A (ja) 直交周波数分割多重受信装置

Legal Events

Date Code Title Description
RD04 Notification of resignation of power of attorney

Free format text: JAPANESE INTERMEDIATE CODE: A7424

Effective date: 20070119

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20081125

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20090203

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20090424

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: 20090519

A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20090529

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20120605

Year of fee payment: 3

R150 Certificate of patent or registration of utility model

Free format text: JAPANESE INTERMEDIATE CODE: R150

Ref document number: 4320026

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20130605

Year of fee payment: 4

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250