Power distribution method based on weighting capacity speedup in millimeter wave multiaerial system
Technical field
The invention belongs to wireless communication technology field, is related to a kind of be used in millimeter wave multiaerial system based on weighting capacity
The power distribution method of speedup, it is a kind of method of the transmission power for multiple terminal distribution base stations.
Background technology
Current cellular radio Communication frequency band is very crowded, and is difficult to the need for meeting growing data service
Ask.Radio wave frequency is improved to the frequency of millimeter wave, can effectively solve the crowded present situation of current frequency band.Also, millimeter wave
Wavelength it is short, the length of antenna can shorten, and base station and terminal can configure a large amount of antennas, improve the utilization of transmission band
Rate.In millimeter wave multiaerial system, base station is configured with hundreds of or thousands of antennas, and each terminal configures more antennas.Due to
Based on direct path, there is correlation in the channel of different terminals for the propagation of signal, cause terminal room interference to be present.Power distribution
It is a kind of effective ways for mitigating terminal room interference.Because the annoyance level of different terminal rooms is different, reasonable distribution base station work(
Rate can reduce interference as far as possible, put forward high-power utilization ratio.
Power distribution in millimeter wave multiaerial system is faced with problem.Power distribution method and base based on geometric programming
Need former power distribution problems approximation, approximate error can cause power distribution knot in gradually convex approximate power distribution method
Fruit is not optimal.Power distribution method maximum power efficiency rather than spectrum efficiency based on fractional programming, are caused interference with
It can not be preferably minimized.
The content of the invention
The purpose of the present invention is the problem of can not distributing base station power in an optimal manner for existing method, that is, to pass through distribution
Base station power is at utmost to reduce the problem of interference, there is provided based on weighting capacity speedup in a kind of millimeter wave multiaerial system
Power distribution method.
For the inventive method based on a communication system comprising a base station and multiple terminals, base-station antenna array is uniform
Rectangle battle array, and include hundreds of or thousands of antennas;Base station reduces signal dimension, the transmission power of base station using beam space method
It is limited.
The inventive method comprises the following steps:
Step 1. initial phase:
The singular value decomposition of beam space channel correlation matrix is calculated, obtains the matrix and singular value group of singular vector composition
Into matrix, calculate the weight of each terminal;General power is distributed in the way of evenly distributing.
The step 2. iteration distribution power stage:
2-1. belongs to the terminal of the set of the numbering for the terminal for needing distribution power for all numberings, calculates weighting
Capacity to the derivative of power, then calculate these derivatives variance and:If the value of the variance sum is more than the threshold value of setting, enter
2-2;If the value of the variance sum is less than or equal to the threshold value of setting, end step 2, terminate algorithm;
2-2. belongs to the terminal of the set for the numbering for gathering the terminal for needing distribution power for all numberings, finds each
The capacity of the weighting of terminal is to the maximum and minimum value of the derivative of power, and the numbering of corresponding terminal;Calculate the two
Second dervative of the capacity of the weighting of terminal to power;According to second dervative, under conditions of ensureing that general power is constant, tool is calculated
There is the power reduction amount of the power incrementss of the terminal of the maximum derivative of weighting and the terminal of the minimum derivative with weighting, make two
The capacity of the weighting of individual terminal is equal after the change to the derivative of power;Calculating has weighting after reducing according to power reduction amount
The power of the terminal of minimum derivative:
If the power of the terminal of the minimum derivative with weighting is more than zero after subtracting decrement, increase according to power incrementss
Add the power of the terminal of the maximum derivative with weighting, reduce the terminal of the minimum derivative with weighting according to power reduction amount
Power;
If the power after decrement that subtracts of the terminal of the minimum derivative with weighting is less than or equal to zero, again by decrement
The power of terminal corresponding to minimum value is arranged to, correspondingly resets incrementss under conditions of general power is constant, by maximum
The power of terminal corresponding to value increases according to incrementss, and the power of terminal corresponding to minimum value is reduced according to decrement, and will
It numbers the set for the numbering that the terminal for not needing distribution power is moved on to from the set of the numbering for the terminal for needing distribution power.
Step 3. return to step 2 judges whether iteration continues.
In millimeter wave multiaerial system, antenna for base station quantity is extremely more.Therefore, the channel vector of the beam space of each terminal
Between there is near orthogonality.So, power system capacity has approximate convex characteristic and monotonicity.Therefore, step 2 reaches each terminal
The derivative approximately equal of the capacity of weighting just can almost reach the maximum of power system capacity, can ignore the mistake brought by approximation
Difference.
The inventive method make use of the approximate convex characteristic and monotonicity of the capacity in millimeter wave multiaerial system to distribute work(
Rate.The inventive method utilizes the capacity of the weighting of each terminal of this property calculation to the derivative of power, that is, the capacity speedup weighted,
Terminal distribution power is given according to the capacity speedup of weighting.In the case of the capacity speedup identical of the weighting of individual terminal, Ke Yishi
Now maximized progressive power system capacity, i.e. antenna for base station quantity tend to be infinite under power system capacity.Because the inventive method utilizes
Be capacity characteristic in millimeter wave multiaerial system, i.e., antenna for base station it is a fairly large number of in the case of capacity characteristic, its approximation is by mistake
Difference is smaller, can obtain more excellent power distribution result.
The innovation key of the inventive method is:Using the approximate convex characteristic and monotonicity of the capacity of millimeter-wave systems, with most
The difference of the derivative of the capacity of the weighting of each terminal of smallization is Target Assignment power, can reach approximately the maximum of power system capacity
Value.The interference of terminal room can be mitigated using the inventive method, improve power system capacity.
Brief description of the drawings
Fig. 1 is the flow chart of the inventive method.
Fig. 2 is the analogous diagram of system spectral efficiency and Between Signal To Noise Ratio in the embodiment of the present invention.
Fig. 3 is system spectral efficiency and the analogous diagram of antenna for base station quantitative relation in the embodiment of the present invention.
Embodiment
To make the object, technical solutions and advantages of the present invention clearer, the present invention is made below in conjunction with the accompanying drawings further
Detailed description.
The application scenarios of the inventive method are introduced first:Base station placed a uniform rectangular battle array for having root antenna, each
There is individual terminal in cell, each terminal has root antenna;Base station constitutes matrix H to the beam space channel of all terminalsb;Hb's
Conjugate transposition isThe power that the terminal is distributed in base station is pk;The transmission power of base station is ρ;The power of distribution is not zero
Terminal numbering composition collection be combined intoThe collection of the numbering composition of remaining terminal is combined into;The weights of the terminal are bk;Tool
The capacity for having the weighting of maximum is k to the numbering of the terminal of the derivative of powerx, corresponding power incrementss are δ (kx), have most
The capacity of small weighting is k to the numbering of the terminal of the derivative of powery, corresponding power reduction amount is δ (ky);fkFor end
The capacity at end;For the terminal capacity to the derivative of power,For weighting the terminal capacity to power
Derivative;σ is the threshold value of iteration.
Referring to Fig. 1, the following operative step of the inventive method is introduced:
(1) initial phase:In the way of evenly distributing, calculated using the weight and base station transmitting power of each terminal
The power distribution of initialization.The operation content of the step is:Empty set will be arranged to, willIt is arranged to include the numbering of all terminals
Set.Calculate the weights b of all terminalskSum, calculate the weights b of transmission power ρ and all terminalskSum ratio, by this
Power of the ratio as each terminal.
(2) the iteration distribution power stage:Calculate sum of the capacity of weighting to the variance of the derivative of power, according to the numerical value and
The magnitude relationship of the threshold value of setting decides whether to enter iteration.If the numerical value is more than the threshold value of setting, continue iteration;Otherwise stop
Only iteration.The capacity of maximum weighted is found to the derivative of power and the capacity of minimum weight to terminal corresponding to the derivative of power.
Calculate the power knots modification of the two terminals and change.Rejudge whether iteration.The operation content of the step is:For all
Numbering belong to setTerminal, calculate derivative of the capacity to power of weightingCalculate the variance of these derivatives
With if this and value are more than threshold value σ, continue follow-up iterative step, otherwise stop iteration.Belong to collection for all numberings
CloseTerminal, find derivative of the capacity to power of weightingMaximum and minimum value corresponding to terminal numbering, most
The numbering of terminal is k corresponding to big valuex, the numbering of terminal is k corresponding to minimum valuey.Calculate the capacity of the weighting of the two terminals
To the derivative and second dervative of power.According to these derivatives, under conditions of ensureing that general power is constant, kth is calculatedxIndividual terminal
Power incrementss δ (kx) and kthyThe power decrement δ (k of individual terminaly), make the capacity of weighting of two terminals to power
Derivative is equal after the change.If kthyThe power of individual terminal subtracts decrement δ (ky) it is more than zero afterwards, by kthxThe work(of individual terminal
Rate increases according to incrementss, by kthyThe power of individual terminal reduces according to decrement;If kthyThe power of individual terminal subtracts reduction
Measure δ (ky) it is less than or equal to zero afterwards, again by δ (ky) it is arranged to kthyThe power of individual terminal, correspondingly in the constant condition of general power
Under reset δ (kx), by kthxThe power of individual terminal increases according to incrementss, by kthyThe power of individual terminal is according to decrement
Reduce, and by kyFrom setMove on to set.
(3) starting for coming back to step (2) judges whether iteration continues.
In order to show the practical performance of the inventive method, implement to test applicant carried out Multi simulation running.L-G simulation test
As a result as shown in Figures 2 and 3, emulated respectively in terms of signal to noise ratio and array antenna quantity two.For intuitively body
Reveal the superiority of the inventive method, by the simulation result of this method and existing channel estimation methods, perfect channel estimation side
Method is contrasted.
As seen from Figure 2, the inventive method significantly enhances the frequency spectrum effect of system in the case where signal to noise ratio is relatively low
Rate, that is to say, that the inventive method can mitigate interference under Low SNR.Fig. 3 shows that the inventive method obtains frequency
Spectrum efficiency improves as array antenna quantity increases, and gradual compared to the gap of the spectrum efficiency of conventional power distribution method
Increase.
The preferred embodiments of the present invention are the foregoing is only, are not intended to limit the invention, all spirit in the present invention
Within principle, any modification, equivalent substitution and improvements done etc., it should be included within the scope of protection of the invention.