Hybrid Beamforming For DFRC System Based On SINR Performance Metric
Hybrid Beamforming For DFRC System Based On SINR Performance Metric
Yuanzhe Dai1 , Kaifeng Han2 , Guo Wei1 , Ying Du2 , Zhiqin Wang2
1
Department of Electronic Engineering and Information Science,
University of Science and Technology of China, Hefei, China
2
China Academy of Information and Communications Technology, Beijing, China
Email: 1 [email protected], 1 [email protected], 2 {hankaifeng, duying1, zhiqin.wang}@caict.ac.cn
978-1-6654-3944-2/21/$31.00
Authorized ©2021
licensed use limited to: INDIAN INSTITUTE TECHNOLOGY DELHI.82
OF IEEE Downloaded on April 05,2023 at 21:29:25 UTC from IEEE Xplore. Restrictions apply.
2021 IEEE/CIC International Conference on Communications in China (ICCC Workshops)
...
...
T
Let A (θ) = ar (θ) at (θ) , the SINR of radar receiver can be
Digital Analog
Beamformer Beamformer UE expressed by
...
...
...
FBB FRF K
α0 wH A (θ0 ) x2
K
NRF
Target of interest γr = " 2 #
PQ
BS receiver E wH αi A (θi ) x + wH w
i=1 (4)
2
Digital α0 wH A (θ0 ) FRF FBB s
Probing signal r MVDR = ,
Filter wH Rin w
w ...
Reflaction signal
Interference where
Nr
Q
X 2 H
Rin = |αi | A (θi ) FRF FBB FH H
BB FRF A (θi ) + I. (5)
Fig. 1. System model of the HAD structure DFRC system
i=1
constraints of analog beamforming matrix. Then the transmit solved through SDR technique by dropping the rank-1 con-
H
beamforming design problem can be expressed as straints, i.e. FBB,k = fBB,k fBB,k , then the problem (P2.1)
Ns can be relaxed as
X
H
FH
max fBB,k RF ΦFRF fBB,k K
{fBB,k },FRF
X
k=1 max tr Φ̃0 FBB,k
{FBB,k }
h FRF fBB,k 2
H
k=1
k
s.t. > Γk , ∀k,
h FRF fBB,i 2 tr H̃k FBB,k
P H
(P1.1) 1+ k (11) X
i6=k (P2.2) s.t. − tr H̃k FBB,i ≥ 1, ∀k,
Γk
i6=k
|FRF (i, j)| = 1, ∀i, j,
K
Ns X
tr FH
RF FBB,k FRF 6 P0 , FBB,k 0, ∀k,
X
H
fBB,i FH
RF FRF fBB,i 6 P0 ,
k=1
i=1 (13)
where Γk is the threshold of the k-th CU’s SINR. where H̃k = h̃k h̃H k . Problem (P2.2) is a convex problem and
This problem involves joint optimization over the digital it can be efficiently solved by CVX toolbox in MATLAB.
beamformer and analog beamformer. However, it is usually The computational complexity of solving this SDR problem
difficult to find their solutions simultaneously. Moreover, the is O K 4 NRF 2.5
[16]. Notice that the rank-1 constraints have
non-convex nature of individual SINR constraints of CUs and been dropped. However, it can be shown that the optimal
constant modulus constraint of analog beamforming matrix solution of problem (P2.2) has the rank-1 property by proving
makes this problem more challenging to solve. To address the following proposition with the theorem proposed in [17].
these challenges, this paper proposes an iterative algorithm Proposition 1: If the problem (P2.2) is feasible and bounded,
for obtaining a sub-optimal solution to (11). then there is always a optimal solution this problem that satisfy
rank(F?BB,k ) = 1, ∀k.
III. P ROPOSED B EAMFORMING D ESIGN A LGORITHM
Proof: See Appendix A.
In this section, we respectively provide digital and analog Therefore the optimal digital beamforming vector for k-th
beamforming design method, and then propose an overall CU can be given by fBB,k ? ?
with fBB,k ?H
fBB,k = F?BB,k .
iterative algorithm to solve the problem (P1.1). Notice that
both digital and analog beamforming design are non-convex B. Analog Beamforming Design
problem and there is no closed-form solution. Thus, we use
SDR technique and FPP-SCA algorithm to optimize the digital For analog beamforming, there are not only SINR con-
and analog beamforming matrices respectively. straints of CUs, but also constant modulus constraints of
analog beamforming matrix, and all of these constrains are
A. Digital Beamforming Design non-convex. Assume that we have obtained a fixed digital
Assuming that we already have a feasible solution of analog beamformer FBB . Then we use the sequential optimization
beamforming matrix FRF . For the sake of conciseness, let and define the constant matrix Φ0 as we did in III-A. The
Φ̃ = FH RF ΦFRF , h̃k = hk FRF , k ∈ {1...K}. Notice that
analog beamforming design problem can be written as
Φ̃ is not a constant matrix, thus the SINR of radar receiver K
(7) is a nonlinear function of the digital beamforming vector
X
tr FH
max RF Φ0 FRF FBB,k
fBB , which makes this problem difficult to solve. Therefore, FRF
k=1
we adopt the sequential optimization algorithm to optimize h FRF fBB,k 2
H
digital beamformer in an iterative fashion, where we fix Φ̃ k
s.t. P H 2 > Γk ,
as a constant matrix Φ̃0 in each iteration process. Then the (P3.1) 1+ hk FRF fBB,i (14)
digital beamforming optimization problem in each iteration i6=k
K
can be written as X
tr FH
Ns RF FBB,k FRF 6 P0 , ∀k,
X k=1
H
max fBB,k Φ̃0 fBB,k (12a)
{FBB,k }
k=1
|FRF (i, j)| = 1,
2
H H
h̃k fBB,k
where FBB,k = fBB,k fBB,k . Due to the SINR constraints of
(P2.1) s.t. P H 2 > Γk , ∀k, (12b) CUs and constant modulus constraints, this is a non-convex
1+ h̃k fBB,i QCQP so it is difficult to solve it directly. Nevertheless, using
i6=k the FPP-SCA method proposed in [16], a suboptimal solution
Ns
X
H
to (P3.1) can be obtained. To adopt FPP-SCA method, first
fBB,i FH
RF FRF fBB,i 6 P0 . (12c) we need to rewrite (P3.1) into the general form in [16]. By
i=1
the property of Kronecker product
This is still a non-convex QCQP problem because (12b) T
tr (ABCD) = vec DT CT ⊗ A vec (B) ,
are non-convex constraints. Nevertheless, it can be efficiently (15)
where A, B, C and D are arbitrary matrices, the problem Algorithm 1 Hybrid Transmit Beamforming Design Algo-
(P3.1) can be equivalently rewritten as rithm
Initialize: {θ0 , θ1 , ..., θQ }, {α0 ,q
α1 , ..., αQ }, [h1 , ..., hK ].
max uH BT ⊗ Φ0 u
(0)
u
Initialization: Let {FBB,k } = NtPN0RF [1, ..., 1]T , and for
T K (0) (0)
F BB,k ⊗ Hk X (i, j)-th element of FRF we let FRF (i, j) = ejϕ(i,j) , where
s.t. uH − FTBB,j ⊗ Hk u
Γk φ(i, j) is the phased of the (i, j)-th element of the conjugate
(P3.2) j6=k (16) transpose of the composite downlink channel, i.e., H =
> 1, ∀k, [h1 , ..., hK ]. Set l = 0.
uH BT ⊗ IN t u 6 P0 ,
Repeat
uH Em u 6 1, uH Em u > 1, ∀m , Set l = l + 1.
(l)
PK 1) Calculate digital beamforming matrix FBB as follows.
where B = i=1 FBB,k , u = vec (FRF ). And Em ∈ • Initialize m = 0 and convergence threshold δ1 .
C Nt NRF ×Nt NRF , where the m-th diagonal element of Em Repeat
is 1, and the remaining elements are 0. Then the following
– Set m = m + 1.
proposition shows the problem relaxed by FPP-SCA. (m−1)
– Calculate Φ̃0 = Φ̃(m−1) and γr according
Proposition 2: The problem (P3.2) can be relaxed
PK into the
T to (8) and (7), respectively.
problem (P3.4) by FPP-SCA, where B̄ = i6=k BB,i ⊗
F
– Calculate the optimal {F?BB,k } according to (13),
Hk , ∀i. Moreover, r, sk , pm and vm are slack variables. λ is
and decompose F?BB,k = fBB,k ? ?H
fBB,k to get the
penalty variable. z ∈ C Nt Nrf is an arbitrary vector, for which ?
optimal fBB,k .
we set z(n) = u(n−1) in the n-th iteration. (m)
Proof: See Appendix B. – Calculate Φ̃(m) and γr according to (8) and
The problem (P3.4) is a convex problem and can be solved (7).
(m) (m−1)
by CVX toolbox in MATLAB. The computational complexity Until |γr − γr | 6 δ1 .
of solving this problem is O([K + Nt Nrf ]3.5 ) [16]. (l)
2) Calculate analog beamforming matrix FRF as follows.
Compared to the traditional SCA method, FPP-SCA adds (l−1)
•Initialize m = 0, z(m) = vec(FRF ), and conver-
slack variables to sustain feasibility thus it does not need
gence threshold δ2 .
a feasible initial point. And the penalty method is used to
Repeat
force the slack variables toward zero, which helps pushing the
iterates towards the feasible point of problem (P3.1). Finally – Set m = m + 1.
(m−1)
the solution will converge to a KKT point of the original – Calculate Φ̃0 = Φ̃(m−1) and γr according
problem (P3.1), and a suboptimal solution to (P3.1) can be to (8) and (7), respectively.
obtained. – Calculate the optimal u(m) according to (17), and
set z(m) = u(m) .
(m)
C. Overall Transmit Beamforming Design Algorithm – Calculate Φ̃(m) and γr according to (8) and
Above all, the overall iteration algorithm of hybrid transmit (7).
(m) (m−1)
beamforming design can be proposed as Algorithm 1. The Until |γr − γr | 6 δ2
digital beamforming optimization and analog beamforming Until convergence.
optimization are alternately repeated until convergence.
Remark 1: Assuming that T0 is average number of outer
iteration between digital and analog beamforming optimiza-
tion, and in each outer iteration, TD is the average number of the reflection coefficient of interference is set as |αi |2 = 30dB,
inner iteration of digital beamforming optimization, and TA is i = 1, 2, 3. The channel vector of each CU is scattering
the average number of inner iteration of analog beamforming channel which is randomly generated.
design. Then the computational complexity of Algorithm 1 is The overall beampatterns for different antenna numbers of
O(T0 [TD (K 4 NRF2.5
)) + TA (K + Nt NRF3.5
)]) the transmitter and receiver is shown in Fig. 2. Assume that
there are 2 CUs and the SINR constraint of each CU is 15dB.
IV. S IMULATION R ESULTS And the number of RF-chains at BS transmitter is set to 2. It
In this section, we evaluate the proposed hybrid transmit is easy to see that the main lobes of the beampatterns point to
beamforming design algorithm via numerical simulations. As- θ0 = 10°, where the target of interest located. And nulls are
sume that both the BS transmitter and receiver are equipped placed in the direction of interferences, i.e., θ1 = −60°, θ2 =
with ULAs with half-wavelength interval between adjacent −30°, θ3 = 40°. As the number of the DFRC BS antennas
antennas. We set the BS transmit power as P0 = 50dBm. increases, the width of the main lobe decreases and the PSLR
A target of interest is located at the spatial angle θ0 = 0°, increases. This suggests that the increasing of the number of
and there are 3 interference located at the spatial angles BS antennas leads to better radar performance.
θ1 = −60°, θ2 = −30°, and θ3 = 40°. The reflection In Fig. 3, the tradeoff between the SINR constraints of CUs
coefficient of target of interest α0 is set as |α0 |2 = 10dB, and and the average SINR of radar is shown from 1,000 Monte-
uH BT ⊗ IN t u 6 P0
(17)
20
41
0
40
-20
-60
38
-80
37
-100
FPP-SCA
-120 36 Codebook
Full Digital
-140
35
-160 15 16 17 18 19 20 21 22 23 24 25
-80 -60 -40 -20 0 20 40 60 80 CU SINR (dB)
Angle(°)
Fig. 3. Tradeoff between the SINR constraints of CUs and the average SINR
Fig. 2. The overall beampatterns for different antenna numbers of DFRC BS, of radar for different design methods, Nt = Nr = 4, K = 2.
Nt = Nr , NRF = 2, K = 2, Γk = 15dB.
where {Cl }, {Aml } and {Xl } are Hermitian matrices (not [5] T. Lin, J. Cong, Y. Zhu, J. Zhang, and K. B. Letaief, “Hybrid beam-
necessarily positive semidefinite), bm ∈ R, Cm ∈ {>, 6, =}, forming for millimeter wave systems using the MMSE criterion,” IEEE
Trans. Commun., vol. 67, no. 5, pp. 3693–3708, May 2019.
∀m. [6] O. El Ayach, S. Rajagopal, S. Abu-Surra, Z. Pi, and R. Heath,
Suppose the above SDP and its dual is solvable. According “Spatially sparse precoding in millimeter wave MIMO systems,” IEEE
to [17], it always has an optimal solution {X?l } such that Trans.Wireless Commun., vol. 13, no. 3, pp. 1499–1513, Mar. 2014.
[7] F. Sohrabi and W. Yu, “Hybrid digital and analog beamforming design
L
X for large-scale antenna arrays,” IEEE J. Sel. Topics Signal Process., vol.
[rank(X?l )]2 6 M . (18) 10, no. 3, pp. 501–513, Apr. 2016.
[8] A. Hassanien and S. A. Vorobyov, “Phased-MIMO radar: A tradeoff
l=1 between phased-array and MIMO radars,” IEEE Trans. Signal Process.,
Then it is easy to prove that there is always an optimal solution vol. 58, no. 6, pp. 3137–3151, Jun. 2010.
[9] D. Wilcox and M. Sellathurai, “On MIMO radar subarrayed trans-
to the problem (P2.2) such that mit beamforming,” IEEE Trans. Signal Process., vol. 60, no. 4, pp.
K 2076–2081, Apr 2012.
X [10] F. Liu, C. Masouros, A. P. Petropulu, H. Griffiths, and L. Hanzo, “Joint
[rank(F?BB,k )]2 6 K + 1. (19) radar and communication design: Applications, State-of-the-art, and the
k=1 road ahead,” IEEE Trans. Commun., vol. 68, no. 6, pp. 3834–3862, Jun.
2020.
Because there is individual SINR constraint for each CU, the [11] F. Liu and C. Masouros, “Hybrid beamforming with sub-arrayed MIMO
optimal solution {F?BB,k }, satisfy rank(F?BB,k ) > 1, ∀k. Then radar: Enabling joint sensing and communication at mmwave band,” in
according to (19), proposition 1 can be proved. Proc. IEEE Int. Conf. Acoust. Speech Signal Process., Brighton, U.K.,
May 2019, pp. 7770–7774.
A PPENDIX B [12] Z. Cheng, B. Liao, Z. He, “Hybrid Transceiver Design for Dual-
Functional Radar-Communication System”, in Proc. 11th IEEE Sensor
P ROOF OF P ROPOSITION 2 Array Multichannel Signal Process. Workshop, Jun. 2020, pp. 8-11.
In order to account for potential infeasibility, the slack [13] Z. Cheng, Z. He and B. Liao, “Hybrid Beamforming for Multi-
carrier Dual-Function Radar-Communication System,” IEEE Trans.
variables can be added to the problem (P3.2), then we obtain Cogn. Commun. Netw., doi: 10.1109/TCCN.2021.3063110.
[14] C. Y. Chen and P. P. Vaidyanathan, “MIMO radar waveform optimization
max uH BT ⊗ Φ0 u − λr − λksk1 − λkpk1 − λkvk1
a,t,r,s,p,v with prior information of the extended target and clutter,” IEEE Trans.
Signal Process., vol. 57, no. 9, pp. 3533–3544, Sep. 2009.
T K [15] G. Cui, H. Li, and M. Rangaswamy, “MIMO radar waveform design
FBB,k ⊗ Hk X
s.t. uH − FTBB,j ⊗ Hk u with constant modulus and similarity constraints,” IEEE Trans. Signal
Γk Process., vol. 62, no. 2, pp. 343–353, Jan. 2014.
j6=k
[16] O. Mehanna, K. Huang, B. Gopalakrishnan, A.Konar, and N. D.
(P3.3) > 1 − sk , sk > 0, ∀k, Sidiropoulos, “Feasible point pursuit and successive approximation of
non-convex QCQPs,” IEEE Signal Process. Lett., vol. 22, no. 7, pp.
uH BT ⊗ IN t u 6 P0 ,
804-808, Jul. 2015.
[17] Y. Huang and D. P. Palomar, “Rank-constrained separable semidefinite
uH Em u 6 1 + pm , uH Em u > 1 + vm , ∀m. programming with applications to optimal beamforming,” IEEE Trans.
(20) Signal Process., vol. 58, no. 2, pp. 664–678, 2010.
Problem (P3.3) is still a non-convex QCQP, and the SCA tech-
nique can be adopted. For any positive semi-definite matrix
C ∈ C n×n and any z, u ∈ C n×1 ,we have (u−z)H C(u−z) >
0. Expanding the above inequality, we can obtain
uH Cu > 2Re zH Cu − zH Cz
(21)
We notice that B, Φ0 , {FBB,k } and {Hk } in (20) are all
positive semi-definite matrices. Thus using the liner restriction
like (21) around the point z, the non-convex problem P(3.3)
can be replaced by convex problem P(3.4). Proposition 2 is
proved.
R EFERENCES
[1] A. Hassanien, M. G. Amin, E. Aboutanios and B. Himed, “Dual-
Function Radar Communication Systems: A Solution to the Spectrum
Congestion Problem,” IEEE Signal Process. Mag., vol. 36, no. 5, pp.
115-126, Sep. 2019
[2] F. Liu, C. Masouros, A. Li, H. Sun, and L. Hanzo, “MU-MIMO commu-
nications with MIMO radar: From co-existence to joint transmission,”
IEEE Trans. Wireless Commun., vol. 17, no. 4, pp. 2755–2770, Apr.
2018.
[3] X. Liu, T. Huang, N. Shlezinger, Y. Liu, J. Zhou, and Y. C. Eldar, “Joint
transmit beamforming for multiuser MIMO communication and MIMO
radar,” IEEE Trans. Signal Process., vol. 55, no. 3, pp. 1213–1226, Jun.
2019.
[4] L. Chen, F. Liu, J. Liu and C. Masouros, “Composite Signalling
for DFRC: Dedicated Probing Signal or Not?” arXiv preprint arXiv:
2009.03528, 2020.