Exercices de traitement numérique du signal

Gabriel Dauphin

1 Cours A : description d’un signal

1.1 Exercices d’application
Exercice 1 (56) On considère un signal temps discret non-périodique défini par xn =
δn − 1.1δn−4 avec fe = 2Hz.
1. Que devient le signal quand on amplifie par un facteur 2 ?
2. Que devient le signal quand on lui ajoute 2 ?
3. Que devient le signal quand on dilate l’échelle des temps par un facteur 2 ?
4. Que devient le signal quand on retarde le signal d’une seconde ?
5. Que devient le signal quand on le quantifie sur 2bits, donnez le résultat graphi-
quement ?
Dans chacun des cas représentez sur une figure ce que devient le signal.

1
1.2 Exercices pour approfondir
Exercice 2 (29) On considère un signal s1(t) = cos(2πt) et s2(t) = | cos(2πt)| où t repré-
sente le temps mesuré en secondes.
1. Représentez s1(t) et s2(t) sur un graphique pour t ∈ [0, 2].
2. Montrez que s1 est périodique de période 1.
3. Proposez une formule à appliquer pour calculer la puissance du signal ?
1+cos(4πt)
4. Démontrez la formule trigonométrique cos2(2πt) = 2
5. Déduisez la puissance de s1.
6. Montrez que s2 est périodique de période 1/2.
7. Proposez une formule à appliquer pour calculer la puissance, si possible la même
que la précédente.
8. Montrez que la puissance de s2 est la même que la puissance de s1.

Exercice 3 (ex28) On considère un robinet qui goutte. On considère que les gouttes d’eau
sont de même taille et ont un volume de 1/20mL. Le débit de la moyen de la fuite est de
0.3Lḣ−1. Expliquez comment ce phénomène peut se modéliser par :
1. un signal temps continu à valeurs réelles,
2. un signal temps continu à valeurs discrètes,
3. un signal temps discret à valeurs réelles,
2
 4. un signal temps discret à valeurs discrètes.
Pour chacun de ces modèles indiquez la période d’échantillonnage et la fréquence d’échan-
tillonnage lorsque cela est nécessaire.

2 Cours B : Echantillonnage d’un signal

2.1 Exercices d’application
Exercice 4 (55) On considère un signal dont les mesures aux instants : t = 0, t = 15s,
t = 30s sont les suivantes 0.5, 0, 1.5.
1. Montrez comment on peut interpréter ces mesures comme celles associées à un
signal temps discret non-périodique. Quelle est la fréquence d’échantillonnage ?
2. Trouvez l’énergie correspondante.
3. Montrez comment on peut interpréter ces mesures comme celles associées à un
signal temps discret périodique. Représentez graphique le signal correspondant.
4. Trouvez la puissance correspondante.
5. Après observation précise de la figure 1, montrez comment au moyen de deux
sinusoïdes et du signal constant égal à 1, on peut interpréter ces mesures comme
celles associées à un signal temps continu périodique.
6. Trouvez la puissance correspondante.
3
F IGURE 1 – Représentation de deux sinusoïdes auquel on a ajouté 1 et de la somme de
4 1. Exercice 4
ces deux sinusoïdes auquel on a encore ajouté
2.2 Exercices pour approfondir
Exercice 5 (33) Un filtre anti-repliement de spectre est souvent placé avant l’échantillon-
nage. A quoi est-ce que cela sert ? Ce filtre est souvent analogique, comment pourrait-on
utiliser un filtre numérique à la place ?

3 Cours C : Série de Fourier, transformée de Fourier

3.1 Exercices d’application
Exercice 6 (51) On considère le signal temps continu et périodique de période 2 défini
par sur [0, 2] par x(t) = 1[0,1](t). Calculez la transformée de Fourier et représentez gra-
phiquement le module de la transformée de Fourier en fonction de la fréquence.

Exercice 7 (53) On considère trois signaux temps continu, x(t), y(t), z(t).
— x(t) est périodique de période 2 et pour t ∈ [0, 2[, il est défini par x(t) = 1[0,1](t).
— y(t) n’est pas périodique et pour t ∈ R, il est défini par x(t) = 1[0,1](t).
— z(t) est périodique de période T et pour t ∈ [0, T [, il est défini par x(t) = 1[0,1](t).
1. Représentez sur un même graphique pour t ∈ [0, 4], x(t), y(t), z(t) avec T = 3
2. Calculez la transformée de Fourier de x(t).
3. Calculez la transformée de Fourier de y(t).
5
 4. Calculez la transformée de Fourier de z(t) en l’exprimant à partir de Yb (f ).
5. Représentez les trois spectres pour f ∈ [−2, 2] avec T = 4.
Exercice 8 (30) On cherche à calculer la transformée de Fourier de s(t) = sin2(2πt) =
1−cos(4πt)
2
1. Représentez sur une même figure les fonctions sin(2πt), cos(2πt), −1/2 cos(4πt) et
sin2(2πt) pour t ∈ [0, 1].
2. Ecrivez sin(2πt) comme une combinaison linéaire d’exponentielles complexes.
3. Montrez que sin(2πt) est périodique de période 1. Déduisez de ceci que la précé-
dente formule est en fait la décomposition en série de Fourier de sin(2πt) en expo-
nentielles complexes. Que valent les coefficients de la série Fourier de sin(2πt) ?
4. Que vaut la transformée de Fourier de sin(2πt) ?
5. En déduire la transformée de Fourier de cos(2πt) = − sin(2π(t − 1/4)) ? (la fonc-
tion cosinus est en avance d’un quart de période par rapport à la fonction sinus,
elle est donc en opposition de phase avec la fonction sinus retardée d’un quart de
période).
6. On observe que la fonction cos(4πt) est une contraction de la fonction cos(2πt),
calculez sa transformée de Fourier ?
7. Quelle est la transformée de Fourier de la fonction constante t 7→ 1 ?
8. En utilisant la formule trigonométrique initiale, quelle est la transformée de Fou-
rier de sin2(2πt) ?
6
 9. Calculez la transformée de Fourier inverse de celle trouvée et retrouvez la formule
trigonométrique initiale.

Exercice 9 (31) On cherche à déterminer la transformée de Fourier de

s(t) = 1[0,1](t) + 1[0,2](t)
1. Représentez le signal s pour t ∈ [0, 2].
2. Calculez la transformée
R ∞ de Fourier de s1(t) = 1[0,1](t) en utilisant la transformée
de Fourier S(f ) = −∞ s(t)e−j2πf tdt, montrez qu’elle se met sous la forme de

−jπf sin(πf )
Ŝ1(f ) = e
πf
3. Expliquez le fait que ce signal ne soit pas à valeurs réelles ?
4. Calculez la transformée de Fourier en f = 0 sans utiliser la formule plus haut.
5. Déduisez la transformée de Fourier de s2(t) = 1[0,2](t)
6. Montrez que la transformée de Fourier de s se met sous la forme suivante :
2 − e−2jπf − e−4jπf
Ŝ(f ) =
2jπf
7. Pour faciliter la représentation du module de la transformée de Fourier, il est en
général souhaitable d’exprimer ce module sous la forme de produit de fonction
7
 simple. Après avoir remarqué que le numérateur s’annule en la fréquence nulle et
effectué une factorisation, montrez que le module de la transformée de Fourier se
met sous la forme suivante :
 
 sin πf  p
|Ŝ(f )| =   5 + 4 cos 2πf
πf 

8. Dessinez à main levée le module de la transformée de Fourier pour f ∈ [−4, 4].

Exercice 10 (6)
Soit le signal défini par x(t) = 0 pour t 6∈] − 1, 3[, x(t) = t pour t ∈]1, 2[, x(t) = 2 − t
pour t ∈]0, 1[ et x(t) = 2 pour t ∈] − 1, 0[ et aussi pour t ∈]2, 3[.
1. Calculez arg(X(f )).
2. Calculez X(0).
R
3. Calculez X(f )df .
R
4. Calculez |X(f )|2df .

3.2 Exercices pour approfondir

Exercice 11 (3)
Donnez le développement en série de Fourier d’un pulse périodique de période T , de
largeur τ et d’amplitude A, centré par rapport à l’origine. En posant K = Tτ , donnez
8
le nombre de raies du lobe principal et des lobes secondaires. Que se passe-t-il pour
K → +∞ en maintenant A/K constant.

Exercice 12 (4)
Donnez la transformée de Fourier d’un pulse de largeur τ et d’amplitude A, centré
autour de l’origine. Donnez la largeur du lobe principale et des lobes secondaires. Que
se passe-t-il pour τ → 0 en maintenant Aτ constant ?

4 Cours D : TFD, TFTD

4.1 Exercices d’application
Exercice 13 (40) On considère deux signaux xn et yn définis par

xn = δn + δn−2 et yn = δn + δn−1 + δn−2 (1)

où δn est la suite nulle sauf en n = 0 où elle vaut 1. On cherche à calculer la transformée

de Fourier. La fréquence d’échantillonnage est notée fe et vaut 1kHz.
1. Dessinez les signaux xn et yn. S’agit-il de signaux à temps discret/temps continu,
s’agit-il de signaux périodiques ou non-périodiques. Quelle transformée de Fou-
rier vous semble adaptée pour de tels signaux ?
2. Calculez la transformée de Fourier de xn, notée X̂(f ).
9
 3. Retrouvez la signal xn en calculant la transformée de Fourier inverse. Pour cela il
est conseillé de traiter séparément les trois cas n = 0, n = 2, n 6∈ {0, 2}.
4. On considère un complexe z, montrez que
3/2
 −3/2 3/2

z z −z
1 + z + z 2 = 1/2 (2)
z z −1/2 − z 1/2
5. Déduisez de (2) que
3
sin( 2 θ)
1 + ejθ + e2jθ =e jθ
(3)
sin( 12 θ)

6. Utilisez (3) pour en déduire la transformée de Fourier de yn, notée Ŷ (f ).

7. Représentez sur f ∈ [−3fe/2, 3fe/2], |Ŷ (f )| en utilisant le fait qu’à basse fré-
quence cela ressemble à un sinus cardinal.

Exercice 14 (45) On considère xn, un signal temps discret périodique de période 4 échan-
tillonné à la fréquence fe = 100Hz. Les premières valeurs de xn sont x0 = x1 = 1 et
x2 = x3 = 0.
Calculez le module de la transformée de Fourier de ce signal. Représentez graphique-
ment le module de la transformée de Fourier en fonction de la fréquence.

Exercice 15 (52)
10
 On considère un signal xn échantillonné à la fréquence fe et défini par

xn = δn + δn−1 + δn−2

On définit yn = xn ∗ xn Calculez yn

4.2 Exercices pour approfondir

Exercice 16 (34)
On considère le signal périodique x1[n] de motif {1, 0, 0, 1} et le signal x2[n] pério-
dique de motif {1, 0, 0, 1, 1, 0, 0, 1}. Calculez les transformées de Fourier discrètes de ces
deux signaux. Montrez comment les deux s’expriment en fonction d’un cosinus et comment
la deuxième aurait pu se déduire de la première.

Exercice 17 (15)
On considère le signal cosinus tel que : x[k] = cos(2πk/6), observé sur une durée
limitée T=N.Te, avec comme fréquence d’échantillonnage f e = 1kHz. On considère 3
cas : N=6, N=12 et N=16.
1. Quelle est la fréquence du signal à temps discret s’il était défini sur une durée
infinie ?
2. Calculez la TFD dans les deux premiers cas. On pourra s’aider de ce que sur
l’ordinateur on trouve les résultats affichés sur la figure 2.
11
3. Le calcul de la TFD dans ces 3 configurations donne les résultats suivants montrés
sur la figure 2. Mettez les bonnes échelles en fréquences pour les trois graphiques.
Confrontez ce résultat à ceux trouvés précédemment. Expliquez pourquoi le troi-
sième cas est différent.
4. Proposez une idée pour atténuer les distorsions dans le 3ème cas ?

F IGURE 2 – s0,se,sa

12
5 Cours E : Repliement de spectre
5.1 Exercices d’application
Exercice 18 (57) On considère le spectre d’un signal défini par
1−r
X(f
b )= (4)
1 − re−j2πf Te
Le module de ce spectre est représenté sur la figure 3 pour une certaine valeur de r
1. Le signal associé à ce spectre est-il temps discret et non-périodique ? Quelle est la
fréquence d’échantillonnage ?
fe
2. À partir de (4) trouvez la valeur du module du spectre en f = 2 ? Dessinez le
graphique associé à ces valeurs en fonction de r ?
3. Quelle est la valeur de r associée à ce graphique, sachant que sur le graphique on
observe que |X(f
b e/2)| = 0.05 ?
4. À partir de la figure 3, trouvez la fréquence de coupure de ce signal, en supposant
qu’on interpréte ce spectre comme la réponse fréquentielle d’un filtre ? S’agit-il
d’un filtre passe-bas/passe-haut/passe-bande/coupe-bande/passe-tout ?

13
F IGURE 3 – Représentation du spectre pour une valeur particulière de r. Exercice 18
14
6 Cours EBis : Filtre et descripteur de signaux
Densité spectrale et autocorrélation
6.1 Exercices d’application
Exercice 19 (41) On considère une suite hn = δn − δn−1 On considère une entrée ayant
les valeurs suivantes
x0 = 1 x1 = 1 x2 = 0 x3 = 0 x4 = 1 x5 = −1
d
Calculez yn = hn ∗ xn Vous pourrez d’abord montrer que
yn = xn − xn−1
Remarquez qu’on a ici calculé la sortie yn d’un filtre de réponse impulsionnelle hn dont
l’entrée est xn.

Exercice 20 (42) On considère une filtre analogique défini par

Z t
y(t) = x(τ ) dτ
t−1

où x(t) est l’entrée et y(t) est la sortie.

1. Calculez y(t) quand x(t) = δ(t) en distinguant le cas où t < 0, t ∈ [0, 1] et t > 1.
On note h(t) le résultat trouve, c’est la réponse impulsionnelle.
15
 2. Tracez la réponse impulsionnelle.
3. Calculez la transformée de Fourier de h(t). On pourra utiliser le fait que
  sin(πf )
TF 1[−1/2,1/2] (f ) =
πf
C’est la réponse fréquentielle notée Ĥ(f ).
4. S’agit-il d’un passe-bas, passe-haut, passe-bande, coupe-bande ou passe-tout ?
Exercice 21 (43) On considère une fréquence d’échantillonnage fe = 100Hz. On consi-
dère un filtre numérique défini par
y = xn−1 (5)
1. On considère une entrée xn = δn. Calculez la sortie yn correspondant à cette
entrée. Cette sortie est notée hn, il s’agit de la réponse impulsionnelle du filtre.
2. Tracez la réponse impulsionnelle
3. Calculez la transformée de Fourier à temps discret de hn. C’est ce qu’on appelle
la réponse fréquentielle notée Ĥ(f ).
4. Tracez le module de la réponse fréquentielle. S’agit-il d’un passe-bas, passe-haut,
passe-bande, coupe-bande ou un passe-tout ?
5. Montrez en utilisant (5) que
Ŷ (f ) = Ĥ(f )X̂(f )
16
Exercice 22 (44) On considère un signal xn que l’on cherche à sur-échantillonner en
doublant la fréquence d’échantillonnage. On suppose que xn = 0 pour n < 0. Le pro-
cédé consiste à d’abord rajouter des échantillons nuls après chaque échantillon, le signal
obtenu est zn
z2n = xn z2n+1 = 0
Puis on applique un filtre au signal zn, la sortie du filtre est notée yn
yn = zn + zn−1
1. On considère le cas de xn défini par
x0 = 2 x1 = 1 x2 = −3 x3 = −2
Tracez sur le même graphique xn, zn et yn.
2. Exprimez z0, . . . z7 en fonction de x0, x1, x2, x3
3. Démontrez les relations suivantes
z2n+1 = xn z2n = xn
Exercice 23 (49) On considère la fréquence d’échantillonnage fe = 100Hz. On considère
le filtre numérique H défini par l’équation aux différences suivante
yn
yn+1 + = xn
2
où xn est l’entrée et xn est la sortie. Calculez la réponse fréquentielle.
17
Exercice 24 (46) On considère un signal temps continue non-périodique x(t) = 1[0,1](t).
Calculez la densité spectrale d’énergie. Représentez graphiquement cette densité spec-
trale d’énergie.

6.2 Exercices pour approfondir

Exercice 25 (54) On considère un signal xn
(−1)n
xn = n
1N[n]
2
échantillonné à la fréquence fe = 2Hz.
1. Sous-échantillonnez ce signal à fe0 = 1Hz en ne conservant qu’un échantillon sur
deux. Calculez le nouveau signal obtenu appelé ya[n].
2. On applique un filtre à xn, ce filtre est défini par la relation entrée sortie :
xn + xn−1
zn = (6)
2
Calculez le signal zn.
3. Sous-échantillonnez le signal zn à fe0 = 1Hz en ne conservant qu’un échantillon
sur deux. Calculez le nouveau signal obtenu appelé yb[n].
4. Tracez les xn, ya[n], yb[n], zn pour les trois premières secondes. Commentez l’inté-
rêt de considérer yb[n].
18
F IGURE 4 – Représentations des spectres X(f
b ), Yba(f ), Z(f
b ), Ybb(f ). Exercice 25
19
 5. Calculez la fréquence de coupure associée au filtre défini par (6).
6. Sur la figure 4, sont représentés les spectres X(f
b ), Z(f
b ), Yba(f ), Ybb(f ). Trouvez
quelle courbe correspond à quel spectre.

7 Cours 1F : Filtres analogiques

7.1 Exercices d’applications
Exercice 26 (5)
Calculez les transformées de Laplace de sin(2πf0t)1R+ (t) et cos(2πf0t)1R+ (t).
1. En utilisant les propriétés de dérivation, reliant le sinus au cosinus, vérifiez que la
TL du cosinus se calcule bien à partir de la TL du sinus.
2. En utilisant les propriétés d’intégration, reliant le cosinus au sinus, vérifiez que la
TL du sinus se calcule bien à partir de la TL du cosinus.

Exercice 27 (11)
1. Calculez la transformée de Laplace de s1(t) = cos(2πf0t)1R+ (t).
2. Calculez la transformée de Fourier de s1(t). Commentez les différences et ressem-
blances entre les deux formules.
3. s2(t) = sin(2πf0t)1R+ (t) s’exprime en fonction de la dérivée de s1(t) en déduire la
transformée de Laplace de s2(t). Commentez la pertinence physique de ce calcul.
20
Exercice 28 (47) On considère un filtre analogique H défini par l’équation différentielle
suivante
dy y
+ =x
dt 2
où x(t) est l’entrée et y(t) est la sortie. Montrez que ce filtre est stable.

7.2 Exercices pour approfondir

Exercice 29 (9)
On considère un filtre de transformée de Laplace H(p) = a p+1
p+b , avec a, b dans R.
1. Pour quelles valeurs de b ce filtre est-il stable ?
2. Calculez la réponse fréquentielle de ce filtre.
3. Donnez la relation entre a et b pour que Ĥ(0) = 1. A quoi sert cette relation ?
4. Représentez le module de la réponse fréquentielle quand b ∈]0, 1[ puis quand b ∈
]1, +∞[. Commentez.
5. Calculez la réponse impulsionnelle de ce filtre. Commentez sur la stabilité du filtre.
6. Reprendre les deux dernières questions en considérant une nouvelle échelle de
temps t0 = 2t appliquée à un nouveau filtre H(p) = p+1
1
.
7. Ecrire la relation entrée-sortie sous la forme d’une équation différentielle.
8. Ecrivez la relation entrée-sortie sous la forme d’une équation intégrale.
21
 R +∞
La transformée de Laplace est définie par H(p) = 0 h(t)e−ptdt
R∞
La transformée de Fourier est définie par Ĥ(f ) = −∞ h(t)e−j2πf tdt
Exercice 30 (39) On considère un filtre de réponse impulsionnelle
h(t) = 1[0,1](t)
L’entrée de ce filtres est notée x(t) et la sortie est notée y(t). Un tel filtre est appelé
moyenneur.
1. Expliquez pourquoi ce filtre est stable ?
2. Montrez que la sortie du filtre s’exprime en fonction de l’entrée
Z t
y(t) = x(τ )dτ
t−1

3. On place en entrée un échelon : x(t) = 1[0,+∞[(t) Calculez la sortie y(t) en distin-

guant le cas t < 0, t ∈ [0, 1] et le cas t > 1.
4. Représentez graphiquement la sortie du filtre y(t).
5. Calculez la fonction de transfert de ce filtre.
6. Expliquez pourquoi ce filtre n’est pas un filtre rationnel ?
7. Pourquoi peut-on dire que ce filtre est stable ?
8. Montrez que Ĥ(0) = 1. Pourquoi est-ce une propriété attendue d’un filtre moyen-
neur.
22
8 Cours 2F : Filtres numériques, MA, AR, ARMA, Trans-
formée en Z
8.1 Exercices d’application

F IGURE 5 – s0,se,sa (exercice 31, (12))

Exercice 31 (12)
On désigne par en et sn respectivement les valeurs de l’entrée et de la sortie du filtre
à l’instant nTe défini dans la figure 5 (p. 23).
1. Montrer que l’algorithme de ce filtre peut s’écrire : sn = aen + bsn−1, (a et b sont
deux coefficients constants).
23
 2. En déduire que la fonction de transfert en z de ce filtre peut s’écrire : T (z) =
a
1−bz −1
.
La transformée en Z s’écrit T Z[hn] = n≥0 hnz −n
P

Exercice 32 (48) On considère le filtre numérique H défini par l’équation aux différences
suivante
yn
yn+1 + = xn
2
où xn est l’entrée et xn est la sortie. Calculez la réponse impulsionnelle.

Exercice 33 (49) On considère la fréquence d’échantillonnage fe = 100Hz. On considère

le filtre numérique H défini par l’équation aux différences suivante
yn
yn+1 + = xn
2
où xn est l’entrée et xn est la sortie. Calculez la réponse fréquentielle.

Exercice 34 (14) On définit deux filtres. Le premier filtre est défini par sa réponse impul-
sionnelle : han = δ[n] + 2δ[n − 1] + δ[n − 2]. Le deuxième filtre est défini par sa fonction
−1
de transfert : H(z) = 1+3z
2−z −1
.
Pour les différents filtres ci-dessus, utilisés avec une fréquence d’échantillonnage de
1MHz, compléter les informations de manière à avoir :
24
 1. Le type de filtre (RII,RIF)
2. La stabilité
3. Le diagramme de pôle et de zéros
4. La réponse impulsionnelle
5. L’allure du module de la réponse fréquentielle.
6. Quelle est l’équation qui lie l’entrée et la sortie.

8.2 Exercices pour approfondir

Exercice 35 (20)
On considère le filtre de fonction de transfert H #(z) = 1 + 3z −1 + 6z −2 + 3z −3 + z −4.
Pourquoi ce filtre est à phase linéaire ? Montrez qu’il se comporte comme un filtre à retard,
quel est ce retard ?

9 Cours G : Synthèse de filtre MA, gabarit, fenêtre, inva-

riant impulsionnel
9.1 Exercices d’applications
Exercice 36 (19)
25
 On cherche à synthétiser avec un filtre à réponse impulsionnelle finie un passe-haut
de fréquence de coupure fc = fe/6, avec fe = 1kHz, à l’aide de la fenêtre triangulaire.
1. On considère une suite temps discret h1n dont la transformée de Fourier est pé-
riodique de période fe et égale à 1[−fc,fc] sur l’intervalle [−fe/2, fe/2]. Quelle
transformée de Fourier utiliser pour calculer han ?
sin(πf )
2. On sait que la transformée de Fourier de 1[−1/2,1/2](t) est πf . Montrez qu’alors
sin(πtfe /3)
la transformée de Fourier inverse de 1[−fc,fc](f ) est πt ?
3. Déduisez que han = sin(πn/3)
πn
4. On note hbn la suite dont la transformée de Fourier vaut 1[−fe/2,−fc] + 1[fc,fe/2] au
sein de l’intervalle [−fe/2, fe/2]. Que vaut hbn.
5. On souhaite maintenant un filtre causal avec une réponse impulsionnelle hcn cau-
sale sur sept termes (i.e. ailleurs que sur ces sept termes la réponse impulsionnelle
est nulle). Que vaut hcn ?
6. On utilise maintenant une fenêtre triangulaire, que vaut le filtre hdn ainsi modifié ?
7. Pourquoi le filtre obtenu est-il à phase linéaire ?
8. Les modules et les phases des filtres recherchés sont représentés sur la figure 6
(p. 27), précisez quelle courbe correspond à quel filtre.

Exercice 37 (35) Les signaux audio stéréo sont numérisés sur 16 bits à la fréquence
d’échantillonnage avec fe = 44.1kHz. log10(2) ≈ 0.3 et log10(3) ≈ 0.5.
26
 27
F IGURE 6 – Module et phase des réponse fréquentielle des différents filtres pour deux
 1. Déterminez le rapport signal sur bruit en décibel pour une sinusoïde à plein
échelle.
2. Le nombre de bits est multiplié par 1.05 (bit d’horloge, correcteur d’erreur, contrôle,
affichage). Quel est le débit du système d’enregistrement ?
3. On peut enregistrer une heure de musique sur un CD. Evaluez le nombre de bits
que l’on peut stocker.
4. A titre de comparaison : un dictionnaire peut posséder jusqu’à 1500 pagesavec
deux colonnes par page, 100 lignes par colonne, 8 mots par ligne, 6 lettres par
mot. Il faut 7 bits pour coder une lettre. Combien de dictionnaires peut on stocker
sur un CD ?

Exercice 38 (37) La fréquence d’échantillonnage est fe = 1kHz. Faites la synthèse d’un

filtre passe-bas de fréquence de coupure fc = 250Hz avec une fenêtre rectangulaire et en
utilisant que 7 termes non-nuls.

10 Exercice pour approfondir

Exercice 39 (36) On étudie l’impact d’une fenêtre sur la synthèse d’un filtre de réponse
impulsionnelle hn. La fréquence d’échantillonnage est fe = 1kHz.

28
 (a) (b)
F IGURE 7 – Transformée de Fourier de wn et wn

29
 (a) (b)
F IGURE 8 – Transformée de Fourier de xn, xnwn (avec N = 20) et xnwn (avec N = 10).

30
 (a) 1
1. Montrez que la transformée de Fourier de wn = N +1 10..N [n]

1 sin πf (N + 1)Te
TF[wn(a)](f ) = e−jπf N Te (7)
N +1 sin πf Te
2. Expliquez en quoi le résultat trouvé se distingue d’un sinus cardinal.
(b) (a) d (a) (b)
3. On considère wn = wn ∗ wn . Calculez wn . Pour cela on pourra d’abord
considérer le cas N = 2, N = 3 puis généraliser en considérant séparément le cas
n ≤ N et le cas n ≥ N .
(b)
4. Comment appele-t-on la fenêtre associée à wn ?
(b) (a) d (a)
wn
5. A partir de la définition de = wn ∗ wn et à partir (7), calculez le module de
(b)
la transformée de Fourier de wn .
6. Indiquez à partir de la figure 7 à quelle courbe 1, 2 ou 3 correspondent les modules
(a) (b)
des transformées de Fourier de wn pour N = 10 et N = 20 et wn pour N = 10 ?
Justifiez votre réponse ?
7. La figure 8 représente les modules des transformées de Fourier de xn = cos(2πf0nTe),
(a) (b)
xnwn et xnwn , avec f0 = 240Hz. Ecrivez (sans faire de véritables calculs) les
modules des transformées de Fourier de ces trois signaux. Indiquez à partir de la
figure 8 à quelle courbe 4,5,6 correspondent les transformées de Fourier de ces
signaux ? Justifiez votre réponse ?
(b) (a)
8. Commentez sur l’intérêt d’utiliser la fenêtre wn par rapport à la fenêtre wn .
31
11 Cours H : Synthèse de filtre AR, gabarit, transformées
bilinéaires
11.1 Exercices d’applications
Exercice 40 (22)
On cherche à synthétiser avec un filtre numérique de Butterworth un passe-bas de
fréquence de coupure fc# = fe/6, avec fe = 1kHz à l’ordre 3.
1. On cherche le gabarit que devrait avoir le filtre analogique pour qu’après ap-
plication de la transformée bilinéaire le filtre numérique obtenu soit celui sou-
haité. L’ingénieur A dit que les fréquences sont transformées par une applica-
tion non-linéaire qui est une tangeante, il propose f # = a tan(bf ), où f désigne
les fréquences du filtre analogique et f # désigne les fréquences du filtre numé-
rique ; quand pour deux fréquences f et f # cette relation est vérifiée alors il y
a égalité entre Ĥ(f ) et Ĥ #(f ). Pourquoi cette relation, telle qu’elle est proposée
par l’ingénieur A n’est pas pertinente ? L’ingénieur B propose une autre rela-
tion f = a tan(bf #). Proposez un raisonnement permettant de fixer b ; puis un
deuxième raisonnement permettant de fixer a en fonction de b en utilisant le fait
qu’à basse fréquence on souhaite que le filtre analogique et le filtre numérique ait
le même comportement. Quel est alors le gabarit ?
2. Sachant que les fonctions de transfert des filtres de Butterworth vérifient H(p)H(−p) =
32
 1 p0
(−1) p2n +1
n montrez que le changement de variable p =
, permet d’approcher le
2πf0
gabarit recherché, calculez pour cela la réponse fréquentielle. Comment choisir
alors f0 ?
3. La table des filtres de Butterworth indique que pour n = 3, H(p) = (p+1)(p12+p+1) .
D’après la façon de construire ces filtres et par suite sans calcul, donnez les pôles
de ces filtres.
4. Calculez H a la fonction de transfert du filtre analogique qui s’approche du gabarit,
ainsi que la réponse fréquentielle.
5. Calculez H b la fonction de transfert du filtre numérique recherché (fonction de
transfert √
et réponse fréquentielle). Les calculs sont simplifiés en montrant d’abord
1−z −1
que p = 3 1+z −1 .
−1 −2 −3
1 + 3z + 3z + z
H b(z) = √ √ √ √
(5 3 + 7) − (7 3 + 3)z + (7 3 − 3)z + (−5 3 + 7)z −3
−1 −2

6. Sans utiliser ce dernier résultat, calculez la réponse fréquentielle du filtre numé-

rique recherché.

Exercice 41 (13)
On cherche à comparer la réponse impulsionnelle d’une filtre analogique avec la ré-
ponse impulsionnelle d’un filtre numérique. La transformée de Laplace du filtre analo-
p
gique étudié est H(p) = p+1
33
 1. Pourquoi H est-il stable ?
2. On considère que la fréquence d’échantillonnage est de f e = 1Hz. Montrez la
transformée bilinéaire H d de H peut se mettre sous la forme H d(z) = a + 1−czb −1 ,
avec a = 2, b = −4/3, c = 1/3.
3. Calculez la réponse impulsionnelle.
4. Du fait de l’utilisation de la transformée bilinéaire, il existe une relation entre
la réponse fréquentielle du filtre analogique et la réponse fréquentielle du filtre
numérique. Quelle est cette relation, redémontrez cette relation et en déduire le
module de la réponse fréquentielle du filtre numérique ? Représentez le module de
la réponse fréquentielle sur [−2fe, 2fe].

Exercice 42 (16)
1
On considère un filtre de transformée de Laplace H(p) = p+1 . L’unité de temps choisie
est la seconde.
1. En quoi cette précision sur l’unité donne un sens plus précis à ce qui précède ?
2. Montrez que la relation entrée-sortie associée est : y 0 + y = u.
3. Ce filtre est-il stable ?
4. On définit un filtre linéaire numérique par sa fonct1ion de transfert H #(z) =
−1
H( T2e 1−z
1+z −1
), où Te est la période d’échantillonnage égale à 1ms. Il s’agit de la
transformée bilinéaire. Quelle est la relation entre la réponse fréquentielle de ce
filre numérique et celle du filtre analogique précédent ?
34
 5. Calculez le filtre numérique associé.
6. Ce filtre est-il stable ?

11.2 Exercices pour approfondir

Exercice 43 (23)
On cherche à construire un filtre numérique passe-bande de fréquences de coupures
fe/8 et 3fe/8, avec une fréquence d’échantillonnage de fe = 1kHz en utilisant les filtres
de Butterworth à l’ordre 1, 2 ou 3. Le gabarit et les filtre numériques synthétisés sont
représentés à droite de la figure 9 (p. 37).
1. Dessinez le gabarit du filtre anal ogique que l’on recherche. On pourra√utiliser
√ √ fe ( 2−1)
que tan(π/8)
√
= 2 − 1 et que tan(3π/8) = 2 + 1. On pose f 1 = π et
fe ( 2+1)
f2 = π
2. On propose le changement de variable
√  0
√ 
f1 f2 p 2π f1f2
p= √ +
f2 − f1 2π f1f2 p0
A quoi doit servir un tel changement de variable ? Quel est√le changement de va-
riable en fréquence correspondant ? Montrez que si f 0 = f1f2 alors f vaut 0.
Montrez que si f 0 = f1 ou si f 0 = f2 alors |f | = 2π
1
. Montrez que si f 0 > f2 ou si
f 0 < f1 alors |f | > 2π
1
. Montrez pourquoi un tel changement de variable permet de
35
 transformer le gabarit d’un filtre de Butterworth en le gabarit du filtre analogique
recherché.
3. Montrez que le filtre numérique s’obtient avec le changement de variable p =
1+z −2
1−z −2
.
1
4. Sur les tables des filtres de Butterworth, on a H(p) = p+1 , H(p) = p2+√12p+1 et
H(p) = (p+1)(p12+p+1) . En déduire trois filtres numériques qui répondent au pro-
1−3z −2 +3z −4 −z −6
blème. Le troisième est H(z) = 6+2z −4
.

Exercice 44 (24)
On cherche à déterminer l’ordre minimal du filtre de Butterworth qui permet de res-
pecter le gabarit d’un filtrepasse-bas.
 On utilise pour le gain une échelle en dB inversée
définie par α = −10 log10 |Ĥ(f )|2 . Cette échelle est commune aux filtres analogiques
et numériques. L’objectif porte sur le filtre numérique à synthétiser, il est que α ∈ [0, α1]
en f # ≤ f1# et α > α2 en f # ≥ f2#.
1. Pourquoi souhaite-t-on que le filtre numérique synthétisé vérifie α = 0 en la fré-
quence nulle ? Pourquoi α ≥ 0 ?
2. Déterminez le gabarit du filtre analogique correspondant, exprimez f1 et f2 en
fonction de f1# et f2#.

36
F IGURE 9 – A gauche : gabarit d’un passe-bas et filtre de Butterworth analogique. A
droite : gabarit du passe-bande et filtres numériques synthétisés en utilisant les filtres de
Butterworth à l’ordre 1, 2 et 3 (exercice 43, (23))

37
 p0
3. On considère le changement de variable p = 2πf0 ,
montrez qu’alors les filtres de
Butterworth se transforment en des filtres de réponses fréquentielles |Ĥ 0(f 0)|2 =
1
 0 2n . Pour cette question et les suivantes, on note p et f pour les variables
f
f0 +1
associées aux filtres de Butterworth définis dans les tables et on note p0 et f 0 pour
les variables associés aux filtres de Butterworth transformés.
4. Montrez que si f 0 = f0 alors α = 3. Les réponses fréquentielles dépendent de f0
et de n, montrez comment f0 et n modifient ces réponses fréquentielles.
5. Montrez que
α2
 
10 10 −1
log10 α1
1 10 10 −1
n=  
2 log10 f2
f1

6. Montrez que par substitution on peut ensuite trouver f0 :

α1
f0 = f1(10 10 − 1)−1/2n

7. Montrez que l’on peut exprimer l’ordre du filtre directement à partir des fréquences

38
du gabarit du filtre numérique.
α2
 
10 10 −1
log10 α2
1 10 10 −1
n=  
2 #
tan(f2 π/fe )
log10 #
tan(f1 π/fe )

39