JP2012222487A - Transmission path characteristic estimation device, transmission path characteristic estimation method, and program - Google Patents

Transmission path characteristic estimation device, transmission path characteristic estimation method, and program Download PDF

Info

Publication number
JP2012222487A
JP2012222487A JP2011084332A JP2011084332A JP2012222487A JP 2012222487 A JP2012222487 A JP 2012222487A JP 2011084332 A JP2011084332 A JP 2011084332A JP 2011084332 A JP2011084332 A JP 2011084332A JP 2012222487 A JP2012222487 A JP 2012222487A
Authority
JP
Japan
Prior art keywords
transmission path
unit
moving average
smoothing
average value
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
JP2011084332A
Other languages
Japanese (ja)
Other versions
JP5610538B2 (en
Inventor
Yoshio Takarada
佳央 寶田
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.)
NEC Network and Sensor Systems Ltd
Original Assignee
NEC Network and Sensor Systems Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by NEC Network and Sensor Systems Ltd filed Critical NEC Network and Sensor Systems Ltd
Priority to JP2011084332A priority Critical patent/JP5610538B2/en
Publication of JP2012222487A publication Critical patent/JP2012222487A/en
Application granted granted Critical
Publication of JP5610538B2 publication Critical patent/JP5610538B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Abstract

PROBLEM TO BE SOLVED: To reduce the influence of noise, and also to improve the time followability of a transmission path characteristic.SOLUTION: Infinite impulse response (IIR) filters 109, 110 calculate the movement average values of transmission path characteristics estimated by a current transmission path characteristic estimation part. A next transmission path characteristic estimation part estimates a next-time transmission path characteristic based on the outputs of the IIR filters 109, 110 and a coefficient stored by a coefficient storage part 114. The coefficient stored by the coefficient storage part 114 is the coefficient which is calculated based on a difference between the outputs of the movement average values of transmission path characteristics previously calculated by the IIR filters 109, 110 with different oblivion coefficients and also a transmission path characteristic estimated by a pilot complex division part 101.

Description

本発明は、連続する複数の受信シンボルから次の受信シンボルの伝送路特性を推定する伝送路特性推定装置、伝送路特性推定方法、及びプログラムに関する。   The present invention relates to a transmission line characteristic estimation apparatus, a transmission line characteristic estimation method, and a program for estimating a transmission line characteristic of a next reception symbol from a plurality of consecutive reception symbols.

連続通信を行うOFDM(直交周波数分割多重:Orthogonal Frequency Division Multiplexing)通信においては、送信装置においてパイロットシンボルを所定の間隔で挿入し、受信装置において当該パイロットシンボルに基づいて伝送路特性を推定して、復調を行う方式が用いられている。なお、パイロットシンボルとは、送信装置と受信装置において既知となっている信号パターンのことである。
伝送路推定の簡易な方法としては、パイロットシンボル間を時間方向で線形に内挿補間することで求める方法が用いられている(例えば、特許文献1を参照)。
In OFDM (Orthogonal Frequency Division Multiplexing) communication in which continuous communication is performed, a transmitter apparatus inserts pilot symbols at predetermined intervals, and a receiver apparatus estimates channel characteristics based on the pilot symbols. A method of performing demodulation is used. A pilot symbol is a signal pattern that is known in the transmission apparatus and the reception apparatus.
As a simple method for channel estimation, a method is used in which pilot symbols are linearly interpolated in the time direction (see, for example, Patent Document 1).

特開2002−374418号公報JP 2002-374418 A

しかしながら、パイロットシンボル間を時間方向で線形に内挿補間する方法は、ノイズの影響が大きく、図7に示すようにまた伝送路特性がI−Q(同相−直交:In phase−Quadrature phase)平面において曲線を描くように変化する場合、推定した伝送路特性と実際の伝送路特性との間でズレが生じてしまうという問題がある。なお図7は、内挿補間による伝送路特性の推定結果を示す図である。
また、この方法を用いる場合、パイロットシンボルとその次のパイロットシンボルを受信するまでその間にあるデータシンボルを復調できないため、受信データを記憶しておくメモリが必要となり、さらに処理の遅延時間が大きくなるという問題がある。
However, the method of linearly interpolating between pilot symbols in the time direction has a large influence of noise, and as shown in FIG. 7, the channel characteristics are IQ (In-phase-Quadrature phase) planes. In the case of changing so as to draw a curve, there is a problem that a deviation occurs between the estimated transmission line characteristic and the actual transmission line characteristic. FIG. 7 is a diagram illustrating the estimation result of the transmission path characteristics by interpolation.
In addition, when this method is used, a data symbol in between cannot be demodulated until a pilot symbol and the next pilot symbol are received, so a memory for storing received data is required, and the processing delay time is further increased. There is a problem.

ノイズの影響を少なくする伝送路特性推定方法として、複数のパイロットシンボルで平均をとることで伝送路特性の推定を行う方法がある。この方法を用いる場合、直線内挿補間による方法と比較して、ノイズの影響を低減することができる。しかし、この方法を用いる場合、復調の対象となるデータシンボルの前後の範囲でパイロットシンボルから求まる伝送路特性の平均をとる必要がある。そのため、直線内挿補間による方法と比べて、さらに伝送路特性の曲線的変動に追従することが難しくなってしまうという問題がある。また、受信データを記憶しておくメモリの容量も直線内挿補間による方法と比較して多くなり、処理の遅延時間が大きくなってしまうという問題がある。
また、メモリ容量の削減と遅延時間の減少のため、復調したいデータシンボルより過去のパイロットシンボルについてのみ平均をとる方法もあるが、この場合は、伝送路特性の時間変動に対する推定値の追従性が悪化してしまうという問題がある。
As a transmission path characteristic estimation method that reduces the influence of noise, there is a method of estimating transmission path characteristics by taking an average of a plurality of pilot symbols. When this method is used, the influence of noise can be reduced as compared with the method using linear interpolation. However, when this method is used, it is necessary to average the transmission path characteristics obtained from the pilot symbols in the range before and after the data symbol to be demodulated. Therefore, there is a problem that it becomes more difficult to follow the curve variation of the transmission path characteristics as compared with the method using linear interpolation. Further, the capacity of the memory for storing the received data is increased as compared with the method using linear interpolation, and there is a problem that the processing delay time is increased.
In order to reduce memory capacity and reduce delay time, there is a method of averaging only pilot symbols past data symbols to be demodulated. In this case, the follow-up of the estimated value with respect to temporal fluctuations in the channel characteristics is possible. There is a problem of getting worse.

本発明は上記の課題を解決するためになされたものであり、連続する複数の受信シンボルから次の受信シンボルの伝送路特性を推定する伝送路特性推定装置であって、前記受信シンボルの受信毎に、現在の伝送路特性を推定する現伝送路特性推定部と、前記受信シンボルの受信毎に、所定の平滑化係数を用いて前記現伝送特性推定部が推定した伝送路特性の移動平均値を算出する第1の平滑化部と、前記受信シンボルの受信毎に、前記第1の平滑化部の平滑化係数と異なる平滑化係数を用いて前記現伝送特性推定部が推定した伝送路特性の移動平均値を算出する第2の平滑化部と、前記現伝送路特性推定部が推定した伝送路特性と前記受信シンボルの前回の受信時に前記第1の平滑化部が算出した移動平均値との差を、前記受信シンボルの前回の受信時に前記第1の平滑化部が算出した移動平均値と前記受信シンボルの前回の受信時に前記第2の平滑化部が算出した移動平均値との差で除算する除算部と、前記第1の平滑化部が算出した移動平均値と前記第2の平滑化部が算出した移動平均値との差に、前記除算部の算出結果を乗算し、当該値に前記第1の平滑部が算出した移動平均値を加算することで、次の受信シンボルの伝送路特性を推定する次伝送路特性推定部とを備えることを特徴とする。   The present invention has been made in order to solve the above-described problem, and is a transmission path characteristic estimation apparatus that estimates a transmission path characteristic of a next reception symbol from a plurality of consecutive reception symbols. A current channel characteristic estimation unit for estimating the current channel characteristic, and a moving average value of the channel characteristic estimated by the current transmission characteristic estimation unit using a predetermined smoothing coefficient for each reception of the received symbol. And a channel characteristic estimated by the current transmission characteristic estimation unit using a smoothing coefficient different from the smoothing coefficient of the first smoothing unit for each reception of the received symbol. A second smoothing unit that calculates a moving average value of the first channel, a transmission channel characteristic estimated by the current transmission channel characteristic estimation unit, and a moving average value calculated by the first smoothing unit at the previous reception of the received symbol The difference between A division unit that divides by the difference between the moving average value calculated by the first smoothing unit at the time of reception and the moving average value calculated by the second smoothing unit at the previous reception of the received symbol; The difference between the moving average value calculated by the first smoothing unit and the moving average value calculated by the second smoothing unit is multiplied by the calculation result of the dividing unit, and the value is calculated by the first smoothing unit. And a next transmission path characteristic estimation unit that estimates the transmission path characteristic of the next received symbol by adding the calculated moving average value.

また、本発明は、連続する複数の受信シンボルから次の受信シンボルの伝送路特性を推定する伝送路特性推定装置を用いた伝送路特性推定方法であって、前記受信シンボルの受信毎に、現伝送路特性推定部は、現在の伝送路特性を推定し、第1の平滑化部は、所定の平滑化係数を用いて前記現伝送特性推定部が推定した伝送路特性の移動平均値を算出し、第2の平滑化部は、前記受信シンボルの受信毎に、前記第1の平滑化部の平滑化係数と異なる平滑化係数を用いて前記現伝送特性推定部が推定した伝送路特性の移動平均値を算出し、除算部は、前記現伝送路特性推定部が推定した伝送路特性と前記受信シンボルの前回の受信時に前記第1の平滑化部が算出した移動平均値との差を、前記受信シンボルの前回の受信時に前記第1の平滑化部が算出した移動平均値と前記受信シンボルの前回の受信時に前記第2の平滑化部が算出した移動平均値との差で除算し、次伝送路特性推定部は、前記第1の平滑化部が算出した移動平均値と前記第2の平滑化部が算出した移動平均値との差に、前記除算部の算出結果を乗算し、当該値に前記第1の平滑部が算出した移動平均値を加算することで、次の受信シンボルの伝送路特性を推定することを特徴とする。   The present invention also provides a transmission path characteristic estimation method using a transmission path characteristic estimation apparatus that estimates a transmission path characteristic of a next reception symbol from a plurality of consecutive reception symbols, and each time a reception symbol is received, The transmission line characteristic estimation unit estimates a current transmission line characteristic, and the first smoothing unit calculates a moving average value of the transmission line characteristic estimated by the current transmission characteristic estimation unit using a predetermined smoothing coefficient. The second smoothing unit uses the smoothing coefficient different from the smoothing coefficient of the first smoothing unit for each reception of the received symbol, to determine the channel characteristic estimated by the current transmission characteristic estimating unit. A moving average value is calculated, and the division unit calculates a difference between the transmission path characteristic estimated by the current transmission path characteristic estimation unit and the moving average value calculated by the first smoothing unit at the previous reception of the received symbol. The first smoothing unit at the previous reception of the received symbol Dividing by the difference between the calculated moving average value and the moving average value calculated by the second smoothing unit at the previous reception of the received symbol, the next transmission path characteristic estimating unit is configured by the first smoothing unit. The difference between the calculated moving average value and the moving average value calculated by the second smoothing unit is multiplied by the calculation result of the dividing unit, and the moving average value calculated by the first smoothing unit is multiplied by the value. The transmission path characteristics of the next received symbol are estimated by addition.

また、本発明は、連続する複数の受信シンボルから次の受信シンボルの伝送路特性を推定する伝送路特性推定装置を、前記受信シンボルの受信毎に現在の伝送路特性を推定する現伝送路特性推定部、前記受信シンボルの受信毎に所定の平滑化係数を用いて前記現伝送特性推定部が推定した伝送路特性の移動平均値を算出する第1の平滑化部、前記受信シンボルの受信毎に前記第1の平滑化部の平滑化係数と異なる平滑化係数を用いて前記現伝送特性推定部が推定した伝送路特性の移動平均値を算出する第2の平滑化部、前記現伝送路特性推定部が推定した伝送路特性と前記受信シンボルの前回の受信時に前記第1の平滑化部が算出した移動平均値との差を、前記受信シンボルの前回の受信時に前記第1の平滑化部が算出した移動平均値と前記受信シンボルの前回の受信時に前記第2の平滑化部が算出した移動平均値との差で除算する除算部、前記第1の平滑化部が算出した移動平均値と前記第2の平滑化部が算出した移動平均値との差に、前記除算部の算出結果を乗算し、当該値に前記第1の平滑部が算出した移動平均値を加算することで、次の受信シンボルの伝送路特性を推定する次伝送路特性推定部として動作させるためのプログラムである。   In addition, the present invention provides a transmission path characteristic estimation device that estimates a transmission path characteristic of a next reception symbol from a plurality of consecutive reception symbols, and a current transmission path characteristic that estimates a current transmission path characteristic for each reception of the reception symbol. An estimation unit, a first smoothing unit that calculates a moving average value of a transmission path characteristic estimated by the current transmission characteristic estimation unit using a predetermined smoothing coefficient for each reception of the reception symbol, and for each reception of the reception symbol A second smoothing unit for calculating a moving average value of the transmission path characteristics estimated by the current transmission characteristic estimation unit using a smoothing coefficient different from the smoothing coefficient of the first smoothing unit, and the current transmission path The difference between the transmission path characteristic estimated by the characteristic estimation unit and the moving average value calculated by the first smoothing unit at the previous reception of the received symbol is expressed as the first smoothing at the previous reception of the received symbol. The moving average value calculated by the A division unit that divides by a difference from the moving average value calculated by the second smoothing unit at the previous reception of a reception symbol, the moving average value calculated by the first smoothing unit, and the second smoothing unit Is multiplied by the calculation result of the division unit, and the moving average value calculated by the first smoothing unit is added to the value, whereby the transmission path characteristic of the next received symbol is calculated. Is a program for operating as a next transmission path characteristic estimation unit.

本発明によれば、除算部の処理によって次の受信シンボルにおける伝送路特性の推定に用いる係数を算出し、次伝送路特性推定部は、当該係数を用いて次の受信シンボルにおける伝送路特性を推定する。これにより、伝送路特性推定装置は、現在の伝送路特性及び過去の伝送路特性の移動平均値に基づいて精度良く次の受信シンボルにおける伝送路特性を推定することができる。   According to the present invention, the coefficient used for estimating the channel characteristic in the next received symbol is calculated by the processing of the division unit, and the next channel characteristic estimating unit uses the coefficient to calculate the channel characteristic in the next received symbol. presume. Thereby, the transmission path characteristic estimation apparatus can estimate the transmission path characteristic in the next received symbol with high accuracy based on the moving average value of the current transmission path characteristic and the past transmission path characteristic.

本発明の一実施形態による伝送路特性推定装置を搭載する復調装置の構成を示す概略ブロック図である。It is a schematic block diagram which shows the structure of the demodulation apparatus carrying the transmission line characteristic estimation apparatus by one Embodiment of this invention. OFDMの送信フレームのパターンの一例を示す図である。It is a figure which shows an example of the pattern of the transmission frame of OFDM. データシンボルに対する復調装置の動作を示すフローチャートである。It is a flowchart which shows operation | movement of the demodulation apparatus with respect to a data symbol. パイロットシンボルに対する復調装置の動作を示すフローチャートである。It is a flowchart which shows operation | movement of the demodulation apparatus with respect to a pilot symbol. 本実施形態による復調装置による伝送路特性の推定結果を示す図である。It is a figure which shows the estimation result of the transmission line characteristic by the demodulation apparatus by this embodiment. OFDMの送信フレームのパターンの他の例を示す図である。It is a figure which shows the other example of the pattern of the transmission frame of OFDM. 内挿補間による伝送路特性の推定結果を示す図である。It is a figure which shows the estimation result of the transmission line characteristic by interpolation.

以下、図面を参照しながら本発明の実施形態について詳しく説明する。
図1は、本発明の一実施形態による伝送路特性推定装置を搭載する復調装置100の構成を示す概略ブロック図である。
復調装置100は、OFDMフレームを受信してFFT(高速フーリエ変換:Fast Fourier Transform)処理を行う受信装置からFFTの出力データ(受信シンボル)の入力を受け付け、各サブキャリアに対する復調処理を行う装置である。復調装置100は、パイロット複素除算部101、伝送路特性記憶部102、位相振幅補償部103、デマップ部104、リマップ部105、複素除算部106、傾き補正部107、スイッチ部108、IIR(無限インパルス応答: Infinite Impulse Response)フィルタ109(第2の平滑化部)、IIRフィルタ110(第1の平滑化部)、移動平均値バッファ111、112、係数算出部113(除算部)、係数記憶部114、次伝送路特性推定部115、伝送路特性バッファ116を備える。
Hereinafter, embodiments of the present invention will be described in detail with reference to the drawings.
FIG. 1 is a schematic block diagram showing a configuration of a demodulation device 100 equipped with a transmission path characteristic estimation device according to an embodiment of the present invention.
The demodulator 100 is an apparatus that receives input of FFT output data (received symbols) from a receiver that receives an OFDM frame and performs FFT (Fast Fourier Transform) processing, and performs demodulation processing on each subcarrier. is there. Demodulator 100 includes pilot complex division unit 101, transmission path characteristic storage unit 102, phase amplitude compensation unit 103, demapping unit 104, remapping unit 105, complex division unit 106, slope correction unit 107, switch unit 108, IIR (infinite impulse). Response: Infinite Impulse Response) filter 109 (second smoothing unit), IIR filter 110 (first smoothing unit), moving average value buffers 111 and 112, coefficient calculation unit 113 (division unit), coefficient storage unit 114 , A next transmission path characteristic estimation unit 115 and a transmission path characteristic buffer 116 are provided.

パイロット複素除算部101は、受信シンボルがパイロットシンボルである場合に、受信したパイロットシンボルを、既知のパイロットシンボルモデルで複素除算することで、伝送路特性を算出する。なお、本実施形態において、受信シンボルに対する複素数計算は、受信シンボルのI成分を実数成分として扱い、Q成分を虚数成分として扱うことで行う。
伝送路特性記憶部102は、パイロット複素除算部101によって算出された伝送路特性を記憶する。なお、伝送路特性記憶部102が記憶する伝送路特性は、パイロット複素除算部101によって伝送路特性が算出されるたびに更新される。
When the received symbol is a pilot symbol, pilot complex division section 101 performs complex division on the received pilot symbol by a known pilot symbol model to calculate transmission path characteristics. In the present embodiment, the complex number calculation for the received symbol is performed by treating the I component of the received symbol as a real component and the Q component as an imaginary component.
The transmission path characteristic storage unit 102 stores the transmission path characteristic calculated by the pilot complex division unit 101. Note that the transmission path characteristics stored in the transmission path characteristics storage unit 102 are updated each time the transmission path characteristics are calculated by the pilot complex division unit 101.

位相振幅補償部103は、受信シンボルがデータシンボルである場合に、受信したデータシンボルを、伝送路特性バッファ116が記憶する伝送路特性で除算することで、送信シンボルを再現する。なお、伝送路特性バッファ116が記憶する伝送路特性は、前回受信シンボルを受信したときに次伝送路特性推定部115によって推定された、受信シンボルの次の受信時における伝送路特性、すなわち現在の伝送路特性である。
デマップ部104は、位相振幅補償部103によって再現された送信シンボルを所定の位相振幅変調方式に従って「0」または「1」のビットデータを出力する。なお、位相振幅変調方式の例としては、QPSK(四位相偏移変調:Quaternary Phase Shift Keying)や16QAM(16値直角位相振幅変調:16 Quadrature Amplitude Modulation)などが挙げられる。
When the received symbol is a data symbol, the phase / amplitude compensation unit 103 reproduces the transmission symbol by dividing the received data symbol by the transmission path characteristic stored in the transmission path characteristic buffer 116. The transmission path characteristic stored in the transmission path characteristic buffer 116 is the transmission path characteristic at the time of the next reception of the received symbol, which is estimated by the next transmission path characteristic estimation unit 115 when the received symbol is received last time, that is, the current Transmission path characteristics.
The demapping unit 104 outputs bit data of “0” or “1” from the transmission symbol reproduced by the phase amplitude compensation unit 103 according to a predetermined phase amplitude modulation method. Examples of the phase amplitude modulation method include QPSK (Quadrature Phase Shift Keying) and 16QAM (16-Quadrature Amplitude Modulation).

リマップ部105は、デマップ部104によって出力されたビットデータから所定の位相振幅変調方式に基づいて送信シンボルを生成する。なお、位相振幅補償部103が再現した送信シンボルには、伝送路特性の推定誤差が含まれるが、リマップ部105によって生成された送信シンボルには誤差が含まれない。
複素除算部106は、受信シンボルがデータシンボルである場合に、リマップ部105によって生成された送信シンボルを受信したデータシンボルで複素除算することで、伝送路特性を算出する。
The remapping unit 105 generates a transmission symbol from the bit data output from the demapping unit 104 based on a predetermined phase amplitude modulation method. The transmission symbol reproduced by the phase amplitude compensation unit 103 includes an estimation error of transmission path characteristics, but the transmission symbol generated by the remapping unit 105 does not include an error.
When the received symbol is a data symbol, the complex division unit 106 calculates transmission path characteristics by performing complex division on the transmission symbol generated by the remapping unit 105 by the received data symbol.

傾き補正部107は、伝送路特性バッファ116が記憶する伝送路特性と、伝送路特性記憶部102が記憶する伝送路特性とに基づいてコンスタレーションが揃うように伝送路特性を補正する。すなわち、傾き補正部107は、受信器の初期動作時に、コンスタレーションが90度、180度、270度の何れかだけ傾いた状態ではまり込んだ場合に、それを訂正するように、伝送路特性を補正する。具体的には、傾き補正部107は、以下に示す式(1)を用いて伝送路特性記憶部102が記憶する伝送路特性の複素共役と、伝送路特性バッファ116が記憶する伝送路特性との複素乗算の時間平均を算出する。   The inclination correction unit 107 corrects the transmission path characteristics so that the constellation is aligned based on the transmission path characteristics stored in the transmission path characteristic buffer 116 and the transmission path characteristics stored in the transmission path characteristic storage unit 102. In other words, the tilt correction unit 107 corrects the transmission line characteristics so that, when the constellation is caught in a state tilted at any one of 90 degrees, 180 degrees, and 270 degrees during the initial operation of the receiver, it is corrected. Correct. Specifically, the inclination correction unit 107 uses the following equation (1) to determine the complex conjugate of the transmission path characteristic stored in the transmission path characteristic storage unit 102 and the transmission path characteristic stored in the transmission path characteristic buffer 116. Calculate the time average of complex multiplication.

Figure 2012222487
Figure 2012222487

但し、h(n)は、受信シンボルのnシンボル目において伝送路特性記憶部102が記憶する伝送路特性である。また、上付き傍線はその値の複素共役を示す記号である。また、H(n)は、受信シンボルのnシンボル目において伝送路特性バッファ116が記憶する伝送路特性である。 Here, h p (n) is a transmission path characteristic stored in the transmission path characteristic storage unit 102 at the nth symbol of the received symbol. The superscript side line is a symbol indicating the complex conjugate of the value. H (n) is a transmission path characteristic stored in the transmission path characteristic buffer 116 at the nth symbol of the received symbol.

そして、傾き補正部107は、上記式(1)によって算出された複素数と実数軸とのなす角度が±45度の範囲内である場合(コンスタレーションの傾きが正しい場合)に、複素除算部106が算出した伝送路特性をスイッチ部108に出力し、±45度の範囲外である場合に、伝送路特性記憶部102が記憶する伝送路特性をスイッチ部108に出力する(伝送路特性を補正する)。
スイッチ部108は、受信シンボルがパイロットシンボルである場合に、パイロット複素除算部101の算出結果を選択してIIRフィルタ109、110に出力し、受信シンボルがデータシンボルである場合に、傾き補正部107の出力を選択してIIRフィルタ109、110に出力する。
なお、パイロット複素除算部101、伝送路特性記憶部102、複素除算部106、傾き補正部107、スイッチ部108は、協働することで、受信シンボルの受信毎に現在の伝送路特性を推定する現伝送路特性推定部として動作する。
Then, when the angle between the complex number calculated by the equation (1) and the real number axis is within a range of ± 45 degrees (when the constellation inclination is correct), the inclination correcting unit 107 Is output to the switch unit 108, and when it is out of the range of ± 45 degrees, the transmission channel characteristic stored in the transmission channel characteristic storage unit 102 is output to the switch unit 108 (correcting the transmission channel characteristic is corrected). To do).
When the received symbol is a pilot symbol, the switch unit 108 selects the calculation result of the pilot complex division unit 101 and outputs it to the IIR filters 109 and 110. When the received symbol is a data symbol, the switch 108 corrects the inclination. Is output to the IIR filters 109 and 110.
Pilot complex division section 101, transmission path characteristic storage section 102, complex division section 106, slope correction section 107, and switch section 108 cooperate to estimate the current transmission path characteristics every time a received symbol is received. It operates as a current transmission path characteristic estimation unit.

IIRフィルタ109、110は、所定の忘却係数を用いて入力値の指数移動平均値を算出する一次IIRフィルタである。なお、IIRフィルタ109とIIRフィルタ110は、それぞれ異なる忘却係数を用いて指数移動平均値を算出する。IIRフィルタ109、110の忘却係数には、1/2、1/4、1/8、1/16、1/32、1/64、1/128のうち異なる2つの値が用いられる。
具体的には、IIRフィルタ109は、以下に示す式(2)によって指数移動平均値を算出し、IIRフィルタ110は、以下に示す式(3)によって指数移動平均値を算出する。
The IIR filters 109 and 110 are first-order IIR filters that calculate an exponential moving average value of input values using a predetermined forgetting factor. The IIR filter 109 and the IIR filter 110 calculate exponential moving average values using different forgetting factors. Two different values of 1/2, 1/4, 1/8, 1/16, 1/32, 1/64, and 1/128 are used as the forgetting coefficients of the IIR filters 109 and 110.
Specifically, the IIR filter 109 calculates an exponential moving average value by the following equation (2), and the IIR filter 110 calculates an exponential moving average value by the following equation (3).

Figure 2012222487
Figure 2012222487
Figure 2012222487
Figure 2012222487

但し、h(n)は、受信シンボルのnシンボル目においてIIRフィルタ109が算出する伝送路特性の指数移動平均値である。また、aは、IIRフィルタ109が用いる忘却係数である。また、hは、スイッチ部108が出力した伝送路特性である。なお、h(n−1)は、受信シンボルのn−1シンボル目においてIIRフィルタ109が算出した伝送路特性の指数移動平均値であり、受信シンボルのnシンボル目において移動平均値バッファ111に記憶されている値である。 Here, h a (n) is an exponential moving average value of transmission path characteristics calculated by the IIR filter 109 at the nth symbol of the received symbol. Further, a is a forgetting factor used by the IIR filter 109. Further, h is a transmission path characteristic output from the switch unit 108. Here, h a (n−1) is an exponential moving average value of the transmission path characteristic calculated by the IIR filter 109 at the n−1th symbol of the received symbol, and is stored in the moving average value buffer 111 at the nth symbol of the received symbol. It is a stored value.

また、h(n)は、受信シンボルのnシンボル目においてIIRフィルタ110が算出する伝送路特性の指数移動平均値である。また、bは、IIRフィルタ110が用いる忘却係数である。なお、h(n−1)は、受信シンボルのn−1シンボル目においてIIRフィルタ110が算出した伝送路特性の指数移動平均値であり、受信シンボルのnシンボル目において移動平均値バッファ112に記憶されている値である。
移動平均値バッファ111、112は、それぞれIIRフィルタ109、110が受信シンボルの前回の受信時に算出した伝送路特性の指数移動平均値を記憶する。
Further, h b (n) is an exponential moving average value of transmission path characteristics calculated by the IIR filter 110 at the nth symbol of the received symbol. B is a forgetting factor used by the IIR filter 110. Note that h b (n−1) is an exponential moving average value of the channel characteristics calculated by the IIR filter 110 at the n−1th symbol of the received symbol, and is stored in the moving average value buffer 112 at the nth symbol of the received symbol. It is a stored value.
Moving average value buffers 111 and 112 store exponential moving average values of transmission path characteristics calculated by IIR filters 109 and 110, respectively, at the previous reception of received symbols.

係数算出部113は、パイロット複素除算部101が算出した伝送路特性h(n)と、移動平均値バッファ111、112が記憶する伝送路特性の指数移動平均値h(n−1)及びh(n−1)とに基づいて、次の受信シンボルにおける伝送路特性の推定に用いる係数αを算出する。具体的には、係数算出部113は、以下に示す式(4)、式(5)の値を算出する。 Coefficient calculation unit 113, a transmission path characteristic h p for the pilot complex division section 101 is calculated (n), exponential moving average of the channel characteristics moving average buffer 111 stores values h a (n-1) and Based on h b (n−1), a coefficient α used to estimate the transmission path characteristics in the next received symbol is calculated. Specifically, the coefficient calculation unit 113 calculates the values of the following expressions (4) and (5).

Figure 2012222487
Figure 2012222487
Figure 2012222487
Figure 2012222487

そして、係数算出部113は、式(4)、式(5)の算出結果の時間平均値を算出し、式(6)に基づいて係数αを算出する。   And the coefficient calculation part 113 calculates the time average value of the calculation result of Formula (4) and Formula (5), and calculates coefficient (alpha) based on Formula (6).

Figure 2012222487
Figure 2012222487

但し、average{}は、カッコ内の値の時間平均値である。なお、時間平均値の算出方法には、単純移動平均、加重移動平均(FIR(有限インパルス応答: Finite Impulse Response)フィルタ平均)、指数移動平均(IIRフィルタ平均)等があるが、例えば単純移動平均を用いる場合、以下に示す式(7)によって係数αを算出することができる。   However, average {} is a time average value of the values in parentheses. Note that the time average value calculation method includes a simple moving average, a weighted moving average (FIR (Finite Impulse Response) filter average), an exponential moving average (IIR filter average), and the like. Is used, the coefficient α can be calculated by the following equation (7).

Figure 2012222487
Figure 2012222487

係数記憶部114は、係数算出部113が算出した係数αを記憶する。なお、係数記憶部114が記憶する係数は、係数算出部113によって係数αが算出されるたびに更新される。
次伝送路特性推定部115は、係数記憶部114が記憶する係数αと、IIRフィルタ109、110が算出した伝送路特性の指数移動平均値h(n)及びh(n)とに基づいて、次の受信シンボルにおける伝送路特性を推定する。具体的には、次伝送路特性推定部115は、以下に示す式(8)を算出することで、次の受信シンボルにおける伝送路特性H(n+1)を推定する。
The coefficient storage unit 114 stores the coefficient α calculated by the coefficient calculation unit 113. The coefficient stored in the coefficient storage unit 114 is updated every time the coefficient α is calculated by the coefficient calculation unit 113.
The next transmission line characteristic estimation unit 115 is based on the coefficient α stored in the coefficient storage unit 114 and the exponential moving average values h a (n) and h b (n) of the transmission line characteristics calculated by the IIR filters 109 and 110. Thus, the transmission path characteristic in the next received symbol is estimated. Specifically, the next transmission path characteristic estimation unit 115 estimates the transmission path characteristic H (n + 1) in the next received symbol by calculating the following equation (8).

Figure 2012222487
Figure 2012222487

伝送路特性バッファ116は、次伝送路特性推定部115によって推定された伝送路特性を記憶する。   The transmission path characteristic buffer 116 stores the transmission path characteristic estimated by the next transmission path characteristic estimation unit 115.

上述した構成を備えることで、復調装置100は、現在の伝送路特性及び過去の伝送路特性の移動平均値に基づいて精度良く次の受信シンボルにおける伝送路特性を推定することができる。   With the above-described configuration, the demodulator 100 can accurately estimate the transmission path characteristic of the next received symbol based on the moving average value of the current transmission path characteristic and the past transmission path characteristic.

ここで、係数算出部113が式(6)を用いて係数αを算出する理由について説明する。
次の受信シンボルにおける伝送路特性H(n+1)を、式(8)によって表す場合、既知のパイロットシンボルから求めた伝送路特性h(n)は、以下に示す式(9)によって表すことができる。
Here, the reason why the coefficient calculation unit 113 calculates the coefficient α using Expression (6) will be described.
When the transmission path characteristic H (n + 1) in the next received symbol is expressed by Expression (8), the transmission path characteristic h p (n) obtained from a known pilot symbol can be expressed by Expression (9) shown below. it can.

Figure 2012222487
Figure 2012222487

ここで、式(9)を等式変形することで、係数αは以下に示す式(10)によって表すことができる。   Here, the coefficient α can be expressed by the following equation (10) by performing equation transformation on the equation (9).

Figure 2012222487
Figure 2012222487

ところで、推定された伝送路特性に含まれるノイズの影響を軽減するためには、時間平均値を用いて係数αを算出することが好ましい。そこで、式(10)の分母が0になるのを防ぐため、分母が正の実数になるように、式(10)の分母と分子にそれぞれ分母の値の複素共役を乗算する。これにより、式(6)が得られる。   Incidentally, in order to reduce the influence of noise included in the estimated transmission path characteristics, it is preferable to calculate the coefficient α using a time average value. Therefore, in order to prevent the denominator of Equation (10) from becoming 0, the denominator and numerator of Equation (10) are each multiplied by a complex conjugate of the denominator value so that the denominator becomes a positive real number. Thereby, Formula (6) is obtained.

以下、復調装置100の動作について説明する。
図2は、OFDMの送信フレームのパターンの一例を示す図である。
図2に示すように、本実施形態において送信装置は、送信フレームを送信する際、データシンボルの12シンボルごとにパイロット信号を1つ挿入する。したがって、復調装置100は、受信装置から入力を受け付けた受信シンボルの12シンボル毎に、パイロットシンボルに基づく伝送路特性推定を行い、それ以外においてデータシンボルの復調処理及び伝送路特性推定を行う。
Hereinafter, the operation of the demodulation device 100 will be described.
FIG. 2 is a diagram illustrating an example of an OFDM transmission frame pattern.
As shown in FIG. 2, in the present embodiment, when transmitting a transmission frame, the transmission apparatus inserts one pilot signal for every 12 data symbols. Therefore, demodulating apparatus 100 performs channel characteristic estimation based on pilot symbols for every 12 received symbols received from the receiving apparatus, and otherwise performs data symbol demodulation processing and channel characteristic estimation.

まず、復調装置100によるデータシンボルの復調処理及び伝送路特性推定の処理について説明する。
図3は、データシンボルに対する復調装置100の動作を示すフローチャートである。
まず、受信装置からデータシンボルの入力を受け付けると、位相振幅補償部103は、入力を受け付けたデータシンボルを伝送路特性バッファ116が記憶する伝送路特性で除算することで、送信シンボルを再現する(ステップS1)。次に、デマップ部104は、位相振幅補償部103によって再現された送信シンボルから、所定の振幅位相変調方式に従って硬判定を行い、ビットデータを生成する(ステップS2)。なお、当該ビットデータは、復調結果として処理装置に出力される。
First, data symbol demodulation processing and transmission path characteristic estimation processing by the demodulator 100 will be described.
FIG. 3 is a flowchart showing the operation of the demodulator 100 for data symbols.
First, when receiving an input of a data symbol from the receiving device, the phase amplitude compensation unit 103 reproduces the transmission symbol by dividing the data symbol that has received the input by the transmission path characteristic stored in the transmission path characteristic buffer 116 ( Step S1). Next, the demapping unit 104 performs a hard decision from the transmission symbol reproduced by the phase amplitude compensation unit 103 according to a predetermined amplitude phase modulation method, and generates bit data (step S2). The bit data is output to the processing device as a demodulation result.

次に、リマップ部105は所定の振幅位相変調方式に基づいて、デマップ部104によって生成されたビットデータから送信シンボルを生成する(ステップS3)。次に、複素除算部106は、入力を受け付けたデータシンボルをリマップ部105が生成した送信シンボルで複素除算することで、伝送路特性を推定する(ステップS4)。次に、傾き補正部107は、伝送路特性記憶部102が記憶する伝送路特性と伝送路特性バッファ116が記憶する伝送路特性とを用いて上記式(1)を演算する(ステップS5)。次に、傾き補正部107は、式(1)の演算結果と実数軸とがなす角の角度が±45度の範囲内であるか否かを判定する(ステップS6)。   Next, the remapping unit 105 generates a transmission symbol from the bit data generated by the demapping unit 104 based on a predetermined amplitude phase modulation method (step S3). Next, the complex division unit 106 estimates transmission path characteristics by performing complex division on the data symbol that has received the input by the transmission symbol generated by the remapping unit 105 (step S4). Next, the inclination correction unit 107 calculates the above equation (1) using the transmission path characteristics stored in the transmission path characteristics storage section 102 and the transmission path characteristics stored in the transmission path characteristics buffer 116 (step S5). Next, the inclination correction unit 107 determines whether or not the angle formed by the calculation result of Expression (1) and the real number axis is within a range of ± 45 degrees (step S6).

式(1)の演算結果と実数軸とがなす角の角度が±45度の範囲内である場合(ステップS6:YES)、傾き補正部107は、複素除算部106が算出した伝送路特性をスイッチ部108に出力する(ステップS7)。他方、式(1)の演算結果と実数軸とがなす角の角度が±45度の範囲外である場合(ステップS6:NO)、傾き補正部107は、伝送路特性記憶部102が記憶する伝送路特性をスイッチ部108に出力する(ステップS8)。   When the angle formed by the calculation result of Expression (1) and the real axis is within a range of ± 45 degrees (step S6: YES), the inclination correction unit 107 calculates the transmission path characteristics calculated by the complex division unit 106. It outputs to the switch part 108 (step S7). On the other hand, when the angle formed by the calculation result of Expression (1) and the real number axis is outside the range of ± 45 degrees (step S6: NO), the inclination correction unit 107 stores the transmission path characteristic storage unit 102. The transmission path characteristics are output to the switch unit 108 (step S8).

スイッチ部108は、ステップS7またはステップS8で伝送路特性の入力を受け付けると、当該伝送路特性をIIRフィルタ109、110に出力する。IIRフィルタ109、110は、それぞれ式(2)、式(3)に基づいて伝送路特性の指数移動平均値を算出する(ステップS9)。このとき、IIRフィルタ109、110は、それぞれ受信シンボルの前回の受信時に算出した伝送路特性の指数移動平均値を移動平均値バッファ111、112に記録する(ステップS10)。
そして、次伝送路特性推定部115は、IIRフィルタ109、110が算出した伝送路特性の指数移動平均値と係数記憶部114が記憶する係数αとを用いて、式(8)に基づいて次の受信シンボルにおける伝送路特性を推定する(ステップS11)。次伝送路特性推定部115は、推定した伝送路特性を伝送路特性バッファ116に記録する(ステップS12)。当該伝送路特性は、次のデータシンボルの受信時に用いられる。
When the switch unit 108 receives the input of the transmission path characteristic in step S7 or step S8, the switching unit 108 outputs the transmission path characteristic to the IIR filters 109 and 110. The IIR filters 109 and 110 calculate the exponential moving average value of the transmission path characteristics based on the expressions (2) and (3), respectively (step S9). At this time, the IIR filters 109 and 110 record the exponential moving average values of the transmission path characteristics calculated at the previous reception of the received symbols in the moving average value buffers 111 and 112, respectively (step S10).
Then, the next transmission path characteristic estimation unit 115 uses the exponential moving average value of the transmission path characteristics calculated by the IIR filters 109 and 110 and the coefficient α stored in the coefficient storage unit 114 based on the equation (8). The transmission path characteristics in the received symbols are estimated (step S11). The next transmission line characteristic estimation unit 115 records the estimated transmission line characteristic in the transmission line characteristic buffer 116 (step S12). The transmission path characteristic is used when the next data symbol is received.

まず、復調装置100によるパイロットシンボルに基づく伝送路特性推定の処理について説明する。
図4は、パイロットシンボルに対する復調装置100の動作を示すフローチャートである。
まず、受信装置からパイロットシンボルの入力を受け付けると、パイロット複素除算部101は、入力を受け付けたパイロットシンボルを予め定められたパイロットシンボルモデルで除算することで、伝送路特性を推定する(ステップS13)。パイロット複素除算部101は、推定した伝送路特性を伝送路特性記憶部102に記録する(ステップS14)。
First, transmission path characteristic estimation processing based on pilot symbols by demodulator 100 will be described.
FIG. 4 is a flowchart showing the operation of the demodulation apparatus 100 for pilot symbols.
First, upon receiving an input of a pilot symbol from the receiving apparatus, pilot complex division section 101 estimates a transmission path characteristic by dividing the pilot symbol that has received the input by a predetermined pilot symbol model (step S13). . The pilot complex division unit 101 records the estimated transmission path characteristic in the transmission path characteristic storage unit 102 (step S14).

また、パイロット複素除算部101は、推定した伝送路特性をスイッチ部108に出力する。スイッチ部108は、伝送路特性の入力を受け付けると、当該伝送路特性をIIRフィルタ109、110に出力する。IIRフィルタ109、110は、それぞれ式(2)、式(3)に基づいて伝送路特性の指数移動平均値を算出する(ステップS15)。このとき、IIRフィルタ109、110は、それぞれ受信シンボルの前回の受信時に算出した伝送路特性の指数移動平均値を移動平均値バッファ111、112に記録する(ステップS16)。   Pilot complex division section 101 also outputs the estimated transmission path characteristics to switch section 108. When receiving an input of the transmission path characteristic, the switch unit 108 outputs the transmission path characteristic to the IIR filters 109 and 110. The IIR filters 109 and 110 calculate the exponential moving average value of the transmission path characteristics based on the equations (2) and (3), respectively (step S15). At this time, the IIR filters 109 and 110 record the exponential moving average values of the transmission path characteristics calculated at the previous reception of the received symbols in the moving average value buffers 111 and 112, respectively (step S16).

また、係数算出部113は、パイロット複素除算部101が推定した伝送路特性と移動平均値バッファ111、112が記憶する受信シンボルの前回の受信時における伝送路特性の指数移動平均値とを用いて、式(6)に基づいて、次の受信シンボルにおける伝送路特性の推定に用いる係数αを算出する(ステップS17)。係数算出部113は、算出した係数αを係数記憶部114に記録する(ステップS18)。   Further, the coefficient calculation unit 113 uses the transmission path characteristic estimated by the pilot complex division unit 101 and the exponential moving average value of the transmission path characteristic at the previous reception time of the received symbol stored in the moving average value buffers 111 and 112. Based on the equation (6), the coefficient α used for estimating the channel characteristics in the next received symbol is calculated (step S17). The coefficient calculation unit 113 records the calculated coefficient α in the coefficient storage unit 114 (step S18).

そして、次伝送路特性推定部115は、IIRフィルタ109、110が算出した伝送路特性の指数移動平均値と係数記憶部114が記憶する係数αとを用いて、式(8)に基づいて次の受信シンボルにおける伝送路特性を推定する(ステップS19)。次伝送路特性推定部115は、推定した伝送路特性を伝送路特性バッファ116に記録する(ステップS20)。当該伝送路特性は、次のデータシンボルの受信時に用いられる。   Then, the next transmission path characteristic estimation unit 115 uses the exponential moving average value of the transmission path characteristics calculated by the IIR filters 109 and 110 and the coefficient α stored in the coefficient storage unit 114 based on the equation (8). The transmission path characteristics in the received symbols are estimated (step S19). The next transmission line characteristic estimation unit 115 records the estimated transmission line characteristic in the transmission line characteristic buffer 116 (step S20). The transmission path characteristic is used when the next data symbol is received.

このように、本実施形態によれば、IIRフィルタ109、110を備えることでノイズの影響を低減する。さらに、本実施形態によれば、忘却係数(平滑化係数)の異なるIIRフィルタ109、110の出力の差から現在の伝送路特性を推定する処理を行うことで、IIRフィルタ109、110による伝送路特性の時間追従性の悪化を改善することができる。   Thus, according to the present embodiment, the influence of noise is reduced by providing the IIR filters 109 and 110. Furthermore, according to the present embodiment, the transmission path by the IIR filters 109 and 110 is performed by performing the process of estimating the current transmission path characteristics from the difference between the outputs of the IIR filters 109 and 110 having different forgetting coefficients (smoothing coefficients). It is possible to improve the deterioration of the time followability of characteristics.

図5は、本実施形態による復調装置100による伝送路特性の推定結果を示す図である。
本実施形態によれば、図5に示すように、伝送路特性が複素平面上で曲線の軌跡を描く場合であっても、係数αの虚数部が適切に調整されることによって、伝送路特性が正確に推定されることとなる。
FIG. 5 is a diagram illustrating the estimation result of the transmission path characteristics by the demodulator 100 according to the present embodiment.
According to the present embodiment, as shown in FIG. 5, even when the transmission path characteristic draws a curved locus on the complex plane, the transmission path characteristic is adjusted by appropriately adjusting the imaginary part of the coefficient α. Is accurately estimated.

また、本実施形態によれば、データシンボル受信時に、既にそのときの伝送路特性が推定されているため、復調装置100は、FFT出力信号のデータを蓄積するための記憶部を備える必要がない。   Further, according to the present embodiment, when the data symbol is received, the transmission path characteristics at that time are already estimated, and thus the demodulator 100 does not need to include a storage unit for accumulating data of the FFT output signal. .

以上、図面を参照してこの発明の一実施形態について詳しく説明してきたが、具体的な構成は上述のものに限られることはなく、この発明の要旨を逸脱しない範囲内において様々な設計変更等をすることが可能である。   As described above, the embodiment of the present invention has been described in detail with reference to the drawings. However, the specific configuration is not limited to the above, and various design changes and the like can be made without departing from the scope of the present invention. It is possible to

図6は、OFDMの送信フレームのパターンの他の例を示す図である。
例えば、本実施形態では、OFDMの送信フレームとして図2に示すように各サブキャリアにおいて同じタイミングでパイロットシンボルとデータシンボルが繰り返されるものを用いる場合を説明したが、これに限られない。例えば、OFDMの送信フレームとしては、図6(A)に示すように、各サブキャリアで異なるタイミングでパイロットシンボルとデータシンボルが繰り返されるものや、図6(B)に示すように、特定のサブキャリアのみにパイロットシンボルが配置されるものが挙げられる。
FIG. 6 is a diagram illustrating another example of an OFDM transmission frame pattern.
For example, in the present embodiment, a case has been described where an OFDM transmission frame in which pilot symbols and data symbols are repeated at the same timing in each subcarrier as shown in FIG. 2 is not limited to this. For example, as an OFDM transmission frame, as shown in FIG. 6A, a pilot symbol and a data symbol are repeated at different timings in each subcarrier, or as shown in FIG. One in which pilot symbols are arranged only on the carrier can be mentioned.

なお、図6(B)に示すように特定のサブキャリアのみにパイロットシンボルが配置される場合は、パイロットシンボルのないサブキャリアは周波数方向または周波数−時間方向で補間すること(図6(B)の○部分)により、パイロットシンボルの代わりとして用いることで、本発明を実施することができる。   When pilot symbols are arranged only on specific subcarriers as shown in FIG. 6B, subcarriers without pilot symbols are interpolated in the frequency direction or frequency-time direction (FIG. 6B). By using as a substitute for a pilot symbol, the present invention can be implemented.

また、本実施形態では、式(6)を用いて係数αを算出する場合について説明したが、これに限られず、例えば、式(6)の分母及び分子に対する時間平均をとらずに係数αを算出するようにしても良い。但し、この場合は、式(6)を用いた場合と比較してノイズの影響を受けやすくなる。
また、分母が0となった場合の例外処理を用意しておくことで、式(6)の代わりに式(10)を用いて係数αを算出することもできる。
In this embodiment, the case where the coefficient α is calculated using the equation (6) has been described. However, the present invention is not limited to this. It may be calculated. However, in this case, it is more susceptible to noise as compared to the case of using Equation (6).
Also, by preparing an exception process when the denominator becomes 0, the coefficient α can be calculated using Expression (10) instead of Expression (6).

また、本実施形態では、係数算出部113が、パイロットシンボルの受信時にのみ係数αを算出する場合を説明したが、これに限られず、データシンボルの受信時に、スイッチ部108の出力を用いて係数αを算出するようにしても良い。但し、パイロットシンボルから推定される伝送路特性は、データシンボルから推定される伝送路特性より精度が高いため、パイロットシンボルの受信時にのみ係数αを算出したほうが、係数αの精度が高い。   Further, in the present embodiment, the case where the coefficient calculation unit 113 calculates the coefficient α only at the time of receiving a pilot symbol has been described. However, the present invention is not limited to this, and the coefficient is calculated using the output of the switch unit 108 when receiving a data symbol. α may be calculated. However, since the channel characteristics estimated from the pilot symbols have higher accuracy than the channel characteristics estimated from the data symbols, the accuracy of the coefficient α is higher when the coefficient α is calculated only when the pilot symbols are received.

上述の復調装置100は内部に、コンピュータシステムを有している。そして、上述した各処理部の動作は、プログラムの形式でコンピュータ読み取り可能な記録媒体に記憶されており、このプログラムをコンピュータが読み出して実行することによって、上記処理が行われる。ここでコンピュータ読み取り可能な記録媒体とは、磁気ディスク、光磁気ディスク、CD−ROM、DVD−ROM、半導体メモリ等をいう。また、このコンピュータプログラムを通信回線によってコンピュータに配信し、この配信を受けたコンピュータが当該プログラムを実行するようにしても良い。   The demodulator 100 described above has a computer system therein. The operation of each processing unit described above is stored in a computer-readable recording medium in the form of a program, and the above processing is performed by the computer reading and executing this program. Here, the computer-readable recording medium means a magnetic disk, a magneto-optical disk, a CD-ROM, a DVD-ROM, a semiconductor memory, or the like. Alternatively, the computer program may be distributed to the computer via a communication line, and the computer that has received the distribution may execute the program.

また、上記プログラムは、前述した機能の一部を実現するためのものであっても良い。さらに、前述した機能をコンピュータシステムにすでに記録されているプログラムとの組み合わせで実現できるもの、いわゆる差分ファイル(差分プログラム)であっても良い。   The program may be for realizing a part of the functions described above. Furthermore, what can implement | achieve the function mentioned above in combination with the program already recorded on the computer system, and what is called a difference file (difference program) may be sufficient.

100…復調装置 101…パイロット複素除算部 102…伝送路特性記憶部 103…位相振幅補償部 104…デマップ部 105…リマップ部 106…複素除算部 107…傾き補正部 108…スイッチ部 109、110…IIRフィルタ 111、112…移動平均値バッファ 113…係数算出部 114…係数記憶部 115…次伝送路特性推定部 116…伝送路特性バッファ   DESCRIPTION OF SYMBOLS 100 ... Demodulator 101 ... Pilot complex division part 102 ... Transmission path characteristic memory | storage part 103 ... Phase amplitude compensation part 104 ... Demap part 105 ... Remap part 106 ... Complex division part 107 ... Inclination correction part 108 ... Switch part 109, 110 ... IIR Filters 111, 112 ... Moving average value buffer 113 ... Coefficient calculation unit 114 ... Coefficient storage unit 115 ... Next transmission line characteristic estimation part 116 ... Transmission line characteristic buffer

Claims (6)

連続する複数の受信シンボルから次の受信シンボルの伝送路特性を推定する伝送路特性推定装置であって、
前記受信シンボルの受信毎に現在の伝送路特性を推定する現伝送路特性推定部と、
前記受信シンボルの受信毎に所定の平滑化係数を用いて前記現伝送特性推定部が推定した伝送路特性の移動平均値を算出する第1の平滑化部と、
前記受信シンボルの受信毎に前記第1の平滑化部の平滑化係数と異なる平滑化係数を用いて前記現伝送特性推定部が推定した伝送路特性の移動平均値を算出する第2の平滑化部と、
前記現伝送路特性推定部が推定した伝送路特性と前記受信シンボルの前回の受信時に前記第1の平滑化部が算出した移動平均値との差を、前記受信シンボルの前回の受信時に前記第1の平滑化部が算出した移動平均値と前記受信シンボルの前回の受信時に前記第2の平滑化部が算出した移動平均値との差で除算する除算部と、
前記第1の平滑化部が算出した移動平均値と前記第2の平滑化部が算出した移動平均値との差に、前記除算部の算出結果を乗算し、当該値に前記第1の平滑部が算出した移動平均値を加算することで、次の受信シンボルの伝送路特性を推定する次伝送路特性推定部と
を備えることを特徴とする伝送路特性推定装置。
A transmission path characteristic estimation device that estimates a transmission path characteristic of a next reception symbol from a plurality of consecutive reception symbols,
A current channel characteristic estimator for estimating a current channel characteristic every time the received symbol is received;
A first smoothing unit that calculates a moving average value of transmission path characteristics estimated by the current transmission characteristic estimation unit using a predetermined smoothing coefficient for each reception of the received symbol;
Second smoothing for calculating a moving average value of the transmission path characteristic estimated by the current transmission characteristic estimation unit using a smoothing coefficient different from the smoothing coefficient of the first smoothing unit for each reception of the received symbol And
The difference between the transmission path characteristic estimated by the current transmission path characteristic estimation unit and the moving average value calculated by the first smoothing unit at the previous reception of the received symbol is expressed as the first time at the previous reception of the reception symbol. A division unit that divides by the difference between the moving average value calculated by one smoothing unit and the moving average value calculated by the second smoothing unit at the previous reception of the received symbol;
The difference between the moving average value calculated by the first smoothing unit and the moving average value calculated by the second smoothing unit is multiplied by the calculation result of the dividing unit, and the value is multiplied by the first smoothing value. A transmission path characteristic estimation device, comprising: a next transmission path characteristic estimation unit that estimates a transmission path characteristic of a next received symbol by adding the moving average values calculated by the unit.
前記除算部は、受信した受信シンボルがパイロットシンボルであるときに計算を行う
ことを特徴とする請求項1に記載の伝送路特性推定装置。
The transmission path characteristic estimation apparatus according to claim 1, wherein the division unit performs calculation when the received symbol received is a pilot symbol.
前記除算部は、前記現伝送路特性推定部が推定した伝送路特性と前記受信シンボルの前回の受信時に前記第1の平滑化部が算出した移動平均値との差に前記受信シンボルの前回の受信時に前記第1の平滑化部が算出した移動平均値と前記受信シンボルの前回の受信時に前記第2の平滑化部が算出した移動平均値との差を乗じた値を、前記受信シンボルの前回の受信時に前記第1の平滑化部が算出した移動平均値と前記受信シンボルの前回の受信時に前記第2の平滑化部が算出した移動平均値との差の絶対値の二乗で除算する
ことを特徴とする請求項1または請求項2に記載の伝送路特性推定装置。
The division unit is configured to calculate a difference between a transmission path characteristic estimated by the current transmission path characteristic estimation unit and a moving average value calculated by the first smoothing unit at the previous reception of the reception symbol, from a previous time of the received symbol. The value obtained by multiplying the difference between the moving average value calculated by the first smoothing unit at the time of reception and the moving average value calculated by the second smoothing unit at the previous reception of the received symbol is the value of the received symbol. Divide by the square of the absolute value of the difference between the moving average value calculated by the first smoothing unit at the previous reception and the moving average value calculated by the second smoothing unit at the previous reception of the received symbol. The transmission path characteristic estimation apparatus according to claim 1 or 2, characterized by the above.
前記第1の平滑化部及び前記第2の平滑化部は、一次IIRフィルタであることを特徴とする請求項1から請求項3の何れか1項に記載の伝送路特性推定装置。   4. The transmission path characteristic estimation apparatus according to claim 1, wherein the first smoothing unit and the second smoothing unit are first-order IIR filters. 5. 連続する複数の受信シンボルから次の受信シンボルの伝送路特性を推定する伝送路特性推定装置を用いた伝送路特性推定方法であって、
前記受信シンボルの受信毎に、
現伝送路特性推定部は、現在の伝送路特性を推定し、
第1の平滑化部は、所定の平滑化係数を用いて前記現伝送特性推定部が推定した伝送路特性の移動平均値を算出し、
第2の平滑化部は、前記受信シンボルの受信毎に、前記第1の平滑化部の平滑化係数と異なる平滑化係数を用いて前記現伝送特性推定部が推定した伝送路特性の移動平均値を算出し、
除算部は、前記現伝送路特性推定部が推定した伝送路特性と前記受信シンボルの前回の受信時に前記第1の平滑化部が算出した移動平均値との差を、前記受信シンボルの前回の受信時に前記第1の平滑化部が算出した移動平均値と前記受信シンボルの前回の受信時に前記第2の平滑化部が算出した移動平均値との差で除算し、
次伝送路特性推定部は、前記第1の平滑化部が算出した移動平均値と前記第2の平滑化部が算出した移動平均値との差に、前記除算部の算出結果を乗算し、当該値に前記第1の平滑部が算出した移動平均値を加算することで、次の受信シンボルの伝送路特性を推定する
ことを特徴とする伝送路特性推定方法。
A transmission path characteristic estimation method using a transmission path characteristic estimation device that estimates a transmission path characteristic of a next reception symbol from a plurality of consecutive reception symbols,
Every time the received symbol is received,
The current channel characteristics estimation unit estimates the current channel characteristics,
The first smoothing unit calculates a moving average value of the transmission path characteristics estimated by the current transmission characteristic estimation unit using a predetermined smoothing coefficient,
The second smoothing unit is a moving average of the transmission path characteristics estimated by the current transmission characteristic estimation unit using a smoothing coefficient different from the smoothing coefficient of the first smoothing unit for each reception of the reception symbol. Calculate the value,
The division unit calculates a difference between the transmission path characteristic estimated by the current transmission path characteristic estimation unit and the moving average value calculated by the first smoothing unit at the previous reception of the received symbol, from the previous reception symbol. Dividing by the difference between the moving average value calculated by the first smoothing unit at the time of reception and the moving average value calculated by the second smoothing unit at the previous reception of the received symbol;
The next transmission path characteristic estimation unit multiplies the difference between the moving average value calculated by the first smoothing unit and the moving average value calculated by the second smoothing unit by the calculation result of the division unit, A transmission path characteristic estimation method, wherein the transmission path characteristic of the next received symbol is estimated by adding the moving average value calculated by the first smoothing unit to the value.
連続する複数の受信シンボルから次の受信シンボルの伝送路特性を推定する伝送路特性推定装置を、
前記受信シンボルの受信毎に現在の伝送路特性を推定する現伝送路特性推定部、
前記受信シンボルの受信毎に所定の平滑化係数を用いて前記現伝送特性推定部が推定した伝送路特性の移動平均値を算出する第1の平滑化部、
前記受信シンボルの受信毎に前記第1の平滑化部の平滑化係数と異なる平滑化係数を用いて前記現伝送特性推定部が推定した伝送路特性の移動平均値を算出する第2の平滑化部、
前記現伝送路特性推定部が推定した伝送路特性と前記受信シンボルの前回の受信時に前記第1の平滑化部が算出した移動平均値との差を、前記受信シンボルの前回の受信時に前記第1の平滑化部が算出した移動平均値と前記受信シンボルの前回の受信時に前記第2の平滑化部が算出した移動平均値との差で除算する除算部、
前記第1の平滑化部が算出した移動平均値と前記第2の平滑化部が算出した移動平均値との差に、前記除算部の算出結果を乗算し、当該値に前記第1の平滑部が算出した移動平均値を加算することで、次の受信シンボルの伝送路特性を推定する次伝送路特性推定部
として動作させるためのプログラム。
A transmission path characteristic estimation device that estimates a transmission path characteristic of a next reception symbol from a plurality of consecutive reception symbols,
A current channel characteristic estimator for estimating a current channel characteristic every time the received symbol is received;
A first smoothing unit that calculates a moving average value of transmission path characteristics estimated by the current transmission characteristic estimation unit using a predetermined smoothing coefficient for each reception of the received symbol;
Second smoothing for calculating a moving average value of the transmission path characteristic estimated by the current transmission characteristic estimation unit using a smoothing coefficient different from the smoothing coefficient of the first smoothing unit for each reception of the received symbol Part,
The difference between the transmission path characteristic estimated by the current transmission path characteristic estimation unit and the moving average value calculated by the first smoothing unit at the previous reception of the received symbol is expressed as the first time at the previous reception of the reception symbol. A division unit that divides by the difference between the moving average value calculated by one smoothing unit and the moving average value calculated by the second smoothing unit at the previous reception of the received symbol;
The difference between the moving average value calculated by the first smoothing unit and the moving average value calculated by the second smoothing unit is multiplied by the calculation result of the dividing unit, and the value is multiplied by the first smoothing value. A program for operating as a next channel characteristic estimation unit that estimates the channel characteristic of the next received symbol by adding the moving average values calculated by the unit.
JP2011084332A 2011-04-06 2011-04-06 Transmission path characteristic estimation device, transmission path characteristic estimation method, and program Active JP5610538B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2011084332A JP5610538B2 (en) 2011-04-06 2011-04-06 Transmission path characteristic estimation device, transmission path characteristic estimation method, and program

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2011084332A JP5610538B2 (en) 2011-04-06 2011-04-06 Transmission path characteristic estimation device, transmission path characteristic estimation method, and program

Publications (2)

Publication Number Publication Date
JP2012222487A true JP2012222487A (en) 2012-11-12
JP5610538B2 JP5610538B2 (en) 2014-10-22

Family

ID=47273579

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2011084332A Active JP5610538B2 (en) 2011-04-06 2011-04-06 Transmission path characteristic estimation device, transmission path characteristic estimation method, and program

Country Status (1)

Country Link
JP (1) JP5610538B2 (en)

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2000082978A (en) * 1998-09-04 2000-03-21 Fujitsu Ltd Demodulator
JP2001237796A (en) * 2000-02-25 2001-08-31 Nippon Telegr & Teleph Corp <Ntt> Multi-carrier signal detector
JP2007318479A (en) * 2006-05-26 2007-12-06 Mitsubishi Electric Corp Ofdm receiver
JP2009027717A (en) * 2007-07-20 2009-02-05 Fujitsu Ltd Carrier interference processing apparatus and method and receiver employing the same
JP2009516484A (en) * 2005-11-15 2009-04-16 クゥアルコム・インコーポレイテッド Method and apparatus for reducing noisy errors by filtering noisy estimates

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2000082978A (en) * 1998-09-04 2000-03-21 Fujitsu Ltd Demodulator
JP2001237796A (en) * 2000-02-25 2001-08-31 Nippon Telegr & Teleph Corp <Ntt> Multi-carrier signal detector
JP2009516484A (en) * 2005-11-15 2009-04-16 クゥアルコム・インコーポレイテッド Method and apparatus for reducing noisy errors by filtering noisy estimates
JP2007318479A (en) * 2006-05-26 2007-12-06 Mitsubishi Electric Corp Ofdm receiver
JP2009027717A (en) * 2007-07-20 2009-02-05 Fujitsu Ltd Carrier interference processing apparatus and method and receiver employing the same

Also Published As

Publication number Publication date
JP5610538B2 (en) 2014-10-22

Similar Documents

Publication Publication Date Title
JP4600559B2 (en) Receiving device, receiving method, and program
JP5278678B2 (en) Receiving method and receiving apparatus
KR101603342B1 (en) Receiving apparatus, receiving method, and recording medium
EP1551120A1 (en) Receiving device, receiving method, and device for measuring transmission channel characteristic
EP2213062A1 (en) Recursive time-frequency channel estimation
JP2010510742A (en) Channel estimation and equalization in OFDM receivers
EP2166721B1 (en) Receiving device, method and program
JP4944106B2 (en) Receiver and propagation path estimation method
JP4545209B2 (en) Orthogonal frequency division multiplexed signal receiving apparatus and receiving method thereof
JP4091804B2 (en) OFDM receiving apparatus and method implemented using polar coordinate system
WO2021047569A1 (en) Method and decoder for suppressing phase noise in an orthogonal frequency division multiplexing signal
JP4524704B2 (en) Demodulation circuit, demodulation method, program, and receiver
JP3870116B2 (en) Equalizer
JP5023007B2 (en) OFDM signal receiving apparatus and relay apparatus
JP5610538B2 (en) Transmission path characteristic estimation device, transmission path characteristic estimation method, and program
JP6028572B2 (en) Receiver
JP2007104574A (en) Multicarrier wireless receiver and receiving method
JP5812827B2 (en) Receiver
JP5099148B2 (en) Receiving device, receiving method, and program
JP5023006B2 (en) OFDM signal receiving apparatus and relay apparatus
JP2019198061A (en) Phase noise compensation device and method, receiver
JP5518261B2 (en) Equalizer and equalization method
JP5716617B2 (en) Signal processing circuit, signal processing method, and reception system
JP5542483B2 (en) OFDM receiver
JP2012019409A (en) Radio communication device

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20131004

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20140729

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

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20140828

R150 Certificate of patent or registration of utility model

Ref document number: 5610538

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150