1.2 Script: clc;close all;clear all; N=81;n=0:1:N-1; Fc=2000; Fs=10000; wc=2*Fc/Fs; h=fir1(N-1,wc,hamming(N)); figure; [H,f]=freqz(h,1,1024,Fs); subplot(211) stem(n,h);grid;xlabel('n');ylabel('h(n)') title('Impulse response using Hamming window') subplot(212) plot(f,abs(H));grid;xlabel('f(Hz)') ylabel('|H(f)|') title('Magnitude response') Impulse response using Hamming window 0.4 0.3 0.2 h(n) 0.1 08 0 10 20 30 40 50 60 70 80 Magnitude response 1 0.8 0.6 14H 0.4 2.1) clc;close all;clear all; Ap=0.5;As=40; fp=4000;fs=8000; Fs=40000; wp=2*fp/Fs;ws=2*fs/Fs; [N wc]=buttord(wp,ws,Ap,As); [b a]=butter(N,wc); figure; freqz(b,a,1024,Fs) -100 Magnitude (dB) -200 -300 -400 0 0.5 1.5 1 Frequency (Hz) 2 x104 -200 Phase (degrees) -400 -600 -800 0 0.5 1.5 2 1 Freque Ap=0.5;As=40; fp=4000;fs=8000; Fs=40000; wp=2*fp/Fs;ws=2*fs/Fs; [N wp]=cheb1ord(wp,ws,Ap,As); [b a]=cheby1(N,Ap,wp); figure; freqz(b,a,1024,Fs) -100 Magnitude (dB) -200 -300 -400 0 0.5 1.5 2 1 Frequency (Hz) *104 0 -100 -200 Phase (degrees) -300 -400 -500 0 0.5 1.5 2 1 Ap=0.5;As=40; fp=4000;fs=8000; Fs=40000; wp=2*fp/Fs;ws=2*fs/Fs; [N ws]=cheb2ord(wp,ws,Ap,As); [b a]=cheby2(N,As,ws); figure; freqz(b,a,1024,Fs) 0 -20 -40 Magnitude (dB) -60 -80 -100 0 0.5 1.5 2 1 Frequency (Hz) *104 -100 Phase (degrees) -200 -300 -400 0 0.5 1.5 2 1 Fre Ap=0.5;As=40; fp=4000;fs=8000; Fs=40000; wp=2*fp/Fs;ws=2*fs/Fs; [N wp]=ellipord(wp,ws,Ap,As); [b a]=ellip(N,Ap,As,wp); figure; freqz(b,a,1024,Fs)