JPWO2007020943A1 - OFDM communication method - Google Patents

OFDM communication method Download PDF

Info

Publication number
JPWO2007020943A1
JPWO2007020943A1 JP2007531011A JP2007531011A JPWO2007020943A1 JP WO2007020943 A1 JPWO2007020943 A1 JP WO2007020943A1 JP 2007531011 A JP2007531011 A JP 2007531011A JP 2007531011 A JP2007531011 A JP 2007531011A JP WO2007020943 A1 JPWO2007020943 A1 JP WO2007020943A1
Authority
JP
Japan
Prior art keywords
pilot sequence
ofdm
communication method
channel
ofdm communication
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.)
Ceased
Application number
JP2007531011A
Other languages
Japanese (ja)
Inventor
▲強▼ ▲呉▼
▲強▼ ▲呉▼
継峰 李
継峰 李
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.)
Panasonic Corp
Panasonic Holdings Corp
Original Assignee
Panasonic Corp
Matsushita Electric Industrial Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Panasonic Corp, Matsushita Electric Industrial Co Ltd filed Critical Panasonic Corp
Publication of JPWO2007020943A1 publication Critical patent/JPWO2007020943A1/en
Ceased legal-status Critical Current

Links

Images

Classifications

    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L27/00Modulated-carrier systems
    • H04L27/26Systems using multi-frequency codes
    • H04L27/2601Multicarrier modulation systems
    • H04L27/2614Peak power aspects
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L27/00Modulated-carrier systems
    • H04L27/26Systems using multi-frequency codes
    • H04L27/2601Multicarrier modulation systems
    • H04L27/2602Signal structure
    • H04L27/261Details of reference signals

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Mobile Radio Communication Systems (AREA)
  • Cable Transmission Systems, Equalization Of Radio And Reduction Of Echo (AREA)

Abstract

OFDM(Orthogonal Frequency Division Multiplexing)信号のPAPR(Peak to Average Power Ratio)を低減しつつ、伝送効率を向上させることができるとともに、受信データの誤り率を低下させることができるOFDM通信方法。この方法において、複数の重み係数のうちから、OFDM信号のPAPRを最小にする最適な重み係数を選択し、前記最適な重み係数に対応するシフト量にてパイロット系列を循環シフトし、循環シフトされた前記パイロット系列を送信データに挿入する。An OFDM communication method capable of improving transmission efficiency and reducing an error rate of received data while reducing PAPR (Peak to Average Power Ratio) of an OFDM (Orthogonal Frequency Division Multiplexing) signal. In this method, an optimum weighting factor that minimizes the PAPR of the OFDM signal is selected from a plurality of weighting factors, and the pilot sequence is cyclically shifted by a shift amount corresponding to the optimum weighting factor, and cyclically shifted. The pilot sequence is inserted into transmission data.

Description

本発明は、OFDM通信方法に関する。   The present invention relates to an OFDM communication method.

OFDM(Orthogonal Frequency Division Multiplexing)は高効率の伝送方式であり、OFDMでは、帯域幅のあるチャネルを複数の直交サブチャネルに分割し、サブチャネルごとに1つのサブキャリアを使用して変調を行い、それぞれのサブキャリアを並列伝送する。このため、全体的なチャネルがフラットではなく周波数選択性があるとしても、各サブチャネルは相対的にフラットであり、また、各サブチャネルの帯域幅はチャネル全体の帯域幅より狭いため、信号波形間の干渉を抑えることができる。   OFDM (Orthogonal Frequency Division Multiplexing) is a highly efficient transmission method. In OFDM, a channel with bandwidth is divided into a plurality of orthogonal subchannels, and modulation is performed using one subcarrier for each subchannel. Each subcarrier is transmitted in parallel. For this reason, even if the overall channel is not flat but frequency selective, each subchannel is relatively flat and the bandwidth of each subchannel is narrower than the overall channel bandwidth, so the signal waveform Interference can be suppressed.

OFDMが一般的なマルチキャリア伝送方式と異なる点は、複数の互いに直交するサブキャリアを周波数軸上において多重することができる点である。各サブキャリアが互いに直交していれば、多重されたサブキャリアから信号を分離することができる。このように、OFDMでは複数の互いに直交するサブキャリアを周波数軸上において多重することができるため、OFDMは、周波数利用効率が非常に良く、伝送効率の高い変調方式として注目されている。   The difference between OFDM and a general multicarrier transmission system is that a plurality of mutually orthogonal subcarriers can be multiplexed on the frequency axis. If the subcarriers are orthogonal to each other, the signal can be separated from the multiplexed subcarriers. In this way, since OFDM can multiplex a plurality of subcarriers orthogonal to each other on the frequency axis, OFDM is attracting attention as a modulation scheme with very good frequency utilization efficiency and high transmission efficiency.

また、OFDMは、マルチパスの影響に対する耐性が高く、シンボル間干渉を防ぐことができ、また、周波数選択性フェージングに強く、チャネルの利用効率が高いという特徴がある。よって、OFDMは、マルチパス伝送及びドップラー周波数シフトのある無線移動チャネルにおける高速データの伝送に好適である。OFDMは、既にヨーロッパDBA、DVB、HIPERLAN及びIEEE802.11への適用が成功している。   In addition, OFDM is characterized by high tolerance to multipath effects, prevention of inter-symbol interference, resistance to frequency selective fading, and high channel utilization efficiency. Therefore, OFDM is suitable for high-speed data transmission in a wireless mobile channel with multipath transmission and Doppler frequency shift. OFDM has already been successfully applied to European DBA, DVB, HIPERLAN, and IEEE 802.11.

複数のサブキャリアからなるOFDM信号では、各サブキャリアが同位相になったときに非常に大きなピーク電力が生じる。このため、OFDM通信では、線形性が高い増幅器を利用しなければ非線形歪みが生じてしまい通信品質が劣化する。しかし、携帯電話機のような小型の移動端末に線形性が高い増幅器を利用することは、部品コストや消費電力の点から好ましくない。よって、送信するOFDM信号自体にピーク電力の発生を回避するような工夫を施す必要がある。このように、OFDMにおいては、変調後の信号(OFDM信号)のPAPR(Peak to Average Power Ratio)を低減する必要がある。   In an OFDM signal composed of a plurality of subcarriers, a very large peak power is generated when the subcarriers have the same phase. For this reason, in OFDM communication, if an amplifier with high linearity is not used, nonlinear distortion occurs and communication quality deteriorates. However, it is not preferable to use an amplifier with high linearity for a small mobile terminal such as a cellular phone in terms of component cost and power consumption. Therefore, it is necessary to devise measures to avoid the generation of peak power in the OFDM signal to be transmitted. Thus, in OFDM, it is necessary to reduce the PAPR (Peak to Average Power Ratio) of the modulated signal (OFDM signal).

これに対し、例えば、Xiaodong Li, Leonard J. Cimini Jr., "Effects of clipping and filtering on the performance of OFDM", IEEE Communications Letters, vol.2, no.5, May 1998 pp.131-133では、帯域幅を制限することでPAPRを低減する方法が提案されている。しかし、帯域幅を制限すると帯域内の干渉及び帯域外の雑音が生じるため、通信品質が低下してしまう。   In contrast, for example, in Xiaodong Li, Leonard J. Cimini Jr., "Effects of clipping and filtering on the performance of OFDM", IEEE Communications Letters, vol.2, no.5, May 1998 pp.131-133, A method for reducing PAPR by limiting the bandwidth has been proposed. However, when the bandwidth is limited, in-band interference and out-of-band noise occur, and communication quality is degraded.

また、Stefan H. MAullerand Johannes B. Huber “A COMPARISON OF PEAK POWER REDUCTION SCHEMES FOR OFDM”, In Proc. of the IEEE Global Telecommunications Conference GLOBECOM '97, pp. 1-5, November 1997, Phoenix, Arizona, USAでは、選択マッピング(SLM)法を用いてPAPRを低減する方法が提案されている。SLM法では、長さがNである任意の(ランダムな)M個の位相シーケンスベクトルP(u=1,…,M)と入力された信号Xとの内積を計算し、この計算により得られたシーケンスに対してIFFT(Inverse Fast Fourier Transform)を行ってM個の時間領域信号を得る。そして、M個の時間領域信号において、PAPRが最も小さい信号を選択して送信に使用する。同時に、使用したランダムな位相シーケンスPを補助情報(side information)として復調のために受信側に送信する。Stefan H. MAullerand Johannes B. Huber “A COMPARISON OF PEAK POWER REDUCTION SCHEMES FOR OFDM”, In Proc. Of the IEEE Global Telecommunications Conference GLOBECOM '97, pp. 1-5, November 1997, Phoenix, Arizona, USA A method of reducing PAPR using a selective mapping (SLM) method has been proposed. In the SLM method, an inner product of an arbitrary (random) M phase sequence vector P u (u = 1,..., M) having a length of N and an input signal X is calculated, and obtained by this calculation. IFFT (Inverse Fast Fourier Transform) is performed on the obtained sequence to obtain M time domain signals. Then, among the M time domain signals, a signal having the smallest PAPR is selected and used for transmission. At the same time, the used random phase sequence P m is transmitted as side information to the receiving side for demodulation.

さらに、Heung-Gyoon Ryu, “A New PAPR Reduction Scheme : SPW(subblock phase weighting)”IEEE Transaction on Consumer Electronics,Vol.48, No.1, pp 81-89, Feb. 2002 では、シーケンスを部分的に伝送する方法(PTS法)が提案されている。この方法の基本的な原理はSLM法と同様であるが、転換ベクトルの構成がSLM法と異なる。PTS法では、まず入力されたデータベクトルをK個の互いに重複しないサブベクトルX,…,Xに分割し、各サブベクトルそれぞれの有効値の数はN/Kである。次いで、すべてのサブベクトルX(j=1,…,K)のサブキャリアにそれぞれ同様の位相要素Pをかける。得られたM通りの線形組み合わせの中から、PAPRが最小となるサブキャリアを選択して送信に使用する。同時に、選択されたサブキャリアに対応する最適な(P,…,P)を補助情報として受信側に送信する。Heung-Gyoon Ryu, “A New PAPR Reduction Scheme: SPW (subblock phase weighting)” IEEE Transaction on Consumer Electronics, Vol. 48, No. 1, pp 81-89, Feb. 2002 A transmission method (PTS method) has been proposed. The basic principle of this method is the same as that of the SLM method, but the configuration of the conversion vector is different from that of the SLM method. In the PTS method, first, an input data vector is divided into K non-overlapping subvectors X 1 ,..., Xk, and the number of valid values for each subvector is N / K. Next, the same phase element P j is applied to the subcarriers of all the subvectors X j (j = 1,..., K). From the obtained M linear combinations, a subcarrier having the smallest PAPR is selected and used for transmission. At the same time, the optimum (P 1 ,..., P k ) corresponding to the selected subcarrier is transmitted to the receiving side as auxiliary information.

上述のように、従来はOFDM信号の受信側に補助情報を送信する必要があるため、その分送信リソースが減り、伝送効率が低下するという問題がある。また送信された補助情報に誤りが生じると、受信データの誤り率が高まるという問題がある。   As described above, conventionally, since it is necessary to transmit auxiliary information to the receiving side of the OFDM signal, there is a problem in that transmission resources are reduced correspondingly and transmission efficiency is lowered. Further, when an error occurs in the transmitted auxiliary information, there is a problem that an error rate of received data increases.

本発明の目的は、OFDM信号のPAPRを低減しつつ、伝送効率を向上させることができるとともに、受信データの誤り率を低下させることができるOFDM通信方法を提供することである。   An object of the present invention is to provide an OFDM communication method capable of improving transmission efficiency and reducing an error rate of received data while reducing PAPR of an OFDM signal.

本発明のOFDM通信方法は、複数の重み係数のうちから、OFDM信号のPAPRを最小にする最適な重み係数を選択する選択ステップと、前記最適な重み係数に対応するシフト量にてパイロット系列を循環シフトするシフトステップと、循環シフトされた前記パイロット系列を送信データに挿入する挿入ステップと、を具備するようにした。   The OFDM communication method of the present invention includes a selection step of selecting an optimum weighting factor that minimizes the PAPR of an OFDM signal from a plurality of weighting factors, and a pilot sequence with a shift amount corresponding to the optimum weighting factor. A shift step for cyclic shift; and an insertion step for inserting the cyclically shifted pilot sequence into transmission data.

本発明によれば、OFDM信号のPAPRを低減させる場合でも、補助情報を送信しないため、伝送効率を向上させることができるとともに、補助情報の誤りによる復号誤りを防止して受信データの誤り率を低下させることができる。   According to the present invention, even when the PAPR of the OFDM signal is reduced, auxiliary information is not transmitted, so that transmission efficiency can be improved and decoding error due to an error in auxiliary information can be prevented to reduce the error rate of received data. Can be reduced.

CPが挿入されたOFDMシンボルを示す図The figure which shows the OFDM symbol in which CP was inserted OFDMシンボルの構成を示す図Diagram showing the structure of an OFDM symbol 本発明の実施の形態1に係るOFDM通信方法の送信側における動作フロー図Operation flow diagram on transmission side of OFDM communication method according to embodiment 1 of the present invention 本発明の実施の形態1に係るOFDM通信方法の受信側における動作フロー図Operation flow diagram on receiving side of OFDM communication method according to embodiment 1 of the present invention 本発明の実施の形態1に係る第一パスの位置により重み係数を確定するステップの概念図The conceptual diagram of the step which determines a weighting coefficient by the position of the 1st path | pass which concerns on Embodiment 1 of this invention. 本発明の実施の形態1に係るチャネル推定を示す模式図Schematic diagram showing channel estimation according to Embodiment 1 of the present invention. 本発明の実施の形態2に係るOFDM通信方法の送信側における動作フロー図Operation flow diagram on transmission side of OFDM communication method according to embodiment 2 of the present invention 本発明の実施の形態2に係るOFDM通信方法の受信側における動作フロー図Operation flow diagram on receiving side of OFDM communication method according to embodiment 2 of the present invention

以下、本発明の実施の形態について図面を参照して詳細に説明する。   Hereinafter, embodiments of the present invention will be described in detail with reference to the drawings.

<モデル定義及びチャネル推定>
次式(1)に示す行列A(n×n)を循環行列と定義する。すなわち、行列Aの各列は1列目を繰り返しシフトすることによって得られるものである。

Figure 2007020943
<Model definition and channel estimation>
A matrix A (n × n) shown in the following equation (1) is defined as a circulant matrix. That is, each column of the matrix A is obtained by repeatedly shifting the first column.
Figure 2007020943

よって、行列Aを式(2)のように書き換えることができる。

Figure 2007020943
Therefore, the matrix A can be rewritten as in equation (2).
Figure 2007020943

式(2)では、行列Aの1列目を(a,a,…,a)と示し、式(1)に基づいて、循環行列Aを構成する。In the equation (2), the first column of the matrix A is denoted as (a 1 , a 2 ,..., A n ), and the circulant matrix A is configured based on the equation (1).

パイロットXに対して、s=IFFT(X)とすると、sをパイロットXの時間領域での表記または時間領域のパイロットと呼ぶ。rを時間領域において受信した信号(CP(Cyclic Prefix)除去後)とすると、SISO−OFDMの時間領域モデルを式(3)で示すことができる。

Figure 2007020943
When s = IFFT (X) for pilot X, s is referred to as pilot X time domain notation or time domain pilot. If r is a signal received in the time domain (after CP (Cyclic Prefix) removal), the SISO-OFDM time domain model can be expressed by Equation (3).
Figure 2007020943

式(3)におけるパイロット系列Sは式(4)に示され、hは式(5)に示され、wは時間領域雑音である。

Figure 2007020943
Figure 2007020943
The pilot sequence S in Equation (3) is shown in Equation (4), h is shown in Equation (5), and w is time domain noise.
Figure 2007020943
Figure 2007020943

そして、式(3)における受信信号rは、sとチャネル行列hとを循環畳み込みした後に雑音wを足したものである。   The received signal r in equation (3) is obtained by adding noise w after circular convolution of s and channel matrix h.

図1はCPが挿入されたOFDMシンボルを示す図である。CP長はチャネルの最大マルチパス遅延長より長いものとする。よって、受信側では、ガード間隔を除去して得られる信号はy(n)であり、そしてFFT(Fast Fourier Transform)処理により、時間領域の信号y(n)を周波数領域の信号Y(k)に変換する。信号Y(k)を次式(6)に示す。ただしFFTのポイント数はNである。

Figure 2007020943
FIG. 1 is a diagram illustrating an OFDM symbol with a CP inserted therein. The CP length is longer than the maximum multipath delay length of the channel. Therefore, on the receiving side, the signal obtained by removing the guard interval is y (n), and the time domain signal y (n) is converted into the frequency domain signal Y (k) by FFT (Fast Fourier Transform) processing. Convert to The signal Y (k) is shown in the following equation (6). However, the number of points of FFT is N.
Figure 2007020943

よって、OFDMシステムの周波数領域における受信モデルは式(7)に示される。

Figure 2007020943
Therefore, the reception model in the frequency domain of the OFDM system is expressed by Equation (7).
Figure 2007020943

式(7)では、XはN×Nの対角行列で、その対角線要素は周波数領域のパイロットシンボルであり、Hは周波数領域でのチャネル応答でN×1のベクトルであり、Wは周波数領域での雑音N×1のベクトルである。   In Equation (7), X is an N × N diagonal matrix, the diagonal elements are frequency domain pilot symbols, H is the channel response in the frequency domain, and is an N × 1 vector, and W is the frequency domain. Is a vector of noise N × 1.

例えば、R. Negi and J. Coiffi, "Pilot tone selection for channel estimation in a mobile OFDM system", IEEE Trans. Consumer Electron., vol.44, pp. 1122-1128, Aug. 1998では、チャネル行列のLS(Least Square)推定(最小平均平方推定)が式(8)と(9)とにより示されている。

Figure 2007020943
Figure 2007020943
For example, in R. Negi and J. Coiffi, “Pilot tone selection for channel estimation in a mobile OFDM system”, IEEE Trans. Consumer Electron., Vol.44, pp. 1122-1128, Aug. 1998, the LS of the channel matrix (Least Square) estimation (minimum mean square estimation) is shown by equations (8) and (9).
Figure 2007020943
Figure 2007020943

よって、チャネル推定を行う際に、まず受信した信号の周波数成分とパイロットの周波数成分を利用し、式(9)に基づいて、チャネルの周波数領域の応答HLSを得る。そして、HLSをIFFT変換し、チャネルの時間領域の応答hLSを得る。hLSを次式(10)で示す。

Figure 2007020943
Therefore, when performing channel estimation, first, the frequency component of the received signal and the frequency component of the pilot are used to obtain a response H LS in the frequency domain of the channel based on Equation (9). Then, the H LS is IFFT transformed to obtain a channel time-domain response h LS . h LS is represented by the following formula (10).
Figure 2007020943

次いで、N個のチャネル推定値の中から、推定値の大きいチャネルを選択し、他のチャネル推定値を「0」にして出力することによって、時間領域での新たなチャネル応答hLSを取得する。例えば、N=64とすると、hLS=[1.2,0.8,0.5,0.001,0.002,0.0011,…,0.0012]が得られる。この場合、選択処理によって得られたチャネル応答はhLS=[1.2,0.8,0.5,0,…,…]となり、すなわち、チャネル数は3個であると考える。Next, a channel with a large estimated value is selected from the N channel estimated values, and another channel estimated value is set to “0” and output, thereby obtaining a new channel response h LS in the time domain. . For example, if N = 64, h LS = [1.2,0.8,0.5,0.001,0.002,0.0011, ..., 0.0012] is obtained. In this case, the channel response obtained by the selection process is h LS = [1.2, 0.8, 0.5, 0,...], That is, the number of channels is considered to be three.

なお、N個のチャネル推定値を閾値と比較する方法により、閾値より大きいチャネル推定値の数をチャネルの個数にしてもよい。この閾値は、予め決められた固定値でも、リアルタイムに計算された数値でもよく、例えば、推定された雑音平均電力のルート値のc倍(cは定数)と設定する閾値である。   Note that the number of channel estimation values larger than the threshold value may be the number of channels by a method of comparing the N channel estimation values with the threshold value. This threshold value may be a predetermined fixed value or a numerical value calculated in real time. For example, the threshold value is set to c times the estimated noise average power route value (c is a constant).

最後に、FFTにより新しいHLS=FFT(hLS)を得る。上述のような処理によって、一部の雑音影響が除去された。Finally, a new H LS = FFT (h LS ) is obtained by FFT. A part of noise influence was removed by the above processing.

<FFTの循環シフト性>
X=fft(x)とすると、xはベクトルの時間領域成分であり(sを行ベクトルとする)、Xはベクトルの周波数領域成分であり、それぞれ次式(12)と(13)に示される。

Figure 2007020943
Figure 2007020943
<Circular shift of FFT>
When X = fft (x), x is a time domain component of the vector (s is a row vector), and X is a frequency domain component of the vector, which are represented by the following equations (12) and (13), respectively. .
Figure 2007020943
Figure 2007020943

xを右方向へnシンボル循環シフトすることをx{(k+n)}と定義し、式(11)と式(12)とに基づいて、x{(k+n)}に対してフーリエー変換を行うと、周波数領域成分がej2πmn/NX[m]となる。すなわち、時間領域において循環シフトすることが、周波数領域での結果として位相シフトとなる。to n symbols cyclic shift to the right direction x is defined as x {(k + n) N}, based on the equations (11) Equations (12), the Fourier chromatography transform to x {(k + n) N} When this is done, the frequency domain component becomes ej2πmn / N X [m]. That is, a cyclic shift in the time domain results in a phase shift as a result in the frequency domain.

式(3)では、雑音を無視すると、s=[s,s,…,s]が得られるため、式(3)を次式(13)に書き換えることができる。

Figure 2007020943
In Expression (3), if noise is ignored, s = [s 1 , s 2 ,..., S N ] is obtained, and therefore Expression (3) can be rewritten into the following Expression (13).
Figure 2007020943

式(13)では、rはhとsとの循環畳み込みを示す。そして、Y=fft(r)、H=fft(h)、X=fft(s)とすると、次式(14)に書き換えることができる。

Figure 2007020943
In Equation (13), r indicates a circular convolution of h and s. When Y = fft (r), H = fft (h), and X = fft (s), the following equation (14) can be rewritten.
Figure 2007020943

すなわち、時間領域の畳み込み計算が周波数領域の掛け算に相当する。よって、式(15)が得られる。

Figure 2007020943
That is, the convolution calculation in the time domain corresponds to multiplication in the frequency domain. Therefore, Formula (15) is obtained.
Figure 2007020943

式(15)より、チャネルの周波数領域値H(m)が得られる。   From the equation (15), the frequency domain value H (m) of the channel is obtained.

送信側では、パイロット系列Sを右方向へnシンボル循環シフトしたもの、つまりs{(k+n)}を送信する。一方、受信側ではsが送信されると想定する場合、次式(16)に示すようになる。

Figure 2007020943
On the transmission side, the pilot sequence S is cyclically shifted by n symbols to the right, that is, s {(k + n) N } is transmitted. On the other hand, when it is assumed that s is transmitted on the receiving side, the following equation (16) is obtained.
Figure 2007020943

すなわち、チャネルを推定する場合、sが送信されたと想定して推定を行うため、次式(17)になる。

Figure 2007020943
That is, when estimating the channel, since estimation is performed assuming that s is transmitted, the following equation (17) is obtained.
Figure 2007020943

式(17)を時間領域に変換すると、式(18)が得られる。

Figure 2007020943
When Expression (17) is converted into the time domain, Expression (18) is obtained.
Figure 2007020943

(実施の形態1)
図2はOFDMシンボルの構成を示す図である。OFDMシンボルでは、Nポイントをチャネル推定用のパイロットとし、その後のMポイントは有効データである。N+Mポイントの間では、チャネルが変化しない(M≧N)と考える。これは、スローフェージングにおいては、チャネルの変化が緩やかであるため、ある一定の時間内ではチャネルが変化しないと見なすことができるからである。
(Embodiment 1)
FIG. 2 is a diagram showing the structure of the OFDM symbol. In the OFDM symbol, N points are pilots for channel estimation, and the subsequent M points are valid data. The channel does not change between N + M points (M ≧ N). This is because in slow fading, the channel changes slowly, so that it can be considered that the channel does not change within a certain period of time.

よって、1つのOFDMデータシンボルを延ばす(Mポイント)ことができる。チャネルを推定する場合、チャネルの長さが一定範囲内のものに対して、小さめのシンボル(Nポイント)を用いて当該チャネルの時間領域応答を推定してもよい。   Therefore, one OFDM data symbol can be extended (M points). When estimating a channel, the time domain response of the channel may be estimated using a smaller symbol (N points) for a channel having a length within a certain range.

図3は本発明の実施の形態1に係るOFDM通信方法の送信側における動作フロー図である。実施の形態1では、基本パイロット系列Sの周波数領域成分をs,s,s,sと定義し、基本パイロット系列Sの時間領域成分をa,a,a,aと定義する。送信側と受信側ともに複数の重み係数を一定の順番で記憶し(S310)、例えば、M個の重み係数P=[p11,p21,…,pk1],P=[p12,p22,…,pk2],…,P=[p1M,p2M,…,pkM]を1〜Mの順番で記憶する。FIG. 3 is an operation flowchart on the transmission side of the OFDM communication method according to Embodiment 1 of the present invention. In the first embodiment, the frequency domain components of the basic pilot sequence S are defined as s 1 , s 2 , s 3 , s 4, and the time domain components of the basic pilot sequence S are a 1 , a 2 , a 3 , a 4. It is defined as A plurality of weighting factors are stored in a certain order on both the transmitting side and the receiving side (S310). For example, M weighting factors P 1 = [p 11 , p 21 ,..., P k1 ], P 2 = [p 12 , p 22 ,..., p k2 ],..., P M = [p 1M , p 2M ,..., p kM ] are stored in the order of 1 to M.

次いで、変調後の信号(OFDM信号)のPAPRを低減するために、OFDM信号のPAPRを最小にする最適な重み係数を選択する(S320)。具体的には、データベクトルをX=[x,x,…,xN−1]とし、このデータベクトルXをk組のサブベクトルに分割する。例えば、N=8,k=2である場合、それぞれ4つの値からなる2組に分割される。xからxを一組として第1サブベクトルXを形成し、X=[x,x,x,x,0,0,0,0]となる。そして、xからxを一組として第2ベクトルXを形成し、X=[0,0,0,0,x,x,x,x]となる。Next, in order to reduce the PAPR of the modulated signal (OFDM signal), an optimal weighting coefficient that minimizes the PAPR of the OFDM signal is selected (S320). Specifically, the data vector is X = [x 0 , x 1 ,..., X N−1 ], and the data vector X is divided into k subvectors. For example, when N = 8 and k = 2, it is divided into two sets each consisting of four values. The first subvector X 1 is formed by taking x 0 to x 3 as a set, and X 1 = [x 0 , x 1 , x 2 , x 3 , 0, 0, 0, 0]. Then, the second vector X 2 is formed with x 4 to x 7 as a set, and X 2 = [0, 0, 0, 0, x 4 , x 5 , x 6 , x 7 ].

そして、2つのサブベクトルを次式(19)に示すように組み合わせる。

Figure 2007020943
Then, the two subvectors are combined as shown in the following equation (19).
Figure 2007020943

ij(1≦i≦k,1≦j≦M)は重み係数であり、Pij=exp(jθij)、且つ、θij∈[0,2π]を満たす。P ij (1 ≦ i ≦ k, 1 ≦ j ≦ M) is a weighting coefficient, and satisfies P ij = exp (jθ ij ) and θ ij ε [0, 2π].

そして、X´ jに対してIDFTを行い、x´ j=IDFT(X´ j)が得られる。式(17)に基づいて、IDFTの線形特性を用いて、次式(20)を得る。

Figure 2007020943
Then, 'performs IDFT with respect to j, x' X j = IDFT (X 'j) is obtained. Based on the equation (17), the following equation (20) is obtained using the linear characteristics of the IDFT.
Figure 2007020943

式(20)により計算されたX´ jのピーク値が最小となるものを時間領域の送信信号とする。この場合、対応する重み係数Pijが最適な重み係数となる。ただし、最適な重み係数は次式(21)を満たす。

Figure 2007020943
A signal having a minimum peak value of X j calculated by Expression (20) is defined as a time domain transmission signal. In this case, the corresponding weight coefficient P ij is the optimum weight coefficient. However, the optimum weight coefficient satisfies the following equation (21).
Figure 2007020943

式(21)により、M個の重み係数から最適な重み係数を選択する。   The optimum weighting factor is selected from the M weighting factors according to the equation (21).

M個の重み係数からQ番目(Q≦M)の重み係数を選択した場合、パイロット系列を(Q−1)×Bだけ循環シフトする(S330)。ただしB≧1(整数)である。そして、シフトされたパイロット系列を時間領域において送信データに挿入する。   When the Qth (Q ≦ M) weighting factor is selected from the M weighting factors, the pilot sequence is cyclically shifted by (Q−1) × B (S330). However, B ≧ 1 (integer). Then, the shifted pilot sequence is inserted into transmission data in the time domain.

例えば、選択した重み係数の順位が2番目である場合、基本パイロット系列を左へ2つシフトして、a,a,a,aとし、これを時間領域において送信データに挿入する。すなわち、挿入したパイロット系列は式(1)に示された行列Aの3列目である。For example, when the rank of the selected weighting factor is the second, the basic pilot sequence is shifted two times to the left to become a 3 , a 4 , a 1 , a 2 , which are inserted into the transmission data in the time domain . That is, the inserted pilot sequence is the third column of the matrix A shown in Equation (1).

図4は、本発明の実施の形態1に係るOFDM通信方法の受信側における動作フロー図である。まず式(9)を用いてチャネル推定を行い(S410)、チャネルの周波数領域における応答HLSを推定する。そして、応答HLSに対してIFFTを行い、チャネルの時間領域における応答hLSを得る。次いで、閾値処理を行って雑音影響を除去する。そして、雑音が除去された時間領域における応答に対してFFTを行い、フィルタ(雑音除去)処理後の周波数領域における応答を得る。FIG. 4 is an operation flowchart on the receiving side of the OFDM communication method according to Embodiment 1 of the present invention. First, channel estimation is performed using equation (9) (S410), and the response H LS in the frequency domain of the channel is estimated. Then, it performs an IFFT on the response H LS, get response h LS in the time domain of the channel. Next, threshold processing is performed to remove noise effects. Then, FFT is performed on the response in the time domain from which noise is removed, and a response in the frequency domain after the filter (noise removal) processing is obtained.

次いで、チャネル応答があった第一パスの位置を確認する(S420)。具体的には、式(17)と(18)より、時間領域のパイロット系列が左へシフトされた分、推定されたチャネル値も同じく左へシフトされていることが分かる。例えば、2パスのチャネルでは、時間領域の基本パイロット系列a,a,a,aが挿入される場合、推定されるチャネル値はh,h,0,0となり、換言すると、第一パスが一番目に位置する。一方、時間領域の基本パイロット系列が左へ1つシフトされてパイロット系列a,a,a,aが挿入される場合、推定されるチャネル値はh,0,0,hとなり、つまり第一パスが二番目に位置する。また、時間領域の基本パイロット系列が左へ2つシフトされてパイロット系列a,a,a,aが挿入される場合、推定されるチャネル値は0,0,h,hとなり、つまり第一パスが三番目に位置する。Next, the position of the first path where there is a channel response is confirmed (S420). Specifically, from equations (17) and (18), it can be seen that the estimated channel value is also shifted to the left as the time-domain pilot sequence is shifted to the left. For example, in a two-path channel, when time domain basic pilot sequences a 1 , a 2 , a 3 , a 4 are inserted, the estimated channel values are h 1 , h 2 , 0, 0, in other words The first pass is located first. On the other hand, when the basic pilot sequence in the time domain is shifted one to the left and the pilot sequences a 2 , a 3 , a 4 , a 1 are inserted, the estimated channel values are h 2 , 0, 0, h 1. That is, the first pass is located second. In addition, when the basic pilot sequence in the time domain is shifted to the left by two and pilot sequences a 3 , a 4 , a 1 , a 2 are inserted, the estimated channel values are 0, 0, h 1 , h 2. That is, the first path is located third.

次いで、確定された第一パスの位置により、送信側で選択された重み係数を確定する(S430)。図5を参照しながら説明する。図5は第一パスの位置により重み係数を確定するステップの概念図である。推定されたチャネル値を右へ順番に、Z,…,Z,0,…,0(後尾がすべて「0」となる、図6を参照)が得られるまでシフトする。本実施の形態では、右へ2つシフトすることで、h,h,0,0が得られるため、受信側は、送信側で2番目の重み係数Pを選択してPAPRを低減したと確定することができる。例えばP=[p12,p22]=[1,−1]である。Next, the weighting factor selected on the transmission side is determined based on the determined position of the first path (S430). This will be described with reference to FIG. FIG. 5 is a conceptual diagram of the steps for determining the weighting coefficient based on the position of the first pass. The estimated channel values are shifted in order to the right until Z,..., Z, 0,..., 0 (the tail is all “0”, see FIG. 6) are obtained. In the present embodiment, h 1 , h 2 , 0, 0 are obtained by shifting two to the right, so that the receiving side selects the second weighting factor P 2 on the transmitting side to reduce PAPR. Can be confirmed. For example, P 2 = [p 12 , p 22 ] = [1, −1].

なお、送信側のパイロット系列及び受信側のチャネル値を1つずつシフトするとは限らず、シフトペースである定数Bを1より大きい整数としてもよい。B>1とすることにより、チャネル推定に誤りがある場合(例えば、雑音の影響により、推定されるマルチパスのチャネルが実際より多い、または、実際より少ない場合)でも、より正確にシフト量を推定することができ、より正確に重み係数を確定することができる。   Note that the pilot sequence on the transmission side and the channel value on the reception side are not necessarily shifted one by one, and the constant B, which is the shift pace, may be an integer greater than one. By setting B> 1, even when there is an error in channel estimation (for example, when there are more or fewer estimated multipath channels due to the influence of noise), the shift amount can be more accurately determined. The weighting factor can be determined more accurately.

その後、受信信号からCPを除去してy=[y,y,…,y]が得られ、yをIFFTすると次式(22)が得られる。

Figure 2007020943
Thereafter, CP is removed from the received signal to obtain y = [y 1 , y 2 ,..., Y N ], and y is IFFT to obtain the following equation (22).
Figure 2007020943

周波数領域のチャネル値をH=[h,h,…,h]とすると、Z=Y/h/p12(i=0,1,2,3)、Z=Y/h/p22(i=4,5,6,7)となる。すなわち、周波数領域の受信信号のチャネル影響を除去して(Y/h)、周波数領域において実際に送信された信号(雑音を含む)が得られる。ここで得られた信号は、有効送信信号(有効データ)に重み係数がかけられた信号(線形組み合わせ信号)であるため、確定された重み係数Pで割る処理を行う。最後に、Zを復調してデータを復元する。If the channel values in the frequency domain are H = [h 1 , h 2 ,..., H N ], Z i = Y i / h i / p 12 (i = 0,1,2,3), Z i = Y i / h i / p 22 (i = 4, 5, 6, 7). That is, the channel effect of the received signal in the frequency domain is removed (Y i / h i ), and a signal (including noise) actually transmitted in the frequency domain is obtained. Signals obtained herein are the effective transmission signal (valid data) weighting factor is applied to the signal (linear combination signal), it performs a process of dividing the determined weighting factor P 2. Finally, Z 1 is demodulated to restore the data.

このように、本実施の形態によれば、基本パイロット系列を定義し、予め複数の重み係数を所定の順番で、それぞれ送信側と受信側とに記憶させる。送信側では、PAPRを低減する過程において使用された最適な重み係数の順位を確定し、その順位に応じて基本パイロット系列をシフトし、シフトされたパイロット系列を時間領域において送信データに挿入する。受信側では、FFTの特性を利用して、受信信号のパイロット系列を用いて推定されたチャネル値の第一パスの位置によって、送信側で使用された最適な重み係数を確定する。このため、本実施の形態によれば、送信側で使用された重み係数を受信側で確定するための補助情報を別途送信する必要がないため、伝送効率を向上することができるとともに、補助情報の誤りによる複合誤りを防止して受信データの誤り率を低下させることができる。   Thus, according to the present embodiment, a basic pilot sequence is defined, and a plurality of weighting factors are stored in advance in a predetermined order on the transmitting side and the receiving side, respectively. On the transmission side, the rank of the optimum weighting factor used in the process of reducing the PAPR is determined, the basic pilot sequence is shifted according to the rank, and the shifted pilot sequence is inserted into the transmission data in the time domain. On the receiving side, the optimum weighting factor used on the transmitting side is determined based on the position of the first path of the channel value estimated using the pilot sequence of the received signal, using the characteristics of the FFT. Therefore, according to the present embodiment, it is not necessary to separately transmit auxiliary information for determining the weighting factor used on the transmission side on the reception side, so that transmission efficiency can be improved and auxiliary information It is possible to reduce the error rate of the received data by preventing the compound error due to the error.

(実施の形態2)
図7は本発明の実施の形態2に係るOFDM通信方法の送信側における動作フロー図である。本実施の形態では、周波数領域のパイロット系列を分割して重み係数にするため、予め送信側と受信側とに複数の重み係数を記憶させる必要がない。
(Embodiment 2)
FIG. 7 is an operation flowchart on the transmission side of the OFDM communication method according to Embodiment 2 of the present invention. In the present embodiment, since the pilot sequence in the frequency domain is divided into weighting factors, it is not necessary to store a plurality of weighting factors in advance on the transmitting side and the receiving side.

上述のように、チャネルを推定する過程において、チャネルのLS推定性能を最良にするため、式(3)のS行列は直交行列である必要がある。すなわち時間領域のパイロット系列がシフト直交である必要がある。以上の条件を満たすブロック状パイロットに対して、周波数領域のパイロット系列が等電力であることが要求される。周波数領域のパイロットシンボルの電力を1とすると、周波数領域のパイロットがe(θは0〜2πの任意値)であるため、上述の条件を満たす。よって、周波数領域のパイロット系列を分割して重み係数にすることができる。As described above, in order to optimize the LS estimation performance of the channel in the process of estimating the channel, the S matrix of Equation (3) needs to be an orthogonal matrix. That is, the time-domain pilot sequence needs to be shift orthogonal. For block pilots that satisfy the above conditions, the pilot sequences in the frequency domain are required to be of equal power. Assuming that the power of the pilot symbol in the frequency domain is 1, the frequency domain pilot is e (θ is an arbitrary value from 0 to 2π), so the above condition is satisfied. Therefore, the frequency domain pilot sequence can be divided into weighting factors.

例えば、N=4とし、周波数領域のパイロット系列をK=2個に分割し、基本パイロット系列をs,s,s,sと定義する。基本パイロット系列の時間領域成分はa,a,a,aである。s,sを1つ目の重み係数とし、s,sを2つ目の重み係数とする(S710)。For example, N = 4, the frequency domain pilot sequence is divided into K = 2, and the basic pilot sequences are defined as s 1 , s 2 , s 3 , s 4 . The time domain components of the basic pilot sequence are a 1 , a 2 , a 3 , a 4 . s 1 and s 2 are set as the first weighting factors, and s 3 and s 4 are set as the second weighting factors (S710).

ついで、実施の形態1と同様方法で、信号のPAPRを最小にする最適な重み係数を選択する(S720)。   Next, an optimum weighting coefficient that minimizes the PAPR of the signal is selected by the same method as in the first embodiment (S720).

選択された最適重み係数の順位に応じて基本パイロット系列をシフトする(S730)。選択された重み係数が2番目の重み係数s,sである場合、基本パイロット系列を左へ2つシフトする。すなわち、挿入されるパイロット系列が式(1)に示された行列Aの3列目a,a,a,aとなる。The basic pilot sequence is shifted according to the order of the selected optimum weight coefficient (S730). If the selected weighting factor is the second weighting factor s 3 , s 4 , the basic pilot sequence is shifted two times to the left. That is, the inserted pilot sequence is the third column a 3 , a 4 , a 1 , a 2 of the matrix A shown in Equation (1).

図8は本発明の実施の形態2に係るOFDM通信方法の受信側における動作フロー図である。図8に示すように、式(9)によりチャネル推定を行い(S810)、式(9)の周波数領域のパイロット系列は基本パイロット系列である。すなわち、式(9)では、X(1)=s、X(2)=s、X(3)=s、X(4)=sである。FIG. 8 is an operation flowchart on the receiving side of the OFDM communication method according to Embodiment 2 of the present invention. As shown in FIG. 8, channel estimation is performed using equation (9) (S810), and the pilot sequence in the frequency domain of equation (9) is a basic pilot sequence. That is, in the formula (9), X (1) = s 1 , X (2) = s 2 , X (3) = s 3 , X (4) = s 4 .

次いで、チャネル応答があった第一パスの位置を確認する(S820)。確定された第一パスの位置により、送信側で選択されたパイロットを確定する(S830)。実施の形態1と同様方法により、チャネル推定値を右へ2つシフトすると、h,h,0,0が得られるため、送信側で選択されたパイロット系列が2番目であると確定することができる。つまり、データ復元に必要な重み係数はs,sであると確定することができる。Next, the position of the first path where there is a channel response is confirmed (S820). Based on the determined position of the first path, the pilot selected on the transmission side is determined (S830). When the channel estimation value is shifted two times to the right by the same method as in the first embodiment, h 1 , h 2 , 0, 0 are obtained, so that the pilot sequence selected on the transmission side is determined to be the second. be able to. That is, it can be determined that the weighting factors necessary for data restoration are s 3 and s 4 .

このように、本実施の形態によれば、送信側では、周波数領域のパイロット系列を分割して重み係数とし、PAPRを低減する過程において選択する重み係数が複数に分割されたパイロット系列での順位を確定し、その順位に応じて基本パイロット系列をシフトして、シフトされたパイロット系列を時間領域において送信データに挿入する。受信側では、FFTの特性を利用して、チャネル応答があった第一パスの位置を確定することによって、送信側で選択された最適な重み係数を確定する。よって、本実施の形態によれば、特定の重み係数を予め送信側と受信側とに記憶させる必要がないため、実施の形態1よりも簡易なOFDM通信方法を提供することができる。   As described above, according to the present embodiment, on the transmitting side, the pilot sequence in the frequency domain is divided into weighting factors, and the order in the pilot sequence in which the weighting factor selected in the process of reducing the PAPR is divided into a plurality of Is determined, the basic pilot sequence is shifted according to the order, and the shifted pilot sequence is inserted into the transmission data in the time domain. On the receiving side, the optimum weighting factor selected on the transmitting side is determined by determining the position of the first path where there was a channel response using the characteristics of the FFT. Therefore, according to the present embodiment, there is no need to store a specific weighting factor in advance on the transmitting side and the receiving side, so that it is possible to provide a simpler OFDM communication method than in the first embodiment.

以上の説明は、典型的な実施の形態について説明したが、本発明の主旨と範囲から逸脱しない限り、種々の変更、置き換え及び追加をすることができる。そのため、本発明は上記の実施の形態に拘らず、請求の範囲及びそれに均等する範囲によって限定されるものである。   Although the above description has been given of typical embodiments, various changes, substitutions and additions can be made without departing from the spirit and scope of the present invention. Therefore, the present invention is not limited to the above-described embodiments, but is limited by the scope of the claims and the scope equivalent thereto.

本明細書は、2005年8月16日出願の中国出願番号200510091685.6に基づくものである。その内容は、全てここに含めておく。   This specification is based on Chinese application No. 2005000915685.6 filed on Aug. 16, 2005. The contents are all included here.

本発明は、送信信号のピーク電力を抑える必要のあるOFDM通信等に好適である。   The present invention is suitable for OFDM communication or the like that needs to suppress the peak power of a transmission signal.

本発明は、OFDM通信方法に関する。   The present invention relates to an OFDM communication method.

OFDM(Orthogonal Frequency Division Multiplexing)は高効率の伝送方式であり、OFDMでは、帯域幅のあるチャネルを複数の直交サブチャネルに分割し、サブチャネルごとに1つのサブキャリアを使用して変調を行い、それぞれのサブキャリアを並列伝送する。このため、全体的なチャネルがフラットではなく周波数選択性があるとしても、各サブチャネルは相対的にフラットであり、また、各サブチャネルの帯域幅はチャネル全体の帯域幅より狭いため、信号波形間の干渉を抑えることができる。   OFDM (Orthogonal Frequency Division Multiplexing) is a highly efficient transmission method. In OFDM, a channel with bandwidth is divided into a plurality of orthogonal subchannels, and modulation is performed using one subcarrier for each subchannel. Each subcarrier is transmitted in parallel. For this reason, even if the overall channel is not flat but frequency selective, each subchannel is relatively flat, and the bandwidth of each subchannel is narrower than the overall channel bandwidth, so the signal waveform Interference can be suppressed.

OFDMが一般的なマルチキャリア伝送方式と異なる点は、複数の互いに直交するサブキャリアを周波数軸上において多重することができる点である。各サブキャリアが互いに直交していれば、多重されたサブキャリアから信号を分離することができる。このように、OFDMでは複数の互いに直交するサブキャリアを周波数軸上において多重することができるため、OFDMは、周波数利用効率が非常に良く、伝送効率の高い変調方式として注目されている。   The difference between OFDM and a general multicarrier transmission system is that a plurality of mutually orthogonal subcarriers can be multiplexed on the frequency axis. If the subcarriers are orthogonal to each other, the signal can be separated from the multiplexed subcarriers. In this way, since OFDM can multiplex a plurality of subcarriers orthogonal to each other on the frequency axis, OFDM is attracting attention as a modulation scheme with very good frequency utilization efficiency and high transmission efficiency.

また、OFDMは、マルチパスの影響に対する耐性が高く、シンボル間干渉を防ぐことができ、また、周波数選択性フェージングに強く、チャネルの利用効率が高いという特徴がある。よって、OFDMは、マルチパス伝送及びドップラー周波数シフトのある無線移動チャネルにおける高速データの伝送に好適である。OFDMは、既にヨーロッパDBA、DVB、HIPERLAN及びIEEE802.11への適用が成功している。   In addition, OFDM is characterized by high tolerance to multipath effects, prevention of inter-symbol interference, resistance to frequency selective fading, and high channel utilization efficiency. Therefore, OFDM is suitable for high-speed data transmission in a wireless mobile channel with multipath transmission and Doppler frequency shift. OFDM has already been successfully applied to European DBA, DVB, HIPERLAN, and IEEE 802.11.

複数のサブキャリアからなるOFDM信号では、各サブキャリアが同位相になったときに非常に大きなピーク電力が生じる。このため、OFDM通信では、線形性が高い増幅器を利用しなければ非線形歪みが生じてしまい通信品質が劣化する。しかし、携帯電話機のような小型の移動端末に線形性が高い増幅器を利用することは、部品コストや消費電力の点から好ましくない。よって、送信するOFDM信号自体にピーク電力の発生を回避するような工夫を施す必要がある。このように、OFDMにおいては、変調後の信号(OFDM信号)のPAPR(Peak to Average Power Ratio)を低減する必要がある。   In an OFDM signal composed of a plurality of subcarriers, a very large peak power is generated when the subcarriers have the same phase. For this reason, in OFDM communication, if an amplifier with high linearity is not used, nonlinear distortion occurs and communication quality deteriorates. However, it is not preferable to use an amplifier with high linearity for a small mobile terminal such as a cellular phone in terms of component cost and power consumption. Therefore, it is necessary to devise measures to avoid the generation of peak power in the OFDM signal to be transmitted. Thus, in OFDM, it is necessary to reduce the PAPR (Peak to Average Power Ratio) of the modulated signal (OFDM signal).

これに対し、例えば、Xiaodong Li, Leonard J. Cimini Jr., "Effects of clipping and filtering on the performance of OFDM", IEEE Communications Letters, vol.2, no.5, May 1998 pp.131-133では、帯域幅を制限することでPAPRを低減する方法が提案されている。しかし、帯域幅を制限すると帯域内の干渉及び帯域外の雑音が生じるため、通信品質が低下してしまう。   In contrast, for example, in Xiaodong Li, Leonard J. Cimini Jr., "Effects of clipping and filtering on the performance of OFDM", IEEE Communications Letters, vol.2, no.5, May 1998 pp.131-133, A method for reducing PAPR by limiting the bandwidth has been proposed. However, when the bandwidth is limited, in-band interference and out-of-band noise occur, and communication quality is degraded.

また、Stefan H. MAullerand Johannes B. Huber “A COMPARISON OF PEAK POWER REDUCTION SCHEMES FOR OFDM”, In Proc. of the IEEE Global Telecommunications Conference GLOBECOM '97, pp. 1-5, November 1997, Phoenix, Arizona, USAでは、選択マッピング(SLM)法を用いてPAPRを低減する方法が提案されている。SLM法では、長さがNである任意の(ランダムな)M個の位相シーケンスベクトルP(u=1,…,M)と入力された信号Xとの内積を計算し、この計算により得られたシーケンスに対してIFFT(Inverse Fast Fourier Transform)を行ってM個の時間領域信号を得る。そして、M個の時間領域信号において、PAPRが最も小さい信号を選択して送信に使用する。同時に、使用したランダムな位相シーケンスPを補助情報(side information)として
復調のために受信側に送信する。
Stefan H. MAullerand Johannes B. Huber “A COMPARISON OF PEAK POWER REDUCTION SCHEMES FOR OFDM”, In Proc. Of the IEEE Global Telecommunications Conference GLOBECOM '97, pp. 1-5, November 1997, Phoenix, Arizona, USA A method of reducing PAPR using a selective mapping (SLM) method has been proposed. In the SLM method, an inner product of an arbitrary (random) M phase sequence vector P u (u = 1,..., M) having a length of N and an input signal X is calculated, and obtained by this calculation. IFFT (Inverse Fast Fourier Transform) is performed on the obtained sequence to obtain M time domain signals. Then, among the M time domain signals, a signal having the smallest PAPR is selected and used for transmission. At the same time, the used random phase sequence P m is transmitted as side information to the receiving side for demodulation.

さらに、Heung-Gyoon Ryu, “A New PAPR Reduction Scheme : SPW(subblock phase weighting)”IEEE Transaction on Consumer Electronics,Vol.48, No.1, pp 81-89, Feb. 2002 では、シーケンスを部分的に伝送する方法(PTS法)が提案されている。この方法の基本的な原理はSLM法と同様であるが、転換ベクトルの構成がSLM法と異なる。PTS法では、まず入力されたデータベクトルをK個の互いに重複しないサブベクトルX,…,Xに分割し、各サブベクトルそれぞれの有効値の数はN/Kである。次いで、すべてのサブベクトルX(j=1,…,K)のサブキャリアにそれぞれ同様の位相要素Pをかける。得られたM通りの線形組み合わせの中から、PAPRが最小となるサブキャリアを選択して送信に使用する。同時に、選択されたサブキャリアに対応する最適な(P,…,P)を補助情報として受信側に送信する。 Heung-Gyoon Ryu, “A New PAPR Reduction Scheme: SPW (subblock phase weighting)” IEEE Transaction on Consumer Electronics, Vol. 48, No. 1, pp 81-89, Feb. 2002 A transmission method (PTS method) has been proposed. The basic principle of this method is the same as that of the SLM method, but the configuration of the conversion vector is different from that of the SLM method. In the PTS method, first, an input data vector is divided into K non-overlapping subvectors X 1 ,..., Xk, and the number of valid values for each subvector is N / K. Next, the same phase element P j is applied to the subcarriers of all the subvectors X j (j = 1,..., K). From the obtained M linear combinations, a subcarrier having the smallest PAPR is selected and used for transmission. At the same time, the optimum (P 1 ,..., P k ) corresponding to the selected subcarrier is transmitted to the receiving side as auxiliary information.

上述のように、従来はOFDM信号の受信側に補助情報を送信する必要があるため、その分送信リソースが減り、伝送効率が低下するという問題がある。また送信された補助情報に誤りが生じると、受信データの誤り率が高まるという問題がある。   As described above, conventionally, since it is necessary to transmit auxiliary information to the receiving side of the OFDM signal, there is a problem in that transmission resources are reduced correspondingly and transmission efficiency is lowered. Further, when an error occurs in the transmitted auxiliary information, there is a problem that an error rate of received data increases.

本発明の目的は、OFDM信号のPAPRを低減しつつ、伝送効率を向上させることができるとともに、受信データの誤り率を低下させることができるOFDM通信方法を提供することである。   An object of the present invention is to provide an OFDM communication method capable of improving transmission efficiency and reducing an error rate of received data while reducing PAPR of an OFDM signal.

本発明のOFDM通信方法は、複数の重み係数のうちから、OFDM信号のPAPRを最小にする最適な重み係数を選択する選択ステップと、前記最適な重み係数に対応するシフト量にてパイロット系列を循環シフトするシフトステップと、循環シフトされた前記パイロット系列を送信データに挿入する挿入ステップと、を具備するようにした。   The OFDM communication method of the present invention includes a selection step of selecting an optimum weighting factor that minimizes the PAPR of an OFDM signal from a plurality of weighting factors, and a pilot sequence with a shift amount corresponding to the optimum weighting factor. A shift step for cyclic shift; and an insertion step for inserting the cyclically shifted pilot sequence into transmission data.

本発明によれば、OFDM信号のPAPRを低減させる場合でも、補助情報を送信しないため、伝送効率を向上させることができるとともに、補助情報の誤りによる復号誤りを防止して受信データの誤り率を低下させることができる。   According to the present invention, even when the PAPR of the OFDM signal is reduced, auxiliary information is not transmitted, so that transmission efficiency can be improved and decoding error due to an error in auxiliary information can be prevented to reduce the error rate of received data. Can be reduced.

以下、本発明の実施の形態について図面を参照して詳細に説明する。   Hereinafter, embodiments of the present invention will be described in detail with reference to the drawings.

<モデル定義及びチャネル推定>
次式(1)に示す行列A(n×n)を循環行列と定義する。すなわち、行列Aの各列は
1列目を繰り返しシフトすることによって得られるものである。

Figure 2007020943
<Model definition and channel estimation>
A matrix A (n × n) shown in the following equation (1) is defined as a circulant matrix. That is, each column of the matrix A is obtained by repeatedly shifting the first column.
Figure 2007020943

よって、行列Aを式(2)のように書き換えることができる。

Figure 2007020943
Therefore, the matrix A can be rewritten as in equation (2).
Figure 2007020943

式(2)では、行列Aの1列目を(a,a,…,a)と示し、式(1)に基づいて、循環行列Aを構成する。 In the equation (2), the first column of the matrix A is denoted as (a 1 , a 2 ,..., A n ), and the circulant matrix A is configured based on the equation (1).

パイロットXに対して、s=IFFT(X)とすると、sをパイロットXの時間領域での表記または時間領域のパイロットと呼ぶ。rを時間領域において受信した信号(CP(Cyclic Prefix)除去後)とすると、SISO−OFDMの時間領域モデルを式(3)で示すことができる。

Figure 2007020943
When s = IFFT (X) for pilot X, s is referred to as pilot X time domain notation or time domain pilot. If r is a signal received in the time domain (after CP (Cyclic Prefix) removal), the SISO-OFDM time domain model can be expressed by Equation (3).
Figure 2007020943

式(3)におけるパイロット系列Sは式(4)に示され、hは式(5)に示され、wは時間領域雑音である。

Figure 2007020943
Figure 2007020943
The pilot sequence S in Equation (3) is shown in Equation (4), h is shown in Equation (5), and w is time domain noise.
Figure 2007020943
Figure 2007020943

そして、式(3)における受信信号rは、sとチャネル行列hとを循環畳み込みした後に雑音wを足したものである。   The received signal r in equation (3) is obtained by adding noise w after circular convolution of s and channel matrix h.

図1はCPが挿入されたOFDMシンボルを示す図である。CP長はチャネルの最大マルチパス遅延長より長いものとする。よって、受信側では、ガード間隔を除去して得られる信号はy(n)であり、そしてFFT(Fast Fourier Transform)処理により、時間領域の信号y(n)を周波数領域の信号Y(k)に変換する。信号Y(k)を次式(6)に示す。ただしFFTのポイント数はNである。

Figure 2007020943
FIG. 1 is a diagram illustrating an OFDM symbol with a CP inserted therein. The CP length is longer than the maximum multipath delay length of the channel. Therefore, on the receiving side, the signal obtained by removing the guard interval is y (n), and the time domain signal y (n) is converted into the frequency domain signal Y (k) by FFT (Fast Fourier Transform) processing. Convert to The signal Y (k) is shown in the following equation (6). However, the number of points of FFT is N.
Figure 2007020943

よって、OFDMシステムの周波数領域における受信モデルは式(7)に示される。

Figure 2007020943
Therefore, the reception model in the frequency domain of the OFDM system is expressed by Equation (7).
Figure 2007020943

式(7)では、XはN×Nの対角行列で、その対角線要素は周波数領域のパイロットシンボルであり、Hは周波数領域でのチャネル応答でN×1のベクトルであり、Wは周波数領域での雑音N×1のベクトルである。   In Equation (7), X is an N × N diagonal matrix, the diagonal elements are frequency domain pilot symbols, H is the channel response in the frequency domain, and is an N × 1 vector, and W is the frequency domain. Is a vector of noise N × 1.

例えば、R. Negi and J. Coiffi, "Pilot tone selection for channel estimation in
a mobile OFDM system", IEEE Trans. Consumer Electron., vol.44, pp. 1122-1128, Aug. 1998では、チャネル行列のLS(Least Square)推定(最小平均平方推定)が式(8)と(9)とにより示されている。

Figure 2007020943
Figure 2007020943
For example, R. Negi and J. Coiffi, "Pilot tone selection for channel estimation in
a mobile OFDM system ", IEEE Trans. Consumer Electron., vol. 44, pp. 1122-1128, Aug. 1998, the LS (Least Square) estimation (minimum mean square estimation) of the channel matrix is 9).
Figure 2007020943
Figure 2007020943

よって、チャネル推定を行う際に、まず受信した信号の周波数成分とパイロットの周波数成分を利用し、式(9)に基づいて、チャネルの周波数領域の応答HLSを得る。そして、HLSをIFFT変換し、チャネルの時間領域の応答hLSを得る。hLSを次式(10)で示す。

Figure 2007020943
Therefore, when performing channel estimation, first, the frequency component of the received signal and the frequency component of the pilot are used to obtain a response H LS in the frequency domain of the channel based on Equation (9). Then, the H LS is IFFT transformed to obtain a channel time-domain response h LS . h LS is represented by the following formula (10).
Figure 2007020943

次いで、N個のチャネル推定値の中から、推定値の大きいチャネルを選択し、他のチャネル推定値を「0」にして出力することによって、時間領域での新たなチャネル応答hLSを取得する。例えば、N=64とすると、hLS=[1.2,0.8,0.5,0.001,0.002,0.0011,…,0.0012]が得られる。この場合、選択処理によって得られたチャネル応答はhLS=[1.2,0.8,0.5,0,…,…]となり、すなわち、チャネル数は3個であると考える。 Next, a channel with a large estimated value is selected from the N channel estimated values, and another channel estimated value is set to “0” and output, thereby obtaining a new channel response h LS in the time domain. . For example, if N = 64, h LS = [1.2,0.8,0.5,0.001,0.002,0.0011, ..., 0.0012] is obtained. In this case, the channel response obtained by the selection process is h LS = [1.2, 0.8, 0.5, 0,...], That is, the number of channels is considered to be three.

なお、N個のチャネル推定値を閾値と比較する方法により、閾値より大きいチャネル推定値の数をチャネルの個数にしてもよい。この閾値は、予め決められた固定値でも、リアルタイムに計算された数値でもよく、例えば、推定された雑音平均電力のルート値のc倍(cは定数)と設定する閾値である。   Note that the number of channel estimation values larger than the threshold value may be the number of channels by a method of comparing the N channel estimation values with the threshold value. This threshold value may be a predetermined fixed value or a numerical value calculated in real time. For example, the threshold value is set to c times the estimated noise average power route value (c is a constant).

最後に、FFTにより新しいHLS=FFT(hLS)を得る。上述のような処理によって、一部の雑音影響が除去された。 Finally, a new H LS = FFT (h LS ) is obtained by FFT. A part of noise influence was removed by the above processing.

<FFTの循環シフト性>
X=fft(x)とすると、xはベクトルの時間領域成分であり(sを行ベクトルとする)、Xはベクトルの周波数領域成分であり、それぞれ次式(12)と(13)に示される。

Figure 2007020943
Figure 2007020943
<Circular shift of FFT>
When X = fft (x), x is a time domain component of the vector (s is a row vector), and X is a frequency domain component of the vector, which are represented by the following equations (12) and (13), respectively. .
Figure 2007020943
Figure 2007020943

xを右方向へnシンボル循環シフトすることをx{(k+n)}と定義し、式(11)と式(12)とに基づいて、x{(k+n)}に対してフーリエー変換を行うと、周波数領域成分がej2πmn/NX[m]となる。すなわち、時間領域において循環シフトすることが、周波数領域での結果として位相シフトとなる。 to n symbols cyclic shift to the right direction x is defined as x {(k + n) N}, based on the equations (11) Equations (12), the Fourier chromatography transform to x {(k + n) N} When this is done, the frequency domain component becomes ej2πmn / N X [m]. That is, a cyclic shift in the time domain results in a phase shift as a result in the frequency domain.

式(3)では、雑音を無視すると、s=[s,s,…,s]が得られるため、式(3)を次式(13)に書き換えることができる。

Figure 2007020943
In Expression (3), if noise is ignored, s = [s 1 , s 2 ,..., S N ] is obtained, and therefore Expression (3) can be rewritten into the following Expression (13).
Figure 2007020943

式(13)では、rはhとsとの循環畳み込みを示す。そして、Y=fft(r)、H=fft(h)、X=fft(s)とすると、次式(14)に書き換えることができる。

Figure 2007020943
In Equation (13), r indicates a circular convolution of h and s. When Y = fft (r), H = fft (h), and X = fft (s), the following equation (14) can be rewritten.
Figure 2007020943

すなわち、時間領域の畳み込み計算が周波数領域の掛け算に相当する。よって、式(15)が得られる。

Figure 2007020943
That is, the convolution calculation in the time domain corresponds to multiplication in the frequency domain. Therefore, Formula (15) is obtained.
Figure 2007020943

式(15)より、チャネルの周波数領域値H(m)が得られる。   From the equation (15), the frequency domain value H (m) of the channel is obtained.

送信側では、パイロット系列Sを右方向へnシンボル循環シフトしたもの、つまりs{(k+n)}を送信する。一方、受信側ではsが送信されると想定する場合、次式(16)に示すようになる。

Figure 2007020943
On the transmission side, the pilot sequence S is cyclically shifted by n symbols to the right, that is, s {(k + n) N } is transmitted. On the other hand, when it is assumed that s is transmitted on the receiving side, the following equation (16) is obtained.
Figure 2007020943

すなわち、チャネルを推定する場合、sが送信されたと想定して推定を行うため、次式(17)になる。

Figure 2007020943
That is, when estimating the channel, since estimation is performed assuming that s is transmitted, the following equation (17) is obtained.
Figure 2007020943

式(17)を時間領域に変換すると、式(18)が得られる。

Figure 2007020943
When Expression (17) is converted into the time domain, Expression (18) is obtained.
Figure 2007020943

(実施の形態1)
図2はOFDMシンボルの構成を示す図である。OFDMシンボルでは、Nポイントをチャネル推定用のパイロットとし、その後のMポイントは有効データである。N+Mポイントの間では、チャネルが変化しない(M≧N)と考える。これは、スローフェージングにおいては、チャネルの変化が緩やかであるため、ある一定の時間内ではチャネルが変化しないと見なすことができるからである。
(Embodiment 1)
FIG. 2 is a diagram showing the structure of the OFDM symbol. In the OFDM symbol, N points are pilots for channel estimation, and the subsequent M points are valid data. The channel does not change between N + M points (M ≧ N). This is because in slow fading, the channel changes slowly, so that it can be considered that the channel does not change within a certain period of time.

よって、1つのOFDMデータシンボルを延ばす(Mポイント)ことができる。チャネルを推定する場合、チャネルの長さが一定範囲内のものに対して、小さめのシンボル(Nポイント)を用いて当該チャネルの時間領域応答を推定してもよい。   Therefore, one OFDM data symbol can be extended (M points). When estimating a channel, the time domain response of the channel may be estimated using a smaller symbol (N points) for a channel having a length within a certain range.

図3は本発明の実施の形態1に係るOFDM通信方法の送信側における動作フロー図である。実施の形態1では、基本パイロット系列Sの周波数領域成分をs,s,s,sと定義し、基本パイロット系列Sの時間領域成分をa,a,a,aと定義する。送信側と受信側ともに複数の重み係数を一定の順番で記憶し(S310)、例えば、M個の重み係数P=[p11,p21,…,pk1],P=[p12,p22,…,pk2],…,P=[p1M,p2M,…,pkM]を1〜Mの順番で記憶する。 FIG. 3 is an operation flowchart on the transmission side of the OFDM communication method according to Embodiment 1 of the present invention. In the first embodiment, the frequency domain components of the basic pilot sequence S are defined as s 1 , s 2 , s 3 , s 4, and the time domain components of the basic pilot sequence S are a 1 , a 2 , a 3 , a 4. It is defined as A plurality of weighting factors are stored in a certain order on both the transmitting side and the receiving side (S310). For example, M weighting factors P 1 = [p 11 , p 21 ,..., P k1 ], P 2 = [p 12 , p 22 ,..., p k2 ],..., P M = [p 1M , p 2M ,..., p kM ] are stored in the order of 1 to M.

次いで、変調後の信号(OFDM信号)のPAPRを低減するために、OFDM信号のPAPRを最小にする最適な重み係数を選択する(S320)。具体的には、データベクトルをX=[x,x,…,xN−1]とし、このデータベクトルXをk組のサブベクトルに分割する。例えば、N=8,k=2である場合、それぞれ4つの値からなる2組に分割される。xからxを一組として第1サブベクトルXを形成し、X=[x,x,x,x,0,0,0,0]となる。そして、xからxを一組として第2ベクトルXを形成し、X=[0,0,0,0,x,x,x,x]となる。 Next, in order to reduce the PAPR of the modulated signal (OFDM signal), an optimal weighting coefficient that minimizes the PAPR of the OFDM signal is selected (S320). Specifically, the data vector is X = [x 0 , x 1 ,..., X N−1 ], and the data vector X is divided into k subvectors. For example, when N = 8 and k = 2, it is divided into two sets each consisting of four values. The first subvector X 1 is formed by taking x 0 to x 3 as a set, and X 1 = [x 0 , x 1 , x 2 , x 3 , 0, 0, 0, 0]. Then, the second vector X 2 is formed with x 4 to x 7 as a set, and X 2 = [0, 0, 0, 0, x 4 , x 5 , x 6 , x 7 ].

そして、2つのサブベクトルを次式(19)に示すように組み合わせる。

Figure 2007020943
Then, the two subvectors are combined as shown in the following equation (19).
Figure 2007020943

ij(1≦i≦k,1≦j≦M)は重み係数であり、Pij=exp(jθij)、且つ、θij∈[0,2π]を満たす。 P ij (1 ≦ i ≦ k, 1 ≦ j ≦ M) is a weighting coefficient, and satisfies P ij = exp (jθ ij ) and θ ij ε [0, 2π].

そして、X´ jに対してIDFTを行い、x´ j=IDFT(X´ j)が得られる。式(17)に基づいて、IDFTの線形特性を用いて、次式(20)を得る。

Figure 2007020943
Then, 'performs IDFT with respect to j, x' X j = IDFT (X 'j) is obtained. Based on the equation (17), the following equation (20) is obtained using the linear characteristics of the IDFT.
Figure 2007020943

式(20)により計算されたX´ jのピーク値が最小となるものを時間領域の送信信号とする。この場合、対応する重み係数Pijが最適な重み係数となる。ただし、最適な重み係数は次式(21)を満たす。

Figure 2007020943
A signal having a minimum peak value of X j calculated by Expression (20) is defined as a time domain transmission signal. In this case, the corresponding weight coefficient P ij is the optimum weight coefficient. However, the optimum weight coefficient satisfies the following equation (21).
Figure 2007020943

式(21)により、M個の重み係数から最適な重み係数を選択する。   The optimum weighting factor is selected from the M weighting factors according to the equation (21).

M個の重み係数からQ番目(Q≦M)の重み係数を選択した場合、パイロット系列を(Q−1)×Bだけ循環シフトする(S330)。ただしB≧1(整数)である。そして、シフトされたパイロット系列を時間領域において送信データに挿入する。   When the Qth (Q ≦ M) weighting factor is selected from the M weighting factors, the pilot sequence is cyclically shifted by (Q−1) × B (S330). However, B ≧ 1 (integer). Then, the shifted pilot sequence is inserted into transmission data in the time domain.

例えば、選択した重み係数の順位が2番目である場合、基本パイロット系列を左へ2つシフトして、a,a,a,aとし、これを時間領域において送信データに挿入する。すなわち、挿入したパイロット系列は式(1)に示された行列Aの3列目である。 For example, when the rank of the selected weighting factor is the second, the basic pilot sequence is shifted two times to the left to become a 3 , a 4 , a 1 , a 2 , which are inserted into the transmission data in the time domain . That is, the inserted pilot sequence is the third column of the matrix A shown in Equation (1).

図4は、本発明の実施の形態1に係るOFDM通信方法の受信側における動作フロー図である。まず式(9)を用いてチャネル推定を行い(S410)、チャネルの周波数領域における応答HLSを推定する。そして、応答HLSに対してIFFTを行い、チャネルの時間領域における応答hLSを得る。次いで、閾値処理を行って雑音影響を除去する。そして、雑音が除去された時間領域における応答に対してFFTを行い、フィルタ(雑音除去)処理後の周波数領域における応答を得る。 FIG. 4 is an operation flowchart on the receiving side of the OFDM communication method according to Embodiment 1 of the present invention. First, channel estimation is performed using equation (9) (S410), and the response H LS in the frequency domain of the channel is estimated. Then, it performs an IFFT on the response H LS, get response h LS in the time domain of the channel. Next, threshold processing is performed to remove noise effects. Then, FFT is performed on the response in the time domain from which noise is removed, and a response in the frequency domain after the filter (noise removal) processing is obtained.

次いで、チャネル応答があった第一パスの位置を確認する(S420)。具体的には、式(17)と(18)より、時間領域のパイロット系列が左へシフトされた分、推定されたチャネル値も同じく左へシフトされていることが分かる。例えば、2パスのチャネルでは、時間領域の基本パイロット系列a,a,a,aが挿入される場合、推定されるチャネル値はh,h,0,0となり、換言すると、第一パスが一番目に位置する。一方、時間領域の基本パイロット系列が左へ1つシフトされてパイロット系列a,a,a,aが挿入される場合、推定されるチャネル値はh,0,0,hとなり、つまり第一パスが二番目に位置する。また、時間領域の基本パイロット系列が左へ2つシフトされてパイロット系列a,a,a,aが挿入される場合、推定されるチャネル値は0,0,h,hとなり、つまり第一パスが三番目に位置する。 Next, the position of the first path where there is a channel response is confirmed (S420). Specifically, from equations (17) and (18), it can be seen that the estimated channel value is also shifted to the left as the time-domain pilot sequence is shifted to the left. For example, in a two-path channel, when time domain basic pilot sequences a 1 , a 2 , a 3 , a 4 are inserted, the estimated channel values are h 1 , h 2 , 0, 0, in other words The first pass is located first. On the other hand, when the basic pilot sequence in the time domain is shifted one to the left and the pilot sequences a 2 , a 3 , a 4 , a 1 are inserted, the estimated channel values are h 2 , 0, 0, h 1. That is, the first pass is located second. In addition, when the basic pilot sequence in the time domain is shifted to the left by two and pilot sequences a 3 , a 4 , a 1 , a 2 are inserted, the estimated channel values are 0, 0, h 1 , h 2. That is, the first path is located third.

次いで、確定された第一パスの位置により、送信側で選択された重み係数を確定する(S430)。図5を参照しながら説明する。図5は第一パスの位置により重み係数を確定するステップの概念図である。推定されたチャネル値を右へ順番に、Z,…,Z,0,…,0(後尾がすべて「0」となる、図6を参照)が得られるまでシフトする。本実施の形態では、右へ2つシフトすることで、h,h,0,0が得られるため、受信側は、送信側で2番目の重み係数Pを選択してPAPRを低減したと確定することができる。例えばP=[p12,p22]=[1,−1]である。 Next, the weighting factor selected on the transmission side is determined based on the determined position of the first path (S430). This will be described with reference to FIG. FIG. 5 is a conceptual diagram of the steps for determining the weighting coefficient based on the position of the first pass. The estimated channel values are shifted in order to the right until Z,..., Z, 0,..., 0 (the tail is all “0”, see FIG. 6) are obtained. In the present embodiment, h 1 , h 2 , 0, 0 are obtained by shifting two to the right, so that the receiving side selects the second weighting factor P 2 on the transmitting side to reduce PAPR. Can be confirmed. For example, P 2 = [p 12 , p 22 ] = [1, −1].

なお、送信側のパイロット系列及び受信側のチャネル値を1つずつシフトするとは限らず、シフトペースである定数Bを1より大きい整数としてもよい。B>1とすることにより、チャネル推定に誤りがある場合(例えば、雑音の影響により、推定されるマルチパスのチャネルが実際より多い、または、実際より少ない場合)でも、より正確にシフト量を推定することができ、より正確に重み係数を確定することができる。   Note that the pilot sequence on the transmission side and the channel value on the reception side are not necessarily shifted one by one, and the constant B, which is the shift pace, may be an integer greater than one. By setting B> 1, even when there is an error in channel estimation (for example, when there are more or fewer estimated multipath channels due to the influence of noise), the shift amount can be more accurately determined. The weighting factor can be determined more accurately.

その後、受信信号からCPを除去してy=[y,y,…,y]が得られ、yをIFFTすると次式(22)が得られる。

Figure 2007020943
Thereafter, CP is removed from the received signal to obtain y = [y 1 , y 2 ,..., Y N ], and y is IFFT to obtain the following equation (22).
Figure 2007020943

周波数領域のチャネル値をH=[h,h,…,h]とすると、Z=Y/h/p12(i=0,1,2,3)、Z=Y/h/p22(i=4,5,6,7)となる。すなわち、周波数領域の受信信号のチャネル影響を除去して(Y/h)、周波数領域において実際に送信された信号(雑音を含む)が得られる。ここで得られた信号は、有効送信信号(有効データ)に重み係数がかけられた信号(線形組み合わせ信号)であるため、確定された重み係数Pで割る処理を行う。最後に、Zを復調してデータを復元する。 If the channel values in the frequency domain are H = [h 1 , h 2 ,..., H N ], Z i = Y i / h i / p 12 (i = 0,1,2,3), Z i = Y i / h i / p 22 (i = 4, 5, 6, 7). That is, the channel effect of the received signal in the frequency domain is removed (Y i / h i ), and a signal (including noise) actually transmitted in the frequency domain is obtained. Signals obtained herein are the effective transmission signal (valid data) weighting factor is applied to the signal (linear combination signal), it performs a process of dividing the determined weighting factor P 2. Finally, Z 1 is demodulated to restore the data.

このように、本実施の形態によれば、基本パイロット系列を定義し、予め複数の重み係数を所定の順番で、それぞれ送信側と受信側とに記憶させる。送信側では、PAPRを低減する過程において使用された最適な重み係数の順位を確定し、その順位に応じて基本パイロット系列をシフトし、シフトされたパイロット系列を時間領域において送信データに挿入する。受信側では、FFTの特性を利用して、受信信号のパイロット系列を用いて推定されたチャネル値の第一パスの位置によって、送信側で使用された最適な重み係数を確定する。このため、本実施の形態によれば、送信側で使用された重み係数を受信側で確定するための補助情報を別途送信する必要がないため、伝送効率を向上することができるとともに、補助情報の誤りによる複合誤りを防止して受信データの誤り率を低下させることができる。   Thus, according to the present embodiment, a basic pilot sequence is defined, and a plurality of weighting factors are stored in advance in a predetermined order on the transmitting side and the receiving side, respectively. On the transmission side, the rank of the optimum weighting factor used in the process of reducing the PAPR is determined, the basic pilot sequence is shifted according to the rank, and the shifted pilot sequence is inserted into the transmission data in the time domain. On the receiving side, the optimum weighting factor used on the transmitting side is determined based on the position of the first path of the channel value estimated using the pilot sequence of the received signal using the characteristics of the FFT. Therefore, according to the present embodiment, it is not necessary to separately transmit auxiliary information for determining the weighting factor used on the transmission side on the reception side, so that transmission efficiency can be improved and auxiliary information It is possible to reduce the error rate of the received data by preventing the compound error due to the error.

(実施の形態2)
図7は本発明の実施の形態2に係るOFDM通信方法の送信側における動作フロー図である。本実施の形態では、周波数領域のパイロット系列を分割して重み係数にするため、予め送信側と受信側とに複数の重み係数を記憶させる必要がない。
(Embodiment 2)
FIG. 7 is an operation flowchart on the transmission side of the OFDM communication method according to Embodiment 2 of the present invention. In the present embodiment, since the pilot sequence in the frequency domain is divided into weighting factors, it is not necessary to store a plurality of weighting factors in advance on the transmitting side and the receiving side.

上述のように、チャネルを推定する過程において、チャネルのLS推定性能を最良にするため、式(3)のS行列は直交行列である必要がある。すなわち時間領域のパイロット系列がシフト直交である必要がある。以上の条件を満たすブロック状パイロットに対して、周波数領域のパイロット系列が等電力であることが要求される。周波数領域のパイロットシンボルの電力を1とすると、周波数領域のパイロットがe(θは0〜2πの任意値)であるため、上述の条件を満たす。よって、周波数領域のパイロット系列を分割して重み係数にすることができる。 As described above, in order to optimize the LS estimation performance of the channel in the process of estimating the channel, the S matrix of Equation (3) needs to be an orthogonal matrix. That is, the time-domain pilot sequence needs to be shift orthogonal. For block pilots that satisfy the above conditions, the pilot sequences in the frequency domain are required to be of equal power. Assuming that the power of the pilot symbol in the frequency domain is 1, the frequency domain pilot is e (θ is an arbitrary value from 0 to 2π), so the above condition is satisfied. Therefore, the frequency domain pilot sequence can be divided into weighting factors.

例えば、N=4とし、周波数領域のパイロット系列をK=2個に分割し、基本パイロット系列をs,s,s,sと定義する。基本パイロット系列の時間領域成分はa,a,a,aである。s,sを1つ目の重み係数とし、s,sを2つ目の重み係数とする(S710)。 For example, N = 4, the frequency domain pilot sequence is divided into K = 2, and the basic pilot sequences are defined as s 1 , s 2 , s 3 , s 4 . The time domain components of the basic pilot sequence are a 1 , a 2 , a 3 , a 4 . s 1 and s 2 are set as the first weighting factors, and s 3 and s 4 are set as the second weighting factors (S710).

ついで、実施の形態1と同様方法で、信号のPAPRを最小にする最適な重み係数を選択する(S720)。   Next, an optimum weighting coefficient that minimizes the PAPR of the signal is selected by the same method as in the first embodiment (S720).

選択された最適重み係数の順位に応じて基本パイロット系列をシフトする(S730)。選択された重み係数が2番目の重み係数s,sである場合、基本パイロット系列を左へ2つシフトする。すなわち、挿入されるパイロット系列が式(1)に示された行列Aの3列目a,a,a,aとなる。 The basic pilot sequence is shifted according to the order of the selected optimum weight coefficient (S730). If the selected weighting factor is the second weighting factor s 3 , s 4 , the basic pilot sequence is shifted two times to the left. That is, the inserted pilot sequence is the third column a 3 , a 4 , a 1 , a 2 of the matrix A shown in Equation (1).

図8は本発明の実施の形態2に係るOFDM通信方法の受信側における動作フロー図で
ある。図8に示すように、式(9)によりチャネル推定を行い(S810)、式(9)の周波数領域のパイロット系列は基本パイロット系列である。すなわち、式(9)では、X(1)=s、X(2)=s、X(3)=s、X(4)=sである。
FIG. 8 is an operation flowchart on the receiving side of the OFDM communication method according to Embodiment 2 of the present invention. As shown in FIG. 8, channel estimation is performed using equation (9) (S810), and the pilot sequence in the frequency domain of equation (9) is a basic pilot sequence. That is, in the formula (9), X (1) = s 1 , X (2) = s 2 , X (3) = s 3 , X (4) = s 4 .

次いで、チャネル応答があった第一パスの位置を確認する(S820)。確定された第一パスの位置により、送信側で選択されたパイロットを確定する(S830)。実施の形態1と同様方法により、チャネル推定値を右へ2つシフトすると、h,h,0,0が得られるため、送信側で選択されたパイロット系列が2番目であると確定することができる。つまり、データ復元に必要な重み係数はs,sであると確定することができる。 Next, the position of the first path where there is a channel response is confirmed (S820). Based on the determined position of the first path, the pilot selected on the transmission side is determined (S830). When the channel estimation value is shifted two times to the right by the same method as in the first embodiment, h 1 , h 2 , 0, 0 are obtained, so that the pilot sequence selected on the transmission side is determined to be the second. be able to. That is, it can be determined that the weighting factors necessary for data restoration are s 3 and s 4 .

このように、本実施の形態によれば、送信側では、周波数領域のパイロット系列を分割して重み係数とし、PAPRを低減する過程において選択する重み係数が複数に分割されたパイロット系列での順位を確定し、その順位に応じて基本パイロット系列をシフトして、シフトされたパイロット系列を時間領域において送信データに挿入する。受信側では、FFTの特性を利用して、チャネル応答があった第一パスの位置を確定することによって、送信側で選択された最適な重み係数を確定する。よって、本実施の形態によれば、特定の重み係数を予め送信側と受信側とに記憶させる必要がないため、実施の形態1よりも簡易なOFDM通信方法を提供することができる。   As described above, according to the present embodiment, on the transmitting side, the pilot sequence in the frequency domain is divided into weighting factors, and the order in the pilot sequence in which the weighting factor selected in the process of reducing the PAPR is divided into a plurality of Is determined, the basic pilot sequence is shifted according to the order, and the shifted pilot sequence is inserted into the transmission data in the time domain. On the receiving side, the optimum weighting factor selected on the transmitting side is determined by determining the position of the first path where there was a channel response using the characteristics of the FFT. Therefore, according to the present embodiment, there is no need to store a specific weighting factor in advance on the transmitting side and the receiving side, so that it is possible to provide a simpler OFDM communication method than in the first embodiment.

以上の説明は、典型的な実施の形態について説明したが、本発明の主旨と範囲から逸脱しない限り、種々の変更、置き換え及び追加をすることができる。そのため、本発明は上記の実施の形態に拘らず、請求の範囲及びそれに均等する範囲によって限定されるものである。   While the above description has been given of exemplary embodiments, various changes, substitutions and additions can be made without departing from the spirit and scope of the present invention. Therefore, the present invention is not limited to the above-described embodiments, but is limited by the scope of the claims and the scope equivalent thereto.

本明細書は、2005年8月16日出願の中国出願番号200510091685.6に基づくものである。その内容は、全てここに含めておく。   This specification is based on Chinese application No. 2005000915685.6 filed on Aug. 16, 2005. The contents are all included here.

本発明は、送信信号のピーク電力を抑える必要のあるOFDM通信等に好適である。   The present invention is suitable for OFDM communication or the like that needs to suppress the peak power of a transmission signal.

CPが挿入されたOFDMシンボルを示す図The figure which shows the OFDM symbol in which CP was inserted OFDMシンボルの構成を示す図Diagram showing the structure of an OFDM symbol 本発明の実施の形態1に係るOFDM通信方法の送信側における動作フロー図Operation flow diagram on transmission side of OFDM communication method according to embodiment 1 of the present invention 本発明の実施の形態1に係るOFDM通信方法の受信側における動作フロー図Operation flow diagram on receiving side of OFDM communication method according to embodiment 1 of the present invention 本発明の実施の形態1に係る第一パスの位置により重み係数を確定するステップの概念図The conceptual diagram of the step which determines a weighting coefficient by the position of the 1st path | pass which concerns on Embodiment 1 of this invention. 本発明の実施の形態1に係るチャネル推定を示す模式図Schematic diagram showing channel estimation according to Embodiment 1 of the present invention. 本発明の実施の形態2に係るOFDM通信方法の送信側における動作フロー図Operation flow diagram on transmission side of OFDM communication method according to embodiment 2 of the present invention 本発明の実施の形態2に係るOFDM通信方法の受信側における動作フロー図Operation flow diagram on receiving side of OFDM communication method according to embodiment 2 of the present invention

Claims (7)

複数の重み係数のうちから、OFDM信号のPAPRを最小にする最適な重み係数を選択する選択ステップと、
前記最適な重み係数に対応するシフト量にてパイロット系列を循環シフトするシフトステップと、
循環シフトされた前記パイロット系列を送信データに挿入する挿入ステップと、
を具備する送信側のOFDM通信方法。
A selection step of selecting an optimum weighting factor that minimizes the PAPR of the OFDM signal from among a plurality of weighting factors;
A shift step of cyclically shifting a pilot sequence by a shift amount corresponding to the optimal weighting factor;
An insertion step of inserting the cyclically shifted pilot sequence into transmission data;
A transmission-side OFDM communication method comprising:
前記選択ステップにおいて、
所定の順番で予め記憶されている前記複数の重み係数のうちから前記最適な重み係数を選択し、
前記シフトステップにおいて、
前記複数の重み係数における前記最適な重み係数の順位に応じたシフト量にて前記パイロット系列を循環シフトする、
請求項1記載のOFDM通信方法。
In the selection step,
Selecting the optimum weighting factor from the plurality of weighting factors stored in advance in a predetermined order;
In the shifting step,
Cyclically shifting the pilot sequence by a shift amount corresponding to the order of the optimum weighting factors in the plurality of weighting factors;
The OFDM communication method according to claim 1.
データベクトルを複数のサブベクトルに分割する分割ステップ、をさらに具備し、
前記選択ステップは、
前記複数の重み係数と前記複数のサブベクトルとを順次乗算して複数の線形組み合わせを得るステップと、
前記複数の線形組み合わせに対して逆フーリエ変換を行って複数の時間領域信号を得るステップと、
前記複数の時間領域信号のうちからPAPRが最小の信号を確定するステップと、
前記確定された信号が対応する線形組み合わせに使用された前記最適な重み係数を選択するとともに、前記複数の重み係数における前記最適な重み係数の順位を確定するステップと、を備える、
請求項2記載のOFDM通信方法。
Further comprising a dividing step of dividing the data vector into a plurality of subvectors;
The selection step includes
Sequentially multiplying the plurality of weighting factors and the plurality of subvectors to obtain a plurality of linear combinations;
Performing inverse Fourier transform on the plurality of linear combinations to obtain a plurality of time domain signals;
Determining a signal with the smallest PAPR from the plurality of time domain signals;
Selecting the optimal weighting factor used for the linear combination to which the determined signal corresponds, and determining the ranking of the optimal weighting factor in the plurality of weighting factors.
The OFDM communication method according to claim 2.
前記シフトステップにおいて、
前記複数の重み係数における前記最適な重み係数の順位Qに応じて、前記パイロット系列を(Q−1)×B(Bは1以上の整数)だけ循環シフトする、
請求項2記載のOFDM通信方法。
In the shifting step,
The pilot sequence is cyclically shifted by (Q−1) × B (B is an integer of 1 or more) according to the order Q of the optimum weighting factors in the plurality of weighting factors.
The OFDM communication method according to claim 2.
前記複数の重み係数は、循環シフト前の前記パイロット系列を複数に分割して生成されるものである、
請求項1記載のOFDM通信方法。
The plurality of weighting factors are generated by dividing the pilot sequence before cyclic shift into a plurality of parts,
The OFDM communication method according to claim 1.
請求項1記載のOFDM通信方法によりシフトされ、前記送信データに挿入された前記パイロット系列を受信する受信ステップと、
受信した前記パイロット系列を用いてチャネル推定を行って第1パスの位置を確定する確定ステップと、
複数の重み係数のうちから前記第1パスの位置により得られるシフト量と対応する重み係数を、請求項1記載のOFDM通信方法により選択された前記最適な重み係数として選択する選択ステップと、
を具備する受信側のOFDM通信方法。
Receiving a pilot sequence shifted by the OFDM communication method according to claim 1 and inserted in the transmission data;
A confirming step of performing channel estimation using the received pilot sequence to determine the position of the first path;
A selection step of selecting a weighting factor corresponding to a shift amount obtained from the position of the first path from among a plurality of weighting factors as the optimum weighting factor selected by the OFDM communication method according to claim 1;
A receiving-side OFDM communication method comprising:
前記選択ステップにおいて選択される前記重み係数が前記パイロット系列の一部である、
請求項6記載のOFDM通信方法。
The weighting factor selected in the selection step is part of the pilot sequence;
The OFDM communication method according to claim 6.
JP2007531011A 2005-08-16 2006-08-16 OFDM communication method Ceased JPWO2007020943A1 (en)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
CN200510091685.6 2005-08-16
CNA2005100916856A CN1917490A (en) 2005-08-16 2005-08-16 Method of reducing peak even ratio of orthoganal frequency division complexing signals
PCT/JP2006/316104 WO2007020943A1 (en) 2005-08-16 2006-08-16 Ofdm communication method

Publications (1)

Publication Number Publication Date
JPWO2007020943A1 true JPWO2007020943A1 (en) 2009-03-26

Family

ID=37738387

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2007531011A Ceased JPWO2007020943A1 (en) 2005-08-16 2006-08-16 OFDM communication method

Country Status (3)

Country Link
JP (1) JPWO2007020943A1 (en)
CN (1) CN1917490A (en)
WO (1) WO2007020943A1 (en)

Families Citing this family (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101843068A (en) * 2008-01-17 2010-09-22 上海贝尔股份有限公司 A method and apparatus for reducing the peak-to-average-power-ratio (papr) of the output symbol of SC-FDMA transmitter
KR101581956B1 (en) * 2008-10-22 2016-01-04 엘지전자 주식회사 Method and apparatus of transmitting signal in wireless communication system
CN101478526B (en) * 2009-01-14 2012-04-25 北京大学深圳研究生院 OFDM system and partial transmission sequence method thereof
WO2012030319A2 (en) * 2009-06-26 2012-03-08 Hypres, Inc. System and method for controlling combined radio signals
CN101945073B (en) * 2009-07-03 2013-02-27 中兴通讯股份有限公司 Pilot frequency-based time offset estimation device and method
CN101867549B (en) * 2010-06-11 2012-10-24 西安电子科技大学 Suppression method for OFDM signal peak-to-average power ratio based on positive semi-definite programming technique
CA2931398A1 (en) * 2013-11-26 2015-06-04 Plusn, Llc System and method for radio frequency carrier aggregation
CN106597384B (en) * 2016-12-21 2019-02-12 中国航空工业集团公司雷华电子技术研究所 A kind of weighing apparatus envelope waveform generation circuit based on OFDM modulation
CN113940040A (en) * 2019-06-13 2022-01-14 株式会社Ntt都科摩 Communication device
EP4074102A4 (en) * 2019-12-13 2023-07-26 ZTE Corporation Time-domain modulation scheme for low peak average power ratio
CN114465860B (en) * 2022-01-11 2023-06-23 珠海中慧微电子有限公司 Peak-average-power-ratio reducing method and device for OFDM (orthogonal frequency division multiplexing) signals and storage medium

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH08274748A (en) * 1995-03-31 1996-10-18 Victor Co Of Japan Ltd Frequency-division multiplex signal generating device
JP2004147126A (en) * 2002-10-24 2004-05-20 Osaka Industrial Promotion Organization Multi-carrier communication method and multi-carrier communication device
JP2005027177A (en) * 2003-07-04 2005-01-27 Mitsubishi Electric Corp Radio transmitter
JP2005524278A (en) * 2002-04-26 2005-08-11 サムスン エレクトロニクス カンパニー リミテッド Apparatus and method for transmitting / receiving additional information of partial transmission sequence in orthogonal frequency division multiplexing communication system
JP2005341056A (en) * 2004-05-25 2005-12-08 Ntt Docomo Inc Transmitter and receiver

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH08274748A (en) * 1995-03-31 1996-10-18 Victor Co Of Japan Ltd Frequency-division multiplex signal generating device
JP2005524278A (en) * 2002-04-26 2005-08-11 サムスン エレクトロニクス カンパニー リミテッド Apparatus and method for transmitting / receiving additional information of partial transmission sequence in orthogonal frequency division multiplexing communication system
JP2004147126A (en) * 2002-10-24 2004-05-20 Osaka Industrial Promotion Organization Multi-carrier communication method and multi-carrier communication device
JP2005027177A (en) * 2003-07-04 2005-01-27 Mitsubishi Electric Corp Radio transmitter
JP2005341056A (en) * 2004-05-25 2005-12-08 Ntt Docomo Inc Transmitter and receiver

Also Published As

Publication number Publication date
WO2007020943A1 (en) 2007-02-22
CN1917490A (en) 2007-02-21

Similar Documents

Publication Publication Date Title
JPWO2007020943A1 (en) OFDM communication method
Muquet et al. Subspace-based blind and semi-blind channel estimation for OFDM systems
JP4398791B2 (en) Transmitter and transmission control method
US7583755B2 (en) Systems, methods, and apparatus for mitigation of nonlinear distortion
US7609616B2 (en) Sparse channel estimation for orthogonal frequency division multiplexed signals
JP4413232B2 (en) Channel evaluator and method for evaluating channel transfer function, and apparatus and method for supplying pilot sequence
JP5632457B2 (en) Reducing peak-to-average power ratio in multicarrier signals
EP1593246A1 (en) Two-dimensional channel estimation for multicarrier multiple input multiple outpout communication systems
JP6026006B2 (en) Transmitting apparatus, receiving apparatus, and communication system
JP5486734B2 (en) Transmission signal generating apparatus and method in single carrier communication system
JP4780419B2 (en) Transmission path estimation apparatus, equalization apparatus, and radio system
Elavarasan et al. Peak-power reduction using improved partial transmit sequence in orthogonal frequency division multiplexing systems
Sudha et al. Low-complexity modified SLM method for PAPR reduction in OFDM systems
US7450490B2 (en) Channel estimation using the guard interval of a multicarrier signal
JP2018007056A (en) Radio communication system, radio transmission device and radio reception device
Soman et al. Pilot based MMSE channel estimation for spatial modulated OFDM systems
JP6110641B2 (en) Method and device for estimating channel impulse response between source and receiver
KR101128287B1 (en) Ofdm receiver being capable of estimating timing error in multi-path fading channel, ofdm system having the same and timing error estimation method thereof
Ce et al. Low-complexity PTS scheme based on phase factor sequences optimization
Wang et al. Channel estimation and equalization for SS-OOFDM system with high mobility
EP4327523A1 (en) Affine frequency division multiplexing waveforms for doubly dispersive channels
Debels et al. Adaptive modulation and coding for BIC-UFMC and BIC-OFDM systems taking CFO into account
US8917585B2 (en) Method for estimating a received signal and corresponding device
Babulkar A Comprehensive Review on Channel Estimation Techniques in MIMO-OFDM
Rao et al. ICI Cancellation Schemes and Selection Criteria in Orthogonal Frequency Division Multiplexing System: A Review

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20090520

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20111025

A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20111220

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

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20120117

A045 Written measure of dismissal of application

Free format text: JAPANESE INTERMEDIATE CODE: A045

Effective date: 20120529