Skip to content

Systèmes d'Ordre 1 et 2

Systèmes de Premier Ordre

Modélisation

Équation différentielle

L'équation différentielle d'un système de premier ordre peut s'exprimer sous la forme suivante :

τdy(t)dt+y(t)=b1dx(t)dt+b0x(t)
  • τ désigne la constante de temps [s].

Fonction de transfert

La fonction de transfert d'un système de premier ordre peut s'exprimer sous la forme normalisée suivante :

H(s)=N(s)τs+1
  • N(s) désigne le numérateur de la fonction de transfert (polynôme de degré inférieur ou égal à 1).

Expression du pôle

Un système de premier ordre possède un unique pôle. Ce pôle s'obtient en cherchant l'unique racine du dénominateur de la fonction de transfert.

τs+1=0p=1τ

Exemples de circuits

Filtre RC Passe-Bas

Circuit : Résistance R en série avec condensateur C. Sortie aux bornes de C.

HLP(s)=T0τs+1

avec τ=RC et T0=1.

Filtre RC Passe-Haut

Circuit : Condensateur C en série avec résistance R. Sortie aux bornes de R.

HHP(s)=Tτsτs+1

avec τ=RC et T=1.

Réponse Temporelle

La solution complète de l'équation différentielle s'exprime sous la forme :

y(t)=yl(t)+yp(t)
  • yl(t): solution libre (régime libre)
yl(t)=λe1τt
  • yp(t): solution particulière (régime forcé). L'expression du régime forcé dépend de l'allure de l'entrée et des coefficients b0 et b1

Exemple : Réponse indicielle d'un passe-bas

Considérons la réponse d'un système de premier ordre à un échelon d'amplitude E, c-à-d x(t)=Eu(t). Le système est supposé initialement au repos (y(0)=0). Comme l'entrée est un échelon, le régime forcé est de la forme yp(t)=αu(t). En remplaçant cette expression dans l'équation différentielle pour t0, nous obtenons yp(t)=T0E. Il en vient que :

y(t)=λe1τt+T0E

La constante d'intégration s'obtient en déterminant une condition initiale. En intégrant l'équation différentielle entre t=0 et t=0+, nous obtenons y(0+)=0. En exploitant cette équation, il en vient que :

y(0+)=λ×1+T0E=0λ=T0E,

Pour un filtre passe-bas de premier ordre, la réponse du système à un échelon d'amplitude E s'exprime finalement sous la forme :

y(t)=T0E(1etτ)
Réponse Indicielle d'un filtre passe-bas de premier ordre
Valeurs Remarquables
  • Valeur initiale : y(0+)=0,
  • Valeur finale : y()=T0E,
  • Valeur maximale: max(y(t))=T0E,
  • Temps de réponse à ±5%: tr=3τ.

Réponse Fréquentielle

La réponse fréquentielle s'obtient en posant s=jω dans l'expression de la fonction de transfert. Nous obtenons :

H(jω)=N(jω)jωωc+1
  • ωc=1τ désigne la pulsation de coupure à -3dB [rad/s].

Passe-Bas

Réponse Fréquentielle d'un filtre passe-bas de premier ordre
  • Module :
|H(jω)|=T0(ωωc)2+1
  • Argument :
arg[H(jω)]=arctan(ωωc)

Passe-Haut

Réponse Fréquentielle d'un filtre passe-haut de premier ordre
  • Module :
|H(jω)|=Tωωc(ωωc)2+1
  • Argument :
arg[H(jω)]=90oarctan(ωωc)

Systèmes de Second Ordre

Modélisation du système

Équation différentielle

L'équation différentielle d'un système de second ordre peut s'exprimer sous la forme suivante :

1ω02d2y(t)dt2+2mω0dy(t)dt+y(t)=b2d2x(t)dt2+b1dx(t)dt+b0x(t)
  • ω0 désigne la pulsation propre [rad/s],
  • m désigne le coefficient d'amortissement,
  • Q=12m désigne le facteur de qualité.

Fonction de transfert

La fonction de transfert d'un système de second ordre peut s'exprimer sous la forme normalisée suivante :

H(s)=N(s)1ω02s2+2mω0s+1
  • N(s) désigne le numérateur de la fonction de transfert (polynôme de degré inférieur ou égal à 2).

Exemples de circuits

Les 4 exemples ci-dessous présentent les fonctions de transfert normalisées (a0=1) de second ordre. À l'exception du passe-haut (m=n), ces fonctions de transfert sont des fonctions rationnelles propres (m<n).

Passe-bas (LP)

HLP(s)=T01ω02s2+2mω0s+1

Exemple de circuit : Filtre RLC série ou structure Sallen-Key passe-bas.

Passe-bande (BP)

HBP(s)=2mTmω0s1ω02s2+2mω0s+1

Exemple de circuit : Circuit RLC série avec sortie aux bornes de R.

Passe-haut (HP)

HHP(s)=Tω02s21ω02s2+2mω0s+1

Exemple de circuit : Structure Sallen-Key passe-haut.

Réjecteur (BR)

HBR(s)=T0(1ω02s2+1)1ω02s2+2mω0s+1

Exemple de circuit : Circuit bouchon RLC.

Analyse du système

Expressions des Pôles

Les pôles s'obtiennent en déterminant les racines du dénominateur de la fonction de transfert c-à-d en déterminant les valeurs de s telles que :

1ω02s2+2mω0s+1=0

Cette équation est une équation du second degré. Le discriminant s'exprime sous la forme suivante:

Δ=(2mω0)24ω02=4ω02(m21)

L'expression du discriminant montre que la valeur de m joue un rôle essentiel dans l'expression des deux racines. Spécifiquement, nous pouvons distinguer 3 cas de figure.

Cas où m>1

Lorsque m>1, le système possède deux pôles réels et distincts. Les pôles s'expriment sous la forme

p1=mω0+ω0m21p2=mω0ω0m21

Mathématiquement, les deux pôles possèdent les propriétés suivantes:

  • le produit des deux pôles est égal à p1p2=m2ω02ω02(m21)=ω02
  • la somme des deux pôles est égale à p1+p2=2mω0.

En utilisant ces 2 propriétés, il est possible d'identifier le coefficient d'amortissement et la pulsation propre via les relations :

ω0=p1p2m=p1+p22ω0

Géométriquement, les deux pôles sont placés sur un cercle de centre mω0 et de rayon ω0m21.

Cas où m=1

Lorsque m=1, le système possède un pôle réel double. Le pôle réel double s'exprime sous la forme

p1=p2=mω0
Cas où m<1

Lorsque m<1, le système possède une paire de pôles complexes-conjugués. Les pôles s'expriment sous la forme

p1=mω0+jω01m2p2=mω0jω01m2

Mathématiquement, il est possible de démontrer que :

  • le module de chaque pôle est égal à |p1|=|p2|=(mω0)2+ω02(1m2)=ω0.
  • l'angle formé entre le pôle p1 et l'axe des réels est donné par cos(θ)=Re(p1)/|p1|=m

En utilisant ces 2 propriétés, il est possible d'identifier le coefficient d'amortissement et la pulsation propre via les relations :

ω0=p1p1=|p1|m=Re(p1)ω0

Équation Différentielle et Réponse Temporelle

En utilisant l'équation différentielle, la réponse à une entrée quelconque s'exprime sous la forme :

y(t)=yl(t)+yp(t)
  • yl(t): solution libre,
  • yp(t): solution particulière.

Cas où m>1

La solution libre s'exprime sous la forme :

yl(t)=λ1e(mm21)ω0t+λ2e(m+m21)ω0t

Nous constatons que la solution libre est donnée par la contribution de deux systèmes de premier ordre ayant pour constantes de temps respectives:

τ1,2=1p1,2=1ω0(m±m21)

Cas où m<1

La solution libre s'exprime sous la forme d'une sinusoïde amortie:

yl(t)=αemω0tcos(ω01m2t+φ)
  • mω0 régit la vitesse de décroissance de l'enveloppe,
  • ω01m2 correspond à la pseudo-pulsation des oscillations [rad/s]
Démonstration

Lorsque m<1, comme la solution libre est réelle, cette solution s'exprime nécessairement sous la forme

yl(t)=λ1ep1t+λ2ep2t=λ1ep1t+λ1ep1t

λ1 est un complexe. Comme cette solution est composée d'un complexe et de son complexe conjugué, il en vient que :

yl(t)=2Re(λ1ep1t)

Cette expression peut également se simplifier en utilisant la décomposition sous forme polaire de la constante d'intégration et la décomposition sous forme algébrique de p1. En effet, nous trouvons :

yl(t)=2Re(|λ1|ejarg[λ1]e(Re(p1)+jIm(p1))t)=2|λ1|eRe(p1)tRe(ej(Im(p1)t+arg[λ1]))=2|λ1|eRe(p1)tcos(Im(p1)t+arg[λ1])

En utilisant l'expression de p1 et en posant α=2|λ1| et φ=arg[λ1], nous obtenons finalement :

yl(t)=αemω0tcos(ω01m2t+φ)
Propriétés
  • Oscillations : pseudo-pulsation et pseudo-période.
ωp=ω01m2Tp=2πω01m2
  • Ratio des amplitudes après une oscillation :
R=yl(t1)yl(t1+Tp)=e2πm1m2>1

Lorsque m1, ωpω0 et Re2πm. Cette propriété peut être utilisée pour identifier rapidement m à partir de deux dépassements consécutifs de même signe.

Cas de la Réponse Indicielle

La réponse indicielle correspond à la réponse du système lorsque l'entrée est un échelon d'amplitude E, c-à-d :

x(t)={Esi t00ailleurs

L'expression de la sortie peut s'obtenir en utilisant l'équation différentielle ou la transformée de Laplace inverse. Spécifiquement, il est possible d'établir que la transformée de Laplace de la sortie est égale à Y(s)=H(s)Es. Dans le domaine temporel, la sortie s'obtient en utilisant une transformée de Laplace inverse. L'expression exacte de y(t) nécessite alors à la fois une décomposition en éléments simples et l'utilisation des tables des transformées de Laplace.

Valeurs limites en 0+ et

Concernant les valeurs aux limites, il est possible d'obtenir rapidement le comportement en t (valeur finale) et à la discontinuité c-à-d t=0+ (valeur initiale) en utilisant les théorèmes de la valeur finale et initiale.

  • Valeur finale: La valeur finale est donnée par
y()=lims0sY(s)=H(0)E
  • Valeur initiale. Pour la valeur initiale, nous nous intéressons à la fois à la valeur initiale de y(t) et sa dérivée dy(t)dt.

    • Concernant y(t), nous obtenons :

      y(0+)=limssY(s)=H()E
    • Concernant la dérivée de y(t), nous obtenons

      y˙(0+)=dy(t)dt|0+=limss2Y(s)=limssH(s)E
Cas du passe-haut

Dans le cas du passe-haut, la valeur de la limite de sH(s)E ne converge pas car la fonction de transfert est une fraction rationnelle impropre. Dans ce cas, il est nécessaire de décomposer la fonction de transfert sous la forme H(s)=H+H1(s)H est un gain et H1(s) est une fonction rationnelle propre. En utilisant cette décomposition et le fait que H n'impacte pas la valeur de y˙(0+) (terme associé à une impulsion de Dirac), il en vient que

y˙(0+)=limssH1(s)E

Cette limite converge vers une valeur finie.

Réponse Fréquentielle

La réponse fréquentielle du système s'obtient en posant s=jω. En utilisant la forme générale d'un second ordre, nous obtenons :

H(jω)=N(jω)1ω2ω02+2jmωω0
  • Expression à la pulsation propre. À la pulsation propre ω=ω0 [rad/s], la réponse fréquentielle s'exprime simplement sous la forme :
H(jω0)=N(jω0)2jm

Expression du module

Le module s'exprime sous la forme générale suivante :

|H(jω)|=|N(jω)|(1Ω2)2+4m2Ω2=|N(jω)|((1Ω2)2+4m2Ω2)1/2

Ω=ωω0 correspond à la pulsation normalisée. En fonction de l'expression du numérateur et de la valeur de m, ce module peut présenter un extremum à la pulsation ωr. Si cet extremum existe, la pulsation ωr s'obtient en annulant la dérivée du module. Lorsque l'extremum correspond à un maximum, la pulsation ωr est appelée pulsation de résonance.

Expression de l'argument

Pour éviter les changements de cadrants liés à l'utilisation de la fonction arctangente, une pré-étape consiste à décomposer la fonction de transfert sous une forme plus simple à manipuler. Spécifiquement, une approche possible consiste à normaliser la partie réelle au dénominateur à 1 pour éviter les indéterminations de phase de la fonction arctangente. En utilisant cette normalisation, nous obtenons :

arg[H(jω)]=arg[N(jω)]arctan(12m(ΩΩ1))

Ω=ωω0 correspond à la pulsation normalisée.

Démonstration

En décomposant la fonction de transfert sous la forme suivante :

H(jω)=N(jω)×11+j2m(ωω0ω0ω)×12jmωω0

Le deuxième terme possède un dénominateur avec une partie réelle positive (pas de changement de cadrant). Nous obtenons alors

arg[H(jω)]=arg[N(jω)]arg(12m(ωω0ω0ω))π2

Formes Normalisées

Passe-bas (LP)

H(s)=T01ω02s2+2mω0s+1
  • gain statique: T0,
  • pulsation propre: ω0 (rad/s),
  • coefficient d'amortissement: m.

Réponse Indicielle

Expression
  • Lorsque m>1, la réponse indicielle s'exprime sous la forme :
y(t)=T0E(1m+m212m21et/τ1mm212m21et/τ2)u(t)

τ1=1/p1 et τ2=1/p2 avec τ1<τ2.

  • Lorsque m<1, la réponse indicielle s'exprime sous la forme :
y(t)=T0E(1emω0t1m2cos(ω01m2tarctan(m1m2)))u(t)
Valeurs Limites
  • Valeur initiale: y(0+)=0 et dy(t)dt|t=0+=0.
  • Valeur finale: y()=T0E
Dépassement (si m<1)

Lorsque m<1, il est possible d'établir que la réponse indicielle présente des extrema en tk=kπω01m2.

Démonstration

Pour obtenir la position des extrema lorsque m<1, il faut dériver l'expression de la réponse indicielle. Pour t0, nous trouvons :

dy(t)dt=ω01m2emω0tsin(ω01m2t)

Cette dérivée s'annule lorsque ω01m2t=kπ avec kR+. Il en vient que :

dy(t)dt|t=tk=0tk=kπω01m2

Il est important de pouvoir quantifier la valeur maximale de la réponse indicielle. En pratique, nous utilisons le plus souvent le premier dépassement relatif qui est une grandeur indépendante de E, T0 et de ω0. Cette grandeur s'exprime sous la forme

Dr(%)=100|max(y(t))y()y()|=100eπm1m2

La valeur du premier dépassement peut s'obtenir rapidement en utilisant une abaque présentant Dr(%)=f(m) (voir abaque).

Temps de réponse

Dans certaines applications, il est nécessaire de quantifier la rapidité du système au moyen du temps de réponse à ±5%. Le temps de réponse tr est défini formellement comme le temps nécessaire pour atteindre définitivement ±5% de la valeur finale y(). Mathématiquement, tr peut être défini de la manière suivante :

tr=inf{t0τt,|y(τ)y()|0.05|y()|}

Il est possible de démontrer que ce temps de réponse dépend uniquement de la pulsation propre ω0 et du coefficient d'amortissement m. Comme cette dépendance est complexe, la valeur du temps de réponse s'obtient généralement en utilisant une abaque présentant ω0tr=g(m) (voir abaque).

Comportement Fréquentiel

Point singulier et comportements asymptotiques
  • Asymptote basse-fréquences: HBF(jω)=T0,
  • Asymptote haute-fréquences: HHF(jω)=T0(ωω0)2,
  • Valeur à la pulsation propre: H(jω0)=T02jm
Phénomène de résonance

Si m<0.7, la dérivée du module s'annule à la pulsation de résonance :

ωr=ω012m2<ω0

À la pulsation de résonance, le module est égal à

|H(jωr)|=T02m1m2>|H(jω0)|

Passe-bande (BP)

H(s)=2mTmω0s1ω02s2+2mω0s+1
  • gain maximum: Tm,
  • pulsation propre: ω0 (rad/s),
  • coefficient d'amortissement: m.

Réponse Indicielle

Valeurs Limites
  • Valeur initiale: y(0+)=0 et dy(t)dt|t=0+=2mω0TmE.
  • Valeur finale: y()=0

Comportement Fréquentiel

Point singulier et comportements asymptotiques
  • Asymptote basse-fréquences: HBF(jω)=2jmTm(ωω0),
  • Asymptote haute-fréquences: HHF(jω)=2jmTm(ωω0)1,
  • Intersection des asymptotes de module en ω=ω0: |Ti|=2m|Tm|
  • Valeur à la pulsation propre: H(jω0)=Tm
Pulsation de coupure à -3dB et bande passante
  • Pulsation de coupure. Le filtre passe-bande possède deux pulsations de coupure à 3dB (|H(jωc)|=|Tm|2). Ces pulsations s'expriment sous la forme :
ωc1=ω0(m+m2+1)ωc2=ω0(m+m2+1)
  • Bande passante à 3 dB. La bande passante à 3dB s'exprime sous la forme :
Δω=ωc2ωc1=2mω0

Passe-haut (HP)

H(s)=Tω02s21ω02s2+2mω0s+1
  • gain haute-fréquence: T,
  • pulsation propre: ω0 (rad/s),
  • coefficient d'amortissement: m.

Réponse Indicielle

Valeurs Limites
  • Valeur initiale: y(0+)=TE et dy(t)dt|t=0+=2mω0TE.
  • Valeur finale: y()=0

Comportement Fréquentiel

Point singulier et comportements asymptotiques
  • Asymptote basse-fréquences: HBF(jω)=T(ωω0)2,
  • Asymptote haute-fréquences: HHF(jω)=T,
  • Valeur à la pulsation propre: H(jω0)=jT2m
Phénomène de résonance

Si m<0.7, la dérivée du module s'annule à la pulsation de résonance :

ωr=ω012m2>ω0

À la pulsation de résonance, le module est égal à

|H(jωr)|=T2m1m2>|H(jω0)|

Réjecteur (Notch)

H(s)=T0(1ω02s2+1)1ω02s2+2mω0s+1
  • gain maximum: T0,
  • pulsation propre: ω0 (rad/s),
  • coefficient d'amortissement: m.

Réponse Indicielle

Valeurs Limites
  • Valeur initiale: y(0+)=T0E et dy(t)dt|t=0+=2mω0T0E.
  • Valeur finale: y()=T0E

Comportement Fréquentiel

Point singulier et comportements asymptotiques
  • Asymptote basse-fréquences: HBF(jω)=T0,
  • Asymptote haute-fréquences: HHF(jω)=T0,
  • Valeur à la pulsation propre: H(jω0)=0
Pulsation de coupure à -3dB et bande rejetée
  • Pulsation de coupure. Le filtre réjecteur possède deux pulsations de coupure à 3dB (|H(jωc)|=|T0|2). Ces pulsations s'expriment sous la forme :
ωc1=ω0(m+m2+1)ωc2=ω0(m+m2+1)
  • Bande rejetée à 3 dB. La bande rejetée à 3dB s'exprime sous la forme :
Δω=ωc2ωc1=2mω0