Notch filter

Notch filter for power line noise rejection:



Matlab code:

fs=1000;
fn=50;
r=0.99;
c=cos(2*pi*fn/fs);
bn=[1 -2*c 1];
an=[1 -2*r*c r*r];
y=filter(bn,an,x);