JP2010016617A - Delay profile estimation device and delay profile estimating method - Google Patents

Delay profile estimation device and delay profile estimating method Download PDF

Info

Publication number
JP2010016617A
JP2010016617A JP2008174584A JP2008174584A JP2010016617A JP 2010016617 A JP2010016617 A JP 2010016617A JP 2008174584 A JP2008174584 A JP 2008174584A JP 2008174584 A JP2008174584 A JP 2008174584A JP 2010016617 A JP2010016617 A JP 2010016617A
Authority
JP
Japan
Prior art keywords
sequence
correlation
delay profile
calculating
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
JP2008174584A
Other languages
Japanese (ja)
Other versions
JP5100536B2 (en
Inventor
Daisuke Shinpo
大介 新保
Naotoshi Maeda
尚利 前田
Eiji Arita
栄治 有田
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.)
Mitsubishi Electric Corp
Original Assignee
Mitsubishi Electric Corp
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 Mitsubishi Electric Corp filed Critical Mitsubishi Electric Corp
Priority to JP2008174584A priority Critical patent/JP5100536B2/en
Priority to CN201210447494.9A priority patent/CN102970097B/en
Priority to CN2009101299556A priority patent/CN101621310B/en
Publication of JP2010016617A publication Critical patent/JP2010016617A/en
Application granted granted Critical
Publication of JP5100536B2 publication Critical patent/JP5100536B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Abstract

<P>PROBLEM TO BE SOLVED: To estimate a delay profile highly precisely by using an autocorrelation characteristic of a pseudo random (PN) sequence in a digital broadcast receiver. <P>SOLUTION: A signal transmitted by a transmitter with a transmission symbol (St) composed of a frame header (Hf) obtained by adding portions of the head (Lpre) and end (Lpost) of a known pseudo random (PN) sequence before and after a PN sequence, and an effective symbol (Se) including information to be transmitted as a transmission unit is received through a transmission line, correlation sequences between two PN sequences (Tf and Tb) and a received signal (Rs) respectively constituting a portion of the frame header (Hf) are calculated, instantaneous power of correlation values to which the respective correlation sequences correspond is compared, and a correlation value with small instantaneous power is outputted as an estimation result of a delay profile. <P>COPYRIGHT: (C)2010,JPO&INPIT

Description

この発明は、伝送路の遅延プロファイルを推定する遅延プロファイル推定装置及び方法に関するものである。   The present invention relates to a delay profile estimation apparatus and method for estimating a delay profile of a transmission line.

地上デジタル放送では、送信機から出力された電波が、建物等の障害物による反射、回折、散乱を受けるため、受信信号に歪が生じる。受信機では、信頼性のある受信性能を実現するために、受信信号から伝送路の遅延プロファイルを推定し、この推定結果を用いて受信信号の歪を補正する必要がある。伝送路の遅延プロファイルを推定する方法として、送信機で擬似ランダム(PN)系列を既知信号として挿入する方法がある。   In terrestrial digital broadcasting, radio waves output from a transmitter are reflected, diffracted, and scattered by an obstacle such as a building, so that a received signal is distorted. In order to realize reliable reception performance, the receiver needs to estimate the delay profile of the transmission path from the received signal and correct the distortion of the received signal using this estimation result. As a method for estimating a delay profile of a transmission line, there is a method of inserting a pseudo random (PN) sequence as a known signal at a transmitter.

中国の地上デジタル放送方式では、送信機は既知の擬似ランダム(PN)系列の先頭から規定長の系列を該PN系列の後に付加し、該PN系列の末尾から規定長の系列を該PN系列の前に付加した構造を持つフレームヘッダを生成し、このフレームヘッダと、送信すべき情報を含む有効シンボルで構成される伝送シンボルを伝送単位として送信している。   In the Chinese terrestrial digital broadcasting system, the transmitter adds a specified length sequence from the beginning of a known pseudo-random (PN) sequence to the end of the PN sequence, and adds a specified length sequence from the end of the PN sequence to the PN sequence. A frame header having a previously added structure is generated, and a transmission symbol composed of the frame header and a valid symbol including information to be transmitted is transmitted as a transmission unit.

受信機は、受信信号を所定のサンプリング周波数で標本化して得られる標本化系列と自己発生させた既知のPN系列との相関をサンプル毎に計算することで得られる相関系列から遅延プロファイルを推定する。受信信号に含まれるPN系列と受信機内で自己発生させたPN系列とのパターンが完全に一致するとき、相関結果は鋭いピークとなる。このピークの値は到来波の受信レベルに比例した値となる。しかし、受信信号に含まれるPN系列の一部と、自己発生させたPN系列の一部のパターンが一致するとき、不要な相関のピークが生じる。この不要な相関のピークは、到来波に応じた相関のピークを中心にして既知のPN系列長だけ離れた位置に存在し、不要な相関のピークの値は到来波の受信レベルから一意に決まる。   The receiver estimates the delay profile from the correlation sequence obtained by calculating the correlation between the sampling sequence obtained by sampling the received signal at a predetermined sampling frequency and the self-generated known PN sequence for each sample. . When the patterns of the PN sequence included in the received signal and the PN sequence self-generated in the receiver completely match, the correlation result has a sharp peak. This peak value is proportional to the incoming wave reception level. However, when a part of the PN sequence included in the received signal matches the pattern of a part of the self-generated PN sequence, an unnecessary correlation peak occurs. The unnecessary correlation peak exists at a position separated by a known PN sequence length around the correlation peak corresponding to the incoming wave, and the value of the unnecessary correlation peak is uniquely determined from the incoming wave reception level. .

受信信号の標本化系列と自己発生させた既知のPN系列との相関系列に含まれる不要な相関のピークを取り除くことで遅延プロファイルが得られる。非特許文献1に、相関系列から最も大きな相関値を検出し、この相関値からPN系列長だけ離れた位置に存在する相関結果から、不要な相関のピークの値を減算する方法が掲載されている。   A delay profile is obtained by removing an unnecessary correlation peak included in a correlation sequence between a sampling sequence of a received signal and a self-generated known PN sequence. Non-Patent Document 1 describes a method of detecting the largest correlation value from a correlation sequence and subtracting the value of an unnecessary correlation peak from a correlation result existing at a position separated from the correlation value by a PN sequence length. Yes.

Guanghui Liu, “ITD−DFE Based Channel Estimation and Equalization in TDS−OFDM Receivers”, IEEE Transactions on Consumer Electronics, Vol.53, No.2, pp.304−309 (第305頁)Guanghui Liu, “ITD-DFE Based Channel Estimation and Evaluation in TDS-OFDM Receivers”, IEEE Transactions on Consumer Electronics, Vol. 53, no. 2, pp. 304-309 (Page 305)

しかしながら、マルチパス伝送路の遅延プロファイルを従来手法で推定する場合、受信信号の標本化系列と自己発生させた既知のPN系列との相関系列から到来波に応じた相関のピークを探索し、この相関のピークから既知のPN系列の長さだけ離れた位置の相関結果から不要な相関のピークの値を減算する処理を到来波の数だけ繰り返す必要があり、計算量が多くなるという問題がある。   However, when estimating the delay profile of the multipath transmission path by the conventional method, the correlation peak corresponding to the incoming wave is searched from the correlation sequence of the sampling sequence of the received signal and the self-generated known PN sequence. It is necessary to repeat the process of subtracting the value of the unnecessary correlation peak from the correlation result at a position separated from the correlation peak by the length of the known PN sequence as many times as the number of incoming waves, which increases the amount of calculation. .

また、マルチパス伝送路において、主波の到来時刻を基準にして遅延時間差が既知のPN系列の長さと等しい先行波、及び遅延波が存在する場合、従来手法で検出する主波に応じた相関のピーク値に、先行波、及び遅延波による不要な相関のピーク値が加算されるため、遅延プロファイルの推定精度が劣化するという問題がある。   Also, in the multipath transmission line, if there is a preceding wave and a delayed wave whose delay time difference is equal to the known PN sequence length with respect to the arrival time of the main wave, the correlation according to the main wave detected by the conventional method Since the peak value of unnecessary correlation due to the preceding wave and the delayed wave is added to the peak value of, the estimation accuracy of the delay profile is degraded.

本発明による遅延プロファイル推定装置は、
送信機で、既知の擬似ランダム(PN)系列における先頭の規定長の系列を該PN系列の後に付加し、該PN系列における末尾の規定長の系列を該PN系列の前に付加した構造を持つフレームヘッダと、送信すべき情報を含む有効シンボルとで構成される伝送シンボルを伝送単位として送信された信号を、伝送路を介して受信し、受信信号を所定のサンプリング周波数で標本化して得られる標本化系列から伝送路の遅延プロファイルを推定する遅延プロファイル推定装置において、
それぞれフレームヘッダの一部を構成する2つの異なるPN系列を生成する第1及び第2のPN系列生成手段と、
前記第1のPN系列生成手段によって生成されたPN系列と受信信号の標本化系列との相関を計算する第1の相関計算手段と、
前記第2のPN系列生成手段によって生成されたPN系列と受信信号の標本化系列との相関を計算する第2の相関計算手段と、
前記第1及び第2の相関計算手段で得られる2つの相関系列から遅延プロファイルを推定するプロファイル値計算手段を備え、
前記プロファイル値計算手段が、前記第1及び第2の相関計算手段から出力される2つの相関系列の互いに対応する相関値の瞬時電力を比較して瞬時電力の小さい相関値を選択し、選択した結果を遅延プロファイルとして出力する
ことを特徴とする。
The delay profile estimation apparatus according to the present invention includes:
The transmitter has a structure in which a first specified length sequence in a known pseudo-random (PN) sequence is added after the PN sequence, and a last specified length sequence in the PN sequence is added before the PN sequence. It is obtained by receiving a signal transmitted using a transmission symbol composed of a frame header and effective symbols including information to be transmitted as a transmission unit, and sampling the received signal at a predetermined sampling frequency. In a delay profile estimation device that estimates a delay profile of a transmission path from a sampling sequence,
First and second PN sequence generating means for generating two different PN sequences, each forming part of the frame header;
First correlation calculating means for calculating a correlation between the PN sequence generated by the first PN sequence generating means and the sampling sequence of the received signal;
Second correlation calculating means for calculating the correlation between the PN sequence generated by the second PN sequence generating means and the sampling sequence of the received signal;
Profile value calculation means for estimating a delay profile from two correlation sequences obtained by the first and second correlation calculation means;
The profile value calculation means selects a correlation value having a small instantaneous power by comparing the instantaneous powers of the correlation values corresponding to each other of the two correlation sequences output from the first and second correlation calculation means The result is output as a delay profile.

本発明によれば、繰返し処理を行う従来手法と比較して遅延プロファイルを推定するために必要な計算量が少なくなるという効果が得られる。   According to the present invention, it is possible to obtain an effect that the amount of calculation required for estimating the delay profile is reduced as compared with the conventional method in which iterative processing is performed.

実施の形態1.
図1はこの発明の実施の形態1に係わる遅延プロファイル推定法推定装置を示す。図示の遅延プロファイル推定装置は、第1のPN系列生成手段1と、第1の相関計算手段2と、第2のPN系列生成手段3と、第2の相関計算手段4と、遅延手段5と、第1の瞬時電力計算手段6と、第2の瞬時電力計算手段7と、最小電力探索手段8と、最小電力選択手段9を備える。
Embodiment 1 FIG.
FIG. 1 shows a delay profile estimation method estimation apparatus according to Embodiment 1 of the present invention. The illustrated delay profile estimation apparatus includes a first PN sequence generation unit 1, a first correlation calculation unit 2, a second PN sequence generation unit 3, a second correlation calculation unit 4, and a delay unit 5. The first instantaneous power calculating means 6, the second instantaneous power calculating means 7, the minimum power searching means 8, and the minimum power selecting means 9 are provided.

図1の遅延プロファイル推定装置には、受信信号Rsを所定のサンプリング周波数で標本化した系列が入力される。本実施の形態では、送信信号を非特許文献1で定められた形式の信号とする。図2(a)にそのような送信信号の構成を示す。送信信号は規定のPN系列で構成されるフレームヘッダHfと送信すべき情報を含む有効シンボルSeからなる伝送シンボルStを伝送単位とする信号である。フレームヘッダHfは、図2(b)に示すように、長さLmサンプルの規定のPN系列Tdの先頭Lpreサンプルを規定のPN系列Tdの後に付加し、規定のPN系列Tdの末尾Lpostサンプルを規定のPN系列Tdの前に付加した構造を持つ。   A sequence obtained by sampling the received signal Rs at a predetermined sampling frequency is input to the delay profile estimation apparatus of FIG. In the present embodiment, the transmission signal is a signal defined in Non-Patent Document 1. FIG. 2A shows the configuration of such a transmission signal. The transmission signal is a signal whose transmission unit is a transmission symbol St composed of a frame header Hf composed of a prescribed PN sequence and an effective symbol Se including information to be transmitted. As shown in FIG. 2B, the frame header Hf adds the first Lpre sample of the specified PN sequence Td of length Lm samples after the specified PN sequence Td, and adds the last Lpost sample of the specified PN sequence Td. It has a structure added before the prescribed PN sequence Td.

図1の遅延プロファイル推定装置における第1のPN系列生成手段1は、図3(a)に示すように、フレームヘッダHfの先頭LmサンプルのPN系列Tfを生成する。
第1の相関計算手段2は、第1のPN系列生成手段1で生成されたPN系列Tfと受信信号Rsの標本化系列との相関を計算することで相関系列Rf(k)を得る。
第2のPN系列生成手段3は、図4(a)に示すように、フレームヘッダHfの末尾LmサンプルのPN系列Tbを生成する。
第2の相関計算手段4は、第2のPN系列生成手段3で生成されたPN系列Tbと受信信号Rsの標本化系列との相関を計算することで相関系列Rb(k)を得る。
The first PN sequence generation means 1 in the delay profile estimation apparatus of FIG. 1 generates a PN sequence Tf of the first Lm samples of the frame header Hf, as shown in FIG.
The first correlation calculating means 2 obtains a correlation sequence Rf (k) by calculating the correlation between the PN sequence Tf generated by the first PN sequence generating means 1 and the sampled sequence of the received signal Rs.
The second PN sequence generation means 3 generates the PN sequence Tb of the last Lm samples of the frame header Hf as shown in FIG.
The second correlation calculation unit 4 obtains a correlation sequence Rb (k) by calculating the correlation between the PN sequence Tb generated by the second PN sequence generation unit 3 and the sampled sequence of the received signal Rs.

遅延手段5は、第1の相関計算手段2から出力される相関系列Rf(k)をLpre+Lpostサンプル間隔だけ遅延させて、遅延相関系列Rf(k−d)を出力する。
第1の瞬時電力計算手段6は、遅延手段5から出力される相関系列Rf(k−d)の各相関値の瞬時電力(相関系列Rf(k−d)の各kの値についての相関値の2乗に比例する値)を計算する。
第2の瞬時電力計算手段7は、第2の相関計算手段4から出力される相関系列Rb(k)の各相関値の瞬時電力(相関系列Rb(k)の各kの値についての相関値の2乗に比例する値)を計算する。
The delay unit 5 delays the correlation sequence Rf (k) output from the first correlation calculation unit 2 by the Lpre + Lpost sample interval, and outputs a delayed correlation sequence Rf (k−d).
The first instantaneous power calculation means 6 is a correlation value for each value of k of the correlation sequence Rf (k−d) of the correlation sequence Rf (k−d) output from the delay means 5. (Proportional to the square of).
The second instantaneous power calculation means 7 calculates the instantaneous power of each correlation value of the correlation sequence Rb (k) output from the second correlation calculation means 4 (correlation value for each k value of the correlation sequence Rb (k)). (Proportional to the square of).

最小電力探索手段8は、第1の瞬時電力計算手段6と第2の瞬時電力計算手段7から出力される瞬時電力のうち小さい方を探索する。
最小電力選択手段9は、最小電力探索手段8における探索結果をもとに、遅延手段5から出力される相関系列Rf(k−d)と第2の相関計算手段4から出力される相関系列Rb(k)から瞬時電力が小さい相関値(相関系列Rf(k−d)及びRb(k)の当該kの値についての相関値のうちの一方)を選択し、推定遅延プロファイル(の当該kの値についての値)として出力する。
The minimum power search means 8 searches for the smaller of the instantaneous powers output from the first instantaneous power calculation means 6 and the second instantaneous power calculation means 7.
The minimum power selection means 9 is based on the search result in the minimum power search means 8 and the correlation sequence Rf (k−d) output from the delay means 5 and the correlation sequence Rb output from the second correlation calculation means 4. From (k), a correlation value (one of correlation values for the k value of the correlation series Rf (k−d) and Rb (k)) is selected and the estimated delay profile Value).

上記のうち、遅延手段5と、第1の瞬時電力計算手段6と、第2の瞬時電力計算手段7と、最小電力探索手段8と、最小電力選択手段9とで、第1の相関手段2で得られる相関系列Rf(k)及び第2の相関計算手段4で得られる相関系列Rb(k)から遅延プロファイルを推定するプロファイル値計算手段10が構成される。   Among the above, the delay means 5, the first instantaneous power calculation means 6, the second instantaneous power calculation means 7, the minimum power search means 8, and the minimum power selection means 9 are the first correlation means 2. The profile value calculation means 10 for estimating the delay profile from the correlation sequence Rf (k) obtained by the above and the correlation series Rb (k) obtained by the second correlation calculation means 4 is configured.

以下、より詳細に説明する。まず、図1における第1の相関計算手段2で得られる相関系列Rf(k)について説明する。
図5(a)に第1の相関計算手段2から出力される相関系列を示す。ここで、伝送路はマルチパスが存在しないとし、フレームヘッダHfの構造は非特許文献1と同様にLm=255、Lpre=82、Lpost=83としている。
Rf(k)は、図3(b)に示すように、PN系列Tfの先頭のサンプルを基準にして受信信号Rsに含まれるフレームヘッダHfの先頭サンプルがkサンプルだけシフトしたときに、これらの系列の相関を計算した結果である。
This will be described in more detail below. First, the correlation sequence Rf (k) obtained by the first correlation calculation means 2 in FIG. 1 will be described.
FIG. 5A shows a correlation sequence output from the first correlation calculation means 2. Here, it is assumed that there is no multipath in the transmission path, and the structure of the frame header Hf is Lm = 255, Lpre = 82, and Lpost = 83, as in Non-Patent Document 1.
As shown in FIG. 3B, Rf (k) is obtained when the first sample of the frame header Hf included in the received signal Rs is shifted by k samples with reference to the first sample of the PN sequence Tf. It is the result of calculating the correlation of the series.

図3(b)においてk=0のとき、受信信号Rsに含まれるPN系列と第1のPN系列生成手段1によって生成されたPN系列Tfが完全に一致するため、図5(a)においてk=0に受信レベルに比例した振幅を持つ相関のピーク値が得られる。
また、図3(b)おいてk=255のとき、受信信号Rsに含まれるPN系列の一部と第1のPN系列生成手段によって生成されたPN系列Tfの一部が一致するため、図5(a)においてk=255に不要な相関のピークが生じる。
When k = 0 in FIG. 3B, the PN sequence included in the received signal Rs and the PN sequence Tf generated by the first PN sequence generating means 1 are completely coincident with each other. A correlation peak value having an amplitude proportional to the reception level is obtained at = 0.
3B, when k = 255, a part of the PN sequence included in the received signal Rs and the part of the PN sequence Tf generated by the first PN sequence generating unit match. In 5 (a), an unnecessary correlation peak occurs at k = 255.

2波モデル伝送路の場合の第1の相関計算手段2から出力される相関系列Rf(k)を図5(b)に示す。ここで、遅延波の遅延時間を200サンプル、D/Uを6[dB]と想定した。図において、k=0に主波に応じた相関のピーク、k=200に遅延波に応じた相関のピークが存在し、それぞれのピークに対して255離れた位置に不要なピークが存在する。   FIG. 5B shows the correlation sequence Rf (k) output from the first correlation calculation means 2 in the case of the two-wave model transmission path. Here, it is assumed that the delay time of the delay wave is 200 samples and D / U is 6 [dB]. In the figure, there is a correlation peak corresponding to the main wave at k = 0, a correlation peak according to the delayed wave at k = 200, and an unnecessary peak at a position 255 away from each peak.

次に、図1における第2の相関計算手段4で得られる相関系列Rb(k)について説明する。図6(a)に第2の相関計算手段4から出力される相関系列Rb(k)を示す。Rb(k)は、図4(b)に示すように、PN系列Tbの先頭のサンプルを基準にして、受信信号Rsに含まれるフレームヘッダHfの先頭サンプルがkサンプルだけシフトしたときに、これらの系列の相関を計算した結果である。図4(b)においてk=165のとき、受信信号Rsに含まれるPN系列と第2のPN系列生成手段3によって生成されたPN系列Tbが完全に一致するため、図6(a)においてk=165に受信レベルに比例した振幅を持つ相関のピーク値が得られる。また、図4(b)におけるk=−90のとき、受信信号Rsに含まれるPN系列の一部と第2のPN系列生成手段3によって生成されたPN系列Tbの一部が一致するため、図6(a)におけるk=−90に不要な相関のピークが生じる。   Next, the correlation sequence Rb (k) obtained by the second correlation calculation means 4 in FIG. 1 will be described. FIG. 6A shows the correlation sequence Rb (k) output from the second correlation calculation means 4. As shown in FIG. 4B, Rb (k) is obtained when the first sample of the frame header Hf included in the received signal Rs is shifted by k samples with reference to the first sample of the PN sequence Tb. It is the result of calculating the correlation of the series. In FIG. 4B, when k = 165, the PN sequence included in the received signal Rs and the PN sequence Tb generated by the second PN sequence generating means 3 completely match. = 165, a correlation peak value having an amplitude proportional to the reception level is obtained. In addition, when k = −90 in FIG. 4B, a part of the PN sequence included in the received signal Rs and a part of the PN sequence Tb generated by the second PN sequence generating means 3 match, An unnecessary correlation peak occurs at k = −90 in FIG.

2波モデル伝送路の場合の第2の相関計算手段4から出力される相関系列Rb(k)を図6(b)に示す。ここで、遅延波の遅延時間を200サンプル、D/Uを6[dB]と想定した。図において、k=165に主波に応じた相関のピーク、k=365に遅延波に応じた相関のピークが存在し、それぞれのピークに対して255離れた位置に不要なピークが存在する。   FIG. 6B shows the correlation sequence Rb (k) output from the second correlation calculation means 4 in the case of the two-wave model transmission path. Here, it is assumed that the delay time of the delay wave is 200 samples and D / U is 6 [dB]. In the figure, there is a correlation peak corresponding to the main wave at k = 165, a correlation peak according to the delayed wave at k = 365, and an unnecessary peak at a position 255 away from each peak.

図5(b)と図6(b)を比較すると、第2の相関計算手段4から出力される相関系列Rb(k)に含まれる主波に応じた相関のピークは、第1の相関計算手段2から出力される相関系列Rf(k)に含まれる主波に応じた相関のピークに比べて165サンプル間隔(即ち、(Lpre+Lpost)サンプル間隔)だけ遅れて出力される。また、遅延波に応じたピークについても同様のことが言える。そこで、図1の遅延プロファイル推定装置における遅延手段5は、2つの相関系列Rf(k)、Rb(k)に含まれる主波と遅延波のピークが存在する位置を互いに対応させる(同時に発生させる)ために、第1の相関計算手段2から出力される系列Rf(k)を165サンプル間隔だけ遅延させ、相関系列Rf(k−165)を得る。その結果、図7(a)及び(b)に示すように、相関系列Rf(k−d)=Rf(k−165)と相関系列Rb(k)のk=165とk=365の値がそれぞれ主波と遅延波に応じた相関のピークとなる。   Comparing FIG. 5B and FIG. 6B, the correlation peak corresponding to the main wave included in the correlation sequence Rb (k) output from the second correlation calculation means 4 is the first correlation calculation. Compared to the correlation peak corresponding to the main wave included in the correlation sequence Rf (k) output from the means 2, the output is delayed by 165 sample intervals (that is, (Lpre + Lpost) sample interval). The same applies to the peak corresponding to the delayed wave. Therefore, the delay means 5 in the delay profile estimation apparatus of FIG. 1 associates the positions where the main wave and the peak of the delay wave included in the two correlation sequences Rf (k) and Rb (k) exist with each other (generate them simultaneously). Therefore, the sequence Rf (k) output from the first correlation calculation means 2 is delayed by 165 sample intervals to obtain a correlation sequence Rf (k-165). As a result, as shown in FIGS. 7A and 7B, the correlation sequence Rf (k−d) = Rf (k−165) and the correlation sequence Rb (k) have values k = 165 and k = 365. The correlation peaks according to the main wave and the delayed wave, respectively.

図7(a)及び(b)に示すように、遅延手段から出力される相関系列Rf(k−165)と第2の相関計算手段4から出力される相関系列Rb(k)に含まれる主波、及び遅延波に応じた相関のピークの存在する位置は対応しいているが、2つの相関系列間で不要な相関のピークが存在する位置は対応していない。つまり、一方の相関系列において不要な相関ピークが存在する位置において、もう一方の相関系列では、非常に小さな相関値となる。   As shown in FIGS. 7A and 7B, the correlation sequence Rf (k-165) output from the delay unit and the correlation sequence Rb (k) output from the second correlation calculation unit 4 are included. The position where the correlation peak corresponding to the wave and the delayed wave exists corresponds, but the position where the unnecessary correlation peak exists between the two correlation sequences does not correspond. That is, at a position where an unnecessary correlation peak exists in one correlation series, the correlation value is very small in the other correlation series.

図1の遅延プロファイル推定装置において、第1の瞬時電力計算手段6と第2の瞬時電力計算手段7は、それぞれ遅延手段5から出力される相関系列Rf(k−165)と第2の相関計算手段4から出力される相関系列Rb(k)の各相関値の瞬時電力((Rf(k−165)の2乗に比例した値Pf(k−165)、Rb(k)の2乗に比例した値Pb(k))を計算する。そして、最小電力探索手段8において第1及び第2の瞬時電力計算手段6、7から出力される瞬時電力のうち小さい方を探索し、この探索結果をもとに、最小電力選択手段9がRf(k−165)及びRb(k)のうち、瞬時電力Pf(k−165),Pb(k)が小さい相関値を選択し出力する。このような選択は、kの各値ごとに行われる。すべてのkについて順次上記の選択を行なうことで、図8に示すような遅延プロファイルを得る。   In the delay profile estimation apparatus of FIG. 1, the first instantaneous power calculation means 6 and the second instantaneous power calculation means 7 are respectively a correlation sequence Rf (k-165) output from the delay means 5 and a second correlation calculation. The instantaneous power of each correlation value of the correlation sequence Rb (k) output from the means 4 (values Pf (k-165) proportional to the square of Rf (k-165), proportional to the square of Rb (k) Then, the minimum power search means 8 searches for the smaller one of the instantaneous powers output from the first and second instantaneous power calculation means 6 and 7, and the search result is obtained as the search result Pb (k)). Basically, the minimum power selection means 9 selects and outputs a correlation value having a small instantaneous power Pf (k-165), Pb (k) among Rf (k-165) and Rb (k). A selection is made for each value of k. By performing the following above selection, obtaining a delay profile as shown in FIG.

以上のように、実施の形態1の遅延プロファイル推定装置は、繰返し処理を行わずに遅延プロファイルを推定する。その結果、従来手法に比べて、遅延プロファイルを推定するために必要な計算量が少なくなるという効果がある。   As described above, the delay profile estimation apparatus according to Embodiment 1 estimates a delay profile without performing iterative processing. As a result, there is an effect that the amount of calculation required for estimating the delay profile is reduced as compared with the conventional method.

実施の形態2.
図9はこの発明の実施の形態2に係わる遅延プロファイル推定装置を示す。図示の遅延プロファイル推定装置は、第1のPN系列生成手段11と、第2のPN系列生成手段12と、第3のPN系列生成手段13と、第1の相関計算手段14と、第2の相関計算手段15と、第3の相関計算手段16と、プロファイル値計算手段17を備える。
Embodiment 2. FIG.
FIG. 9 shows a delay profile estimation apparatus according to Embodiment 2 of the present invention. The illustrated delay profile estimation apparatus includes a first PN sequence generation unit 11, a second PN sequence generation unit 12, a third PN sequence generation unit 13, a first correlation calculation unit 14, Correlation calculation means 15, third correlation calculation means 16, and profile value calculation means 17 are provided.

図9の遅延プロファイル推定装置に入力される信号は受信信号Rsを所定のサンプリング周波数で標本化した系列であるとする。送信信号とフレームヘッダは、実施の形態1と同様の形式の信号とする。   Assume that the signal input to the delay profile estimation apparatus in FIG. 9 is a sequence obtained by sampling the received signal Rs at a predetermined sampling frequency. The transmission signal and the frame header are signals in the same format as in the first embodiment.

図9における第1、第2及び第3のPN系列生成手段11、12、13で生成されるPN系列の一例T1、T2、T3を図10(a)〜(d)に、フレームヘッダHfとともに示す。
第1のPN系列生成手段11は、図10(d)に示すように、図10(a)のフレームヘッダHfの先頭から数えて第1の所定数C1a番目、例えば15サンプル目から、第2の所定数C1b番目、例えば305サンプル目までの長さ(サンプル数で数えた長さ)S1b、例えば290サンプルのPN系列T1を生成する。PN系列T1の後端からフレームヘッダHfの後端までの長さ(サンプル数で数えた長さ)をS1cで表す。S1cは、
S1c=Lm+Lpre+Lpost−C1a−S1b …(1a)
で与えられる。
Examples of PN sequences T1, T2, and T3 generated by the first, second, and third PN sequence generating units 11, 12, and 13 in FIG. 9 are shown in FIGS. 10 (a) to 10 (d) together with the frame header Hf. Show.
As shown in FIG. 10 (d), the first PN sequence generation means 11 starts from the first predetermined number C1a-th, for example, the 15th sample from the top of the frame header Hf in FIG. A predetermined number C1b, for example, a length up to the 305th sample (length counted by the number of samples) S1b, for example, a PN sequence T1 of 290 samples is generated. The length from the rear end of the PN sequence T1 to the rear end of the frame header Hf (the length counted by the number of samples) is represented by S1c. S1c
S1c = Lm + Lpre + Lpost−C1a−S1b (1a)
Given in.

第2のPN系列生成手段12は、図10(c)に示すように、フレームヘッダHfの先頭から数えて第3の所定数C2a番目、例えば55サンプル目から、第4の所定数C2b番目、例えば、305サンプル目までの長さ(サンプル数で数えた長さ)S2b、例えば250サンプルのPN系列T2を生成する。PN系列T2の後端からフレームヘッダHfの後端までの長さ(サンプル数で数えた長さ)をS2cで表す。S2cは、
S2c=Lm+Lpre+Lpost−C2a−S2b …(1b)
で与えられる。
As shown in FIG. 10C, the second PN sequence generation means 12 counts from the head of the frame header Hf to the third predetermined number C2a, for example, from the 55th sample, the fourth predetermined number C2b. For example, a PN sequence T2 having a length up to the 305th sample (a length counted by the number of samples) S2b, for example, 250 samples, is generated. The length from the rear end of the PN sequence T2 to the rear end of the frame header Hf (the length counted by the number of samples) is represented by S2c. S2c
S2c = Lm + Lpre + Lpost−C2a−S2b (1b)
Given in.

第3のPN系列生成手段13は、図10(b)に示すように、フレームヘッダHfの先頭から数えて第5の所定数C3a番目、例えば80サンプル目から、第6の所定数C3b番目、例えば、280サンプル目までの長さ(サンプル数で数えた長さ)S3b、例えば200サンプルのPN系列T3を生成する。PN系列T3の後端からフレームヘッダHfの後端までの長さ(サンプル数で数えた長さ)をS3cで表す。S3cは、
S3c=Lm+Lpre+Lpost−C3a−S3b …(1c)
で与えられる。
As shown in FIG. 10 (b), the third PN sequence generation means 13 is a fifth predetermined number C3a-th counting from the head of the frame header Hf, for example, from the 80th sample, the sixth predetermined number C3b-th, For example, a length up to the 280th sample (a length counted by the number of samples) S3b, for example, a PN sequence T3 of 200 samples is generated. The length from the rear end of the PN sequence T3 to the rear end of the frame header Hf (the length counted by the number of samples) is represented by S3c. S3c
S3c = Lm + Lpre + Lpost−C3a−S3b (1c)
Given in.

C1a、C2a、C3b、S1b、S2b、S3bは、下記の式(2a)〜(4c)を満たすように定められる。
S1b>(Lm−Lpre−Lpost) …(2a)
S2b>(Lm−Lpre−Lpost) …(2b)
S3b>(Lm−Lpre−Lpost) …(2c)
C1a+S1b>Lm …(3a)
C2a+S2b>Lm …(3b)
C3a+S3b>Lm …(3c)
S1b+S1c>Lm …(4a)
S2b+S2c>Lm …(4b)
S3b+S3c>Lm …(4c)
C1a, C2a, C3b, S1b, S2b, and S3b are determined so as to satisfy the following expressions (2a) to (4c).
S1b> (Lm-Lpre-Lpost) (2a)
S2b> (Lm-Lpre-Lpost) (2b)
S3b> (Lm-Lpre-Lpost) (2c)
C1a + S1b> Lm (3a)
C2a + S2b> Lm (3b)
C3a + S3b> Lm (3c)
S1b + S1c> Lm (4a)
S2b + S2c> Lm (4b)
S3b + S3c> Lm (4c)

受信信号Rsの標本化系列と、第1、第2及び第3のPN系列生成手段11、12、13で生成されたPN系列T1、T2、T3との相関系列R1(k)、R2(k)、R3(k)をそれぞれ図11(a)乃至(c)に示す。ここで、マルチパスが存在しない伝送路を想定している。
図11(a)乃至(c)から、到来波に応じた相関のピークR1a、R2a、R3aが存在する位置k=0を中心に、255サンプル離れた位置k=±255に不要な相関のピークR1b、R1c、R2b、R2c、R3b、R3cが生じることがわかる。
不要な相関のピークR1c、R2c、R3cは、フレームヘッダHfの構成の際、PN系列の先頭から規定長Lpreの系列を該PN系列の後に付加しているために生じるものであり、
不要な相関のピークR1b、R2b、R3bは、フレームヘッダHfの構成の際、PN系列の末尾から規定長Lpostの系列を該PN系列の前に付加しているために生じるものである。
Correlation sequences R1 (k), R2 (k) between the sampling sequence of the received signal Rs and the PN sequences T1, T2, T3 generated by the first, second and third PN sequence generation means 11, 12, 13 ) And R3 (k) are shown in FIGS. 11 (a) to 11 (c), respectively. Here, it is assumed that the transmission path has no multipath.
11 (a) to 11 (c), an unnecessary correlation peak at a position k = ± 255, which is 255 samples away, centering on the position k = 0 where the correlation peaks R1a, R2a, R3a exist according to the incoming wave. It can be seen that R1b, R1c, R2b, R2c, R3b, and R3c are generated.
Unnecessary correlation peaks R1c, R2c, and R3c are caused by adding a sequence of the specified length Lpre from the beginning of the PN sequence after the PN sequence when the frame header Hf is configured.
Unnecessary correlation peaks R1b, R2b, and R3b are generated because a sequence of the specified length Lpost is added before the PN sequence from the end of the PN sequence when the frame header Hf is configured.

相関のピークR1a、R2a、R3a、R1b、R1c、R2b、R2c、R3b、R3cは、それぞれ以下の式(5a)〜(7c)で与えられる。
R1a=S1b …(5a)
R1b=Lpre+Lpost−C1a …(5b)
R1c=Lpre+Lpost−S1c …(5c)
R2a=S2b …(6a)
R2b=Lpre+Lpost−C2a …(6b)
R2c=Lpre+Lpost−S2c …(6c)
R3a=S3b …(7a)
R3b=Lpre+Lpost−C3a …(7b)
R3c=Lpre+Lpost−S3c …(7c)
The correlation peaks R1a, R2a, R3a, R1b, R1c, R2b, R2c, R3b, and R3c are given by the following equations (5a) to (7c), respectively.
R1a = S1b (5a)
R1b = Lpre + Lpost−C1a (5b)
R1c = Lpre + Lpost−S1c (5c)
R2a = S2b (6a)
R2b = Lpre + Lpost−C2a (6b)
R2c = Lpre + Lpost−S2c (6c)
R3a = S3b (7a)
R3b = Lpre + Lpost−C3a (7b)
R3c = Lpre + Lpost−S3c (7c)

図11(a)に示すように、受信信号RsとPN系列T1との相関系列R1(k)における到来波に応じた相関のピークの値R1aは290、不要なピークR1c、R1bの値は50と150である。
また、図11(b)に示すように、受信信号RsとPN系列T2との相関系列R2(k)における到来波に応じたピークの値R2aは250、不要なピークの値R2c、R2bは50と110であり、
図11(c)に示すように、受信信号RsとPN系列T3との相関系列R3(k)における到来波に応じたピークの値R3aは200、不要なピークの値R3c、R3bは25と85である。
As shown in FIG. 11A, the correlation peak value R1a corresponding to the incoming wave in the correlation sequence R1 (k) between the received signal Rs and the PN sequence T1 is 290, and the unnecessary peaks R1c and R1b are 50 values. And 150.
Further, as shown in FIG. 11B, the peak value R2a corresponding to the incoming wave in the correlation sequence R2 (k) between the received signal Rs and the PN sequence T2 is 250, and the unnecessary peak values R2c and R2b are 50. And 110,
As shown in FIG. 11 (c), the peak value R3a corresponding to the incoming wave in the correlation sequence R3 (k) between the received signal Rs and the PN sequence T3 is 200, and unnecessary peak values R3c and R3b are 25 and 85. It is.

マルチパス伝送路の場合、複数の到来波を合わせて受信する。この受信信号RsとPN系列との相関の計算式は次式(8)で表される。   In the case of a multipath transmission path, a plurality of incoming waves are received together. The equation for calculating the correlation between the received signal Rs and the PN sequence is expressed by the following equation (8).

Figure 2010016617
Figure 2010016617

ここで、r(i)は受信信号の標本化系列、r(i)はj番目の到来波の標本化系列、pn(i)はPN系列、MはPN系列の長さ、Lは到来波の数とする。 Here, r (i) is the sampling sequence of the received signal, r j (i) is the sampling sequence of the jth arrival wave, pn (i) is the PN sequence, M is the length of the PN sequence, and L is the arrival The number of waves.

式(8)の最後の行は、複数の到来波を合わせて受信した信号とPN系列との相関値が、各到来波とPN系列との相関結果をすべての到来波に関して足し合わせた値と等しいことを意味する。   The last line of Equation (8) shows that the correlation value between the signal received by combining a plurality of incoming waves and the PN sequence is the sum of the correlation results of each incoming wave and the PN sequence for all the incoming waves. Means equal.

図12(a)に受信信号Rsと第1のPN系列生成手段11で生成されるPN系列T1との相関系列R1(k)を示し、図13(a)に受信信号Rsと第2のPN系列生成手段12で生成されるPN系列T2との相関系列R2(k)を示し、図14(a)における受信信号Rsと第3のPN系列生成手段13で生成されるPN系列T3との相関系列R3(k)を示す。
ここで、3波モデル伝送路を想定し、受信する到来波として、主波と、主波に対して規定のPN系列Tdの長さ(255サンプル)だけ遅れて受信する第1の遅延波と、主波に対して規定のPN系列Tdの長さの2倍(510サンプル)だけ遅れて受信する第2の遅延波が存在する場合を想定する。このような場合を想定するのは、発明が解決しようとする課題の欄で述べたように、主波の到来時刻を基準にして遅延時間差が既知のPN系列の長さと等しい先行波、及び遅延波が存在する場合、従来手法で検出する主波に応じた相関のピーク値に、先行波、及び遅延波による不要な相関のピーク値が加算されるため、遅延プロファイルの推定精度が劣化するという問題があったが、本発明では、そのような問題が解決できることを示すためである。
主波、第1の遅延波、及び第2の遅延波の受信レベルをそれぞれα、β、γとする。
FIG. 12 (a) shows a correlation sequence R1 (k) between the received signal Rs and the PN sequence T1 generated by the first PN sequence generating means 11, and FIG. 13 (a) shows the received signal Rs and the second PN. 14 shows a correlation sequence R2 (k) with the PN sequence T2 generated by the sequence generation means 12, and the correlation between the received signal Rs in FIG. 14 (a) and the PN sequence T3 generated by the third PN sequence generation means 13. The series R3 (k) is shown.
Here, assuming a three-wave model transmission path, as the incoming wave to be received, a main wave, and a first delayed wave that is received with a delay of a specified PN sequence Td length (255 samples) from the main wave, Assume that there is a second delayed wave that is received with a delay of twice the length of the prescribed PN sequence Td (510 samples) with respect to the main wave. Such a case is assumed because, as described in the section of the problem to be solved by the invention, a preceding wave whose delay time difference is equal to the length of a known PN sequence with reference to the arrival time of the main wave, and a delay If there is a wave, an unnecessary correlation peak value due to the preceding wave and the delayed wave is added to the correlation peak value corresponding to the main wave detected by the conventional method, so that the estimation accuracy of the delay profile deteriorates. This is to show that the present invention can solve such a problem.
The reception levels of the main wave, the first delay wave, and the second delay wave are α, β, and γ, respectively.

式(8)より、図12(a)に示した相関系列R1(k)は、図12(b)〜(d)に示すような3つの到来波(主波、第1の遅延波、第2の遅延波)とPN系列T1とのそれぞれの相関系列の和で表される。   From equation (8), the correlation sequence R1 (k) shown in FIG. 12 (a) has three incoming waves (the main wave, the first delay wave, and the first delay wave as shown in FIGS. 12 (b) to (d). 2 delayed waves) and the sum of the correlation sequences of the PN sequence T1.

図12(a)において、例えばk=255(この255は、規定のPN系列Tdの長さ(サンプル数で表す)に等しく、且つ主波に対する第1の遅延波の遅延時間(サンプル数で表す)、及び第1の遅延波に対する第2の遅延波の遅延時間(サンプル数で表す)に等しい)の相関値R1(255)は、第1の遅延波に応じた相関のピーク値290βと、主波、及び第2の遅延波による不要な相関のピーク値150α、50γの和(150α+290β+50γ)で表される。
また、k=0の相関値R1(0)は、主波に応じた相関のピーク値290αと、第1の遅延波による不要な相関のピーク値50βの和(200α+50β)で表され、k=510の相関値R1(510)は、第2の遅延波に応じた相関のピーク値290γと第1の遅延波による不要な相関のピーク値150βとの和(150β+290γ)で表される。
In FIG. 12A, for example, k = 255 (255 is equal to the length of the prescribed PN sequence Td (expressed by the number of samples), and the delay time of the first delayed wave with respect to the main wave (expressed by the number of samples). ) And the delay time (expressed in number of samples) of the second delay wave with respect to the first delay wave) is equal to the correlation peak value 290β corresponding to the first delay wave, It is represented by the sum (150α + 290β + 50γ) of peak values 150α and 50γ of unnecessary correlation due to the main wave and the second delayed wave.
The correlation value R1 (0) at k = 0 is represented by the sum (200α + 50β) of the correlation peak value 290α corresponding to the main wave and the unnecessary correlation peak value 50β due to the first delayed wave, and k = The correlation value R1 (510) of 510 is represented by the sum (150β + 290γ) of the correlation peak value 290γ according to the second delayed wave and the unnecessary correlation peak value 150β due to the first delayed wave.

同様にして、図13(a)に示した相関系列R2(k)は、図13(b)〜(d)に示すような3つの到来波(主波、第1の遅延波、第2の遅延波)とPN系列T2とのそれぞれの相関系列の和で表される。
同様にして、図14(a)に示した相関系列R3(k)は、図14(b)〜(d)に示すような3つの到来波(主波、第1の遅延波、第2の遅延波)とPN系列T3とのそれぞれの相関系列の和で表される。
Similarly, the correlation sequence R2 (k) shown in FIG. 13A has three incoming waves (the main wave, the first delay wave, and the second wave as shown in FIGS. 13B to 13D). (Delayed wave) and the sum of correlation sequences of the PN sequence T2.
Similarly, the correlation sequence R3 (k) shown in FIG. 14A has three incoming waves (main wave, first delay wave, second wave as shown in FIGS. 14B to 14D). (Delayed wave) and the sum of the correlation sequences of the PN sequence T3.

図12(a)、図13(a)および図14(a)に示すように、受信信号RsとPN系列T1との相関系列R1(k)の各相関値は、到来波に応じたピーク値と、他の到来波による不要な相関のピーク値との和である可能性がある。そこで、任意のkにおける受信信号とPN系列T1との相関結果R1(k)を3つの変数x、y、zを用いて次式(9a)のように表すことができる。   As shown in FIGS. 12 (a), 13 (a) and 14 (a), each correlation value of the correlation sequence R1 (k) between the received signal Rs and the PN sequence T1 is a peak value corresponding to the incoming wave. And the sum of unnecessary correlation peak values due to other incoming waves. Therefore, the correlation result R1 (k) between the received signal at an arbitrary k and the PN sequence T1 can be expressed as the following equation (9a) using three variables x, y, and z.

150x+290y+50z=R1(k) …(9a)   150x + 290y + 50z = R1 (k) (9a)

ここで、式(9a)において、左辺の第2項の「290y」は、R1(k)に含まれる到来波に応じた相関のピーク値を示しており、左辺の第1項の「150x」は、上記到来波を基準に255サンプルだけ早く到来した信号による不要な相関のピーク値、左辺の第3項の「50z」は、上記到来波を基準に255サンプルだけ遅れて到来する信号による不要な相関のピーク値を示す。即ち、yに乗算されている値(290)は、図11(a)に示す相関系列に含まれる到来波に応じた相関のピークR1aであり、x、zに乗算されている値(150、50)は、図11(a)に示す不要なピークの値R1b、R1cである。   Here, in Expression (9a), “290y” in the second term on the left side indicates the peak value of the correlation according to the incoming wave included in R1 (k), and “150x” in the first term on the left side. Is an unnecessary correlation peak value due to a signal arriving earlier by 255 samples with respect to the incoming wave, and the third term “50z” on the left side is unnecessary due to a signal arriving with a delay of 255 samples with respect to the incoming wave. The peak value of the correlation is shown. That is, the value (290) multiplied by y is the correlation peak R1a corresponding to the incoming wave included in the correlation sequence shown in FIG. 11 (a), and the values (150, 50) are unnecessary peak values R1b and R1c shown in FIG.

受信信号RsとPN系列T2、T3との相関結果R2(k)、R3(k)も、受信信号RsとPN系列T1との相関結果R1(k)と同様に次式(9b)、(9c)で表される。   Correlation results R2 (k) and R3 (k) between the received signal Rs and the PN sequences T2 and T3 are also expressed by the following equations (9b) and (9c), similarly to the correlation result R1 (k) between the received signal Rs and the PN sequence T1. ).

110x+250y+50z=R2(k) …(9b)
85x+200y+25z=R3(k) …(9c)
110x + 250y + 50z = R2 (k) (9b)
85x + 200y + 25z = R3 (k) (9c)

ここで、x、zは不要な相関ピークの原因となる到来波の受信レベルを示しており、yは所望の到来波の受信レベルを示している。即ち、式(9b)におけるyに乗算されている値(250)は、図11(b)に示す相関系列に含まれる到来波に応じた相関のピークR2aであり、x、zに乗算されている値(110、50)は、図11(b)に示す不要なピークの値R2b、R2cであり、式(9c)におけるyに乗算されている値(200)は、図11(c)に示す相関系列に含まれる到来波に応じた相関のピークR3aであり、式(9c)におけるx、zに乗算されている値(85、25)は、図11(c)に示す不要なピークの値R3b、R3cである。
式(9a)、(9b)、(9c)は、x、y、zを未知の変数とする連立方程式として扱える。
Here, x and z indicate the reception level of an incoming wave that causes an unnecessary correlation peak, and y indicates the reception level of a desired incoming wave. That is, the value (250) multiplied by y in the equation (9b) is a correlation peak R2a corresponding to the incoming wave included in the correlation sequence shown in FIG. 11 (b), and is multiplied by x and z. The values (110, 50) are unnecessary peak values R2b, R2c shown in FIG. 11 (b). The value (200) multiplied by y in the equation (9c) is shown in FIG. 11 (c). The correlation peak R3a corresponding to the incoming wave included in the correlation sequence shown, and the values (85, 25) multiplied by x and z in the equation (9c) are the unnecessary peaks shown in FIG. Values R3b and R3c.
Expressions (9a), (9b), and (9c) can be treated as simultaneous equations with x, y, and z as unknown variables.

図11(a)〜(c)に示すピークの値を同図に示す記号で置き換えて、式(9a)〜(9c)を一般化すれば下記の式(10a)〜(10c)のようになる。
R1b・x+R1a・y+R1c・z=R1(k) …(10a)
R2b・x+R2a・y+R2c・z=R2(k) …(10b)
R3b・x+R3a・y+R3c・z=R3(k) …(10c)
When the peak values shown in FIGS. 11A to 11C are replaced with the symbols shown in FIG. 11 and the equations (9a) to (9c) are generalized, the following equations (10a) to (10c) are obtained. Become.
R1b.x + R1a.y + R1c.z = R1 (k) (10a)
R2b.x + R2a.y + R2c.z = R2 (k) (10b)
R3b.x + R3a.y + R3c.z = R3 (k) (10c)

次式に式(9a)〜(9c)を連立で解いたときの解を示す。
x={15×R1(k)−11×R2(k)−8×R3(k)}/360 …(11a)
y={−3×R1(k)+R2(k)+4×R3(k)}/180 …(11b)
z={−15×R1(k)+107×R2(k)−112×R3(k)}/1800
…(11c)
The following equation shows solutions when equations (9a) to (9c) are solved simultaneously.
x = {15 * R1 (k) -11 * R2 (k) -8 * R3 (k)} / 360 (11a)
y = {− 3 × R1 (k) + R2 (k) + 4 × R3 (k)} / 180 (11b)
z = {− 15 × R1 (k) + 107 × R2 (k) −112 × R3 (k)} / 1800
... (11c)

図9におけるプロファイル値計算手段17は、第1、第2及び第3の相関計算手段14、15、16から出力される相関値R1(k)、R2(k)、R3(k)を式(11a)、(11b)、(11c)に代入し、サンプル毎に連立方程式(9a)、(9b)、(9c)の解を求める。そして、3つの解x、y、zのうち、yを遅延プロファイルの推定結果として出力する。   9 calculates the correlation values R1 (k), R2 (k), and R3 (k) output from the first, second, and third correlation calculation units 14, 15, and 16 by the formula ( 11a), (11b), and (11c) are substituted, and solutions of simultaneous equations (9a), (9b), and (9c) are obtained for each sample. Of the three solutions x, y, and z, y is output as a delay profile estimation result.

例えば、図12(b)〜(d)に示すような3つの到来波(主波、第1の遅延波、第2の遅延波)が存在する場合、図12(a)、図13(a)、図14(a)に示される、第1、第2及び第3の相関計算手段14、15、16から出力される相関系列R1(k)、R2(k)、R3(k)におけるk=0の相関結果R1(0)(=290α+50β)、R2(0)(=250α+50β)、R3(0)(=200α+25β)を式(11a)、(11b)、(11c)に代入すると、x=0、y=α、z=βとなり、yが主波の受信レベルとなる。
また、k=255における相関結果R1(255)、R2(255)、R3(255)を式(11a)、(11b)、(11c)に代入すると、x=α、y=β、z=γとなり、yが第1の遅延波の受信レベルとなる。
k=510における相関結果R1(510)、R2(510)、R3(510)を式(11a)、(11b)、(11c)に代入すると、x=β、y=γ、z=0となり、yが第2の遅延波の受信レベルとなる。
また、相関のピークが存在しない範囲では、相関結果が非常に小さな値となるので、x、y、zはすべて0近傍の小さな値となる。
For example, when there are three incoming waves (main wave, first delay wave, and second delay wave) as shown in FIGS. 12B to 12D, FIGS. 12A and 13A are used. ), K in the correlation sequences R1 (k), R2 (k), and R3 (k) output from the first, second, and third correlation calculating means 14, 15, and 16 shown in FIG. = 0 correlation results R1 (0) (= 290α + 50β), R2 (0) (= 250α + 50β), R3 (0) (= 200α + 25β) are substituted into equations (11a), (11b), and (11c), x = 0, y = α, z = β, and y is the reception level of the main wave.
Further, when the correlation results R1 (255), R2 (255), and R3 (255) at k = 255 are substituted into the equations (11a), (11b), and (11c), x = α, y = β, z = γ And y is the reception level of the first delayed wave.
Substituting correlation results R1 (510), R2 (510), and R3 (510) at k = 510 into equations (11a), (11b), and (11c) yields x = β, y = γ, and z = 0. y is the reception level of the second delayed wave.
In the range where no correlation peak exists, the correlation result has a very small value, so x, y, and z all have small values near zero.

以上、相関値のピークR1a、R1b、R1c、R2a、R2b、R2c、R3a、R3b、R3bが図11(a)〜(c)に示す値である場合について説明したが、一般的には、式(10a)〜(10c)をyについて連立で解いて、ピーク値R1a〜R3c、並びに相関系列R1(k)、R2(k)、R3(k)を変数とする、yを表す式(式(11b)と同様の式)を作成しておき、そのようにして作成された式に、任意のkに対する相関系列R1(k)、R2(k)、R3(k)の値を代入することで、当該kの遅延プロファイルの値が求められる。すべてのkについて求めた値の系列が遅延プロファイルy(k)となる。   The correlation value peaks R1a, R1b, R1c, R2a, R2b, R2c, R3a, R3b, and R3b have been described with reference to the values shown in FIGS. (10a) to (10c) are solved simultaneously for y, and the peak value R1a to R3c and the correlation series R1 (k), R2 (k), and R3 (k) are used as variables. 11b), and the values of the correlation sequences R1 (k), R2 (k), and R3 (k) for an arbitrary k are substituted into the formula thus created. , The value of k delay profile is obtained. A series of values obtained for all k is a delay profile y (k).

以上のように、実施の形態2の遅延プロファイル推定装置は、主波に対して遅延時間差が既定のPN系列長と等しい遅延波が存在する場合でも、不要な相関のピークを除去し、遅延プロファイルを正確に推定することができる。言い換えると、不要な相関のピークを除去し、到来波に応じた相関のピークのみを出力することができる。また、繰返し処理を行わずに遅延プロファイルを推定するため、従来手法に比べて計算量が少ないという効果がある。   As described above, the delay profile estimation apparatus according to the second embodiment removes unnecessary correlation peaks even when a delay wave having a delay time difference equal to a predetermined PN sequence length exists with respect to the main wave, and delay profile Can be estimated accurately. In other words, unnecessary correlation peaks can be removed, and only correlation peaks according to the incoming wave can be output. In addition, since the delay profile is estimated without performing iterative processing, there is an effect that the amount of calculation is smaller than that of the conventional method.

さらに、本発明の実施の形態1及び実施の形態2のいずれにおいても、計算量を少なくすることができるので、装置におけるエネルギー消費を削減することができる。   Furthermore, in both Embodiment 1 and Embodiment 2 of the present invention, the amount of calculation can be reduced, so that energy consumption in the apparatus can be reduced.

この発明の実施の形態1の遅延プロファイル推定装置を示すブロック図である。It is a block diagram which shows the delay profile estimation apparatus of Embodiment 1 of this invention. (a)及び(b)は、この発明の実施の形態1及び実施の形態2における送信信号の構成を示す図及びフレームヘッダHfを示す図である。(A) And (b) is the figure which shows the structure of the transmission signal in Embodiment 1 and Embodiment 2 of this invention, and the figure which shows the frame header Hf. (a)及び(b)は、この発明の実施の形態1で使用するPN系列TfをフレームヘッダHfとともに示す図である。(A) And (b) is a figure which shows the PN series Tf used in Embodiment 1 of this invention with the frame header Hf. (a)及び(b)は、この発明の実施の形態1で使用するPN系列TbをフレームヘッダHfとともに示す図である。(A) And (b) is a figure which shows the PN series Tb used in Embodiment 1 of this invention with the frame header Hf. (a)及び(b)は、この発明の実施の形態1における受信信号RsとPN系列Tfとの相関系列Rf(k)を示す図である。(A) And (b) is a figure which shows correlation series Rf (k) of the received signal Rs and PN series Tf in Embodiment 1 of this invention. (a)及び(b)は、この発明の実施の形態1における受信信号RsとPN系列Tbとの相関系列Rb(k)を示す図である。(A) And (b) is a figure which shows correlation series Rb (k) of the received signal Rs and PN series Tb in Embodiment 1 of this invention. (a)及び(b)は、この発明の実施の形態1で遅延プロファイルを推定するための相関系列Rf(k−165)及びRb(k)を示す図である。(A) And (b) is a figure which shows correlation series Rf (k-165) and Rb (k) for estimating a delay profile in Embodiment 1 of this invention. この発明の実施の形態1における遅延プロファイルの推定結果を示す図である。It is a figure which shows the estimation result of the delay profile in Embodiment 1 of this invention. この発明の実施の形態2の遅延プロファイル推定装置を示すブロック図である。It is a block diagram which shows the delay profile estimation apparatus of Embodiment 2 of this invention. (a)〜(d)は、この発明の実施の形態2で使用する3つのPN系列T1、T2、T3をフレームヘッダHfとともに示す図である。(A)-(d) is a figure which shows the three PN series T1, T2, T3 used with Embodiment 2 of this invention with the frame header Hf. (a)〜(c)は、この発明の実施の形態2で使用するPN系列と受信信号Rsとの相関結果を示す図である。(A)-(c) is a figure which shows the correlation result of the PN series and reception signal Rs which are used in Embodiment 2 of this invention. (a)〜(d)は、この発明の実施の形態2における3波モデル伝送路での受信信号Rsと第1のPN系列生成手段11で生成されたPN系列T1との相関結果を示す示す図である。(A)-(d) shows the correlation result of the received signal Rs on the three-wave model transmission line and the PN sequence T1 generated by the first PN sequence generating means 11 in the second embodiment of the present invention. FIG. (a)〜(d)は、この発明の実施の形態2における3波モデル伝送路での受信信号Rsと第2のPN系列生成手段12で生成されたPN系列T2との相関結果を示す示す図である。(A)-(d) shows the correlation result of the received signal Rs on the three-wave model transmission line and the PN sequence T2 generated by the second PN sequence generating means 12 in the second embodiment of the present invention. FIG. (a)〜(d)は、この発明の実施の形態2における3波モデル伝送路での受信信号Rsと第3のPN系列生成手段13で生成されたPN系列T3との相関結果を示す示す図である。(A)-(d) shows the correlation result of the received signal Rs on the three-wave model transmission line and the PN sequence T3 generated by the third PN sequence generating means 13 in the second embodiment of the present invention. FIG.

符号の説明Explanation of symbols

1 第1のPN系列生成手段、 2 第1の相関計算手段、 3 第2のPN系列生成手段、 4 第2の相関計算手段、 5 遅延手段、 6 第1の瞬時電力計算手段、 7 第2の瞬時電力計算手段、 8 最小電力探索手段、 9 最小電力選択手段、 10 プロファイル値計算手段、 11 第1のPN系列生成手段、 12 第2のPN系列生成手段、 13 第3のPN系列生成手段、 14 第1の相関計算手段、 15 第2の相関計算手段、 16 第3の相関計算手段、 17 プロファイル値計算手段。
DESCRIPTION OF SYMBOLS 1 1st PN sequence production | generation means, 2 1st correlation calculation means, 3 2nd PN sequence production | generation means, 4 2nd correlation calculation means, 5 Delay means, 6 1st instantaneous electric power calculation means, 7 2nd Instantaneous power calculation means, 8 minimum power search means, 9 minimum power selection means, 10 profile value calculation means, 11 first PN sequence generation means, 12 second PN sequence generation means, 13 third PN sequence generation means , 14 first correlation calculation means, 15 second correlation calculation means, 16 third correlation calculation means, 17 profile value calculation means.

Claims (8)

送信機で、既知の擬似ランダム(PN)系列における先頭の規定長の系列を該PN系列の後に付加し、該PN系列における末尾の規定長の系列を該PN系列の前に付加した構造を持つフレームヘッダと、送信すべき情報を含む有効シンボルとで構成される伝送シンボルを伝送単位として送信された信号を、伝送路を介して受信し、受信信号を所定のサンプリング周波数で標本化して得られる標本化系列から伝送路の遅延プロファイルを推定する遅延プロファイル推定装置において、
それぞれフレームヘッダの一部を構成する2つの異なるPN系列を生成する第1及び第2のPN系列生成手段と、
前記第1のPN系列生成手段によって生成されたPN系列と受信信号の標本化系列との相関を計算する第1の相関計算手段と、
前記第2のPN系列生成手段によって生成されたPN系列と受信信号の標本化系列との相関を計算する第2の相関計算手段と、
前記第1及び第2の相関計算手段で得られる2つの相関系列から遅延プロファイルを推定するプロファイル値計算手段を備え、
前記プロファイル値計算手段が、前記第1及び第2の相関計算手段から出力される2つの相関系列の互いに対応する相関値の瞬時電力を比較して瞬時電力の小さい相関値を選択し、選択した結果を遅延プロファイルとして出力する
ことを特徴とする遅延プロファイル推定装置。
The transmitter has a structure in which a first specified length sequence in a known pseudo-random (PN) sequence is added after the PN sequence, and a last specified length sequence in the PN sequence is added before the PN sequence. It is obtained by receiving a signal transmitted using a transmission symbol composed of a frame header and effective symbols including information to be transmitted as a transmission unit, and sampling the received signal at a predetermined sampling frequency. In a delay profile estimation device that estimates a delay profile of a transmission path from a sampling sequence,
First and second PN sequence generating means for generating two different PN sequences, each forming part of the frame header;
First correlation calculating means for calculating a correlation between the PN sequence generated by the first PN sequence generating means and the sampling sequence of the received signal;
Second correlation calculating means for calculating the correlation between the PN sequence generated by the second PN sequence generating means and the sampling sequence of the received signal;
Profile value calculation means for estimating a delay profile from two correlation sequences obtained by the first and second correlation calculation means;
The profile value calculation means selects a correlation value having a small instantaneous power by comparing the instantaneous powers of the correlation values corresponding to each other of the two correlation sequences output from the first and second correlation calculation means A delay profile estimation apparatus, characterized in that the result is output as a delay profile.
前記第1及び第2のPN系列生成手段が生成するPN系列がともに前記規定長の長さを有することを特徴とする請求項1に記載の遅延プロファイル推定装置。   2. The delay profile estimation apparatus according to claim 1, wherein both of the PN sequences generated by the first and second PN sequence generation units have the specified length. 前記第1のPN系列生成手段が生成するPN系列が、フレームヘッダの先頭から前記規定長分の長さのPN系列と同じものであり、
前記第2のPN系列生成手段が生成するPN系列が、フレームヘッダの後端から前記規定長分の長さのPN系列と同じものである
ことを特徴とする請求項1に記載の遅延プロファイル推定装置。
The PN sequence generated by the first PN sequence generation means is the same as the PN sequence having a length corresponding to the specified length from the beginning of a frame header,
2. The delay profile estimation according to claim 1, wherein the PN sequence generated by the second PN sequence generation unit is the same as the PN sequence having a length corresponding to the specified length from a rear end of a frame header. apparatus.
送信機で、既知の擬似ランダム(PN)系列における先頭の規定長の系列を該PN系列の後に付加し、該PN系列における末尾の規定長の系列を該PN系列の前に付加した構造を持つフレームヘッダと、送信すべき情報を含む有効シンボルとで構成される伝送シンボルを伝送単位として送信された信号を、伝送路を介して受信し、受信信号を所定のサンプリング周波数で標本化して得られる標本化系列から伝送路の遅延プロファイルを推定する遅延プロファイル推定装置において、
それぞれフレームヘッダの一部を構成する3つの異なるPN系列を生成する第1、第2及び第3のPN系列生成手段と、
前記第1のPN系列生成手段で生成されるPN系列と受信信号の標本化系列との相関を計算する第1の相関計算手段と、
前記第2のPN系列生成手段で生成されるPN系列と受信信号の標本化系列との相関を計算する第2の相関計算手段と、
前記第3のPN系列生成手段で生成されるPN系列と受信信号の標本化系列との相関を計算する第3の相関計算手段と、
前記第1、第2及び第3の相関計算手段で得られる3つの相関系列から遅延プロファイルを推定するプロファイル値計算手段を備え、
前記プロファイル値計算手段が、前記第1、第2及び第3の相関計算手段から出力される3つの相関系列の互いに対応する相関値をもとに遅延プロファイルを推定する
ことを特徴とする遅延プロファイル推定装置。
The transmitter has a structure in which a first specified length sequence in a known pseudo-random (PN) sequence is added after the PN sequence, and a last specified length sequence in the PN sequence is added before the PN sequence. It is obtained by receiving a signal transmitted using a transmission symbol composed of a frame header and effective symbols including information to be transmitted as a transmission unit, and sampling the received signal at a predetermined sampling frequency. In a delay profile estimation device that estimates a delay profile of a transmission path from a sampling sequence,
First, second and third PN sequence generating means for generating three different PN sequences, each forming part of a frame header;
First correlation calculating means for calculating a correlation between the PN sequence generated by the first PN sequence generating means and the sampling sequence of the received signal;
Second correlation calculation means for calculating the correlation between the PN sequence generated by the second PN sequence generation means and the sampling sequence of the received signal;
Third correlation calculating means for calculating the correlation between the PN sequence generated by the third PN sequence generating means and the sampling sequence of the received signal;
Profile value calculation means for estimating a delay profile from three correlation sequences obtained by the first, second and third correlation calculation means;
The profile value calculating means estimates a delay profile based on correlation values corresponding to each other of the three correlation sequences output from the first, second and third correlation calculating means. Estimating device.
前記プロファイル値計算手段は、
前記第1の相関計算手段から出力される相関系列の到来波に応じた相関のピーク及び第1及び第2の不要な相関のピークと、
前記第2の相関計算手段から出力される相関系列の到来波に応じた相関のピーク及び第1及び第2の不要な相関のピークと、
前記第3の相関計算手段から出力される相関系列の到来波に応じた相関のピーク及び第1及び第2の不要な相関のピークと
に基づいて前記遅延プロファイルを推定するものであり、
前記第1の不要な相関のピークは、前記フレームヘッダの構成に際し、前記PN系列における末尾の規定長の系列を該PN系列の前に付加しているために生じるものであり、
前記第2の不要な相関のピークは、前記フレームヘッダの構成に際し、前記PN系列における先頭の規定長の系列を該PN系列の後に付加しているために生じるものである
ことを特徴とする請求項4に記載の遅延プロファイル推定装置。
The profile value calculation means includes
A correlation peak corresponding to an incoming wave of a correlation sequence output from the first correlation calculation means and first and second unnecessary correlation peaks;
A correlation peak according to an incoming wave of a correlation sequence output from the second correlation calculation means and first and second unnecessary correlation peaks;
Estimating the delay profile based on a correlation peak corresponding to an incoming wave of a correlation sequence output from the third correlation calculation means and first and second unnecessary correlation peaks;
The first unnecessary correlation peak is generated when a frame having a specified length at the end of the PN sequence is added before the PN sequence when the frame header is configured.
The second unnecessary correlation peak is caused by adding a sequence having a specified length at the beginning of the PN sequence after the PN sequence in the configuration of the frame header. Item 5. The delay profile estimation device according to Item 4.
前記プロファイル値計算手段は、
前記第1の相関計算手段から出力される相関系列R1(k)の到来波に応じた相関のピークR1a及び第1及び第2の不要な相関のピークR1b、R1cと、
前記第2の相関計算手段から出力される相関系列R2(k)の到来波に応じた相関のピークR2a及び第1及び第2の不要な相関のピークR2b、R2cと、
前記第3の相関計算手段から出力される相関系列R3(k)の到来波に応じた相関のピークR3a及び第1及び第2の不要な相関のピークR3b、R3cとを関係付ける下記の式、
R1b・x+R1a・y+R1c・z=R1(k)
R2b・x+R2a・y+R2c・z=R2(k)
R3b・x+R3a・y+R3c・z=R3(k)
を連立で解いたときのyの値を、前記遅延プロファイルの値として算出する
ことを特徴とする請求項4に記載の遅延プロファイル推定装置。
The profile value calculation means includes
A correlation peak R1a and first and second unnecessary correlation peaks R1b, R1c according to an incoming wave of the correlation sequence R1 (k) output from the first correlation calculation means;
A correlation peak R2a corresponding to an incoming wave of the correlation sequence R2 (k) output from the second correlation calculating means, and first and second unnecessary correlation peaks R2b, R2c;
The following equation relating the correlation peak R3a according to the incoming wave of the correlation sequence R3 (k) output from the third correlation calculation means and the first and second unnecessary correlation peaks R3b, R3c:
R1b.x + R1a.y + R1c.z = R1 (k)
R2b.x + R2a.y + R2c.z = R2 (k)
R3b.x + R3a.y + R3c.z = R3 (k)
The delay profile estimation apparatus according to claim 4, wherein a value of y when the two are solved simultaneously is calculated as the value of the delay profile.
送信機で、既知の擬似ランダム(PN)系列における先頭の規定長の系列を該PN系列の後に付加し、該PN系列における末尾の規定長の系列を該PN系列の前に付加した構造を持つフレームヘッダと、送信すべき情報を含む有効シンボルとで構成される伝送シンボルを伝送単位として送信された信号を、伝送路を介して受信し、受信信号を所定のサンプリング周波数で標本化して得られる標本化系列から伝送路の遅延プロファイルを推定する遅延プロファイル推定方法において、
それぞれフレームヘッダの一部を構成する2つの異なるPN系列を生成する第1及び第2のPN系列生成ステップと、
前記第1のPN系列生成ステップによって生成されたPN系列と受信信号の標本化系列との相関を計算する第1の相関計算ステップと、
前記第2のPN系列生成ステップによって生成されたPN系列と受信信号の標本化系列との相関を計算する第2の相関計算ステップと、
前記第1及び第2の相関計算ステップで得られる2つの相関系列から遅延プロファイルを推定するプロファイル値計算ステップを備え、
前記プロファイル値計算ステップが、前記第1及び第2の相関計算ステップから出力される2つの相関系列の互いに対応する相関値の瞬時電力を比較して瞬時電力の小さい相関値を選択し、選択した結果を遅延プロファイルとして出力する
ことを特徴とする遅延プロファイル推定方法。
The transmitter has a structure in which a first specified length sequence in a known pseudo-random (PN) sequence is added after the PN sequence, and a last specified length sequence in the PN sequence is added before the PN sequence. It is obtained by receiving a signal transmitted using a transmission symbol composed of a frame header and effective symbols including information to be transmitted as a transmission unit, and sampling the received signal at a predetermined sampling frequency. In a delay profile estimation method for estimating a delay profile of a transmission path from a sampling sequence,
First and second PN sequence generation steps for generating two different PN sequences each forming part of the frame header;
A first correlation calculating step of calculating a correlation between the PN sequence generated by the first PN sequence generating step and the sampling sequence of the received signal;
A second correlation calculating step of calculating a correlation between the PN sequence generated by the second PN sequence generating step and the sampling sequence of the received signal;
A profile value calculation step for estimating a delay profile from two correlation sequences obtained in the first and second correlation calculation steps;
The profile value calculating step selects a correlation value having a small instantaneous power by comparing the instantaneous powers of the correlation values corresponding to each other of the two correlation sequences output from the first and second correlation calculating steps. A delay profile estimation method characterized in that the result is output as a delay profile.
送信機で、既知の擬似ランダム(PN)系列における先頭の規定長の系列を該PN系列の後に付加し、該PN系列における末尾の規定長の系列を該PN系列の前に付加した構造を持つフレームヘッダと、送信すべき情報を含む有効シンボルとで構成される伝送シンボルを伝送単位として送信された信号を、伝送路を介して受信し、受信信号を所定のサンプリング周波数で標本化して得られる標本化系列から伝送路の遅延プロファイルを推定する遅延プロファイル推定方法において、
それぞれフレームヘッダの一部を構成する3つの異なるPN系列を生成する第1、第2及び第3のPN系列生成ステップと、
前記第1のPN系列生成ステップで生成されるPN系列と受信信号の標本化系列との相関を計算する第1の相関計算ステップと、
前記第2のPN系列生成ステップで生成されるPN系列と受信信号の標本化系列との相関を計算する第2の相関計算ステップと、
前記第3のPN系列生成ステップで生成されるPN系列と受信信号の標本化系列との相関を計算する第3の相関計算ステップと、
前記第1、第2及び第3の相関計算ステップで得られる3つの相関系列から遅延プロファイルを推定するプロファイル値計算ステップを備え、
前記プロファイル値計算ステップが、前記第1、第2及び第3の相関計算ステップから出力される3つの相関系列の互いに対応する相関値をもとに遅延プロファイルを推定する
ことを特徴とする遅延プロファイル推定方法。
The transmitter has a structure in which a first specified length sequence in a known pseudo-random (PN) sequence is added after the PN sequence, and a last specified length sequence in the PN sequence is added before the PN sequence. It is obtained by receiving a signal transmitted using a transmission symbol composed of a frame header and effective symbols including information to be transmitted as a transmission unit, and sampling the received signal at a predetermined sampling frequency. In a delay profile estimation method for estimating a delay profile of a transmission path from a sampling sequence,
First, second and third PN sequence generation steps for generating three different PN sequences each forming part of the frame header;
A first correlation calculating step of calculating a correlation between the PN sequence generated in the first PN sequence generating step and the sampling sequence of the received signal;
A second correlation calculating step of calculating a correlation between the PN sequence generated in the second PN sequence generating step and the sampling sequence of the received signal;
A third correlation calculating step of calculating a correlation between the PN sequence generated in the third PN sequence generating step and the sampling sequence of the received signal;
A profile value calculating step for estimating a delay profile from the three correlation sequences obtained in the first, second and third correlation calculating steps;
The delay profile is characterized in that the profile value calculation step estimates a delay profile based on correlation values corresponding to each other of the three correlation sequences output from the first, second and third correlation calculation steps. Estimation method.
JP2008174584A 2008-07-03 2008-07-03 Delay profile estimation apparatus and delay profile estimation method Active JP5100536B2 (en)

Priority Applications (3)

Application Number Priority Date Filing Date Title
JP2008174584A JP5100536B2 (en) 2008-07-03 2008-07-03 Delay profile estimation apparatus and delay profile estimation method
CN201210447494.9A CN102970097B (en) 2008-07-03 2009-04-10 Delay distribution estimation device and delay distribution estimation method
CN2009101299556A CN101621310B (en) 2008-07-03 2009-04-10 Delay distribution estimation device and delay distribution estimation method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2008174584A JP5100536B2 (en) 2008-07-03 2008-07-03 Delay profile estimation apparatus and delay profile estimation method

Publications (2)

Publication Number Publication Date
JP2010016617A true JP2010016617A (en) 2010-01-21
JP5100536B2 JP5100536B2 (en) 2012-12-19

Family

ID=41514392

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2008174584A Active JP5100536B2 (en) 2008-07-03 2008-07-03 Delay profile estimation apparatus and delay profile estimation method

Country Status (2)

Country Link
JP (1) JP5100536B2 (en)
CN (2) CN101621310B (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2010226161A (en) * 2009-03-19 2010-10-07 Mitsubishi Electric Corp Delay profile estimating device and image display

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH0983485A (en) * 1995-09-13 1997-03-28 Canon Inc Code division multiplex communication method and device
JPH10178368A (en) * 1996-12-19 1998-06-30 Fujitsu Ltd Rake mixing circuit
JP2009065673A (en) * 2007-09-05 2009-03-26 Samsung Electronics Co Ltd Cfo synchronization method and apparatus of ofdm receiver
JP2009201083A (en) * 2008-01-22 2009-09-03 Nippon Telegr & Teleph Corp <Ntt> Radio transmitting apparatus, radio receiving apparatus, radio communication system, radio communication method, radio transmission method, and radio reception method
JP2009284036A (en) * 2008-05-20 2009-12-03 Mitsubishi Electric Corp Transmission path estimating device and method

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
SE516143C2 (en) * 1999-05-10 2001-11-26 Ericsson Telefon Ab L M Estimate of pure delay
US7142586B2 (en) * 2002-09-18 2006-11-28 Telefonaktiebolaget Lm Ericsson (Publ) Robust delay estimation architecture
WO2004052039A1 (en) * 2002-12-02 2004-06-17 Nokia Corporation Estimation of a signal delay
CN1719819A (en) * 2005-07-08 2006-01-11 广州海格通信有限公司 A kind of improved OFDM time synchronization method based on the PN sequence
CN100561997C (en) * 2005-11-29 2009-11-18 华为技术有限公司 Sequence generation system and method, transmission and method for synchronous and power Adjustment System

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH0983485A (en) * 1995-09-13 1997-03-28 Canon Inc Code division multiplex communication method and device
JPH10178368A (en) * 1996-12-19 1998-06-30 Fujitsu Ltd Rake mixing circuit
JP2009065673A (en) * 2007-09-05 2009-03-26 Samsung Electronics Co Ltd Cfo synchronization method and apparatus of ofdm receiver
JP2009201083A (en) * 2008-01-22 2009-09-03 Nippon Telegr & Teleph Corp <Ntt> Radio transmitting apparatus, radio receiving apparatus, radio communication system, radio communication method, radio transmission method, and radio reception method
JP2009284036A (en) * 2008-05-20 2009-12-03 Mitsubishi Electric Corp Transmission path estimating device and method

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2010226161A (en) * 2009-03-19 2010-10-07 Mitsubishi Electric Corp Delay profile estimating device and image display

Also Published As

Publication number Publication date
JP5100536B2 (en) 2012-12-19
CN102970097A (en) 2013-03-13
CN102970097B (en) 2015-05-27
CN101621310B (en) 2013-06-12
CN101621310A (en) 2010-01-06

Similar Documents

Publication Publication Date Title
JP4783481B1 (en) Ultrasonic measurement method and ultrasonic measurement apparatus
US7336647B2 (en) System and method for ranging
KR101092557B1 (en) Apparatus for detecting synchronization and VSB receiver using the same and method thereof
RU2733419C1 (en) Transmitter and receiver and corresponding methods
JP2009545754A5 (en)
US7616723B2 (en) Method for symbol timing synchronization and apparatus thereof
JP2005505966A5 (en)
KR20130104289A (en) Apparatus and method for estimating offset value, receiving apparatus, and method for processing signal in the same
EP3110095A1 (en) Ambiguity resolution in positioning measurements
US9516616B2 (en) Method and apparatus for estimating frequency errors
CN102298097A (en) Method for estimating thunder impulse signal Time Difference of Arrival (TDOA)
JP5042118B2 (en) Transmission path estimation apparatus and method
US7728763B2 (en) High resolution ranging apparatus and method using UWB
CN101611606B (en) Method for symbol synchronization of received digital signal and digital signal receiver using the same method
KR101421156B1 (en) Apparatus and Method for acquiring GPS Signal using Multi-stage Partial Cross-Correlator
CN104035109A (en) Weak signal capturing method based on 1/5 bit differential circulation coherent integration
JP5100536B2 (en) Delay profile estimation apparatus and delay profile estimation method
JP5918351B2 (en) Signal search method, signal search program, signal search device, GNSS signal reception device, and information equipment terminal
CN105656511B (en) Differential correlation acquisition method suitable for environment with frequency offset and low signal-to-noise ratio
JP3998023B2 (en) Ranging and position measuring method using spread spectrum signal and apparatus for performing the method
CN105681240B (en) High-precision frequency offset estimation method suitable for low signal-to-noise ratio environment
JP3795885B2 (en) Reception device and reception control method
US8238408B2 (en) Method for processing a sampled rectified ultra wide band signal
KR101644560B1 (en) 2-STEP FDOA/FDOA estimation Method and Apparatus
US20120242528A1 (en) Receiving apparatus and receiving method

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20110210

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20120615

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20120703

A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20120803

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

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

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

Free format text: PAYMENT UNTIL: 20151005

Year of fee payment: 3

R150 Certificate of patent or registration of utility model

Ref document number: 5100536

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

Free format text: JAPANESE INTERMEDIATE CODE: R150

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