Detailed Description
In order to ensure that the aviation broadband OFDM mobile communication system can well suppress noise in different channel environments, it is necessary to dynamically track changes in the channel impulse response length. The invention provides a method for estimating the channel impulse response length, which can adaptively estimate the size of the channel impulse response length in real time, thereby ensuring the channel estimation performance under the complicated and changeable aviation broadband channel condition.
Fig. 1 is a flowchart illustrating a method for estimating a channel impulse response length according to an embodiment of the present invention, as shown in fig. 1:
101. and performing Least Square (LS) channel estimation on the frequency domain data after synchronization and frequency offset compensation through a Least Square (LS) criterion at the pilot frequency position to obtain channel frequency domain response at the pilot frequency position.
After receiving data, a receiver needs to synchronize first, and after synchronization is successful, further frequency offset estimation and compensation are generally needed to be performed on signals. For example, when the base station transmits data to the user terminal, the receiver is a receiver on the user terminal side, and when the user terminal transmits data to the base station, the receiver is a receiver on the base station side.
Assuming that the frequency domain data matrix at the pilot frequency position after synchronization and frequency offset compensation is Y (diagonal matrix), the frequency domain data matrix at the pilot frequency position of the transmitting end stored locally is X (diagonal matrix), HLSUsing the least squares criterion for the channel frequency domain response truth matrix at the pilot under the LS criterion, we obtain:
minf(HLS)=min{(Y-XHLS)H(Y-XHLS)}
from the calculus theorem, the function minf (H) isLS) Taking the minimum value, it is necessary to make the function pair HLSHas a partial derivative of 0, i.e.Obtaining an estimated value matrix of channel frequency domain response at the pilot frequency under the LS criterion:
102. and interpolating the channel frequency domain response at the pilot frequency to obtain the channel frequency domain response at the data.
It should be noted that the interpolation algorithm may be the prior art, and generally, the algorithm used for interpolation may select different interpolation algorithms according to the system requirements and the pilot pattern used by the system, and it is assumed that the present embodiment uses a first-order linear interpolation algorithm, where the first-order linear interpolation algorithm uses two adjacent pilot symbols to perform linear interpolation.
Assuming that the interpolation in the time axis direction of the two-dimensional pilot pattern is performed first, the interpolation formula is:
wherein k represents the kth subcarrier, i represents the ith OFDM symbol, p and q are intermediate variables, have no specific physical significance and are convenient for expression of a formula, and p is more than or equal to 1 and is less than or equal to deltat-1,qΔt<i<(q+1)Δt,q=0,1,2....;ΔtThe distance between two adjacent pilots on the time axis.
If linear interpolation is needed in the direction of the frequency axis, the interpolation formula is as follows:
wherein k represents the kth subcarrier, i represents the ith OFDM symbol, p and q are intermediate variables, have no specific physical significance and are convenient for expression of a formula, and q is deltaf<k<(q+1)Δf,1≤p≤Δf-1,ΔfIs the spacing between two adjacent pilots on the frequency axis.
103. And estimating the channel impulse response length of the channel frequency domain response according to the channel frequency domain response at the pilot frequency and the time domain channel response corresponding to the channel frequency domain response at the data.
Assuming that, under an ideal condition, the channel impulse response coefficient is within the channel impulse response length, the gaussian white noise is uniformly distributed in all the channel samples, and the true value of the channel time domain impulse response of the ith OFDM symbol on the kth subcarrier is:
step 103 may specifically be:
1031. the cumulative mean square error function of the channel impulse response is derived from the definition of the cumulative mean square error.
Wherein the cumulative mean square error function is:
wherein,is the estimated value of the channel time domain impulse response of the ith OFDM symbol on the kth subcarrier, and contains noise components.
The cumulative mean square error function can further be derived from equation (a): the function is a concave function and has a minimum value, and the value L corresponding to the minimum value of the cumulative mean square error function e (L) is determined as the channel impulse response length of the ith OFDM symbol;
wherein,for estimated noise error, hk,iIs the channel time domain response true value at the data,the estimated value of the time domain response of the channel at the data position is L, the length of the impulse response of the channel is expressed, N represents the total number of subcarriers, k is more than or equal to 0 and less than or equal to N, M represents the number of OFDM symbols of the orthogonal frequency division multiplexing technology, i represents the ith OFDM symbol of the orthogonal frequency division multiplexing technology, and i is more than or equal to 1 and less than or equal to M;
therefore, in a scene that the channel environment is changeable in the airborne mobile communication system, even if the channel frequency domain response at the pilot frequency and the channel frequency domain response at the data change with the change of the channel environment, the embodiment can change the estimated channel impulse response length in real time according to the time domain channel response corresponding to the channel frequency domain response at the pilot frequency and the channel frequency domain response at the data along with the change of the channel environment.
1032. And taking the maximum value of the channel impulse response length of the M OFDM symbols as the channel impulse response length of the channel frequency domain response.
In this embodiment, the length of the channel impulse response of M OFDM symbols is maximized, so that the more the number of zero forcing elements included in the time domain channel response is, the more the noise is removed, and the more accurate the channel estimation is.
Fig. 2 is a schematic flow chart of a channel estimation method of a civil aviation air-ground broadband wireless communication system applied in the embodiment of the present invention, as shown in fig. 2:
and performing Least Square (LS) channel estimation on the frequency domain data after synchronization and frequency offset compensation through a Least Square (LS) criterion at the pilot frequency position to obtain channel frequency domain response at the pilot frequency position. And interpolating the channel frequency domain response at the pilot frequency to obtain the channel frequency domain response at the data.
Assuming that, under an ideal condition, the channel impulse response coefficient is within the channel impulse response length, and the gaussian white noise is uniformly distributed in all the channel samples, the channel frequency domain impulse response of the ith OFDM symbol on the kth subcarrier is:
j mathematically represents the imaginary part of the complex number;
rewriting the formula into
Wherein N represents the total number of subcarriers, and k is more than or equal to 0 and less than or equal to N; l denotes a channel impulse response length.
The frequency domain response vector of each OFDM symbol isWherein the matrix F is a Fourier transform matrix, NDFor the number of useful sub-carriers, N is not less than NDWhen N includes useful subcarriers and virtual subcarriers, N is greater than ND,. When only useful subcarriers are contained in N, NDAnd N are both equal.
Carrying out QR decomposition on the matrix F to obtain a matrix:
wherein,is a normalized matrix of the number of bits in the matrix,is an upper triangular matrix, where NDThe elements of the L rows are zero.
Carrying out QR decomposition on the channel frequency domain response H of the data obtained by interpolation in the step 102 to obtainTo: further written as:wherein, F is a Fourier transform matrix, and H is a time domain response matrix corresponding to the frequency domain response H.
In mathematical theory, there is H = Fh, but in practical applications, there is noise in the channel estimation, and therefore, it is necessary to superimpose the noise(ii) a Where F = QR is a mathematical matrix decomposition method, called QR decomposition of a matrix.
The transposition of channel frequency domain response H multiplied by Q, namely the left and right sides of the expression are respectively multiplied by QHTo obtain a vector
Since R is an upper triangular matrix, there is ND-L lines of zero elements, it can be seen that the first item on the right side has only the first L elements, the last NDThe L elements are zero and the second term on the right of the equation is the noise term.
Suppose thatWherein the first L elements mainly comprise channel impulse response and the last N elementsD-L elements contain mainlyNoise. Therefore, the first L elements are useful data needed by the user, and in order to further reduce the influence of noise on the estimation performance, the last N elements can be usedDL noise elements zero forcing yields:
the final channel frequency response is:
the embodiment of the invention obtains the channel frequency domain response of the pilot frequency position by carrying out LS channel estimation on the frequency domain data after synchronization and frequency offset compensation through the least square criterion of the pilot frequency position; interpolating the channel frequency domain response at the pilot frequency position to obtain the channel frequency domain response at the data position; and estimating the channel impulse response length of the channel frequency domain response according to the channel frequency domain response at the pilot frequency and the time domain channel response corresponding to the channel frequency domain response at the data. Thereby enabling: in a scene that the channel environment is changeable in the airborne mobile communication system, the embodiment can estimate the channel impulse response length according to the channel frequency domain response at the pilot frequency and the time domain channel response corresponding to the channel frequency domain response at the data, and the channel impulse response length changes in real time along with the change of the channel environment, so that the method can be applied to the airborne broadband wireless communication system with complicated and changeable channel environment and can improve the channel estimation performance.
Fig. 3 is a schematic structural diagram of an apparatus for estimating a channel impulse response length according to another embodiment of the present invention, which is applied to a channel estimation scenario of a civil aviation air-ground broadband wireless communication system, and as shown in fig. 3, the apparatus includes:
the first processing module 31 is configured to perform Least Square (LS) channel estimation on the frequency domain data after synchronization and frequency offset compensation through a pilot frequency to obtain a channel frequency domain response at the pilot frequency;
a second processing module 32, configured to interpolate the channel frequency domain response at the pilot frequency to obtain a channel frequency domain response at the data;
a third processing module 33, configured to estimate a channel impulse response length of the channel frequency domain response according to the channel frequency domain response at the pilot frequency and the time domain channel response corresponding to the channel frequency domain response at the data.
For example, the third processing module 33 is specifically configured to:
using the cumulative mean square error function: determining the value L corresponding to the minimum value of the cumulative mean square error function e (L) as the channel impulse response length of the ith OFDM symbol;
wherein,for estimated noise error, hk,iIs the channel time domain response true value at the data,the estimated value of the time domain response of the channel at the data position is L, the length of the impulse response of the channel is expressed, N represents the total number of subcarriers, k is more than or equal to 0 and less than or equal to N, M represents the number of OFDM symbols of the orthogonal frequency division multiplexing technology, i represents the number of OFDM symbols of the ith orthogonal frequency division multiplexing technology, and i is more than or equal to 1 and less than or equal to M;
and taking the maximum value of the channel impulse response length of the M OFDM symbols as the channel impulse response length of the channel frequency domain response.
The embodiment of the invention obtains the channel frequency domain response of the pilot frequency position by carrying out LS channel estimation on the frequency domain data after synchronization and frequency offset compensation through the least square criterion of the pilot frequency position; interpolating the channel frequency domain response at the pilot frequency position to obtain the channel frequency domain response at the data position; and estimating the channel impulse response length of the channel frequency domain response according to the channel frequency domain response at the pilot frequency and the time domain channel response corresponding to the channel frequency domain response at the data. Thereby enabling: in a scenario where the channel environment is changeable in the airborne mobile communication system, even if the channel frequency domain response at the pilot frequency and the channel frequency domain response at the data change with the change of the channel environment, the present embodiment may change the estimated channel impulse response length in real time according to the time domain channel response corresponding to the channel frequency domain response at the pilot frequency and the channel frequency domain response at the data according to the change of the channel environment, and therefore, may be applicable to an airborne broadband wireless communication system where the channel environment is complicated and changeable, and may improve the channel estimation performance in a scenario where the channel environment is changeable in the airborne mobile communication system, the present embodiment may estimate the channel impulse response length according to the channel frequency domain response at the pilot frequency and the time domain channel response corresponding to the channel frequency domain response at the data according to the change of the channel environment, and the channel impulse response length changes in real time with the change of the channel environment, therefore, the method can be applied to aviation broadband wireless communication systems with complex and variable channel environments, and can improve the channel estimation performance.
Finally, it should be noted that: the above examples are only intended to illustrate the technical solution of the present invention, but not to limit it; although the present invention has been described in detail with reference to the foregoing embodiments, it will be understood by those of ordinary skill in the art that: the technical solutions described in the foregoing embodiments may still be modified, or some technical features may be equivalently replaced; and the modifications or the substitutions do not make the essence of the corresponding technical solutions depart from the scope of the technical solutions of the embodiments of the present invention.