Noms et prénoms
UNIVERSITE DJILLALI LIABES FACULTE DE GENIE ELECTRIQUE …………………………………………………..
DEPARTEMENT DE TELECOMMUNICATION
LICENCE EN TELECOMMUNICATION
TRAVAUX PRATIQUE TRAITEMENT DU SIGNAL …………………………………………………..
…………………………………………………..
TP N° 5
ANALYSE ET SYNTHESE DES FILTRES
ANALOGIQUE
Objectif : Le but de ce TP est d’effectuer une synthèse et analyse de filtres analogiques de types
butterworth
Connaissances requises : Chapitre 3(Filtrage analogique).
Matériels :
- Ordinateur de bureau.
- Système d’exploitation Windows.
- Logiciel MATLAB.
Travail demandé
A) synthèse d'un filtre passe-bas de butterworth
On se donne les objectifs suivants:
type de filtre : passe-bas
bande passante allant de 0 à fp=200 Hz
atténuation maximale tolérée Amax=2 dB
bande affaiblie ou atténuée allant de fa=350 à l'infini
atténuation minimale exigée Amin=40 dB
1) Tracer le gabarit du filtre sur la figure suivante:
50 dB
figure 1 500 Hz
2) En utilisant la commande "buttord", on calcule l'ordre nécessaire pour satisfaire ce gabarit:
fp=200 ;
fa=350 ;
%pour le filtre analogique sur MATLAB on introduit les pulsations
wp=2*pi*fp;
wa=2*pi*fa;
rp=2; %attenuation maximale en bande passante 2dB
ra=40; %attenuation minimale en bande coupée 40dB
[n, wn] = buttord(wp, wa, rp, ra, 's');
%n ordre du filtre necessaire pouvant satisfaire le gabarit exigé
-1-
%wn frequence de coupure naturelle à -3dB en rad/s
3) En utilisant la commande "butter" concevoir le filtre passe-bas d'ordre n et de pulsation de
coupure wn, ensuite tracer la fonction de transfert et l'atténuation en fonction de la fréquence.
[b,a]=butter(n,wn,'s'); %b vecteur contenant les coefficients du numérateur
%a vecteur contenant les coefficients du dénominateur
[h,w]=freqs(b,a);%calcul de la fonction de transfert du filtre
hh=abs(h);%modulte de la fonction de transfert
hh_db=20*log10(hh);%module en dB
aa=abs(1./h);%atténuation en module
aa_db=20*log10(aa);%atténuation en dB
f=w/(2*pi);%axe de frequence en Hz à la place de rad/s
subplot(211);plot(f,hh_db);grid;
title('module de la fonction de transfert du filtre en dB')
xlabel('fréquence en Hz');
subplot(212);plot(f,aa_db);grid;
title('Atténuation du filtre en dB');
grid;axis([0 500 0 60]);
4) Dessiner l'atténuation du filtre en dB sur la même figure1
5) Que remarquez-vous ?
………………………………………………………………………………………………………………………………………………………..
……………………………………………………………………………………………………………………………………………………….
B) Analyse du filtre conçut
1) Construire un signal composé de deux fréquences sinusoïdales f1= 100 et f2=400 :
fe=2000;
f1=100;
f2=400;
t=0:1/fe:1;
x1=cos(2*pi*f1*t);
x2=cos(2*pi*f2*t);
x=x1+x2;
2) Tracer le spectre de x(t) :
spx=(abs(fft(x)))/fe;
fr=linspace(0,fe,length(t));
figure;subplot(211);plot(t,x);grid;axis([0 5/f1 -1.5 1.5])
xlabel('temps en s');title('signal x(t)')
subplot(212);plot(fr,spx);grid;axis([0 1000 0 1]);
xlabel('frequence en Hz');
title('spectre du signal x(t)');
3) Opération de filtrage du signal x(t) par le filtre passe-bas conçut
syst=tf(b,a);%création du filtre
s=lsim(syst,x,t);%filtrage de x par le filtre
figure;subplot(211),plot(t,s);grid;axis([5/f1 10/f1 -1.5 1.5])
title('signal filtré');
4) Observer le spectre de la sortie du filtre , que remarquez-vous
…………………………………………………………………………………………………………………………………………………
…………………………………………………………………………………………………………………………………………………
spy=(abs(fft(s)))/fe;%spectre du signal de sortie du filtre
subplot(212);plot(fr,spy);grid;axis([0 1000 0 1]);
xlabel('fréquence en Hz');
title('spectre du signal filtré');
1) Conclusion
-2-
