CN101577692A - Channel estimating method of orthogonal frequency division multiplexing system and device thereof - Google Patents

Channel estimating method of orthogonal frequency division multiplexing system and device thereof Download PDF

Info

Publication number
CN101577692A
CN101577692A CNA2009100872659A CN200910087265A CN101577692A CN 101577692 A CN101577692 A CN 101577692A CN A2009100872659 A CNA2009100872659 A CN A2009100872659A CN 200910087265 A CN200910087265 A CN 200910087265A CN 101577692 A CN101577692 A CN 101577692A
Authority
CN
China
Prior art keywords
time domain
impulse response
domain impulse
channel
effective
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CNA2009100872659A
Other languages
Chinese (zh)
Other versions
CN101577692B (en
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.)
Beijing University of Posts and Telecommunications
Original Assignee
Beijing University of Posts and Telecommunications
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 Beijing University of Posts and Telecommunications filed Critical Beijing University of Posts and Telecommunications
Priority to CN2009100872659A priority Critical patent/CN101577692B/en
Publication of CN101577692A publication Critical patent/CN101577692A/en
Application granted granted Critical
Publication of CN101577692B publication Critical patent/CN101577692B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Radio Transmission System (AREA)
  • Mobile Radio Communication Systems (AREA)

Abstract

The invention discloses a channel estimating method of an orthogonal frequency division multiplexing system and a device thereof. The method comprises the following steps of: obtaining initial channel time domain impulse response; obtaining channel time domain impulse response corresponding to an effective multidiameter delay position in the initial channel time domain impulse response and generating an interference matrix according to the effective multidiameter delay position; multiplying the channel time domain impulse response corresponding to the multidiameter delay position by the inverse of the interference matrix to obtain the channel time domain impulse response after modification; and carrying out Fourier transformation to the channel time domain impulse response after modification to obtain channel frequency domain response estimation. In the method and the device, by constructing the multidiameter interference, the inverse of the interference matrix is multiplied by the effective channel time domain impulse response so as to restrict the energy leakage, therefore, the problem of the energy leakage caused by the incoming of virtual subcarrier as well as the floor effect can be effectively eliminated, the estimation precision of the channel is improved, and the system complexity can be controlled by rapid matrix algorithm.

Description

A kind of channel estimation methods of ofdm system and device
Technical field
The present invention relates to the wireless mobile telecommunication technology field, relate in particular to a kind of channel estimation methods and device of ofdm system.
Background technology
OFDM (OFDM) technology is multiplexing on the different orthogonal subcarrier after with the high-speed data-flow serial to parallel conversion, has improved spectrum efficiency when disturbing between data reducing.And the subchannel of each subcarrier correspondence is the frequency-flat fading channel, has effectively reduced the intersymbol interference that multipath fading causes.At present, the OFDM technology is able to extensive use in WLAN (wireless local area network), enhancement mode 3-G (Generation Three mobile communication system) etc.
For guaranteeing the superperformance of ofdm system in the wireless mobile communications environment, receiving terminal need provide to be estimated as far as possible accurately to wireless fading channel.Therefore, channel estimating is the important step that can not lack in the ofdm system frenulum signal processing.
Fig. 1 illustrates the system schematic that data send and receive in the ofdm system.As shown in Figure 1, at transmitting terminal, data are through chnnel coding, modulation, Space Time Coding, and the signal of generation and pilot tone signal carry out the resource mapping together, launch after handling through IFFT variation, framing again.At receiving terminal, after separating through FFT conversion, pilot data, the pilot signal of separation is used for channel estimating, deciphers when according to the result of channel estimating isolated data-signal being carried out sky, the data that send through operation reduction such as demodulation, channel decodings again.
In the ofdm system of reality, for avoiding sending the influence that signal is subjected to low pass filter, the parton carrier wave of reserving band edge that system distributes usually is as the virtual subnet carrier wave.Because the pilot frequency sequence that falls in the virtual carrier scope can't provide this frequency domain scope interior channel information, if the channel frequency domain response estimation value of lost part frequency domain information is carried out inverse fourier transform, the time domain impulse response that obtains will be subjected to the influence of energy leakage and the phase mutual interference causes the decline of precision of channel estimation.Existing channel estimation technique is not considered the influence of energy leakage mostly, in time domain limited sampled point is intercepted, and has lost effective time domain impulse response, floor effect occurs.Therefore, demand a kind of method of effective elimination energy leakage influence urgently, improve the precision of channel estimation in the real system.
Summary of the invention
The technical problem that the present invention will solve provides a kind of channel estimation methods of ofdm system, the precision of channel estimation in can the raising system.
The invention provides a kind of channel estimation methods of ofdm system, comprising: the patterns available that is positioned at effective subcarrier scope is carried out channel frequency domain response estimation, carry out inverse Fourier transform, obtain initial channel time domain impulse response; With the corresponding channel time domain impulse response in effective multidiameter delay position, generate interference matrix in the acquisition initial channel time domain impulse response according to effective multidiameter delay position; Contrary and the channel time domain impulse response corresponding with effective multidiameter delay position of interference matrix multiplied each other, obtain revised channel time domain impulse response; Revised channel time domain impulse response is carried out Fourier transform, obtain channel frequency domain response estimation.
According to an embodiment of channel estimation methods of the present invention, this method is also determined described effective multidiameter delay position according to default effective diameter decision threshold.Wherein, this default effective diameter decision threshold Thr is:
Thr = ρP 2 N av / ξ
Wherein, ρ=J 0(2 π f DΔ mT s), T sBe OFDM symbol period, f DBe maximum doppler frequency, ξ is the signal-to-noise ratio (SNR) estimation value, and P is a frequency pilot sign power, N AvBe the OFDM symbolic number that to add up in order to calculating decision threshold Thr.
According to an embodiment of channel estimation methods of the present invention, generate interference matrix according to described effective multidiameter delay position by following formula:
Figure A20091008726500072
Wherein δ i , j ( l p , l q ) = W N β l q , l p sin ( π ( l q - l p + α i ) N e / N ) sin ( π ( l q - l p + α i ) D f / N ) , W N β l q , l p = exp ( - j 2 π β l q , l p / N ) , l 0, l 1..., l L '-1Be described effective multidiameter delay position, Be the phase rotation coefficient of different transmitting antenna time domain impulse responses, α iBe the shift factor of different transmitting antenna time domain impulse responses, N is an IFFT length, N eBe effective sub-carrier number, D fBe the frequency pilot sign frequency domain interval.
Further, for the ofdm system of frequency division multi-antenna pilot, β l p , l q = l q i - ( l q - l p ) D f / 2 , α i=0; Perhaps, divide the ofdm system of multi-antenna pilot for sign indicating number, β l p , l q = ( l q - l p + α i ) / 2 , α i=(i-1) K 0, D f=1, K wherein 0The spacing of time domain impulse response on time shaft for each transmitting antenna.
According to an embodiment of channel estimation methods of the present invention, the step of revised channel time domain impulse response being carried out Fourier transform acquisition channel frequency domain response estimation comprises: revised channel time domain impulse response is mended 0 make number of samples extend to the FFT/IFFT of system length N; Channel time domain impulse response sequence after the expansion is carried out the FFT conversion, obtain the channel frequency domain response on each subcarrier.
Another technical problem that the present invention will solve provides a kind of channel estimating apparatus of ofdm system, the precision of channel estimation in can the raising system.
The invention provides a kind of channel estimating apparatus of ofdm system, comprise: initial time domain impulse response acquisition module, be used for the patterns available that is positioned at effective subcarrier scope is carried out channel frequency domain response estimation, described channel frequency domain response estimation is carried out inverse Fourier transform obtain initial channel time domain impulse response, send described initial channel time domain impulse response; Effective time domain impulse response acquisition module, be used to receive initial channel time domain impulse response from described initial time domain impulse response acquisition module, obtain in the described initial channel time domain impulse response and the effective corresponding channel time domain impulse response in multidiameter delay position, send and the corresponding channel time domain impulse response in described effective multidiameter delay position; The interference matrix acquisition module is used for generating interference matrix according to described effective multidiameter delay position, sends described interference matrix; Time domain impulse response correcting module, be used to receive from described effective time domain impulse response acquisition module with the corresponding channel time domain impulse response in described effective multidiameter delay position, reception is from the described interference matrix of described interference matrix acquisition module, contrary and the channel time domain impulse response corresponding with described multidiameter delay position of described interference matrix multiplied each other, obtain revised channel time domain impulse response, send described revised channel time domain impulse response; The frequency domain response estimation acquisition module, be used to receive described revised channel time domain impulse response from described time domain impulse response correcting module, described revised channel time domain impulse response is carried out Fourier transform, obtain revised channel frequency domain response estimation.
Further, this channel estimating apparatus also comprises: effective multidiameter delay position acquisition module, be used for determining described effective multidiameter delay position according to default effective diameter decision threshold, and send described effective multidiameter delay position; Described effective time domain impulse response acquisition module also is used to receive the described effective multidiameter delay position from described effective multidiameter delay position acquisition module; Described interference matrix acquisition module also is used to receive the described effective multidiameter delay position from described effective multidiameter delay position acquisition module.
Further, above-mentioned frequency domain response estimation acquisition module makes number of samples extend to the FFT/IFFT of system length N to described revised channel time domain impulse response benefit 0; Channel time domain impulse response sequence after the described expansion is carried out the FFT conversion, obtain the revised channel frequency domain response on each subcarrier.
Channel estimation methods provided by the invention and device by the interference of structure multipath, multiply by the inhibition energy leakage with the contrary of interference matrix with efficient channel time domain impulse response mutually, eliminate floor effect, improve channel estimated accuracy.
Description of drawings
Fig. 1 illustrates the system schematic that data send and receive in the ofdm system;
Fig. 2 illustrates the flow chart of an embodiment of the channel estimation methods of ofdm system of the present invention;
Fig. 3 illustrates the flow chart of another embodiment of the channel estimation methods of ofdm system of the present invention;
Fig. 4 illustrates the block diagram of an embodiment of the channel estimating apparatus of ofdm system of the present invention;
Fig. 5 illustrates the block diagram of another embodiment of the channel estimating apparatus of ofdm system of the present invention.
Embodiment
With reference to the accompanying drawings the present invention is described more fully, exemplary embodiment of the present invention wherein is described.In the accompanying drawings, identical label is represented identical or similar assembly or element.
Fig. 2 illustrates the flow chart of an embodiment of the channel estimation methods of ofdm system of the present invention.
As shown in Figure 2,, the patterns available that is positioned at effective subcarrier scope is carried out channel frequency domain response estimation, carry out inverse Fourier transform then, obtain initial channel time domain impulse response in step 202.For example, the channel frequency domain is estimated and can be carried out based on least square or minimum mean square error criterion by carrying out on each pilot sub-carrier.
In step 204, obtain in the initial channel time domain impulse response and the effective corresponding channel time domain impulse response in multidiameter delay position, and according to effective multidiameter delay position generation interference matrix.Effectively multi-path location can be by determining according to predefined effective diameter decision threshold; For the known system of channel multi-path time delay, can directly determine effective multidiameter delay position according to known multidiameter delay information.For the system of channel multi-path time delay the unknown, the method that keeps multipath comprises also preceding Lmax point is all remained that wherein Lmax is a maximum multipath time delay; Perhaps, by amplitude from big to small, keep the individual point of L ', L ' is an empirical value, and is different in the different channels environment.
In step 206, the contrary and channel time domain impulse response corresponding with effective multidiameter delay position of interference matrix multiplied each other, obtain revised channel time domain impulse response;
In step 208, revised channel time domain impulse response is carried out Fourier transform, obtain revised channel frequency domain response estimation.
Fig. 3 illustrates the flow chart of another embodiment of the channel estimation methods of ofdm system of the present invention.
As shown in Figure 3, in step 302, received signal is carried out the FFT conversion with the signal that receives, and will in the received signal through the FFT conversion frequency pilot sign be separated.
In step 304, try to achieve channel frequency domain response estimation according to least square or minimum mean square error criterion, behind IFFT, obtain initial time domain impulse response estimated value
Figure A20091008726500101
I, j represent i root transmitting antenna and j root reception antenna respectively.
In step 306,, seek effective multidiameter delay position l according to predefined effective diameter decision threshold 0, l 1..., l L '-1, keep and the corresponding efficient channel time domain impulse response in effective multidiameter delay position h ~ i , j = [ h ^ i , j ( l 0 ) , h ^ i , j ( l 1 ) , . . . h ^ i , j ( l L ′ - 1 ) ] T .
In step 308, according to the interference δ of each multidiameter delay correspondence of effective multidiameter delay position reconstruct of determining I, j, and structure interference matrix A I, j:
Wherein δ i , j ( l p , l q ) = W N β l q , l p sin ( π ( l q - l p + α i ) N e / N ) sin ( π ( l q - l p + α i ) D f / N ) , W N β l q , l p = exp ( - j 2 π β l q , l p / N ) ,
Figure A20091008726500115
Be the phase rotation coefficient of different transmitting antenna time domain impulse responses, α iShift factor for different transmitting antenna time domain impulse responses.According to different pilot tone mode (for example, frequency division multi-antenna pilot, sign indicating number divide multi-antenna pilot), these two parameters are different.
Use interference matrix A I, jThe contrary efficient channel time domain impulse response that multiply by h ~ i , j = [ h ^ i , j ( l 0 ) , h ^ i , j ( l 1 ) , . . . h ^ i , j ( l L ′ - 1 ) ] T , The time domain impulse response h that obtains revising I, j:
h ‾ i , j = A i , j - 1 h ~ i , j - - - ( 2 )
In step 310, with the time domain impulse response h that revises I, jAdd 0 to N point, carry out FFT conversion, the channel frequency domain response estimation that obtains revising then.
According to one embodiment of present invention,, can obtain predefined effective decision threshold Thr according to following formula for the system of channel multi-path time delay the unknown:
Thr = ρP 2 N av / ξ - - - ( 3 )
Wherein, ρ=J 0(2 π f DΔ mT s), T sBe OFDM symbol period, f DBe maximum doppler frequency, ξ is the signal-to-noise ratio (SNR) estimation value, and P is a frequency pilot sign power, N AvBe the OFDM symbolic number that to add up in order to calculating decision threshold Thr.According to formula (3) initial channel time domain impulse response being carried out effective diameter keeps:
h ~ i , j ( n ) = h ^ i , j ( n ) , | h ^ i , j ( n ) | 2 ≥ Thr 0 , else - - - ( 4 )
According to one embodiment of present invention, for the system of channel multi-path time delay the unknown,, seek effective multidiameter delay position according to predefined decision threshold.And a plurality of dual-mode antennas can be united statistics to reaching a plurality of frequency pilot signs, to determine effective multidiameter delay position l 0, l 1..., l L '-1
According to one embodiment of present invention, for channel multi-path time delay system, can determine effective multidiameter delay position l directly according to known multidiameter delay information 0, l 1..., l L '-1, and the time domain impulse response of effective multidiameter delay position correspondence kept:
h ~ i , j = [ h ^ i , j ( l 0 ) , h ^ i , j ( l 1 ) , . . . h ^ i , j ( l L ′ - 1 ) ] T .
According to one embodiment of present invention, can be with interference matrix A I, jBe modified to the complex conjugate symmetrical matrix, find the solution more precise channels time domain impulse response, thereby reach the effect that reduces computation complexity with Qiao Lisiji (Cholesky) decomposition method.
It may be noted that method of the present invention is applicable to that the single-shot list is received, single-shot is overcharged or the ofdm communication system of MIMO.
Below by a concrete application examples method of the present invention is further detailed.
In this application examples, the FFT/IFFT length of supposing ofdm system is N, and effectively sub-carrier number is N e, the virtual subnet carrier number is N-N eThe frequency-region signal Y that receives on the receiver j root antenna j[k] is:
Y j[k]=H j,i[k]P i[k]+V j[k],k∈Ω uc (5)
Wherein, H J, iThe channel frequency domain response of [k] expression i root transmitting antenna to the j root reception antenna, P i[k] is the frequency pilot sign of i root transmitting antenna, Ω UcBe effectively subcarrier set.
The first step: on each pilot sub-carrier, carry out estimating, and carry out the IFFT conversion and obtain initial channel time domain impulse response based on the channel frequency domain of least square or minimum mean square error criterion
Figure A20091008726500123
For the ofdm system of frequency division multi-antenna pilot method, promptly on the subcarrier of each transmitting antenna correspondence divided by separately pilot frequency sequence, obtain the noisy channel frequency domain response estimation value of each transmitting antenna:
H ^ j , i [ k ] = Y j [ k ] / P i [ k ] = H j , i [ k ] + V j [ k ] / P i [ k ] , k ∈ Ω uc - - - ( 6 )
To each transmitting antenna Carry out N/D respectively fPoint IFFT conversion obtains the initial channel time domain impulse response of each transmitting antenna
Figure A20091008726500133
Wherein N is an IFFT length, D fBe the frequency pilot sign frequency domain interval.
Dividing the ofdm system of multi-antenna pilot method for sign indicating number, is example with the exponential type pilot tone, and the conjugation that is about to first transmitting antenna frequency pilot sign multiply by frequency-region signal Y j[k] obtains the adding up of noisy channel frequency domain response estimation value of each transmitting antenna:
H ^ [ k ] = Y j [ k ] P i * [ k ] / | P i [ k ] | 2 = Σ i = 1 T x H j , i [ k ] e - j 2 πk ( i - 1 ) K 0 / N + V ^ j [ k ] , k ∈ Ω uc - - - ( 7 )
To frequency-region signal
Figure A20091008726500135
Carry out N point IFFT conversion, obtain the initial channel time domain impulse response of all transmitting antennas
Figure A20091008726500136
Time domain impulse response non-overlapping copies on time shaft of each transmitting antenna, its spacing is K 0
Second step: set effective diameter decision threshold Thr according to above-mentioned formula (3), according to the initial channel time domain impulse response of above-mentioned formula (4) from first step gained Middle effective multidiameter delay position l that seeks 0, l 1..., l L '-1, and keep the time domain impulse response of described effective multidiameter delay position correspondence: h ~ i , j = [ h ^ i , j ( l 0 ) , h ^ i , j ( l 1 ) , . . . h ^ i , j ( l L ′ - 1 ) ] T . For improving the accuracy that multipath detects, can determine effective multidiameter delay position by the statistics of uniting to a plurality of transmitting antennas and a plurality of frequency pilot signs.
The 3rd step: according to the interference of definite each multidiameter delay correspondence of effective multidiameter delay position reconstruct of above-mentioned steps δ i , j ( l p , l q ) = W N β l q , l p sin ( π ( l q - l p + α i ) N e / N ) / sin ( π ( l q - l p + α i ) D f / N ) .
For different pilot tone modes, disturb δ I, jDifferent.For example, for the ofdm system of frequency division multi-antenna pilot method, β l p , l q = l q i - ( l q - l p ) D f / 2 , α i=0。Divide the ofdm system of multi-antenna pilot method for sign indicating number, β l p , l q = ( l q - l p + α i ) / 2 , α i=(i-1)K 0,D f=1。
Then according to δ I, jFormula (1) structure interference matrix A by the front I, j
With the time domain impulse response that keeps in second step
Figure A200910087265001312
With interference matrix A I, jAnti-phase take advantage of (referring to top formula (2)), to eliminate energy leakage that the virtual subnet carrier wave introduces, obtain channel time domain impulse response h more accurately to the influence that channel estimating causes I, j
Preferably, A I, jBe deformed into the complex conjugate symmetrical matrix.In frequency division multi-antenna pilot mode is example, will
Figure A20091008726500141
With
Figure A20091008726500142
Merge new interference matrix
Figure A20091008726500143
(p, q) individual element is δ ~ i ( l p , l q ) = W N ( ( l q - l p ) D f / 2 ) sin ( π ( l q - l p ) N e / N ) / sin ( π ( l q - l p ) D f / N ) , This interference matrix
Figure A20091008726500145
Be the complex conjugate symmetrical matrix.Can utilize Cholesky to decompose the solving equation group and replace matrix inversion, thereby reach the effect that reduces computation complexity.
The 4th step: the 3rd step gained channel time domain impulse response is estimated h ‾ i , j = [ h ‾ i , j ( l 0 ) , h ^ i , j ( l 1 ) , . . . , h ^ i , j ( l L ′ - 1 ) ] T Carry out Fourier transform after adding 0 to N point, obtain revised channel frequency domain response estimation.
Fig. 4 illustrates the block diagram of an embodiment of the channel estimating apparatus of ofdm system of the present invention.As shown in Figure 4, this channel estimating apparatus comprises initial time domain impulse response acquisition module 41, effective time domain impulse response acquisition module 42, interference matrix acquisition module 43, time domain impulse response correcting module 44 and frequency domain response estimation acquisition module 45.Wherein, initial time domain impulse response acquisition module 41 is used for the patterns available that is positioned at effective subcarrier scope is carried out channel frequency domain response estimation, channel frequency domain response estimation is carried out inverse Fourier transform, obtain initial channel time domain impulse response, send the initial channel time domain impulse response that obtains and arrive effective time domain impulse response acquisition module 42.Effective time domain impulse response acquisition module 42, be used to receive initial channel time domain impulse response from initial time domain impulse response acquisition module 41, with the corresponding channel time domain impulse response in effective multidiameter delay position, send and the corresponding channel time domain impulse response in effective multidiameter delay position in the acquisition initial channel time domain impulse response.Interference matrix acquisition module 43 is used for generating interference matrix according to effective multidiameter delay position, sends interference matrix.Time domain impulse response correcting module 44, be used to receive the channel time domain impulse response corresponding with effective multidiameter delay position from effective time domain impulse response acquisition module 42, reception is from the interference matrix of interference matrix acquisition module 43, contrary and the channel time domain impulse response corresponding with effective multidiameter delay position of interference matrix multiplied each other, obtain revised channel time domain impulse response, send revised channel time domain impulse response to frequency domain response estimation acquisition module 45.Frequency domain response estimation acquisition module 45 is used to receive the revised channel time domain impulse response from time domain impulse response correcting module 44, and revised channel time domain impulse response is carried out Fourier transform, obtains revised channel frequency domain response estimation.
Fig. 5 illustrates the block diagram of another embodiment of the channel estimating apparatus of ofdm system of the present invention.The channel estimating apparatus of Fig. 5 is compared with Fig. 4, increased effective multidiameter delay position acquisition module 56, be used for determining effective multidiameter delay position, send effective multidiameter delay position to effective time domain impulse response acquisition module 42 and interference matrix acquisition module 43 according to default effective diameter decision threshold.Effectively time domain impulse response acquisition module 43 receives effective multidiameter delay position from effective multidiameter delay position acquisition module 56.Interference matrix acquisition module 43 is from effective multidiameter delay position of effective multidiameter delay position acquisition module 56.Other module is described can be referring to the description of respective modules among Fig. 4, for for purpose of brevity, this no longer auspicious stating.
Wherein, the default effective diameter decision threshold Thr that effective multidiameter delay position acquisition module 56 adopts among Fig. 5 is Thr = ρP 2 N av / ξ . Wherein, ρ=J 0(2 π f DΔ mT s), T sBe OFDM symbol period, f DBe maximum doppler frequency, ξ is the signal-to-noise ratio (SNR) estimation value, and P is a frequency pilot sign power, N AvBe the OFDM symbolic number that to add up in order to calculating decision threshold Thr.
According to an embodiment of channel estimating apparatus of the present invention, the effective diameter decision threshold that effective multidiameter delay position acquisition module basis is preset is by adding up definite described effective multidiameter delay position to a plurality of dual-mode antennas to reaching uniting of a plurality of frequency pilot signs.In addition, the interference matrix acquisition module generates interference matrix A by top formula (1) according to effective multidiameter delay position I, j
According to an embodiment of channel estimating apparatus of the present invention, time domain impulse response correcting module is modified to the complex conjugate symmetrical matrix with interference matrix, calculates the contrary of interference matrix by the fast algorithm that adopts matrix inversion, to reduce computation complexity.
According to an embodiment of channel estimating apparatus of the present invention, the frequency domain response estimation acquisition module mends 0 to revised channel time domain impulse response makes number of samples extend to the FFT/IFFT of system length N; Channel time domain impulse response sequence after the expansion is carried out the FFT conversion, obtain the revised channel frequency domain response on each subcarrier.
Method and apparatus of the present invention is by the interference of structure multipath, the contrary of interference matrix be multiply by the inhibition energy leakage mutually with efficient channel time domain impulse response, be applicable to single antenna and many antennas ofdm system, can be applied to flexibly or the system of Unknown Channel multidiameter delay information.The method can effectively be eliminated the Energy Leaking that the virtual subnet carrier wave is introduced, and eliminates floor effect, improves precision of channel estimation, and can be by quick matrix algorithm control system complexity.
Description of the invention provides for example with for the purpose of describing, and is not exhaustively or limit the invention to disclosed form.Many modifications and variations are obvious for the ordinary skill in the art.Selecting and describing embodiment is for better explanation principle of the present invention and practical application, thereby and makes those of ordinary skill in the art can understand the various embodiment that have various modifications that the present invention's design is suitable for special-purpose.

Claims (11)

1. the channel estimation methods of an ofdm system is characterized in that, comprising:
The patterns available that is positioned at effective subcarrier scope is carried out channel frequency domain response estimation, carry out inverse Fourier transform, obtain initial channel time domain impulse response;
Obtain in the described initial channel time domain impulse response and the effective corresponding channel time domain impulse response in multidiameter delay position, generate interference matrix according to described effective multidiameter delay position;
Contrary and the channel time domain impulse response corresponding with described effective multidiameter delay position of described interference matrix multiplied each other, obtain revised channel time domain impulse response;
Described revised channel time domain impulse response is carried out Fourier transform, obtain channel frequency domain response estimation.
2. the channel estimation methods of ofdm system according to claim 1 is characterized in that, also comprises step:
Determine described effective multidiameter delay position according to default effective diameter decision threshold.
3. the channel estimation methods of ofdm system according to claim 2 is characterized in that, described default effective diameter decision threshold Thr is:
Thr = ρP 2 N av / ξ
Wherein, ρ=J 0(2 π f DΔ mT s), Ts is the OFDM symbol period, and fD is a maximum doppler frequency, and ξ is the signal-to-noise ratio (SNR) estimation value, and P is a frequency pilot sign power, the OFDM symbolic number of Nav for adding up in order to calculating decision threshold Thr.
4. according to the channel estimation methods of claim 2 or 3 described ofdm systems, it is characterized in that the default effective diameter decision threshold of described basis determines that the step of described effective multidiameter delay position comprises:
According to default effective diameter decision threshold by to a plurality of dual-mode antennas to and the statistics of uniting of a plurality of frequency pilot signs determine described effective multidiameter delay position.
5. the channel estimation methods of ofdm system according to claim 1 is characterized in that, generates interference matrix by following formula according to described effective multidiameter delay position:
Figure A2009100872650003C1
Wherein δ i , j ( l p , l q ) = W N β l q , l p sin ( π ( l q - l p + α i ) N e / N ) sin ( π ( l q - l p + α i ) D f / N ) , W N β l q , l p = exp ( - j 2 π β l q , l p / N ) , l 0, l 1..., l L '-1Be described effective multidiameter delay position,
Figure A2009100872650003C4
Be the phase rotation coefficient of different transmitting antenna time domain impulse responses, α iBe the shift factor of different transmitting antenna time domain impulse responses, N is an IFFT length, N eBe effective sub-carrier number, D fBe the frequency pilot sign frequency domain interval.
6, the channel estimation methods of ofdm system according to claim 1 is characterized in that, describedly described revised channel time domain impulse response is carried out the step that Fourier transform obtains channel frequency domain response estimation comprises:
Described revised channel time domain impulse response is mended 0 makes number of samples extend to the FFT/IFFT of system length N;
Channel time domain impulse response sequence after the described expansion is carried out the FFT conversion, obtain the channel frequency domain response on each subcarrier.
7. the channel estimating apparatus of an ofdm system is characterized in that, comprising:
Initial time domain impulse response acquisition module, be used for the patterns available that is positioned at effective subcarrier scope is carried out channel frequency domain response estimation, described channel frequency domain response estimation is carried out inverse Fourier transform obtain initial channel time domain impulse response, send described initial channel time domain impulse response;
Effective time domain impulse response acquisition module, be used to receive initial channel time domain impulse response from described initial time domain impulse response acquisition module, obtain in the described initial channel time domain impulse response and the effective corresponding channel time domain impulse response in multidiameter delay position, send and the corresponding channel time domain impulse response in described effective multidiameter delay position;
The interference matrix acquisition module is used for generating interference matrix according to described effective multidiameter delay position, sends described interference matrix;
Time domain impulse response correcting module, be used to receive from described effective time domain impulse response acquisition module with the corresponding channel time domain impulse response in described effective multidiameter delay position, reception is from the described interference matrix of described interference matrix acquisition module, contrary and the channel time domain impulse response corresponding with described effective multidiameter delay position of described interference matrix multiplied each other, obtain revised channel time domain impulse response, send described revised channel time domain impulse response;
The frequency domain response estimation acquisition module, be used to receive described revised channel time domain impulse response from described time domain impulse response correcting module, described revised channel time domain impulse response is carried out Fourier transform, obtain revised channel frequency domain response estimation.
8. the channel estimating apparatus of ofdm system according to claim 7 is characterized in that, also comprises:
Effectively multidiameter delay position acquisition module is used for determining described effective multidiameter delay position according to default effective diameter decision threshold, sends described effective multidiameter delay position;
Described effective time domain impulse response acquisition module also is used to receive the described effective multidiameter delay position from described effective multidiameter delay position acquisition module; Described interference matrix acquisition module also is used to receive the described effective multidiameter delay position from described effective multidiameter delay position acquisition module.
9. the channel estimating apparatus of ofdm system according to claim 8, it is characterized in that the effective diameter decision threshold that acquisition module basis in described effective multidiameter delay position is preset is by determining described effective multidiameter delay position to a plurality of dual-mode antennas to the statistics of uniting that reaches a plurality of frequency pilot signs.
10. the channel estimating apparatus of ofdm system according to claim 7 is characterized in that, described interference matrix acquisition module generates interference matrix by following formula according to described effective multidiameter delay position:
Figure A2009100872650004C1
Wherein δ i , j ( l p , l q ) = W N β l q , l p sin ( π ( l q - l p + α i ) N e / N ) sin ( π ( l q - l p + α i ) D f / N ) , W N β l q , l p = exp ( - j 2 π β l q , l p / N ) , l 0, l 1..., l L '-1Be described effective multidiameter delay position, Be the phase rotation coefficient of different transmitting antenna time domain impulse responses, α iBe the shift factor of different transmitting antenna time domain impulse responses, N is an IFFT length, N eBe effective sub-carrier number, D fBe the frequency pilot sign frequency domain interval.
11, the channel estimating apparatus of ofdm system according to claim 7 is characterized in that, described frequency domain response estimation acquisition module mends 0 to described revised channel time domain impulse response makes number of samples extend to the FFT/IFFT of system length N; Channel time domain impulse response sequence after the described expansion is carried out the FFT conversion, obtain the revised channel frequency domain response on each subcarrier.
CN2009100872659A 2009-06-15 2009-06-15 Channel estimating method of orthogonal frequency division multiplexing system and device thereof Expired - Fee Related CN101577692B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN2009100872659A CN101577692B (en) 2009-06-15 2009-06-15 Channel estimating method of orthogonal frequency division multiplexing system and device thereof

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN2009100872659A CN101577692B (en) 2009-06-15 2009-06-15 Channel estimating method of orthogonal frequency division multiplexing system and device thereof

Publications (2)

Publication Number Publication Date
CN101577692A true CN101577692A (en) 2009-11-11
CN101577692B CN101577692B (en) 2012-06-13

Family

ID=41272481

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2009100872659A Expired - Fee Related CN101577692B (en) 2009-06-15 2009-06-15 Channel estimating method of orthogonal frequency division multiplexing system and device thereof

Country Status (1)

Country Link
CN (1) CN101577692B (en)

Cited By (17)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102143098A (en) * 2010-12-03 2011-08-03 无锡物联网产业研究院 Channel estimation method and channel estimation device in orthogonal frequency division multiplexing system
CN102790737A (en) * 2011-05-17 2012-11-21 中兴通讯股份有限公司 Synchronization method and device of system
CN103067134A (en) * 2012-03-14 2013-04-24 王晓安 Signal channel feedback method in orthogonal frequency division multiplexing system
CN103701736A (en) * 2013-12-31 2014-04-02 大唐移动通信设备有限公司 Interference inhibiting method and device in LTE (long term evolution) system
CN103905350A (en) * 2014-04-08 2014-07-02 电子科技大学 S-IFFT-based channel estimation method for OFDM system
WO2015061953A1 (en) * 2013-10-28 2015-05-07 华为技术有限公司 Chanel estimation method and apparatus
CN104935534A (en) * 2015-06-23 2015-09-23 北京北方烽火科技有限公司 Channel estimation method and device based on demodulation reference signal
CN106559362A (en) * 2015-09-24 2017-04-05 联芯科技有限公司 The combined channel and data estimation method and system of fast time variant OFDM channels
CN103813344B (en) * 2012-11-12 2017-06-20 中国移动通信集团浙江有限公司 A kind of interference matrix establishment method and system
CN107302392A (en) * 2016-04-14 2017-10-27 中兴通讯股份有限公司 The method and device that antenna data merges under a kind of time delay channel
CN108599914A (en) * 2018-03-12 2018-09-28 西安电子科技大学 Co-channel full duplex bi-directional relaying transmission method while based on orthogonal frequency division multiplexing
CN108768567A (en) * 2018-06-01 2018-11-06 北京邮电大学 A kind of multipath cluster-dividing method, device, electronic equipment and readable storage medium storing program for executing
CN109150783A (en) * 2017-06-28 2019-01-04 大唐移动通信设备有限公司 A kind of channel estimation methods and device
WO2020216056A1 (en) * 2019-04-26 2020-10-29 中兴通讯股份有限公司 Channel estimation method and apparatus
CN112235223A (en) * 2020-09-11 2021-01-15 北京邮电大学 Anti-interference communication method and electronic equipment
CN113037590A (en) * 2021-03-02 2021-06-25 成都爱瑞无线科技有限公司 Time delay estimation method and device used in communication system
CN115118358A (en) * 2021-03-17 2022-09-27 大唐移动通信设备有限公司 Method, device, equipment and storage medium for estimating angle of interference source

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101064571B (en) * 2006-04-29 2010-09-29 上海贝尔阿尔卡特股份有限公司 Apparatus for enhancing channel evaluation in OFDM receiver and its method
CN101127753B (en) * 2007-09-29 2011-12-14 北京邮电大学 A channel estimation method applicable to multi-carrier system
CN101227440B (en) * 2008-02-14 2010-07-21 北京创毅视讯科技有限公司 Apparatus and method for eliminating zero-frequency interference in channel estimation

Cited By (26)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102143098A (en) * 2010-12-03 2011-08-03 无锡物联网产业研究院 Channel estimation method and channel estimation device in orthogonal frequency division multiplexing system
CN102143098B (en) * 2010-12-03 2014-03-05 无锡物联网产业研究院 Channel estimation method and channel estimation device in orthogonal frequency division multiplexing system
CN102790737A (en) * 2011-05-17 2012-11-21 中兴通讯股份有限公司 Synchronization method and device of system
CN103067134A (en) * 2012-03-14 2013-04-24 王晓安 Signal channel feedback method in orthogonal frequency division multiplexing system
CN103067134B (en) * 2012-03-14 2016-01-20 王晓安 Channel feedback method in ofdm system
CN103813344B (en) * 2012-11-12 2017-06-20 中国移动通信集团浙江有限公司 A kind of interference matrix establishment method and system
WO2015061953A1 (en) * 2013-10-28 2015-05-07 华为技术有限公司 Chanel estimation method and apparatus
CN105264849A (en) * 2013-10-28 2016-01-20 华为技术有限公司 Chanel estimation method and apparatus
CN103701736A (en) * 2013-12-31 2014-04-02 大唐移动通信设备有限公司 Interference inhibiting method and device in LTE (long term evolution) system
CN103701736B (en) * 2013-12-31 2017-08-25 大唐移动通信设备有限公司 Disturbance restraining method and device in LTE system
CN103905350A (en) * 2014-04-08 2014-07-02 电子科技大学 S-IFFT-based channel estimation method for OFDM system
CN103905350B (en) * 2014-04-08 2017-01-04 电子科技大学 Channel estimation methods based on S IFFT in a kind of ofdm system
CN104935534B (en) * 2015-06-23 2018-04-27 北京北方烽火科技有限公司 A kind of channel estimation methods and device based on demodulated reference signal
CN104935534A (en) * 2015-06-23 2015-09-23 北京北方烽火科技有限公司 Channel estimation method and device based on demodulation reference signal
CN106559362A (en) * 2015-09-24 2017-04-05 联芯科技有限公司 The combined channel and data estimation method and system of fast time variant OFDM channels
CN106559362B (en) * 2015-09-24 2019-09-20 联芯科技有限公司 The combined channel sum number of fast time variant OFDM channel method and system according to estimates
CN107302392A (en) * 2016-04-14 2017-10-27 中兴通讯股份有限公司 The method and device that antenna data merges under a kind of time delay channel
CN107302392B (en) * 2016-04-14 2020-11-27 中兴通讯股份有限公司 Method and device for combining antenna data under time delay channel
CN109150783A (en) * 2017-06-28 2019-01-04 大唐移动通信设备有限公司 A kind of channel estimation methods and device
CN108599914A (en) * 2018-03-12 2018-09-28 西安电子科技大学 Co-channel full duplex bi-directional relaying transmission method while based on orthogonal frequency division multiplexing
CN108768567A (en) * 2018-06-01 2018-11-06 北京邮电大学 A kind of multipath cluster-dividing method, device, electronic equipment and readable storage medium storing program for executing
WO2020216056A1 (en) * 2019-04-26 2020-10-29 中兴通讯股份有限公司 Channel estimation method and apparatus
CN112235223A (en) * 2020-09-11 2021-01-15 北京邮电大学 Anti-interference communication method and electronic equipment
CN113037590A (en) * 2021-03-02 2021-06-25 成都爱瑞无线科技有限公司 Time delay estimation method and device used in communication system
CN115118358A (en) * 2021-03-17 2022-09-27 大唐移动通信设备有限公司 Method, device, equipment and storage medium for estimating angle of interference source
CN115118358B (en) * 2021-03-17 2023-10-27 大唐移动通信设备有限公司 Method, device, equipment and storage medium for estimating interference source angle

Also Published As

Publication number Publication date
CN101577692B (en) 2012-06-13

Similar Documents

Publication Publication Date Title
CN101577692B (en) Channel estimating method of orthogonal frequency division multiplexing system and device thereof
CN102687441B (en) Channel estimation for OFDM systems
US8320507B2 (en) Mobile communication system, receiving device, and method
EP1473862B1 (en) Apparatus and method for transmitting training symbol groups in an OFDM communications system using multiple antennas
KR100729726B1 (en) System and Method for Timing Acquisition and Carrier Frequency Offset Estimation in Wireless Communication Based on OFDM
US8588204B2 (en) Efficient channel estimation method using superimposed training for equalization in uplink OFDMA systems
CN101958865A (en) Method for generating demodulating reference signal and device thereof
CN102143101A (en) Mirror-extended frequency domain windowing orthogonal frequency division multiple access channel estimation method
CN101291311B (en) Synchronization implementing method and device for multi-input multi-output orthogonal frequency division multiplexing system
CN102045285A (en) Channel estimation method and device and communication system
CN101958866B (en) Pilot frequency insertion method and module
Tureli et al. Multicarrier synchronization with diversity
Yin et al. Design and performance analysis of AFDM with multiple antennas in doubly selective channels
Kalwar et al. Analysis of effect of carrier frequency offset on performance of LTE uplink
CN104205694A (en) Channel estimation method and receiver
US8223865B2 (en) Method for the blind estimation of OFDM signal parameters by adapted filtering
CN102223336B (en) Wireless communication method and equipment
US7324435B2 (en) Methods and apparatus for receiving OFDM signals
Jhingan et al. Performance analysis of OFDM system augmented with SC diversity combining technique in presence of CFO
Khlifi et al. Hybrid LS-LMMSE channel estimation technique for LTE downlink systems
Yuan et al. Research on channel estimation for OFDM receiver based on IEEE 802.11 a
Tureli et al. Performance analysis of OFDM carrier synchronization with diversity
Zhang et al. A low complexity CFO estimation method for UFMC systems
KR20070113563A (en) Method for channel estimation, channel estimator, mobile station and base station
Prakash et al. A new blind carrier frequency offset estimation scheme for 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
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20120613

Termination date: 20180615

CF01 Termination of patent right due to non-payment of annual fee