IIR unterschiedliche Koeffizienten +AVR

Mit Scilab (www.scilab.org) und dem Befehl eqiir lassen sich auch 
IIR-Filter berechnen, sie werden schon in Teilfilter 2.Ordnung zerlegt. 
Ist das die alte oder neuere Appnote? in der alten war auch noch ein IIR 
6.Ordnung beschrieben, das schwingt sehr zuverlässig, außerdem ein paar 
Schreibfehler im Programm.

Ein Beispiel für ein elliptisches Stopbandfilter, das ich mit dem 
ATMega48 gebaut habe:

// Stopband-Filter für 1750 Hz, Samplerate =18029Hz
// Eckfrequenzen 1400Hz, 1700Hz, 1800 Hz, 2100Hz:
omega=[2*%pi*(1400/18029),2*%pi*(1700/18029),2*%pi*(1800/18029),2*%pi*(2 
100/18029)];
// Maximale Durchgangswelligkeit 0,5 dB
deltapass = (1-10**(-0.5/20)) ;
// minimale Sperrband-Dämpfung 50 dB
deltastop = 10**(-50/20);
// IIR-Filter berechnen:
[sos,gain,zeroes,poles] = eqiir('sb','el',omega,deltapass,deltastop);
// Second order sections anzeigen:
sos
gain

                            2                             2 
2
          1 - 1.6219938z + z            1 - 1.6396794z + z           1 - 
1.656616z + z
sos = --------------------------    ------------------------   
-------------------------
                                2                             2 
2
      0.9383702 - 1.5100165z + z     0.8466245 - 1.513936z + z 
0.9485511 - 1.663563z + z

gain = 0.8683583

die Brüche müssen noch in Zähler und Nenner durch z^2 dividiert werden, 
damit die übliche Schreibweise mit z^(-2) und z^(-1) herauskommt.
( hoffentlich passt das in eine Zeile, sonst wird es unleserlich).

Ich hatte es befürchtet, in Scilab kommt es lesbarer heraus, je nach 
Fensterbreite.

http://www.mikrocontroller.net/articles/Digitalfilter_mit_ATmega
Ich habe mal endlich meinen Artikel zum Thema ins Wiki eingetragen. Noch 
nicht komplett, aber das wesentliche ist gesagt. Kommentare, am besten 
unter der dortigen Rubrik "Diskussion" sehr erwünscht.