CN110602006A - Channel estimation method under LTE system - Google Patents

Channel estimation method under LTE system Download PDF

Info

Publication number
CN110602006A
CN110602006A CN201910812011.2A CN201910812011A CN110602006A CN 110602006 A CN110602006 A CN 110602006A CN 201910812011 A CN201910812011 A CN 201910812011A CN 110602006 A CN110602006 A CN 110602006A
Authority
CN
China
Prior art keywords
channel estimation
channel
frequency domain
mean square
square error
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.)
Pending
Application number
CN201910812011.2A
Other languages
Chinese (zh)
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.)
SHENZHEN HIPAD COMMUNICATION TECHNOLOGY Co Ltd
Original Assignee
SHENZHEN HIPAD COMMUNICATION TECHNOLOGY Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by SHENZHEN HIPAD COMMUNICATION TECHNOLOGY Co Ltd filed Critical SHENZHEN HIPAD COMMUNICATION TECHNOLOGY Co Ltd
Priority to CN201910812011.2A priority Critical patent/CN110602006A/en
Publication of CN110602006A publication Critical patent/CN110602006A/en
Pending legal-status Critical Current

Links

Classifications

    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L25/00Baseband systems
    • H04L25/02Details ; arrangements for supplying electrical power along data transmission lines
    • H04L25/0202Channel estimation
    • H04L25/024Channel estimation channel estimation algorithms
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L25/00Baseband systems
    • H04L25/02Details ; arrangements for supplying electrical power along data transmission lines
    • H04L25/0202Channel estimation
    • H04L25/024Channel estimation channel estimation algorithms
    • H04L25/0256Channel estimation using minimum mean square error criteria

Abstract

The invention provides a channel estimation method under an LTE system, which comprises the steps of firstly carrying out LS least square method channel estimation on a frequency domain at a receiver side to obtain a preliminary channel estimation result, then carrying out MMSE minimum mean square error channel estimation to improve the accuracy of channel estimation, then carrying out IFFT (inverse fast Fourier transform) on the frequency domain channel estimation result, utilizing the time limitation of channel impulse response to carry out truncation on signals except for channel duration, forcibly setting zero, reducing the energy of noise, and finally carrying out FFT (fast Fourier transform) to return to the frequency domain, thereby realizing the improvement of the accuracy of channel estimation. The invention has the beneficial effects that: the mean square error of the channel estimation is greatly reduced, the precision of the channel estimation is improved, and a good effect is obtained.

Description

Channel estimation method under LTE system
Technical Field
The present invention relates to a channel estimation method, and in particular, to a channel estimation method in an LTE system.
Background
In the current information society, the development of mobile communication technology is receiving more and more attention. The ultimate goal of the development of mobile communication technology is to enable anyone to communicate with anyone or an object of any kind at any time and any place.
The development of communication technology has passed through the continuous development of the first, second, third and current fourth generation LTE, and even the continuous evolution of the future 5G. The development of communication technology has taught that the most core technology of communication is coding, decoding, modulating and demodulating. However, correct demodulation and decoding of data information first relies on accurate estimation of the radio channel.
The channel estimation can be performed in a time domain or a frequency domain, and in an LTE (Long Term Evolution) system, the channel estimation is generally performed in the frequency domain in consideration of the wideband characteristics and the corresponding complexity of LTE.
Therefore, how to improve the accuracy of channel estimation is an urgent technical problem to be solved by those skilled in the art.
Disclosure of Invention
In order to solve the problems in the prior art, the invention provides a channel estimation method in an LTE system.
The invention provides a channel estimation method under an LTE system, which comprises the steps of firstly carrying out LS least square method channel estimation on a frequency domain at a receiver side to obtain a preliminary channel estimation result, then carrying out MMSE minimum mean square error channel estimation to improve the accuracy of channel estimation, then carrying out IFFT (inverse fast Fourier transform) on the frequency domain channel estimation result, utilizing the time limitation of channel impulse response to carry out truncation on signals except for channel duration, forcibly setting zero, reducing the energy of noise, and finally carrying out FFT (fast Fourier transform) to return to the frequency domain, thereby realizing the improvement of the accuracy of channel estimation.
As a further improvement of the present invention, the channel estimation method comprises the steps of:
s1, frequency domain equalization;
s2, LS least square method channel estimation;
s3, MMSE minimum mean square error channel estimation;
s4, noise estimation;
s5, time domain truncation.
As a further improvement of the present invention, in step S1, the system equation of the channel estimation is expressed in the form of discrete convolution:
wherein the content of the first and second substances,in order to input the signal, the signal is,the response of the channel impulse is determined,in order to observe the signal(s),for input signalsLength of (d);
expressed in matrix form as follows:
wherein the content of the first and second substances,is composed of
Is the impulse response of the channel;
is an observed signal;
is complex gaussian noise;
according to the property of discrete Fourier transform, multiplying time domain circular convolution by a frequency domain;
to apply frequency domain multiplication, the relationship of the above expressions is first converted to the form of a circular convolution by the following methodCut into lengths ofThen add up:
the observation signal and noise are also treated in the same way:
the signals thus processed are related to the convolution:
expressed in a matrix, then:
andis formed byAnda column vector formed by the arrangement of the column vectors,
within the matrix theory, the matrices are represented in upper case bold and the vectors are represented in lower case bold, and in order to comply with this rule, the expression is used hereinRepresenting frequency domain signals to avoid confusion, if at allArranged as vectorsHandle barArranged as vectorsThen the discrete fourier transform can be represented by a matrix:
inverse transformation:
whereinA discrete Fourier transform matrix of N points;
note the book
Equation of pairsDFT conversion is carried out on two sides:
from this, the equation in the frequency domain is expressed as:
wherein the content of the first and second substances,
according to the convolution theorem, the time domain circular convolution corresponds to the frequency domain product,must be a diagonal matrix:
as a further improvement of the present invention, in step S2,
LS least squares channel estimation does not take noise effects into account, so equations in the frequency domainThe method comprises the following steps:
constructing an error vector:
where, e is the error vector and,is a received frequency-domain signal that is,is a frequency-domain signal that is transmitted,
is a frequency domain LS least squares method channel estimation,presentation pairEstimating a signal vector;
to minimize the mean square error of the signal estimate, the condition needs to be satisfied:
transforming the above equation:
the above J-pair matrixCalculating the partial derivatives, and making the partial derivatives zero, then:
namely:
therefore, the method comprises the following steps:
thereby obtaining a frequency domain LS least squares channel estimate
As a further improvement of the present invention, in step S3,
MMSE minimum mean square error channel estimation fully considers the influence of noise, further carries out minimum mean square error channel estimation on the basis of completing the preliminary estimation of the channel by an LS least square method,
when MMSE minimum mean square error channel estimation is adopted, orderFor equation in frequency domainConstructing an error vector:
where, e is the error vector and,is a characteristic of the frequency domain of the channel,is channel frequency domain MMSE minimum mean square error estimation;
to minimize the mean square error of the signal estimate, the condition needs to be satisfied:
transforming the above equation:
the above J-pair matrixCalculating the partial derivatives, and making the partial derivatives zero, then:
namely:
therefore, the method comprises the following steps:
order toRepresents the noise power, then
Therefore, there are
Wherein, in order to simplify the calculation, the actual engineering realization is considered inIn the calculation of (1) inApproximate substitutionThereby obtaining a frequency domain MMSE minimum mean square error channel estimation
As a further improvement of the present invention, in step S4,
the upper typeIn which there is also a noise termNeed to estimate, in order to estimate, the system noiseThe following calculations were performed:
suppose inOn one OFDM symbolThe received signal on the subcarriers is represented as:
whereinIs the number of sub-carriers,represents the m th OFDM symbolThe channel frequency response on the sub-carriers,represents the m th OFDM symbolThe symbols transmitted on the individual sub-carriers,represents the m th OFDM symbolAdditive Gaussian stationary noise on subcarriers with mean of 0 and variance of
Since the following steps are performed for each subcarrier k, for convenience of description, the subscript k is omitted in the following equation, and the pilot symbol positions are taken as 3 and 10 as examples, that is, the uplink LTE pilot symbol positions:
1) calculating channel frequency response by LS channel estimation
2) Calculating channel frequency domain response difference;
3) Calculating a state compensation quantity;
4) ComputingWherein
5) ByCalculating
6) Computing
7) ComputingAnd thus a noise variance estimate.
As a further improvement of the present invention, in step S5,
obtained by adopting frequency domain MMSE minimum mean square error channel estimation methodThen, more accurate channel estimation is obtained through a time domain truncation processing method;
in a wireless propagation environment, the delay spread of the channel, i.e. the coherence bandwidth, is substantially constant, which can be considered as a known quantity, and is further obtained by the already obtained LS channel estimate:
1) calculating an autocorrelation function of frequency domain channel response according to LS channel estimation;
2) determining a 3db bandwidth of the autocorrelation function;
3) calculating frequency domain correlation bandwidth
4) Determining root mean square delay spread
5) According to the Gaussian normal probability distribution, the maximum time delay expansion probability is less thanGet it
Assuming that the maximum delay spread of the channel has been obtained by channel estimationAccording to an inverse Fourier transformFor frequency domain MMSE minimum mean square error channel estimationAnd (3) carrying out inverse Fourier transform:
wherein the content of the first and second substances,
frequency domain MMSE minimum mean square error channel estimation;
for the purpose of the inverse fourier transformation,
estimating a time domain expression for an MMSE minimum mean square error channel, wherein an N-point discrete sequence is arranged on a time domain and is recorded as follows:
the prior information of the time delay expansion length of the channel is utilized to exceedSetting the information response outside the length to be zero, and performing time domain truncation to obtain time domain channel estimation after time domain truncation:
,0,0,,0
Time domain channel estimation after time domain truncationBy performing fourier transform, more accurate channel estimation can be achieved:
wherein
The high-precision channel estimation is finally obtained through frequency domain equalization, LS least square method channel estimation, MMSE minimum mean square error channel estimation, inverse Fourier transform, time domain truncation and Fourier transform.
The invention has the beneficial effects that: by the scheme, the mean square error of channel estimation is greatly reduced, the precision of the channel estimation is improved, and a good effect is achieved.
Drawings
Fig. 1 is a flowchart of a channel estimation method in an LTE system according to the present invention.
Fig. 2 is a mean square error comparison diagram of the LS channel estimation, MMSE + time domain truncation channel estimation algorithm.
Detailed Description
The invention is further described with reference to the following description and embodiments in conjunction with the accompanying drawings.
As shown in fig. 1, in a channel estimation method in an LTE system, at a receiver, LS least square method channel estimation is performed in a frequency domain to obtain a preliminary channel estimation result, MMSE minimum mean square error channel estimation is performed to improve accuracy of the channel estimation, IFFT inverse transformation is performed on the frequency domain channel estimation result, signals (i.e., noise) outside a channel duration are truncated and forcibly zeroed by using time limitation of channel impulse response, energy of the noise is reduced, and finally FFT transformation is performed to return to the frequency domain, thereby improving accuracy of the channel estimation. Simulation shows that the mean square error of channel estimation is greatly reduced, the precision of channel estimation is improved, and a good effect is achieved.
As shown in fig. 1, a channel estimation method in an LTE system includes the following specific implementation steps:
frequency domain equalization
For a narrow-band system, in a channel estimation algorithm, a large number of zero elements exist in a correlation matrix, and the operation amount of time domain equalization is not large, so that the time domain equalization is feasible. However, in the LTE broadband system, the maximum bandwidth of LTE can reach 20MHz, and in this case, the channel delay is large, and the computation amount of time domain equalization is too large, so the channel estimation of the LTE system generally adopts a frequency domain equalization technology.
In general, the system equation for channel estimation, in the form of a discrete convolution, can be expressed as:
wherein the content of the first and second substances,in order to input the signal, the signal is,the response of the channel impulse is determined,in order to observe the signal(s),for input signalsLength of (d).
In matrix form, the following can be expressed:
wherein the content of the first and second substances,is composed of
Is the impulse response of the channel;
is an observed signal;
is complex gaussian noise.
According to the property of discrete Fourier transform, the convolution of time domain circles corresponds to the multiplication of frequency domains, the convolution operation of the time domain is more complex, and the conversion to the frequency domain is simpler.
To apply frequency domain multiplication, the relationship of the above expressions is first converted to a circular convolution form byCut into lengths ofThen add up:
the observation signal and noise are also treated in the same way:
the signals thus processed are related to the convolution:
expressed in a matrix as
Andis formed byAnda column vector formed by the arrangement of the column vectors,
in the matrix theory, the matrix is generally represented by upper bold and the vector is represented by lower bold, and in order to comply with the rule, the algorithm expression of the invention is expressed byThe frequency domain signal is represented so as not to be confused. If it is at handArranged as vectorsHandle barArranged as vectorsThen the discrete fourier transform can be represented by a matrix:
inverse transform to
WhereinIs a discrete fourier transform matrix of N points.
Note the book
Equation of pairsDFT conversion is carried out on two sides:
from this, the equation in the frequency domain can be expressed as:
wherein
According to the convolution theorem, the time domain circular convolution corresponds to the frequency domain product,must be a diagonal matrix:
(II) LS least squares channel estimation
LS least squares channel estimation does not take noise effects into account, so the above equation in the frequency domainThe method comprises the following steps:
constructing an error vector:
where, e is the error vector and,is a received frequency-domain signal that is,is a frequency-domain signal that is transmitted,is a frequency domain LS least squares method channel estimation,presentation pairEstimating a signal vector;
to minimize the mean square error of the signal estimate, the condition needs to be satisfied:
transforming the above equation:
the above J-pair matrixCalculating the partial derivatives, and making the partial derivatives zero, then:
namely:
therefore, the method comprises the following steps:
thereby obtaining a frequency domain LS least squares channel estimate
(III) MMSE minimum mean square error channel estimation
MMSE minimum mean square error channel estimation fully considers the influence of noise, and further carries out minimum mean square error channel estimation on the basis of completing preliminary estimation of a channel by an LS least square method
When MMSE minimum mean square error channel estimation is adopted, orderFor the above equation in the frequency domainConstructing an error vector:
where, e is the error vector and,is a characteristic of the frequency domain of the channel,is channel frequency domain MMSE minimum mean square error estimation;
to minimize the mean square error of the signal estimate, the condition needs to be satisfied:
transforming the above equation:
the above J-pair matrixCalculating the partial derivatives, and making the partial derivatives zero, then:
namely:
therefore, the method comprises the following steps:
order toRepresents the noise power, then
Therefore, there are
Wherein, in order to simplify the calculation, the actual engineering realization is considered inIn the calculation of (1) inApproximate substitutionThereby obtaining a frequency domain MMSE minimum mean square error channel estimation
(IV) noise estimation
Note the above formulaIn which there is also a noise termNeed to estimate, in order to estimate, the system noiseWe can do the following calculations.
Suppose inOn one OFDM symbolThe received signal on the subcarriers is represented as:
whereinIs the number of sub-carriers,represents the m th OFDM symbolThe channel frequency response on the sub-carriers,represents the m th OFDM symbolThe symbols transmitted on the individual sub-carriers,represents the m th OFDM symbolAdditive Gaussian stationary noise on subcarriers with mean of 0 and variance of
Since the following steps are performed for each subcarrier k, for convenience of description, the subscript k is omitted in the following equation, and the pilot symbol positions are taken as 3 and 10 as examples (uplink LTE pilot symbol positions):
1) calculating channel frequency response by LS channel estimation
2) Calculating channel frequency domain response difference;
3) Calculating a state compensation quantity;
4) ComputingWherein
5) ByCalculating
6) Computing
7) ComputingAnd thus a noise variance estimate.
(V) time domain truncation
Obtained by adopting frequency domain MMSE minimum mean square error channel estimation methodAnd then, more accurate channel estimation can be obtained by a time domain truncation processing method.
Note that there is actually also information that can be used to reduce the impact of burst noise, and in a wireless propagation environment, the delay spread of the channel, i.e. the coherence bandwidth, is substantially constant, which can be considered as a known quantity, and is further obtained by the already obtained LS channel estimation:
1) computation of autocorrelation function of frequency domain channel response from LS channel estimation
2) Determining the 3db bandwidth of an autocorrelation function
3) Calculating frequency domain correlation bandwidth
4) Determining root mean square delay spread
5) According to the Gaussian normal probability distribution, the maximum time delay expansion probability is less thanGet it
Assuming that the maximum delay spread of the channel has been obtained by channel estimationAccording to an inverse Fourier transformFor frequency domain MMSE minimum mean square error channel estimationAnd (3) carrying out inverse Fourier transform:
wherein the content of the first and second substances,
minimum mean square error channel estimation for frequency domain MMSE
Is an inverse Fourier transform
Estimating a time domain expression for an MMSE minimum mean square error channel, wherein an N-point discrete sequence is arranged on a time domain and is recorded as follows:
the prior information of the time delay expansion length of the channel is utilized to exceedSetting the information response outside the length to be zero, and performing time domain truncation to obtain time domain channel estimation after time domain truncation:
,0,0,,0
Time domain channel estimation after time domain truncationBy performing fourier transform, more accurate channel estimation can be achieved:
wherein
The high-precision channel estimation is finally obtained through frequency domain equalization, LS least square method channel estimation, MMSE minimum mean square error channel estimation, inverse Fourier transform, time domain truncation and Fourier transform.
As shown in fig. 1, it is a block diagram of a specific implementation of the method of the present invention, and the specific implementation is as follows:
(one) LS least squares channel estimation
According toAnd performing LS channel estimation, wherein the obtained LS channel estimation result is used for noise estimation, frequency domain channel response autocorrelation calculation and coefficient generation of MMSE channel estimation.
(II) noise estimation
And according to the received signal and the LS channel estimation result, applying the noise estimation method to estimate the noise variance for MMSE coefficient generation.
MMSE estimation module
And generating an MMSE coefficient for MMSE channel estimation according to the LS channel estimation result, the frequency domain channel response correlation function calculated by the LS channel estimation result and the noise variance estimation.
(IV) time domain truncation
And performing inverse Fourier transform on the frequency domain channel estimation result obtained by MMSE channel estimation to obtain time domain channel estimation.
Determining the bandwidth of the correlation function 3db according to the frequency domain channel response correlation function, and calculating the root mean square delay spreadDetermining the maximum delay spreadAnd performing time domain truncation on the obtained time domain channel estimation.
And performing Fourier transform on the time domain channel after time domain truncation to finally obtain high-precision frequency domain channel estimation.
As shown in fig. 2, the mean square error comparison of the LS channel estimation, MMSE + time domain truncation channel estimation algorithm is also consistent with the expected result, the LS channel estimation does not consider the influence of noise factors, and therefore the mean square error of the channel estimation is maximum; the MMSE channel estimation fully considers the influence of noise factors, the mean square error of the channel estimation is smaller, but the time delay expansion of the channel is not considered in the MMSE channel estimation, so the channel estimation is still not accurate enough, and the MMSE + time domain truncation channel estimation algorithm fully utilizes the time delay expansion information of the channel and carries out time domain truncation processing on the basis of the MMSE channel estimation, so the accuracy of the channel estimation is further improved, and the high-precision frequency domain channel estimation is obtained.
Simulation shows that the mean square error of channel estimation is greatly reduced, the precision of channel estimation is improved, and a good effect is achieved.
The foregoing is a more detailed description of the invention in connection with specific preferred embodiments and it is not intended that the invention be limited to these specific details. For those skilled in the art to which the invention pertains, several simple deductions or substitutions can be made without departing from the spirit of the invention, and all shall be considered as belonging to the protection scope of the invention.

Claims (7)

1. A channel estimation method under an LTE system is characterized in that: at the receiver side, firstly, LS least square method channel estimation is carried out on a frequency domain to obtain a primary channel estimation result, then MMSE minimum mean square error channel estimation is carried out to improve the accuracy of channel estimation, then IFFT inverse transformation is carried out on the frequency domain channel estimation result, signals beyond the channel duration are truncated by utilizing the time limitation of channel impulse response, zero is forcibly set, the energy of noise is reduced, and finally FFT transformation is carried out to return to the frequency domain, so that the accuracy of channel estimation is improved.
2. The channel estimation method under the LTE system according to claim 1, wherein: the channel estimation method comprises the following steps:
s1, frequency domain equalization;
s2, LS least square method channel estimation;
s3, MMSE minimum mean square error channel estimation;
s4, noise estimation;
s5, time domain truncation.
3. The channel estimation method under the LTE system according to claim 2, wherein: in step S1, the system equation for channel estimation is expressed in the form of discrete convolution:
wherein the content of the first and second substances,in order to input the signal, the signal is,the response of the channel impulse is determined,in order to observe the signal(s),for input signalsLength of (d);
expressed in matrix form as follows:
wherein the content of the first and second substances,is composed of
Is the impulse response of the channel;
is an observed signal;
is complex gaussian noise;
according to the property of discrete Fourier transform, multiplying time domain circular convolution by a frequency domain;
to apply frequency domain multiplication, the relationship of the above expressions is first converted to the form of a circular convolution by the following methodCut into lengths ofThen add up:
the observation signal and noise are also treated in the same way:
the signals thus processed are related to the convolution:
expressed in a matrix, then:
andis formed byAnda column vector formed by the arrangement of the column vectors,
within the matrix theory, the matrices are represented in upper case bold and the vectors are represented in lower case bold, and in order to comply with this rule, the expression is used hereinRepresenting frequency domain signals to avoid confusion, if at allArranged as vectorsHandle barArranged as vectorsThen the discrete fourier transform can be represented by a matrix:
inverse transformation:
whereinA discrete Fourier transform matrix of N points;
note the book
Equation of pairsDFT conversion is carried out on two sides:
from this, the equation in the frequency domain is expressed as:
wherein the content of the first and second substances,
according to the convolution theorem, the time domain circular convolution corresponds to the frequency domain product,must be a diagonal matrix:
4. the channel estimation method under the LTE system according to claim 3, wherein: in the step S2, in step S2,
LS least squares channel estimation does not take noise effects into account, so equations in the frequency domainThe method comprises the following steps:
constructing an error vector:
where, e is the error vector and,is a received frequency-domain signal that is,is a frequency-domain signal that is transmitted,
is a frequency domain LS least squares method channel estimation,presentation pairEstimating a signal vector;
to minimize the mean square error of the signal estimate, the condition needs to be satisfied:
transforming the above equation:
the above J-pair matrixCalculating the partial derivatives, and making the partial derivatives zero, then:
namely:
therefore, the method comprises the following steps:
thereby obtaining a frequency domain LS least squares channel estimate
5. The channel estimation method under LTE system according to claim 4, characterized in that: in the step S3, in step S3,
MMSE minimum mean square error channel estimation fully considers the influence of noise, further carries out minimum mean square error channel estimation on the basis of completing the preliminary estimation of the channel by an LS least square method,
when MMSE minimum mean square error channel estimation is adopted, orderFor equation in frequency domainConstructing an error vector:
where, e is the error vector and,is a characteristic of the frequency domain of the channel,is channel frequency domain MMSE minimum mean square error estimation;
to minimize the mean square error of the signal estimate, the condition needs to be satisfied:
transforming the above equation:
the above J-pair matrixCalculating the partial derivatives, and making the partial derivatives zero, then:
namely:
therefore, the method comprises the following steps:
order toRepresents the noise power, then
Therefore, there are
Wherein, in order to simplify the calculation, the actual engineering realization is considered inIn the calculation of (1) inApproximate substitutionThereby obtaining a frequency domain MMSE minimum mean square error channel estimation
6. The channel estimation method under the LTE system according to claim 5, wherein: in the step S4, in step S4,
the upper typeIn which there is also a noise termNeed to estimate, in order to estimate, the system noiseThe following calculations were performed:
suppose inOn one OFDM symbolThe received signal on the subcarriers is represented as:
whereinIs the number of sub-carriers,represents the m th OFDM symbolThe channel frequency response on the sub-carriers,represents the m th OFDM symbolThe symbols transmitted on the individual sub-carriers,represents the m < th >On OFDM symbolAdditive Gaussian stationary noise on subcarriers with mean of 0 and variance of
Since the following steps are performed for each subcarrier k, for convenience of description, the subscript k is omitted in the following equation, and the pilot symbol positions are taken as 3 and 10 as examples, that is, the uplink LTE pilot symbol positions:
1) calculating channel frequency response by LS channel estimation
2) Calculating channel frequency domain response difference;
3) Calculating a state compensation quantity;
4) ComputingWherein
5) ByCalculating
6) Computing
7) ComputingAnd thus a noise variance estimate.
7. The channel estimation method under the LTE system according to claim 6, wherein: in the step S5, in step S5,
obtained by adopting frequency domain MMSE minimum mean square error channel estimation methodThen, more accurate channel estimation is obtained through a time domain truncation processing method;
in a wireless propagation environment, the delay spread of the channel, i.e. the coherence bandwidth, is substantially constant, which can be considered as a known quantity, and is further obtained by the already obtained LS channel estimate:
1) calculating an autocorrelation function of frequency domain channel response according to LS channel estimation;
2) determining a 3db bandwidth of the autocorrelation function;
3) calculating frequency domain correlation bandwidth
4) Determining root mean square delay spread
5) According to the Gaussian normal probability distribution, the maximum time delay expansion probability is less thanGet it
Assuming that the maximum delay spread of the channel has been obtained by channel estimationAccording to an inverse Fourier transformFor frequency domain MMSE minimum mean square error channel estimationAnd (3) carrying out inverse Fourier transform:
wherein the content of the first and second substances,
frequency domain MMSE minimum mean square error channel estimation;
for the purpose of the inverse fourier transformation,
estimating a time domain expression for an MMSE minimum mean square error channel, wherein an N-point discrete sequence is arranged on a time domain and is recorded as follows:
the prior information of the time delay expansion length of the channel is utilized to exceedSetting the information response outside the length to be zero, and performing time domain truncation to obtain time domain channel estimation after time domain truncation:
,0,0,,0
Time domain channel estimation after time domain truncationBy performing fourier transform, more accurate channel estimation can be achieved:
wherein
The high-precision channel estimation is finally obtained through frequency domain equalization, LS least square method channel estimation, MMSE minimum mean square error channel estimation, inverse Fourier transform, time domain truncation and Fourier transform.
CN201910812011.2A 2019-08-30 2019-08-30 Channel estimation method under LTE system Pending CN110602006A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910812011.2A CN110602006A (en) 2019-08-30 2019-08-30 Channel estimation method under LTE system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910812011.2A CN110602006A (en) 2019-08-30 2019-08-30 Channel estimation method under LTE system

Publications (1)

Publication Number Publication Date
CN110602006A true CN110602006A (en) 2019-12-20

Family

ID=68856790

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910812011.2A Pending CN110602006A (en) 2019-08-30 2019-08-30 Channel estimation method under LTE system

Country Status (1)

Country Link
CN (1) CN110602006A (en)

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
TW200915751A (en) * 2007-09-28 2009-04-01 Univ Southern Taiwan Tech Orthogonal frequency division multiplexing system and its channel estimating appliance and method
US20120327991A1 (en) * 2010-03-04 2012-12-27 Tomasz Hrycak Method for channel estimation
CN104486266A (en) * 2014-12-12 2015-04-01 中国科学院自动化研究所 Method and device for estimating channel based on MIMO-OFDM system
CN106453165A (en) * 2016-12-05 2017-02-22 珠海全志科技股份有限公司 DFT-BASE channel estimation method and device
CN108418770A (en) * 2018-01-22 2018-08-17 南京邮电大学 Frequency domain channel reciprocity compensation method based on channel estimation errors in extensive MIMO
CN108881077A (en) * 2018-07-03 2018-11-23 重庆大学 A kind of Wiener filtering channel estimation methods based on superimposed pilot

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
TW200915751A (en) * 2007-09-28 2009-04-01 Univ Southern Taiwan Tech Orthogonal frequency division multiplexing system and its channel estimating appliance and method
US20120327991A1 (en) * 2010-03-04 2012-12-27 Tomasz Hrycak Method for channel estimation
CN104486266A (en) * 2014-12-12 2015-04-01 中国科学院自动化研究所 Method and device for estimating channel based on MIMO-OFDM system
CN106453165A (en) * 2016-12-05 2017-02-22 珠海全志科技股份有限公司 DFT-BASE channel estimation method and device
CN108418770A (en) * 2018-01-22 2018-08-17 南京邮电大学 Frequency domain channel reciprocity compensation method based on channel estimation errors in extensive MIMO
CN108881077A (en) * 2018-07-03 2018-11-23 重庆大学 A kind of Wiener filtering channel estimation methods based on superimposed pilot

Similar Documents

Publication Publication Date Title
US20040005010A1 (en) Channel estimator and equalizer for OFDM systems
CN103107969B (en) Incremental iterative time-varying channel evaluation and inter carrier interference (ICI) elimination method of fast orthogonal frequency division multiplexing (OFDM) system
CN101064571B (en) Apparatus for enhancing channel evaluation in OFDM receiver and its method
Ma et al. A low complexity MMSE for OFDM systems over frequency-selective fading channels
US8259864B2 (en) Methods and systems for fourier-quadratic basis channel estimation in OFDMA systems
CN107332797B (en) Channel estimation method in power line OFDM communication system
US8824534B2 (en) Method for channel estimation
CN110677361B (en) Signal equalization method, equalizer and storage medium for orthogonal time-frequency space system
CN101808053A (en) Channel estimation method and device based on OFDM (Orthogonal Frequency Division Multiplexing)
CN107222438B (en) Simplified BEM channel estimation method of high-speed mobile SC-FDMA system
CN113890796A (en) High-speed channel estimation device and method based on OTFS system modulation and demodulation
CN103873422A (en) Method for eliminating multipath interference in system symbol in underwater sound orthogonal frequency-division multiplexing system
Wu et al. DFT-spread orthogonal time frequency space system with superimposed pilots for terahertz integrated sensing and communication
CN103873406A (en) Inter-frame interference elimination method used for underwater sound orthogonal frequency-division multiplexing communication system
Rottenberg et al. Efficient equalization of time-varying channels in MIMO OFDM systems
EP2615784B1 (en) A method of equalizing an OFDM signal
CN102780656A (en) Method and device for eliminating multi-symbol subcarrier jamming and performing channel estimation jointly
CN102984113B (en) Fast time varying channel equalization method of OFDM system based on time varying
CN102113285A (en) A simplified equalizationscheme for distributed resource allocation in multi-carrier systems
JP6220844B2 (en) MIMO system testing apparatus and testing method
Yu et al. Iterative estimation of doubly selective underwater acoustic channel using basis expansion models
CN115412416B (en) Low-complexity OTFS signal detection method for high-speed moving scene
CN1984109A (en) Channel estimater and channel estimating method in telecommunication system
CN110602006A (en) Channel estimation method under LTE system
KR20200005807A (en) Method and apparatus for performing blind selected mapping

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20191220

WD01 Invention patent application deemed withdrawn after publication