Pan Tompkins Bandpassfilter

[Hans (Gast)]

Hallo Leute,

ich bin gerade dabei den Bandpassfilter (5-11Hz) des 
Pan/Tompkins-Algorithmus in C zu implementieren. Leider liefert er mir 
nicht das gewünschte Ergebnis, sondern im Ausgangssignal ist immer noch 
eine deutliche 50Hz-Schwingung zu erkennen und überhaupt kein 
EKG-Signal. Allerdings erkenne ich in meiner Implementierung keinen 
Fehler. Außerdem hatte ich den kompletten Algorithmus vorher schon mit 
einem IIR-Butterworth-Bandpass 3.Ordnung am laufen, also kann es auch 
nicht an dem Rest des Codes liegen.
Ich wollte nur den IIR-Bandpass durch den originalen Bandpass aus dem 
Pan-Tompkins-Algorithmus ersetzen um einen linearen Phasengang zu 
erhalten und so das Ergebnis etwas zu verbessern.

Der Bandpass besteht aus einem Hoch- und einem Tiefpass deren Funktionen 
folgendermaßen gegeben sind:

Hochpass:
y(nT) = 32x(nT-16T) - 1/32[y(nT-T) + x(nT) - x(nT-32T)]

Tiefpass:
y(nT) = 2y(nT-T) - y(nT-2T) + x(nT) - 2x(nT-6T) + x(nT-12T)

Habe ich vielleicht doch einen Fehler übersehen?

-----------------Filterfunktionen---------------------
signed long lowpass(signed long *x, signed long *y)
{  
  y[0] = 2*y[1] - y[2] + ((x[0] - 2*x[6] + x[12]));

  return y[0];

  y[1] = y[2];
  y[0] = y[1];
}


signed long highpass(signed long *x, signed long *y)
{  
  y[0] = x[16] - (y[1] >> 5) - (x[0] >> 5) + (x[32] >> 5);

  return y[0];

  y[0] = y[1];
}

-------------------Aufruf in der Main----------------------
x_lowpass[0] = data;
data = lowpass(x_lowpass, y_lowpass);
//alle Werte um eine Stelle weiterschieben
for(l = 11; l > 0; l--)
{
  x_lowpass[l] = x_lowpass[l+1];
}
x_lowpass[0] = x_lowpass[1];

x_highpass[0] = data;
data = highpass(x_highpass, y_highpass);
//alle Werte um eine Stelle weiterschieben
for(l = 31; l > 0; l--)
{
  x_highpass[l] = x_highpass[l+1];
}
x_highpass[0] = x_highpass[1];

[Paul B. (Gast)]

Woher stammen denn Deine Eingangsdaten? Ich kenne mich mit EKG nicht so 
aus, aber 50Hz klingen mir verdächtig nach ungewollten Störungen. Sind 
die empfangen worden oder simuliert?