Inductive-pass integrated hybrid beam forming method based on one-bit digital-to-analog converter

Abstract

The invention relates to an inductance-pass integrated hybrid beam forming method based on a one-bit digital-to-analog converter. The method comprises the following steps: designing an analog precoder and a one-bit digital-to-analog converter quantized signal passing through a digital precoder based on a trade-off optimization problem between radar and communication performance; a transmitting beam vector is designed according to the designed analog precoder and the quantized signal of the one-bit digital-to-analog converter, an optimization expression of the transmitting beam vector is solved, and an optimal transmitting beam vector is obtained; and designing an optimization problem of quantized signals of the analog precoder and the one-bit digital-to-analog converter according to the optimal transmitting beam vector, solving the optimization problem by using a Matlab toolbox and a BTSA-GA algorithm, and generating the sensing and communication integrated hybrid beam according to the optimized quantized signals of the analog precoder and the one-bit digital-to-analog converter. By adopting the method, the complexity and the power consumption of the large-scale MIMO system can be reduced.

Description

Technical Field

[0001] The present application relates to the field of communication technology, and in particular to a hybrid beamforming method based on a one-bit digital-to-analog converter and integrated sensing. Background Art

[0002] In recent years, the demand for integrated sensing and communication (ISAC) smart devices in the commercial field has grown rapidly. With the continuous development of radar and wireless communication technology, and people's high expectations for the multifunctionality and intelligence of embedded small devices, ISAC systems have entered the public eye and aroused extensive discussions in industry and academia. The two systems use the same platform, share the same hardware and spectrum resources, and form a specific type of RF technology. Despite its good communication and radar detection performance, the ISAC system with millimeter wave and massive MIMO as the core technology requires a large number of hardware components, such as RF chain, power amplifier, digital-to-analog converter, etc., which will undoubtedly lead to a significant increase in the total power consumption and hardware cost of the base station. Pure digital beamforming directly designs the quantized transmission signal vector, and each antenna maps an RF chain. It is suitable for systems with a small number of antennas. For massive MIMO systems, in order to reduce system complexity and power consumption, a hybrid beamforming design is required. Summary of the invention

[0003] Based on this, it is necessary to provide a hybrid beamforming method based on a one-bit digital-to-analog converter integrated sensing that can reduce the complexity and power consumption of large-scale MIMO systems in order to address the above technical problems.

[0004] A hybrid beamforming method based on a one-bit digital-to-analog converter integrated sensing and communication, the method comprising:

[0005] The signal sent by the base station is obtained, and the signal sent by the base station is modeled to obtain a signal model; the signal model includes an analog precoder to be solved; based on the trade-off optimization problem between radar and communication performance, an analog precoder and a one-bit digital-to-analog converter quantized signal after the digital precoder are designed;

[0006] Designing a transmit beam vector according to the designed analog precoder and a one-bit digital-to-analog converter quantized signal, and establishing an optimized expression for the transmit beam vector;

[0007] Solve the optimization expression of the transmit beam vector, expand the cost function of the optimization expression, express the expanded expression in Lagrangian form, and obtain the optimality condition of the expanded expression based on the Lagrangian form; solve the optimality condition to obtain the optimal transmit beam vector;

[0008] According to the optimal transmit beam vector, the optimization problem of designing the analog precoder and the one-bit digital-to-analog converter quantization signal is solved by using Matlab toolbox and BTSA-GA algorithm respectively to obtain the optimized analog precoder and the one-bit digital-to-analog converter quantization signal; based on the optimized analog precoder and the one-bit digital-to-analog converter quantization signal, a sensing-communication integrated hybrid beam is generated.

[0009] The above hybrid beamforming method based on the integration of one-bit digital-to-analog converter sensing and communication is different from the design idea of ​​directly designing the quantized transmission signal vector in the digital beamforming scheme. This application designs an analog precoder and a one-bit digital-to-analog converter quantized signal after the digital precoder. Given the optimal communication beamformer and the desired radar beam pattern, the analog precoder and the one-bit digital-to-analog converter quantized signal after the digital precoder are designed. Under the constraints of the system transmit power, the constant modulus constraint of the analog precoder, and the discrete non-convex constraint of the quantized signal, the weighted sum of the communication beamforming error and the radar beamforming error is minimized. The ISAC system can obtain a better peak-to-peak ratio of the radar sensing waveform and superior communication bit error rate performance. Then, according to the optimal transmit beam vector, the optimization problem of the analog precoder and the one-bit digital-to-analog converter quantization signal is designed. The Matlab toolbox and the BTSA-GA algorithm are used to alternately solve the analog precoder under the constant envelope constraint of the phase-shifting network and the one-bit digital-to-analog converter quantization signal after digital precoding. A global optimization algorithm, namely the BTSA-GA algorithm, is designed to solve the discrete constrained optimization problem, avoiding the performance loss caused by relaxing the non-convex symbol constraints in some classic nonlinear precoding algorithms. According to the optimized analog precoder and the one-bit digital-to-analog converter quantization signal, a sensing-integrated hybrid beam is generated, and the detection and perception functions of the radar are realized while the base station communicates with multiple downlink users. The scarce spectrum resources are fully utilized to achieve larger channel capacity, higher spectrum efficiency and better target detection performance, which can greatly reduce the hardware cost, system complexity and power consumption of each RF link. BRIEF DESCRIPTION OF THE DRAWINGS

[0010] Figure 1 It is a flowchart of a hybrid beamforming method based on a one-bit digital-to-analog converter integrated sensing in one embodiment;

[0011] Figure 2 A diagram of a massive MIMO ISAC system model based on a one-bit digital-to-analog converter in one embodiment;

[0012] Figure 3 The figure is a flowchart of the overall process of the algorithm of the present application in one embodiment.

[0013] Figure 4 is a flowchart of a BTSA-GA algorithm in one embodiment;

[0014] Figure 5 A diagram showing the system reachability and rate comparison of the method of the present application in one embodiment;

[0015] Figure 6 A comparison diagram of bit error rates of the method of the present application in one embodiment;

[0016] Figure 7 FIG. 4 is a radar beam pattern diagram of the method of the present application in one embodiment. DETAILED DESCRIPTION

[0017] In order to make the purpose, technical solution and advantages of the present application more clearly understood, the present application is further described in detail below in conjunction with the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are only used to explain the present application and are not used to limit the present application.

[0018] In one embodiment, Figure 1 As shown, a hybrid beamforming method based on a one-bit digital-to-analog converter and sensing integration is provided, which can be applied to Figure 2 The massive MIMO ISAC system based on a one-bit digital-to-analog converter includes the following steps:

[0019] Step 102, obtain the signal sent by the base station, model the signal sent by the base station, and obtain a signal model; the signal model includes an analog precoder to be solved; based on the trade-off optimization problem between radar and communication performance, the analog precoder and the one-bit digital-to-analog converter quantized signal after the digital precoder are designed.

[0020] Narrowband massive MIMO system, where dual-function base stations are equipped with N T Transmitting antennas and N r receiving antennas, and N RF The base station sends data streams to downlink users. At the same time, it supports the detection and perception of k targets in the same time-frequency resources.

[0021] set up is the downlink channel matrix from the base station to the user equipment, then the received signal at the user end can be modeled as

[0022]

[0023] In the formula, is the power control parameter, To simulate the precoder, is the digital precoder, is a one-bit quantization, n is zero mean, and variance is σ 2 The observation noise consists of independent complex Gaussian components. The data flow vector satisfies Where P s is the total power of the unprecoded data, N s ×N s The identity matrix of .

[0024] The mmWave channel is modeled as having N c Scattering clusters and each scattering cluster contains N p The propagation environment of the scattering path. In this model, the channel matrix is ​​as follows

[0025]

[0026] Among them, α lr is the small-scale fading coefficient, satisfying the standard complex Gaussian distribution φ lr and θ lr They are the angle of arrival (AoA) and the angle of departure (AoD). and is the array steering vector that depends only on the transmit and receive antenna array structures.

[0027]

[0028] A single-bit DAC senses the signal The elements in are:

[0029]

[0030] A typical MIMO radar beamformer can be constructed by combining a H (θ)C Rad a(θ) is matched with the desired transmit beam pattern φ(θ) to solve, where C Rad is the covariance matrix of the detection signal. Mathematically, solving C Rad The constrained least squares problem is expressed as:

[0031]

[0032] Among them, θ i is a fine grid of points covering the target direction of arrival (DOA), N is the number of points, and δ is the scale parameter. This is a convex optimization problem that can be solved directly using the Matlab toolbox CVX.

[0033] The trade-off optimization problem between radar and communication performance can be described as follows:

[0034]

[0035] Among them, F FDis the optimal unconstrained precoding matrix, which consists of the first n columns of the right singular matrix of H. ρ∈[0,1] is a weighting factor that balances radar performance and communication performance.

[0036] Step 104, designing a transmit beam vector according to the designed analog precoder and a one-bit digital-to-analog converter quantized signal, and establishing an optimized expression of the transmit beam vector.

[0037] First, design the transmit beam vector The optimization problem of y can be equivalent to:

[0038]

[0039] in,

[0040] The transmit beam vector is designed according to the designed analog precoder and a one-bit digital-to-analog converter quantized signal, and an optimized expression of the transmit beam vector is established to ensure that it can meet the dual requirements of radar and communication.

[0041] Step 106, solving the optimization expression of the transmit beam vector, expanding the cost function of the optimization expression, expressing the expanded expression in Lagrangian form, and obtaining the optimality condition of the expanded expression based on the Lagrangian form; solving the optimality condition to obtain the optimal transmit beam vector.

[0042] Define Q = D H D, G = D H B, the cost function of (8) is further expanded as:

[0043]

[0044] The optimization problem (9) is a matrix version of the trust region subproblem (TRS) that satisfies strong duality, with a duality gap of zero. Its Lagrangian form is expressed as:

[0045]

[0046] where λ is the dual variable associated with the equality constraint. The optimality condition of (9) is

[0047]

[0048]

[0049] in, λ opt is the optimal dual variable, y opt It is the optimal solution to the dual problem and also the optimal solution to the original optimization problem. (12) and (13) guarantee the feasibility of the original problem and the feasibility of the dual problem respectively.

[0050] based on Obviously, It is the matrix M H The eigenvalue of +M, and the corresponding eigenvector satisfying (12) and (13), is exactly the optimal solution y opt .

[0051] The cost function is expanded into Lagrangian form to obtain the optimality condition and then solve the optimal transmit beam vector.

[0052] Design an analog precoder, and a one-bit digital-to-analog converter (1-bit digital-to-analog converter) to sense the quantized signal after the digital precoder. Given the optimal communication beamformer and the desired radar beam pattern, design an analog precoder and a one-bit digital-to-analog converter to sense the quantized signal after the digital precoder, and minimize the weighted sum of the communication beamforming error and the radar beamforming error under the constraints of the system transmit power, the constant modulus constraint of the analog precoder, and the discrete non-convex constraint of the quantized signal. The ISAC system can obtain a better peak-to-side ratio (PSLR) of the radar sensing waveform and superior communication bit error rate performance.

[0053] Step 108, design the optimization problem of the analog precoder and the one-bit digital-to-analog converter quantization signal according to the optimal transmit beam vector, use the Matlab toolbox and the BTSA-GA algorithm to solve the optimization problem respectively, and obtain the optimized analog precoder and the one-bit digital-to-analog converter quantization signal; generate an integrated hybrid beam based on the optimized analog precoder and the one-bit digital-to-analog converter quantization signal.

[0054] fixed F A The optimization problem can be written as:

[0055]

[0056] The optimization problem (14) is a typical Riemannian manifold optimization problem, which can be solved using the Matlab toolbox Manopt. The direction in which the cost function in the Riemannian manifold decreases the most at x is its Riemannian gradient:

[0057]

[0058] The cost function is its Euclidean gradient, described in detail as follows:

[0059]

[0060] Fixed F A , The optimization problem can be written as:

[0061]

[0062] By quantizing the real part and the imaginary part by one bit respectively, the above optimization problem can be further rewritten as:

[0063]

[0064] in, Represents the set {1,2,...,2N RF},

[0065]

[0066] The optimization problem (18) is solved by binary tree species-genetic optimization algorithm (BTSA-GA). The specific process of the algorithm is as follows.

[0067] In the initialization phase of the BTSA-GA algorithm, a batch of trees are generated:

[0068] T i,j = l j +r i,j (h j -l j ), (twenty two)

[0069] Among them, l j is the lower bound of the search space -1, h j is the upper bound of the search space 1, r i,j It is a random number generated for each dimension at each position in the range [0,1].

[0070] To T i,j Binary discretization:

[0071]

[0072] Find the current best solution by:

[0073] S=minf(T i ),i∈1,2,…,2N RF , (twenty four)

[0074]

[0075] Trees generate seeds, and the positions of trees and seeds in the n-dimensional search space can be viewed as possible solutions to the optimization problem. The positions of seeds are obtained through two search equations:

[0076]

[0077] Among them, S i,j is the jth dimension of the i-th seed generated by the i-th tree, Ti,j is the jth dimension of the i-th tree, B j is the jth dimension of the best tree position obtained so far, T r,j is the jth dimension of the rth tree randomly selected from the population, β i,j is a randomly generated scaling factor in the range [-1,1]. The first equation in (22) focuses on global search to prevent the algorithm from entering the local optimal solution, and the second equation focuses on local search to facilitate rapid convergence of the algorithm. γ is a random number generated for each position and dimension in the range (0,1), and ST is a selection control parameter.

[0078] In order to avoid the tendency of iteration moving to the local optimum, the threshold H is set. i The number of generations of generating seeds exceeds H and there is no better seed to update T i When T is updated through crossover and mutation operations in the genetic algorithm i , expand the search space, and thus reach the global optimal solution.

[0079] According to the optimal transmit beam vector, an optimization problem of designing an analog precoder and a one-bit digital-to-analog converter quantization signal is proposed. The optimization problem is solved using the Matlab toolbox and the BTSA-GA algorithm respectively to obtain the optimized analog precoder and the one-bit digital-to-analog converter quantization signal. The optimized analog precoder and the one-bit digital-to-analog converter quantization signal can generate a more accurate integrated hybrid beam, thereby providing a larger channel capacity, higher spectrum efficiency and better target detection performance. At the same time, the one-bit digital-to-analog converter is deployed in a large-scale MIMO ISAC system, which greatly reduces the hardware cost, system complexity and power consumption of each RF link. The non-convex subproblem of the one-bit digital-to-analog converter quantization signal is solved. The BTSA-GA algorithm proposed in this application has lower performance loss than the linear quantization precoder.

[0080] In the above-mentioned hybrid beamforming method based on one-bit digital-to-analog converter sensing integration, different from the design idea of ​​directly designing the quantized transmission signal vector in the digital beamforming scheme, the present application designs an analog precoder, and a one-bit digital-to-analog converter quantized signal after the digital precoder. Given the optimal communication beamformer and the desired radar beam pattern, the analog precoder and the one-bit digital-to-analog converter quantized signal after the digital precoder are designed. Under the constraints of the system transmit power, the constant modulus constraint of the analog precoder, and the discrete non-convex constraint of the quantized signal, the weighted sum of the communication beamforming error and the radar beamforming error is minimized. The ISAC system can obtain a better peak-to-peak ratio of the radar sensing waveform and superior communication bit error rate performance. Then, according to the optimal transmit beam vector, the optimization problem of the analog precoder and the one-bit digital-to-analog converter quantization signal is designed. The Matlab toolbox and the BTSA-GA algorithm are used to alternately solve the analog precoder under the constant envelope constraint of the phase-shifting network and the one-bit digital-to-analog converter quantization signal after digital precoding. A global optimization algorithm, namely the BTSA-GA algorithm, is designed to solve the discrete constrained optimization problem, avoiding the performance loss caused by relaxing the non-convex symbol constraints in some classic nonlinear precoding algorithms. According to the optimized analog precoder and the one-bit digital-to-analog converter quantization signal, a sensing-integrated hybrid beam is generated, and the detection and perception functions of the radar are realized while the base station communicates with multiple downlink users. The scarce spectrum resources are fully utilized to achieve larger channel capacity, higher spectrum efficiency and better target detection performance, which can greatly reduce the hardware cost, system complexity and power consumption of each RF link.

[0081] In one embodiment, an analog precoder and a one-bit digital-to-analog converter quantized signal after a digital precoder are designed based on a trade-off optimization problem between radar and communication performance, including:

[0082] The analog precoder and the one-bit DAC quantized signal after the digital precoder are modeled as the objective function of the trade-off optimization problem between radar and communication performance. The trade-off optimization problem model between radar and communication performance is established under the constraints of system transmit power constraint, constant modulus constraint of analog precoder and discrete non-convex constraint of quantized signal.

[0083]

[0084] Among them, F A represents the analog precoder, represents a one-bit DAC quantized signal, ρ∈[0,1] is a weighting factor that weighs radar performance and communication performance, the superscript H represents conjugate transpose, and F FD is the optimal unconstrained precoding matrix, C Rad represents the radar beam covariance matrix, represents an element in a one-bit DAC quantized signal, represents the quantized value set of a one-bit digital-to-analog converter, represents the set {1,2,...,N RF},N T represents the number of transmitting antennas of the dual-function base station, n represents the column subscript index of the element in the mth row and nth column in the analog precoder, and m represents the row subscript index of the element in the mth row and nth column in the quasi-precoder.

[0085] In one embodiment, a transmit beam vector is designed according to the designed analog precoder and a one-bit digital-to-analog converter quantized signal, and an optimized expression of the transmit beam vector is established, including:

[0086] According to the designed analog precoder and one-bit DAC quantized signal, the transmit beam vector is designed as: The optimal expression for establishing the transmit beam vector is:

[0087]

[0088] in, F A represents the analog precoder, represents a one-bit DAC quantized signal, ρ∈[0,1] is a weighting factor for weighing radar performance and communication performance, the superscript H represents the transpose operation, and F FD is the optimal unconstrained precoding matrix, C Rad represents the radar beam covariance matrix, P t It represents the total transmit power of the dual-function base station, and the superscript T represents transposition.

[0089] In one embodiment, the cost function of the optimization expression is expanded, including:

[0090] Define Q = D H D, G = D H B. Expand the cost function of the optimization expression to

[0091]

[0092] Where y represents the transmit beam vector, ρ∈[0,1] is a weighting factor that weighs radar performance and communication performance. The superscript H represents the transposition operation. FD is the optimal unconstrained precoding matrix, C Rad represents the radar beam covariance matrix, P t It represents the total transmit power of the dual-function base station, and the superscript T represents transposition.

[0093] In one embodiment, the expanded expression is expressed in Lagrangian form, including:

[0094] The expanded expression is expressed in Lagrangian form as

[0095]

[0096] where λ is the dual variable associated with the equality constraint.

[0097] In one embodiment, obtaining the optimality condition of the expanded expression based on the Lagrangian form includes:

[0098] The optimality condition for the expanded expression based on the Lagrangian form is:

[0099]

[0100] Where Q = D H D, G = D H B, λ opt is the optimal dual variable, y opt is a candidate solution for the optimal transmit beam vector, ρ∈[0,1] is a weighting factor that weighs radar performance and communication performance, the superscript H represents the transpose operation, and F FD is the optimal unconstrained precoding matrix, C Rad represents the radar beam covariance matrix, P t Indicates the total transmit power of the dual-function base station, Indicates that the dimension is 2N T ×2N T The superscript T indicates the transpose.

[0101] In one embodiment, the optimization problem of designing an analog precoder and a one-bit digital-to-analog converter quantized signal according to an optimal transmit beam vector includes:

[0102] The optimization problem of designing the analog precoder based on the optimal transmit beam vector by fixing the one-bit DAC quantization signal is:

[0103]

[0104] Among them, F A represents the analog precoder, represents a one-bit DAC quantized signal, y opt is the candidate solution for the optimal transmit beam vector, N T represents the number of transmitting antennas of the dual-function base station, and n represents the analog precoder F A The element in the mth row and nth column [F A ]m,n The column subscript index of m represents the analog precoder F A The element in the mth row and nth column [F A ] m,n The row subscript index of .

[0105] In one embodiment, the analog precoder is fixed, and the optimization problem of designing a one-bit digital-to-analog converter quantized signal according to the optimal transmit beam vector is:

[0106]

[0107] in, represents an element in a one-bit DAC quantized signal, represents the quantized value set of a one-bit digital-to-analog converter, represents the set {1,2,...,N RF}.

[0108] The real part and the imaginary part are quantized by one bit respectively. The optimization problem of the above one-bit DAC quantized signal is further rewritten as:

[0109]

[0110] in, Represents the set {1,2,...,2N RF},

[0111]

[0112] It represents the operator for taking the real part of a complex number. Represents the operator for taking the imaginary part of a complex number.

[0113] In one embodiment, the optimization problem is solved using the Matlab toolbox and the BTSA-GA algorithm to obtain an optimized analog precoder and a one-bit digital-to-analog converter quantized signal, including:

[0114] The optimization problem of the analog precoder is solved using the Matlab toolbox Manopt to obtain the optimized analog precoder.

[0115] In one embodiment, the BTSA-GA algorithm is used to solve the optimization problem of a one-bit digital-to-analog converter quantized signal. In the initialization stage of the BTSA-GA algorithm, a batch of trees are generated as follows:

[0116] T i,j = l j +r i,j (h j -l j ),

[0117] Among them, l j is the lower bound of the search space -1, h j is the upper bound of the search space 1, r i,j It is a random number generated for each dimension at each position in the range [0,1];

[0118] To T i,j Binary discretization:

[0119]

[0120] Find the current best solution by:

[0121] S = min f (T i ),i∈1,2,…,2N RF ,

[0122]

[0123] Trees generate seeds, and the positions of trees and seeds in the n-dimensional search space are regarded as possible solutions to the optimization problem. The positions of seeds are obtained through two search equations:

[0124]

[0125] Among them, S i,j is the jth dimension of the i-th seed generated by the i-th tree, T i,j is the jth dimension of the i-th tree, B j is the jth dimension of the best tree position obtained so far, T r,j is the jth dimension of the rth tree randomly selected from the population, β i,j is a randomly generated scaling factor in the range [-1,1], γ is a random number generated for each position and dimension in the range (0,1), and ST is a selection control parameter;

[0126] Set the threshold H, when T i The number of generations of generating seeds exceeds H and there is no better seed to update T i When T is updated through crossover and mutation operations in the genetic algorithm i , and obtain the global optimal solution, that is, the optimized one-bit DAC quantized signal.

[0127] In a specific embodiment, in combination Figure 3 As shown, based on the embodiment of the present application, a hybrid beamforming method for integrating sensing and communication of a digital-to-analog converter is proposed in the present application. A two-variable alternating minimization algorithm is firstly designed to transmit the beam y of the MIMO ISAC receiver. opt Then, according to the designed transmit beam, the analog precoder F under the constant envelope constraint of the phase shift network is solved alternately.A and the 1-bit quantized signal after digital precoding because The solution needs to obtain the quantized signal after the digital precoder through exhaustive search Second, when there are many RF chains, the solution process is too complicated. To address this problem, the present application designs a global optimization algorithm, namely the improved binary tree seed genetic optimization algorithm (BTSA-GA) algorithm to solve the discrete constrained optimization problem. The algorithm flow chart is as follows: Figure 4 shown.

[0128] like Figure 5 As shown in FIG. 1 , this application compares the achievable communication and rate of the proposed one-bit hybrid beamforming design method with other linear one-bit digital beamforming design methods based on Monte Carlo experiments. Consider a medium-sized 32×4 MIMO system, from Figure 5 It can be clearly seen that as the signal-to-noise ratio increases, the achievable communication and rate also increases accordingly. Among the precoding methods, the zero forcing algorithm (ZF) is a linear precoding method with lower complexity. Obviously, the achievable communication and rate of the nonlinear precoder with a one-bit digital-to-analog converter designed in this application is better than the linear precoding method, and as the signal-to-noise ratio increases, the gap in achievable and rate also increases accordingly. When the signal-to-noise ratio is 10dB, the achievable communication and rate of the nonlinear precoder is almost twice that of the linear precoding method. The algorithm proposed in this application has a achievable communication and rate performance close to that of the nonlinear digital beamforming design method, and has significant advantages in achievable and rate performance compared to the linear digital beamforming design method.

[0129] Figure 6 This is a comparison chart of the bit error rate performance of the one-bit hybrid beamforming design method proposed in this application and the bit error rate of other precoding schemes. Compared with the one-bit linear digital beamforming design scheme, such as ZF (zero forcing)-1bit, the one-bit hybrid beamforming design method proposed in this application has significant advantages in bit error rate performance.

[0130] In order to evaluate the radar detection and sensing performance, this application sets three target directions of interest, and the beam tracking angles are [θ1 = -30°, θ2 = 0°, θ1 = 30°], and the radar beam pattern is normalized and measured in dB. Figure 7 As shown, the three peaks correspond to beams pointing in the directions of three targets of interest, indicating that the sensory-communication integrated hybrid beamforming method proposed in this application meets the detection and perception functions of the radar. Figure 7 It indicates that the quantization error makes the beam pattern quality slightly worse based on a one-bit DAC, and a performance loss of less than 1dB is tolerable.

[0131] It should be understood that although Figure 1 The steps in the flowchart are shown in sequence as indicated by the arrows, but these steps are not necessarily executed in the order indicated by the arrows. Unless otherwise specified in this document, there is no strict order restriction for the execution of these steps, and these steps can be executed in other orders. Moreover, Figure 1 At least part of the steps may include multiple sub-steps or multiple stages. These sub-steps or stages are not necessarily executed at the same time, but can be executed at different times. The execution order of these sub-steps or stages is not necessarily sequential, but can be executed in turn or alternately with other steps or at least part of the sub-steps or stages of other steps.

[0132] The technical features of the above embodiments may be combined arbitrarily. To make the description concise, not all possible combinations of the technical features in the above embodiments are described. However, as long as there is no contradiction in the combination of these technical features, they should be considered to be within the scope of this specification.

[0133] The above-mentioned embodiments only express several implementation methods of the present application, and the descriptions thereof are relatively specific and detailed, but they cannot be understood as limiting the scope of the invention patent. It should be pointed out that, for a person of ordinary skill in the art, several variations and improvements can be made without departing from the concept of the present application, and these all belong to the protection scope of the present application. Therefore, the protection scope of the patent of the present application shall be subject to the attached claims.

Claims

1. A hybrid beamforming method based on a one-bit digital-to-analog converter integrated sensing, the method is applied to a narrowband massive MIMO system including a dual-function base station, characterized in that: The method comprises: Acquire a signal sent by a base station, model the signal sent by the base station, and obtain a signal model; the signal model includes an analog precoder to be solved; design the analog precoder and a one-bit digital-to-analog converter quantized signal after the digital precoder based on a trade-off optimization problem between radar and communication performance; Designing a transmit beam vector according to the designed analog precoder and a one-bit digital-to-analog converter quantized signal, and establishing an optimized expression for the transmit beam vector; Solving the optimization expression of the transmit beam vector, expanding the cost function of the optimization expression, expressing the expanded expression in Lagrangian form, and obtaining the optimality condition of the expanded expression based on the Lagrangian form; solving the optimality condition to obtain the optimal transmit beam vector; According to the optimal transmit beam vector, the optimization problem of designing the analog precoder and the one-bit digital-to-analog converter quantization signal is solved respectively by using Matlab toolbox and BTSA-GA algorithm to obtain the optimized analog precoder and the one-bit digital-to-analog converter quantization signal; and a sense-communication integrated hybrid beam is generated according to the optimized analog precoder and the one-bit digital-to-analog converter quantization signal.

2. The method according to claim 1, characterized in that: The analog precoder and the one-bit digital-to-analog converter quantized signal after the digital precoder are designed based on the trade-off optimization problem between radar and communication performance, including: The analog precoder and the one-bit DAC quantized signal after the digital precoder are modeled as the objective function of the trade-off optimization problem between radar and communication performance. The trade-off optimization problem model between radar and communication performance is established under the constraints of system transmit power constraint, constant modulus constraint of the analog precoder and discrete non-convex constraint of the quantized signal: Among them, F A represents the analog precoder, represents a one-bit DAC quantized signal, ρ∈[0,1] is a weighting factor that weighs radar performance and communication performance, the superscript H represents conjugate transpose, and F FD is the optimal unconstrained precoding matrix, C Rad represents the radar beam covariance matrix, P t Indicates the total transmit power of the dual-function base station, represents an element in a one-bit DAC quantized signal, represents a one-bit DAC quantization set, represents the set {1,2,...,N RF },N RF Indicates the number of RF chains of dual-function base stations, N T represents the number of transmitting antennas of the dual-function base station, and n represents the analog precoder F A The element in the mth row and nth column [F A ] m,n The column subscript index of m represents the analog precoder F A The element in the mth row and nth column [F A ] m,n The row subscript index of .

3. The method according to claim 1, characterized in that A transmit beam vector is designed according to the analog precoder and a one-bit digital-to-analog converter quantized signal obtained by the design, and an optimization expression of the transmit beam vector is established, including: According to the designed analog precoder and one-bit DAC quantized signal, the transmit beam vector is designed as: The optimal expression for establishing the transmit beam vector is: in, F A represents the analog precoder, represents a one-bit DAC quantized signal, ρ∈[0,1] is a weighting factor for weighing radar performance and communication performance, the superscript H represents the transpose operation, and F FD is the optimal unconstrained precoding matrix, C Rad represents the radar beam covariance matrix, P t It represents the total transmit power of the dual-function base station, and the superscript T represents transposition.

4. The method according to claim 1, characterized in that: Expand the cost function of the optimization expression, including: Define Q = D H D, G = D H B. Expand the cost function of the optimization expression to Where y represents the transmit beam vector, ρ∈[0,1] is a weighting factor that weighs radar performance and communication performance. The superscript H represents the transposition operation. FD is the optimal unconstrained precoding matrix, C Rad represents the radar beam covariance matrix, P t It represents the total transmit power of the dual-function base station, and the superscript T represents transposition.

5. The method according to claim 4, characterized in that The expanded expression is expressed in Lagrangian form, including: The expanded expression is expressed in Lagrangian form as where λ is the dual variable associated with the equality constraint.

6. The method according to claim 1, characterized in that The optimality conditions of the expanded expression are obtained based on the Lagrangian form, including: The optimality condition for the expanded expression based on the Lagrangian form is Where Q = D H D, G = D H B, λ opt is the optimal dual variable, y opt is a candidate solution for the optimal transmit beam vector, ρ∈[0,1] is a weighting factor that weighs radar performance and communication performance, the superscript H represents the transpose operation, and F FD is the optimal unconstrained precoding matrix, C Rad represents the radar beam covariance matrix, P t Indicates the total transmit power of the dual-function base station, Indicates that the dimension is 2N T ×2N T The superscript T indicates the transpose.

7. The method according to claim 1, characterized in that The optimization problem of designing the analog precoder and the one-bit digital-to-analog converter quantized signal according to the optimal transmit beam vector includes: The optimization problem of designing the analog precoder according to the optimal transmit beam vector by fixing a one-bit digital-to-analog converter quantized signal is: Among them, F A represents the analog precoder, represents a one-bit DAC quantized signal, y opt is the candidate solution for the optimal transmit beam vector, N T represents the number of transmitting antennas of the dual-function base station, n represents the column subscript index of the element in the mth row and nth column in the analog precoder, and m represents the row subscript index of the element in the mth row and nth column in the analog precoder.

8. The method according to claim 7, characterized in that The method further comprises: The optimization problem of designing a one-bit digital-to-analog converter quantized signal according to the optimal transmit beam vector with a fixed analog precoder is: in, represents an element in a one-bit DAC quantized signal, represents a set of quantized values ​​of a one-bit digital-to-analog converter, represents the set {1,2,...,N RF }. The real part and the imaginary part are quantized by one bit respectively. The optimization problem of the above one-bit DAC quantized signal is further rewritten as: in, Represents the set {1,2,...,2N RF }, It represents the operator for taking the real part of a complex number. Represents the operator for taking the imaginary part of a complex number.

9. The method according to claim 1, characterized in that: The optimization problem is solved by using Matlab toolbox and BTSA-GA algorithm respectively, and an optimized analog precoder and a one-bit digital-to-analog converter quantization signal are obtained, including: The optimization problem of the analog precoder is solved by using Matlab toolbox Manopt to obtain an optimized analog precoder.

10. The method according to claim 9, characterized in that The method further comprises: The BTSA-GA algorithm is used to solve the optimization problem of the quantized signal of the one-bit digital-to-analog converter. In the initialization stage of the BTSA-GA algorithm, a batch of trees are generated as follows: T i,j =l j +r i,j (h j -l j ), Among them, l j is the lower bound of the search space -1, h j is the upper bound of the search space 1, r i,j It is a random number generated for each dimension at each position in the range [0,1]; To T i,j Binary discretization: Find the current best solution by: S=min f(T i ),i∈1,2,…,2N RF , Trees generate seeds, and the positions of trees and seeds in the n-dimensional search space are regarded as possible solutions to the optimization problem. The positions of seeds are obtained through two search equations: Among them, S i,j is the jth dimension of the i-th seed generated by the i-th tree, T i,j is the jth dimension of the i-th tree, B j is the jth dimension of the best tree position obtained so far, T r,j is the jth dimension of the rth tree randomly selected from the population, β i,j is a randomly generated scaling factor in the range [-1,1], γ is a random number generated for each position and dimension in the range (0,1), and ST is a selection control parameter; Set the threshold H, when T i The number of generations of generating seeds exceeds H and there is no better seed to update T i When T is updated through crossover and mutation operations in the genetic algorithm i , and obtain the global optimal solution, that is, the optimized one-bit DAC quantized signal.

Images