- Research
- Open access
- Published:
Synergistic optimization strategy for beamforming and power allocation in dual-functional radar-communication systems
EURASIP Journal on Advances in Signal Processing volume 2024, Article number: 89 (2024)
Abstract
The dual-functional radar-communication (DFRC) integrated system presents an ideal solution to address the challenge of spectrum resource congestion in future networks. This paper explores an adaptive power allocation technique based on beamforming to enhance the word error probability (WEP) performance of the DFRC system. Initially, a joint optimization model is developed to minimize the WEP while adhering to constraints on radar signal-to-interference-noise ratio (SINR), peak-to-average-power ratio, sidelobe level, and total transmit power. This model incorporates dual-function transmit beam, radar, and communication receive beam patterns. Subsequently, the proposed subproblem convex relaxation alternating update (SCRAU) algorithm is introduced to achieve a locally optimal solution for multi-carrier power allocation. This algorithm decomposes the original non-convex optimization problem into three sub-problems with lower complexity and iteratively optimizes them. Simulation experiments validate that the SCRAU algorithm can simultaneously fulfill radar and communication functions. The SCRAU algorithm demonstrates superior WEP performance compared to current advanced algorithms.
1 Introduction
The advancement of 5G and the anticipated development of 6G technologies are escalating the demand for wireless spectrum within the global communication industry. To tap into additional spectrum resources, communication spectrum is being extended to higher frequency bands, leading to overlaps with radar spectrum [1, 2]. The dual-functional radar-communication (DFRC) integrated system, which achieves both radar and communication functions concurrently by utilizing shared hardware and spectrum resources, is considered a paradigm for the next generation of wireless systems and networks [3,4,5].
Currently, the predominant dual-function implementation schemes involve fast-time modulation for information embedding (IE) and the sidelobe control information embedding (SCIE) scheme [6]. In the DFRC integrated system, the primary drawback of fast information embedding schemes lies in the potential for inefficient bandwidth utilization, increased system complexity, and adverse impacts on communication performance [7, 8]. Conversely, the SCIE scheme in the DFRC integrated system offers a significant advantage in reducing interference, particularly in mitigating the interference from the radar system on the communication component. This contributes to enhancing the performance of the communication system, ensuring the reliability and stability of the communication link [9].
The SCIE scheme comes in various forms, including amplitude modulation (AM) [10], phase modulation (PM) [11] and quadrature amplitude modulation (QAM) [6]. A general framework for information embedding in frequency-hopped (FH) multiple-input-multiple-output (MIMO) DFRC system is introduced in [10]. A unified formula is defined to accommodate various existing modulation methods. The proposed hybrid modulation strategy significantly increases the data rate, but the scheme still has obvious defects in range sidelobe suppression and target detection. In [11], an adaptive DFRC system based on PM iteratively updates target scene parameters by learning from the radar scene at the receiver and dynamically reallocating transmit power. The goal is to adapt to the time-varying characteristics of radar targets and surroundings. A two-part waveform optimization method is proposed that improves target detection without affecting communication function. The optimization problem of DFRC transmit beamforming and receive beamforming is studied in [12]. The alternating direction sequential relaxation programming (ADSRP) algorithm is proposed, optimizing the detection performance of DFRC while simultaneously achieving the dual functions. In the integrated radar and communication system discussed in [13], an orthogonal frequency division multiplexing (OFDM) integrated signal with QAM is employed. To address the challenge of high peak-to-average power ratio (PAPR), an active constellation expansion algorithm is applied to control the PAPR of the integrated signal. This algorithm enhances target detection performance with minimal impact on the word error probability (WEP).
WEP measures the proportion of incorrectly interpreted words at the communication receiving end. A lower WEP indicates fewer data misinterpretations, helping to reduce the need for redundant information, improve system efficiency, and ensure better communication quality under various conditions [12]. PAPR has significant implications for the integration of radar and communication systems. Higher PAPR may lead to nonlinear distortions in power amplifiers, impacting the system's dynamic range [14]. In DFRC integrated systems, optimizing PAPR handling is crucial for system performance, reliability, flexibility, and adaptability to various communication environments. Optimizing power allocation enables efficient system operation, enhances communication reliability, minimizes overall energy consumption, and ensures stable system performance in complex environments [15]. Meanwhile, precise control of beamforming can further minimize radiation in non-target directions, reducing additional energy consumption. By synergistically optimizing beamforming and power allocation, the DFRC system can achieve better signal transmission and target detection performance, and improve the overall effectiveness of the system.
Increasing transmit power can improve radar SINR and communication WEP, but blindly increasing transmit power may consume resources and increase the risk of interception [16]. For the current research on multi-carrier power allocation in DFRC systems, the problem of optimizing communication WEP while considering total transmit power, PAPR and SINR is still not sufficiently concerned. To address this problem, we propose a multi-carrier power allocation problem for DFRC signals based on beamforming with sidelobe control information embedded in it. To ensure the Quality of Server (QoS), we comprehensively investigate transmit/receive beamforming in DFRC integrated systems. Since DFRC systems share transmission signals and equipment, the radar and communication performance are highly coupled. Relying only on radar or communication channel states for power allocation can severely degrade performance. Therefore, limited power allocation while realizing radar and communication functions has become an urgent problem.
This paper makes the following key contributions:
-
(1)
Beamforming-Based Multi-Carrier Power Allocation: We introduce a beamforming-based approach to present a multi-carrier power allocation problem for integrated radar-communication signals. The objective is to minimize the word error probability (WEP) while considering constraints on radar signal-to-interference-noise ratio (SINR), sidelobe level (SLL), peak-to-average-power ratio (PAPR), and total transmit power.
-
(2)
SCRAU Algorithm for Non-convex Optimization: To address the non-convex optimization problem, we propose the subproblem convex relaxation alternating update (SCRAU) algorithm. This algorithm decomposes the original non-convex problem into three sub-problems, achieving a locally optimal solution through alternating iterative optimization.
The paper is structured in the following manner: Section II introduces the system model and signal model; Section III elucidates the proposed algorithm; Section IV provides simulation results and analysis, and Section V presents the conclusions.
Notations: \(\left( \cdot \right)^{T}\),\(\left( \cdot \right)^{\dag }\),\(\left( \cdot \right)^{*}\),\(\left| \cdot \right|\),\(\left\| \cdot \right\|\),\(E\left\{ \cdot \right\}\),\(C\),\(C^{M \times 1}\),\(C^{M \times M}\) and \(I^{M \times M}\) denote transpose, conjugate transpose, conjugate, absolute value, the Euclidean norm, average value, complex ensemble, the set of \(M \times 1\) vectors with complex entries, the set of \(M \times M\) matrices with complex entries, and \(M \times M\) identity matrix, respectively. \(\mathcal{C}\mathcal{N}\left( { \cdot { ,} \cdot } \right)\) denotes the complex Gaussian distribution. \({\text{Re}} \left\{ \cdot \right\}\) denotes the real part of the number in parentheses.
2 System model and signal model
Figure 1 illustrates the DFRC integrated system, which consists of dual-function transmitting co-located array, communication receiving array and radar receiving array. The dual-function transmitting co-located array contains \(M_{T}\) array elements, and the radar receiving array contains \(M_{R,r}\) array elements, and communication receiving array contains \(M_{R,c}\) array elements. The integrated signals transmitted by the dual-function transmitting array can realize both radar and communication functions. We consider the scenario where all arrays are assumed to be uniform linear arrays, and the spacing between array elements is set to half the wavelength. The radar and communication receiving arrays have the same number of array elements, that is \(M_{R,r} = M_{R,c} = M_{T}\). In this model, the radar transmitting and the receiving arrays are positioned in close proximity. Looking from the viewpoint of the target in Fig. 1, the DFRC transmitting array and the radar receiving array are assumed to be positioned at the same angle. Assuming that we have information about each angle and the channel state information (CSI) for both the radar and communication [12].
\(\left\{ {s_{n} } \right\}_{n = 1}^{N}\) is defined as an integrated multi-carrier signal with \(N\) subcarriers. Let \(\Delta f\) be the frequency interval of the subcarriers and assume that \(\Delta f\) is large enough, then the subcarriers are independent of each other. The power of the \(n{ - }th\) subcarrier is \(E\left\{ {\left| {s_{n} } \right|^{2} } \right\} = P_{n}\). In the DFRC integrated system, the radar and communication functions are implemented by the \(s_{n}\) of the dual-function transmitter. The main lobe of the beam pointing to the detection target is used to detect the information of the target, and the sidelobe of the beam pointing to the communication user transmits the communication symbol. Radar performance is affected by the source of interference, and from the perspective of the communication, the reflected signal from both the target and the interference source can affect the communication performance.
\({\varvec{q}}_{n} = \left[ {q_{1,n} ,q_{2,,n} ,...,q_{{M_{T} ,n}} } \right]^{\dag } \in C^{{M_{T} \times \, 1}}\) is defined as the transmit beamforming vector of \(M_{T} \times 1\) on the \(n{ - }th\) subcarrier, where \(\left\| {{\varvec{q}}_{n} } \right\| = 1\). At the dual-function transmitter, \(s_{n}\) and \({\varvec{q}}_{n}\) are multiplied to obtain the integrated transmission signal \({\varvec{q}}_{n} s_{n}\) of \(M_{T}\) elements on the \(n{ - }th\) subcarrier. Target detection is achieved through the analysis of echo signals received by the radar receiving array. Simultaneously, the communication receiving array processes the received signals to facilitate the transmission of information.
\({\varvec{v}}_{n} = \left[ {v_{1,n} ,v_{2,n} ...,v_{{M_{R} ,n}} } \right]^{\dag } \in C^{{M_{R} \times 1}}\) is defined as the radar receive beamforming vector on the \(n{ - }th\) subcarrier which is of size \(M_{R} \times 1\). Assume that the channel between the dual-function transmitting co-located array and the target (the interference source) is \(H \in C^{{M_{T} \times M_{T} }}\)(\(H_{I} \in C^{{M_{T} \times M_{T} }}\)), respectively. The target and the interference source are located at angles \(\omega_{0}\) and \(\omega_{1}\), respectively. Then the received signal on the \(n{ - }th\) subcarrier at the radar receiving array is
where \({\varvec{H}}_{n} = a_{0,n} {\varvec{\beta}}(\omega_{0} ){\varvec{\alpha}}(\omega_{0} )^{T}\) and \({\varvec{H}}_{I,n} = a_{1,n} {\varvec{\beta}}(\omega_{1} ){\varvec{\alpha}}(\omega_{1} )^{T}\).\({\varvec{\alpha}}(\omega ) \in C^{{M_{T} }}\) and \({\varvec{\beta}}(\omega ) \in C^{{M_{R} }}\) are the transmitting and receiving steering vectors pointing to the angle \(\omega\), respectively. \(a_{0,n} \in C\) and \(a_{1,n} \in C\) are the channel coefficients on the corresponding \(n{ - }th\) subcarrier and obey complex Gaussian distributions \(\mathcal{C}\mathcal{N}\left( {0,\,\,\sigma_{a,0,n}^{2} } \right)\) and \(\mathcal{C}\mathcal{N}\left( {0,\,\,\sigma_{a,1,n}^{2} } \right)\), respectively. \(\varpi_{r,n} \in C^{{M_{R} \times 1}}\) represents the noise vector of the radar on the \(n{ - }th\) subcarrier and obeys Gaussian distribution \(\mathcal{C}\mathcal{N}\left( {0,\sigma_{r,n}^{2} {\varvec{I}}_{{M_{R} \times M_{R} }} } \right)\).
Assume that \({\varvec{u}}_{n} = {\varvec{q}}_{n} s_{n}\), the received signal \(y_{r,n}\) is filtered by \({\varvec{v}}_{n}\) to improve the radar SINR of the \(n{ - }th\) subcarrier at the radar receiving array.
where \({{\varvec{\uppsi}}}_{n} = E\left\{ {{\varvec{H}}_{I,n} {\varvec{u}}_{n} {\varvec{u}}_{n}^{\dag } {\varvec{H}}_{I,n}^{\dag } } \right\}\).
Assume the communication receiving array is positioned at azimuth angle \(\partial\), the baseband communication transmission signal on the \(n{ - }th\) subcarrier is
where each value of \({\varvec{d}}_{n}\) corresponds to only one particular communication symbol. In order to transmit \(\log_{2} \left( L \right)\)-bit symbols (one pulse duration) on each subcarrier, \(L\) \({\varvec{d}}_{n}\) different transmit beams need to be designed and stored in the code word dictionary \({\mathcal{D}} = \left\{ {{\varvec{d}}_{1,n} ,...,{\varvec{d}}_{L,n} } \right\}\). The transmission of communication symbols is achieved by controlling the amplitude of the expression (3), i.e., the sidelobe level of the transmitted beam in the direction of the azimuth angle. Therefore, during the duration of the pulse, a \({\varvec{d}}_{l,n}\) is selected from \({\mathcal{D}} = \left\{ {{\varvec{d}}_{1,n} ,...,{\varvec{d}}_{L,n} } \right\}\). At the communication receiving end, the corresponding symbol is determined by judging the \({\varvec{d}}_{l,n}\), and then the symbol transmission is completed.
\({\varvec{w}}_{n} = \left[ {w_{1,n} ,w_{2,n} ,...,w_{{M_{R} ,n}} } \right]^{\dag } \in C^{{M_{R} \times 1}}\) is defined as the communication receive beamforming vector on the \(n{ - }th\) subcarrier, which is of size \(M_{R} \times 1\).The received signal at the communication receiving array comprises both the desired signal and interference signal. The desired signal comes from the dual-function transmitting array. The interference signal comes from the reflection of the target or from the reflection of the interference source. The received signal on the \(n{ - }th\) subcarrier at the communication receiving array is
where \(\partial^{ * }\),\(\vartheta_{0}\) and \(\vartheta_{1}\) are the angles of arrival of the transmitted signal (reception angle of the communication array), the angle of the desired communication signal and the angle of the communication interference signal, respectively. \(\gamma_{n}\) is the communication channel coefficient on the \(n{ - }th\) subcarrier and obeys \(\mathcal{C}\mathcal{N}\left( {0,\,\,\sigma_{\gamma ,n}^{2} } \right)\). \(b_{0,n} \in C\)(\(b_{1,n} \in C\)) is the channel coefficient of the \(n{ - }th\) subcarrier from the transmitter to the reflectors (target and interference source) and then to the communication user.\(b_{0,n} \in C\) (\(b_{1,n} \in C\)) obeys \(\mathcal{C}\mathcal{N}\left( {0,\,\,\sigma_{b,0,n}^{2} } \right)\)(\(\mathcal{C}\mathcal{N}\left( {0,\,\,\sigma_{b,1,n}^{2} } \right)\)). \(\varpi_{c,n} \in C^{{M_{R} \times 1}}\) is the communication noise vector on the \(n{ - }th\) subcarrier and obeys Gaussian distribution \(\mathcal{C}\mathcal{N}\left( {0,\,\,\sigma_{c,n}^{2} {\varvec{I}}_{{M_{R} \times M_{R} }} } \right)\).
The received signal \(y_{c,n}\) is filtered by \({\varvec{w}}_{n}\) to improve the communication SINR on the \(n - th\) subcarrier at the communication receiving array.
where \({\mathbf{\mathcal{M}}}_{n} = E\left\{ {{\varvec{f}}_{n} {\varvec{u}}_{n} {\varvec{u}}_{n}^{\dag } {\varvec{f}}_{n}^{\dag } } \right\}\) and \({\varvec{f}}_{n} = a_{0,k} {\varvec{\beta}}(\vartheta_{0} ){\varvec{\alpha}}(\omega_{0} )^{T} + a_{1,k} {\varvec{\beta}}(\vartheta_{1} ){\varvec{\alpha}}(\omega_{1} )^{T}\).
The amplitude of the transmitted sidelobes for each pulse is selected from a pre-designed codebook \({\varvec{d}}_{l} = \left( {d_{l,1} ,...,d_{l,N} } \right)^{T} \in C^{N}\). This codebook consists of \(L\) code words of length \(N\), i.e., \({\mathcal{D}} = \left\{ {{\varvec{d}}_{1} ,...,{\varvec{d}}_{L} } \right\}\). The identification of the communication symbol for the communication user is based on the comparison of the amplitude of the received signal \(\left| {y_{c,n} } \right|\) on the \(n{ - }th\) subcarrier with \(L - 1\) thresholds. We assume that \({\varvec{d}}_{\ell }\) is the communication transmit symbol, during the pulse duration, the expression for the communication WEP is [17].
where \({\varvec{d}}_{l,n}\) is a code word different from \({\varvec{d}}_{s,n}\).
3 Multi-carrier power allocation algorithm for integrated signals
The objective is to minimize WEP in a DFRC integrated system with beamforming, considering constraints on radar SINR, PAPR, and total transmit power. The multi-carrier power allocation scheme \({\varvec{P}} = [\left\| {{\varvec{u}}_{1} } \right\|^{2} ,\left\| {{\varvec{u}}_{2} } \right\|^{2} ,...,\left\| {{\varvec{u}}_{N} } \right\|^{2} ]^{T}\) can be derived through the optimization problem presented as follows [18]
where C1 indicates that the sum of radar SINR should not exceed \(\Re\); C2 indicates that the transmit power cannot exceed \(P_{tot}\); C3 controls the amplitude of the side flap pointing in azimuth to accomplish communication symbol transmission; C4 indicates that the PAPR cannot exceed \(\delta_{p}\). In addition, it can be seen from constraint C3 that in order to transmit \(\log_{2} L\) bit symbols/pulses, optimization problem (7) needs to be solved \(L\) times. The PAPR is denoted as [19]
Since the optimization problem (7) is high-dimensional and non-convex, it is hard to solve the closed-form solution directly. Drawing inspiration from the algorithm introduced in [20, 21], the divisible structure of the objective function allows for its expression as follows:
To obtain a suboptimal solution of problem (9), i.e., \(\left\{ {{\varvec{u}}_{n}^{*} ,{\varvec{v}}_{n}^{*} ,{\varvec{w}}_{n}^{*} } \right\}\), an iterative SCRAU algorithm is introduced, aiming to tackle the three low-complexity sub-problems. Let \(\left\{ {{\varvec{u}}_{n}^{t} ,{\varvec{v}}_{n}^{t} ,{\varvec{w}}_{n}^{t} } \right\}\) denote the solution in the \(t{ - }th\) iteration, \(t \ge 1\). Firstly, a feasible initial value \(\left\{ {{\varvec{u}}_{n}^{0} ,{\varvec{v}}_{n}^{0} ,{\varvec{w}}_{n}^{0} } \right\}\) is set and then iteratively solve sub-problems \({\varvec{u}}_{n}\),\({\varvec{v}}_{n}\) and \({\varvec{w}}_{n}\). To attain the solution by the \(\left( {t + 1} \right){ - }th\) iterative step, the following iterative optimization procedure is given:
-
(1)
update \({\varvec{u}}_{n}^{t + 1}\) by fixing \({\varvec{v}}_{n}^{t}\) and \({\varvec{w}}_{n}^{t}\);
-
(2)
update \({\varvec{v}}_{n}^{t + 1}\) by fixing \({\varvec{u}}_{n}^{t + 1}\) and \({\varvec{w}}_{n}^{t}\);
-
(3)
update \({\varvec{w}}_{n}^{t + 1}\) by fixing \({\varvec{u}}_{n}^{t + 1}\) and \({\varvec{v}}_{n}^{t + 1}\).
The iteration process stops when the difference between the optimization results in the \((t + 1){ - }th\) iteration and the \(t{ - }th\) iteration is not more than \(\varepsilon\).
where \(W_{l}^{t} \left( {{\varvec{u}}_{n} ,{\varvec{v}}_{n} ,{\varvec{w}}_{n} } \right)\) denotes the minimum value of WEP obtained in the \(t{ - }th\) iteration, and \(\varepsilon\) is a predefined constant used as the stopping criterion for iterations.
A. Radar Transmit Beamformer Update
When \({\varvec{v}}_{n}^{t}\) and \({\varvec{w}}_{n}^{t}\) are fixed to solve out \({\varvec{u}}_{n}^{t + 1}\), (7a) takes the following shape
Therefore, when \({\varvec{v}}_{n}^{t}\) and \({\varvec{w}}_{n}^{t}\) are fixed, \({\varvec{u}}_{n}^{t + 1}\) can be obtained by solving \({\mathbf{\mathcal{P}}}_{1}\).
Because of the non-convexity of C1 and C4, the optimization problem (12) is non-convex. We can transform problem (12) into a relaxed convex problem and subsequently solve it.
1) Relaxation of the objective function in (12)
Since \(SINR_{c,n}\) is non-convex in the objective function (12), we perform a convex transformation on \(SINR_{c,n}\). Firstly, \(SINR_{c,n}\) can be rewritten as
where \({\varvec{\zeta}}_{n}^{t} = \sigma_{r,n}^{2} \left| {\user2{(w}_{n}^{t} )^{\dag } {\varvec{\beta}}\left( {\partial^{ * } } \right)} \right|^{2}\), \(\xi_{n}^{t} = \sigma_{c,n}^{2} \left( {{\varvec{w}}_{n}^{t} } \right)^{\dag } {\varvec{w}}_{n}^{t}\) and \(\user2{\mathcal{G}}_{n}^{t} = {\varvec{f}}_{n}^{t} {\varvec{w}}_{n}^{t} \left( {{\varvec{w}}_{n}^{t} } \right)^{\dag } \left( {{\varvec{f}}_{n}^{t} } \right)^{\dag }\). In accordance with the linear approximation from the first-order Taylor series expansion [22]
where the equation holds true when \(x = x^{0}\). Subsequently, according to (14), it can be derived that
where \({\varvec{u}}_{n}^{t}\) denotes the transmit beam formation vector obtained in the previous iteration. The objective function (11) is equivalently replaced by (15), i.e.,
where \(\eta_{l} = \left| {{\varvec{d}}_{l,n} - {\varvec{d}}_{s,n} } \right|^{2}\).
2) Relaxation of C1 in (12)
To solve out \({\varvec{u}}_{n}^{t + 1}\) by fixing \({\varvec{v}}_{n}^{t}\) and \({\varvec{w}}_{n}^{t}\), we rewrite C1 as
where \({\mathbf{\mathcal{T}}}_{n} = E\left\{ {{\varvec{H}}_{I,n} {\varvec{v}}_{n}^{t} \left( {{\varvec{v}}_{n}^{t} } \right)^{\dag } {\varvec{H}}_{I,n}^{T} } \right\}\). The following Eq. (18) is obtained by shifting the terms left and right on both sides of the inequality.
3) Relaxation of C4 in (12)
From [19], constraint C4 can be recast as
where
We use the first-order expansion condition to obtain
According to (19) and (21), the constraint C4 can be equivalently replaced as
Problem (12) can be formulated as the relaxation problem presented below
The convexity of the optimization problem (23) enables the acquisition of its optimal solution through the utilization of the CVX toolbox [23].
B. Radar Receive Beamformer Update
When \({\varvec{u}}_{n}^{t + 1}\) and \({\varvec{w}}_{n}^{t}\) are fixed to solve out \({\varvec{v}}_{n}^{t + 1}\), \({\varvec{v}}_{n}^{t + 1}\) can be obtained by (24)
The optimization problem (24) can be equated to the minimum variance distortionless response (MVDR) beamforming problem [24]:
Therefore, the local analytic solution of \({\varvec{v}}_{n}^{t + 1}\) is
C. Communication Receive Beamformer Update
When \({\varvec{u}}_{n}^{t + 1}\) and \({\varvec{v}}_{n}^{t + 1}\) are fixed to solve out \({\varvec{w}}_{n}^{t + 1}\), \({\varvec{w}}_{n}^{t + 1}\) can be obtained by (27)
The optimization problem (27) can be equated to the MVDR beamforming problem:
Therefore, the local analytic solution of \({\varvec{w}}_{n}^{t + 1}\) is
D. Algorithm Complexity and Convergence Analysis
1) Algorithm Complexity Analysis
The complexity of the proposed SCRAU is related to \(N\), the dimension of \({\varvec{u}}_{n}\),\({\varvec{v}}_{n}\),\({\varvec{w}}_{n}\) and the number of iterations. The relaxed convex problem (23) is solved using the interior point method (IPM) (via the CVX toolbox) during the \(t{ - }th\) iteration to get \({\varvec{u}}_{n}^{t}\). Its computational complexity is \(O\left( {NM_{T}^{3.5} } \right)\)[25, 26]. \({\varvec{v}}_{n}^{t}\) and \({\varvec{w}}_{n}^{t}\) are obtained from problems (25) and (28), respectively, and the computational complexity of both is \(O\left( {NM_{R}^{3} } \right)\). Therefore, for each iteration, the computational complexity of SCRAU is \(O\left( {NM_{T}^{3.5} + 2NM_{R}^{3} } \right)\).
2) Convergence Analysis
During the \(t{ - }th\) iteration, the WEP is non-increasing when optimizing \({\varvec{u}}_{n}\), i.e., \(W_{l} \left( {\left. {{\varvec{u}}_{n}^{t + 1} } \right|{\varvec{v}}_{n}^{t} ,{\varvec{w}}_{n}^{t} } \right) \le W_{l} \left( {\left. {{\varvec{u}}_{n}^{t} } \right|{\varvec{v}}_{n}^{t} ,{\varvec{w}}_{n}^{t} } \right)\); when optimizing \({\varvec{v}}_{n}\), the WEP is non-increasing, i.e., \(W_{l} \left( {\left. {{\varvec{u}}_{n}^{t} } \right|{\varvec{v}}_{n}^{t} ,{\varvec{w}}_{n}^{t} } \right) \le W_{l} \left( {\left. {{\varvec{v}}_{n}^{t + 1} } \right|{\varvec{u}}_{n}^{t + 1} ,{\varvec{w}}_{n}^{t} } \right)\); since \({\varvec{w}}_{n}\) is not an independent variable of WEP when \({\varvec{u}}_{n}\) and \({\varvec{v}}_{n}\) are fixed, then \(W_{l} \left( {\left. {{\varvec{w}}_{n}^{t + 1} } \right|{\varvec{u}}_{n}^{t + 1} ,{\varvec{v}}_{n}^{t + 1} } \right) = W_{l} \left( {\left. {{\varvec{v}}_{n}^{t + 1} } \right|{\varvec{u}}_{n}^{t + 1} ,{\varvec{w}}_{n}^{t} } \right)\). Thus, it can be obtained that \(W_{l} \left( {{\varvec{u}}_{n} ,{\varvec{v}}_{n} ,{\varvec{w}}_{n} } \right)\) is non-increasing during the iteration process, i.e., \(W_{l} \left( {{\varvec{u}}_{n}^{t + 1} ,{\varvec{v}}_{n}^{t + 1} ,{\varvec{w}}_{n}^{t + 1} } \right) \le W_{l} \left( {{\varvec{u}}_{n}^{t} ,{\varvec{v}}_{n}^{t} ,{\varvec{w}}_{n}^{t} } \right)\).
Figure 2 illustrates the relationship between WEP and the number of iterations for the SCRAU algorithm under various radar SINR values and transmit power conditions. In Fig. 2, with an increasing number of iterations, WEP stabilizes at a fixed value in all three cases. According to the monotonic convergence theorem, the SCRAU algorithm proposed in this paper is convergent [27].
4 Simulation results and analysis
The efficacy of the SCRAU algorithm is verified through simulation in this section. The simulation scenario is shown in Fig. 3. Assume that there is an interference source, the target and the interference source are located at angles \(\omega_{0} = 0^{ \circ }\) and \(\omega_{1} = 60^{ \circ }\), respectively. Assume that the number of the total subcarriers is \(N = 16\), the numbers of transmitting and receiving elements are \(M_{T} = M_{R,c} = M_{R,r} = 16\). The communication receiving array is located at angle \(\partial = 40^{ \circ }\), and \(L = 2\) types of sidelobe are pre-designed, i.e., −20 dB and −25 dB. From the perspective of the communication receiving array, it receives desired signals from the communication arrival angle \(\partial^{ * } = 0^{ \circ }\) and interference signals from angles \(\vartheta_{0} = 50^{ \circ }\) and \(\vartheta_{1} = 70^{ \circ }\). Without loss of generality, we assume a noise variance of 1 for both radar and communication, i.e., \(\sigma_{r,n}^{2} = \sigma_{c,n}^{2} = 1\). All other channel coefficients are set to 0.01, i.e., \(\sigma_{a,0,n}^{2} = \sigma_{a,j,n}^{2} = \sigma_{b,0,n}^{2} = \sigma_{b,1,n}^{2} = 0.01\), meanwhile, the variance of the communication channel coefficients is set to 1, i.e.,\(\sigma_{\gamma ,n}^{2} = 1\). Assume sufficient frequency spacing between shared subcarriers for independence of observation noise and signal reflections. The iteration stopping constant is set as \(\varepsilon = 10^{ - 5}\).
To better understand the intensity distribution of the transmitted signal in different directions, Fig. 4 shows the relationship between the normalized transmit power distribution patterns and spatial angles for \(L = 2\) transmit beam patterns. Figure 4 shows that the two transmit beams have the same amplitude of the main lobe, this implies that these two beams have comparable detection ranges. The angle of the communication receiving array is \(\partial = 40^{ \circ }\), and the communication SLLs are −20 dB and −25 dB, respectively. During the pulse duration, the communication receiving array determines the corresponding communication symbol by comparing the amplitude of received signal from the side-lope during the pulse duration. In addition, the presence of a deep null at azimuth \(\omega_{1} = 60^{ \circ }\), this means that the system can effectively suppress or minimize interference from that direction, thus improving the system's performance and immunity to interference (Table 1).
Figure 5 shows the communication receiving beam when the main lobe direction is \(\partial^{ * } = 0^{ \circ }\), and the interference source direction is located in \(\vartheta_{0} = 50^{ \circ }\) and \(\vartheta_{1} = 70^{ \circ }\). From Fig. 5, it can be seen that there are deep nulls in both azimuths \(\vartheta_{0} = 50^{ \circ }\) and \(\vartheta_{1} = 70^{ \circ }\). This shows that the proposed algorithm can reduce the beamforming of interference, so that it has less impact on communication. The direction diagram of the communication receive beam is mainly determined by the azimuths \(\partial^{ * }\),\(\vartheta_{0}\) and \(\vartheta_{1}\).
In order to further delve into the power allocation algorithm of beamforming-based DFRC integrated system, we investigate the adaptive power allocation of all \(N\) subcarriers under specific channel conditions. A traditional communication system (TCS) without considering the radar SINR constraint and WEP minimization is used as a comparison simulation to verify the performance of the DFRC integrated system. Assuming the state of the communication channel is predetermined. The normalized communication signal-to-noise ratio (SNR) [12, 28] is as shown in Fig. 6a, i.e., \(10\log_{10} \left( {\sigma_{\gamma ,n}^{2} /\sigma_{c,n}^{2} } \right)\), where it can be seen that subcarriers 2 and 9 are favorable to the communication system. Figure 6b shows the channel strength of the desired radar signal, i.e., \(10\log_{10} \left( {\sigma_{a,0,n}^{2} /\sigma_{r,n}^{2} } \right)\); Fig. 6c shows the channel strength of the radar interference signal, i.e., \(10\log_{10} \left( {\sigma_{a,j,n}^{2} /\sigma_{r,n}^{2} } \right)\), where subcarriers 3, 4 and 5 interfere more strongly with the radar; Fig. 6d shows the channel strength of the communication interference signal, i.e., \(10\log_{10} \left( {\sigma_{b,1,n}^{2} /\sigma_{c,n}^{2} } \right)\).
Figure 7 illustrates the multi-carrier waveform power allocation for the DFRC integrated system under TCS, SINR = 5, and SINR = 10 for the specific channel states described in Fig. 6. When \(\delta_{p} = 2.5\), the power allocation results of \(P_{tot} = 10kw\) and \(P_{tot} = 30kw\) are shown in Fig. 7a and b, respectively. It is evident from Fig. 7a and b that, in the case of TCS, the highest power allocation is assigned to subcarrier 2. This is because the TCS only considers the communication channel state and does not consider radar SINR constraints. With radar SINR constraints, the DFRC integrated system considers both radar and communication system channel conditions. When the SINR = 10, the subcarrier 9 is allocated with the most power because the subcarrier 9 has the worst radar channel state and communication channel is better, thus more power is needed to meet the constraint of SINR = 10. In addition, higher radar SINR requirements result in a more 'uniform' power allocation across subcarriers.
In general, under low radar SINR, the allocation of power to subcarriers is predominantly influenced by the communication channel state. Conversely, under high radar SINR constrains, the subcarrier power allocation is more significantly impacted by the radar channel state. By comparing Fig. 7a and b, it is clear that increasing total transmitted power leads to a reduction in power allocation discrepancies between the three scenarios. Analysis of Fig. 7a and c shows that when the PAPR constraint is lower, the gap in the distribution of the transmit power over each subcarrier is smaller. In addition, from the deep nulls in Figs. 4 and 5, it is evident that the interference channel states of radar and communication have almost negligible effects on the power allocation.
The proposed SCRAU algorithm is founded upon the sidelobe control information embedding scheme. To validate the efficacy of the SCRAU algorithm, this study compares its performance with that of the QAM-IE algorithm. The QAM-IE algorithm solution process is detailed in [29]. The fundamental principle of the QAM-IE algorithm is to control the amplitude and phase of the sidelobe directed toward the communication user while preserving the main lobe gain of the radar. The communication receiver demodulates the received signals to determine the corresponding communication symbols.
WEP performance of SCRAU and QAM-IE for \(\delta_{p} = 2.5\) with TCS, radar SINR = 5 and 10 is shown in Fig. 8. Figure 8 show that TCS achieves the lowest WEP in both SCRAU and QAM-IE algorithms. TCS, not considering radar performance, allocates transmit power based on communication channel conditions to minimize WEP. From Fig. 8, it is evident that an increase in SINR corresponds to an elevated value of WEP, i.e., the performance of communication accuracy decreases. Increasing the total transmit power leads to a gradual decrease in WEP under the three conditions. This indicates that a judicious increase in transmit power can effectively reduce WEP, enhancing the reliability of communication transmission. Figure 8 show that the SCRAU algorithm, with the same and transmit power and radar SINR, outperforms the QAM-IE algorithm in combating noise, resulting in a significantly improved WEP performance.
The WEP performance of SCRAU and QAM-IE for radar SINR = 10 with \(\delta_{p} = 2.5\) and \(\delta_{p} = 2.1\) is visualized in Fig. 9. It is clear from Fig. 9 that the increase in total transmitted power is associated with a decrease in WEP. When \(\delta_{p} = 2.5\), the performance of WEP is better. This is due to the fact that the over-restriction on PAPR allows for a more uniform allocation of subcarrier power in specific channel states as shown in Fig. 7. This restriction prevents the DFRC system from adequately considering and accurately achieving the goal of minimizing the communication WEP, which leads to an increase in the WEP. In all cases, the WEP values obtained by the SCRAU algorithm are lower than those of the QAM-IE algorithm. The reason for this difference is that the SCRAU algorithm performs adaptive power allocation based on the channel state conditions with the aim of minimizing the communication WEP while meeting the radar performance constraints. The SCRAU algorithm improves the reliability of the communication transmission and ultimately optimizes the overall performance of the DFRC system. On the communication side, reducing the WEP helps to ensure the integrity and accuracy of data transmission. On the other hand, the QAM-IE algorithm uses a single carrier waveform and does not consider the radar communication coupling of the DFRC system.
The realization of this result can be attributed to several key factors: (1) the proposed method, by modulating in specific directions, mitigates interference from other directions, thereby enhancing the system's performance in complex environments; (2) the method in this paper results in higher achievable rates as well as larger distances between symbols in the symbol space on the subcarriers, which reduces the communication word error rate; (3) the modulation method in this paper is easier to integrate with the radar system for radar and communication integration; (4) QAM-IE optimizes only the transmit beamforming vectors while higher gain is obtained by optimizing the transmit and receive beamforming vectors; (5) the SCRAU algorithm adaptively allocates power according to the specific channel state.
5 Conclusion
This paper thoroughly investigates the adaptive power allocation problem in dual-functional radar-communication (DFRC) systems. We formulate the optimization criterion as a minimization of word error probability (WEP) under constraints of radar signal-to-interference-noise ratio (SINR), total transmit power, peak-to-average power ratio (PAPR), and sidelobe information. To address the non-convex nature of this problem, we introduce the subproblem convex relaxation alternating update (SCRAU) algorithm, which decomposes the original problem into three sub-problems: (1) dual-function transmit beamforming vector optimization subproblem; (2) design of the radar receive beamforming vector; and (3) design of the communication receive beamforming vector. We prove the decrease and boundedness of WEP, ensuring the convergence of the algorithm. Finally, we compare the SCRAU algorithm with the quadrature amplitude modulation information embedding (QAM-IE) algorithm under various radar SINR, total transmit power, and PAPR conditions. Numerical results demonstrate that the SCRAU algorithm effectively improves WEP performance, enabling adaptive power allocation by jointly considering radar and communication channel states.
A future research direction to consider: in our simulations and analyses, the power limit is set to be fixed within a single pulse duration. Future work will further explore optimization strategies for power allocation under different PRI conditions to ensure system performance under various operating conditions.
Availability of data and materials
Not applicable.
References
L. Zheng, M. Lops, Y.C. Eldar, X. Wang, Radar and communication coexistence: an overview: a review of recent methods. IEEE Signal Process. Mag. 36, 85–99 (2019). https://doi.org/10.1109/MSP.2019.2907329
S.K. Agrawal, A. Samant, S.K. Yadav, Spectrum sensing in cognitive radio networks and metacognition for dynamic spectrum sharing between radar and communication system: a review. Phys. Commun. 52, 101673 (2020). https://doi.org/10.1016/j.phycom.2022.101673
C. Shi, Y. Wang, F. Wang, S. Salous, J. Zhou, Joint optimization scheme for subcarrier selection and power allocation in multicarrier dual-function radar-communication system. IEEE Syst. J. 15, 947–958 (2020). https://doi.org/10.1109/JSYST.2020.2984637
L. Chen, Z. Wang, Y. Du, Y. Chen, F.R. Yu, Generalized transceiver beamforming for DFRC with MIMO radar and MU-MIMO communication. IEEE J. Sel. Areas Commun. 40, 1795–1808 (2022). https://doi.org/10.1109/JSAC.2022.3155515
J. Zeng, B. Liao, Transmit and receive hybrid beamforming design for OFDM dual-function radar-communication systems. EURASIP J. Adv. Signal Process. 2023, 37 (2023). https://doi.org/10.1186/s13634-023-00996-2
T. Tian, G. Li, H. Deng, J. Lu, Adaptive bit/power allocation with beamforming for dual-function radar-communication. IEEE Wirel. Commun. Lett. 11, 1186–1190 (2022). https://doi.org/10.1109/LWC.2022.3160674
A. Hassanien, M.G. Amin, E. Aboutanios, B. Himed, Dual-function radar communication systems: a solution to the spectrum congestion problem. IEEE Signal Process. Mag. 36, 115–126 (2019). https://doi.org/10.1109/MSP.2019.2900571
A. Hassanien, M.G. Amin, Y.D. Zhang, F. Ahmad, Signaling strategies for dual-function radar communications: an overview. IEEE Aerosp. Electron. Syst. Mag. 31, 36–45 (2016). https://doi.org/10.1109/MAES.2016.150225
A. Hassanien, M.G. Amin, Y.D. Zhang, F. Ahmad, Dual-function radar-communications: information embedding using sidelobe control and waveform diversity. IEEE Trans. Signal Process. 64, 2168–2181 (2016). https://doi.org/10.1109/TSP.2015.2505667
W. Baxter, E. Aboutanios, A. Hassanien, Joint radar and communications for frequency-hopped MIMO systems. IEEE Trans. Signal Process. 70, 729–742 (2022). https://doi.org/10.1109/TSP.2022.3142909
Y. Yao, P. Miao, Z. Chen, Cognitive waveform optimization for phase-modulation-based joint radar-communications system. IEEE Access. 8, 33276–33288 (2020). https://doi.org/10.1109/ACCESS.2020.2974787
T. Tian, T. Zhang, L. Kong, Y. Deng, Transmit/receive beamforming for MIMO-OFDM based dual-function radar and communication. IEEE Trans. Veh. Technol. 70, 4693–4708 (2021). https://doi.org/10.1109/TVT.2021.3072094
Y. Yang, J. Mei, D. Hu et al., Research on reducing PAPR of QAM-OFDM radar-communication integration sharing signal. J Eng 2019, 8042–8046 (2019). https://doi.org/10.1049/joe.2019.0726
Y. Huang, S. Hu, S. Ma, Z. Liu, M. Xiao, Designing low-PAPR waveform for OFDM-based RadCom systems. IEEE Trans. Wireless Commun. 21, 6979–6993 (2022). https://doi.org/10.1109/TWC.2022.3153606
Z. Li, J. Xie, H. Zhang, J. Ge, C. Qi, A robust power allocation strategy based on benefit–cost ratio for multiple target guidance in the C-MIMO radar system under blanket jamming. EURASIP J. Adv. Signal Process. 2022, 48 (2022). https://doi.org/10.1186/s13634-022-00880-5
C. Shi, S. Salous, F. Wang, J. Zhou, Low probability of intercept-based adaptive radar waveform optimization in signal-dependent clutter for joint radar and cellular communication systems. EURASIP J. Adv. Signal Process. 2016, 111 (2016). https://doi.org/10.1186/s13634-016-0411-6
J.G. Proakis, Digital communications, 4th edn. (McGraw Hill, New York USA, 2001), p.444
A.D. Maio, Y. Huang, M. Piezzo, S. Zhang, A. Farina, Design of optimized radar codes with a peak to average power ratio constraint. IEEE Trans. Signal Process. 59, 2683–2697 (2011). https://doi.org/10.1109/TSP.2011.2128313
Z. Cheng, Z. He, B. Liao, M. Fang, MIMO radar waveform design with PAPR and similarity constraints. IEEE Trans. Signal Process. 66, 968–981 (2017). https://doi.org/10.1109/TSP.2017.2780052
A. Tkacenko, P.P. Vaidyanathan, Iterative greedy algorithm for solving the FIR paraunitary approximation problem. IEEE Trans. Signal Process. 54, 146–160 (2005). https://doi.org/10.1109/TSP.2005.861054
M.N. Omidvar, M. Yang, Y. Mei, X. Li, X. Yao, DG2: a faster and more accurate differential grouping for large-scale black-box optimization. IEEE Trans. Evol. Comput. 21, 929–942 (2017). https://doi.org/10.1109/TEVC.2017.2694221
Y. Sun, P. Babu, D.P. Palomar, Majorization-minimization algorithms in signal processing, communications, and machine learning. IEEE Trans. Signal Process. 65, 794–816 (2016). https://doi.org/10.1109/TSP.2016.2601299
M. Grant, S. Boyd, CVX: matlab software for disciplined convex programming, version 2.1. (2014). Accessed 28 Feb 2024. http://cvxr.com/cvx
J. Capon, High-resolution frequency-wavenumber spectrum analysis. Proc. IEEE 57, 1408–1418 (1969). https://doi.org/10.1109/PROC.1969.7278
X. Yu, K. Alhujaili, G. Cui, V. Monga, MIMO radar waveform design in the presence of multiple targets and practical constraints. IEEE Trans. Signal Process. 68, 1974–1989 (2020). https://doi.org/10.1109/TSP.2020.2979602
Z. Cheng, B. Liao, Z. He, J. Li, J. Xie, Joint design of the transmit and receive beamforming in MIMO radar systems. IEEE Trans. Veh. Technol. 68, 7919–7930 (2019). https://doi.org/10.1109/TVT.2019.2927045
A. Aubry, A. DeMaio, A. Farina, M. Wicks, Knowledge-aided (potentially cognitive) transmit signal and receive filter design in signal-dependent clutter. IEEE Trans. Aerosp. Electron. Syst. 49, 93–117 (2013). https://doi.org/10.1109/TAES.2013.6404093
F. Wang, H. Li, M.A. Govoni, Power allocation and co-design of multicarrier communication and radar systems for spectral coexistence. IEEE Trans. Signal Process. 67, 3818–3831 (2019). https://doi.org/10.1109/TSP.2019.2920598
A. Ahmed, Y.D. Zhang, Y. Gu, Dual-function radar-communications using QAM-based sidelobe modulation. Digital Signal Process. 82, 166–174 (2018). https://doi.org/10.1016/j.dsp.2018.06.018
Funding
This work was supported by the National Natural Science Foundation of China under Grant Nots. 61871203.
Author information
Authors and Affiliations
Contributions
Each author made a significant contribution to the article. JT was responsible for the simulation design and implementation and the original draft of the article. ZK Z participated in the simulation design and review and revision of the article. All authors reviewed and approved the final version of the article.
Corresponding author
Ethics declarations
Ethics approval and consent to participate
All procedures performed in this paper were in accordance with the ethical standards of research community.
Competing interests
We declare that we have no financial and personal relationships with other people or organizations that can inappropriately influence our work and there is no professional or other personal interest of any nature or kind in any product, service and/or company that could be construed as influencing the position presented in or the review of the manuscript entitled.
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License, which permits any non-commercial use, sharing, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if you modified the licensed material. You do not have permission under this licence to share adapted material derived from this article or parts of it. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by-nc-nd/4.0/.
About this article
Cite this article
Tao, J., Zhang, Z. Synergistic optimization strategy for beamforming and power allocation in dual-functional radar-communication systems. EURASIP J. Adv. Signal Process. 2024, 89 (2024). https://doi.org/10.1186/s13634-024-01184-6
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s13634-024-01184-6