Etude des circuits RLC

Projet de P6

Cedric Dangeard Sami Benjelloun Vincent Bachelot

Amer Tahle Clara Gainon de Forsan de Gabriac
Table des matières

I Présentation du projet 3
1 Introduction 4

2 Oragnisation du projet et répartition du travail 5

II Rapport 6
3 Les circuits RLC série 7
3.1 Charge d’un condensateur . . . . . . . . . . . . . . . . . . . . . . 7
3.1.1 Cas du régime pseudo-périodique (m < 1) . . . . . . . . . 8
3.1.2 Cas du régime apériodique (m > 1) . . . . . . . . . . . . . 9
3.1.3 Cas du régime critique (m = 1) . . . . . . . . . . . . . . . 10
3.2 Décharge d’un condensateur . . . . . . . . . . . . . . . . . . . . . 12
3.2.1 Cas du Régime pseudo périodique (m > 1) . . . . . . . . 12
3.2.2 Cas du régime apériodique (m > 1) . . . . . . . . . . . . . 13
3.2.3 Cas du régime critique, m = 0 . . . . . . . . . . . . . . . 13

4 Résolution d’équations de circuits avec la méthode de Laplace 14

4.1 Introduction à la méthode de Laplace . . . . . . . . . . . . . . . 14
4.2 Transformée de Laplace d’une fonction échelon dans un circuit RC 15
4.3 Résolution détaillée d’une équation différentielle du second ordre
d’un circuit RLC avec Ve(t) fonction échelon . . . . . . . . . . . 16

5 Moédlisation et application 21
5.1 Résolution d’équations de circuits avec la méthode de Laplace
sous mapple et affichage des graphiques associés . . . . . . . . . . 21
5.1.1 Réponse d’un circuit RC à un échelon de tension . . . . . 21
5.1.2 Réponse d’un circuit RLC à un échelon de tension . . . . 22
5.1.3 Réponse d’un circuit RC à un signal créneau . . . . . . . 23
5.1.4 Conclusion sur notre travail Maple . . . . . . . . . . . . . 24
5.2 Conception d’un récepteur radio AM grâce aux circuits RLC . . 25
5.2.1 L’émetteur : . . . . . . . . . . . . . . . . . . . . . . . . . . 25
5.2.2 Le récepteur : . . . . . . . . . . . . . . . . . . . . . . . . . 26

6 Conclusion 29

1
III Annexes 30
7 Méthode de Laplace 31
7.1 Résolution détaillée d’une équation différentielle du second ordre
d’un circuit RLC avec Ve(t) fonction d’impulsion unitaire . . . . 31
7.2 Résolution détaillée d’une équation différentielle du second ordre
d’un circuit RLC avec Ve(t) fonction sinusoı̈dale . . . . . . . . . 32

8 Méthode de Laplace appliquée sur maple (code mapple) 35

8.1 fonction échelon d’un circuit RC . . . . . . . . . . . . . . . . . . 35
8.2 fonction crénau d’un circuit RC . . . . . . . . . . . . . . . . . . . 38
8.3 fonction échelon d’un circuit RLC . . . . . . . . . . . . . . . . . . 49

2
 Première partie

Présentation du projet

3
Chapitre 1

Introduction

De nos jours, le circuit RLC est l’un des systèmes électriques les plus répandus
dans la vie quotidienne. On le retrouve en effet dans les radios de nos voitures,
dans la plupart de nos appareils électroménagers et même dans les disposi-
tifs antivols des magasins. Ce projet de physique nous offre donc l’opportu-
nité de maı̂triser un dispositif très courant dans l’industrie et que l’on sera
nécessairement amené à rencontrer au cours de notre carrière. La théorie oc-
cupe une place prépondérante dans cette étude. Le circuit RLC est en effet
régi par une équation différentielle générale que nous détaillerons par la suite.
Néanmoins, en fonction de la tension en entrée qui lui est appliquée, il en
résulte des réponses indicielles et donc des modes de fonctionnement radica-
lement distincts. C’est la raison pour laquelle l’un des objectifs de ce projet
fut de répertorier les différents régimes transitoires du circuit RLC. Ce travail
théorique s’est poursuivi ensuite sur le logiciel Maple. En nous appuyant sur
ce travail préliminaire, nous avons simulé le fonctionnement d’un circuit sous
plusieurs régimes et confronté la théorie avec des résultats obtenus lors d’une
séance de travaux pratiques réalisée en cours d’année dernière. Cependant pour
la résolution de ces équations différentielles, il est à préciser que le logiciel Maple
utilise une méthode qui nous était jusque-là totalement inconnue : la méthode
dite de Laplace. Se basant sur l’utilisation d’intégrales, cette dernière permet
d’extraire le problème du domaine du temps et simplifie par la même occa-
sion les calculs. C’est pourquoi un pan entier du projet consista en l’analyse et
l’explication de cette méthode efficace, de sorte que nous puissions comprendre
parfaitement les outils que nous avons manipulés et notamment Maple. Enfin
nous reviendrons sur l’aboutissement de ce projet : la réalisation d’un récepteur
radio. Il nous a permis de mettre en application et concrétiser toutes les connais-
sances acquises au cours de ce semestre.

4
Chapitre 2

Oragnisation du projet et
répartition du travail

5
Deuxième partie

Rapport

6
Chapitre 3

Les circuits RLC série

Dans cette partie, nous analyserons les différents comportements possibles

du circuit RLC en série. Nous étudierons le cas de la charge du condensateur,
c’est à dire quand celui-ci est soumis à un échelon de tension, de 0V à E puis
celui de la décharge du condensateur, pour un échelon de tension de E à 0V.
Rappelons au passage que c’est justement la charge du condensateur qui est
l’objet du TP, fil conducteur de ce projet. Ce travail préliminaire permettra
donc de démontrer les formules qui y figurent et ainsi de préparer au mieux la
simulation sous Maple.

3.1 Charge d’un condensateur

En utilisant la loi d’Ohm et la loi des mailles on obtient :

q(t) di
Rxi(t) + )+L× =E (3.1)
C dt

7
 Néanmoins, la formule suivante relie l’intensité du circuit à la charge et donc
à la tension aux bornes du condensateur :

dq(t) dUc
i(t) = =C× (3.2)
dt dt
Grâce à ce résultat que l’on  injecte  dans (1), on obtient l’équation
différentielle suivante qui régit le fonctionnement du circuit RLC en charge :

d2 Uc (t) R dUc (t)

 
Uc (t) E
+ + = (3.3)
dt L dt LC LC

Cette équation différentielle est du second ordre. La méthode de résolution

consistera, comme toujours, à résoudre l’équation différentielle sans second membre.
Uc0 (t) = E est une solution particulière évidente. La solution générale de
l’équation différentielle s’exprimera donc comme une somme de la solution à
l’équation différentielle sans second membre et de Uc0 . Ecrivons l’équation différentielle
sans second membre :

d2 Uc (t)
 
dUc (t)
+ 2λ + ω02 (Uc (t) = 0 (3.4)
dt dt
R 1
avec : λ = 2L et ω0 = √LC
On définit le coefficient
q d’amortissement m du système de la manière sui-
λ R C
vante : m = ω02
= 2 L
L’équation caractéristique est la suivante : X 2 + 2λX + ω02 X = 0, de discri-
minant ∆ = 4(λ2 − ω02 ). Trois comportements en régime transitoire découleront
du signe de ∆.

3.1.1 Cas du régime pseudo-périodique (m < 1)

C’est le cas dans lequel on a :
r
2 L
Delta < 0 ⇔ λ < ω02 ⇔R<2 ⇔m<1 (3.5)
C
(
s1 = −λ − jω
L’équation caractéristique admet deux racines complexes conjuguées :
s2 = −λ + jω
p
en posant ω = 2 − λ2
0

La solution de l’équation différentielle devient donc :

Uc (t) = E + M es1 t + N es2 t Uc (t) = E + e−λ (M ejωt + N e−jωt ) (3.6)

En revenant à la définition des exponentielles complexes, on peut aussi écrire

Uc (t) de la manière suivante :

Uc (t) = E + e−λt (M 0 cos(ωt) + N 0 sin(ωt)) ⇔ Uc (t) = −Be−λt (cos(ωt + φ) + E

(3.7)

8
 De plus : i(t) = C dUdt
c (t)
= −λBCe−λt cos(φ + ωt) − ωBCe−λt sin(φ + ωt)
En utilisant les conditions initiales, on obtient les deux équations suivantes :
i(t = 0) = 0 ⇔ −λBCcos(φ) − ωBCsin(φ) = 0 (3.8)

Uc (t = 0) = 0 ⇔ −Bcos(φ) + E = 0 (3.9)
q
λ2
p
Ainsi φ = tan− 1( ωλ et B = cos(φ)
E
= E 1 + (tan2 (φ) = E 1 + ω2 Dans

cette configuration, on parle de régime pseudo-périodique. En effet avec la

courbe de Uc (t) suivante on constate que l’allure sinusoı̈dale est modulée par un
terme d’exponentielle d’amortissement (ici e−λt ). L’amortissement des oscilla-
tions est d’autant plus faible que le coefficient est faible.

Ce signal présente une pseudo-période T, supérieure à la période propre du

circuit T0 = 2Π
ω0 :

2π 2π
T = =p 2 (3.10)
ω ω0 − λ2

3.1.2 Cas du régime apériodique (m > 1)

C’est le cas dans lequel on a :
r
2 L
∆>0⇔λ > ω02 ⇔R>2 ⇔m>1 (3.11)
C
L’équation caractéristique admet deux racines réelles négatives :
q q
s1 = λ + λ2 − ω02 s2 = λ − λ2 − ω02 (3.12)
Ainsi la solution à l’équation différentielle sans second membre est de la
forme :
Uc1 (t) = M.es1 .t + N.es2 .t (3.13)
Donc la solution générale est du type :
Uc (t) = Uc0 + Uc1 = E + M.es1 .t + N.es2 .t (3.14)

dUc (t)
Deplus : i(t) = C × = CM s1 × es1 .t + CN s2 × es2 .t (3.15)
dt

9
 On est en mesure de déterminer les constantes M et N ; en utilisant les
conditions initiales on obtient deux équations suivantes :

i(t = 0) = 0 ⇔ M s1 + N s2 = 0Uc (t = 0) = 0 ⇔ E + M + N = 0 (3.16)

d’où
E.s2 −E.s1
M= N= (3.17)
s1 − s2 s1 − s2

Ce régime transitoire est dit apériodique par opposition au précédent par

opposition au précédent. La tension au cours du temps s’exprimant comme une
somme d’exponentielles décroissantes, le signal ne présente aucune oscillation
et donc aucune période. Le régime permanent est atteint d’autant plus vite que
est grand. Le courant dans le circuit tend à s’annuler, tout comme les tensions
aux bornes de la bobine et de la resistance. Toute la tension de la source est
contenue à terme aux bornes du condensateur d’où :

lim Uc (t) = E (3.18)

t→∞

Figure 3.1 – Régime apériodique

3.1.3 Cas du régime critique (m = 1)

C’est le cas dans lequel on a :
r
2 L
∆=0⇔λ = ω02 ⇔R=2 ⇔m=1 (3.19)
C
L’équation caractéristique admet une racine réelle double :
R
s = −λ = − (3.20)
2L
Ainsi la solution de l’équation différentielle est de la forme :

Uc (t) = E + (M + N t)e−λt (3.21)

10
 Ainsi on obtient

i(t) = e−λt (−CM λ + CN − λN t) (3.22)

De nouveau, en réutilisant les conditions initiales on trouve :

M =E N = −λE (3.23)

L’expression de Uc(t) devient :

Uc (t) = −E(1 + λt)e−λt + E (3.24)

Ici il s’agit du régime critique. Uc(t) ne présente pas d’oscillation et ce

régime ne diffère pas qualitativement du régime du régime apériodique décrit
précédemment. Dans cette situation, le régime permanent est atteint le plus
rapidement.

Figure 3.2 – Régime critique

11
3.2 Décharge d’un condensateur
En utilisant la loi d’Ohm et la loi des mailles, on peut écrire :

q(t) di(t)
R × i(t) + +L× =0 (3.25)
C dt
Néanmoins, la formule suivante relie l’intensité du circuit à la charge et donc
à la tension aux bornes du condensateur :
dq(t) dUc (t)
i(t) = =C× (3.26)
dt dt
Grâce à ce résultat que l’on  injecte  dans (1), on obtient l’équation
différentielle suivante qui régit le fonctionnement du circuit RLC en décharge :

d2 Uc (t) R dUc (t) Uc (t)

+ × + =0 (3.27)
dt L dt LC
Ici, du fait de la tension nulle appliquée au circuit, le second membre est
nul. La solution de cette équation différentielle coı̈ncide donc avec la solution
de l’équation différentielle sans second membre. Pour chacun des régimes, on
réutilisera donc les résultats précédents. Les constantes seront déterminées à
partir des nouvelles conditions initiales (à savoir Uc(t=0)=E et i(t)=0A).

3.2.1 Cas du Régime pseudo périodique (m > 1)

La solution à l’équation différentielle est la suivante

Uc (t) = −Be−λt × cos(ωt + φ) (3.28)

On en déduit l’expression de l’intensité :

dUc (t)
i(t) = C = −λBCe−λt × cos(ωt + φ) − ωBCe−λt × sin(ωt + φ) (3.29)
dt
On détermine les constantes à l’aide des conditions initiales

i(t = 0) = 0 ⇔ −λBCcos(φ) − ωBCsin(φ) (3.30)

Uc (t = 0) = E ⇔ −Bcos(φ) = E (3.31)
Ainsi :
λ
φ = tan−1 (− ) (3.32)
ω
r
E p λ2
B=− = −E. 1 + (tanφ)2 = −E 1+ (3.33)
cos(φ) ω2
En décharge comme pour la charge, Uc(t) oscille autour de sa valeur initiale
(ici 0V). La pseudo-fréquence est toujours inférieure à la fréquence propre. Plus
le coefficient d’amortissement m est faible, plus les oscillations sont faibles.

12
3.2.2 Cas du régime apériodique (m > 1)
L’équation de la solution différentielle est la suivante :

Uc (t) = M es1 t + N es2 t (3.34)

avec q q
s1 = λ + λ2 − ω02 s2 = λ − λ2 − ω02 (3.35)

On en déduit l’expression de l’intensité :

duc (t)
i(t) = C × = CM s1 × es1 t + CN s2 × es2 t (3.36)
dt
On détermine les ocnstantes M et N d’après les conditions initiales

i(t = 0) = 0 ⇔ M s1 + N s2 = 0 (3.37)

Uc (t = 0) = E ⇔ M + N = E (3.38)
Ce qui donne :
E.s2 −E.s1
M= N= (3.39)
s1 − s2 s1 − s2

En décharge, le condensateur atteint la valeur 0, le tout sans oscillation. La

décroissance de E à 0V est d’autant plus rapide que est grand.Le courant dans
le circuit tend à s’annuler, tout comme les tensions aux bornes de la bobine et
de la resistance. Toute la tension de la source est contenue à terme aux bornes
du condensateur d’où :

lim Uc (t) = E (3.40)

t→∞

3.2.3 Cas du régime critique, m = 0

La solution de l’équation différentielle est de la forme :

Uc (t) = (M + N t)e−λt (3.41)

Ainsi on obtient

i(t) = e−λt (−CM λ + CN − λN t) (3.42)

On détermine les constantes à l’aide des conditions initiales : M = E et

N = Eλ
L’expression de Uc(t) devient :

Uc (t) = E(1 + λt)e−λt (3.43)

13
Chapitre 4

Résolution d’équations de
circuits avec la méthode de
Laplace

4.1 Introduction à la méthode de Laplace

La caractérisation des circuits eĺectriques se fait souvent par détermination
puis résolution d’équations différentielles. C’est la raison pour laquelle nous nous
sommes intéressé à une méthode de résolution, qui est la méthode de Laplace.Il
faut savoir que la méthode de Laplace permet de résoudre des équations intégro-
différentielle dans le domaine temporel.
Cette méthode est très efficace quant à la résolution des équations caractérisant
les réseaux, puisqu’elle consiste en la simplification d’ équations dans le domaine
temporel via les transformées de Laplace. On passe donc à une équation sans
dérivées partielles dans le domaine de Laplace puis à la résolution de l’équation
obtenu et enfin en appliquant l’inverse de la transformée de Laplace á la solu-
tion trouvée, pour revenir dans le domaine temporel, on retrouve la solution de
l’équation initiale.
La définition de la transformée de Laplace pour une fonction f(t) est la
suivante : Z ∞
F (p) = e−pt f (t)dt (4.1)
0

notée L(f(t))=F(p), avec p la variable de Laplace (p est un complexe).

Note : la transformée de Laplace existe si et seulement si l’intégrale impropre
(1) converge.
De même, on définie l’inverse de la transformée de Laplace de la façon
suivante :
Z γ+i∞
1
f (t) = L−1 (F (p)) = ept F (p)dp (4.2)
2Πi γ−i∞

où γ est choisie pour que l’intégrale soit convergente.

14
4.2 Transformée de Laplace d’une fonction échelon
dans un circuit RC
pour introduire la transoformé de Laplace, on va prendre un exemple assez
simple pour se familiariser avec cette nouvelle méthode. On a alors choisit la
fonction échelon unité d’un circuit RC que nous avons trouvés dans le TP réalisé
en premier cycle préparatoire. L’équation différentiel dans le domaine temporel
donnée dans le TP pour ce type de circuit était :

dVs (t)
τ + Vs (t) = Ve (t) (4.3)
dt
On
 pose donc les transformés de Laplace de notre équation :
L{Ve (t)} = Ve (p)

L{Vs (t)} = Vs (p)
 dVs (t)

L{ dt } = pVs (p)
puis, on obtient en les remplaçant dans l’équation une nouvelle expression
qui est dans le domaine de Laplace :

τ vs (p)p + vs (p) = ve (p) (4.4)

on pose désormais τ = RC et on isole vs puisqu’on veut exprimer la tension
de sortie qui est en faite la tension aux bornes du condensateur :

ve (p)
vs (p) = (4.5)
1 + RCp
Ensuite, on sait d’après le tableau des transformé de Laplace, que dans le
domaine de Laplace,comme ve est une fonction échellon :
1
ve (p) = (4.6)
p

en remplaçant ve(p) on trouve alors :

1
vs (p) = (4.7)
p(1 + RCp)

Ensuite il nous allons décomposer cette fraction en éléments simples pour

trouver deux racines distincte et pouvoir appliquer la transformé de Laplace
inverse :

A B A(1 + RCp) + Bp A + p(ARC + B)

vs (p) = + = = (4.8)
p 1 + RCp p(1 + RCp) p(1 + RCp)

Par identification on trouve A et B :

(
A=1
(4.9)
ARC + B = 0
donc :

15
 (
A=1
(4.10)
B = −RC
On a alors :
1 RC 1 τ
vs (p) = − = − 1 (4.11)
p 1 + RCp p τ ( τ + p)
d’où :
1 1
vs (p) = − 1 (4.12)
p ( τ + p)
En passant alors par la transformé de Laplace inverse, d’après le tableau on
a:
(
L−1 ( p1 ) = EΓ(t)
t (4.13)
L−1 ( 1 1+p) = e− τ
τ

on obtient alors l’équation de vs(t) dans le domaine temporel :

t
vs (t) = (EΓ(t) − EΓ(t)e− τ ) (4.14)
Finalement :
t
vs (t) = EΓ(t)(1 − e− τ ) (4.15)

4.3 Résolution détaillée d’une équation différentielle

du second ordre d’un circuit RLC avec Ve(t)
fonction échelon
Avant de Résoudre l’équation différentille nous l’avons établie grâce au schéma
électrique donné dans le TP (en annexe) réalisé au premier cycle préparatoire :

16
 d’après le schéma on a l’équation de ve tel que :

di(t)
Ve (t) = Vs (t) + Ri(t) + L (4.16)
dt
or Z
1
Vc = idt (4.17)
C
donc
dVc
i(t) = C (4.18)
dt
de plus la tension aux bornes du condensateur est ici la tension de sortie
donc on notera Vc =Vs , on aura donc :
dVs
i=c (4.19)
dt
on aura donc en remplaçant ces notations dans l’équation précédente :

dV s(t) d dVs (t)

Ve (t) = Vs (t) + RC + L (c ) (4.20)
dt dt dt
donc :

dV s(t) d2 V s(t)
Ve (t) = Vs (t) + RC + LC (4.21)
dt dt
Or dans l’ équation donnée dans le TP est sous la forme :

d2 V s(t) dVs (t)

+ 2mω0 + ω02 Vs (t) = ω02 Ve (t) (4.22)
dt2 dt
On divise alors les deux membres de l’équation (4.21) par LC :

Ve (t) Vs (t) R dV s(t) d2 V s(t)

= + + (4.23)
LC LC L dt dt
en posant :
(q
1
LC = ω0
(4.24)
R
2mω0 = L

donc :
q
1

LC = ω0
√ q (4.25)
m = R R LC R C
2Lω0 = 2L = 2 L

On obtient alors l’équation (4.22)

Maintenant à partir de l’équation (4.22) on doit appliquer la méthode de
Laplace pour obtenir sa solution. Avant tout on pose les transformés de Laplace
suivante :

17
 

 L(Ve (t)) = ve (p)

L(V (t)) = v (p)
s s
dVs (t) (4.26)
L( ) = pvs (p)
 d2dt



vs (t) 2
L( dt2 ) = p vs (p)
On passe maintenant l’équation (4.22) dans le modèle de Laplace :

p2 vs (p) + 2ω0 mpvs (p) + ω02 vs (p) = ω02 ve (p) (4.27)

On isole vs (p) :

ω02 ve (p)
vs (p) = (4.28)
p2 + 2mω0 p + ω02
or dans ce cas nous devons résoudre l’équation avec ve une fonction échelon
or d’après le tableau des transformés de Laplace on a :
1
L(Ve (t)) = p = ve (p)
D’où :

ω02
vs (p) = (4.29)
p(p2 + 2mω0 p + ω02 )
Il nous faut alors décomposer cette racine (décomposition en éléments simples)
pour trouver deux racines distincte et pouvoir appliquer la transformé de La-
place inverse :
A B C
vs (p) = ω02 ( + + ) (4.30)
p p − p1 p − p2
calculons les coefficient A,B,C :
ω02

A = limpvs (p) = ω2 = 1

 0
ω 2 (p−p )
1 ω02
B = lim (p − p1 )vs (p) = lim p(p−p
0
=
p1 (p1 −p2 )
1 )(p−p2 )
C = lim (p − p )v (p) = lim ω02 (p−p2 ) = ω02


2 s p(p−p1 )(p−p2 ) p2 (p2 −p1 )
A présent,nous devons calculer les racines p1 et p2 du polynôme correspon-
dant au dénominateur de la formule de vs(p) :
nous avons le polynôme du 2éme degrè suivant :

p2 + 2mω0 p + ω02 = 0 (4.31)

donc :

∆ = (2mω0 )2 − 4ω02 = 4m2 ω02 − 4ω02 = 4ω02 (m2 − 1) = −4ω02 (1 − m2 ) (4.32)

or d’après l’énnoncé du Tp en annexe m est comprit entre 0 et 1 donc :

√ p
∆ < 0 ⇒ ∆ = 2iω0 1 − m2 (4.33)
donc on obtient :
( √
p1 = ω0 (−m + i 1 − m2 ) = ε1 ω0
√ (4.34)
p2 = ω0 (−m − i 1 − m2 ) = ε2 ω0

18
 avec :
( √
ε1 = −m + i 1 − m2
√ (4.35)
ε2 = −m − i 1 − m2
finalement en remplaçant les valeurs de A,B et C dans vs et en continuant
à utiliser p1 et p2 pour ne pas charger les calculs on obtient :

1 1 1
vs (p) = ω02 ( + √ − √ ) (4.36)
p (2iε1 1 − m )(p − p1 ) (2iε2 1 − m2 )(p − p2 )
2

Appliquons désormais la transformation inverse de Laplace :

1 1
Vs (t) = ω02 (Γ(t) + √ e p1 t − √ e p2 t (4.37)
2iε1 1 − m2 2iε2 1 − m2
Nous pouvons alors poser que ε2 est le conjugué de ε1 :

ω2 e p1 t ep2 t
Vs (t) = ω02 Γ(t) + √ 0 ( − ∗ ) (4.38)
2i 1 − m2 ε1 ε1
∗
ω2 eε1 ω0 t eε1 ω0 t
Vs (t) = ω02 Γ(t) + √ 0 ( − ) (4.39)
2i 1 − m2 ε1 ε∗1

√ √
1−m2 1−m2
ω2 ei ω0 t e−i ω0 t
Vs (t) = ω02 Γ(t) + √ 0 e−mω0 t ( − ) (4.40)
2i 1 − m2 ε1 ε∗1

√ √
1−m2 ω0 t 1−m2 ω0 t
ω2 ε∗ ei − ε1 e−i
Vs (t) = ω02 Γ(t) + √ 0 e−mω0 t ( 1 2 ) (4.41)
2i 1 − m2 |ε1 |
Maintenant pour simplifier les calculs, posons f(t) tel que :
√
1−m2 ω0 t
f (t) = ε∗1 ei (4.42)
on aura alors :

√ p √ √ p √
1−m2 ω0 t 2 2 2
f (t) = −mei −i 1 − m2 ei 1−m ω0 t +me−i 1−m ω0 t −i 1 − m2 e−i 1−m ω0 t
(4.43)

√ √ p √ √
1−m2 ω0 t 1−m2 ω0 t 2 2
f (t) = −m(ei − e−i ) − i 1 − m2 (ei 1−m ω0 t − e−i 1−m ω0 t )
(4.44)
On remarque alors que,diviser par 2i, on obtient la formule d’Euler pour le
sin, on a donc à multiplier par 2i pour équilibrer. Ce qui nous donne :

p p p
f (t) = −m2isin( 1 − m2 ω0 t) − i 1 − m2 2isin( 1 − m2 ω0 t) (4.45)

p p p
f (t) = 2i(−m − i 1 − m2 )sin( 1 − m2 ω0 t) = 2iε∗1 sin( 1 − m2 ω0 t) (4.46)

19
 en reprenant l’expression de Vs dans le domaine temporel on obtient :

ω 2 e−mω0 t p
Vs (t) = ω02 Γ(t) + √0 2 2iε∗1 sin( 1 − m2 ω0 t) (4.47)
2i 1 − m2 |ε1 |

finalement on aura Vs(t) tel que :

√
Γ(t) + e−mω0 t ε∗1 sin( 1 − m2 ω0 t)
Vs (t) = ω02 ( 2√ ) (4.48)
|ε1 | 1 − m2

En vue de la contrainte de page imposé dans le cahier des charges du projet,

nous n’avons pas intégrer deux autres cas de résolution d’équation différentielle
avec la méthode de Laplace. En effet, vous trouverez en annexe la résolution
d’une équation différentielle d’un circuit RLC avec Ve(t) fonction d’impulsion
unitaire ainsi qu’avec ve(t) fonction sinusoı̈dale. Toutes ces fonctions étudiées
dans cette partie (et son complément en annexe) présentent tous les cas possible
de résolution pour un circuit RLC.Je vous invite donc vivement de les regarder.

20
Chapitre 5

Moédlisation et application

5.1 Résolution d’équations de circuits avec la

méthode de Laplace sous mapple et affichage
des graphiques associés
5.1.1 Réponse d’un circuit RC à un échelon de tension
Le premier problème que nous avons réussi à traiter correctement avec Maple
fut la simulation de la tension aux bornes du condensateur dans un circuit RC.
L’équation différentielle à résoudre étant du 1er ordre, nous avons pu utiliser
les fonctions de base de Maple pour obtenir un résultat conforme à ce que nous
avions prévu.
Pour obtenir un graph intéressant nous avons pris les valeurs d’un TP de
P3 (R = 300Ω; C = 10−9 F ) Puis on a crée un signal d’entrée en utilisant
la fonction Heaviside. Pour un échelon de tension de 0 à 5 Volts on obtient
l’équation suivante :
E(t) = 5.Heaviside(t) (5.1)
Dans l’optique de la résolution de notre problème par Maple on entre les condi-
tions initiales ainsi que l’équation différentielles qui régis ce circuit.

∂(V s(t)) V S(t) E

+ = (5.2)
∂t τ τ
Puis via l’utilisation de la fonction Maple dsolve en lui spécifiant d’utiliser la
méthode de Laplace on obtient l’équation de la tension au borne du condensateur
pour t > 0. On obtient ainsi le graph ci-dessous.

21
5.1.2 Réponse d’un circuit RLC à un échelon de tension
On a ensuite cherché à simuler grâce à Maple la réponse d’un circuit RLC
à un échelon de tension. Cette fois ci nous avons simulé la décharge du conden-
sateur, la tension d’entrée sera donc défini par :

E(t) = E(1 − Heaviside(t)) (5.3)

Nous avons cette fois du procéder par étape pour résoudre l’équation de
second ordre suivante :
∂ 2 (U c(t)) R ∂(U c(t)) U c(t)
+ × + =0 (5.4)
∂t2 L ∂t LC
Nous avons aussi posé et résolus l’équation caractéristique pour obtenir les
résultats intermédiaires et confirmer les calculs de la machine.

22
5.1.3 Réponse d’un circuit RC à un signal créneau
Nous nous sommes aussi attaqués au problème de la résolution d’un système
avec une entrée créneau. Le premier problème fut de modéliser le signal d’entrée,
n’existant pas, à notre connaissance, de fonction mathématiques donnant un
signal créneau, nous avons recherché de nombreuses solutions, l’idée d’utiliser
des séries de Fourrier nous est évidement venu à l’esprit mais le manque de
connaissance que nous avions sur le sujet nous a assez rapidement écarté de
cette piste, nous nous sommes donc rabattu sur l’utilisation d’une somme de
piecewise ce qui nous permettait d’obtenir un faible nombre de créneau, mais
de façon peu couteuse.

L’équation différentielles bien qu’étant la même que pour l’échelon de tension

apparait sur plusieurs page, ceci est dû à la définition en piecewise de la tension
d’entrée. Malgré le caractère illisible de l’équation différentielle ainsi que de la
solution trouvé par Maple grâce à la méthode de Laplace, on a tout de même pu
obtenir un graph bien plus compréhensible qui correspond au résultat attendu
puisque que l’on observe la charge, puis les cycles charge/décharge d’un circuit
RLC soumis à une tension créneau.

23
5.1.4 Conclusion sur notre travail Maple
N’ayant jamais auparavant travaillé sur Maple, nous avons dû consacrer une
grande partie de notre temps à la recherche, essayer de trouver des travaux
de résolution d’équation différentielles sur ce logiciel pouvant nous indiquer une
voix, nous en avons trouvé quelques-uns, que ce soit sur internet ou sur les nom-
breux documents que Mr Gleyse a mis à notre disposition. Notre problème ma-
jeur a été de tout d’abords comprendre les commandes de base Maple, d’éviter
les erreurs syntaxique, et de lire la doc des commandes les plus récurrentes,
mais une fois les bases du logiciel domptées, nous avons pu progresser dans nos
simulations. La simulation des entrées en échelons a été plutôt facile à modéliser
une fois les équations différentielles recalculées. Comme dit précédemment, la
simulation d’un circuit en entrée échelons, nous a donné bien plus de mal. Les
majeurs problèmes, auront été de tout d’abords trouvé une modélisation du
signal d’entrée compatible avec la fonction Maple  dsolve , nous avons es-
sayé une série de commande présente dans un des documents fourni par notre
professeur, qui produisait une composée de fonction Heaviside, mais la fonction
résultante n’était (à notre grand étonnement) pas considéré comme définie aux
points de discontinuité par maple. Nous avons aussi travaillé sur la simulation
d’un circuit RLC à signal d’entrée en créneau, mais aucun des programmes
résultants ne nous a donné des résultats suffisamment satisfaisants pour être
évoqués dans ce rapport.

24
5.2 Conception d’un récepteur radio AM grâce
aux circuits RLC
Pour bien comprendre les applications des circuits RLC, nous avons suivit 6
séances avec M. BENJELLOUN professeur d’électronique a l’ESIGELEC.
Suite aux cours suivit nous avons pu voir que les circuits RLC sont souvent
utilisés dans la création et le traitement de signaux (ondes électromagnétiques).
Il existe différentes applications de circuit RLC dans les radio AM la première
est le filtre et la seconde l’oscillateur. On sépare le fonctionnement des radios
en 2 parties : l’émetteur et le récepteur.

5.2.1 L’émetteur :
L’émetteur est composé de 4 parties principales :
- Le signal à transmettre de basse fréquence [F1 ; F2 ] (20Hz à 25kHz fréquences
des sons audible par l’homme).
- Le multiplicateur ou modulateur : ce composant multiplie le signal de mo-
dulation par la porteuse (un signal pure de haute fréquence f0 ) pour avoir un
signal de haute fréquence qui contient l’information à transmettre de fréquence
centrale f0 , c’est le signal modulé.
- Un passe bande qui ne laisse passer qu’une bande (-25kHz a 25kHz) centré sur
la fréquence centrale.
- L’antenne qui convertit le signal modulé en ondes életromagnétiques.

Le fonctionnement :
Le signal à transmettre (un son) est ajouté a une composante continue (pour
éviter la sur-modulation), on obtient alors le signal de modulation. Ensuite on
multiplie ce signal avec la porteuse pour obtenir un signal modulé, c’est la
modulation. Le modulateur fait varier l’amplitude de la porteuse en fonction du
signal de modulation. La fréquence de la porteuse est fixe pour chaque radio et
elle est générée par un circuit RLC en parallèle. Le signale modulé est filtré par
un passe bande qui élimine les harmoniques de la porteuse (2f0 , 3f0 , 4f0 , etc.),
comme tout filtre c’est un circuit RLC.
In1 : Signal de modulation, il est de la forme :

Ue (t) = A + Bcos(Ωt) avec Ω = 2πF (5.5)

F ∈ [F1 ; F2 ] tel que F1 ≥ 20Hz et F2 ≤ 25kHz (5.6)

In2 : signal de la porteuse, elle est de la forme :
1
U0 (t) = Ccos(ω0 t) avec ω0 = 2πf0 et f0 = √ (5.7)
2π LC
Out1 : Le signal module, il est de la forme :

Out1(t) = Ue (t) × U0 (t) (5.8)

Out2 : Le signal module apres filtrage des harmonique

25
 Figure 5.1 – Schéma synoptique d’un emetteur

Figure 5.2 – Signal de modulation

5.2.2 Le récepteur :
Le récepteur est composé de 5 parties :
-L’antenne de reception qui capte toutes les ondes émisent dans l’atmosphère et
les convertit sous la forme d’un signal electrique.
-Un préselecteur (circuit LC) qui élimine les fréquences images ainsi que les
parasites.
-Un multiplicateur ou mélangeur, ce mutiplicateur a pour objetif de mutltiplier
tout les signaux reçus par une certaine frequence fp réglable. Un des signaux
emis de fréquence centrale f0 sera centré sur une frequence intermediaire, FI.
-Un filtre intermediaire, ce filtre passe bande ne laisse passer qu’une bande (-
25kHz à 25kHz) centré sur la frequence intermediaire fixe pour tout les radio
(455kHz pour les radio AM).
-Un demodulateur, il a la fonction inverse du modulateur. Il permet de lire
l’information contenu dans le signal modulé et on obtient une copie du signal
audio original.

Le fonctionnement :
L’antenne reçoit et convertit les ondes de l’atmosphere en signal électrique
par l’intermedaire d’un circuit RLC. Une premiere étape est de purifier ce signal
en éléminant les fréquences images (ou interferences) et les parasites. Pour per-
mettre au recepteur d’avoir qu’un seul filtre, le mélangeur permet de choisir la
porteuse qu’on désire démoduler. Pour éliminer tous les autres signaux reçu, le
signal crée par le mélangeur passe par le filtre intermediaire qui élimine tout les
signaux non centré sur FI, donc il nous reste un seul signal puisque FI est fixe
selon la fréquence fp choisit on a pas le meme signal de fréquences centrale f0 .

26
 Figure 5.3 – Signal de la porteuse

Figure 5.4 – Signal modulé

Ce signal est ensuite démodulé pour enfin obtenir le signal originale. Ce dernier
peut-être utilisé pour le brancher à un haut parleur.

Figure 5.5 – Schéma synoptique d’un recepteur

In3 : Signal module, ce signal est composé de tout les signaux modulés emis
(Out1), pour extraire l’information qu’on désire il faut séparer la porteuse du
reste, il est de la forme :

In3(t) = Out1(t) + Out10 + Out100 + Out1(3) + Out1(4) + ... (5.9)

Le signal de fréquence fp est de la meme forme que la porteuse de In2.

Out3 : Ce signal est de la forme :

Out3(t) = In3(t) × Up (5.10)

Out4 : C’est le signal modulé qu’on a choisit, suite au filtre, il a pour
fréquence centrale FI selon l’equation suivant :

F I = |fp − f0 | (5.11)
Out5 : Le signal original (on l’obtient avec la demodulation de Out4).

Remarque : Pour bien comprendre le systeme de modulation, les multiplicateurs

sont basés sur une diode qui permet de faire la multiplication selon l’equation :

27
 1
cosΩt × cosω0 t = [cos (ω0 t + Ωt) + cos (ω0 t − Ωt)] (5.12)
2
En multipliant les signaux on obtient un signal particulier. En regardant la
transformé de Fourier on remarque qu’on a une fréquence centrale f0 ainsi que
les fréquences f0 +F et f0 − F . C’est là que l’information est stockée dans la
porteuse.

28
Chapitre 6

Conclusion

Ce projet nous aura permis à tous de mieux comprendre non pas tant le
fonctionnnement des circuits RLC, puisque nous les avions déjà étudié en STPI,
que leur modes de résolution et leur utilité. En effet, la méthode de résolution
de Laplace et les simulations Maple nous ont permis d(avoir une vision plus
globale de ces circuits, puisque la simplification de leur résolution nous a donné
du recul en nous permettant d’en étudié un plus grand nombre que l’année
dernière. C’est justement cette facilité à les étudier et à les simuler qui nous a
donné envie d’aller plus loin, et d’appliquer nos connaissances. C’est la raison
pour laquelle nous avons pris des cours supplémentaires, de façon à pouvoir
exploiter au mieux notre projet, et simuler un émetteur/récepteur radio. C’est
cette partie non-obligatoire qui nous a permis à tous de comprendre l’intéret
et l’utilisation de logiciels comme Maple ou d’outils mathématiques comme les
transformées de Laplace.

29
Troisième partie

Annexes

30
Chapitre 7

Méthode de Laplace

7.1 Résolution détaillée d’une équation différentielle

du second ordre d’un circuit RLC avec Ve(t)
fonction d’impulsion unitaire
On a l’équation différentielle suivante à résoudre :

d2 V s(t) dVs (t)

+ 2mω0 + ω02 Vs (t) = ω02 Ve (t) (7.1)
dt2 dt
or on sait que :


 L(Ve (t)) = ve (p)

L(V (t)) = v (p)
s s
dVs (t) (7.2)
L( ) = pvs (p)
 d2dt



vs (t) 2
L( dt2 ) = p vs (p)
on obtient donc grâce aux transformés de Lapalce :

ω02
vs (p) = (7.3)
p2 + 2mω0 p + ω02
d’après le tableau des transformés inverses de Laplace en Annexe on sait
que :

sin(ωp t)e−mωp t
 
−1 1
L 2 2
= (7.4)
p + 2ζωn p + ωn ωp
donc dans notre cas on aura :

sin(ω0 t)e−mω0 t
 
1
L−1 2 2 = (7.5)
p + 2mω0 p + ω0 ω0
Donc Vs (t) est sous la forme :

sin(ω0 te−mω0 t
Vs (t) = (7.6)
ω0

31
7.2 Résolution détaillée d’une équation différentielle
du second ordre d’un circuit RLC avec Ve(t)
fonction sinusoı̈dale
Nous démarons toujours avec la même équation différentielle qui est celle
donnée dans le TP :

d2 V s(t) dVs (t)

2
+ 2mω0 + ω02 Vs (t) = ω02 Ve (t) (7.7)
dt dt
or on sait ,d’près le tableau des transformés de Laplace en Annexe, que :


 L(Ve (t)) = ve (p)

L(V (t)) = v (p)
s s
dVs (t) (7.8)
L( ) = pvs (p)
 d2dt



vs (t) 2
L( dt2 ) = p vs (p)
en les remplaçants dans l’équation différentielle précédente :

p2 vs (p) + 2ω0 mpvs (p) + ω02 vs (p) = ω02 ve (p) (7.9)

en isolant vs :

ω02 ve (p)
vs (p) = (7.10)
p2 + 2mω0 p + ω02
or Ve(t) est une fonction sinusoı̈dale tel que :
ω
Ve (t) = U0 sin(ωt) ⇒ ve (p) = U0 (7.11)
p2 + ω 2
en remplaçant dans vs(p) :

ω02 ωU0
vs (p) = (7.12)
(p2 + 2
ω )(2mω0 p + p2 + ω02 )
décomposons maintenant cette expression en éléments simples :
Ap + B Cp + D
vs (p) = 2 + 2 (7.13)
p2 + 2mω0 p + ω0 p + ω2
en développant :

Ap3 + Aω 2 p + Bp2 + Bω 2 + Cp3 + Dp2 + 2mω0 Cp2 + 2mω0 Dp + ω02 Cp + ω02 D

vs (p) =
(p2 + ω 2 )(2mω0 p + p2 + ω02 )
(7.14)
d’où :

p3 (A + C) + p2 (B + D + 2mω0 C) + p(Aω 2 + 2mω0 D + ω02 C) + (Bω 2 + ω02 D)

vs (p) =
(p2 + ω 2 )(2mω0 p + p2 + ω02 )
(7.15)
soit, par identification :

32
 

 A+C =0

B + D + 2mω C = 0
0
2
(7.16)


 Aω + 2mω 0 D + ω02 C = 0
 2
Bω + ω02 D = 1


 C = −A

B + D − 2mω A = 0
0
(7.17)


 Aω + 2mω0 D − ω02 A = 0
2
 2
Bω + ω02 D = 1


C = −A
B = A(2mω − omega20 −ω2 )


0 2mω0
ω02 −ω 2 (7.18)


D = A 2mω0

 2 ω 2 −ω 2
Bω + ω0 02m A = 1
On a donc deux expressions de A en fonction de B :
ω02 −ω 2
(
A(−2mω0 + 2mω )+B =0
2 ω0
0 (7.19)
2
A(ω0 − ω ) 2m + Bω 2 = 1
donc :  
0 1 

1 ω2 
A =   (7.20)
ω02 −ω 2 
 −2ω0 2mω0 

 2 ω0
(ω0 − ω 2 ) 2m 2 
ω
donc :
−1
A= (ω02 −ω 2 )2
(7.21)
−2mω0 ω 2 − 4m2
on obtient alors les valeurs de B,C et D grâce à A :
4m2

A = 8m3 ω2 ω0 −ω 4 4 2 2
0 −ω +2ω0 ω



B = − 2m
 −4m2 −ω 2
ω0 8m3 ω 2 ω0 −ω04 −ω 4 +2ω02 ω 2
(7.22)


 C = −A
ω02 −ω 2

D = A 2mω

0

d’où l’expression de vs(p) :

ω02 −ω 2 ω02 −ω 2
A(p + 2mω0 − 2mω0 )ω0 ω (−Ap + A 2mω0 )
vs (p) = + ωω0 (7.23)
p2 + 2mω0 p + ω02 p2 + ω2

en factorisant par Aω0 ω :

33
  
ω02 −ω 2 ω02 −ω 2
p 2mω0 p 2mω0 
vs (p) = Aω0 ω  − 2 − 2 +
p2 + 2mω0 p + ω02 p + 2mω0 p + ω02 p + ω2 p2 + ω 2
(7.24)
Appliquons alors les transformés inverse de Laplace (d’après le tableau) :

p
√
√ −1 e−mω0 t sin(ω0 1 − m2 t − θ)
 −
 L 1( p2 +2mω 0 p+ω 2) =
1−m 2
 0

 ω02 −ω 2
ω02 −ω 2

L 1( p2 +2mω0 p+ω2 ) = 2mω2 sin(ω0 t)e−mω0 t
− 2mω0





 0 0



p
 L− 1( p2 +ω 2 ) = cos(ωt)
(7.25)





 ω02 −ω 2
ω02 −ω 2

 − 2mω0


 L 1( p2 +ω2 ) = 2mω 0ω
sin(ωt)



avec le theta de l’équation (1) du système qui vaut :

θ = cos−1 (m) (7.26)

Il ne reste plus qu’à remplacer ces transformés inverses pour retrouver l’ex-
pression de Vs(t) (domaine temporel) :

Vs (t) =
 √ ω02 −ω 2 ω02 −ω 2

ω0 ωA √ −1 e−mω0 t sin(ω0 1− m2 t − θ) − sin(ω0 t)e−mω0 t − cos(ωt) +
1−m2 2mω02 2mω0 ω sin(ωt)

34
Chapitre 8

Méthode de Laplace
appliquée sur maple (code
mapple)

8.1 fonction échelon d’un circuit RC

35
>

>

(1)

>
(2)
>

>

(3)
 (3)
>

(4)
>

(5)
>

(6)
>

>
8.2 fonction crénau d’un circuit RC

38
>

(1)

>

>

(2)

>

(3)

>
>

(4)

>
>
 (5)
>

(6)
(5)

(6)
>
>
 (7)

>
>
>
8.3 fonction échelon d’un circuit RLC

49
>
>

(1)
>
(2)
>
>

>

(3)
>

(4)
>
(5)
>
(6)
>

(7)
>

(8)

>