CN1901518A - Method for channel evaluating of weak energy parallel PN sequence - Google Patents
Method for channel evaluating of weak energy parallel PN sequence Download PDFInfo
- Publication number
- CN1901518A CN1901518A CN 200510035892 CN200510035892A CN1901518A CN 1901518 A CN1901518 A CN 1901518A CN 200510035892 CN200510035892 CN 200510035892 CN 200510035892 A CN200510035892 A CN 200510035892A CN 1901518 A CN1901518 A CN 1901518A
- Authority
- CN
- China
- Prior art keywords
- matrix
- sequence
- channel
- channel estimating
- signal
- 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
Links
Images
Landscapes
- Cable Transmission Systems, Equalization Of Radio And Reduction Of Echo (AREA)
Abstract
This invention provides a method and a device for using parallel PN sequence of weak energy to estimate channels, which utilizes matrix Q determined by the channel S/N ratio by the applied training sequence and the measured channel S/N ratio and matrix z determined by a receiving signal matrix y and a mean matrix A to get the pulse response estiomation value characterizing that the pulse response estimation value is hIMMSEE=Qz, which is a channel estimation technology between the biased and the MMSE estimations and can estimate channel pulse response accurately.
Description
Technical field:
The present invention relates to the method that a kind of weak energy parallel PN sequence carries out channel estimating, belong to the digital mobile communication field of using the electromagnetic wave technology, particularly digital terrestrial television, single carrier ofdm communication system, multi-carrier OFDM communication system, wireless lan (wlan) etc. adopt the channel estimation technique in the OFDM modulated digital communication system.
Background technology:
As everyone knows, the high frequency band radio wave will be through multiple transmission environment reflection such as house, vehicle, high mountain or diffraction in transmission course.That is: the primary signal that receives of big line amplitude difference not only, and phase place also there are differences.The result of these signal amplitude combinations will cause the violent fluctuation of signal amplitude, promptly so-called multipath fading.
In ofdm communication system, this multipath fading wireless channel be frequency selectivity with time selectivity, so, before the ofdm signal demodulation, dynamic channel is estimated it is very important.
Channel estimating can adopt to be inserted training sequence with certain cycle and realizes on the subcarrier of OFDM symbol, also can adopt the mode of inserting training sequence on each OFDM symbol.The former is called block-type pilot channel estimation, is mainly used in the estimation to slow fading channel.The latter is called the Comb Pilot channel estimating, is mainly used in the estimation of fast fading channel.These two kinds of methods all can be further divided into least square (Least Square, LS) method or least mean-square error (Minimum Mean-Square-Error, MMSE) method etc.Yet these methods have all taken useful bandwidth and have reduced data transmission efficiency.
It then is another kind of thinking that training sequence is superimposed upon the method for carrying out channel estimating on the time-domain signal.The bandwidth resources that this method biggest advantage is exactly a preciousness of no use exchange performance for estimating channel for.
The thought that overlying training sequence is used for channel estimating originates from nineteen sixty-five, referring to document Van Duc Nguyen, ChristianHansen and Hans-Peter Kuchenbecher, " Performance of channel estimation using pilot symbols forcoherent OFDM system ", WPMC ' 00 November12-15 2000, Bangkok, Thailand, pp:842-847.
Nineteen ninety-five, this thought was used to digital communication, see document C.E.Kastenholz and W.P.Birkemeier, A simultaneousinformation transfer ans channel sounding modulation technique for wide-band channels, IEEETrans.on Communication Technology, June 1965, pp:162-165.
Thereafter different channel estimation methods also are suggested gradually, see following document:
[1]B.Farhang-Boroujeny,Pilot-based?channel?identification:proposal?for?semi-blind?identificationof?communication?channels,IEE?Electronics?Letters,vol.31,no.13,June?1995,pp:1044-104.
[2]G.T.Zhou,M.Viberg,and?T.Mckelvey,Superimposed?periodic?pilots?for?blind?channelestimation,Proc.35
th?Asilomar?Conference?on?Signals,Systems,and?Computers,PacificGrove,CA,Nov.2001,pp:653-657.
[3]N.Chen?and?G.T.Zhou,A?superimposed?periodic?pilot?scheme?for?semi-blind?channel?estimationof?OFDM?systems,Proc.10
th?IEEE?DSP?Workshop,Pine?Mountain,GA,Oct.2002,pp:362-365.
[4]J.K.Tugnait?and?W.Luo,On?channel?estimation?using?superimposed?training?and?first-orderstatistics,Proc.IEEE?Int.Conf.Acoust,Speech,Signal?Processing,Hong?Kong,China,Apr.2003,pp:624-627.
[5]G.T.Zhong?and?N.Chen,Superimposed?training?for?doubly?selective?channels,Proc.IEEEStatistical?Signal?Processing?Workshop,St.Louis,MO,Sept.2003,pp:73-76.
[6]X.Meng?and?J.K.Tugnait,Semi-blind?channel?estimation?and?detection?using?superimposedtraining,ICASSP,2004,pp:417-420.
[7]Ning?Chen,G.Tong?Zhou,Superimposed?training?for?OFDM:a?peak-to-average?power?ratioanalysis,IEEE?Trans.on?Signal?Processing,submitted,2004.
Though the method performance that these employing parallel PN sequences carry out channel estimating is pretty good,, its computation complexity is still bigger.
Summary of the invention:
The objective of the invention is: propose a kind of channel estimation methods and device that utilizes weak energy parallel PN sequence to carry out orthogonal FDM communication system, can estimate channel parameter effectively, and reduce computation complexity with respect to the channel estimation methods of prior art.
For achieving the above object, according to an aspect of the present invention, a kind of method of utilizing weak energy parallel PN sequence to carry out channel estimating has been proposed.Its technical scheme is: at transmitting terminal, and the PN sequence of the weak energy of stack on a symbol of time domain OFDM modulation signal; With the matrix Q that the channel signal to noise ratio that can survey determines channel estimating is reduced to linear operation by this PN sequence in the receiving terminal utilization, directly obtains channel state parameter.The impulse response estimated value that it is characterized in that channel
Wherein
z=Ay
Y is the received signal matrix.
According to another aspect of the present invention, a kind of device that utilizes weak energy parallel PN sequence to carry out channel estimating has been proposed.Utilization obtains the impulse response estimated value of channel by the matrix Q of this PN sequence and the channel signal to noise ratio decision that can survey, by received signal matrix y and the definite matrix z of average matrix A.The impulse response estimated value that it is characterized in that channel
The method and apparatus that utilizes weak energy parallel PN sequence to carry out channel estimating that the present invention proposes is a kind of between the channel estimation technique that does not have between estimation partially and the MMSE estimation, can estimate channel state parameter effectively and be used for demodulation, thereby effectively improve systematic function.
Channel estimating model of the present invention and principle are described in detail as follows:
For a sub-carrier number is N
cThe system of OFDM modulation, X
i(k) frequency-region signal of k symbol of i piece of expression (block).Suppose that in system protection is longer than channel delay at interval, then do not have inter-block-interference (inter-block-interference, IBI).In order to simplify, can omit piece mark i.Time domain OFDM signal { x (n) } can be expressed as:
The training sequence p (n) of the weak energy of stack on signal has:
Then, length is N
gProtection be filled in the information symbol front at interval, have:
In the formula, G is the gain of power amplifier.General hypothesis G=1.
For frequency selective fading channels, can represent by constant finite impulse response filter h (n) of time spent.So received signal can be expressed as:
Wherein,
In the formula, * represents linear convolution, and h (n) is the impulse response of transmitter filter, frequency-selective channel, filter for receiver, and w (n) is the zero-mean additive noise, and r is the total number in channel fading footpath, h
iBe the multiple impulse response in i bar footpath, f
DiBe the Doppler frequency shift in i bar footpath, λ postpones the expansion label, and T is the sampling period, τ
iIt is the time delay of carrying out normalized i bar footpath by the sampling time.
At receiving terminal, received signal by A/D and low pass filter after, protection is at interval deleted, has:
After eliminating Cyclic Prefix, linear convolution just becomes circular convolution (representing with ), at this moment, has:
(7)
Wherein,
q(n)=p(n)h(n) (8)
u(n)=s(n)h(n)+w(n) (9)
With matrix (7) formula of rewriting be:
y=Ph+u (10)
Wherein, P is a N
c* L circular matrix, it first classifies [p[0], p[1 as] ..., p[N
c-1]]
T, y=[y[0], y[1] ..., y[N
c-1]]
T, h=[h[0], h[1] ..., h[L-1]]
T, u=[u[0], u[1] ..., u[N
c-1]]
T, T represents transposition herein.
For the time frequency-selective channel that becomes, selection cycle is the PN sequence p[n of P] carry out channel estimating.In order to estimate the channel that length is L, require P 〉=L.Simultaneously, for further simplification, suppose
It is an integer.At this moment, size of structure is P * N
cMean matrix
I wherein
PIt is the unit matrix of a P * P.Can get by (10):
Ay=APh+Au (11)
At this moment, suppose z=Ay, then z be P * 1 to y[n] matrix that averages, its i element is:
In like manner, establish u=Au, u be P * 1 to y[n] matrix that averages, its i element is:
Because p[n] be that one-period is the sequence of P, suppose that first column element of the matrix P correspondence in (11) formula is:
Then have:
So (11) formula of can rewriting is:
z=Ay=P
0h+ u (15)
Least mean-square error estimates that (MMSEE) is based on a kind of statistical method that Bayesan estimates, supposes that wherein h is a random vector.Suppose that h obeys the Gaussian Profile of zero-mean, and uncorrelated with u, at this moment, have:
Wherein:
D=P
0 HP
0 (18)
C
h=E{hh
H} (19)
The corresponding channel frequency response of MMSEE is:
Wherein G is a matrix, and its element can be expressed as:
[G]
n,k=e
-j2πnk/N |n|≤N
α,0≤k≤L-1 (21)
That is:
Wherein:
Know that by formula (16) the MMSEE channel estimation methods needs the prior information C of channel
h, the accurate acquisition of this information is difficult for; But also must be to C
hInvert, when channel impulse response length was longer, this amount of calculation also was higher.
Propose a kind of linear minimum mean-squared error algorithm below and estimate channel, reduce computation complexity.
Suppose E{ (hh
H)
-1}=E{|1/h
k|
2I, wherein I is a unit matrix.Definition average signal-to-noise ratio SNR is E{|x
k|
2/ σ
2], can obtain by (16) formula:
Wherein,
B=E{|x
k|
2}·E{|1/h
k|
2} (25)
In (25) formula, B is a constant relevant with signal constellation (in digital modulation) and channel conditions.
In practice, in order to simplify, B can directly be defined as a constant c.At this moment, because X no longer is a factor of matrix computations, so P
0 HP
0+ (c/SNR) the contrary of I need not recomputate along with the variation of transmission data X.In addition, because P
0 HP
0Known in advance, SNR can survey in receiving end, at this moment
Be the function of SNR, in real system,
Can directly obtain by look-up table.So formula (24) just is reduced to a multiplication of matrices computing, has:
Wherein:
At this moment, amount of calculation is all littler than LS method.Because the LS method relates to matrix inversion operation, and Q directly obtains by look-up table in (26) formula, does not have computation complexity; (26) calculating of formula just is mainly reflected in the multiplying of matrix Q and matrix z, and this operand also is indispensable for the LS method, so the LMMSEE method computation complexity after the simplification is all lower than LS method.Certainly, comparing with the amount of calculation of MMSEE method, then is to have significantly reduced.
The present invention is applicable to and adopts OFDM to carry out communication system for modulation, especially is applicable to European digital terrestrial television (DVB-T) system.Though technology of the present invention is primarily aimed at the system of frequency division multiple access, but any engineer with knowledge background such as signal processing, communications, can design corresponding channel estimating apparatus at code division multiple access and time division multiple access according to the present invention, it all should be included in thought of the present invention and the scope.
In the present invention,, can adopt the mode of look-up table to realize, so the present invention compares with existing method has advantage on computation complexity because do not need to calculate inverse matrix.
Simultaneously, prove also that through the ofdm communication system link simulation compare with other conventional methods, the present invention can accurately estimate channel parameter values.
Description of drawings:
Fig. 1 is at the stack PN of different-energy sequence LS method and the art of this patent LMMSEE method channel estimating performance comparison diagram.Fig. 1 is N
c=252 o'clock, LS method, LMMSEE method channel estimating performance were relatively.Wherein, the energy factors that overlying training sequence distributed in LMMSEE method and the LS method is fixed as β=0.3 and 0.5, c=1.0.As seen from the figure, LMMSEE method performance when low signal-to-noise ratio is better than the LS method slightly, and when signal to noise ratio was higher, the LS method was more excellent.
Fig. 2 is at stack PN sequence energy fixedly the time, LS method and the art of this patent LMMSEE method channel estimating performance comparison diagram.Fig. 2 is LS method and LMMSEE method channel estimating MSE comparison diagram, β=0.5 wherein, c=1.0.As seen, two kinds of method performances are close.
Fig. 3 is during to stack PN sequence allocation different-energy, the art of this patent LMMSEE method channel estimating performance comparison diagram.
Wherein, N
c=252, B=1.As seen from the figure, be lower than 5dB and β=0.3 o'clock in signal to noise ratio, systematic function and channel estimating performance and ideal situation are approaching.
Embodiment:
Below by concrete enforcement technical scheme of the present invention is further described.
Apply the present invention to ofdm system.This system emulation major parameter is: channel model h (l)=0.6912+2.5340i, 1.0603+0.4829i ,-0.0014-0.5519i ,-0.6523-0.0743i, 0,0.2098+0.1348i ,-0.0287-0.0617i; Pilot tone p (n)=σ
pExp{j π/7}exp{-j π (n
2+ 1)/and 7}, cycle P=7, its DFT are P (k)=2 σ
pExp{j π (k
2+ 1)/7} is a chirp.Our fixed transmission power
The 4QAM modulation.
Concrete steps are:
1, transmitting terminal is sent into the OFDM base band signal modulated.
2, selection cycle is that the PN sequence of P is carried out channel estimating.In order to estimate the channel that length is L, require P 〉=L.The PN sequence of the weak energy of stack on the unique therein information symbol.
3, be N with each information symbol tail length
gData be replicated in the information symbol front at interval as protection, through handling, generation transmits.
4, at receiving terminal, received signal by A/D and low pass filter after, protection is at interval deleted, obtains the received signal matrix:
y=Ph+u
5, suppose
It is an integer.At this moment, size of structure is P * N
cMean matrix
I wherein
PIt is the unit matrix of a P * P.Calculate matrix z, wherein z=Ay.
6, establish P
0=AP calculates P
0
7, actual measurement receiving end signal to noise ratio snr is set invariant c, obtains matrix Q by look-up table.
Claims (7)
1, the present invention relates to a kind of method and apparatus that utilizes weak energy parallel PN sequence to carry out channel estimating, it is characterized in that comprising the steps:
Step 1 transmitting terminal is sent into the OFDM base band signal modulated.
Step 2 selection cycle is that the PN sequence of P is carried out channel estimating.In order to estimate the channel that length is L, require P 〉=L.The PN sequence of the weak energy of stack on the unique therein information symbol.
Step 3 is N with each information symbol tail length.Data be replicated in the information symbol front at interval as protection, through handling, generation transmits.
Step 4 is at receiving terminal, received signal by A/D and low pass filter after, protection is deleted at interval, obtains the received signal matrix:
y=Ph+u
Step 5 hypothesis
It is an integer.At this moment, size of structure is P * N
cMean matrix
I wherein
PIt is the unit matrix of a P * P.Calculate matrix z, wherein:
z=Ay
Step 6 is established P
0=AP calculates P
0
Step 7 actual measurement receiving end signal to noise ratio snr is set invariant c, obtains matrix Q by look-up table.
2, said as claim 1, a kind of method of utilizing weak energy parallel PN sequence to carry out channel estimating is characterized in that: c=1.
3, the method for utilizing weak energy parallel PN sequence to carry out channel estimating as claimed in claim 1 is characterized in that: under the situation of having determined training sequence, and matrix
Change along with the signal to noise ratio of actual measurement, so can obtain matrix Q by signal to noise ratio by the mode of look-up table; While z=Ay, matrix A just becomes a constant matrices thereupon after having determined training sequence, so matrix z does linear change with received signal matrix y, at this moment, also can obtain matrix z by the mode of look-up table; According to
Just can be by the very simple impulse response estimated value that obtains channel of calculating.
4, the device that utilizes weak energy parallel PN sequence to carry out channel estimating as claimed in claim 1 is characterized in that: utilize the impulse response estimated value of being determined channel by the training sequence that is adopted, the received signal y that measures and channel signal to noise ratio snr
5, as the right 4 described devices that utilize weak energy parallel PN sequence to carry out channel estimating, it is characterized in that: after training sequence is determined, matrix Q and Z can calculate according to selected training sequence in advance, be stored in the memory device, after measuring received signal y, signal to noise ratio snr, obtain Q and Z by look-up table, thereby calculate the impulse response estimated value of channel
6, said as claim 1, a kind of method and apparatus that utilizes weak energy parallel PN sequence to carry out channel estimating is characterized in that: can be used for employing OFDM such as single carrier, multicarrier and protocol of wireless local area network and DTB Digital Terrestrial Broadcasting TV agreement and carry out channel estimating in the communication system for modulation.
7, as described in the right 1, a kind of method and apparatus that utilizes weak energy parallel PN sequence to carry out channel estimating is characterized in that: the present invention is applicable to and adopts OFDM to carry out communication system for modulation, especially is applicable to European digital terrestrial television (DVB-T) system.Though technology of the present invention is primarily aimed at the system of frequency division multiple access, but any engineer with knowledge background such as signal processing, communications, can design corresponding channel estimating apparatus at code division multiple access and time division multiple access according to the present invention, these all should be included in thought of the present invention and the scope.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB200510035892XA CN100493052C (en) | 2005-07-20 | 2005-07-20 | Method for channel evaluating of weak energy parallel PN sequence |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB200510035892XA CN100493052C (en) | 2005-07-20 | 2005-07-20 | Method for channel evaluating of weak energy parallel PN sequence |
Publications (2)
Publication Number | Publication Date |
---|---|
CN1901518A true CN1901518A (en) | 2007-01-24 |
CN100493052C CN100493052C (en) | 2009-05-27 |
Family
ID=37657261
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CNB200510035892XA Expired - Fee Related CN100493052C (en) | 2005-07-20 | 2005-07-20 | Method for channel evaluating of weak energy parallel PN sequence |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN100493052C (en) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101883065A (en) * | 2009-05-09 | 2010-11-10 | 电子科技大学中山学院 | Nonlinear companding method and device for superimposed training sequence |
CN101425987B (en) * | 2007-10-30 | 2011-05-04 | 华为技术有限公司 | Channel estimation method and apparatus |
CN101547019B (en) * | 2008-03-25 | 2012-10-03 | 卓胜微电子(上海)有限公司 | Method for evaluating channel in DTMB system and device realizing same |
-
2005
- 2005-07-20 CN CNB200510035892XA patent/CN100493052C/en not_active Expired - Fee Related
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101425987B (en) * | 2007-10-30 | 2011-05-04 | 华为技术有限公司 | Channel estimation method and apparatus |
CN101547019B (en) * | 2008-03-25 | 2012-10-03 | 卓胜微电子(上海)有限公司 | Method for evaluating channel in DTMB system and device realizing same |
CN101883065A (en) * | 2009-05-09 | 2010-11-10 | 电子科技大学中山学院 | Nonlinear companding method and device for superimposed training sequence |
Also Published As
Publication number | Publication date |
---|---|
CN100493052C (en) | 2009-05-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP4419957B2 (en) | WIRELESS COMMUNICATION SYSTEM, RECEIVING APPARATUS AND RECEIVING METHOD FOR Performing Multicarrier Transmission | |
TWI241810B (en) | Channel estimation in OFDM systems | |
JP4046515B2 (en) | Iterative maximum likelihood channel estimation and signal detection system and method for OFDM systems | |
CN108924943B (en) | Maximum correlation estimation detection method based on narrowband Internet of things random access channel | |
CN1642159A (en) | Apparatus and method for estimating interference and noise in a communication system | |
CN102365833B (en) | For the channel estimating of control channel in ofdm system | |
CN104410590A (en) | Short-wave OFDM (Orthogonal Frequency Division Multiplexing) interference suppression joint channel estimation method based on compressed sensing | |
CN1890935A (en) | Channel estimation for OFDM systems | |
CN101141425A (en) | Time-division pilot based channel estimation method of mobile communication system | |
CN104735014B (en) | A kind of time synchronization method related based on leading symbol difference | |
CN1756248A (en) | MIMO OFDM mobile communication system and channel estimating method | |
US20050190800A1 (en) | Method and apparatus for estimating noise power per subcarrier in a multicarrier system | |
CN1703034A (en) | A MIMO-OFDM system based channel estimation method | |
CN101309243A (en) | Novel OFDM parameterized channel estimator | |
CN1753395A (en) | Symbol timing method for multi-antenna wireless communication system | |
CN107124194B (en) | Marine digital intermediate frequency emits equipment | |
Gupta et al. | Channel estimation for wavelet based OFDM system | |
CN1909527A (en) | Frequency domain channel estimation method of crossing frequency division multiplexing system with time-domain enveloping weighting | |
CN101635598A (en) | Method and device for estimating noise power | |
CN101252560A (en) | High-performance OFDM frame synchronization algorithm | |
EP1994701A1 (en) | Identification of maximum cyclic deley in an ofdm system based on coherence bandwith of the channel | |
CN1901518A (en) | Method for channel evaluating of weak energy parallel PN sequence | |
US9369329B1 (en) | Low-complexity non-data-aided estimation of symbol time offset in OFDM systems | |
CN102487364B (en) | Channel estimation method and apparatus thereof | |
US8731109B2 (en) | Methods and systems for effective channel estimation in OFDM systems |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
C17 | Cessation of patent right | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20090527 Termination date: 20100720 |