26.11.2014 Views

Exercices et quelques corrigés

Exercices et quelques corrigés

Exercices et quelques corrigés

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

E.N.S.I.L Deuxième Année ELT<br />

EXERCICE<br />

TRAITEMENT DU SIGNAL (FILTRAGE)<br />

1- Nous devons concevoir un filtre pour séparer le signal audio du signal ADSL. Supposons que le signal audio<br />

est un signal passe bas limité jusqu'à 3200 Hz <strong>et</strong> que le spectre du signal ADSL commence à partir de 20<br />

kHz. Nous tolérons une atténuation maximale de 1 dB pour le signal audio <strong>et</strong> souhaitons avoir une<br />

atténuation minimale de 50 dB sur le signal ADSL. Donner le H(p) du filtre à réaliser pour les cas suivant :<br />

Filtre Butterworth<br />

- Calculer le n (3,5)<br />

- Tracer le gabarit du filtre prototype<br />

- Donner le Hp(P) du filtre prototype<br />

- Dénormaliser pour obtenir H(p)<br />

- Donner le circuit LC du filtre pour les résistances de source <strong>et</strong> de charge égale à 50 ohm.<br />

Filtre Chebyshev<br />

- Calculer le n (2,83)<br />

- Tracer le gabarit du filtre prototype<br />

- Donner le Hp(P) du filtre prototype<br />

- Dénormaliser pour obtenir H(p)<br />

- Donner le circuit LC du filtre pour les résistances de source <strong>et</strong> de charge égale à 50 ohm.<br />

Solution<br />

Filtre Butterworth : n=3,5 => 4<br />

A l’aide de la table, le résultat suivant est obtenue :<br />

A p (p) = p 4 + 2,613p 3 + 3,414p 2 + 2,613p + 1<br />

A(p) = A p (p) avec p =<br />

= 4,2*10 -5 p<br />

Où ε = = 0,509 => ε 1/4 = 0,84<br />

A(p) = (4,2*10 -5 ) 4 p 4 + 2,613*(4,2*10 -5 ) 3 p 3 + 3,414(4,2*10 -5 ) 2 p 2 + 2,613(4,2*10 - 5)p + 1<br />

Filtre de Chebyshev :<br />

Amin = 50, Amax = 1, Ω p =1, Ω a = = 6,25 n=2,83 => 3<br />

1


E.N.S.I.L Deuxième Année ELT<br />

AdB<br />

50<br />

1<br />

1<br />

6,25<br />

Ω<br />

Ω =<br />

ε = 1 dB => table : A(p) = 2 n-1 ε (p 3 + 0,988p 2 + 1,238p + 0,491)<br />

H(p) = H p (p) où p= =<br />

A(p) = 0,25*10 -12 p 3 + 0,5*10 -8 p 2 + 0,125*10 -3 p + 1<br />

2- Nous allons concevoir un filtre Chebyshev inversé avec les paramètres suivants :<br />

N=2, Ω a =2, Ω p =1, ε=1<br />

Donnez le H(p)<br />

3- Nous avons un signal électrique qui contrôle la sortie d'une source ultrason. Ce signal est contaminé par un<br />

bruit 50 Hz qui se trouve dans l'environnement du travail. Il nous faut donc un filtre coupe bande afin<br />

d'éliminer c<strong>et</strong>te interférence indésirable.<br />

Il nous faut donc au moins 19 dB d'atténuation autours de 50 Hz avec une largeur de bande coupée de 10 Hz<br />

(45-55). Ensuite, dans les bandes passantes, il ne faut pas plus de 2 dB d'atténuation pour les fréquences<br />

inférieures à 35 <strong>et</strong> supérieur à 100 Hz.<br />

• Tracer le gabarit du filtre désiré.<br />

• Transformer ce gabarit en gabarit d'un filtre prototype<br />

• Calculer l'ordre du filtre pour une approximation Butterworth<br />

• Utilisant le tableau donné en cours, donner la fonction de transfert du filtre prototype<br />

• Dénormaliser c<strong>et</strong>te fonction de transfert pour obtenir la fonction de transfert du filtre réel<br />

• Identifier les éléments du circuit ci-dessous pour réaliser le filtre passe bas prototype<br />

L<br />

Vin<br />

C<br />

R=<br />

50<br />

Vout<br />

• Transformer ce filtre en filtre équivalent coupe-bande réel.<br />

4- Calculer la fonction de transfert des circuits ci-dessous :<br />

2


E.N.S.I.L Deuxième Année ELT<br />

V1<br />

R 1<br />

C<br />

V1<br />

R1<br />

R2<br />

V2<br />

Vout<br />

V2<br />

Vout<br />

R 2<br />

R1<br />

Utilisant ces circuits, donner le circuit réalisant les fonctions ci-dessous :<br />

1<br />

1<br />

H ( p)<br />

= <strong>et</strong> H ( p)<br />

=<br />

2<br />

p<br />

1+<br />

1+<br />

ap + bp<br />

a<br />

5- Utilisant la méthode à base d’intégrateurs, donner le circuit actif qui synthétise la fonction de transfert cidessous<br />

:<br />

H ( p)<br />

=<br />

2<br />

p + 2 p + 2<br />

2<br />

p + 3p<br />

− 5<br />

6- Ci-contre le diagramme des pôles <strong>et</strong> des zéros du filtre<br />

2<br />

a0 p + a1 p + a2<br />

H ( p)<br />

=<br />

est donné. Quel est le type de ce filtre (passe<br />

3 2<br />

b0 p + b1 p + b2 p + b3<br />

bas, passe haut, passe bande, …) ? Justifier.<br />

×<br />

×<br />

×<br />

jω<br />

2<br />

σ<br />

Quelle est la phase pour f = infinie ? Justifier.<br />

7- Un filtre passe haut est à concevoir. Le gabarit de ce filtre est donné sur la figure ci-dessous.<br />

A(jω)| dB<br />

20 dB<br />

1 dB<br />

500 1500<br />

f(Hz)<br />

a- Normalisez ce gabarit pour obtenir un filtre passe bas prototype. Donner le gabarit du filtre prototype.<br />

b- Nous voulons implanter ce filtre par la méthode Butterworth. Quel sera l'ordre de ce filtre ?<br />

c- Où se trouvent les pôles de la fonction de transfert résultant ?<br />

d- Utilisant le tableau, quelle est la fonction de transfert H(p) de ce filtre prototype ?<br />

e- Dénormalisez le filtre <strong>et</strong> calculez la fonction de transfert du filtre passe haut réel.<br />

3


E.N.S.I.L Deuxième Année ELT<br />

8- (Problème posé à l'examen de l'année 2002-2003) Supposons H(z)=1-0.5z -1 la fonction de transfert d'un<br />

filtre numérique.<br />

1. Donner h(n) (formes mathématique <strong>et</strong> graphique). Est-ce un filtre causal ?<br />

2. Présenter les pôles <strong>et</strong> les zéros de c<strong>et</strong>te fonction sur le plan Z. En déduire, de manière qualitative, le<br />

comportement du système dans le domaine des fréquences (passe bas, passe bande, …).<br />

3. Calculer |H(e jω )| <strong>et</strong> le tracer.<br />

4. R<strong>et</strong>racer la même fonction de transfert si on suppose une fréquence d'échantillonnage à 10 kHz<br />

(abscisse est la fréquence réelle)<br />

5. Afin de supprimer les bases fréquences, sur le plan Z, où est-ce que l'on doit ajouter un zéro (ou<br />

éventuellement un pôle) à c<strong>et</strong>te fonction de transfert ? Donner la nouvelle fonction de transfert H(z) (il<br />

nous faut toujours un système causal). Est-ce que c'est normal que la somme des coefficients de H(z)<br />

fasse zéro ? Justifier.<br />

6. Donner le nouveau h(n) (formes mathématique <strong>et</strong> graphique).<br />

9- Dans le problème qui suit, nous cherchons à synthétiser un filtre numérique équivalent à un filtre analogique<br />

passe-bas du premier ordre de constante de temps τ = 1 ms. Nous allons utiliser la méthode d'invariance<br />

impulsionnel.<br />

I) Filtre analogique :<br />

1) Donner l’expression de la fonction de transfert H a (jω) = S(jω)/E(jω). Tracer le diagramme de bode de<br />

H a (jω).<br />

2) Tracer la réponse impulsionnelle h a (t) de c<strong>et</strong>te structure.<br />

C1<br />

R1<br />

R1 = R2 = 1 KΩ<br />

C1 = 1 µF<br />

e(t)<br />

R2<br />

s(t)<br />

II) Synthèse par invariance impulsionnelle.<br />

Nous remplaçons le filtre analogique par un filtre numérique. Celui-ci est contenu dans la chaîne suivante :<br />

fe<br />

e (t) Filtre anti<br />

échantillonner<br />

e (n) FILTRE s CNA<br />

s (t)<br />

(n)<br />

repliement<br />

CAN<br />

numérique<br />

bloquer<br />

Enlève toutes les<br />

composantes spectrales<br />

>= à Fe/2<br />

Nous cherchons une fonction de transfert en Z qui aura la même réponse impulsionnelle que le filtre analogique<br />

à AOP ci dessus.<br />

4


E.N.S.I.L Deuxième Année ELT<br />

Filtre à réponse impulsionnelle finie RII.<br />

La période d’échantillonnage est Te = 1 ms.<br />

1) Calculer le h(n)<br />

2) Donner la TZ notée H(z) de h(n).<br />

3) Donner H(e jωTe ) <strong>et</strong> calculer H(e j0 ). Que constatez-vous ?<br />

4) Est-ce qu'on pouvait utiliser la transformation de c<strong>et</strong>te méthode directement sur le H a (p)?<br />

Si la période d’échantillonnage est Te = 0,5 ms.<br />

5) Donner la TZ notée H(z) de h(n).<br />

6) Calculer <strong>et</strong> tracer H(e jωTe )<br />

Comment faut-il choisir Fe pour que le filtre numérique ressemble au filtre analogique ? Expliquer d’où vient la<br />

différence des fonctions de transfert entre les filtres numériques <strong>et</strong> le filtre analogique.<br />

Est-ce que le bloqueur présenté sur le schéma modifie la réponse fréquentielle du système ?<br />

Utilisation MATLAB<br />

Pour tracer le diagramme de bode sous matlab :<br />

Sys=tf(1000,[1 1000]; bode (sys); (trace le digramme de bode de 1000/(p+1000) )<br />

0<br />

Bode Diagram<br />

-5<br />

-10<br />

Magnitude (dB)<br />

-15<br />

-20<br />

-25<br />

-30<br />

-35<br />

-40<br />

0<br />

Phase (deg)<br />

-45<br />

-90<br />

10 1 10 2 10 3 10 4 10 5<br />

5


E.N.S.I.L Deuxième Année ELT<br />

fe=10000; B=1000; A=[1 1000];<br />

[Ha,Fa]=freqs(B,A);<br />

[Bz,Az]=impinvar(B,A,fe);<br />

[Hz,Fz]=freqz(Bz,Az);<br />

plot(Fz*fe/pi/2,20*log10(abs(Hz)),'r',Fa/pi/2,20*log10(abs(Ha))); grid;<br />

5<br />

0<br />

-5<br />

-10<br />

-15<br />

-20<br />

-25<br />

0 500 1000 1500 2000 2500<br />

10- Considérer la fonction de transfert passe bas H(p) suivant :<br />

6<br />

H a<br />

( p)<br />

=<br />

2<br />

p + 5p<br />

+ 6<br />

• Utiliser la méthode invariance impulsionnelle pour transformer ce filtre en numérique.<br />

6


E.N.S.I.L Deuxième Année ELT<br />

H a<br />

6 1 6 −6<br />

( p)<br />

=<br />

= +<br />

p + 2 p + 3 p + 2 p + 3<br />

h ( t)<br />

= 6e<br />

a<br />

h(<br />

n)<br />

= Th<br />

−2t<br />

a<br />

u(<br />

t)<br />

− 6e<br />

−3t<br />

( nT ) = 6Te<br />

u(<br />

t)<br />

−2nT<br />

u(<br />

n)<br />

− 6Te<br />

−3nT<br />

u(<br />

n)<br />

H<br />

6T<br />

1−<br />

e z<br />

−6T<br />

1 − e z<br />

( z)<br />

= +<br />

−2T<br />

−1<br />

−3T<br />

− 1<br />

• La méthode invariance impulsionnelle conserve la réponse impulsionnelle. Essayer de conserver la réponse<br />

à échelon (la réponse indicielle) <strong>et</strong> calculer le H(z).<br />

t<br />

−3t<br />

−2t<br />

−3nT<br />

−2nT<br />

La réponse à un échelon = s(<br />

t)<br />

= ∫ h(<br />

τ ) dτ<br />

= ( 1 + 2e<br />

− 3e<br />

) u(<br />

t)<br />

s(<br />

n)<br />

= ( 1 + 2e<br />

− 3e<br />

) u(<br />

n)<br />

−∞<br />

1 2 3<br />

S ( z)<br />

= + −<br />

−1<br />

−3T<br />

−1<br />

−2T<br />

1<br />

1 − z 1 − e z 1 − e z<br />

−<br />

Pour calculer H(z), on utilise la relation h(n)=s(n)-s(n-1) H(z)=S(z)(1-z -1 ).<br />

−1<br />

−1<br />

2(1 − z ) 3(1 − z )<br />

H ( z)<br />

= 1 + −<br />

−3T<br />

−1<br />

−2T<br />

−1<br />

1 − e z 1 − e z<br />

Réponse à un échelon<br />

Réponse impulsionnelle<br />

Réponse fréquentielle du filtre analogique <strong>et</strong> du filtre numérique.<br />

11- Supposons que nous disposons d’un signal analogique que nous voulons filtrer par un filtre passe bas<br />

numérique. On veut garder les composantes en dessous de 2 Hz (à 1 dB près) <strong>et</strong> atténuer les composantes en<br />

dessus de 4 Hz d’au moins 10 dB.<br />

jω<br />

• Tracer le gabarit du filtre analogique désiré |H(jΩ| <strong>et</strong> le gabarit recherché numérique H ( e ) supposant<br />

une fréquence d’échantillonnage égale à 10 Hz.<br />

dB<br />

7


E.N.S.I.L Deuxième Année ELT<br />

Nous procédons de la conception de ce gabarit résultant dans le domaine numérique, en repassant dans le<br />

domaine analogique.<br />

• Transformer ce gabarit dans le domaine analogique <strong>et</strong> calculer le filtre Chebyshev analogique correspondant<br />

(H a (p)) pour des méthodes suivantes (vous pouvez utiliser les tables pour obtenir le polynôme<br />

correspondant). Pour chaque cas, calculer le H(z) du filtre numérique.<br />

- invariance impulsionnelle<br />

- bilinéaire<br />

• Donner la structure d’implantation du filtre obtenu par la méthode bilinéaire.<br />

Le système devient alors le suivant :<br />

Filtre anti<br />

repliement<br />

CAN<br />

Filtre<br />

numérique<br />

CNA<br />

• Quel est le rôle du filtre anti-repliement ?<br />

• Si on utilisait la méthode Kaiser pour calculer directement le filtre numérique, quel serait l’ordre du filtre.<br />

Solution :<br />

Le gabarit du filtre analogique est le suivant :<br />

Puisqu’il y a une relation linéaire entre le système analogique <strong>et</strong> le système échantillonné, on peut<br />

ω Ω Ω Ω<br />

écrire = = = T ⇒ ω = ΩT<br />

. C’est à dire : ω<br />

p<br />

= 2 π *2 /10 = 2 π / 5 , ωa<br />

= 2 π *4 /10 = 4 π / 5<br />

2 π Ω 2 π f 2 π<br />

e<br />

e<br />

Le gabarit numérique à satisfaire <strong>et</strong> le suivant :<br />

Avec la méthode invariance impulsionnelle le gabarit numérique doit être transformé en analogique. Il est plus<br />

simple d’obtenir un filtre analogique prototype. Ainsi on pourra utiliser directement les tables. Pour cela, on fait<br />

2 π / 5<br />

4 π / 5<br />

de telle sorte que Ω’ p =1. Ω ′<br />

p<br />

= = 1 ⇒ T = 2 π / 5 . Alors Ω ′<br />

a<br />

= = 2 . Avec ces valeurs on calcule le<br />

T<br />

2 π / 5<br />

n :<br />

0.1Amin<br />

−1<br />

10 −1<br />

0.1Amax<br />

10 −1<br />

−1<br />

Ωa<br />

Ω<br />

p<br />

cosh<br />

n =<br />

cosh /<br />

=1.86 n=2<br />

8


E.N.S.I.L Deuxième Année ELT<br />

Utilisant la table pour le filtre de Chebyshev à 1 dB on a :<br />

2<br />

A( p) = p + 1.098 p + 1.103<br />

1 1<br />

0.1Amax<br />

Alors : H ( p)<br />

= ≈<br />

car ε = 10 −1 ≈ 0.5 .<br />

n−1 2 2<br />

2 ε ( p + 1.098 p + 1.103) p + 1.098 p + 1.103<br />

Il faut maintenant calculer H(z) sachant que<br />

1 T<br />

⇒<br />

p − p 1−<br />

e z −<br />

0<br />

p0T<br />

1<br />

0 1 0<br />

1<br />

( ) A A<br />

A T AT<br />

H p = + ⇒ H ( z)<br />

= + avec T = 2 π / 5 .<br />

a 0 1<br />

1 1<br />

p p0 p p1<br />

1 e<br />

p T −<br />

z 1 e p T −<br />

− − − − z<br />

En remplaçant, on aura :<br />

−1<br />

0.63z<br />

H ( z)<br />

=<br />

1 2<br />

1 −<br />

−<br />

−<br />

0.43z<br />

+ 0.25z<br />

Méthode bilinéaire<br />

Pout obtenir un filtre prototype analogique, on fixe le Ω<br />

p<br />

= 1 <strong>et</strong> on calcule la valeur de T.<br />

2 ω<br />

p<br />

0.4π<br />

Ω<br />

p<br />

= tan = 1 ⇒ T = 2 tan = 1.45<br />

T 2<br />

2<br />

De la même manière Ω<br />

a<br />

= 4.23 . Avec ces valeurs on peut calculer le n.<br />

On obtient le même H(p) :<br />

C<strong>et</strong>te fois ci :<br />

0.1Amin<br />

−1<br />

10 −1<br />

0.1Amax<br />

10 −1<br />

−1<br />

Ωa<br />

Ω<br />

p<br />

cosh<br />

n =<br />

cosh /<br />

=1.159 2<br />

1 1<br />

H ( p)<br />

= ≈<br />

n−1 2 2<br />

2 ε ( p + 1.098 p + 1.103) p + 1.098 p + 1.103<br />

0.218 + 0.436 + 0.218<br />

H ( z) = H ( p)<br />

=<br />

−1 −2<br />

z z<br />

−1<br />

2 1−<br />

z<br />

p=<br />

−1 −2<br />

T 1<br />

1 z<br />

−<br />

+ 1− 0.35z<br />

+ 0.33z<br />

Remarque : c<strong>et</strong>te expression peut être obtenue avec matlab utilisant la commande :<br />

[A,B]=bilinear(1,[1 1.098 1.103],1/1.45) ;<br />

Remarque : Utilisant la commande freqz de matlab, dessiner les deux fonctions de transfert. Pourquoi le filtre<br />

résultant de la méthode invariance impulsionnelle ne respecte pas le gabarit souhaité alors que la méthode<br />

bilinéaire le respecte ?<br />

Structure directe II<br />

b = 0.218, b = 0.436, b = 0.218, a = 0.35, a = − 0.33<br />

0 1 2 1 2<br />

9


E.N.S.I.L Deuxième Année ELT<br />

Le rôle du filtre anti-repliement est de supprimer les composantes fréquentielles éventuelles du signal analogique<br />

en entrée au-dessus de la fréquence f e /2. Ces fréquences ne respectant pas le critère de Shannon, allaient créer du<br />

bruit dû au recouvrement de spectre.<br />

Utilisation de la méthode Kaiser.<br />

−0.5<br />

20logδ 2<br />

= − 10 ⇒ δ 2<br />

= 10 = 0.31<br />

20log(1 − δ ) = − 1 ⇒ δ1 = 0.1<br />

δ = min{ δ1, δ2} = 0.1⇒ A = 20<br />

Utilisant la fenêtre Kaiser ⇒ β = 0 ⇒ Fenêtre rectangulaire<br />

1<br />

4 − 2 20 −8<br />

∆ ω = 2 π <strong>et</strong> M = = 4.18 ⇒ 5<br />

10<br />

2.285∆ω<br />

Avec Matlab : [n,wn,b<strong>et</strong>a]=kaiserord( [0.4*pi 0.6*pi] , [1 0] , [0.31 0.1] , 2*pi) ; qui donnera n=5, wn=0.6,<br />

b<strong>et</strong>a=0.<br />

12- Un filtre numérique est à implanter. La fréquence d'échantillonnage est fixée à 8kHz. Idéalement, on veut<br />

supprimer la bande fréquentielle de 1kHz à 2 kHz.<br />

a- Tracez, pour ce filtre idéal, le module de H d (e jω ) (ω varie de -π à π).<br />

b- Si on utilise la méthode d'échantillonnage en fréquence en prenant N points sur c<strong>et</strong>te réponse idéale,<br />

quel serait, approximativement, |H(e jω )| dB (pas de calcul, pas de formule) ? Expliquez-vous.<br />

c- Quelle procédure peut-on envisager pour augmenter l'atténuation dans la bande coupée ?<br />

d- On veut utiliser la méthode Kaiser pour la conception de ce filtre (prenez δ=0.01, largeur de la bande de<br />

transition pour chaque discontinuité 500 Hz). Calculez la taille du filtre ? Calculez le h d (n). Donnez<br />

l'expression de h(n).<br />

Solution :<br />

1<br />

1<br />

ω = ΩT,<br />

alors ω1<br />

= 2π1000 = π / 4 , ω2<br />

= 2π 2000 = π / 2<br />

8000<br />

8000<br />

Si on utilise la méthode d’échantillonnage fréquentielle, on prend <strong>quelques</strong> points sur la réponse idéal <strong>et</strong> on<br />

calcule la transformé inverse de Fourier. Ce sui va nous donner un h(n). H(e jω ) la TF de ce H(n), va<br />

effectivement passer par ces points mais en dehors de ces points le comportement du filtre peut dévier du filtre<br />

idéal. Pour améliorer, il faudra prendre plus de points, d’où un filtre plus long.<br />

10


E.N.S.I.L Deuxième Année ELT<br />

Pour augmenter l’atténuation dans la bande coupée, on peut prendre un point dans la bande de transition qui est<br />

différent de 0 ou 1. La valeur de ce point est à optimiser pour maximiser l’atténuation dans la bande coupée. En<br />

contre partie, la bande de transition sera plus large. Voir figure ci-dessous<br />

A=-20log 10 δ=40, ensuite<br />

⎧0.1102(<br />

A −8.7)<br />

A > 50<br />

⎪<br />

1 π<br />

0.4<br />

β = ⎨0.5842(<br />

A − 21) + 0.07886( A − 21) 21≤<br />

A ≤ 50 alors, β = 3.4 , ∆ ω = 2π<br />

500 = , 8000 8<br />

⎪<br />

⎩0.0<br />

A < 21<br />

A −8<br />

M = = 35.66 M=36. La fenêtre Kaiser est connue maintenait :<br />

2.285∆ω<br />

2<br />

I ⎡<br />

0<br />

β (1 [( n a) / a]<br />

⎤<br />

⎢<br />

− −<br />

⎥<br />

w( n) =<br />

⎣<br />

⎦<br />

0 ≤ n ≤ M<br />

I<br />

0<br />

[ β ]<br />

Il faudra juste multiplier c<strong>et</strong>te fenêtre par la réponse idéale du filtre :<br />

1 π<br />

jω<br />

jωn<br />

sin(<br />

hd<br />

( n) = H<br />

d<br />

( e ) e dω<br />

2π<br />

∫ = n π / 4) − sin( n π / 2)<br />

− π<br />

nπ<br />

Alors le h(n) est le suivant<br />

h( n) = h ( n − 18) w( n)<br />

d<br />

13- (Problème posé à l'examen de l'année 2001-2002) Nous avons à notre disposition un système à<br />

microprocesseur qui nous perm<strong>et</strong> de réaliser des filtres numériques. Nous voulons donc implanter le système<br />

suivant :<br />

X(Ω)<br />

CAN H(e jω ) CNA<br />

idéal<br />

Y(Ω)<br />

- Où faut-il ajouter un filtre anti-repliement <strong>et</strong> quel est son rôle ?<br />

Nous cherchons à satisfaire le gabarit ci-dessous :<br />

20 dB<br />

20log|A(j2πf)|<br />

1 dB<br />

1 5<br />

f(kHz)<br />

a- Sachant que la fréquence d'échantillonnage est de 20 kHz, tracer le gabarit du filtre numérique à satisfaire<br />

b- Utilisant la méthode bilinéaire pour calculer le filtre numérique, r<strong>et</strong>ransformer le gabarit numérique obtenu<br />

en gabarit analogique<br />

c- Calculer la fonction de transfert H(p) du filtre Butterworth prototype satisfaisant ce gabarit analogique<br />

d- Dénormaliser ce H(p)<br />

11


E.N.S.I.L Deuxième Année ELT<br />

e- Donner l'expression du filtre numérique résultant (calcul exact de H(z) n'est pas nécessaire)<br />

14- Pour le filtre suivant, calculer H(z).<br />

x(n)<br />

a 1<br />

b 0<br />

z -1<br />

b 1<br />

z -1<br />

y(n)<br />

a 2<br />

b 2<br />

15- Pour un filtre RIF anti-symétrique de taille N, nous avons h(n)=-h(N-1-n). Quelle est la phase de ce filtre<br />

pour 0


E.N.S.I.L Deuxième Année ELT<br />

CNA idéal<br />

y(n)<br />

Numérique<br />

vers delta<br />

y T (t)<br />

Filtre passe bas idéal<br />

f c =1/(2T)<br />

y a (t)<br />

∑ ∞<br />

n=<br />

−∞<br />

y ( t)<br />

= y(<br />

n)<br />

δ ( t − nT )<br />

T<br />

Quelle est la relation entre Y a (Ω) <strong>et</strong> Y(e jω )?<br />

20- Considérer le système ci-dessous :<br />

x a (t)<br />

CAN<br />

idéal<br />

x(n)=<br />

x a (nT)<br />

H(e jω )<br />

y(n)<br />

CNA<br />

idéal<br />

y a (t)<br />

T<br />

T<br />

Quelle la fonction de transfert équivalente du système?<br />

21- Supposons qu'en absence d'un CNA idéal le système ci-dessous est utilisé :<br />

y(n)<br />

Numérique<br />

vers delta<br />

y T (t)<br />

h(t)<br />

1<br />

T<br />

y a (t)<br />

∑ ∞<br />

n=<br />

−∞<br />

y ( t)<br />

= y(<br />

n)<br />

δ ( t − nT )<br />

T<br />

Quelle la fonction de transfert du système ci-dessous :<br />

x a (t)<br />

CAN<br />

idéal<br />

x(n)=<br />

x a (nT)<br />

H(e jω )<br />

y(n)<br />

CNA non<br />

idéal<br />

y a (t)<br />

T<br />

T<br />

−1<br />

−2<br />

1 + 2z<br />

+ z<br />

22- Soit H ( z)<br />

=<br />

,<br />

−1<br />

−2<br />

1−<br />

0.75z<br />

+ 0.125z<br />

- écrire l'équation aux différences,<br />

- donner structures de la réalisation de la forme directe I <strong>et</strong> II.<br />

23- Supposons que H(z) est un filtre passe bas. Que peut-on dire de H(1/z)? Est-ce que le filtre résultant est<br />

stable?<br />

24- Pour la structure ci-dessous calculer la fonction de transfert. Présenter les pôles <strong>et</strong> les zéros du système<br />

correspondant sur le plan Z.<br />

13


E.N.S.I.L Deuxième Année ELT<br />

25- Calculer le module <strong>et</strong> la phase de H(e jω −1<br />

z − 0.54<br />

) de H ( z)<br />

= .<br />

−1<br />

1 − 0.54z<br />

Si c<strong>et</strong>te fonction est réalisée par la structure ci-contre, identifier<br />

les coefficients a, b <strong>et</strong> c.<br />

26- Pour les systèmes présentés ci-dessous :<br />

• calculer la fonction de transfert<br />

• dessiner la structure transposée<br />

• vérifier si la fonction de transfert du système transposé est équivalente à celle du système initial<br />

27- (examen septembre 2004)<br />

a) Expliquer le théorème de la transposition.<br />

b) Dessiner une structure d'implantation de la fonction de transfert<br />

−1<br />

1−<br />

0.5z<br />

−1<br />

H ( z)<br />

= + (1 − 2z<br />

+ z<br />

−1<br />

−2<br />

1+<br />

0.25z<br />

+ 0.1z<br />

−2<br />

)<br />

c) Utilisant le théorème de la transposition dessiner la structure transposée du circuit ci-dessus.<br />

28- (Examen 2005-2006) Un filtre passe bas est souhaité avec au maximum 1 dB d'atténuation jusqu'à 3 kHz, <strong>et</strong><br />

au moins 40 dB d'atténuation dans la bande coupée (à partir de 4 kHz). Ce filtre va être implanté par un<br />

microprocesseur avec une fréquence d'échantillonnage à 24 kHz. Combien de multiplications par seconde<br />

sont nécessaires pour chacune des méthodes suivantes :<br />

a- filtre RII avec la méthode invariance impulsionnelle (le filtre analogique abstrait est de type<br />

Butterworth) (Structure directe II)<br />

b- filtre RII avec la méthode bilinéaire (le filtre analogique abstrait est de type Butterworth) (Structure<br />

directe II)<br />

c- Filtre RIF de type Kaiser (Structure transversale <strong>et</strong> profitant de la symétrie)<br />

d- Filtre RIF avec l'approximation equi-ondulation (Structure transversale <strong>et</strong> profitant de la symétrie)<br />

14


E.N.S.I.L Deuxième Année ELT<br />

solution<br />

A max = 1dB ; A min = 40dB ; ω p = 3000*2π ; ω a = 4000*2π ; k = 24 kHz<br />

a) Une transformation linéaire :<br />

ω = Ω T ,<br />

2π<br />

3000 π<br />

ω<br />

p<br />

= = ,<br />

24000 4<br />

La méthode invariance impulsionnelle ;<br />

2π<br />

4000 π<br />

ω<br />

a<br />

= =<br />

24000 3<br />

Butterworth : = 18,35 => n = 19<br />

Transformation invariance impulsionnelle :<br />

H a (p) =<br />

Directe II :<br />

19*2=38 multiplications par d’échantillons en sortie <br />

38*fe=38*24000=912000 mult/sec<br />

b) On calcule le gabarit du filtre analogique abstrait :<br />

H(Z) =<br />

T=1 => Ω = => Ω a = = = 1,1547<br />

nc = 15,9 => nc = 16<br />

Ω p = = = 0,8284<br />

H(p) = => H(Z) = H(p) où p=2<br />

H(Z) =<br />

Directe II :<br />

17 + 16=33 multiplications par d’échantillons en sortie<br />

33*fe=38*24000=792000 mult/sec<br />

c) Filtre kaiser :<br />

-20log(1- δ 1 ) = A min => δ 1 = 1-10 -1/20 = 1-0,89 = 0,11<br />

-20logδ 2 = 40 => δ 2 = 0,01 ; ∆ = min(δ 1 , δ 2 ) = 0,01 => A=40<br />

M =<br />

H(Z) = h 0 + h 1 Z -1 + … + h 54 Z -54<br />

Filtre à phase linéaire : il faut 54/2 + 1 = 27 + 1 = 28 mult/éch<br />

28*24000 = 672000 mult/sec<br />

d) Filtre equi-ondulation<br />

N =<br />

Filtre à phase linéaire : 32/2 + 1 = 17 mult/éch<br />

17*24000 = 408000 mult/sec<br />

15


E.N.S.I.L Deuxième Année ELT<br />

29- (Examen 2005- 2006) Le comportement fréquentiel idéal d'un filtre analogique est présenté sur la figure. On<br />

nous a donné un H(P) qui présente approximativement le même comportement fréquentiel. On remplace P<br />

−1<br />

1−<br />

z<br />

par P = 20 pour obtenir un H(z). Tracez H(e jω ) entre (-π,π).<br />

1<br />

−1<br />

+ z<br />

1<br />

0.8<br />

0.5<br />

|H(jΩ)|<br />

2 5 6<br />

Ω<br />

30- La fonction de transfert d’un filtre Bessel d’ordre 3, donnant un r<strong>et</strong>ard de T 0 =1 sec, est la suivante :<br />

15<br />

H ( p)<br />

=<br />

( T p) + 6( T p) + 15( T p) + 15<br />

3 2<br />

0 0 0<br />

Si on trace H ( jω ) , on constate que c’est un filtre passe bas avec ω − 3<br />

≃ 1.75 . On souhaite utiliser un filtre de<br />

Bessel d’ordre 3 mais avec f − 3dB<br />

= 3000 kHz. Quelle est la fonction de transfert correspondant ? Quel sera le<br />

r<strong>et</strong>ard de ce filtre ?<br />

Le circuit correspondant à la fonction de transfert ci-dessus est le suivant :<br />

dB<br />

Donner le circuit correspondant pour des charges 50 ohm..<br />

Solution :<br />

C’est une transformation passe-bas passe-bas ω − 3dB<br />

≃ 1.75 doit être transformé en 2π<br />

3000 . Alors<br />

2π<br />

× 3000<br />

ωr<br />

= = 10771<br />

1.75<br />

15<br />

H ( p) = H1( p ) =<br />

p=<br />

p/<br />

ω<br />

3 2<br />

r<br />

⎛ T0 p ⎞ ⎛ T0 p ⎞ ⎛ T0<br />

p ⎞<br />

⎜ 6 15 15<br />

10771<br />

⎟ + ⎜ + +<br />

10771<br />

⎟ ⎜<br />

10771<br />

⎟<br />

⎝ ⎠ ⎝ ⎠ ⎝ ⎠<br />

1<br />

Où T<br />

0<br />

= 1 . Le délai de groupe devient τ<br />

g<br />

= = 92.8µs.<br />

10771<br />

Le circuit est le suivant :<br />

16


E.N.S.I.L Deuxième Année ELT<br />

50 L<br />

V in<br />

C 50<br />

1 C 2<br />

1.255 1<br />

0.192 1<br />

c<br />

1<br />

= × = 2.33µF<br />

, c<br />

1<br />

= × = 0.357µF<br />

,<br />

10771 50<br />

10771 50<br />

0.553 1<br />

L = × = 2.57mH<br />

10771 50<br />

31- Utilisant des intégrateurs, synthétiser la fonction de transfert ci-dessous :<br />

2<br />

p + 2 p + 2<br />

H ( p)<br />

=<br />

2<br />

p + 3p<br />

− 5<br />

Solution : On vérifie d’abord si les pôles sont sur le demi-plan gauche (stabilité). Puis :<br />

2 2<br />

Y p + 2 p + 2 1+ 2 / P + 2 / P<br />

H ( p)<br />

= = =<br />

2 2<br />

X p + 3p − 5 1+ 3 / p − 5 / P<br />

1 ⎧<br />

⎨<br />

1<br />

p ⎩ p<br />

On peut donc écrire : Y = X + 2X − 3Y + [ 2X + 5Y<br />

]<br />

1 U 2X 5Y<br />

p<br />

1<br />

= 2 − 3 + , <strong>et</strong> Y = X + W .<br />

p<br />

= [ + ] <strong>et</strong> W { X Y U }<br />

⎫<br />

⎬ . On prend<br />

⎭<br />

Avec ce circuit, on a V<br />

out<br />

1 ⎡ X Y ⎤<br />

= − ⎢ + ⎥ . Si on prend C=100µF, R 1 =5 k <strong>et</strong> R 2 =2 k, on obtiendra V out =-U.<br />

Cp ⎣ R1 R2<br />

⎦<br />

1 W = 2 X − 3 Y + U<br />

p<br />

On peut aussi construire W : { }<br />

Avec R 1 =5k, R 2 =3.33 k <strong>et</strong> R 3 =10 k. Pour construire Y, on aura donc<br />

On pourra donc connecter tous les signaux manquants.<br />

17


E.N.S.I.L Deuxième Année ELT<br />

18


E.N.S.I.L Deuxième Année ELT<br />

Corrigé<br />

19


E.N.S.I.L Deuxième Année ELT<br />

20


E.N.S.I.L Deuxième Année ELT<br />

21


E.N.S.I.L Deuxième Année ELT<br />

Diagramme pôle <strong>et</strong> zéro<br />

22


E.N.S.I.L Deuxième Année ELT<br />

Filtre passe haut<br />

SimulationElectronics workbench<br />

23


E.N.S.I.L Deuxième Année ELT<br />

24

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!