Untitled

b=[0.0048,0.1224,0.0408];
a=[1.00,-1.2978,0.7875,-0.1632];
Fs=2000;
figure(1)
freqz(b,a,512,fs)
title('Magnitude and phase response of IRR filter');
f1=100;f2=900;
ts=1/Fs;
l=2000;
n=0:1-1;
y=sin(2*pi*fi/Fs*n)+0.1 sin(2*pi*f2/Fs*n);
figure(2)
subplot(2,1,1);
plot(n(1:100)*ts*1000,y(1:100));
title('sine wave with noise ');
xlabel('Time in ms ');