Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
Next Article in Journal
A Modified Artificial Bee Colony Algorithm Based on the Self-Learning Mechanism
Previous Article in Journal
PHEFT: Pessimistic Image Processing Workflow Scheduling for DSP Clusters
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Novel Design of Sparse Prototype Filter for Nearly Perfect Reconstruction Cosine-Modulated Filter Banks

1
School of Electronics and Information Engineering, Tianjin Polytechnic University, Tianjin 300387, China
2
Tianjin Key Laboratory of Optoelectronic Detection Technology and System, Tianjin 300387, China
3
College of Electronic Information and Optical Engineering, Nankai University, Tianjin 300071, China
*
Author to whom correspondence should be addressed.
Algorithms 2018, 11(5), 77; https://doi.org/10.3390/a11050077
Submission received: 15 April 2018 / Revised: 8 May 2018 / Accepted: 21 May 2018 / Published: 22 May 2018

Abstract

:
Cosine-modulated filter banks play a major role in digital signal processing. Sparse FIR filter banks have lower implementation complexity than full filter banks, while keeping a good performance level. This paper presents a fast design paradigm for sparse nearly perfect-reconstruction (NPR) cosine-modulated filter banks. First, an approximation function is introduced to reduce the non-convex quadratically constrained optimization problem to a linearly constrained optimization problem. Then, the desired sparse linear phase FIR prototype filter is derived through the orthogonal matching pursuit (OMP) performed under the weighted l 2 norm. The simulation results demonstrate that the proposed scheme is an effective paradigm to design sparse NPR cosine-modulated filter banks.

1. Introduction

Maximally decimated, uniform, cosine-modulated filter banks play a major role in applications where signals must be processed in sub-bands [1,2,3,4,5]. There are two main methods obtaining a cosine-modulated filter bank. (1) Design analysis and synthesis filter banks respectively under the certain conditions; (2) Design one or two prototype filters and get the analysis and synthesis filter banks by cosine modulating. In this paper, we choose the later one because the design is straightforward, e.g., this method focus on designing a single prototype filter.
There are several methods to design the prototype filter of cosine-modulated filter banks [6,7,8]. The weighted constrained least square (WCLS) algorithm proposed in reference [6] casts the design problem as a linear minimization of prototype filter coefficients, which makes the value at ω = π / 2 M is 0.707 . However, the order of the prototype filter is very high. It is generally true that group delay increases with the order of the filter. The quantum-behaved particle swarm optimization (QBPSO) algorithm proposed in reference [7] appropriately relaxes the limitation of perfect-reconstruction conditions. This method designs the prototype filter by making an unconstrained optimization to adjust the pass band cutoff frequency of objective function, which minimizes the cost function to meet the reconstruction conditions. Although the amplitude distortion is guaranteed, the aliasing distortion cannot meet the perfect-reconstruction condition of filter banks, and the attenuation of prototype filter is restricted. The design method based on gradient information (GI) in reference [8] uses the gradient vector of objective function to design prototype filter. However, the order of prototype coefficient is high, the aliasing distortion and the amplitude distortion of filter banks are not controlled well. The order of NPR cosine-modulated filter bank designed by these methods is high and not sparse, which makes an increasing number of the multipliers and adders in practical application, and results in a high hardware cost and power consumption.
In this paper, we present a fast design paradigm for sparse NPR cosine-modulated filter banks. A sparse FIR filter bank contains a large number of zero coefficients, such that the multipliers and adders corresponding to these zero coefficients are not required, which results in a lower hardware cost and power consumption with allowed accuracy of distortion. The design procedure can be divided into two stages. In the first stage, an approximation function f ( ω ) is introduced to reduce the non-convex quadratically constrained optimization problem to a linearly constrained optimization problem. The desired sparse linear phase FIR prototype filter can be derived through the orthogonal matching pursuit [9] performed under the weighted l 2 norm in the second stage. The simulation results demonstrate that the proposed scheme is an effective paradigm to design sparse NPR cosine-modulated filter banks.
The paper organization is as follows: In Section 2, we briefly summarize the design framework of perfect-reconstruction cosine-modulated filter banks. The design of sparse linear phase FIR prototype filter of NPR cosine-modulated filter banks is showed in Section 3. The representative experimental results and comparison data are discussed in Section 4. The conclusion of proposed scheme is showed in Section 5.

2. Perfect-Reconstruction Cosine-Modulated Filter Banks

A M-channel maximally decimated filter bank as illustrated in Figure 1 is composed of M finite impulse response duration analysis filters with impulse responses h k ( n ) , the down- and up- sampling by factor M of sub-band signals, and synthesis FIR filters g k ( n ) . The impulse response of filter bank is obtained by cosine modulation with a low pass filter:
h k ( n ) = 2 h ( n ) cos [ π M ( k + 1 2 ) ( n V 2 ) + ( 1 ) k π M ] ,
g k ( n ) = 2 h ( n ) cos [ π M ( k + 1 2 ) ( n V 2 ) ( 1 ) k π M ]
with 0 k M 1 and 0 n N 1 . h ( n ) is the impulse response of the FIR prototype filter, V denotes the system delay and N is the order of the prototype filter. The input-output relation of the system in the z domain is given by
X ^ ( z ) = A 0 ( z ) X ( z ) + l = 1 M 1 A l ( z ) X ( z e j 2 π l / M )
where
A 0 ( z ) = 1 M k = 0 M 1 H k ( z ) G k ( z )
is the distortion transfer function which determining the distortion of non-aliased component X ( z ) caused by the system, and
A l ( z ) = 1 M k = 0 M 1 H k ( z e j 2 π l / M ) G k ( z )
are the aliasing transfer functions which determine the attenuation of aliased components X ( z e j 2 π l / M ) for l = 1 , , M 1 . It follows that the filter bank holds the perfect-reconstruction property if and only if
A 0 ( z ) = z V ,
A l ( z ) = 0 , l = 1 , , M 1 .

3. The Design of Sparse Linear Phase FIR Prototype Filter of NPR Cosine-Modulated Filter Banks

Without loss of generality, the linear phase FIR filter H ( e j ω ) is assumed to be the Type I I which suggests that the length N of H ( e j ω ) is even and its tap weights satisfying h ( n ) = h ( N 1 n ) for all 0 n N 1 . For the other types of the linear phase filters, an argument similar to the one developed in this paper can be followed. The frequency response corresponding to an N-th order Type I I linear phase FIR prototype filter can be represented as
H ( e j ω ) = e j N 1 2 ω H 0 ( ω ) = e j N 1 2 ω 2 n = 0 N 2 1 h ( n ) cos ( ( N 1 2 n ) ω ) .
where H 0 ( ω ) is the real-valued amplitude response of prototype filter. The perfect-reconstruction conditions Equations (6) and (7) for linear phase cosine-modulated filter banks can be reduced to
| H 0 ( e j ω ) | 2 + | H 0 ( e j ( ω π M ) ) | 2 = 1 , ω [ 0 , π M ] ,
| H 0 ( e j ω ) | = 0 , ω ( π M , π ] .
Thus, the problem of sparse linear phase FIR prototype filter design for cosine-modulated filter banks can be formulated as
min h 0
s . t . | c ( ω ) h | 2 + | c ( ω π M ) h | 2 = 1 , ω [ 0 , π M ] ,
| c ( ω ) h | = 0 , ω ( π M , π ]
where we have
c ( ω ) = [ cos N 1 2 ω , , cos ( N 1 2 n ) ω , , cos 1 2 ω ] ,
h = 2 [ h ( 0 ) , h ( 1 ) , , h ( n ) , , h ( N 2 1 ) ] T .
Define a approximation function as f ( ω ) , which is required to satisfy the following equations:
f ( ω 0 ) = 1 ,
f 2 ( ω ) + f 2 ( ω π M ) = 1 , ω ( ω 0 , π M ω 0 ) ,
f ( π M ω 0 ) = 0
where ω 0 = α π / 2 M and 0 α 1 . There are lots of methods to obtain the function f ( ω ) . In this paper, we establish a mathematical model like
f ω = cos κ ω ϕ .
Substituting Equation (14a), Equation (14c) and f π / 2 M = 2 / 2 into Equation (15), we can obtain
κ = M 2 ( 1 α ) , ϕ = M ω 0 2 ( 1 α ) .
Utilizing the approximation function f ( ω ) , the problem of designing sparse NPR cosine-modulated filter banks can be expressed as a constrained optimization problem for the coefficients of prototype filter in the following form:
min h 0
s . t . | c ( ω ) h 1 | δ p , ω [ 0 , ω 0 ) ,
| c ( ω ) h f ( ω ) | δ t , ω [ ω 0 , π M ω 0 ] ,
| c ( ω ) h | δ s , ω ( π M ω 0 , π ] ,
c ( π 2 M ) h = 2 2
where δ p δ t and δ s represent the ripples of pass band, transition band and stop band respectively. To compute a solution of problem Equation (17), we follow the standard discretized procedure and replace the continuous parameter ω with discretized samples distributed in the frequency set [ 0 , π ] . The L × N 2 discretized sampling matrix B can be constructed as
B = [ B p ; B t ; B s ] ,
B p = c ( ω 1 ( p ) ) c ( ω L p ( p ) ) , [ ω 1 ( p ) , ω 2 ( p ) , , ω L p ( p ) ] 0 , ω 0 , B t = c ( ω 1 ( t ) ) c ( ω L t ( t ) ) , [ ω 1 ( t ) , ω 2 ( t ) , , ω L t ( t ) ] ω 0 , π M ω 0 , B s = c ( ω 1 ( s ) ) c ( ω L s ( s ) ) , [ ω 1 ( s ) , ω 2 ( s ) , , ω L s ( s ) ] π M ω 0 , π
with L = L p + L t + L s . Define L × 1 vector d and e as
d = [ 1 1 f ( ω 1 ( t ) ) f ( ω L t ( t ) ) 0 0 ] T ,
e = [ δ p δ p δ t δ t δ s δ s ] T .
The discretized and normalized formulation of problem Equation (17) can be expressed as:
min h 0
s . t . Bh d e
c ( π 2 M ) h = 2 2
Choose the initial L × 1 weight vector w ( 1 ) = [ 1 , , 1 ] T . At the k iteration with k 1 , the procedure of computing the sparse FIR prototype filter proceeds as follows.
Step 1:
Normalize column vectors b 0 , b 1 , , b N 2 1 of B under the weight vector w ( k ) as
B ( k ) = [ b 0 ( k ) , b 1 ( k ) , , b n ( k ) , , b N 2 1 ( k ) ] ,
b n ( k ) = b n / [ i = 0 L w ( k ) ( l ) | b n ( l ) | 2 ] 1 2 , 0 n N 2 1 .
Step 2:
The OMP algorithm [9] is employed to solve the following problem
min h ( k ) B ( k ) h ( k ) d 2 2
s . t . h ( k ) 0 k .
Let k × 1 vector s ( k ) denote the final solution of Equation (22), and the final L × 1 residual vector is given as
r ( k ) = Φ ( k ) s ( k ) d .
where Φ ( k ) = [ b n 1 ( k ) , b n 2 ( k ) , , b n k ( k ) ] is the matrix of chosen atoms and Λ ( k ) = { n 1 , n 2 , , n k } is the nonzero index set.
Step 3:
Solve the following linear programming problem:
min h ( k ) , μ μ
s . t . | B ( k ) h ( k ) d | e + μ · 1 L × 1 ,
c ( π 2 M ) h ( k ) = 2 2
h n ( k ) = 0 , n { 0 , 1 , , N 2 1 } Λ ( k ) .
If μ 0 , then the computed filter h ^ ( k ) is the sparse FIR prototype filter of cosine-modulated filter banks which meeting our design specifications.
Step 4:
Compute the new weight vector by formula from reference [10]
w l ( k + 1 ) = 1 + r l ( k ) max { h ^ ( k ) 100 } 2 1 4 , 1 l L
where w l ( k + 1 ) and r l ( k ) are the i-th entries of the new weight vector w ( k + 1 ) and the residual vector r ( k ) respectively. Replace w ( k ) with w ( k + 1 ) and repeat the procedure from Steps 1 to 4.

4. Examples

In this section, the performance of the proposed scheme is evaluated by several design examples on a Intel Core i7-6700 3.40GHz PC with MATLAB, and compares to other NPR cosine-modulated filter banks design algorithms.
Example 1.
This example is concerned with the design of NPR cosine-modulated filter banks with M = 4 . The filter order of the FIR linear phase prototype filter is N = 140 . The numbers of discretized samples for pass band, transition band and stop band are L p = 6 , L t = 20 and L s = 800 , respectively. Choose δ p = δ t = δ s = 1 × 10 3 . The amplitude distortion function magnitude e a m ( ω ) = 1 A 0 ( e j ω ) , the total aliasing transfer function magnitude e a ( ω ) = l = 1 M 1 A l e j ω 2 , and the amplitude response of the sparse prototype filter derived from the proposed scheme are shown in Figure 2. The amplitude distortion e a m = max ω e a m ( ω ) and total aliasing error e a = max ω e a ( ω ) of the corresponding cosine-modulated filter bank are e a m = 1.70 × 10 3 and e a = 5.06 × 10 7 , respectively. It takes 31 iterations for proposed method to converge, with the CPU time equals to 0.2716 seconds in each iteration.
We compare the performance of the proposed scheme with the WCLS [6] scheme and QBPSO [7] scheme. At the similar reconstruction performances, the number of nonzero tap weights of the prototype filter yielded from the proposed method is 18 less than that of the WCLS scheme and the QBPSO scheme. In Table 1, we compare the filter order, the number of nonzero tap weights, the amplitude distortion e a m and the total aliasing error e a corresponding to the proposed scheme, the WCLS scheme and the QBPSO scheme.
Example 2.
This example is concerned with the design of NPR cosine-modulated filter banks with M = 8. The filter order of the FIR linear phase prototype filter is N = 160. The numbers of discretized samples for pass band, transition band and stop band are L p = 6, L t = 93 and L s = 91, respectively. Choose δ p = δ t = δ s = 1 × 10 3 . The amplitude distortion function magnitude e a m ( ω ) , the total aliasing transfer function magnitude e a ( ω ) , and the amplitude response of the sparse prototype filter derived from the proposed scheme are shown in Figure 3. The amplitude distortion e a m and total aliasing error e a of the corresponding cosine-modulated filter bank are e a m = 5.2 × 10 4 and e a = 9.85 × 10 6 . It takes 49 iterations for proposal method to converge, with the CPU time equals to 0.0591 seconds in each iteration.
We compare the performance of the proposed scheme with the WCLS [6] scheme and GI [8] scheme. At the similar reconstruction performances, the number of nonzero tap weights of the prototype filter yielded from the proposed method is 46 less than that of the WCLS scheme, the number of nonzero tap weights of the prototype filter yielded from the proposed method is 34 less than that of the GI scheme. In Table 1, we compare the filter order, the number of nonzero tap weights, the amplitude distortion e a m and the total aliasing error e a corresponding to the proposed scheme, the WCLS scheme and the GI scheme.
Example 3.
This example is concerned with the design of NPR cosine-modulated filter banks with M = 16 . The filter order of the FIR linear phase prototype filter is N = 254 . The numbers of discretized samples for pass band, transition band and stop band are L p = 4 , L t = 7 and L s = 94 , respectively. Choose δ p = δ t = δ s = 1 × 10 3 . The amplitude distortion function magnitude e a m ( ω ) , the total aliasing transfer function magnitude e a ( ω ) , and the amplitude response of the sparse prototype filter derived from the proposed scheme are shown in Figure 4. The amplitude distortion e a m and total aliasing error e a of the corresponding cosine-modulated filter bank are e a m = 5.21 × 10 4 and e a = 4.24 × 10 6 . It takes 83 iterations for proposed method to converge, with the CPU time equals to 0.1108 seconds in each iteration.
We compare the performance of the proposed scheme with the WCLS [6] scheme and QBPSO [7] scheme. At the similar reconstruction performances, the number of nonzero tap weights of the prototype filter yielded from the proposed method is 58 less than that of the WCLS scheme, the number of nonzero tap weights of the prototype filter yielded from the proposed method is 90 less than that of the QBPSO scheme. In Table 1, we compare the filter order, the number of nonzero tap weights, the amplitude distortion e a m and the total aliasing error e a corresponding to the proposed scheme, the WCLS scheme and the QBPSO scheme.

5. Conclusions

In this paper, we present a fast design paradigm for sparse NPR cosine-modulated filter banks. Different to traditional method of designing filter banks, the proposed algorithm aims to add a sparsity into filter banks while others have not, and ensures the stop-band attenuation at least 60 dB while most other methods cannot guarantee those with the same performance of amplitude distortion and aliasing distortion. The desired sparse linear phase FIR prototype filter can be derived through the orthogonal matching pursuit [9] performed under the weighted l 2 norm. The simulation results demonstrate that the proposed scheme is an effective paradigm to design sparse NPR cosine-modulated filter banks.

Author Contributions

W.X. and Y.L. contributed to the overall idea and algorithm, coded the algorithm in MATLAB and contributed to the detailed writing of the manuscript. J.M., J.Z. and X.G. contributed to the revision, preparation, and publishing of the paper. All authors have read and approved the final manuscript.

Acknowledgments

This work was supported by the Nature Science Foundation under Grant 61501324 and the Nature Science Foundation of Tianjin under Grant 14JCYBJC16100.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
NPRNearly Perfect-Reconstruction
WCLSWeighted Constrained Least Square
QBPSOQuantum-Behaved Particle Swarm Optimization
GIgradient information
OMPorthogonal matching pursuit

References

  1. Hang, S.; Tian, L.; Ma, X.; Wei, Y. A reconfigurable sound wave decomposition filter bank for hearing aids based on nonlinear transformation. IEEE Tran. Biomed. Circuits Syst. 2016, 10, 487–496. [Google Scholar] [CrossRef] [PubMed]
  2. Ling, B.W.; Ho, C.Y.; Teo, K.; Siu, W.; Cao, J.; Dai, Q. Optimal design of cosine modulated nonuniform linear phase FIR filter bank via both stretching and shifting frequency response of single prototype filter. IEEE Trans. Signal Proc. 2014, 62, 2517–2530. [Google Scholar] [CrossRef]
  3. Kyochi, S.; Ikehara, M. A class of near shift-invariant and orientation-selective transform based on delay-less oversampled even-stacked cosine-modulated filter banks. IEICE Trans. Fundam. Electron. Commun. Comput. Sci. 2010, E93.A, 724–733. [Google Scholar] [CrossRef]
  4. Zhao, N.; Pu, F.; Xu, X.; Chen, N. Cognitive wideband spectrum sensing using cosine-modulated filter banks. Int. J. Electron. 2015, 102, 1890–1901. [Google Scholar] [CrossRef]
  5. Sakiyama, A.; Watanabe, K.; Tanaka, Y. Spectral graph wavelets and filter banks with low approximation error. IEEE Trans. Signal Inf. Proc. Over Netw. 2016, 2, 230–245. [Google Scholar] [CrossRef]
  6. Kumar, A.; Singh, G.K.; Anand, R.S. A simple design method for the cosine-modulated filter banks using weighted constrained least square technique. J. Frankl. Inst. Geophys. 2011, 348, 606–621. [Google Scholar] [CrossRef]
  7. Tan, F.; Zhang, T.; Gao, C.; Huang, L. Optimal design of cosine modulated filter banks using quantum-behaved particle swarm optimization algorithm. Image Signal Proc. 2011, 5, 2280–2284. [Google Scholar] [CrossRef]
  8. Zhang, Z.J. Efficient design of cosine modulated filter banks based on gradient information. IEEE Signal Proc. Lett. 2007, 14, 940–943. [Google Scholar] [CrossRef]
  9. 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]
  10. Bube, K.P. Hybrid l1/l2 minimization with applications to tomography. Geophysics 1997, 62, 1183. [Google Scholar] [CrossRef]
Figure 1. M-channel maximally decimated filter bank.
Figure 1. M-channel maximally decimated filter bank.
Algorithms 11 00077 g001
Figure 2. M = 4 channel cosine-modulated filter banks in example 1. (a) The amplitude distortion function magnitude e a m ( ω ) ; (b) The total aliasing transfer function magnitude e a ( ω ) ; (c) Amplitude response of sparse prototype filter.
Figure 2. M = 4 channel cosine-modulated filter banks in example 1. (a) The amplitude distortion function magnitude e a m ( ω ) ; (b) The total aliasing transfer function magnitude e a ( ω ) ; (c) Amplitude response of sparse prototype filter.
Algorithms 11 00077 g002
Figure 3. M = 8 channel cosine-modulated filter banks in example 2. (a) The amplitude distortion function magnitude e a m ( ω ) ; (b) The total aliasing transfer function magnitude e a ( ω ) ; (c) Amplitude response of sparse prototype filter.
Figure 3. M = 8 channel cosine-modulated filter banks in example 2. (a) The amplitude distortion function magnitude e a m ( ω ) ; (b) The total aliasing transfer function magnitude e a ( ω ) ; (c) Amplitude response of sparse prototype filter.
Algorithms 11 00077 g003
Figure 4. M = 16 channel cosine-modulated filter banks in Example 3. (a) The amplitude distortion function magnitude e a m ( ω ) ; (b) The total aliasing transfer function magnitude e a ( ω ) ; (c) Amplitude response of sparse prototype filter.
Figure 4. M = 16 channel cosine-modulated filter banks in Example 3. (a) The amplitude distortion function magnitude e a m ( ω ) ; (b) The total aliasing transfer function magnitude e a ( ω ) ; (c) Amplitude response of sparse prototype filter.
Algorithms 11 00077 g004
Table 1. Comparisons for Several Parameters of the Designed Filters Derived from Different Methods.
Table 1. Comparisons for Several Parameters of the Designed Filters Derived from Different Methods.
MDesign MethodFilter OrderNumber of Nonzero TapsAmplitude DistortionTotal Aliasing Error
WCLS8080 3.10 × 10 3 1.06 × 10 6
4QBPSO8080 7.01 × 10 4 8.49 × 10 8
Proposed Method14062 1.70 × 10 3 5.06 × 10 7
WCLS144144 1.80 × 10 3 3.90 × 10 7
8GI132132 2.14 × 10 5 1.30 × 10 5
Proposed Method16098 5.20 × 10 4 9.85 × 10 6
WCLS224224 1.30 × 10 3 9.03 × 10 7
16QBPSO256256 4.73 × 10 4 2.34 × 10 6
Proposed Method254166 5.21 × 10 4 4.24 × 10 6

Share and Cite

MDPI and ACS Style

Xu, W.; Li, Y.; Miao, J.; Zhao, J.; Gao, X. A Novel Design of Sparse Prototype Filter for Nearly Perfect Reconstruction Cosine-Modulated Filter Banks. Algorithms 2018, 11, 77. https://doi.org/10.3390/a11050077

AMA Style

Xu W, Li Y, Miao J, Zhao J, Gao X. A Novel Design of Sparse Prototype Filter for Nearly Perfect Reconstruction Cosine-Modulated Filter Banks. Algorithms. 2018; 11(5):77. https://doi.org/10.3390/a11050077

Chicago/Turabian Style

Xu, Wei, Yi Li, Jinghong Miao, Jiaxiang Zhao, and Xin Gao. 2018. "A Novel Design of Sparse Prototype Filter for Nearly Perfect Reconstruction Cosine-Modulated Filter Banks" Algorithms 11, no. 5: 77. https://doi.org/10.3390/a11050077

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