Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
Next Article in Journal
Compatible Biomass Model with Measurement Error Using Airborne LiDAR Data
Next Article in Special Issue
Joint Design of Complementary Sequence and Receiving Filter with High Doppler Tolerance for Simultaneously Polarimetric Radar
Previous Article in Journal
Assimilating GNSS TEC with an LETKF over Yunnan, China
Previous Article in Special Issue
MIMO DFRC Signal Design in Signal-Dependent Clutter
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Technical Note

Joint Design of Transmitting Waveform and Receiving Filter via Novel Riemannian Idea for DFRC System

1
School of Communication Engineering, Hangzhou Dianzi University, Hangzhou 310018, China
2
School of Information Science and Engineering, Harbin Institute of Technology, Weihai 264209, China
*
Author to whom correspondence should be addressed.
Remote Sens. 2023, 15(14), 3548; https://doi.org/10.3390/rs15143548
Submission received: 22 April 2023 / Revised: 3 July 2023 / Accepted: 11 July 2023 / Published: 14 July 2023
(This article belongs to the Special Issue Theory and Applications of MIMO Radar)

Abstract

:
Recently, the problem of target detection in noisy environments for the Dual-Functional Radar Communication (DFRC) integration system has been a hot topic. In this paper, to suppress the noise and further enhance the target detection performance, a novel manifold Riemannian Improved Armijo Search Conjugate Gradient algorithm (RIASCG) framework has been proposed which jointly optimizes the integrated transmitting waveform and receiving filter. Therein, the reference waveform is first designed to achieve excellent pattern matching of radar beamforming. Furthermore, to ensure the quality of system information transmission, the energy of multi-user interference (MUI) of communication signals is incorporated as the constraint. Additionally, the typical similarity constraint is introduced to ensure the transmitting waveform with a good ambiguity function. Finally, simulation results demonstrate that the designed waveform not only enhances the system’s target detection performance in noisy environments but also achieves a relatively good multi-user communication ability when compared with other prevalent waveforms.

1. Introduction

In the past decades, with the rapid development of commercial wireless communication and remote sensing image processing [1,2,3], the demand for available frequency bands has also increased dramatically where most of the frequency resources are allocated to radar requirements and also limit the communication throughput. Therein, the integration of sensing and communication to improve spectrum utilization has become a research hotspot [4,5,6].
At present, the radar communication integration system can be divided into two categories according to whether the spectrum is shared: radar communication coexistence (RCC) [6] and dual functional radar communication (DFRC) [7]. For RCC, it means that radar and communication work independently on the same platform through some resource diversity, including frequency division multiplexing, power multiplexing and time division multiplexing [8]. Authors in [9] presented the power allocation and subcarrier selection scheme to minimize the transmission power while ensuring the presence of mutual information between radar and communication. Authors in [10] discussed the integration of radar communication by orthogonalizing radar and communication signals. Unlike RCC, DFRC integrates radar and communication functions through a single integrated waveform and allocates its power to a specific spatial area to detect targets while transmitting the user communication signals. Early research focused on the integrated waveform design where the digital information was embedded into the radar waveform by modifying the traditional radar waveform and controlling its sidelobe in the direction of the objective user [11,12]. Typically, when the communication user is in the mainlobe of the radar waveform, the communication rate would drop greatly, and the communication function would fail. Moreover, the communication symbols can also modulate the radar waveform to achieve coexistence for dual functions [13,14,15]. Authors in [16,17] tried to use the existing communication waveforms to achieve radar’s tasks. From the communication view, DFRC signals need high-quality communication performance, such as a high communication rate, which can be improved by minimizing multi-user interference (MUI).
Considering the scenarios involving multiple communication users and multiple radar sensing targets, the DFRC waveform has faced more challenges [18]. For communication, the signal-to-interference-plus-noise ratio (SINR) of each user should be considered. For radar, the SINR should also be considered. For DFRC systems, the performance of communication and sensing is coupled together, which means that any improvement of communication may deteriorate the radar performance, and vice versa. Therefore, in fact, when DFRC systems work in multi-user and multi-target scenarios, they inevitably face multiple performance tradeoffs between multi-users, multi-targets and also communication and perception. To solve this problem, authors in [19] designed the transmitting waveform by minimizing the joint least squares of weighted squared error and total MUI, where the weighting factor is used to balance two systems. Moreover, in this integration, most communication signals are modulated by multiple carriers, which will inevitably lead to a high peak average power ratio (PAPR) and incur distortion at the RF end. To improve the power efficiency of transmitters, PAPR constraints or constant modulus (CM) are widely used. Authors in [20] tried to jointly optimize the MUI and radar SINR by alternating the minimization and gradient projection frameworks. Authors in [21,22] designed the radar transmitting waveform under a given PAPR and similarity constraints. Note that radar waveforms as well as DFRC waveforms with constant modulus are also in need [23,24]. Authors in [25] designed a CM-integrated waveform by synthesizing different signals in the direction of communication and radar and also proposed an iterative optimization amplitude weighting method. No matter the PAPR/CM constraint or similarity constraint, these prior works would have to tackle the non-convex optimization problem with a heavy computation burden. How to design DFRC waveforms within the non-convex framework has been a hot topic. Furthermore, authors in [26] designed the CM DFRC waveform to minimize MUI and maximize the similarity between integrated signals and reference radar waveforms. Authors in [27] considered the joint design of the receiving filter and transmitting waveform with the maximum signal-to-noise ratio ( SNR ) and proposed a novel algorithm based on manifold ideas which give us lots of inspiration.
In this paper, to improve the detection performance of the integration system in noisy environments, the joint design of the system’s transmitting waveform and receiving filter has been proposed to enhance the output SNR , which is based on the Riemannian Improved Armijo Search Conjugate Gradient algorithm (RIASCG) framework. This framework could transform the non-convex optimization problem into a novel convex one within Riemannian manifold space. Firstly, the MIMO radar waveform, with the constant modulus and similarity (CM&S) constraint as well as good directivity, was designed. Furthermore, the waveform with the minimum MUI was also considered. By using the manifold principle, the CM&S constraint was transformed into an unconstrained Riemannian space. Particularly in the Riemannian space, the final solution can be obtained through the iterative closed-form. Finally, we compare their performance with several existing ones.
The organization of this paper is as follows. The system model and problem formulation are presented in Section 2 and Section 3. The DFRC waveform design by the novel algorithm is proposed in Section 4. Section 5 presents the numerical results. Finally, conclusions are drawn in Section 6.
Notation: Lower-case letters x and upper-case letters X denote vectors and matrices, respectively. The symbols ( ) T , ( ) H and ( ) * stand for the transpose, the conjugate transpose and the conjugate operators, respectively. The set of N × N complex matrices and the set of n -dimensional complex numbers vectors are denoted by N × N and n , respectively. The l 2 norm is denoted by the symbol | | | | 2 , and the Frobenius norm is represented by the symbol | | | | F . I N stands for the identity matrix of size N × N . Finally, the notations Re x and Im x are denoted as the real and imaginary part of x , E represents the expectation operator, the symbol denotes the Kronecker product and the symbol represents the Hadamard product.

2. System Model

Cognitive radar adjusts its waveforms via artificial intelligence or machine learning as shown in Figure 1, which is regarded as a closed-loop feedback cycle. This adaptive system makes it more intelligent and offers higher robustness in waveform optimization compared with the traditional one. This paper focuses on the design of transmitting waveforms and receiving filters in DFRC systems where the integrated waveform is suitable for both target detection and information transmission. Namely, to achieve this, the optimization problem should satisfy the transmitted beampattern favorable for target detection while generating minimal MUI to multiple users in the downlink. The DFRC system needs to perform the following two tasks simultaneously: (i) target detection and (ii) communication with single-antenna users in the downlink.

2.1. Communication Model

Assume that the channel between the bifunctional base-station and the communication user is flat Rayleigh fading and the channel characteristics remain unchanged for a certain period of time. The integrated MIMO radar-antenna system has N t transmitting antennas and N r receiving antennas, and the frame length of the dual-functional waveform is assumed to be N . Suppose that there are M communication users subjected to the interference from K irrelevant signals when detecting the target. Specifically, the signals received by M communication users can be represented as
Y = H X + W
where H = [ h 1 , h 2 , , h M ] T M × N t is the channel matrix. The transmission signal matrix is X = [ x 1 , x 2 , , x N ] N t × N , where x i N t × 1 denotes the i -th transmitted symbol vector, and W = w 1 , w 2 , , w N M × N is the white Gaussian noise matrix of the receiver.
Furthermore, as the integrated system would transmit a signal matrix S M × N to M users, (1) can be rewritten as
Y = S + H X S MUI + W
where the second term in (2) represents the multi-user interference [19], and its energy can be expressed as φ X = H X S F 2 with separable property H X S F 2 = i = 1 N H x i s i 2 . According to [20], the achievable sum-rate of the users can be defined as
ϑ m = 1 M log 2 1 + γ m
where γ m represents the SINR t of per-frame received by the m -th user, i.e.,
γ m = E s m , i 2 E h m T x i s m , i 2 + N 0
where s m , i represents the i -th code-unit for the m -th user, N 0 is the power of the received noise, and the energy of the m , i term of MUI can be expressed as E ( | h m T x i s m , i | ) . By minimizing the energy of MUI, the achievable rate of the system can be maximized, which is equivalent to minimizing φ X .

2.2. Detection Model

This section primarily aims to synthesize waveforms to achieve beampattern matching while maximizing SINR r . The SINR r of the MIMO system is determined by the transmitted waveform and its covariance matrix R . By optimizing the transmitting waveform vector, the quality of the radar-output signal can also be enhanced, which in turn improves the SNR and anti-interference ability of radar and also improves the sensitivity and accuracy of target detection. Assuming that a point target exists in the direction θ 0 , and K independent interference sources are located at θ k θ k θ 0 , k = 1 , 2 , , K , then the received signals at the i -th frame ( i = 1 , 2 , , N ) is formulated by
y i = α 0 a r θ 0 a t T θ 0 x i + k = 1 K α k a r θ k a t T θ k x i + v i
where α 0 , α 1 , , α K are the amplitudes of target and interference sources, and v i denotes the receiver noise. The a r θ N r and a t θ N t are the propagation vector and steering vector for the direction θ where the transmit and receive arrays are assumed to be linear uniform ones with a t θ = 1 N t e j π 0 sin θ , , e j π N t 1 sin θ T and a r θ = 1 N r e j π 0 sin θ , , e j π N r 1 sin θ T . To simplify the expression, the N vectors corresponding to y i in Equation (5) can be represented as
y ˜ = α 0 A 0 x ˜ + k = 1 K α k A k x ˜ + v ˜
where A k = I N a r θ k a t T θ k , y ˜ = [ y 1 T , , y N T ] T , x ˜ = [ x 1 T , , x N T ] T , v ˜ = [ v 1 T , , v N T ] T .
To improve the detection performance, it is required to process the received signal. Regarding the receiving filter w , the resulting output can be expressed as
r = w H y ˜ = α 0 w H A 0 x ˜ Target + w H k = 1 K α k A k x ˜ Interference + w H v ˜ Noise
In Equation (7), the first term represents the desired signal, the second term represents the interference signal, and the third term represents the noise. Thus, the SINR r of the filter output can be expressed as
SINR r = σ 0 2 w H A 0 x ˜ 2 w H k = 1 K σ k 2 A k x ˜ x ˜ H A k H w + σ v 2 w H w
where E [ | α k | 2 ] = σ k 2 represents the complex amplitude of α k . Generally, in a Gaussian noise environment, the larger the SINR r , the better the detection performance would be. Equation (8) can be transformed into a convex problem for a fixed x ˜ .
min w   w H k = 1 K σ ¨ k A k x ˜ x ˜ H A k H + I w s . t . w H A 0 x ˜ = 1
where σ ¨ k = σ k 2 / σ 0 2 . Furthermore, to accurately obtain some information of target or interference in the environment, waveforms in all directions should be equipped with some low sidelobes so as to reduce mutual interference. The beampattern power located at θ can be expressed as
P θ = a t H θ R a t θ
Considering that each array has the same emission energy, i.e., unitary power, the covariance matrix R can be designed as
R bb = 1 N t , b = 1 , , N t
where R bb represents the b , b -th element of covariance matrix R . ϕ θ denotes an expected transmit beampattern where { θ k } k = 1 K is a fine grid covering the points of interest. Assuming that there are K ˜ expected target locations, the objective is to detect the target at locations { θ ˜ k } k = 1 K ˜ , which can be confirmed by calculating the Capon spatial spectrum or generalized likelihood ratio test (GLRT) [28]. Here, the dominant peak position has been calculated by the GLRT pseudo-spectrum so as to form the desired beampattern { θ ˜ k } k = 1 K , and K is the resulting estimate of K ˜ , i.e.,
ϕ θ = 1 , θ θ ˜ k Δ 2 , θ ˜ k + Δ 2 , k 1 , , K 0 , o t h e r s
where Δ is the beamwidth selected by each target.
Based on the aforementioned discussions, it is necessary to design a matrix R that minimizes the least squares error between the transmitted beampattern P θ and the expected beampattern ϕ θ , while also minimizing the cross-correlation terms from different backscattered signals. Consequently, the corresponding problem could be expressed as
min R , α J R , α = 1 K k = 1 K ω k a t H θ k R a t θ k α ϕ θ k 2 + 2 ω c K ˜ K ˜ 1 p = 1 K ˜ 1 q = p + 1 K ˜ a t H θ ˜ p R a t θ ˜ q 2
where ω k represents the weight factor of the k -th source, ω c represents the weight factor of the cross-correlation term, and α is a scaling factor that needs to be optimized.
Moreover, some appropriate constraints are further imposed on the covariance matrix R , which also considers the need of the low cross-correlation beampattern. The designed R must be positive semidefinite, and all diagonal elements of R must be equal to the uniform antenna power while satisfying the uniform basic-power constraint. The covariance optimization problem can be formulated as
min R , α J R , α     s . t . R 0       R bb = 1 N t , b = 1 , , N t       α > 0
To solve the optimization problem of variables R , α , the convex optimization toolbox CVX can be employed [29]. Note that when designing MIMO radar systems with low cross-correlation beampatterns, it is crucial to make a balance between detection and communication performance. It is necessary to define the optimization objectives and constraints based on the specific application requirements. This ensures that the design effectively meets the desired performance criteria.

3. Optimization Modeling of Radar and Communication Integrated System

In this section, the optimization objective is to maximize SINR r while considering the signal-dependent interference. The design of the covariance matrix and transmitting waveform of MIMO radar will be addressed simultaneously. In practical radar applications, to improve the detection performance while maintaining the ability of multi-user communication, good ambiguity function and range resolution characteristics, it is necessary to introduce similarity constraints, such as
1 N N t x ˜ x ˜ 0 ξ
where x ˜ 0 is the reference waveform and ξ is the similarity coefficient. To enhance detection performance and facilitate multi-user communication, the permissible range of ξ is usually set as 0 ξ 2 c m , where c m = 1 / N N t . It should be noted that the communication performance may be severely degraded under the above constraints. Once the output result is obtained as the integrated referenced waveform X 0 , then it is equivalent to add the communication information into the integrated signal. Therefore, the optimization problem for the design of directional beamforming can be formulated as
min X 0 H X 0 S F 2   s . t . 1 N X 0 X 0 H = R
According to [19], the problem described in (16) can be characterized as an Orthogonal Procrustes problem (OPP) with a closed-form solution X 0 = N F U I N t × N V H , where R = F F H represents the Cholesky decomposition or other valid square-root decomposition, while U Σ V H = F H H H S represents the singular value decomposition (SVD). Consequently, the benchmark radar waveform X 0 is uncorrelated with the power of the desired constellation S .
Next, a compromise constraint, i.e., ρ H X S F 2 + 1 ρ X X 0 F 2 Υ , is proposed, where Υ defines the maximum permissible level for the communication performance metric and radar waveform similarity error. By taking the limited transmission energy into account, the optimization problem can be formulated as
min x ˜   w H k = 1 K σ ¨ k A k x ˜ x ˜ H A k H + I w   s . t . ρ H X S F 2 + 1 ρ X X 0 F 2 Υ     1 N N t x ˜ x ˜ 0 ξ   1 N X F 2 = P T
where 0 ρ 1 denote the weight factor, and P T represents the total power of all N t antennas per symbol. The constraint conditions can be transformed into the following composite form [19], i.e.,
    ρ H X S F 2 + 1 ρ X X 0 F 2 = [ ρ H T , 1 ρ I N t ] T X [ ρ S T , 1 ρ X 0 T ] T F 2
Denoting C = [ ρ H T , 1 ρ I N t ] T and D = [ ρ S T , 1 ρ X 0 T ] T , Equation (18) can be rewritten and extended to
C X D F 2 = tr C X D H C X D = tr X H C H C X tr X H C H D tr D H C X + tr D H D
Further defining Q = C H C and G = C H D , Equation (19) can be rewritten as
tr X H Q X 2 Re tr X H G
where Q is a Hermitian matrix, Equation (20) can be written in the form of a Lagrange multiplier with respect to the total power as follows
L X , λ = tr X H Q X 2 Re tr X H G + λ X F 2 N P T
where λ is the dual variable of the equality constraint. Defining X ˙ and λ ˙ as the optimal point and the dual optimal point with zero duality gap, according to the trust-region subproblem (TRS) optimality conditions [30], the following conclusions show
X ˙ = Q + λ ˙ I N t G Q + λ ˙ I N t G F 2 = V Λ + λ ˙ I N t 1 V H G F 2 = N P T λ ˙ λ m i n
where the notation refers to the Moore–Penrose pseudoinverse of the matrix. Furthermore, the matrix Q can be decomposed into Q = V Λ V H where λ m i n denotes the minimum eigenvalue of Q . It can be further proved that formula (22) has a unique solution, i.e.,
P λ = V Λ + λ I N t 1 V H G F 2 = n = 1 N t j = 1 N V H G n , j 2 λ + λ n 2
From deduction of (23), it is noted that when λ ˙ λ m i n , the function P λ is strictly decreasing and convex. Therefore, the golden-section search method can be employed to determine the optimal solution for λ ˙ .

4. Waveform Optimization Algorithm

In this section, a novel RIASCG algorithm is proposed to optimize the objective function. The Riemann gradient of function h x ˜ is defined as grad h x ˜ , which can be obtained by projecting the gradient on the Euclidean space. Here, Grad h x ˜ represents the Euclidean gradient of h x ˜ , and the contraction operator Retr maps the vector on the tangent space T x ˜ M s at the vicinity of manifold x ˜ M s . The next iteration point x ˜ l + 1 is considered when the objective value satisfies the descent condition. In each descent process, a more accurate step size d l needs to be selected to ensure faster convergence. To achieve this, an improved Riemannian manifold conjugate gradient algorithm based on the Armijo back-tracking line-search idea is proposed, offering several advantages over the first-order conjugate gradient algorithm [27]:
(1)
Faster convergence speed: The second-order conjugate gradient algorithm, utilizing second-order derivative information, could more accurately determine the search direction and step size compared with the first-order conjugate gradient one, resulting in better results in the same number of iterations.
(2)
More effective optimization for high-dimensional data: The first-order conjugate gradient algorithm may have a slow convergence speed when optimizing high-dimensional data, while the second-order conjugate gradient algorithm can better overcome this problem.
(3)
Stronger numerical stability: The second-order conjugate gradient algorithm can better avoid numerical instability, which is particularly prominent in optimizing high-dimensional data.
(4)
Fewer iterations: Due to faster convergence, the second-order conjugate gradient algorithm typically requires fewer iterations to achieve the same optimization effect, which is particularly important for optimizing large-scale data.
Suppose that the above constraints can be denoted as M s . For any z N N t , its projection operator of the sequence x ˜ l M s can be expressed as
Proj x ˜ l M s z = z Re z * x ˜ l x ˜ l
Once the input z N N t is given, it is possible to apply a universal compression function to effectively address common constraints, i.e.,
Retr z = arg   min x ˜ M s x ˜ z 2
The closed-form solution of w in Equation (17) can be derived, i.e.,
w = k = 1 K σ ¨ k A k x ˜ x ˜ H A k H + I 1 A 0 x ˜ x ˜ H A 0 H k = 1 K σ ¨ k A k x ˜ x ˜ H A k H + I 1 A 0 x ˜
As a result, a subproblem of the optimization problem has been obtained, i.e.,
min x ˜ x ˜ H A 0 H k = 1 K σ ¨ k A k x ˜ x ˜ H A k H + I 1 A 0 x ˜    s . t . x ˜ ϵ M s
The subproblem can be further expressed as
h x ˜ = x ˜ H A 0 H k = 1 K σ ¨ k A k x ˜ x ˜ H A k H + I 1 A 0 x ˜
The Euclidean gradient of the smooth extension h ˙ x ˜ can be denoted as Grad h ˙ x ˜ , i.e.,
Grad h ˙ x ˜ = 2 A 0 H k = 1 K σ ¨ k A k x ˜ x ˜ H A k H + I 1 A 0 x ˜ x ˜ H x ˜ A 0 H k = 1 K σ ¨ k A k x ˜ x ˜ H A k H + I 1 A 0 x ˜ = 2 A 0 H k = 1 K σ ¨ k A k x ˜ x ˜ H A k H + I 1 A 0 x ˜ 1 N N t x ˜ H A 0 H k = 1 K σ ¨ k A k x ˜ x ˜ H A k H + I 1 × k = 1 K σ ¨ k I N N t A k x ˜ x ˜ H x ˜ 1 x ˜ x ˜ H x ˜ N N t T A k H × k = 1 K σ ¨ k A k x ˜ x ˜ H A k H + I 1 A 0 x ˜
The steps of manifold RIASCG for DFRC waveform design can be summarized in Algorithm 1.
Algorithm 1: The Manifold RIASCG for DFRC Waveform Design.
Input:   l = 0 , x ˜ 0 , β , d , ϖ , η , X 0 , H , S , P T , σ 1 0 , 1 , σ 2 0 , 1 , weight factor 0 ρ 1 .
Output:  x ˜ l , w .
While  h x ˜ i + 1 h x ˜ i J  do
 1. Compute C = [ ρ H T , 1 ρ I N t ] T , D = [ ρ S T , 1 ρ X 0 T ] T , Q = C H C   and   G = C H D .
 2. Compute the eigenvalue decomposition Q = V Λ V H of Q , set the searching interval as [ λ m i n , b ], where b 0 is a searching upper-bound.
 3. Find the optimal solution λ ˙ to (22) using golden-section search.
 4. Compute X = Q + λ ˙ I N G .
 5. Once the global minimizer X is obtained, given its separability property, it can be employed as the reference waveform for the similarity constraint, denoted as x ˜ 0 .
 6. Compute Grad h ˙ x ˜ i according to (29).
 7. Compute Proj x ˜ i Grad h x ˜ i according to (24).
 8. Compute the improved Armijo back-tracking line-search parameter d i = ϖ η β , β is the smallest non-negative integer satisfying
    h x ˜ l h x ˜ l ϖ η β Grad h x ˜ l σ 1 d l Proj x ˜ l Grad h x ˜ l 2 2
                      + σ 2 d l 2 Proj x ˜ l Grad h x ˜ l 2 2
 9. Perform the projection step x ˜ ˙ l + 1 = x ˜ l d l Proj x ˜ l Grad h x ˜ l .
 10. Obtain   x ˜ l + 1 = Retr x ˜ l   according to (25).
 11. l = l + 1 .
End while
Compute w according to (26).

5. Numerical Result

To evaluate the performance of the designed waveforms, we first formulated the simulation scenarios. A uniform linear array configuration is assumed for both the transmitting and receiving arrays. Units in these arrays have an element spacing of half a wavelength. The number of transmitting antennas is N t = 10 , the number of receiving antennas is N r = 10 , and the target echo power is 10   dB . Additionally, the interference power is set to 20   dB , the noise power is 0   dB and the code-length of the waveform is N = 8 , while ϖ and η are drawn from [ 0 , 0.5 ] and [ 0 , 1 ] , respectively. For convenience, P T = 1 and assume that each entry of the channel matrix H is modeled as flat fading one and also is independently and identically distributed with a standard complex Gaussian distribution h i , j C N 0 , 1 . The constellation selected for the communication users is the unit-power QPSK alphabet, and the threshold value of h x ˜ i + 1 h x ˜ i is set as 10 4 .
By employing Equation (14), the objective is to construct the covariance matrix R for problem (16), and then its performance would be compared with the omni-directional waveform. Assume that there are four users, and the direction of arrival (DoA) information for K ˜ = 3 targets with unit complex amplitude is approximately { 50 , 0 , 50 } , which can be obtained by the Capon or GLRT method. Three symmetrical beampatterns of interest are denoted as θ ˜ 1 = 50 , θ ˜ 2 = 0 , θ ˜ 3 = 50 , with a beampattern width of Δ = 20 . Next, the performance of the omni-directional and directional beampattern would be discussed, and a radar-communication compromise waveform would be designed to achieve some flexible trade-off between radar and communication for practical needs. Considering the trade-off design for radar communication, the Pareto weight factor ρ = 0.2 is introduced, and ‘Omni’ and ‘Directional’ are denoted as the omnidirectional and directional beam-patterns. Further, the waveforms with strict equality constraints are denoted as ‘Strict’, while the trade-off designs are denoted as ‘Tradeoff’. In Figure 2, it is evident that the proposed method exhibits significantly lower sidelobe levels compared to the method of [19].
Moreover, to further assess the robustness of different algorithms, the proposed RIASCG would be compared with RCG [31], MM [32], MM-SQUAREM [33] and RCG-Armijo algorithms. To enhance the comparability, two additional constraints are also incorporated: the constant modulus constraint [34] and the e-uncertainty constant modulus constraint [35]. As shown in Figure 3, in comparison to the first-order RCG-Armijo algorithm, the second-order conjugate RIASCG algorithm used in this paper demonstrates faster convergence. The first-order conjugate gradient algorithm needs some precise line searching, which could increase computational costs and lead to some direction inconsistency in certain cases. Particularly in ill-conditioned problems, this algorithm might encounter direction loss within iterations and fail to converge to the minimum value. In such scenarios, the second-order conjugate gradient algorithm appears to be more advantageous than the first-order one, as it utilizes more information to determine the search direction and reduces the direction inconsistency. The RIASCG algorithm has demonstrated fewer iterations and a faster convergence rate compared with other algorithms.
To evaluate the scalability and robustness of algorithms, the computation runtime comparisons were also conducted for different N in Table 1 considering the CM&S constraint. Obviously, the practicality and scalability of RIASCG has effectively demonstrated that the comparative runtime performance has outperformed other prior works, such as RCG, RCG-Armijo, MM and MM-SQUAREM, especially for larger N . These comparisons could underscore the effectiveness and applicability of RIASCG.
Next, the effect of the similarity constraint has been discussed. As shown in Figure 4, when the similarity coefficient between two users is lower, it indicates a lower similarity between their signals, which can result in a greater degree of interference between them. When multiple users transmit signals simultaneously, the interference will affect the quality of the received signal and also result in a lower SINR r . In Figure 5, waveforms optimized by RIASCG algorithm as ξ decreasing have shown different shapes of ambiguity function. The similarity between the designed waveform and the reference waveform would be gradually increased, which would result in a better-formed ambiguity function. However, this also leads to a reduction in the degrees of freedom for waveform design.
Furthermore, given a specified radar beampattern, the design of dual-function waveforms is further analyzed. As demonstrated in Figure 6, we define the transmit SNR as SNR = P T / N 0 , and then examine the relationship between the transmit SNR of the communication signal and the average achievable sum-rate. Here, ϵ represents the stopping criterion for the golden-section search iteration. As the transmit SNR increases, the effect of the signal synthesis error on the average achievable sum-rate increases. Therefore, the sum-rate of the synthesized signal is slightly lower than that of a single communication signal.
Moreover, in Figure 7, the relationship between the average achievable rate per-user vs. the detection probability is illustrated when the receive SNR has 6   dB and the false-alarm probability has P F A = 10 7 . With an increase in the number of users, the detection probability would be decreased for a fixed average achievable rate per-user, which suggests that the increasing degrees of freedom could further minimize MUI energy. In Figure 8, the relationship between the average achievable rate per-user vs. the number of iterations has been demonstrated. Once the number of communication users increases, the feasible solution space decreases, which leads to a decline in the average achievable rate.

6. Conclusions

In this paper, we aimed to design an integrated waveform for DFRC using an RIASCG framework, which offers a flexible trade-off between radar and communication performance. To accomplish this, the manifold principle was leveraged to transform the constrained CM&S problem into unconstrained Riemann spaces, ensuring radar beampattern constraints and their trade-offs. Simulations have demonstrated the convergence performance and superiority of RIASCG when compared with other existing algorithms. Moreover, by adjusting the similarity coefficient, the designed waveform exhibited desirable properties in terms of the ambiguity function. As this paper demonstrated, we mainly focused on presenting a novel optimizing idea to tackle the joint optimization of integrated transmitting waveforms and receiving filters, which ignores the relative motion or Doppler shift. Next, in our future research, we will investigate the influence of relative motion and the Doppler shift and try to assess the robustness and adaptability of the proposed techniques in realistic scenarios. Additionally, we will also investigate the impact of varying channel conditions and different fading models, such as frequency-selective fading or multipath fading.

Author Contributions

Conceptualization, Y.Z. and X.F.; data curation, Z.Z. (Zhongqing Zhao) and F.T.; formal analysis, X.F. and P.S.; funding acquisition, X.F. and Y.Z.; investigation, Y.Z., X.F. and Z.Z. (Zhanfeng Zhao); methodology, X.F., Z.Z. (Zhanfeng Zhao) and P.S.; project administration, Y.Z., X.F. and Z.Z. (Zhongqing Zhao); resources, Z.Z. (Zhanfeng Zhao) and P.S.; supervision, X.F.; validation, Y.Z. and F.T.; writing—original draft, X.F. and Z.Z. (Zhongqing Zhao); writing—review and editing, X.F. and F.T. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported in part by the major scientific and technological innovation projects of Shandong Province of China (Grant No. 2021ZLGX05 and 2022ZLGX04).

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Liu, R.; Li, M.; Liu, Y.; Wu, Q.; Liu, Q. Joint Transmit Waveform and Passive Beamforming Design for RIS-Aided DFRC Systems. IEEE J. Sel. Top. Signal Process. 2022, 16, 995–1010. [Google Scholar] [CrossRef]
  2. Qian, X.; Zeng, Y.; Wang, W.; Zhang, Q. Co-Saliency Detection Guided by Group Weakly Supervised Learning. IEEE Trans. Multimed. 2023, 25, 1810–1818. [Google Scholar] [CrossRef]
  3. Lin, S.; Zhang, M.; Cheng, X.; Wang, L.; Xu, M.; Wang, H. Hyperspectral Anomaly Detection via Dual Dictionaries Construction Guided by Two-Stage Complementary Decision. Remote Sens. 2022, 14, 1784. [Google Scholar] [CrossRef]
  4. Yuan, W.; Wei, Z.; Li, S.; Yuan, J.; Ng, D.W.K. Integrated Sensing and Communication-Assisted Orthogonal Time Frequency Space Transmission for Vehicular Networks. IEEE J. Sel. Top. Signal Process. 2021, 15, 1515–1528. [Google Scholar] [CrossRef]
  5. Liu, R.; Li, M.; Liu, Q.; Swindlehurst, A.L. Joint Waveform and Filter Designs for STAP-SLP-Based MIMO-DFRC Systems. IEEE J. Sel. Areas Commun. 2022, 40, 1918–1931. [Google Scholar] [CrossRef]
  6. Labib, M.; Marojevic, V.; Martone, A.F.; Reed, J.H.; Zaghloui, A.I. Coexistence between Communications and Radar Systems: A Survey. URSI Radio Sci. Bull. 2017, 2017, 74–82. [Google Scholar] [CrossRef]
  7. Ma, D.; Shlezinger, N.; Huang, T.; Liu, Y.; Eldar, Y.C. Joint Radar-Communication Strategies for Autonomous Vehicles: Combining Two Key Automotive Technologies. IEEE Signal Process. Mag. 2020, 37, 85–97. [Google Scholar] [CrossRef]
  8. Chen, Y.; Gu, X. Time Allocation for Integrated Bi-Static Radar and Communication Systems. IEEE Commun. Lett. 2021, 25, 1033–1036. [Google Scholar] [CrossRef]
  9. Tian, T.; Li, G.; Zhou, T. Power Distribution for an OFDM-Based Dual-Function Radar-Communication Sensor. IEEE Sens. Lett. 2020, 4, 1–4. [Google Scholar] [CrossRef]
  10. Yu, X.; Yao, X.; Yang, J.; Zhang, L.; Kong, L.; Cui, G. Integrated Waveform Design for MIMO Radar and Communication via Spatio-Spectral Modulation. IEEE Trans. Signal Process. 2022, 70, 2293–2305. [Google Scholar] [CrossRef]
  11. Hassanien, A.; Amin, M.G.; Zhang, Y.D.; Ahmad, F. Dual-Function Radar-Communications: Information Embedding Using Sidelobe Control and Waveform Diversity. IEEE Trans. Signal Process. 2016, 64, 2168–2181. [Google Scholar] [CrossRef]
  12. Cheng, Z.; Wu, L.; Wang, B.; Shankar, M.R.B.; Ottersten, B. Double-Phase-Shifter Based Hybrid Beamforming for MmWave DFRC in the Presence of Extended Target and Clutters. IEEE Trans. Wirel. Commun. 2023, 22, 3671–3686. [Google Scholar] [CrossRef]
  13. Ahmed, A.; Zhang, Y.D.; Gu, Y. Dual-Function Radar-Communications Using QAM-Based Sidelobe Modulation. Digit. Signal Process. 2018, 82, 166–174. [Google Scholar] [CrossRef]
  14. Huang, T.; Xu, X.; Liu, Y.; Shlezinger, N.; Eldar, Y.C. A Dual-Function Radar Communication System Using Index Modulation. In Proceedings of the 2019 IEEE 20th International Workshop on Signal Processing Advances in Wireless Communications (SPAWC), Cannes, France, 2–5 July 2019; pp. 1–5. [Google Scholar]
  15. Yang, J.; Cui, G.; Yu, X.; Kong, L. Dual-Use Signal Design for Radar and Communication via Ambiguity Function Sidelobe Control. IEEE Trans. Veh. Technol. 2020, 69, 9781–9794. [Google Scholar] [CrossRef]
  16. Sen, S. PAPR-Constrained Pareto-Optimal Waveform Design for OFDM-STAP Radar. IEEE Trans. Geosci. Remote Sens. 2014, 52, 3658–3669. [Google Scholar] [CrossRef]
  17. Zhou, S.; Liang, X.; Yu, Y.; Liu, H. Joint Radar-Communications Co-Use Waveform Design Using Optimized Phase Perturbation. IEEE Trans. Aerosp. Electron. Syst. 2019, 55, 1227–1240. [Google Scholar] [CrossRef]
  18. Hassanien, A.; Amin, M.G.; Aboutanios, E.; Himed, B. Dual-Function Radar Communication Systems: A Solution to the Spectrum Congestion Problem. IEEE Signal Process. Mag. 2019, 36, 115–126. [Google Scholar] [CrossRef]
  19. Liu, F.; Zhou, L.; Masouros, C.; Li, A.; Luo, W.; Petropulu, A. Toward Dual-Functional Radar-Communication Systems: Optimal Waveform Design. IEEE Trans. Signal Process. 2018, 66, 4264–4279. [Google Scholar] [CrossRef] [Green Version]
  20. Tsinos, C.G.; Arora, A.; Chatzinotas, S.; Ottersten, B. Joint Transmit Waveform and Receive Filter Design for Dual-Function Radar-Communication Systems. IEEE J. Sel. Top. Signal Process. 2021, 15, 1378–1392. [Google Scholar] [CrossRef]
  21. De Maio, A.; Huang, Y.; Piezzo, M.; Zhang, S.; Farina, A. Design of Optimized Radar Codes With a Peak to Average Power Ratio Constraint. IEEE Trans. Signal Process. 2011, 59, 2683–2697. [Google Scholar] [CrossRef]
  22. Cheng, Z.; He, Z.; Liao, B.; Fang, M. MIMO Radar Waveform Design With PAPR and Similarity Constraints. IEEE Trans. Signal Process. 2018, 66, 968–981. [Google Scholar] [CrossRef]
  23. Cui, G.; Li, H.; Rangaswamy, M. MIMO Radar Waveform Design With Constant Modulus and Similarity Constraints. IEEE Trans. Signal Process. 2014, 62, 343–353. [Google Scholar] [CrossRef]
  24. BouDaher, E.; Hassanien, A.; Aboutanios, E.; Amin, M.G. Towards a Dual-Function MIMO Radar-Communication System. In Proceedings of the 2016 IEEE Radar Conference (RadarConf), Philadelphia, PA, USA, 2–6 May 2016; pp. 1–6. [Google Scholar]
  25. Jiang, M.; Liao, G.; Yang, Z.; Liu, Y.; Chen, Y. Integrated Radar and Communication Waveform Design Based on a Shared Array. Signal Process. 2021, 182, 107956. [Google Scholar] [CrossRef]
  26. Shi, S.; Wang, Z.; He, Z.; Cheng, Z. Constrained Waveform Design for Dual-Functional MIMO Radar-Communication System. Signal Process. 2020, 171, 107530. [Google Scholar] [CrossRef]
  27. Liu, F.; Masouros, C.; Li, A.; Sun, H.; Hanzo, L. MU-MIMO Communications With MIMO Radar: From Co-Existence to Joint Transmission. IEEE Trans. Wirel. Commun. 2018, 17, 2755–2770. [Google Scholar] [CrossRef]
  28. Bose, A.; Khobahi, S.; Soltanalian, M. Efficient Waveform Covariance Matrix Design and Antenna Selection for MIMO Radar. Signal Process. 2021, 183, 107985. [Google Scholar] [CrossRef]
  29. Grant, M.; SP, B. CVX: MATLAB Software for Disciplined Convex Programming. 2014. Available online: http://cvxr.com/cvx (accessed on 19 March 2023).
  30. Rendl, F.; Wolkowicz, H. A Semidefinite Framework for Trust Region Subproblems with Applications to Large Scale Minimization. Math. Program. 1997, 77, 273–299. [Google Scholar] [CrossRef]
  31. Alhujaili, K.; Monga, V.; Rangaswamy, M. Transmit MIMO Radar Beampattern Design via Optimization on the Complex Circle Manifold. IEEE Trans. Signal Process. 2019, 67, 3561–3575. [Google Scholar] [CrossRef] [Green Version]
  32. Hunter, D.R.; Lange, K. A Tutorial on MM Algorithms. Am. Stat. 2004, 58, 30–37. [Google Scholar] [CrossRef]
  33. Song, J.; Babu, P.; Palomar, D. Sequence Design to Minimize the Weighted Integrated and Peak Sidelobe Levels. IEEE Trans. Signal Process. 2015, 64, 2051–2064. [Google Scholar] [CrossRef] [Green Version]
  34. He, H.; Li, J.; Stoica, P. Waveform Design for Active Sensing Systems: A Computational Approach; Cambridge University Press: Cambridge, UK; New York, NY, USA, 2012. [Google Scholar]
  35. Zhao, L.; Song, J.; Babu, P.; Palomar, D.P. A Unified Framework for Low Autocorrelation Sequence Design via Majorization–Minimization. IEEE Trans. Signal Process. 2017, 65, 438–453. [Google Scholar] [CrossRef]
Figure 1. MIMO DFRC system based on cognitive radar idea.
Figure 1. MIMO DFRC system based on cognitive radar idea.
Remotesensing 15 03548 g001
Figure 2. Comparison of radar beampattern corresponding to different cases.
Figure 2. Comparison of radar beampattern corresponding to different cases.
Remotesensing 15 03548 g002
Figure 3. Comparison of convergence rates under different constraints. (a) Unconstrained; (b) Similarity constraints; (c) Constant modulus constraint; (d) e-uncertainty constant modulus constraint.
Figure 3. Comparison of convergence rates under different constraints. (a) Unconstrained; (b) Similarity constraints; (c) Constant modulus constraint; (d) e-uncertainty constant modulus constraint.
Remotesensing 15 03548 g003
Figure 4. SINR r vs. iteration number for different ξ .
Figure 4. SINR r vs. iteration number for different ξ .
Remotesensing 15 03548 g004
Figure 5. Ambiguity functions of the presented waveforms as ξ decreasing.
Figure 5. Ambiguity functions of the presented waveforms as ξ decreasing.
Remotesensing 15 03548 g005
Figure 6. SNR comparison of different stop conditions for golden-section search iterations.
Figure 6. SNR comparison of different stop conditions for golden-section search iterations.
Remotesensing 15 03548 g006
Figure 7. The average achievable rate per-user vs. radar detection probability.
Figure 7. The average achievable rate per-user vs. radar detection probability.
Remotesensing 15 03548 g007
Figure 8. The average achievable rate per-user vs. number of iterations.
Figure 8. The average achievable rate per-user vs. number of iterations.
Remotesensing 15 03548 g008
Table 1. Runtime comparisons of different methods under CM&S constraint.
Table 1. Runtime comparisons of different methods under CM&S constraint.
Algorithm N = 4 N = 8 N = 16 N = 32
RIASCG0.6875 s2.7031 s14.1719 s100.0785 s
RCG4.3751 s9.9843 s188.2811 s709.3284 s
RCG-Armijo0.8281 s3.6718 s40.1406 s232.0167 s
MM6.7968 s14.4375 s328.3755 s2562.1734 s
MM-SQUAREM1.5937 s5.5156 s113.8751 s630.3911 s
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Zhao, Y.; Zhao, Z.; Tong, F.; Sun, P.; Feng, X.; Zhao, Z. Joint Design of Transmitting Waveform and Receiving Filter via Novel Riemannian Idea for DFRC System. Remote Sens. 2023, 15, 3548. https://doi.org/10.3390/rs15143548

AMA Style

Zhao Y, Zhao Z, Tong F, Sun P, Feng X, Zhao Z. Joint Design of Transmitting Waveform and Receiving Filter via Novel Riemannian Idea for DFRC System. Remote Sensing. 2023; 15(14):3548. https://doi.org/10.3390/rs15143548

Chicago/Turabian Style

Zhao, Yinan, Zhongqing Zhao, Fangqiu Tong, Ping Sun, Xiang Feng, and Zhanfeng Zhao. 2023. "Joint Design of Transmitting Waveform and Receiving Filter via Novel Riemannian Idea for DFRC System" Remote Sensing 15, no. 14: 3548. https://doi.org/10.3390/rs15143548

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop