FPGA Central - World's 1st FPGA / CPLD Portal

FPGA Central

World's 1st FPGA Portal

 

Go Back   FPGA Groups > NewsGroup > DSP

DSP comp.dsp newsgroup, mailing list

Reply
 
LinkBack Thread Tools Display Modes
  #1 (permalink)  
Old 04-03-2007, 03:40 AM
robert bristow-johnson
Guest
 
Posts: n/a
Default Re: Cascade form IIR Filters

On Apr 2, 7:31 pm, "robert bristow-johnson"
<r...@audioimagination.com> wrote:
>
> we were just talking about this at the music-dsp mailing list. do you
> want explicit closed form formulae for biquad coefs for Butterworth.
> they exist. if you want i can get my email copies (and correct them,
> 'cuz i did copulate 'em up a bit at first, but i think we fixed it).
>


use this in reference to the Audio EQ cookbook that presently (and
forever, i hope) lives at:

http://www.musicdsp.org/files/Audio-EQ-Cookbook.txt

for an Nth order Butterworth (N even or odd) you will have N/2 biquad
sections ((N-1)/2 for odd N or floor(N/2) for either even or odd N),
each will have the same resonant frequency w0 and will have Q:

Q = 1/( 2*sin((pi/N)*(n + 1/2)) )

where 0 <= n < (N-1)/2

then, use the cookbook LPF or HPF to get the coefs for each biquad
section.

if it's an odd order Butterworth, you'll need an additional 1st order
section:

H(z) = (b0 + b1*z^-1)/(a0 + a1*z^-1)


1st order LPF: H(s) = 1/(s+1)

b0 = sin(w0)
b1 = sin(w0)
a0 = cos(w0) + sin(w0) + 1
a1 = sin(w0) - cos(w0) - 1


1st order HPF: H(s) = s/(s+1)

b0 = cos(w0) + 1
b1 = -(cos(w0) + 1)
a0 = cos(w0) + sin(w0) + 1
a1 = sin(w0) - cos(w0) - 1


credits to Peter Schoffhauzer <[email protected]> for
catching and fixing errors.

r b-j


Reply With Quote
Reply

Bookmarks

Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Re: Cascade form IIR Filters Jerry Avins DSP 0 04-03-2007 02:34 AM
The Cascade Algorithm - Recursion, Wavelets moosedude DSP 0 07-25-2006 05:23 PM
Implementing cascade IIR filters kmap DSP 2 06-08-2006 08:26 AM
Cascade Form IIR Filters bngguy DSP 4 04-20-2004 10:59 PM
Convert cascade to parallel Jon Harris DSP 3 03-04-2004 07:58 AM


All times are GMT +1. The time now is 01:56 AM.


Powered by vBulletin® Version 3.8.0
Copyright ©2000 - 2012, Jelsoft Enterprises Ltd.
Search Engine Friendly URLs by vBSEO 3.2.0
Copyright 2008 @ FPGA Central. All rights reserved