Leçon
La théorie — Laboratoire de transformée de Laplace
La transformée de Laplace échange une fonction du temps contre une fonction d’une variable complexe s, changeant l’analyse en t en algèbre en s. Le catalogue ci-dessus est un dictionnaire de huit paires, et à côté de chacune la page affiche quatre faits qui se lisent sur F(s) seul, sans jamais revenir au domaine temporel.
Ce que signifie chaque symbole
F(s)- la transformée. Un quotient de polynômes pour chaque entrée de ce catalogue — c’est la seule raison d’être des trois lignes suivantes.
poles- les valeurs de
squi annulent le dénominateur. Tout le reste de la page découle de leur position. ROC- la région de convergence — le demi-plan où l’intégrale de définition converge réellement. L’étape 3 montre le peu d’information propre que porte cette ligne.
f(∞)- la valeur finale, affichée comme un nombre, comme
∞, ou comme un refus.
D’où vient la formule
- Lisez la paire. Pour l’exponentielle par défaut,
e^(−t)devient1/(s+1): un taux de décroissance dans le temps devient la position d’un pôle unique. - Trouvez les pôles en annulant le dénominateur.
s+1=0donne s = −1, ce qu’affiche la ligne des pôles. - Comparez maintenant la ligne des pôles et celle du ROC, et continuez à comparer en parcourant le catalogue. Le ROC est la partie réelle du pôle le plus à droite. Exponentielle : pôle en −1, ROC
Re(s) > −1. Échelon : pôle en 0, ROCRe(s) > 0. Sinus : pôles en±j2, dont la partie réelle est nulle, donc ROCRe(s) > 0à nouveau. Huit lignes, une règle. - Les deux lignes de valeurs viennent des deux extrémités du même
F(s):f(0⁺) = lim(s→∞) sF(s)etf(∞) = lim(s→0) sF(s). D’où l’accord de l’échelon et du cosinus surf(0⁺) = 1et leur désaccord sur la fin.
Comment lire ce que vous voyez
Parcourez les huit entrées et regardez deux colonnes apparemment indépendantes se comporter comme une seule. Échelon, rampe et puissance ont leurs pôles à l’origine et affichent tous Re(s) > 0. L’exponentielle et les deux sinusoïdes amorties ont un pôle de partie réelle −1 et affichent toutes Re(s) > −1. Sinus et cosinus semblent faire exception — leurs pôles sont imaginaires — et n’en font pas : la partie réelle de ±j2 est nulle, donc ils affichent Re(s) > 0 comme ceux de l’origine. Regardez ensuite les fins, que les pôles décident aussi : une partie réelle −1 donne 0, un pôle simple à l’origine donne une constante, un pôle double à l’origine donne ∞, et une paire purement imaginaire ne donne aucune fin.
- Suppose
- Que
f(t) = 0avantt = 0— c’est la transformée unilatérale, ce qui explique que chaque entrée soit écrite contre un échelon unité. Que l’amplitudeAmetF(s)à l’échelle sans déplacer de pôle, donc modifie la valeur initiale et jamais le ROC ni la fin. Et que chaque entrée est un quotient de polynômes : une transformée avec coupure ou singularité essentielle n’a pas de liste de pôles à lire, et trois des quatre lignes n’auraient rien à dire. - Ne tient plus quand
- Le théorème de la valeur finale est celui qu’on cite sans ses conditions, et ce catalogue contient le contre-exemple. Appliqué aveuglément au sinus,
lim(s→0) s·ω/(s²+ω²)vaut0: net, définitif, et faux —sin(ωt)n’a pas de valeur finale. Le théorème ne vaut que si les pôles desF(s)sont strictement dans le demi-plan gauche, et ceux du sinus sont sur l’axe imaginaire. La page refuse au lieu de rapporter la limite : c’est pourquoi cette ligne indique n’existe pas — oscille et non le 0 que la formule produit volontiers. La rampe échoue au même test par l’autre côté et reçoit∞. Deux des huit entrées brisent la règle qu’on allait généraliser à partir des six autres.
Problèmes entièrement résolus
-
Valeurs aux extrémités et décroissance de f(t) = e −1,5t 6 étapes
f(t) = e−1,5t pour t ≥ 0. Obtenez F(s) à partir de la définition, lisez le pôle et la région de convergence au fil du même calcul, et déterminez les deux valeurs limites sans jamais inverser de transformée. Répondez ensuite à la question à laquelle aucun des deux théorèmes des valeurs limites ne répond : combien de temps la décroissance prend-elle réellement ?
-
La transformée est une intégrale unique, et tout le reste sur le panneau est un sous-produit de son évaluation. Les deux exponentielles partageant la même base, elles se regroupent en un seul exposant, −(s + 1,5)t.
-
C'est à la borne supérieure que réside la difficulté. e−(s+1,5)t tend vers 0 lorsque t augmente uniquement quand la partie réelle de s + 1,5 est positive, et cette unique inégalité vous fournit les deux éléments restants : le domaine de convergence Re(s) > −1,5 et le pôle en s = −1,5 situé exactement sur sa frontière. Le pôle n'est pas une propriété isolée à mémoriser — c'est la valeur de s pour laquelle l'intégrale cesse d'exister.
-
Démontrez le théorème de la valeur initiale plutôt que de simplement le citer. Une valeur élevée de s pondère e−st si fortement vers les petites valeurs de t que sF(s) devient une sonde de f au voisinage de l'origine. Ici, s/(s + 1,5) tend vers 1, et e0 = 1 le confirme — vous avez lu f(0+) à partir d'une transformée que vous n'avez jamais inversée.
-
Faites tendre s dans l'autre sens pour la valeur finale. Cette limite est assortie d'une condition d'application : elle ne constitue un théorème que si chaque pôle de sF(s) se trouve strictly à gauche de l'axe imaginaire, et l'unique pôle ici se situe en −1,5. Le cas oscillant ci-dessous est précisément ce dont cette condition vous protège.
-
Aucun de ces deux théorèmes ne comporte d'horloge. L'un décrit t = 0 et l'autre t = ∞, et entre les deux, la courbe peut mettre une microseconde ou un siècle. Résolvez f(τ) = f(0)/e pour trouver l'échelle manquante.
-
Laissez maintenant le panneau de côté et demandez-vous quand le signal est réduit à 1 % de sa valeur initiale.
Réponse
F(s) = 1/(s + 1,5), un pôle en −1,5 et τ = 0,667 s — le panneau affiche le pôle et la constante de temps, et vous avez dérivé les deux à partir d'une seule intégrale. Conservez plutôt la dernière ligne. Puisque a s'élimine de ln(100)/a divisé par 1/a, atteindre 1 % prend toujours 4,61 constantes de temps, quel que soit le taux de décroissance : 4,61 est une propriété du nombre 100, non de ce signal. Tout ce que contrôle le pôle est la durée d'un τ. Doublez a pour le porter à 3 et chaque durée sur cette page est divisée par deux. Faites-le glisser dans l'autre sens, vers l'axe imaginaire, et τ tend vers l'infini — et une fois que le pôle atteint l'axe, la valeur finale cesse complètement d'exister, ce qui fait l'objet du problème suivant.
-
-
La limite de la valeur finale pour f(t) = sin(4t) qui donne 0 5 étapes
f(t) = sin(4t) pour t ≥ 0. La limite de la valeur finale s'effectue sans aucun problème ici et donne 0. Rien n'explose et rien n'est indéterminé. Expliquez pourquoi 0 est néanmoins la mauvaise réponse, et ce que le théorème vous demandait réellement de vérifier.
-
Construisez cette transformée à partir du résultat du problème précédent plutôt qu'à l'aide d'une table. Euler décompose un sinus en une paire d'exponentielles complexes, et chacune obéit à la règle que vous avez déjà établie à partir de l'intégrale.
-
Réduisez au même dénominateur. Le 2j se simplifie avec j4 − (−j4) au numérateur, de sorte qu'une construction complexe donne une fonction réelle de s. Le panneau écrit le dénominateur sous la forme s² + 4², ce qui revient au même.
-
Localisez les pôles. Le dénominateur s'annule en s = ±j4, et leur partie réelle est 0 — ni négative, ni positive, mais exactement sur l'axe imaginaire. C'est précisément le cas que la condition ci-dessus visait à exclure.
-
Calculez quand même la limite de la valeur finale, car rien ne vous en empêche. Le numérateur tend vers 0, le dénominateur tend vers 16, et le quotient vaut 0. L'arithmétique ne vous avertira pas que vous n'en aviez pas le droit.
-
Comparez cela avec la fonction elle-même. sin(4t) a pour période 2π/4 = 1,5708 s, ce que le panneau arrondit à 1,57 s, et au quart de chacune de ces périodes, elle revient à +1. Une fonction qui continue de revenir indéfiniment à 1 ne peut pas converger vers 0, ni vers quoi que ce soit. Le panneau n'affiche aucun nombre dans la ligne de la valeur finale et indique plutôt que le signal oscille.
Réponse
Pôles en ±j4, période 1,57 s, et aucune valeur finale — bien que la limite s'évalue proprement à 0. La leçon n'est pas que l'arithmétique est fragile ; elle est que l'arithmétique n'est pas le théorème. Comparez les deux sF(s) côte à côte et la différence réside dans la position d'un seul pôle. Ce qui rend 0 intéressant plutôt que simplement faux, c'est qu'il s'agit de la réponse correcte à chaque question voisine : amortissez le sinus de n'importe quel a > 0 et les pôles se déplacent en −a ± j4, le théorème s'applique, et la valeur finale est réellement 0 — pour tout a, aussi petit soit-il. La limite donne la réponse à tous les problèmes voisins et à aucun de celui-ci. Et un 0 prouvable peut tout de même être inutile : pour a = 0,1, l'horloge du problème 1 donne ln(100)/0,1 = 46 s pour descendre à 1 %, soit environ 29 cycles complets d'oscillation à cette fréquence.
-
Références (3)
- When the final-value theorem applies, and when the table must say it does not: E. Gluskin, "Let us teach this generalization of the final-value theorem." European Journal of Physics 24(6), 591–597, 2003.
- Pole location and stability, at length: Katsuhiko Ogata, Modern Control Engineering, 5th edition. Prentice Hall, 2010. ISBN 978-0-13-615673-4.
- Where the region of convergence is derived from the pole positions: A. V. Oppenheim & A. S. Willsky, Signals and Systems, 2nd edition. Prentice Hall, 1997 — chapter 9, on the Laplace transform and its region of convergence. ISBN 978-0-13-814757-0.