"Hima" <himabindu81@gmail.com> wrote in message
news:V9ydna-3xs4YuknenZ2dnUVZ_tadnZ2d@giganews.com...

> Hi
> I need to design Smith-Barnwell filters for N=4,8,16,32,64 in Matlab.
> I need to take these steps
>
> (1) design an equiripple half-band filter using the Remez exchange
>
> (2) make the filter frequency response positive-definite
>
> (3) spectral factorization
>
> I am new to filter design , can anyone guide as to where do I start .
> It would be great.
>
>
> Where can I get remex exchange algorithm to understand and read...
>

fred harris at UCSD has a method for using the Parks-McClellan program
directly to design a half-band filter.
1) Design a filter that is one-half the frequency response as follows:
Passband 0.,.4
Stopband, 0.5,0.5 note the stopband is zero width - but that's OK.
2) Put zero-valued coefficients between each of the coefficients found.
This doubles the sample rate of the existing filter so that the frequency
response repeats twice (and, it happens to mirror image in sign).
3) Multiply the coefficients by 0.5 for scaling
4) Replace the center coefficient (zero) with 0.5.
This shifts the frequency response upward by 0.5 so that it is all positive.
I have the P-M program .exe posted at:
ftp.mission-systems-inc.com/outgoing/Parks_McClellan/ParksMcClellan.exe
And a half-band program specifically at:
ftp://ftp.mission-systems-inc.com/outgoing/Halfband/
You might look up Hermann & Schussler's paper on the rest.
Here's something you might try:
Design the half-band filter.
Once you get the design, note the stopband peak ripple value.
Now, redesign the filter using the stopband peak ripple value for the
stopband criterion instead of zero.
Another method would be to add the stopband peak ripple value to the center
coefficient - thus shifting the filter upward by that amount. This makes
the negative stopband ripples peaks into double zeros and the entire filter
frequency response is then positive ... and I don't know about the
"definite" part....
Hermann & Schussler would then take the square root of the filter response -
a type of factorization. There are many choices here including minimum
phase, maximum phase, or a mixture that's somewhere in between depending on
which zero pairs you choose and whether they are inside or outside the unit
circle.
Smith Barnwell may be different.... but these are some nice ideas for filter
design.
Fred

Reply by mlimber●January 23, 20062006-01-23

Hima wrote:

> I need to design Smith-Barnwell filters for N=4,8,16,32,64 in Matlab.

[snip]
You might consider posting in a MATLAB group to get more specific help
than you've gotten here so far, but you should do a little more
searching on your own first.

> Where can I get remex exchange algorithm to understand and read...

Hi
I need to design Smith-Barnwell filters for N=4,8,16,32,64 in Matlab.
I need to take these steps
(1) design an equiripple half-band filter using the Remez exchange
(2) make the filter frequency response positive-definite
(3) spectral factorization
I am new to filter design , can anyone guide as to where do I start .
It would be great.
Where can I get remex exchange algorithm to understand and read...
Thankyou
Hima