Filter Koeffizienten auf XMega rechnen

[Tobias B. (horschtx)]

Hallo Zusammen,

ich versuche gerade Filter Koeffizienten für den DSP ADAU1452 von Analog 
zu berechnen.

Als Controller habe ich einen XMega128A1. Aber ich denke der Typ sollte 
doch relativ egal sein.

Das sind die Formel wie sie Analog vorgibt:

ω0 = 2*pi*f0/Fs
alpha = sin(ω0) / 2.0 * 1/sqrt(2)

a0 =   1 + alpha
a1 =  -2*cos(ω0)
a2 =   1 - alpha
b0 = -(1 + cos(ω0)) * gainLinear / 2
b1 = -(1 + cos(ω0)) * gainLinear
b2 = -(1 + cos(ω0)) * gainLinear / 2

Ich wollte jetzt erstmal den Koeffizienten a2 rechnen, der scheinte mir 
der einfachste zu sein.

Hier ist mal mein Codeausschnitt den ich bis jetzt so habe:

  double alpha,w0,a2;
  int32_t param227,param223;
  
  w0 = (2*M_PI*filter_hp_top_berechnen)/192000;      //W0 berechnen
  alpha = sin(w0*(M_PI/180))/2.0*(1/sqrt(2));        //alpha berechnen
  a2=1-alpha;                        //a2 berechnen
    
  param223 = a2 * 16777216;                //auf das 8.24 Format bringen
  param227 = param223 + 2147483648;
  
  werte[3] = (uint8_t)param227;              //Auf die bytes verteilen
  werte[2] = (uint8_t)(param227>>8);
  werte[1] = (uint8_t)(param227>>16);
  werte[0] = (uint8_t)(param227>>24);
  werte[0] = werte[0] ^ 0x80;                //Vorzeichen Anpassen

Die Aufbereitung der Daten auf das 8.24 Format funktioniert. Das hab ich 
mal mit festen Werten getestet.

Aber wenn ich jetzt den Controller die Koeffizienten Rechnen lasse 
scheint es zu ungenau zu sein. Die unteren 2 Bytes unterscheiden sich 
doch sehr von den Rechnung aus dem Sigma Studio.

Gibt es eine Möglichket die Genauigkeit zu erhöhen?
Oder hab ich für die Anwenung einfach den falschen Controller genommen.

Die Funktion mit der Filterberechnung ist jetzt nicht ultra wichtig, 
wäre nur ein nice to have. Momentan nehme ich noch fest vorgegebene 
Werte für die Filter.

Schon einmal Danke im voraus

[Karl M. (Gast)]

Naja welchen Compiler verwendest du ?

bei AVR Gcc ist double = float (32Bit)

[Tobias B. (horschtx)]

jo ist AVR Gcc

[Wolfgang (Gast)]

Tobias B. schrieb:
> Aber wenn ich jetzt den Controller die Koeffizienten Rechnen lasse
> scheint es zu ungenau zu sein.

Dann solltest du dir erstmal überlegen, welche Genauigkeit der 
Koeffizienten erforderlich ist, ohne dass deine Filtercharakteristik zu 
sehr abweicht.

Tobias B. schrieb:
> jo ist AVR Gcc

Da ist dem Compiler das "double" ziemlich egal.