exp6

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(1.8*pi/tr_width);
wn=(rectwin(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('Rectangular 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-->')

% ______________________________________________________________

% 6b
%
% clc;
% clear all;
% wp =input('Enter pass band frequency: ') * pi;
% ws =input('Enter stop band frequency: ') * pi;
% tr_width =abs(ws-wp);
% wc = (wp+ws)/2;
% M = ceil(6.1*pi/tr_width);
% wn=(triang(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('Triangular 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-->');