CN113423112A - RIS assisted multi-carrier NOMA transmission system parameter optimization method - Google Patents

RIS assisted multi-carrier NOMA transmission system parameter optimization method Download PDF

Info

Publication number
CN113423112A
CN113423112A CN202110675019.6A CN202110675019A CN113423112A CN 113423112 A CN113423112 A CN 113423112A CN 202110675019 A CN202110675019 A CN 202110675019A CN 113423112 A CN113423112 A CN 113423112A
Authority
CN
China
Prior art keywords
user
model
ris
subcarrier
reflection coefficient
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
CN202110675019.6A
Other languages
Chinese (zh)
Other versions
CN113423112B (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.)
Southeast University
Original Assignee
Southeast 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 Southeast University filed Critical Southeast University
Priority to CN202110675019.6A priority Critical patent/CN113423112B/en
Publication of CN113423112A publication Critical patent/CN113423112A/en
Application granted granted Critical
Publication of CN113423112B publication Critical patent/CN113423112B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04W—WIRELESS COMMUNICATION NETWORKS
    • H04W24/00—Supervisory, monitoring or testing arrangements
    • H04W24/02—Arrangements for optimising operational condition
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L41/00—Arrangements for maintenance, administration or management of data switching networks, e.g. of packet switching networks
    • H04L41/14—Network analysis or design
    • H04L41/145—Network analysis or design involving simulating, designing, planning or modelling of a network
    • Y—GENERAL 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
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02D—CLIMATE 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/00—Reducing energy consumption in communication networks
    • Y02D30/70—Reducing energy consumption in communication networks in wireless communication networks

Abstract

The invention discloses a parameter optimization method of an RIS auxiliary multi-carrier NOMA transmission system, which takes the maximum system energy efficiency as a target and respectively establishes a subcarrier distribution model, a user decoding sequence model and a user power and RIS reflection coefficient model; aiming at a subcarrier distribution model, a matching method based on channel gain sequencing is provided; secondly, aiming at a user decoding sequence model in a subcarrier, a low-complexity decoding sequence optimization method is provided; and finally, aiming at the user power and RIS reflection coefficient model, a combined optimization method of alternating iteration of the user power and the RIS reflection coefficient is used. The method of the invention can effectively improve the channel quality and improve the system energy efficiency.

Description

RIS assisted multi-carrier NOMA transmission system parameter optimization method
Technical Field
The invention relates to the technical field of NOMA communication, in particular to a parameter optimization method of a RIS auxiliary multi-carrier NOMA transmission system.
Background
At present, in the application of maximizing the energy efficiency of multi-carrier NOMA downlink transmission assisted by RIS, a joint optimization method for maximizing the energy efficiency aiming at special scenes of two users and without considering a decoding sequence exists, but the research and the application of optimizing subcarrier allocation, multi-user power allocation and RIS reflection coefficient control are less while considering the decoding sequence.
Disclosure of Invention
In view of the above, the present invention provides a method for optimizing parameters of an RIS assisted multi-carrier NOMA transmission system, so as to solve the problems mentioned in the background art. Aiming at an RIS-assisted multi-carrier NOMA downlink transmission system, the invention researches the joint optimization problem of subcarrier allocation, decoding sequence, power allocation and RIS reflection coefficient by taking the maximum system energy efficiency as an objective function and considering the constraint of the maximum transmission power and the RIS reflection coefficient of a base station.
In order to solve the problems, the invention adopts the technical scheme that:
a parameter optimization method for RIS auxiliary multi-carrier NOMA transmission system includes the following steps:
step S1, establishing an energy efficiency maximization model, wherein the energy efficiency maximization model A1 is expressed as:
Figure BDA0003120635130000011
constraint conditions are as follows:
C1:
Figure BDA0003120635130000012
C2:
Figure BDA0003120635130000013
C3:
Figure BDA0003120635130000014
C4:
Figure BDA0003120635130000015
C5:
Figure BDA0003120635130000016
C6:θm∈[0,2π),m＝1,…,M
in the formula, the number of RIS reflecting units is M, and the phase of the mth reflecting unit is thetam,m＝0,1,…,M，
Figure BDA0003120635130000021
Is the reflection coefficient matrix of the RIS, diag [ ·]Representing vector diagonalization; k is the number of the users,
Figure BDA0003120635130000022
representing a set of users; n is the number of sub-carriers,
Figure BDA0003120635130000023
representing a set of subcarriers, specifying that each user can be and can only be allocated to one subcarrier; deltan,kAn indicator variable, delta, for subcarrier allocationn,kWith 1 is meant that subcarrier n is allocated to user k, δn,k0 denotes no partition, δ ═ δ { δ ═ δ1,1,…,δ1,K,…,δn,1,…,δn,K,…,δN,1,…,δN,KAllocating a state set for the sub-carriers; a isn,kDenotes the power allocated to user k on subcarrier n, a ═ a1,1,…,a1,K,…,an,1,…,an,K,…,aN,1,…,aN,KThe user power set is used as the power set; v. ofn,k，fnAnd gn,kChannel matrices representing subcarriers n to user k, subcarriers n to RIS and RIS to user k, respectively, (-)HRepresents a conjugate transpose; on(k) Indicating the decoding order, o, of user k on subcarrier nn(k) J denotes that user k is the j-th decoded user, and o ═ { o ═1(1),…,o1(K),…,on(1),…,on(K),…,oN(1),…,oN(K) Is a decoding order set; sigma2Represents the variance of additive white gaussian noise with a mean of 0; pcTotal circuit power consumed for the system; user k on subcarrier n has rate Rn,kIs shown as
Figure BDA0003120635130000024
| · | represents modulo a complex number; suppose on(k)＜on(j) Rate R at which the signal of user k is decoded at user jn,j→kIs shown as
Figure BDA0003120635130000025
PmaxIs the maximum transmit power at the base station;
step S2, solving the energy efficiency maximization model, and the specific steps are as follows:
step S201, the mth element on the diagonal of the reflection coefficient matrix theta of the initial RIS is:
Figure BDA0003120635130000026
channel matrix of subcarrier n to user k is vn,kThe corresponding channel gain is | vn,kI, n ', k' are
Figure BDA0003120635130000027
The minimum subcarrier number and the user number, and the angle (·) represents the angle corresponding to the complex number; channel gain of user k is
Figure BDA0003120635130000028
User k selection
Figure BDA0003120635130000029
Maximum subcarrier
Figure BDA00031206351300000210
Order to
Figure BDA00031206351300000211
δn,k＝0,n≠n*(ii) a After all users select the subcarrier, obtaining a subcarrier distribution state set delta;
step S202, obtaining the value of the sub-carrier distribution state set delta by step S201, and enabling | v of the user k not distributed to the sub-carrier nn,k0, |; on subcarrier n, according to the method for all users
Figure BDA00031206351300000212
Is sorted from small to large, and the | v of the user k is recordedn,kThe | value is ordered as inThen on(k)＝in(ii) a After all the subcarriers are sequenced, a decoding sequence set o is obtained;
step S203, after the subcarrier distribution state set delta and the decoding sequence set o are obtained, simplifying a constraint C3 of a model A1; the energy efficiency maximization model B1, which transforms model a1 into a reflection coefficient matrix Θ for user power set a and RIS, is:
Figure BDA0003120635130000031
constraint conditions are as follows:
C1:
Figure BDA0003120635130000032
C2:
Figure BDA0003120635130000033
C3:
Figure BDA0003120635130000034
C4:θm∈[0,2π),m＝1,…,M
step S204, iteratively solving the energy efficiency maximization model B1, wherein in the ith iteration, the solving process is as follows:
step S2041, splitting an energy efficiency maximization model B1 to respectively obtain a model for optimizing a user power set a and a model for optimizing a reflection coefficient matrix theta of the RIS; the model B2 for user power set a is expressed as:
Figure BDA0003120635130000035
constraint conditions are as follows:
C1:
Figure BDA0003120635130000036
C2:
Figure BDA0003120635130000037
Θ(i-1)the matrix theta of the reflection coefficient of the RIS is obtained for the i-1 iteration; solving the model B2 by using a Dinkelbach method to obtain a user power set a of the ith iteration(i)；
Step S2042,
Figure BDA0003120635130000038
Is a(i)The n × k-th element of (1); the maximization model B3 for the reflectance matrix Θ for RIS is expressed as:
Figure BDA0003120635130000041
constraint conditions are as follows:
C1:
Figure BDA0003120635130000042
C2:θm∈[0,2π),m＝1,…,M
step S2043, introducing a relaxation variable:
Figure BDA0003120635130000043
η＝[η1,1,…,η1,K,…,ηn,1,…,ηn,K,…,ηN,1,…,ηN,K]；
order to
Figure BDA0003120635130000044
t＝[t1,1,…,t1,K,…,tn,1,…,tn,K,…,tN,1,…,tN,K]，
Figure BDA0003120635130000045
(·)*Represents a conjugate operation, then
Figure BDA0003120635130000046
Introducing vector beta ═ beta1,1,…,β1,K,…,βn,1,…,βn,K,…,βN,1,…,βN,K]Model B3 is transformed into model B4, which model B4 is represented as:
Figure BDA0003120635130000047
constraint conditions are as follows:
C1:|sm|＝1,m＝1,…,M
C2:
Figure BDA0003120635130000048
C3:
Figure BDA0003120635130000049
C4:
Figure BDA00031206351300000410
and step S2044, carrying out iterative solution on the model B4 by using a continuous convex approximation method and a penalty function method, wherein in the u-th iteration, the model B4 is converted into the model B5 as follows:
Figure BDA0003120635130000051
constraint conditions are as follows:
C1:|sm|≤1,m＝1,…,M
C2:
Figure BDA0003120635130000052
C3:
Figure BDA0003120635130000053
C4:
Figure BDA0003120635130000054
wherein
Figure BDA0003120635130000055
And s(u-1)Is the value of the u-1 th iteration,
Figure BDA0003120635130000056
is s(u-1)The mth element of (1); r (·) represents taking a real part;
step S2045, substituting variable values S, eta and beta obtained in the u-th iteration into a model B5 to calculate an objective function, and when the value of the objective function is not changed any more or the absolute value of the change is smaller than a threshold value epsilon1Then, the iteration method is terminated, s is output, and the s is diagonalized to obtain a RIS reflection coefficient matrix theta; otherwise, carrying out the next iteration;
step S205, substituting the user power set a obtained in step S2041 and the RIS reflection coefficient matrix theta obtained in step S2045 in the ith iteration into the model B1 to calculate the value of the objective function, and when the value of the objective function is not changed any more or the absolute value of the change is smaller than the threshold epsilon2Then, the iteration method is terminated, and a user power set a and an RIS reflection coefficient matrix theta are output; otherwise, the next iteration is carried out.
Further, the given threshold value epsilon in step S20451Is 10-4。
Further, a given threshold value ε in step S2052Is 10-4
The invention has the beneficial effects that:
aiming at a multi-carrier NOMA downlink transmission system assisted by RIS, the invention takes the maximum system energy efficiency as an objective function, and performs joint optimization on subcarrier allocation, decoding sequence of users in subcarriers, user power and a reflection coefficient matrix of RIS. A corresponding mathematical optimization model is established for the problem, and the problem is decomposed into three sub-problems of subcarrier allocation, user decoding order optimization and user power and RIS reflection coefficient alternative optimization: aiming at the solution of the subcarrier allocation sub-problem, a matching method based on channel gain sequencing is provided; aiming at the solution of the decoding sequence of the user in the subcarrier, a decoding sequence optimization method with low complexity is provided; aiming at the sub-problem of the joint optimization of the user power and the RIS reflection coefficient, a Dinkelbach method is used for optimizing the user power, a continuous convex approximation method and a penalty function method are used for optimizing the RIS reflection coefficient, and a joint optimization method of alternating iteration of the user power and the RIS reflection coefficient is provided. By adopting the method, the energy efficiency of the system can be effectively and simply improved.
Detailed Description
The technical solutions in the embodiments of the present invention will be clearly and completely described below, and it is obvious that the described embodiments are some, but not all, embodiments of the present invention. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
Example 1
The embodiment provides a parameter optimization method for an RIS-assisted multi-carrier NOMA transmission system, which is used for researching the joint optimization problem of subcarrier allocation, decoding sequence, power allocation and RIS reflection coefficient by taking the maximum system energy efficiency as an objective function and considering the constraint of the maximum transmission power and the RIS reflection coefficient of a base station aiming at an RIS-assisted multi-carrier NOMA downlink transmission system. In order to solve the corresponding non-convex problem, the original problem is decomposed into four sub-problems, the four sub-problems are solved through three steps, firstly, a matching method based on the principle of channel gain maximization is provided, secondly, a user decoding sequence optimization method in sub-carriers with lower complexity is provided, and finally, when the sub-carrier distribution and decoding sequence are given, a Dinkelbach method is used for optimizing user power, and a continuous convex approximation and penalty function method is used for optimizing the RIS reflection coefficient.
Specifically, in this embodiment, in an RIS-assisted multi-carrier NOMA downlink transmission system, an RIS is placed between a base station and users, the base station transmits superimposed signals to K users through a direct channel and a reflected channel via the RIS by using a power domain multiplexing technique according to the NOMA principle, and the users are equipped with single antennas. Note the book
Figure BDA0003120635130000061
Representing a set of users, the number of subcarriers N,
Figure BDA0003120635130000062
representing a set of sub-carriers, each definedA user can and can only be allocated to one subcarrier.
First step, an,kRepresenting the power, s, allocated to user k on subcarrier nn,kThe transmission data representing user k on subcarrier n is an independent random variable with a mean of 0 and a variance of 1. Deltan,kAn indicator variable, delta, for subcarrier allocationn,kWith 1 is meant that subcarrier n is allocated to user k, δn,k0 means no allocation. Transmission signal x of subcarrier nnCan be expressed as:
Figure BDA0003120635130000063
secondly, the number of the reflection units of the RIS is M, and the phase position of the mth reflection unit is recorded as thetamE [0,2 pi)), M1, 2, … M, diagonal matrix
Figure BDA0003120635130000064
Matrix of reflection coefficients, diag [. cndot.) representing the RIS]Representing the vector diagonalization. By the symbol vn,k，
Figure BDA0003120635130000065
And
Figure BDA0003120635130000066
channel matrices representing subcarriers n to user k, subcarriers n to RIS and RIS to user k, respectively, (-)HRepresenting the conjugate transpose of the matrix. c represents that the user receives white Gaussian noise, the mean value of which is 0 and the variance of which is sigma2. The base station is assumed to perfectly acquire the channel state information of all channels, which are assumed to be flat-fading rayleigh channels. The signal y received by user k on subcarrier n, after passing through the direct channel and the reflected channeln,kCan be expressed as:
Figure BDA0003120635130000071
third step, with on(k) Representing on subcarrier nDecoding order of user k, on(k) J indicates that user k is the jth decoded user. While demodulating the signal of user k on subcarrier n, the other non-demodulated users l, on(l)＞on(k) The signal of (b) is treated as interference, the signal to interference plus noise ratio y of the kth user on subcarrier nn,kThe calculation is as follows:
Figure BDA0003120635130000072
where | represents the rate R of the user modulo a complex numbern,kCan be represented as Rn,k＝log2(1+γn,k)。
Fourth, assume on(k)＜on(j) Calculating the velocity of user k at user j is expressed as:
Figure BDA0003120635130000073
fifth step, PcFor the total circuit power consumed by the system, the total power consumed by the RIS assisted NOMA system is calculated as:
Figure BDA0003120635130000074
sixthly, calculating the energy efficiency of the system as follows:
Figure BDA0003120635130000075
to obtain subcarrier allocation, user decoding order within subcarriers, user power and RIS reflection coefficient and apply to the above system, let δ be { δ ═ δ1,1,…,δ1,K,…,δn,1,…,δn,K,…,δN,1,…,δN,KDenotes a subcarrier allocation status set, o ═ o1(1),…,o1(K),…,on(1),…,on(K),…,oN(1),…,oN(K) Is a decoding order set, a ═ a1,1,…,a1,K,…,an,1,…,an,K,…,aN,1,…,aN,KEstablishing an energy efficiency maximization model for a user power set, wherein the energy efficiency maximization model is expressed as follows:
Figure BDA0003120635130000081
constraint conditions are as follows:
C1:
Figure BDA0003120635130000082
C2:
Figure BDA0003120635130000083
C3:
Figure BDA0003120635130000084
C4:
Figure BDA0003120635130000085
C5:
Figure BDA0003120635130000086
C6:θm∈[0,2π),m＝1,…,M
solving the energy efficiency maximization model, which comprises the following specific steps:
step 1: the mth element on the diagonal of the reflection coefficient matrix Θ of the initial RIS is:
Figure BDA0003120635130000087
channel matrix of subcarrier n to user k is vn,kThe corresponding channel gain is | vn,kI, n ', k' are
Figure BDA0003120635130000088
The smallest subcarrier number and user number, angle (·) represents the angle corresponding to the complex number. Channel gain of user k is
Figure BDA0003120635130000089
User k selection
Figure BDA00031206351300000810
Maximum subcarrier
Figure BDA00031206351300000811
Order to
Figure BDA00031206351300000812
And after all users select the subcarrier, obtaining a subcarrier distribution state set delta.
Step 2: obtaining the value of a subcarrier allocation state set delta by the step 1, and enabling users k | v not allocated to the subcarrier nn,kAnd | ═ 0. On subcarrier n, according to the method for all users
Figure BDA00031206351300000813
Is sorted from small to large, and the | v of the user k is recordedn,kThe | value is ordered as inThen on(k)＝in(ii) a And obtaining a decoding sequence set o after sequencing all the subcarriers.
And step 3: after the subcarrier allocation state set delta and the decoding sequence set o are obtained, the constraint C3 of the model A1 is simplified. The energy efficiency maximization model B1, which transforms model a1 into a reflection coefficient matrix Θ for user power set a and RIS, is:
Figure BDA0003120635130000091
constraint conditions are as follows:
C1:
Figure BDA0003120635130000092
C2:
Figure BDA0003120635130000093
C3:
Figure BDA0003120635130000094
C4:θm∈[0,2π),m＝1,…,M
and 4, step 4: and (3) iteratively solving the energy efficiency maximization model B1, wherein in the ith iteration, the solving process is as follows:
step 4-1: and splitting the energy efficiency maximization model B1 to respectively obtain a model for optimizing the user power set a and a model for optimizing the reflection coefficient matrix theta of the RIS. The model B2 for user power set a is expressed as:
Figure BDA0003120635130000095
constraint conditions are as follows:
C1:
Figure BDA0003120635130000096
C2:
Figure BDA0003120635130000097
Θ(i-1)the matrix theta of the reflection coefficient of the RIS obtained for the i-1 th iteration. Solving the model B2 by using a Dinkelbach method to obtain a user power set a of the ith iteration(i)。
Step 4-2:
Figure BDA0003120635130000098
is a(i)The nth × k elements of (a), the maximization model B3 of the reflection coefficient matrix Θ with respect to RIS is expressed as:
Figure BDA0003120635130000099
constraint conditions are as follows:
C1:
Figure BDA00031206351300000910
C2:θm∈[0,2π),m＝1,…,M
step 4-3: introducing a relaxation variable:
Figure BDA0003120635130000101
η＝[η1,1,…,η1,K,…,ηn,1,…,ηn,K,…,ηN,1,…,ηN,K]；
order to
Figure BDA0003120635130000102
t＝[t1,1,…,t1,K,…,tn,1,…,tn,K,…,tN,1,…,tN,K]，
Figure BDA0003120635130000103
(·)*Represents a conjugate operation, then
Figure BDA0003120635130000104
Introducing vectors
β＝[β1,1,…,β1,K,…,βn,1,…,βn,K,…,βN,1,…,βN,K]Model B3 is transformed into model B4, which model B4 is represented as:
Figure BDA0003120635130000105
constraint conditions are as follows:
C1:|sm|＝1,m＝1,…,M
C2:
Figure BDA0003120635130000106
C3:
Figure BDA0003120635130000107
C4:
Figure BDA0003120635130000108
step 4-4: the solution was iterated using the successive convex approximation method and the penalty function method for model B4, and in the u-th iteration model B4 was transformed into model B5 as follows:
Figure BDA0003120635130000109
constraint conditions are as follows:
C1:|sm|≤1,m＝1,…,M
C2:
Figure BDA00031206351300001010
C3:
Figure BDA00031206351300001011
C4:
Figure BDA00031206351300001012
wherein
Figure BDA00031206351300001013
And s(u-1)Is the value of the u-1 th iteration,
Figure BDA00031206351300001014
is s(u-1)R (-) represents the real part.
And 4-5: substituting variable values s, eta and beta obtained in the u-th iteration into the model B5 to calculate the target function, and when the value of the target function is not changed any more or the absolute value of the change is less than the threshold value 10-4Then, the iteration method is terminated, s is output, and the s is diagonalized to obtain a RIS reflection coefficient matrix theta; otherwise, the next iteration is carried out.
And 5: substituting the user power set a obtained in the step 4-1 and the RIS reflection coefficient matrix theta obtained in the step 4-5 in the ith iteration into a model B1 to calculate the value of the objective function, and when the value of the objective function is not changed any more or the absolute value of the change is less than a threshold value 10-4Then, the iteration method is terminated, and a user power set a and an RIS reflection coefficient matrix theta are output; otherwise, the next iteration is carried out.
The invention is not described in detail, but is well known to those skilled in the art.
The foregoing detailed description of the preferred embodiments of the invention has been presented. It should be understood that numerous modifications and variations could be devised by those skilled in the art in light of the present teachings without departing from the inventive concepts. Therefore, the technical solutions available to those skilled in the art through logic analysis, reasoning and limited experiments based on the prior art according to the concept of the present invention should be within the scope of protection defined by the claims.

Claims (3)

1. A parameter optimization method for RIS assisted multi-carrier NOMA transmission system is characterized by comprising the following steps:
step S1, establishing an energy efficiency maximization model, wherein the energy efficiency maximization model A1 is expressed as:
Figure FDA0003120635120000011
constraint conditions are as follows:
Figure FDA0003120635120000012
Figure FDA0003120635120000013
Figure FDA0003120635120000014
Figure FDA0003120635120000015
Figure FDA0003120635120000016
C6：θm∈[0，2π)，m＝1，…，M
in the formula, the number of RIS reflecting units is M, and the phase of the mth reflecting unit is thetam，m＝0，1，…，M，
Figure FDA0003120635120000017
Is the reflection coefficient matrix of the RIS, diag [ ·]Representing vector diagonalization; k is the number of the users,
Figure FDA0003120635120000018
representing a set of users; n is the number of sub-carriers,
Figure FDA0003120635120000019
representing a set of subcarriers, specifying that each user can be and can only be allocated to one subcarrier; deltan，kAn indicator variable, delta, for subcarrier allocationn，kWith 1 is meant that subcarrier n is allocated to user k, δn，k0 denotes no partition, δ ═ δ { δ ═ δ1，1，…，δ1，K，…，δn，1，…，δn，K，…，δN，1，…，δN，KAllocating a state set for the sub-carriers; a isn，kDenotes the power allocated to user k on subcarrier n, a ═ a1，1，…，a1，K，…，an，1，…，an，K，…，aN，1，…，aN，KThe user power set is used as the power set; v. ofn，k，fnAnd gn，kRespectively representing subcarriers n to user k, subcarriers n to RIS and RIS to userk channel matrix, (.)HRepresents a conjugate transpose; on(k) Indicating the decoding order, o, of user k on subcarrier nn(k) J denotes that user k is the j-th decoded user, and o ═ { o ═1(1)，…，o1(K)，…，on(1)，…，on(K)，…，oN(1)，…，oN(K) Is a decoding order set; sigma2Represents the variance of additive white gaussian noise with a mean of 0; pcTotal circuit power consumed for the system; user k on subcarrier n has rate Rn，kIs shown as
Figure FDA0003120635120000021
| · | represents modulo a complex number; suppose on(k)＜on(j) Rate R at which the signal of user k is decoded at user jn，j→kIs shown as
Figure FDA0003120635120000022
PmaxIs the maximum transmit power at the base station;
step S2, solving the energy efficiency maximization model, and the specific steps are as follows:
step S201, the mth element on the diagonal of the reflection coefficient matrix theta of the initial RIS is:
Figure FDA0003120635120000023
channel matrix of subcarrier n to user k is vn，kThe corresponding channel gain is | vn，kI, n ', k' are
Figure FDA0003120635120000024
The minimum subcarrier number and the user number, and the angle (·) represents the angle corresponding to the complex number; channel gain of user k is
Figure FDA0003120635120000025
User k selection
Figure FDA0003120635120000026
Maximum subcarrier n*，
Figure FDA0003120635120000027
Order to
Figure FDA0003120635120000028
δn，k＝0，n≠n*(ii) a After all users select the subcarrier, obtaining a subcarrier distribution state set delta;
step S202, obtaining the value of the sub-carrier distribution state set delta by step S201, and enabling | v of the user k not distributed to the sub-carrier nn，k0, |; on subcarrier n, according to the method for all users
Figure FDA0003120635120000029
Is sorted from small to large, and the | v of the user k is recordedn，kThe | value is ordered as inThen on(k)＝in(ii) a After all the subcarriers are sequenced, a decoding sequence set o is obtained;
step S203, after the subcarrier distribution state set delta and the decoding sequence set o are obtained, simplifying a constraint C3 of a model A1; the energy efficiency maximization model B1, which transforms model a1 into a reflection coefficient matrix Θ for user power set a and RIS, is:
Figure FDA0003120635120000031
constraint conditions are as follows:
Figure FDA0003120635120000032
Figure FDA0003120635120000033
Figure FDA0003120635120000034
C4：θm∈[0，2π)，m＝1，…，M
step S204, iteratively solving the energy efficiency maximization model B1, wherein in the ith iteration, the solving process is as follows:
step S2041, splitting an energy efficiency maximization model B1 to respectively obtain a model for optimizing a user power set a and a model for optimizing a reflection coefficient matrix theta of the RIS; the model B2 for user power set a is expressed as:
Figure FDA0003120635120000035
constraint conditions are as follows:
Figure FDA0003120635120000036
Figure FDA0003120635120000037
Θ(i-1)the matrix theta of the reflection coefficient of the RIS is obtained for the i-1 iteration; solving the model B2 by using a Dinkelbach method to obtain a user power set a of the ith iteration(i)；
Step S2042,
Figure FDA0003120635120000038
Is a(i)The n × k-th element of (1); the maximization model B3 for the reflectance matrix Θ for RIS is expressed as:
Figure FDA0003120635120000041
constraint conditions are as follows:
Figure FDA0003120635120000042
C2：θm∈[0，2π)，m＝1，…，M
step S2043, introduce a relaxation variable
Figure FDA0003120635120000043
η＝[η1，1，…，η1，K，…，ηn，1，…，ηn，K，…，ηN，1，…，ηN，K]；
Order to
Figure FDA0003120635120000044
t＝[t1，1，…，t1，K，…，tn，1，…，tn，K，…，tN，1，…，tN，x]，
Figure FDA0003120635120000045
(·)*Represents a conjugate operation, then
Figure FDA0003120635120000046
Introducing vector beta ═ beta1，1，…，β1，K，…，βn，1，…，βn，K，…，βN，1，…，βN，K]Model B3 is transformed into model B4, which model B4 is represented as:
Figure FDA0003120635120000047
constraint conditions are as follows:
C1：|sm|＝1，m＝1，…，M
Figure FDA0003120635120000048
Figure FDA0003120635120000049
Figure FDA00031206351200000410
and step S2044, carrying out iterative solution on the model B4 by using a continuous convex approximation method and a penalty function method, wherein in the u-th iteration, the model B4 is converted into the model B5 as follows:
Figure FDA0003120635120000051
constraint conditions are as follows:
C1：|sm|≤1，m＝1，…，M
Figure FDA0003120635120000052
Figure FDA0003120635120000053
Figure FDA0003120635120000054
wherein
Figure FDA0003120635120000055
And s(u-1)Is the value of the u-1 th iteration,
Figure FDA0003120635120000056
is s(u-1)The mth element of (1); r (·) represents taking a real part;
step S2045, obtaining the variable value of the u-th iterations, eta, beta are substituted into the model B5 to calculate the objective function when the value of the objective function is no longer changed or the absolute value of the change is less than the threshold value epsilon1Then, the iteration method is terminated, s is output, and the s is diagonalized to obtain a RIS reflection coefficient matrix theta; otherwise, carrying out the next iteration;
step S205, substituting the user power set a obtained in step S2041 and the RIS reflection coefficient matrix theta obtained in step S2045 in the ith iteration into the model B1 to calculate the value of the objective function, and when the value of the objective function is not changed any more or the absolute value of the change is smaller than the threshold epsilon2Then, the iteration method is terminated, and a user power set a and an RIS reflection coefficient matrix theta are output; otherwise, the next iteration is carried out.
2. A RIS assisted multi-carrier NOMA transmission system parameter optimization method according to claim 1, characterized in that, the given threshold value e in step S20451Is 10-4。
3. A RIS assisted multi-carrier NOMA transport system parameter optimization method according to claim 1, characterized by the given threshold value ε in step S2052Is 10-4。
CN202110675019.6A 2021-06-18 2021-06-18 RIS assisted multi-carrier NOMA transmission system parameter optimization method Active CN113423112B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110675019.6A CN113423112B (en) 2021-06-18 2021-06-18 RIS assisted multi-carrier NOMA transmission system parameter optimization method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110675019.6A CN113423112B (en) 2021-06-18 2021-06-18 RIS assisted multi-carrier NOMA transmission system parameter optimization method

Publications (2)

Publication Number Publication Date
CN113423112A true CN113423112A (en) 2021-09-21
CN113423112B CN113423112B (en) 2023-02-10

Family

ID=77788970

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110675019.6A Active CN113423112B (en) 2021-06-18 2021-06-18 RIS assisted multi-carrier NOMA transmission system parameter optimization method

Country Status (1)

Country Link
CN (1) CN113423112B (en)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114466388A (en) * 2022-02-16 2022-05-10 北京航空航天大学 Intelligent super-surface-assisted wireless energy-carrying communication method
CN114499602A (en) * 2021-12-28 2022-05-13 南京邮电大学 Method for designing transmitting precoding and phase shift matrix in RIS auxiliary MIMO-NOMA communication system
CN114614925A (en) * 2022-03-10 2022-06-10 南京航空航天大学 Reconfigurable intelligent surface assisted energy efficiency optimization method in millimeter wave non-orthogonal multiple access system
CN116471613A (en) * 2023-04-28 2023-07-21 南京邮电大学 Energy efficiency optimization method and device for multi-user MIMO system
CN113923767B (en) * 2021-09-23 2023-10-13 怀化建南电子科技有限公司 Energy efficiency maximization method for multi-carrier cooperation non-orthogonal multiple access system

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111953391A (en) * 2020-07-09 2020-11-17 东南大学 Intelligent reflector assisted multi-user MIMO uplink energy efficiency and spectrum efficiency combined optimization method
CN112272384A (en) * 2020-11-03 2021-01-26 广东工业大学 Communication system throughput optimization method based on reconfigurable intelligent surface
CN112601242A (en) * 2020-12-17 2021-04-02 南京邮电大学 Intelligent reflecting surface assisted two-cell NOMA uplink low-power-consumption transmission method
CN112822703A (en) * 2021-02-03 2021-05-18 广东工业大学 Intelligent reflecting surface assisted performance gain optimization method for non-orthogonal multiple access system
CN112929068A (en) * 2021-02-04 2021-06-08 重庆邮电大学 SDR-based IRS-NOMA system beam forming optimization method

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111953391A (en) * 2020-07-09 2020-11-17 东南大学 Intelligent reflector assisted multi-user MIMO uplink energy efficiency and spectrum efficiency combined optimization method
CN112272384A (en) * 2020-11-03 2021-01-26 广东工业大学 Communication system throughput optimization method based on reconfigurable intelligent surface
CN112601242A (en) * 2020-12-17 2021-04-02 南京邮电大学 Intelligent reflecting surface assisted two-cell NOMA uplink low-power-consumption transmission method
CN112822703A (en) * 2021-02-03 2021-05-18 广东工业大学 Intelligent reflecting surface assisted performance gain optimization method for non-orthogonal multiple access system
CN112929068A (en) * 2021-02-04 2021-06-08 重庆邮电大学 SDR-based IRS-NOMA system beam forming optimization method

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
YUN XU等: "Fair Non-Orthogonal Multiple Access Communication Systems with Reconfigurable Intelligent Surface", 《2020 IEEE 31ST ANNUAL INTERNATIONAL SYMPOSIUM ON PERSONAL, INDOOR AND MOBILE RADIO COMMUNICATIONS》 *

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113923767B (en) * 2021-09-23 2023-10-13 怀化建南电子科技有限公司 Energy efficiency maximization method for multi-carrier cooperation non-orthogonal multiple access system
CN114499602A (en) * 2021-12-28 2022-05-13 南京邮电大学 Method for designing transmitting precoding and phase shift matrix in RIS auxiliary MIMO-NOMA communication system
CN114499602B (en) * 2021-12-28 2023-03-07 南京邮电大学 Transmission precoding and phase shift matrix design method in RIS auxiliary MIMO-NOMA communication system
CN114466388A (en) * 2022-02-16 2022-05-10 北京航空航天大学 Intelligent super-surface-assisted wireless energy-carrying communication method
CN114466388B (en) * 2022-02-16 2023-08-08 北京航空航天大学 Intelligent super-surface-assisted wireless energy-carrying communication method
CN114614925A (en) * 2022-03-10 2022-06-10 南京航空航天大学 Reconfigurable intelligent surface assisted energy efficiency optimization method in millimeter wave non-orthogonal multiple access system
CN114614925B (en) * 2022-03-10 2023-11-24 南京航空航天大学 Energy efficiency optimization method in reconfigurable intelligent-surface-assisted millimeter wave non-orthogonal multiple access system
CN116471613A (en) * 2023-04-28 2023-07-21 南京邮电大学 Energy efficiency optimization method and device for multi-user MIMO system
CN116471613B (en) * 2023-04-28 2023-10-10 南京邮电大学 Energy efficiency optimization method and device for multi-user MIMO system

Also Published As

Publication number Publication date
CN113423112B (en) 2023-02-10

Similar Documents

Publication Publication Date Title
CN113423112B (en) RIS assisted multi-carrier NOMA transmission system parameter optimization method
CN112865893B (en) Intelligent reflector assisted SM-NOMA system resource allocation method
He et al. Model-driven deep learning for physical layer communications
CN1983910B (en) Method for detecting signal in multi-antenna digital wireless telecommunication system
CN110505020B (en) Optical nonlinear equalization method based on transfer learning
CN111786724B (en) Multi-wavelength LED underwater visible light communication modulation method based on deep learning
CN101827060B (en) Adaptive modulation-demodulation method base on fractional order Fourier transform
CN101674160A (en) Signal detection method and device for multiple-input-multiple-output wireless communication system
CN109412993A (en) A kind of having based on deep learning inhibits the full duplex cognition water sound communication signal detection method of self-interference
CN109039534A (en) A kind of sparse CDMA signals detection method based on deep neural network
CN112953653B (en) Single-carrier multi-user underwater acoustic communication method
CN110098870B (en) Optical generalized spatial modulation method based on OB-MMSE detection algorithm
CN112215335A (en) System detection method based on deep learning
CN113825159A (en) Wireless energy-carrying communication system robust resource allocation method based on intelligent reflector
CN116318311A (en) Transmission method based on reconfigurable intelligent surface anti-phase index modulation
Hassan et al. Polynomial GSVD beamforming for two-user frequency-selective MIMO channels
CN101340224A (en) Bit distribution method and apparatus
Alaca et al. CNN-based signal detector for IM-OFDMA
CN114124168A (en) MIMO-NOMA system signal detection method and system based on deep learning
CN112543056B (en) Power and grouping combined optimization method for PD-NOMA-VLC system
CN112564830B (en) Deep learning-based dual-mode orthogonal frequency division multiplexing index modulation detection method and device
CN117220740A (en) Beam forming method of communication and interference integrated system in non-cooperative scene
CN104618297B (en) Signal transmitting and method of reseptance based on OFDM-DCSK communication systems
CN107181705A (en) A kind of half-blind channel estimating method and system
Li et al. Model-driven deep learning scheme for adaptive transmission in MIMO-SCFDE system

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