*Article* **Cramér-Rao Bound of Joint DOA-Range Estimation for Coprime Frequency Diverse Arrays**

**Zihuan Mao 1,2, Shengheng Liu 1,2,\*, Si Qin <sup>3</sup> and Yongming Huang 1,2**


**Abstract:** Frequency diverse array (FDA) produces a beampattern with controllable direction and range by slightly shifting the carrier frequencies across the elements, which is attractive in many applications. By further incorporating coprime array structure and coprime frequency offsets, improved degrees-of-freedom and spatial/range resolutions have been achieved. For such a relatively new array configuration, theoretical performance analyses are essential to explore the potentials and to facilitate practical implementation. In this work, we consider coprime-FDA-based joint/separate angle-range estimation of far-field targets that exhibit two different types of Swerling fluctuation behavior, which are respectively modelled as deterministic and stochastic sources. Analytical expressions of the Cramér–Rao bounds (CRB) and numerical simulations for both cases are provided. The results reveal that the relationship between CRB and coprime FDA parameters is not simply monotonic. As shown in the numerical simulations, the CRB of coprime FDA outperforms that of uniform FDA-MIMO for more than 60% under commonly-adopted coprime patterns. The presented results can be used as a guideline for optimal design of coprime FDA.

**Keywords:** Cramér–Rao bound (CRB); direction of arrival (DOA); parameter estimation; coprime array; frequency diverse array (FDA)

#### **1. Introduction**

The concept of frequency diverse array (FDA) was first proposed by Antonik et al. in 2006 [1,2]. In the pioneer work, a small and progressive frequency offset upon the carrier frequency is applied across the radiating elements. The use of element-to-element frequency offset generates a time-range-angle-dependent beampattern. Such a pattern is highly attractive for numerous applications including radar target localization [3–5], synthetic aperture radar (SAR) imaging [6,7], and wireless communications [8], since it enhances the flexibility of beam scanning and offers the potential to resist interference in a specific spatial region [9]. Nevertheless, the FDA beampattern generally exhibits periodicity in angle, range, and time [10,11]. Additionally, a standard FDA, i.e., a uniform linear FDA employing linearly increasing frequency offsets, yields coupling angle and range responses in the far-field beampattern. These unfavorable characteristics conspire to a spatiotemporal periodic S-shaped energy distribution in the angle-range plane, which further leads to ambiguity in target localization.

Owing to the fact that the FDA beampattern distribution can be controlled by tuning the frequency offsets, various solutions, such as time-dependent [12], logarithmic [13,14], and random [15] offsets, have been suggested to achieve a spatial-focusing and also preferably a time-invariant beampattern. Some sophisticated selection mechanisms have also been designed, which includes a transmit subarray strategy [16], adaptive selection scheme [17], and a multi-carrier transmission scheme [18,19], to name but a few. Additionally, by imposing certain evaluation metrics and penalty functions, the problem of frequency

**Citation:** Mao, Z.; Liu, S.; Qin, S.; Haung, Y. Cramér-Rao Bound of Joint DOA-Range Estimation for Coprime Frequency Diverse Arrays. *Remote Sens.* **2022**, *14*, 583. https://doi.org/ 10.3390/rs14030583

Academic Editors: Zhihuo Xu, Jianping Wang and Yongwei Zhang

Received: 8 December 2021 Accepted: 25 January 2022 Published: 26 January 2022

**Publisher's Note:** MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

**Copyright:** © 2022 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https:// creativecommons.org/licenses/by/ 4.0/).

<sup>1</sup> School of Information Science and Engineering, Southeast University, Nanjing 210096, China; mzh@seu.edu.cn (Z.M.); huangym@seu.edu.cn (Y.H.)

offset selection can be recast as an optimization problem [19]; different algorithms [20–22] have been employed to determine the optimal frequency offset.

The FDA angle-range estimation can also be decoupled from an array structure perspective. For example, many efforts have been made to incorporate the multiple-input multiple-output (MIMO) concept into FDA [4,23,24] to form a range-dependent-only beampattern. For MIMO enjoys the merits of increased degrees-of-freedom (DoF) and spatial diversity gain [25], it has attracted considerable interests. On the other hand, several types of sparse arrays have emerged to provide an enlarged array aperture and to overcome the DoF limitation for a given number of physical sensors. Among these sparse arrays, nested and co-prime configurations are compelling because they are more flexible and permit systematical design guided by explicit close-form expression [26,27]. But the latter is usually preferred, primarily due to the reduced mutual coupling [28,29]. FDA-MIMO radars exploit the DoF in the angle-range domain to jointly estimate the direction and range parameters of the targets [23]. In order to further achieve an improved DoF and spatial/range resolutions, a coprime FDA is proposed [30]. On the one hand, by incorporating the idea of coprime sampling to circumvent the limitation of physical sampling and using the virtual difference coarray concept, coprime FDA shows prominent advantages in terms of DOA-range resolution, localization accuracy, and the number of resolvable targets. Under the space-frequency difference equivalence of a coprime FDA, the unqiue non-negative lags of coarray produce extensive DoFs compared to the physical array [31]. On the other hand, the robustness of the coprime configuration in spatially filtering the interference [32] shows huge potential of adaptive beamforming using coprime FDA.

Objective performance analysis and optimal subarray design strategy are of paramount importance for arrays in practice, and they are often interconnected. In general, the Cramér– Rao bound (CRB) serves as a statistical benchmark to evaluate the degree of suboptimality that a certain unbiased parameter estimator exhibits. CRB also provides insights into the optimization of the array structure (cf., e.g., [20,33–35]). In the field of array signal processing, conventional research with respect to CRB mostly focuses on DOA estimation [36–38]. Motivated by the potentials of coprime FDAs and given the fact that their CRBs with respect to joint DOA-range estimation and the dependence between joint and separate estimations remain as yet uninvestigated, the technical contributions of this work are threefold, which are summarized as follows.


This paper is organized as follows. In Section 2, the configuration of coprime FDA and the transmit-receive steering vectors in deterministic case and stochastic cases are formulated. Based on the steering vectors, the array signal model is given. In Section 3, the explicit analytical form CRBs of joint estimation and separate estimation in deterministic signal cases are derived in terms of Gaussian distribution probability density function (PDF). When it comes to the stochastic signal case, the CRB is derived on the basis of the matrix form solution [42] in Section 4. Moreover, Section 5 presents some numerical simulations of CRBs in deterministic and stochastic source cases for different uniform FDA structures and coprime FDA structures. The comparisons of CRBs in different scenarios, i.e., two kinds of statistical characteristics and joint/separate estimation are analyzed. Ultimately, we make our conclusion in Section 6.

Notations: Lower (upper)-case bold characters are used to denote vectors (matrices). Vectors are by default in column orientation. (·)<sup>T</sup> and (·)<sup>H</sup> respectively represent the transpose and Hermitian operators of a matrix. tr(·) returns the trace of a matrix. Symbol stands for the Hardmard-product. <sup>E</sup>[·] returns the expected value of a discrete random variable. ∇ is the Nabla symbol. returns the real part of a complex value. Operators (·)<sup>⊥</sup> represents the ortho-complement of a projector matrix.

#### **2. Signal Model**

In this work we consider a coprime FDA that consists of two collocated uniform linear subarrays. The transmit and receive subarrays respectively have *N* and *M* elements, where *N* and *M* are coprime integers. The array configuration of the coprime FDA is illustrated in Figure 1.

**Figure 1.** Configuration of the coprime FDA under investigation.

We assume that the unit inter-element spacing *d* is *λ*/2, where *λ* denotes the wavelength. For a coprime FDA, a frequency increment as a multiple of Δ*f* is applied across the array element. As such, the carrier frequency of the *n*th transmit element can be expressed as

$$f\_n = f\_0 + (n-1) \cdot M \cdot \Delta f,\ n = 1, 2, \cdots, \prime, N,\tag{1}$$

where *f*<sup>0</sup> represents the FDA carrier frequency.

In array signal processing, the modelling of signal sources commonly falls into two categories [43], i.e., the deterministic and stochastic models. The former model category is adopted to describe the Swerling 0 type radar targets [39,40] and also in wireless communications [44]. The latter is related to Swerling I type targets and generic interferences [45] in radar applications as well as acoustic signals [46], where the signals are assumed to be driven by a Gaussian random process. In the following, we investigate radar detection performance for Swerling 0 and I targets, where the receive signals are respectively represented by deterministic and stochastic models.

#### *2.1. Deterministic Signal Model*

We first establish a deterministic signal model in this subsection. As stated above, we assume that a *priori* knowledge on the scattered waveform from the target is accessible to the radar operators. Consider a far-field Swerling 0 target at position (*θ*,*r*), where *r* ∈ [0, c/2Δ*f* ] denotes the range from the far-field target to the first element and *θ* ∈ [−*π*/2, *π*/2] denotes the angle between arrival direction and array normal. Then, in the transmit stage, the phase of the signal radiated from the first element is

$$p\_1 = \frac{2\pi f\_1}{\mathbf{c}}r.\tag{2}$$

Likewise, we can express the phase of the signal transmitted from the second element as

$$\varphi\_2 = \frac{2\pi f\_2}{\mathfrak{c}} (r - Md\sin\theta) = \frac{2\pi (f\_0 + M\Delta f)}{\mathfrak{c}} (r - Md\sin\theta). \tag{3}$$

The phase difference between the first and the second elements caused by the path length difference is

$$
\Delta q \varphi = q \varphi - \varphi\_1 = \frac{2\pi}{\mathfrak{c}} (-f\_0 Md \cdot \sin \theta + M \Delta f \cdot r - M^2 \Delta f d \cdot \sin \theta). \tag{4}
$$

Similarly, the phase of the signal transmitted from the *n*th transmit element is

$$\varphi\_n = \frac{2\pi f\_n}{\mathfrak{c}}(r - (n - 1)Md\sin\theta) = \frac{2\pi (f\_0 + M(n - 1)\Delta f)}{\mathfrak{c}}(r - (n - 1)Md\sin\theta). \tag{5}$$

The phase difference Δ*ϕ<sup>n</sup>* between the first and the *n*th transmit element is

$$
\Delta \boldsymbol{\varrho}\_{\text{fl}} = \boldsymbol{\varrho}\_{\text{l}} - \boldsymbol{\varrho}\_{1} = \frac{2\pi}{\mathbf{c}} \left[ -f\_{0}(\mathbf{n} - 1)M\boldsymbol{d} \cdot \sin\theta + (n - 1)M\Delta f \cdot \boldsymbol{r} - M^{2}(n - 1)^{2}\Delta f \boldsymbol{d} \cdot \sin\theta \right]. \tag{6}
$$

By taking the first transmit element as the reference, the transmit steering vector **a***t*(*θ*,*r*) is given by

$$\mathbf{a}\_{l}(\theta, r) = \begin{bmatrix} a\_{l1}(\theta, r), \ a\_{l2}(\theta, r), \ \cdots, \ a\_{lN}(\theta, r) \ \end{bmatrix}^{\mathrm{T}},\tag{7}$$

where *atn*(*θ*,*r*) = *<sup>e</sup>j*Δ*ϕ<sup>n</sup>* , *<sup>n</sup>* <sup>=</sup> 1, ··· , *<sup>N</sup>*, and *<sup>j</sup>* <sup>=</sup> √−1.

In the receive stage, the receive steering vector corresponding to the carrier frequency *fn* = *f*<sup>0</sup> + (*n* − 1)*M*Δ*f* can be derived as

$$\begin{array}{l} \mathbf{a}\_{\mathrm{rn}}(\theta,r) = \exp\left(\frac{j2\pi \cdot (n-1)M\Delta f \cdot r}{\mathbb{c}}\right) . \\\ \left[ \begin{array}{l} 1, \quad \exp\left(-j\frac{2\pi}{\mathsf{c}}(f\_{0}+(n-1)M\Delta f)Nd\sin\theta\right), \cdots \\\ \exp\left(-j\frac{2\pi}{\mathsf{c}}(f\_{0}+(n-1)M\Delta f)(M-1)Nd\sin\theta\right) \end{array} \right] \end{array} \tag{8}$$

The first term denotes the phase shift caused by carrier frequency *fn*, and the vector in the steering vector corresponds to the sparse receive array. To concentrate on the impact of coprime FDA structure, we consider an ideal case without mutual coupling and other hardware impairments in this work. Therefore, the overall transmit-receive steering vector **<sup>a</sup>**(*θ*,*r*) <sup>∈</sup> <sup>C</sup>*MN*×<sup>1</sup> which captures the phase difference generated in the transmit-receive stages can be written as

$$\mathbf{a}(\theta, r) = \left[ a\_{t1}(\theta, r) \cdot \mathbf{a}\_{r1\prime}^{\mathrm{T}} \cdot \cdots \cdot a\_{tN}(\theta, r) \cdot \mathbf{a}\_{rN}^{\mathrm{T}} \right]^{\mathrm{T}}.\tag{9}$$

More concretely, the entry in **a**(*θ*,*r*) that corresponds to the *n*th transmit element and the *m*th receive element is

$$a\_{\rm u,m}(\theta,r) = \exp\left(-j\frac{2\pi}{\mathfrak{c}}\left(f\_0M(n-1)d\sin\theta - 2\Delta f M(n-1)r + \Delta f M^2(n-1)^2d\sin\theta\right)\right)$$

$$\cdot \exp\left(-j\frac{2\pi}{\mathfrak{c}}\left(f\_0 + M(n-1)\Delta f\right)(m-1)Nd\sin\theta\right). \tag{10}$$

Overall, the output of matched filter corresponding to the *n*th transmit element and the *m*th receive element can be expressed as

$$y\_{n,m}(k) = a\_{n,m}(\theta, r)s\_{\text{det}}(k) + n\_{n,m}(k), \quad k = 1, \cdots, K,\tag{11}$$

where *s*det(*k*), *k* = 1, ··· , *K* denotes a signal with non-zero mean and zero variance [45]. Furthermore, *nn*,*m*(*k*) is the additive white Gaussian noise. Stacking *yn*,*m*(*k*) for all *n* =1, ··· , *N* and *m* = 1, ··· , *M* yields an *MN* × 1 vector and we have the deterministic receive signal vector

$$\mathbf{y}\_{\text{det}}(k) = \mathbf{a}(\theta, r)\mathbf{s}\_{\text{det}}(k) + \mathbf{n}(k), \quad k = 1, \cdots, \mathbb{K}. \tag{12}$$

#### *2.2. Stochastic Signal Model*

In this subsection we discuss the stochastic signal model that is used in Swerling I type target detection. Concretely, the statistics are governed by a chi-squared probability density function with two DoFs, and the echoed signal is modelled as a complex Gaussian process with zero mean and known variance. Suppose that we have a target at (*ψ*, *l*). By following a similar derivation to (12), we obtain the stochastic receive signal vector below:

$$\mathbf{y}\_{\text{sto}}(k) = \mathbf{a}(\psi, l)\mathbf{s}\_{\text{sto}}(k) + \mathbf{n}(k), \quad k = 1, \dots, K,\tag{13}$$

where *s*sto(*k*), *k* = 1, ··· , *K* denotes a stochastic signal driven by a complex Gaussian process N (0, *<sup>δ</sup>*<sup>2</sup> *<sup>s</sup>* ). We so far have established two different signal models for coprime FDA, which serve as the foundation for the following analyses of CRBs for joint DOA-range estimation in the following sections.

#### **3. CRB of Deterministic Signal for Coprime FDA**

#### *3.1. Deterministic Signal Model and CRB Derivation*

The deterministic CRB is the lower bound of the square error for target parameter estimation, which precisely defines the performance boundary of a specific array system, and is commonly used as the benchmark for the evaluation of estimation algorithms [43]. In this section, we derive the CRB of deterministic signal model.

Recall that the far-field target is located at (*θ*,*r*). The matched filter at output of the receive antenna array can be expressed as (12). The covariance matrix of the zero-mean complex Gaussian white noise vector **n**(*k*) is

$$\mathbf{R\_n} = \mathbb{E}\left[\mathbf{n}(k)\mathbf{n}^{\mathrm{H}}(k)\right] = \sigma\_n^2 \mathbf{I}\_{MN\_{\prime}} \tag{14}$$

where *σ*<sup>2</sup> *<sup>n</sup>* is the noise power, and **I***MN* denotes the *MN* × *MN* identity matrix.

For simplicity of derivation, the subscript of **y**det(*k*) is omitted and the estimation parameters (*θ*,*r*) are denoted as *α*. We first consider the single-snapshot case, in which the observed vector follows a Gaussian distribution, i.e., **<sup>y</sup>** ∼ N (**my**(*α*), *<sup>σ</sup>*<sup>2</sup> *<sup>n</sup>***I***MN*), where **my**(*α*) represents the mean value of the observed vector **y**. As such, the PDF of **y** is

$$P\_{\mathbf{y}|\mathbf{a}}(\mathbf{y}) = \frac{1}{\det(\pi \mathbf{R\_n})} \cdot \exp\left\{- (\mathbf{y} - \mathbf{m\_y}(\mathbf{a})) \mathbf{R\_n}^{-1} (\mathbf{y} - \mathbf{m\_y}(\mathbf{a}))^H \right\},\tag{15}$$

where det(·) denotes the matrix determinant. Then, the log-likelihood function of **y** is expressed as

$$\mathcal{L}\_{\mathbf{y}}(\mathfrak{a}) \stackrel{\triangle}{=} \ln P\_{\mathbf{y}|\mathfrak{a}}(\mathbf{y}) = -\ln \det(\pi \mathbf{R}\_{\mathbf{n}}) - \left\{ (\mathbf{y} - \mathbf{m}\_{\mathbf{Y}}(\mathfrak{a})) \mathbf{R}\_{\mathbf{n}}^{-1} (\mathbf{y} - \mathbf{m}\_{\mathbf{Y}}(\mathfrak{a}))^{H} \right\}. \tag{16}$$

The CRB gives the lower bound of the variance of unbiased estimation for parameter set *α*, which is denoted as **a**ˆ. That is to say, the covariance matrix **C**(*α*) - <sup>E</sup>[(*α*<sup>ˆ</sup> <sup>−</sup> *<sup>α</sup>*)(*α*<sup>ˆ</sup> <sup>−</sup> *<sup>α</sup>*)T] satisfies the following inequality:

$$\mathbf{C}(\mathfrak{a}) \succeq \mathbf{C} \mathbf{c} \mathbf{a}(\mathfrak{a}) \triangleq \mathbf{J}^{-1}, \tag{17}$$

where **J** is the Fisher information matrix (FIM) giving by

$$\mathbf{J} = -\mathbb{E}\left[\nabla\_{\mathbf{a}}\nabla\_{\mathbf{a}}\mathbf{T}^{\mathrm{T}}\mathcal{L}\_{\mathbf{y}}(\mathbf{a})\right],\tag{18}$$

and ∇*<sup>α</sup>* = *<sup>∂</sup> ∂θ ∂ ∂r* T . Furthermore, for deterministic signals,

$$J\_{i, \mathbf{j}} = \frac{\partial^2}{\partial \alpha\_i \partial \alpha\_j} [\mathcal{L}\_\mathbf{y}(\mathfrak{a})] = 2 \cdot \Re \left\{ \frac{\partial \mathbf{m}\_\mathbf{y}^\mathbf{H}(\mathfrak{a})}{\partial \alpha\_i} \cdot \mathbf{R}\_\mathbf{n}^{-1} \cdot \frac{\partial \mathbf{m}\_\mathbf{y}(\mathfrak{a})}{\partial \alpha\_j} \right\},\tag{19}$$

where *Ji*,*<sup>j</sup>* represents the element of the *i*th row and the *j*th column, *α<sup>i</sup>* represents the *i*th parameter in *α*. The detailed proof of (19) can be found in Appendix A. Therefore, the CRBs for parameter estimation are obtained as follows

$$\text{CRB}(\mathfrak{a}) = \begin{bmatrix} \text{CRB}(\theta) \\\\ \text{CRB}(r) \end{bmatrix} = \text{diag}(\mathbf{J}^{-1} \odot \mathbf{I}), \tag{20}$$

where diag(·) returns the diagonal elements of the matrix as a vector.

#### *3.2. CRB of Joint DOA-Range Estimation*

As mentioned above, the steering vector of coprime FDA is a function of both range and angle which leads to the capacity of joint DOA-range estimation. For the situation where the DOA and the range of the far-field target are both unknown and need to be estimated jointly, the corresponding FIM is expressed as

$$\mathbf{J} = \begin{bmatrix} J\_{\theta\theta} & J\_{\theta r} \\ & & \\ J\_{r\theta} & J\_{rr} \end{bmatrix}. \tag{21}$$

Based on the deterministic signal model established in Section 2.1, we first derive the analytical expression of Fisher information with respect to *θ*. Formula (19) is first rewritten as

$$J\_{1,1} = \frac{\partial^2}{\partial \theta^2} [\mathcal{L}\_\mathbf{y}(\mathbf{a})] = 2 \cdot \frac{\partial \mathbf{m}\_\mathbf{y}^\mathrm{H}(\mathbf{a})}{\partial \theta} \cdot \mathbf{R}\_\mathbf{n}^{-1} \cdot \frac{\partial \mathbf{m}\_\mathbf{y}(\mathbf{a})}{\partial \theta} = 2 \cdot \frac{A\_\mathbf{y}^2}{\sigma\_\mathbf{n}^2} \cdot \frac{\partial \mathbf{a}^\mathrm{H}(\theta, r)}{\partial \theta} \cdot \frac{\partial \mathbf{a}(\theta, r)}{\partial \theta},\tag{22}$$

where *A*<sup>2</sup> **<sup>y</sup>** = |**my**(*α*)| <sup>2</sup> represents the power of the deterministic signal. For each element *an*,*m*(*θ*,*r*) in the steering vector **a**(*θ*,*r*), we readily obtain

$$\begin{split} \frac{\partial a\_{n,m}(\theta,r)}{\partial \theta} &= a\_{n,m}(\theta,r) \cdot (-j\eta) \cdot \left\{ f\_0 \cdot \left[ M(n-1) + N(m-1) \right] \right. \\ &+ \Delta f \cdot \left[ M^2(n-1)^2 + MN(m-1)(n-1) \right] \Big\}\_{\prime} \end{split} \tag{23}$$

where *<sup>η</sup>* = <sup>2</sup>*π<sup>d</sup>* cos *<sup>θ</sup>*/c. Substituting (23) into (22) and considering the property <sup>Δ</sup>*<sup>f</sup>* <sup>2</sup> *<sup>f</sup>* <sup>2</sup> 0 , we have the following approximation

$$\begin{split} &\frac{\partial\mathbf{a}^{H}(\theta,r)}{\partial\theta} \cdot \frac{\partial\mathbf{a}(\theta,r)}{\partial\theta} = \eta^{2}f\_{0}^{2} \sum\_{m=0}^{M-1} \sum\_{n=0}^{N-1} (Mn+Nm)^{2} \\ &+ \eta^{2}\Delta f^{2} \sum\_{m=0}^{M-1} \sum\_{n=0}^{N-1} (M^{2}n^{2}+MNmn)^{2} + 2\eta^{2}\Delta f f\_{0} \sum\_{m=0}^{M-1} \sum\_{n=0}^{N-1} (M^{3}n^{3}+2M^{2}Nmn^{2}+MN^{2}m^{2}n) \\ &\approx \eta^{2}f\_{0}^{2} \left(M^{3} \sum\_{n=0}^{N-1} n^{2} + N^{3} \sum\_{m=0}^{M-1} m^{2} + 2MN \sum\_{m=0}^{M-1} m \sum\_{n=0}^{N-1} n \right) \\ &+ 2\eta^{2}\Delta f f\_{0} \left(M^{4} \sum\_{n=0}^{N-1} n^{3} + 2M^{2}N \sum\_{m=0}^{M-1} m \sum\_{n=0}^{N-1} n^{2} + MN^{2} \sum\_{m=0}^{M-1} m^{2} \sum\_{n=0}^{N-1} n \right). \end{split} \tag{24}$$

Accordingly, the Fisher information with respect to *θ* for the coprime FDA can be expressed as (25),

$$\begin{split} \mathcal{I}\_{\theta\theta} &= f\_{1,1} \\ &= 2 \cdot \text{SNR} \cdot \frac{4\pi^{2}d^{2}\cos^{2}\theta}{\mathfrak{c}^{2}} \cdot \left\{ f\_{0}^{2} \cdot \left( M^{3} \sum\_{n=0}^{N-1} n^{2} + N^{3} \sum\_{m=0}^{M-1} m^{2} + 2MN \sum\_{m=0}^{M-1} m \sum\_{n=0}^{N-1} n \right) \right. \\ &\left. + 2 \cdot \Delta f f\_{0} \cdot \left( M^{4} \sum\_{n=0}^{N-1} n^{3} + 2M^{2}N \sum\_{m=0}^{M-1} m \sum\_{n=0}^{N-1} n^{2} + MN^{2} \sum\_{m=0}^{M-1} m^{2} \sum\_{n=0}^{N-1} n \right) \right\} \\ &= 2 \cdot \text{SNR} \cdot \frac{4\pi^{2}d^{2}\cos^{2}\theta f\_{0}}{\mathfrak{c}^{2}} \cdot \left\{ f\_{0} \cdot \frac{\left( 7M^{3} - 6M^{2} + M^{3} \right)N^{3} + \left( -6M^{3} + 3M^{2} \right)N^{2} + M^{3}N}{6} \right. \\ &\left. + \Delta f \cdot \frac{\left( 9M^{4} - 7M^{3} + M^{2} \right)N^{4} + \left( -14M^{4} + 9M^{3} - M^{2} \right)N^{3} + \left( 5M^{4} - 2M^{3} \right)N^{2}}{6} \right\}. \tag{25} \end{split}$$

We next consider the Fisher information with respect to *r*. Similar to the derivation of (22), we obtain the Fisher information with respect to *r* as follows:

$$J\_{2,2} = \frac{\partial^2}{\partial r^2} [\mathcal{L}\_\mathbf{y}(\mathfrak{a})] = 2 \cdot \frac{A\_\mathbf{y}^2}{\sigma\_\mathbf{n}^2} \cdot \frac{\partial \mathbf{a}^\mathrm{H}(\theta, r)}{\partial r} \cdot \frac{\partial \mathbf{a}(\theta, r)}{\partial r}. \tag{26}$$

We also get

$$\frac{\partial a\_{n,m}(\theta,r)}{\partial r} = a\_{n,m}(\theta,r) \cdot \frac{j2\pi\Delta f}{\mathfrak{c}} \cdot (n-1)2M\_{\prime} \tag{27}$$

and

$$\frac{\partial \mathbf{a}^{\rm H}(\theta, r)}{\partial r} \cdot \frac{\partial \mathbf{a}(\theta, r)}{\partial r} = \left(\frac{2\pi \Delta f}{\mathbf{c}}\right)^2 \cdot 4M^3 \cdot \sum\_{n=0}^{N-1} n^2. \tag{28}$$

Thus, the Fisher information with respect to *r* is expressed as

$$\begin{split} \mathcal{J}\_{rr} &= \mathcal{J}\_{2,2} = 2 \text{SNR} \cdot \frac{4\pi^2 \Delta f^2}{\text{c}^2} \cdot 4M^3 \cdot \sum\_{n=0}^{N-1} n^2 \\ &= 2 \cdot \text{SNR} \cdot \frac{4\pi^2 \Delta f^2}{\text{c}^2} \cdot \left(\frac{8M^3 N^3 - 12M^3 N^2 + 4M^3 N}{6}\right) . \end{split} \tag{29}$$

Thus far, *Jθθ* and *Jrr* have been derived. In addition, we know that *Jθ<sup>r</sup>* = *Jr<sup>θ</sup>* and we have

$$J\_{1,2} = 2 \cdot \Re \left\{ \frac{\partial \mathbf{m}\_\mathbf{y}^\mathrm{H}(\mathbf{a})}{\partial \theta} \cdot \mathbf{R}\_\mathbf{n}^{-1} \cdot \frac{\partial \mathbf{m}\_\mathbf{y}(\mathbf{a})}{\partial r} \right\} = 2 \cdot \frac{A\_\mathbf{y}^2}{\sigma\_\mathbf{n}^2} \cdot \Re \left\{ \frac{\partial \mathbf{a}^\mathrm{H}(\theta, r)}{\partial \theta} \cdot \frac{\partial \mathbf{a}(\theta, r)}{\partial r} \right\}. \tag{30}$$

We further obtain

$$\begin{split} \frac{\partial \mathbf{a}^{\mathrm{H}}(\theta,r)}{\partial \theta} \cdot \frac{\partial \mathbf{a}(\theta,r)}{\partial r} &= \frac{8\pi^{2}d\cos\theta \Delta f}{c^{2}} \cdot \\ \left\{ f\_{0}(\mathcal{M}^{3} \sum\_{n=0}^{N-1} n^{2} + \mathcal{M}N \sum\_{m=0}^{M-1} m \sum\_{n=0}^{N-1} n) + \Delta f(\mathcal{M}^{4} \sum\_{n=0}^{N-1} n^{3} + \mathcal{M}^{2}N \sum\_{m=0}^{M-1} m \sum\_{n=0}^{N-1} n^{2}) \right\}. \end{split} \tag{31}$$

Hence, Fisher information with respect to *θ* and *r* is computed as (32),

$$\begin{split} I\_{\theta r} &= I\_{1,2} = 2 \cdot \text{SNR} \cdot \frac{8\pi^2 d \cos\theta \Delta f}{c^2} \\ &\quad \left\{ f\_0 \left( M^3 \sum\_{n=0}^{N-1} n^2 + M \sum\_{m=0}^{M-1} m \sum\_{n=0}^{N-1} n \right) + \Delta f \left( M^4 \sum\_{n=0}^{N-1} n^3 + M^2 N \sum\_{m=0}^{M-1} m \sum\_{n=0}^{N-1} n^2 \right) \right\} \\ &= 2 \cdot \text{SNR} \cdot \frac{8\pi^2 d \cos\theta \Delta f}{c^2} \cdot \left\{ f\_0 \cdot \frac{(7M^3 - 3M^2)N^3 + (-9M^3 + 3M^2)N^2 + 2M^3 N}{12} \right. \\ &\left. + \Delta f \cdot \frac{(5M^4 - 2M^3)N^4 + (-9M^4 + 3M^3)N^3 + (4M^4 - M^3)N^2}{12} \right\}. \end{split} \tag{32}$$

Substituting (25), (29) and (32) into (20) yields the analytical expression of the CRB for the joint DOA-range estimation, given as

$$\text{CRB}\_{\theta} = \left[ \mathbf{J}^{-1} \right]\_{1,1} = \frac{1}{J\_{\theta\theta} - \frac{\text{J}xLrt}{\text{J}rr}},\tag{33}$$

$$\text{CRB}\_r = \left[\text{J}^{-1}\right]\_{2,2} = \frac{1}{l\_{rr} - \frac{l\_{\theta r} l\_{xt}}{l\_{\theta \theta}}},\tag{34}$$

where · *<sup>i</sup>*,*<sup>j</sup>* denotes the element of the *i*th row and the *j*th column. For the *K*-snapshot case, the Fisher information are scaled up by a factor of *K* on the basis of (25), (29) and (32). As a result, the CRBs are scaled down by a factor of *K* on the basis of (33) and (34).

According to above results, we notice that all Fisher information increase with the coprime FDA parameters, i.e., coprime integers *M*, *N* and the carrier increment Δ*f* . Clearly, the Fisher information *Jθθ* of an FDA is generally greater than or equal to that of a phased array and, in particular, the equality holds when Δ*f* = 0. Nevertheless, this does not lead to the conclusion that CRB*<sup>θ</sup>* decreases with the increase of Δ*f* . The reason is that the CRB of joint DOA-range estimation is also influenced by the fisher information *Jθr*, which increases with Δ*f* . As a result, the mutual dependence of the DOA and the range estimations are stronger in the joint estimation case. As will be shown in the numerical results of Section 5, the relationship between CRB*<sup>θ</sup>* and Δ*f* is nonmonotonic.

#### *3.3. CRB of Separate Estimation*

If either parameter of *α* = (*θ*,*r*) is known, the problem of joint parameter estimation degrades to a separate parameter estimation problem. In this subsection, we assume a certain one of the target parameters is known from a previous detection [47,48]. Based on the Fisher information with respect to *θ* and *r* established in Section 2.2, we derive the analytical expression of CRB of DOA (range) estimation with coprime FDA parameters while the range (DOA) is known. The CRBs of separate estimation, which are denoted as SCRB, are given by

$$\text{SCRB}\_{\theta} = \frac{1}{f\_{\theta\theta}} \quad \text{and} \quad \text{SCRB}\_{\mathcal{I}} = \frac{1}{f\_{rr}}.\tag{35}$$

We readily observe from (35) that, SCRB*<sup>θ</sup>* decreases with the increase of the coprime FDA parameters, but the impact of Δ*f* to SCRB*<sup>θ</sup>* is very limited, for coprime FDA is narrowband in nature. This result agrees with our common sense that the performance of DOA

estimation is fundamentally determined by the array aperture. Likewise, SCRB*r* decreases with the increase of the coprime FDA parameters.

#### **4. CRB of Stochastic Signal for Coprime FDA**

*4.1. Stochastic Signal Model and CRB Derivation*

In the following, we examine the stochastic signal case and derive the CRB for parameter estimation in coprime FDA. For notational convenience, we also simplify *s*sto(*k*) and the estimation parameters (*ψ*, *l*) as *s*(*k*) and *ξ*, respectively. The mean and variance of the observed vectors **y**sto(*k*) are respectively assumed to be

$$\mathbf{m}\_{\mathbf{y}}(\xi) = 0,\tag{36}$$

$$\mathbf{R}\_{\mathbf{y}}(\xi) = \mathbf{a}\sigma\_s^2 \mathbf{a}^H + \sigma\_n^2 \mathbf{I}\_{MN}. \tag{37}$$

The log-likelihood function of **y** is expressed as

$$\mathcal{L}\_{\mathbf{y}}(\boldsymbol{\xi}) \stackrel{\triangle}{=} \ln P\_{\mathbf{y}|\boldsymbol{\xi}}(\mathbf{y}) = -\ln \det(\pi \mathbf{R}\_{\mathbf{y}}) - \mathbf{y} \mathbf{R}\_{\mathbf{y}}^{-1} \mathbf{y}^{H},\tag{38}$$

where det(·) represents the matrix determinant. Similar to the deterministic signal case, according to the matrix form solution [42], the FIM of stochastic signal can be written as

$$\mathbf{J} = -\mathbb{E}\left[\nabla\_{\xi}\nabla\_{\xi}\prescript{\mathbf{T}}{}{\mathbf{L}}\_{\mathbf{Y}}(\xi)\right] = \frac{2}{\sigma\_{\mathbf{n}}^{2}}\Re\{\left(\sigma\_{s}^{2}\,\mathbf{a}\mathbf{R}\_{\mathbf{y}}^{-1}\,\mathbf{a}\sigma\_{s}\right)\cdot\left(\nabla\_{\xi}\mathbf{a}^{\mathrm{H}}\mathbf{P}\_{\mathbf{a}}^{\perp}\,\nabla\_{\xi}\mathbf{a}\right)\},\tag{39}$$

where ∇*<sup>ξ</sup>* = *<sup>∂</sup> ∂ψ ∂ ∂l* T and **P**⊥ **<sup>a</sup>** = **<sup>I</sup>** − **<sup>a</sup>**(**a**H**a**)-1**a**H. In order to facilitate the numerical computing, the former factor in (39) is expressed as

$$\sigma\_s \mathbf{a}^{\text{H}} \mathbf{R}\_\mathbf{y}^{-1} \mathbf{a} \sigma\_s^2 = [\mathbf{I} - \mathbf{a}^{\text{H}} \mathbf{a} \Sigma (\mathbf{a}^{\text{H}} \mathbf{a} \Sigma + \mathbf{I})^{-1}] \mathbf{a}^{\text{H}} \mathbf{a} \Sigma = \frac{\text{MN} \cdot \text{SNR}}{1 + \text{MN} \cdot \text{SNR}} \cdot \sigma\_s^2,\tag{40}$$

where Σ = *σ*<sup>2</sup> *<sup>s</sup>* /*σ*<sup>2</sup> *<sup>n</sup>* = SNR. The readers are referred to the proof given in Appendix B for more details. The latter factor in (39) can be reformulated as

$$\nabla\_{\xi} \mathbf{a}^{\mathsf{H}} \mathbf{P}\_{\mathsf{A}}{}^{\bot} \nabla\_{\xi} \mathbf{a} = \nabla\_{\xi} \mathbf{a}^{\mathsf{H}} [\mathbf{I} - \mathbf{a} (\mathbf{a}^{\mathsf{H}} \mathbf{a})^{-1} \mathbf{a}^{\mathsf{H}}] \nabla\_{\xi} \mathbf{a} = \nabla\_{\xi} \mathbf{a}^{\mathsf{H}} \cdot \nabla\_{\xi} \mathbf{a} - \frac{1}{MN} \cdot \nabla\_{\xi} \mathbf{a}^{\mathsf{H}} \cdot \mathbf{a} \mathbf{a}^{\mathsf{H}} \cdot \nabla\_{\xi} \mathbf{a}. \tag{41}$$

Hence, the FIM of stochastic signal can be derived as

$$\mathbf{J} = 2 \cdot \frac{\text{MN} \cdot \text{SNR}^2}{1 + \text{MN} \cdot \text{SNR}} \cdot \left(\nabla\_{\xi} \mathbf{a}^{\text{H}} \cdot \nabla\_{\xi} \mathbf{a} - \frac{1}{\text{MN}} \cdot \nabla\_{\xi} \mathbf{a}^{\text{H}} \cdot \mathbf{a} \mathbf{a}^{\text{H}} \cdot \nabla\_{\xi} \mathbf{a}\right). \tag{42}$$

#### *4.2. CRB of Joint DOA-Range Estimation*

Similar to the analysis of the deterministic signals, in deriving the analytical form CRB of joint DOA-range estimation, the corresponding joint FIM is expressed as

$$\mathbf{J} = \begin{bmatrix} J\_{\Psi\Psi} & J\_{\Psi I} \\ & & \\ J\_{I\Psi} & J\_{II} \end{bmatrix}. \tag{43}$$

According to Formula (42), the Fisher information with respect to *ψ* is

$$J\_{1,1} = 2 \cdot \frac{\text{MN} \cdot \text{SNR}^2}{1 + \text{MN} \cdot \text{SNR}} \cdot \left(\frac{\partial \mathbf{a}^{\text{H}}(\boldsymbol{\psi}, l)}{\partial \boldsymbol{\psi}} \cdot \frac{\partial \mathbf{a}(\boldsymbol{\psi}, l)}{\partial \boldsymbol{\psi}} - \frac{1}{\text{MN}} \cdot \frac{\partial \mathbf{a}^{\text{H}}(\boldsymbol{\psi}, l)}{\partial \boldsymbol{\psi}} \cdot \mathbf{a} \mathbf{a}^{\text{H}} \cdot \frac{\partial \mathbf{a}(\boldsymbol{\psi}, l)}{\partial \boldsymbol{\psi}}\right). \tag{44}$$

For arbitrary element *an*,*m*(*ψ*, *l*) in the steering vector **a**(*ψ*, *l*), which corresponds to the *n*th transmit element and the *m*th receive element, we have

$$\frac{\partial a\_{\rm n,m}(\boldsymbol{\psi},l)}{\partial \boldsymbol{\psi}} = \frac{\partial a\_{\rm n,m}(\boldsymbol{\theta},r)}{\partial \boldsymbol{\theta}}|\_{\boldsymbol{\theta}=\boldsymbol{\psi},r=l}.\tag{45}$$

Since <sup>Δ</sup>*<sup>f</sup>* <sup>2</sup> *<sup>f</sup>* <sup>2</sup> <sup>0</sup> , we further have

$$\frac{\partial \mathbf{a}^{\mathrm{H}}(\boldsymbol{\psi},l)}{\partial \boldsymbol{\psi}} \cdot \frac{\partial \mathbf{a}(\boldsymbol{\psi},l)}{\partial \boldsymbol{\psi}} = \frac{\partial \mathbf{a}^{\mathrm{H}}(\boldsymbol{\theta},r)}{\partial \boldsymbol{\theta}} \cdot \frac{\partial \mathbf{a}(\boldsymbol{\theta},r)}{\partial \boldsymbol{\theta}}|\_{\boldsymbol{\theta}=\boldsymbol{\psi},r=l\prime} \tag{46}$$

and

$$\begin{aligned} \frac{\partial \mathbf{a}^{\mathrm{H}}(\boldsymbol{\psi},l)}{\partial \boldsymbol{\psi}} \cdot \mathbf{a} \mathbf{a}^{\mathrm{H}} \cdot \frac{\partial \mathbf{a}(\boldsymbol{\psi},l)}{\partial \boldsymbol{\psi}} &= \left\| \left| \frac{\partial \mathbf{a}^{\mathrm{H}}(\boldsymbol{\psi},l)}{\partial \boldsymbol{\psi}} \cdot \mathbf{a} \right| \right\|^{2} \\ \boldsymbol{\phi} &= \theta^{2} \left\{ f\_{0} \left( M^{2} \sum\_{n=0}^{N-1} n + N^{2} \sum\_{m=0}^{M-1} m \right) + \Delta f \left( M^{3} \sum\_{n=0}^{N-1} n^{2} + MN \sum\_{m=0}^{M-1} m \sum\_{n=0}^{N-1} n \right) \right\}^{2}, \end{aligned} \tag{47}$$

where *ϑ* = 2*πd* cos *ψ*/c. Substituting (47) into (44), the Fisher information with respect to *ψ* is derived as (48),

$$\begin{split} \mathcal{J}\_{\theta\theta} &= 2 \cdot \frac{\text{MN} \cdot \text{SNR}^2}{1 + \text{MN} \cdot \text{SNR}} \cdot \frac{4\pi^2 d^2 \cos^2 \theta f\_0}{\text{c}^2} \cdot \left\{ f \bullet \left( M^3 \sum\_{n=0}^{N-1} n^2 + N^3 \sum\_{n=0}^{M-1} M^2 + 2\text{MN} \sum\_{n=0}^{M-1} M \sum\_{n=0}^{N-1} \right) \right. \\ &+ 2 \cdot \Delta f \cdot \left( M^3 \sum\_{n=0}^{N-1} n^3 + 2M^2 N \sum\_{m=0}^{M-1} m \sum\_{n=0}^{N-1} n^2 + MN^2 \sum\_{n=0}^{M-1} m^2 \sum\_{n=0}^{N-1} n \right) - \frac{f\_0}{MN} \cdot \left( M^2 \sum\_{n=0}^{N-1} n + N^2 \sum\_{n=0}^{M-1} m \right)^2 \\ &- 2 \cdot \frac{\Delta f}{MN} \cdot \left( M^2 \sum\_{n=0}^{N-1} n + N^2 \sum\_{m=0}^{M-1} m \right) \left( M^3 \sum\_{n=0}^{N-1} n^2 + MN \sum\_{m=0}^{M-1} m \sum\_{n=0}^{N-1} n \right) \\ &= 2 \cdot \frac{MN \cdot \text{SNR}^2}{1 + MN \cdot \text{SNR}} \cdot \frac{4\pi^2 d^2 \cos^2 \psi f\_0}{\text{c}^2} \cdot \left\{ f\_0 \cdot \frac{(2M^3 - M)N^3 - M^3 N}{12} \right. \\ &+ \Delta f \cdot \frac{\left( 4M^4 - M^3 - M^2 \right) N^4 + \left( -3M^4 + M^2 \right) N^3 + \left( -3M^4 + M^3 \right) N^2 + 2M^4 N}$$

In the same way, the Fisher information with respect to *l* is

$$J\_{2,2} = 2 \cdot \frac{MN \cdot \text{SNR}^2}{1 + MN \cdot \text{SNR}} \cdot \left(\frac{\partial \mathbf{a}^H(\boldsymbol{\psi}, l)}{\partial l} \cdot \frac{\partial \mathbf{a}(\boldsymbol{\psi}, l)}{\partial l} - \frac{1}{MN} \cdot \frac{\partial \mathbf{a}^H(\boldsymbol{\psi}, l)}{\partial l} \cdot \mathbf{a} \mathbf{a}^H \cdot \frac{\partial \mathbf{a}(\boldsymbol{\psi}, l)}{\partial l}\right). \tag{49}$$

For arbitrary element *an*,*m*(*ψ*, *l*), we can get

$$\frac{\partial a\_{n,m}(\boldsymbol{\psi},l)}{\partial l} = \frac{\partial a\_{n,m}(\theta,r)}{\partial r}\Big|\_{\theta=\boldsymbol{\psi},r=l}.\tag{50}$$

Likewise, the first term in between the parentheses in (49) can be reformulated as

$$\frac{\partial \mathbf{a}^{\rm H}(\boldsymbol{\psi}, l)}{\partial l} \cdot \frac{\partial \mathbf{a}(\boldsymbol{\psi}, l)}{\partial l} = \frac{\partial \mathbf{a}^{\rm H}(\boldsymbol{\theta}, r)}{\partial r} \cdot \frac{\partial \mathbf{a}(\boldsymbol{\theta}, r)}{\partial r}|\_{\boldsymbol{\theta} = \boldsymbol{\psi}, r = l\prime} \tag{51}$$

and the latter term can be derived as

$$\frac{\partial \mathbf{a}^{\mathrm{H}}(\boldsymbol{\psi},l)}{\partial l} \cdot \mathbf{a} \mathbf{a}^{\mathrm{H}} \cdot \frac{\partial \mathbf{a}(\boldsymbol{\psi},l)}{\partial l} = \left| \left| \frac{\partial \mathbf{a}^{\mathrm{H}}(\boldsymbol{\psi},l)}{\partial l} \cdot \mathbf{a} \right| \right|^{2} = \frac{4\pi^{2}\Delta f^{2}}{\mathrm{c}^{2}} \left\{ 2M^{2} \sum\_{n=0}^{N-1} n \right\}^{2}.\tag{52}$$

Then, the Fisher information with respect to *l* is expressed as

$$J\_{ll} = 2 \cdot \frac{MN \cdot \text{SNR}^2}{1 + MN \cdot \text{SNR}} \cdot \frac{4\pi^2 \Delta f^2}{\text{c}^2} \cdot 4M^3 \cdot \left\{ \sum\_{n=0}^{N-1} n^2 - \frac{1}{N} \left( \sum\_{n=0}^{N-1} n \right)^2 \right\}$$

$$= 2 \cdot \frac{MN \cdot \text{SNR}^2}{1 + MN \cdot \text{SNR}} \cdot \frac{4\pi^2 \Delta f^2}{\text{c}^2} \cdot \frac{M^3 N^3 - M^3 N}{3} \cdot \tag{53}$$

In the scenarios where the DOA and the range of the stochastic signal source are both unknown and need to be estimated jointly, the Fisher information with respect to *ψ* and *l* is expressed as

$$J\_{\mathbb{M}} = J\_{1,2} = 2 \cdot \frac{MN \cdot \text{SNR}^2}{1 + MN \cdot \text{SNR}} \cdot \left(\frac{\partial \mathbf{a}^H(\boldsymbol{\upphi}, l)}{\partial \boldsymbol{\upphi}} \cdot \frac{\partial \mathbf{a}(\boldsymbol{\upphi}, l)}{\partial l} - \frac{1}{MN} \cdot \frac{\partial \mathbf{a}^H(\boldsymbol{\upphi}, l)}{\partial \boldsymbol{\upphi}} \cdot \mathbf{a} \mathbf{a}^H \cdot \frac{\partial \mathbf{a}(\boldsymbol{\upphi}, l)}{\partial l}\right). \tag{54}$$

The two terms in between the parentheses in (54) are respectively further developed as

$$\frac{\partial \mathbf{a}^{\mathrm{H}}(\boldsymbol{\psi},l)}{\partial \boldsymbol{\psi}} \cdot \frac{\partial \mathbf{a}(\boldsymbol{\psi},l)}{\partial l} = \frac{\partial \mathbf{a}^{\mathrm{H}}(\boldsymbol{\theta},r)}{\partial \boldsymbol{\theta}} \cdot \frac{\partial \mathbf{a}(\boldsymbol{\theta},r)}{\partial r}|\_{\boldsymbol{\theta}=\boldsymbol{\psi},r=l\prime} \tag{55}$$

and

$$\frac{\partial \mathbf{a}^{\mathrm{H}}(\boldsymbol{\psi},l)}{\partial \boldsymbol{\psi}} \cdot \mathbf{a} \mathbf{a}^{\mathrm{H}} \cdot \frac{\partial \mathbf{a}(\boldsymbol{\psi},l)}{\partial l} = \frac{8\pi^{2}d \cos\psi \Delta f}{c^{2}} \cdot \left\{ f\_{0} \left( M^{2} \sum\_{n=0}^{N-1} n + N^{2} \sum\_{m=0}^{M-1} m \right) \right. \\ \left. + \Delta f \left( M^{3} \sum\_{n=0}^{N-1} m \sum\_{n=0}^{N-1} n \right) \right\} \cdot M^{2} \sum\_{n=0}^{N-1} n. \tag{56}$$
 
$$\Delta f = \Delta f \left( M^{3} \sum\_{n=0}^{N-1} n^{2} + MN \sum\_{m=0}^{M-1} m \sum\_{n=0}^{N-1} n \right) \cdot M^{2} \sum\_{n=0}^{N-1} n. \tag{57}$$

As such, the Fisher information with respect to *ψ* and *l* is derived as (57),

$$\begin{split} \mathcal{J}\_{\mathsf{M}} &= 2 \cdot \frac{MN \cdot \text{SNR}^{2}}{1 + MM \cdot \text{SNR}} \cdot \frac{8\pi^{2}d \cos\theta \Delta f}{c^{2}} \\ &\quad \left\{ f\_{0} \cdot \left( M^{3} \sum\_{n=0}^{N-1} n^{2} + MN \sum\_{n=0}^{M-1} m \sum\_{n=0}^{N-1} n \right) + \Delta f \cdot \left( M^{4} \sum\_{n=0}^{N-1} n^{3} + M^{2} N \sum\_{n=0}^{M-1} m \sum\_{n=0}^{N-1} n^{2} \right) \\ &\quad - \frac{1}{MN} \left[ f\_{0} \cdot \left( M^{2} \sum\_{n=0}^{N-1} n + N^{2} \sum\_{n=0}^{M-1} m \right) - \Delta f \cdot \left( M^{3} \sum\_{n=0}^{N-1} n^{2} + MN \sum\_{n=0}^{M-1} m \sum\_{n=0}^{N-1} n \right) \right] \cdot \left( M^{2} \sum\_{n=0}^{N-1} n \right) \right\} \\ &= 2 \cdot \frac{MN \cdot \text{SNR}^{2}}{1 + MN \cdot \text{SNR}} \cdot \frac{8\pi^{2}d \cos\theta \Delta f}{c^{2}} \cdot \left\{ f\_{0} \cdot \frac{M^{3}N^{3} - M^{3}N}{12} + \Delta f \cdot \frac{(3M^{4} - M^{3})N^{4} - 5M^{4}N^{3} + M^{3}N^{2} + 2M^{4}N}{24} \right\}. \end{split}$$

Substituting (48), (53) and (57) into (43) yields the analytical expression of the CRB for the joint DOA-range estimation in stochastic signal case, i.e.,

$$\text{CRB}\_{\Psi} = \left[ \mathbf{J}^{-1} \right]\_{1,1} = \frac{1}{J\_{\Psi\Psi} - \frac{I\_{\Psi} J\_{l\Psi}}{J\_{\text{fl}}}} \, \tag{58}$$

$$\text{CRB}\_{l} = \left[\mathbf{J}^{-1}\right]\_{2,2} = \frac{1}{J\_{ll} - \frac{J\_{\Psi} \eta l \eta}{J\_{\Psi\Psi}}}.\tag{59}$$

Note that, the *<sup>K</sup>*-snapshot PDF *<sup>P</sup>***y**1,**y**2,··· ,**y***K*|*α*(**y**) is the product of *<sup>K</sup>* single-snapshot PDFs and the result resembles the aforementioned situation of the deterministic signal.

#### *4.3. CRB of Separate Estimation*

Similar to the deduction in the deterministic signal case, the analytical expression of CRB of DOA (range) estimation while range (DOA) is known can be expressed as

$$\text{SCRB}\_{\psi} = \frac{1}{f\_{\psi\psi}} \quad \text{and} \quad \text{SCRB}\_{l} = \frac{1}{f\_{ll}}.\tag{60}$$

We can see from (58), (59) and (60) that, the dependence of CRB*<sup>ψ</sup>* (SCRB*ψ*) and the coprime FDA parameters *M*, *N*, and Δ*f* is quite similar to that in the deterministic signal case. On the other hand, CRB*<sup>l</sup>* (SCRB*l*) decreases with the increase of the coprime FDA parameters. The detailed analyses will be provided in Section 5.

#### **5. Numerical Simulations and Analyses**

In this section, numerical simulations are presented to exemplify the CRBs in different scenarios. We assume that the coprime FDA under investigation operates at the X-band with a carrier frequency of *f*<sup>0</sup> = 10 GHz. Both deterministic and stochastic signal cases are investigated in the simulations. The location of the deterministic and stochastic signal sources are set as *α* = *ξ*, and the azimuth angle is 30◦ while the range is 2.8 km. To observe and analyze the CRB of coprime FDA, the following commonly-used coprime designs are simulated as examples: (a) *M* = 5, *N* = 3, Δ*f* = 15 kHz; (b) *M* = 7, *N* = 3, Δ*f* = 15 kHz; (c) *M* = 5, *N* = 4, Δ*f* = 15 kHz; (d) *M* = 5, *N* = 3, Δ*f* = 30 kHz; (e) *M* = 7, *N* = 3, Δ*f* = 30 kHz; (f) *M* = 5, *N* = 4, Δ*f* = 30 kHz. Note that the CRB results can be easily extended to other parameter settings using the closed-form solutions. We also include the CRBs of uniform FDA-MIMO [49] with the same number of array sensors in the numerical analyses to demonstrate the benefits of coprime configuration. The markers on the curves are used to distinguish between the CRBs of coprime and uniform FDAs.

Figure 2 shows the results of deterministic CRB*<sup>θ</sup>* and stochastic CRB*<sup>ψ</sup>* for joint estimation versus SNR. The zoomed-in view in Figure 2a,b indicates that the impact of the frequency increment values on the CRBs of DOA estimation is trivial. Figure 3 presents the curves of the deterministic CRB*<sup>r</sup>* and stochastic CRB*<sup>l</sup>* versus SNR. We can see from Figure 3 that, in both deterministic and stochastic cases, the CRBs of range estimation decrease with the increase of the coprime FDA parameters *M*, *N*, and Δ*f* . Furthermore, in contrast to the DOA estimation, the range estimation performance is significantly influenced by the unit frequency increment Δ*f* . As for the comparison between the coprime FDA and the uniform FDA-MIMO, we infer from the numerical results that the CRB of coprime FDA is more than 60% which is lower than that of uniform FDA-MIMO owing to the coprime configuration. The CRB performance for the detection of signals with different statistical characteristics is evaluated in Figure 4. Three groups of array design are simulated for deterministic and stochastic signal models. The simulation results suggest that the deterministic CRB outperforms the stochastic one, which tallies with the common knowledge in array signal processing.

In the following, we investigate the variation of angle CRB with the increase of frequency increment Δ*f* . A variable *∂*CRB is used to quantify the magnitude of the CRB variation, which is defined as

$$
\partial\_{\rm CRB} \stackrel{\Delta}{=} \frac{\rm CRB\_B - \rm CRB\_0}{\rm CRB\_0},
\tag{61}
$$

where CRB0 and CRB*<sup>B</sup>* denote the CRB of DOA when bandwidth is 0 and *B*, respectively. We let *B* and *κ* = *B*/ *f*<sup>0</sup> respectively denote the overall and relative bandwidth of the coprime FDA. To ensure that the narrow-band assumption of coprime FDA is met, we assume that *κ* is confined to ≤0.1. As depicted in Figure 5a, CRBs of DOA estimation in stochastic signal case decline with the increase of Δ*f* . The curves in Figure 5a also implies a non-monotonic dependence of the deterministic CRB*<sup>θ</sup>* on the bandwidth. From (32) we notice that *Jθ<sup>r</sup>* increases with the increase of Δ*f* . Hence, the dependence between

parameters in the joint estimation cases is stronger. As a result, deterministic CRB*<sup>θ</sup>* first decreases and then increases when the bandwidth reaches a certain level. Nevertheless, in the case of joint estimation of stochastic signals as shown in Figure 5b, this phenomenon is inconspicuous.

**Figure 2.** CRB of DOA in joint DOA-range estimation. (**a**) Deterministic signal case. (**b**) Stochastic signal case.

**Figure 3.** CRB of range in joint DOA-range estimation. (**a**) Deterministic signal case. (**b**) Stochastic signal case.

**Figure 4.** Comparison between two types of CRB. (**a**) CRB of DOA estimation in deterministic case and stochastic case. (**b**) CRB of range estimation in deterministic case and stochastic case.

**Figure 5.** Influence of Δ*f* on the DOA estimation performance. (**a**) *∂*CRB*<sup>θ</sup>* and *∂*CRB*<sup>ψ</sup>* versus *κ*. (**b**) *∂*SCRB*<sup>θ</sup>* and *∂*SCRB*<sup>ψ</sup>* versus *κ*.

We compare the CRBs of separate and joint estimations under both deterministic and stochastic signal models, which are respectively given in Figures 6 and 7. In addition, we examine the root mean square errors (RMSEs) of the classical spatial spectrum estimation algorithms including the two-dimensional MUltiple SIgnal Classification (2- D MUSIC) [50] and two-dimensional minimum variance distortionless response (2-D MVDR) [51], and compare them with the CRBs. We observe from Figure 8 that both the RMSEs of 2-D MUSIC and 2-D MVDR can approach the CRB gradually as SNR increases.

**Figure 6.** CRBs of joint/separate DOA-range estimations in deterministic signal case. (**a**) CRB of DOA estimation. (**b**) CRB of range estimation.

**Figure 7.** CRBs of joint/separate DOA-range estimations in stochastic signal case. (**a**) CRB of DOA estimation. (**b**) CRB of range estimation.

**Figure 8.** RMSEs of 2-D MUSIC and 2-D MVDR algorithm. The coprime FDA parameters are set as *M* = 5, *N* = 3 and Δ*f* = 15 kHz. The number of snapshots is set as 100 and the result is obtained through 500 Monte Carlo trials. (**a**) RMSE of DOA estimation. (**b**) RMSE of range estimation.

In addition to the analyses of the numerical results, from the analytical form expressions given in (33)–(35) and (58)–(60), we arrive at the following remarks.


Furthermore, we observe from Figure 2 that the CRBs of both DOA and range estimations decrease with the increase of coprime integers *M*, *N*, which leads to the expansion of array aperture and bandwidth. This result tallies with our common sense. Accordingly, we can use the largest *M*, *N* possible to design a theoretically best-performance coprime FDA. Nevertheless, in reality, array designs are often subject to resource constraints such as the total number of antenna elements and/or bandwidth. Therefore, a rational design with limited resources is important in practice. From the analytical CRB results given in the above sections, we notice that a larger total number of antenna elements *M* + *N* is not equivalent to superior CRB. The above conclusion is drawn in Figure 3a. More concretely, the CRB of range estimation for the coprime FDA choosing *M* = 5, *N* = 4, Δ*f* = 15 kHz, which uses a smaller number of elements, performs better than the coprime FDA choosing *M* = 7, *N* = 3, Δ*f* = 15 kHz. Based on this conclusion, we propose to use the CRB of coprime FDA to choose the optimal solution of coprime integers *M* and *N* to assist the array design.

Consider a practical scenario where the total number of antennas is *L* and the upper limit of bandwidth is *B*˜. By further denoting the coprime integers set as S*c*, we formulate the coprime FDA array design as the following optimization problem:

$$\begin{aligned} \min\_{M,N} & \quad \text{CRB}\_{\theta} + \gamma \cdot \text{CRB}\_{r}, \\ & \text{s.t.} \quad M + N \le L, (M - 1)N \cdot \Delta f \le \mathcal{B}, \\ & \quad \{M, N\} \in \mathbb{S}\_{c}. \end{aligned} \tag{62}$$

where *γ* represents the trade-off coefficient that controls the compromise between the CRB of DOA (range) estimation. This problem can be readily solved via exhaustive search. In doing so, we are able to design a theoretically optimal coprime FDA configuration with constrained physical resources.

#### **6. Conclusions**

In this work, the analytical expressions of the CRBs in deterministic/stochastic signal cases for joint DOA-range estimation using coprime FDA are derived. Extensive numerical simulations are also provided and analyzed. We are able to arrive at the following conclusions from the simulation results: (i) The CRB of DOA estimation is weakly related with the frequency increment; (ii) The difference between the deterministic and stochastic CRBs is distinct; (iii) The relation between frequency increment and the CRB of DOA for joint estimation is nonmonotonic in the deterministic signal case; (iv) The CRB of joint DOA-range estimation is slightly worse than the separate parameter estimation in both the deterministic and stochastic signal cases; and (v) The CRB of coprime FDA outperforms that of uniform FDA-MIMO for more than 60% under commonly-adopted coprime patterns. The difference between the CRBs of joint and separate DOA-range estimation methods is analyzed and we also propose an intuitive method for coprime FDA design based on CRB minimization.

**Author Contributions:** Conceptualization, Z.M. and S.L.; Data curation, Z.M.; Formal analysis, Z.M.; Funding acquisition, S.L. and Y.H.; Investigation, Z.M.; Methodology, Z.M.; Project administration, S.L.; Resources, S.L. and Y.H.; Software, Z.M.; Supervision, S.L. and Y.H.; Validation, S.L. and S.Q.; Visualization, Z.M.; Writing—original draft, Z.M. and S.L.; Writing—review & editing, S.L., S.Q. and Y.H. All authors have read and agreed to the published version of the manuscript.

**Funding:** This research was supported in part by the National Natural Science Foundationof China under Grant Nos. 62001103, U1936201 and the Basic Research Program of Jiangsu Province under Grant No. BK20190338.

**Institutional Review Board Statement:** Not applicable.

**Informed Consent Statement:** Not applicable.

**Data Availability Statement:** The data and code used in this study are available upon request to the corresponding author.

**Conflicts of Interest:** The authors declare that they have no conflict of interest regarding the publication of this manuscript.

#### **Appendix A**

The log-likelihood function L**y**(*α*) can be written as

<sup>−</sup> *Ji*,*<sup>j</sup>* <sup>=</sup> *<sup>∂</sup>*<sup>2</sup> *∂αi∂α<sup>j</sup>* [L**y**(*α*)] = tr <sup>−</sup>**R**−<sup>1</sup> **<sup>y</sup>** (*α*) *∂***Ry**(*α*) *∂α<sup>j</sup>* **R**−<sup>1</sup> **<sup>y</sup>** (*α*) *∂***Ry**(*α*) *∂α<sup>i</sup>* + **R**−<sup>1</sup> **<sup>y</sup>** (*α*) *∂*2**Ry**(*α*) *αiα<sup>j</sup>* − **<sup>y</sup>**<sup>H</sup> <sup>−</sup> **<sup>m</sup>**<sup>H</sup> **<sup>y</sup>** (*α*) ( <sup>−</sup> **<sup>R</sup>**−<sup>1</sup> **<sup>y</sup>** (*α*) *∂***Ry**(*α*) *∂α<sup>j</sup>* **R**−<sup>1</sup> **<sup>y</sup>** (*α*) *∂***Ry**(*α*) *∂α<sup>i</sup>* **R**−<sup>1</sup> **<sup>y</sup>** (*α*) + **R**−<sup>1</sup> **<sup>y</sup>** (*α*) *∂*2**Ry**(*α*) *αiα<sup>j</sup>* **R**−<sup>1</sup> **<sup>y</sup>** (*α*) + **Ry** <sup>−</sup>1(*α*) *∂***Ry**(*α*) *∂α<sup>j</sup>* **R**−<sup>1</sup> **<sup>y</sup>** (*α*) *∂***Ry**(*α*) *∂α<sup>i</sup>* **R**−<sup>1</sup> **<sup>y</sup>** (*α*)) **y** − **my**(*α*) − 2 4 <sup>−</sup> *<sup>∂</sup>*2**m**<sup>H</sup> **<sup>y</sup>** (*α*) *αiα<sup>j</sup>* **R**−<sup>1</sup> **<sup>y</sup>** (*α*) **y** − **my**(*α*) + *∂***m**<sup>H</sup> **<sup>y</sup>** (*α*) *∂α<sup>i</sup>* **R**−<sup>1</sup> **<sup>y</sup>** (*α*) *∂***Ry**(*α*) *∂α<sup>j</sup>* **R**−<sup>1</sup> **<sup>y</sup>** (*α*) **y** − **my**(*α*) + *∂***m**<sup>H</sup> **<sup>y</sup>** (*α*) *∂α<sup>j</sup>* **R**−<sup>1</sup> **<sup>y</sup>** (*α*) *∂***Ry**(*α*) *∂α<sup>i</sup>* **R**−<sup>1</sup> **<sup>y</sup>** (*α*) **y** − **my**(*α*) + *∂***m**<sup>H</sup> **<sup>y</sup>** (*α*) *∂α<sup>i</sup>* **R**−<sup>1</sup> **<sup>y</sup>** (*α*) *∂***my**(*α*) *∂α<sup>j</sup>* 5 . (A1)

where tr(·) returns the trace of the matrix. We notice that the first three terms in between the curly brackets are zero, and (60) can be simplified as

*Ji*,*<sup>j</sup>* <sup>=</sup> <sup>−</sup>tr4 <sup>−</sup> **<sup>R</sup>**−<sup>1</sup> **<sup>y</sup>** (*α*) *∂***Ry**(*α*) *∂α<sup>j</sup>* **R**−<sup>1</sup> **<sup>y</sup>** (*α*) *∂***Ry**(*α*) *∂α<sup>i</sup>* + **R**−<sup>1</sup> **<sup>y</sup>** (*α*) *∂*2**Ry**(*α*) *αiα<sup>j</sup>* + **R**−<sup>1</sup> **<sup>y</sup>** (*α*) *∂***Ry**(*α*) *∂α<sup>j</sup>* **R**−<sup>1</sup> **<sup>y</sup>** (*α*) *∂***Ry**(*α*) *∂α<sup>i</sup>* <sup>−</sup> **<sup>R</sup>**−<sup>1</sup> **<sup>y</sup>** (*α*) *∂*2**Ry**(*α*) *αiα<sup>j</sup>* <sup>−</sup> **<sup>R</sup>**−<sup>1</sup> **<sup>y</sup>** (*α*) *∂***Ry**(*α*) *∂α<sup>i</sup>* **R**−<sup>1</sup> **<sup>y</sup>** (*α*) *∂***Ry**(*α*) *∂α<sup>j</sup>* 5 + 2 4 *∂***m**<sup>H</sup> **<sup>y</sup>** (*α*) *∂α<sup>i</sup>* **R**−<sup>1</sup> **<sup>y</sup>** (*α*) *∂***my**(*α*) *∂α<sup>j</sup>* 5 <sup>=</sup> tr **R**−<sup>1</sup> **<sup>y</sup>** (*α*) *∂***Ry**(*α*) *∂α<sup>i</sup>* **R**−<sup>1</sup> **<sup>y</sup>** (*α*) *∂***Ry**(*α*) *∂α<sup>j</sup>* + 2 4 *∂***m**<sup>H</sup> **<sup>y</sup>** (*α*) *∂α<sup>i</sup>* **R**−<sup>1</sup> **<sup>y</sup>** (*α*) *∂***my**(*α*) *∂α<sup>j</sup>* 5 . (A2)

In the deterministic signal case, we have *<sup>∂</sup>***Ry**(*α*) *<sup>∂</sup>α<sup>i</sup>* <sup>=</sup> *<sup>∂</sup>***Ry**(*α*) *<sup>∂</sup>α<sup>j</sup>* = 0. As a result, (19) is obtained.

#### **Appendix B**

For the stochastic signal case, according to the theorem of Woodbury matrix identity **R**−<sup>1</sup> **<sup>y</sup>** = <sup>1</sup> *σ*2 *n* [**<sup>I</sup>** <sup>−</sup> *<sup>σ</sup>*<sup>2</sup> *s σ*2 *n* **<sup>a</sup>**( *<sup>σ</sup>*<sup>2</sup> *s σ*2 *n* **a**H**a** + **I**)−1**a**H] and denoting *σ*<sup>2</sup> *<sup>s</sup>* /*σ*<sup>2</sup> *<sup>n</sup>* (SNR) as Σ, we have

$$\begin{split} \mathbf{a}^{\operatorname{H}} \mathbf{R}\_{\mathbf{y}}^{-1} \mathbf{a} \sigma\_{s}^{2} &= \mathbf{a}^{\operatorname{H}} [\mathbf{I} - \mathbf{a} \Sigma (\mathbf{a}^{\operatorname{H}} \mathbf{a} \Sigma + \mathbf{I})^{-1} \mathbf{a}^{\operatorname{H}}] \mathbf{a} \Sigma \\ &= \mathbf{a}^{\operatorname{H}} \mathbf{a} \Sigma - \mathbf{a}^{\operatorname{H}} \mathbf{a} \Sigma (\mathbf{a} \mathbf{a}^{\operatorname{H}} \Sigma + \mathbf{I})^{-1} \mathbf{a}^{\operatorname{H}} \mathbf{a} \Sigma \\ &= [\mathbf{I} - \mathbf{a}^{\operatorname{H}} \mathbf{a} \Sigma (\mathbf{a}^{\operatorname{H}} \mathbf{a} \Sigma + \mathbf{I})^{-1}] \mathbf{a}^{\operatorname{H}} \mathbf{a} \Sigma \\ &= [(\mathbf{a}^{\operatorname{H}} \mathbf{a} \Sigma + \mathbf{I})(\mathbf{a}^{\operatorname{H}} \mathbf{a} \Sigma + \mathbf{I})^{-1} - \mathbf{a}^{\operatorname{H}} \mathbf{a} \Sigma (\mathbf{a}^{\operatorname{H}} \mathbf{a} \Sigma + \mathbf{I})^{-1}] \mathbf{a}^{\operatorname{H}} \mathbf{a} \Sigma \\ &= (\mathbf{a}^{\operatorname{H}} \mathbf{a} \Sigma + \mathbf{I})^{-1} \mathbf{a}^{\operatorname{H}} \mathbf{a} \Sigma = \frac{\mathbf{M} \mathcal{N} \cdot \text{SNR}}{1 + \mathbf{M} \mathcal{N} \cdot \text{SNR}}. \end{split} \tag{A3}$$

Thus, Formula (40) is obtained.

#### **References**


MDPI St. Alban-Anlage 66 4052 Basel Switzerland Tel. +41 61 683 77 34 Fax +41 61 302 89 18 www.mdpi.com

*Remote Sensing* Editorial Office E-mail: remotesensing@mdpi.com www.mdpi.com/journal/remotesensing

Academic Open Access Publishing

www.mdpi.com ISBN 978-3-0365-7955-9