CN109905917A - Based on wireless resource allocation methods in the NOMA communication system for wirelessly taking energy - Google Patents

Based on wireless resource allocation methods in the NOMA communication system for wirelessly taking energy Download PDF

Info

Publication number
CN109905917A
CN109905917A CN201910065711.XA CN201910065711A CN109905917A CN 109905917 A CN109905917 A CN 109905917A CN 201910065711 A CN201910065711 A CN 201910065711A CN 109905917 A CN109905917 A CN 109905917A
Authority
CN
China
Prior art keywords
user
cluster
power
iteration
channel
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201910065711.XA
Other languages
Chinese (zh)
Other versions
CN109905917B (en
Inventor
张军
刘晓光
蔡曙
蔡艳
王海荣
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Nanjing Post and Telecommunication University
Original Assignee
Nanjing Post and Telecommunication University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Nanjing Post and Telecommunication University filed Critical Nanjing Post and Telecommunication University
Priority to CN201910065711.XA priority Critical patent/CN109905917B/en
Publication of CN109905917A publication Critical patent/CN109905917A/en
Application granted granted Critical
Publication of CN109905917B publication Critical patent/CN109905917B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02DCLIMATE CHANGE MITIGATION TECHNOLOGIES IN INFORMATION AND COMMUNICATION TECHNOLOGIES [ICT], I.E. INFORMATION AND COMMUNICATION TECHNOLOGIES AIMING AT THE REDUCTION OF THEIR OWN ENERGY USE
    • Y02D30/00Reducing energy consumption in communication networks
    • Y02D30/70Reducing energy consumption in communication networks in wireless communication networks

Landscapes

  • Mobile Radio Communication Systems (AREA)

Abstract

The invention discloses propose to include the following steps: 1) to calculate separately user uplink pilot frequency sequence length, regularization parameter, power division coefficient and power partition coefficient using iterative method based on " system spectral efficiency maximization " principle based on wireless resource allocation methods in the NOMA communication system for wirelessly taking energy;2) the system spectral efficiency absolute value of the difference after current and preceding an iteration is calculated, and judges whether the value is less than or equal to restrain decision threshold;If so, stopping calculating, iteration result at this time is required optimal value, otherwise continues iteration;3) user realizes the optimum allocation of system wireless resource using optimal uplink pilot sequence length, optimal regularization parameter, optimal power allocation coefficient and optimal power division coefficient.This method has the advantages that each user can be made using the information decoding at optimal power division coefficient coordination reception end and collection of energy.

Description

Based on wireless resource allocation methods in the NOMA communication system for wirelessly taking energy
Technical field
The present invention relates to wireless communication technology fields, more particularly to based on wirelessly take can NOMA communication system in wirelessly provide Source distribution method.
Background technique
With the fast development of the 5th third-generation mobile communication technology, it is impossible to meet futures for traditional orthogonal multiple access techniques The connection of 5G network mass users, therefore non-orthogonal multiple (Non-orthogonal Multiple Access, NOMA) technology The main reason for coming into being, using NOMA technology is that it can provide clothes in same time slot and frequency resource for multiple users Business, it is again effective while support can largely connect user to save frequency spectrum resource, however, will be brought using NOMA more dry It disturbs, while the fairness of user must be taken into consideration.Another common-denominator target is exactly to improve efficiency to the maximum extent in the following 5G network, In traditional data transmission procedure, the wireless energy signal of Base Transmitter is all taken as useless energy, causes the waste of resource.Closely Occurs a kind of new technology for capableing of simultaneous transmission wireless messages and energy year a bit, i.e., wirelessly taking can communicate.Wirelessly takes and can communicate It is one of the key technology for realizing green communications.It wirelessly takes to communicate and is not only adapted to low power applications, be also applied for high-power disappear Consume scene.But research shows that energy with receiving end energy collected in information simultaneous transmission is that this disappears with its achievable rate The relationship of that length, therefore, how reasonable distributing radio resource is decoded with the collection of energy and information at coordination reception end to Guan Chong It wants.
Summary of the invention
In view of the above-mentioned problems, the invention proposes based on radio resource allocation side in the NOMA communication system for wirelessly taking energy Method, this method are based on " system spectral efficiency maximization " principle, are ensuring system user fairness and effective use user's collection energy Under the constraint of amount, optimal uplink pilot sequence length, regularization parameter, power division coefficient and power partition coefficient are calculated, With the optimum allocation of radio resource in realization system, the specific implementation steps are as follows:
Step (101): it is based on " system spectral efficiency maximization " principle, calculates separately user uplink pilot frequency sequence using iterative method Column length, regularization parameter, power division coefficient and power partition coefficient, and the ascending pilot frequency sequence gone out according to each iteration Column length, regularization parameter, power division coefficient and power partition coefficient calculate the system spectral efficiency of current iteration number;
Step (102): the system spectral efficiency absolute value of the difference after calculating current and preceding an iteration, and whether judge the value If the determination result is YES then stopping calculating less than or equal to convergence decision threshold, iteration result at this time is required optimal value, If judging result is no, return step (101) continuation iteration;
Step (103): utilizing above-mentioned optimal pilot sequence length, and user sends ascending pilot frequency, and base station is according to receiving Pilot signal carries out channel estimation, obtains the imperfect channel state information of all users, and then combine optimal regularization parameter Regularization force zero pre-coding matrix is calculated, base station is based on above-mentioned optimal power allocation coefficient and sends a signal to each user, each The a part for receiving signal is used for collection of energy based on optimal power division coefficient by user, and another part is decoded for information.
Further, above-mentioned based on wireless resource allocation methods in the NOMA communication system for wirelessly taking energy, in step (101) in, iterative algorithm the following steps are included:
Step (201): initialization system user uplink pilot sequence length, regularization parameter, power division coefficient and Power partition coefficient, respectively symbolically areβ(0)Concurrently set convergence decision threshold ∈=10-8, repeatedly Generation number i=0;
Step (202): it utilizesβ(0)Computing system spectrum efficiency initial value;
Step (203): being based on " system spectral efficiency maximization " principle, calculates the system spectral efficiency of i+1 time iteration
Step (204): judgementIt is whether true, if the determination result is YES, (205) are thened follow the steps, If judging result be it is no, then follow the steps (206);
Step (205): the uplink pilot sequence length that i+1 time iteration goes outRegularization parameter β(i+1), power point Cut coefficientAnd power partition coefficientAs required optimal value, is wirelessly transferred according to above-mentioned optimal value;
Step (206): executing i=i+1, and return step (203) continues iteration.
Further, above-mentioned based on wireless resource allocation methods in the NOMA communication system for wirelessly taking energy, in step (202) in, steps are as follows for the specific calculating of system spectral efficiency initial value:
Step (301): using system user uplink pilot sequence length, regularization parameter, power partitioning parameters it is initial Value, the channel calculated at this time estimates note parameter, to obtain the imperfect channel state information of all users of system;
Step (302): it using the non-ideal estimation channel matrix and regularization parameter initial value of strong user, calculates at this time Send pre-coding matrix;
Step (303): pre-coding matrix is sent using user's real channel vector sum in cluster, user's is effective in calculating cluster It is interfered between power and cluster;
Step (304): interference and uplink pilot sequence length, canonical between the effective power and cluster of user in cluster are utilized Change the initial value of parameter, power partitioning parameters, power distribution parameter, computing system spectrum efficiency initial value.
Further, above-mentioned based on wireless resource allocation methods in the NOMA communication system for wirelessly taking energy, in step (301) in, the imperfect channel state information of all users of acquisition system is implemented as follows:
Based on energy a part that in the single cell multi-user NOMA downlink wireless communication system for wirelessly taking energy, user is collected It is consumed for the circuit of uplink channel estimation and channel estimation phase, another part is used for transmission the circuit of information phase Consumption;Assuming that base station is furnished with N root antenna, user's single antenna, total number of users is K (K > N), wherein > indicate to be greater than;There is N in system A cluster, there are two users for each cluster, and wherein channel gain is preferably strong user, indicate that channel gain is poor to be with (n, 1) Weak user indicates that n ∈ [1,2 ..., N], ∈ expression belong to (n, 2);The uplink that each user sends to base station is assumed simultaneously Pilot number is Tt, channel coherence interval is T, and T is fixed value;Channel usage amount approximately constant, TtChannel usage amount for Downlink channels estimation, T-TtChannel usage amount be used for down link data information transmission, multiple Gauss noise power be σ2; Utilize formula EN, i=(T-Tt)ω(1-psN, i ρCalculate the ENERGY E that i-th of user collects in n-th of clusterN, i, wherein ω ∈ [0, 1] energy conversion efficiency, p are indicatedsIndicate that the power division coefficient of each user, ρ are DL SNR ratio;It will Substitute into calculating formula EN, iIn, it calculatesWhen user collect energyUtilize formula The uplink transmission power for calculating the i-th user in n-th of cluster isWhereinFor with The circuit at family consumes, andAnd p0For constant;It utilizesCalculate the channel estimation parameter of the i-th user in n-th of cluster τN, i, willIt is updated to calculating formula τN, i, calculateWhen channel estimation parameterTo ask ? The estimation channel gain of Shi Suoyou user, wherein subscript x(i)Indicate result of the x Jing Guo i-th iteration;
Assuming that the real channel gain of i-th of user is modeled as in n-th of clusterWherein [1,2] i ∈,Indicate the real channel gain in base station and n-th of cluster between i-th of user, vector magnitude is 1 × N;βN, iFor constant, Indicate the large-scale fading factor of i-th of user in n-th of cluster,Indicate that the multiple Gauss random vector of 1 × N, element all take From 0 mean value, varianceIndependent same distribution;The estimation channel gain of user is available in base station, and all users send to base station Row pilot tone, base station receives ascending pilot frequency and carries out channel estimation, by channel estimation parameterSubstitute into formula Obtain the estimation channel gain of the i-th user in n-th of clusterVector magnitude is 1 ×N;Wherein TN, iIt is the certainty nonnegative definite matrix of N × N, indicates the transmission correlation matrix of antenna for base station, and TN, i=IN, INIt is N The unit matrix of × N,Indicate that the multiple Gauss random vector of 1 × N, element all obey 0 mean value, varianceIt is independent with point Cloth, τN, iFor the channel estimation parameter of the i-th user in n-th of cluster, the accuracy of channel estimation, τ are indicatedN, i∈ [0,1];Subscript (·)HThe conjugate transposition operation of representing matrix,Indicate arithmetic square root,The square root calculation of representing matrix, ()2Table Show quadratic power;And it is modeled with the rapid fading real channel gain of strong user in cluster and weak user as follows:
Wherein n ∈ [1,2 ..., N],WithVector magnitude be all 1 × N,Indicate that the multiple Gauss of 1 × N is random Vector, element all obey 0 mean value, varianceIndependent same distribution,WithIndependently of each other, θnFor constant, and θn∈ [0, 1], error coefficient is indicated.
Further, above-mentioned based on wireless resource allocation methods in the NOMA communication system for wirelessly taking energy, in step (302) it in, calculates transmission pre-coding matrix and is implemented as follows:
Using the estimation channel gain of all strong users, by all strong subscriber channel Gain generating channel estimate matrixs Matrix size is N × N;Transmitting terminal uses the precoding of regularization force zero, and regularization parameter β is utilizedAnd β, it calculates Wherein ()-1The inverse operation of representing matrix;Utilize formulaIt calculates to send and prelist Code matrix G, wherein G meets tr (GGH)≤NP, P > 0, tr () representing matrix ask mark operation.Utilize formula tr (GGH)≤NP, It calculates Wherein ε indicates that the normalized parameter for meeting base station transmitting power constraint, N are antenna for base station number, and P is The general power of Base Transmitter;Enable β=β(0), acquire β=β(0)When pre-coding matrix G.
Further, above-mentioned based on wireless resource allocation methods in the NOMA communication system for wirelessly taking energy, in step (303) it in, interferes and is implemented as follows between the effective power and cluster of user in calculating cluster:
For based on wirelessly take can NOMA communication system,Remove strong user in n-th of cluster and estimates channel vector? It arrivesMatrix size is (N-1) × N;Utilize formulaCalculate strong user and weak user in n-th of cluster Real channel gain be respectivelyWithAccording to G,AndUtilize formula Calculate the effective power U of the strong user of n-th of clusters, utilize formulaCalculate the weak user of n-th of cluster Effective power Uw;Utilize formulaOther clusters are to n-th of cluster in computing system In strong user interference, utilizeOther clusters are to n-th of cluster in computing system The interference of weak user.
Further, above-mentioned based on wireless resource allocation methods in the NOMA communication system for wirelessly taking energy, in step (304) in, computing system spectrum efficiency initial value is implemented as follows:
Utilize following formula
Computing system spectrum efficiency Rsp;Wherein αN, 1For user power distribution coefficient strong in n-th of cluster, constant μ is greater than 0, Ex{f (x) } it indicates to ask expectation to f (x) about variable x, | |2Indicate that the square operation of vector mould, ∑ () indicate summation operation, Log () indicates logarithm operation;It willIt is updated to system spectral efficiency calculating formula RspIn, calculate system Spectrum efficiency initial value
Further, above-mentioned based on wireless resource allocation methods in the NOMA communication system for wirelessly taking energy, in step (203) in, the system spectral efficiency of i+1 time iteration is calculatedIt is implemented as follows: first in i-th iteration result Row pilot sequence lengthRegularization parameter β(i), power division coefficientUtilize following formulas:
Calculate strong user power distribution coefficient in n-th of cluster of i+1 time iterationWherein ρ is DL SNR Than RminIndicate the minimum channel capacity of weak user, and RminFor constant;Secondly, utilizing parameter beta(i)Calculate system Unite spectrum efficiency Rsp, and calculated using linear search optimalThen, parameter is utilizedCalculate system Unite spectrum efficiency Rsp, and β is calculated using linear search(i+1);Then, parameter is utilizedComputing system spectrum effect Rate Rsp, and calculated using linear searchFinally, with computing system spectrum efficiency initial valueEqually, will Substitute into calculating formula RspIn, calculate the system spectral efficiency of i+1 time iteration
Further, above-mentioned based on wireless resource allocation methods in the NOMA communication system for wirelessly taking energy, in step (102) in, maximum number of iterations is set, first judges whether current iteration number is greater than maximum number of iterations, if judging result is It is then to stop calculating;If judging result is no, after calculating current and preceding an iteration system spectral efficiency absolute value of the difference, And judging whether the value is less than or equal to convergence decision threshold and if the determination result is YES then stops calculating, current iteration goes out at this time Uplink pilot sequence length, regularization parameter, power division coefficient and power partition coefficient are required optimal value, if sentencing Disconnected result be it is no, then return step (101) continues iteration, until meeting currently and the system spectral efficiency after preceding an iteration is poor Absolute value is less than or equal to convergence decision threshold.
The invention adopts the above technical scheme compared with prior art, has following technical effect that
(1) this method uses RZF precoding, and compared with ZF precoding, performance is significantly improved, and establishes same cluster In model between strong user and the rapid fading real channel gain of weak user, pass through and adjust error coefficient θn, lifting system performance Computation complexity is reduced simultaneously;
(2) this method is based on " system spectral efficiency maximization " principle, is ensuring system user fairness and effective use use Family is collected under the constraint of energy, is calculated optimal uplink pilot sequence length, regularization parameter, power using iterative method and is divided Coefficient and power partition coefficient make simultaneously so that interference is preferably minimized and efficiently uses the energy of user's collection between system cluster Each user significantly improves system spectral using the information decoding at optimal power division coefficient coordination reception end and collection of energy Efficiency.
Detailed description of the invention
Fig. 1 is that the process of the present invention based on wireless resource allocation methods in the NOMA communication system for wirelessly taking energy is shown It is intended to.
Fig. 2 is the specific implementation step flow chart of iterative algorithm in Fig. 1 step (101).
Fig. 3 is the specific implementation step flow chart that system spectral efficiency initial value is calculated in Fig. 2 step (202).
Specific embodiment
To make the objectives, technical solutions, and advantages of the present invention clearer, below in conjunction with the accompanying drawings and the specific embodiments Present invention is further described in detail.
As shown in Figure 1, of the present invention based on wireless resource allocation methods in the NOMA communication system for wirelessly taking energy, base In the NOMA communication system for wirelessly taking energy, including a base station and multiple users, base station is furnished with multiple antennas, each user Dan Tian Line;This method is ensuring system user fairness and effective use user's collection energy according to " system spectral efficiency maximization " principle Under the constraint of amount, optimal uplink pilot sequence length, regularization parameter, power division coefficient and power distribution system are calculated Number, thus realize the optimum allocation of radio resource, specifically includes the following steps:
Step (101): it is based on " system spectral efficiency maximization " principle, calculates separately user uplink pilot frequency sequence using iterative method Column length, regularization parameter, power division coefficient and power partition coefficient, and the ascending pilot frequency sequence gone out according to each iteration Column length, regularization parameter, power division coefficient and power partition coefficient calculate the system spectral efficiency of current iteration number;
Step (102): the system spectral efficiency absolute value of the difference after calculating current and preceding an iteration, and whether judge the value Less than or equal to convergence decision threshold;If so, stopping calculating, iteration result at this time is required optimal value, otherwise continues to change Generation;
Step (103): utilizing above-mentioned optimal pilot sequence length, and user sends ascending pilot frequency, and base station is according to receiving Pilot signal carries out channel estimation, obtains the imperfect channel state information of all users, and then combine optimal regularization parameter Regularization force zero pre-coding matrix is calculated, base station is based on above-mentioned optimal power allocation coefficient and sends a signal to each user, each The a part for receiving signal is used for collection of energy based on optimal power division coefficient by user, and another part is decoded for information.
Specifically, as shown in Fig. 2, iterative algorithm carries out as follows in step (101):
Step (201): initialization system user uplink pilot sequence length, regularization parameter, power division coefficient and Power partition coefficient, respectively symbolically areConcurrently set convergence decision threshold ∈=10-8、 The number of iterations i=0;
Step (202): it utilizesComputing system spectrum efficiency initial value;Wherein, system spectral is imitated The specific calculating step of rate initial value is as shown in Figure 3:
Step (301): using system user uplink pilot sequence length, regularization parameter, power partitioning parameters it is initial Value, calculates channel estimation parameter at this time, to obtain the imperfect channel state information of all users of system;Specific implementation is such as Under:
Based in the single cell multi-user NOMA downlink wireless communication system for wirelessly taking energy, it is assumed that base station is furnished with N root antenna, User's single antenna, total number of users are K (K > N), wherein > indicate to be greater than;There is N number of cluster in system, there are two users for each cluster, wherein Channel gain is preferably strong user, is indicated with (n, 1), and channel gain is poor for weak user, is indicated with (n, 2), n ∈ [1, 2 ..., N], ∈ expression belongs to;Assume that the ascending pilot frequency quantity that each user sends to base station is T simultaneouslyt, channel coherence interval It is fixed value for T, T;Channel usage amount approximately constant, TtChannel usage amount be used for uplink channel estimation, T-TtChannel Usage amount is used for the transmission of down link data information, and multiple Gauss noise power is σ2;Utilize formula EN, i=(T-Tt)ω(1- psN, iρ calculates the ENERGY E that i-th of user collects in n-th of clusterN, i, wherein ω ∈ [0,1] indicates energy conversion efficiency, psTable Show that the power division coefficient of each user, ρ are DL SNR ratio;It willSubstitute into calculating formula EN, i In, it calculatesWhen user collect energyUtilize formula The uplink transmission power for calculating the i-th user in n-th of cluster isWhereinIt is consumed for the circuit of user, andAnd p0It is normal Amount;It utilizesCalculate the channel estimation parameter τ of the i-th user in n-th of clusterN, i, willIt is updated to calculating Formula τN, i, calculateWhen channel estimation parameterTo acquireShi Suoyou user's Estimate channel gain, wherein subscript x(i)Indicate result of the x Jing Guo i-th iteration;
Assuming that the real channel gain of i-th of user is modeled as in n-th of clusterWherein [1,2] i ∈,Indicate the real channel gain in base station and n-th of cluster between i-th of user, vector magnitude is 1 × N, βN, iFor constant, Indicate the large-scale fading factor of i-th of user in n-th of cluster,Indicate that the multiple Gauss random vector of 1 × N, element all take From 0 mean value, varianceIndependent same distribution;The estimation channel gain of user is available in base station, and all users send to base station Row pilot tone, base station receives ascending pilot frequency and carries out channel estimation, by channel estimation parameterSubstitute into formula Obtain the estimation channel gain of the i-th user in n-th of clusterVector magnitude is 1 ×N;Wherein TN, iIt is the certainty nonnegative definite matrix of N × N, indicates the transmission correlation matrix of antenna for base station, and TN, i=IN, INIt is N The unit matrix of × N,Indicate that the multiple Gauss random vector of 1 × N, element all obey 0 mean value, varianceIt is independent with point Cloth, τN, iFor the channel estimation parameter of the i-th user in n-th of cluster, the accuracy of channel estimation, τ are indicatedN, i∈ [0,1], subscript (·)HThe conjugate transposition operation of representing matrix,Indicate arithmetic square root,The square root calculation of representing matrix, ()2Table Show quadratic power;Rapid fading real channel gain modeling with strong user in cluster and weak user is as follows:
WhereinWithVector magnitude be all 1 × N,Indicate that the multiple Gauss random vector of 1 × N, element all take From 0 mean value, varianceIndependent same distribution,WithIndependently of each other;θnFor constant, and θn∈ [0,1] indicates error coefficient;
Step (302): it using the non-ideal estimation channel matrix and regularization parameter initial value of strong user, calculates at this time Send pre-coding matrix;Wherein, transmission pre-coding matrix is calculated to be implemented as follows:
Using the estimation channel gain of all strong users, by all strong subscriber channel Gain generating channel estimate matrixs Matrix size is N × N, and transmitting terminal uses the precoding of regularization force zero, and regularization parameter β is utilizedAnd β, it calculates Wherein ()-1The inverse operation of representing matrix;Utilize formulaIt calculates to send and prelist Code matrix G, wherein G meets tr (GGH)≤NP, P > 0, tr () representing matrix ask mark operation;Utilize formula tr (GGH)≤NP, It calculates Wherein ε indicates that the normalized parameter for meeting base station transmitting power constraint, N are antenna for base station number, and P is The general power of Base Transmitter;Enable β=β(0), acquire β=β(0)When pre-coding matrix G;
Step (303): pre-coding matrix is sent using user's real channel vector sum in cluster, user's is effective in calculating cluster It is interfered between power and cluster;Wherein, it calculates to interfere between the effective power and cluster of user in cluster and be implemented as follows:
For based on wirelessly take can NOMA communication system,Remove strong user in n-th of cluster and estimates channel vector? It arrivesMatrix size is (N-1) × N;Utilize formulaCalculate strong user and weak user in n-th of cluster Real channel gain be respectivelyWithAccording to G,AndUtilize formula Calculate the effective power U of the strong user of n-th of clusters, utilize formulaCalculate the weak user of n-th of cluster Effective power Uw;Utilize formulaOther clusters are to n-th of cluster in computing system In strong user interference, utilizeOther clusters are to n-th of cluster in computing system The interference of weak user;
Step (304): interference and uplink pilot sequence length, canonical between the effective power and cluster of user in cluster are utilized Change the initial value of parameter, power partitioning parameters, power distribution parameter, computing system spectrum efficiency initial value;Wherein, computing system is composed Efficiency initial value is implemented as follows:
Utilize the user uplink pilot frequency sequence length of i-th iterationRegularization parameter β(i), power division coefficient And power partition coefficientUtilize following formula
The system spectral efficiency of i-th iteration can be calculatedWherein αN, 1It is distributed for user power strong in n-th of cluster and is Number, constant μ are greater than 0, Ex{ f (x) } indicates to ask expectation to f (x) about variable x, | |2Indicate the square operation of vector mould, ∑ () indicates summation operation, and log () indicates logarithm operation;Therefore, willIt is updated to system spectral effect Rate calculating formula RspIn, calculate system spectral efficiency initial value
Step (203): it is based on " system spectral efficiency maximization " principle, first according to i-th iteration result ascending pilot frequency sequence Column lengthRegularization parameter β(i), power division coefficientUtilize following formulas:
Calculate strong user power distribution coefficient in n-th of cluster of i+1 time iterationWherein ρ is DL SNR Than RminIndicate the minimum channel capacity of weak user, and RminFor constant;Secondly, utilizing parameter beta(i)Calculate system Unite spectrum efficiency Rsp, and calculated using linear search optimalThen, parameter is utilizedCalculate system Unite spectrum efficiency Rsp, and β is calculated using linear search(i+1);Then, parameter is utilizedComputing system spectrum effect Rate Rsp, and calculated using linear searchFinally, with computing system spectrum efficiency initial valueEqually, will Substitute into calculating formula RspIn, calculate the system spectral efficiency of i+1 time iteration
Step (204): judgementIt is whether true, if the determination result is YES, (205) are thened follow the steps, If judging result be it is no, then follow the steps (206);
Step (205): the uplink pilot sequence length that i+1 time iteration goes outRegularization parameter β(i+1), power point Cut coefficientAnd power partition coefficientAs required optimal value, is wirelessly transferred according to above-mentioned optimal value;
Step (206): executing i=i+1, and return step (203) continues iteration.
Specifically, in step (102), maximum number of iterations S is set, enables current iteration number s=i+1, judgement is current Whether the number of iterations s is greater than maximum number of iterations S, if the determination result is YES, then stops calculating;If judging result be it is no, count System spectral efficiency absolute value of the difference after calculating current and preceding an iteration, and judge whether the value is less than or equal to restrain decision gate Limit, if the determination result is YES, then stops calculating, the uplink pilot sequence length that current iteration goes out at this timeRegularization parameter β(i+1), power division coefficientAnd in n-th of cluster strong user power partition coefficientAs required optimal value, If judging result is no, return step (101) continuation iteration, until meeting the system spectral efficiency after current and preceding an iteration Absolute value of the difference is less than or equal to convergence decision threshold;
Specifically, in step (103), the calculated system user optimal pilot sequence length of step (102), Yong Hufa Ascending pilot frequency is sent, base station carries out channel estimation according to the pilot signal received, gets the non-ideal communication channel shape of all users State information, and then optimal regularization parameter is combined to calculate regularization force zero pre-coding matrix, base station is based on above-mentioned optimal power point Distribution coefficient sends a signal to each user, and a part for receiving signal is used for energy based on optimal power division coefficient by each user Amount is collected, and another part is decoded for information, to realize the optimum allocation of system wireless resource.
The above description is merely a specific embodiment, but scope of protection of the present invention is not limited thereto, any In the technical scope disclosed by the present invention, any changes or substitutions that can be easily thought of by those familiar with the art, all answers It is included within the scope of protection of the present invention.

Claims (9)

1. based on wireless resource allocation methods in the NOMA communication system for wirelessly taking energy, it is characterised in that: specific step is as follows:
Step (101): it is based on " system spectral efficiency maximization " principle, it is long to calculate separately user uplink pilot frequency sequence using iterative method Degree, regularization parameter, power division coefficient and power partition coefficient;And it is long according to the uplink pilot sequence that each iteration goes out Degree, regularization parameter, power division coefficient and power partition coefficient calculate the system spectral efficiency of current iteration number;
Step (102): the system spectral efficiency absolute value of the difference after calculating current and preceding an iteration, and judge whether the value is less than If the determination result is YES then stop calculating equal to convergence decision threshold, iteration result at this time is required optimal value, if sentencing Disconnected result be no, then return step (101) continuation iteration;
Step (103): utilizing above-mentioned optimal pilot sequence length, and user sends ascending pilot frequency, and base station is according to the pilot tone received Signal carries out channel estimation, obtains the imperfect channel state information of all users, and then optimal regularization parameter is combined to calculate Regularization force zero pre-coding matrix, base station are based on above-mentioned optimal power allocation coefficient and send a signal to each user, each user The a part for receiving signal is used for collection of energy based on optimal power division coefficient, another part is decoded for information.
2. according to claim 1 based on wireless resource allocation methods in the NOMA communication system for wirelessly taking energy, feature Be: in step (101), iterative algorithm the following steps are included:
Step (201): initialization system user uplink pilot sequence length, regularization parameter, power division coefficient and power Distribution coefficient, respectively symbolically areβ(0)Concurrently set convergence decision threshold ∈=10-8, iteration time Number i=0;
Step (202): it utilizesβ(0)Computing system spectrum efficiency initial value;
Step (203): being based on " system spectral efficiency maximization " principle, calculates the system spectral efficiency of i+1 time iteration
Step (204): judgementIt is whether true, if the determination result is YES, (205) are thened follow the steps, if sentencing Disconnected result be it is no, then follow the steps (206);
Step (205): the uplink pilot sequence length that i+1 time iteration goes outRegularization parameter β(i+1), power segmentation system NumberAnd power partition coefficientAs required optimal value, is wirelessly transferred according to above-mentioned optimal value;
Step (206): executing i=i+1, and return step (203) continues iteration.
3. according to claim 2 based on wireless resource allocation methods in the NOMA communication system for wirelessly taking energy, feature Be: in step (202), steps are as follows for the specific calculating of system spectral efficiency initial value:
Step (301): using system user uplink pilot sequence length, regularization parameter, power partitioning parameters initial value, meter Channel estimation parameter at this time is calculated, to obtain the imperfect channel state information of all users of system;
Step (302): using the non-ideal estimation channel matrix and regularization parameter initial value of strong user, transmission at this time is calculated Pre-coding matrix;
Step (303): pre-coding matrix is sent using user's real channel vector sum in cluster, calculates the effective power of user in cluster It is interfered between cluster;
Step (304): joined using interference between the effective power and cluster of user in cluster and uplink pilot sequence length, regularization The initial value of number, power partitioning parameters, power distribution parameter, computing system spectrum efficiency initial value.
4. according to claim 3 based on wireless resource allocation methods in the NOMA communication system for wirelessly taking energy, feature Be: in step (301), the imperfect channel state information of all users of acquisition system is implemented as follows:
Based in the single cell multi-user NOMA downlink wireless communication system for wirelessly taking energy, energy a part that user collects is used for The consumption of the circuit of uplink channel estimation and channel estimation phase, the circuit that another part is used for transmission information phase disappear Consumption;Assuming that base station is furnished with N root antenna, user's single antenna, total number of users is K (K > N), wherein > indicate to be greater than;Have in system N number of Cluster, there are two users for each cluster, and wherein channel gain is preferably strong user, indicate that poor channel gain is weak with (n, 1) User indicates that n ∈ [1,2 ..., N], ∈ expression belong to (n, 2);Assume that each user leads to the uplink that base station is sent simultaneously Frequency amount is Tt, channel coherence interval is T, and T is fixed value;Channel usage amount approximately constant, TtChannel usage amount be used for uplink Link channel estimation, T-TtChannel usage amount be used for down link data information transmission, multiple Gauss noise power be σ2;Benefit With formula EN, i=(T-Tt)ω(1-psN, iρ calculates the ENERGY E that i-th of user collects in n-th of clusterN, i, wherein [0,1] ω ∈ Indicate energy conversion efficiency, psIndicate that the power division coefficient of each user, ρ are DL SNR ratio;It will Substitute into calculating formula EN, iIn, it calculatesWhen user collect energyUtilize formula The uplink transmission power for calculating the i-th user in n-th of cluster isWhereinFor with The circuit at family consumes, andAnd p0For constant;It utilizesCalculate the channel estimation ginseng of the i-th user in n-th of cluster Number τN, i, willIt is updated to calculating formula τN, i, calculateWhen channel estimation parameterTo It acquires The estimation channel gain of Shi Suoyou user, wherein subscript x(i)Indicate result of the x Jing Guo i-th iteration;
Assuming that the real channel gain of i-th of user is modeled as in n-th of clusterWherein [1,2] i ∈, Indicate the real channel gain in base station and n-th of cluster between i-th of user, vector magnitude is 1 × N;βN, iFor constant, indicate The large-scale fading factor of i-th of user in n-th of cluster,Indicate that the multiple Gauss random vector of 1 × N, element all obey 0 Mean value, varianceIndependent same distribution;The estimation channel gain of user is available in base station, and all users send uplink to base station and lead Frequently, base station receives ascending pilot frequency and carries out channel estimation, by channel estimation parameterSubstitute into formula Obtain the estimation channel gain of the i-th user in n-th of clusterVector magnitude is 1 ×N;Wherein TN, iIt is the certainty nonnegative definite matrix of N × N, indicates the transmission correlation matrix of antenna for base station, and TN, i=IN, INIt is N The unit matrix of × N,Indicate that the multiple Gauss random vector of 1 × N, element all obey 0 mean value, varianceIt is independent with point Cloth, τN, iFor the channel estimation parameter of the i-th user in n-th of cluster, the accuracy of channel estimation, τ are indicatedN, i∈ [0,1];Subscript (·)HThe conjugate transposition operation of representing matrix,Indicate arithmetic square root,The square root calculation of representing matrix, ()2Table Show quadratic power;And it is modeled with the rapid fading real channel gain of strong user in cluster and weak user as follows:
Wherein n ∈ [1,2 ..., N],WithVector magnitude be all 1 × N,Indicate 1 × N multiple Gauss at random to Amount, element all obey 0 mean value, varianceIndependent same distribution,WithIndependently of each other, θnFor constant, and θn∈ [0,1], Indicate error coefficient.
5. according to claim 3 based on wireless resource allocation methods in the NOMA communication system for wirelessly taking energy, feature It is: in step (302), calculates transmission pre-coding matrix and be implemented as follows:
Using the estimation channel gain of all strong users, by all strong subscriber channel Gain generating channel estimate matrixsMatrix is big Small is N × N;Transmitting terminal uses the precoding of regularization force zero, and regularization parameter β is utilizedAnd β, it calculates Wherein ()-1The inverse operation of representing matrix;Utilize formulaIt calculates and sends pre-coding matrix G, wherein G Meet tr (GGH)≤NP, P > 0, tr () representing matrix ask mark operation;Utilize formula tr (GGH)≤NP is calculated Wherein ε indicates that the normalized parameter for meeting base station transmitting power constraint, N are antenna for base station number, and P is Base Transmitter General power;Enable β=β(0), acquire β=β(0)When pre-coding matrix G.
6. according to claim 3 based on wireless resource allocation methods in the NOMA communication system for wirelessly taking energy, feature It is: in step (303), interferes and be implemented as follows between the effective power and cluster of user in calculating cluster:
For based on wirelessly take can NOMA communication system,Remove strong user in n-th of cluster and estimates channel vectorIt obtainsMatrix size is (N-1) × N;Utilize formulaCalculate the true of strong user and weak user in n-th of cluster Real channel gain is respectivelyWithAccording to G,AndUtilize formulaIt calculates The effective power U of the strong user of n-th of clusters, utilize formulaCalculate having for the weak user of n-th of cluster Imitate power Uw;Utilize formulaOther clusters are to strong in n-th of cluster in computing system The interference of user utilizesOther clusters are to the weak use of n-th of cluster in computing system The interference at family.
7. according to claim 3 based on wireless resource allocation methods in the NOMA communication system for wirelessly taking energy, feature Be: in step (304), computing system spectrum efficiency initial value is implemented as follows:
Utilize following formula
Computing system spectrum efficiency Rsp;Wherein αN, 1For user power distribution coefficient strong in n-th of cluster, constant μ is greater than 0, Ex{f(x)} It indicates to ask expectation to f (x) about variable x, | |2Indicate that the square operation of vector mould, ∑ () indicate summation operation, log () indicates logarithm operation;It willβ(0)It is updated to system spectral efficiency calculating formula RspIn, calculate system spectral effect Rate initial value
8. according to claim 2 based on wireless resource allocation methods in the NOMA communication system for wirelessly taking energy, feature It is: in step (203), calculates the system spectral efficiency of i+1 time iterationIt is implemented as follows: first with i-th Iteration result uplink pilot sequence lengthRegularization parameter β(i), power division coefficientUtilize following formulas:
Calculate strong user power distribution coefficient in n-th of cluster of i+1 time iterationWherein ρ is DL SNR ratio, RminIndicate the minimum channel capacity of weak user, and RminFor constant;Secondly, utilizing parameter beta(i)Computing system spectrum Efficiency Rsp, and calculated using linear search optimalThen, parameter is utilizedComputing system spectrum Efficiency Rsp, and β is calculated using linear search(i+1);Then, parameter is utilizedβ(i+1),Computing system spectrum efficiency Rsp, And it is calculated using linear searchFinally, with computing system spectrum efficiency initial valueEqually, willβ(i+1) Substitute into calculating formula RspIn, calculate the system spectral efficiency of i+1 time iteration
9. according to claim 1 based on wireless resource allocation methods in the NOMA communication system for wirelessly taking energy, feature It is: in step (102), sets maximum number of iterations, first judge whether current iteration number is greater than maximum number of iterations, if sentencing Disconnected result be it is yes, then stop calculating;If judging result be it is no, calculate currently and preceding an iteration after system spectral efficiency difference Absolute value, and judge whether the value is less than or equal to convergence decision threshold and if the determination result is YES then stops calculating, currently at this time Uplink pilot sequence length, regularization parameter, power division coefficient and the power partition coefficient that iteration goes out are required optimal Value, if judging result is no, return step (101) continuation iteration, until meeting the system spectral after current and preceding an iteration Inefficient absolute value is less than or equal to convergence decision threshold.
CN201910065711.XA 2019-01-23 2019-01-23 Wireless energy-carrying NOMA communication system wireless resource allocation method Active CN109905917B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910065711.XA CN109905917B (en) 2019-01-23 2019-01-23 Wireless energy-carrying NOMA communication system wireless resource allocation method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910065711.XA CN109905917B (en) 2019-01-23 2019-01-23 Wireless energy-carrying NOMA communication system wireless resource allocation method

Publications (2)

Publication Number Publication Date
CN109905917A true CN109905917A (en) 2019-06-18
CN109905917B CN109905917B (en) 2022-12-02

Family

ID=66944124

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910065711.XA Active CN109905917B (en) 2019-01-23 2019-01-23 Wireless energy-carrying NOMA communication system wireless resource allocation method

Country Status (1)

Country Link
CN (1) CN109905917B (en)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110430576A (en) * 2019-07-17 2019-11-08 大连理工大学 The non-orthogonal multiple access safety in network of unmanned plane auxiliary is taken can communication means
CN110881010A (en) * 2019-11-15 2020-03-13 南京邮电大学 Statistical CSI-assisted multi-user NOMA downlink transmission method
CN113014304A (en) * 2021-02-22 2021-06-22 南京邮电大学 Position design method for unmanned aerial vehicle-assisted relay multi-user wireless communication
CN117395626A (en) * 2023-12-11 2024-01-12 厦门大学 Underwater acoustic network water quality monitoring data collection method based on meta learning and NOMA

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2013000169A1 (en) * 2011-06-29 2013-01-03 中国人民解放军理工大学 Resource allocation method for maximizing throughput in cooperative cognitive simo network
CN108234101A (en) * 2018-01-11 2018-06-29 郑州航空工业管理学院 Efficiency maximizes pilot signal design method and large-scale multi-antenna system
CN108430040A (en) * 2017-02-14 2018-08-21 王晋良 Federated user is grouped with power distribution method and uses the base station of the method
CN109005551A (en) * 2018-07-10 2018-12-14 南京邮电大学 A kind of multi-user's NOMA downlink power distributing method of imperfect channel state information

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2013000169A1 (en) * 2011-06-29 2013-01-03 中国人民解放军理工大学 Resource allocation method for maximizing throughput in cooperative cognitive simo network
CN108430040A (en) * 2017-02-14 2018-08-21 王晋良 Federated user is grouped with power distribution method and uses the base station of the method
CN108234101A (en) * 2018-01-11 2018-06-29 郑州航空工业管理学院 Efficiency maximizes pilot signal design method and large-scale multi-antenna system
CN109005551A (en) * 2018-07-10 2018-12-14 南京邮电大学 A kind of multi-user's NOMA downlink power distributing method of imperfect channel state information

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110430576A (en) * 2019-07-17 2019-11-08 大连理工大学 The non-orthogonal multiple access safety in network of unmanned plane auxiliary is taken can communication means
CN110430576B (en) * 2019-07-17 2023-04-11 大连理工大学 Safe energy-carrying communication method in non-orthogonal multiple access network assisted by unmanned aerial vehicle
CN110881010A (en) * 2019-11-15 2020-03-13 南京邮电大学 Statistical CSI-assisted multi-user NOMA downlink transmission method
CN110881010B (en) * 2019-11-15 2022-01-11 南京邮电大学 Statistical CSI-assisted multi-user NOMA downlink transmission method
CN113014304A (en) * 2021-02-22 2021-06-22 南京邮电大学 Position design method for unmanned aerial vehicle-assisted relay multi-user wireless communication
CN117395626A (en) * 2023-12-11 2024-01-12 厦门大学 Underwater acoustic network water quality monitoring data collection method based on meta learning and NOMA
CN117395626B (en) * 2023-12-11 2024-02-09 厦门大学 Underwater acoustic network water quality monitoring data collection method based on meta learning and NOMA

Also Published As

Publication number Publication date
CN109905917B (en) 2022-12-02

Similar Documents

Publication Publication Date Title
CN109905917A (en) Based on wireless resource allocation methods in the NOMA communication system for wirelessly taking energy
CN109005551B (en) Multi-user NOMA downlink power distribution method of non-ideal channel state information
CN107613567B (en) Wireless sensor network resource allocation method based on wireless energy transmission
CN108770007A (en) Wireless portable communications system Multipurpose Optimal Method based on NOMA
CN111405596B (en) Resource optimization method for large-scale antenna wireless energy-carrying communication system under Rice channel
CN105680920A (en) Method for optimizing throughput of multiuser multi-antenna digital-energy integrated communication network
CN108832977A (en) The sparse nonopiate access implementing method in the extensive airspace MIMO
CN104796900A (en) Cellular network D2D (device-to-device) communication resource distributing method based on auction theory
CN113905443A (en) Spectrum efficiency-based resource allocation method for multi-carrier NOMA system under imperfect CSI
CN106028456A (en) Power allocation method of virtual cell in 5G high density network
CN105721127A (en) Power line channel communication resource allocation method based on minimum required rates of users
CN110881010A (en) Statistical CSI-assisted multi-user NOMA downlink transmission method
CN108462996A (en) A kind of non-orthogonal multiple network resource allocation method
CN106059734B (en) Massive mimo system pilot distribution method based on edge customer interference metric values
Huang et al. Energy efficient resource allocation algorithm in multi-carrier NOMA systems
Li et al. Pilot assignment based on weighted-count for cell-free massive MIMO systems
CN106851726A (en) A kind of cross-layer resource allocation method based on minimum speed limit constraint
Zhang et al. Large system analysis of downlink MISO-NOMA system via regularized zero-forcing precoding with imperfect CSIT
CN108650710B (en) Resource allocation method of wireless energy supply communication network based on hybrid multiple access
CN107426775B (en) Distributed multi-user access method for high-energy-efficiency heterogeneous network
Sun et al. Spectral efficiency and scalability analysis for multi-level cooperative cell-free networks
CN107733488B (en) Water injection power distribution improvement method and system in large-scale MIMO system
Zeng et al. Downlink power allocation optimization in pattern division multiple access
CN115622617A (en) NOMA-based cognitive backscatter communication robust transmission method
Wei et al. Joint user association and resource optimization for multi-cell NOMA networks

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant