Summary of the invention
The object of this invention is to provide a kind of reasonable distribution and utilize limited pilot resources, support the technical method of dense distribution formula MIMO mobile communication.
The invention provides the pilot multiplex method in a kind of Dense nodes configuration-system, it is characterized in that:
(1) according to the user number K of access and the statistical channel state information of user, the interference matrix G between structuring user's, wherein G
i,jrepresent that i-th user (UE) is to the interference between a jth UE, the diagonal entry of G is set to 0.
(2) system designs one group of orthogonal guide frequency in advance, and available orthogonal guide frequency number is N, and pilot tone sequence number is expressed as 1,2 ..., N;
(3) set pilot frequency distribution vector as ξ=[ξ
1... ξ
i... ξ
k], ξ
i=n (n=1...N) represents that i-th UE uses the n-th orthogonal guide frequency.If the UE number using the n-th orthogonal guide frequency is K
n, then have
(4) carry out pilot frequency distribution according to interference matrix to user, pilot frequency distribution is obtained by following iterative process:
Steps A, initialization pilot frequency distribution vector, ξ
i=1 (i=1 ..., K), all users use No. 1 pilot tone; If I
maxfor use same pilot UE between maximum interference, be initially
n
tempfor the interim orthogonal guide frequency sequence number in iteration, be initially n
temp=1.If ξ '=(ξ
1' ... ξ
i' ... ξ '
k) for recording the pilot frequency distribution vector of current iteration process, be initially ξ '=ξ.
Step B, carries out iteration,
Successively pilot frequency distribution is carried out to user.When to user i (i=1,2 ..., K) when distributing pilot tone, traversal pilot tone 1 to pilot tone N, if
so
current pilot n is distributed to user i simultaneously, continue to perform ergodic process.
Step C, judges pilot frequency distribution result.If ξ ' is identical with ξ, then terminate pilot frequency distribution process, obtain pilot frequency distribution result; Otherwise, make ξ '=ξ, continue step B.
For the pilot multiplex method in above-mentioned a kind of Dense nodes configuration-system, the building method of a kind of interference matrix G is:
Wherein d
i,jrepresent the distance between i-th UE and a jth UE.
The present invention adopts technique scheme, has following beneficial effect:
1, when ensureing that CSI is accurately estimated in base station, can the good orthogonal pilot frequency sequence of designed in advance, and pilot number can be fixing, which greatly enhances uniformity and the operability of real system, the pilot-frequency expense of system can be reduced simultaneously, improve the information rate of system;
2, when orthogonal guide frequency number is given, can distribute simply and effectively and utilize existing pilot resources, reducing the pilot pollution between user, improving the accuracy that CSI is estimated in base station;
3, system is adaptive, change at intra-cell users number, intra-cell users change in location, during the change of intra-cell users clustered pattern, the pilot multiplex method that this invents provides still can be stablized and reasonably distribute pilot tone, minimized by pilot pollution between user, guarantee system normally works;
4, this to invent the pilot multiplex method complexity that provides low, but reliability and feasibility is good, system is always disturbed little, and total capacity is much larger than random pilot allocative decision;
5, this invents the pilot multiplex method that provides and has taken into account fairness, ensures that the Pilot Interference residing for each user minimizes, number of users and pilot number ratio comparatively rationally time, there is not the possibility that can not communicate under individual user is in strong jamming.
Embodiment
The present invention program is understood better in order to make those skilled in the art person, below in conjunction with the accompanying drawing in the embodiment of the present invention, clear, complete description is carried out to the technical scheme in the embodiment of the present invention, obviously, described embodiment is only the present invention's part embodiment, instead of whole embodiments.Based on the embodiment in the present invention, those of ordinary skill in the art, not making the every other embodiment obtained under creative work prerequisite, should belong to the scope of protection of the invention.
Suppose total K user in community, available orthogonal guide frequency number is N.
By pilot tone by user grouping, each group the multiplexing same pilot tone of user, be divided into N (1,2 ..., N) and group.Such as, the user of first group by pilot tone 1, i-th group use pilot tone i, N group pilot tone N.User's number symbol n of each group
1, n
2... n
n(often organizing number of users not necessarily identical) represents.
The pilot signal received is:
[formula 1]
Represent and sue for peace to the pilot signal of the i-th group of user the received user of useful pilot tone i (namely), user's number is n
i.G
i,jbe the matrix of MN × 1, represent the channel matrix of i-th group of jth user to M RAUs
Wherein, Λ
i,jrepresent that corresponding large scale declines, h
i,jrepresent multipath fading, n
p,irepresent additive Gaussian noise, variance is γ
p.The MMSE channel estimating of [formula 1] can be expressed as
Definition channel estimation errors is,
Be defined as follows matrix:
G=[G
1… G
N],
Represent ideal communication channel matrix.Correspondingly, the definition matrix that obtains of channel estimating and channel estimation errors matrix as follows,
Then ideal communication channel matrix can be expressed as
So at BPU, utilize the channel matrix that channel estimating obtains, received signal vector and the relation between transmitting can be expressed as
Wherein y is the received signal vector of MN × 1,
S
lbe the transmission signal vector of l group user W × 1, n is covariance matrix is ε (nn
h)=γ
uLi
mNthe multiple additive white Gaussian noise in MN × 1.Further, the covariance matrix of channel estimating interference plus noise can be calculated by following formula,
The lower bound of the uplink channel capacity of known system can be determined by following formula
Due to det (I+A
-1b)=det (A+B)-det (A), can obtain
Wherein,
use law of great number, can obtain,
Wherein
Definition matrix
so can obtain
Definition matrix
for block diagonal matrix
Ξ=diag([Ξ
1…Ξ
n…Ξ
N])
Known
So the limiting value that can obtain lower bound capacity is
The specific embodiments of pilot frequency multiplexing, as shown in Figure 3:
Step 1: system designs one group of orthogonal guide frequency in advance, available orthogonal guide frequency number is N, and pilot tone sequence number is expressed as 1,2 ..., N;
If pilot frequency distribution vector is ξ=[ξ
1... ξ
i... ξ
k], ξ
i=n (n=1...N) represents that i-th UE uses the n-th orthogonal guide frequency.If the UE number using the n-th orthogonal guide frequency is K
n, then have
Step 2: according to the user number K of access and the location distribution of user, the interference matrix G between structuring user's, wherein G
i,jrepresent that i-th user (UE) is to the interference between a jth UE, the diagonal entry of G is set to 0,
Wherein d
i,jrepresent the distance between i-th UE and a jth UE.
Step 3: pilot frequency distribution is carried out to user according to interference matrix, pilot frequency distribution is obtained by following iterative process:
(1) initialization pilot frequency distribution vector, ξ
i=1 (i=1 ..., K), all users use No. 1 pilot tone; If I
maxfor use same pilot UE between maximum interference, be initially
n
tempfor the interim orthogonal guide frequency sequence number in iteration, be initially n
temp=1.If ξ '=(ξ
1' ... ξ
i' ... ξ '
k) for recording the pilot frequency distribution vector of current iteration process, be initially ξ '=ξ.
(2) iteration is carried out
Successively pilot frequency distribution is carried out to user.When to user i (i=1,2 ..., K) and when distributing pilot tone, carry out the traversal of pilot tone 1 to pilot tone N, if
so
current pilot n is distributed to user i simultaneously, continue traversal.
(3) pilot frequency distribution result is judged.If ξ ' is identical with ξ, then terminate pilot frequency distribution process, obtain pilot frequency distribution result; Otherwise, make ξ '=ξ, continue the iterative process of (2).
At the end of iterative process, just obtain pilot frequency distribution result.At this moment, carry out performance evaluation, comparison system capacity, reference object is random pilot allocative decision.The capacity comparison of change user uplink power, as shown in Figure 1.As seen from Figure 1, the solution of the present invention is better than random pilot scheme when ascending power is greater than 5dBm, along with the increase of power, finally have an appointment 32% gain.The capacity comparison of change orthogonal guide frequency number, as shown in Figure 2.As seen from Figure 2, when the solution of the present invention changes pilot number at fixing ascending power, be better than random pilot scheme all the time, have an appointment at most 32% gain.
In the embodiment that the application provides, should be understood that, disclosed method, not exceeding in the spirit and scope of the application, can realize in other way.Current embodiment is a kind of exemplary example, should as restriction, and given particular content should in no way limit the object of the application.Such as, multiple unit or assembly can in conjunction with or another system can be integrated into, or some features can be ignored, or do not perform.