Next Article in Journal
A Novel Single-Sample Retinal Vessel Segmentation Method Based on Grey Relational Analysis
Previous Article in Journal
LMI-Based MPC Design Applied to the Single-Phase PWM Inverter with LC Filter under Uncertain Parameters
Previous Article in Special Issue
A Novel Generalized Nested Array MIMO Radar for DOA Estimation with Increased Degrees of Freedom and Low Mutual Coupling
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Communication

Efficient Aperture Fill Time Correction for Wideband Sparse Array Using Improved Variable Fractional Delay Filters

1
National Key Laboratory of Electromagnetic Space Security, Chengdu 610036, China
2
College of Information Science and Engineering, Hohai University, Nanjing 210098, China
*
Author to whom correspondence should be addressed.
Sensors 2024, 24(13), 4327; https://doi.org/10.3390/s24134327
Submission received: 31 May 2024 / Revised: 20 June 2024 / Accepted: 2 July 2024 / Published: 3 July 2024

Abstract

:
To solve the problem of aperture fill time (AFT) for wideband sparse arrays, variable fractional delay (VFD) FIR filters are applied to eliminate linear coupling between spatial and time domains. However, the large dimensions of the filter coefficient matrix result in high system complexity. To alleviate the computational burden of solving VFD filter coefficients, a novel multi–regultion minimax (MRMM) model utilizing the sparse representation technique has been presented. The error function is constrained by the introduction of L2–norm and L1–norm regularizations within the minimax criterion. The L2–norm effectively resolves the problems of overfitting and non–unique solutions that arise in the sparse optimization of traditional minimax (MM) models. Meanwhile, the use of multiple L1–norms enables the optimal design of the smallest sub–filter number and order of the VFD filter. To solve the established nonconvex model, an improved sequential–alternating direction method of multipliers (S–ADMM) algorithm for filter coefficients is proposed, which utilizes sequential alternation to iteratively update multiple soft–thresholding problems. The experimental results show that the optimized VFD filter reduces system complexity significantly and corrects AFT effectively in a wideband sparse array.

1. Introduction

A wideband sparse array has the advantages of improving spatial resolution, reducing hardware complexity, and improving target detection and recognition performance. However, the aperture fill time (AFT) phenomenon, inherent to wideband arrays, significantly diminishes the performance of beamforming. Variable fractional delay (VFD) filters involve the implementation of precise variable delay in digital systems and are widely used for aperture correction in wideband arrays [1,2]. High–precision delay filters are usually represented as a sum of integer and fractional terms. While cascaded registers allow for the implementation of integer terms, fractional terms require more sophisticated techniques [3,4]. The pioneering Farrow structure [5], known for its flexible and efficient delay correction, modifies and extends the traditional structure, effectively reducing implementation complexity and improving design accuracy.
In general, the design of a Farrow structure VFD (FS–VFD) filter primarily involves the minimax method [6,7] and the weighted least squares (WLS) method [8,9]. The minimax method is implemented by minimizing the maximum magnitude value between the actual frequency response and the ideal frequency response, ensuring accuracy within the approximation error range [7]. Some studies have transformed the solution problem into convex optimization problems in order to ensure the optimal value of the coefficient solution; it still requires numerous iterations and updates, leading to high computational complexity. The WLS design is implemented by minimizing the mean square error of the two frequencies, yielding an optimal closed solution [8]. However, the higher–order filter design leads to a cascading increase in computational complexity for both design criterion, necessitating additional storage resources [9]. High computational complexity and unstable values affect the efficiency of FS–VFD filters, so an ideal design approach needs to ensure accuracy while reducing computational complexity.
FS–VFD FIR filters comprise several sub–filters with a large number of coefficients, distinguishing them from the numerical scalars in digital FIR filters [10]. To reduce computation, the literature [11] reduces the number of filter coefficients using the WLS design criterion and the (inverse) symmetry of the coefficients but does not take into account the change in system complexity. Filter coefficient sparse processing is an optimization technique aimed at reducing system complexity and improving processing efficiency by reducing the number of non–zero coefficients in a filter [12,13]. Through sparse optimization, we can obtain a more compact and efficient filter structure, which is especially important for resource–constrained environments or application scenarios requiring real–time performance. The literature [14] studies VFDs with sparse coefficients for reducing computation. Leveraging the polynomial properties of the coefficients, a two–stage technique is proposed for a VFD FIR filter, aimed at obtaining a least squares solution under sparse target coefficients. To reduce system complexity and computation of the Farrow structure, the literature [15] considers sparse–optimized Farrow structure VFD (SFS–VFD) filters based on least squares and coefficient (inverse) symmetry to ensure high accuracy and efficient correction of aperture effects. It is noteworthy that while much research focuses on implementing FS–VFD filters with the WLS criterion, there remains a gap in related research based on the minimax criterion [16].
Therefore, this paper proposes a new multi–regulation minimax sparse model for VFD design. The error function is constrained by introducing L2–norm and L1–norm regularizations under the minimax criterion. The L2–norm solves the overfitting and non–unique solution problems of the traditional MM model in sparse optimization, and the L1–norm achieves the optimal design of the smallest sub–filter number and order for VFD filters. An improved sequential ADMM algorithm is proposed for solving this nonconvex model. Subsequently, the proposed MRMM–VFD filters are applied to a wideband sparse array to correct the AFT. Finally, simulation results verify the effectiveness and feasibility of MRMM–VFD filters.

2. Problem Formulation

The ideal frequency response of a VFD filter is expressed as follows:
H d ( ω , p ) = e j ω p  
where ω [ 0 , α π ] is the passband range, and 0 < α 1 , p [ 0.5 , 0.5 ] is the fractional delay parameter range. The block diagram of the VFD filter structure based on the Farrow structure is shown in Figure 1. Therefore, the actual filter frequency response can be expressed as follows:
H ( ω , p ) = n = N N + 1 m = 0 M a ( n , m ) p m e j ω n
where a ( n , m ) are the filter coefficients; M is the number of sub–filters; and N is the sub–filter order. Since VFD filters have coefficient (anti)symmetry, taking an odd number N as an example, the coefficients of the even sub–filter (with symmetry) and the odd sub–filter (with anti–symmetry) are split to obtain the following form:
H ( ω , p ) = m = 0 M e c m T b e m p 2 m j m = 1 M o s m T b o m p 2 m 1
where
c m T = [ cos ( 1 2 ω ) cos ( 3 2 ω ) cos ( ( N + 1 2 ) ω ) ] t s m T = [ sin ( 1 2 ω )   sin ( 3 2 ω ) sin ( ( N + 1 2 ) ω ) ] t
Convert it into the following form:
H ( ω , p ) = f T b e j g T b o
where
f m T = c m T p 2 m , g m T = s m T p 2 m 1
f T = f 0 T f 1 T     f I T g T = g 1 T g 2 T     g I + 1 T
The minimax design of a VFD filter refers to minimizing the peak error between the ideal frequency response and the actual frequency response. The solution problem for the filter coefficients of the minimax model is given as follows:
minimize   max | H ( ω , p ) H d ( ω , p ) |
The error function is defined as follows:
e H ( ω , p ) = H ( ω , p ) H d ( ω , p ) = e R ( ω , p ) j e I ( ω , p )
where
e R ( ω , p ) = f b e cos ( ω p ) e I ( ω , p ) = g b o sin ( ω p )
Then, problem (8) can be transformed into the following form:
minimize   f b e cos ( ω p ) + j ( g b o sin ( ω p ) )
where | |   ·   | | denotes the L norm.

3. Improved Sparse–Constrained Design for the VFD FIR Filter

3.1. The Cost Function Based on the Multi–Regulation Minimax Criterion

Building upon the traditional MM model established in the above section, SDP and SOCP algorithms are usually used to optimize the solution [17] but can easily lead to high system complexity. Owing to the characteristics of VFD filters with more near–zero coefficients, sparse constraint theory is introduced into the optimization of filter weight coefficients, and an improved MM model, the multi–regulation minimax (MRMM) model, is proposed in this paper. The model is introduced with L2 and L1 regularization to constrain the model coefficients. The L2–norm solves the overfitting and non–unique solution problems of the traditional MM model in sparse optimization [18]. The L1–norm realizes the design of the minimum number of sub–filters and orders of the VFD filter. The MRMM model is formulated as follows:
minimize Φ x ψ + y 2 2 + λ 1 n = 0 ( N + 1 ) ( M + 1 ) 1 k = 0 n x ( N + 1 ) ( M + 1 ) k 1 + λ 2 n = 1 N + 1 m = 0 M k = 1 n x ( N + 1 k ) + M ( N + 1 ) 1
where the filter coefficient vector x = [ b e , b o ] T and the parameters Φ and ψ in the term | |   ·   | | are related to Equation (11) in Section 2. According to reference [19], the diagonal elements of H are associated with the symmetry of the coefficients of the filter. y is the vector associated with the real impulse response a ( n , m ) of the filter. According to coefficient symmetry b ( n , m ) = 2 a ( n , m ) ,   for   m = 0 , 1 , , M , we define that the diagonal elements of H are h 1 = h 2 = = h ( N + 1 ) ( M + 1 ) = 1 / 2 , and the diagonal matrix H can be reduced to a multiple of the unit matrix I ( N + 1 ) ( M + 1 ) . In the L term, we give
Φ = f 0 0 g ,   ψ = cos ( ω p ) sin ( ω p )
The first and second terms in problem (12) denote the minimax error of the VFD filter; the third term denotes the optimal filter order constraint; and the fourth term denotes the optimal number of sub–filters constraint. The last two sparse terms should essentially be formulated as L0–norm problems. However, they cannot be solved directly due to NP Hard and a highly nonconvex nature. Therefore, the last two sparse terms are converted to L1–norm form after relaxation [20,21]. It can be seen that the orders of all sub–filters are optimized simultaneously by the second L1–norm term, which has a complex nested summation process and cannot be solved directly. Based on the properties of paradigm triangular inequality | | i b i | | i | | b i | | [22], problem (12) can be approximated as follows:
minimize Φ x ψ + x T H x + λ 1 n = 0 ( N + 1 ) ( M + 1 ) 1 k = 0 n x ( N + 1 ) ( M + 1 ) k + λ 2 n = 1 N + 1 k = 1 n x ( N + 1 k ) + + x ( ( N + 1 k ) + M ( N + 1 ) )
A further derivation gives the following:
minimize Φ x ψ + x T H x + λ 1 n = 0 ( N + 1 ) ( M + 1 ) 1 ( ( N + 1 ) ( M + 1 ) n ) x ( N + 1 ) ( M + 1 ) n + λ 2 n = 1 ( M + 1 ) ( N + 1 ) r e m [ n ( N + 1 ) ] x n
The last two items are integrated and reformulated in L1–norm form as follows:
minimize Φ x ψ + x T H x + λ 1 Q x 1 + λ 2 S x 1
Q and S in the formula are, respectively, as follows:
Q = d i a g ( 1 , 2 , , ( N + 1 ) ( M + 1 ) ) s i = d i a g ( 1 , 2 , , N + 1 ) S = d i a g ( s 1 , s 2 , , s M + 1 )

3.2. S–ADMM Algorithm for Computing VFD Filter Coefficients

Since problem (16) needs to solve the optimization problem of the minimum sum of different paradigms, and the solution variables are multi–dimensional and complex, it is difficult for general optimization algorithms or tools to find the optimal solution. The ADMM algorithm combines the decomposability of the dual ascent method with the superior convergence of the multiplier method, making it suitable for solving large–scale distributed optimization problems.
The core idea of the ADMM algorithm is to decompose the original high–dimensional optimization problem into several low–dimensional subproblems, and gradually approach the optimal solution by solving these subproblems with alternating iterations. This decomposition and iteration strategy enables the algorithm to effectively deal with high–dimensional data and has better convergence and solving efficiency. However, the standard ADMM algorithm is only applicable to the problem of minimizing the sum of two norms. Therefore, in order to solve problem (3), the ADMM algorithm was improved and a sequential ADMM algorithm is proposed.
By introducing the auxiliary variables z 1 , z 2 , z 3 , problem (16) is expressed as a minimum optimization problem with constraints as follows:
minimize   z 1 ψ + 1 2 x T H x + λ 1 z 2 1 + λ 2 z 3 1 s u b j e c t   t o   Φ Q S x + I 0 0 z 1 + 0 I 0 z 2 + 0 0 I z 3 = 0
The variable x of the original problem is decomposed by the S–ADMM algorithm into multiple variables for sequential alternating iterations. Through constructing the augmented Lagrange function and updating the independent variables and the Lagrange multipliers, the iterations are repeated and stopped until the residuals reach a predetermined accuracy, thereby determining the optimal values of the variables. By defining the Lagrange multiplier vector u = [ u 1 , u 2 , u 3 ] , the following iterative solution formula is obtained using the dual gradient descent method:
x ¯ k + 1 = arg min x L ρ ( x , z 1 k , z 2 k , z 3 k , u k ) = arg min x   1 2 x T H x + ρ 1 2 Φ x ψ z 1 k + u 1 k 2 2 + ρ 2 2 Q x z 2 k + u 2 k 2 2 + ρ 3 2 S x z 3 k + u 3 k 2 2 x k + 1 = β x ¯ k + ( 1 β x k ) z 1 k + 1 = arg min z 1 L ρ ( x k , z 1 , u 1 k ) = arg min z 1   z 1 ψ + ρ 1 2 Φ x k z 1 + u 1 k 2 2 z 2 k + 1 = arg min z 2 L ρ ( x k , z 2 , u 2 k ) = arg min x   λ 1 z 2 1 + ρ 2 2 Q x k z 2 + u 2 k 2 2 z 3 k + 1 = arg min z 3 L ρ ( x k , z 3 , u 3 k ) = arg min x   λ 2 z 3 1 + ρ 3 2 S x k z 3 + u 3 k 2 2
where β takes the value in the range of (0, 1]. It is the core correction parameter that ensures the convergence of the algorithm [23]. It is not difficult to find that solving for x is a quadratic optimization problem, while solving for the auxiliary variable z 1 , z 2 , z 3 is a multi–soft–threshold optimization problem. Multiple variables are iteratively updated in an alternating sequence and constraints are taken into account using the Lagrange multiplier u. The iterative formula for the Lagrange multiplier vector u is as follows:
u 1 k + 1 = u 1 k + ρ 1 ( Φ x k z 1 k + 1 ) u 2 k + 1 = u 2 k + ρ 2 ( Q x k z 2 k + 1 ) u 3 k + 1 = u 3 k + ρ 3 ( S x k z 3 k + 1 )
The algorithm sets a predefined termination condition x k + 1 x k 2 ( N + 1 ) ( M + 1 ) μ and determines whether the convergence criterion is met, where μ is the threshold value. If it is satisfied, the iteration is stopped; the optimal solution is found; and the final coefficients are obtained by using the symmetry of the coefficients; otherwise, the iteration is continued.
The coefficient vector x, i.e., the corresponding even sub–filter coefficients b e and odd sub–filter coefficients b o , is computed according to the S–ADMM algorithm. The complete filter coefficient matrix is obtained using the following coefficient (inverse) symmetry formula:
b ( n , m ) = 2 a ( n , m ) ,   for   m = 0 ,   1 , ,   M

4. Wideband Sparse Array AFT Correction Using MRMM–VFD

Sparse array design strives to optimally deploy sensors to achieve desirable beamforming characteristics, reduce system hardware costs, and lower computational complexity. In wideband signal models, sparse arrays can avoid conflicts caused by frequency spread of the signal. Reference [24] frames the design problem as optimally selecting P sensors out of N possible equally spaced locations (P << N), and a wideband sparse array design method based on the tapped delay line (TDL) structure was proposed.
Assume that all of the sensors are omnidirectional with the same response, and the signals impinge upon the array from the far field. The beamformer output y [ n ] is a linear combination of differently delayed versions of the received array signals x m [ n ] , m = 0 , , M 1 . The distance from the zeroth sensor to the subsequent sensor is denoted by d m for m = 0 , , M 1 , with d 0 = 0 . And the incident signal arrives at an angle θ .
The steering vector of the array as a function of the normalized frequency Ω = w T s and the arrival angle θ is as follows:
minimize Φ x ψ + x T H x + λ 1 n = 0 ( N + 1 ) ( M + 1 ) 1 k = 0 n x ( N + 1 ) ( M + 1 ) k + λ 2 n = 1 N + 1 k = 1 n x ( N + 1 k ) + + x ( ( N + 1 k ) + M ( N + 1 ) )
where μ m = d m c T s for m = 0 ,   1 , ,   M 1 , and T indicates a transpose operation. The response of the array is then given by the following:
P ( Ω , θ ) = w H s ( Ω , θ )
where w H is the Hermitian transpose of the weight vector of the array, which is given by the following:
w = w 0 T w 1 T w M 1 T T
w m = w m , 0 w m , 1 w m , J 1 T
Then, a compressive sensing sparse algorithm is used to match the array response with the desired/reference response, resulting in a sparse wideband array distribution [24]. For the obtained wideband sparse array, the MRMM–VFD filter is employed for aperture crossing correction, with the specific operational diagram shown in Figure 2.
The MRMM–VFD filter proposed in this paper has coefficient invariance, indicating that only the delay input needs to be altered while the filter coefficients remain unchanged when processing with different time delay parameters. In the aperture fill time correction of wideband sparse arrays, the MRMM–VFD greatly reduces computational complexity and system complexity.

5. Experimental Analysis

To verify the effectiveness of the designed MRMM sparse algorithm, the error results for different parameter settings are given in this section. Regarding the filter order N, a higher value results in a steeper frequency response of the filter, allowing for more accurate delay performance. However, excessively high orders may induce system instability, leading to increased noise and distortion; here, the value of N was chosen to be 30. The number of sub–filters M was set to 5, 7, 9, and 11. Regarding the value of the number of sub–filters, the amplitude response error decreases with an increase in the number of sub–filters. However, when it increases to a certain amount, the error reduction is no longer obvious. Additionally, an excessively large number of sub–filters can lead to unnecessary group phase delay errors [11,25]. The cutoff frequency parameter α = 1 , i.e., the full–pass bandwidth range was chosen for the experiment. In the MRMM algorithm, the penalty factors λ 1 and λ 2 are two important regularization parameters, which control the trade–off between data fitting and sparsity. We kept adjusting these two parameters with a sampling interval of 0.01, and the designed MRMM–VFD can achieve better performance when λ 1 = 0.11 and λ 2 = 0.27 . To evaluate the error accuracy of the filter, the maximum magnitude error, the normalized mean square error (NMSE), and the maximum group delay error were used as evaluation metrics [25], and the expressions are shown, respectively, as follows:
ε 1 = max { 20 log H ( ω , p ) H d ( ω , p ) }
ε 2 = 0 α π 0.5 0.5 | H ( ω , p ) H d ( ω , p ) | 2 d p d ω 0 α π 0.5 0.5 | H d ( ω , p ) | 2 d p d ω 1 2
ε p = max { | τ ( ω , p ) p | }
Firstly, Figure 3a shows the convergence of primal and pairwise residuals in the S–ADMM algorithm, where the dotted lines are the set threshold; both of which can converge with an accuracy of 10−5. The algorithm sets the maximum number of iterations to 100 but as can be seen in Figure 3b, the minimum of function (16) is actually found after 20 iterations, with a clear and fast trend. The results of the standard MM algorithm and the weighted least squares sparse algorithm are compared to demonstrate the effectiveness of the MRMM sparse model in this paper, and the error results of different algorithms are given in Table 1. When M = 5, the maximum magnitude error and the normalized mean square error of the MRMM sparse algorithm are larger than the standard MM algorithm, but as the number of sub–filters increases, the maximum magnitude error gradually approaches the standard MM algorithm. The normalized mean square error of the MRMM algorithm is smaller compared to the S–WLS sparse algorithm for different numbers of sub–filters. In addition, the more the number of sub–filters, the smaller the filter error, but too many sub–filters will make the system redundant, so M = 9 is chosen for practical applications. The amplitude–frequency response error plots of the VFD filters using the MRMM algorithm, the standard MM algorithm, and the S–WLS sparse algorithm for the N = 30, M = 9 condition are given as shown in Figure 4. In comparison, the MRMM algorithm is flatter than the S–WLS sparse algorithm. The maximum magnitude error is −19.81 dB for the MRMM algorithm and −18.71 dB for the S–WLS sparse algorithm. The amplitude error is reduced by 1.0 dB compared to the un–sparse–processed MM algorithm. Therefore, it can be concluded that the proposed MRMM sparse algorithm is highly effective with high error accuracy.
Next, the amplitude–frequency response error and phase–frequency response error of the designed VFD filter were analyzed. Experimental results are given for the number of sub–filters M = 9 and the initial order of sub–filters N = 30. The magnitude–frequency response plots of different algorithms are given as shown in Figure 5, where it can be seen that the normalized magnitude of the un–sparse standard MM algorithm is close to one between the band range [0,0.9], and the normalized magnitude of the sparse–processed S–WLS algorithm and the MRMM algorithm fluctuates around one between the band range [0,0.9], but the range of fluctuation is within acceptable limits. In the bandwidth range [0.9,1], the maximum amplitude fluctuation of the MRMM algorithm is only 0.05, while the maximum amplitude fluctuation of the S–WLS algorithm in this range is 0.12. Next, Figure 6 shows the phase–frequency curves of the different algorithms where fractional delay p = 0.3; the maximum phase fluctuation of the S–WLS algorithm is 0.004; and the maximum phase fluctuation of the MRMM algorithm is only 0.002. Therefore, it can be proved that the MRMM algorithm proposed in this paper has more stable amplitude–frequency characteristics than the existing sparse algorithms.
To evaluate the sparsity of the MRMM algorithm and the effect of the sparse design on system complexity, the sparsity rate is defined as η = P s P × 100 % where Ps represents the number of zero coefficients after being sparse and P = 2 ( N + 1 ) ( M + 1 ) represents the number of filter coefficients under the initial setting. The experimental parameters are M = 9, N = 30, and P = 620. The results of the sparsity rate η and operation time T c comparison of different algorithms are given in Table 2. The data in the table show that at M = 5 and M = 7, the sparsity rate of the MRMM algorithm is 66.67% and 75.40%, while the sparsity rate of the S–WLS algorithm is 69.89% and 75.80%, respectively. It can be seen that the sparsity of the S–WLS algorithm is a little bit better when the number of sub–filters is small, and the MRMM algorithm is better when the number of sub–filters is large. However, the sparsity of the sparse algorithms MRMM and S–WLS is much greater than that of the standard MM algorithm, regardless of the experimental conditions. In terms of computation time, we can see that the MRMM algorithm exhibits higher computational efficiency, especially when dealing with more sub–filters, and the MRMM algorithm is able to complete the computation faster. Therefore, it can be shown the sparsity of the MRMM algorithm designed in this paper fully satisfies our desired design.
Next, the system complexity of the VFD filter designed by the sparse MRMM algorithm is evaluated. The presence of zero coefficients at lower orders saves the use of corresponding multipliers, and the presence of zero coefficients at higher orders saves the corresponding multipliers and adders. By statistically analyzing the coefficient matrix elements, it can be determined that the filter designed by the standard MM algorithm requires 844 multipliers and 421 adders, which results in extremely high complexity for the system. However, the multipliers and adders required by the sparse MRMM algorithm are 240 and 119, corresponding to a reduction of 71.56% and 71.73%, respectively. It can be concluded that the use of multipliers and adders is significantly reduced by the MRMM sparse design algorithm, which effectively reduces system complexity.
The above experiment analysis proves that the improved VFD filter based on the MRMM model can achieve our purpose well. We then used it to correct the aperture effect in a wideband array. Assume that the angle of incidence of the signal is 10°; the center frequency is 1.2 GHz; and the f H : f L = 3 : 1 . The uniformly distributed array element spacing is 0.5 λ min , and the number of full arrays is 30. Then, the number of arrays after sparse reconstruction is 20, as shown in Figure 7, which is designed according to [24].
Compared with the conventional wideband aperture effect correction for a TDL structure [24], our proposed MRMM–VFD–based wideband aperture effect correction method can keep the filter coefficients unchanged under different scan angles, and only the delay parameter needs to be adjusted. Compared with the LCMM–VFD filter, our proposed MRMM–VFD can significantly reduce the use of multipliers and adders and decrease system complexity. Figure 8a shows the aperture crossing phenomenon of a wideband sparse array, and it can be seen that serious pointing deviations occur in the received signals of different arrays. Figure 8b shows that the received signals of each array are aggregated to the same distance unit after aperture correction using the improved VFD filter, which completely solves beam dispersion. The antenna orientation diagrams of the wideband sparse array with different frequency signals are given in Figure 9. It can be seen above that it is effective to use the proposed VFD filter for wideband sparse arrays during aperture correction.

6. Conclusions

In order to reduce the system complexity of the VFD FIR filter, a new MRMM sparse algorithm under a great–minimal criterion is proposed. The error function is constrained by introducing L2 and L1 regularizations. To solve this nonconvex model, an improved S–ADMM algorithm is proposed, which uses sequential alternation to iteratively update multiple soft–threshold problems. Experimental results show that the proposed MRMM algorithm completely satisfies the desired error value and sparsity. The use of multipliers and adders is reduced by 71.56% and 71.73%, respectively, leading to a significant decrease in system complexity. Simulation results confirm that the MRMM–VFD effectively corrects AFT in a wideband sparse array.

Author Contributions

Conceptualization, J.G. and M.X.; methodology, J.G., M.X. and M.S.; software, M.S. and W.Z.; validation, J.G., M.X., M.S. and W.Z.; formal analysis, J.G. and M.X.; investigation, J.G. and M.X.; resources, J.G., M.X. and M.S.; data curation, J.G. and M.X.; writing—original draft preparation, J.G.; writing—review and editing, M.X., M.S. and W.Z.; visualization, J.G.; supervision, J.G. and M.X.; project administration, J.G. and M.X.; funding acquisition, J.G., M.X., M.S. and W.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Natural Science Foundation of China [No.62271190], and Science and Technology on Electronic Information Control Laboratory [No.20220003].

Data Availability Statement

The data used to support the findings of this study are available from the corresponding author upon request.

Acknowledgments

We sincerely thank all the reviewers for their patient reviews and valuable suggestions on this article.

Conflicts of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

References

  1. Edussooriya, C.U.S.; Bruton, L.T.; Naeini, M.A.; Agathoklis, P. Using 1-D Variable Fractional-Delay Filters to Reduce the Computational Complexity of 3-D Broadband Multibeam Beamformers. IEEE Trans. Circuits Syst. II Express Briefs 2014, 61, 278–283. [Google Scholar] [CrossRef]
  2. Canese, L.; Cardarilli, G.C.; Di Nunzio, L.; Fazzolari, R.; Giardino, D.; Re, M.; Spanó, S. Efficient Digital Implementation of a Multirate-Based Variable Fractional Delay Filter for Wideband Beamforming. IEEE Trans. Circuits Syst. II Express Briefs 2023, 70, 2231–2235. [Google Scholar] [CrossRef]
  3. Kokila, R.; Chithra, K.; Dhilsha, R. Wideband Beamforming Using Modified Farrow Structure FIR Filtering Method for Sonar Applications. In Proceedings of the 2019 International Symposium on Ocean Technology (SYMPOL), Ernakulam, India, 11–13 December 2019. [Google Scholar]
  4. Wang, J.; Cai, D.; Lei, P. Aperture effect influence and analysis of stepped frequency wideband phased array radar. In Proceedings of the IEEE CIE International Conference on Radar, Chengdu, China, 24–27 October 2011. [Google Scholar]
  5. Selva, J. An Efficient Structure for the Design of Variable Fractional Delay Filters Based on the Windowing Method. IEEE Trans. Signal Process. 2008, 56, 3370–3775. [Google Scholar] [CrossRef]
  6. Deng, T.-B. Decoupling Minimax Design of Low-Complexity Variable Fractional-Delay FIR Digital Filters. IEEE Trans. Circuits Syst. I Regul. Pap. 2011, 58, 2398–2408. [Google Scholar] [CrossRef]
  7. Deng, T.-B. Minimax Design of Low-Complexity Even-Order Variable Fractional-Delay Filters Using Second-Order Cone Programming. IEEE Trans. Circuits Syst. II Express Briefs 2011, 58, 692–696. [Google Scholar] [CrossRef]
  8. Huang, Y.-D.; Pei, S.-C.; Shyu, J. WLS Design of Variable Fractional-Delay FIR Filters Using Coefficient Relationship. IEEE Trans. Circuits Syst. II Express Briefs 2009, 56, 220–224. [Google Scholar] [CrossRef]
  9. Deng, T.-B.; Qin, W.-B.; Qin, W. Coefficient relation-based minimax design and low-complexity structure of variable fractional-delay digital filters. Signal Process. 2012, 93, 923–932. [Google Scholar] [CrossRef]
  10. Farrow, C.W. A continuously variable digital delay element. In Proceedings of the IEEE International Symposium on Circuits and Systems, Espoo, Finland, 7–9 June 1988. [Google Scholar]
  11. Deng, T.-B. Symmetric Structures for Odd-Order Maximally Flat and Weighted-Least-Squares Variable Fractional-Delay Filters. IEEE Trans. Circuits Syst. I Regul. Pap. 2007, 54, 2718–2732. [Google Scholar] [CrossRef]
  12. Jiang, A.M.; Kwan, H.K.; Zhu, Y.P.; Liu, X.F.; Xu, N.; Tang, Y.B. Design of sparse FIR filters with joint optimization of sparsity and filter order. IEEE Trans. Circuits Syst. I Regul. Pap. 2015, 62, 195–204. [Google Scholar] [CrossRef]
  13. Chen, W.Q.; Huang, M.; Lou, X. Design of Sparse FIR Filters With Reduced Effective Length. IEEE Trans. Circuits Syst. I Regul. Pap. 2019, 66, 1496–1506. [Google Scholar] [CrossRef]
  14. Lu, W.-S.; Hinamoto, T. Variable fractional delay FIR filters with sparse coefficients. In Proceedings of the 2012 IEEE International Symposium on Circuits and Systems (ISCAS), Seoul, Republic of Korea, 20–23 May 2012. [Google Scholar]
  15. Zhou, W.J.; Shen, M.W.; Xu, M.; Han, G.; Zhang, Y. Sparsity-optimised farrow structure variable fractional delay filter for wideband array. IET Signal Process. 2023, 17, e12228. [Google Scholar] [CrossRef]
  16. Shentu, X.; Lai, X.; Wang, T.; Cao, J. Efficient ADMM-Based Algorithm for Regularized Minimax Approximation. IEEE Signal Process. Lett. 2023, 30, 210–214. [Google Scholar] [CrossRef]
  17. Chiu, W.-Y.; Chen, B.-S.; Yang, C.-Y. Robust Relative Location Estimation in Wireless Sensor Networks with Inexact Position Problems. IEEE Trans. Mob. Comput. 2012, 11, 935–946. [Google Scholar] [CrossRef]
  18. Nakamoto, M.; Aikawa, N. Minimax Design of Sparse IIR Filters Using Sparse Linear Programming. IEICE Trans. Fundam. Electron. Commun. Comput. Sci. 2012, 104, 1006–1018. [Google Scholar] [CrossRef]
  19. Nordebo, S.; Claesson, I. Minimum norm design of two-dimensional weighted Chebyshev FIR filters. IEEE Trans. Circuits Syst. II Analog Digit. Signal Process. 1997, 44, 251–253. [Google Scholar] [CrossRef]
  20. Lu, W.-S.; Hinamoto, T. Digital filters with sparse coefficients. In Proceedings of the 2010 IEEE International Symposium on Circuits and Systems, Paris, France, 30 May–2 June 2010. [Google Scholar]
  21. Baran, T.; Wei, D.; Oppenheim, A.V. Linear Programming Algorithms for Sparse Filter Design. IEEE Trans. Signal Process. 2010, 58, 1605–1617. [Google Scholar] [CrossRef]
  22. Matsuoka, R.; Kyochi, S.; Ono, S.; Okuda, M. Joint Sparsity and Order Optimization Based on ADMM With Non-Uniform Group Hard Thresholding. IEEE Trans. Circuits Syst. I Regul. Pap. 2018, 65, 1602–1613. [Google Scholar] [CrossRef]
  23. Yanamshetti, R.; Bharatkar, S.S.; Chatterjee, D.; Ganguli, A.K. A hybrid fuzzy based loss minimization technique for fast efficiency optimization for variable speed induction machine. In Proceedings of the 2009 6th International Conference on Electrical Engineering/Electronics, Computer, Telecommunications and Information Technology, Chonburi, Thailand, 6–9 May 2009. [Google Scholar]
  24. Hawes, M.B.; Liu, W. Sparse array design for wideband beamforming with reduced complexity in tapped delay-lines. IEEE/ACM Trans. Audio Speech Lang. Process. 2014, 22, 1236–1247. [Google Scholar] [CrossRef]
  25. Zhao, R.; Hong, X. Matrix-Based Algorithms for the Optimal Design of Variable Fractional Delay FIR Filters. IEEE Trans. Circuits Syst. I Regul. Pap. 2019, 66, 4740–4752. [Google Scholar] [CrossRef]
Figure 1. The Farrow structure the VFD FIR filter.
Figure 1. The Farrow structure the VFD FIR filter.
Sensors 24 04327 g001
Figure 2. Aperture crossing correction of wideband sparse array.
Figure 2. Aperture crossing correction of wideband sparse array.
Sensors 24 04327 g002
Figure 3. Convergence of the S−ADMM algorithm: (a) shows the convergence of r 2 and s 2 ; and (b) shows the convergence of the cost function given by Equation (16).
Figure 3. Convergence of the S−ADMM algorithm: (a) shows the convergence of r 2 and s 2 ; and (b) shows the convergence of the cost function given by Equation (16).
Sensors 24 04327 g003
Figure 4. Amplitude–frequency response error plots for different algorithms. (a) is the MM algorithm; (b) is the S–WLS algorithm; and (c) is the MRMM algorithm.
Figure 4. Amplitude–frequency response error plots for different algorithms. (a) is the MM algorithm; (b) is the S–WLS algorithm; and (c) is the MRMM algorithm.
Sensors 24 04327 g004aSensors 24 04327 g004b
Figure 5. Amplitude–frequency response plots of different algorithms. (a) is the MM algorithm; (b) is the S–WLS algorithm; and (c) is the MRMM algorithm.
Figure 5. Amplitude–frequency response plots of different algorithms. (a) is the MM algorithm; (b) is the S–WLS algorithm; and (c) is the MRMM algorithm.
Sensors 24 04327 g005
Figure 6. Phase–frequency plot for different algorithms at a fractional delay p = 0.3 .
Figure 6. Phase–frequency plot for different algorithms at a fractional delay p = 0.3 .
Sensors 24 04327 g006
Figure 7. Sparse array location distribution.
Figure 7. Sparse array location distribution.
Sensors 24 04327 g007
Figure 8. Wideband sparse array antenna orientation diagram. (a) is the result of aperture crossing; and (b) is the result of aperture correction.
Figure 8. Wideband sparse array antenna orientation diagram. (a) is the result of aperture crossing; and (b) is the result of aperture correction.
Sensors 24 04327 g008
Figure 9. Wideband sparse array antenna orientation diagram.
Figure 9. Wideband sparse array antenna orientation diagram.
Sensors 24 04327 g009
Table 1. Comparison of error performance of different algorithms.
Table 1. Comparison of error performance of different algorithms.
M = 5M = 7M = 9M = 11
ε 1 ε 2 ε 1 ε 2 ε 1 ε 2 ε 1 ε 2
MM−19.830.0067−20.010.0061−18.690.0085−23.240.0047
S–WLS−17.740.01−20.780.0054−18.710.0072−20.440.0061
MRMM−18.370.0097−19.450.0080−19.810.0066−21.980.0058
Table 2. Comparison of sparsity and operation time of different algorithms.
Table 2. Comparison of sparsity and operation time of different algorithms.
M = 5M = 7M = 9M = 11
η Tc η Tc η Tc η Tc
MRMM66.67%
(248)
0.16 s75.40%
(374)
0.29 s82.64%
(512)
0.61 s85.75%
(638)
0.73 s
LCMM37.97%
(139)
0.62 s21.77%
(108)
2.05 s31.93%
(198)
2.40 s32.52%
(242)
3.67 s
S–WLS69.89%
(260)
0.14 s75.80%
(376)
0.26 s81.29%
(504)
0.77 s85.21%
(634)
0.92 s
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Gu, J.; Xu, M.; Zhou, W.; Shen, M. Efficient Aperture Fill Time Correction for Wideband Sparse Array Using Improved Variable Fractional Delay Filters. Sensors 2024, 24, 4327. https://doi.org/10.3390/s24134327

AMA Style

Gu J, Xu M, Zhou W, Shen M. Efficient Aperture Fill Time Correction for Wideband Sparse Array Using Improved Variable Fractional Delay Filters. Sensors. 2024; 24(13):4327. https://doi.org/10.3390/s24134327

Chicago/Turabian Style

Gu, Jie, Min Xu, Wenjing Zhou, and Mingwei Shen. 2024. "Efficient Aperture Fill Time Correction for Wideband Sparse Array Using Improved Variable Fractional Delay Filters" Sensors 24, no. 13: 4327. https://doi.org/10.3390/s24134327

APA Style

Gu, J., Xu, M., Zhou, W., & Shen, M. (2024). Efficient Aperture Fill Time Correction for Wideband Sparse Array Using Improved Variable Fractional Delay Filters. Sensors, 24(13), 4327. https://doi.org/10.3390/s24134327

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