[go: up one dir, main page]
More Web Proxy on the site http://driver.im/
Next Article in Journal
A Generative Adversarial Network to Synthesize 3D Magnetohydrodynamic Distortions for Electrocardiogram Analyses Applied to Cardiac Magnetic Resonance Imaging
Previous Article in Journal
Smart Farm Security by Combining IoT Sensor Network and Virtualized Mycelium Network
You seem to have javascript disabled. Please note that many of the page functionalities won't work as expected without javascript enabled.
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Sensors
Volume 23
Issue 21
10.3390/s23218690
sensors-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

High-Precision DOA Estimation Based on Synthetic Aperture and Sparse Reconstruction

by
Yang Fang
1,*,
Xiaolong Wei
2 and
Jianjun Ma
3
1
State Key Laboratory for Strength and Vibration of Mechanical Structures, Shaanxi Engineering Research Centre of NDT and Structural Integrity Evaluation, School of Aerospace Engineering, Xi’an Jiaotong University, Xi’an 710049, China
2
Science and Technology on Plasma Dynamics Laboratory, Air Force Engineering University, Xi’an 710000, China
3
National Laboratory of Radar Signal Processing, Xidian University, Xi’an 710071, China
*
Author to whom correspondence should be addressed.
Sensors 2023, 23(21), 8690; https://doi.org/10.3390/s23218690
Submission received: 9 September 2023 / Revised: 17 October 2023 / Accepted: 18 October 2023 / Published: 24 October 2023
(This article belongs to the Section Radar Sensors)
Browse Figures
Figure 1
<p>Radar target space geometry.</p> "> Figure 2
<p>Spatial synthetic aperture geometry.</p> "> Figure 3
<p>Diagram representing the flow of the proposed algorithm.</p> "> Figure 4
<p>Initial rough estimation result. (The area pointed by the red arrow in the figure is an enlargement of the target area.)</p> "> Figure 5
<p>OMP estimation result of refined dictionary.</p> "> Figure 6
<p>OMP-RELAX estimation result.</p> "> Figure 7
<p>SBL estimation result.</p> "> Figure 8
<p>Comparison of angle estimation for flickering and non-flickering targets.</p> "> Figure 9
<p>Reconstruction error of a non-flickering target.</p> "> Figure 10
<p>RMSE results across different heading angles.</p> "> Figure 11
<p>Mean correlation coefficients of the dictionary across varied spatial sampling points.</p> "> Figure 12
<p>Angle measurement precision for targets across different distances.</p> "> Figure 13
<p>Angle estimation results of the real data by the proposed method.</p> ">
Versions Notes

Abstract

:
The direction-of-arrival (DOA) estimation is predominantly influenced by the antenna’s aperture size. However, space constraints on flight platforms often necessitate the use of antennas with smaller apertures and fewer array elements. This inevitably imposes limitations on the DOA estimation’s resolution and degrees of freedom. To address these precision constraints, we introduce an accurate DOA estimation method based on spatial synthetic aperture model. This method adopts a two-stage strategy to ensure both efficiency and precision in DOA estimation. Initially, the orthogonal matching pursuit (OMP) reconstruction algorithm processes the original aperture data, providing a rough estimate of target angles that guides the aircraft’s flight direction. Subsequently, the early estimations merge with the aircraft’s motion space samples, forming equivalent spatially synthesized array samples. The refined angle estimation then employs the OMP-RELAX algorithm. Moreover, with the off-grid issue in mind, we devise an estimation method integrating Bayesian parameter estimation with dictionary sequence refinement. The proposed technique harnesses the spatial synthetic aperture for pinpoint estimation, effectively addressing the challenges of atomic orthogonality and angular off-grid on estimation accuracy. Importantly, the efficiency of deploying sparse reconstruction for angle estimation is bolstered by our phased strategy, eliminating the necessity for fine grid analysis across the entire observation scene. Moreover, the poor estimation accuracy caused by coherent source targets and angular-flickering targets is improved by sparse reconstruction. Through simulation and experiment, we affirm the proposed method’s efficacy in angle estimation. The results indicate that target angle estimation errors are limited to within 1°. Furthermore, we assess the impact of variables such as target state, heading angle, spatial sampling points, and target distance on the estimation accuracy of our method, showcasing its resilience and adaptability.

1. Introduction

Enhancing the Direction of Arrival (DOA) estimation is one of the ways to improve radar detection accuracy [1,2]. However, the spatial constraints on the platform mean that only small aperture radars can be accommodated. With radar beam-width inversely related to antenna aperture, achieving precise angles via traditional DOA estimation techniques poses challenges [3,4,5]. This tension between the demand for precise DOA estimation—especially for coherent source DOA—and radar resource constraints underscores the limitations of current DOA methods. Angle glint is another obstacle: echo power may diminish at certain angles or times due to adversarial radar interference. This makes angle estimation particularly challenging. Thus, there is an urgent need to investigate high-resolution DOA estimation methods adaptable to radar constraints. Synthetic aperture techniques for DOA estimation, which may offer improved accuracy, present an intriguing area of research.
Numerous DOA estimation methods have emerged over the past few decades. Spatial spectrum estimation techniques such as MUSIC [6,7,8] and ESPRIT [9,10] are commonly utilized for incoherent source angle estimation. These methods depend on statistical analyses, harnessing covariance from multiple snapshots and employing eigenvalue analysis to realize super-resolution DOA estimation. However, for coherent sources, steer vectors often converge to the same eigenvalue subspace, rendering this method less effective [11,12,13]. Additionally, these methods suffer from low efficiency due to eigenvalue decomposition, making them ill-suited for time-sensitive DOA estimation [14,15]. Alternative techniques like the iterative adaptive (IAA) method [16] have been suggested. These methods utilize weighted least squares (WLS) for super-resolution spectral estimation and the Toeplitz-FFT method to improve computational efficiency [17,18]. However, even with Toeplitz-FFT integration, they might not adequately meet the efficiency standards for radar [19]. For coherent source DOA estimation, temporal smoothing techniques combined with MUSIC and ESPRIT have been proposed [20,21]. In reference [22], a dual-polarization synthetic sparse nested array consisting of three subarrays is proposed. Based on the sparse nested array, a DOA estimation method is proposed to improve the estimation accuracy, which is first rough estimation and then precise estimation. Recently, methods harnessing compressive sensing have been advanced for both coherent and incoherent DOA estimation [23,24]. These methods capitalize on sparsity in the angle domain, with [25] integrating sequential snapshots with sparsity priors for Bayesian inference-based sequential DOA estimation, and [26,27] applying multiple snapshots to devise the MMV-OMP sparsity reconstruction method for coherent DOA estimation. Although these methods can provide super-resolution estimation, they may experience performance degradation in radar applications due to the small antenna aperture.
To address these challenges, this paper introduces a two-stage DOA estimation technique. Initially, the original aperture facilitates a rough angle estimation, which then determines the refinement measurement grid and sets the flight trajectory. Subsequently, merging the initial target angles with spatial samples from the moving target, these samples are treated as equivalent spatial synthetic array samples. A measurement dictionary, rooted in spatial samples, is established, and sparse reconstruction is applied for angle estimation. To bolster DOA estimation accuracy, a sequential estimation method, amalgamating Bayesian proportion parameter estimation with dictionary refinement, is designed. The method presented holds three distinct advantages:
(1)
Enhanced Precision: The spatial domain’s synthetic array, with its expanded aperture and heightened resolution, combined with sequential estimation that integrates Bayesian proportion parameter estimation and dictionary refinement, delivers a more precise DOA.
(2)
Augmented Efficiency: Rough angle estimation results serve as the foundation for subsequent estimations, significantly diminishing the algorithm’s computational complexity. And the DOA estimation process uses the greedy method of sparsity reconstruction, catering to the need for rapid estimations.
(3)
Mitigation of angle glint issues: The spatial synthesis array sample and target angle reconstruction can be considered equivalent to a spatial matched filter. Even though the target’s power may be weak over a certain period, DOA estimation can still be accomplished using spatial accumulation over an extended duration.
The remainder of this paper is organized as follows: Section 2 outlines the signal model and method, while Section 3 discusses comparative tests and performance analysis. Lastly, conclusions on the proposed method are summarized in the final section.

2. Method Description

2.1. Synthesis Array DOA Estimation Signal Model

We assume that our radar system is tasked with estimating three coherent sources, as depicted in Figure 1. To facilitate this, a uniform circular array is incorporated into the radar system. The origin point signifies the initial measurement position, and the axis direction is illustrated in Figure 1. The azimuth angle φ is the angle formed between the target and the X-axis. Concurrently, the pitch angle θ represents the angle between the target and the XOY plane. Targets are configured in an equilateral triangle and are labeled as P 1 , P 2 , and P 3 , with their respective coordinates and angles designated as P 1 ( x 1 , y 1 , z 1 ) , P 2 ( x 2 , y 2 , z 2 ) , P 3 ( x 3 , y 3 , z 3 ) and P 1 ( φ 1 , θ 1 ) , P 2 ( φ 2 , θ 2 ) , P 3 ( φ 3 , θ 3 ) . respectively.
The coordinates of the three targets relate as follows:
[ x 2 − x 1   y 2 − y 1 ] T = Φ [ x 3 − x 1   y 3 − y 1 ] T [ x 3 − x 2   y 3 − y 2 ] T = Φ [ x 3 − x 1   y 3 − y 1 ] T ,
where Φ is the rotation matrix:
Φ = cos ( 60 )   − sin ( 60 ) sin ( 60 )   cos ( 60 ) .
Angles can be deduced from the target’s coordinates as
θ i = arctan ( x i y i ) ,   φ i = arctan ( z i x i 2 + y i 2 ) ,   i = 1 , 2 , 3 .
Let us denote the target steering vector as a ( φ i , θ i ) and the signal envelope as s ( t ) . If we assume noise n ( t ) in every snapshot to be independent and identically distributed Gaussian white noise with a noise power σ 2 , the resultant signal can be expressed as
y ( t ) = ∑ i = 1 I A i a ( φ i , θ i ) s i ( t ) + n ( t ) ,
where A i is the reflection coefficient of the i t h target, and y(t) denotes all the reflected signals in the observation scene received by the antennas.
To enlarge the array aperture, the movement of the platform is employed to design the spatial synthesis array, depicted in Figure 2.
We suppose the platform’s heading angle is α , its velocity is v , and the spatial sampling time is T s p . If the target’s initial azimuth is φ 0 and its coordinates are given by p ( x 0 , y 0 , z 0 ) , the vertical distance from the target to the airline can be expressed as
d = y 0 − tan ( α − φ 0 ) x 0 1 + tan 2 ( α − φ 0 ) .
The center point of the scene is defined as the projection of the target in the flight direction, so distance w between scene center o ′ and the platform is
w = d tan ( α − φ ) .
Consequently, the distance from the i t h spatial sample to the scene center is
l i = ( w − i T s p v ) 2 + d 2 + h 2 ,
where h is the elevation of the platform from the ground. The phase difference from the i t h spatial sample of the platform to the original spatial sample is
Δ i = 2 π λ ( l i − l 0 ) .
Under these conditions, the comprehensive spatial sample steering vector becomes
a s p = [ exp ( j Δ 0 ) , exp ( j Δ 1 ) … exp ( j Δ M ) ] T .
After accumulating M spatial samples, the spatial echo across arrays can be described as
y s p = a s p ⊙ [ s ( 0 • T s p ) , s ( 1 • T s p ) … s ( M • T s p ) ] T ,
y s p = a s p ⊙ [ s 0 ⋅ T s p , s 1 ⋅ T s p , … , s M ⋅ T s p ]
where a 1 , a 2 … a M T ⊙ b 1 , b 2 … b M T = a 1 b 1 , a 2 b 2 , … , a M b M T [ a 1 , a 2 , … a M ] T ⊙ [ b 1 , b 2 , … b M ] T = [ a 1 b 1 , a 2 b 2 , … a M b M ] T , the radar equivalent aperture is L = v M T s p . With an adequate number of samples, the synthesized beam width is sufficiently narrow to discern target angles.

2.2. Rough DOA Estimation

Initially, we employ the maximum coherence method to provide a rough estimation of target angles using the original radar array. The angle dictionary is formulated as
A ( φ , θ ) = [ a ( φ j , θ i ) ] , i = 1 , 2 …. I , j = 1 , 2 …. J ,
where A consists of spatial steering vectors, and the angle grid is partitioned according to the beam width. Given that multiple targets fall within a single beam, the maximum coherence between several snapshots and the angle dictionary can be harnessed to ascertain the rough angle. This estimation can be represented by the following equation:
[ φ ^ r o u g h , θ ^ r o u g h ] = m e a n ( max ( φ t , θ t ) s ( t ) T A ( φ , θ ) ) , t = 0 , T s , 2 T s …. ,
where T s denotes the snapshot sampling time. Each snapshot provides a rough angle estimation through the metric of maximum coherence. At this time, the angle grid in the dictionary atom can be divided according to the beam width of the array. The estimation of the rough angle can be used to plan the flight direction of the aircraft. Due to the limitation of the aperture, when the angle interval of the target is too small, it is unable to distinguish multiple targets, so the rough estimation only needs to obtain the approximate angle of the target. By averaging across numerous snapshots, a more refined rough estimation is acquired. Subsequently, to differentiate between targets, the spatial sample synthesis array is leveraged to expand the radar’s aperture, facilitating a more accurate angle estimation.

2.3. Refined DOA Estimation

In order to both minimize the count of spatial samples and enhance the efficiency of angle estimation, a sparse reconstruction framework is proposed. This serves as an effective method for super-resolution estimation of target angles. Let us denote the precise angle grid as Δ φ p r e c i o u s and Δ θ p r e c i o u s . (It is worth noting that the angle division here builds upon the preceding step, utilizing the roughly estimated angle from that phase as the foundation for this more refined estimation.) The complete spatial domain’s precise steering vector, denoted as A p r e c i o u s , is derived; refer to Equations (5)–(9). Consequently, our signal model can be expressed as
y s p = A p r e c i o u s x + n ,
where A p r e c i o u s = [ a ( Δ φ p r e c i o u s , i , Δ θ p r e c i o u s , j ) ] , and x stands for the echo intensity vector associated with a given angle. x is sparse, signifying that the majority of its elements are zero, with a much smaller number being non-zero. As such, x can be estimated via the following optimization function:
min ‖ y s p − A p r e c i o u s x ‖ 2 2 , s . t . ‖ x ‖ 0 < ε .
This presents a non-convex challenge. However, it can be addressed through greedy algorithms like MP (Matching Pursuit) and OMP (Orthogonal Matching Pursuit). To incorporate multiple snapshot data, the MMV-OMP (Multiple Measurement Vector-OMP) approach can be employed to resolve Equation (14).

2.4. Dynamic Off-Grid Estimation Based on Bayesian Theory

The method previously discussed allows for more accurate estimation of adjacent target angles. However, with spatial sampling and estimation, two pressing issues remain:
(1)
Off-Grid Estimation: Our dictionary dictates that angles can only be estimated based on its specific grid. Unfortunately, actual angles often deviate from this grid, leading to potential errors.
(2)
Non-Orthogonal Atoms: We utilize the OMP method for angle estimation, but it is crucial to recognize that the atoms in our angle dictionary are not orthogonal. Given that amplitude estimation relies on orthogonal projection, the estimated angle typically falls between two adjacent atoms.
Beyond post-sampling estimation, a dynamic estimation process could be crafted to conduct spatial sampling and angle estimation concurrently. This real-time process can utilize error information as it emerges, enhancing estimation precision. To tackle these issues and refine precision, we propose a DOA estimation method based on sequential off-grid Bayesian reference, which refines the measurement dictionary.
To address non-orthogonal atoms, the OMP-RELAX method offers a potential solution to boost estimation performance [28]. The RELAX method, a super-resolution spectral estimation technique, uses iterative re-projection to diminish atom correlation introduced by non-orthogonal atom side-lobes. Thus, RELAX is deployed post each OMP estimation.
For the off-grid challenge, we present a dynamic Bayesian off-grid method, integrated with sequential sampling. As spatial sampling grows, so does angle estimation resolution. Given there are N samples during k t h estimation, re-estimating the angle post each N T sample becomes feasible. Focusing on the azimuth angle synthesis aperture, the azimuth beam-width during k t h estimation can be gauged as
φ b w = 1 cos ( φ 0 ) 50.8 λ N v T s p ,
where φ 0 defines the target azimuth, and λ signifies the wavelength. From this, we can derive the beam width during the ( k + 1 ) t h estimation and delineate the angle grid of the dictionary for subsequent estimations using the beam width. The angle can be re-estimated after each N T sample, and the beam width can be calculated during the ( k + 1 ) t h estimation as
Δ φ k + 1 = 1 R • 1 cos ( φ 0 ) 50.8 λ ( N + N T ) v T s p .
∆ φ k + 1 = 1 R ⋅ 1 cos ⁡ φ 0 50.8 λ ( N + N T ) v T s p
Considering the k t h estimation, the measurement dictionary can be refined according to the resolution near the angle x ^ k . Combining the angle resolution from the aforementioned equation, the real angle lies within
x ∈ [ x ^ k − Δ φ k 2 , x ˜ k + 1 + Δ φ k + 1 2 ] .
We posit x ^ k < x ˜ k + 1 and introduce proportion parameters α ( 0 ≤ α ≤ 1 ) to describe the real angle as
x = α • ( x ^ k − Δ φ k 2 ) + ( 1 − α ) • ( x ˜ k + 1 + Δ φ k + 1 2 ) .
x = α ⋅ x ^ k − ∆ φ k 2 + ( 1 − α ) ⋅ ( x ~ k + 1 + ∆ φ k + 1 2 )
It is pertinent to note that the real angle measurement atom evolves as a function of the angle from Equations (5)–(9). Using the first-order Taylor approximation of Equation (5), we can express
A ( x ) = A ( x ˜ k + 1 ) + A ′ ( x ˜ k + 1 ) ( x − x ˜ k + 1 ) + O ( x − x ˜ k + 1 ) .
Integrating the above into the derived Equation (18) and overlooking higher-order terms, the atom becomes
A ( x ) ≈ [ A ( x ˜ k + 1 ) + Δ φ k + 1 2 A ′ ( x ˜ k + 1 ) ] + α • [ A ′ ( x ˜ k + 1 ) ( x ^ k − x ˜ k + 1 − Δ φ k + 1 2 − Δ φ k 2 ) ] .
A x ≈ A x ~ k + 1 + ∆ φ k + 1 2 A ′ x ~ k + 1 + α A ′ x ~ k + 1 ( x ^ k − x ~ k + 1 − ∆ φ k + 1 2 − ∆ φ k 2 )
Estimating α enables pinpointing the exact angle. The subsequent section provides the estimation of α within the Bayesian context.
The relationship between the measurement signal and real parameters in the ( k + 1 ) t h spatial sampling can be expressed as
y k + 1 = A ( x ) x x + n .
Substituting Equation (20) into Equation (21) results in an accurate measurement model based on estimation parameters, as follows:
y k + 1 = [ P ( x ˜ k + 1 ) + α • Q ( x ˜ k + 1 , x ^ k ) ] x x + n ,
y k + 1 = P x ~ k + 1 + α ⋅ Q ( x ~ k + 1 , x ^ k ) x x + n
where
P ( x ˜ k + 1 ) = [ A ( x ˜ k + 1 ) + Δ φ k + 1 2 A ′ ( x ˜ k + 1 ) ] Q ( x ˜ k + 1 , x ^ k ) = [ A ′ ( x ˜ k + 1 ) ( x ^ k − x ˜ k + 1 − Δ φ k + 1 2 − Δ φ k 2 ) ] .
Assuming the noise is Gaussian white noise, the probability distribution is n ~ N ( 0 , σ 2 n ) , and the amplitude distribution is x x ~ N ( x x ˜ k + 1 , σ 2 x ) , where x x ˜ k + 1 is the estimated amplitude after OMP-RELAX refined estimation. Consequently, the posterior distribution of amplitude x x can be represented by a complex Gaussian using Bayesian formulation as
p ( x x | y k + 1 ; σ n , σ x ) = p ( y k + 1 | x x ; σ n , σ x ) p ( x x | σ x ) ∫ p ( y k + 1 | x x ; σ n , σ x ) p ( x x | σ x ) d x x .
The complex Gaussian distribution is p ( x x | y k + 1 , σ n , σ x ) ~ C N ( μ k + 1 , Σ k + 1 ) , where
μ k + 1 = Σ k + 1 ( β A ( x ) H y k + 1 + Σ x − 1 x ˜ k + 1 ) Σ k + 1 = ( β A ( x ) H A ( x ) + Σ x − 1 ) − 1 .
Referencing [29,30], the proportion parameter α can be estimated using the maximum posterior with the EM iteration method. In the E-step, parameter x x is treated as a hidden variable, and given the fixed parameters y k + 1 , σ n , and σ x in the i t h iteration, the E-step expectation can be presented as
E x x | y k + 1 ; σ n ( i ) , σ x ( i ) [ ln ( p ( y k + 1 , x x   |   σ n , σ x ) ) ] ,
where p ( y k + 1 , x x | σ n , σ x ) can be deduced using Bayesian formulation as
p ( y k + 1 , x x   |   σ n , σ x ) = p ( y k + 1   |   x x ; σ n , σ x ) p ( x x   |   σ x ) .
In this case, Equation (26) can be calculated. The next step, the M-step, involves estimating the parameters by maximizing Expectation (26):
σ n i + 1 , σ x i + 1 = arg max σ n , σ x E x x   |   y k + 1 ; σ n ( i ) , σ x ( i ) [ ln ( p ( y k + 1 , x x   |   σ n , σ x ) ) ] .
Following a similar approach as in [31], the proportion parameter α can be estimated concurrently with σ x i + 1 , and the amplitude estimation can use the expectation of the posterior distribution, that is,
α ^ l i + 1 = ℜ ( μ l * q l H ( y k + 1 − P μ ) ) − ℜ ( q l H P γ i ) p l 2 2 μ l + γ i , i p l 2 2 x ^ x i + 1 = μ i + 1 .
After several iterations, proportion parameter α converges near the real value, allowing the estimation in the ( k + 1 ) t h spatial sampling to be expressed as
x ^ k + 1 = α ^ k + 1 x ^ k + ( 1 − α ^ k + 1 ) x ˜ k + 1 .

2.5. Overall Procedure of the Proposed Method

Figure 3 illustrates the comprehensive flow of the proposed method, which can be broken down into four main steps:
Step 1: Initial Angle Estimation using a Single Snapshot Signal
  • Initially, using the beam width, the angle interval is coarsely divided.
  • The dictionary matrix is then assembled using Equation (11).
  • For a singular snapshot dataset, the OMP algorithm is used to solve Equation (12) for singular peak reconstruction, which provides the primary rough estimation of the target angle.
Step 2: Angle Estimation via the Synthesized Array Signal
  • Based on the previously estimated result [ φ ^ r o u g h , θ ^ r o u g h ] , the aircraft’s heading angle α is planned.
  • The number Q of sampling points is predetermined.
  • Leveraging multiple snapshot datasets, a locally refined angle comprehensive guide vector matrix A p r e c i o u s is fashioned using Equations (5)–(9).
  • The MMV-OMP algorithm is then deployed to provide a rough angle estimation.
Step 3: Correction for Non-Orthogonal Dictionary Atomic Errors
  • The OMP-RELAX algorithm rectifies errors stemming from non-orthogonal atoms Φ .
  • The angle estimation derived from the previous stage serves as the input.
  • The angle undergoes refinement by Φ r e f i n e = Φ − δ , Φ − δ + ( Δ φ , Δ θ ) η , ⋯ , Φ + δ , where Δ φ and Δ θ correspond to the azimuth and elevation angle interval divisions during the initial estimation, respectively.
  • The refinement parameter η = ( η φ , η θ ) is determined based on the desired resolution precision.
  • The angle constraint δ = ( δ φ , δ θ ) is set in accordance with the current angle interval division, typically half of the present angle interval.
  • Utilizing this, the revised dictionary matrix A ( Φ r e f i n e ) is acquired. Subsequently, the OMP algorithm procures K revised target angle estimations.
Step 4: Off-Grid Rectification
  • Drawing from the previous step’s estimated value and coupled with spatial sequence sampling, the sequential Bayesian dynamic estimation technique is employed to address positional deviations due to off-grid anomalies.
  • The initial values ( α 0 , σ n 0 , σ x 0 ) of the undisclosed parameter are set, with the iteration number fixed at K.
  • Using Equation (28), the matrices Σ k + 1 and μ k + 1 are refreshed.
  • Equations (29) and (30) facilitate updates to ( σ n 0 , σ x 0 ) and α 0 , respectively.
  • Through iterative solutions, the final precise target angle estimation is derived.

3. Simulation and Experimental Verification and Analysis

3.1. Validation of Effectiveness

To verify the effectiveness of the method proposed in this paper, simulation data were used. The radar and target scene settings are depicted in Figure 2. The circular array is made up of seven omnidirectional array elements with an array radius of 8 cm. The radar emits a linear frequency-modulated signal at a carrier frequency of 3.75 GHz. The carrier’s altitude is 1.5 km, the distance from the radar to the scene’s center is 5 km, and the spatial sampling interval measures 0.04 m. The angular coordinates of the three targets in the scene are P 1 ( 20.1 ° , − 17.4 ° ) , P 2 ( 23.7 ° , − 16.2 ° ) , and P 3 ( 26.1 ° , − 17.6 ° ) .
Given the proximity of the targets, their azimuth and elevation angles are notably similar. During the initial rough estimation, both azimuth and elevation intervals are set at 1°, with search ranges designated at (0, 90°) for the azimuth and at (−90°, 0°) for elevation angles. Figure 4 presents the initial rough estimation result, highlighting an estimated azimuth angle of 24° and an elevation angle of −17°. The area pointed by the red arrow in the figure is an enlargement of the target area. Due to the pronounced correlation among adjacent atoms in the dictionary (with a correlation coefficient nearing one), the search results in a singular target point. Since rough estimation involves azimuth and elevation angles, a two-dimensional angle estimation is provided. In the first step of the rough estimation, the three targets set are located at the angles corresponding to the same atom. This angle is used to offer a rough location of the three adjacent targets.
Using the initial rough estimation ( φ r o u g h = 24 ° ), we plan the trajectory. The half-power beam width corresponding to the synthetic aperture is φ h f = 0.015 ° . Consequently, the angular search range for the azimuth angle becomes ( φ r o u g h − φ h f ) : Δ φ r o u g h / 4 : ( φ r o u g h + φ h f ) , where Δ φ r o u g h represents the azimuth angle search grid during the initial estimation. Since the target elevation angle remains relatively consistent throughout the search phase, the search range for the elevation angle is set at ( θ r o u g h − 1 ° ) : Δ θ r o u g h / 4 : ( θ r o u g h + 1 ° ) . We then construct a dictionary matrix corresponding to this refined angle based on the established azimuth and elevation angle grid. Figure 5 displays the refined estimation result. It is evident that all three targets have an estimated elevation angle of −17°, and their estimated azimuth angles are 19°, 23°, and 25.5°, respectively. While the exact azimuth angles of the targets are not pinpointed, the possibility of estimating the three targets within the same beam range became feasible after multiple spatial samplings and refined angle divisions. As there is a greater focus on azimuth angle information during the search state, the rough estimation of the elevation angle is used for a refined azimuth angle estimation. Compared to the first step’s rough estimation, it is clear that the three targets can be distinguished, but the angle values are still not accurate.
To further enhance estimation precision, we employ the RELAX algorithm to fine-tune the target azimuth estimation further. Following the RELAX algorithm’s methodology, when estimating the k t h target, the center of the previous estimation result φ k is utilized to construct the refined angle search grid ( φ k − Δ φ r o u g h / 4 ) : Δ φ r o u g h : ( φ k + Δ φ r o u g h / 4 ) for the k t h target. In our study, we define Δ φ r o u g h = 0.1 ° , with the resulting estimations illustrated in Figure 6. The azimuth angles for the three targets are determined to be 19.5°, 23.8°, and 26.6°, which signifies a marked improvement in estimation accuracy compared to the results of the OMP algorithm. Figure 6, based on the elevation angle estimation result of Figure 4 and the azimuth angle estimation result of Figure 5, re-divides the azimuth grid. The results indicate that after processing with the OMP-RELAX algorithm, the estimation accuracy is further improved.
Acknowledging that a target might not align with the precise grid during sparse reconstruction, this study introduces a sequential SBL technique to address the off-grid issue. The outcome of this method is presented in Figure 7. In comparison to the actual target value; the error margin is a mere 0.1°. Collectively, these experiments underscore that our proposed method boasts commendable estimation accuracy. Figure 7, based on the angle estimation of Figure 6, displays the estimation results processed by the SBL algorithm. It is evident from the figure that the estimation result is almost consistent with the target angle.

3.2. Robustness Analysis

In this section, we delve into the robustness of the proposed method by assessing its efficacy under varying conditions such as target state (flickering and non-flickering), signal-to-noise ratio (SNR), heading angle, spatial sampling, and target distance. The angle measurement performance for flickering and non-flickering targets is gauged via 1000 Monte Carlo simulations, utilizing the root mean square error (RMSE) of both target azimuth angle and amplitude to measure accuracy. RMSE is defined as
R M S E D O A = 1 M K ∑ k = 1 K ∑ m = 1 M x k − x ^ k , m 2 2 R M S E A m p = 1 M K ∑ k = 1 K ∑ m = 1 M Θ k − Θ ^ k , m 2 2 ,
where M denotes the number of Monte Carlo trials and K represents the number of targets.

3.2.1. Algorithm’s Adaptability to Various Target States

Figure 8 presents a comparative analysis of angle estimation results for diverse target types using the proposed method and other algorithms over a range of SNRs (0~20 dB). In the experiment, evaluations are conducted at intervals of 20 points, totaling 50 evaluations. Observations indicate that our method exhibits superior robustness, especially at lower SNRs, in comparison to alternative algorithms. Although flickering targets can be estimated; they demand a higher SNR than non-flickering targets. This is attributed to the inability of flickering targets to fully accumulate energy during spatial sampling, which subsequently leads to a loss of SNR at specific spatial sampling points. Moreover, the method’s precision in measuring both target types is notable, stemming from three pivotal reasons:
  • The OMP-RELAX process minimizes non-orthogonal atom interference.
  • The introduced algorithm mitigates the off-grid problem.
  • Spatial sampling dynamically refines the dictionary, superseding the original fixed angle grid.
The amplitude reconstruction error is depicted in Figure 9. We focus our amplitude estimation exclusively on non-flickering targets for two primary reasons. Firstly, the estimation performance for non-flickering targets surpasses that of flickering targets. Secondly, and more crucially, the inherent unpredictability of target flickering hampers the acquisition of precise amplitude estimation results. The proposed method demonstrates superior amplitude estimation accuracy compared to using OMP alone. This enhanced accuracy can be attributed to two main factors:
  • The sequential Bayesian method, as introduced in this study, thoroughly integrates prior information concerning noise and target amplitude.
  • Building upon OMP-RELAX, the off-grid Bayesian estimation addresses the issue of position deviation arising from the true angle not aligning with the grid point, thereby bolstering estimation precision.
Figure 9 presents the reconstruction amplitude error for non-flickering targets. The error remains consistent whether the data undergo three (OMP- > MMV-OMP- > OMP-RELAX) or four (OMP- > MMV-OMP- > OMP-RELAX- > SBL) processing steps. This consistency is attributed to the minor angle estimation differences between the three and four-step processes for non-flickering targets. Both are in close proximity to the actual target position, resulting in similar amplitude values.

3.2.2. Impact of Heading Angle on DOA Estimation Accuracy

In this section, our evaluation focuses on non-flickering targets. We consider a heading angle range of − 50 ° : 5 : 50 ° , with an SNR set at 20 dB. We conduct 200 Monte Carlo simulations for each specified heading angle, the results of which are depicted in Figure 10.
As can be seen from the figure, under different heading conditions, the proposed method has higher angle measurement accuracy than the OMP algorithm. The target azimuth angle is ( φ 1 , φ 2 , φ 3 ) ∈ [ 20 ° , 30 ° ] for the proposed method, and the OMP algorithm in the heading angle range [ 10 ° , 40 ° ] has a significantly lower angle measurement accuracy compared to other heading angles, with RMSE reaching its peak at heading angles of 20° and 25°. The reason for this phenomenon is that when the heading angle is close to the target angle, the spatial sampling of the aircraft cannot form an effective aperture, so the angle resolution does not improve.

3.2.3. Impact of the Number of Sampling Points on DOA Estimation Results

Sparse reconstruction theory posits that the mutual correlation between dictionary atoms plays a pivotal role in shaping the reconstruction results. Specifically, a diminished mutual correlation is conducive to enhanced reconstruction results. In line with the dictionary design framework proposed in this paper, it is observed that the number of sampling points directly influences the mutual correlation within the dictionary. Hence, we assess the impact of varied sampling point numbers on measurement results by examining the mutual correlation coefficients associated with dictionaries of different sampling sizes.
For the analysis, the aircraft’s maximum spatial sampling points are capped at 1000, maintaining an SNR of 20 dB and a heading angle of −60°. Figure 11 visually represents the mean correlation coefficients of the dictionary matrix across distinct sampling points.
From the figure, a discernible trend emerges. As we increment the number of spatial sampling points, there is a concomitant decline in the dictionary matrix’s average correlation coefficient. This suggests an ongoing enhancement in reconstruction performance. Put another way, with the surge in the count of spatial sampling points, there is an expansion in the spatial synthetic aperture, a contraction in the beam width, and a bolstering of the angular resolution prowess.

3.2.4. Impact of Target Distance on DOA Estimation Accuracy

To assess how the distance between the aircraft and the target affects angle measurement precision, we examine two distinct scenarios: those involving flickering targets and those with non-flickering targets. Figure 12 illustrates the angle measurement precision for targets at varying distances.
From the data presented in the figure, it is evident that even as the distance to the target extends, there is negligible variation in the angle measurement accuracy for both categories of targets—flickering and non-flickering. This observation underscores the robustness of the proposed algorithm with respect to varying target distances.

3.3. Field Experiment

To validate the efficacy of the algorithm, we test the proposed method using real- data. The parameters of the radar system align with those used in the simulation. Additionally, the environment surrounding the radar is characterized by a hilly landscape interspersed with rocks and vegetation, meeting the conditions of a complex terrain.
In our configuration, we set the aircraft’s altitude at 1 km and positioned the radar 5 km from the scene’s center. We arrange three corner reflectors in a triangular pattern on a flat terrain. The aircraft’s actual position is determined using its GPS. Through calculations, we derive the true azimuth angle of the corner reflectors relative to the aircraft. As shown in Figure 13, the star-shaped line represents the true angle values of the three corner reflectors in relation to the aircraft, while the triangular line displays the results of 100 actual measurements. It is clear that during initial sampling, due to a coarse angle estimate, all three reflectors are within the same beam range, exhibiting identical azimuth angles. However, as we increase the measurement count, and consequently the number of spatial sampling points, the spatial synthetic aperture expands, enhancing the accuracy of measurements. By the time we conduct around 40 measurements, we accurately estimate the corner reflectors’ precise azimuth information.

3.4. Computational Complexity Analysis

When the original aperture is used to plan the flight direction, the time complexity of the OMP algorithm is O(1), a total of four multiplication operations of complex numbers, in which the main computational cost is the inverse operation of M × 1 data, where M denotes the number of elements. Subsequently, in the process of precise estimation, according to the target number, the MMV-OMP algorithm needs to cycle for K times, and so does the OMP-RELAX algorithm, with the time complexity of these two algorithms O(K) together. The main amount of computation in this stage is the same as that of the OMP algorithm, which comes from the inverse operation of M × 1 data. For the off-grid algorithm, assuming that the number of iterations is Q and the number of atoms in the dictionary is P, the time complexity is O(QP), a total of (2Q + 3P) multiplication operations of complex numbers along with the inverse operation of the T × T matrix as the main computational cost.
Moreover, the single operation time of the OMP algorithm is 0.054 s, the single operation of the MMV-OMP algorithm is 0.006 s, the single operation of the OMP-RLEAX is 0.015 s, the single operation of the off-SBL algorithm is 1.221 s. In the calculation, only one OMP rough estimation is contained in the entire algorithm; then, the precise estimation is carried out on the basis of the rough estimation. On the whole, the calculation time of the proposed algorithm is 1.4 s per process. What is more, the hardware configuration of the computer is CPU·of Intel(R) Core (TM) 7-4770HOCPU @ 2.2 GHz:2.19 GHz, and the memory is 8 GB. In conclusion, although three different algorithms are designed in this paper, only a small number of cycles are required in the calculation process of each algorithm, and high-dimensional data operation is not involved. Thus, the proposed method can ensure the requirement of fast operation. Furthermore, according to the actual measurement result, the estimation accuracy is significantly improved compared with the single goniometric algorithm.

4. Conclusions

This paper addresses the issue of limited array aperture DOA estimation performance, which falls short in satisfying the requirements for precise detection. Our approach unfolds as follows:
  • Spatial Synthetic Aperture Model Introduction: We present a uniform motion spatial synthetic aperture model. Under this paradigm, the radar enhances DOA estimation accuracy by forming a spatial synthetic aperture through spatial sampling.
  • Algorithm Design for High-Precision DOA Estimation: To fulfill both accuracy and real-time demands, we introduce an algorithm that initially employs OMP (Orthogonal Matching Pursuit) for rough target angle estimation and trajectory planning. This algorithm combines the roughly determined target angle with the aircraft’s spatial samples. A measurement dictionary, built upon these spatial samples, then utilizes both OMP-RELAX and sequential off-grid Bayesian techniques for a precise estimation of the spatial sampling data.
Our method offers three distinct advantages over conventional techniques:
  • Enhanced Precision: The spatially synthesized array boasts a substantially enlarged aperture and heightened resolution. Our sequential estimation process, paired with Bayesian proportion parameter estimation and dictionary refinement, yields more precise DOA data.
  • Increased Efficiency: By employing the greedy method of sparsity reconstruction in DOA estimation, our approach ensures rapid estimations.
  • Angle Glint Issues Resolution: The spatial synthetic array sampling and target angle reconstruction can be viewed as a spatial matched filter. Even if the target’s power is diminished over a specific duration, DOA estimation remains feasible through spatial accumulation across an extended period.
In conclusion, the efficacy of our proposed algorithm is substantiated using computer simulation data. The experimental design comprehensively considers variables like target state (both flickering and non-flickering), heading angle, target distance, and evaluates their influence on algorithm precision. While this study focuses on a uniform linear motion model, we do not explore the potential effects of flight trajectory errors on algorithmic accuracy—an avenue we aim to delve into in future research.

Author Contributions

Conceptualization, Y.F.; methodology, Y.F.; investigation, X.W.; formal analysis, X.W.; writing—original draft preparation, Y.F. and J.M.; writing—review and editing, Y.F. and J.M.; funding acquisition, Y.F. All authors have read and agreed to the published version of the manuscript.

Funding

This work is supported in part by the National Natural Science Foundation of China under Grant No.52107009, Key R&D Project in Shaanxi Province under Grant No.2023-YBGY-411.

Data Availability Statement

The data presented in this study are available on request from the corresponding author after obtaining permission of an authorized person.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Ma, J.J.; Ma, H.; Liu, H.W.; Liu, W.; Cheng, X.Z. A novel DOA estimation for low-elevation target method based on multiscattering center equivalent model. IEEE Geosci. Remote Sens. Lett. 2023, 20, 3501605. [Google Scholar] [CrossRef]
  2. Dai, X.R.; Zhang, X.F.; Wang, Y.F. Extended DOA-matrix method for DOA estimation via two parallel linear arrays. IEEE Commun. Lett. 2019, 23, 1981–1984. [Google Scholar] [CrossRef]
  3. Liu, Y.; Liu, H.W.; Xia, X.G.; Zhang, L.; Jiu, B. Projection techniques for altitude estimation over complex multipath condition-based VHF radar. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2018, 11, 2362–2375. [Google Scholar] [CrossRef]
  4. Shi, Z.; Zhang, X.F.; Xu, L. DOA estimation of multiple sources for a moving array in the presence of phase noise. IET Signal Process. 2019, 13, 29–35. [Google Scholar] [CrossRef]
  5. Liu, G.H.; Chen, H.; Sun, X.Y.; Qiu, R.C. Modified MUSIC algorithm for DOA estimation with Nyström approximation. IEEE Sens. J. 2016, 16, 4673–4674. [Google Scholar] [CrossRef]
  6. Stoica, P.; Nehorai, A. Music, maximum likelihood and Cramer-Rao bound. IEEE Trans. Acoust. Speech Signal Process. 1989, 37, 720–741. [Google Scholar] [CrossRef]
  7. Lee, J.; Park, L.; Chun, J. Weighted two-dimensional root MUSIC for joint angle-Doppler estimation with MIMO radar. IEEE Trans. Aerosp. Electron. Syst. 2019, 55, 1474–1482. [Google Scholar] [CrossRef]
  8. Liu, Z.Y.; Wu, J.Y.; Yang, S.Y.; Lu, W. DOA estimation method based on EMD and MUSIC for mutual interference in FMCW automotive radars. IEEE Geosci. Remote Sens. Lett. 2022, 19, 3504005. [Google Scholar] [CrossRef]
  9. Roy, R.; Kailath, T. ESPRIT-estimation of signal parameters via rotational invariance techniques. Opt. Eng. 1990, 29, 296–313. [Google Scholar]
  10. Li, J.F.; Jiang, D.F.; Zhang, X.F. DOA estimation based on combined unitary ESPRIT for coprime MIMO radar. IEEE Commun. Lett. 2017, 21, 96–99. [Google Scholar] [CrossRef]
  11. Wang, S.H.; Cao, Y.H.; Liu, Y.T. A method of Robust low-angle target height and compound reflection coefficient joint estimation. J. Syst. Eng. Electron. 2022, 33, 322–329. [Google Scholar] [CrossRef]
  12. Chen, S.; Zhao, Y.B.; Hu, Y.L.; Cao, C.H.; Pan, X.J. Target height and multipath attenuation joint estimation with complex scenarios for very high frequency radar. Front. Inf. Technol. Electron. Eng. 2022, 23, 937–949. [Google Scholar] [CrossRef]
  13. Li, J.; Sadler, B.; Viberg, M. Sensor array and multichannel signal processing. IEEE Signal Process. Mag. 2011, 28, 157–158. [Google Scholar] [CrossRef]
  14. Zoltowski, M.; Haber, F. A vector space approach to direction finding in a coherent multipath environment. IEEE Trans. Antennas Propag. 1986, 34, 1069–1079. [Google Scholar] [CrossRef]
  15. Zhang, X.; Liu, X.M.; Yu, H.X. Improved MUSIC algorithm for DOA estimation of coherent signals via Toeplitz and fourth-order-cumulants. Int. J. Control Autom. 2015, 8, 261–272. [Google Scholar] [CrossRef]
  16. Barcelo, M.; Vicario, J.L.; Seco-Granados, G. A reduced complexity approach to IAA beamforming for efficient DOA estimation of coherent sources. Eurasip J. Adv. Signal Process. 2011, 2011, 521265. [Google Scholar] [CrossRef]
  17. Yuan, X. Coherent sources direction finding and polarization estimation with various compositions of spatially spread polarized antenna arrays. Signal Process. 2014, 102, 265–281. [Google Scholar] [CrossRef]
  18. Wang, Q.; Chen, H.; Zhao, G.H.; Chen, B.; Wang, P.C. An improved direction finding algorithm based on Toeplitz approximation. Sensors 2013, 13, 746–757. [Google Scholar] [CrossRef] [PubMed]
  19. Hyder, M.M.; Mahata, K. Direction-of-Arrival estimation using a mixed L2 norm approximation. IEEE Trans. Signal Process. 2010, 58, 4646–4655. [Google Scholar] [CrossRef]
  20. Wen, J.; Liao, B.; Guo, C.T. Spatial smoothing based methods for direction-of-arrival estimation of coherent signals in nonuniform noise. Digit. Signal Process. 2017, 67, 116–122. [Google Scholar] [CrossRef]
  21. Qi, C.Y.; Wang, Y.L.; Zhang, Y.S.; Han, Y. Spatial difference smoothing for DOA estimation of coherent signals. IEEE Signal Process. Lett. 2005, 12, 800–802. [Google Scholar]
  22. Jiang, G.J.; Yang, Y.L. Synthetic sparse nested array with extended aperture and reduced mutual coupling for DOA estimation. Signal Process. 2023, 204, 108849. [Google Scholar] [CrossRef]
  23. Gorodnitsky, I.F.; Rao, B.D. Sparse signal reconstruction from limited data using FOCUSS: A re-weighted minimum norm algorithm. IEEE Trans. Signal Process. 1997, 45, 600–616. [Google Scholar] [CrossRef]
  24. Wang, J.; Kwon, S.; Shim, B. Generalized orthogonal matching pursuit. IEEE Trans. Signal Process. 2012, 60, 6202–6216. [Google Scholar] [CrossRef]
  25. Gerstoft, P.; Mecklenbräuker, C.; Xenaki, A.; Nannuru, S. Multisnapshot sparse Bayesian learning for DOA. IEEE Signal Process. Lett. 2016, 23, 1469–1473. [Google Scholar] [CrossRef]
  26. Chen, J.; Huo, X.M. Theoretical results on sparse representations of multiple-measurement vectors. IEEE Trans. Signal Process. 2006, 54, 4634–4643. [Google Scholar] [CrossRef]
  27. Zhuang, S.Y.; Zhao, W.; Wang, Q.; Wang, Z.; Chen, L.; Huang, S.L. A high-resolution algorithm for supraharmonic analysis based on multiple measurement vectors and Bayesian compressive sensing. Energies 2019, 12, 2559. [Google Scholar] [CrossRef]
  28. Herzet, C.; Drémeau, A.; Soussen, C. Relaxed recovery conditions for OMP/OLS by exploiting both coherence and decay. IEEE Trans. Inf. Theory 2016, 62, 459–470. [Google Scholar] [CrossRef]
  29. Yang, Z.; Xie, L.H.; Zhang, C.S. Off-grid direction of arrival estimation using sparse Bayesian inference. IEEE Trans. Signal Process. 2012, 61, 38–43. [Google Scholar] [CrossRef]
  30. Chen, F.H.; Dai, J.S.; Hu, N.; Ye, Z.F. Sparse Bayesian learning for off-grid DOA estimation with nested arrays. Digit. Signal Process. 2018, 82, 187–193. [Google Scholar] [CrossRef]
  31. Lin, J.C.; Ma, X.C.; Yan, S.F.; Jiang, L. Pattern synthesis of sparse linear array by off-grid Bayesian compressive sampling. Electron. Lett. 2015, 51, 2141–2143. [Google Scholar] [CrossRef]
Sensors 23 08690 g001
Figure 1. Radar target space geometry.
Figure 1. Radar target space geometry.
Sensors 23 08690 g001
Sensors 23 08690 g002
Figure 2. Spatial synthetic aperture geometry.
Figure 2. Spatial synthetic aperture geometry.
Sensors 23 08690 g002
Sensors 23 08690 g003
Figure 3. Diagram representing the flow of the proposed algorithm.
Figure 3. Diagram representing the flow of the proposed algorithm.
Sensors 23 08690 g003
Sensors 23 08690 g004
Figure 4. Initial rough estimation result. (The area pointed by the red arrow in the figure is an enlargement of the target area.)
Figure 4. Initial rough estimation result. (The area pointed by the red arrow in the figure is an enlargement of the target area.)
Sensors 23 08690 g004
Sensors 23 08690 g005
Figure 5. OMP estimation result of refined dictionary.
Figure 5. OMP estimation result of refined dictionary.
Sensors 23 08690 g005
Sensors 23 08690 g006
Figure 6. OMP-RELAX estimation result.
Figure 6. OMP-RELAX estimation result.
Sensors 23 08690 g006
Sensors 23 08690 g007
Figure 7. SBL estimation result.
Figure 7. SBL estimation result.
Sensors 23 08690 g007
Sensors 23 08690 g008
Figure 8. Comparison of angle estimation for flickering and non-flickering targets.
Figure 8. Comparison of angle estimation for flickering and non-flickering targets.
Sensors 23 08690 g008
Sensors 23 08690 g009
Figure 9. Reconstruction error of a non-flickering target.
Figure 9. Reconstruction error of a non-flickering target.
Sensors 23 08690 g009
Sensors 23 08690 g010
Figure 10. RMSE results across different heading angles.
Figure 10. RMSE results across different heading angles.
Sensors 23 08690 g010
Sensors 23 08690 g011
Figure 11. Mean correlation coefficients of the dictionary across varied spatial sampling points.
Figure 11. Mean correlation coefficients of the dictionary across varied spatial sampling points.
Sensors 23 08690 g011
Sensors 23 08690 g012
Figure 12. Angle measurement precision for targets across different distances.
Figure 12. Angle measurement precision for targets across different distances.
Sensors 23 08690 g012
Sensors 23 08690 g013
Figure 13. Angle estimation results of the real data by the proposed method.
Figure 13. Angle estimation results of the real data by the proposed method.
Sensors 23 08690 g013
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Fang, Y.; Wei, X.; Ma, J. High-Precision DOA Estimation Based on Synthetic Aperture and Sparse Reconstruction. Sensors 2023, 23, 8690. https://doi.org/10.3390/s23218690

AMA Style

Fang Y, Wei X, Ma J. High-Precision DOA Estimation Based on Synthetic Aperture and Sparse Reconstruction. Sensors. 2023; 23(21):8690. https://doi.org/10.3390/s23218690

Chicago/Turabian Style

Fang, Yang, Xiaolong Wei, and Jianjun Ma. 2023. "High-Precision DOA Estimation Based on Synthetic Aperture and Sparse Reconstruction" Sensors 23, no. 21: 8690. https://doi.org/10.3390/s23218690

APA Style

Fang, Y., Wei, X., & Ma, J. (2023). High-Precision DOA Estimation Based on Synthetic Aperture and Sparse Reconstruction. Sensors, 23(21), 8690. https://doi.org/10.3390/s23218690

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

Article Metrics

Back to TopTop