Method for channel estimation by using sequential pilot frequency sequence
Technical Field
The invention belongs to the technical field of large-scale multi-antenna wireless communication, and particularly relates to a method for performing channel estimation by using a sequential pilot frequency sequence.
Background
In the transmission scheme, a large number of low-power antennas are arranged at a base station end, and the number of the antennas far exceeds the number of single-antenna users scheduled at the same time. Massive MIMO technology can enable wireless communication systems to achieve very high throughput. In the conventional channel estimation method, users transmit pilot sequences simultaneously, and a base station performs channel estimation on all users at the same time. The conventional method requires the design of complex orthogonal pilot sequences, and the design of the pilot sequences depends on the number of users.
Disclosure of Invention
In view of the above-mentioned deficiencies or needs in the art, the present invention provides a method for channel estimation and data communication with multiple users transmitting pilots sequentially. The channel of the user can be estimated without complex pilot orthogonal sequences and the system performance is hardly lost.
In order to achieve the above object, the present invention provides a method for channel estimation using sequential pilot sequences, comprising:
step 1, K users in a cell send pilot frequencies to a base station in a sequential manner: k users in a cell simultaneously send signals to a base station on the same channel, and the base station simultaneously provides data transmission service for the K users by adopting a large-scale MIMO technology; the K users sequentially send pilot frequency to the base station, only one user sends the pilot frequency to the base station at one time point, and other K-1 users simultaneously send uplink data to the base station;
step 2, the base station carries out channel estimation in a sequential pilot frequency sending mode: a base station provides service for K users on one channel simultaneously; at a certain time point, only user K sends pilot frequency to the base station, the rest K-1 users all send uplink data to the base station, and the channels of the rest K-1 users are obtained in the previous pilot frequency sending and channel estimation stages; after receiving the signals sent by the K users, the base station estimates the channel of the K users.
In an embodiment of the present invention, in step 2, after receiving the signals sent by the K users, the base station estimates a channel of the user K, specifically, the base station estimates the channel by using a projection method:
computingThen Y' is used to channel g of user kkCarrying out direct estimation;
wherein: y' is Y atIs a space perpendicular to the space spanned by the column vectors ofThe component of (a) to (b),is Y atThe column vector of (a),is composed ofA projection matrix of the space spanned by the column vectors of (a),a channel matrix consisting of channels of other K-1 users which have been estimated by the base station, giIs an N x 1 dimensional channel vector from user i to the base station,k, Y is a vector obtained after a base station receives a signal matrix Y sent by K users and performs a vectorization operator, and N is the number of antennas of the base station.
In an embodiment of the present invention, in step 2, after receiving the signals sent by the K users, the base station estimates a channel of the user K, specifically, performs channel estimation by using a linear minimum mean square error method, where the estimated user channel is:
wherein: skPilot sequences of length L, p, sent for a userdPower at which data is transmitted for a user, σ2Is the power of white gaussian noise and,a channel matrix formed by the channels of other K-1 users which are estimated by the base station,is channel gkSecond moment of (I)N×NIs an N × N dimensional identity matrix, IL×LThe unit matrix is L multiplied by L dimension, Y is a vector obtained after the base station receives a signal matrix Y sent by K users and carries out a vectorization operator, and N is the number of antennas of the base station.
In an embodiment of the present invention, the sequential pilot sending mode is a packet sequential pilot sending mode, and specifically includes: grouping users, different groups of users in different time slotsSending pilot frequency, wherein K is the number of groups of users, the length of a time slot used for sending the pilot frequency in each group is L symbol periods, and the correlation time of a channel is T; the users in the same group use simple orthogonal pilot to transmit pilot in the same time slot, i.e. the pilot sequences transmitted by users i and k in the same group And satisfy the orthogonal relation <math>
<mrow>
<munderover>
<mi>Σ</mi>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<mrow>
<mi>L</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</munderover>
<msubsup>
<mi>s</mi>
<mi>j</mi>
<mi>i</mi>
</msubsup>
<msup>
<mrow>
<mo>(</mo>
<msubsup>
<mi>s</mi>
<mi>j</mi>
<mi>k</mi>
</msubsup>
<mo>)</mo>
</mrow>
<mo>*</mo>
</msup>
<mo>=</mo>
<mn>0</mn>
<mo>,</mo>
</mrow>
</math> Where L is the length of the pilot sequence.
In an embodiment of the present invention, the sequential pilot sending method is a user grouping sequential pilot sending method in multiple cells, and specifically includes: grouping users of adjacent cells, wherein each group comprises a plurality of users from different cells, the users of different groups send pilot frequencies in different time slots, the users of the same group send pilot frequencies in the same time slot, and the pilot frequency sequences of the users of the same group are mutually orthogonal.
In general, compared with the prior art, the technical solution of the present invention can effectively solve the problem that the pilot sequence in the conventional channel estimation method needs to be designed in a complex orthogonal manner, and can significantly reduce the length of the required pilot sequence and improve the throughput of the system. The method is characterized in that:
1. solving complex orthogonal design problems in pilot design
In the prior art, users simultaneously transmit pilots to a base station, and in order to perform good estimation on the channels of the users, the pilot sequences transmitted by the users need to be kept orthogonal, so that complex orthogonal design is needed when designing the pilots. The technical scheme of the invention does not need to carry out orthogonal design on the pilot frequency of each user.
2. Reducing pilot length per user
In the prior art, the length of the pilot frequency required for ensuring the orthogonality of the pilot frequency sequence of each user is not less than the number of the users, so the pilot frequency length is increased along with the increase of the number of the users. The pilot frequency length in the technical scheme of the invention can be very short and cannot be increased along with the increase of the number of users.
3. Improving the throughput of a system
Although the pilot signal transmitted by the user and the data signal transmitted by other users do not have orthogonality, the base station can still obtain a good estimated channel of the user in a massive MIMO system. The technical scheme of the invention reduces the length of the pilot frequency, increases the time frequency resources for sending the data signals, thereby improving the utilization rate of the time frequency resources, and can support more users and improve the throughput of the system to a certain extent.
Drawings
Fig. 1 is a schematic diagram of a sequential pilot sequence transmission scheme in the present invention.
Detailed Description
In order to make the objects, technical solutions and advantages of the present invention more apparent, the present invention is described in further detail below with reference to the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are merely illustrative of the invention and are not intended to limit the invention. In addition, the technical features involved in the embodiments of the present invention described below may be combined with each other as long as they do not conflict with each other.
As shown in fig. 1, in the massive MIMO communication system proposed by the present invention, the method for performing channel estimation and data communication by sequentially transmitting pilots by multiple users is as follows:
multiple users in a cell simultaneously transmit signals to a base station on the same channel (frequency), and the base station simultaneously provides data transmission services for the multiple users by adopting a large-scale MIMO technology. A plurality of users send pilot frequency to the base station sequentially, only one user sends the pilot frequency to the base station at one time point, and other users can send uplink data to the base station at the same time.
Suppose a base station has N antennas and serves K users simultaneously on one channel. At a certain time point, only user K sends pilot frequency to the base station, the rest K-1 users all send uplink data to the base station, and the channels of the rest K-1 users are obtained in the previous pilot frequency sending and channel estimation stages. After receiving the signals sent by the K users, the base station may use a projection method, an LMMSE (linear minimum mean square error) method, to suppress interference caused by data signals from K-1 users when estimating channels of user K.
Let the channel matrix be G ═ G1 g2...gK]Wherein g isiIs an N x 1 dimensional channel vector from user i to the base station. Assume that the pilot length is L. The signal matrix sent by the user is X ═ X1 x2...xK]TWherein x isi=[xi,1 xi,2...xi,L]TSignal sequence, x, transmitted in L symbol periods for the ith userk=sk=[si,1 si,2...si,L]TThe signal matrix received by the base station is Y ═ GX + W, and W is zero-mean gaussian white noise in the channel. Before estimating the channel of user K, the channels of other K-1 users are estimated, and the channel matrix formed by the channels isThe matrix formed by the real channels of other K-1 users isThe K-1 users transmit data matrix ofThe implementation of the projection method and LMMSE method is described in detail below.
Projection method: when a certain user is estimated, a method of projecting a signal received by a base station to a vertical space of a known user channel space for estimation is adopted.The projection matrix of the space spanned by the column vectors ofThus Y is inProjected onto the spanned space isThus y is inIs formed as a space perpendicular to the spaceComponent ofThe channel of user k is then directly estimated using Y'.
LMMSE method: the known observation vector y ═ Vec (Y) is used to determine the unknown parameter g by using LMMSEkThe formula for the estimation isWhere Vec (·) is the vectorization operator. Suppose a known channel gkSecond moment of From E { gkSubstituting 0 for E { y } into 0
By calculating:
let w ═ vec (w),
wherein,is a Kronecker product operation, pdPower at which data is transmitted for a user, σ2Is the power of gaussian white noise.
Thereby to obtain
<math>
<mrow>
<msub>
<mover>
<mi>g</mi>
<mo>^</mo>
</mover>
<mi>k</mi>
</msub>
<mo>=</mo>
<msub>
<mi>R</mi>
<mrow>
<msub>
<mi>g</mi>
<mi>k</mi>
</msub>
<msub>
<mi>g</mi>
<mi>k</mi>
</msub>
</mrow>
</msub>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>s</mi>
<mi>k</mi>
</msub>
<mo>⊗</mo>
<msub>
<mi>I</mi>
<mrow>
<mi>N</mi>
<mo>×</mo>
<mi>N</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mi>H</mi>
</msup>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>s</mi>
<mi>k</mi>
</msub>
<mo>⊗</mo>
<msub>
<mi>I</mi>
<mrow>
<mi>N</mi>
<mo>×</mo>
<mi>N</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<msub>
<mi>R</mi>
<mrow>
<msub>
<mi>g</mi>
<mi>k</mi>
</msub>
<msub>
<mi>g</mi>
<mi>k</mi>
</msub>
</mrow>
</msub>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>s</mi>
<mi>k</mi>
</msub>
<mo>⊗</mo>
<msub>
<mi>I</mi>
<mrow>
<mi>N</mi>
<mo>×</mo>
<mi>N</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mi>H</mi>
</msup>
<mo>+</mo>
<msub>
<mi>p</mi>
<mi>d</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>I</mi>
<mrow>
<mi>L</mi>
<mo>×</mo>
<mi>L</mi>
</mrow>
</msub>
<mo>⊗</mo>
<msub>
<mover>
<mi>G</mi>
<mo>·</mo>
</mover>
<mi>k</mi>
</msub>
<mo>)</mo>
</mrow>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>I</mi>
<mrow>
<mi>L</mi>
<mo>×</mo>
<mi>L</mi>
</mrow>
</msub>
<mo>⊗</mo>
<msub>
<mover>
<mi>G</mi>
<mo>·</mo>
</mover>
<mi>k</mi>
</msub>
<mo>)</mo>
</mrow>
<mi>H</mi>
</msup>
<mo>+</mo>
<msup>
<mi>σ</mi>
<mn>2</mn>
</msup>
<msub>
<mi>I</mi>
<mrow>
<mi>NL</mi>
<mo>×</mo>
<mi>NL</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mi>y</mi>
<mo>.</mo>
</mrow>
</math>
Since in practice the base station cannot know the channel matrixThus using the estimated valueInstead. Thus, it is possible to prevent the occurrence of,
in an embodiment of the present invention, the sequential pilot sending mode is a packet sequential pilot sending mode, and specifically includes: the users are grouped, and the users of different groups send the pilot frequency in different time slots. If the users are divided into K groups, the length of the time slot for transmitting the pilot frequency in each group is L symbol periods, and the correlation time of the channel is T, it is required to satisfy K × L ≤ T. The users in the same group use simple orthogonal pilot to transmit pilot in the same time slot, i.e. the pilot sequences transmitted by users i and k in the same group And satisfy the orthogonal relation <math>
<mrow>
<munderover>
<mi>Σ</mi>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<mrow>
<mi>L</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</munderover>
<msubsup>
<mi>s</mi>
<mi>j</mi>
<mi>i</mi>
</msubsup>
<msup>
<mrow>
<mo>(</mo>
<msubsup>
<mi>s</mi>
<mi>j</mi>
<mi>k</mi>
</msubsup>
<mo>)</mo>
</mrow>
<mo>*</mo>
</msup>
<mo>=</mo>
<mn>0</mn>
<mo>,</mo>
</mrow>
</math> Where L is the length of the pilot sequence.
In another embodiment of the present invention, the sequential pilot sending method is a user grouping sequential pilot sending method in multiple cells, and specifically includes: grouping users of adjacent cells, wherein each group comprises a plurality of users from different cells, the users of different groups send pilot frequencies in different time slots, the users of the same group send pilot frequencies in the same time slot, and the pilot frequency sequences of the users of the same group are mutually orthogonal.
The following describes the method of the present invention by taking the sequential pilot transmission of multiple users in the same cell as an example, which is specifically as follows.
Example (b):
in this specific embodiment, 10 users in a cell simultaneously transmit signals to a base station with 128 antennas on the same channel (frequency), and the base station simultaneously provides data transmission services for these multiple users by using massive MIMO technology, assuming that the pilot length of each user is L ═ 1.
At a certain point in time, only user k (k ═ 1,2, …,10) transmits a pilot to the base station, the remaining 9 users all transmit uplink data to the base station, and the channels of the remaining 9 users have been obtained in the previous pilot transmission and channel estimation phase. After receiving the signals sent by the 10 users, the base station may suppress interference caused by data signals from other 9 users when estimating the channel of user k by using a projection method or an LMMSE (linear minimum mean square error) method, and estimate the channel of user k. Suppose that the signal vector received by the base station is Y, and the channel matrix formed by the known channel vectors of the users is Y
Projection method: when estimating the channel of user k, it uses the method that projects the signal received by base station to the vertical space of the known user channel space to estimate. Y is atIs formed as a space perpendicular to the spaceComponent ofAnd then Y' is used for carrying out channel estimation on the user k.
LMMSE method:
wherein sigma2Is the power of the channel Gaussian white noise, skPilot sequences of length L, p, sent for a userdThe power at which the data is transmitted for the user,is channel gkSecond order moment of (a).
It will be understood by those skilled in the art that the foregoing is only a preferred embodiment of the present invention, and is not intended to limit the invention, and that any modification, equivalent replacement, or improvement made within the spirit and principle of the present invention should be included in the scope of the present invention.