% filtro ricorsivo reiezione 50 Hz
%recursive band rejection filter 50 Hz
function [cMA,cAR]=rico(z,B,fc,T);

% input parameters
%	z  (0.01) attenuazione minima alla frequenza fc
%   minimum attenuation at the center frequency frequency band

%	B  (2) larghezza di banda corrispondente alla attenuazione 0.707
%	Bandwidth which corresponds to attenuation 0.707

% 	fc (50) frequenza di centro banda 
%   Central frequency

% 	T  intervallo di campionamento 
%   Sampling interval 

% output parameters
%	cMA  filter coefficients (MA part)
%	cAR  filter coefficients (AR part)


b=pi*B*T;
a=b*z;
c1=-2*(1-a)*cos(2*pi*fc*T);
c2=(1-a)^2;
c3=2*(1-b)*cos(2*pi*fc*T);
c4=-(1-b)^2;
cMA=[1 c1 c2];
cAR=[1 -c3 -c4];

end