Angular domain channel estimation method based on matrix reconstruction for symmetrical nonuniform array

Abstract

An angular domain channel estimation method based on matrix reconstruction for a symmetrical nonuniform array is a combined two-stage channel estimation and channel equalization scheme provided based on an SNLA model. In a first stage, a matrix reconstruction method is used to estimate a path AOA, and compared with traditional channel estimation based on ULA, the matrix reconstruction method achieves a higher resolution ratio. In a second stage, an LS method is used to obtain a path gain. According to the angular domain channel estimation method, a mean square error of channel estimation, a bit error of data transmission and complexity of a traditional scheme are significantly reduced. A simulation result indicates that compared with the traditional method, the angular domain channel estimation method can achieve a lower MSE and BER.

Claims

1. An angular domain channel estimation method based on matrix reconstruction for a symmetrical nonuniform array, defining that a system comprises a base station with M antennas and a user with a single antenna, wherein the M antennas form a symmetrical nonuniform linear array and the symmetrical nonuniform linear array is divided into a dense symmetrical uniform linear subarray, a first sparse uniform linear subarray and a second sparse uniform linear subarray; the dense symmetrical uniform linear subarray has 2M.sub.1+1 array elements; each of the array elements has a spacing d, and d=/2, being a half of a wavelength; each of the first sparse uniform linear subarray and the second sparse uniform linear subarray comprises M.sub.2 array elements, each of the array elements has a spacing (M.sub.1+1)d, and M=2(M.sub.1+M.sub.2)+1; the first sparse uniform linear subarray and the second sparse uniform linear subarray are respectively deployed on two sides of the dense symmetrical uniform linear subarray; an array element in a middle of the dense symmetrical uniform linear subarray is selected as a reference array element, and rest array elements are symmetrically distributed by taking the reference array element as a center; since a wireless channel experiences limited scattering propagation, the channel has a sparse multi-path structure and the user is defined to have L scattering paths; and the channel estimation method comprises: performing path angle estimation based on a matrix reconstruction method, specifically as follows: enabling a user side to send a training signal s t at a time t, and in all snapshots, enabling |s.sub.t|=1, then a receiving signal at a position of a base station antenna being:
y.sub.t=h.sub.ts.sub.t+n.sub.t=Ag.sub.ts.sub.t+n.sub.t wherein h.sub.t is a user uplink channel, n.sub.t is an additive white Gaussian noise obeying complex Gaussian distribution CN(0,.sup.2I), .sup.2I is a variance of the additive white Gaussian noise, and Ag.sub.t is a form of matrix multiplication of the channel h.sub.t:
g.sub.t=[g.sub.1,t, . . . ,g.sub.l,t].sup.TC.sup.L1
A=[a(.sub.1), . . . ,a(.sub.L))]C.sup.ML g.sub.l,t being a channel gain of the user at the time t and at an l.sup.th scattering path and obeying complex Gaussian distribution g.sub.i,tCN(0,1), .sub.l representing an angle-of-arrivals of an l.sup.th path of the user, and a vector a(.sub.l)C.sup.M1 representing an array manifold vector, l=1, . . . L; enabling x.sub.t=g.sub.ts.sub.tC.sup.L1 to obtain a receiving signal covariance matrix:
R.sub.y=Ey.sub.ty.sub.t.sup.H=R.sub.h+.sup.2I=AR.sub.xA.sup.H+.sup.2I wherein R.sub.h=Eh.sub.th.sub.t.sup.H and R.sub.x=Ex.sub.tx.sub.t.sup.H; vectorizing the covariance matrix R.sub.y to obtain a vector z:
z=vec(R.sub.y)=p+.sub.m.sup.2 wherein =A*AC.sup.|M|.sup.2.sup.L, p=[g.sub.1.sup.2.sub.1.sup.2, . . . , g.sub.L.sup.2.sub.L.sup.2].sup.T, g.sub.l.sup.2 and .sub.l.sup.2 respectively represent a transmission signal power and a path gain power, 1lL, .sub.m.sup.2 is a noise power, =[e.sub.1.sup.T, e.sub.2.sup.T, . . . , e.sub.M.sup.T].sup.T, e.sub.i is a column vector, except that the i.sup.th position is 1, the rest are 0, the vector z is equivalent to receiving data with an array manifold matrix (A*A), and array element positions of the vector z are given by a set D=d.sub.id.sub.j, i,j=1, 2, . . . , M, d.sub.i represents a distance from an i.sup.th array element to the reference array element, repeated elements in the set D are deleted to obtain a set B, integer elements of the set B correspond to positions of virtual array elements, the repeated data in the receiving data z are removed and corresponding rows are rearranged to cause the rows to correspond to the positions of the virtual array to obtain a new vector:
{tilde over (z)}=A.sub.Bp+.sub.m.sup.2e.sub.0 wherein {tilde over (z)}C.sup.|B|1 is a receiving signal of the virtual array, and A.sub.BC.sup.|B|L is an array manifold matrix corresponding to the virtual array, |B|=M+2(M.sub.1+(M.sub.1+1)M.sub.2), e.sub.0C.sup.|B|1, and except that a central term is 1, the rest are 0; reconstructing the received data {tilde over (z)} into a covariance matrix ${\tilde{R}}_y{}^{\left(\frac{.Math.B.Math.+1}{2}\right)\left(\frac{.Math.B.Math.+1}{2}\right)}$ of the virtual array, the matrix {tilde over (R)}.sub.y having a toeplitz matrix property, that is, elements on the same diagonal line being the same, so during construction of the matrix {tilde over (R)}.sub.y, only constructing data in a first column and a first row, constructing previous $\frac{.Math.B.Math.+1}{2}$ data in the vector {tilde over (z)} into a first column of the matrix {tilde over (R)}.sub.y, constructing last $\frac{.Math.B.Math.-1}{2}$ data in me vector {tilde over (z)} into a first row of the matrix {tilde over (R)}.sub.y, and then complementing {tilde over (R)}.sub.y according to the property that the elements on the same diagonal line of {tilde over (R)}.sub.y are the same; based on feature value decomposition of {tilde over (R)}.sub.y, representing: ${\tilde{R}}_y=\left[\begin{array}{cc}{U}_S& U_N\end{array}\right]\left[\begin{array}{cc}{.Math.}_S& 0\\ 0& {.Math.}_N\end{array}\right]\left[\begin{array}{c}{U}_S^H\\ U_N^H\end{array}\right]$ wherein U.sub.S is a signal subspace formed by a feature vector corresponding to a large feature value, and U.sub.N is a noise subspace formed by a feature vector corresponding to a small feature value; multiplying both sides of the matrix by U.sub.N to obtain:
{tilde over (R)}.sub.yU.sub.N=(A.sub.1R.sub.xA.sub.1.sup.H+.sub.m.sup.2I)U.sub.N=.sub.m.sup.2U.sub.n, wherein A.sub.1C.sup.(|B|+1)/2L represents an array manifold matrix corresponding to the virtual array and meets:
A.sub.1R.sub.xA.sub.1.sup.HU.sub.N=0 since a column vector of A.sub.1 corresponds to a signal transmitting direction, a direction of a signal source is estimated by the characteristic; based on a multiple signal classification algorithm, defining a spatial spectrum signal P.sub.music({circumflex over ()}.sub.1) as: $P_{\mathrm{music}}\left({\stackrel{}{}}_l\right)=\frac{1}{{\tilde{a}\left({}_l\right)}^H{U}_N{U}_N^H\stackrel{}{a}\left({}_l\right)}$ wherein in a case that a denominator ({circumflex over ()}.sub.l).sup.HU.sub.NU.sub.N.sup.H({circumflex over ()}.sub.l) reaches a minimum value, ({circumflex over ()}.sub.l) is an l.sup.th column vector of the matrix A.sub.1, P.sub.music({circumflex over ()}.sub.l) reaches a maximum value, a direction-of-arrivals .sub.l is estimated according to a peak value of P.sub.music({circumflex over ()}.sub.l), thereby obtaining all path angle information {circumflex over ()}=[{circumflex over ()}.sub.1, {circumflex over ()}.sub.2, . . . , {circumflex over ()}.sub.L]; and performing path gain estimation, specifically as follows: obtaining an array manifold matrix based on the obtained {tilde over ()}, sending a pilot signal u.sub.t, estimating path gains in different time blocks based on the obtained , and constructing a cost function: $J\left({\stackrel{}{g}}_t\right)={.Math.\frac{1}{u_t}{y}_t-\hat{A}{\hat{g}}_t.Math.}^2={\left(\frac{1}{u_t}{y}_t-\hat{A}{\hat{g}}_t\right)}^H\left(\frac{1}{u_t}{y}_t-\hat{A}{\hat{g}}_t\right)={\left(\frac{1}{u_t}\right)}^2{y}_t^H{y}_t-\frac{1}{u_t}{y}_t^H\hat{A}{\hat{g}}_t-\frac{1}{u_t}{\stackrel{}{g}}_t^H{\hat{A}}^H{y}_t+{\stackrel{}{g}}_t^H{\hat{A}}^H\hat{A}{\hat{g}}_t$ minimizing the cost function to obtain a channel gain estimation .sub.t, specifically, by calculating a partial derivative of the cost function relative to .sub.t, obtaining: $\frac{J\left({\stackrel{}{g}}_t\right)}{{\stackrel{}{g}}_t}=-\frac{1}{u_t}{y}_t^H\hat{A}+{\hat{g}}_t^H{\hat{A}}^H\hat{A}$ in a case that ${\stackrel{}{g}}_t^H{\hat{A}}^H\hat{A}=\frac{1}{u_t}{y}_t^H\hat{A},$ a solution of the channel gain being: ${\stackrel{}{g}}_t=\frac{1}{u_t}{\left({\hat{A}}^H\hat{A}\right)}^{-1}{\hat{A}}^H{y}_t,t=1,.Math.,T$ wherein T is a time block, then within one time block, the whole channel estimation result expression is:
.sub.t=.sub.t,t=1, . . . ,T.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

(1) FIG. 1 is a sparse channel estimation system model based on SNLA;

(2) FIG. 2 is a symmetrical nonuniform array geometric structure;

(3) FIG. 3 is a two-stage channel estimation signal processing process;

(4) FIG. 4 is a spatially normalized frequency spectrum of SNLA and ULA;

(5) FIG. 5 is an NMSE and SNR curve graph of LRSCR, SOMP, MUSIC and MR algorithms of a single-user system; and

(6) FIG. 6 is a schematic structural diagram of comparing with an ideal CSI, an SNLA CSI estimated by the provided method, a ULA CSI estimated by an MUSIC method and a bit error ratio estimated by an LS method, respectively.

DETAILED DESCRIPTION OF THE EMBODIMENTS

(7) The technical solution of the present invention will be clearly and completely described below with reference to the accompanying drawings and simulation examples, so that those skilled in the art can better understand the present invention.

(8) The present invention considers a single-user communication system model, as shown in FIG. 1. The communication system is composed of a M-antenna base station and a single-antenna user, where M=2(M.sub.1+M.sub.2)+1. At the BS side, an SNLA is designed and includes a dense symmetrical uniform linear subarray 1 and two sparse uniform linear subarrays, where the subarray 2 and the subarray 3 are respectively located on a left side and a right side of the subarray 1, as shown in FIG. 2. The dense subarray 1 totally has 2M.sub.1+1 elements, and each array has a spacing d, where d=/2, and is a half of a wavelength. Both the sparse subarray 2 and the sparse subarray 3 include M.sub.2 elements, and each array has a spacing (M.sub.1+1)d. For channel estimation, a matrix reconstruction (MR) is provided to estimate a channel path angle based on a Vandermonde structure of SNLA.

(9) Since a wireless channel will experience limited scattering propagation to cause the channel to have a sparse multi-path structure, it is assumed that the user has L scattering paths. Then, the channel may be described by a geometric model with L(L<M) scattering bodies, where each path is represented by a path angle and a path gain. In this channel modeling, the angle of the scattering path is kept unchanged within a relatively long time, but the channel coefficient is changed very quickly, so the user uplink channel h.sub.t may be represented as:

(10) 0$\begin{array}{cc}{h}_t={.Math.}_{l=1}^L{g}_{l,t}a\left({}_l\right),t,k=1,.Math.,K& \left(1\right)\end{array}$ where g.sub.l,t is defined as the channel gain of the user at the time t and at an l.sup.th scattering path, obeying complex Gaussian distribution g.sub.l,tCN(0,1). .sub.l represents the (DOA) of the l.sup.th path, and the vector a(0.sub.l)C.sup.M1 represents the array manifold vector, having the following forms:
a(.sub.l)=[e.sup.jd.sup.1, . . . ,e.sup.js.sup.M].sup.T(2) where

(11) $=\frac{2}{}\sin \left({}_{k,l}\right),d_i\left(1iM\right)$
represents the distance from the i.sup.th array element to a reference element. It can be clearly seen from FIG. 2 that, from left to right, the index information of the array is respectively M.sub.1M.sub.2, . . . , M.sub.11, M.sub.1, . . . , 0, . . . , M.sub.1, M.sub.1+1, . . . , M.sub.1+M.sub.2, the array is placed at the position d.sub.i, then d.sub.i=(M.sub.1M.sub.2(M.sub.1+1))d, . . . , (2M.sub.1+1)d, M.sub.1d, . . . , 0, . . . , M.sub.1d, (2M.sub.1+1)d, . . . , (M.sub.1+M.sub.2(M.sub.1+1))d. The channel is represented as a form of matrix multiplication, and the expression h.sub.t in (1) may be represented as:
h.sub.t=Ag.sub.t(3) where g.sub.t=[g.sub.1,t, . . . , g.sub.L,t].sup.TC.sup.L1, A=[a(.sub.1), . . . , a(.sub.L))]C.sup.ML.

(12) In the uplink channel estimation, the user side sends a training signal s.sub.t, and ensures |s.sub.t|=1 in all snapshots. The receiving signal at the position of the base station antenna may be represented as:
y.sub.t=h.sub.ts.sub.t+n.sub.t=Ag.sub.ts.sub.t+n.sub.t(4) where n.sub.t is defined as an additive white Gaussian noise obeying complex Gaussian distribution CN(0,.sup.2I). According to (4), at the time t, the covariance matrix of the receiving signal may be represented as:
R.sub.y=E{y.sub.ty.sub.t.sup.H}=R.sub.h+.sup.2I=AR.sub.xA.sup.H+.sup.2I(5) where R.sub.h=E{h.sub.th.sub.t.sup.H}, R.sub.x=E{x.sub.tx.sub.t.sup.H}. According to (5), the channel path angle information may be estimated based on the MR method. It is assumed that in a block fading channel, the path angle is changed slowly and is kept constant in a block, but the path gain is changed very quickly. Therefore, the present invention designs a signal processing framework of two-stage channel estimation considering the difference of the path angle and the path gain changed with time, as shown in FIG. 3. At a first stage, the channel covariance matrix R.sub.y, is Obtained by the array signal processing related method, and the path angle information is retrieved based on the MR method. During the estimation of the channel angle information, it is only necessary to know the covariance matrix of the receiving signal. Therefore, at the first stage, the channel angle information can be acquired without sending the pilot signal. At a second stage, the path gain is obtained by the least square method by sending pilot information in different blocks. Therefore, in signal processing, CSI estimation is divided into two sub-problems: path angle estimation and path gain estimation.

(13) The present invention converts the channel estimation into the DOA estimation problem and the path gain estimation problem. In more detail, first, the MR method based on the SNAL structure is proposed to estimate DOA. Then, the channel path gain is Obtained by an LS method. The M-dimensional channel covariance matrix may be expanded as (M+2(M.sub.1+(M.sub.1+1)M.sub.2)+1)/2-dimensional virtual covariance matrix by the MR method, thereby increasing the degree of spatial freedom and improving the accuracy of the channel estimation.

(14) The path angle estimation based on the matrix reconstruction is: according to the receiving signal expression in the formula (4), x.sub.t=g.sub.ts.sub.tC.sup.L1, the receiving signal may be further represented as:
y.sub.t=Ag.sub.ts.sub.t+n.sub.t=Ax.sub.t+n.sub.t(6)

(15) Since the channel path gain information is changed quickly within a block time, it is assumed that the channel has L paths, and after N snapshots, XC.sup.LT may be represented as:

(16) $\begin{array}{cc}X=\left(\begin{array}{ccc}{g}_{11}{s}_1& .Math.& g_{1N}{s}_N\\ g_{21}{s}_1& .Math.& g_{21}{s}_N\\ .Math.& & .Math.\\ g_{L1}{s}_1& .Math.& g_{L1}{s}_N\end{array}\right)& \left(7\right)\end{array}$

(17) In the formula (7), each row of X is independent and irrelevant. Therefore, the autocorrelation matrix

(18) of x.sub.t is a diagonal matrix. Q.sub.i=g.sub.i.sup.2.sub.i.sup.2(1iL) is defined, where g.sub.i.sup.2 and .sub.i.sup.2(1lL) respectively represent a transmission signal power and a path gain power, and the obtained covariance matrix of the receiving signal is:

(19) $\begin{array}{cc}{R}_y={\mathrm{AR}}_x{A}^H+{}_m^{2}I=A\left(\begin{array}{cccc}{Q}_1& & & \\ & Q_2& & \\ & & & \\ & & & Q_L\end{array}\right)A^H+{}_m^{2}I& \left(8\right)\end{array}$

(20) The element in the covariance matrix [R.sub.y].sub.m,n(1mM, 1nM) may be regarded as the receiving data at the position coordinate d.sub.md.sub.n array element, so one array element can be virtually created at the position where no array element is originally through the known physical array element. The covariance matrix R.sub.y is vectorized to obtain the following vector:
z=vec(R.sub.y)=vec(AR.sub.xA)+.sub.m.sup.2vec(I)=(A*A)p+.sub.m.sup.2vec(I)=p+.sub.m.sup.2(9) where =A*AC.sup.|M|.sup.2.sup.L, p=[g.sub.1.sup.2.sub.1.sup.2, . . . , g.sub.L.sup.2.sub.L.sup.2].sup.T, .sub.m.sup.2 is a noise power, =[e.sub.1.sup.T, e.sub.2.sup.T, . . . , e.sub.M.sup.T].sup.T, e.sub.i is a column vector, and except that an i.sup.th position is 1, the rest are 0. Similar to the formula (6), the vector z is equivalent to the receiving data with the array manifold matrix (A*A) and corresponds to a larger array, and the array element position is given by a set D={d.sub.id.sub.j}, i,j=1,2, . . . , M.

(21) Since the difference value will be the same when the difference value between any pair of original physical array elements is calculated, that is, the same array element is virtually created, so the vector z is redundant. The repeated elements in the set D are deleted to obtain a set B, the integer elements of the set B correspond to the position of the virtual array element, the repeated data in the receiving data z is removed, and the corresponding rows are rearranged to cause the row to correspond to the virtual array position to obtain a new vector:
{tilde over (z)}=A.sub.Bp+.sub.m.sup.2e.sub.0(10) where {tilde over (z)}C.sup.|B|1 is a receiving signal of the virtual array, and A.sub.BC.sup.|B|L is an array manifold matrix corresponding to the virtual array, |B|=M+2(M.sub.1+(M.sub.1+1)M.sub.2), e.sub.0C.sup.|B|1, and except that a central term is 1, the rest are 0.

(22) The received data {tilde over (z)} is reconstructed into the covariance matrix

(23) ${\tilde{R}}_y{}^{\left(\frac{.Math.B.Math.+1}{2}\right)\left(\frac{.Math.B.Math.+1}{2}\right)}$
of the virtual array, where the matrix {tilde over (R)}.sub.y has the toeplitz matrix property, that is, the elements on the same diagonal line are the same. Therefore, during the construction of the matrix {tilde over (R)}.sub.y, only data in the first column and the first row is required to be constructed; constructing the previous

(24) $\frac{.Math.B.Math.+1}{2}$
data in me vector {tilde over (z)} into me first column of the matrix {tilde over (R)}.sub.y, the last

(25) $\frac{.Math.B.Math.-1}{2}$
data in the vector {tilde over (z)} is constructed into the first row of the matrix {tilde over (R)}.sub.y, and then {tilde over (R)}.sub.y is completed according to the property that the elements on the same diagonal line of {tilde over (R)}.sub.y are the same. Based on feature value decomposition of {tilde over (R)}.sub.y, it may be represented as:

(26) $\begin{array}{cc}{\tilde{R}}_y=\left[U_S{U}_N\right]\left[\begin{array}{cc}{.Math.}_S& 0\\ 0& {.Math.}_N\end{array}\right]\left[\begin{array}{c}{U}_S^H\\ U_N^H\end{array}\right]& \left(11\right)\end{array}$ where U.sub.S is a signal subspace formed by a feature vector corresponding to a large feature value and U.sub.N is a noise subspace formed by a feature vector corresponding to a small feature value. Both sides of the matrix are multiplied by U.sub.N to obtain:
{tilde over (R)}.sub.yU.sub.N=(A.sub.1R.sub.xA.sub.1.sup.H+.sub.m.sup.2I)U.sub.N=.sub.m.sup.2U.sub.n(12) where A.sub.1C.sup.(|B|+1)/2L represents an array manifold matrix corresponding to the virtual array and meets:
A.sub.1R.sub.xA.sub.1.sup.HU.sub.N=0(13)

(27) Since the column vector of A.sub.1 corresponds to a signal transmitting direction, a direction of a signal source may be estimated by the characteristic. Due to the influence of the noise, the general signal subspace and the noise subspace cannot be completely orthogonal, and based on a multiple signal classification (MUSIC) algorithm, a spatial spectrum signal P.sub.music({circumflex over ()}.sub.l) is defined as:

(28) $\begin{array}{cc}{P}_{\mathrm{music}}\left({\hat{}}_l\right)=\frac{1}{{\tilde{a}\left({}_l\right)}^H{U}_N{U}_N^H\tilde{a}\left({}_l\right)}& \left(14\right)\end{array}$ where in a case that the denominator ({circumflex over ()}.sub.l).sup.HU.sub.NU.sub.N.sup.H({circumflex over ()}.sub.l) reaches a minimum value, ({circumflex over ()}.sub.l) is the l.sup.th column of vector of the matrix A.sub.1, P.sub.music({circumflex over ()}.sub.l) reaches a maximum value, the direction-of-arrivals .sub.l may be estimated according to the peak valine of P.sub.music({circumflex over ()}.sub.l), thereby obtaining all path angle information {circumflex over ()}=[{circumflex over ()}.sub.1, {circumflex over ()}.sub.2, . . . , {circumflex over ()}.sub.L].

(29) Based on the expression R.sub.y in the formula (8), during the estimation of the path angle information, it is only necessary to know the signal sending statistical information. Therefore, the first stage of the channel estimation does not require the pilot signal, thereby greatly reducing the pilot overhead. In addition, for the ULA with M antennas, the maximum antenna array aperture is (M1)d.sup.[3]. The maximum virtual array aperture of the designed SNLA structure can be increased to ((M+2(M.sub.1+(M.sub.1+1)M.sub.2)+1)/2)d. Compared with the ULA scheme, the angle estimation accuracy of the SNLA is further improved.

(30) The path gain estimation is: once the path angle is estimated by the MR-based algorithm, the array manifold matrix can be obtained. The pilot signal u.sub.t is sent at this stage, so that the path gains in different blocks can be estimated based on the obtained . To obtain the channel gain estimation g.sub.t, it is necessary to minimize the following cost function:

(31) $\begin{array}{cc}J\left({\hat{g}}_t\right)={.Math.\frac{1}{u_t}{y}_t-\hat{A}{\hat{g}}_t.Math.}^2={\left(\frac{1}{u_t}{y}_t-\hat{A}{\hat{g}}_t\right)}^H\left(\frac{1}{u_t}{y}_t-\hat{A}{\hat{g}}_t\right)={\left(\frac{1}{u_t}\right)}^2{y}_t^H{y}_t-\frac{1}{u_t}{y}_t^H\hat{A}{\hat{g}}_t-\frac{1}{u_t}{\hat{g}}_t^H{\hat{A}}^H{y}_t+{\hat{g}}_t^H{\hat{A}}^H\hat{A}{\hat{g}}_t& \left(15\right)\end{array}$

(32) By calculating a partial derivative of the cost function relative to .sub.t, the following can be obtained:

(33) 0$\begin{array}{cc}\frac{J\left({\stackrel{}{g}}_t\right)}{{\stackrel{}{g}}_t}=-\frac{1}{u_t}{y}_t^H\hat{A}+{\hat{g}}_t^H{\hat{A}}^H\hat{A}& \left(16\right)\end{array}$

(34) In a case that

(35) ${\hat{g}}_t^H{\hat{A}}^H\hat{A}=\frac{1}{u_t}{y}_t^H\hat{A},$
a solution of the channel gain is:

(36) $\begin{array}{cc}{\stackrel{}{g}}_t=\frac{1}{u_t}{\left({\hat{A}}^H\hat{A}\right)}^{-1}{\hat{A}}^H{y}_t,t=1,.Math.,T& \left(17\right)\end{array}$

(37) Finally, within a time block, the whole channel estimation result expression is:
.sub.t=.sub.t,t=1, . . . ,T(18)

(38) Based on the estimated channel matrix, the communication symbol can be equalized by a ZF algorithm.

Simulation Example

(39) A base station deploys an SNLA, in each Monte Carlo simulation, the channel path DOA is randomly distributed at (90, 90). The channel estimation property is described by a normalized mean square error (NMSE), that is, NMSE=Eh.sub.2.sup.2/Eh.sub.2.sup.2, where =[.sub.1, . . . ,.sub.T]. The signal to noise ratio (SNR) is defined as

(40) $10{\log }^{\frac{P_s}{{}^2}},$
where P.sub.s is the normalized signal power which is fixed at 1. In a single-user communication system, parameters are set as follows: M=15, K=1, L=4, T=16. To verify the angular resolution of the SNLA and ULA, this example simulates a set of normalized frequency spectrum with a dense incident wave direction, where the path angle is specified as (10, 4, 4, 10). A simulation result is shown in FIG. 4. Note: dots represent the true angle. It is observed that compared with the ULA scheme, the frequency spectrum of the SNLA is sharper in the incident wave direction. Therefore, the angle estimation accuracy based on the SNLA is higher than the ULA.

(41) To verify the accuracy of the provided channel estimation algorithm, the relationship curve of NMSE and SNR under different conditions is drawn in FIG. 5. In simulation, the base station antennas are respectively arranged as the SNLA, the ILIA and a random array (RAND). The path angle is set as: (16, 4, 4,16). Therefore, the NMSE of the method (that is, MR-SNLA) of the present invention is much smaller than those of the existing simultaneous orthogonal matching pursuit (SOMP), low rank structured covariance reconstruction (LRSCR) and multiple signal classification (MUSIC) algorithms. In addition, by changing the antenna distribution at the base station, it is observed that the SNLA is more excellent than the ULA and the RAND. This is because the method of the present invention not only utilizes the sparsity of the channel, but also utilizes the geometric structure of the array, thereby improving the channel estimation accuracy.

(42) To further describe the advantages of the provided method, FIG. 6 compares the BER performance of different channel estimation algorithms modulated by quadrature phase shift keying (QPSK). For different channel estimation methods, the signal is modulated by the same equalizer (that is, 7T). To maintain the fairness of comparison, the system model and parameters required by simulation are the same as those in FIG. 5. In FIG. 6, the bit error ratio of the channel estimation method provided by the present invention under the SNLA is closer to the perfect CSI condition, and has a certain advantage compared with the traditional channel estimation method.