Cours et Travaux Dirigs de

Traitement du Signal Dterministe

Benot Decoux ([email protected])

- Exercices -

1re partie : "Notions de base et tudes temporelles"

1
Bases du traitement de signal

Exercice

Calculer lamplitude de la drive dun signal sinusodal damplitude gale 1 et de frquence 2

Hertz.

Rponse

La drive du signal sinusodal dfini par exemple par :

s( t ) = A cos(t + )
est dfinie par :
s( t ) = A sin(t + )
donc lamplitude du signal driv est A. Lapplication numrique donne :
A = 2 2 = 4

Exercice

Exprimer la fonction chelon unit sous forme dune fonction signe damplitude judicieusement
choisie et dune constante.

Rponse
1 1
u (t ) = + sgn( t )
2 2

Exercice

Exprimer la fonction rectangulaire x ( t ) = A.rect[t T ] laide de 2 signaux chelons.

Rponse
x ( t ) = A.u ( t + T / 2) A.u ( t T / 2)

Exercice

1) Calculer la valeur moyenne et la valeur efficace dun signal carr, compris entre 0 et 5V, de rapport
cyclique 1/2.

2) Mme chose pour un rapport cyclique 1/3.

3) Calculer la valeur moyenne dun signal sinusodal damplitude A, dfini par :

s( t ) = A cos(t + )

4) Calculer la valeur efficace de ce signal.

Solutions

1) Soit s(t) ce signal. Comme il est priodique, sa valeur moyenne est dfinie par :
T T/2 T/2
1 1 1 5 T/2 5 T
Smoy = s( t )dt = s( t )dt = 5dt = [t ]0 = = 2,5V
T0 T 0 T 0 T T 2
Sa valeur efficace est dfinie par :

2
 T T/2 T/2
1 1 1 25 25 T
= s 2 ( t )dt = 0 s ( t)dt = T 25dt = T [t ]
2 2 T/2
S eff 0 = = 12,5V 2
T0 T 0
T 2
Soit
Seff 3,5V
2) Valeur moyenne :
T T/3 T/3
1 1 1 5 5 T
= s( t )dt = 0 s( t )dt = T 5dt = T [t ]
T/3
Smoy 0 = 1,66V
T0 T 0
T 3
Valeur efficace :
T T /3 T/3
1 1 1 25 25 T
= s 2 ( t )dt = 0 s ( t )dt = T 25dt = T [t ]
2 2 T /3
S eff 0 = = 8,33V 2
T0 T 0
T 3
Soit
Seff 2,9V
3)
T
1 t 0 +T 1 T 1 A
S moy = A cos(t + )dt = A cos(t + )dt = sin(t + )
T 0t T 0 T 0
A A A
= [sin(T + ) sin(0 + )] = [sin( 2 + ) sin()] = [sin() sin()] = 0
T T T

4)
1 T 2 A2 T

T 0
2 2
Seff = A cos ( t + ) dt = cos 2 (t + )dt
T 0

On utilise la formule de trigonomtrie :

1
cos 2 a = (1 + cos 2a )
2
do

{ dt + cos(2t + )dt}= A2T [t] + sin(22t + )

T
A2 T A2 T T
2

0 1 + cos(2t + )dt = 2T
2 T
S eff = 0
2T 0 0
0
A
2
sin( 2T + ) sin( ) A sin( ) sin() A
2 2
= T + = T + = 2
2T 2 2T 2
Soit :
A
Seff =
2
Les lectroniciens connaissent bien ce rsultat.

Exercice

Soit x(t) un signal carr logique TTL (tat bas : 0V ; tat haut : 5V) de rapport cyclique 1/2 et de priode
T=0,1s.
1) Calculer son nergie sur une priode. En dduire son nergie totale.
2) Calculer sa puissance totale et sa puissance moyenne.
3) En dduire sa valeur efficace.

Rponses

1) Son nergie sur une priode est dfinie par :

T T/2 T/2
25 T
E T = x 2 ( t )dt = x ( t )dt = 25dt = 25[t ]
2 T/2
0 = = 12,5 T = 1,25 Joule
0 0 0
2
3
Son nergie totale est gale :

x ( t )dt = 25[t ] = 25[ + ] =

2
ET =

2) La puissance moyenne totale est identique la puissance calcule sur une priode, dfinie par :
T T/2 T/2
1 1 1 25 T / 2 25
PT = x 2 ( t )dt = x 2 ( t )dt = 25dt = [t ]0 = = 12,5 Watts
T0 T 0 T 0 T 2
3) La valeur efficace est la racine carre de la puissance (calcule sur une priode, ou totale) :
X eff = 12,5 = 3,53 Volt

Exercice

Calculer lnergie et la puissance totales des signaux suivants (on prendra T=1 quand ncessaire pour
les applications numriques) :
Echelon de Heaviside
Fonction porte de largeur T et de hauteur 1/T, centre sur 0

Rponse

1) Echelon de Heaviside.

Energie :
+ + +
E = s ( t )dt = s ( t )dt = 1.dt = [t ]0 = + 0 = +
2 2 +

0 0
Puissance totale :
T/2
1 1 T/2 1 T 1
s 2 ( t )dt = lim [t ]0 = lim = Watt
T T
P = lim
T T T T 2
T / 2 2

2) Fonction porte de largeur T et de hauteur 1/T, centre sur 0.

Energie :
+ T/2 T/2
1 1 T/2 1 T T
E = s ( t )dt = s ( t )dt =
2 2
1.dt = [t ]T / 2 = ( + ) = 1 Joule
T / 2
T T / 2 T T 2 2
Puissance totale :
E
lim = 0
T T

Convolution-Rponse impulsionnelle

Exercice

On considre le produit de convolution entre 2 signaux x(t) et y() :

+
x ( t ) * y( t ) = x ().y( t ).d

Par un changement de variable adquat, montrer que le produit de convolution est commutatif.

Solution

On cherche dmontrer que :

4
 x ( t ) * y( t ) = y ( t ) * x ( t )
Appelons s(t)=x(t)*y(t). Si lon effectue le changement de variable =t- , on obtient :

s( t ) = x ( t ' ).y(' ).d'
+
soit
+
s( t ) = x ( t ' ).y(' ).d'

que lon peut r-crire
+
s( t ) = x ( t ).y().d = y( t ) * x ( t )

Ce qui dmontre que le produit de convolution est commutatif.

Exercice

1) Simplifier les intgrales suivantes :

s( t )(t )dt ; s( t )( t + 1)dt

o s(t) est un signal quelconque, causal puis non causal.
2) Calculer la valeur numrique des intgrales suivantes :

r( t )( t 1)dt o r(t) est la fonction rampe

0

Solution
1)

s(t )( t )dt = s(0)( t)dt = s(0) (t )dt =s(0)

De mme :

s( t )( t + 1)dt =s(1)

2)

r( t )( t 1)dt = r(1)(t 1)dt = r (1) (t 1)dt = r(1)

0 0 0

Exercice

Montrer que la convolution dun signal e(t) avec la fonction rectangle dfinie par :
1 t t0
h ( t ) = rect
T T
(centre sur t0, damplitude 1/T et de largeur T), avec t0=-T/2, correspond un filtrage de type
moyenneur.

Solution

La dfinition de la convolution donne :

+ + t
1 t T / 2 1
s( t ) = h ( t )e()d = rect e()d = e()d

T T T t T
qui est la dfinition de la moyenne mobile.

5
Exercice

1) Montrer que lopration de moyenne mobile (ou glissante) est une convolution avec la fonction
rectangulaire.
2) Exprimer la rponse impulsionnelle correspondante.

Solution

1)
t + +
1 1 t T / 2
s( t ) = e()d = rect e()d = h ( t )e()d
T t T T T
2)
1 t + T / 2
h (t ) = rect
T T

Exercice

1) Dterminer la rponse indicielle (rponse un signal chelon de Heaviside) dun circuit RC dont la
rponse impulsionnelle est dfinie par :
1 t
h(t) = exp
RC RC
avec t0 (0 pour t<0).

2) Reprsenter cette rponse impulsionnelle ainsi que la rponse du circuit.

Solution

1) Cette rponse est dfinie par :

S( u ( t )) = u ( ) h ( t )d

1 t ( t ) / RC
RC 0
= e d
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

2)
1/RC 1

h(t)=(1/RC)e t /RC
t t

6
Exercice

Calculer la rponse dun circuit RC une rampe de pente 1, partir de sa rponse impulsionnelle.

Solution
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 lintgration 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 e
/ RC
d = e e [
0 e ]
d
0 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 )

Exercice

On applique lentre dun filtre passe-bas une impulsion damplitude 10V et de dure 0,00001s. Sa
rponse observe en sortie est dfinie par y(t)=e-3000t.

1) Calculer laire dfinie par limpulsion dentre et laxe des abscisses. Lexprimer en fonction de
limpulsion de Dirac sous la forme c.(t).
2) Reprsenter y(t), en prcisant sa pente lorigine.
3) En dduire la "vraie" rponse impulsionnelle du systme, que lon notera h(t).
4) Dterminer lexpression de la rponse de ce systme une entre chelon unit.
5) Mme question pour un signal rampe de pente 1 pour t0 et nul pour t<0.

Solution

1) La produit de la dure de limpulsion dentre et son amplitude est gal 10-510=10-4V.s Lentre
peut tre assimile une impulsion de Dirac pondre : 10-4 (t).

2) y(t) est une exponentielle commenant au point (0,1) et dcroissant avec une pente initiale de 3000.

3) Puisque le signal dentre ne reprsente quune fraction de limpulsion de Dirac, on peut considrer
que la vraie rponse impulsionnelle sobtient en pondrant la rponse limpulsion donne dans lnonc
de la manire suivante :
10-4 (t) e-3000t
(t) 104e-3000t =h(t)

4) La rponse du filtre tout autre entre est obtenue par la convolution entre cette entre et sa rponse
impulsionnelle. La rponse lchelon u(t) est donc donne par :
t
y( t ) = u ()h ( t )d
0

7
 t t
10 4 3000 t
= h ( t )d = 10 4
e d =
3000
3000
e 0 =
3
e [3
]
e = (1 e 3000 t )
10 3000 t 0 10
[ ]
0 0
5) Soit r(t) lexpression de cette rampe. On remplace lexpression de u(t) par celle de r(t) dans lintgrale
de convolution :

y( t ) = r ()h ( t )d = .h ( t )d
0
On restreint lintervalle dintgration [0,t] :
t t t
y( t ) = .h ( t )d = 10 4 .e 3000 ( t ) d =10 4 e 3000 t .e 3000 d
0 0 0
On doit utiliser lintgration par parties :
(uv)' = u ' v + uv' uv = u ' v + uv' u' v = uv uv'
-3000t
En prenant u=e et v=, il vient :
1 t
10 4 e 3000 t 3000 t t 3000
y( t ) = 10 4 e 3000 t
3000
[
t
e 3000 0 ]
1
3000 e 3000
d =
3000
e [0 e ] d
0 0
10 4 e 3000 t 3000 t 10 4 e 3000 t 3000 t
=
3000
[
e 0 +
1
3000
] e [
3000 t
0 = ]3000
te [ ]
+
1
3000
[
e3000 t 1 ]

10 1 1 3000 t
= t + e
3 900 900

Exercice

1) Soient e et h deux squences de valeurs provenant respectivement de lchantillonnage dun signal et

de la rponse impulsionnelle dun systme :
e={0,0,0,1,1,1,1,0,0,0} et h={1,-1}
Calculer la squence s rsultant de la convolution numrique e*h.
2) Interprter ces rsultats du point de vue des plages de frquences limines et conserves.
3) Quel est le signal qui permettrait de connatre la squence h ?
4) Proposer une squence h permettant de raliser un moyennage du signal e. Mme question quen 2)
5) Dmontrer laide de cet exemple que le produit de convolution est bien commutatif.

Rponse

1) On utilise la formule de la convolution discrte :

N 1
s k = e k * h k = e k i .h i avec k=0,1,,M+N-2
i =0
N est le nombre dlments de h et M celui de e : N=2 et M=10. Elle possde un nombre dlments gal
la somme de ceux des deux squences, moins 1, soit 10+2-1=11. On obtient la squence suivante :
s={0,0,0,1,0,0,0,-1,0,0,0}
2) On peut dire que les basses frquences sont limines puisque les longues suites de 1, assimilables
des frquences basses, sont limines.
3) La squence h correspond la rponse impulsionnelle dun systme. Il suffit donc dappliquer une
impulsion de Kronecker (lquivalent numrique de limpulsion de Dirac) en entre de ce systme. La
squence utiliser est :
e={1,0,0,0,0,0,0,0,0,0}
On peut facilement dmontrer que le rsultat de la convolution entre e et h est s={1,-1,0,0,0,0,0,0,0,0,0},
donc h={1,-1} est bien la rponse impulsionnelle du systme.
Ce rsultat dmontre galement que limpulsion de Kronecker est llement neurtre de la convolution.

8
4) Par exemple h={1/3,1/3,1/3}. On obtient :
s={0,0,0,1/3,2/3,1,1,2/3,1/3,0,0,0}
Ici ce sont plutt les hautes frquences qui sont limines.
5) On calcule e k * h k . Cest dire que concrtement, on retourne la squence {en} et on la dcale. On
constate que lon obtient le mme rsultat, ce qui vrifie que le produit de convolution est bien
commutatif.

Exercice

Soit limage suivante (les cases vide contiennent des 0) :

1 1 1
1 1 1
1 1 1
1 1 1

et le filtre de Prewitt suivant :

1 0 -1
1 0 -1
1 0 -1

On rappelle lexpression de la convolution bi-dimensionnelle discrte :

M / 2 1 M / 2 1
s(i, j) = h (k + M / 2, l + M / 2).e(i + k, j + l) , i,j=0,N-1, k entier
k = M / 2 l= M / 2
o N est la taille de limage en pixels par ct (carre), et M celle du filtre.
1) Donner la valeur du produit de convolution pour le pixel de coordonnes (2,2).
2) Calculer limage rsultant de la convolution de toute limage avec le filtre.

Solution

1)
1 1
s(2,2) = h (k + 1, l + 1).e(2 + k,2 + l)
k = 1l= 1
= h (0,0).e(1,1) + h (1,0).e(2,1) + h (2,0).e(3,1)
+ h (0,1).e(1,1) + h (1,1).e(2,2) + h (2,1).e(3,2)
+ h (0,2).e(1,2) + h (1,2).e(2,3) + h (2,3).e(3,3)
= 2

2) Le rsultat pour toute limage est :

9
 ? ? ?
? 0 -1 0 1
? 0 -2 0 2
0 -3 0 3
0 -3 0 3

Corrlation

Exercice

1) Calculer la fonction dautocorrlation dun signal porte dfini par

x ( t ) = rect[( t T / 2) T ]
2) Conclure sur les proprits de la corrlation utiles pour la mesure de ressemblance.

3) Dterminer lnergie de ce signal partir de sa fonction dautocorrlation.

4) Calculer la fonction dautocorrlation du signal carr obtenu par rptition de la fonction porte tous
les instants kT, avec T=2T et k entier.

5) Montrer quand dans ce dernier cas, la borne infrieure de lintgrale peut tre diffrente de la valeur
choisie dans la question prcdente.

Solution

1) La fonction possde une largeur gale T, une amplitude gale 1 et est centre sur T/2.

La simplification de cette intgrale va se ramener dterminer les bornes dintgration, selon diffrents
cas, car la fonction rect va tre remplace par 1. On peut distinguer 4 cas :
- 1er cas : si T+<0, soit <-T, le produit des 2 fonctions est nul, donc la fonction dautocorrlation
galement.
- 2e cas : -T<<0 : lintervalle de valeurs de pour lequel le produit nest pas nul est [0, +T]. Sur cet
invervalle, ce produit vaut 1. On a donc :
+ T
C xx () = 1.dt = [t ]0
+ T
= +T
t =0
- 3e cas : 0<<T : lintervalle sur lequel le produit nest pas nul est [, T]. Sur cet invervalle, ce produit
vaut 1. On a donc :
C xx () = 1.dt = [t ] = T
T T
t =
- 4e cas : >T :le produit des fonctions est nouveau nul.

2) La fonction dautocorrlation est maximale quand la concidence entre le signal et lui-mme est
maximale. Elle traduit donc la ressemblance entre les 2 signaux, permettant de dterminer le dcalage
pour lequel cette ressemblance est maximale.

3) On utilise directement la proprit selon laquelle la fonction dautocorrlation en 0 est gale lnergie
du signal, soit ici T.
10
4) Pour un signal priodique, lexpression de la fonction dautocorrlation est :
1 T
C xx () = x ( t ).x ( t )dt
T t =0
Ici, T reprsente la priode du signal, gale 2 fois la dure de la fonction porte. On prendra T pour ne
pas confondre avec le T dsignant la dure du signal porte.
1 T'
C xx () = x ( t ).x ( t )dt
T ' t =0
Pour dterminer cette fonction, on peut utiliser directement les rsultats du signal porte. En effet, la
fonction obtenue tait comprise entre T et +T, et nulle en dehors de cet intervalle. Ici, la mme fonction
va rapparatre aprs une priode du signal et ainsi de suite, le dcalage augmentant ( ) indfiniment. De
mme pour <0. La fonction dautocorrlation est donc pdiodique, de priode gale celle du signal.

5) On peut remarquer graphiquement que si lon avait effectu lintgration sur lintervalle [-T/2, T/2]
plutt que [0,T], le rsultat aurait t le mme.

Exercice

Soit un signal dfini par :

A
x(t) = t[u ( t ) u ( t T )]
T
1) Reprsenter ce signal.

2) Calculer sa fonction dautocorrlation, et la reprsenter.

3) Dterminer son nergie partir de sa fonction dautocorrlation.

Solution

1) On peut galement exprimer le signal sous la forme :

A
t pour t [0, T]
x(t) = T
0 ailleurs
La forme de ce signal est une dent de scie de largeur T : en t=0, sa valeur est 0 ; en T, elle vaut A.

2) Sa fonction dautocorrlation est dfinie par :

+ + A A
C xx () = x ( t ).x ( t )dt = t.[u ( t ) u ( t T)]. ( t ).[u ( t ) u ( t T)]dt
t = t = T T
Les diffrents cas tudier sont les mmes quavec un signal porte. Les 2 seuls cas pour lesquels le
produit des fonctions nest pas nul sont :
- 2e cas : -T<<0 : lintervalle de valeurs de pour lequel le produit nest pas nul est [0, +T]. On a donc :
+T A A
C xx () = t. ( t ).dt
t =0 T T
+ T
A2 A2 t3 t2 A 2 ( + T)3 ( + T) 2
Tt=0
= 2 t
+ T
t (
.dt 2
= )
T 2 3 2 0
=
T 2 3

2


2 2
A
( A
) (
= 2 2( + T)3 3( + T) 2 = 2 2( + T)( 2 + 2T + T 2 ) 3( 2 + 2T + T 2 )
6T 6T
)
2
A
(
= 2 2(3 + 2 2 T + T 2 + 2 T + 2T 2 + T 3 ) 3(3 + 2 2 T + T 2 )
6T
)
11
 A2
= 2 (2(3 + 3 2 T + 3T 2 + T 3 ) 3(3 + 2 2 T + T 2 ) )
6T
A2
= 2 ( 3 + 3T 2 + 2T 3 )
6T
e
- 3 cas : 0<<T : lintervalle sur lequel le produit nest pas nul est [, T]. On a donc :
T A A
C xx () = t. ( t ).dt
t = T T
T
A2 t3 t2 A2 T3 T 2 3 3 A 2
= 2 = 2
T 3 2 T 3

2
+ =
3 2 6T 2
(
2T 3 3T 2 + 3 )
On remarque que lon a le mme rsultat que dans le cas prcdent, mais avec un changement de signe
pour . On peut donc crire ces 2 rsultats sous la forme dun seul :
A2
( 3
)
2 2T 3 3T 2 + si < 0
C xx () = 6T
0 si > 0
Cette fonction est un polynme dordre 3. Elle est symtrique par rapport 0. Quand 0 par valeurs
ngatives, la fonction se comporte comme a+b.

3) Lnergie du signal est la valeur de sa fonction dautocorrlation pour =0. On remplace donc par 0
A 2T
dans lexpression prcdente. On obtient C xx () = .
3

Exercice

Calculer la fonction dautocorrlation du signal sinusodal.

Solution

Dans le cas simplifi dune amplitude unit et dune phase nulle, la fonction sinusodale est dfinie par :
s( t ) = sin(0 t )
pour une pulsation 0.
On applique lexpression de lautocorrlation :
1 T 1 T
C xx () = x ( t ).x ( t )dt = sin(0 t ). sin(0 ( t ))dt
T t =0 T t =0
On utilise la formule de trigonomtrie :
cos(a b) cos(a + b)
sin(a ). sin( b) =
2
donc :
1 T 1 T
C xx () =
2T t = 0
cos(0 )dt
2T t =0
cos(0 (2 t ))dt
cos(0 ) T 1 T
C xx () =
2T 1dt
t =0 2T t =0
cos(0 (2 t ))dt
cos(0 ) 1
= + [sin(0 (2T )) sin(0 )]
2 2 0 T
cos(0 )
C xx () =
2

12
Exercice : Corrlation entre un signal long et dun signal court

1) On considre le signal long suivant, sous la forme dune squence {xn}, o n reprsente les indices des
lments dans la squence, commenant 0 :
signal long : 1 0 1 2 1 2 2 0 0
Et le signal court sous la forme dune squence {yn} :
signal court : 0 1 2
On cherche dtecter la prsence du signal court dans le signal long. Pour cela, on dfinit le produit de
corrlation de la manire suivante :
1 N1
C xy (k ) = x i + k .y i
N i =0
o N est le nombre dlments du signal court : N=3.

1) Calculer le produit de convolution pour k=0, k=1 et k=4.

2) Par interprtation de ces rsultats, indiquer si la dtection de la ressemblance entre les 2 signaux est
bien effectue de cette manire.

3) Recommencer avec les versions centres des signaux (pour chacun, la moyenne de ses lments est
retranche de chaque lment) :

signal long : 0 -1 0 1 0 1 1 -1 -1
signal long : -1 0 1

4) Mme question que 2) pour les rsultats du 3).

5) Conclure sur lutilit de centrer les signaux pour rechercher des motifs dans un signal par corrlation,
et r-crire lexpression de la fonction de corrlation prenant en compte ce centrage.

Solution

1) Cxy(0)=2 ; Cxy(1)=5 ; Cxy(4)=6

2) Ces calculs ne permettent pas de dtecter la ressemblance maximale, car la valeur maximale ne
correspond pas celle-ci.

3) Cxy(0)=0 ; Cxy(1)=2 ; Cxy(4)=1

4) Ici, la valeur maximale correspond la ressemblance maximale entre le signal court et la zone
compare du signal long.

5) Le centrage des signaux permet donc de transformer ce calcul en un moyen de dtecter la ressemblance
entre 2 signaux. On peut r-crire la fonction de corrlation de la faon suivante :
1 N 1
C xy (k ) = ( x i+ k x ).( y i y)
N i =0
o x et y reprsentent respectivement les moyennes des squences {xn} et {yn}.

13
Equations diffrentielles

Exercice

On considre lquation diffrentielle suivante :

ay'+ y = 1
1) Dterminer lensemble des solutions de lquation homogne associe.
2) Dterminer une solution particulire de lquation complte.
3) En dduire la solution gnrale de lquation complte.
4) En dduire la forme gnrale de la rponse dun circuit RC (=tension de sortie) un chelon de
tension E.
5) Prciser cette rponse sachant que le condensateur est dcharg t=0.

Indications

Pour la question 2, on pourra penser une solution constante.

Pour la question 5, il faudra utiliser une considration lectrique : la tension aux bornes dun
condensateur ne peut tre discontinue. Donc t=0, on aura s(0-)=s(0+).

Solution

1) Lensemble des solutions de lquation homogne est dfini par :

x

e a

rel quelconque.

2) Prenons une solution constante K. K=0 donc K=1.

3) La solution gnrale de lquation complte est donc :

x

y = 1 + e a

4) Cette solution nous permet de connatre la solution dans le cas dun circuit RC dfini par lquation :
ds( t )
RC + s( t ) = e( t )
dt
La solution gnrale de lquation sans second membre (=rgime libre) est donc :
t

e RC

Pour la solution gnrale de lquation complte, il faut prciser e(t). Si lon sintresse la rponse
indicielle (=rponse un chelon), lquation est dfinie par :
ds( t )
RC + s( t ) = E
dt
et on va chercher une solution particulire de cette quation sous la forme s(t)=cte. On a ds(t)/dt=0, donc
cette constante est E. La solution gnrale de lquation complte est donc :
t

s ( t ) = E + e RC

Il reste dterminer la constante , qui dpend des conditions initiales. Celles-ci sont dtermines par
des considrations lectriques : si on suppose le condensateur initialement dcharg, la tension ses
bornes est nulle. On a donc :
s(0-)=0
mais aussi, en remplaant t par 0 dans lexpression de la sortie :
s(0+)=E+
En posant :
14
 s(0-)=s(0+)
on obtient :
=-E
Donc, finalement :
t

s( t ) = E(1 e RC
)

Remarque : on aurait pu considrer dautres conditions initiales. Par exemple, avec C initialement
charg :
s(0-)=U0.
avec toujours
s(0+)=E+
on a
= U0-E
et finalement :
t

s ( t ) = E + ( U 0 E )e RC

Transforme de Laplace

Exercice

Calculer la Transforme de Laplace dun signal carr dfini par x(t)=A pour t dans [kT, kT+], k
entier.

Solution

On utilise la proprit ci-dessus. On calcule donc dabord la transforme dune priode. Soit le
rapport cyclique du signal, la transforme dune priode est dfinie par :

X 0 (p) = A.e pt .dt
0

=
A
p
[
1 e p ]
La transforme du signal carr est donc :
A 1 e p
X ( p) =
p 1 e pT

Exercice : Dtermination de la rponse impulsionnelle dun circuit RC

Dterminer la rponse impulsionnelle dun circuit RC. On supposera le condensateur initialement

dcharg.

Indications

On dtermine lquation diffrentielle du circuit, on en dduit sa transmittance de Laplace et donc sa

rponse impulsionnelle.

Solution

On a vu que lquation diffrentielle qui rgissait ce circuit tait :

15
 ds( t )
RC + s( t ) = e( t )
dt
Par application de la transforme de Laplace, on obtient :
RCpS(p) + S(p) = E(p)
car les conditions initiales (cest dire s(0)=0) sont nulles.
La transmittance du systme est donc :
S(p) 1
T ( p) = =
E(p) 1 + RCp
Par utilisation de la transforme
A
s( t ) = A.ea t L
S( p) =
p+a
(avec a>0) et en r-crivant :
S(p) 1 1
T ( p) = =
E(p) RC 1
p+
RC
on en dduit la rponse impulsionnelle du systme :
1 t / RC
s( t ) = .e
RC

Exercice

Soit une fraction rationnelle dfinie par :

1
F(p) = 2
p + 3p + 2
1) Dterminer sa transforme de Laplace inverse.
2) En dduire la rponse impulsionnelle dun systme possdant F(p) pour transmittance.
3) Calculer la rponse de ce systme un signal chelon u(t).

Solution

1) On peut dmontrer facilement quon peut factoriser le dnominateur sous la forme :

1
F(p) =
(p + 1)(p + 2)
Par la mthode des rsidus, on la dcompose en lments simples :
1 1
F(p) = +
(p + 1) (p + 2)
et, par utilisation des transformes de Laplace lmentaires (voir table des transformes), on dduit
directement par transformation inverse :
f ( t ) = (e t e 2 t ).u ( t )
2) Si F(p) est la transmittance dun systme, f(t) est sa rponse impulsionnelle puisque la transforme
inverse dune transmittance est une rponse impulsionnelle.

3) La transforme dun chelon u(t) tant 1/p, on a :

S(p) 1
S(p) = F(p)E(p) = =
p p(p + 1)(p + 2)
En utilisant la dcomposition prcdente :
A A2 A3
F(p) = 1 + +
p (p + 1) (p + 2)
On calcule les coefficients Ai laide de la formule des rsidus :
16
 1
A1 = [F(p).p]p=0 =
2
A 2 = [F(p).(p + 1)]p=1 = 1
1
A 3 = [F(p).(p + 2)]p=2 =
2
On a donc :
0,5 1 0,5
F(p) = + +
p (p + 1) (p + 2)
do la rponse temporelle recherche :
f ( t ) = 0,5 e t + 0,5e 2 t

Exercice

Effectuer la dcomposition en lments simples de la fonction :

3p 3p
F(p) = =
(p + 1)(p + 2p + 1) (p + 1) 3
2

Solution

La dcomposition en lments simples donne :

A0 A1 A2
F(p) = 3
+ 2
+
(p + 1) (p + 1) (p + 1)
avec
1
[ ]
A 0 = F(p).(p + 1) 3 p=1 = [3p ]p=1 = 3
0!
1 d (F(p)(p + 1) 3 ) d (3p)
A1 = = = [3]p= 1 = 3
1! dp p=1 dp p=1
1 d 2 (F(p)(p + 1) 3 )
A2 = = [0]p=1 = 0
2! dp 2 p = 1
do :
3 3
F(p) = 2
+
(p + 1) (p + 1) 3

Exercice

On considre un circuit RLC dont le fonctionnement est rgit par lquation diffrentielle suivante :
d 2 s( t ) ds( t )
2
+ 2m0 + 02 s( t ) = e( t )
dt dt
avec
1 R C
0 = et m=
LC 2 L

On supposera les conditions initiales nulles (cest dire le condensateur initialement dcharg), ce
qui se traduit par :
ds(0 )
s (0 ) = =0
dt
Dterminer sa rponse impulsionnelle, avec les valeurs numriques suivantes :
17
 R=1k, C=100F, L=1mH

Solution

En utilisant la transforme de Laplace, lquation devient :

p 2S(p) + 2m0 pS(p) + 02S(p) = E(p)

[ ]
S(p) p 2 + 2m0 p + 02 = E(p)

S(p) 1
= 2
E (p) p + 2m0 p + 02

On factorise le dnominateur :
S(p) 1
=
E(p) (p p1 )(p p 2 )
puis on dcompose en lments simples :
S(p) A1 A2
= +
E (p) (p p1 ) (p p 2 )
avec
1
A1 = et A 2 = A1
(p1 p 2 )
Les ples dpendent du signe de , donc de la valeur de m :
R C 100 10 4
m= = = 50 0,01 = 50 0,1 = 5
2 L 2 10 2
m>1 donc les 2 ples sont rels, dfinis par :
p1 = 0 (m + m 2 1) et p 2 = 0 ( m m 2 1)
De mme,
1 1 1
0 = = = 3 = 1000rad / s
LC 2
10 10 4 10
donc les valeurs des ples sont :
p1 = 1000(5 + 24 ) 101 et p 2 = 1000(5 24 ) 990
On a donc :
1 1 1
A1 = 1,12 10 3 = A 2
(p1 p 2 ) 101 990 101 990
La rsultat
S(p) 1 1
= 1,12 10 3
E ( p) (p 101) (p 990)

La rponse temporelle correspondante est donc :

s( t ) = 1,12 10 3 (e 101t e 990 t )

18