1234567891011121314151617181920212223242526272829303132 |
- % 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
|