JP2000252969A - Method for detecting symbol timing - Google Patents

Method for detecting symbol timing

Info

Publication number
JP2000252969A
JP2000252969A JP11052060A JP5206099A JP2000252969A JP 2000252969 A JP2000252969 A JP 2000252969A JP 11052060 A JP11052060 A JP 11052060A JP 5206099 A JP5206099 A JP 5206099A JP 2000252969 A JP2000252969 A JP 2000252969A
Authority
JP
Japan
Prior art keywords
symbol
absolute value
sequence
correlation sequence
correlation
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
JP11052060A
Other languages
Japanese (ja)
Other versions
JP3619384B2 (en
Inventor
Kinichi Higure
欽一 日暮
Takehiko Kobayashi
岳彦 小林
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.)
Hitachi Denshi KK
Original Assignee
Hitachi Denshi KK
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 Hitachi Denshi KK filed Critical Hitachi Denshi KK
Priority to JP05206099A priority Critical patent/JP3619384B2/en
Publication of JP2000252969A publication Critical patent/JP2000252969A/en
Application granted granted Critical
Publication of JP3619384B2 publication Critical patent/JP3619384B2/en
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Abstract

PROBLEM TO BE SOLVED: To reduce the influence caused by eliminating a non-principal wave component by selecting a numerical series that has the least numerical amplitude among plural distributed numerical series and defining the timing containing the selected numerical series as the symbol timing. SOLUTION: The N sets of numerical series R0, R1...RN-1 which are stored in a two-dimensional register 107 of (N sets × F stages) are read out by N sets of arithmetic circuits 108 which calculate the rates Ti (i=0, 1...N-1) between two main waves showing the largest value and the second largest value of received signals and the absolute value of a correlative series for those numerical series respectively. A largest value retriever 109 retrieves the largest value of the rate Ti calculated between two main waves and the absolute value of the correlative series for each of series R0 to RN-1. Then (i) having the largest Ti is defined as i0 with a sample point where the remainder is i0 when N is divided by (n) defined as the symbol timing respectively, and i0 is outputted via an output terminal 110. Thus, the timing that contains a numerical series having the least numerical amplitude is defined as the symbol timing.

Description

【発明の詳細な説明】DETAILED DESCRIPTION OF THE INVENTION

【0001】[0001]

【発明の属する技術分野】本発明は、ディジタル通信に
おいて用いられる受信機で、受信信号から伝送信号を取
り出すときに使用されるシンボルタイミングの検出方法
に関するものである。
BACKGROUND OF THE INVENTION 1. Field of the Invention The present invention relates to a method for detecting a symbol timing used for extracting a transmission signal from a received signal in a receiver used in digital communication.

【0002】[0002]

【従来の技術】文字、音声、画像などのメディアを統合
したマルチメディアの利用が進展するに伴い、ディジタ
ル通信においても動画像の伝送が行なわれるようにな
り、動画像の伝送を行なうディジタル無線通信では、高
速での情報伝送が要求されるようになってきている。こ
のディジタル無線通信に用いられる受信機では、受信信
号の中から正しく文字、音声、画像などからなる伝送信
号を取り出すためにシンボルタイミングの検出を正確に
行なうことが重要であるが、ディジタル無線通信の一つ
の手段として使用される移動無線通信は、受信におい
て、特に、遅延波の影響を受けやすいという問題があ
る。送信機から送信された電波は、さまざまな経路を伝
搬し、それぞれ異なった遅延時間を伴って受信機で受信
される。この遅延波の影響を改善するため、受信機に波
形等化器が使用されるが、波形等化器は受信信号に含ま
れるシンボル信号を検出するタイミングにより性能が左
右されるため、適正なシンボルタイミングの検出が必要
となる。シンボルタイミングを検出する方法としては、
送信側から、送信信号の各フレームの先頭に既知のシン
ボル系列を配置して送信し、受信側において、受信信号
と既知のシンボル系列との相関系列を算出し、算出した
相関系列が最大値となる点をシンボルタイミングとする
シンボルタイミング検出方法がある。
2. Description of the Related Art As the use of multimedia that integrates media such as characters, voices, images, and the like progresses, moving images are also transmitted in digital communication, and digital wireless communication for transmitting moving images is performed. Then, high-speed information transmission is required. In a receiver used in this digital wireless communication, it is important to accurately detect symbol timing in order to correctly extract a transmission signal including characters, voices, images, etc. from a received signal. The mobile radio communication used as one means has a problem in reception that it is particularly susceptible to delayed waves. Radio waves transmitted from the transmitter propagate through various paths and are received by the receiver with different delay times. In order to improve the effect of the delayed wave, a waveform equalizer is used in the receiver, but the performance of the waveform equalizer depends on the timing of detecting the symbol signal included in the received signal. Timing detection is required. As a method of detecting the symbol timing,
From the transmitting side, a known symbol sequence is arranged and transmitted at the beginning of each frame of the transmission signal, and at the receiving side, a correlation sequence between the received signal and the known symbol sequence is calculated, and the calculated correlation sequence is the maximum value. There is a symbol timing detection method using a certain point as a symbol timing.

【0003】図4に示すシンボルタイミング検出回路を
使用して、従来のシンボルタイミング検出方法を説明す
る。なお、以下の説明において使用されるNは、1シン
ボルあたりのサンプル数、Mは、既知なシンボル系列の
シンボル数、Fは、1フレームあたりのシンボル数であ
る。受信機で受信された受信信号が処理され、所要の受
信信号が入力端子401を介して同期検波器402に入
力される。同期検波器402は、入力端子401を介し
入力される受信信号を同期検波してベースバンド信号4
03を取り出し、同期検波して取り出したベースバンド
信号403を低域通過フィルタ405へ出力する。低域
通過フィルタ405は、同期検波器402から入力され
る同期検波して取り出したベースバンド信号403の高
周波成分を除去し、高周波成分を除去したベースバンド
信号をA/D変換器406へ出力する。A/D変換器4
06は、低域通過フィルタ405から入力される高周波
成分を除去したベースバンド信号をシンボル間隔の1/
Nの時間間隔毎にサンプリングしてディジタル信号に変
換し、変換したディジタル信号を相関演算回路407へ
出力する。なお、以下の説明で、「シンボル間隔の1/
Nの時間間隔毎にサンプリング」を「1シンボルあたり
N回オーバーサンプリング」と呼ぶ。
A conventional symbol timing detection method will be described using the symbol timing detection circuit shown in FIG. In the following description, N is the number of samples per symbol, M is the number of symbols in a known symbol sequence, and F is the number of symbols per frame. A received signal received by the receiver is processed, and a required received signal is input to the synchronous detector 402 via the input terminal 401. The synchronous detector 402 performs synchronous detection on the received signal input via the input terminal 401 and performs baseband signal 4
03 is extracted, and the baseband signal 403 extracted by synchronous detection is output to the low-pass filter 405. The low-pass filter 405 removes the high-frequency component of the baseband signal 403 extracted by synchronous detection input from the synchronous detector 402 and outputs the baseband signal from which the high-frequency component has been removed to the A / D converter 406. . A / D converter 4
06 is the baseband signal from which the high-frequency component input from the low-pass filter 405 has been removed, that is,
The digital signal is sampled and converted into a digital signal at every N time intervals, and the converted digital signal is output to the correlation operation circuit 407. In the following description, "1/3 of the symbol interval"
"Sampling every N time intervals" is called "N times oversampling per symbol".

【0004】相関演算回路407は、A/D変換器40
6から入力されるディジタル信号と、別途、入力される
それぞれの長さMシンボルの既知シンボル系列{u
(0),u(1),・・・u(M−1)}との相関系列
{r(n)}を算出し、算出した相関系列{r(n)}
を絶対値演算回路408へ出力する。絶対値演算回路4
08は、相関演算回路407から入力される相関系列
{r(n)}の絶対値{|r(n)|}を算出し、算出し
た絶対値{|r(n)|}をレジスタ409へ出力する。
[0004] The correlation operation circuit 407 includes an A / D converter 40.
6 and a separately input known symbol sequence {u} of length M symbols.
(0), u (1),..., U (M−1)} are calculated, and the calculated correlation sequence {r (n)} is calculated.
Is output to the absolute value calculation circuit 408. Absolute value calculation circuit 4
08 calculates the absolute value {| r (n) |} of the correlation sequence {r (n)} input from the correlation operation circuit 407 and sends the calculated absolute value {| r (n) |} to the register 409. Output.

【0005】ここで、前記相関演算回路407の詳細な
ブロック図を図3に示し説明する。相関演算回路407
は、入力端子301に接続しているM−1段の縦続接続
された遅延素子302と、M−1段の縦続接続された遅
延素子302の第1段目の遅延素子302の入力側およ
びM−1段の遅延素子302の各出力側に、一方の乗算
入力が接続され、他方の乗算入力に既知シンボル系列
{u}が入力するようにそれぞれ接続されたM個の乗算
回路303と、M個の乗算回路303の各出力信号をつ
ぎつぎに積算していくように接続されたM−1個の加算
回路304とで構成され、最後の加算回路304の積算
出力が出力端子305に接続した構成となっている。
Here, a detailed block diagram of the correlation operation circuit 407 will be described with reference to FIG. Correlation operation circuit 407
Are the M-1 cascaded delay elements 302 connected to the input terminal 301, the input side of the first-stage delay element 302 of the M-1 cascaded delay elements 302, and M One multiplication input is connected to each output side of the -1 stage delay element 302, and M multiplication circuits 303 each connected so that the known symbol sequence {u} is input to the other multiplication input; A configuration in which M-1 addition circuits 304 are connected so as to successively integrate the output signals of the multiplication circuits 303, and the integration output of the last addition circuit 304 is connected to an output terminal 305. It has become.

【0006】入力端子301にA/D変換器406(図
4参照)から入力されるディジタル信号は、縦続接続さ
れたM−1段の遅延素子302によりそれぞれ所要時間
遅延されつつシフトする。M−1段の遅延素子302で
それぞれ所要時間遅延された各ディジタル信号は、M個
の乗算回路303の2個目からM個目までの乗算回路3
03(1個目の乗算回路303には、入力端子301か
ら遅延されていないディジタル信号が入力される)へ出
力される。このディジタル信号は、A/D変換器406
において1シンボルあたりN回オーバーサンプリングさ
れており、M−1段の遅延素子302のそれぞれにおい
て、Nサンプル、すなわち、1シンボル分づつ遅延す
る。
A digital signal input from an A / D converter 406 (see FIG. 4) to an input terminal 301 is shifted while being delayed for a required time by an M-1 stage delay element 302 connected in cascade. Each digital signal delayed for a required time by the delay element 302 of the (M-1) -th stage is multiplied by the second to M-th multiplication circuits 3 of the M multiplication circuits 303.
03 (a digital signal that has not been delayed from the input terminal 301 is input to the first multiplication circuit 303). This digital signal is supplied to an A / D converter 406
Are oversampled N times per symbol, and each of the delay elements 302 in the (M-1) -th stage delays by N samples, that is, by one symbol.

【0007】ここで、入力端子301から入力される信
号系列を{y}として、時刻nにおけるM−1段目の遅
延素子302の出力信号をy(n)とするために、入力
端子301に入力される入力信号をy(n+(M−1)
N)とする。この状態において、あるサンプルタイミン
グで、入力端子301に入力された入力信号{y
(n),y(n+N),・・・y(n+(M−1)
N)}と既知シンボル系列{u}との相関系列が算出さ
れるとすると、つぎのサンプルタイミングでは入力信号
{y(n+1),y(n+N+1),・・・y(n+
(M−1)N+1)}と既知シンボル系列{u}との相
関系列が算出される。なお、相関演算は、最初に入力端
子301から入力されたディジタル信号が、M−1段の
遅延素子302をつぎつぎにシフトされ、M−1段目の
遅延素子302から出力される時点から開始し、M−1
段目の遅延素子302にy(n)が出力された時点での
相関系列をr(n)とする。この相関系列r(n)は、
M個の乗算回路303それぞれの入力端子306に既知
シンボル系列{u(0),u(1),・・・u(M−
1)}に近いパターンの系列が入力されるとき、絶対値
が最大値を与える。
Here, in order to set the signal sequence input from the input terminal 301 to {y} and to set the output signal of the (M−1) th stage delay element 302 at time n to y (n), the input terminal 301 The input signal to be input is y (n + (M-1)
N). In this state, the input signal {y input to the input terminal 301 at a certain sample timing
(N), y (n + N), ... y (n + (M-1)
Assuming that a correlation sequence between (N)} and the known symbol sequence {u} is calculated, input signals {y (n + 1), y (n + N + 1),.
A correlation sequence between (M−1) N + 1)} and the known symbol sequence {u} is calculated. Note that the correlation operation starts from the point in time when the digital signal first input from the input terminal 301 is successively shifted through the delay element 302 of the (M-1) th stage and output from the delay element 302 of the (M-1) th stage. , M-1
The correlation sequence at the time when y (n) is output to the delay element 302 at the stage is r (n). This correlation sequence r (n) is
A known symbol sequence {u (0), u (1),... U (M−
1) When a series of patterns close to} is input, the absolute value gives the maximum value.

【0008】さらに図4において、絶対値演算回路40
8から算出した絶対値{|r(n)|}が入力されるレジ
スタ409は、1フレームあたりのサンプル数FN個
(Fは、1フレームあたりのシンボル数、Nは、1シン
ボルあたりのサンプル数)のレジスタである。レジスタ
409は、1フレーム分の相関系列の絶対値{|r
(0)|,|r(1)|,・・・|r(FN−1)|}を記
憶する。最大値検索器410は、レジスタ409から1
フレーム分の相関系列の絶対値を読み出し、最大値を検
索し、最大値を与える点をn0 とする。最大値を与える
点n0 を、1シンボルあたりのサンプル数Nで割った余
りn0modNを算出し、n0 をNで割った余りがn0
modNと等しくなるサンプル点をシンボルタイミング
とし、n0 modNを出力端子411へ出力する。
[0008] Further, in FIG.
The register 409, to which the absolute value {| r (n) |} calculated from Eq. 8 is input, stores the number of samples FN per frame (F is the number of symbols per frame, and N is the number of samples per symbol). ). The register 409 stores the absolute value of the correlation sequence for one frame {| r
(0) |, | r (1) |,... | R (FN-1) |} are stored. The maximum value searcher 410 reads 1
The absolute value of the correlation sequence for the frame is read, the maximum value is searched, and the point giving the maximum value is defined as n 0 . The remainder n 0 mod N is calculated by dividing the point n 0 giving the maximum value by the number N of samples per symbol, and the remainder obtained by dividing n 0 by N is n 0.
A sample point equal to modN is set as a symbol timing, and n 0 modN is output to the output terminal 411.

【0009】図5に示す受信信号の相関系列は、図7に
1フレームのシンボル構成を示すように、1フレームが
192シンボル(F=192)の先頭に16シンボル
(M=16)の長さの基準とする既知シンボル系列を配
置し、π/4シフトQPSKで信号点配置して送信し、
遅延時間が1シンボル周期の遅延波(D/U(所要信号
対不要信号の比)=0dB)1波が加わった場合の相関
系列(N=4サンプル/シンボル)例である。また、図
6は、図5に示す受信信号の相関系列を最大値付近で拡
大表示した図である。図5および図6に示すように、サ
ンプル数がn=60で相関系列の絶対値 |r(n)| が
最大値となるため、最大値検索器410でn0 =60が
選択され、n0 modN=60mod4=0となり、n
を4で割った余りが0となるサンプル点(n=0,4,
8,・・・)をシンボルタイミングとし、0を出力端子
411へ出力する。
The correlation sequence of the received signal shown in FIG. 5 has a length of 16 symbols (M = 16) at the beginning of 192 symbols (F = 192) as shown in FIG. A known symbol sequence as a reference is arranged, and signal points are arranged and transmitted using π / 4 shift QPSK.
This is an example of a correlation sequence (N = 4 samples / symbol) when one delay wave (D / U (ratio of required signal to unnecessary signal) = 0 dB) with a delay time of one symbol period is added. FIG. 6 is an enlarged view of the correlation sequence of the received signal shown in FIG. 5 near the maximum value. As shown in FIGS. 5 and 6, the absolute value | r (n) | of the correlation sequence becomes the maximum value when the number of samples is n = 60, so that n 0 = 60 is selected by the maximum value searcher 410, and n 0 mod N = 60 mod 4 = 0, and n
Points at which the remainder of dividing by 4 is 0 (n = 0, 4,
8,...) Is set as the symbol timing, and 0 is output to the output terminal 411.

【0010】シンボルタイミングの検出方法において
は、一般的に受信信号は帯域制限されているので、送信
信号系列がランダムな信号系列であると考え、遅延波が
存在しないとすると、受信信号と既知シンボル系列との
相関系列r(n)は、図8に示す直接波と既知シンボル
系列との相関系列801のようになる。また、遅延時間
が1シンボル、D/U=0dB、直接波と同位相の遅延
波が存在するとすると、遅延波と既知シンボル系列との
相関系列r(n)は、図8に破線で示す遅延波と既知シ
ンボル系列との相関系列802のようになる。また、直
接波と遅延波とを合成した受信信号と既知シンボル系列
との相関系列r(n)は、直接波と既知シンボル系列と
の相関系列801と遅延波と既知シンボル系列との相関
系列802との和となる合成した受信信号と既知シンボ
ル系列との相関系列803のようになる。
In the symbol timing detection method, the received signal is generally band-limited, so that the transmitted signal sequence is considered to be a random signal sequence, and if there is no delayed wave, the received signal and the known symbol The correlation sequence r (n) with the sequence is like the correlation sequence 801 between the direct wave and the known symbol sequence shown in FIG. Assuming that a delay time is one symbol, D / U = 0 dB, and a delayed wave having the same phase as the direct wave exists, a correlation sequence r (n) between the delayed wave and the known symbol sequence is represented by a delay indicated by a broken line in FIG. It becomes like a correlation sequence 802 between a wave and a known symbol sequence. The correlation sequence r (n) between the received signal obtained by combining the direct wave and the delayed wave and the known symbol sequence is a correlation sequence 801 between the direct wave and the known symbol sequence and a correlation sequence 802 between the delayed wave and the known symbol sequence. And a correlation sequence 803 between the combined received signal and the known symbol sequence, which is the sum of

【0011】さらに、図9(a)に、前述の合成した受
信信号と既知シンボル系列との相関系列803の絶対値
を算出した絶対値|r(n)|901を図示する。従来技
術によるシンボルタイミング検出方法では、合成した受
信信号と既知シンボル系列との相関系列の絶対値901
が最大値となる点nMAX 902が検出され、図9(b)
に示す最大値となる点nMAX 902を基準にシンボル周
期で抜き出した相関系列903がシンボルタイミングと
なる。
FIG. 9A shows an absolute value | r (n) | 901 obtained by calculating the absolute value of the correlation sequence 803 between the synthesized received signal and the known symbol sequence. In the symbol timing detection method according to the related art, an absolute value 901 of a correlation sequence between a combined received signal and a known symbol sequence is used.
The point n MAX 902 where is the maximum value is detected, and FIG.
The correlation sequence 903 extracted at the symbol period based on the point n MAX 902 having the maximum value shown in FIG.

【0012】遅延波の影響を改善する波形等化器は、大
きく分類すると判定帰還型等化器(DFE)と最尤系列
推定器(ビタビ等化器)の2つに分けることができる
が、ここでは後者の最尤系列推定器を等化器として使用
する場合の従来技術によるシンボルタイミング検出方法
の問題について考える。図10に、最尤系列推定器の概
念を仮想回路とし、ブロック図化したものを示す。伝搬
路推定器1006で、伝送路1002を伝送されてきた
受信信号y(t)1005から、伝送路1002の伝搬
路特性h(t)の推定値h’(t)を求める。ある適当
な相関系列{an }を、伝送路特性h(t)の推定値
h’(t)の推定伝送路1007に通したレプリカy’
(t;{an })1008を算出する。受信信号y
(t)1005とレプリカy’(t;{an })100
8とを比較回路1009で比較し、受信信号y(t)1
005とレプリカy’(t;{an})1008との信
号間距離d({an })が最小となる相関系列{an
を検出回路1010で検索する。信号間距離d({a
n })が最小となる相関系列を検索するためには、相関
系列{an }のすべてについて検索しなければならない
が、ビタビアルゴリズムを使用することにより検索を効
率的に行なうことができる。
The waveform equalizer for improving the influence of the delay wave can be roughly classified into two types, a decision feedback equalizer (DFE) and a maximum likelihood sequence estimator (Viterbi equalizer). Here, the problem of the symbol timing detection method according to the related art when the latter maximum likelihood sequence estimator is used as an equalizer will be considered. FIG. 10 is a block diagram showing a concept of a maximum likelihood sequence estimator as a virtual circuit. The propagation path estimator 1006 obtains an estimated value h ′ (t) of the propagation path characteristic h (t) of the transmission path 1002 from the received signal y (t) 1005 transmitted through the transmission path 1002. There suitable correlation sequence {a n}, and the transmission path characteristic h (t) of the estimated value h 'replicas through estimated CIR 1007 (t) y'
Calculating a; ({a n} t) 1008. Received signal y
(T) 1005 and the replica y '(t; {a n }) 100
8 is compared by a comparison circuit 1009, and the received signal y (t) 1
005 and replica y '(t; {a n }) between the signals of 1008 distance d ({a n}) correlation sequence which is minimum {a n}
Is searched by the detection circuit 1010. Signal distance d ({a
In order to search for a correlation sequence that minimizes n {), it is necessary to search for all of the correlation sequences {a n }, but the search can be performed efficiently by using the Viterbi algorithm.

【0013】[0013]

【発明が解決しようとする課題】従来技術によるシンボ
ルタイミング検出方法においては、シンボルタイミング
を決定した場合、図9(b)に相関系列903として示
すように無限に続き、非主力波成分の除去を行なうと、
非主力波成分を除去したことの影響が大きいという問題
がある。本発明は、前記問題を解決し、図9(c)に相
関系列904として示すように、相関系列の主力波2波
の割合が最大となるようにシンボルタイミングを決定
し、非主力波成分を除去したことの影響が少なくなるよ
うにしたシンボルタイミング検出方法を提供することを
目的とする。
In the symbol timing detection method according to the prior art, when the symbol timing is determined, the non-main component is removed infinitely, as shown as a correlation sequence 903 in FIG. 9B. When you do
There is a problem that the effect of removing the non-mainstream wave component is great. The present invention solves the above problem, and determines the symbol timing so that the ratio of two main waves in the correlation sequence is maximized, as shown as a correlation sequence 904 in FIG. An object of the present invention is to provide a symbol timing detection method in which the influence of the removal is reduced.

【0014】[0014]

【課題を解決するための手段】前記目的を達成するため
に、本発明のシンボルタイミング検出方法は、ディジタ
ル通信に用いる受信機に使用するシンボルタイミング検
出方法であって、複数のシンボルからなる各フレームの
先頭に長さMシンボルの既知なシンボル系列を配置した
フレーム構成の送信信号を前記受信機で受信して受信信
号とし、該受信信号を同期検波してベースバンド信号と
してから1シンボルあたりN回(Nは自然数)オーバー
サンプリングしてディジタル信号とし、該ディジタル信
号の前記既知なシンボル系列との相関を演算して相関系
列を算出し、該相関系列の絶対値からシンボルタイミン
グを検出するシンボルタイミング検出方法において、前
記相関系列の絶対値を1シンボルあたりのサンプル数N
の周期で、N組の数値系列に振り分け、該振り分けたN
組の数値系列のそれぞれについて絶対値が大きい2つの
絶対値が前記相関系列の絶対値に占める割合を算出し、
該割合が最も大きくなる相関系列を検索し、検索した相
関系列が存在するタイミングをシンボルタイミングとす
るものである。
In order to achieve the above object, a symbol timing detecting method according to the present invention is a symbol timing detecting method used for a receiver used in digital communication, wherein each symbol comprises a plurality of symbols. The receiver receives a transmission signal having a frame structure in which a known symbol sequence having a length of M symbols is arranged at the head of the reception signal as a reception signal, and synchronously detects the reception signal to form a baseband signal, and then N times per symbol. (N is a natural number) Oversampling into a digital signal, calculating a correlation of the digital signal with the known symbol sequence, calculating a correlation sequence, and detecting a symbol timing from an absolute value of the correlation sequence. In the method, the absolute value of the correlation sequence is calculated as the number of samples per symbol N
In a cycle of N, the set is divided into N sets of numerical sequences, and the sorted N
Calculating a ratio of two absolute values having a large absolute value to the absolute value of the correlation sequence for each of the set of numerical sequences;
The correlation sequence having the largest ratio is searched, and the timing at which the searched correlation sequence exists is used as the symbol timing.

【0015】さらに、詳細には、本発明のシンボルタイ
ミング検出方法は、振り分けたN組の数値系列のそれぞ
れについて絶対値が大きい2つの絶対値が相関系列の絶
対値に占める割合の算出は、前記振り分けたN組の数値
系列のそれぞれについて絶対値の最大値を算出し、該絶
対値の最大値を与えるサンプル点から前後1シンボルの
範囲内で、絶対値の最大値および2番目に大きい値を検
索し、該絶対値の最大値および2番目に大きい値が前記
相関系列の絶対値に占める割合を算出し、N組の数値系
列のそれぞれから前記割合が最も大きくなる相関系列を
検索し、検索した相関系列が存在するタイミングをシン
ボルタイミングとするものである。
More specifically, in the symbol timing detection method according to the present invention, the ratio of two absolute values having a large absolute value to the absolute value of the correlation sequence is calculated for each of the sorted N sets of numerical sequences. The maximum value of the absolute value is calculated for each of the allocated N sets of numerical sequences, and the maximum value of the absolute value and the second largest value within the range of one symbol before and after the sample point giving the maximum value of the absolute value are calculated. Calculating the ratio of the absolute value of the maximum value and the second largest value of the absolute value to the absolute value of the correlation sequence, searching for the correlation sequence having the highest ratio from each of the N sets of numerical sequences; The timing at which the calculated correlation sequence exists is used as the symbol timing.

【0016】また、本発明のシンボルタイミング検出方
法は、絶対値の最大値および2番目に大きい値が相関系
列の絶対値に占める割合の算出は、前記絶対値の最大値
を与えるサンプル点から前後(M−1)シンボルの範囲
内で、絶対値の最大値および2番目に大きい値が相関系
列の絶対値に占める割合の算出し、N組の数値系列のそ
れぞれから前記割合が最も大きくなる相関系列を検索
し、検索した相関系列が存在するタイミングをシンボル
タイミングとするものである。
In the symbol timing detecting method according to the present invention, the maximum absolute value and the ratio of the second largest value to the absolute value of the correlation sequence are calculated before and after the sample point giving the maximum absolute value. (M-1) Within the range of the symbol, the maximum value of the absolute value and the ratio of the second largest value to the absolute value of the correlation sequence are calculated, and the correlation in which the ratio is the largest from each of the N sets of numerical sequences is calculated. A sequence is searched, and the timing at which the searched correlation sequence exists is used as the symbol timing.

【0017】また、本発明のシンボルタイミング検出方
法は、ディジタル通信に用いる受信機に使用するシンボ
ルタイミング検出方法であって、所定のビット数のデー
タの複数のシンボルからなる各フレームの先頭に長さM
シンボルの既知なシンボル系列を配置したフレーム構成
の送信信号を受信して受信信号とし、該受信信号を同期
検波器で検波してベースバンド信号とし、該ベースバン
ド信号をA/D変換器で1シンボルあたりN回(Nは自
然数)オーバーサンプリングしてディジタル信号とし、
該ディジタル信号を相関演算回路で前記既知なシンボル
系列との相関を演算して相関系列を算出し、該相関系列
を絶対値演算回路で演算して相関系列の絶対値を算出
し、該相関系列の絶対値の少なくとも1フレームをレジ
スタに記憶し、該レジスタから前記相関系列の絶対値の
少なくとも1フレームを読み出して最大値検索器で最大
値を検索し、該最大値を与えるサンプル点をシンボルタ
イミングとするシンボルタイミング検出方法において、
前記絶対値演算回路で算出した前記相関系列の絶対値
を、循環切換器で順次に、循環して1シンボルあたりの
サンプル数Nの周期で、N組の数値系列に振り分け、該
振り分けたN組の数値系列をN組×F段の2次元レジス
タに記憶し、該N組×F段の2次元レジスタからそれぞ
れの組に接続されているN組の演算回路に読み出して、
前記N組の数値系列のそれぞれについて絶対値の最大値
および2番目に大きい値が前記相関系列の絶対値に占め
る割合を算出し、前記N組の数値系列のそれぞれの割合
から最大値検索器で前記割合が最も大きくなる相関系列
を検索し、検索した前記割合が最も大きくなる相関系列
が存在するタイミングをシンボルタイミングとするもの
である。
The symbol timing detection method according to the present invention is a symbol timing detection method used for a receiver used for digital communication, wherein a symbol having a predetermined number of bits has a length at the beginning of each frame. M
A transmission signal having a frame configuration in which a known symbol sequence of symbols is arranged is received and used as a reception signal, and the reception signal is detected by a synchronous detector to generate a baseband signal, and the baseband signal is converted into one by an A / D converter. Oversampling N times (N is a natural number) per symbol into a digital signal,
The digital signal is correlated with the known symbol sequence by a correlation calculation circuit to calculate a correlation sequence, and the correlation sequence is calculated by an absolute value calculation circuit to calculate the absolute value of the correlation sequence. Is stored in a register, at least one frame of the absolute value of the correlation sequence is read out from the register, the maximum value is searched by a maximum value searcher, and a sample point giving the maximum value is determined by symbol timing. In the symbol timing detection method,
The absolute values of the correlation sequences calculated by the absolute value calculation circuit are sequentially circulated by a circulation switcher and distributed to N sets of numerical sequences at a period of N samples per symbol, and the N sets of the allocated sets are arranged. Is stored in an N-set × F-stage two-dimensional register, and is read out from the N-set × F-stage two-dimensional register to N-set arithmetic circuits connected to the respective sets.
The maximum value of the absolute value and the ratio of the second largest value to the absolute value of the correlation sequence are calculated for each of the N sets of numerical sequences. The correlation sequence with the largest ratio is searched, and the timing at which the searched correlation sequence with the largest ratio exists is defined as the symbol timing.

【0018】[0018]

【発明の実施の形態】本発明の実施の形態を説明する前
に、まず、本発明のシンボルタイミング検出方法につい
て基本的な説明をする。図1に示す、シンボルタイミン
グを検出する基本的なブロック図を使用する。送信側か
らは、図7に示すように、1つのフレームごとに、長さ
Mシンボルの基準とする既知のシンボル系列{u(k)
|k=0,1,・・・M−1}を先頭に配置したフレー
ム構成の送信信号が送信される。受信側では、受信され
た受信信号が所要の処理をされて、既知のシンボル系列
{u(k)|k=0,1,・・・M−1}を先頭に配置
したフレーム構成のディジタル信号(図1には記載して
いないが、A/D変換器において、1シンボルあたりN
回オーバーサンプリングされている)となり、相関演算
回路102の入力端子101を介し、縦続接続されてい
るM−1段のNサンプル遅延素子103の第1段の遅延
素子103およびM個の乗算回路104の第1の乗算回
路104に入力される。縦続接続されているM−1段の
Nサンプル遅延素子103は、入力される1シンボルあ
たりN回オーバーサンプリングされたディジタル信号に
それぞれNサンプル分の遅延を与えるため、例えば、第
1段の遅延素子103に信号系列y(n)が入力された
とすると、第1段の遅延素子103の出力はNサンプル
前のy(n−N)、第2段の遅延素子103の出力は2
Nサンプル前のy(n−2N)、最後の第M−1段の遅
延素子103の出力は(M−1)Nサンプル前のy(n
−(M−1)N)となる。
DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS Before describing an embodiment of the present invention, a basic description will first be given of a symbol timing detecting method of the present invention. The basic block diagram for detecting the symbol timing shown in FIG. 1 is used. From the transmitting side, as shown in FIG. 7, for each frame, a known symbol sequence と す る u (k) as a reference for a length M symbol
| K = 0, 1,..., M−1} are transmitted at the beginning. On the receiving side, the received signal is processed as required, and a known symbol sequence {u (k) | k = 0, 1,... (Although not shown in FIG. 1, in the A / D converter, N
Times over-sampled), the first delay element 103 of the M-1 N-sample delay elements 103 and the M multiplication circuits 104 connected in cascade through the input terminal 101 of the correlation operation circuit 102. Is input to the first multiplication circuit 104. The M-stage N-sample delay element 103 connected in cascade gives a delay of N samples to the input digital signal oversampled N times per symbol. Assuming that the signal sequence y (n) is input to 103, the output of the first-stage delay element 103 is y (n-N) N samples before, and the output of the second-stage delay element 103 is 2
The output of y (n−2N) before N samples and the output of the last M−1 stage delay element 103 are y (n) before (M−1) N samples.
− (M−1) N).

【0019】以下の説明では、最後の第M−1段の遅延
素子103の出力をy(n)とするために、図1にも記
載しているように、第1の遅延素子103の入力をy
(n+(M−1)N)、出力をy(n+(M−2)
N)、第2の遅延素子103の出力をy(n+(M−
3)N)(以下、同様)とする。第1の遅延素子103
の入力および各遅延素子103の出力にはそれぞれ乗算
回路104が接続されており、別途それぞれ入力される
既知シンボル系列{u}の各係数との積が算出され、M
−1個の加算回路105に入力され、M−1個の加算回
路105それぞれにおいて、各遅延素子103の出力と
既知シンボル系列{u}の各係数との積が積算され、入
力信号y(n)と既知シンボル系列{u}との相関系列
r(n)が算出される。このような算出を計算式で示す
と、下記(1)式となる。
In the following description, in order to set the output of the last (M-1) th delay element 103 to y (n), as shown in FIG. To y
(N + (M-1) N) and output is y (n + (M-2)
N), the output of the second delay element 103 is y (n + (M−
3) N) (the same applies hereinafter). First delay element 103
, And the output of each delay element 103 are connected to a multiplication circuit 104, and the product of each of the coefficients of the known symbol sequence {u} separately input is calculated.
-1 adder circuits 105, and in each of the M-1 adder circuits 105, the product of the output of each delay element 103 and each coefficient of the known symbol sequence {u} is integrated, and the input signal y (n ) And the known symbol sequence {u} are calculated as a correlation sequence r (n). This calculation is represented by the following equation (1).

【数1】 (Equation 1)

【0020】相関演算回路102において、(1)式に
より算出した相関系列r(n)を絶対値演算回路111
に入力し、絶対値演算回路111で算出した絶対値 |r
(n)| をNの周期で循環切換器106によりN組(周
期Nと合致させた組数)の数値系列 R0 ={|r(0)|,|r(N)|,・・・|r((F−1)N)|}, R1 ={|r(1)|,|r(N+1)|,・・・|r((F−1)N+1)|}, ・・・・・・ RN-1 ={|r(N−1)|,|r(2N−1),・・・|r(FN−1)|} に振り分け、N組×F段(1フレームあたりのシンボル
数Fと合致させた段数で、1つの組の段数)の2次元レ
ジスタ107に記憶する。例えば、N=4の場合、n=
5ならば2次元レジスタ107の1組目の2段目に、n
=12ならば4組目の3段目にそれぞれ記憶する。循環
切換器106によって2次元レジスタ107に振り分け
られたN組の数値系列R0 ,R1 ,・・・RN-1 のそれ
ぞれについて、主力波2波の割合を算出する。
In the correlation operation circuit 102, the correlation sequence r (n) calculated by the equation (1) is used as an absolute value operation circuit 111.
And the absolute value | r calculated by the absolute value calculation circuit 111
(N) | is set to N series (the number of sets matched with the cycle N) by the circulation switching unit 106 at N cycles. R 0 = {| r (0) |, | r (N) |,... | r ((F−1) N) |}, R 1 = {| r (1) |, | r (N + 1) |,... | r ((F−1) N + 1) |},. ... R N-1 = {| r (N-1) |, | r (2N-1), ... | r (FN-1) |}, and N sets × F stages (per frame) Is stored in the two-dimensional register 107 of the number of stages matched with the number of symbols F of (1). For example, if N = 4, n =
If it is 5, the second stage of the first set of the two-dimensional register 107 has n
If = 12, they are stored in the third row of the fourth set. For each of the N sets of numerical sequences R 0 , R 1 ,... R N−1 distributed to the two-dimensional register 107 by the circulation switch 106, the ratio of the two main force waves is calculated.

【0021】この主力波2波の場合の割合算出を、図1
2を使用して説明する。図12に示す相関系列の各絶対
値1201は、循環切換器106(図1参照)が絶対値
演算回路111から入力される絶対値|r(n)|を振り
分けた相関系列の絶対値である。まず、循環切換器10
6によって振り分けられる前の相関系列r(n)(n=
0,1,・・・FN−1)が最大となるサンプル点n
max 1202の前後1シンボルの範囲1203で、相関
系列の各絶対値1201の中の最大値rmax,i 1204
と2番目に大きい値r2nd,i 1205(主力波2波)と
を検索し、これら最大値と2番目に大きい値とが存在す
るサンプル点をそれぞれmmax,i 、m2nd,iとする。つ
ぎに、mmax,i を中心に前後(M−1)シンボルの範囲
1206で相関系列の絶対値1201の和Si を算出
し、さらに、主力波2波1204、1205の和(r
max,i +r2nd,i )と相関系列の絶対値1201の和S
i との比、Ti =(rmax,i +r2nd,i )/Si を算出
する。各数値系列R0 ,R1 ,・・・RN-1 のそれぞれ
について主力波2波と相関系列の絶対値1201との割
合Ti を計算した後、この割合Ti の最大値を最大値検
索器109(図1参照)により検索し、そのときのiを
0 とし、nをNで割った余りがi0 となるサンプル点
をシンボルタイミングとし、i0 を出力端子110から
出力する。
The calculation of the ratio in the case of two main force waves is shown in FIG.
2 will be described. Each absolute value 1201 of the correlation sequence shown in FIG. 12 is an absolute value of the correlation sequence obtained by the circulation switch 106 (see FIG. 1) sorting the absolute value | r (n) | input from the absolute value calculation circuit 111. . First, the circulation switch 10
6, the correlation sequence r (n) (n =
0, 1,... FN-1) is the maximum sample point n
Within a range 1203 of one symbol before and after max 1202, the maximum value r max, i 1204 of each absolute value 1201 of the correlation sequence
And the second largest value r 2nd, i 1205 (two main waves) are searched, and the sample points where these maximum value and the second largest value exist are defined as m max, i and m 2nd, i , respectively. . Next, the sum S i of the absolute value 1201 of the correlation sequence is calculated in a range 1206 of (M−1) symbols before and after m max, i , and further, the sum (r
max, i + r 2nd, i ) and the sum S of the absolute value 1201 of the correlation sequence
The ratio of the i, to calculate the T i = (r max, i + r 2nd, i) / S i. After calculating the ratio T i between the two main waves and the absolute value 1201 of the correlation sequence for each of the numerical sequences R 0 , R 1 ,... R N−1 , the maximum value of this ratio T i is set to the maximum value. searching by the search unit 109 (see FIG. 1), the i at that time and i 0, the sample point remainder of dividing n by n becomes i 0 to the symbol timing, and outputs the i 0 from the output terminal 110.

【0022】つぎに、上述の基本原理にもとづく、本発
明のシンボルタイミング検出方法の実施の形態を説明す
る。図2は、受信機で使用される、本発明のシンボルタ
イミング検出方法を使用したシンボルタイミング検出回
路の一実施例を示す。送信機からは、図7に示すよう
に、1つのフレームごとに、基準とする既知シンボル系
列を先頭に配置したフレーム構成の送信信号が送信され
る。受信機では、受信された受信信号が所要の処理をさ
れ、1つのフレーム(Fシンボルで構成)ごとに基準と
なる既知シンボル系列(Mシンボルで構成)を先頭に配
置したフレーム構成の受信信号となり、図2に示す本発
明のシンボルタイミング検出方法を使用したシンボルタ
イミング検出回路に入力される。
Next, an embodiment of the symbol timing detecting method according to the present invention based on the above basic principle will be described. FIG. 2 shows an embodiment of a symbol timing detection circuit used in a receiver and using the symbol timing detection method of the present invention. As shown in FIG. 7, the transmitter transmits a transmission signal having a frame structure in which a reference known symbol sequence is arranged at the head for each frame. In the receiver, the received signal received is subjected to required processing to become a received signal having a frame configuration in which a known symbol sequence (constituted by M symbols) serving as a reference is placed at the head of each frame (constituted by F symbols). Are input to a symbol timing detection circuit using the symbol timing detection method of the present invention shown in FIG.

【0023】受信信号は、入力端子201を介し、同期
検波器202に入力される。同期検波器202は、入力
端子201を介し入力される受信信号からベースバンド
信号203を同期検波して取り出し、同期検波して取り
出したベースバンド信号203を低域通過フィルタ20
5へ出力する。低域通過フィルタ205は、同期検波器
202から入力される同期検波して取り出したベースバ
ンド信号203の高周波成分を除去し、高周波成分を除
去したベースバンド信号をA/D変換器206へ出力す
る。A/D変換器206は、低域通過フィルタ205か
ら入力される高周波成分を除去したベースバンド信号を
シンボル間隔の1/Nの時間間隔毎にサンプリングして
ディジタル信号に変換し、変換したディジタル信号を相
関演算回路102へ出力する。
The received signal is input to the synchronous detector 202 via the input terminal 201. The synchronous detector 202 synchronously detects and extracts the baseband signal 203 from the received signal input via the input terminal 201, and synchronously detects and extracts the extracted baseband signal 203 from the low-pass filter 20.
Output to 5 The low-pass filter 205 removes the high-frequency component of the baseband signal 203 extracted by synchronous detection input from the synchronous detector 202 and outputs the baseband signal from which the high-frequency component has been removed to the A / D converter 206. . The A / D converter 206 samples the baseband signal from which the high-frequency component input from the low-pass filter 205 has been removed at every 1 / N of a symbol interval, converts it into a digital signal, and converts the converted digital signal. Is output to the correlation operation circuit 102.

【0024】相関演算回路102は、A/D変換器20
6から入力されるディジタル信号と、別途、入力される
長さMシンボルの基準とする既知シンボル系列{u
(0),u(1),・・・u(M−1)}との相関系列
{r(n)}を算出し、算出した相関系列{r(n)}
を絶対値演算回路(例えば、2乗回路)111へ出力す
る。絶対値演算回路111は、相関演算回路102から
入力される相関系列{r(n)}の絶対値{|r(n)
|}を算出し、算出した絶対値{|r(n)|}を循環切
換器106へ出力する。循環切換器106は、絶対値演
算回路111から入力される絶対値 {|r(n)|} を
Nの周期で、nをNで割った余りによりN組(周期Nと
合致させた組数)の数値系列に振り分け、N組×F段
(1フレームあたりのシンボル数Fと合致させた段数
で、1つの組の段数)の2次元レジスタ107へ出力す
る。2次元レジスタ107は、循環切換器106から入
力されるN組に振り分けられた数値系列 R0 ={r(0),r(N),・・・r((F−1)N)}, R1 ={r(1),r(N+1),・・・r((F−1)N+1)}, ・・・・・ RN-1 ={r(N−1),r(2N−1),・・・r(FN−1)} を記憶する。
The correlation operation circuit 102 includes an A / D converter 20
6 and a known symbol sequence {u} separately used as a reference for a length M symbol that is input.
(0), u (1),..., U (M−1)} are calculated, and the calculated correlation sequence {r (n)} is calculated.
Is output to an absolute value calculation circuit (for example, a squaring circuit) 111. The absolute value calculation circuit 111 calculates the absolute value {| r (n) of the correlation sequence {r (n)} input from the correlation calculation circuit 102.
Is calculated, and the calculated absolute value {| r (n) |} is output to the circulation switch 106. The circulation switching unit 106 converts the absolute value {| r (n) |} input from the absolute value calculation circuit 111 into N cycles and modulo n by N to obtain N sets (the number of sets matched with the cycle N). ), And outputs to the two-dimensional register 107 of N sets × F stages (the number of stages in one set with the number of stages matched with the number of symbols F per frame). The two-dimensional register 107 stores a numerical sequence R 0 = {r (0), r (N),... R ((F−1) N)}, which is distributed to N sets input from the circulation switch 106. R 1 = {r (1), r (N + 1),... R ((F−1) N + 1)},..., R N−1 = {r (N−1), r (2N−) 1),... R (FN-1)} are stored.

【0025】N組×F段の2次元レジスタ107に記憶
されたN組の数値系列R0 ,R1 ,・・・RN-1 は、N
組×F段の2次元レジスタ107のそれぞれの組に接続
されているN組の演算回路108にそれぞれ読み出され
る。N組の演算回路108は、2次元レジスタ107か
ら読み出されるN組の数値系列R0 ,R1 ,・・・R
N-1 のそれぞれについて、受信信号の最大値と2番目に
大きい値を示す主力波2波と相関系列の絶対値との割合
Ti(i=0,1,・・・N−1)を算出する。主力波
2波と相関系列の絶対値との割合Tiの算出は、図12
に示すように、循環切換器106によってN組に振り分
けられる前の相関系列r(n)(n=0,1,・・・F
N−1)が最大となる点nmax 1202の前後1シンボ
ルの範囲1203で、相関系列の絶対値1201の中の
主力波である最大値rmax,i 1204と2番目に大きい
値r2nd,i 1205とを検索し、これら最大値と2番目
に大きい値とが存在する位置をそれぞれmmax,i 、m
2nd,i とする。つぎに、mmax,i を中心に前後(M−
1)シンボルの範囲1206で相関系列の絶対値120
1の和Si を算出し、さらに、主力波2波1204、1
205の和(rmax,i +r2nd,i )と相関系列の絶対値
1201の和Si との比Ti=(rmax,i +r2nd,i
/Si を算出する。N組の各数値系列R0 ,R1 ,・・
・RN-1 のそれぞれについて算出された主力波2波と相
関系列の絶対値1201との割合Tiの最大値を最大値
検索器109で検索し、Tiが最大となるiをi0
し、nをNで割った余りがi0 となるサンプル点をシン
ボルタイミングとし、i0 を出力端子110から出力す
る。
The N sets of numerical sequences R 0 , R 1 ,..., R N-1 stored in the N sets × F-stage two-dimensional register 107 are N
The data is read out to N sets of arithmetic circuits 108 connected to each set of the set × F-stage two-dimensional registers 107. N sets of arithmetic circuits 108 read out from the two-dimensional register 107 and have N sets of numerical sequences R 0 , R 1 ,.
For each of N-1 , the ratio Ti (i = 0, 1,... N-1) between the maximum value of the received signal and the two main force waves indicating the second largest value and the absolute value of the correlation sequence is calculated. I do. The calculation of the ratio Ti between the two main waves and the absolute value of the correlation sequence is shown in FIG.
, The correlation sequence r (n) (n = 0, 1,... F) before being allocated to N sets by the circulation switch 106.
In a range 1203 of one symbol before and after the point n max 1202 at which N−1) becomes the maximum, the maximum value r max, i 1204 which is the main wave in the absolute value 1201 of the correlation sequence and the second largest value r 2nd, i 1205, and the positions where these maximum value and the second largest value exist are defined as m max, i , m
2nd, i . Next, before and after centering on m max, i (M−
1) The absolute value 120 of the correlation sequence in the symbol range 1206
Calculating a first sum S i, further major wave 2 wave 1204,1
The ratio Ti = (r max, i + r 2nd, i ) between the sum (r max, i + r 2nd, i ) of 205 and the sum S i of the absolute value 1201 of the correlation sequence
/ Si is calculated. N sets of each numerical series R 0 , R 1 ,.
The maximum value search unit 109 searches for the maximum value of the ratio Ti between the two main waves calculated for each of R N-1 and the absolute value 1201 of the correlation sequence, and sets i at which Ti is maximum to i 0 , The sampling point at which the remainder of dividing n by N becomes i 0 is set as the symbol timing, and i 0 is output from the output terminal 110.

【0026】図5および図6に示す相関系列の図は、上
述したように1フレームあたりのシンボル数F=192
シンボル、既知シンボルの長さM=16シンボルで、1
シンボルあたりN=4回オーバーサンプリングされてお
り、遅延波と直接波とで振幅、位相が同じで、遅延時間
が1シンボルの主力波2波のモデルの例である。循環切
換器106により、N=4の周期で、4つの相関系列R
0 〜R3 に振り分けると、図13〜16の黒丸の点に示
す振り分けとなる。なお、図中の破線は、振り分けられ
る前の図6に示す相関系列である。図13に示す相関系
列のR0 は、nをN=4で割った余りが0になる |r
(n)| を抽出した相関系列(以下同様)である。
The diagrams of the correlation sequences shown in FIGS. 5 and 6 show the number of symbols per frame F = 192 as described above.
Symbol, known symbol length M = 16 symbols, 1
This is an example of a model of two main force waves each of which is oversampled N = 4 times per symbol, has the same amplitude and phase for the delayed wave and the direct wave, and has a delay time of one symbol. By the circulation switch 106, four correlation sequences R in N = 4 periods
If distributed to 0 to R 3, the distribution shown by the black point in FIG. 13-16. The broken line in the figure is the correlation sequence shown in FIG. 6 before being sorted. R 0 of the correlation sequence shown in FIG. 13 has a remainder of 0 obtained by dividing n by N = 4.
(N) | is a correlation sequence extracted (the same applies hereinafter).

【0027】振り分けられる前の相関系列が最大値を与
える点nmax 1202の前後1シンボルの範囲1203
内で、図13〜16に示す相関系列の主力波2波を検索
すると、図中に示すように図13の相関系列R0 につい
てはrmax,0 =57.36(mmax,0 =60)、r
2nd,0 =7.56(m2nd,0 =64)、図14の相関系
列R1 についてはrmax,1 =51.65(mmax,1 =6
1)、r2nd,1 =19.10(m2nd,1 =57)、図1
5の相関系列R2 についてはrmax,2 =36.62(m
max,2 =62)、r2nd,2 =35.78(m2nd,2 =5
8)、図16の相関系列R3 についてはrmax,3 =4
9.82(mmax,3 =59)、r2nd,3 =18.21
(m2nd,3 =63)となる。
A range 1203 of one symbol before and after a point n max 1202 where the correlation sequence before being sorted has the maximum value.
When two main waves of the correlation sequence shown in FIGS. 13 to 16 are searched, r max, 0 = 57.36 (m max, 0 = 60) for the correlation sequence R 0 shown in FIG. ), R
2nd, 0 = 7.56 (m 2nd, 0 = 64), and r max, 1 = 51.65 (m max, 1 = 6 for the correlation sequence R 1 in FIG. 14)
1), r 2nd, 1 = 19.10 (m 2nd, 1 = 57), FIG.
For the correlation sequence R 2 of 5, r max, 2 = 36.62 (m
max, 2 = 62), r 2nd, 2 = 35.78 (m 2nd, 2 = 5)
8), r max, 3 = 4 for correlation sequence R 3 in FIG. 16
9.82 (m max, 3 = 59), r 2nd, 3 = 18.21
(M 2nd, 3 = 63).

【0028】つぎに、mmax,i (i=0,1,2,3)
を中心に、前後M−1=9シンボルの範囲内で図13〜
16に示す相関系列の和を算出すると、図13の相関系
列R0 についてはS0 =575.5、図14の相関系列
1 についてはS1 =577.6、図15の相関系列R
2 についてはS2 =572.9、図16の相関系列R3
についてはS3 =583.0となる。したがって、主力
波2波の割合Tiは、図17に示すように、図13の相
関系列R0 についてはT0 =0.1128、図14の相
関系列R1 についてはT1 =0.1225、図15の相
関系列R2 についてはT2 =0.1264、図16の相
関系列R3 についてはT3 =0.1167となり、T2
が最大となるため、nを4で割った余りが2となるサン
プル点(n=2,6,10,・・・)をシンボルタイミ
ングとし、2を出力端子110から出力する。
Next, m max, i (i = 0,1,2,3)
13 to 13 in the range of M-1 = 9 symbols before and after
When the sum of the correlation sequences shown in FIG. 16 is calculated, S 0 = 575.5 for the correlation sequence R 0 in FIG. 13, S 1 = 577.6 for the correlation sequence R 1 in FIG. 14, and the correlation sequence R in FIG.
2 , S 2 = 572.9, and the correlation sequence R 3 in FIG.
Is S 3 = 583.0. Therefore, as shown in FIG. 17, the ratio Ti of the two main waves is T 0 = 0.1128 for the correlation sequence R 0 in FIG. 13, T 1 = 0.1225 for the correlation sequence R 1 in FIG. T 2 = .1264 for correlation sequence R 2 in FIG. 15, T 3 = 0.1167 next for correlation sequence R 3 in FIG. 16, T 2
Is the maximum, so the sampling point (n = 2, 6, 10,...) Where the remainder of dividing n by 4 becomes 2 is set as the symbol timing, and 2 is output from the output terminal 110.

【0029】最大補償遅延時間が1シンボルで、主力波
2波を除く非主力波を除去するビタビ等化器を使用した
場合、上述のように、主力波2波の割合Tiの値が大き
いほど非主力波の割合が小さいため、非主力波成分除去
の影響が少なくなる。従来技術によって検出した場合の
シンボルタイミングはn=0,4,8・・・であり、4
で割った余りが0なので、T0 を選択した場合のシンボ
ルタイミングと同じである。一方、本発明のシンボルタ
イミング検出方法による上述の例では、T0 より大きい
2 (最大値)を選択したため、従来技術のシンボルタ
イミング検出方法で検出したシンボルタイミングよりも
非主力波成分除去の影響が少ないシンボルタイミングを
検出している。本発明のシンボルタイミング検出方法に
より、主力波2波の割合が最も大きくなるシンボルタイ
ミングが検出され、主力波2波以外の非主力波成分を除
去するビタビ等化器を使用する場合、非主力波成分の除
去の影響が最も少なくなる。
When a maximum compensation delay time is one symbol and a Viterbi equalizer that removes non-main wave except the two main wave is used, as described above, as the ratio Ti of the two main wave becomes larger, Since the ratio of the non-mainstream wave is small, the influence of the removal of the non-mainstream wave component is reduced. The symbol timings detected by the prior art are n = 0, 4, 8,.
The remainder when divided by 0 is the same as the symbol timing when T 0 is selected. On the other hand, in the above-described example according to the symbol timing detection method of the present invention, since T 2 (maximum value) larger than T 0 is selected, the influence of removal of the non-mainstream wave component is smaller than that of the symbol timing detected by the conventional symbol timing detection method. Are detected with less symbol timing. According to the symbol timing detection method of the present invention, the symbol timing at which the ratio of the two main wave components is maximized is detected, and when using a Viterbi equalizer that removes non-main component wave components other than the two main force components, The effect of component removal is minimal.

【0030】[0030]

【発明の効果】本発明によれば、相関系列の主力波2波
の割合が最大となるようにシンボルタイミングを決定
し、非主力波成分を除去したことの影響が少なくなるよ
うにしたシンボルタイミング検出方法を提供することが
できる。
According to the present invention, the symbol timing is determined so that the ratio of the two main waves in the correlation sequence is maximized, and the influence of the removal of the non-main wave component is reduced. A detection method can be provided.

【図面の簡単な説明】[Brief description of the drawings]

【図1】本発明のシンボルタイミング検出方法によりシ
ンボルタイミングを検出する基本部分のブロック図。
FIG. 1 is a block diagram of a basic part for detecting a symbol timing by a symbol timing detection method according to the present invention.

【図2】本発明のシンボルタイミング検出方法を使用し
たシンボルタイミング検出回路のブロック図。
FIG. 2 is a block diagram of a symbol timing detection circuit using the symbol timing detection method of the present invention.

【図3】相関演算回路の詳細ブロック図。FIG. 3 is a detailed block diagram of a correlation operation circuit.

【図4】従来技術によるシンボルタイミング検出方法を
使用したシンボルタイミング検出回路のブロック図。
FIG. 4 is a block diagram of a symbol timing detection circuit using a conventional symbol timing detection method.

【図5】受信信号の相関系列(主力波2波モデル)。FIG. 5 shows a correlation sequence of a received signal (two main wave models).

【図6】図5に示す受信信号の相関系列の最大値付近の
拡大図。
FIG. 6 is an enlarged view of the vicinity of the maximum value of the correlation sequence of the received signal shown in FIG. 5;

【図7】伝送信号の1フレームのシンボル構成を示す
図。
FIG. 7 is a diagram showing a symbol configuration of one frame of a transmission signal.

【図8】受信信号の理想的な相関系列(主力波2波モデ
ル:遅延1シンボル)。
FIG. 8 shows an ideal correlation sequence of a received signal (two main wave models: one symbol delay).

【図9】サンプルタイミングによる相関系列の違いを説
明する図(主力波2波モデル:遅延1シンボル)。
FIG. 9 is a view for explaining a difference in a correlation sequence depending on a sample timing (a main wave two-wave model: one delay symbol).

【図10】最尤系列推定器の概念を仮想回路としたブロ
ック図。
FIG. 10 is a block diagram showing a concept of a maximum likelihood sequence estimator as a virtual circuit.

【図11】等化器の最大補償遅延時間の制限および非主
力波の除去を説明する図。
FIG. 11 is a view for explaining the limitation of the maximum compensation delay time of the equalizer and the removal of non-mainstream waves.

【図12】主力波2波の割合算出を説明する図。FIG. 12 is a view for explaining the calculation of the ratio of two main force waves.

【図13】図6の相関系列を循環切換器で振り分けた相
関系列R0 を示す図。
FIG. 13 is a diagram showing a correlation sequence R 0 obtained by allocating the correlation sequence of FIG. 6 by a circulation switch.

【図14】図6の相関系列を循環切換器で振り分けた相
関系列R1 を示す図。
FIG. 14 is a diagram showing a correlation sequence R 1 obtained by allocating the correlation sequence of FIG. 6 by a circulation switch.

【図15】図6の相関系列を循環切換器で振り分けた相
関系列R2 を示す図。
FIG. 15 is a diagram showing a correlation sequence R 2 obtained by allocating the correlation sequence of FIG. 6 by a circulation switch.

【図16】図6の相関系列を循環切換器で振り分けた相
関系列R3 を示す図。
FIG. 16 is a diagram showing a correlation sequence R 3 obtained by allocating the correlation sequence of FIG. 6 by a circulation switch.

【図17】図13〜16の相関系列R0 、R1 、R2
3 について算出した主力波2波の割合Tiを示す図。
FIG. 17 shows correlation sequences R 0 , R 1 , R 2 ,
Shows the percentage Ti main wave 2 wave calculated for R 3.

【符号の説明】[Explanation of symbols]

101、201、301、401…入力端子、102、
407…相関演算回路、103、302…遅延素子、1
04、303…乗算回路、105、304…加算回路、
106…循環切換器、107…2次元レジスタ、108
…演算回路、109、410…最大値検索器、110、
305、411…出力端子、111、408…絶対値演
算回路、202、402…同期検波器、205、405
…低域通過フィルタ、206、406…A/D変換器、
409…レジスタ。
101, 201, 301, 401 ... input terminal, 102,
407: correlation operation circuit, 103, 302 ... delay element, 1
04, 303 ... multiplication circuit, 105, 304 ... addition circuit,
106: circulation switching device, 107: two-dimensional register, 108
... arithmetic circuit, 109, 410 ... maximum value searcher, 110,
305, 411: output terminal, 111, 408: absolute value calculation circuit, 202, 402: synchronous detector, 205, 405
... low-pass filters, 206 and 406 ... A / D converters
409 ... register.

Claims (5)

【特許請求の範囲】[Claims] 【請求項1】 ディジタル通信に用いる受信機に使用す
るシンボルタイミング検出方法であって、 前記受信機側で既知なシンボル系列を、フレームのあら
かじめ定めた位置に配置した構成の送信信号を前記受信
機で受信して受信信号とし、 該受信信号と前記既知なシンボル系列との相関系列を算
出し、 該相関系列をシンボル周期で複数の数値系列に振り分
け、 該振り分けられた複数の数値系列の中から、数値の振幅
が最小となる数値系列を選択し、 該選択した数値系列が存在するタイミングをシンボルタ
イミングとすることを特徴とするシンボルタイミング検
出方法。
1. A method for detecting a symbol timing used in a receiver used for digital communication, comprising: transmitting a transmission signal having a configuration in which a symbol sequence known on the receiver side is arranged at a predetermined position in a frame. And calculates a correlation sequence between the received signal and the known symbol sequence. The correlation sequence is divided into a plurality of numerical sequences in a symbol cycle, and from among the plurality of distributed numerical sequences. A method for selecting a numerical sequence having the smallest numerical value amplitude, and using a timing at which the selected numerical sequence exists as a symbol timing.
【請求項2】 ディジタル通信に用いる受信機に使用す
るシンボルタイミング検出方法であって、 複数のシンボルからなる各フレームの先頭に長さMシン
ボルの既知なシンボル系列を配置したフレーム構成の送
信信号を前記受信機で受信して受信信号とし、 該受信信号を同期検波してベースバンド信号としてから
1シンボルあたりN回(Nは自然数)オーバーサンプリ
ングしてディジタル信号とし、 該ディジタル信号の前記既知なシンボル系列との相関を
演算して相関系列を算出し、 該相関系列の絶対値からシンボルタイミングを検出する
シンボルタイミング検出方法において、 前記相関系列の絶対値を1シンボルあたりのサンプル数
Nの周期で、N組の数値系列に振り分け、 該振り分けたN組の数値系列のそれぞれについて絶対値
が大きい2つの絶対値が前記相関系列の絶対値に占める
割合を算出し、 該割合が最も大きくなる相関系列を検索し、 検索した相関系列が存在するタイミングをシンボルタイ
ミングとすることを特徴とするシンボルタイミング検出
方法。
2. A symbol timing detection method used in a receiver used for digital communication, comprising: a transmission signal having a frame configuration in which a known symbol sequence having a length of M symbols is arranged at the head of each frame including a plurality of symbols. The received signal is received by the receiver to obtain a received signal. The received signal is synchronously detected to be a baseband signal, and then oversampled N times (N is a natural number) per symbol to obtain a digital signal. The known symbol of the digital signal In a symbol timing detection method for calculating a correlation sequence by calculating a correlation with a sequence and detecting a symbol timing from an absolute value of the correlation sequence, the absolute value of the correlation sequence is calculated at a period of the number N of samples per symbol. Sorted into N sets of numerical series, and the absolute value of each of the sorted N sets of numerical series is large. A symbol sequence for calculating a ratio of two absolute values to the absolute value of the correlation sequence, searching for a correlation sequence having the largest ratio, and setting a timing at which the searched correlation sequence exists as a symbol timing. Timing detection method.
【請求項3】 請求項1記載のシンボルタイミング検出
方法において、 振り分けたN組の数値系列のそれぞれについて絶対値が
大きい2つの絶対値が相関系列の絶対値に占める割合の
算出は、 前記振り分けたN組の数値系列のそれぞれについて絶対
値の最大値を算出し、 該絶対値の最大値を与えるサンプル点から前後1シンボ
ルの範囲内で、絶対値の最大値および2番目に大きい値
を検索し、 該絶対値の最大値および2番目に大きい値が前記相関系
列の絶対値に占める割合を算出し、 N組の数値系列のそれぞれから前記割合が最も大きくな
る相関系列を検索し、 検索した相関系列が存在するタイミングをシンボルタイ
ミングとすることを特徴とするシンボルタイミング検出
方法。
3. The symbol timing detection method according to claim 1, wherein the ratio of two absolute values having a large absolute value to the absolute value of the correlation sequence is calculated for each of the allocated N sets of numerical sequences. The maximum value of the absolute value is calculated for each of the N sets of numerical sequences, and the maximum value of the absolute value and the second largest value are searched within one symbol before and after the sample point that gives the maximum value of the absolute value. Calculating the proportion of the absolute value of the maximum value and the second largest value of the absolute value to the absolute value of the correlation sequence, searching for the correlation sequence having the largest ratio from each of the N sets of numerical sequences, A symbol timing detection method, wherein timing at which a sequence exists is used as symbol timing.
【請求項4】 請求項2記載のシンボルタイミング検出
方法において、 絶対値の最大値および2番目に大きい値が相関系列の絶
対値に占める割合の算出は、 前記絶対値の最大値を与えるサンプル点から前後(M−
1)シンボルの範囲内で、絶対値の最大値および2番目
に大きい値が相関系列の絶対値に占める割合を算出し、 N組の数値系列のそれぞれから前記割合が最も大きくな
る相関系列を検索し、 検索した相関系列が存在するタイミングをシンボルタイ
ミングとすることを特徴とするシンボルタイミング検出
方法。
4. The symbol timing detecting method according to claim 2, wherein calculating the maximum value of the absolute value and the ratio of the second largest value to the absolute value of the correlation sequence comprises: Before and after (M-
1) Calculate the ratio of the maximum absolute value and the second largest value to the absolute value of the correlation sequence within the range of the symbol, and search the N sets of numerical sequences for the correlation sequence with the highest ratio. A timing at which the searched correlation sequence exists is used as a symbol timing.
【請求項5】 ディジタル通信に用いる受信機に使用す
るシンボルタイミング検出方法であって、 所定のビット数のデータの複数のシンボルからなる各フ
レームの先頭に長さMシンボルの既知なシンボル系列を
配置したフレーム構成の送信信号を受信して受信信号と
し、 該受信信号を同期検波器で検波してベースバンド信号と
し、 該ベースバンド信号をA/D変換器で1シンボルあたり
N回(Nは自然数)オーバーサンプリングしてディジタ
ル信号とし、 該ディジタル信号を相関演算回路で前記既知なシンボル
系列との相関を演算して相関系列を算出し、 該相関系列を絶対値演算回路で演算して相関系列の絶対
値を算出し、 該相関系列の絶対値の少なくとも1フレームをレジスタ
に記憶し、 該レジスタから前記相関系列の絶対値の少なくとも1フ
レームを読み出して最大値検索器で最大値を検索し、該
最大値を与えるサンプル点をシンボルタイミングとする
シンボルタイミング検出方法において、 前記絶対値演算回路で算出した前記相関系列の絶対値
を、循環切換器で順次に、循環して1シンボルあたりの
サンプル数Nの周期で、N組の数値系列に振り分け、 該振り分けたN組の数値系列をN組×F段の2次元レジ
スタに記憶し、 該N組×F段の2次元レジスタからそれぞれの組に接続
されているN組の演算回路に読み出して、前記N組の数
値系列のそれぞれについて絶対値の最大値および2番目
に大きい値が前記相関系列の絶対値に占める割合を算出
し、 前記N組の数値系列のそれぞれの割合から最大値検索器
で前記割合が最も大きくなる相関系列を検索し、 検索した前記割合が最も大きくなる相関系列が存在する
タイミングをシンボルタイミングとすることを特徴とす
るシンボルタイミング検出方法。
5. A method for detecting a symbol timing used in a receiver used for digital communication, wherein a known symbol sequence having a length of M symbols is arranged at the head of each frame composed of a plurality of symbols of data having a predetermined number of bits. A transmission signal having a frame configuration is received and converted into a reception signal. The reception signal is detected by a synchronous detector to generate a baseband signal. The baseband signal is N times per symbol by an A / D converter (N is a natural number). ) Oversampling into a digital signal, calculating the correlation of the digital signal with the known symbol sequence in a correlation operation circuit to calculate a correlation sequence, and calculating the correlation sequence in an absolute value operation circuit to calculate the correlation sequence. Calculating an absolute value; storing at least one frame of the absolute value of the correlation sequence in a register; Also reads one frame, searches for a maximum value with a maximum value searcher, and uses a sample point giving the maximum value as a symbol timing, wherein the absolute value of the correlation sequence calculated by the absolute value calculation circuit is calculated. Sequentially circulated by the circulation switcher and distributed to N sets of numerical sequences at a period of N samples per symbol, and the allocated N sets of numerical sequences are stored in N sets × F-stage two-dimensional registers. The N sets of arithmetic circuits connected to each set are read out from the N sets × F stages of two-dimensional registers, and the maximum value and the second largest value of the absolute value are obtained for each of the N sets of numerical sequences. Calculates the percentage of the absolute value of the correlation sequence in the absolute value, and searches the maximum value searcher for the correlation sequence in which the ratio is the largest from the respective ratios of the N sets of numerical sequences. A symbol timing detection method, wherein a timing at which a correlation sequence having the largest combination exists is defined as a symbol timing.
JP05206099A 1999-03-01 1999-03-01 Symbol timing detection method Expired - Fee Related JP3619384B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP05206099A JP3619384B2 (en) 1999-03-01 1999-03-01 Symbol timing detection method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP05206099A JP3619384B2 (en) 1999-03-01 1999-03-01 Symbol timing detection method

Publications (2)

Publication Number Publication Date
JP2000252969A true JP2000252969A (en) 2000-09-14
JP3619384B2 JP3619384B2 (en) 2005-02-09

Family

ID=12904282

Family Applications (1)

Application Number Title Priority Date Filing Date
JP05206099A Expired - Fee Related JP3619384B2 (en) 1999-03-01 1999-03-01 Symbol timing detection method

Country Status (1)

Country Link
JP (1) JP3619384B2 (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2007537661A (en) * 2004-05-12 2007-12-20 トムソン ライセンシング Correction of symbol timing ambiguity
JP2008236206A (en) * 2007-03-19 2008-10-02 New Japan Radio Co Ltd Bit synchronization circuit

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2007537661A (en) * 2004-05-12 2007-12-20 トムソン ライセンシング Correction of symbol timing ambiguity
JP2008236206A (en) * 2007-03-19 2008-10-02 New Japan Radio Co Ltd Bit synchronization circuit

Also Published As

Publication number Publication date
JP3619384B2 (en) 2005-02-09

Similar Documents

Publication Publication Date Title
US7966360B2 (en) Finite impulse response filter and digital signal receiving apparatus
CN111694027B (en) Ultra-large dynamic spread spectrum signal capturing method and device
JP2000504514A (en) Digital architecture for reproducing NRZ / NRZI data
TWI299948B (en) Apparatus and method of searching for known sequences
JPH06197083A (en) Digital data demodulator
CA2035774C (en) Frame phase estimation method and circuit
JP2000504539A (en) Method and receiver for calculating impulse response
US7206365B2 (en) Decision sequence generating method and associated receiver with a decision feedback equalizer
KR100534592B1 (en) Apparatus and method for digital communication receiver
JP3619384B2 (en) Symbol timing detection method
JP2001211102A (en) Rake receiver
KR100395497B1 (en) Method for Canceling Interference for Parallel Type for CDMA Receiver
JPH10200503A (en) Adaptive equalizing circuit for spread spectrum communication
US6928105B2 (en) Vector tree correlator for variable spreading rates
EP1458097A1 (en) Arbitrary sampling rate conversion
US6304591B1 (en) Match filter architecture based upon parallel I/O
JPH118575A (en) Maximum likelihood system estimate receiver
JP2801782B2 (en) Frame phase estimation method and circuit
JP3244722B2 (en) Maximum likelihood sequence estimator
KR20050021491A (en) Fsk receiver having a variable threshold slicer stage and corresponding method
JP2000091965A (en) Equalizer and equalizing method
JPH0897810A (en) Digital transmission system able to be synchronized in initializing sequence
KR100289404B1 (en) Apparatus and method for reducing pattern jitter by using quasi locally symmetric wave signal
JP2002084267A (en) Timing detection device and timing detection method
JP2001103044A (en) Method for synchronizing symbol timing

Legal Events

Date Code Title Description
A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20040715

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20040720

A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20040916

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

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20041112

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

Year of fee payment: 3

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

Free format text: PAYMENT UNTIL: 20081119

Year of fee payment: 4

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

Free format text: PAYMENT UNTIL: 20091119

Year of fee payment: 5

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

Free format text: PAYMENT UNTIL: 20091119

Year of fee payment: 5

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

Free format text: PAYMENT UNTIL: 20101119

Year of fee payment: 6

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

Free format text: PAYMENT UNTIL: 20111119

Year of fee payment: 7

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

Free format text: PAYMENT UNTIL: 20111119

Year of fee payment: 7

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

Free format text: PAYMENT UNTIL: 20121119

Year of fee payment: 8

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

Free format text: PAYMENT UNTIL: 20121119

Year of fee payment: 8

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

Free format text: PAYMENT UNTIL: 20131119

Year of fee payment: 9

LAPS Cancellation because of no payment of annual fees