exp7

clc;
clear all;
%DEFINING PARAMETERES
wp =0.65 * pi;
ws =0.5 * pi;

% tr_width = abs(ws-wp);
% wc = (wp+ws)/2;
%
% %DEFINING WINDOW
% M = ceil(6.6*pi/tr_width);
% wn=(hamming(M));
% n = -(M-1)/2:(M-1)/2;
%
% fc = wc/(2*pi);
% hd = -2*fc*(sinc(2*fc*n));
%
% h=hd.*wn';
% [HW,WW] = freqz(h,1);
%
% subplot(2,1,1);
% stem(n,wn);
% title('Hamming Window 101915122');
% xlabel('n-->');
% ylabel('W[n]-->');
%
% subplot(2,1,2);
% plot(WW./pi,abs(HW));
% title('Low Pass Filter 101915122');
% xlabel('Normalised Freq-->');
% ylabel('Magnitude-->')

% ______________________________________________________________
%
% 7b
%
% clc;
% clear all;
% %DEFINING PARAMETERES
% % wp =input('Enter pass band frequency: ') * pi;
% % ws =input('Enter stop band frequency: ') * pi;
%
tr_width =abs(ws-wp);
wc = (wp+ws)/2;

%DEFINING WINDOW
M = ceil(6.2*pi/tr_width);
wn=(hanning(M));
n = -(M-1)/2:(M-1)/2;

fc = wc/(2*pi);
hd = -2*fc*(sinc(2*fc*n));

h=hd.*wn';
[HW,WW] = freqz(h,1);

subplot(2,1,1);
stem(n,wn);
title('Hanning Window 101915122');
xlabel('n-->');
ylabel('W[n]-->');

subplot(2,1,2);
plot(WW./pi,abs(HW));
title('Low Pass Filter 101915122');
xlabel('Normalised Freq-->');
ylabel('Magnitude-->')