Mode selection and transceiver design method in NOMA-assisted ISAC system

Abstract

The invention discloses a mode selection and transceiver design method of a non-orthogonal multiple access (NOMA) auxiliary sensing and communication integration (ISAC) in a network auxiliary full-duplex cellular-free large-scale multiple-input multiple-output (MIMO) system, and relates to a mode selection and transceiver design method of a non-orthogonal multiple access (NOMA) auxiliary sensing and communication integration (ISAC) system of a network auxiliary full-duplex cellular-free large-scale multiple-input multiple-output (MIMO) system, and the mode selection and transceiver design method of the non-orthogonal multiple access (NOMA) auxiliary sensing and communication integration (ISAC) system of the network auxiliary full-duplex cellular-free MIMO system of the network auxiliary full-duplex cellular-free MIMO system of the network auxiliary full-duplex cellular-free MIMO system of the network auxiliary full-duplex cellular-free MIMO system. Comprising the following steps of: converting a joint mode selection and transceiver design problem of NOMA-ISAC in a cellular-free large-scale MIMO system into an optimization problem of maximizing a signal-to-interference-noise ratio (SINR) of a radar, and decomposing the optimization problem into two sub-problems, and respectively adopting a simulated annealing algorithm and a continuous convex approximation algorithm to obtain the optimal working mode, precoding vector and uplink power. Simulation results show that the sensing performance of the system can be effectively improved under the condition that rate constraints of communication users are met.

Description

Technical Field

[0001] The present invention relates to the field of wireless communication technologies, and in particular to a mode selection and transceiver design method in a NOMA-assisted ISAC system. Background Art

[0002] Integrated sensing and communication (ISAC) is one of the six major application scenarios of the sixth-generation mobile communication system. ISAC enables devices to perceive the surrounding environment while communicating, which plays an important role in promoting emerging industries including the Internet of Things and autonomous driving.

[0003] However, integrating communication and perception functions on the same device will also bring strong interference problems, which poses a challenge to the actual deployment of ISAC. To solve this problem, more advanced interference management technologies, such as non-orthogonal multiple access (NOMA), are used to assist the ISAC system. At the same time, in order to give full play to the potential coordinated perception capabilities of multiple base stations, we integrate NOMA-ISAC into the network-assisted full-duplex non-cellular massive MIMO system for research. The network-assisted full-duplex (NAFD) that has emerged in recent years can complete the transmission and reception of data on the same time-frequency resource block. Compared with the traditional TDD / FDD mode, it can double the spectrum efficiency of the mobile communication link. At the same time, compared with the CCFD mode, it can effectively alleviate the self-interference problem caused by full-duplex, thereby effectively improving the system performance.

[0004] There are still many problems to be solved in implementing NOMA-ISAC in network-assisted full-duplex cell-free massive MIMO systems. Reasonable RAU working mode and transceiver design can effectively improve system performance. Therefore, it is very important to jointly optimize the working mode and transceiver design. Summary of the invention

[0005] The purpose of the present invention is to provide a mode selection and transceiver design method in a NOMA-assisted ISAC system to effectively improve system performance.

[0006] To achieve this object, the present invention adopts the following technical solutions:

[0007] A method for mode selection and transceiver design in a NOMA-assisted ISAC system, comprising the following steps:

[0008] Step 1, establish the optimal receiver expression of the NOMA-assisted ISAC in the network-assisted full-duplex non-cellular massive MIMO system, and then obtain the expression of the system-perceived SINR, and then jointly optimize the duplex mode selection, downlink precoding and uplink power allocation to maximize the system-perceived SINR and model it as a mathematical optimization problem; since the optimization problem is non-convex, the block coordinate descent method is used to decompose the optimization problem into two sub-problems of mode selection and transceiver design, and the optimal solution is obtained based on the simulated annealing algorithm and the continuous convex approximation algorithm respectively;

[0009] Step 2: When the number of iterations of the simulated annealing algorithm reaches the upper limit, the optimal working mode is output; when the number of iterations of the continuous convex approximation algorithm reaches the upper limit, the transceiver design scheme is output.

[0010] Preferably: the optimization problem model for maximizing the system perceived SINR is:

[0011]

[0012] Where w = {w k |0≤k≤K D} and p={p j |1≤j≤K U} respectively represent the downlink precoding w k and the set of uplink power, T m =diag[0 (m-1)N , 1 N , 0 (M-m)N ], M represents the number of remote antenna units, N represents the number of antennas, K U Indicates the number of uplink users, K D represents the number of downlink users, T represents the total number of radars, q D ∈{0, 1} M×1 Select the downlink RAU working mode vector, q D (m) indicates the downlink RAU working mode corresponding to the mth remote antenna unit. If it is 1, it means that the mth remote antenna unit is in the downlink transmission mode, otherwise it is in the uplink reception mode. j The power of uplink user j, γ D，k represents the SINR of downlink user k, which is obtained by decoding the signal s k When γ , downlink user k regards the signals from other downlink users and all uplink users as interference; U，j represents the SINR of uplink user j, which is obtained by CPU analyzing the signal d from uplink user j. jWhen decoding, the signals from the uplink users with higher priority have been decoded and completely eliminated from the received signal, and the signals from the uplink users with lower priority and the signals of all downlink users are treated as interference; γ D，min , γ U，min , P U，max , P D，max They represent the minimum SINR required by downlink users, the minimum SINR required by uplink users, the maximum transmit power of uplink users, and the maximum transmit power of downlink remote antenna units respectively; s，t represents the perceived SINR corresponding to the radar target t.

[0013] Preferably: in step 1, the block coordinate descent method is used to decompose the optimization problem into two sub-problems of mode selection and transceiver design, and the optimal solution is obtained based on the simulated annealing algorithm and the continuous convex approximation algorithm respectively. The specific steps include:

[0014] Step 1.1: Initialize q D , randomly generate a 0-1 vector, the i-th element represents the working mode of the i-th RAU, 1 means it works in downlink transmission mode, and 0 means it works in uplink reception mode;

[0015] Step 1.2: Under a given working mode, the optimal downlink precoding and uplink power are solved based on continuous convex approximation, and the corresponding perceived SINR is obtained;

[0016] Step 1.3: Set the maximum number of iterations L of the simulated annealing algorithm out , the control factor θ is a pre-set parameter used to control the population convergence speed. If the current number of iterations does not reach the set upper limit, the evolution operation is performed; during the evolution process, q D The mth element q in D (m) is changed to 1-q D (m), by solving the optimization problem (1), we can obtain the optimal downlink precoding and uplink power corresponding to the new duplex mode, as well as the corresponding perceived SINR; the perceived SINR before and after the modification is respectively denoted as and

[0017] Step 1.4: If Then change qD(m) to 1-qD(m), if Assume q D (m) The probability of still having P is changed to 1-q D (m), that is, if q D If (m) is 1, change it to 0; otherwise, change it to 1. The expression of P is

[0018]

[0019] When the current number of iterations reaches the set upper limit, go to step 2.

[0020] Preferably, step 1.2 specifically includes:

[0021] Defining the Matrix 0≤k≤K D And ignoring the rank constraint, we get The SINR of downlink user k is expressed as

[0022]

[0023] in, represents the effective channel vector from downlink user k to all RAUs, represents the downlink RAU mode selection matrix, and the uplink RAU mode selection matrix is h D，k represents the channel vector from downlink user k to all RAUs, h IUI，j，k represents the channel between downlink user k and uplink user j, Indicates the noise power at the downlink user;

[0024] The downlink RAU power constraint in the optimization problem shown in equation (1) is rewritten as

[0025]

[0026] On this basis, under a given duplex mode, the sub-problems of downlink precoding and uplink power optimization are remodeled as

[0027]

[0028] Based on the continuous convex approximation algorithm, the non-convex problem shown in the above formula (5) is transformed into the following convex problem for solution

[0029]

[0030] Among them, n s ={n s，t |1≤t≤T} and m s ={m s，t |1≤t≤T} is the temporary variable n introduced s，t and m s，t A collection of Respectively represent m s，t and n s，t The value at the nth iteration, and Respectively represent Ψ t and Φ jThe value at the nth iteration, represents the echo signal matrix reflected by radar target t, h e，D，t represents the channel from radar target t to all downlink RAUs, Uplink RAU mode selection matrix for providing radar target perception services, represents the effective channel vector from uplink user j to all RAU antennas, h U，j represents the channel vector from uplink user j to all RAU antennas, represents the uplink noise power, η is the radar cross section, E represents the disturbance matrix; the function g(.), f t (.) and k j (.) are defined as

[0031]

[0032]

[0033] in, h e，U，t represents the channel from radar target t to all uplink RAUs; the optimal downlink precoding and uplink power in a given working mode are obtained by iteratively calculating formula (6) multiple times, and the corresponding perceived SINR is obtained.

[0034] Preferably, in step 2, check whether the current iteration number i has reached the maximum number of iterations L out , if i<L out , then return to step 1, otherwise, output the current optimal duplex mode selection and transceiver design solution.

[0035] Beneficial effects: The method of the present invention converts the RAU duplex mode and transceiver design of NOMA-ISAC in the network-assisted full-duplex non-cellular massive MIMO system into vectors to be optimized, converts the joint optimization problem of mode selection and transceiver design into an optimization problem of maximizing the system perceived SINR, and proposes an algorithm based on simulated annealing and continuous convex approximation to solve the problem. The mode selection method based on simulated annealing and the transceiver design method based on continuous convex approximation provided by the present invention can effectively improve system performance compared with the MRT precoding and random mode selection method. BRIEF DESCRIPTION OF THE DRAWINGS

[0036] Figure 1 It is an algorithm flow chart of a specific implementation mode of the present invention;

[0037] Figure 2 It is a schematic diagram of a network-assisted full-duplex cell-free massive MIMO system according to a specific embodiment of the present invention;

[0038] Figure 3 It is a performance comparison diagram of the method according to the specific implementation mode of the present invention and the MRT precoding and random mode selection method. DETAILED DESCRIPTION

[0039] The present invention will be further described below in conjunction with the accompanying drawings.

[0040] First, let’s make the following explanation: Consider a typical single-cell network-assisted full-duplex non-cellular massive MIMO system with M remote antenna units (RAUs). All RAUs are distributed in the cell according to Poisson distribution and connected to the central processing unit (CPU) through a low-latency backhaul link. Each RAU is equipped with N antennas. There are K U uplink users and K D The total number of downlink users and radars is denoted as K and T respectively. The channel model and signal transmission model of the system are briefly introduced as follows.

[0041] The channel vector from downlink user k to all RAU antennas is defined as but

[0042]

[0043] In the formula, represents the large-scale fading matrix, λ m，k is the large-scale fading factor from RAU m to downlink user k, represents the small-scale fading vector.

[0044] Similarly, the channel vector from uplink user j to all antennas on the RAU is defined as The channel from uplink user j to downlink user k is The self-interference matrix between RAUs is: The echo signal matrix reflected by radar target t is H e，t It can be expressed as

[0045]

[0046] In the network-assisted full-duplex cell-free massive MIMO system, for the convenience of expression, the binary vector q is defined U ∈{0, 1} M×1 and q D ∈{0, 1} M×1 As the uplink and downlink RAU working mode selection vectors, respectively. Specifically, if RAU i works in the uplink receiving mode, then q U (i) = 1 and q D(i) = 0; otherwise, if RAU i works in downlink transmission mode, then q D (i) = 1 and q U (i)=0, it is easy to get q U +q D =1 M×1 In addition, the RAU mode selection matrix of the system is further defined as and

[0047] For downlink transmission, the signals sent by M RAUs It can be expressed as

[0048]

[0049] in, and They represent the symbols sent to downlink user i and for target sensing, respectively. and represent the corresponding precoding vectors respectively. Therefore, the signal received by downlink user k is It can be expressed as

[0050]

[0051] in, represents the signal sent by uplink user j, then the power of uplink user j can be expressed as p j =E{|d j | 2}. Additionally, is additive Gaussian white noise, satisfying

[0052] Similarly, for uplink transmission, the signal received by all RAUs It can be expressed as

[0053]

[0054] in, is an independent and identically distributed additive Gaussian white noise, satisfying

[0055] Since supporting NOMA requires relatively complex hardware, it is often difficult to implement for single-antenna users. Therefore, it is assumed that the NOMA strategy is only used in uplink communications to assist the ISAC system. Without loss of generality, it is assumed that uplink users are arranged in increasing order of channel quality, that is, For downlink user k, after decoding s k When , the signals from other downlink users and all uplink users are regarded as interference, and the corresponding communication SINR can be obtained as

[0056]

[0057] in, represents the effective channel vector from downlink user k to all RAUs. When the CPU receives the signal d from uplink user j j When decoding is performed, the signals from the higher priority uplink users have already been decoded and are completely eliminated from the received signal. Define the receive vector To decode d j , then the corresponding SINR is

[0058]

[0059] in, represents the effective channel vector from uplink user j to all RAU antennas, represents the effective interference channel caused by the radar target echo and the downlink transmission. In addition, the covariance matrix of the downlink ISAC signal can be written as After decoding and eliminating all signals from uplink users, define To capture the reflected signal of radar target t. In addition, it is assumed that only the L uplink RAUs closest to the sensing area provide sensing services for the radar target. Definition and As the corresponding binary uplink RAU selection vector and matrix, it is necessary to satisfy On this basis, the perceived SINR corresponding to the radar target t is:

[0060]

[0061] in, and The corresponding optimal filters are

[0062]

[0063]

[0064] in, Similarly, we can get E represents the disturbance matrix. Substituting the above optimal filters into the corresponding SINR expressions, we can get

[0065]

[0066]

[0067] Duplex mode selection, downlink precoding and uplink power allocation have a direct impact on the performance of NOMA-assisted ISAC in network-assisted full-duplex non-cellular massive MIMO system. Reasonable duplex mode selection, downlink precoding and uplink power allocation are crucial to improving the performance of the system. To this end, the present invention maximizes the perceived SINR of the system under the constraint of satisfying the communication user rate by optimizing the duplex mode selection, downlink precoding and uplink power allocation. The optimization problem can be abstracted as the following mathematical problem

[0068]

[0069] Where w = {w k |0≤k≤K D} and p={p j |1≤j≤K U} respectively represent the set of downlink precoding and uplink power. At the same time, T m =diag[0 (m-1 ) N , 1 N , 0 (M-m)N ].

[0070] Based on the analysis of communication and radar perception SINR in the above system, the working mode selection, downlink precoding and uplink power allocation method based on simulated annealing algorithm and continuous convex approximation algorithm provided by the present invention can be divided into the following steps:

[0071] Step 1: Initialize q D , randomly generate a 0-1 vector, the i-th element represents the working mode of the i-th RAU, 1 means it works in the downlink transmission mode, and 0 means it works in the uplink reception mode.

[0072] Step 2: Under a given working mode, the optimal downlink precoding and uplink power are solved based on continuous convex approximation. Inspired by the semi-positive definite relaxation algorithm, the matrix is ​​defined as 0≤k≤K D And ignore the rank constraint. Then we can get The SINR of downlink user k can be expressed as

[0073]

[0074] in, In addition, the downlink RAU power constraint in the optimization problem shown in equation (13) can be rewritten as

[0075]

[0076] On this basis, under a given duplex mode, the sub-problems of downlink precoding and uplink power optimization can be remodeled as

[0077]

[0078] The optimization problem shown in formula (16) is difficult to solve because the objective function and some constraints are non-convex. Based on the continuous convex approximation algorithm, the above non-convex problem is transformed into the following convex problem for solution

[0079]

[0080] Among them, n s ={n s，t |1≤t≤T} and m s ={m s，t |1≤t≤T} is the set of temporary variables introduced. Functions g(.), f t (.) and k j (.) are defined as

[0081]

[0082]

[0083]

[0084] in, and Respectively represent Ψ t and Φ j The value at the nth iteration. The optimal downlink precoding and uplink power in a given working mode are obtained by calculating formula (17) multiple times through iterations, and the corresponding perceived SINR is obtained.

[0085] Step 3: Set the maximum iteration number L of the simulated annealing algorithm out , control factor θ, if the current number of iterations does not reach the upper limit, the evolution operation is performed. During the evolution process, according to step 2, q D The mth element q in D (m) is changed to 1-q D The optimal downlink precoding and uplink power corresponding to the duplex mode of (m), and the corresponding perceived SINR. The perceived SINR before and after the modification is respectively denoted as and

[0086] Step 4: If Then change qD(m) to 1-qD(m), if Different from the greedy strategy, assuming q D (m) The probability of still having P is changed to 1-q D (m), the expression of P is

[0087]

[0088] In this way, we can avoid falling into the local optimal solution. When the current number of iterations reaches the set upper limit, go to step 5;

[0089] Step 5: Output the current optimal duplex mode selection and transceiver design solution.

[0090] The above is only a preferred embodiment of the present invention. It should be pointed out that for ordinary technicians in this technical field, several improvements and modifications can be made without departing from the principle of the present invention. These improvements and modifications should also be regarded as the scope of protection of the present invention.

Claims

1. A mode selection and transceiver design method in a NOMA-assisted ISAC system, characterized by: The following steps are involved: Step 1, establish the optimal receiver expression of the NOMA-assisted ISAC in the network-assisted full-duplex non-cellular massive MIMO system, and then obtain the expression of the system-perceived SINR, and then jointly optimize the duplex mode selection, downlink precoding and uplink power allocation to maximize the system-perceived SINR and model it as a mathematical optimization problem; since the optimization problem is non-convex, the block coordinate descent method is used to decompose the optimization problem into two sub-problems of mode selection and transceiver design, and the optimal solution is obtained based on the simulated annealing algorithm and the continuous convex approximation algorithm respectively; Step 2: When the number of iterations of the simulated annealing algorithm reaches the upper limit, the optimal working mode is output; When the number of iterations of the continuous convex approximation algorithm reaches an upper limit, the transceiver design solution is output.

2. The method for mode selection and transceiver design in a NOMA-assisted ISAC system according to claim 1, characterized in that: The optimization problem model for maximizing the system perceived SINR is: Where w = {w k |0≤k≤K D } and p={p j |1≤j≤K U } respectively represent the downlink precoding w k and the set of uplink power, T m =diag[0 (m-1)N , 1 N , 0 (M-m)N ], M represents the number of remote antenna units, N represents the number of antennas, J U Indicates the number of uplink users, K D represents the number of downlink users, T represents the total number of radars, q D ∈{0, 1} M×1 Select the downlink RAU working mode vector, q D (m) indicates the downlink RAU working mode corresponding to the mth remote antenna unit. If it is 1, it means that the mth remote antenna unit is in the downlink transmission mode, otherwise it is in the uplink reception mode. j The power of uplink user j, γ D,k represents the SINR of downlink user k, which is obtained by decoding the signal s k When γ , downlink user k regards the signals from other downlink users and all uplink users as interference; U,j represents the SINR of uplink user j, which is obtained by CPU analyzing the signal d from uplink user j. j When decoding, the signals from the uplink users with higher priority have been decoded and completely eliminated from the received signal, and the signals from the uplink users with lower priority and the signals of all downlink users are treated as interference; γ D,min , γ U,min , P U,max , P D,max They represent the minimum SINR required by downlink users, the minimum SINR required by uplink users, the maximum transmit power of uplink users, and the maximum transmit power of downlink remote antenna units respectively; s,t represents the perceived SINR corresponding to the radar target t.

3. The method for mode selection and transceiver design in a NOMA-assisted ISAC system according to claim 2, characterized in that: In step 1, the block coordinate descent method is used to decompose the optimization problem into two sub-problems: mode selection and transceiver design, and the optimal solution is obtained based on the simulated annealing algorithm and the continuous convex approximation algorithm respectively. The specific steps include: Step 1.1: Initialize q D , randomly generate a 0-1 vector, the i-th element represents the working mode of the i-th RAU, 1 means it works in downlink transmission mode, and 0 means it works in uplink reception mode; Step 1.2: Under a given working mode, the optimal downlink precoding and uplink power are solved based on continuous convex approximation, and the corresponding perceived SINR is obtained; Step 1.3: Set the maximum number of iterations L of the simulated annealing algorithm out , the control factor θ is a pre-set parameter used to control the population convergence speed. If the current number of iterations does not reach the set upper limit, the evolution operation is performed; during the evolution process, q D The mth element q in D (m) is changed to 1-q D (m), by solving the optimization problem (1), we can obtain the optimal downlink precoding and uplink power corresponding to the new duplex mode, as well as the corresponding perceived SINR; the perceived SINR before and after the modification is respectively denoted as and Step 1.4: If Then q D (m) is changed to 1-q D (m) if Assume q D (m) The probability of still having P is changed to 1-q D (m), that is, if q D If (m) is 1, change it to 0; otherwise, change it to 1. The expression of P is When the current number of iterations reaches the set upper limit, go to step 2.

4. The method for mode selection and transceiver design in a NOMA-assisted ISAC system according to claim 3, characterized in that: Step 1.2 specifically includes: Defining the Matrix And ignoring the rank constraint, we get The SINR of downlink user k is expressed as in, represents the effective channel vector from downlink user k to all RAUs, represents the downlink RAU mode selection matrix, and the uplink RAU mode selection matrix is h D，k represents the channel vector from downlink user k to all RAUs, h IUI，j，k represents the channel between downlink user k and uplink user j, Indicates the noise power at the downlink user; The downlink RAU power constraint in the optimization problem shown in equation (1) is rewritten as On this basis, under a given duplex mode, the sub-problems of downlink precoding and uplink power optimization are remodeled as Based on the continuous convex approximation algorithm, the non-convex problem shown in the above formula (5) is transformed into the following convex problem for solution Among them, n s ={n s，t |1≤t≤T} and m s ={m s，t |1≤t≤T} is the temporary variable n introduced s,t and m s,t A collection of Respectively represent m s,t and n s,t The value at the nth iteration, and Respectively represent Ψ t and Φ j The value at the nth iteration, represents the echo signal matrix reflected by radar target t, h e，D，t represents the channel from radar target t to all downlink RAUs, Uplink RAU mode selection matrix for providing radar target perception services, represents the effective channel vector from uplink user j to all RAU antennas, h U，j represents the channel vector from uplink user j to all RAU antennas, represents the uplink noise power, η is the radar cross section, E represents the disturbance matrix; the function g(.), f t (.) and k j (.) are defined as in, h e，U，t represents the channel from radar target t to all uplink RAUs; the optimal downlink precoding and uplink power in a given working mode are obtained by iteratively calculating formula (6) multiple times, and the corresponding perceived SINR is obtained.

5. The method for mode selection and transceiver design in a NOMA-assisted ISAC system according to claim 1, characterized in that: In step 2, check whether the current iteration number i has reached the maximum number of iterations L out , if i<L out , then return to step 1, otherwise, output the current optimal duplex mode selection and transceiver design solution.

Images