NEAR-FIELD BROADBAND UPLINK MIMO TRANSMISSION METHOD ASSISTED BY DYNAMIC METASURFACE ANTENNA

Abstract

The present invention discloses a near-field broadband uplink MIMO transmission method assisted by a dynamic metasurface antenna. The method includes: Broadband signals sent by a plurality of users distributed in a near-field region are processed with a large-size dynamic metasurface antenna as a receive antenna on a base station side, which can reduce system hardware costs and power consumption; and compared with the current hybrid beamforming based on a phase shifter and a conventional antenna, hybrid beamforming based on the dynamic metasurface antenna can effectively improve transmission performance. The present invention proposes an algorithm framework jointly designing a dynamic metasurface antenna and a baseband beamformer and including method such as matrix-weighted mean square error sum (MWMSE) minimization, alternate optimization, matrix vectorization, and MM. The present invention implements near-field broadband large-scale MIMO uplink transmission assisted by a dynamic metasurface antenna with low algorithm complexity and good convergence.

Claims

1. A near-field broadband uplink MIMO transmission method assisted by a dynamic metasurface antenna, wherein in the method, a sum rate maximization problem is constructed based on a broadband large-scale MIMO uplink single-cell system and a channel model that considers a near-field effect, frequency-selective fading, and a spatial broadband effect, wherein a dynamic metasurface antenna array is used on a base station side, and the sum rate maximization problem is solved in a manner of jointly designing a baseband beamforming matrix and a weight matrix of a dynamic metasurface antenna, to maximize a near-field broadband large-scale MIMO uplink sum rate; and the transmission method comprises: step S1, giving a weight matrix of a dynamic metasurface antenna, and solving, based on MWMSE transformation, a baseband beamforming matrix on each subcarrier according to a system sum rate maximization criterion; step S2, giving a baseband beamforming matrix, and solving, based on matrix vectorization, an MM method, and a convex optimization method, a weight matrix of a dynamic metasurface antenna according to the system sum rate maximization criterion; and step S3, cyclically performing step S1 and step S2 until a difference between two adjacent system sum rates is less than a given threshold, wherein in a moving process of users, as state information of a channel from each user to the base station side is changed, the near-field broadband uplink MIMO transmission method assisted by a dynamic metasurface antenna is dynamically implemented.

2. The near-field broadband uplink MIMO transmission method assisted by a dynamic metasurface antenna according to claim 1, wherein a specific expression of the channel model that considers the near-field effect, the frequency-selective fading, and the spatial broadband effect is: $\begin{array}{cc}{g}_u\left(f\right)=\underset{p=0}{\overset{P_u}{.Math.}}{a}_{u,p}\left(f\right)?b_{u,p}\left(f\right)& \left(1\right)\end{array}$ in the expression (1), $\begin{array}{cc}{a}_{u,p}\left(f\right)=\left[\begin{array}{c}\sqrt{{?}_{1,1,u,p}}{A}_{1,1,u,p}\left(f\right),\\ \sqrt{{?}_{1,2,u,p}}{A}_{1,2,u,p}\left(f\right),.Math.,\sqrt{{?}_{M,L,u,p}}{A}_{M,L,u,p}\left(f\right)\end{array}\right]& \left(2\right)\end{array}$ $\begin{array}{cc}{b}_{u,p}\left(f\right)=\left[\begin{array}{c}{e}^{-\mathrm{?2?}\left(f+f_c\right)\frac{.Math.p_{u,p}-p_{1,1}.Math.}{2}},\\ e^{-\mathrm{?2?}\left(f+f_c\right)\frac{.Math.p_{u,p}-p_{1,2}.Math.}{2}},.Math.,e^{-\mathrm{?2?}\left(f+f_c\right)\frac{.Math.p_{u,p}-p_{M,L}.Math.}{2}}\end{array}\right]& \left(3\right)\end{array}$ wherein a.sub.u,p(f) and b.sub.u,p(f) represent a channel gain that considers a near-field effect, frequency selectivity, and a spatial broadband effect and a response matrix of an antenna array respectively; ?.sub.m,l,u,p and A.sub.m,l,u,p(f) represent a large-scale fading factor of a (p).sup.th transmission path between a (l).sup.th metamaterial on a (m).sup.th microstrip of a base station antenna and a user u and a channel gain coefficient respectively, p.sub.u,p and p.sub.m,l represent a scatterer position of the (p).sup.th transmission path between the user u and the base station and a position of the (l).sup.th metamaterial on the (m).sup.th microstrip of the base station antenna respectively, f and f.sub.c represent a frequency and a center frequency respectively, and c represents a signal transmission speed equal to 3?10.sup.8; the channel gain coefficient is expressed as: $\begin{array}{cc}{A}_{m,l,u,p}\left(f\right)=.Math.{?}_{u,p}\left(f\right).Math.\sqrt{F\left({?}_{m,l,u,p}\right)}\frac{c}{4?\left(f+f_c\right){.Math.p_{u,p}-p_{m.l}.Math.}_F}& \left(4\right)\end{array}$ in the expression (4), ?.sub.m,l,u,p=(?.sub.m,l,u,p?.sub.m,l,u,p) represents height-azimuth of a signal reflected from the user u by a (p).sup.th reflector and reaching a (l).sup.th antenna unit on the (m).sup.th microstrip of the base station antenna, and F(?.sub.m,l,u,p) is expressed as: $\begin{array}{cc}F\left({?}_{m,l,u,p}\right)=\{\begin{array}{cc}6{\cos }^2\left({?}_{m,l,u,p}\right),& {?}_{m,l,u,p}?\left[0,\frac{?}{2}\right]\\ 0,& \mathrm{otherwise}\end{array}& \left(5\right)\end{array}$ ?.sub.u,p(f) refers to a reflection coefficient of a reflector on a (p).sup.th path of the user u, and is expressed as: $\begin{array}{cc}{?}_{u,p}\left(f\right)=\{\begin{array}{cc}\frac{\cos {?}_{1,u,p}-n_t\cos {?}_{t,u,p}}{\cos {?}_{1,u,p}+n_t\cos {?}_{t,u,p}}{e}^{-\left(\frac{8{{?}^2\left(f+f_c\right)}^2{?}_{\mathrm{rough}}^2{\cos }^2{?}_{1,u,p}}{c^2}\right)},& p=1,2,.Math.,P\\ 1,& p=0\end{array}& \left(6\right)\end{array}$ in the expression (6), n.sub.t is a refractive index, ?.sub.rough.sup.2 is a roughness coefficient of a reflection surface, and cos ?.sub.i,u,p and cos ?.sub.t,u,p are an incident angle and a reflection angle of the signal of the user u on a (p).sup.th reflection object respectively.

3. The near-field broadband uplink MIMO transmission method assisted by a dynamic metasurface antenna according to claim 2, wherein the sum rate maximization problem is defined as a first optimization problem, and a specific expression of the problem is: $\begin{array}{cc}{?}_1:\underset{Q,W_s,?s?\left\{1,2,.Math.,S\right\}}{\max }\underset{s=0}{\overset{S-1}{.Math.}}{?}_B\log .Math.I_U+\frac{P_t}{{?}_B{?}^2}{W}_s^H{\mathrm{QH}}_s{G}_s{G}_s^H{H}_s^H{Q}^H{W_s\left(W_s^H{\mathrm{QH}}_s{H}_s^H{Q}^H{W}_s\right)}^{-1}.Math.& \left(7\right)\end{array}$ $s.t.{\left(Q\right)}_{m_1\left(m_2-1\right)L+l}=\{\begin{array}{cc}{q}_{m_1,l},& m_1=m_2\\ 0,& m_1?m_2\end{array},$ $q_{m_1,l}?,?m_1,l$ in the expression (7), a target function is a sum rate of a broadband large-scale MIMO uplink system based on a dynamic metasurface antenna, S represents a quantity of subcarriers, and ?.sub.B represents a subcarrier spacing and is expressed as a ratio of a bandwidth B to the quantity of subcarriers S; I.sub.U is an identity matrix of U?U, ?.sup.2 is a variance of noise, P.sub.t represents a transmit power, and U is a quantity of users in a cell; G.sub.s=[g.sub.1,s,g.sub.2,s, . . . , g.sub.U,s]?.sup.N.sup.R.sup.?U represents a channel matrix of an (s).sup.th subcarrier, W.sub.s?.sup.M?U represents a baseband beamformer of the (s).sup.th subcarrier, H.sub.s?.sup.N.sup.R.sup.?N.sup.R describes a frequency-selective effect of a signal propagated on a microwave transmission band, and Q?.sup.M?N.sup.R represents a weight matrix of a dynamic metasurface antenna; and log is a logarithm operation, and |.Math.| is a matrix determinant obtaining operation.

4. The near-field broadband uplink MIMO transmission method assisted by a dynamic metasurface antenna according to claim 3, wherein the step S1 specifically comprises: step S101, obtaining an equivalent mean square error sum minimization problem of sum rate maximization with MWMSE transformation: $\begin{array}{cc}\underset{M_s,W_s,?s?\left\{1,2,...,S\right\}}{\min }\underset{s=0}{\overset{S-1}{.Math.}}\mathrm{tr}\left\{M_s{E}_s\left(Q,W_s\right)\right\}-{\log }_2.Math.M_s.Math.& \left(8\right)\end{array}$ in the expression (8), E.sub.s(Q, W.sub.s) is a mean square error sum matrix, whose specific expression is: $\begin{array}{cc}{E}_s\left(Q,W_s\right)=P_t\left(W_s^H{\mathrm{QH}}_s{G}_s-I_U\right){\left(W_s^H{\mathrm{QH}}_s{G}_s-I_U\right)}^H+{?}_B{?}^2{W}_s^H{\mathrm{QH}}_s{H}_s^H{Q}^H{W}_s& \left(9\right)\end{array}$ M.sub.s is a weighted auxiliary matrix and is obtained from an expression (10), whose expression is: $\begin{array}{cc}{M}_s^{\mathrm{opt}}={E_s\left(Q,W_s\right)}^{-1}& \left(10\right)\end{array}$ step S102, optimizing a baseband beamforming matrix on each subcarrier when a weight matrix of a dynamic metasurface antenna is given, to obtain a system sum rate optimization problem based on MWMSE transformation, wherein the problem is defined as a second optimization problem, and a specific expression of the second optimization problem is: $\begin{array}{cc}{?}_2:\underset{W_s,?s?\left\{1,2,.Math.,S\right\}}{\min }\underset{s=0}{\overset{S-1}{.Math.}}\mathrm{tr}\left\{M_s{E}_s\left(Q,W_s\right)\right\}-{\log }_2.Math.M_s.Math.& \left(11\right)\end{array}$ when E.sub.s(Q, W.sub.s) and M.sub.s are given, W.sub.s is given by the following expression, wherein ?s?{1, 2, . . . , S}: $\begin{array}{cc}{W}_s^{\mathrm{opt}}={\left(P_t{\mathrm{QH}}_s{G}_s{G}_s^H{H}_s^H{Q}^H+{?}_B{?}^2{\mathrm{QH}}_s{H}_s^H{Q}^H\right)}^{-1}{\mathrm{QH}}_s{G}_s& \left(12\right)\end{array}$ step S103, iteratively updating the mean square error sum matrix E.sub.s(Q, W.sub.s), the weighted auxiliary matrix M.sub.s, and the baseband beamforming matrix W.sub.s, wherein ?s?{1, 2, . . . , S}, to obtain a baseband beamforming matrix meeting that a difference between adjacent two iteration results is less than a given threshold as a solution of a baseband beamforming matrix whose system sum rate is maximized when a weight matrix of a dynamic metasurface antenna is given.

5. The near-field broadband uplink MIMO transmission method assisted by a dynamic metasurface antenna according to claim 4, wherein the step S2 specifically comprises: step S201, obtaining, when a baseband beamforming matrix is given, a problem of optimizing a weight matrix of a dynamic metasurface antenna to maximize a system sum rate, and obtaining an equivalent mean square error sum minimization problem with MWMSE transformation, wherein the problem is defined as a third optimization problem, whose specific expression is: $\begin{array}{cc}{?}_3:\underset{Q}{\min }\underset{s=0}{\overset{S-1}{.Math.}}\mathrm{tr}\left\{P_t{M}_s{W}_s^H{\mathrm{QH}}_s{G}_s{G}_s^H{H}_s^H{Q}^H{W}_s\right\}-\mathrm{tr}\left\{P_t{M}_s{W}_s^H{\mathrm{QH}}_s{G}_s\right\}-\mathrm{tr}\left\{P_t{M}_s{G}_s^H{H}_s^H{Q}^H{W}_s\right\}+\mathrm{tr}\left\{{?}_B{?}^2{M}_s{W}_s^H{\mathrm{QH}}_s{H}_s^H{Q}^H{W}_s\right\}& \left(13a\right)\end{array}$ $\begin{array}{cc}s.t.{\left(Q\right)}_{m_1,\left(m_2-1\right)L+l}=\{\begin{array}{cc}{q}_{m_1,l},& m_1=m_2\\ 0,& m_1?m_2\end{array},& \left(13b\right)\end{array}$ $\begin{array}{cc}{q}_{m_1,l}??,?m_1,l& \left(13c\right)\end{array}$ step S202, removing a physical structure constraint of a dynamic metasurface antenna, that is, the expression (13b) with a matrix vectorization method, to obtain a problem of maximizing an equivalent expression of a system sum rate having only a weight feasible domain constraint of the dynamic metasurface antenna, wherein the problem is defined as a fourth optimization problem; step S203, designing a weight of the dynamic metasurface antenna with an MM algorithm and a convex optimization algorithm in consideration of four weight feasible domain constraints; and step S204, iteratively updating the mean square error sum matrix E.sub.s(Q, W.sub.s), the weighted auxiliary matrix M.sub.s, and the weight matrix Q of the dynamic metasurface antenna, to obtain a weight matrix meeting that a difference between adjacent two iteration results is less than a given threshold as a solution of the weight matrix of the dynamic metasurface antenna whose system sum rate is maximized when a baseband beamforming matrix is given.

6. The near-field broadband uplink MIMO transmission method assisted by a dynamic metasurface antenna according to claim 5, wherein the step S202 specifically comprises: step S2021, pulling the matrix Q into q=[q.sub.1,1,q.sub.1,2, . . . ,q.sub.m,(m?1)L+l, . . . q.sub.M,ML].sup.T, wherein q.sub.m,(m?1)L+l represents an element of a (m).sup.th row and a (l).sup.th column of the matrix Q; step S2022, obtaining the following with a matrix vectorization rule: $\begin{array}{cc}\mathrm{tr}\left\{{\mathrm{QC}}_s\right\}=q^T{c}_s,\mathrm{and}\mathrm{tr}\left\{Q^H{C}_s^H\right\}=c_s^{H}q& \left(14a\right)\end{array}$ $\begin{array}{cc}\mathrm{tr}\left\{Q^H{B}_s{\mathrm{QA}}_s\right\}=q^H\left(B_s.Math.N_L\right)?A_s^{T}q& \left(14b\right)\end{array}$ $\begin{array}{cc}\mathrm{tr}\left\{Q^H{B}_s{\mathrm{QD}}_s\right\}=q^H\left({\tilde{B}}_s.Math.I_L\right)?D_{s}q& \left(14c\right)\end{array}$ in the expressions (14a) to (14c), A.sub.s=H.sub.sG.sub.sG.sub.s.sup.HH.sub.s.sup.H, B.sub.s=W.sub.sM.sub.sW.sub.s.sup.H, C.sub.s=H.sub.sG.sub.sM.sub.sW.sub.s.sup.H, and D.sub.s=H.sub.sH.sub.s.sup.H; L represents a quantity of metasurface units on each microwave transmission band, N.sub.L is an all 1's matrix of L?L, I.sub.L is an identity matrix of L?L, {tilde over (B)}.sub.s is a diagonal matrix, and ({tilde over (B)}.sub.s).sub.m,m=(B.sub.s).sub.m,m; and c.sub.s=[(C.sub.s).sub.1,1,(C.sub.s).sub.2,1, . . . ,(C.sub.s).sub.(m?1)L+l,m, . . . , (C.sub.s).sub.ML,M].sup.T; and step S2023, performing matrix vectorization transformation, wherein the fourth optimization problem is specifically expressed as: $\begin{array}{cc}{?}_4:\underset{q}{\min }{q}^H\mathrm{Sq}-2\mathrm{Re}\left\{q^H{c}^{*}\right\}& \left(15a\right)\end{array}$ $\begin{array}{cc}s.t.q_{m_1,l}??,?m_1,l& \left(15b\right)\end{array}$ in the expression (15a), $S=\underset{s=0}{\overset{S-1}{.Math.}}{P}_t\left(B_s.Math.N_L\right)?A_s^T+{?}_B{?}^2\left({\tilde{B}}_s.Math.I_L\right)?D_s,\mathrm{and}c=\underset{s=0}{\overset{S-1}{.Math.}}{P}_t{c}_s.$

7. The near-field broadband uplink MIMO transmission method assisted by a dynamic metasurface antenna according to claim 6, wherein the step S203 specifically comprises: step S2031, considering four weight feasible domain constraints, comprising: an unconstrained weight, an amplitude weight, a binary amplitude weight, and a Lorentzian constraint phase weight; step S2032, for the unconstrained weight and the amplitude weight, expressing problems of optimizing a weight vector of the dynamic metasurface antenna to maximize the system sum rate as and respectively, wherein the problem and the problem are specifically expressed as: $\begin{array}{cc}:\underset{q}{\min }{q}^H\mathrm{Sq}-2\mathrm{Re}\left\{c^{T}q\right\}& \left(16\right)\end{array}$ $\begin{array}{cc}:\underset{q}{\min }{q}^T\mathrm{Sq}-2\mathrm{Re}\left\{c^T\right\}q& \left(17a\right)\end{array}$ $\begin{array}{cc}s.t.q_{m_1,l}?\left[a,b\right],b>a>0,?m_1,l& \left(17b\right)\end{array}$ wherein the problem and the problem are convex problems and solved through the convex optimization algorithm; step S2033, for the binary amplitude weight, expressing a problem of optimizing a weight vector of the dynamic metasurface antenna to maximize the system sum rate as , wherein the problem is specifically expressed as: $\begin{array}{cc}{?}_7:\underset{q}{\min }{q}^T\mathrm{Sq}-2\mathrm{Re}\left\{c^T\right\}q& \left(18a\right)\end{array}$ $\begin{array}{cc}s.t.q_{m_1,l}?c.Math.\left\{0,1\right\},c>0,?m_1,l& \left(18b\right)\end{array}$ wherein the problem is solved through a brute-force search method; and step S2034, for the Lorentzian constraint phase weight, expressing a problem of optimizing a weight vector of the dynamic metasurface antenna to maximize the system sum rate as , wherein the problem is specifically expressed as: $\begin{array}{cc}{?}_8:\underset{q}{\min }{q}^H\mathrm{Sq}-2\mathrm{Re}\left\{c^{T}q\right\}& \left(19a\right)\end{array}$ $\begin{array}{cc}s.t.q_{m_1,l}?\left\{\frac{?+e^{\mathrm{??}}}{2}|??\left[0,2?\right]\right\},?m_1,l& \left(19b\right)\end{array}$ wherein represents an imaginary unit; and solving the problem through the MM algorithm.

8. The near-field broadband uplink MIMO transmission method assisted by a dynamic metasurface antenna according to claim 7, wherein the solving the problem through the MM algorithm specifically comprises the following steps: step S20341, expressing the weight vector of the dynamic metasurface antenna as $q=\frac{1}{2}\left({\mathrm{?1}}_{N_R}+p\right),$ wherein 1.sub.N.sub.R is an all 1's vector, and $p={\left[e^{{\mathrm{??}}_1},e^{{\mathrm{??}}_2},.Math.,e^{{\mathrm{??}}_{N_R}}\right]}^T;$ and step S20342, finding a tractable asymptotic function, whose expression is: $\begin{array}{cc}{p}^H\mathrm{Sp}?p^H\mathrm{Tp}-2\mathrm{Re}\left\{p^H\left(T-S\right)p^{\left(?\right)}\right\}+{\left(p^{\left(?\right)}\right)}^H\left(T-S\right)p^{\left(?\right)}& \left(20\right)\end{array}$ T=?.sub.maxI, wherein ?.sub.max is a maximum eigenvalue of S; and therefore, the problem is transformed into: $\begin{array}{cc}{?}_9^{\left(?\right)}\underset{p}{\max }\mathrm{Re}\left\{p^H{a}^{\left(?\right)}\right\}& \left(21a\right)\end{array}$ $\begin{array}{cc}s.t.p_n?\left\{e^{\mathrm{??}}|??\left[0,2?\right]\right\},?n& \left(21b\right)\end{array}$ wherein =(?.sub.maxI.sub.N.sub.R?S)+2c*?; and solving the problem by alternately optimizing the vectors a and p.

9. The near-field broadband uplink MIMO transmission method assisted by a dynamic metasurface antenna according to claim 8, wherein the solving the problem by alternately optimizing the vectors a and p specifically comprises: first, initializing and , and setting an iteration index .sub.1=0 and a threshold ?.sub.1; then, giving , and calculating =and ?n; then, giving , and calculating =(?.sub.maxI.sub.N.sub.R?S)+2c*?; and finally, calculating the system sum rate , and if a difference between the (.sub.1).sup.th system sum rate and the (.sub.1+1).sup.th system sum rate is less than the given threshold ?.sub.1, jumping out of the loop, using as a solution meeting the Lorentzian constraint phase weight under the system sum rate maximization criterion when the baseband beamforming matrix is given: otherwise .sub.1=.sub.1+1, and performing the previous three steps again.

10. The near-field broadband uplink MIMO transmission method assisted by a dynamic metasurface antenna according to claim 9, wherein the step S3 specifically comprises: step S301, initializing a baseband beamforming matrix W.sub.s(0), wherein ?s?{1, 2, . . . , S}, a weight matrix Q.sup.(0) of the dynamic metasurface antenna, a weighted auxiliary matrix M.sup.(0), and a system sum rate R.sub.S.sup.(0), wherein a quantity of iteration times is .sub.2=0, and a threshold is ?.sub.2; step S302, giving a weight matrix of the dynamic metasurface antenna, and solving a baseband beamforming matrix according to the expression (12), wherein ?s?{1, 2, . . . , S}; step S303, giving and the baseband beamforming matrix , wherein ?s?{1, 2, . . . , S}, and calculating a mean square error sum matrix E.sub.s(Q, W.sub.s according to the expression (9), wherein s?{1, 2, . . . , S}; step S304, giving the mean square error sum matrix E.sub.s(Q, W.sub.s, wherein ?s?{1, 2, . . . , S}, and calculating a weighted auxiliary matrix according to the expression (10), wherein ?s?{1, 2, . . . , S}; step S305, giving and , wherein ?s?{1, 2, . . . , S}, and solving weight matrices of four dynamic metasurface antennas according to the expressions (16) to (21) respectively; and step S306, calculating a system sum rate , and if |?|??.sub.2 holds, jumping out of the loop, and using (,), wherein ?s?{1,2, . . . , S}, as a solution meeting the baseband beamforming matrix under the system sum rate maximization criterion and the weight matrix of the dynamic metasurface antenna; otherwise .sub.2=.sub.2+1, and performing step S302 to step S306 again.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

[0054] FIG. 1 is a schematic diagram of near-field broadband large-scale MIMO uplink communication assisted by a dynamic metasurface antenna according to Embodiment 1.

[0055] FIG. 2 is a flowchart of an algorithm of near-field broadband uplink transmission assisted by a dynamic metasurface antenna based on an alternate optimization method and with system sum rate maximization as a criterion according to Embodiment 1.

[0056] FIG. 3 is a flowchart of a Lorentzian constraint phase weight algorithm based on an MM method according to Embodiment 1.

DETAILED DESCRIPTION

[0057] In order to make the objectives, technical solutions, and advantages of the embodiments of the present invention clearer, the technical solutions in the embodiments of the present invention will be clearly and completely described with reference to the accompanying drawings in the embodiments of the present invention. Obviously, it is a part of the embodiments of the present invention, but not all the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by a person of ordinary skill in the art without creative efforts shall fall within the protection scope of the present invention.

Embodiment 1

[0058] Referring to FIG. 1 to FIG. 3, this embodiment provides a near-field broadband uplink MIMO transmission method assisted by a dynamic metasurface antenna. In the method, the dynamic metasurface antenna is arranged on a base station side and configured to observe and capture signals from a channel. In the uplink transmission method, with system sum rate maximization as a criterion, based on methods such as MWMSE minimization, alternate optimization, matrix vectorization and MM, a baseband beamformer and a weight matrix of the dynamic metasurface antenna are determined, and formation of signal beams is completed, thereby improving transmission performance.

[0059] As channel state information in a communication system is changed, the base station side repeats the foregoing steps according to updated channel state information, to perform the near-field broadband uplink MIMO transmission method assisted by a dynamic metasurface antenna, thereby dynamically updating the transmission process, so as to ensure transmission performance.

[0060] Specifically, in this embodiment, the foregoing MWMSE minimization specifically includes:

[0061] A process of solving a sum rate expression of a near-field broadband uplink transmission system is relatively complex, so that in this embodiment, an original sum rate maximization problem is transformed into a mean square error sum minimization problem with MWMSE transformation, and a weighted auxiliary matrix is introduced based on original variables to reduce complexity of transmission optimization.

[0062] Specifically, in this embodiment, the foregoing alternate optimization method includes: for the weight matrix of the dynamic metasurface antenna and the weighted auxiliary matrix that are given, a baseband beamforming matrix is obtained with a closed-form solution; for the weight matrix of the dynamic metasurface antenna and the baseband beamforming matrix that are given, the weighted auxiliary matrix is obtained with a closed-form solution:

[0063] for the baseband beamforming matrix and the weighted auxiliary matrix that are given, the weight matrix of the dynamic metasurface antenna is designed with mean square error sum minimization as a criterion and with methods such as matrix vectorization and MM; and joint optimization of the baseband beamforming matrix and the weight matrix of the dynamic metasurface antenna is alternately implemented until a difference between two adjacent system sum rates is less than a given threshold.

[0064] More specifically, in this embodiment, the foregoing designing the weight matrix of the dynamic metasurface antenna with mean square error sum minimization as a criterion and with methods such as matrix vectorization and MM specifically includes the following steps:

[0065] Terms in the optimization problem that are unrelated to the weight matrix and that may be considered as constants are neglected, to obtain a simplified mean square error sum minimization problem: [0066] a target function is transformed into a matrix tracing form through MWMSE transformation, and the matrix tracing form may be transformed into a vector multiplication form with the matrix vectorization method, where the transformation may eliminate a block structure of the weight matrix of the dynamic metasurface antenna in the mean square error sum minimization problem, to reduce problem solving complexity; and [0067] solving of a weight vector of the dynamic metasurface antenna in four weight feasible domains is considered, and the four feasible domains are an unconstrained (complex plane) weight, an amplitude weight, a binary amplitude weight and a Lorentzian phase constraint weight respectively, where [0068] for such two feasible domains as the unconstrained (complex plane) weight and the amplitude weight, the weight vector is solved with a common convex optimization algorithm; [0069] for such a feasible domain as the binary amplitude weight, the weight vector is solved with brute-force search; and [0070] for such a feasible domain as the Lorentzian phase constraint weight, the weight vector is solved with an MM method.

[0071] More specifically, in this embodiment, the foregoing solving the weight vector with an MM method specifically includes: [0072] through the matrix vectorization method, an optimization variable is simplified from the weight matrix of the dynamic metasurface antenna into the weight vector, and the target function is simplified from matrix optimization into vector optimization; [0073] when a weighted auxiliary variable introduced by MWMSE transformation and the baseband beamforming matrix are considered as constants to solve the weight vector of the dynamic metasurface antenna, the target function is a non-convex function of the weight vector and iteratively solved with the MM method; [0074] in each time of iteration, a target function is replaced with its upper bounding function, a closed expression of an upper bounding problem is given, a target function in the next time of iteration is updated with this solution, a value of the original target function is calculated, and the iteration terminates when a difference between target functions in adjacent two times of iteration is less than a given threshold; and after the termination, the weight vector is changed again into a matrix as a solution to the mean square error sum minimization problem when the baseband beamforming matrix and the weighted auxiliary variable are given.

[0075] In this embodiment, to describe the transmission method more clearly and in more detail, the transmission method is specifically described with a specific application scenario and includes:

[0076] (1) A sum rate maximization problem is constructed based on a broadband large-scale MIMO uplink single-cell system and a channel model that considers a near-field effect, frequency-selective fading, and a spatial broadband effect, and the problem is defined as a first optimization problem, where a dynamic metasurface antenna array is used on a base station side of the broadband large-scale MIMO uplink single-cell system, and the first optimization problem is solved in a manner of jointly designing a baseband beamforming matrix and a weight matrix of a dynamic metasurface antenna, to maximize a near-field broadband large-scale MIMO uplink sum rate, where the step specifically includes:

[0077] As shown in FIG. 1, the method is based on a broadband large-scale MIMO uplink single-cell system, the system includes a plurality of single-antenna users and one base station, and a dynamic metasurface antenna array is used on a base station side as a signal receive antenna. The array is formed by M microwave transmission bands, and L metasurface units are mounted on each microwave transmission band. Therefore, the dynamic metasurface antenna array is formed by a total of N.sub.RML, a cell includes U single-antenna users, a set of user is {1,2, . . . , U}, and N.sub.u antennas are configured for each user.

[0078] Q?.sup.M?N.sup.R represents a weight matrix of a dynamic metasurface antenna, whose expression is:

[00022]$\begin{array}{cc}{\left(Q\right)}_{m_1,\left(m_2-1\right)L+l}=\{\begin{array}{cc}{q}_{m_1.l},& m_1=m_2\\ 0,& m_1?m_2\end{array}& \left(1\right)\end{array}$

[0079] In the expression (1), m.sub.1?{1, 2, . . . , M}, m.sub.2?{1,2, . . . , M}, l?{1, 2, . . . , L}, and q.sub.m.sub.1.sub.l represents a gain of a (l).sup.th antenna unit on a (m).sup.th microwave transmission band for a signal, that is, a change for a signal amplitude or phase. Specifically, metamaterial units may be considered as a resonance circuit, a change of the units for a signal may be modeled into a weight multiplier of an amplitude, a binary amplitude or a Lorentzian constraint phase, and a specific expression is: [0080] amplitude: q?=[a,b],b>a>0; and binary amplitude: q?=c.Math.{0,1},c>0; and Lorentzian constraint phase:

[00023]$q??=\left\{\frac{?+e^{\mathrm{??}}}{2}|??\left[0,2?\right]\right\};$

and where represents an imaginary unit.

[0081] Specifically, the near-field broadband uplink transmission system assisted by the dynamic metasurface antenna has characteristics such as large base station antenna array aperture, high signal carrier frequency, and large transmission bandwidth, these characteristics cause wireless communication to possibly occur in a near-field region of the base station, and meanwhile signal transmission is affected by frequency-selective fading and the spatial broadband effect. Therefore, the following channel model is introduced into this embodiment, and a specific expression is:

[00024]$\begin{array}{cc}{g}_u\left(f\right)=\underset{p=0}{\overset{P_u}{.Math.}}{a}_{u,p}\left(f\right)?b_{u,p}\left(f\right)& \left(2\right)\end{array}$

in the expression (2),

[00025]$\begin{array}{cc}{a}_{u,p}\left(f\right)=\left[\sqrt{{?}_{1,1,u,p}}{A}_{1,1,u,p}\left(f\right),\sqrt{{?}_{1,2,u,p}}{A}_{1,2,u,p}\left(f\right),.Math.,\sqrt{{?}_{M,L,u,p}}{A}_{M,L,u,p}\left(f\right)\right]& \left(3a\right)\end{array}$$\begin{array}{cc}{b}_{u,p}\left(f\right)=\left[\begin{array}{c}{e}^{-\mathrm{?2?}\left(f+f_c\right)\frac{.Math.p_{u,p}-p_{1,1}.Math.}{c}},\\ e^{-\mathrm{?2?}\left(f+f_c\right)\frac{.Math.p_{u,p}-p_{1,2}.Math.}{c}},.Math.,e^{-2?\left(f+f_c\right)\frac{.Math.p_{u,p}-p_{M,L}.Math.}{c}}\end{array}\right]& \left(3b\right)\end{array}$

where a.sub.u,p(f) and b.sub.u,p(f) represent a channel gain that considers a near-field effect, frequency selectivity, and a spatial broadband effect and a response matrix of an antenna array respectively: ?.sub.m,l,u,p and A.sub.m,l,u,p(f) represent a large-scale fading factor of a (p).sup.th transmission path between a (l).sup.th metamaterial on a (m).sup.th microstrip of a base station antenna and a user u and a channel gain coefficient respectively, p.sub.u,p and P.sub.m,l represent a scatterer position of the (p).sup.th transmission path between the user u and the base station and a position of the (l).sup.th metamaterial on the (m).sup.th microstrip of the base station antenna respectively, f and f.sub.c represent a frequency and a center frequency respectively, and c represents a signal transmission speed equal to 3?10.sup.8;

[0082] Specifically, a specific expression of the channel gain coefficient is:

[00026]$\begin{array}{cc}{A}_{m,l,u,p}\left(f\right)=.Math.{?}_{u,p}\left(f\right).Math.\sqrt{F\left({?}_{m,l,u,p}\right)}\frac{c}{4?\left(f+f_c\right){.Math.p_{u,p}-p_{m,l}.Math.}_F}& \left(4\right)\end{array}$

in the expression (4), ?.sub.m,l,u,p=(?.sub.m,l,u,p?.sub.m,l,u,p) represents height-azimuth of a signal reflected from the user u by a (p).sup.th reflector and reaching a (l).sup.th antenna unit on the (m).sup.th microstrip of the base station antenna, and a specific expression of F(?.sub.m,l,u,p) is:

[00027]$\begin{array}{cc}F\left({?}_{m,l,u,p}\right)=\{\begin{array}{cc}6{\cos }^2\left({?}_{m,l,u,p}\right),& {?}_{m,l,u,p}?\left[0,\frac{?}{2}\right]\\ 0,& \mathrm{otherwise}\end{array}& \left(5\right)\end{array}$

?.sub.u,p(f) refers to a reflection coefficient of a reflector on a (p).sup.th path of the user u, and a specific expression thereof is:

[00028]$\begin{array}{cc}{?}_{u,p}\left(f\right)=\{\begin{array}{cc}\frac{\begin{array}{c}\cos {?}_{i,u,p}-\\ n_t\cos {?}_{t,u,p}\end{array}}{\begin{array}{c}\cos {?}_{i,u,p}+\\ n_t\cos {?}_{t,u,p}\end{array}}{e}^{-\left(\frac{8{{?}^2\left(f+f_c\right)}^2{?}_{\mathrm{rough}}^2{\cos }^2{?}_{i,u,p}}{c^2}\right)},& p=1,2,.Math.,P\\ 1,& p=0\end{array}& \left(6\right)\end{array}$

in the expression (6), n.sub.t is a refractive index, ?.sub.rough.sup.2 is a roughness coefficient of a reflection surface, and cos ?.sub.i,u,p and cos ?.sub.t,u,p are an incident angle and a reflection angle of the signal of the user u on a (p).sup.th reflection object respectively.

[0083] To sum up, the sum rate of the system may be expressed as:

[00029]$\begin{array}{cc}{R}_S=\underset{s=0}{\overset{S-1}{.Math.}}{?}_B\log .Math.\begin{array}{c}{I}_U+\\ \frac{P_t}{{?}_B{?}^2}{W}_s^{H}Q{H}_s{G}_s{G}_s^H{H}_s^H{Q}^H{W_s\left(W_s^{H}Q{H}_s{H}_s^H{Q}^H{W}_s\right)}^{-1}\end{array}.Math.& \left(7\right)\end{array}$

in the expression (7), S represents a quantity of subcarriers, and ?.sub.B represents a subcarrier spacing and is expressed as a ratio of a bandwidth B to the quantity of subcarriers S, that is,

[00030]${?}_B=\frac{B}{S};I_U$

is an identity matrix of U?U, ?.sup.2 is a variance of noise, P.sub.t represents a transmit power, and U is a quantity of users in a cell; G.sub.s=[g.sub.1,s,g.sub.2,s, . . . , g.sub.U,s]?.sup.N.sup.R.sup.?U represents a channel matrix of an (s).sup.th subcarrier, W.sub.s?.sup.M?U represents a baseband beamformer of the (s).sup.th subcarrier, H.sub.s?.sup.N.sup.R.sup.?N.sup.R describes a frequency-selective effect of a signal propagated on a microwave transmission band, and Q?.sup.M?N.sup.R represents a weight matrix of a dynamic metasurface antenna; and log is a logarithm operation, and |.Math.| is a matrix determinant obtaining operation.

[0084] Specifically, a baseband beamforming matrix and a weight matrix of a dynamic metasurface antenna are jointly designed, to maximize a near-field broadband large-scale MIMO uplink sum rate, and a specific expression of the foregoing first optimization problem is:

[00031]$\begin{array}{cc}{?}_1\text{:}\underset{Q,W_s,?s?\left\{1,2,.Math.,S\right\}}{\min }{R}_{S}s.t.{\left(Q\right)}_{m_1,\left(m_2-1\right)L+l}=\{\begin{array}{cc}{q}_{m_1,l},& m_1=m_2\\ 0,& m_1?m_2\end{array},q_{m_1,l}??,?m_1,l& \left(8\right)\end{array}$

[0085] In this problem, calculation complexity of the target function is very high, constraints are complex, and a plurality of target matrices need to be jointly optimized.

[0086] Therefore, this embodiment proposes a near-field broadband uplink MIMO transmission method assisted by a dynamic metasurface antenna, including methods such as MWMSE transformation, alternate optimization, matrix vectorization, and MM. Involved algorithms are described in detail below with reference to the foregoing optimization problem model.

[0087] (2) The first optimization problem in step (1) is equivalent to a mean square error sum minimization problem, the problem is defined as a second optimization problem, and then the second optimization problem is solved with the alternate optimization method where system sum rate maximization is used as a criterion, where when the second optimization problem is solved, a weighted auxiliary matrix is introduced based on an original variable to reduce complexity of transmission optimization; and the step (2) specifically includes: [0088] the first optimization problem posed in the step (1) is a typical sum rate maximization problem and is equivalent to a matrix-weighted mean square error sum minimization problem , the problem is defined as the second optimization problem in this embodiment, and a specific expression thereof is:

[00032]$\begin{array}{cc}{?}_2\text{:}\underset{M_s,Q,W_s,?s?\left\{1,2,.Math.,S\right\}}{\min }\underset{s=0}{\overset{S-1}{.Math.}}\mathrm{tr}\left\{M_s{E}_s\left(Q,W_s\right)\right\}-{\log }_2.Math.M_s.Math.& \left(9a\right)\end{array}$$\begin{array}{cc}s.t.{\left(Q\right)}_{m_1,\left(m_2-1\right)L+l}=\{\begin{array}{cc}{q}_{m_1,l},& m_1=m_2\\ 0,& m_1?m_2\end{array}& \left(9b\right)\end{array}$$\begin{array}{cc}{q}_{m_1,l}??,?m_1,l& \left(9c\right)\end{array}$

[0089] In the foregoing expression, M.sub.s is a weighted auxiliary matrix, and E.sub.s(Q, W.sub.s) is a mean square error sum matrix, whose specific expression is:

[00033]$\begin{array}{cc}{E}_s\left(Q,W_s\right)=P_t\left(W_s^H{\mathrm{QH}}_s{G}_s-I_U\right){\left(W_s^H{\mathrm{QH}}_s{G}_s-I_U\right)}^H+{?}_B{?}^2{W}_s^H{\mathrm{QH}}_s{H}_s^H{Q}^H{W}_s& \left(10\right)\end{array}$

[0090] Specifically, In this embodiment, the second optimization problem is solved through alternate optimization, specifically including:

[0091] When Q and W.sub.s are given, where ?s?{1, 2, . . . , S}, M.sub.s may be obtained from the following expression, where ?s?{1, 2, . . . , S}:

[00034]$\begin{array}{cc}{M}_s^{\mathrm{opt}}={E_s\left(Q,W_s\right)}^{-1}& \left(11\right)\end{array}$

[0092] When Q and M.sub.s are given, where ?s?{1, 2, . . . , S}, W.sub.s may be given by the following expression, where ?s?{1, 2, . . . , S}:

[00035]$\begin{array}{cc}{W}_s^{\mathrm{opt}}={\left(P_t{\mathrm{QH}}_s{G}_s{G}_s^H{H}_s^H{Q}^H+{?}_B{?}^2{\mathrm{QH}}_s{H}_s^H{Q}^H\right)}^{-1}{\mathrm{QH}}_s{G}_s& \left(12\right)\end{array}$

[0093] When W.sub.s and M.sub.s are given, where ?s?{1, 2, . . . , S}, Q may be obtained by solving a third optimization problem , and a specific expression of the third optimization problem is:

[00036]$\begin{array}{cc}{?}_3\text{:}\underset{s=0}{\overset{S-1}{.Math.}}\mathrm{tr}\left\{P_t{M}_s{W}_s^H{\mathrm{QH}}_s{G}_s{G}_s^H{H}_s^H{Q}^H{W}_s\right\}-\mathrm{tr}\left\{P_t{M}_s{W}_s^H{\mathrm{QH}}_s{G}_s\right\}-\mathrm{tr}\left\{P_t{M}_s{G}_s^H{H}_s^H{Q}^H{W}_s\right\}+\mathrm{tr}\left\{{?}_B{?}^2{M}_s{W}_s^H{\mathrm{QH}}_s{H}_s^H{Q}^H{W}_s\right\}& \left(13a\right)\end{array}$$\begin{array}{cc}s.t.{\left(Q\right)}_{m_1,\left(m_2-1\right)L+1}=\{\begin{array}{cc}{q}_{m_1,l},& m_1=m_2\\ 0,& m_1?m_2\end{array},& \left(13b\right)\end{array}$$\begin{array}{cc}{q}_{m_1,l}??,?m_1,l& \left(13c\right)\end{array}$

[0094] The third optimization problem is obtained by substituting E.sub.s(Q, W.sub.s) into the second optimization problem , where ?s?{1, 2, . . . , S} and leaving out terms unrelated to Q.

[0095] Specifically, as shown in FIG. 2, this embodiment gives a flow of an algorithm of near-field broadband uplink transmission assisted by a dynamic metasurface antenna based on an alternate optimization method and with system sum rate maximization as a criterion, and a detailed process of the algorithm is as follows: [0096] Step 1, initialize a baseband beamforming matrix W.sub.s.sup.(0), where ?s?{1, 2, . . . , S}, a weight matrix Q.sub.(0) of the dynamic metasurface antenna, a weighted auxiliary matrix M.sub.(0), and a system sum rate R.sub.S.sup.(0), where a quantity of iteration times is .sub.2=0, and a threshold is ?.sub.2. [0097] Step 2, give a weight matrix of the dynamic metasurface antenna, and solve a baseband beamforming matrix according to the expression (12), where ?s?{1, 2, . . . , S}. [0098] Step 3, give and the baseband beamforming matrix , where ?s?{1, 2, . . . , S}, and calculate a mean square error sum matrix E.sub.s(Q, W.sub.s according to the expression (10), where ?s?{1, 2, . . . , S}. [0099] Step 4, give the mean square error sum matrix E.sub.s(Q, W.sub.s, where ?s?{1, 2, . . . , S}, and solve a weighted auxiliary matrix according to the expression (11), where ?s?{1, 2, . . . , S}. [0100] Step 5, give and , where ?s?{1, 2, . . . , S}, and solve a problem to obtain a weight matrix of a dynamic metasurface antenna. [0101] Step 6, calculate a system sum rate , and if |?|??.sub.2 holds, jump out of the loop, and using (,), where ?s?{1, 2, . . . , S}, as a solution meeting the baseband beamforming matrix under the system sum rate maximization criterion and the weight matrix of the dynamic metasurface antenna: otherwise .sub.2=.sub.2+1, and perform step 2 to step 6 again.

[0102] (3) Solve the weight matrix of the dynamic metasurface antenna based on the matrix vectorization method

[0103] Specifically, in this embodiment, for the third optimization problem , the following method is taken:

[0104] For the block structure (13b) of the matrix Q, with the matrix vectorization method, a matrix tracing form in the target function is transformed into a vector multiplication form. Specifically, the matrix Q is pulled into q=[q.sub.1,1,q.sub.1,2, . . . ,q.sub.m,(m?1)L+l, . . . ,q.sub.M,ML].sup.T, where q.sub.m,(m?1)L+l represents an element of a (m).sup.th row and a (l).sup.th column of the matrix Q. Moreover, a matrix vectorization rule is used for obtaining:

[00037]$\begin{array}{cc}\mathrm{tr}\left\{{\mathrm{QC}}_s\right\}=q^T{c}_s,\mathrm{and}\mathrm{tr}\left\{Q^H{C}_s^H\right\}=c_s^{H}q& \left(14a\right)\end{array}$$\begin{array}{cc}\mathrm{tr}\{Q^H{B}_s{\mathrm{QA}}_s=q^H\left(B_s.Math.N_L\right)?A_s^{T}q& \left(14b\right)\end{array}$$\begin{array}{cc}\mathrm{tr}\{Q^H{B}_s{\mathrm{QD}}_s=q^H\left(B_s.Math.I_L\right)?D_{s}q& \left(14c\right)\end{array}$ [0105] where A.sub.s=H.sub.sG.sub.sG.sub.s.sup.HH.sub.s.sup.H, B.sub.s=W.sub.sM.sub.sW.sub.s.sup.H, C.sub.s=H.sub.sG.sub.sM.sub.sW.sub.s.sup.H, and D.sub.s=H.sub.sH.sub.s.sup.H; [0106] L represents a quantity of metasurface units on each microwave transmission band, N.sub.L is an all 1's matrix of L?L, I.sub.L is an identity matrix of L?L, {tilde over (B)}.sub.s is a diagonal matrix, and ({tilde over (B)}.sub.s).sub.m,m=(B.sub.s).sub.m,m; and [0107] c.sub.s=[(C.sub.s).sub.1,1,(C.sub.s).sub.2,1, . . . ,(C.sub.s).sub.(m?1)L+l,m, . . . ,(C.sub.s).sub.ML,M].sup.T.

[0108] The expression (14) is substituted into the target function of the problem , to obtain:

[00038]$f\left(q\right)=\underset{s=0}{\overset{S-1}{.Math.}}{q}^H\left[P_t\left(B_s.Math.N_L\right)?A_s^T+{?}_B{?}^2\left({\tilde{B}}_s.Math.I_L\right)?D_s\right]q-2\mathrm{Re}\left\{P_t{q}^H{c}_s^{*}\right\}$

[0109] Therefore, a problem of optimizing a weight vector of the dynamic metasurface antenna to maximize the system sum rate may be expressed as:

[00039]$\begin{array}{cc}{?}_4:\underset{q}{\min }{q}^H\mathrm{Sq}-2\mathrm{Re}\left\{q^H{c}^{*}\right\}& \left(15a\right)\end{array}$$\begin{array}{cc}s.t.q_{m_1,l}??,?m_1,l& \left(15b\right)\end{array}$$\mathrm{where}S=\underset{s=0}{\overset{S-1}{.Math.}}{P}_t\left(B_s.Math.N_L\right)?A_s^T+{?}_B{?}^2\left({\tilde{B}}_s.Math.I_L\right)?D_s,c=\underset{s=0}{\overset{S-1}{.Math.}}{P}_t{c}_s.$

[0110] For the constraint (15b), four weight feasible domains are considered, are an unconstrained feasible domain, an amplitude feasible domain, a binary amplitude feasible domain and a Lorentzian constraint phase feasible domain respectively, and are specifically as follows: [0111] a. For the unconstrained feasible domain, a problem of optimizing a weight vector of the dynamic metasurface antenna to maximize the system sum rate may be expressed as , whose specific expression is:

[00040]$\begin{array}{cc}{?}_5:\underset{q}{\min }{q}^H\mathrm{Sq}-2\mathrm{Re}\left\{c^{T}q\right\}& \left(16\right)\end{array}$ [0112] is a convex problem, and can be solved with a conventional convex optimization method. [0113] b. For the amplitude feasible domain, a problem of optimizing a weight vector of the dynamic metasurface antenna to maximize the system sum rate may be expressed as , whose specific expression is:

[00041]$\begin{array}{cc}{?}_6:\underset{q}{\min }{q}^T\mathrm{Sq}-2\mathrm{Re}\left\{c^T\right\}q& \left(17a\right)\end{array}$$\begin{array}{cc}s.t.q_{m_1,l}?\left[a,b\right],b>a>0,?m_1,l& \left(17b\right)\end{array}$

[0114] Similarly, is a convex problem, and can be solved with a conventional convex optimization method. [0115] c. For the binary amplitude feasible domain, a problem of optimizing a weight vector of the dynamic metasurface antenna to maximize the system sum rate may be expressed as , whose specific expression is:

[00042]$\begin{array}{cc}{?}_7:\underset{q}{\min }{q}^T\mathrm{Sq}-2\mathrm{Re}\left\{c^T\right\}q& \left(18a\right)\end{array}$$\begin{array}{cc}s.t.q_{m_1,l}?c.Math.\left\{0,1\right\},c>0,?m_1,l& \left(18b\right)\end{array}$ [0116] can be solved with a brute-force search method. [0117] d. For the Lorentzian constraint phase weight, a problem of optimizing a weight vector of the dynamic metasurface antenna to maximize the system sum rate may be expressed as , whose specific expression is:

[00043]$\begin{array}{cc}{?}_8:\underset{q}{\min }{q}^H\mathrm{Sq}-2\mathrm{Re}\left\{c^{T}q\right\}& \left(19a\right)\end{array}$$\begin{array}{cc}s.t.q_{m_1,l}?\left\{\frac{?+e^{\mathrm{??}}}{2}|??\left[0,2?\right]\right\},?m_1,l& \left(19b\right)\end{array}$

where represents an imaginary unit.

[0118] Specifically, in this embodiment, can be solved with an MM method (an ordered convex optimization method), specifically including: First, a tractable valid upper bound is found, the problem is replaced with a problem about an upper bound replacement function, and then a weight vector of a dynamic metasurface antenna is obtained with an alternate optimization method. An algorithm of solving a Lorentzian constraint phase weight with an MM method is described in detail below.

[0119] More specifically, the foregoing solving a Lorentzian constraint phase weight with an MM method specifically includes the following steps:

[0120] First, the weight vector of the dynamic metasurface antenna is expressed as

[00044]$q=\frac{1}{2}\left({\mathrm{?I}}_{N_R}+p\right),$

where 1.sub.N.sub.R is an all 1's vector, and

[00045]$p={\left[e^{?{?}_1},e^{{\mathrm{??}}_2},.Math.,e^{?{?}_{N_R}}\right]}^T.$

Therefore, the Lorentzian constraint phase is simplified into a modulo-1 constraint phase.

[0121] In this case, the expression (19a) may be written as:

[00046]$f\left(p\right)=\frac{1}{4}{p}^H\mathrm{Sp}+\mathrm{Re}\left\{\frac{?}{2}{p}^{H}S{1}_{N_R}-p^H{c}^{*}\right\}+\mathrm{Re}\left\{{\mathrm{?1}}_{N_R}^T{c}^{*}\right\}+\frac{1}{4}{1}_{N_R}^{T}S{1}_{N_R}$

[0122] The function f (p) is a non-convex quadratic function about P, and with the MM method that is a sequential convex optimization method, a compact upper bounding function of f(p) may be obtained. First, a tractable compact upper bounding function is found and expressed as:

[00047]$\begin{array}{cc}{p}^H\mathrm{Sp}?p^H\mathrm{Tp}-2\mathrm{Re}\left\{p^H\left(T-S\right)p^{\left(?\right)}\right\}+{\left(p^{\left(?\right)}\right)}^H)\left(T-S\right)p^{\left(?\right)}& \left(20\right)\end{array}$

T=?.sub.maxI, ?.sub.max is a maximum eigenvalue of S, and the expression (20) is substituted into the problem to obtain a problem , whose specific expression is:

[00048]$\begin{array}{cc}{?}_9^{\left(?\right)}\underset{p}{\max }\mathrm{Re}\left\{p^H{a}^{\left(?\right)}\right\}& \left(21a\right)\end{array}$$\begin{array}{cc}s.t.p_n?\left\{e^{\mathrm{??}}|??\left[0,2?\right]\right\},?n& \left(21b\right)\end{array}$

where =(?.sub.max I.sub.N.sub.R?S)+2c*?. Additionally, in the target function of the problem , terms unrelated to the variable p are left out. The problem , can be solved by alternately optimizing the vectors a and p. Moreover, in each time of iteration, the vectors a and p may be obtained through a closed-form solution.

[0123] FIG. 3 gives a flow of an algorithm of solving a Lorentzian constraint phase weight based on an MM method, and a detailed process of the algorithm is specifically: [0124] Step 1, initialize and , and set an iteration index .sub.1=0 and a threshold ?.sub.1. [0125] Step 2, give , and calculate =and ?n. [0126] Step 3, give , and calculate =(?.sub.maxI.sub.N.sub.R?S)+2c*?. [0127] Step 4, calculate the system sum rate , and if a difference between the (.sub.1).sup.th system sum rate and the (.sub.1+1).sup.th system sum rate is less than the given threshold ?.sub.1, jump out of the loop, using as a solution meeting the Lorentzian constraint phase weight under the system sum rate maximization criterion when the baseband beamforming matrix is given: otherwise .sub.1=.sub.1+1, and perform steps 2 to 4 again.

[0128] As channel state information in a communication system is changed, the base station side dynamically implements, according to updated channel state information, near-field broadband large-scale MIMO uplink transmission assisted by a dynamic metasurface antenna with system sum rate maximization as a criterion, thereby dynamically updating transmission, to ensure transmission performance.

[0129] Any content not described in detail in the present invention is a technology publicly known by a person skilled in the art. The specific exemplary embodiments of the present invention are described in detail above. It should be understood that, a person of ordinary skill in the art may make many modifications and changes according to the idea of the present invention without creative effort. Therefore, any technical solution that may be obtained by a person skilled in the art through logic analysis, reasoning or limited experiments based on the existing technology according to the idea of the present invention should fall within the protection scope determining by the claims.