Next Article in Journal
Resource Allocation of UAV-Assisted IoT Node Secure Communication System
Previous Article in Journal
Computer Vision and Image Processing in Structural Health Monitoring: Overview of Recent Applications
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Nonlinear Optimization Design Algorithm for Nearly Linear-Phase 2D IIR Digital Filters

1
Department of Electrical and Computer Engineering, University of Victoria, Victoria, BC V8W 3P6, Canada
2
Department of Computer Science, Camosun College, Victoria, BC V8P 5J2, Canada
*
Author to whom correspondence should be addressed.
Signals 2023, 4(3), 575-590; https://doi.org/10.3390/signals4030030
Submission received: 7 May 2023 / Revised: 27 June 2023 / Accepted: 20 July 2023 / Published: 2 August 2023

Abstract

:
In this paper, a new optimization method for the design of nearly linear-phase two-dimensional infinite impulse (2D IIR) digital filters with a separable denominator is proposed. A design framework for 2D IIR digital filters is formulated as a nonlinear constrained optimization problem where the group delay deviation in the passband is minimized under prescribed soft magnitude constraints and hard stability requirements. To achieve this goal, sub-level sets of the group delay deviations are utilized to generate a sequence of filters, from which the one with the best performance is selected. The quality of the obtained filter is evaluated using three quality factors, namely, the passband magnitude quality factor  Q h  and the group delay deviation quality factor  Q τ , while the third one is a new quality factor  Q s  that assesses the performance in the stopband relative to the minimum filter gain in the passband. The proposed framework is implemented using the interior-point (IP) method in a MATLAB environment, and the experimental results show that filters designed using the proposed method have good magnitude response and low group delay deviation. The performance of the resulting filters is compared with the results of other methods.

1. Introduction

Two-dimensional (2D) digital filters have found many signal processing applications, such as image processing, video signal filtering, beamforming systems, and remote sensing. Digital filters are broadly classified into two main categories, namely, finite impulse response (FIR) filters and infinite impulse response (IIR) filters [1]. An FIR filter is one whose impulse response is of finite duration. The output of such a filter is calculated solely from the current and previous input values. This type of filter is hence said to be non-recursive. On the other hand, an IIR filter is one whose impulse response (theoretically) continues forever in time. They are also termed as recursive filters. In addition to the current and previous input values, the current output of such a filter depends on the previous output values.
FIR filters are generally easier to implement, as they are non-recursive, always stable (by definition), and can be designed to have a linear phase characteristic. They are, however, not well-suited for very sharp (narrow transition band) frequency responses. IIR filters, on the other hand, are much more difficult to design due to their non-linear phase and their stability constraints. However, they can effectively accommodate very narrow transition band frequency responses, which makes them suitable for a broad range of applications. Moreover, selective IIR filters are much more efficient than FIR filters, as the latter require higher orders (i.e., more multiplications) [1,2,3,4,5,6,7,8,9,10]. The efficiency of IIR filters stems from their ability to satisfy the required filter specifications with a much lower filter order, thereby reducing the computational requirements. In addition, IIR filters can have a much smaller group delay compared with FIR filters [1,11,12,13]. Due to these advantages, IIR filters find extensive applications in the domains of de-noising of digital images, biomedical imaging and digital mammography, X-ray image enhancement, and seismic data processing.
Since many signal processing applications are very sensitive to phase non-linearity, the design of IIR filters needs to take into consideration not only the amplitude requirements but also the phase constraints. Therefore, for most of these applications, the designed IIR filters should have a nearly linear-phase response. Moreover, hard constraints on the location of the poles of IIR filters must be enforced to ensure the stability of the obtained filters [14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29]. These and other design requirements usually lead to highly nonlinear constrained optimization problems that require highly sophisticated optimization methods [1,12,29,30,31,32]. Research on 2D IIR filters has produced many constrained and non-constrained optimization problems, depending on the researchers’ design objectives. Moreover, each researcher has designed specialized algorithms and techniques to handle their optimization problems. Most of these methods can be classified into linear and nonlinear techniques. In linear design techniques, the nonlinear optimization problem is approximated by a sequence of linear optimization problems [1,11], which are then solved using linear programming. Starting from an initial filter that satisfies most of the design specifications, the algorithm linearizes the nonlinear optimization around this initial filter to produce a linear optimization problem. Solving this problem leads to an improved filter, around which the linearization is carried out. This iterative process is repeated until an acceptable filter is achieved. Although the linear methods can be very fast, and the designer might be able to find an initial filter satisfying most of the design specifications, there is no guarantee that the will converge or produce accurate solutions. As such, nonlinear techniques are preferred. Some researchers have attempted to develop design methods based on modern optimization algorithms.
In [29], the Sequential Partial Optimization algorithm was proposed to convert the minmax design problem into a sequence of smaller sub-problems, each updating only a single second-order denominator factor.
In [33], a modern optimization methodology known as semi-definite programming (SDP) is utilized. This SDP method includes minmax and weighted least-squares designs of FIR filters as well as minmax formulations of stable separable-denominator IIR filters. In [33] describes how SDP has the potential to serve as a unified optimization engine for designing a wide range of 2D digital filters.
The Quasi-Newton method has been applied to design 2D IIR digital filters [7], but it may result in unstable filters. Neural networks were used in [34] to design 2D IIR filters, which can guarantee the stability of the designed filters with low computational time. Some researchers have employed the genetic algorithm (GA) to design 2D IIR digital filters [32,35,36,37,38]. Others have proposed hybrid methods that combine different evolutionary algorithms, such as biogeography-based optimization and particle swarm optimization in [39], to tackle the nonlinear 2D IIR optimization design.
In this paper, a robust nonlinear optimization method for the design of nearly linear-phase 2D IIR digital filters with separable-denominators is proposed. This design problem is formulated as a nonlinear constrained optimization problem where the group delay deviation is minimized under prescribed soft magnitude constraints and hard stability requirements. We also rely on sub-level sets of the group delay deviations to obtain a sequence of filters, from which we choose the most optimal one. The performance of the obtained 2D IIR filters is evaluated using three quality factors. These first two quality factors are the usual passband magnitude quality factor  Q h  and the group delay deviation quality factor  Q τ , while the third one is a new quality factor  Q s  that assesses the performance of the filter within the stopband relative to the minimum filter gain in the passbands, and that is because amplitude response anomalies in transition regions are undesirable [40]. The obtained results show that these quality factors are all necessary to judge the optimality of the filter. The proposed framework is implemented using the interior-point (IP) method in a MATLAB environment, and the results are very promising.
The remainder of the paper is organized as follows. The constrained optimization problem for 2D IIR digital filters is formulated in Section 2. Section 3 is for the stability requirements for the 2D IIR digital filters with a separable denominator. Section 4 presents the constrained optimization design problem for a 2D IIR digital filters. Section 5 is dedicated to the quality of the design, and Section 6 to the experimental results. Section 7 concludes with a brief summary and future work directions.

2. Formulation of the Design Problem

Much attention has been paid to the optimal design of 2D filters with separable denominators [18,29,30,33,41,42,43]. The stability of such 2D filters can be easily guaranteed by ensuring the roots of both denominators are within the unit bi-disc independently. As such, the internal stability triangles for 1D second-order factors are used as stability constraints [44]. In addition, separable denominators reduce the complexity of the 2D filter and the induced optimization problem and also result in fewer coefficients [29,45]. Although the phase of the optimized 2D IIR filter is only nearly linear, one can see that the separable denominators 2D IIR filter offers good selectivity and computation efficiency and reduced system delay without compromising the required amplitude specifications when compared with the corresponding other 2D filters [45]. As such, we will assume in this paper that the filters of interest have separable denominators. The desired 2D IIR filter can then be represented by the following transfer function [46,47,48]:
H ( a , b , z 1 , z 2 ) = N ( z 1 , z 2 , b ) D ( z 1 , z 2 , a ) = H 0 N ( z 1 , z 2 , b ) D 1 ( z 1 , a 1 ) D 2 ( z 2 , a 2 ) ,
where b is a vector of the form  [ b 00 b 01 b N 1 N 2 ] , and  N ( z 1 , z 2 , b )  is a multivariate polynomial of the form:
N ( z 1 , z 2 , b ) = j = 0 N 2 i = 0 N 1 b i j z 1 i z 2 j ,
D 1 ( z 1 , a 1 )  and  D 2 ( z 2 , a 2 )  are polynomials in  z 1  and  z 2 , respectively, of the form
D 1 ( z 1 , a 1 ) = n = 1 J 1 ( 1 + a 1 n , 1 z 1 1 + a 1 n , 2 z 1 2 ) ,
D 2 ( z 2 , a 2 ) = n = 1 J 2 ( 1 + a 2 n , 1 z 2 1 + a 2 n , 2 z 2 2 ) ,
where  J 1 = N 1 / 2  and  J 2 = N 2 / 2  and  N 1  and  N 2  are even.
The vectors  a 1  and  a 2  are given by:
a 1 = [ a 11 , 1 a 11 , 2 a 12 , 1 a 12 , 2 a 1 J 1 , 1 a 1 J 1 , 2 ] ,   and
a 2 = [ a 21 , 1 a 21 , 2 a 22 , 1 a 22 , 2 a 2 J 2 , 1 a 2 J 2 , 2 ] .
We then construct a vector c by concatenating the vectors b a 1 , and  a 2  as follows [44]:
c = [ b a 1 a 2 H 0 ] T .
Using the same notation symbol, the transfer function can be written as follows:
H ( c , z 1 , z 2 ) = N ( z 1 , z 2 , b ) D ( z 1 , z 2 , a ) .
Now, the main design problem is to determine the coefficients in the numerator and denominator of Equation (7) in such a way that the obtained transfer function  H ( z 1 , z 2 )  satisfies all the desired specifications as explained in the following sections.

2.1. Passband Group Delay Deviations

Let  ϕ ( c , w 1 , w 2 )  and  ψ ( c , w 1 , w 2 )  be the phase angles of  N ( e j w 1 , e j w 2 , b )  and  D ( e j w 1 , e j w 2 , a ) , respectively. The group delay of the transfer function  H ( c , z 1 , z 2 )  is given as a partial derivative with respect to  w 1  and  w 2  as follows:
τ 1 ( c , w 1 , w 2 ) = ϕ ( c , w 1 , w 2 ) w 1 ψ ( c , w 1 , w 2 ) w 1 ,   and τ 2 ( c , w 1 , w 2 ) = ϕ ( c , w 1 , w 2 ) w 2 ψ ( c , w 1 , w 2 ) w 2 .
Letting  τ 10  and  τ 20  be the intended average passband group delays of the initial filter in each dimension over a domain  Ω p b , the deviation can be computed as follows:
e g i ( x , w 1 , w 2 ) = τ i ( x , w 1 , w 2 ) τ i 0 , for   i = 1 , 2 , ( w 1 , w 2 ) Ω p b ,
where  x = [ c T τ 10 τ 20 ] T , and  Ω p b  is the passband region [44].
At the start of the optimization,  τ 10  and  τ 20  can be assigned values equal to the average group delays of the initial filter. As the optimization progresses, their values can be adjusted along with the coefficients of the filter so as to minimize the group delay deviation. The end result will be that  τ 10  and  τ 20  will track the average value of the group delays in each dimension. As such, we have added them to the vector of variables to be optimized.

2.2. Passband Amplitude Error

For the passband error, if we consider  H d ( w 1 , w 2 )  to be the desired frequency response of the filter and that  c k  is the value of vector c at the start of the  k t h  iteration, a passband error function at frequencies  w 1  and  w 2  can be defined as
e h ( x , w 1 , w 2 ) | H ( c k , w 1 , w 2 ) | | H d ( w 1 , w 2 ) | .
Assuming that the desired amplitude response is unity in the passband, therefore, the passband error function becomes:
e h ( x , w 1 , w 2 ) | H ( c k , w 1 , w 2 ) | 1 , ( w 1 , w 2 ) Ω p b .
The expression above can be used as is in the case of the nonlinear approach.
For the transition and stopband regions, the magnitude response is subject to prescribed levels of transition-band and stopband gains, as explained in Section 4.

3. Filter Stability

One of the major problems underlying the design task is satisfying the stability criterion for the 2D IIR filter transfer function. For the stability constraints of the 2D IIR filter, since the filter transfer function is separable in the denominator, the stability analysis is similar to the 1D case used in [11,44,49], where all the poles of  D 1 ( z 1 , a 1 )  and  D 2 ( z 2 , a 2 )  have to be located inside the unit circle U of the z-plane for both  z 1  and  z 2 .

4. The Constrained Optimization Problem

To solve the design problem, the group delay deviation is minimized under hard and soft constraints. Hard constraints are the requirements that the filter must satisfy. These requirements must include at least the stability constraints of the filter. On the other hand, soft constraints are the requirements that the filter should try to satisfy. These can, for example, include the passband ripple constraint, transition-band gain. and stopband attenuation.
The general constrained optimization problem at hand can be stated as follows:
min x max w 1 , w 2 Ω p b { | e g 1 ( x , w 1 , w 2 ) | , | e g 2 ( x , w 1 , w 2 ) | }
subject to : passband error Γ p b , transition - band gain Γ t b , stopband gain Γ s b , stability margin ε s ,
where  Γ p b , Γ t b , and  Γ s b  are the maximum prescribed levels for the passband error, transition-band gain, and stopband gain, respectively. The stability margin,  ε s , is defined as  ε s = 1 r m a x ( p )  where  r m a x ( p )  is the maximum pole radius allowed [11]. The above general optimization problem for the nearly linear-phase 2D IIR digital filters can now be expressed as follows:
min x max w 1 , w 2 Ω p b { | τ 1 ( x , w 1 , w 2 ) τ 10 | , | τ 2 ( x , w 1 , w 2 ) τ 20 | }
subject to : | e h ( x , w 1 , w 2 ) | Γ p b , ( w 1 , w 2 ) Ω p b , | H ( x , w 1 , w 2 ) | Γ t b , ( w 1 , w 2 ) Ω t b , | H ( x , w 1 , w 2 ) | Γ s b , ( w 1 , w 2 ) Ω s b , | τ 1 τ 10 | Γ g 1 , ( w 1 , w 2 ) Ω p b , | τ 2 τ 20 | Γ g 2 , ( w 1 , w 2 ) Ω p b , B 1 x 1 γ , B 2 x 1 γ ,
where  Γ g 1  and  Γ g 2  are the maximum errors between the group delays of the actual and the initial filters.  γ = 1 ( 1 ε s ) 2 B 1  and  B 2  are the stability matrices for  D 1 ( z 1 , a 1 )  and  D 2 ( z 2 , a 2 ) , respectively [11,44]. The 2D IIR filter design problem can be readily solved using interior-point method. This is because the interior-point method is efficient from a computational point of view and also highly efficient in practice [50,51]. To make our implementation easier, we used  Γ g = Γ g 1 = Γ g 2 . By varying  Γ g , and therefore  Γ g 1  and  Γ g 2 , we obtain different filters from which we choose the best filter by looking at the quality factors discussed in the next section. The proposed framework is implemented using the interior-point (IP) method in a MATLAB environment and the flowchart of the optimization process is shown in Figure 1.

5. Quality of the Design

Design efficiency and the performance of the proposed method are evaluated in terms of the quality of the group delay characteristic  Q τ , the quality of the magnitude response in the passband  Q h , and a new quality factor  Q s  for the magnitude attenuation in the stopband. We introduced this new factor,  Q s , to measure the quality of the magnitude in the stopband. The quality of the group delay characteristic of the 2D filter will be measured by using the normalized maximum variation of the filter group delay over the passband as a percentage as follows:
Q τ = max τ 1 m a x τ 1 m i n τ 1 m a x + τ 1 m i n , τ 2 m a x τ 2 m i n τ 2 m a x + τ 2 m i n × 100 % ,
To measure the quality of the magnitude in the passband, we use:
Q h = H m a x H m i n H m a x + H m i n × 100 % ,
where  H m a x  and  H m i n  are the maximum and minimum magnitudes in the passband, respectively. We observe that when we design the filter, we improve the passband ripple and the group delays, but sometimes the stopband attenuation is not good enough. Based on that observation, we introduced a new quality factor for the stopband attenuation, which is measured in terms of the minimum variation of the passband ripple and the maximum variation of the stopband attenuation. To measure the quality of the magnitude response in the stopband, we use:
Q s = H m a x s H m i n H m a x s ,
where  H m a x s  is the maximum magnitude in the stopband region, and  H m i n  is the minimum magnitude in the passband. This quality factor,  Q s , is valid for  H m i n > H m a x s , i.e., negative values are excluded. Note that these quality factors are inverse measures, and therefore a good filter is one where these quality factors are small.

6. Experimental Results

This section is dedicated to implementing Equations (12) and (13) to optimize the design of 2D IIR filters with nearly linear phase. Our proposed approach will be evaluated against other methods described in the literature. Specifically, we will compare the designs obtained using the proposed method to the ones presented in [8]. Some of these examples are discussed in [52,53] as well. The desired frequency response will be defined in each example, and the group delays  τ 1  and  τ 2  along the  ω 1  and  ω 2  directions are the same, i.e.,  τ 1 = τ 2 = τ . In each of the designed examples, the starting values for  τ 10  and  τ 20  are set as the mean values of the passband group delays of the initial filter, and the stability margin  ε s  is set to 0.02, unless otherwise mentioned in the design example. The MATLAB environment is used to solve the constrained optimization problems. The MATLAB Symbolic Math Toolbox and the interior-point method (IP) were used to implement the constrained optimization problem.

6.1. 2D Highpass Filter

In this example, it is required to design a 2D highpass filter with the ideal filter specifications shown in Figure 2, where  θ i = w i T . In addition to the given design specifications, we are including the transition band in our design, which aligns with the filter obtained in [8]. To assess the performance of our filter and the one proposed by [8], we have calculated the maximum transfer function magnitude in the transition band  H max t  for both filters, as presented in Table 1.
In this example, the max. iterations = 140,  ϵ = 1 × 10 7 ε s = 0.1 Γ p b = 0.01 Γ t b = 1.5 Γ s b = 0.1 H 0 = 1 . We implement the proposed algorithm for various values of  Γ g  starting from  0.1  to  4.0 . Using the quality factors defined in Section 5, the filter with  Γ g = 1.65  is the most optimal filter among the filters obtained. Figure 3 shows the magnitude response of the designed filter using the proposed method, and the filter presented in [8]. Figure 4 and Figure 5 show the group delays  τ 1  and  τ 2  of the obtained 2D IIR highpass filter and the filter in [8], respectively. Table 1 compares the results of the highpass filter designed by the proposed algorithm with the highpass filter presented in [8].
The quality measure of the magnitude response in the passband, defined by Equation (15), of the obtained 2D IIR highpass filter is  Q h = 50.76 , where  H m i n = 0.437  and  H m a x = 1.338 . In the same manner, the quality measure of the magnitude response in the passband of the filter presented in [8] is  Q h = 60.54 , where  H m i n = 0.262  and  H m a x = 1.066 . The quality of the magnitude response in the stopband defined by Equation (16) is  Q s = 8.104  for the filter obtained by the proposed method and for the filter presented in [8],  Q s  is negative, indicating too low stopband attenuation compared with the minimum gain in the passband. The quality of the group delay characteristic, defined by Equation (14), of the highpass filter designed by the proposed method is better than the one designed in [8]. It can be seen that  Q τ = 64.09 , where  τ m i n = 0.962  and  τ m a x = 4.397  for the filter obtained, and  Q τ = 99.84 , whereas  τ m i n = 0.009  and  τ m a x = 11.614  for the filter designed in [8].

6.2. 2D Lowpass Filter

The proposed method is used to design a 2D IIR lowpass filter with the ideal design specifications shown in Figure 6. The filter was presented in [8], and its order was selected as (4,4). The transition band constraint is included in this example, and  H m a x t  is computed for both filters, where the results are shown in Table 2.
This 2D lowpass filter is designed by the proposed method using the following values: Max. iterations = 140,  ϵ = 1 × 10 7 ε s = 0.1 Γ p b = 0.0001 Γ t b = 1 Γ s b = 0.01 H 0 = 1 τ 1 = τ 2 = 2 . The algorithm is implemented for different values of  Γ g , starting from 0.1 to 4. The filter with  Γ g = 2.35  has the best quality in this example.
Figure 7 shows the impulse and magnitude responses of the filter obtained by using the proposed method and the filter presented in [8]. Figure 8 and Figure 9 show the group delays,  τ 1  and  τ 2 , of the obtained 2D lowpass filter for  Γ g = 2.35 . The group delays,  τ 1  and  τ 2 , of the 2D lowpass filter presented in [8] are shown in Figure 10 and Figure 11, respectively. Table 2 compares the results of the filter designed by the proposed method and the filter obtained in [8]. The quality measure of the magnitude response in the passband, defined by Equation (15) of the 2D lowpass filter obtained by the proposed method is  Q h = 14.25 , where  H m i n = 0.843  and  H m a x = 1.120 , and for the filter presented in [8] is  Q h = 25.26 , where  H m i n = 0.642  and  H m a x = 1.069 . The quality of the magnitude response in the stopband defined by Equation (16) is  Q s = 0.916  for the filter designed by the proposed method, and  Q s = 2.275  for the filter presented in [8]. The quality of the group delay characteristic, defined by Equation (14), for the 2D lowpass filter obtained by the proposed method is  Q τ = 47.650 , where  τ m i n = 2.217  and  τ m a x = 6.864 . The same quality measure for the 2D lowpass filter presented in [8] is found to be  Q τ = 93.774 , where  τ m i n = 1.738  and  τ m a x = 54.091 . We can see that the filter designed with the proposed method is better than the one designed in [8], especially with respect to the passband phase linearity.

6.3. 2D Bandpass Filter

In this example, a 2D IIR bandpass filter with the design specifications shown in Figure 12 is presented. It is also required that the order of this filter remain low, specifically (4,4) as in [8]. In addition to these design specifications and to assess the performance of our filter and the one proposed in [8], we have included the transition band constraint in our design. As a result, we calculated the maximum transfer function magnitude  H max t  for both filters, as presented in Table 3. The algorithm is implemented with various values of  Γ g , starting from 0.1 to 4, and the resulting filter with  Γ g = 1.6  has the best quality.
The magnitude response of the obtained 2D IIR bandpass filter and the filter designed in [8] are shown in Figure 13. The group delays  τ 1  and  τ 2  of the designed 2D bandpass filter are illustrated in Figure 14. Figure 15 shows the group delays of the same filter in the first quadrant. The group delays of the 2D bandpass filter presented in [8] are illustrated in Figure 16. The first quadrant plot of the group delays  τ 1  and  τ 2  of the same filter are shown in Figure 17. Table 3 compares the results of the bandpass filter designed by the proposed algorithm and the bandpass filter obtained in [8]. The quality measure of the magnitude response in the passband, defined by Equation (15), of the 2D IIR bandpass filter obtained by the proposed method is  Q h = 19.09 , where  H m i n = 0.809  and  H m a x = 1.191 . In the same manner, the quality measure of the magnitude in the passband region of the 2D IIR bandpass filter presented in [8] is  Q h = 29.667 , where  H m i n = 0.633  and  H m a x = 1.167 . The quality of the magnitude response in the stopband defined by Equation (16) is  Q s = 0.668  for the filter obtained by the proposed method, and  Q s = 0.734  for the filter presented in [8]. The quality of the group delay characteristic, defined by Equation (14), of the 2D IIR bandpass filter obtained by the proposed method is  Q τ = 47.150 , where  τ m i n = 1.822  and  τ m a x = 5.073 . The same measure is applied to the 2D IIR bandpass filter designed in [8], which gives  Q τ = 57.315 , where  τ m i n = 1.828  and  τ m a x = 6.737 .
In this example, the filter obtained with the proposed algorithm has a better passband quality factor,  Q h , stopband quality factor,  Q s , and group delay quality factor,  Q τ , than the filter presented in [8].

7. Conclusions

A new framework for a nonlinear optimization design method for nearly linear-phase 2D IIR digital filters with separable denominators is presented. The method is used to solve the design problem, expressed as a nonlinear constrained optimization problem. In our implementation, we distinguish two kinds of constraints: soft and hard constraints. Meeting the soft constraints is not compulsory and, therefore, not decisive in rejecting the filter. The hard constraints, however, must be satisfied to accept the filter. The designed filters are evaluated with three quality factors, namely passband amplitude and group delay quality factors, and a new quality factor for the stopband magnitude response. compared with other 2D IIR filters designed with other methods, the results show that the framework produces very good results for different filter specifications.
Evaluating the resulted filters, this method improves the linearity of the phase in the passband with magnitude improvements, generally without significant deterioration in the magnitude response. The method is also flexible enough to be used after any other 2D IIR design algorithm by taking the resulting filter as the initial filter and running this method to obtain possible further improvements.
In future work, we intend to improve our nonlinear optimization method by directly optimizing the quality factors instead of the group delay deviations only.

Author Contributions

Conceptualization, A.O. and P.A.; methodology, A.O., P.A. and D.S.; software, A.O., P.A. and B.M.; validation, A.O., P.A. and D.S.; formal analysis, A.O., P.A. and D.S.; investigation, A.O., P.A. and D.S.; resources, A.O., P.A. and B.M.; data curation, A.O. and P.A.; writing—original draft preparation, A.O.; writing—review and editing, A.O., P.A. and D.S.; visualization, A.O., P.A. and D.S.; supervision, P.A. and D.S.; project administration, A.O., P.A. and D.S. All authors have read and agreed to the published version of the manuscript.

Funding

Funding provided by NSERC.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
2DTwo Dimensional
1DOne Dimensional
IIRInfinite Impulse Response
FIRFinite Impulse Response
IPInterior Point
SDPSemi Definite Programming
GAGenetic Algorithm
MinMinimum
MaxMaximum
UUnit Circle

References

  1. Antoniou, A. Digital Signal Processing: Signals, Systems, and Filters; McGraw-Hill: New York, NY, USA, 2005. [Google Scholar]
  2. Matei, R. Analytic design of directional and square-shaped 2D IIR filters based on digital prototypes. Multidimens. Syst. Signal Process. 2019, 30, 2021–2043. [Google Scholar] [CrossRef]
  3. Chan, T.S.; Kumar, A. Reliable ear identification using 2-D quadrature filters. Pattern Recognit. Lett. 2012, 33, 1870–1881. [Google Scholar] [CrossRef]
  4. Wu, C.W. Bit-level pipelined 2-D digital filters for real-time image processing. IEEE Trans. Circuits Syst. Video Technol. 1991, 1, 22–34. [Google Scholar] [CrossRef]
  5. Wijesekara, R.T.; Edussooriya, C.U.; Bruton, L.T.; Agathoklis, P. A low-complexity 2-D spatially-interpolated FIR trapezoidal filter for enhancing broadband plane waves. In Proceedings of the 2017 10th International Workshop on Multidimensional (nD) Systems (nDS), IEEE, Zielona Gora, Poland, 13–15 September 2017; pp. 1–6. [Google Scholar]
  6. Hua, J.; Kuang, W.; Gao, Z.; Meng, L.; Xu, Z. Image denoising using 2-D FIR filters designed with DEPSO. Multimed. Tools Appl. 2014, 69, 157–169. [Google Scholar] [CrossRef]
  7. Lu, W.S. Two-Dimensional Digital Filters; CRC Press: Boca Raton, FL, USA, 1992; Volume 80. [Google Scholar]
  8. Chen, W.K. The Circuits and Filters Handbook; CRC Press: Boca Raton, FL, USA, 2018. [Google Scholar]
  9. Oppenheim, A.V.; Schafer, R.W. Discrete-Time Signal Processing; Pearson Education: Atlanta, GA, USA, 2014. [Google Scholar]
  10. Dudgeon, D.E.; Mersereau, R.M. Multidimensional Digital Signal Processing; Prentice-Hall Signal Processing Series; Prentice-Hall: Englewood Cliffs, NJ, USA, 1984. [Google Scholar]
  11. Nongpiur, R.C.; Shpak, D.J.; Antoniou, A. Improved design method for nearly linear-phase IIR filters using constrained optimization. IEEE Trans. Signal Process. 2013, 61, 895–906. [Google Scholar] [CrossRef] [Green Version]
  12. Holford, S.; Agathoklis, P. The use of model reduction techniques for designing IIR filters with linear phase in the passband. IEEE Trans. Signal Process. 1996, 44, 2396–2404. [Google Scholar] [CrossRef]
  13. Lertniphonphun, W.; McClellan, J.H. Unified design algorithm for complex FIR and IIR filters. In Proceedings of the IEEE International Conference on Acoustics, Speech, and Signal Processing, Proceedings (Cat. No. 01CH37221), Salt Lake City, UT, USA, 7–11 May 2001; Volume 6, pp. 3801–3804. [Google Scholar]
  14. O’connor, B.; Huang, T. Stability of general two-dimensional recursive digital filters. IEEE Trans. Acoust. Speech Signal Process. 1978, 26, 550–560. [Google Scholar] [CrossRef]
  15. Katayama, T. Subspace Methods for System Identification; Springer Science & Business Media: Berlin/Heidelberg, Germany, 2006. [Google Scholar]
  16. Boyd, S.; El Ghaoui, L.; Feron, E.; Balakrishnan, V. Linear Matrix Inequalities in System and Control Theory; SIAM: Philadelphia, PA, USA, 1994; Volume 15. [Google Scholar]
  17. Lim, J.S. Two-Dimensional Signal and Image Processing; Prentice Hall: Englewood Cliffs, NJ, USA, 1990; 710p. [Google Scholar]
  18. Lu, W.S.; Hinamoto, T. Optimal design of IIR digital filters with robust stability using conic-quadratic-programming updates. IEEE Trans. Signal Process. 2003, 51, 1581–1592. [Google Scholar]
  19. Lai, X.; Lin, Z. Minimax design of IIR digital filters using a sequential constrained least-squares method. IEEE Trans. Signal Process. 2010, 58, 3901–3906. [Google Scholar] [CrossRef]
  20. Kockanat, S.; Karaboga, N. The design approaches of two-dimensional digital filters based on metaheuristic optimization algorithms: A review of the literature. Artif. Intell. Rev. 2015, 44, 265–287. [Google Scholar] [CrossRef]
  21. Jiang, A.; Kwan, H.K. IIR digital filter design with new stability constraint based on argument principle. IEEE Trans. Circuits Syst. I Regul. Pap. 2009, 56, 583–593. [Google Scholar] [CrossRef]
  22. Anderson, B.; Agathoklis, P.; Jury, E.; Mansour, M. Stability and the matrix Lyapunov equation for discrete 2-dimensional systems. IEEE Trans. Circuits Syst. 1986, 33, 261–267. [Google Scholar] [CrossRef]
  23. Lodge, J.; Fahmy, M. Stability and overflow oscillations in 2-D state-space digital filters. IEEE Trans. Acoust. Speech Signal Process. 1981, 29, 1161–1171. [Google Scholar] [CrossRef]
  24. Lu, W.S.; Lee, E. Stability analysis for two-dimensional systems via a Lyapunov approach. IEEE Trans. Circuits Syst. 1985, 32, 61–68. [Google Scholar] [CrossRef]
  25. Reddy, H.; Rajan, P. A comprehensive study of two-variable Hurwitz polynomials. IEEE Trans. Educ. 1989, 32, 198–209. [Google Scholar] [CrossRef]
  26. Paraskevopoulos, P.; Panagopoulos, P.; Vaitsis, G.; Varoufakis, S.; Antoniou, G. Model reduction of 2-D systems via orthogonal series. Multidimens. Syst. Signal Process. 1991, 2, 69–83. [Google Scholar] [CrossRef]
  27. Guo, T.Y.; Hwang, C.; Shieh, L.S.; Chen, C.H. Reduced-order models of 2-D linear discrete separable-denominator system using bilinear Routh approximations. In IEE Proceedings G (Electronic Circuits and Systems); IET: London, UK, 1992; Volume 139, pp. 45–56. [Google Scholar]
  28. Knorn, S.; Middleton, R.H. Stability of two-dimensional linear systems with singularities on the stability boundary using LMIs. IEEE Trans. Autom. Control 2013, 58, 2579–2590. [Google Scholar] [CrossRef]
  29. Lai, X.; Meng, H.; Cao, J.; Lin, Z. A sequential partial optimization algorithm for minimax design of separable-denominator 2-D IIR filters. IEEE Trans. Signal Process. 2017, 65, 876–887. [Google Scholar] [CrossRef]
  30. Dumitrescu, B. Optimization of two-dimensional IIR filters with nonseparable and separable denominator. IEEE Trans. Signal Process. 2005, 53, 1768–1777. [Google Scholar] [CrossRef]
  31. Gorinevsky, D.; Boyd, S. Optimization-based design and implementation of multidimensional zero-phase IIR filters. IEEE Trans. Circuits Syst. I Regul. Pap. 2006, 53, 372–383. [Google Scholar] [CrossRef] [Green Version]
  32. Wysocka-Schillak, F. Design of separable 2-D IIR Filters with approximately linear phase in the passband using genetic algorithm. In Proceedings of the 2008 Conference on Human System Interactions, IEEE, Krakow, Poland, 25–27 May 2008; pp. 66–70. [Google Scholar]
  33. Lu, W.S. A unified approach for the design of 2-D digital filters via semidefinite programming. IEEE Trans. Circuits Syst. Fundam. Theory Appl. 2002, 49, 814–826. [Google Scholar] [CrossRef]
  34. Mladenov, V.M.; Mastorakis, N.E. Design of two-dimensional recursive filters by using neural networks. IEEE Trans. Neural Netw. 2001, 12, 585–590. [Google Scholar] [CrossRef] [PubMed]
  35. Mastorakis, N.E.; Gonos, I.F.; Swamy, M. Design of two-dimensional recursive filters using genetic algorithms. IEEE Trans. Circuits Syst. Fundam. Theory Appl. 2003, 50, 634–639. [Google Scholar] [CrossRef] [Green Version]
  36. Tsai, J.T.; Ho, W.H.; Chou, J.H. Design of two-dimensional IIR digital structure-specified filters by using an improved genetic algorithm. Expert Syst. Appl. 2009, 36, 6928–6934. [Google Scholar] [CrossRef]
  37. Liang, L.; Ahmadi, M.; Sid-Ahmed, M. Design of 2D IIR filters with canonical signed-digit coefficients using genetic algorithm. In Proceedings of the 2003 46th Midwest Symposium on Circuits and Systems, IEEE, Cairo, Egypt, 27–30 December 2003; Volume 2, pp. 633–635. [Google Scholar]
  38. Kawamata, M.; Imakubo, J.; Higuchi, T. Optimal design method of 2-D IIR digital filters based on a simple genetic algorithm. In Proceedings of the Proceedings of 1st International Conference on Image Processing, IEEE, Austin, TX, USA, 13–16 November 1994; Volume 1, pp. 780–784. [Google Scholar]
  39. Lv, C.; Yan, S.; Cheng, G.; Xu, L.; Tian, X. Design of two-dimensional IIR digital filters by using a novel hybrid optimization algorithm. Multidimens. Syst. Signal Process. 2017, 28, 1267–1281. [Google Scholar] [CrossRef]
  40. Shpak, D.J.; Antoniou, A. A generalized Remez method for the design of FIR digital filters. IEEE Trans. Circuits Syst. 1990, 37, 161–174. [Google Scholar] [CrossRef]
  41. Charalambous, C. Design of 2-dimensional circularly-symmetric digital filters. In IEE Proceedings G (Electronic Circuits and Systems); IET: London, UK, 1982; Volume 129, pp. 47–54. [Google Scholar]
  42. Lu, W.S.; Pei, S.C.; Tseng, C.C. A weighted least-squares method for the design of stable 1-D and 2-D IIR digital filters. IEEE Trans. Signal Process. 1998, 46, 1–10. [Google Scholar] [CrossRef]
  43. Miyata, T.; Aikawa, N.; Sugita, Y.; Yoshikawa, T. A design method for separable-denominator 2D IIR filters with a necessary and sufficient stability check. IEICE Trans. Fundam. Electron. Commun. Comput. Sci. 2009, 92, 307–310. [Google Scholar] [CrossRef]
  44. Omar, A.; Shpak, D.; Agathoklis, P. Improved Design Method for Nearly Linear-Phase IIR Filters Using Constrained Optimization. J. Circuits Syst. Comput. 2021, 30, 2150207. [Google Scholar] [CrossRef]
  45. Omar, A.; Shpak, D.; Agathoklis, P. Nearly Linear-Phase 2-D Recursive Digital Filters Design Using Balanced Realization Model Reduction. J. Circuits Syst. Comput. 2022; submitted in October. [Google Scholar]
  46. Kim, K.; Kim, J.; Nam, S. Design of computationally efficient 2D FIR filters using sampling-kernel-based interpolation and frequency transformation. Electron. Lett. 2015, 51, 1326–1328. [Google Scholar] [CrossRef]
  47. Aly, S.; Fahmy, M. Design of two-dimensional recursive digital filters with specified magnitude and group delay characteristics. IEEE Trans. Circuits Syst. 1978, 25, 908–916. [Google Scholar] [CrossRef]
  48. Antoniou, A.; Lu, W.S. Design of two-dimensional digital filters by using the singular value decomposition. IEEE Trans. Circuits Syst. 1987, 34, 1191–1198. [Google Scholar] [CrossRef]
  49. Guindon, D.; Shpak, D.J.; Antoniou, A. Design methodology for nearly linear-phase recursive digital filters by constrained optimization. IEEE Trans. Circuits Syst. I Regul. Pap. 2009, 57, 1719–1731. [Google Scholar] [CrossRef]
  50. Andersen, E.D.; Gondzio, J.; Mészáros, C.; Xu, X. Implementation of Interior Point Methods for Large Scale Linear Programming; HEC/Université de Geneve: Dordrecht, The Netherlands, 1996. [Google Scholar]
  51. Liu, L.; Fan, L. The complexity analysis of an efficient interior-point algorithm for linear optimization. In Proceedings of the 2010 Third International Joint Conference on Computational Science and Optimization, IEEE, Huangshan, China, 28–31 May 2010; Volume 2, pp. 21–24. [Google Scholar]
  52. Wirski, R. Synthesis and Realization of Two-Dimensional Separable Denominator Orthogonal Systems via Decomposition Into 1-D Systems. IEEE Trans. Circuits Syst. I Regul. Pap. 2019, 66, 4309–4322. [Google Scholar] [CrossRef]
  53. Reddy, H.C.; Khoo, I.H.; Rajan, P. 2-D symmetry: Theory and filter design applications. IEEE Circuits Syst. Mag. 2003, 3, 4–33. [Google Scholar] [CrossRef]
Figure 1. Flowchart of the optimization process.
Figure 1. Flowchart of the optimization process.
Signals 04 00030 g001
Figure 2. Ideal specifications for 2D highpass filter in [8].
Figure 2. Ideal specifications for 2D highpass filter in [8].
Signals 04 00030 g002
Figure 3. Magnitude response of the 2D highpass filters [8].
Figure 3. Magnitude response of the 2D highpass filters [8].
Signals 04 00030 g003
Figure 4. Group delays in the passband of the obtained highpass filter ( Γ g = 1.65 ).
Figure 4. Group delays in the passband of the obtained highpass filter ( Γ g = 1.65 ).
Signals 04 00030 g004
Figure 5. Group delays in the passband of the highpass filter in [8].
Figure 5. Group delays in the passband of the highpass filter in [8].
Signals 04 00030 g005
Figure 6. Ideal specifications for 2D lowpass filter in [8].
Figure 6. Ideal specifications for 2D lowpass filter in [8].
Signals 04 00030 g006
Figure 7. Magnitude response of the 2D lowpass filter [8].
Figure 7. Magnitude response of the 2D lowpass filter [8].
Signals 04 00030 g007
Figure 8. Group delays in the passband of the obtained lowpass filter ( Γ g = 2.35 ).
Figure 8. Group delays in the passband of the obtained lowpass filter ( Γ g = 2.35 ).
Signals 04 00030 g008
Figure 9. First quadrant group delays for lowpass ( Γ g = 2.35 ).
Figure 9. First quadrant group delays for lowpass ( Γ g = 2.35 ).
Signals 04 00030 g009
Figure 10. Group delays in the passband of the lowpass filter in [8].
Figure 10. Group delays in the passband of the lowpass filter in [8].
Signals 04 00030 g010
Figure 11. First quadrant group delays of the lowpass in [8].
Figure 11. First quadrant group delays of the lowpass in [8].
Signals 04 00030 g011
Figure 12. Ideal specifications for a 2D bandpass filter in [8].
Figure 12. Ideal specifications for a 2D bandpass filter in [8].
Signals 04 00030 g012
Figure 13. Magnitude response of the bandpass filter [8].
Figure 13. Magnitude response of the bandpass filter [8].
Signals 04 00030 g013
Figure 14. Group delays in the passband of the obtained bandpass filter.
Figure 14. Group delays in the passband of the obtained bandpass filter.
Signals 04 00030 g014
Figure 15. Group delays of the obtained bandpass filter in the first quadrant.
Figure 15. Group delays of the obtained bandpass filter in the first quadrant.
Signals 04 00030 g015
Figure 16. Group delays in the passband of the bandpass filter in [8].
Figure 16. Group delays in the passband of the bandpass filter in [8].
Signals 04 00030 g016
Figure 17. Group delays of the bandpass filter in [8] in the first quadrant.
Figure 17. Group delays of the bandpass filter in [8] in the first quadrant.
Signals 04 00030 g017
Table 1. Comparison results for highpass filter.
Table 1. Comparison results for highpass filter.
Proposed MethodMethod in [8]
Orders(4,4)(4,4)
H min 0.4370.262
H max 1.3381.066
Q h 50.7660.54
H max t 0.6640.849
H max s 0.3890.426
Q s 8.104 2.598
τ min 0.9620.009
τ max 4.39711.614
Q τ 64.0999.84
Table 2. Comparison results for lowpass filter.
Table 2. Comparison results for lowpass filter.
Proposed MethodMethod in [8]
Orders(4,4)(4,4)
H min 0.8430.642
H max 1.1201.069
Q h 14.2525.26
H max t 1.00791.020
H max s 0.4030.446
Q s 0.9162.275
τ min 2.2171.738
τ max 6.86454.091
Q τ 47.65093.774
Table 3. Comparison results for bandpass filter.
Table 3. Comparison results for bandpass filter.
Proposed MethodMethod in [8]
Orders(4,4)(4,4)
H min 0.8090.633
H max 1.1911.167
Q h 19.10029.667
H max t 0.8550.798
H max s 0.3240.268
Q s 0.6680.734
τ min 1.8221.828
τ max 5.0736.737
Q τ 47.15057.315
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

Omar, A.; Shpak, D.; Agathoklis, P.; Moa, B. A Nonlinear Optimization Design Algorithm for Nearly Linear-Phase 2D IIR Digital Filters. Signals 2023, 4, 575-590. https://doi.org/10.3390/signals4030030

AMA Style

Omar A, Shpak D, Agathoklis P, Moa B. A Nonlinear Optimization Design Algorithm for Nearly Linear-Phase 2D IIR Digital Filters. Signals. 2023; 4(3):575-590. https://doi.org/10.3390/signals4030030

Chicago/Turabian Style

Omar, Abdussalam, Dale Shpak, Panajotis Agathoklis, and Belaid Moa. 2023. "A Nonlinear Optimization Design Algorithm for Nearly Linear-Phase 2D IIR Digital Filters" Signals 4, no. 3: 575-590. https://doi.org/10.3390/signals4030030

APA Style

Omar, A., Shpak, D., Agathoklis, P., & Moa, B. (2023). A Nonlinear Optimization Design Algorithm for Nearly Linear-Phase 2D IIR Digital Filters. Signals, 4(3), 575-590. https://doi.org/10.3390/signals4030030

Article Metrics

Back to TopTop