Filter Design for Orthogonal Two-Channel
Filter Bank - Documentation
Go back to previous
page.
A program for the design of a minimum phase lowpass FIR
analysis filter for a two-channel orthogonal filter bank
(wavelet system).
The filter has a specified ripple size and a specified
degree of flatness.
The filters produced by this program achieve a trade-off
between the classical Daubechies (Herrmann) filter and an
equi-ripple solution.
This permits a trade-off between frequency selectivity and
the number of vanishing moments.
The Daubechies (maximally flat) filter is obtained as a
special case by using L/2==N below.
program name : fir_orthog.m
subprogram : oref.m, local_max.m, choose.m
________________________________________________________
help on fir_orthog _____________________________________
[h,h2,h2] = fir_orthog(L,N,del)
Design of a minimum phase lowpass filter for a two
channel orthogonal filter bank (wavelet system).
by : Ivan Selesnick, Rice University.
input
L : length of filter (must be even)
N : degree of flatness
del : ripple size (in magnitude squared)
need L/2-N even and L/2 >= N
output
h : minimum phase wavelet filter
h = conv(h2,h2); h2 contains the roots at z=-1,
h2 contains all the other roots.
subprograms:
oref.m, local_max.m, choose.m
% EXAMPLE
L = 14;
N = 5;
del = 0.01;
[h,h2,h2] = fir_orthog(L,N,del);
________________________________________________________
Please report any bugs or send comments regarding the programs
to selesi@ece.rice.edu