在MatLab中实现50 Hz的陷波滤波器

Question

我试图在MatLab中实现50 Hz的陷波滤波器。这是代码：

     subplot(2,1,2);
     %Fse = 220;
     time = 0:1/f:sB-1/f;
     x = buf;
     wo = 60 / (1000/2);
     bw = wo / 50;
     [b,a] = iirnotch(wo,bw);
     y = filter(b,a,x);
     ydft = fft(y);
     xdft = fft(x);
     plot(x,y);

buf的值为buf = zeros([fse*sB,numel(name)]);，采样频率为220 Hz 而且我没有得到预期的结果。你能告诉我我哪里错了吗？ 目前的输出如下：

所需的输出是这样的：

我正在谈论第二张图。第一张图很好。