CN103391266A - Frequency domain channel response obtaining method and device - Google Patents

Frequency domain channel response obtaining method and device Download PDF

Info

Publication number
CN103391266A
CN103391266A CN2012101478204A CN201210147820A CN103391266A CN 103391266 A CN103391266 A CN 103391266A CN 2012101478204 A CN2012101478204 A CN 2012101478204A CN 201210147820 A CN201210147820 A CN 201210147820A CN 103391266 A CN103391266 A CN 103391266A
Authority
CN
China
Prior art keywords
threshold value
noise amplitude
channel response
amplitude threshold
domain channel
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
CN2012101478204A
Other languages
Chinese (zh)
Other versions
CN103391266B (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.)
China Academy of Telecommunications Technology CATT
Datang Mobile Communications Equipment Co Ltd
Original Assignee
China Academy of Telecommunications Technology CATT
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 China Academy of Telecommunications Technology CATT filed Critical China Academy of Telecommunications Technology CATT
Priority to CN201210147820.4A priority Critical patent/CN103391266B/en
Publication of CN103391266A publication Critical patent/CN103391266A/en
Application granted granted Critical
Publication of CN103391266B publication Critical patent/CN103391266B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Noise Elimination (AREA)

Abstract

The invention discloses a frequency domain channel response obtaining method and device. The method includes calculating a first noise amplitude threshold value g by using time domain channel response; calculating a second noise amplitude threshold value G by using the time domain channel response and the first noise amplitude threshold value g; calculating a third noise amplitude threshold value G' by using the second noise amplitude threshold value G; obtaining noise-restrained frequency domain channel response by using the third noise amplitude threshold value G'. By means of the method and device, channel estimation performance can be improved, and especially channel estimation performance under small bandwidth can be improved effectively.

Description

A kind of acquisition methods of domain channel response and equipment
Technical field
The present invention relates to communication technical field, especially related to a kind of acquisition methods and equipment of domain channel response.
Background technology
Channel estimating refers to the time-frequency position that utilizes known pilot signal transmitted and pilot tone to send, and the data-signal that receives in corresponding time-frequency position, obtains the process of spatial channel information, namely obtains domain channel response.Channel estimation process schematic diagram as shown in Figure 1, with LTE(Long Term Evolution, Long Term Evolution) downlink system is example, and it is S that known descending pilot frequency sends sequence, and the signal that receives is Y, and the spatial frequency domain channel is H+n, and H represents the impulse response of wireless fading channel, n represents white Gaussian noise, Y=(H+n) S, thus can estimate domain channel response
Figure BDA00001631131300011
In order further to improve the precision of channel estimating, and estimate simultaneously noise power-value, can also adopt time domain to get the method for window; Its detailed process comprises: first frequency domain channel is changed to time domain, again according to CP(Cyclic Prefix, Cyclic Prefix) length or prior information define the position with the signal time delay maximum diameter, then take out the noise in time domain footpath that does not contain useful signal power, the calculating noise performance number, finally noise path is removed conversion and return frequency domain, with the channel estimation value after the acquisition noise suppression.
The method that noise window is set according to CP length is: with the time elongatedness of CP length as time delay maximum diameter in system, the position of calculating noise window; Wherein, stipulated the Ts(sampled point of CP under different bandwidth in the LTE system) number N CP, FFT(Fast Fourier Transformation, fast Fourier transform) and points N FFT, sub-carrier number N SC, the time-domain position of time delay maximum diameter can be expressed as:
Figure BDA00001631131300012
If adopt IDFT(Inverse Discrete Fourier Transform, inverse discrete Fourier transform) frequency domain is transformed to time domain, the footpath number of time domain channel response is same OFDM(Orthogonal Frequency Division Multiplexing, OFDM) pilot tone point number N on symbol Pliot, residue footpath number is N Pliot-N τ, consider that the signal power leakage can cause afterbody to comprise the problem of available signal power, actual spendable noise path still less, as is got original noise window { N τ+ 1, N τ+ 2 ..., N Pliot-0.5N τ; If adopt mirror image IDFT that frequency domain is transformed to time domain, the footpath number of mirror image time domain channel response is
Figure BDA00001631131300021
The position in useful signal footpath can symmetrical being present in
Figure BDA00001631131300022
Scope in, desirable original noise window is
{2N τ+1,2N τ+2,…,N pliot,N pliot,+2,N pliot+3,…,2N pliot+1-2N τ}
In realizing process of the present invention, the inventor finds in prior art to exist at least following problem:
Support the various bandwidth allocation plan in the LTE system, wherein minimum bandwidth is configured to 1.4MHz, only take 6 PRB (PhVsical Resource Block, Physical Resource Block), because the pilot number on each PRB is consistent, therefore for little bandwidth, its pilot number is relatively less, the amount of information that can use when carrying out channel estimating is less, thereby can cause noise to estimate the measurement result mistake, and channel estimating performance descends.
Concrete, while according to existing mode, carrying out the calculating of original noise window, window is long very little under little bandwidth, and as under the 1.4M bandwidth, this value is If adopt IDFT that frequency domain is transformed to time domain, the footpath number of time domain channel response is 12, residue footpath number is 6, consider the power leakage problem of afterbody useful signal, actual spendable window is long can be still less, if adopt mirror image IDFT that frequency domain is transformed to time domain, the footpath number of mirror image time domain channel response is 24, and the position in useful signal footpath can symmetrical being present in
Figure BDA00001631131300024
Scope in, remaining noise path is approximately 1, the footpath number is very few, thereby it is inaccurate to cause noise to be estimated.In addition, under the long time delay scene, the power in useful signal footpath can be leaked to footpath, noise window edge, when window internal diameter number is less, larger on the estimated value impact of noise power, thus the noise suppression effect of Signal estimation is exerted an influence.
Summary of the invention
The embodiment of the present invention provides a kind of acquisition methods and equipment of domain channel response, to improve the performance of channel estimating.In order to achieve the above object, the embodiment of the present invention provides following technical scheme:
The embodiment of the present invention provides a kind of acquisition methods of domain channel response, comprising:
Utilize time domain channel response to calculate the first noise amplitude threshold value g;
Utilize described time domain channel response and described the first noise amplitude threshold value g to calculate the second noise amplitude threshold value G;
Utilize described the second noise amplitude threshold value G to calculate the 3rd noise amplitude threshold value G ';
Utilize the domain channel response after described the 3rd noise amplitude threshold value G ' obtains noise suppression.
The embodiment of the present invention provides a kind of equipment that obtains of domain channel response, comprising:
The first computing module, be used for utilizing time domain channel response to calculate the first noise amplitude threshold value g;
The second computing module, be used for utilizing described time domain channel response and described the first noise amplitude threshold value g to calculate the second noise amplitude threshold value G;
The 3rd computing module, be used for utilizing described the second noise amplitude threshold value G to calculate the 3rd noise amplitude threshold value G ';
Obtain module, be used for utilizing the domain channel response after described the 3rd noise amplitude threshold value G ' obtains noise suppression.
Compared with prior art, the embodiment of the present invention has the following advantages at least: in the embodiment of the present invention, can improve the performance of channel estimating, especially can effectively improve the channel estimating performance under little bandwidth.
Description of drawings
In order to be illustrated more clearly in technical scheme of the present invention, in below describing embodiment, the accompanying drawing of required use is briefly described, apparently, accompanying drawing in the following describes is only some embodiments of the present invention, for those of ordinary skills, under the prerequisite of not paying creative work, can also obtain according to these accompanying drawings other accompanying drawing.
Fig. 1 is channel estimation process schematic diagram of the prior art;
Fig. 2 is the acquisition methods schematic flow sheet of a kind of domain channel response of providing of the embodiment of the present invention one;
Fig. 3 be a kind of domain channel response of providing of the embodiment of the present invention two obtain the device structure schematic diagram.
Embodiment
Below in conjunction with the accompanying drawing in the present invention, the technical scheme in the present invention is clearly and completely described, obviously, described embodiment is only a part of embodiment of the present invention, rather than whole embodiment.Based on the embodiment in the present invention, those of ordinary skills, not making under the creative work prerequisite the every other embodiment that obtains, belong to the scope of protection of the invention.
Embodiment one
The embodiment of the present invention one provides a kind of acquisition methods of domain channel response, and the method can be applied to LTE and LTE-A(LTE-Advanced, senior LTE) etc. system; And the method can be applied to the noise of little bandwidth (being that bandwidth is less than preset bandwidth value) and estimate and channel estimating (namely obtaining domain channel response) process, little bandwidth in the embodiment of the present invention mainly refers to: the time domain footpath number in noise window is less than the band width configuration of a threshold value, and this threshold value can be determined by emulation; And the method can be applied to the up-downgoing channel estimation process simultaneously; As shown in Figure 2, the acquisition methods of this domain channel response comprises the following steps:
Step 201, obtain the time domain channel response of each each frequency pilot sign of port of each reception antenna.
In the embodiment of the present invention, obtain the time domain channel response of each each frequency pilot sign of port of each reception antenna, specifically comprise: the channel estimation in frequency domain value of obtaining pilot frequency locations on each frequency pilot sign of each port of each reception antenna according to following formula:
Figure BDA00001631131300041
Adopt IDFT mode or mirror image IDFT mode that the channel estimation in frequency domain value transform of pilot frequency locations on each frequency pilot sign of each port of each reception antenna is arrived time domain, obtain the time domain channel response of each each frequency pilot sign of port of each reception antenna
Figure BDA00001631131300042
Wherein, r represents reception antenna, and p represents port, and l represents frequency pilot sign,
Figure BDA00001631131300043
The channel estimation in frequency domain value that represents pilot frequency locations on l frequency pilot sign of p port on r reception antenna, For each known pilot frequency sequence,
Figure BDA00001631131300045
The reception signal that takes out in the pilot frequency locations of each port for each reception antenna, and
Figure BDA00001631131300046
The time domain channel response that represents l frequency pilot sign of p port on r reception antenna.
Further, adopt the IDFT mode with each frequency pilot sign of each port of each reception antenna on the channel estimation in frequency domain value transform of pilot frequency locations in the process of time domain, obtain the time domain channel response of each each frequency pilot sign of port of each reception antenna
Figure BDA00001631131300047
Mode for by following formula:
Figure BDA00001631131300048
Adopt mirror image IDFT mode with each frequency pilot sign of each port of each reception antenna on the channel estimation in frequency domain value transform of pilot frequency locations in the process of time domain, obtain the time domain channel response of each each frequency pilot sign of port of each reception antenna
Figure BDA00001631131300051
Mode for by following formula:
h l r , p = IDFT ( H ′ ) , H ′ = H l r , p ( i ) , i = 1 , . . . , N pilot H l r , p ( 2 N pilot + 1 - i ) , i = N pilot + 1 , . . . , 2 N pilot ;
Wherein, above-mentioned N PliotFor the footpath number of time domain channel response is pilot tone point number on same orthogonal frequency division multiplex OFDM symbol.
It should be noted that for each frequency pilot sign of each port of each reception antenna, need to carry out respectively follow-up correlation step.
Step 202, utilize time domain channel response to calculate the first noise amplitude threshold value g.
In the embodiment of the present invention, utilize time domain channel response to calculate the process of the first noise amplitude threshold value g, specifically include but not limited to following mode:
Mode one, utilize following formula to calculate the first noise amplitude threshold value g:g=ρ a; Wherein, ρ is weight coefficient, can be determined by emulation, and a is all signal footpath amplitude mean values, and
Figure BDA00001631131300053
Mean represents the time domain channel response after taking absolute value
Figure BDA00001631131300054
Average.
Mode two, to utilize following formula to calculate the first noise amplitude threshold value g(g this moment be the signal footpath amplitude mean value that improves in noise window):
Figure BDA00001631131300055
J={j|j ∈ win}; Wherein, mean represents the time domain channel response after taking absolute value
Figure BDA00001631131300056
Average.
Need to prove, the difference that improves noise window and original noise window is, improves in noise window and can comprise part useful signal footpath, and namely the noise window original position need to not be as the criterion with the time delay maximum diameter, but need to as much as possible noise path be included.
In addition, the window method of getting of IDFT mode and mirror image IDFT mode is different, when adopting the IDFT mode to get window, and can be according to win={2N 1+ 1,2N 1+ 2 ..., N Pliot-N 1Method get window, when adopting mirror image IDFT mode to get window, can be according to win={N 2+ 1, N 2+ 2 ..., N Pliot, N Pliot+ 2, N Pliot+ 3 ..., 2N Pliot+ 1-N 2Method get window; Wherein, N 1<0.5N τ, N 2<2N τ, value can be determined by emulation, and N τFor the time-domain position of time delay maximum diameter, N PliotFor the footpath number of time domain channel response is pilot tone point number on same OFDM symbol.
Step 203, utilize time domain channel response and the first noise amplitude threshold value g to calculate the second noise amplitude threshold value G.
In the embodiment of the present invention, utilize time domain channel response and the first noise amplitude threshold value g to calculate the second noise amplitude threshold value G, comprise: utilize time domain channel response and the first noise amplitude threshold value g calculating noise amplitude mean value, and utilize noise amplitude mean value calculation the second noise amplitude threshold value G.
Utilize time domain channel response and the first noise amplitude threshold value g calculating noise amplitude mean value, comprise: be useful on all the signal footpath amplitudes and the first noise amplitude threshold value g that calculate the first noise amplitude threshold value g are compared, take out the signal footpath of range value lower than the first noise amplitude threshold value g, thereby calculate noise amplitude mean value; Concrete, can utilize following formula calculating noise amplitude mean value
Figure BDA00001631131300061
A l r , p = mean ( | h l r , p ( k ) | ) , And k = { k | | h l r , p ( k ) | < g } ; k = { k | | h l r , p ( k ) | < g } Expression is with the time domain channel response after the taking absolute value of each each frequency pilot sign of port of each reception antenna
Figure BDA00001631131300065
All signal footpath amplitudes and the first noise amplitude threshold value g compare, the signal footpath range value of taking-up is lower than the signal footpath of the first noise amplitude threshold value g, mean represents the time domain channel response after taking absolute value
Figure BDA00001631131300066
Average.
Utilize noise amplitude mean value calculation the second noise amplitude threshold value G, comprising: utilize following formula to calculate the second noise amplitude threshold value G:
Figure BDA00001631131300067
Expression is with the noise amplitude mean value on each frequency pilot sign of each port of each reception antenna
Figure BDA00001631131300068
Multiply by 2 after averaging, obtain the second noise amplitude threshold value G.
In the embodiment of the present invention,, based on two noise amplitude threshold value G obtained above, can also carry out noise and estimate, to estimate accurately the noise power under little bandwidth, thereby, according to the noise of noise power-value filtering pilot tone point channel estimating, improve the detection performance.
In the embodiment of the present invention, can take out
Figure DEST_PATH_GDA00002592097200018
In all range values lower than the signal footpath of the second noise amplitude threshold value G as noise path, and calculating noise footpath power average value is the noise power estimation value; Concrete, can utilize following formula calculating noise power estimated value
Figure DEST_PATH_GDA00002592097200019
And m = { m | | h l r , p ( m ) | < G } ; m = { m | | h l r , p ( m ) | < G } Expression is with the time domain channel response after the taking absolute value of each each frequency pilot sign of port of each reception antenna
Figure DEST_PATH_GDA000025920972000112
In all range values lower than the signal footpath of the second noise amplitude threshold value G as noise path, mean represents the time domain channel response after taking absolute value
Figure DEST_PATH_GDA000025920972000113
Square average.
It should be noted that when adopting mirror image IDFT, due to the signal footpath
Figure BDA00001631131300071
The range value perseverance be 0, therefore, if adopt mirror image IDFT, also need the signal footpath
Figure BDA00001631131300072
Remove.
Step 204, utilize the second noise amplitude threshold value G to calculate the 3rd noise amplitude threshold value G '.
In the embodiment of the present invention, utilize the second noise amplitude threshold value G to calculate the 3rd noise amplitude threshold value G ', comprising: according to following formula, calculate the 3rd noise amplitude threshold value G ': G '=min (λ 1H max, λ 2G); Wherein, min represents to get λ 1Multiply by h maxAnd λ 2Multiply by the little value of the second noise amplitude threshold value G, λ 1And λ 2, for default value, can determine by emulation h maxFor
Figure DEST_PATH_GDA00002592097200023
Each footpath in the maximum of amplitude.
Step 205, utilize the domain channel response after the 3rd noise amplitude threshold value G ' obtains noise suppression.In the embodiment of the present invention, utilize the domain channel response after the 3rd noise amplitude threshold value G ' obtains noise suppression, comprising: determine In all range values higher than the 3rd noise amplitude threshold value G ', and the signal footpath in original noise window is not initial useful signal footpath; Determine that the signal footpath in initial useful signal footpath and its both sides Δ scope is the useful signal footpath, Δ is default value, can determine by emulation; , with other footpath zero setting except the useful signal footpath in all signal footpaths, obtain the time-domain signal response after noise suppression
Figure BDA00001631131300076
Time-domain signal response transform after noise suppression to frequency domain, is obtained the domain channel response after noise suppression
Figure BDA00001631131300077
Further, in employing IDFT obtains the process of domain channel response, can realize by following formula: In employing mirror image IDFT obtains the process of domain channel response, can realize by following formula: H &prime; = DFT ( h ~ l r , p ) , H l r , p ( i ) = 0.5 H &prime; ( i ) + 0.5 H &prime; ( 2 N pilot + 1 - i ) .
In sum, in the embodiment of the present invention, can estimate comparatively accurately the noise power under little bandwidth, improve and detect performance; And can effectively improve channel estimating performance under little bandwidth, especially under the long time delay scene, performance boost is obvious.
Embodiment two
Based on the inventive concept same with said method, a kind of equipment that obtains of domain channel response also is provided in the embodiment of the present invention, as shown in Figure 3, this equipment comprises:
The first computing module 11, be used for utilizing time domain channel response to calculate the first noise amplitude threshold value g;
The second computing module 12, be used for utilizing described time domain channel response and described the first noise amplitude threshold value g to calculate the second noise amplitude threshold value G;
The 3rd computing module 13, be used for utilizing described the second noise amplitude threshold value G to calculate the 3rd noise amplitude threshold value G ';
Obtain module 14, be used for utilizing the domain channel response after described the 3rd noise amplitude threshold value G ' obtains noise suppression.
Described acquisition module 14, also for obtain the channel estimation in frequency domain value of pilot frequency locations on each frequency pilot sign of each port of each reception antenna according to following formula:
Figure BDA00001631131300081
Wherein,
Figure BDA00001631131300082
The channel estimation in frequency domain value that represents pilot frequency locations on l frequency pilot sign of p port on r reception antenna,
Figure BDA00001631131300083
For known pilot frequency sequence,
Figure BDA00001631131300084
The reception signal that takes out in the pilot frequency locations of each port for each reception antenna;
Adopt inverse discrete Fourier transform IDFT mode or mirror image IDFT mode that the channel estimation in frequency domain value transform of pilot frequency locations on each frequency pilot sign of each port of each reception antenna is arrived time domain, obtain the time domain channel response of each each frequency pilot sign of port of each reception antenna
Figure BDA00001631131300085
Described the first computing module 11, concrete for utilizing following formula to calculate the first noise amplitude threshold value g:g=ρ a; Wherein, ρ is weight coefficient, The time domain channel response that represents l frequency pilot sign of p port on r reception antenna, mean represent the time domain channel response after taking absolute value Average; Perhaps,
Utilize following formula to calculate the first noise amplitude threshold value g: J={j|j ∈ win}; Wherein,
Figure BDA000016311313000810
The time domain channel response that represents l frequency pilot sign of p port on r reception antenna, mean represent the time domain channel response after taking absolute value Average; When adopting the IDFT mode to get window, win={2N 1+ 1,2N 1+ 2 ..., N Pliot-N 1, when adopting mirror image IDFT mode to get window, win={N 2+ 1, N 2+ 2 ..., N Pliot, N Pliot+ 2, N Pliot+ 3 ..., 2N Pliot+ 1-N 2; And N 1<0.5N τ, N 2<2N τ, and N τFor the time-domain position of time delay maximum diameter, N PliotFor the footpath number of time domain channel response is pilot tone point number on same orthogonal frequency division multiplex OFDM symbol.
Described the second computing module 12, concrete being used for utilized described time domain channel response and described the first noise amplitude threshold value g calculating noise amplitude mean value, and utilizes described the second noise amplitude threshold value G of described noise amplitude mean value calculation.
Described the second computing module 12, be further used for utilizing following formula to calculate described noise amplitude mean value
Figure BDA00001631131300092
And
Figure BDA00001631131300093
Wherein,
Figure BDA00001631131300094
The time domain channel response that represents l frequency pilot sign of p port on r reception antenna,
Figure BDA00001631131300095
Expression is with the time domain channel response after the taking absolute value of each each frequency pilot sign of port of each reception antenna
Figure BDA00001631131300096
All signal footpath amplitudes and the first noise amplitude threshold value g compare, the signal footpath range value of taking-up is lower than the signal footpath of the first noise amplitude threshold value g, mean represents the time domain channel response after taking absolute value
Figure BDA00001631131300097
Average.
Described the second computing module 12 is further used for utilizing following formula to calculate described the second noise amplitude threshold value G:
Figure BDA00001631131300098
Its expression is with the noise amplitude mean value on each frequency pilot sign of each port of each reception antenna
Figure BDA00001631131300099
Multiply by 2 after averaging, obtain described the second noise amplitude threshold value G.
This equipment also comprises: the 4th computing module 15 is used for utilizing following formula calculating noise power estimated value
Figure BDA000016311313000910
( &sigma; l r , p ) 2 = mean ( | h l r , p ( m ) | 2 ) , And m = { m | | h l r , p ( m ) | < G } ; Wherein,
Figure BDA000016311313000913
The time domain channel response that represents l frequency pilot sign of p port on r reception antenna,
Figure BDA000016311313000914
Expression is with the time domain channel response after the taking absolute value of each each frequency pilot sign of port of each reception antenna
Figure BDA000016311313000915
In all range values lower than the signal footpath of the second noise amplitude threshold value G as noise path, mean represents the time domain channel response after taking absolute value
Figure BDA000016311313000916
Square average.
Described the 3rd computing module 13, concrete being used for calculated the 3rd noise amplitude threshold value G ' according to following formula: Wherein, min represents to get λ 1Multiply by h maxAnd λ 2Multiply by the little value of the second noise amplitude threshold value G, γ 1And λ 2For default value, h maxFor
Figure BDA000016311313000918
Each footpath in the maximum of amplitude,
Figure BDA000016311313000919
The time domain channel response that represents l frequency pilot sign of p port on r reception antenna.
Described acquisition module 14, concrete being used for determined In all range values higher than the 3rd noise amplitude threshold value G ', and the signal footpath in original noise window is not initial useful signal footpath, and is described
Figure BDA000016311313000921
The time domain channel response that represents l frequency pilot sign of p port on r reception antenna; Determine that the signal footpath in described initial useful signal footpath and its both sides Δ scope is the useful signal footpath, described Δ is default value; , with other footpath zero setting except described useful signal footpath in all signal footpaths, obtain the time-domain signal response after noise suppression; Time-domain signal response transform after described noise suppression to frequency domain, is obtained the domain channel response after noise suppression.
Wherein, the modules of apparatus of the present invention can be integrated in one, and also can separate deployment.Above-mentioned module can be merged into a module, also can further split into a plurality of submodules.
Through the above description of the embodiments, those skilled in the art can be well understood to the present invention and can realize by the mode that software adds essential general hardware platform, can certainly pass through hardware, but in a lot of situation, the former is better execution mode.Based on such understanding, the part that technical scheme of the present invention contributes to prior art in essence in other words can embody with the form of software product, this computer software product is stored in a storage medium, comprise that some instructions are with so that a computer equipment (can be personal computer, server, the perhaps network equipment etc.) carry out the described method of each embodiment of the present invention.
It will be appreciated by those skilled in the art that accompanying drawing is the schematic diagram of a preferred embodiment, the module in accompanying drawing or flow process might not be that enforcement the present invention is necessary.
It will be appreciated by those skilled in the art that the module in the device in embodiment can be distributed in the device of embodiment according to the embodiment description, also can carry out respective change and be arranged in the one or more devices that are different from the present embodiment.The module of above-described embodiment can be merged into a module, also can further split into a plurality of submodules.
The invention described above embodiment sequence number, just to describing, does not represent the quality of embodiment.
Above disclosed be only several specific embodiment of the present invention, still, the present invention is not limited thereto, the changes that any person skilled in the art can think of all should fall into protection scope of the present invention.

Claims (18)

1. the acquisition methods of a domain channel response, is characterized in that, comprising:
Utilize time domain channel response to calculate the first noise amplitude threshold value g;
Utilize described time domain channel response and described the first noise amplitude threshold value g to calculate the second noise amplitude threshold value G;
Utilize described the second noise amplitude threshold value G to calculate the 3rd noise amplitude threshold value G ';
Utilize the domain channel response after described the 3rd noise amplitude threshold value G ' obtains noise suppression.
2. the method for claim 1, is characterized in that, the described time domain channel response that utilizes calculates the first noise amplitude threshold value g, also comprises before:
Obtain the channel estimation in frequency domain value of pilot frequency locations on each frequency pilot sign of each port of each reception antenna according to following formula:
Figure FDA00001631131200011
Wherein,
Figure FDA00001631131200012
The channel estimation in frequency domain value that represents pilot frequency locations on l frequency pilot sign of p port on r reception antenna,
Figure FDA00001631131200013
For known pilot frequency sequence, The reception signal that takes out in the pilot frequency locations of each port for each reception antenna;
Adopt inverse discrete Fourier transform IDFT mode or mirror image IDFT mode that the channel estimation in frequency domain value transform of pilot frequency locations on each frequency pilot sign of each port of each reception antenna is arrived time domain, obtain the time domain channel response of each each frequency pilot sign of port of each reception antenna
Figure FDA00001631131200015
3. the method for claim 1, is characterized in that, the described time domain channel response that utilizes calculates the first noise amplitude threshold value g, comprising:
Utilize following formula to calculate the first noise amplitude threshold value g:g=ρ a; Wherein, ρ is weight coefficient,
Figure FDA00001631131200017
The time domain channel response that represents l frequency pilot sign of p port on r reception antenna, mean represent the time domain channel response after taking absolute value
Figure FDA00001631131200018
Average; Perhaps,
Utilize following formula to calculate the first noise amplitude threshold value g:
Figure FDA00001631131200019
J={j|j ∈ win}; Wherein, The time domain channel response that represents l frequency pilot sign of p port on r reception antenna, mean represent the time domain channel response after taking absolute value
Figure FDA000016311312000111
Average; When adopting the IDFT mode to get window, win={2N 1+ 1,2N 1+ 2 ..., N Pliot-N 1, when adopting mirror image IDFT mode to get window, win={N 2+ 1, N 2+ 2 ..., N Plion, N Pliot+ 2, N Pliot+ 3 ..., 2N Pliot+ 1-N 2; And N 1<0.5N τ, N 2<2N τ, and N τFor the time-domain position of time delay maximum diameter, N PliotFor the footpath number of time domain channel response is pilot tone point number on same orthogonal frequency division multiplex OFDM symbol.
4. the method for claim 1, is characterized in that, utilizes described time domain channel response and described the first noise amplitude threshold value g to calculate the second noise amplitude threshold value G, comprising:
Utilize described time domain channel response and described the first noise amplitude threshold value g calculating noise amplitude mean value, and utilize described the second noise amplitude threshold value G of described noise amplitude mean value calculation.
5. method as claimed in claim 4, is characterized in that, utilizes described time domain channel response and described the first noise amplitude threshold value g calculating noise amplitude mean value, comprising:
Utilize following formula to calculate described noise amplitude mean value
Figure FDA00001631131200021
Figure FDA00001631131200022
And
Figure FDA00001631131200023
Wherein,
Figure FDA00001631131200024
The time domain channel response that represents l frequency pilot sign of p port on r reception antenna,
Figure FDA00001631131200025
Expression is with the time domain channel response after the taking absolute value of each each frequency pilot sign of port of each reception antenna
Figure FDA00001631131200026
All signal footpath amplitudes and the first noise amplitude threshold value g compare, the signal footpath range value of taking-up is lower than the signal footpath of the first noise amplitude threshold value g, mean represents the time domain channel response after taking absolute value
Figure FDA00001631131200027
Average.
6. method as claimed in claim 4, is characterized in that, utilizes described the second noise amplitude threshold value G of described noise amplitude mean value calculation, comprising:
Utilize following formula to calculate described the second noise amplitude threshold value G:
Figure FDA00001631131200028
Its expression is with the noise amplitude mean value on each frequency pilot sign of each port of each reception antenna
Figure FDA00001631131200029
Multiply by 2 after averaging, obtain described the second noise amplitude threshold value G.
7. the method for claim 1, is characterized in that, utilizes described time domain channel response and described the first noise amplitude threshold value g to calculate the second noise amplitude threshold value G, also comprises afterwards:
Utilize following formula calculating noise power estimated value
Figure FDA000016311312000210
Figure FDA000016311312000211
And
Figure FDA000016311312000212
Wherein, The time domain channel response that represents l frequency pilot sign of p port on r reception antenna, Expression is with the time domain channel response after the taking absolute value of each each frequency pilot sign of port of each reception antenna
Figure FDA00001631131200032
In all range values lower than the signal footpath of the second noise amplitude threshold value G as noise path, mean represents the time domain channel response after taking absolute value
Figure FDA00001631131200033
Square average.
8. the method for claim 1, is characterized in that, utilizes described the second noise amplitude threshold value G to calculate the 3rd noise amplitude threshold value G ', comprising:
Calculate the 3rd noise amplitude threshold value G ': G '=min (λ according to following formula 1H max, λ 2G); Wherein, min represents to get λ 1Multiply by h maxAnd λ 2Multiply by the little value of the second noise amplitude threshold value G, λ 1And λ 2For default value, h maxFor
Figure RE-FDA00002592097100014
Each footpath in the maximum of amplitude,
Figure RE-FDA00002592097100015
The time domain channel response that represents l frequency pilot sign of p port on r reception antenna.
9. the method for claim 1, is characterized in that, utilizes the domain channel response after described the 3rd noise amplitude threshold value G ' obtains noise suppression, comprising:
Determine
Figure FDA00001631131200037
In all range values higher than the 3rd noise amplitude threshold value G ', and the signal footpath in original noise window is not initial useful signal footpath, and is described
Figure FDA00001631131200038
The time domain channel response that represents l frequency pilot sign of p port on r reception antenna;
Determine that the signal footpath in described initial useful signal footpath and its both sides Δ scope is the useful signal footpath, described Δ is default value;
, with other footpath zero setting except described useful signal footpath in all signal footpaths, obtain the time-domain signal response after noise suppression;
Time-domain signal response transform after described noise suppression to frequency domain, is obtained the domain channel response after noise suppression.
10. the equipment that obtains of a domain channel response, is characterized in that, comprising:
The first computing module, be used for utilizing time domain channel response to calculate the first noise amplitude threshold value g;
The second computing module, be used for utilizing described time domain channel response and described the first noise amplitude threshold value g to calculate the second noise amplitude threshold value G;
The 3rd computing module, be used for utilizing described the second noise amplitude threshold value G to calculate the 3rd noise amplitude threshold value G ';
Obtain module, be used for utilizing the domain channel response after described the 3rd noise amplitude threshold value G ' obtains noise suppression.
11. equipment as claimed in claim 10, is characterized in that,
Described acquisition module, also for obtain the channel estimation in frequency domain value of pilot frequency locations on each frequency pilot sign of each port of each reception antenna according to following formula:
Figure FDA00001631131200041
Wherein, The channel estimation in frequency domain value that represents pilot frequency locations on l frequency pilot sign of p port on r reception antenna,
Figure FDA00001631131200043
For known pilot frequency sequence,
Figure FDA00001631131200044
The reception signal that takes out in the pilot frequency locations of each port for each reception antenna;
Adopt inverse discrete Fourier transform IDFT mode or mirror image IDFT mode that the channel estimation in frequency domain value transform of pilot frequency locations on each frequency pilot sign of each port of each reception antenna is arrived time domain, obtain the time domain channel response of each each frequency pilot sign of port of each reception antenna
Figure FDA00001631131200045
12. equipment as claimed in claim 10, is characterized in that,
Described the first computing module, concrete for utilizing following formula to calculate the first noise amplitude threshold value g:g=ρ a; Wherein, ρ is weight coefficient,
Figure FDA00001631131200046
The time domain channel response that represents l frequency pilot sign of p port on r reception antenna, mean represent the time domain channel response after taking absolute value
Figure FDA00001631131200048
Average; Perhaps,
Utilize following formula to calculate the first noise amplitude threshold value g:
Figure FDA00001631131200049
J={j|j ∈ win}; Wherein,
Figure FDA000016311312000410
The time domain channel response that represents l frequency pilot sign of p port on r reception antenna, mean represent the time domain channel response after taking absolute value
Figure FDA000016311312000411
Average; When adopting the IDFT mode to get window, win={2N 1+ 1,2N 1+ 2 ..., N Pliot-N 1, when adopting mirror image IDFT mode to get window, win={N 2+ 1, N 2+ 2 ..., N Pliot, N Pliot+ 2, N Pliot+ 3 ..., 2N Pliot+ 1-N 2; And N 1<0.5N τ, N 2<2N τ, and N τFor the time-domain position of time delay maximum diameter, N PliotFor the footpath number of time domain channel response is pilot tone point number on same orthogonal frequency division multiplex OFDM symbol.
13. equipment as claimed in claim 10, is characterized in that,
Described the second computing module, concrete being used for utilized described time domain channel response and described the first noise amplitude threshold value g calculating noise amplitude mean value, and utilizes described the second noise amplitude threshold value G of described noise amplitude mean value calculation.
14. equipment as claimed in claim 13, is characterized in that,
Described the second computing module, be further used for utilizing following formula to calculate described noise amplitude mean value
Figure FDA00001631131200051
Figure FDA00001631131200052
And
Figure FDA00001631131200053
Wherein,
Figure FDA00001631131200054
The time domain channel response that represents l frequency pilot sign of p port on r reception antenna,
Figure FDA00001631131200055
Expression is with the time domain channel response after the taking absolute value of each each frequency pilot sign of port of each reception antenna
Figure FDA00001631131200056
All signal footpath amplitudes and the first noise amplitude threshold value g compare, the signal footpath range value of taking-up is lower than the signal footpath of the first noise amplitude threshold value g, mean represents the time domain channel response after taking absolute value
Figure FDA00001631131200057
Average.
15. equipment as claimed in claim 13, is characterized in that,
Described the second computing module is further used for utilizing following formula to calculate described the second noise amplitude threshold value G: Its expression is with the noise amplitude mean value on each frequency pilot sign of each port of each reception antenna
Figure FDA00001631131200059
Multiply by 2 after averaging, obtain described the second noise amplitude threshold value G.
16. equipment as claimed in claim 10, is characterized in that, also comprises:
The 4th computing module, be used for utilizing following formula calculating noise power estimated value
Figure FDA000016311312000510
Figure FDA000016311312000511
And
Figure FDA000016311312000512
Wherein,
Figure FDA000016311312000513
The time domain channel response that represents l frequency pilot sign of p port on r reception antenna,
Figure FDA000016311312000514
Expression is with the time domain channel response after the taking absolute value of each each frequency pilot sign of port of each reception antenna
Figure FDA000016311312000515
In all range values lower than the signal footpath of the second noise amplitude threshold value G as noise path, mean represents the time domain channel response after taking absolute value
Figure FDA000016311312000516
Square average.
17. equipment as claimed in claim 10, is characterized in that,
Described the 3rd computing module, concrete being used for calculated the 3rd noise amplitude threshold value G ' according to following formula:
Figure FDA000016311312000517
Wherein, min represents to get λ 1Multiply by h maxAnd λ 2Multiply by the little value of the second noise amplitude threshold value G, λ 1And λ 2For default value, h maxFor
Figure FDA00001631131200061
Each footpath in the maximum of amplitude,
Figure FDA00001631131200062
The time domain channel response that represents l frequency pilot sign of p port on r reception antenna.
18. equipment as claimed in claim 10, is characterized in that,
Described acquisition module, concrete being used for determined
Figure FDA00001631131200063
In all range values higher than the 3rd noise amplitude threshold value G ', and the signal footpath in original noise window is not initial useful signal footpath, and is described The time domain channel response that represents l frequency pilot sign of p port on r reception antenna;
Determine that the signal footpath in described initial useful signal footpath and its both sides Δ scope is the useful signal footpath, described Δ is default value;
, with other footpath zero setting except described useful signal footpath in all signal footpaths, obtain the time-domain signal response after noise suppression;
Time-domain signal response transform after described noise suppression to frequency domain, is obtained the domain channel response after noise suppression.
CN201210147820.4A 2012-05-11 2012-05-11 The acquisition methods of a kind of domain channel response and equipment Active CN103391266B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201210147820.4A CN103391266B (en) 2012-05-11 2012-05-11 The acquisition methods of a kind of domain channel response and equipment

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201210147820.4A CN103391266B (en) 2012-05-11 2012-05-11 The acquisition methods of a kind of domain channel response and equipment

Publications (2)

Publication Number Publication Date
CN103391266A true CN103391266A (en) 2013-11-13
CN103391266B CN103391266B (en) 2016-12-14

Family

ID=49535426

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201210147820.4A Active CN103391266B (en) 2012-05-11 2012-05-11 The acquisition methods of a kind of domain channel response and equipment

Country Status (1)

Country Link
CN (1) CN103391266B (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110445733A (en) * 2019-06-27 2019-11-12 熊军 Iteration self-adapting channel denoising method and iteration self-adapting channel denoise device
CN111371436A (en) * 2020-03-16 2020-07-03 中国人民解放军空军研究院战略预警研究所 Radar antenna scanning period measuring method, device, equipment and storage medium

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101127745A (en) * 2006-08-16 2008-02-20 大唐移动通信设备有限公司 A chancel estimation method and device
CN101267422A (en) * 2008-03-10 2008-09-17 电子科技大学 A frequency domain channel estimation method for OFDM multiplex system
EP2120359A1 (en) * 2007-02-07 2009-11-18 NEC Corporation Channel estimation device, equalization device, and equalization method in the estimation
CN101707574A (en) * 2009-11-27 2010-05-12 北京邮电大学 Channel estimation method and device
CN101742629A (en) * 2008-11-26 2010-06-16 中兴通讯股份有限公司 Channel estimation method and device of TD-SCDMA system terminal

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101127745A (en) * 2006-08-16 2008-02-20 大唐移动通信设备有限公司 A chancel estimation method and device
EP2120359A1 (en) * 2007-02-07 2009-11-18 NEC Corporation Channel estimation device, equalization device, and equalization method in the estimation
CN101267422A (en) * 2008-03-10 2008-09-17 电子科技大学 A frequency domain channel estimation method for OFDM multiplex system
CN101742629A (en) * 2008-11-26 2010-06-16 中兴通讯股份有限公司 Channel estimation method and device of TD-SCDMA system terminal
CN101707574A (en) * 2009-11-27 2010-05-12 北京邮电大学 Channel estimation method and device

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110445733A (en) * 2019-06-27 2019-11-12 熊军 Iteration self-adapting channel denoising method and iteration self-adapting channel denoise device
CN111371436A (en) * 2020-03-16 2020-07-03 中国人民解放军空军研究院战略预警研究所 Radar antenna scanning period measuring method, device, equipment and storage medium

Also Published As

Publication number Publication date
CN103391266B (en) 2016-12-14

Similar Documents

Publication Publication Date Title
CN101534266B (en) Channel estimation method for Orthogonal Frequency Division Multiplexing system and device
CN101917363B (en) Method and device for estimating Doppler frequency shift
WO2015139590A1 (en) Frequency offset estimation and compensation method and device
JP2020519162A (en) Method and apparatus for determining uplink synchronization timing deviation
EP2515459A1 (en) Frequency offset compensation and equalization method and device in an orthogonal frequency division multiplexing system
WO2015165354A1 (en) Power delay profile (pdp) estimation method and device
WO2016128728A3 (en) Obtaining channel state information in a multicarrier wireless communication network
JP5592421B2 (en) MMSE channel estimation apparatus and method considering offset compensation in wireless communication system
CN104052555B (en) A kind of method of radio channel multi-path parameter Estimation under ofdm system
CN101635598B (en) Method and device for estimating noise power
EP3028392B1 (en) Methods, systems and computer readable media for simulating per user equipment (ue) doppler shifts for testing air interface devices
CN105072063B (en) A kind of method and apparatus for inhibiting interference signal
EP2924936B1 (en) Channel estimation method and device
CN103391266A (en) Frequency domain channel response obtaining method and device
CN102113285A (en) A simplified equalizationscheme for distributed resource allocation in multi-carrier systems
CN104243124A (en) Computing method and device of sounding reference signal channel coefficient matrix
CN102104433B (en) Symbol detection method and equipment
CN106302286A (en) A kind of signal to noise ratio blind estimating method of ofdm system
CN104753834B (en) A kind of channel estimation methods and device
EP2991298B1 (en) Channel equalization method and device, and receiver
CN102957641A (en) Channel estimation method of TDS (Total Dissolved Solids)-OFDM (Orthogonal Frequency Division Multiplexing) system
CN104243365A (en) Estimation method for noise variance
KR100733596B1 (en) System and method for estimating channel state information used in adaptive wireless transceiver
Wypich et al. 5G NR Based Radar-on-Demand Using Channel Impulse Response Estimate
CN102281170B (en) Speed measuring method and device based on channel estimation

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
CP01 Change in the name or title of a patent holder

Address after: 100191 No. 40, Haidian District, Beijing, Xueyuan Road

Patentee after: CHINA ACADEMY OF TELECOMMUNICATIONS TECHNOLOGY

Address before: 100191 No. 40, Haidian District, Beijing, Xueyuan Road

Patentee before: CHINA ACADEMY OF TELECOMMUNICATIONS TECHNOLOGY

CP01 Change in the name or title of a patent holder
TR01 Transfer of patent right

Effective date of registration: 20210604

Address after: 100085 1st floor, building 1, yard 5, Shangdi East Road, Haidian District, Beijing

Patentee after: DATANG MOBILE COMMUNICATIONS EQUIPMENT Co.,Ltd.

Address before: 100191 No. 40, Haidian District, Beijing, Xueyuan Road

Patentee before: CHINA ACADEMY OF TELECOMMUNICATIONS TECHNOLOGY

TR01 Transfer of patent right