CN106851833A - MIMO power distribution methods and system based on high specific transfer pre-coding - Google Patents

Abstract

A kind of MIMO power distribution methods and system based on high specific transfer pre-coding are the embodiment of the invention provides, the method includes：The amplitude and phase of the radio circuit gain of the terminal and base station of the multi-input multi-output system of high specific transfer pre-coding are obtained, up channel model and down channel model is determined.According to up channel model and down channel model, it is determined that approximate and speed.According to approximate and speed, it is the nonconvex property majorized function of target to set up to maximize approximate and speed.By the way that function to be converted to the formula of convexity by nonconvex property, nonconvex property majorized function is converted into convexity majorized function.When approximate and speed is maximum, the optimal solution of convexity majorized function is determined, the descending power of multi-input multi-output system is configured according to optimal solution.By the MIMO power distribution methods based on high specific transfer pre-coding of the invention, realize improving the reasonability of telecommunication system resources configuration, save the communication resource of system.

Description

MIMO power distribution methods and system based on high specific transfer pre-coding

Technical field

The present invention relates to wireless communication technology field, more particularly to the MIMO power based on high specific transfer pre-coding point Method of completing the square and system.

Background technology

Due to MIMO (Multiple-Input Multiple-Output, multiple-input and multiple-output) systems and a single aerial system Compare, with the superior function such as diversity and spatial multiplexing gain, speed higher and reliability, and be considered as cellular communication system one Plant important technology.In recent years, extensive mimo system, i.e. base station install substantial amounts of antenna, by antenna increase to hundreds of with To service tens users simultaneously, to obtain speed and reliability higher.Extensive mimo system is considered as 5G communication systems Unified kind of important key technology.

Base station carries out precoding using descending CSI (Channel State Information, channel condition information) can be with Preferably service more users.For the system of FDD (Frequency Division Duplexing, FDD), it is Downlink channel condition information is estimated, base station should be at least equal to the antenna number of base station to the descending pilot frequency number that user launches.Cause This, when the scale of system become it is very big, i.e. when the antenna number of base station increases to hundreds of, descending pilot frequency expense and descending CSI's Feed back unacceptable by what is become.And for TDD (Time Division Duplexing, time division duplex) system, due to up and Downlink in same frequency range, when sending between interval less than channel coherency time when, up-downgoing channel model will experience identical Physical Attenuation, that is, meet the reciprocity of channel, that is, down channel is the transposition of up channel.Therefore, down channel Estimation can be obtained by uplink channel estimation.As long as and the ascending pilot frequency quantity of user's transmitting ensures more than or equal to number of users just The estimation of up channel can be completed, it is clear that the order of magnitude of user is more much smaller than antenna for base station number.Pilot tone can so be greatly reduced The feedback of expense and CSI.

In the prior art, base station obtains each upward signal in up channel respectively, according to the reciprocity of channel, it is determined that The power of each downstream signal, and the downstream signal that will be determined power configuration on corresponding antenna.But, actual communication System not only includes radio propagation channel, RF (Radio Frequency, the radio frequency) circuit also including both link ends transceiver Part.Generally, RF circuits include blender, A/D and D/A converter, power amplifier etc., and RF circuits by the temperature of external environment condition The influence such as degree and humidity is very big.So, even if the change at random of transceiver RF circuits will cause uplink and downlink signals in same frequency Duan Shang, it is also difficult to keep the reciprocity of this channel.Therefore in the prior art directly according to upward signal model configurating downlink work( Rate, can cause the unreasonable of telecommunication system resources configuration, cause the waste of the system communication resource.

The content of the invention

The purpose of the embodiment of the present invention is to provide a kind of multi-input multi-output system based on high specific transfer pre-coding Downlink power distributing method and system, to realize improving the reasonability that telecommunication system resources are configured, save the communication resource of system. Concrete technical scheme is as follows：

A kind of MIMO power distribution methods based on high specific transfer pre-coding, are applied to MRT high specific transfer pre-codings Extensive MIMO multi-input multi-output systems radio frequency mismatch case, including：

Obtain the width of the radio circuit gain of the terminal and base station of the multi-input multi-output system of high specific transfer pre-coding Degree and phase, determine up channel model and down channel model；

According to the up channel model and the down channel model, it is determined that approximate and speed, wherein, it is described approximate and Speed is obtained by default approximate formula, the corresponding overall transmission rate of the down channel model；

According to the approximate and speed, it is the nonconvex property majorized function of target to set up to maximize the approximate and speed；

By the way that function to be converted to the formula of convexity by nonconvex property, the nonconvex property majorized function is converted into convexity optimization Function；

When the approximate and speed is maximum, the optimal solution of the convexity majorized function is determined, matched somebody with somebody according to the optimal solution Put the descending power of the multi-input multi-output system.

A kind of MIMO power distribution systems based on high specific transfer pre-coding, are applied to MRT high specific transfer pre-codings Extensive MIMO multi-input multi-output systems radio frequency mismatch case, including：

Channel model builds module, terminal and base for obtaining the multi-input multi-output system of high specific transfer pre-coding The amplitude and phase of the radio circuit gain stood, determine up channel model and down channel model；

First computing module, for according to the up channel model and the down channel model, it is determined that approximate and speed Rate, wherein, it is described it is approximate obtained by default approximate formula with speed, the corresponding total transmission of the down channel model Speed；

Second computing module, is target for according to the approximate and speed, setting up to maximize the approximate and speed Nonconvex property majorized function；

3rd computing module, for the formula by the way that function to be converted to convexity by nonconvex property, by nonconvex property optimization Function is converted into convexity majorized function；

Power configuration module, for when the approximate and speed is maximum, determining the optimal solution of the convexity majorized function, The descending power of the multi-input multi-output system is configured according to the optimal solution.

MIMO power distribution methods and system based on high specific transfer pre-coding provided in an embodiment of the present invention, according to end The amplitude and phase of the radio circuit gain of end and base station, determine up channel model and down channel model, according to up letter Road model and down channel model, it is determined that approximate and speed, when approximate and speed is maximum, determines the optimal solution of descending power. According to optimal solution allocation of downlink power, it is possible to achieve improve the reasonability of telecommunication system resources configuration, the communication money of system is saved Source.Certainly, implementing any product of the invention or method must be not necessarily required to while reaching all the above advantage.

Brief description of the drawings

In order to illustrate more clearly about the embodiment of the present invention or technical scheme of the prior art, below will be to embodiment or existing The accompanying drawing to be used needed for having technology description is briefly described, it should be apparent that, drawings in the following description are only this Some embodiments of invention, for those of ordinary skill in the art, on the premise of not paying creative work, can be with Other accompanying drawings are obtained according to these accompanying drawings.

Fig. 1 shows for a kind of flow of the MIMO power distribution methods based on high specific transfer pre-coding of the embodiment of the present invention It is intended to；

Fig. 2 is another flow of the MIMO power distribution methods based on high specific transfer pre-coding of the embodiment of the present invention Schematic diagram；

Fig. 3 is the schematic flow sheet that optimal solution is determined by alternative manner of the embodiment of the present invention；

The curve map that Fig. 4 changes for the system and speed of the embodiment of the present invention with antenna for base station number；

The curve map that Fig. 5 changes for the system and speed of the embodiment of the present invention with base station transmitting power；

Fig. 6 mismatches the curve map of variance change for the system and speed radio frequency circuit of the embodiment of the present invention；

The curve map that Fig. 7 changes for the rate gain of the embodiment of the present invention with base station transmitting power；

Fig. 8 mismatches the curve map of variance change for the rate gain radio frequency circuit of the embodiment of the present invention；

Fig. 9 is the schematic diagram of the MIMO power distribution systems based on high specific transfer pre-coding of the embodiment of the present invention.

Specific embodiment

Below in conjunction with the accompanying drawing in the embodiment of the present invention, the technical scheme in the embodiment of the present invention is carried out clear, complete Site preparation is described, it is clear that described embodiment is only a part of embodiment of the invention, rather than whole embodiments.It is based on Embodiment in the present invention, it is every other that those of ordinary skill in the art are obtained under the premise of creative work is not made Embodiment, belongs to the scope of protection of the invention.

In actual MIMO communication system, not only including radio propagation channel, also including the RF of both link ends transceiver Circuit part.Generally, the transmission circuit of RF circuits is due to the influence of external environment condition, such as temperature and humidity, causes transceiver RF electricity The gain on road produces change.Even if this causes that the TDD system that up channel and down channel are operated in same frequency range can not Enough meet the reciprocity of channel.This mismatch of RF circuits obviously can deteriorate the performance of system, be unfavorable for the precoding of base station With the distribution of power.Even if carrying out reciprocity calibration research in terms of hardware circuit and software algorithm two, it is also difficult to accomplish ideal Eliminate the mismatch problem of RF circuits.

Therefore, it is necessary to the performance when channel has non-matching to system is evaluated, and mismatched based on channel When resource allocation problem, such as base station power assignment problem；It is necessary research under the conditions of channel is unmatched, base station is using each System obtains performance when planting method for precoding, and then studies the influence that each mismatch parameter is brought to system, and according to system How the obtained performance of system carries out resource allocation, such as power distribution, so that it is determined that the best configuration of system, reduces RF and mismatch Influence to system.

By taking the scene of single cell as an example, M (M is positive integer) root antennas are installed while servicing K (K in base station to the embodiment of the present invention It is positive integer) individual user, meet M ＞ ＞ K, and system operatio is in tdd mode.First, on the basis of RF circuit gain models On, channel model when RF is mismatched is set up, use MRT (Max Ratio Transmisson, high specific transmission) for base station Precoding, analyzes the approximate and speed obtained when extensive mimo system channel is mismatched.Then, according to base station transmitting power Constraints, it is the optimization problem of target to set up to maximize approximate and speed.Due to the target letter of the optimization problem of foundation Number is nonconvex property, will be optimized by default logarithm lower bound inequality (lower bound will become very tight near specific value) The object function of problem is converted into and optimizes its lower bound.Because the power that the function after conversion distributes to each user is still non-convex , and then a kind of exponential transform is used, the problem is converted into the difference of a linear function and logarithmic function so that the problem becomes Into convex problem.Finally, by the parameter in iteration more new lower bound inequality, with Step wise approximation optimal solution.

Referring to Fig. 1, Fig. 1 is the one of the MIMO power distribution methods based on high specific transfer pre-coding of the embodiment of the present invention Schematic flow sheet is planted, extensive multi-input multi-output system radio frequency is applied to and is mismatched environment, including：

S101, obtains the radio circuit gain of the terminal and base station of the multi-input multi-output system of high specific transfer pre-coding Amplitude and phase, determine up channel model and down channel model.

Investigation end side (user side) and the amplitude and phase model of base station side RF circuit gains, then set up RF circuits not Up channel model and down channel model during matching.

S102, according to up channel model and down channel model, it is determined that approximate and speed, wherein, approximate and speed is Obtained by default approximate formula, the corresponding overall transmission rate of down channel model.

The corresponding overall transmission rate of down channel model is calculated, as system and speed, by approximate formula by system Approximate and speed is converted to speed.

S103, according to approximate and speed, it is the nonconvex property majorized function of target to set up to maximize approximate and speed.

It is determined that after approximate and speed, in order that mimo system transmission rate is maximized, setting up to maximize approximate and speed Rate is the nonconvex property majorized function of target.

S104, by the way that function to be converted to the formula of convexity by nonconvex property, is converted into convexity excellent by nonconvex property majorized function Change function.

Nonconvex property function cannot solve extreme value, need for nonconvex property majorized function to be converted to convexity majorized function.

S105, when approximate and speed is maximum, determines the optimal solution of convexity majorized function, and multi input is configured according to optimal solution The descending power of multiple output system.

(when i.e. approximate and speed is maximum), the solution of Convexity Functions, as MIMO when determining that convexity majorized function takes maximum The foundation of system allocation of downlink power.

MIMO power distribution methods based on high specific transfer pre-coding provided in an embodiment of the present invention, according to terminal and base The amplitude and phase of the radio circuit gain stood, determine up channel model and down channel model, according to up channel model With down channel model, it is determined that approximate and speed, in approximate and speed maximum, determines the optimal solution of descending power.According to most It is excellent to de-assign descending power, it is possible to achieve to improve the reasonability of telecommunication system resources configuration, save the communication resource of system.

Optionally, set up radio frequency and mismatch channel model, including：

Set up down channel model when radio frequency is mismatched：

Wherein, H_DIt is down channel model, U_rIt is the radio circuit gain matrix that terminal is received, and U_r=diag { u_r,1, u_r,2,…,u_r,k,…,u_r,K, u_r,kIt is the radio circuit gain that k-th terminal is received, k ∈ [1, K], K is positive integer, and Amplitude obeys logarithm normal distribution δ_u,rIt is default parameter, the amplitude of the radio circuit gain received according to k-th terminal is calculated,Phase is obeyed and is uniformly distributed For k-th terminal is received The maximum of the radio circuit gain-phase distortion arrived,It is common Rayleigh channel, andIn element obey average be 0th, variance is 1 multiple Gauss variable, B_tIt is radio circuit gain matrix that base station sends, and B_t=diag { b_t,1, b_t,2,…,b_t,m,…,b_t,M, b_t,mBe base station send m-th radio circuit gain, m ∈ [1, M], andAmplitude obeys logarithm normal distribution δ_b,tIt is default parameter, the amplitude of the m-th radio circuit gain sent according to base station is calculated,Phase is obeyed and is uniformly distributed For the m that base station sends The maximum of individual radio circuit gain-phase distortion.

Set up up channel model when radio frequency is mismatched：

Wherein, H_UIt is up channel model, B_rIt is the gain matrix of the radio circuit that base station receives, and B_r=diag {b_r,1,b_r,2,…,b_r,m,…,b_r,M, b_r,mIt is m-th radio circuit gain that base station receives, andAmplitude obeys logarithm normal distribution δ_b,tIt is default parameter, the amplitude of the m-th radio circuit gain received according to base station is calculated,Phase is obeyed and is uniformly distributed For base station receive The m maximum of radio circuit gain-phase distortion,ForTransposition, U_tIt is the gain square of the radio circuit that terminal sends Battle array, and U_t=diag { u_t,1,u_t,2,…,u_t,k,…,u_t,K, u_t,kIt is k-th radio circuit gain of terminal transmission, andAmplitude obeys logarithm normal distributionδ_u,rFor Default parameter, the amplitude according to k-th radio circuit gain of terminal transmission is calculated,Phase Obedience is uniformly distributed For k-th radio circuit gain-phase distortion of terminal transmission most Big value.

Due to the mismatch of RF circuits, have It is H_UTransposition.

In embodiments of the present invention, up channel model and down channel model that radio frequency mismatches channel are established, is The follow-up approximate and speed for determining system provides technical support.

Optionally, according to up channel model and down channel model, it is determined that approximate and speed, including：

Step one, according to up channel model and down channel model, by power constraints, determines that high specific is transmitted The precoding coefficients of precoding.

Step 2, according to precoding coefficients, by up channel model and down channel model, determines that the signal of terminal is done Disturb noise ratio.

Step 3, according to Signal Interference and Noise Ratio, passes sequentially through aromatic formula and default approximate formula, it is determined that it is approximate and Speed.

Described approximate formula be can be converted to system and speed herein, the approximate and speed for being easy to solve it is any Approximate formula, for example

In embodiments of the present invention, by up channel model and down channel model, the approximate and speed of system is determined, To determine that corresponding optimal solution provides technical support during approximate and speed maximum.

The requirement equal in order to meet transmission signal power before and after MRT precodings, being calculated precoding coefficients is

Optionally, according to up channel model and down channel model, it is determined that approximate and speed, including：

Obtain the output signal of high specific transmission method precoding in down channel model：Wherein, X is The signal of high specific transfer pre-coding output, α is precoding coefficients, and W is that base station uses high specific transmitting pre-encoding matrix, andD is large scale fading matrix, and D=diag { β₁,β₂,…,β_k,…,β_K, β_kIt is k-th big chi of terminal Degree fading coefficients, P is the transmission power matrix of base station, P=diag { P₁,P₂,…,P_k,…,P_K, P_kFor base station is distributed to k-th The transmission power of terminal, the random signal vector that m ∈ [1, M], S send for terminal, and s=[s₁,s₂,…,s_k,…,s_K], and Meet E (ss^H)=I_K, E (ss^H) it is ss^HExpectation, I_KIt is K rank unit matrixs, s_kIt is k-th random signal vector of terminal, s^H It is the conjugate transposition of S.

Can be obtained by the property of normThereforeCan WillSubstitute into, determine precoding coefficients：

Wherein, power constraints are E (| | x | |²)=p_max, p_maxIt is the maximum transmission power of base station, E (| | x | |²) be | |x||²Expectation, | | x | | represent x norm, M be base station install antenna number, Tr () be trace function, W^HFor the conjugation of W turns Put,It is H_UConjugate matrices.

According to precoding coefficients, up channel model and down channel model, k-th signal of terminal reception is determined：

Wherein, y_kIt is k-th signal of terminal reception, p_kFor base station is k-th power of terminal distribution,h _kForKth OK,It is B_rConjugate matrices, n_kIt is k-th noise of terminal received signals.

Forh _kConjugate transposition,Forh _jConjugate transposition,h _jForJth row, j ∈ [1, K].p_kMeetWherein, p_maxIt is the maximum transmission power of base station,h _kForRow k.

According to k-th signal of terminal reception, k-th Signal Interference and Noise Ratio of terminal is determined：

Wherein, γ_kIt is k-th Signal Interference and Noise Ratio of terminal.

According to Signal Interference and Noise Ratio, by aromatic formula, system and speed are determined：

Wherein, R is system and speed, identifies the sum of each terminal transmission rate,ForExpectation, system and speed are the corresponding total transmission rate of down channel model.

By approximate formulaSystem and speed are converted to closely Sihe speed：

Work as definition When, then the approximate and speed of system can write a Chinese character in simplified form into

Wherein,It is approximate and speed.

In embodiments of the present invention, by up channel model and down channel model, the approximate and speed of system is determined, To determine that approximate and velocity maximum provides technical support.

Optionally, according to approximate and speed, it is the nonconvex property majorized function of target to set up to maximize approximate and speed, bag Include：

According to approximate and speed, as constraints, it is target to set up to maximize approximate and speed to the power with base station Nonconvex property majorized function：

p_k>=0, k=1,2 ..., K

Wherein,For base station is k-th optimal power of terminal distribution,It is the set of optimal power,It is approximate and speed, p_kFor base station is k-th power of terminal distribution, p_maxFor the emission maximum work(that base station can provide Rate, K is positive integer, represents the number of terminal, and s.t. (subject to) refers to the condition for meeting regulation.

In embodiments of the present invention, it is the nonconvex property majorized function of target to set up to maximize approximate and speed, to determine Corresponding optimal solution provides technical support during approximate and velocity maximum.

Optionally, by the way that function to be converted to the formula of convexity by nonconvex property, nonconvex property majorized function is converted into convexity Majorized function, including：

Step one, according to nonconvex property majorized function, by default logarithm lower bound inequality, nonconvex property majorized function is turned It is changed to optimization lower limit function.

Step 2, according to optimization lower limit function, by default exponential transform, determines convexity majorized function.

In embodiments of the present invention, nonconvex property majorized function is converted into convexity majorized function, can calculate it is approximate and The maximum of speed, to determine that corresponding optimal solution provides technical support when approximate and speed takes maximum.

Optionally, by the way that function to be converted to the formula of convexity by nonconvex property, nonconvex property majorized function is converted into convexity Majorized function, including：

According to nonconvex property majorized function, by logarithm lower bound inequality log (1+z) >=λ logz+ μ, nonconvex property is optimized into letter Number is converted to optimization lower limit function：

Wherein,It is optimization lower limit function, λ, μ and z are default parameter, (z₀It is default numerical value) when, in z=z₀Near, inequality log (1+z) >=λ log z+ μ are very tight by what is become, i.e. the lower bound of approximate and speed, will be with optimization lower limit function near specific value Neutralize the value of speed closely.

Therefore the optimization problem can be switched to be converted into the lower bound for optimizing the problem：

p_k>=0, k=1,2 ..., K

In specific value, the optimal solution of this optimization lower limit function is by with the optimal solution of former problem closely.

According to optimization lower limit function, pass throughCan obtain

BecauseAnd μ_kThe parameter unrelated with optimized variable, thus optimization problem convexity majorized function It is equivalent to：

Wherein,ForOptimal solution set.

The function obtained by conversion above is convex function, becauseBe it is linear,It is on { p_k, k=1,2 ..., K joint convex functions.

In embodiments of the present invention, nonconvex property majorized function is converted into convexity majorized function, can calculate it is approximate and The maximum of speed, to determine that corresponding optimal solution provides technical support when approximate and speed takes maximum.

Optionally, when approximate and speed is maximum, the optimal solution of convexity majorized function is determined, including：

Step A, obtains and according to the initial value of base station power allocation result, calculates parameter in convexity majorized function Value.

Step B, the value of parameter substitutes into convexity majorized function in the convexity majorized function that will be calculated, and determines that convexity optimizes letter Several most new explanations.

Step C, according to the most new explanation of convexity majorized function, by exponential transform, determines the most new explanation after exponential transform.

Step D, determines whether the most new explanation after exponential transform meets default stop condition, if newest after exponential transform Solution meets stop condition, then using exponential transform after most new explanation as optimal solution, if the most new explanation after exponential transform is unsatisfactory for stopping Only condition, then the numerical value of the initial value of base station power allocation result is updated to the numerical value of the most new explanation after exponential transform, is returned Step A is continued executing with, until the most new explanation of convexity majorized function meets stop condition.

Set the initial value of base station power allocation resultPower is divided equally in such as base station for K userAccording toWhereinCalculate initial valueSolution is above-mentioned Formula (1) obtains most new explanationThen, according to variation relation formulaCalculateAgain According to what is obtainedUndated parameter λ₁,λ₂,...,λ_K, until meeting stop condition.

Stop condition during the present invention is implemented is any condition for meeting the embodiment of the present invention, including but not It is limited to：Reach the iterative steps specified or meet the condition of convergence | | pⁿ⁺¹-pⁿ| | ＜ ε, wherein, ε is the limits of error,

In embodiments of the present invention, optimal solution is determined by iterative method, the optimal solution can make approximate and speed maximum, profit With the optimal descending power for de-assigning mimo system, it is possible to increase the reasonability of telecommunication system resources configuration, system is saved The communication resource.

The mismatch of RF circuits is degrading the performance of system, while being also that the precoding of base station and the distribution of power bring New challenge.For TDD system, even if carrying out reciprocity calibration research in terms of hardware circuit and software algorithm two to ensure letter The symmetry in road, also how difficulty accomplishes the preferable mismatch problem for eliminating RF circuits, and reciprocity calibration reliability is low.Therefore, have very much Resource allocation problem when necessity is evaluated and mismatched based on channel when channel has non-matching to the performance of system, Such as base station power assignment problem.Need to be calibrated using reciprocity.

In embodiments of the present invention, first, mismatched in RF on the basis of the channel model set up, analyze MRT precodings The approximate and speed that system is obtained.Then, according to the constraints of base station transmitting power, establish to maximize whole system speed Rate sum is the optimization problem of target.Due to the nonconvex property of object function, the embodiment of the present invention employs a kind of logarithm lower bound not Equation, the lower bound of logarithm lower bound inequality will become very tight near specific value, thus the targeted transformation of optimization is optimization Its lower bound.Because the power that object function after conversion distributes to each user is still non-convex, and then a kind of exponential transform is used, The problem is converted into the difference of a linear function and logarithmic function so that the problem becomes convex problem；Finally by iteration Update, with Step wise approximation optimal solution.

Referring to Fig. 2, Fig. 2 is descending for the multi-input multi-output system based on high specific transfer pre-coding of the embodiment of the present invention Another schematic flow sheet of power distribution method, including：

S201, according to the user side and base station side RF radio circuit gain models of investigation, sets up and includes wireless channel and RF Circuit is in interior whole communication channel model.

Set up the unmatched up channel model of RF circuits and the unmatched down channel model of RF circuits.

S202, based on communication channel model, does according to the signal that MRT method for precoding sets up down channel user's receiving terminal Noise ratio is disturbed, and then goes out the approximate and speed of system according to aromatic formula analysis.

First, determine the precoding coefficients of MRT precodings to meet base station transmitting power requirement.Secondly, according to precoding Coefficient, writes out k-th signal of user's reception, so that it is determined that k-th SINR of user (Signal to Interference Plus Noise Ratio, Signal Interference and Noise Ratio).Then, according to k-th SINR of user, using aromatic formula system System is traveled through and speed.Finally by approximate formulaDetermine that system is obtained approximate And speed.

S203, the approximate and speed based on system obtained above, with base station power as constraints, sets up optimization problem It is as follows：

p_k>=0, k=1,2 ..., K

S204, by logarithm lower bound inequality and exponential transform, convex asking is converted into by above-mentioned optimization problem by non-convex problem Topic is solved.

By logarithm lower bound inequality log (1+z) >=λ log z+ μ, can obtain

By exponential transformCan obtain

The convexity majorized function of optimization problem is equivalent to：

S205, undated parameter is iterated, until reaching the condition of convergence.

According to the result of S204, iterative parameter is updatedUntil reaching the condition of convergence：As pre-set Iterative steps meet the condition of convergence | | pⁿ⁺¹-pⁿ| | ＜ ε, wherein, ε is the limits of error. AndThe upper right corner N when representing nth iteration, the parameter of acquisition, such asWhen representing nth iteration, the p being calculated_K。

MIMO power distribution methods based on high specific transfer pre-coding provided in an embodiment of the present invention, according to terminal and base The amplitude and phase of the radio circuit gain stood, determine up channel model and down channel model, according to up channel model With down channel model, it is determined that approximate and speed, in approximate and speed maximum, determines the optimal solution of descending power.According to most It is excellent to de-assign descending power, it is possible to achieve to improve the reasonability of telecommunication system resources configuration, save the communication resource of system.

Referring to Fig. 3, Fig. 3 is the schematic flow sheet that optimal solution is determined by alternative manner of the embodiment of the present invention, including：

S301, initiation parameter.

Initialize the power distribution result of base stationAccording toMeter Calculate initial valueThen basisInitial value is calculated againStep-up error limits ε, and iterative steps n:=1.

S302, the convexity majorized function after solution conversion.

The power distribution result of the base station according to S301, calls CVX kits solution formula (1) problem, obtains optimal solution

S303, seeks the solution of former problem.

According to mapping relationsAnd the solution in S302, calculate the solution of former problem

S304, updates iterative parameter.

According toIn S303As a result, undated parameterThen basisUndated parameter

S305, judges stopping criterion for iteration.

Whether iteration result changes twice before and after judging, if being unsatisfactory for the condition of convergence：||pⁿ⁺¹-pⁿ| | ＜ ε, return S302 continues iteration；If being unsatisfactory for the condition of convergence, terminate iteration and perform S306.

S306, output result.

The optimal power distribution result of output, i.e. optimal solutionFor base station precoding.

In embodiments of the present invention, by the method for iteration, optimal solution is determined, the numerical value of optimal solution is more accurate.

The curve map that Fig. 4 changes for the system and speed of the embodiment of the present invention with antenna for base station number, in Base Transmitter total work When rate is under the conditions of p=20dB, to give channel perfect match and channel mismatch (i.e. user side and base station side is all mismatched) Two kinds of situations, base station is distributed using constant power and uses the MIMO based on high specific transfer pre-coding provided in an embodiment of the present invention Power distribution method carries out the variation diagram of system and speed during power distribution with antenna for base station number.

Wherein, using the system and rate curve of constant power distribution when curve 1 is the perfect match for calculating, curve 2 is Using the system and rate curve of the inventive method during the perfect match for calculating, during the perfect match that curve 3 is obtained for emulation The system and rate curve distributed using constant power, be using the inventive method during the perfect match that curve 4 is obtained for emulation System and rate curve, curve 5 are that the radio circuit for calculating mismatches the system and speed distributed using constant power when variance is 0.3 Rate curve, curve 6 is that system and speed that the radio circuit for calculating is mismatched when variance is 0.3 using the inventive method are bent Line, curve 7 mismatches the system and rate curve distributed using constant power when variance is 0.3 for the radio circuit that emulation is obtained, Curve 8 is the system and rate curve when the radio circuit mismatch variance that emulation is obtained is 0.3 using the inventive method.

Can be obtained by Fig. 4, in the range of whole analogue system, system velocity when channel is preferable will be mismatched better than channel When system velocity, that is, the mismatch of channel can deteriorate the performance of system.And either in channel perfect matchWhen, or it is 0.3 to occur mismatch radio circuit mismatching variance in channel When, when the method provided using the embodiment of the present invention is allocated to the power of base station, system is obtained and rate capability Being better than base station carries out constant power distribution is obtained and speed.In addition, with the increase of antenna for base station number, being carried using the present invention The system and speed value added that the method for confession is obtained also are being improved constantly.When antenna for base station number is smaller, when increasing antenna number, make It is more obvious that the method provided with the embodiment of the present invention obtains speed value added.

The curve map that Fig. 5 changes for the system and speed of the embodiment of the present invention with base station transmitting power, in antenna for base station number M Under the conditions of=100, (i.e. user side and base station side is all mismatched) two kinds of feelings when giving channel perfect match and channel mismatch Condition, the speed that system is obtained when base station is distributed using constant power and carries out power distribution using the algorithm that the present invention is provided is with base station The variation diagram of transmission power.

Wherein, the system and rate curve distributed using constant power when curve 11 is the perfect match for calculating, curve 12 For calculate perfect match when using the inventive method system and rate curve, curve 13 is the perfect match that obtains of emulation The system and rate curve of Shi Caiyong constant powers distribution, using the inventive method during the perfect match that curve 14 is obtained for emulation System and rate curve, curve 15 are that the radio circuit for calculating mismatches the system distributed using constant power when variance is 0.3 And rate curve, curve 16 is that the radio circuit for calculating mismatches the system and speed using the inventive method when variance is 0.3 Rate curve, curve 17 mismatches the system and speed distributed using constant power when variance is 0.3 for the radio circuit that emulation is obtained Curve, curve 18 is bent using the system and speed of the inventive method when the radio circuit mismatch variance that emulation is obtained is 0.3 Line.

Can be obtained by Fig. 5, in the range of whole simulation base station transmission power, the system and speed when channel is preferable are still all excellent The mismatch of system velocity when channel is mismatched, i.e. channel can deteriorate the performance of system.And in channel perfect match and letter Road is mismatchedIn the case of two kinds, the power of base station is divided using method provided by the present invention Timing, performance that system is obtained and speed is better than base station carries out constant power distribution is obtained and speed.In addition, with base The increase of transmission power of standing, speed value added is obtained using method provided by the present invention, and what is also become is increasing.But work as base Transmission power of standing inherently than it is larger when, be further added by transmission power, the rate capability for being obtained of system is not almost just further added by. This be due to when base station transmission power than it is larger when, due to being MRT precodings, the interference between user also begins to become big, institute It is unobvious with the speed increase of system.Therefore, when base station transmitting power than it is larger when, because the interference between user becomes big, use The speed value added that the method for the present invention is obtained also will not always become big, and be held in certain fixed value.

Fig. 6 mismatches the curve map of variance change for the system and speed radio frequency circuit of the embodiment of the present invention, in base station Transmission power be p=20dB under the conditions of, investigated antenna for base station number two kinds of situations of M=100 and M=200, when channel mismatch When, the speed that system is obtained when the method provided using the present invention carries out power distribution with user side/base station side RF circuits not With varianceVariation diagram.

Can be obtained by Fig. 6, when RF circuits mismatch variance become larger when, system and speed constantly reduce.I.e. channel is got over Mismatch, system and rate capability it is poorer.This is because channel gets over mismatch, uncertain bigger, the MRT precodings of channel Effect is poorer, thus the performance of system is poorer.Additionally, in the case of two kinds of M=100 and M=200, using of the invention real When applying the method that example provided the power of base station being allocated, the rate capability that system is obtained is better than constant power distribution and obtains Performance.However, when RF circuits mismatch variance and gradually increase, being entered to the power of base station using algorithm provided by the present invention Row distribution and constant power distribute obtained poor performance and are then gradually reduced.Meanwhile, when RF circuits mismatch variance is larger, M= In the case of two kinds of 100 and M=200, the speed that system is obtained also becomes closer to, i.e., when RF circuits mismatch variance is larger, Increasing the advantage of antenna for base station number will also tend to disappearing.This be due to when RF circuits mismatch variance it is larger when, channel it is not true It is qualitative also bigger, when channel becomes not knowing completely, the power of base station is allocated and is increased the work that antenna number rises With all will be very limited.

The curve map that Fig. 7 changes for the rate gain of the embodiment of the present invention with base station transmitting power, in antenna for base station number M= When 100, when having investigated channel perfect match and channel mismatchFour kinds of situations, base station is adopted The rate gain that system is obtained when being distributed with constant power and carrying out power distribution using the method that the present invention is provided is with Base Transmitter The variation diagram of power.

Wherein, rate gain is defined as：

Can be obtained by Fig. 7, when base station power is smaller, increase transmission power, obtain change in gain obvious.Work as base station Transmission power than it is larger when, then increase transmission power and obtain gain and then tend towards stability.Reason has been carried in analyzing Fig. 5 Arrive, when base station power than it is larger when, the interference of user also than larger, be now further added by transmission power system and speed increase then What is become is very limited, thus gain increase starts to slow down.In addition, when RF circuits mismatch variance ratio is smaller, the speed of system Gain is more than rate gain during perfect match channel, but when mismatching variance and being very big, the rate gain of system is small on the contrary Rate gain when perfect match channel.This be due to when there is serious mismatch in channel, system and rate capability open Begin to deteriorate, and because the uncertain change of channel is big, even if will also become using the speed value added that method proposed by the present invention is obtained Very little, thus gain diminishes.

Fig. 8 mismatches the curve map of variance change for the rate gain radio frequency circuit of the embodiment of the present invention, in base station hair Power is penetrated under the conditions of p=20dB, to have investigated antenna for base station number M=100, in the case of 200,300 three kinds, is provided using the present invention Method carry out the rate gain that system during power distribution obtains and mismatch variance with user side/base station side RF circuitsVariation diagram.

Can be obtained by Fig. 8, in M=100, under 200,300 three kinds of values, rate gain all mismatches variance and first becomes with RF circuits Reduce after big.Namely when RF circuits mismatch variance ratio is smaller, the rate gain obtained using method proposed by the present invention It is obvious.But when there is serious mismatch in channel, that is, not when not knowing completely of channel change, the side proposed with we Method is compared, and constant power distribution method can also obtain good performance.Additionally, when antenna for base station number than it is larger when, obtain optimal speed RF circuits corresponding during rate gain mismatch variance yields also than larger.

Referring to Fig. 9, Fig. 9 shows for the MIMO power distribution systems based on high specific transfer pre-coding of the embodiment of the present invention It is intended to, including：

Channel model builds module 901, the terminal of the multi-input multi-output system for obtaining high specific transfer pre-coding With the amplitude and phase of the radio circuit gain of base station, up channel model and down channel model are determined.

First computing module 902, for according to up channel model and down channel model, it is determined that approximate and speed, its In, approximate and speed is obtained by default approximate formula, the corresponding overall transmission rate of down channel model.

Second computing module 903, is the non-convex of target for according to approximate and speed, setting up to maximize approximate and speed Property majorized function.

3rd computing module 904, for the formula by the way that function to be converted to convexity by nonconvex property, letter is optimized by nonconvex property Number is converted into convexity majorized function.

Power configuration module 905, for when approximate and speed is maximum, determining the optimal solution of convexity majorized function, according to Optimal solution configures the descending power of multi-input multi-output system.

MIMO power distribution systems based on high specific transfer pre-coding provided in an embodiment of the present invention, according to terminal and base The amplitude and phase of the radio circuit gain stood, determine up channel model and down channel model, according to up channel model With down channel model, it is determined that approximate and speed, in approximate and speed maximum, determines the optimal solution of descending power.According to most It is excellent to de-assign descending power, it is possible to achieve to improve the reasonability of telecommunication system resources configuration, save the communication resource of system.

It should be noted that the system of the embodiment of the present invention is to be applied to the above-mentioned MIMO based on high specific transfer pre-coding The system of power distribution method, then all embodiments of the above-mentioned MIMO power distribution methods based on high specific transfer pre-coding are equal Suitable for the device, and can reach same or analogous beneficial effect.

Optionally, channel model builds module 901, including：

Down channel model buildings submodule, for setting up down channel model when radio frequency is mismatched：

Wherein, H_DIt is down channel model, U_rIt is the radio circuit gain matrix that terminal is received, and U_r=diag { u_r,1, u_r,2,…,u_r,k,…,u_r,K, ur, k are k-th radio circuit gain that terminal is received, k ∈ [1, K], and K is positive integer, and Amplitude obeys logarithm normal distribution Phase is obeyed and is uniformly distributed For

The maximum of the radio circuit gain-phase distortion that k-th terminal is received,It is common Rayleigh channel, andIn element to obey average be multiple Gauss variable that 0, variance is 1, B_tIt is the radio circuit gain matrix that base station sends, B_t=diag { b_t,1,b_t,2,…,b_t,m,…,b_t,M, b_t,mBe base station send m-th radio circuit gain, m ∈ [1, M], andAmplitude obeys logarithm normal distributionPhase is obeyed and is uniformly distributed For m-th maximum of radio circuit gain-phase distortion that base station sends.

Up channel model buildings submodule, for setting up up channel model when radio frequency is mismatched：

Wherein, H_UIt is up channel model, B_rIt is the gain matrix of the radio circuit that base station receives, and B_r= diag{b_r,1,b_r,2,…,b_r,m,…,b_r,M, b_r,mFor m-th radio circuit that base station receives increases Benefit, andAmplitude obeys logarithm normal distributionPhase is obeyed and is uniformly distributed It is m-th maximum of radio circuit gain-phase distortion that base station receives,ForTransposition, U_tFor terminal sends Radio circuit gain matrix, and U_t=diag { u_t,1,u_t,2,…,u_t,k,…,u_t,K, u_t,kIt is penetrating for k-th terminal transmission Frequency circuit gain, andAmplitude obeys logarithm normal distributionPhase is obeyed and is uniformly distributed It is k-th maximum of the radio circuit gain-phase distortion of terminal transmission.

In embodiments of the present invention, up channel model and down channel model that radio frequency mismatches channel are established, is The follow-up approximate and speed for determining system provides technical support.

Optionally, the first computing module 902, including：

Precoding coefficients determination sub-module, for according to up channel model and down channel model, by power constraint Condition, determines the precoding coefficients of high specific transfer pre-coding.

Signal to Interference plus Noise Ratio determination sub-module, for according to precoding coefficients, by up channel model and down channel model, Determine the Signal Interference and Noise Ratio of terminal.

Approximate and speed determination sub-module, for according to Signal Interference and Noise Ratio, passing sequentially through aromatic formula and default Approximate formula, it is determined that approximate and speed.

In embodiments of the present invention, by up channel model and down channel model, the approximate and speed of system is determined, To determine that corresponding optimal solution provides technical support during approximate and speed maximum.

Optionally, precoding coefficients determination sub-module specifically for：

Obtain the output signal of high specific transmission method precoding in down channel model：Wherein, X is The signal of high specific transfer pre-coding output, α is precoding coefficients, and W is that base station uses high specific transmitting pre-encoding matrix, andD is large scale fading matrix, and D=diag { β₁,β₂,…,β_k,…,β_K, β_kIt is k-th large scale of terminal Fading coefficients, P is the transmission power matrix of base station, P=diag { P₁,P₂,…,P_k,…,P_K, P_kFor base station distributes to k-th eventually The transmission power at end, the random signal vector that S sends for terminal, and s=[s₁,s₂,…,s_k,…,s_K], and meet E (ss^H)= I_K, E (ss^H) it is ss^HExpectation, I_KIt is K rank unit matrixs, s_kIt is k-th random signal vector of terminal.

According to output signal, by power constraints and down channel model, precoding coefficients are determined：

Wherein, power constraints are E (| | x | |²)=p_max, E (| | x | |²) it is | | x | |²Expectation, M be base station install Antenna number.

Optionally, Signal to Interference plus Noise Ratio determination sub-module specifically for：

According to precoding coefficients, up channel model and down channel model, k-th signal of terminal reception is determined：

Wherein, y_kIt is k-th signal of terminal reception, p_kFor base station is k-th power of terminal distribution,h _kForKth OK,It is B_rConjugate matrices, n_kIt is k-th noise of terminal received signals.

According to k-th signal of terminal reception, k-th Signal Interference and Noise Ratio of terminal is determined：

Wherein, γ_kIt is k-th Signal Interference and Noise Ratio of terminal.

Optionally, approximate and speed determination sub-module specifically for：

According to Signal Interference and Noise Ratio, by aromatic formula, system and speed are determined：

Wherein, R is system and speed, identifies the sum of each terminal transmission rate,ForExpectation；

By approximate formulaSystem and speed are converted to closely Sihe speed：

Wherein,It is approximate and speed, b and c is the coefficient of setting, and

In embodiments of the present invention, by up channel model and down channel model, the approximate and speed of system is determined, To determine that corresponding optimal solution provides technical support during approximate and speed maximum.

Optionally, the second computing module 903, including：

Nonconvex property majorized function setting up submodule, for according to approximate and speed, the power with base station to be built as constraints Vertical is the nonconvex property majorized function of target to maximize approximate and speed：

p_k>=0, k=1,2 ..., K

Wherein,For base station is k-th optimal power of terminal distribution,It is the set of optimal power,It is approximate and speed, p_kFor base station is k-th power of terminal distribution, p_maxFor the emission maximum work(that base station can provide Rate, K is positive integer, represents the number of terminal.

In embodiments of the present invention, it is the nonconvex property majorized function of target to set up to maximize approximate and speed, to determine Corresponding optimal solution provides technical support during approximate and velocity maximum.

Optionally, the 3rd computing module 904, including：

Optimization lower limit function setting up submodule, for according to nonconvex property majorized function, by default logarithm lower bound Formula, optimization lower limit function is converted to by nonconvex property majorized function.

Convexity majorized function setting up submodule, for according to optimization lower limit function, by default exponential transform, determining convex Property majorized function.

In embodiments of the present invention, nonconvex property majorized function is converted into convexity majorized function, can calculate it is approximate and The maximum of speed, to determine that corresponding optimal solution provides technical support when approximate and speed takes maximum.

Optionally, optimization lower limit function setting up submodule specifically for：

According to nonconvex property majorized function, by logarithm lower bound inequality log (1+z) >=λ log z+ μ, nonconvex property is optimized Function is converted to optimization lower limit function：

Wherein,It is optimization lower limit function, λ, μ and z are default parameter,

Optionally, convexity majorized function setting up submodule specifically for：

According to optimization lower limit function, pass throughDetermine convexity majorized function：

Wherein,ForOptimal solution set.

In embodiments of the present invention, nonconvex property majorized function is converted into convexity majorized function, can calculate it is approximate and The maximum of speed, to determine that corresponding optimal solution provides technical support when approximate and speed takes maximum.

Optionally, power configuration module 905, including：

Parameter computation module, for obtaining and according to the initial value of base station power allocation result, calculates convexity optimization The value of parameter in function.

First most new explanation calculating sub module, the value for parameter in the convexity majorized function that will calculate substitutes into convexity optimization Function, determines the most new explanation of convexity majorized function.

Second most new explanation calculating sub module, for the most new explanation according to convexity majorized function, by exponential transform, it is determined that referring to Most new explanation after transformation of variables.

Optimal solution output sub-module, for determining whether the most new explanation after exponential transform meets default stop condition, if Most new explanation after exponential transform meets stop condition, then using exponential transform after most new explanation as optimal solution, if after exponential transform Most new explanation be unsatisfactory for stop condition, then the numerical value of the initial value of base station power allocation result is updated to after exponential transform most The numerical value of new explanation, return parameters calculating sub module is continued executing with, until the most new explanation of convexity majorized function meets stop condition.

In embodiments of the present invention, optimal solution is determined by iterative method, the optimal solution can make approximate and speed maximum, profit With the optimal descending power for de-assigning mimo system, it is possible to increase the reasonability of telecommunication system resources configuration, system is saved The communication resource.

It should be noted that herein, such as first and second or the like relational terms are used merely to a reality Body or operation make a distinction with another entity or operation, and not necessarily require or imply these entities or deposited between operating In any this actual relation or order.And, term " including ", "comprising" or its any other variant be intended to Nonexcludability is included, so that process, method, article or equipment including a series of key elements not only will including those Element, but also other key elements including being not expressly set out, or also include being this process, method, article or equipment Intrinsic key element.In the absence of more restrictions, the key element limited by sentence "including a ...", it is not excluded that Also there is other identical element in process, method, article or equipment including the key element.

Each embodiment in this specification is described by the way of correlation, identical similar portion between each embodiment Divide mutually referring to what each embodiment was stressed is the difference with other embodiment.Especially for system reality Apply for example, because it is substantially similar to embodiment of the method, so description is fairly simple, related part is referring to embodiment of the method Part explanation.

Presently preferred embodiments of the present invention is the foregoing is only, is not intended to limit the scope of the present invention.It is all Any modification, equivalent substitution and improvements made within the spirit and principles in the present invention etc., are all contained in protection scope of the present invention It is interior.

Claims (10)

1. a kind of MIMO power distribution methods based on high specific transfer pre-coding, are applied to MRT high specific transfer pre-codings Extensive MIMO multi-input multi-output system radio frequency mismatch cases, it is characterised in that including：

Obtain high specific transfer pre-coding multi-input multi-output system terminal and base station radio circuit gain amplitude and Phase, determines up channel model and down channel model；

According to the up channel model and the down channel model, it is determined that approximate and speed, wherein, the approximate and speed Obtained by default approximate formula, the corresponding overall transmission rate of the down channel model；

According to the approximate and speed, it is the nonconvex property majorized function of target to set up to maximize the approximate and speed；

By the way that function to be converted to the formula of convexity by nonconvex property, the nonconvex property majorized function is converted into convexity optimization letter Number；

When the approximate and speed is maximum, the optimal solution of the convexity majorized function is determined, institute is configured according to the optimal solution State the descending power of multi-input multi-output system.

2. method according to claim 1, it is characterised in that the determination up channel model and down channel model, Including：

Set up down channel model when radio frequency is mismatched：

$H_D=U_r\tilde{H}{B}_t$

Wherein, the H_DIt is the down channel model, the U_rIt is the radio circuit gain matrix that the terminal is received, andThe u_r,kIt is the radio circuit gain that k-th terminal is received, K ∈ [1, K], K are positive integer, andIt is describedIt is right that amplitude is obeyed Number normal distributionIt is describedPhase is obeyed and is uniformly distributedIt is describedIt is the maximum of the radio circuit gain-phase distortion that k-th terminal is received, It is describedIt is common Rayleigh channel, and it is describedIn element to obey average be multiple Gauss variable that 0, variance is 1, institute State B_tIt is the radio circuit gain matrix that the base station sends, B_t=diag { b_t,1,b_t,2,…,b_t,m,…,b_t,M, it is described b_t,mIt is m-th radio circuit gain that the base station sends, m ∈ [1, M], M is positive integer, andIt is describedAmplitude obeys logarithm normal distributionIt is describedPhase is obeyed and is uniformly distributedIt is describedFor m-th maximum of radio circuit gain-phase distortion that the base station sends Value；

Set up up channel model when radio frequency is mismatched：

$H_U=B_r{\tilde{H}}^T{U}_t$

Wherein, the H_UIt is the up channel model, the B_rIt is the gain matrix of the radio circuit that the base station receives, and B_r=diag { b_r,1,b_r,2,…,b_r,m,…,b_r,M, the b_r,mIt is m-th radio circuit gain that the base station receives, andIt is describedAmplitude obeys logarithm normal distribution It is describedPhase is obeyed and is uniformly distributedIt is describedIt is the base M-th maximum of radio circuit gain-phase distortion that station receives, it is describedFor describedTransposition, the U_tFor institute The gain matrix of the radio circuit of terminal transmission is stated, andThe u_t,kFor K-th radio circuit gain of terminal transmission, andIt is describedIt is right that amplitude is obeyed Number normal distributionIt is describedPhase is obeyed and is uniformly distributedIt is describedIt is k-th maximum of the radio circuit gain-phase distortion of terminal transmission.

3. method according to claim 1, it is characterised in that described according to the up channel model and the descending letter Road model, it is determined that approximate and speed, including：

According to the up channel model and the down channel model, by power constraints, determine that the high specific is passed The precoding coefficients of defeated precoding；

According to the precoding coefficients, by the up channel model and the down channel model, the terminal is determined Signal Interference and Noise Ratio；

According to the Signal Interference and Noise Ratio, aromatic formula and the approximate formula are passed sequentially through, determine the approximate and speed.

4. method according to claim 2, it is characterised in that described according to the up channel model and the descending letter Road model, it is determined that approximate and speed, including：

Obtain the output signal of high specific transmission method precoding in the down channel model：Wherein, the X It is the signal of high specific transfer pre-coding output, the α is the precoding coefficients, and the W is the base station using most Greatly than transmitting pre-encoding matrix, andThe D is large scale fading matrix, and D=diag { β₁,β₂,…, β_k,…,β_K, the β_kIt is k-th large scale fading coefficients of terminal, the P is the transmission power matrix of the base station, P= diag{P₁,P₂,…,P_k,…,P_K, the P_kFor k-th transmission power of terminal is distributed in the base station, the s is the end Hold the random signal vector for sending, and s=[s₁,s₂,…,s_k,…,s_K], and meet E (ss^H)=I_K, E (ss^H) it is ss^HPhase Hope, the I_KIt is K rank unit matrixs, the s_kIt is k-th random signal vector of terminal；

According to the output signal, by power constraints and the down channel model, the precoding coefficients are determined：

$\mathrm{\α}=\frac{1}{\sqrt{{\mathrm{Me}}^{2{\mathrm{\δ}}_{u,t}^2+2{\mathrm{\δ}}_{b,r}^2}{\mathrm{\Σ}}_{k=1}^K{\mathrm{\β}}_k}}$

Wherein, the power constraints are E (| | x | |²)=p_max, the p_maxIt is the maximum transmission power of the base station, it is described E(||x||²) | | the x | | that is described²Expectation, the M is the antenna number that the base station is installed；

According to the precoding coefficients, the up channel model and the down channel model, k-th terminal reception is determined Signal：

$y_k=\mathrm{\α}\sqrt{p_k}{\mathrm{\β}}_k{u}_{r,k}{u}_{t,k}{\underset{\‾}{h}}_k{B}_t{B}_r^{*}{\underset{\‾}{h}}_k^H{s}_k+\mathrm{\α}\sqrt{{\mathrm{\β}}_k}\underset{j=1,j\≠k}{\overset{K}{\mathrm{\Σ}}}\sqrt{{\mathrm{\β}}_j{p}_j}{u}_{r,k}{u}_{t,j}{\underset{\‾}{h}}_k{B}_t{B}_r^{*}{\underset{\‾}{h}}_j^H{s}_j+n_k$

Wherein, the y_kIt is the signal that k-th terminal is received, the p_kFor the base station is k-th work(of terminal distribution Rate, it is describedh _kFor describedRow k, it is describedIt is the B_rConjugate matrices, the n_kIt is k-th terminal received signals Noise；

According to the signal that k-th terminal is received, k-th Signal Interference and Noise Ratio of terminal is determined：

${\mathrm{\γ}}_k=\frac{p_k{\mathrm{\β}}_k^2\left|u_{r,k}{|}^2\right|u_{t,k}{|}^2|{\underset{\‾}{h}}_k{B}_t{B}_r^{*}{\underset{\‾}{h}}_k^H{|}^2}{{\mathrm{\β}}_k\underset{j=1,j\≠k}{\overset{K}{\Σ}}{\mathrm{\β}}_j{p}_j\left|u_{r,k}{|}^2\right|u_{t,j}{|}^2|{\underset{\‾}{h}}_k{B}_t{B}_r^{*}{\underset{\‾}{h}}_j^H{|}^2+{\mathrm{\α}}^{-2}}$

Wherein, the γ_kIt is the Signal Interference and Noise Ratio of k-th terminal；

According to the Signal Interference and Noise Ratio, by the aromatic formula, the system and speed are determined：

$R=E\[\underset{k=1}{\overset{K}{\Σ}}{\log }_2(1+{\mathrm{\γ}}_k)\]$

Wherein, the R is the system and speed, identifies the sum of the transmission rate of each terminal, describedExpectation；

By the approximate formulaThe system and speed are converted into institute State approximate and speed：

${\tilde{R}}^{MRT}\≈\underset{k=1}{\overset{K}{\Σ}}{\log }_2(1+\frac{{\mathrm{bp}}_k{\mathrm{\β}}_k^2}{{\mathrm{\β}}_k{\mathrm{\Σ}}_{j=1,j\≠k}^K{\mathrm{\β}}_j{p}_j+c})$

Wherein, it is describedIt is the approximate and speed, the b and c is the coefficient of setting, and

5. method according to claim 4, it is characterised in that described according to the approximate and speed, sets up to maximize The approximate and speed is the nonconvex property majorized function of target, including：

According to the approximate and speed, the power with the base station is set up to maximize the approximate and speed as constraints It is the nonconvex property majorized function of target：

$\left\{p_k^{opt}\right\}=\arg max{\tilde{R}}^{MRT}$

$\begin{array}{cc}s.t.& {\mathrm{\Σ}}_{k=1}^K{p}_k\≤p_{max}\end{array}$

p_k>=0, k=1,2 ..., K

Wherein, it is describedIt is described for the base station is k-th optimal power of terminal distributionIt is the optimal work( The set of rate, it is describedIt is the approximate and speed, the p_kIt is described for the base station is k-th power of terminal distribution p_maxIt is the maximum transmission power that the base station can provide, the K is positive integer, represents the number of the terminal.

6. method according to claim 1, it is characterised in that the public affairs by the way that function to be converted to convexity by nonconvex property Formula, convexity majorized function is converted into by the nonconvex property majorized function, including：

According to the nonconvex property majorized function, by default logarithm lower bound inequality, by nonconvex property majorized function conversion It is optimization lower limit function；

According to the optimization lower limit function, by default exponential transform, the convexity majorized function is determined.

7. method according to claim 5, it is characterised in that the public affairs by the way that function to be converted to convexity by nonconvex property Formula, convexity majorized function is converted into by the nonconvex property majorized function, including：

It is by logarithm lower bound inequality log (1+z) >=λ logz+ μ, the nonconvex property is excellent according to the nonconvex property majorized function Change function and be converted to optimization lower limit function：

${\tilde{R}}^{MRT}\≥{\tilde{R}}_{LB}^{MRT}=\underset{k=1}{\overset{K}{\mathrm{\Σ}}}\left\{{\mathrm{\λ}}_k\[\log \left({\mathrm{b\β}}_k^2{p}_k\right)-\log \left({\mathrm{\β}}_k{\mathrm{\Σ}}_{j=1,j\≠k}^K{\mathrm{\β}}_j{p}_j+c\right)\]+{\mathrm{\μ}}_k\right\}$

Wherein, it is describedIt is the optimization lower limit function, the λ, the μ and the z are default parameter,

$z_k=\frac{{\mathrm{bp}}_k{\mathrm{\β}}_k^2}{{\mathrm{\β}}_k{\mathrm{\Σ}}_{j=1,j\≠k}^K{\mathrm{\β}}_j{p}_j+c};$

According to the optimization lower limit function, pass throughDetermine the convexity majorized function：

$\left\{{\tilde{p}}_k^{opt}\right\}=\arg max\underset{k=1}{\overset{K}{\mathrm{\Σ}}}\[{\mathrm{\λ}}_k{\tilde{p}}_k-{\mathrm{\λ}}_{k}ln({\mathrm{\β}}_k{\mathrm{\Σ}}_{j=1,j\≠k}^K{\mathrm{\β}}_j{e}^{{\tilde{p}}_j}+c)\]$

$\begin{array}{cc}s.t.& {\mathrm{\Σ}}_{k=1}^K{e}^{{\tilde{p}}_k}\≤p_{max}\end{array}$

Wherein, it is describedForOptimal solution set.

8. method according to claim 1, it is characterised in that in the approximate and speed maximum, determine the convexity The optimal solution of majorized function, including：

Step A, obtains and according to the initial value of base station power allocation result, calculates parameter in the convexity majorized function Value；

Step B, the value of parameter substitutes into the convexity majorized function in the convexity majorized function that will be calculated, and determines described convex The most new explanation of property majorized function；

Step C, according to the most new explanation of the convexity majorized function, by exponential transform, determines the most new explanation after exponential transform；

Step D, determines whether the most new explanation after the exponential transform meets default stop condition, if after the exponential transform Most new explanation meets the stop condition, then using the exponential transform after most new explanation as the optimal solution, if the index becomes Most new explanation after changing is unsatisfactory for the stop condition, then the numerical value of the initial value of the base station power allocation result is updated into institute The numerical value of the most new explanation after exponential transform is stated, return to step A is continued executing with, until the most new explanation of the convexity majorized function meets The stop condition.

9. a kind of MIMO power distribution systems based on high specific transfer pre-coding, are applied to MRT high specific transfer pre-codings Extensive MIMO multi-input multi-output system radio frequency mismatch cases, it is characterised in that including：

Channel model builds module, terminal and base station for obtaining the multi-input multi-output system of high specific transfer pre-coding The amplitude and phase of radio circuit gain, determine up channel model and down channel model；

First computing module, for according to the up channel model and the down channel model, it is determined that approximate and speed, its In, it is described it is approximate obtained by default approximate formula with speed, the corresponding overall transmission rate of the down channel model；

Second computing module, it is non-as target to maximize the approximate and speed for according to described approximate and speed, setting up Convexity majorized function；

3rd computing module, for the formula by the way that function to be converted to convexity by nonconvex property, by the nonconvex property majorized function It is converted into convexity majorized function；

Power configuration module, for when the approximate and speed is maximum, determining the optimal solution of the convexity majorized function, according to The optimal solution configures the descending power of the multi-input multi-output system.

10. system according to claim 9, it is characterised in that first computing module, including：

Precoding coefficients determination sub-module, for according to the up channel model and the down channel model, by power Constraints, determines the precoding coefficients of the high specific transfer pre-coding；

Signal to Interference plus Noise Ratio determination sub-module, for according to the precoding coefficients, by the up channel model and described descending Channel model, determines the Signal Interference and Noise Ratio of the terminal；

Approximate and speed determination sub-module, for according to the Signal Interference and Noise Ratio, pass sequentially through aromatic formula and it is described closely Like formula, the approximate and speed is determined.