clc
clear all
close all
% Ce programme réalise la conception d'un filtre passe-bande 
% par la méthode des fenêtres.
% Définition des paramètres 
N=50; % Nombre des coefficients du filtre 
NB=20000; % Nombre des échantillons de la réponse en fréquence
Fe=20000; % Fréquence d'échantillonnage
Te=1/Fe; % Période d'échantillonnage
fc1=3000; % Fréquence de coupure basse
fc2=7000; % Fréquence de coupure haute
n=0:N-1; % Indices du temps de la réponse impulsionnelle
n1=n-(N-1)/2; % Translation pour rendre la réponse impulsionnelle causale
% réponse impulsionnelle causale
h=(2*fc2*Te)*(sinc(2*fc2*n1*Te))-(2*fc1*Te)*(sinc(2*fc1*n1*Te));
h=h.*(rectwin(length(h)))'; % troncature par une fenêtre
% Ici vous devez changer, chaque fois, la fenêtre de troncature
% Tapez "help rectwin" dans l'espace de travail pour trouver les autres fenêtres
subplot 211
stem(n,h); % Représentation de la réponse impulsionnelle
f=0:(NB/2-1); % Indices des fréquences de la réponse en fréquence
a=[1];
H=freqz(h,a,f,Fe); % Calcul de la réponse en fréquence
subplot 212
plot(f,abs(H))