← Back to Patent Explorer
US5436940A Expired / Lapsed US ⏱️ 0.0 years remaining

Quadrature mirror filter banks and method

Assignee: Massachusetts Institute of Technology
•
Filing: 1992-06-11
•
Granted: 1995-07-25
•
Expires: 2012-06-11
📄 Original PDF 🌐 Google Patents 🔄 Sync Now

Abstract

Abstract M-channel pseudo-quadrature-mirror-filter (QMF) banks using a prototype filter having a linear-phase spectral-factor of a 2M th band filter. The overall transfer function of the analysis filter/synthesis filter system is a delay, and the aliasing cancellation has all the significant aliasing terms canceled. Consequently, the aliasing level at the output of the pseudo-QMF banks is comparable to the stopband attenuation of the prototype filter, with the error at the output of the analysis filter/synthesis filter system approximately equal to the aliasing error at the level of the stopband attenuation. The pseudo-QMF banks have the stopband attenuation of the analysis filters and thus synthesis filters of -100 dB. The resulting reconstruction error is also on the order of -100 dB. Optimization of the pseudo-QMF banks by a quadratic-constrained least-squares formulation converges very fast as both a cost function and constraints are quadratic functions with respect to unknown parameters, providing a much higher stopband attenuation compared to previous filter banks.

Inventors

  • 👤 Truong Q. Nguyen

Bibliographic Details

Application Number —
Jurisdiction US
Data Source google_patents
Last Synchronized 2026-11-01 09:03
CPC Classifications
H; H04; H04B; H04B1/00; H04B1/66; H04B1/667

Claims Analysis

Total Claims: 16
Claim 1
Independent Claim 510 words
A pseudo-quadrature-mirror filter bank for near-perfect-reconstruction pseudo-quadrature-mirror filtering of an input signal, comprising: a plurality of analysis filters, each of said plurality of analysis filters including, a first delay chain, operatively coupled to the input signal, forming a set of 2M parallel paths for buffering the input signal; a first cascade of 2M polyphase components of an impulse transfer function H(z) of impulse response, h(n), operatively coupled to said first delay chain; and means, operatively coupled to said first cascade of 2M polyphase components, for generating a 2M-point Discrete Fourier Transform (DFT) to implement a 2M-point Discrete Cosine Transform (DCT) of the input signal, with each analysis filter having an impulse response, h k (n), of a k th analysis filter, where M is the number of subband signals, obtained by cosine-modulating an impulse response, h(n), of a prototype filter with linear phase, according to: ##EQU65## and N is the length of the impulse response, h(n), of the prototype filter; a plurality of synthesis filters, each of said plurality of synthesis filters including, a second delay chain, operatively coupled to the input signal, forming a set of 2M parallel paths for buffering the input signal; a second cascade of 2M polyphase components of an impulse transfer function H(z) of impulse response, h(n), operatively coupled to said second delay chain; and means, operatively coupled to said second cascade of 2M polyphase components, for generating a 2M-point Discrete Fourier Transform (DFT) to implement a 2-point Discrete Cosine Transform (DCT) of the input signal, with each of said plurality of synthesis filters operatively coupled to a respective one of said plurality of analysis filters, each synthesis filter having an impulse response, f k (n), of a k th synthesis filter, obtained by cosine-modulating the impulse response, h(n), of the prototype filter according to: ##EQU66## and N is the length of the impulse response, h(n), of the prototype filter; and wherein each impulse response, h(n), is found in accordance with: ##STR5## for even N, where n=2M(m-l)+2m 1 -1 and x is the greatest integer less than x, for x equal to (m+1)/2, and in accordance with: ##STR6## for odd N, where n=2M(m-f)+2m 1 -1 and x is the greatest integer less than x, for x equal to any of 1+m/2 and m/2, where J is an inverse identity matrix, matrix V is defined to be: ##EQU67## wherein each impulse response, h(n), is found to minimize the stopband error: ##EQU68## where P is a real, symmetric and positive definite matrix, with the elements, using a notation P k,l for denoting a (k,l) th element of matrix P, ##EQU69## where N is even, and where P is a real, symmetric and positive definite matrix, with the elements ##EQU70## where N is odd, and where K is the number of stopbands of H(e j ω), β i are their relative weights and ω i ,1 and ω i ,2 are the bandedges of these stopbands, and ##EQU71## and wherein the filter H k (z) is optimized by finding a least squares optimization h opt such that: ##EQU72##
Claim 2
Dependent Claim 510 words
A pseudo-quadrature-mirror filter bank for near-perfect-reconstruction pseudo-quadrature-mirror filtering of an input signal, comprising: a plurality of analysis filters, each of said plurality of analysis filters including, a first delay chain, operatively coupled to the input signal, forming a set of 2M parallel paths for buffering the input signal; a first cascade of 2M polyphase components of an impulse transfer function H(z) of impulse response, h(n), operatively coupled to said first delay chain; and means, operatively coupled to said first cascade of 2M polyphase components, for generating a 2M-point Discrete Fourier Transform (DFT) to implement a 2M-point Discrete Cosine Transform (DCT) of the input signal, with each analysis filter having an impulse response, h k (n), of a k th analysis filter, where M is the number of subband signals, obtained by cosine-modulating an impulse response, h(n), of a prototype filter with linear phase, according to: ##EQU65## and N is the length of the impulse response, h(n), of the prototype filter; a plurality of synthesis filters, each of said plurality of synthesis filters including, a second delay chain, operatively coupled to the input signal, forming a set of 2M parallel paths for buffering the input signal; a second cascade of 2M polyphase components of an impulse transfer function H(z) of impulse response, h(n), operatively coupled to said second delay chain; and means, operatively coupled to said second cascade of 2M polyphase components, for generating a 2M-point Discrete Fourier Transform (DFT) to implement a 2-point Discrete Cosine Transform (DCT) of the input signal, with each of said plurality of synthesis filters operatively coupled to a respective one of said plurality of analysis filters, each synthesis filter having an impulse response, f k (n), of a k th synthesis filter, obtained by cosine-modulating the impulse response, h(n), of the prototype filter according to: ##EQU66## and N is the length of the impulse response, h(n), of the prototype filter; and wherein each impulse response, h(n), is found in accordance with: ##STR5## for even N, where n=2M(m-l)+2m 1 -1 and x is the greatest integer less than x, for x equal to (m+1)/2, and in accordance with: ##STR6## for odd N, where n=2M(m-f)+2m 1 -1 and x is the greatest integer less than x, for x equal to any of 1+m/2 and m/2, where J is an inverse identity matrix, matrix V is defined to be: ##EQU67## wherein each impulse response, h(n), is found to minimize the stopband error: ##EQU68## where P is a real, symmetric and positive definite matrix, with the elements, using a notation P k,l for denoting a (k,l) th element of matrix P, ##EQU69## where N is even, and where P is a real, symmetric and positive definite matrix, with the elements ##EQU70## where N is odd, and where K is the number of stopbands of H(e j ω), β i are their relative weights and ω i ,1 and ω i ,2 are the bandedges of these stopbands, and ##EQU71## and wherein the filter H k (z) is optimized by finding a least squares optimization h opt such that: ##EQU72##
Claim 3
Dependent Claim 42 words
The pseudo-quadrature-mirror filter bank set forth in claim 1, wherein the impulse response, h(n), provides plurality of analysis filters and the plurality of synthesis filters with a stopband attenuation less than -100 dB and with a reconstruction error less than -100 dB.
Claim 4
Dependent Claim 26 words
The pseudo-quadrature-mirror filter bank set forth in claim 1, wherein the stopband error h t Ph is minimized by subroutine DNOONF of the IMSL Math Library.
Claim 5
Dependent Claim 319 words
A pseudo-quadrature-mirror-filter bank for near-perfect-reconstruction pseudo-quadrature-mirror filtering an input signal, constructed by a process comprising the steps of: finding an impulse response, h(n), of a prototype filter in accordance with: ##EQU73## for even N, where n=2M(m-l)+2m 1 -1 and x is the greatest integer less than x, for x equal to (m+1)/2, and in accordance with: ##EQU74## for odd N, where n=2M(m-l)+2m 1 and x is the greatest integer less than x, for x equal to any of 1+m/2 and m/2, where J is an inverse identity matrix, matrix V is defined to be: ##EQU75## finding the impulse response, h(n), minimizing the stopband error: ##EQU76## wherein P is a real, symmetric, and positive definite matrix, with the elements, using a notation P k,l for denoting a (k,l) th element of matrix P, ##EQU77## where N is even, and wherein P is a real, symmetric and positive definite matrix, with the elements ##EQU78## where N is odd, with K is the number of stopbands of H(e j ω), β i are their relative weights and ω i ,1 and ω i ,2 are the bandedges of these stopbands, and ##EQU79## optimizing, by least squares optimization, H k (z) such that: ##EQU80## generating a plurality of analysis filters from the impulse response, h(n), each analysis filters having an impulse response, h k (n), of a k th analysis filter generated by cosine-modulating the impulse response, h(n), with linear phase, according to: ##EQU81## and N is the length of the impulse response, h(n), of the prototype filter; generating a plurality of synthesis filters, from the impulse response, h(n), each synthesis filter having an impulse response, f k (n), of a k th synthesis filter, of the prototype filter according to: ##EQU82## and N is the length of the impulse response, h(n), of the prototype filter; and coupling the output of each k th analysis filter with each k th synthesis filter, respectively.
Claim 6
Dependent Claim 319 words
A pseudo-quadrature-mirror-filter bank for near-perfect-reconstruction pseudo-quadrature-mirror filtering an input signal, constructed by a process comprising the steps of: finding an impulse response, h(n), of a prototype filter in accordance with: ##EQU73## for even N, where n=2M(m-l)+2m 1 -1 and x is the greatest integer less than x, for x equal to (m+1)/2, and in accordance with: ##EQU74## for odd N, where n=2M(m-l)+2m 1 and x is the greatest integer less than x, for x equal to any of 1+m/2 and m/2, where J is an inverse identity matrix, matrix V is defined to be: ##EQU75## finding the impulse response, h(n), minimizing the stopband error: ##EQU76## wherein P is a real, symmetric, and positive definite matrix, with the elements, using a notation P k,l for denoting a (k,l) th element of matrix P, ##EQU77## where N is even, and wherein P is a real, symmetric and positive definite matrix, with the elements ##EQU78## where N is odd, with K is the number of stopbands of H(e j ω), β i are their relative weights and ω i ,1 and ω i ,2 are the bandedges of these stopbands, and ##EQU79## optimizing, by least squares optimization, H k (z) such that: ##EQU80## generating a plurality of analysis filters from the impulse response, h(n), each analysis filters having an impulse response, h k (n), of a k th analysis filter generated by cosine-modulating the impulse response, h(n), with linear phase, according to: ##EQU81## and N is the length of the impulse response, h(n), of the prototype filter; generating a plurality of synthesis filters, from the impulse response, h(n), each synthesis filter having an impulse response, f k (n), of a k th synthesis filter, of the prototype filter according to: ##EQU82## and N is the length of the impulse response, h(n), of the prototype filter; and coupling the output of each k th analysis filter with each k th synthesis filter, respectively.
Claim 7
Dependent Claim 49 words
The pseudo-quadrature-mirror-filter bank constructed by the process set forth in claim 4, wherein the step of finding the impulse response, h(n), minimizing the stopband error: ##EQU83## includes the steps of: computing the stopband error h t Ph and the gradient of h t Ph; and minimizing h t Ph.
Claim 8
Dependent Claim 37 words
The pseudo-quadrature-mirror-filter bank constructed by the process set forth in claim 4, wherein the step of optimizing by least squares optimization includes the steps of: linearizing a set of quadratic constraints; and minimizing a cost function Φ.
Claim 9
Dependent Claim 62 words
The Pseudo-quadrature-mirror-filter bank constructed by the process set forth in claim 4, further comprising the steps of: buffering, using a delay chain, the input signal to form a set of 2M parallel paths; cascading the buffered input signal using a cascade of 2M polyphase components of H(z); and implementing a 2M-point Discrete Cosine Transform (DCT) using a 2M-point Discrete Fourier Transform (DFT).
Claim 10
Independent Claim 318 words
A method, using a pseudo-quadrature-mirror-filter bank, for near-perfect-reconstruction pseudo-quadrature-mirror filtering an input signal, comprising the steps of: finding an impulse response, h(n), of a prototype filter in accordance with: ##STR7## for even N, where n=2M(m-l)+2m 1 -1 and x is the greatest integer less than x, for x equal to (m+1)/2, and in accordance with: ##STR8## for odd N, where n=2M(m-l)+2m 1 and x is the greatest integer less than x, for x equal to any of 1+m/2 and m/2, where J is an inverse identity matrix, matrix V is defined to be: and ##EQU84## finding the impulse response, h(n), minimizing the stopband error: ##EQU85## wherein P is a real, symmetric and positive definite matrix, with the elements, using a notation P k,l for denoting a (k,l) th element of matrix P, ##EQU86## where N is even, and wherein P is a real, symmetric and positive definite matrix, with the elements ##EQU87## where N is odd, with K is the number of stopbands of H(e j ω), β i are their relative weights and ω i ,1 and ω i ,2 are the bandedges of these stopbands, and optimizing, by least squares optimization, H k (z) such that: ##EQU88## generating a plurality of analysis filters from the impulse response, h(n), each analysis filters having an impulse response, h k (n), of a k th analysis filter generated by cosine-modulating the impulse response, h(n), with linear phase, according to: ##EQU89## and N is the length of the impulse response, h(n), of the prototype filter; generating a plurality of synthesis filters, from the impulse response, h(n), each synthesis filter having an impulse response, f k (n), of a k th synthesis filter, of the prototype filter according to: ##EQU90## and N is the length of the impulse response, h(n), of the prototype filter; and coupling the output of each k th analysis filter with each k th synthesis filter, respectively.
Claim 11
Independent Claim 318 words
A method, using a pseudo-quadrature-mirror-filter bank, for near-perfect-reconstruction pseudo-quadrature-mirror filtering an input signal, comprising the steps of: finding an impulse response, h(n), of a prototype filter in accordance with: ##STR7## for even N, where n=2M(m-l)+2m 1 -1 and x is the greatest integer less than x, for x equal to (m+1)/2, and in accordance with: ##STR8## for odd N, where n=2M(m-l)+2m 1 and x is the greatest integer less than x, for x equal to any of 1+m/2 and m/2, where J is an inverse identity matrix, matrix V is defined to be: and ##EQU84## finding the impulse response, h(n), minimizing the stopband error: ##EQU85## wherein P is a real, symmetric and positive definite matrix, with the elements, using a notation P k,l for denoting a (k,l) th element of matrix P, ##EQU86## where N is even, and wherein P is a real, symmetric and positive definite matrix, with the elements ##EQU87## where N is odd, with K is the number of stopbands of H(e j ω), β i are their relative weights and ω i ,1 and ω i ,2 are the bandedges of these stopbands, and optimizing, by least squares optimization, H k (z) such that: ##EQU88## generating a plurality of analysis filters from the impulse response, h(n), each analysis filters having an impulse response, h k (n), of a k th analysis filter generated by cosine-modulating the impulse response, h(n), with linear phase, according to: ##EQU89## and N is the length of the impulse response, h(n), of the prototype filter; generating a plurality of synthesis filters, from the impulse response, h(n), each synthesis filter having an impulse response, f k (n), of a k th synthesis filter, of the prototype filter according to: ##EQU90## and N is the length of the impulse response, h(n), of the prototype filter; and coupling the output of each k th analysis filter with each k th synthesis filter, respectively.
Claim 12
Dependent Claim 44 words
The method set forth in claim 8, wherein the step of finding the impulse response, h(n), minimizing the stopband error: ##EQU91## includes the steps of: computing the stopband error h t Ph and the gradient of h t Ph; and minimizing h t Ph.
Claim 13
Dependent Claim 32 words
The method set forth in claim 8, wherein the step of optimizing by least squares optimization includes the steps of: linearizing a set of quadratic constraints; and minimizing a cost function Φ.
Claim 14
Dependent Claim 57 words
The method set forth in claim 8, further comprising the steps of: buffering, using a delay chain, the input signal to form a set of 2M parallel paths; cascading the buffered input signal using a cascade of 2M polyphase components of H(z); and implementing a 2M-point Discrete Cosine Transform (DCT) using a 2M-point Discrete Fourier Transform (DFT).
Claim 15
Independent Claim 159 words
A method, using a pseudo-quadrature-mirror-filter bank, for near-perfect-reconstruction pseudo-quadrature-mirror filtering an input signal, comprising the steps of: finding an impulse response, h(n), of a prototype filter minimizing the stopband error: ##EQU92## optimizing the impulse response, h(n), of the prototype filter; generating a plurality of analysis filters from the impulse response, h(n), each analysis filters having an impulse response, h k (n), of a k th analysis filter generated by cosine-modulating the impulse response, h(n), with linear phase, according to: ##EQU93## and N is the length of the impulse response, h(n), of the prototype filter; generating a plurality of synthesis filters, from the impulse response, h(n), each synthesis filter having an impulse response, f k (n), of a k th synthesis filter, of the prototype filter according to: ##EQU94## and N is the length of the impulse response, h(n), of the prototype filter; and coupling the output of each k th analysis filter with each k th synthesis filter, respectively.
Claim 16
Independent Claim 159 words
A method, using a pseudo-quadrature-mirror-filter bank, for near-perfect-reconstruction pseudo-quadrature-mirror filtering an input signal, comprising the steps of: finding an impulse response, h(n), of a prototype filter minimizing the stopband error: ##EQU92## optimizing the impulse response, h(n), of the prototype filter; generating a plurality of analysis filters from the impulse response, h(n), each analysis filters having an impulse response, h k (n), of a k th analysis filter generated by cosine-modulating the impulse response, h(n), with linear phase, according to: ##EQU93## and N is the length of the impulse response, h(n), of the prototype filter; generating a plurality of synthesis filters, from the impulse response, h(n), each synthesis filter having an impulse response, f k (n), of a k th synthesis filter, of the prototype filter according to: ##EQU94## and N is the length of the impulse response, h(n), of the prototype filter; and coupling the output of each k th analysis filter with each k th synthesis filter, respectively.

References Cited (Backward Citations)

US5214678A Cited
US5109419A Cited

USPTO Annuity & Maintenance Windows

Annuity Stage Window Opens Payment Deadline Surcharge Period Until Estimated Fee Status
3.5 Years (Year 4) 1998-07-24 1999-01-23 1999-07-25 $2000.00 Window Closed / Lapsed
7.5 Years (Year 8) 2002-07-24 2003-01-23 2003-07-25 $3760.00 Window Closed / Lapsed
11.5 Years (Year 12) 2006-07-24 2007-01-23 2007-07-25 $7700.00 Window Closed / Lapsed

Assigned Categories & Tags

No categories assigned to this patent yet.

Assign Category Tag