Einfacher IIR-Filter in C implementieren

[Thomas (Gast)]

Hallo zusammen,

im Rahmen meines Studiums soll ich einen einfachen IIR-Filter mit C-Code 
implementieren. Leider hat mir die Suchfunktion nicht wirklich 
weitergeholfen.
Der Filter hat folgende Differenzengleichung:

y(k) = c * {x(k) – x(k-2)} – y (k-1)

Mein C-Code würde folgendermaßen aussehen:

void IIR(float *x, float *y)
{
  int k;
  
  y[0] = c * x[0];

  for (k=1; k < NUM_SAMPLES; k++)
  {
    if(k < 2)
    {
      y[k] = c * x[k] - y [k-1];
    }
    else
    {
      y[k] = c * (x[k] - x[k-2]) - y[k-1];
    }
  }
}

Ist die Implementierung so möglich?

Ich habe auch schon mal was von der Verwendung eines Ringpuffers 
gelesen.  Wie würde hier die Implementierung aussehen?

Vielen Dank
Thomas

[Matthias (Gast)]

In deinem Fall berechnest du offline die gefilterten Daten.
Wenn du einen Filter für online-Verwendung brauchst, könnte man zb 
folgendes schreiben:

#define NTAPS 3

// Buffer
float X[NTAPS];
float Y[NTAPS];

// Koeffizienten
float a[]={c,0,-c};
float b[]={0,1,0};

float IIR(float* x, float* y, float* a, float* b)
{
  int k;

    for(k = 1; y[0] = a[0]*x[0]; k < NTAPS; k++)
    {
        y[0] = y[0] + a[k]*x[k] - b[k]*y[k];
    }

    for(k = NTAPS-1; k > 0; k--) // verschiebe Werte um eins nach hinten
    {
       y[k] = y[k-1];
       x[k] = x[k-1];
    }

    return y[0]; // Filter-Ausgang
}

Das wäre der allgemeine Fall und ist nicht als Ringbuffer gebaut.
y[0] stellt neuesten Wert dar. X[0] kriegt immer den neuen Messwert 
zugewiesen, usw.

[Besucher (Gast)]

Google mkfilter

der erste Treffer. Dort findest Du alles, Sourcecode für Filter sowie 
Berechnung der Koeffizienten

[Gebhard R. (Firma: Raich Gerätebau & Entwicklung) (geb)]

Ich hab mal sowas hier hereingestellt:

Beitrag "IIR Filter in Integer"

Grüsse