Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
Next Article in Journal
A Dual-Mode Surface Acoustic Wave Delay Line for the Detection of Ice on 64°-Rotated Y-Cut Lithium Niobate
Previous Article in Journal
Fusion of Multimodal Imaging and 3D Digitization Using Photogrammetry
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

An Improved OMP Algorithm for Enhancing the Anti-Interference Performance of Array Antennas

1
College of Big Data and Intelligent Engineering, Southwest Forestry University, Kunming 650224, China
2
College of Mathematics and Physics, Southwest Forestry University, Kunming 650224, China
*
Author to whom correspondence should be addressed.
Sensors 2024, 24(7), 2291; https://doi.org/10.3390/s24072291
Submission received: 20 February 2024 / Revised: 30 March 2024 / Accepted: 31 March 2024 / Published: 4 April 2024
(This article belongs to the Section Radar Sensors)

Abstract

:
The demand for precise positioning in noisy environments has propelled the development of research on array antenna radar systems. Although the orthogonal matching pursuit (OMP) algorithm demonstrates superior performance in signal reconstruction, its application efficacy in noisy settings faces challenges. Consequently, this paper introduces an innovative OMP algorithm, DTM_OMP_ICA (a dual-threshold mask OMP algorithm based on independent component analysis), which optimizes the OMP signal reconstruction framework by utilizing two different observation bases in conjunction with independent component analysis (ICA). By implementing a mean mask strategy, it effectively denoises signals received by array antennas in noisy environments. Simulation results reveal that compared to traditional OMP algorithms, the DTM_OMP_ICA algorithm shows significant advantages in noise suppression capability and algorithm stability. Under optimal conditions, this algorithm achieves a noise suppression rate of up to 96.8%, with its stability also reaching as high as 99%. Furthermore, DTM_OMP_ICA surpasses traditional denoising algorithms in practical denoising applications, proving its effectiveness in reconstructing array antenna signals in noisy settings. This presents an efficient method for accurately reconstructing array antenna signals against a noisy backdrop.

1. Introduction

The frequency-modulated continuous-wave (FMCW) radar positioning system based on an antenna array achieves target localization by transmitting and receiving echo signals with varying frequencies due to object reflections. The determination of frequency components in the signals is a critical aspect of this technology [1,2,3]. However, in practical scenarios, the received echo signals by the antenna array are often accompanied by noise interference, making it challenging to accurately analyze the frequency and phase from the mixer [4,5]. This ultimately affects the localization accuracy of the FMCW radar positioning system [6,7,8].
Therefore, effectively reducing the interference of noise in the echo signals to achieve accurate frequency and phase have become hot research topics. Traditional methods are intermediate frequency signal filtering [9] or windowing [10] to address this issue. The discrete-time signal processing theory proposed by Oppenheim, Schafer, and Buck laid the foundation for signal sampling and reconstruction [11]. Parks and Burrus introduced methods for designing digital filters [12] to cope with the complexity of signal processing, but the complexity of filter design also increased accordingly. Hayes’s statistical digital signal processing approach emphasizes signal analysis in the presence of random interference [13] and uses probability models to describe the relationship between signals and noise. While this method can effectively reduce the interference of noise to some extent, its accuracy is not ideal. In order to achieve effective denoising of the echo signals received by the antenna array in a complex noise interference environment, more advanced techniques are needed compared to traditional methods.
Compressed sensing (CS) [14,15,16,17,18] exploits the sparsity of signals or their sparsity in a certain transform domain, breaking the limitations of the traditional Nyquist sampling theorem. By compressing signal sampling, CS allows for perfect reconstruction of signals even at sampling frequencies far below twice the highest Nyquist frequency. In the context of FMCW radar positioning systems based on antenna arrays, the sparsity of the target object relative to the positioning background [19] made compressed sensing suitable. The frequency count of the intermediate signals, formed by the echoes received by the antenna array after reflecting off the target object, is only related to the number of target objects. This aligns with the sparsity of the frequency of intermediate signals in FMCW radar systems based on antenna arrays, which is equivalent to the sparsity of the target relative to its background. This similarity makes the application of compressed sensing technology essential for signal processing in antenna array FMCW radar systems [20,21]. The literature [22] introduces a radar signal processing method based on compressed sensing, providing new approaches and technical support for reconstructing radar signals in complex environments. However, it is found that this method yields unstable results in reconstructing noisy radar signals.
The OMP algorithm [23], as a classic method in compressed sensing, has been widely studied due to its advantages, such as low algorithm complexity, fast computation speed, and high reconstruction accuracy for sparse signals. However, it is sensitive to noise, particularly in scenarios where the reconstructed signal contains substantial noise interference, leading to unstable reconstruction results. Additionally, the algorithm’s performance in reconstructing sparse signals depends on prior knowledge of signal sparsity. When applying the OMP algorithm to reconstruct the intermediate signals formed by the noisy signals received by antenna arrays, a key challenge arises in enhancing the algorithm’s resistance to noise interference and determining the sparsity of intermediate signals to optimize the construction of the sensing matrix [24]. Thus, addressing the challenging task of setting prior knowledge to determine the sparsity of intermediate signals and eliminating noise interference to achieve accurate reconstruction of noisy signals using the OMP algorithm becomes crucial.
This paper proposes an improved algorithm, DTM_OMP_ICA, for the orthogonal matching pursuit (OMP) algorithm under a masking strategy [25], based on independent component analysis (ICA) [26,27]. This algorithm enhances the traditional OMP framework. It does so by employing OMP algorithms based on two distinct observation bases to derive two different sparse representations of the input signal. These representations are then used as dual inputs to the ICA algorithm, aimed at bolstering the robustness of the signal’s sparse representation. Furthermore, the algorithm integrates a mean mask strategy to accomplish denoising reconstruction of the input signal. The goal of this novel signal processing approach is to ensure the robustness and stability of the algorithm while accurately reconstructing the intermediate frequency signals constituted by noisy echo signals received by array antennas. The organization of the paper is as follows: Section 1 introduces the research motivation, significance, and content. Section 2 provides an overview of the process strategy of the orthogonal matching pursuit algorithm. In Section 3, the improved algorithm based on the masking strategy for the orthogonal matching pursuit algorithm (DTM_OMP_ICA) is proposed, and its theoretical feasibility is validated through mathematical formula derivation. Section 4 analyzes the frequency characteristics of signals reconstructed by the DTM_OMP_ICA algorithm under different noise environments, verifying the accuracy and stability of the algorithm. Additionally, the DTM_OMP_ICA algorithm is applied to denoise and reconstruct simulated signals received by continuous modulation wave radar array antennas in noisy environments to explore its performance in practical applications. Section 5 summarizes the proposed algorithm and provides prospects for future research.

2. The Principle of CS

CS exploits the sparsity of a signal in a certain transformed domain to achieve signal reconstruction with only a small number of compressed projection data. The core idea of its reconstruction algorithm is to project the signal into a lower-dimensional space and then solve an optimization problem to reconstruct the original signal. The objective of this optimization problem is to find the sparsest representation of the signal, minimizing the difference between the projected data and the original signal. In this way, compressed sensing can significantly reduce sampling and storage requirements while maintaining a high level of reconstruction quality.

2.1. The Sampling Process of Compressed Sensing

According to compressed sensing theory, signal compression encoding does not directly measure the signal itself. Instead, it involves projecting the signal onto a set of measurement matrices to obtain observed values. A one-dimensional signal, denoted as x (with a length of N ), can be represented using an orthogonal transform basis (sparse basis) ψ and sparse representation θ , where the sparsity level is K , as shown in (1):
x = ( ψ × θ ) T
where ψ N × N , θ N × 1 , and “ ( · ) T ” represent the matrix transpose operation.
Assuming an observation basis φ , which is uncorrelated with the transform basis ψ , allows for the compression of the signal x . This compression is achieved by applying φ , resulting in a one-dimensional observation signal y of length M . Consequently, this process compresses signal x from its original length N to M , as outlined in Equation (2).
y = φ × x
where φ M × N , x N × 1 , y M × 1 .
To ensure precise signal reconstruction, the design of the observation basis φ requires that during the transformation of the signal from x to y , the M observed values obtained should not destroy the essential information of the original signal x . This means that the number of observed values M should satisfy the following equation, and it should hold that K M N :
M K × log N K
Combining (1) and (2), we can obtain (4):
y = φ × x = φ × ψ × θ = A × θ
where A represents the sensing matrix. A = φ × ψ , A ϵ M × N .
Due to the fact that the number of equations represented by (4) is much smaller than the number of unknowns, the solutions are non-unique, affecting the precise reconstruction of the signal x . However, because the signal x is k -sparse in its orthogonal transform basis ψ , if the φ in (4) satisfies the restricted isometry property (RIP) [28,29], then the signal x can be accurately reconstructed by the optimal solution from M observed values. The equivalent condition for RIP property is that the observation basis φ and the sparse basis ψ are uncorrelated.

2.2. Compressed Sensing Signal Reconstruction Algorithm

When the observation basis φ satisfies the RIP criterion, the compressed sensing signal reconstruction algorithm can obtain the sparse representation θ of the signal through solving the inverse problem of (4). Subsequently, it can precisely reconstruct the length- N signal x with K sparsity in the orthogonal transform basis ψ from the length- M observation vector y . Its mathematical expression is:
m i n θ θ 0     s . t .     y = A × θ
Currently, to effectively solve for θ to achieve signal reconstruction, algorithms are mainly divided into three major categories: convex optimization algorithms [30,31,32], greedy matching pursuit methods, and combination algorithms [33].
Among them, the greedy matching pursuit method iteratively searches for the support set of the sparse vector and reconstructs the signal using a constrained support least-squares estimate. Due to its low complexity, this algorithm can rapidly reconstruct signals and is suitable for real-time computing scenarios. Such algorithms mainly include matching pursuit (MP) [34], orthogonal matching pursuit (OMP), the compressive sampling matching pursuit (CoSaMP) algorithm [35], and the regularized OMP (ROMP) [36,37] algorithm, among others. However, these algorithms are sensitive to noise, and their noise resistance performance can be improved by combining them with other methods, such as preprocessing with wavelet transforms [38] or improving the sparse representation of signals through various dictionary learning methods [39].

2.3. Orthogonal Matching Pursuit (OMP)

The core challenge of OMP lies in determining which column a i of the matrix A contribute significantly to the vector y in (6). The distance between two vectors can be measured using the inner product of vectors. Therefore, vector y can be projected onto the column vectors of matrix A , and the numerical value of the inner product is used to assess the contribution of the column vectors of matrix A to vector y .
C o n t r i b u t i o n y , a i = y , a i a i
where a i represents the i -th column of matrix A . θ in (4), as the projection of signal x onto its orthogonal transformation basis ψ , is K-sparse; therefore, the algorithm only needs to iterate K times to find the K column vectors in matrix A that contribute the most to y , along with their positions. This allows for solving the current optimal solution θ to complete the reconstruction of signal x . The above describes the workflow of the orthogonal matching pursuit algorithm, and the pseudo-code for the algorithm is shown as follows:
In the pseudo-code shown in Algorithm 1, the input A represents the sensing matrix, y is a one-dimensional observation vector of length M , and K denotes the sparsity level of the original signal x (typically, the choice of K is related to the energy distribution of the one-dimensional signal x in its orthogonal transform basis ψ ). The output θ represents the sparse representation of the signal x .
Algorithm 1: Orthogonal matching pursuit (OMP)
Input: A , y , K
Output: θ
1 Initialization r 0 = θ , Λ 0 = ;
2 Normalize all columns of A to unit L2 norm;
3 for k = 1,2, …, K do
4            Step 1:       λ k = j Λ k 1 arg max a j , r κ 1 ;
5            Step 2:       Λ k = Λ K 1 λ k ;
6            Step 3:       x k Λ k = x arg min A Λ k x y 2 ;
7            Step 4:       θ k ~ = A x k ;
8            Step 5:       γ k = θ θ k ~ ;
9 end
The above algorithm illustrates that the composition of the sensing matrix A is crucial for the reconstruction performance of the orthogonal matching pursuit algorithm [23]. It has been proven [40] that independently and identically distributed Gaussian random measurement matrices can serve as universal compressed sensing observation basis. Therefore, random Gaussian matrices are commonly used as observation basis. In addition to random Gaussian matrices, other commonly used observation basis include random Bernoulli matrices [41], random orthogonal matrices [42], Toeplitz matrices [43], and sparse random matrices [44], among others.
The theory of compressed sensing (CS) and its applications in denoising and signal reconstruction continue to garner interest among researchers. Particularly, due to its efficacy in processing sparse signals, the orthogonal matching pursuit (OMP) algorithm has emerged as a focal point of study. In environments characterized by the presence of noise, the robust recovery algorithm proposed by V Meena and G Abhilash, which improves upon the OMP algorithm, has demonstrated the effectiveness of compressed sensing under conditions of high signal-to-noise ratios [45]. Similarly, TJ Thomas and colleagues developed a new algorithm that not only enhances the recovery capabilities of OMP but also offers a new direction for the denoising of ECG signals [46].
Further, the comprehensive review by L Li et al. evaluates the potential applications of the OMP algorithm and its variants in noisy settings, providing valuable perspectives on the feasibility of compressed sensing across various application scenarios [47]. Moreover, the research by C Cheng and D Lin, applying the OMP algorithm in the process of image restoration, showcases the potential and challenges of compressed sensing theory in the field of image processing [48].
These developments indicate that while compressed sensing and the OMP algorithm have proven effective in multiple application scenarios, signal reconstruction in high-noise environments remains a challenge. To address these challenges effectively, the DTM_OMP_ICA algorithm introduced in this paper differs from traditional OMP algorithms in that it does not rely on a priori assumptions about signal sparsity. The algorithm can adaptively determine the sparsity level of the input signal, enhancing its flexibility and efficiency in processing signals with unknown or variable sparsity. Furthermore, by employing a threshold masking technique, the DTM_OMP_ICA algorithm significantly improves the accuracy and stability of signal reconstruction in noisy conditions.

3. Dual-Threshold Mask OMP Based on ICA

3.1. Principle of ICA

Independent component analysis (ICA) [24,25] is a multivariate statistical analysis method that aims to separate mutually independent components from mixed multidimensional data based on the independence assumption. Its mathematical expression is shown in (7).
S = W × X
where S represents the data matrix of independent components, X is the original observation matrix, and W is the mixing matrix used to transform the original data X into the independent components S .
In the ICA algorithm model, the core objective is to estimate the mixing matrix W . This objective is achieved by maximizing the independence between components. Typically, the independence between components is measured through information entropy or higher-order statistics, serving as a guide for extracting independent components. In practical scenarios, noisy signals are often a mixture of source signals and additive noise. In other words, the original observation matrix X in (7) can be represented as follows:
X = α 1 x 1 + α 2 x 2
where X represents the noisy signal, x 1 is the pure signal, x 2 is the additive noise, α 1 and α 2 represent the energy proportions of the pure signal and additive noise in the noisy signal X , and their values are determined by the signal-to-noise ratio of the noisy signal X .
Assuming that the clean signal x 1 and additive noise x 2 come from mutually independent, different sources, which satisfy the independence assumption of the ICA algorithm, this allows for the ICA algorithm to effectively capture the independence information between the original signal and additive noise in the noisy signal by estimating the mixing matrix W in (8). This provides strong support for signal denoising processing.

3.2. DTM_OMP_ICA

When using the OMP algorithm to reconstruct the noisy signal X in (8), its observed that signal y can be expressed as follows:
y = φ × X = φ × α 1 x 1 + α 2 x 2 = A × θ
where φ , α 1 , α 2 , A , θ , x 1 , and x 2 have the same meanings as defined in (2), (4), and (8).
In the theory of compressed sensing, random matrices are commonly chosen as the observation basis φ , and the Fourier transform matrix is used as the orthogonal transform basis ψ . Moreover, for the observation basis φ , the number of rows M and the number of columns N satisfy (3). This implies that the row vectors of the observation basis φ are usually linearly independent, making the observation basis φ typically full-rank in rows. Additionally, the Fourier transform basis, being the Fourier transform of the identity matrix E , is also a full-rank matrix. Therefore,
r φ = M ,             r ψ = N
At the same time, the sensing matrix A = φ × ψ . Therefore, it is inevitable that the sensing matrix A is a full-rank matrix in rows, i.e.,
r A = r φ = M
For sensing matrix A , there must exist a pseudo-inverse matrix A + such that (12) and (13) hold as follows:
A × A + × A = A
A + × A × A + = A +
Therefore, (9) can be rewritten as
θ = A + × y = A + × φ × α 1 x 1 + α 2 x 2 = C × α 1 x 1 + α 2 x 2
where C = A + × φ .
Since the term C in (14) involves only matrix correlations, the sparse representations ( θ R G M , θ S R M ) obtained by OMP algorithms based on two different random matrices as observation bases ( φ R G M , φ S R M ) for the noisy signal X can be considered as different linear transformations of the noisy signal X , as shown in (15) and (16):
θ R G M = A R G M + × φ R G M × α 1 x 1 + α 2 x 2 = C R G M × α 1 x 1 + α 2 x 2
θ S R M = A S R M + × φ S R M × α 1 x 1 + α 2 x 2 = C S R M × α 1 x 1 + α 2 x 2
These distinct sparse representations of the noisy signal X correspond to different weighted mixtures of the clean signal x 1 and the additive noise x 2 . Simultaneously, x 1 and x 2 originate from independent sources, satisfying the independence assumption of the ICA algorithm. Therefore, θ R G M and θ S R M can be used as inputs to the ICA algorithm. Leveraging the ICA algorithm allows for obtaining sparse representations, namely, θ 1 related solely to the clean signal x 1 and θ 2 related exclusively to the additive noise x 2 , from the noisy signal X . Subsequently, θ 1 will be utilized for denoising and signal reconstruction of the noisy signal X .
Based on the above theory, the present paper proposes the DTM_OMP_ICA algorithm to eliminate or reduce the impact of noise on signal reconstruction. The flowchart of the DTM_OMP_ICA algorithm is shown in Figure 1.
In Figure 1, when denoising and reconstructing the noisy signal X_Noise, the proposed DTM_OMP_ICA algorithm employs OMP algorithms based on two different observation bases (RGM-OMP and SRM-OMP) to construct a dual-path process. This results in obtaining distinct sparse representations, namely, θ R G M and θ S R M , for the noisy signal X_Noise. The algorithm combines the usage of the ICA algorithm (Fast ICA) to decouple and obtain sparse representations most correlated with the X_Clear component in the noisy signal X_Noise, achieving denoising and reconstruction of X_Noise. To ensure that the reconstructed signal retains the frequencies and corresponding amplitudes of X_Clear from X_Noise as much as possible, a mean mask technique is applied based on the frequency domain mean of the reconstructed signal. This technique generates a mask matrix suitable for the frequency domain of X_Noise, multiplies it with the spectrum of X_Noise, and employs inverse Fourier transform for denoising and reconstructing the noisy signal X_Noise.
The pseudo-code for DTM_OMP_ICA is shown as follows:
In the pseudo-code shown in Algorithm 2, “path 1” and “path 2” represent the dual-path process constructed using OMP algorithms with different observation bases. “ O r t h ( r a n d ( · ) ) ” is utilized to generate a random Gaussian matrix. “ S p a r s e R a n d o m M t x ( · ) ” is employed to generate a sparse random matrix. “ O M P ( · ) ” indicates the use of the OMP algorithm.
Algorithm 2: Dual-threshold mask OMP based on ICA (DTM_OMP_ICA)
Input:  X _ N o i s e
Output:  R e c o n s t r u c t e d   X _ C l e a r
1 Initialization K = 0 ,   M = 0 ,   N = l e n ( X _ N o i s e ) ;
2 do
3      Step 1:   K = O u t L i n e _ E x t r a ( X _ N o i s e ) ;
4      Step 2:   M = K × l o g ( N / K ) ;
5      Step 3:   ψ = 1 N × F F T ( I ) ;
6      path1 do
7          Step 1:   φ R G M = O r t h ( r a n d ( M , N ) ) ;
8          Step 2:   y R G M = φ R G M × X _ N o i s e ; A R G M = φ R G M × ψ ;
9          Step 3:   θ R G M = O M P ( A R G M , y R G M , K ) ;
10     end path1
11     path2 do
12          Step 1:   φ S R M = S p a r s e R a n d o m M t x ( M , N , d = 4 ) ;
13          Step 2:   y S R M = φ S R M × X _ N o i s e ; A S R M = φ S R M × ψ ;
14          Step 3:   θ S R M = O M P ( A S R M , y S R M , K ) ;
15    end path2
16      Step 4:   θ 1 = c o r r ( F a s t I C A ( θ R M G , θ S R M ) ) ;
17      Step 5:   R e c o n s t r u c t e d _ X = ψ T × θ 1 T ;
18      Step 6:   M a s k = M e a n M a s k ( F F T ( R e c o n s t r u c t e d _ X ) ) ;
19      Step 7:   F X = d o t ( F F T ( X _ N o i s e ) ,   M a s k ) ;
20      Step 8:   R e c o n s t r u c t e d   X _ C l e a r = i F F T ( F X ) ;
21 end
The pseudo-code function O u t L i n e _ E x t r a represents the outlier detection algorithm [49], used to identify anomalies in the signal spectrum and determine the signal’s sparsity level. Its definition is as follows:
O u t L i n e _ E x t r a 1 N i = 1 N X i m e d i a n X
where N represents the size of the data to be examined X , X i denotes the i -th sample point in the data X .
The function M e a n M a s k represents the generation of a mask matrix using the frequency domain mean of the input signal. It is defined as follows:
M e a n M a s k i 0 ,             F x i m e a n ( F x ) < 0 1 ,             F x i m e a n ( F x ) 0
where “ M a s k ” represents the generated mask matrix; “ F x ” is used to denote the spectrum of the noisy signal x after Fourier transformation, and “ F x ” and “ M a s k ” have the same length; “ m e a n ( · ) ” indicates the mean calculation; and “ i ” represents the i -th spectral line.

4. Experiments and Results

This paper conducted two sets of simulation experiments. The first set evaluated the reconstruction performance of a noisy multi-frequency signal, focusing on accuracy, stability, and processing time. The second set involved the reconstruction experiment of a noisy array antenna. The comparative algorithms for both sets were the RGM-OMP and SRM-OMP algorithms. Each experiment category consisted of 100 rounds of independent repetitions.

4.1. Multi-Frequency Signals with Noise Experiment

4.1.1. Data

The experiment simulated a multi-frequency signal with four different frequencies and added noise with varying SNR. This was carried out to validate the stability and accuracy of the proposed DTM_OMP_ICA algorithm in denoising and reconstructing multi-frequency signals in the presence of noise.
The multi-frequency signal used in this paper was represented as:
x = cos 2 π f 1 t + 0.7 cos 2 π 2 f 2 t + 1.6 cos 2 π f 3 t + 1.2 cos 2 π f 4 t
where f 1 = 500   H z ,   f 2 = 1000   H z ,   f 3 = 1500   H z ,   f 4 = 2000   H z .
The main parameter settings of the algorithm were shown in Table 1.
The noisy signal used in this paper is represented as
x = x + n o i s e
where n o i s e represents Gaussian white noise. Four sets of noisy signals with different SNRs (0 dB, 3 dB, 5 dB, and 9 dB) were used to compare in the experiments. The time-domain waveforms of all multi-frequency signals used in this paper are shown in Figure 2.
In Figure 2, the time-domain waveform graphs from left to right represent the clean multi-frequency signal x and the noisy multi-frequency signal x at SNRs of 0 dB, 3 dB, 5 dB, and 9 dB, respectively.

4.1.2. Results of Multi-Frequency Signal with Noise

Spectral Analysis of Reconstructed Signals

After 100 rounds of independent repeated experiments, the spectra of the reconstructed signals S using three different algorithms for the signals at the mentioned five different SNRs were shown in Figure 3.
Figure 3 presents a comparative analysis of the signal reconstruction capabilities of RGM-OMP, SRM-OMP, and the DTM_OMP_ICA algorithm proposed in this paper under various noise conditions. The first row of five spectral graphs (Figure 3(a(1)–a(5))) displays the spectrum of the multi-frequency signal x and its noisy counterparts x under different signal-to-noise ratios (SNRs of 0 dB, 3 dB, 5 dB, and 9 dB, respectively). Following these, the next three rows correspond to the spectra of the reconstructed signal S after processing by the RGM-OMP, SRM-OMP, and DTM_OMP_ICA algorithms, respectively, for the aforementioned multi-frequency signals.
Figure 3 clearly demonstrates that in an ideal noise-free environment, all three algorithms exhibit excellent reconstruction performance, successfully retaining the original frequency components of the signal in the spectrum of the reconstructed signal S (as shown in Figure 3(b(1),c(1),d(1))). However, as the noise level increases, the reconstructed signal spectra of the RGM-OMP and SRM-OMP algorithms begin to show evident frequency components arising from noise interference (as depicted in Figure 3(b(2)–b(5),c(2)–c(5))), and the number of interference frequency components appearing in the spectra of the signals reconstructed by the RGM-OMP and SRM-OMP algorithms increases with the noise, largely due to the randomness of the observation basis they employ. In contrast, the signal reconstructed by the DTM_OMP_ICA algorithm not only retains the integrity of the original pure signal x frequency components in its spectrum from the noisy multi-frequency signal x but also significantly reduces the frequency components arising from noise interference (as shown in Figure 3(d(2)–d(5))). The experimental results indicate that the DTM_OMP_ICA algorithm surpasses traditional OMP algorithms in robustness when reconstructing noisy multi-frequency signals.

Evaluations of Reconstructed Signals

To objectively assess the noise-resistant reconstruction stability of the algorithm, this paper proposed the following stability evaluation.
(i)
Evaluation of Signal Reconstruction Accuracy
This paper employs statistical parameters ((21)–(23)) to perform algorithm stability analysis.
A.
The accuracy of reconstructed signal frequencies, denoted as A c c F :
A c c F = A c c F N u m T o t a l N u m × 100 %
where A c c F N u m represents the number of times the correct frequencies appear in the reconstructed signal and T o t a l N u m is the total number of frequency components in the reconstructed signal.
B.
The error rate of reconstructed signal frequencies, denoted as E r r F :
E r r F = E r r F N u m T o t a l N u m × 100 %
where E r r F N u m represents the number of times incorrect frequencies appear in the reconstructed signal.
C.
The accuracy of both the frequency and amplitude of the reconstructed signal, denoted as A c c F A :
A c c F A = A c c F A N u m A c c F N u m × 100 %
where A c c F A N u m represents the number of times both the correct frequency and amplitude appear in the reconstructed signal.
After 100 rounds of independent repeated experiments, the line charts illustrating the frequency statistical parameters of the reconstructed signals using three different algorithms are shown in Figure 4, Figure 5 and Figure 6, respectively.
Figure 4, Figure 5 and Figure 6 depicted line charts of frequency statistical parameters, including A c c F , E r r F , and A c c F A , for the signals reconstructed by three different algorithms in 100 rounds of independent repeated experiments. Subplots (a–d) in Figure 4, Figure 5 and Figure 6 corresponded to the frequency statistical parameters of the reconstructed signals for different SNRs of the noisy multi-frequency signals x (SNR of 0 dB, 3 dB, 5 dB, and 9 dB, respectively) using each algorithm. In each subplot, the orange dashed line represents the proposed DTM_OMP_ICA algorithm, while the purple and blue solid lines correspond to the RGM-OMP and SRM-OMP algorithms, respectively. To simplify and visualize the frequency statistical parameters of the reconstructed signals, the values of the parameters shown in (18), (19), and (20) were averaged every 10 rounds of independent repeated experiments.
Figure 4, Figure 5 and Figure 6 indicated that the reconstruction performance of the traditional OMP algorithm was affected by changes in the SNR values of the noisy multi-frequency signals. In comparison, the algorithm proposed in this paper exhibited a relatively stable denoising and reconstruction capability for noisy multi-frequency signals, and it had a higher probability of reconstructing the frequency components of the pure signals within the noisy multi-frequency signals.
After 100 independent repeated experiments, this paper conducted an overall analysis of the frequency statistical parameters of the reconstructed signals using three different algorithms at different SNRs, as shown in Figure 7.
In Figure 7, the x-axis represents three different frequency statistical parameters, and the y-axis represents the mean values of frequency statistical parameters of signals reconstructed by different algorithms after 100 rounds of independent repeated experiments with four different signal-to-noise ratios (SNRs) for noisy multi-frequency signals. The purple, blue, and orange bars represent the RGM-OMP, SRM-OMP, and DTM_OMP_ICA algorithms, respectively.
The results shown in Figure 7 indicate that the proposed DTM_OMP_ICA algorithm excelled in the accuracy of reconstructing signal frequencies ( A c c F ), achieving an average accuracy of 96.8%, significantly surpassing RGM-OMP with 45.01% and SRM-OMP with 44.38%. In terms of the error rate of reconstructed signal frequencies ( E r r F ) , DTM_OMP_ICA was at 3.18%, much lower than RGM-OMP at 54.99% and SRM-OMP at 55.62%. Additionally, DTM_OMP_ICA attained an accuracy of 79.88% in reconstructing the frequencies and amplitudes of signals ( A c c F A ), surpassing other algorithms by 15.94% and 15%, respectively.
(ii)
Evaluation of Signal Reconstruction Stability
The mean and variance of the reconstructed frequencies implemented by the algorithm were statistically analyzed using 100 experiments. First, k -means clustering was used to cluster the reconstructed frequencies obtained by the algorithm. Since the signal X was a linear combination of four frequencies, the number of clusters was set to 4. The statistical results of the three algorithms mentioned above are shown in Figure 8.
Figure 8 meticulously illustrates the frequency clustering results of the reconstructed signal S after processing the noisy multi-frequency signal x through RGM-OMP, SRM-OMP, and the DTM_OMP_ICA algorithm proposed in this paper under different signal-to-noise ratio (SNR) environments (0 dB, 3 dB, 5 dB, and 9 dB). Each row of four subfigures corresponds to the same algorithm, with the three rows of subfigures from top to bottom, respectively, showcasing the frequency clustering results of the reconstructed signal S after using the RGM-OMP, SRM-OMP, and DTM_OMP_ICA algorithms (for example: Figure 8(c(1)) displays the frequency clustering result of the reconstructed signal S after processing through the DTM_OMP_ICA algorithm in a noise environment with an SNR of 0 dB).
A comparative analysis of the results shown in each subfigure of Figure 8 clearly reveals that after reconstructing the noisy multi-frequency signal x using traditional OMP signal reconstruction algorithms based on two different observation bases, the frequency clustering results of the reconstructed signal S demonstrate that the average value of each frequency cluster deviates from the frequency of the original noise-free signal. Additionally, the variance within each frequency cluster in the clustering results is considerable, especially more pronounced in environments with stronger noise interference. For example, in a noise environment with an SNR of 0 dB, the frequency distribution variance within the four frequency clusters of the signal reconstructed by the RGM-OMP algorithm are 32,355, 24,389, 10,043, and 7251, as shown in Figure 8(a(1)), where each frequency cluster is surrounded by numerous outlier frequencies with a wide distribution range. Similarly, in the frequency clustering results of the signal reconstructed by the SRM-OMP algorithm shown in Figure 8(b(1)), each frequency cluster is also surrounded by outlier frequencies with a large distribution range, with variances within each cluster of 10,587, 12,967, 19,081, and 4697. These experimental results indicate that the stability of the frequency distribution of the reconstructed signal weakens as noise interference increases after reconstructing the noisy multi-frequency signal using traditional OMP algorithms.
In contrast, the frequency clustering results of the signal S reconstructed by the DTM_OMP_ICA algorithm proposed in this paper show that the average value of each frequency cluster is closer to the frequency of the original noise-free signal, and the variance within the four frequency clusters is near zero. Even in a strong noise interference environment with an SNR of 0 dB, the frequency clustering results of the signal reconstructed by the DTM_OMP_ICA algorithm demonstrate that the average values of the four frequency clusters are closer to the frequencies of the original noise-free signal, with the variance within each cluster being significantly low (0.07, 0, 0.0042, 0.005), as illustrated in Figure 8(c(1)). The occurrence of outlier frequencies within each clustered frequency is greatly reduced compared to Figure 8(a(1),b(1)). Similarly, under noise interference of 3 dB, 5 dB, and 9 dB, the frequency clustering results of the signal reconstructed by the DTM_OMP_ICA algorithm show that the average value within each frequency cluster matches the original noise-free signal frequency exactly, and the variance is closer to zero. The experimental results indicate that the DTM_OMP_ICA algorithm’s stability in reconstructing noisy multi-frequency signals surpasses that of traditional OMP algorithms.
Based on the observations, it can be concluded that the DTW_OMP_ICA algorithm proposed in this paper exhibited superior noise resistance and stability compared to traditional OMP algorithms. In noisy environments, it excelled at reconstructing the original signal frequencies more accurately.
(iii)
Time Consumption
In order to exclude the effect of computer hardware chance on the running time of the algorithm, the experiment was conducted with 10 sets of 100 rounds of independent repetitive experiments with the same settings as before. The computer was configured with CPU: AMD Ryzen 7 5800H (8 cores, 16 threads, 3.2 GHz); RAM: DDR4-3200 16 GB; and the final algorithmic runtime results were averaged over the 10 groups of recorded durations, and counted using the durations of a single round, ten rounds, fifty rounds, and one hundred rounds of different algorithmic-based runs, the results of which are shown in Table 2.
Table 2 shows the statistics of the length of time the algorithms were run independently, where the first column shows the three different algorithms used, and columns 2, 3, 4, and 5 show the time consumed by each algorithm in different rounds of independently repeated experiments.
Table 2 indicates that the RGM-OMP algorithm had the lowest time cost for reconstructing noisy multi-frequency signals, followed by the SRM-OMP algorithm, but the difference is minimal. The proposed algorithm, which employs both RGM-OMP and SRM-OMP to construct a dual-path process and combines them with ICA to find the most relevant sparse representation of the clean signal in the noisy signal, incurred the highest time cost. However, it was still less than the sum of the time taken by the RGM-OMP and SRM-OMP algorithms.
While the DTM_OMP_ICA algorithm incurred a time cost approximately 1.7 times that of the traditional OMP algorithm, considering its overall performance in denoising and reconstructing noisy multi-frequency signals, the DTM_OMP_ICA algorithm still maintained a certain advantage over the traditional OMP algorithm.

4.2. Array Antenna Signal Denoising Experiment

To verify the denoising and reconstruction performance of the proposed DTM_OMP_ICA in practical applications, this paper applies it to the denoising process of signals received by an array antenna. Furthermore, to comprehensively evaluate the performance of our algorithm in real-world scenarios, this paper also conducts a comparative analysis using two orthogonal matching pursuit algorithms based on different observation bases—RGM-OMP and SRM-OMP—as well as two classic signal denoising techniques—the singular value decomposition denoising algorithm [50] and the wavelet threshold denoising algorithm [51].

4.2.1. Data

The array antenna signal data used in this paper were derived from the simulated data provided in Problem A of the 2022 China Postgraduate Mathematical Modeling Competition. The data were generated from the echo signals received by a simulated FMCW (frequency-modulated continuous-wave) radar with 86 array antennas in a noisy environment. The specific parameters of the signals received by the array antenna are shown in Table 3.

4.2.2. Results of Simulated Radar Signal Processing

This paper employed the proposed DTM_OMP_ICA algorithm, along with the RGM-OMP algorithm, SRM-OMP algorithm, SVD decomposition denoising algorithm, and wavelet threshold denoising algorithm used for comparative experiments, and applied them to denoise and reconstruct the noisy echo signals received by 86 array antennas in a simulated FMCW radar system. The reconstructed signals were then subjected to Fourier transform for frequency domain analysis, and the results are presented in Figure 9.
Figure 9a displays the spectral graph of the noisy echo signal S1 received by 86 array antennas. Subsequently, Figure 9b–f each present the spectral graphs of signals processed through various denoising algorithms, including RMG-OMP, SRM-OMP, DTM_OMP_ICA, SVD decomposition denoising, and wavelet threshold denoising. The signals reconstructed by these algorithms are, respectively, denoted as S3, S4, S2, S5, and S6.
Figure 9 indicates that, after denoising and reconstructing the noisy signals received by the array antennas using the DTM_OMP_ICA algorithm proposed in this paper, the reconstructed signals effectively reduced the majority of interference from noise while preserving the original frequency characteristics of the signals. Since the positioning performance of FMCW radar systems relies on the effective frequencies of received signals, the DTM_OMP_ICA algorithm could provide assistance in achieving precise localization of target objects for radar systems in noise-disturbed environments.

5. Conclusions and Discussion

This paper proposed the DTM_OMP_ICA algorithm to validate and eliminate the potential negative impact resulting from the strong randomness of the observation matrix used in the OMP algorithm. Unlike traditional OMP algorithms, the DTM_OMP_ICA algorithm does not directly reconstruct the signal. Instead, it generates a denoising mask matrix to perform denoising on the noisy signal. This approach effectively addresses the limitations of traditional OMP algorithms in noise resistance. The reconstructed signal, while preserving the original signal’s characteristics, reduces interference from most of the noise, making this algorithm significantly superior when reconstructing noisy signals.
This paper conducted comparative experiments by designing the same simulation experiments and contrasting them with the OMP reconstruction algorithm based on two different observation matrices. The experimental results demonstrate that the proposed method exhibits outstanding noise resistance, with a maximum noise resistance performance of up to 98%. Additionally, this algorithm showcases excellent stability, with a maximum stability performance of up to 99%, which is significantly better than traditional OMP algorithms. Admittedly, the above superior performance comes at the cost of time: the time consumption of the DTM_OMP_ICA algorithm is about 1.7 times that of two traditional algorithms. In practical applications, this method provides an innovative approach to the problem of reconstructing noisy signals and holds promising prospects in the relevant fields.
Additionally, the experiments found that the threshold setting for the mean mask in the proposed algorithm has an impact on the experimental results. When the threshold is set from 5 times to 30 times the mean energy of the reconstructed signal’s frequencies for constructing the mask matrix, the algorithm demonstrates good reconstruction performance. However, setting a threshold exceeding 50 times the mean energy of the reconstructed signal’s frequencies may result in the loss of some frequency features of the signal itself. On the other hand, setting a threshold less than five times the mean energy of the reconstructed signal’s frequencies for the mask matrix may lead to reduced noise resistance when reconstructing noisy signals with a high SNR. Therefore, effectively setting the mask threshold will be studied in our future work.

Author Contributions

Conceptualization, M.G.; methodology, M.G., D.L. and Y.Z.; software, M.G.; validation, M.G. and Y.Y.; formal analysis, W.L.; investigation, R.X.; resources, L.G.; data curation, Z.W.; writing—original draft preparation, M.G.; writing—review and editing, D.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Agricultural Joint Fund of Yunnan Province, China under grant no. 202301BD070001-086; the Scientific Research Foundation of the Education Department of Yunnan Province, China under grant no. 2022J0495; the National Natural Science Foundation of China under grant no. 31860332; the National Natural Science Foundation of China under grant no. 32360388; and Research on the Application of Multi-Target Swarm Intelligence Algorithms with the Multi-Modal in Biological Data.

Data Availability Statement

The data used in this paper were sourced from the simulated intermediate-frequency (IF) noisy signals of the 2022 China Graduate Mathematical Modeling Competition, Problem A.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Ausherman, D.A.; Kozma, A.; Walker, J.L.; Jones, H.M.; Poggio, E.C. Developments in radar imaging. IEEE Trans. Aerosp. Electron. Syst. 1984, AES-20, 363–400. [Google Scholar]
  2. Bahl, P.; Padmanabhan, V.N. RADAR: An in-building RF-based user location and tracking system. In Proceedings of the IEEE INFOCOM 2000. Conference on Computer Communications. Nineteenth Annual Joint Conference of the IEEE Computer and Communications Societies (Cat. No. 00CH37064), Tel Aviv, Israel, 26–30 March 2000; Volume 2, pp. 775–784. [Google Scholar]
  3. Griffiths, H.D.; Baker, C.J. Passive coherent location radar systems. Part 1: Performance prediction. IEE Proc.-Radar Sonar Navig. 2005, 152, 153–159. [Google Scholar] [CrossRef]
  4. Chen, D.; Chen, B.; Qin, G. Angle estimation using ESPRIT in MIMO radar. Electron. Lett. 2008, 44, 1. [Google Scholar]
  5. Chen, J.; Hong, G.; Su, W. Angle estimation using ESPRIT without pairing in MIMO radar. Electron. Lett. 2008, 44, 1422–1423. [Google Scholar]
  6. Stove, A.G. Linear FMCW radar techniques IEE Proceedings F (Radar and Signal Processing). IET Digit. Libr. 1992, 139, 343–350. [Google Scholar]
  7. Li, Y.A.; Hung, M.H.; Huang, S.J.; Lee, J. A fully integrated 77 GHz FMCW radar system in 65 nm CMOS. In Proceedings of the 2010 IEEE International Solid-State Circuits Conference-(ISSCC), San Francisco, CA, USA, 7–11 February 2010; pp. 216–217. [Google Scholar]
  8. Brooker, G.M. Understanding millimetre wave FMCW radars. In Proceedings of the 1st International Conference on Sensing Technology, Palmerston North, New Zealand, 21–23 November 2005. [Google Scholar]
  9. Reindl, L.; Ruppel, C.C.; Berek, S.; Knauer, U.; Vossiek, M.; Heide, P.; Oréans, L. Design, fabrication, and application of precise SAW delay lines used in an FMCW radar system. IEEE Trans. Microw. Theory Tech. 2001, 49, 787–794. [Google Scholar] [CrossRef]
  10. Piper, S.O. Receiver frequency resolution for range resolution in homodyne FMCW radar. In Proceedings of the Conference Proceedings National Telesystems Conference 1993, Atlanta, GA, USA, 16–17 June 1993; pp. 169–173. [Google Scholar]
  11. Oppenheim, A.V.; Schafer, R.W.; Buck, J.R. Discrete-Time Signal Processing, 2nd ed.; Prentice Hall: Upper On River City, NJ, USA, 1999. [Google Scholar]
  12. Parks, T.W.; Burrus, C.S. Digital Filter Design; John Wiley & Sons: Hoboken, NJ, USA, 1987. [Google Scholar]
  13. Hayes, M.H. Statistical Digital Signal Processing and Modeling; John Wiley & Sons: Hoboken, NJ, USA, 1982. [Google Scholar]
  14. Candès, E.J. Compressive sampling. In Proceedings of the International Congress of Mathematicians, Madrid, Spain, 22–30 August 2006; pp. 1433–1452. [Google Scholar]
  15. Candès, E.J.; Romberg, J.; Tao, T. Robust uncertainty principles: Ex act signal reconstruction from highly incomplete frequency information. IEEE Trans. Inf. Theory 2006, 52, 489–509. [Google Scholar] [CrossRef]
  16. Donoho, D.L. Compressed sensing. IEEE Trans. Inf. Theory 2006, 52, 1289–1306. [Google Scholar] [CrossRef]
  17. Candès, E.J.; Romberg, J. Signal recovery from random projections. In Proceedings of the SPIE 5674, Computational Imaging III, San Jose, CA, USA, 11 March 2005; pp. 56–74. [Google Scholar]
  18. Donoho, D.L.; Tsaig, Y. Extensions of compressed sensing. Signal Process. 2006, 86, 533–548. [Google Scholar]
  19. Schwartz, J.; Steinberg, B. Ultrasparse, ultrawide band arrays. IEEE Trans. Ultrason. Ferro Electr. Freq. Control. 1998, 45, 376–393. [Google Scholar] [CrossRef]
  20. Wang, Y.; Wang, W.; Zhou, M.; Ren, A.; Tian, Z. Remote monitoring of human vital signs based on 77-GHz mm-wave FMCW radar. Sensors 2020, 20, 2999. [Google Scholar] [CrossRef]
  21. Wang, G.; Gu, C.; Inoue, T.; Li, C. A hybrid FMCW-interferometry radar for indoor precise positioning and versatile life activity monitoring. IEEE Trans. Microw. Theory Tech. 2014, 62, 2812–2822. [Google Scholar] [CrossRef]
  22. Herman, M.A.; Strohmer, T. High-resolution radar via compressed sensing. IEEE Trans. Signal Process. 2009, 57, 2275–2284. [Google Scholar] [CrossRef]
  23. Tropp, J.A.; Gilbert, A.C. Signal recovery from random measurements via orthogonal matching pursuit. IEEE Trans. Inf. Theory 2007, 53, 4655–4666. [Google Scholar] [CrossRef]
  24. Duarte-Carvajalino, J.M.; Sapiro, G. Learning to sense sparse signals: Simultaneous sensing matrix and sparsifying dictionary optimization. IEEE Trans. Image Process. 2009, 18, 1395–1408. [Google Scholar] [CrossRef] [PubMed]
  25. Erbe, C.; Reichmuth, C.; Cunningham, K.; Lucke, K.; Dooling, R. Communication masking in marine mammals: A review and research strategy. Mar. Pollut. Bull. 2016, 103, 15–38. [Google Scholar] [CrossRef] [PubMed]
  26. Hyvärinen, A.; Oja, E. Independent component analysis: Algorithms and applications. Neural Netw. 2000, 13, 411–430. [Google Scholar] [CrossRef]
  27. Zi, J.; Lv, D.; Liu, J.; Huang, X.; Yao, W.; Gao, M.; Xi, R.; Zhang, Y. Improved Swarm Intelligent Blind Source Separation Based on Signal Cross-Correlation. Sensors 2021, 22, 118. [Google Scholar] [CrossRef] [PubMed]
  28. Candes, E.J. The restricted isometry property and its implications for compressed sensing. Comptes Rendus. Math. 2008, 346, 589–592. [Google Scholar] [CrossRef]
  29. Baraniuk, R.; Davenport, M.; DeVore, R.; Wakin, M. A simple proof of the restricted isometry property for random matrices. Constr. Approx. 2008, 28, 253–263. [Google Scholar] [CrossRef]
  30. Rosen, J.B. The gradient projection method for nonlinear programming. Part I. Linear constraints. J. Soc. Ind. Appl. Math. 1960, 8, 181–217. [Google Scholar] [CrossRef]
  31. Rosen, J.B. The gradient projection method for nonlinear programming. Part II. Nonlinear constraints. J. Soc. Ind. Appl. Math. 1961, 9, 514–532. [Google Scholar] [CrossRef]
  32. Li, S.; Zhao, G.; Zhang, W.; Qiu, Q.; Sun, H. ISAR imaging by two-dimensional convex optimization-based compressive sensing. IEEE Sens. J. 2016, 16, 7088–7093. [Google Scholar] [CrossRef]
  33. Wang, R.; Zhang, J.; Ren, S.; Li, Q. A reducing iteration orthogonal matching pursuit algorithm for compressive sensing. Tsinghua Sci. Technol. 2016, 21, 71–79. [Google Scholar] [CrossRef]
  34. Wu, D.; Zhu, W.P.; Swamy, M.N.S. The theory of compressive sensing matching pursuit considering time-domain noise with application to speech enhancement. IEEE/ACM Trans. Audio Speech Lang. Process. 2014, 22, 682–696. [Google Scholar] [CrossRef]
  35. Do, T.T.; Gan, L.; Nguyen, N.; Tran, T.D. Sparsity adaptive matching pursuit algorithm for practical compressed sensing. In Proceedings of the 2008 42nd Asilomar Conference on Signals, Systems and Computers, Pacific Grove, CA, USA, 26–29 October 2008; pp. 581–587. [Google Scholar]
  36. Needell, D.; Vershynin, R. Signal recovery from incomplete and inaccurate measurements via regularized orthogonal matching pursuit. IEEE J. Sel. Top. Signal Process. 2010, 4, 310–316. [Google Scholar] [CrossRef]
  37. Needell, D.; Vershynin, R. Uniform uncertainty principle and signal recovery via regularized orthogonal matching pursuit. Found. Comput. Math. 2009, 9, 317–334. [Google Scholar]
  38. Yang, G.; Zhang, Q.; Que, P.W. Matching-pursuit-based adaptive wavelet-packet atomic decomposition applied in ultrasonic inspection. Russ. J. Nondestruct. Test. 2007, 43, 62–68. [Google Scholar] [CrossRef]
  39. Li, S.; Fang, L. Signal denoising with random refined orthogonal matching pursuit. IEEE Trans. Instrum. Meas. 2011, 61, 26–34. [Google Scholar] [CrossRef]
  40. Candes, E.J.; Tao, T. Near-optimal signal recovery from random projections: Universal encoding strategies? IEEE Trans. Inf. Theory 2006, 52, 5406–5425. [Google Scholar] [CrossRef]
  41. Zhang, G.; Jiao, S.; Xu, X.; Wang, L. Compressed sensing and reconstruction with bernoulli matrices. In Proceedings of the 2010 IEEE International Conference on Information and Automation, Harbin, China, 20–23 June 2010; pp. 455–460. [Google Scholar]
  42. Stewart, G.W. The efficient generation of random orthogonal matrices with an application to condition estimators. SIAM J. Numer. Anal. 1980, 17, 403–409. [Google Scholar] [CrossRef]
  43. Haupt, J.; Bajwa, W.U.; Raz, G.; Nowak, R. Toeplitz compressed sensing matrices with applications to sparse channel estimation. IEEE Trans. Inf. Theory 2010, 56, 5862–5875. [Google Scholar] [CrossRef]
  44. Indyk, P. Sparse Recovery Using Sparse Random Matrices. In Proceedings of the LATIN 2010: Theoretical Informatics, 9th Latin American Symposium, Oaxaca, Mexico, 19–23 April 2010; p. 157. [Google Scholar]
  45. Meena, V.; Abhilash, G. Robust recovery algorithm for compressed sensing in the presence of noise. IET Signal Process. 2016, 10, 227–236. [Google Scholar] [CrossRef]
  46. Thomas, T.J.; Arun, A.; Sheeba, R.J. Compressed Sensing Recovery using Modified Newton Gradient Pursuit Algorithm and its Application to ECG with Denoising. In Proceedings of the 2018 International Conference on Signal Processing and Communications (SPCOM), Bangalore, India, 16–19 July 2018; IEEE: Piscataway, NJ, USA, 2018; pp. 412–416. [Google Scholar]
  47. Li, L.; Fang, Y.; Liu, L.; Peng, H.; Kurths, J.; Yang, Y. Overview of compressed sensing: Sensing model, reconstruction algorithm, and its applications. Appl. Sci. 2020, 10, 5909. [Google Scholar] [CrossRef]
  48. Cheng, C.; Lin, D. Based on compressed sensing of orthogonal matching pursuit algorithm image recovery. J. Internet Things 2020, 2, 37–45. [Google Scholar] [CrossRef]
  49. Leys, C.; Ley, C.; Klein, O.; Bernard, P.; Licata, L. Detecting outliers: Do not use standard deviation around the mean, use absolute deviation around the median. J. Exp. Soc. Psychol. 2013, 49, 764–766. [Google Scholar] [CrossRef]
  50. Jha, S.K.; Yadava, R.D.S. Denoising by singular value decomposition and its application to electronic nose data processing. IEEE Sens. J. 2010, 11, 35–44. [Google Scholar] [CrossRef]
  51. Chang, S.G.; Yu, B.; Vetterli, M. Adaptive wavelet thresholding for image denoising and compression. IEEE Trans. Image Process. 2000, 9, 1532–1546. [Google Scholar] [CrossRef]
Figure 1. DTM_OMP_ICA.
Figure 1. DTM_OMP_ICA.
Sensors 24 02291 g001
Figure 2. Time−domain waveform plots of the multi−frequency signal. (a) Original clean multi−frequency signal x . (b) Noisy multi−frequency signal x at 0 dB SNR. (c) Noisy multi−frequency signal x at 3 dB SNR. (d) Noisy multi−frequency signal x at 5 dB SNR. (e) Noisy multi−frequency signal x at 9 dB SNR.
Figure 2. Time−domain waveform plots of the multi−frequency signal. (a) Original clean multi−frequency signal x . (b) Noisy multi−frequency signal x at 0 dB SNR. (c) Noisy multi−frequency signal x at 3 dB SNR. (d) Noisy multi−frequency signal x at 5 dB SNR. (e) Noisy multi−frequency signal x at 9 dB SNR.
Sensors 24 02291 g002
Figure 3. The spectrum of signals reconstructed by different algorithms. (a(1)a(5)) Spectrum of the original multi-frequency signal x and its noisy counterparts at SNRs of 0 dB, 3 dB, 5 dB, and 9 dB, respectively. (b(1)b(5)) Spectrum of the signal reconstructed by RGM-OMP under the same SNRs. (c(1)c(5)) Spectrum of the signal reconstructed by SRM-OMP under the same SNRs. (d(1)d(5)) Spectrum of the signal reconstructed by DTM_OMP_ICA under the same SNRs.
Figure 3. The spectrum of signals reconstructed by different algorithms. (a(1)a(5)) Spectrum of the original multi-frequency signal x and its noisy counterparts at SNRs of 0 dB, 3 dB, 5 dB, and 9 dB, respectively. (b(1)b(5)) Spectrum of the signal reconstructed by RGM-OMP under the same SNRs. (c(1)c(5)) Spectrum of the signal reconstructed by SRM-OMP under the same SNRs. (d(1)d(5)) Spectrum of the signal reconstructed by DTM_OMP_ICA under the same SNRs.
Sensors 24 02291 g003
Figure 4. Reconstructed frequency statistical parameter A c c F . (a) SNR = 0. (b) SNR = 3. (c) SNR = 5. (d) SNR = 9.
Figure 4. Reconstructed frequency statistical parameter A c c F . (a) SNR = 0. (b) SNR = 3. (c) SNR = 5. (d) SNR = 9.
Sensors 24 02291 g004
Figure 5. Reconstructed frequency statistical parameter E r r F . (a) SNR = 0. (b) SNR = 3. (c) SNR = 5. (d) SNR = 9.
Figure 5. Reconstructed frequency statistical parameter E r r F . (a) SNR = 0. (b) SNR = 3. (c) SNR = 5. (d) SNR = 9.
Sensors 24 02291 g005
Figure 6. Reconstructed frequency statistical parameter A c c F A . (a) SNR = 0. (b) SNR = 3. (c) SNR = 5. (d) SNR = 9.
Figure 6. Reconstructed frequency statistical parameter A c c F A . (a) SNR = 0. (b) SNR = 3. (c) SNR = 5. (d) SNR = 9.
Sensors 24 02291 g006
Figure 7. Overall analysis of reconstructed frequency statistical parameters for three algorithms.
Figure 7. Overall analysis of reconstructed frequency statistical parameters for three algorithms.
Sensors 24 02291 g007
Figure 8. Frequency clustering statistical results. (a(1)a(4)) Frequency clustering results of the reconstructed signal S using RGM-OMP at SNRs of 0 dB, 3 dB, 5 dB, and 9 dB, respectively. (b(1)b(4)) Frequency clustering results of the reconstructed signal S using SRM-OMP at the same SNRs. (c(1)c(4)) Frequency clustering results of the reconstructed signal S using DTM_OMP_ICA at the corresponding SNRs.
Figure 8. Frequency clustering statistical results. (a(1)a(4)) Frequency clustering results of the reconstructed signal S using RGM-OMP at SNRs of 0 dB, 3 dB, 5 dB, and 9 dB, respectively. (b(1)b(4)) Frequency clustering results of the reconstructed signal S using SRM-OMP at the same SNRs. (c(1)c(4)) Frequency clustering results of the reconstructed signal S using DTM_OMP_ICA at the corresponding SNRs.
Sensors 24 02291 g008
Figure 9. Spectrum of signals before and after reconstruction. (a) Original signal. (b) RGM-OMP. (c) SRM-OMP. (d) DTM_OMP_ICA. (e) SVD denoising. (f) Wavelet thresholding denoising.
Figure 9. Spectrum of signals before and after reconstruction. (a) Original signal. (b) RGM-OMP. (c) SRM-OMP. (d) DTM_OMP_ICA. (e) SVD denoising. (f) Wavelet thresholding denoising.
Sensors 24 02291 g009
Table 1. Initial parameter settings for the experiment.
Table 1. Initial parameter settings for the experiment.
ParameterValueIntroduction
K8The sparsity of the signal in its orthogonal transformation basis (since the Fourier transform is two-sided, the sparsity is 4 × 2).
N1024The original number of sampled points in the signal.
t128 (unit: ms)Duration of the Signal.
fs8000 (unit: Hz)The signal’s sampling frequency based on the Nyquist sampling theorem will be used for subsequent frequency domain analysis of the reconstructed signal.
Table 2. Algorithm runtime summary.
Table 2. Algorithm runtime summary.
RoundRound 1Round 10Round 50Round 100
Algorithm
RGM-OMP45.41 (ms)227.62 (ms)880.43 (ms)1839.89 (ms)
SRM-OMP46.43 (ms)277.01 (ms)926.09 (ms)1915.21 (ms)
DTM_OMP_ICA91.14 (ms)414.41 (ms)1446.20 (ms)2785.59 (ms)
Table 3. Parameters and significance of signals received by array antennas.
Table 3. Parameters and significance of signals received by array antennas.
VariablesSignificance
T = 3.2 × 10 5   ( s ) Represents the chirp period, pulse time
T s = 1.25 × 10 7   ( s ) Sampling interval
N f = 32 N f represents the number of chirp cycles
L = 0.0815   ( m ) The aperture of the array antennas
γ = 78.986 × 10 12   ( H z / s ) γ represents the frequency modulation slope
N a = 86 The number of equivalent virtual antenna arrays
C = 3 × 10 8   ( m / s ) Speed of light
f 0 = 78.8 × 10 9   ( H z ) Represents the carrier frequency
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

Gao, M.; Zhang, Y.; Yu, Y.; Lv, D.; Xi, R.; Li, W.; Gu, L.; Wang, Z. An Improved OMP Algorithm for Enhancing the Anti-Interference Performance of Array Antennas. Sensors 2024, 24, 2291. https://doi.org/10.3390/s24072291

AMA Style

Gao M, Zhang Y, Yu Y, Lv D, Xi R, Li W, Gu L, Wang Z. An Improved OMP Algorithm for Enhancing the Anti-Interference Performance of Array Antennas. Sensors. 2024; 24(7):2291. https://doi.org/10.3390/s24072291

Chicago/Turabian Style

Gao, Mingyuan, Yan Zhang, Yueyun Yu, Danju Lv, Rui Xi, Wei Li, Lianglian Gu, and Ziqian Wang. 2024. "An Improved OMP Algorithm for Enhancing the Anti-Interference Performance of Array Antennas" Sensors 24, no. 7: 2291. https://doi.org/10.3390/s24072291

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