Travaux Dirigés de
Traitement du Signal
Sujet n°1 : "Etude temporelle des signaux"
- corrigé partiel : exercices 5 et 6 -
Exercice 5 : Réponse d’un système linéaire et invariant à une entrée quelconque
On considère un circuit RC dont la réponse impulsionnelle est définie par :
1 t
h (t ) = exp( − )
RC RC
1) Représenter cette réponse. (*)
2) Déterminer sa réponse indicielle (réponse à un signal échelon de Heaviside) avec t≥0 (0
pour t<0) et la représenter. (*)
3) En déduire sa réponse à un signal porte défini par :
T
x(t) = rectT(t − ) (**)
2
avec T>>RC.
4) Même chose pour un signal rampe définie par x(t)=t entre 0 et 1, x(t)=0 ailleurs. (**)
Rappels de cours
Dans le domaine temporel, la relation entrée-sortie d’un système linéaire et invariant
(SLI) est un produit de convolution entre l’entrée et la réponse impulsionnelle :
+∞ +∞
s(t) = e(t) * h(t) = ∫ h(t − τ)e(τ)dτ = ∫ h(τ)e(t − τ)dτ
−∞ −∞
Dans le cas d’un système causal et d’un signal d’entrée également causal, le domaine
d’intégration se restreint à [0, t] :
t t
s(t) = ∫ h(t − τ)e(τ)dτ = ∫ h(τ)e(t − τ)dτ
0 0
Pour la borne inférieure (0), l’explication est simple : le signal d’entrée e(t) étant causal,
il n’est pas défini pour t<0, ce qui provoque une intégration nulle sur ]-∞, t].
Pour la borne supérieure, il faut bien comprendre ce que représente physiquement
l’opération de convolution sur les signaux. Prenons l’exemple d’un circuit RC série, dont la
réponse impulsionnelle est :
1/RC
h(t)=(1/RC)e –t /RC
t
La variable d’intégration étant ici τ (différente de t, puisque la valeur de sortie à chaque
instant t est le résultat d’une intégrale), il faut considérer h(τ) :
h(τ)
1 τ
On considère un décalage vers la droite, égal à t (>0) :
h(τ-t)
τ
t
On effectue un changement de variable τ’=τ-t :
h(τ’)
τ'
-t 0
On inverse la variable :
h(-τ’)
τ
-t 0
On fait ré-apparaître la variable initiale (τ) :
h(t-τ)
τ
0 t
Cette fonction est alors multipliée par le signal d’entrée (on prend l’exemple d’un échelon
unité) :
e(τ)
1
τ
0
La multiplication donne :
h(t-τ)
τ
t
2
C’est cette fonction qui est intégrée. On voit que du fait que la réponse impulsionnelle est
causale, l’intégration sera nulle sur l’intervalle ]t,+∞[. On peut donc limiter l’intervalle
d’intégration entre 0 et t.
L’intégration consiste à mesurer la surface entre la courbe et l’axe des abscisses, ce qui
donne un point (pour t) du signal de sortie :
τ
t
τ
t
En appliquant cette procédure pour toutes les valeurs de t>0, on obtient la réponse du
système à un échelon, appelée réponse indicielle :
Solution
1)
1/RC
h(t)=(1/RC)e –t /RC
t
2) Pour obtenir la réponse indicielle, on utilise le produit de convolution comme pour
n’importe quel signal. Pour un signal d’entrée e(t), cette réponse est définie par :
+∞ +∞
s(t) = (e * h)(t) = ∫ h(t − τ)e(τ)dτ = ∫ h(τ)e(t − τ)dτ
−∞ −∞
où h(t) est la réponse impulsionnelle du système, définie par :
+∞ +∞
h(t) = ∫ δ(t − τ)e(τ)dτ = ∫ e(τ)δ(t − τ)dτ
−∞ −∞
Pour e(t)=u(t), on a :
∞
s(t) = ∫ u(τ)h(t − τ)dτ
−∞
1 t −( t −τ ) / RC
RC ∫0
= e dτ
3
1 t −t / RC τ / RC
RC ∫0
= e e dτ
e − t / RC t τ / RC
RC ∫0
= e dτ
[
= e − t / RC e τ / RC 0 ]
t
[
= e − t / RC e t / RC − 1 ]
− t / RC
= 1− e
La représentation graphique de cette fonction est la suivante :
3) Il faut distinguer plusieurs cas :
- 1er cas : t<0 ; le produit entre h(t-τ) et e(τ) est nul, donc l’intégrale est nulle ainsi que le
signal de sortie (comme précédemment).
- 2e cas : 0<t<T ; on retombe dans le cas où le signal d’entrée est l’échelon de Heaviside (voir
ci-dessus), donc :
y(t) = 1 − e−t / RC
- 3e cas : t>T ; il est alors plus simple d’effectuer un décalage temporel avant le calcul : t’=t-T.
Pour simplifier, on peut considérer que la valeur du signal à T est égale à 1. Le signal d’entrée
est donc égal à 1-u(t) :
1-u(τ)
h(t-τ)
τ
0 t
Le cas qui nous intéresse ici est alors τ >0. Le produit des 2 fonctions 1-u(t) et h(t- τ) étant nul
sur cet intervalle, il ne faut plus considérer les signaux comme causals (en effet, si on le faisait
le produit de convolution serait nul sur cet intervalle). L’intervalle d’intégration, initialement
]-∞,+∞[., se limite à ]-∞,0[. On a donc :
∞
s(t) = ∫ (1 − u(τ))h(t − τ)dτ
−∞
0
= ∫ (1 − u(τ))h(t − τ)dτ
−∞
1 0 − (t − τ) / RC
RC ∫− ∞
= e dτ
[
= e−t / RC eτ / RC − ∞ ]
0
[
= e−t / RC e0 − e−∞ ]
4
= e−t / RC
On obtient une exponentielle décroissante, correspondant physiquement à la décharge du
condensateur dans la résistance.
Remarque : si le signal à T=0 n’avait pas été égal à 1 mais à K (K<1), cela aurait équivalent à
avoir un signal d’entrée égal à K(1-u(t)), et le résultat aurait été :
y(t) = Ke−t / RC
4)
-1er cas : identique au cas du signal porte.
-2e cas :
t t
y( t ) = ∫ r (τ)h ( t − τ)dτ = ∫ τ.h ( t − τ)dτ
0 0
t t t
1 1 − t / RC
y( t ) = ∫ τ.h ( t − τ)dτ =
∫ τ.e −( t −τ ) / RC dτ = e ∫ τ.e τ / RC dτ
0
RC 0
RC 0
On doit utiliser l’intégration par parties :
(uv)' = u ' v + uv' ↔ uv = ∫ u ' v + ∫ uv' ↔ ∫ u' v = uv − ∫ uv'
En prenant u’=eτ/RC et v=τ, on a : u=RCeτ/RC et v’=1. Donc :
1 − t / RC − t / RC τ / RC t t τ / RC
t
y( t ) =
RC
e RC τe [
τ / RC t
0 − RC ]
∫0 e τ / RC
d τ =e τe [ 0 − ∫e ] dτ
0
{[ ] [ ]} [ ] [ ]
= e − t / RC τe τ / RC 0 − RC e τ / RC 0 = e − t / RC { te t / RC − RC e t / RC − 1 }
t t
= t − RC(1 − e − t / RC )
3e cas : identique à celui de la porte.
Exercice 6 : Convolution avec réponse impulsionnelle (***)
1) Le calcul de la moyenne mobile (ou glissante) d’un signal sur une durée T consiste à
calculer, à chaque instant t, l’intégrale d’un signal de t-T à T. Montrer que ce calcul est une
convolution avec la fonction rectangulaire centrée sur T/2, définie par :
1 T
h(t) = rectT(t − )
T 2
2) Estimer sans calcul le résultat de la convolution de ce signal avec le signal défini par :
e(t) = rect 2T(t − T)
puis le démontrer par le calcul.
Solution
1) La moyenne mobile est obtenue par :
t
1
s(t) = ∫ e(τ)dτ
T t −T
Cette expression est équivalente à une intégration sur l’intervalle ]-∞,+∞[, du signal e(τ)
multiplié avec une fonction porte. Il faut donc d’abord exprimer correctement cette dernière.
5
Une fonction porte de largeur T et centrée sur 0 est définie par (en fonction de la variable
d’intégration τ) :
rectT(τ)
La même fonction centrée sur t>0 est définie par :
rectT(τ − t)
Décalée de T vers la gauche, son expression devient :
T
rectT(τ − t + )
2
L’intégrale peut donc s’écrire :
+∞
1 T
s(t) = ∫ rectT(τ − t + )e(τ)dτ
T −∞ 2
Il reste à montrer que celle fonction porte est une réponse impulsionnelle. Pour cela on
cherche à identifier cette intégrale avec celle du produit de convolution :
+∞
s(t) = ∫ h(t − τ)e(τ)dτ
−∞
Par identification :
1
rect T(τ − t + T / 2) = h(t − τ)
T
soit, avec le changement de variable t’=t-τ :
1
h(t' ) = rect T(−t'+T / 2)
T
e
Avec un 2 changement de variable : t’’=t’-T/2 ;
1 1
h(t' ' ) = rectT(−t' ' ) = rect T(t' ' )
T T
car la fonction est paire.
Si l’on effectue le changement de variable inverse (t’=t’’+T/2) pour revenir à t’ :
1
h(t' ) = rect T(t'−T / 2)
T
que l’on peut noter également :
1
h(t) = rect T (t − T / 2)
T
on obtient bien une fonction causale puisque de largeur T et décalée vers la droite de T/2
(donc nulle pour t<0).
2) Il s’agit d’une fonction porte de hauteur 1, centrée sur T et de largeur 2T. Puisque la
convolution avec la fonction h(t) précédente constitue un calcul de moyenne mobile, on peut
déduire la forme du signal de sortie sans le calculer :
- de -∞ à 0 ; le produit entre h(t-τ) et e(τ) est nul, donc l’intégrale est nulle ainsi que le signal
de sortie ;
- de 0 à T ; la moyenne du signal d’entrée donne une droite de pente 1, donc s(T)=1 ;
- de T à 2T ; la moyenne du signal d’entrée est constante, égale à 1 ;
- de 2T à 3T ; la moyenne du signal d’entrée est une droite de pente –1, donc s(3T)=0 ;
- de 3T à +∞ ; la moyenne du signal d’entrée est nulle.
rect2T(t-T) (1/T)rectT(t-T/2) s(t)
1 1
* 1/T =
0 2T 0 T 6 0 T 2T 3T t
