Open AccessArticle

Direct Position Determination of Non-Gaussian Sources for Multiple Nested Arrays: Discrete Fourier Transform and Taylor Compensation Algorithm

Hao Hu

^1,2,3,

Meng Yang

^1,2,

Qi Yuan

^1,2

Mingyi You

^4,5,

Xinlei Shi

^1,2,*

and

Yuxin Sun

^1,2

College of Electronic Information Engineering, Nanjing University of Aeronautics and Astronautics, Nanjing 211106, China

Key Laboratory of Dynamic Cognitive System of Electromagnetic Spectrum Space, Nanjing University of Aeronautics and Astronautics, Ministry of Industry and Information Technology, Nanjing 211106, China

Jiangsu Automation Research Institute, Lianyungang 222000, China

⁴

China Electronics Technology Group Corporation No.36 Research Institute, Jiaxing 314033, China

⁵

National Key Laboratory of Electromagnetic Space Security, Jiaxing 314033, China

Author to whom correspondence should be addressed.

Sensors 2024, 24(12), 3801; https://doi.org/10.3390/s24123801

Submission received: 16 May 2024 / Revised: 4 June 2024 / Accepted: 10 June 2024 / Published: 12 June 2024

(This article belongs to the Special Issue Robust Parameter Estimation with Sensor Arrays in Complex Electromagnetic Environments)

Download

Browse Figures

Versions Notes

Abstract

This paper delves into the problem of direct position determination (DPD) for non-Gaussian sources. Existing DPD algorithms are hindered by their high computational complexity from exhaustive grid searches and a disregard for the received signal characteristics by multiple nested arrays (MNAs). To address these issues, the paper proposes a novel DPD algorithm for non-Gaussian sources with MNAs: the Discrete Fourier Transform (DFT) and Taylor compensation algorithm. Initially, the fourth-order cumulant matrix of the received signal is computed, and the vectorizing method is applied. Subsequently, a computationally efficient DPD cost function is proposed by leveraging a normalized DFT matrix to reduce complexity. Finally, first-order Taylor compensation is utilized to enhance the accuracy of the localization results. The superiority of the proposed algorithm is demonstrated through numerical simulation results.

Keywords:

non-Gaussian signal; multiple nested arrays; direct position determination; Discrete Fourier Transform; Taylor compensation

1. Introduction

The advancement of wireless positioning technology has led to its widespread application in various sectors, including military defense, emergency rescue [1,2], resource exploration, intelligent transportation, and more [3,4,5]. Precisely locating enemy radiation sources is crucial for victory in battlefield scenarios. Wireless localization technology is categorized into traditional two-step positioning and direct position determination (DPD) technology [6,7,8].

The two-step positioning framework can be further classified into angle of arrival (AOA) [9], time of arrival (TOA) [10], time difference of arrival (TDOA) [11], received signal strength (RSS) [12], and other categories based on intermediate parameters estimation. While this technology estimates intermediate parameters from received signals and then solves spatial geometry problems to determine radiation source positions, it suffers from information loss between steps and suboptimal accuracy due to errors in parameter matching processes [13].

A DPD algorithm is suggested in [14] as a solution to address the challenges mentioned. This algorithm has garnered attention for its ability to achieve superior localization accuracy compared to the two-step algorithm, particularly in low signal-to-noise ratio (SNR) environments. The DPD algorithm directly determines the source position from the original received data, eliminating the need for estimating signal parameters and thereby avoiding estimation errors associated with intermediate parameters in the two-step approach. Moreover, the DPD algorithm does not require a matching procedure. While no intermediate parameters are essential in the DPD framework, certain signal parameters must still be taken into account in the algorithm model. The maximum likelihood (ML) DPD estimator in [15] incorporates joint information from angle of arrival (AOA) and time difference of arrival (TDOA) to achieve high localization accuracy through exhaustive search, albeit at the cost of increased complexity in scenarios with multiple sources. To mitigate this complexity, the subspace data fusion (SDF) technology [16] based on multiple signal classification (MUSIC) [17] is proposed. Furthermore, the Capon DPD algorithm, which avoids eigenvalue decomposition (EVD), is introduced for situations with multiple sources [18]. However, all the aforementioned algorithms treat the source signal as arbitrary, whereas studies have proved that better positioning accuracy can be achieved when the property of the source signal is a priori known [19,20]. A localization algorithm tailored for orthogonal frequency division multiplexing (OFDM) signals is detailed in [21]. Due to the high computational demands of ML algorithms, an extended SDF DPD algorithm is proposed in [22]. Additionally, features of non-circular signals are harnessed to enhance positioning accuracy and increase the degree of freedom (DOF) as evidenced by studies such as [23,24,25].

In the realm of localization research, it is commonly presumed that a signal adheres to a Gaussian distribution, with the second-order cumulant (SOC) [26] being utilized to derive a probability density function that encapsulates all signal information. However, practical scenarios often involve signals that deviate from Gaussian distribution, rendering lower-order cumulants insufficient in capturing all signal details. Therefore, the fourth-order cumulant (FOC) introduced in [27] is employed for signal analysis. Unlike SOC, FOC is adept at disregarding Gaussian noise and expanding array elements, leading to enhanced parameter estimation precision as evidenced in [28,29,30]. Moreover, prevalent position estimation techniques rely on uniform linear arrays (ULAs), which suffer from densely positioned elements, inducing heightened mutual coupling. Proposals to mitigate this challenge include sparse arrays with larger apertures and reduced mutual coupling, exemplified by classic coprime arrays and nested arrays (NAs) as put forth in [31,32,33].

In this paper, we propose a novel DPD algorithm utilizing multiple nested arrays (MNAs) for non-Gaussian sources. The main contributions are summarized as follows:

The property of non-Gaussian sources is fully exploited to suppress Gaussian noise and augment the virtual array aperture, which benefits the available degrees of freedom (DOFs).
We propose a novel low-complexity DPD algorithm of non-Gaussian sources utilizing MNAs. We deploy the Discrete Fourier Transform (DFT) method to construct a computationally efficient DPD cost function to reduce the high computational complexity caused by exhaustive grid search.
We utilize the Taylor compensation method to improve the localization accuracy at the expense of calculating the position estimation bias. It should be emphasized that even when the source position does not fall on the preset grid, the proposed algorithm can still estimate the position of sources accurately.
Complexity analysis and extensive numerical results are presented to verify the superiority of the proposed algorithm in terms of location accuracy, resolution capability, and computational complexity.

The following structure of the article is, the MNAs localization model is introduced in Section 2. The proposed algorithm is derived in Section 3. Performance analysis, numerical results, and conclusions are drawn in Section 4, Section 5 and Section 6, respectively.

Notations: Vectors and matrices are lower-case bold and upper-case bold, respectively.

{(\cdot)}^{T}

{(\cdot)}^{*}

, and

{(\cdot)}^{H}

denote transposition, conjugate, and Hermitian transpose, respectively.

x (n)

extracts the nth element of vector

x

;

E (\cdot)

represents the expectation operator; ⊗ is the Kronecker product;

vec (\cdot)

stacks the columns of a matrix into a vector;

\partial (x) / \partial (y)

denotes the derivation of x with respect to y;

diag (x)

turns vector

x

into a diagonal matrix; and

∥ \cdot ∥

denotes the

l_{2}

norm.

2. Model Formulation

In the context depicted in Figure 1, we analyze a scenario where K incoherent non-Gaussian sources emit narrow-band stationary signals in a two-dimensional localization setting, with unknown positions. It is assumed that the quantity of sources K is predetermined, and various established techniques can be employed to ascertain it [34,35]. To pinpoint these unidentified sources, we deploy L spatially distributed sensor arrays, each furnished with a NA comprising M array elements. The positions of the sources and sensor arrays are represented as

p_{k} = {[x_{k}, y_{k}]}^{T} (k = 1, \dots, K)

and

u_{l} = {[x_{l}, y_{l}]}^{T} (l = 1, \dots, L)

, respectively. The signal received by the lth sensor array at the tth

(t = 1, \dots, T)

snapshot can be formulated as per [36]:

x_{l} (t) = A_{l} (p) s (t) + n_{l} (t)

(1)

with the following notational definitions:

$p = {[p_{1}^{T}, \dots, p_{K}^{T}]}^{T} \in C^{2 K \times 1}$ is the position vector of K sources;
$A_{l} (p) = [a_{l} (p_{1}), a_{l} (p_{2}), \dots, a_{l} (p_{K})] \in C^{M \times K}$ is the array manifold of the lth NA, and $a_{l} (p_{k}) = {[e^{j 2 π d_{1} \frac{sin φ_{l, k}}{λ}}, e^{j 2 π d_{2} \frac{sin φ_{l, k}}{λ}}, \dots, e^{j 2 π d_{M} \frac{sin φ_{l, k}}{λ}}]}^{T} \in C^{M \times 1}$ denotes the array steering vector, where $sin φ_{l, k} = \frac{p_{k} (1) - u_{l} (1)}{∥u_{l} - p_{k}∥}$ . Array element locations $d_{1}, d_{2}, \dots, d_{M}$ are given by the set

$S = \{d, 2 d, \dots, N_{1} d, (N_{1} + 1) d, 2 (N_{1} + 1) d, \dots, N_{2} (N_{1} + 1) d\}$

(2)

where $M = N_{1} + N_{2}$ , and d denotes the space between array elements, which is usually set as half of wavelength $λ / 2$ ;
$s (t) = {[s_{1} (t), s_{2} (t), \dots, s_{K} (t)]}^{T} \in C^{K \times 1}$ is the signal vector transmitted from the kth source at time t, where $k = 1, 2, \dots, K$ , $t = 1, 2, \dots, T$ ;
$n_{l} (t)$ denotes the independent additive white Gaussian noise vector of the lth sensor array.

3. Proposed Algorithm

Define the FOC of a stationary stochastic process x as [37]:

C u m (x_{k_{1}}, x_{k_{2}}, x_{k_{3}}^{*}, x_{k_{4}}^{*}) = E (x_{k_{1}} x_{k_{2}} x_{k_{3}}^{*} x_{k_{4}}^{*}) - E (x_{k_{1}} x_{k_{2}}) E (x_{k_{3}}^{*} x_{k_{4}}^{*}) - E (x_{k_{1}} x_{k_{3}}^{*}) E (x_{k_{2}} x_{k_{4}}^{*}) - E (x_{k_{1}} x_{k_{4}}^{*}) E (x_{k_{2}} x_{k_{3}}^{*})

(3)

where

k_{1}

k_{2}

k_{3}

k_{4}

1 \leq k_{1}, k_{2}, k_{3}, k_{4} \leq T

are integers. And the FOC matrix of the received signal

x_{l} (t)

can be expressed as [38,39]:

R_{l, 4} = \sum_{k = 1}^{K} c_{s_{k}, 4} [a_{l} (p_{k}) \otimes a_{l}^{*} (p_{k})] {[a_{l} (p_{k}) \otimes a_{l}^{*} (p_{k})]}^{H} = \sum_{k = 1}^{K} c_{s_{k}, 4} a_{l, 4 x} (p_{k}) a_{l, 4 x}^{H} (p_{k})

(4)

where

a_{l, 4 x} (p_{k}) = a_{l} (p_{k}) \otimes a_{l}^{*} (p_{k}), 1 \leq k \leq K

c_{s_{k}, 4} = C u m (s_{k} (t), s_{k} (t), s_{k}^{*} (t), s_{l, k}^{*} (t))

denotes the FOC of the kth signal and is calculated using (3). Notably, the elements of

a_{l, 4 x} (p_{k})

can be constructed with

e^{j 2 π (d_{i} - d_{j}) \frac{sin φ_{l, k}}{λ}}

d_{i}, d_{j} \in S

. In practice,

a_{l, 4 x} (p_{k})

can be considered the steering vector of

R_{l, 4}

, which expands from

a_{l} (p_{k})

. To further improve the accuracy of localization, we vectorize

R_{l, 4}

by column [40]:

{\tilde{z}}_{l} = vec (R_{l, 4}) = {\tilde{A}}_{l} (p) c

(5)

where

{\tilde{A}}_{l} (p)

takes the form:

{\tilde{A}}_{l} (p) = [a_{l, 4 x}^{*} (p_{1}) \otimes a_{l, 4 x} (p_{1}), \dots, a_{l, 4 x}^{*} (p_{K}) \otimes a_{l, 4 x} (p_{K})] = [a_{l, v e c} (p_{1}), \dots, a_{l, v e c} (p_{K})]

(6)

and

c = {[c_{s_{1}, 4}, c_{s_{2}, 4}, \dots, c_{s_{K}, 4}]}^{T}

(7)

Without loss of generality, suppose that the number of sensors of two subarrays of each NA is equal, i.e.,

N_{1} = N_{2} = M / 2

. Then, the closed-form expression of the

a_{l, v e c} (p_{k})

location can be stacked as [40]

S_{v e c} = \{- M_{v e c} d, - (M_{v e c} - 1) d, \dots, (M_{v e c} - 1) d, M_{v e c} d |M_{v e c} = 2 N_{1} (N_{1} + 1) - 2\}

(8)

with many redundant elements. According to [41], matrix

A_{l, s o r t} (p) \in C^{(2 M_{v e c} + 1) \times K}

is formed by removing the redundancy from

a_{l, v e c} (p_{k})

. Meanwhile, the corresponding virtual signal vector

z_{l, s o r t}

is given by

z_{l, s o r t} = A_{l, s o r t} (p) c

(9)

Then, according to the geometry relationship showed in Figure 1, for a random grid

p = {[x, y]}^{T}

, we have

sin φ_{l} = \frac{x - x_{l}}{\sqrt{{(x - x_{l})}^{2} + {(y - y_{l})}^{2}}}

(10)

Let

q_{l} = \frac{J}{2} sin φ_{l}

, where

J = 2 M_{v e c} + 1

, yielding

q_{l} = \frac{J (x - x_{l})}{2 \sqrt{{(x - x_{l})}^{2} + {(y - y_{l})}^{2}}}

(11)

Define DFT vector

v_{l, p} \in C^{J \times 1}

, the

τ

th element of which is

v_{l, p} (τ) = {(e^{- j 2 π q_{l} / J})}^{τ}

(12)

Substitute (11) into (12), yielding

v_{l, p} (τ) = e^{\frac{- j π τ (x - x_{l})}{\sqrt{{(x - x_{l})}^{2} + {(y - y_{l})}^{2}}}}

(13)

where

1 \leq τ \leq J

. Combining the information of MNAs, we can construct the following DPD problem:

{\hat{p}}_{k}^{ini} = \arg max_{p} \sum_{l = 1}^{L} \sum_{τ = 1}^{J} v_{l, p} (τ) z_{l, s o r t} (τ)

(14)

where the estimated location

{\hat{p}}_{k}^{ini} = {[{\hat{x}}_{k}^{ini}, {\hat{y}}_{k}^{ini}]}^{T}

k = 1, 2, \dots, K

is based on the assumption that the position of the source falls on the preset grid. When the assumption is not met, the off-grid error will always exist. To overcome this issue, we apply the Taylor compensation method to the initial position estimation.

Define vectors

p_{x} = {[x_{1}, x_{2}, \dots, x_{K}]}^{T}

(15)

p_{y} = {[y_{1}, y_{2}, \dots, y_{K}]}^{T}

(16)

{\hat{p}}_{x}^{ini} = {[{\hat{x}}_{1}^{ini}, {\hat{x}}_{2}^{ini}, \dots, {\hat{x}}_{K}^{ini}]}^{T}

(17)

{\hat{p}}_{y}^{ini} = {[{\hat{y}}_{1}^{ini}, {\hat{y}}_{2}^{ini}, \dots, {\hat{y}}_{K}^{ini}]}^{T}

(18)

z_{s o r t} = {[z_{1, s o r t}^{T}, z_{2, s o r t}^{T}, \dots, z_{L, s o r t}^{T}]}^{T}

(19)

and matrix

A_{s o r t} (p) = {[A_{1, s o r t}^{T} (p), A_{2, s o r t}^{T} (p), \dots, A_{L, s o r t}^{T} (p)]}^{T}

(20)

According to (9), we have

z_{s o r t} = A_{s o r t} (p) c

(21)

Performing the first-order Taylor expansion of

A_{s o r t} (p)

{\hat{p}}^{ini} = {[{({\hat{p}}_{1}^{ini})}^{T}, \dots, {({\hat{p}}_{K}^{ini})}^{T}]}^{T}

, and ignoring the second-order and higher-order terms, we have:

A_{s o r t} (p) \approx A_{s o r t} ({\hat{p}}^{ini}) + \sum_{i = 1}^{2 K} {(p (i) - {\hat{p}}^{ini} (i)) \frac{\partial A_{s o r t} (p)}{\partial p (i)}|}_{p = {\hat{p}}^{ini}}

(22)

where

p (i)

and

{\hat{p}}^{ini} (i)

extract the ith element of

p

and

{\hat{p}}^{ini}

, respectively. Substitute (15)–(18) into (22), yielding

A_{s o r t} (p) \approx A_{s o r t} ({\hat{p}}^{ini}) + {\frac{\partial A_{s o r t} (p)}{\partial p_{x}^{T}}|}_{p_{x} = {\hat{p}}_{x}^{ini}} Λ_{x} + {\frac{\partial A_{s o r t} (p)}{\partial p_{y}^{T}}|}_{p_{y} = {\hat{p}}_{y}^{ini}} Λ_{y}

(23)

where

Λ_{x} = diag (δ_{x})

Λ_{y} = diag (δ_{y})

with

δ_{x} = p_{x} - {\hat{p}}_{x}^{ini}

δ_{y} = p_{y} - {\hat{p}}_{y}^{ini}

. Substitute (23) into (21), yielding

{\hat{z}}_{s o r t} \approx (A_{s o r t} ({\hat{p}}^{ini}) + {\frac{\partial A_{s o r t} (p)}{\partial p_{x}^{T}}|}_{p_{x} = {\hat{p}}_{x}^{ini}} Λ_{x} + {\frac{\partial A_{s o r t} (p)}{\partial p_{y}^{T}}|}_{p_{y} = {\hat{p}}_{y}^{ini}} Λ_{y}) c

(24)

Rewrite (24) as

{\hat{z}}_{s o r t} \approx [A_{s o r t} ({\hat{p}}^{ini}), {\frac{\partial A_{s o r t} (p)}{\partial p_{x}^{T}}|}_{p_{x} = {\hat{p}}_{x}^{ini}}, {\frac{\partial A_{s o r t} (p)}{\partial p_{y}^{T}}|}_{p_{y} = {\hat{p}}_{y}^{ini}}] [\begin{matrix} c \\ ω_{x} \\ ω_{y} \end{matrix}]

(25)

where

ω_{x} = Λ_{x} c

and

ω_{y} = Λ_{y} c

. Let

B = [A_{s o r t} ({\hat{p}}^{ini}), {\frac{\partial A_{s o r t} (p)}{\partial p_{x}^{T}}|}_{p_{x} = {\hat{p}}_{x}^{ini}}, {\frac{\partial A_{s o r t} (p)}{\partial p_{y}^{T}}|}_{p_{y} = {\hat{p}}_{y}^{ini}}]

y = [\begin{matrix} c \\ ω_{x} \\ ω_{y} \end{matrix}]

, the least square solution of

y

is readily given by

\hat{y} = {(B^{H} B)}^{- 1} B^{H} {\hat{z}}_{s o r t}

(26)

It should be noted that the estimated value

\hat{c}

is composed of the first to Kth elements of

\hat{y}

{\hat{ω}}_{x}

is composed of the

(K + 1)

th to

2 K

th elements of

\hat{y}

, and

{\hat{ω}}_{y}

is composed of the

(2 K + 1)

th to

3 K

th elements of

\hat{y}

. Therefore, the accurate position estimation after Taylor compensation is given by

{\hat{p}}_{x} = {\hat{p}}_{x}^{ini} + {\hat{ω}}_{x} . / \hat{c}

(27)

{\hat{p}}_{y} = {\hat{p}}_{y}^{ini} + {\hat{ω}}_{y} . / \hat{c}

(28)

The main steps of the proposed algorithm are summarized as in Algorithm 1.

Algorithm 1: Main Steps of the FOC-DFT-Taylor Algorithm
Input: ${\{x_{l} (t)\}}_{t = 1}^{T}$ , K, ${\{u_{l}\}}_{l = 1}^{L}$ , $λ$ , d, M, $N_{1}$ , $N_{2}$ .
	Calculate the FOC matrix ${\hat{R}}_{l, 4}$ according to (4); Vectorize the matrix ${\hat{R}}_{l, 4}$ and remove the redundancy to obtain the virtual signal vector $z_{l, s o r t}$ utilizing (5) and (8); Construct the normalized DFT vector according to (12) and obtain initial estimated values of the source position ${\hat{p}}_{k}^{ini} = {[{\hat{x}}_{k}, {\hat{y}}_{k}]}^{T}, k = 1, 2, \dots, K$ using (14); Perform the first-order Taylor expansion of $A_{s o r t} (p)$ at ${\hat{p}}^{ini}$ and construct the least square constraint utilizing (23) and (26); Obtain accurate position estimation after compensation according to (27) and (28).
Output: ${\hat{p}}_{x}$ , ${\hat{p}}_{y}$

4. Performance Discussion

4.1. Complexity

Here, the main computational complexity of different algorithms is compared in terms of complex number multiplication times. The parameters used are listed as follows: M is the number of array elements; K denotes the number of sources; L is the number of NAs; T represents the number of snapshot;

J = 2 M_{v e c} + 1

is the length of virtual signal vector

z_{l, s o r t}

; and we denote the number of search grids for x and y directions as recorded as

D_{x}

and

D_{y}

. The main complexity of the proposed algorithm lies in the calculation of the FOC matrix,

O (L T M^{4})

; grid search,

O (L J D_{x} D_{y})

; complex differentiation,

O (2 L K J)

; and compensation,

O (27 K^{3} + 9 K^{2} + 9 K^{2} L J + 3 K L J)

. For the sake of comparison, the complexity of the proposed algorithm (termed as FOC-DFT-Taylor), the classic SDF DPD algorithm with spatial smoothing technology (termed as FOC-SS-SDF-DPD) [22,42], the Capon DPD algorithm with spatial smoothing technology (termed as FOC-SS-Capon-DPD) [18], and the proposed DFT DPD algorithm with SOC (termed as SOC-DFT-DPD) are given in Table 1, where

N_{a} = (J + 1) / 2

denotes the length of the smoothing window.

Figure 2 shows the computational complexity of four different algorithms with the parameters setting as below:

M = 8

K = 12

N_{1} = N_{2} = 4

T = 500

J = 77

N_{a} = 39

D_{x} = D_{y} = 400

, and the number of NAs is changed from 3 to 10. As we can see, compared to the FOC-SS-SDF-DPD and FOC-SS-Capon-DPD algorithm, the proposed FOC-DFT-DPD and FOC-DFT-Taylor algorithms reduce the computational complexity thanks to the computationally efficient cost function. The complexity of the proposed FOC-DFT-Taylor is a bit heavier than that of the FOC-DFT-DPD and SOC-DFT-DPD algorithms, while it exhibits higher localization accuracy, which will be shown in Section 5.

Figure 3 shows the computational complexity of four different algorithms with the parameters setting as below:

K = 8

N_{1} = N_{2} = M / 2

T = 300

D_{x} = D_{y} = 500

, and the number of array elements is changed from 6 to 20. As we can see, compared to the FOC-SS-SDF-DPD and FOC-SS-Capon-DPD algorithm, the proposed FOC-DFT-DPD and FOC-DFT-Taylor algorithms reduce the computational complexity, as they are less sensitive to the number of array elements M. The complexity of the proposed FOC-DFT-Taylor algorithm is a bit heavier than that of FOC-DFT-DPD and SOC-DFT-DPD algorithms.

4.2. Advantages

Low Complexity: Compared to the FOC-SS-SDF-DPD and FOC-SS-Capon-DPD algorithms, the proposed algorithm reduces much of the computational complexity thanks to the computationally efficient cost function.
High DOF: Compared to the FOC-SS-SDF-DPD and FOC-SS-Capon-DPD algorithms, the proposed algorithm does not need SS technology for decoherence, which is caused by vectorizing the FOC matrix. Thus, the proposed algorithm can estimate more sources.
Suppress Gaussian Noise: The proposed algorithm takes full advantage of the characteristics of non-Gaussian signals, and the Gaussian noise is suppressed during the process of calculating the FOC matrix.
High Accuracy: When the sources do not fall on the preset grids, i.e., the off-grid error exists, the proposed algorithm can still estimate the position of sources accurately thanks to the Taylor compensation.
High-Resolution Capability: Compared to FOC-SS-SDF-DPD and FOC-SS-Capon-DPD methods, the proposed method has higher resolution capability.

5. Numerical Analysis

We conduct numerical simulations in this section to evaluate various algorithms and demonstrate the effectiveness of the proposed approach. The superiority of the proposed method is established through comparisons. Our evaluation includes assessing localization accuracy using root mean square error (RMSE) across different algorithms, which is defined as

RMSE = \frac{1}{K} \sum_{k = 1}^{K} \sqrt{\sum_{m_{c} = 1}^{M_{c}} {∥ {\hat{p}}_{k, m_{c}} - p_{k} ∥}^{2} / M_{c}}

(29)

here,

M_{c}

denotes the number of Monte Carlo runs, and

{\hat{p}}_{k, m_{c}}

is the estimated position of the kth source in the

m_{c}

th Monte Carlo experiment. The resolution capability of the algorithms is assessed by the definition of the effective estimate rate (EER), which is calculated by

EER = \frac{1}{M_{c}} \sum_{m_{c}}^{M_{c}} α_{m_{c}}

(30)

with

α_{m_{c}} = \{\begin{matrix} 1, if \frac{1}{K} \sum_{k = 1}^{K} ∥ {\hat{p}}_{k, m_{c}} - p_{k} ∥ < ε_{error} \\ 0, otherwise \end{matrix}

(31)

where

ε_{error}

represents the error threshold.

5.1. Effectiveness Analysis

Figure 4 shows scatter diagrams of the proposed FOC-DFT-Taylor algorithm, where the number of Monte Carlo runs is 300,

K = 2

, and

L = 4

. The location of sources and SAs are set as

u_{1} = {[- 900 m, 1000 m]}^{T}

u_{2} = {[- 300 m, 500 m]}^{T}

u_{3} = {[300 m, - 180 m]}^{T}

u_{4} = {[- 800 m, 900 m]}^{T}

p_{1} = {[100.5 m, 900.2 m]}^{T}

p_{2} = {[900.6 m, 200.7 m]}^{T}

. The results presented in Figure 4a demonstrate the successful estimation of sources by the FOC-DFT-DPD algorithm. Furthermore, analysis of Figure 4b reveals that the estimated positions closely align with the true positions following Taylor compensation, providing a compelling clue for the efficacy of the proposed FOC-DFT-Taylor algorithm.

5.2. RMSE Results

Figure 5 depicts the RMSE of four algorithms versus SNR when the number of Monte Carlo runs is 600,

M = 6

N_{1} = N_{2} = 3

T = 500

D_{x} = D_{y} = 600

L = 4

K = 2

u_{1} = {[- 900 m, - 1200 m]}^{T}

u_{2} = {[- 300 m, - 1100 m]}^{T}

u_{3} = {[300 m, - 1000 m]}^{T}

u_{4} = {[- 800 m, - 900 m]}^{T}

p_{1} = {[200.5 m, 623.5 m]}^{T}

p_{2} = {[920.5 m, - 174.5 m]}^{T}

, SNR is changed from

- 5 dB

30 dB

, and the gap of search grids is 1 m. Observe that the RMSE values of the four algorithms exhibit close proximity except for the SOC-DFT-DPD algorithm when SNR is below 10 dB, thanks to the property of non-Gaussian sources. The proposed FOC-DFT-Taylor algorithm demonstrates superior performance over the others at SNR levels exceeding 10 dB, showcasing its exceptional capabilities. Moreover, the RMSE associated with the proposed FOC-DFT-Taylor algorithm consistently remains below that of the FOC-DFT-DPD algorithm at the expense of performing Taylor compensation.

Figure 6 shows the RMSE of four algorithms versus T when the number of Monte Carlo runs is 600,

M = 6

N_{1} = N_{2} = 3

SNR = 15 dB

D_{x} = D_{y} = 600

L = 4

K = 2

u_{1} = {[- 900 m, - 1200 m]}^{T}

u_{2} = {[- 300 m, - 1100 m]}^{T}

u_{3} = {[300 m, - 1000 m]}^{T}

u_{4} = {[- 800 m, - 900 m]}^{T}

p_{1} = {[100.5 m, 815.5 m]}^{T}

p_{2} = {[970.5 m, 54.5 m]}^{T}

, the number of snapshots is changed from 200 to 900, and the gap of the search grids is 1 m. The localization performance of the FOC-DFT-DPD and SOC-DFT-DPD remains consistent, as they do not operate as super resolution algorithms and are constrained by M. Thanks to the property of non-Gaussian sources, the localization accuracy of the FOC-DFT-DPD is higher than that of SOC-DFT-DPD. Moreover, the RMSE values of the FOC-DFT-Taylor algorithm consistently outperform the other three algorithms due to its compensation method, showcasing its superior accuracy in localization.

Figure 7 shows the RMSE performance of different algorithms versus the number of array elements M over 500 Monte Carlo runs, when

N_{1} = N_{2} = M / 2

D_{x} = D_{y} = 600

L = 4

K = 2

u_{1} = {[- 900 m, - 1200 m]}^{T}

u_{2} = {[- 300 m, - 1100 m]}^{T}

u_{3} = {[300 m, - 1000 m]}^{T}

u_{4} = {[- 800 m, - 900 m]}^{T}

p_{1} = {[100.5 m, 815.5 m]}^{T}

p_{2} = {[970.5 m, 54.5 m]}^{T}

T = 500

SNR = 20 dB

, the number of array elements is changed from 4 to 18, and the gap of the search grids is 1 m. Observe that when the array size exceeds 8, the RMSEs of FOC-SS-SDF-DPD, FOC-SS-Capon-DPD, and FOC-DFT-DPD algorithms stabilize around 0.7 m due to their inability to address off-grid errors. In comparison, the RMSE values of the proposed FOC-DFT-Taylor algorithm consistently outperform the other algorithms thanks to its compensation method, demonstrating superior accuracy in localization. Unfortunately, the SOC-DFT-DPD performs poorly, as it fails to take full advantage of non-Gaussian sources. Furthermore, the RMSE values of FOC-DFT-DPD decrease as the array size increases, indicating that localization accuracy is constrained by the number of array elements.

5.3. EER Analysis

Figure 8 shows the EER of four algorithms versus SNR when the number of Monte Carlo runs is 600,

M = 8

N_{1} = N_{2} = 4

T = 500

D_{x} = D_{y} = 600

L = 4

K = 2

u_{1} = {[- 900 m, - 1200 m]}^{T}

u_{2} = {[- 300 m, - 1100 m]}^{T}

u_{3} = {[300 m, - 1000 m]}^{T}

u_{4} = {[- 800 m, - 900 m]}^{T}

p_{1} = {[100.5 m, 900.2 m]}^{T}

p_{2} = {[900.6 m, 200.7 m]}^{T}

, SNR is changed from 6 dB to 20 dB, the gap of search grids is 1 m, and the error threshold

ε_{error}

is set as 0.75 m. The EER achieved by the proposed FOC-DFT-Taylor algorithm surpasses 90%, contrasting with the EERs below 90% for the other four algorithms when SNR reaches a level of 14 dB. Unfortunately, the SOC-DFT-DPD performs poorly, as it fails to take full advantage of non-Gaussian sources. Notably, the proposed FOC-DFT-Taylor algorithm exhibits superior resolution capabilities compared to the remaining algorithms in the study.

Figure 9 shows the EER of four algorithms versus T, when the number of Monte Carlo runs is 600,

M = 6

N_{1} = N_{2} = 3

SNR = 20 dB

D_{x} = D_{y} = 600

L = 4

K = 2

u_{1} = {[- 900 m, - 1200 m]}^{T}

u_{2} = {[- 300 m, - 1100 m]}^{T}

u_{3} = {[300 m, - 1000 m]}^{T}

u_{4} = {[- 800 m, - 900 m]}^{T}

p_{1} = {[200.5 m, 623.5 m]}^{T}

p_{2} = {[920.5 m, - 174.5 m]}^{T}

, the number of snapshots is changed from 20 to 600, the gap of search grids is 1 m, and the error threshold

ε_{error}

is set as 1 m. The EER of the proposed FOC-DFT-Taylor algorithm shows a slight increase compared to the FOC-SS-SDF-DPD algorithm, yet with a significant reduction in complexity. Unfortunately, the SOC-DFT-DPD performs poorly, as it fails to take full advantage of non-Gaussian sources. Additionally, when compared to the FOC-SS-Capon-DPD and FOC-DFT-DPD algorithms, the EER of the FOC-DFT-Taylor algorithm is higher, underscoring its superior resolution capabilities.

6. Conclusions

In this article, the focus is on discussing the direct position determination (DPD) of non-Gaussian sources and proposing a new algorithm called the FOC-DFT-Taylor algorithm with multiple nested arrays (MNAs). This novel algorithm is designed to create an efficient DPD cost function that significantly reduces complexity when compared to conventional FOC-SS-SDF-DPD and FOC-SS-Capon-DPD algorithms. By leveraging the characteristics of non-Gaussian signals, this proposed algorithm offers a greater number of achievable degrees of freedom (DOFs). Additionally, due to the utilization of Taylor compensation, this new algorithm demonstrates improved localization accuracy and resolution capability over the traditional FOC-SS-SDF-DPD and FOC-SS-Capon-DPD algorithms.

Author Contributions

Conceptualization, H.H. and M.Y. (Meng Yang); methodology, Q.Y. and M.Y. (Mingyi You); software, X.S.; validation, H.H., M.Y. (Meng Yang), Q.Y., M.Y. (Mingyi You), X.S. and Y.S.; formal analysis, Y.S.; investigation, H.H.; resources, M.Y. (Meng Yang); data curation, Q.Y.; writing—original draft preparation, M.Y. (Mingyi You); writing—review and editing, H.H.; visualization, M.Y. (Meng Yang); supervision, Q.Y.; project administration, X.S.; funding acquisition, X.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National science foundation of China under Grants 62371225 and 62371227.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data used to support the findings of this study are included within the article.

Conflicts of Interest

The authors declare no conflicts of interest.

Abbreviations

The following abbreviations are used in this manuscript:

DPD	Direct Position Determination
MNA	Multiple Nested Array
DFT	Discrete Fourier Transform
AOA	Angle of Arrival
TDOA	Time Difference of Arrival
TOA	Time of Arrival
RSS	Received Signal Strength
ML	Maximum Likelihood
SNR	Signal-to-noise Ratio
SDF	Subspace Data Fusion
MUSIC	Multiple Signal Classification
EVD	Eigenvalue Decomposition
OFDM	Orthogonal Frequency Division Multiplexing
SOC	Second-order Cumulant
FOC	Fourth-order Cumulant
ULA	Uniform Linear Array
NA	Nested Array
DOF	Degree of Freedom
SS	Spatial Smoothing
RMSE	Root Mean Square Error
EER	Effective Estimate Rate

References

Coluccia, A.; Fascista, A. A review of advanced localization techniques for crowdsensing wireless sensor networks. Sensors 2019, 19, 988. [Google Scholar] [CrossRef]
Sibanyoni, S.V.; Ramotsoela, D.T.; Silva, B.J.; Hancke, G.P. A 2-D acoustic source localization system for drones in search and rescue missions. IEEE Sens. J. 2018, 19, 332–341. [Google Scholar] [CrossRef]
Zhao, K.; Zhao, T.; Zheng, Z.; Yu, C.; Ma, D.; Rabie, K.; Kharel, R. Optimization of time synchronization and algorithms with TDOA based indoor positioning technique for Internet of Things. Sensors 2020, 20, 6513. [Google Scholar] [CrossRef]
Zhang, R.; Cheng, L.; Wang, S.; Lou, Y.; Gao, Y.; Wu, W.; Ng, D.W.K. Integrated Sensing and Communication with Massive MIMO: A Unified Tensor Approach for Channel and Target Parameter Estimation. IEEE Trans. Wirel. Commun. 2024; early access. [Google Scholar]
Ma, Y.; Miao, C.; Long, W.; Zhang, R.; Chen, Q.; Zhang, J.; Wu, W. Time-Modulated Arrays in Scanning Mode Using Wideband Signals for Range-Doppler Estimation with Time-Frequency Filtering and Fusion. IEEE Trans. Aerosp. Electron. Syst. 2023, 60, 980–990. [Google Scholar] [CrossRef]
Wang, Y.; Ho, K. Unified near-field and far-field localization for AOA and hybrid AOA-TDOA positionings. IEEE Trans. Wirel. Commun. 2017, 17, 1242–1254. [Google Scholar] [CrossRef]
Wu, M.; Xiong, N.; Tan, L. Adaptive range-based target localization using diffusion Gauss–Newton method in industrial environments. IEEE Trans. Ind. Inform. 2019, 15, 5919–5930. [Google Scholar] [CrossRef]
Weiss, A.J.; Amar, A. Direct position determination of multiple radio signals. EURASIP J. Adv. Signal Process. 2005, 2005, 653549. [Google Scholar] [CrossRef]
Lo, T.K.; Leung, H.; Litva, J. Artificial neural network for AOA estimation in a multipath environment over the sea. IEEE J. Ocean. Eng. 1994, 19, 555–562. [Google Scholar] [CrossRef]
Yu, K.; Guo, Y.J. Statistical NLOS identification based on AOA, TOA, and signal strength. IEEE Trans. Veh. Technol. 2008, 58, 274–286. [Google Scholar] [CrossRef]
Wen, F.; Shi, J.; Zhang, Z. Joint 2D-DOD, 2D-DOA, and polarization angles estimation for bistatic EMVS-MIMO radar via PARAFAC analysis. IEEE Trans. Veh. Technol. 2019, 69, 1626–1638. [Google Scholar] [CrossRef]
Kulas, L. RSS-based DOA estimation using ESPAR antennas and interpolated radiation patterns. IEEE Antennas Wirel. Propag. Lett. 2017, 17, 25–28. [Google Scholar] [CrossRef]
Zeng, H.; Shi, S.; Yue, H.; Li, J. Coherent signal localization with multiple arrays: A successive perspective. Multidimens. Syst. Signal Process. 2023, 34, 321–340. [Google Scholar] [CrossRef]
Weiss, A.J. Direct position determination of narrowband radio frequency transmitters. IEEE Signal Process. Lett. 2004, 11, 513–516. [Google Scholar] [CrossRef]
Tirer, T.; Weiss, A.J. High resolution direct position determination of radio frequency sources. IEEE Signal Process. Lett. 2015, 23, 192–196. [Google Scholar] [CrossRef]
Yin, J.; Wu, Y.; Wang, D. Direct position determination of multiple noncircular sources with a moving array. Circuits Syst. Signal Process. 2017, 36, 4050–4076. [Google Scholar] [CrossRef]
Vallet, P.; Mestre, X.; Loubaton, P. Performance analysis of an improved MUSIC DOA estimator. IEEE Trans. Signal Process. 2015, 63, 6407–6422. [Google Scholar] [CrossRef]
Tzafri, L.; Weiss, A.J. Application of capon method to direct position determination. ICT Express 2016, 2, 5–9. [Google Scholar] [CrossRef]
Shu, T.; Li, L.; He, J. Near-field localization for non-circular sources in the presence of sensor phase uncertainties. IEEE Wirel. Commun. Lett. 2020, 10, 562–566. [Google Scholar] [CrossRef]
Wang, Q.; Wang, X.; Chen, H.; Zhu, X.; Liu, W.; Yan, W.; Wang, L. An effective localization method for mixed far-field and near-field strictly non-circular sources. Digit. Signal Process. 2019, 94, 125–136. [Google Scholar] [CrossRef]
Lin, Y.; Jin, S.; Matthaiou, M.; You, X. Channel estimation and user localization for IRS-assisted MIMO-OFDM systems. IEEE Trans. Wirel. Commun. 2021, 21, 2320–2335. [Google Scholar] [CrossRef]
Lu, Z.; Wang, J.; Ba, B.; Wang, D. A novel direct position determination algorithm for orthogonal frequency division multiplexing signals based on the time and angle of arrival. IEEE Access 2017, 5, 25312–25321. [Google Scholar] [CrossRef]
Xie, J.; Tao, H.; Rao, X.; Su, J. Efficient method of passive localization for near-field noncircular sources. IEEE Antennas Wirel. Propag. Lett. 2015, 14, 1223–1226. [Google Scholar] [CrossRef]
Zhang, Y.K.; Xu, H.Y.; Ba, B.; Wang, D.M.; Geng, W. Direct position determination of non-circular sources based on a Doppler-extended aperture with a moving coprime array. IEEE Access 2018, 6, 61014–61021. [Google Scholar] [CrossRef]
Shi, X.; Zhang, X.; Qian, Y.; Ahmad, M. Closed-form rectilinear emitters localization with multiple sensor arrays: Phase alignment and optimal weighted least-square method. IEEE Sens. J. 2023, 23, 7266–7278. [Google Scholar] [CrossRef]
Wang, H.; Zhao, J.; Qian, L. A study based on fourth-second order normalized cumulant in TDOA estimation. In Proceedings of the IEEE 10th International Conference on Signal Processing Proceedings, Beijing, China, 24–28 October 2010; pp. 215–218. [Google Scholar]
Hou, C.; Wang, D. The Fourth Order Cumulant-Based MUSIC Algorithm for Harmonic Retrieval in Colored Noise. In Proceedings of the 2009 WRI International Conference on Communications and Mobile Computing, Kunming, China, 6–8 January 2009; Volume 2, pp. 564–567. [Google Scholar]
He, Q.; Cheng, Z.; Wang, Z.; He, Z. Mixed far-field and near-field source separation and localization based on FOC matrix differencing. Digit. Signal Process. 2022, 131, 103753. [Google Scholar] [CrossRef]
Liu, C.; Ye, Z.; Zhang, Y. DOA estimation based on fourth-order cumulants with unknown mutual coupling. Signal Process. 2009, 89, 1839–1843. [Google Scholar] [CrossRef]
An, L.; Shu, W. Propagator method for DOA estimation using fourth-order cumulant. In Proceedings of the 2011 7th International Conference on Wireless Communications, Networking and Mobile Computing, Wuhan, China, 23–25 September 2011; pp. 1–4. [Google Scholar]
Alinezhad, P.; Seydnejad, S.R.; Abbasi-Moghadam, D. DOA estimation in conformal arrays based on the nested array principles. Digit. Signal Process. 2016, 50, 191–202. [Google Scholar] [CrossRef]
Peng, Z.; Ding, Y.; Ren, S.; Wu, H.; Wang, W. Coprime nested arrays for DOA estimation: Exploiting the nesting property of coprime array. IEEE Signal Process. Lett. 2021, 29, 444–448. [Google Scholar] [CrossRef]
Qin, G.; Zhang, Y.D.; Amin, M.G. DOA estimation exploiting moving dilated nested arrays. IEEE Signal Process. Lett. 2019, 26, 490–494. [Google Scholar] [CrossRef]
Han, K.; Nehorai, A. Improved source number detection and direction estimation with nested arrays and ULAs using jackknifing. IEEE Trans. Signal Process. 2013, 61, 6118–6128. [Google Scholar] [CrossRef]
Chen, Y.; Wang, W.; Wang, Z.; Xia, B. A source counting method using acoustic vector sensor based on sparse modeling of DOA histogram. IEEE Signal Process. Lett. 2018, 26, 69–73. [Google Scholar] [CrossRef]
Li, J.; He, Y.; Zhang, X.; Wu, Q. Simultaneous localization of multiple unknown emitters based on UAV monitoring big data. IEEE Trans. Ind. Inform. 2021, 17, 6303–6313. [Google Scholar] [CrossRef]
Ye, Z.; Zhang, Y. DOA estimation for non-Gaussian signals using fourth-order cumulants. IET Microwaves Antennas Propag. 2009, 3, 1069–1078. [Google Scholar] [CrossRef]
Dogan, M.C.; Mendel, J.M. Applications of cumulants to array processing. I. Aperture extension and array calibration. IEEE Trans. Signal Process. 1995, 43, 1200–1216. [Google Scholar] [CrossRef]
Yuen, N.; Friedlander, B. DOA estimation in multipath: An approach using fourth-order cumulants. IEEE Trans. Signal Process. 1997, 45, 1253–1263. [Google Scholar] [CrossRef]
Ye, C.; Chen, W.; Zhu, B.; Tang, L. Direction of arrival estimation of non-Gaussian signals for nested arrays: Applying fourth-order difference co-array and the successive method. ETRI J. 2021, 43, 869–880. [Google Scholar] [CrossRef]
Pal, P.; Vaidyanathan, P.P. Coprime sampling and the MUSIC algorithm. In Proceedings of the 2011 Digital Signal Processing and Signal Processing Education Meeting (DSP/SPE), Sedona, AZ, USA, 4–7 January 2011; pp. 289–294. [Google Scholar]
Pannert, W. Spatial smoothing for localized correlated sources–Its effect on different localization methods in the nearfield. Appl. Acoust. 2011, 72, 873–883. [Google Scholar] [CrossRef]

Figure 1. Geometry of multiple nested arrays localization.

Figure 2. Comparison of complexity versus L.

Figure 3. Comparison of complexity versus M.

Figure 4. Scatter diagrams of the proposed algorithm when

M = 4

N_{1} = N_{2} = 2

T = 300

D_{x} = D_{y} = 600

u_{1} = {[- 900 m, 1000 m]}^{T}

u_{2} = {[- 300 m, 500 m]}^{T}

u_{3} = {[300 m, - 180 m]}^{T}

u_{4} = {[- 800 m, 900 m]}^{T}

p_{1} = {[100.5 m, 900.2 m]}^{T}

p_{2} = {[900.6 m, 200.7 m]}^{T}

SNR = 5 dB

. (a) Initial estimation. (b) After Taylor compensation.

Figure 4. Scatter diagrams of the proposed algorithm when

M = 4

N_{1} = N_{2} = 2

T = 300

D_{x} = D_{y} = 600

u_{1} = {[- 900 m, 1000 m]}^{T}

u_{2} = {[- 300 m, 500 m]}^{T}

u_{3} = {[300 m, - 180 m]}^{T}

u_{4} = {[- 800 m, 900 m]}^{T}

p_{1} = {[100.5 m, 900.2 m]}^{T}

p_{2} = {[900.6 m, 200.7 m]}^{T}

SNR = 5 dB

. (a) Initial estimation. (b) After Taylor compensation.

Figure 5. Comparison of RMSE versus SNR.

Figure 6. Comparison of RMSE versus T.

Figure 7. Comparison of RMSE versus M.

Figure 8. Comparison of EER versus SNR.

Figure 9. Comparison of EER versus T.

Table 1. Comparison of computational complexity.

Methods	Computational Complexity
SOC-DFT-DPD	$O (L T M^{2} + L M_{v e c} D_{x} D_{y})$
FOC-DFT-DPD	$O (L T M^{4} + L J D_{x} D_{y})$
FOC-DFT-Taylor	$O (L T M^{4} + L J D_{x} D_{y} + 27 K^{3} + 9 K^{2} + 9 K^{2} L J + 5 K L J)$
FOC-SS-SDF-DPD [22]	$O (L T M^{4} + 2 L N_{a}^{3} + L D_{x} D_{y} N_{a} (N_{a} - K))$
FOC-SS-Capon-DPD [18]	$O (L T M^{4} + L N_{a}^{3} + L D_{x} D_{y} (N_{a}^{2} + N_{a}))$

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.

© 2024 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/).

Share and Cite

MDPI and ACS Style

Hu, H.; Yang, M.; Yuan, Q.; You, M.; Shi, X.; Sun, Y. Direct Position Determination of Non-Gaussian Sources for Multiple Nested Arrays: Discrete Fourier Transform and Taylor Compensation Algorithm. Sensors 2024, 24, 3801. https://doi.org/10.3390/s24123801

AMA Style

Hu H, Yang M, Yuan Q, You M, Shi X, Sun Y. Direct Position Determination of Non-Gaussian Sources for Multiple Nested Arrays: Discrete Fourier Transform and Taylor Compensation Algorithm. Sensors. 2024; 24(12):3801. https://doi.org/10.3390/s24123801

Chicago/Turabian Style

Hu, Hao, Meng Yang, Qi Yuan, Mingyi You, Xinlei Shi, and Yuxin Sun. 2024. "Direct Position Determination of Non-Gaussian Sources for Multiple Nested Arrays: Discrete Fourier Transform and Taylor Compensation Algorithm" Sensors 24, no. 12: 3801. https://doi.org/10.3390/s24123801

APA Style

Hu, H., Yang, M., Yuan, Q., You, M., Shi, X., & Sun, Y. (2024). Direct Position Determination of Non-Gaussian Sources for Multiple Nested Arrays: Discrete Fourier Transform and Taylor Compensation Algorithm. Sensors, 24(12), 3801. https://doi.org/10.3390/s24123801

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

Article Menu

Direct Position Determination of Non-Gaussian Sources for Multiple Nested Arrays: Discrete Fourier Transform and Taylor Compensation Algorithm

Abstract

1. Introduction

2. Model Formulation

3. Proposed Algorithm

4. Performance Discussion

4.1. Complexity

4.2. Advantages

5. Numerical Analysis

5.1. Effectiveness Analysis

5.2. RMSE Results

5.3. EER Analysis

6. Conclusions

Author Contributions

Funding

Institutional Review Board Statement

Informed Consent Statement

Data Availability Statement

Conflicts of Interest

Abbreviations

References

Share and Cite

Article Metrics

Article Access Statistics

Further Information

Guidelines

MDPI Initiatives

Follow MDPI