function [y,B,A]=filt50hz(x,tau);
% Фильтр-пробка на 50 Гц - вырезаетс? спектр с шириной 1 Гц
%
%    [y,[B,A]]=filt50hz(x,tau);

df=0.5;     % полуширина вырезаемой полосы

f0=1/tau;   % Частота опроса
n=1;		% номер гармоники сетевой помехи
f50=abs(n*f0-n*50); % Частота просачивани? помехи 50 Гц
[B,A]=cheby1(7,0.01,[f50-df f50+df]/(0.5*f0),'stop');

y=filtfilt(B,A,x);
