Dans la fenêtre d’impression, choisissez « Enregistrer au format PDF » pour garder une copie sur votre appareil.
1Formules et méthodes
Convolution et systèmes LTI
Convolution discrète
y[n] = Σ x[k]·h[n − k]
Convolution continue
y(t) = ∫ x(τ)·h(t − τ) dτ
Longueur et somme
Ly = Lx + Lh − 1, Σy = Σx·Σh
Début, fin et aire
t1 + t3 … t2 + t4, ∫y = ∫x·∫h
Échelon dans e^(−at)u(t)
y(t) = (A/a)(1 − e−at)
Système du premier ordre
y′ + a·y = b·x → h(t) = b·e−atu(t)
Récurrence du premier ordre
y[n] = a·y[n − 1] + b·x[n] → h[n] = b·aⁿ
Sinusoïde traversant h(t)
|H(jω)|, ∠H(jω); H(jω) = b/(a + jω)
Amplitudes de la série de Fourier
carré ±A : 4A/(kπ), k impair ; triangle ±A : 8A/(π²k²), k impair ; train d’impulsions : A·D et (2A/(kπ))|sin kπD|
Méthode
Esquissez les deux signaux et notez où chacun est non nul.
Retournez h et faites-la glisser : pour chaque instant de sortie, trouvez la plage où les deux signaux se recouvrent.
Multipliez et additionnez (cas discret) ou multipliez et intégrez (cas continu) sur ce seul recouvrement.
Vérifiez le résultat : la sortie commence à la somme des instants de début, se termine à la somme des instants de fin, et son aire ou sa somme est le produit des aires ou des sommes des entrées.
2Exercices
1.Convolution et systèmes LTIMoyen
L’entrée x(t) = 10·e−5tu(t) est appliquée à un système de réponse h(t) = e−2tu(t). Trouvez la sortie y(1,5) à t = 1,5 s.
x(t) = 10·e−5tu(t)h(t) = e−2tu(t)t = 1,5 s
y(1.5) =
2.Convolution et systèmes LTIDifficile
L’entrée x(t) = 10·cos(20t) est appliquée à un système de réponse h(t) = 2·e−10tu(t). Trouvez la phase φ de la sortie en régime permanent par rapport à l’entrée, en degrés.
h(t) = 2e−10tu(t)x(t) = 10cos(20t)
φ = °
3.Transformées de LaplaceMoyen
Trouvez f(t) pour F(s) = (2s + 30)/(s² + 8s + 25) en complétant le carré, et calculez f(0,4).
F(s) = (2s + 30)/(s² + 8s + 25)t = 0,4 s
f(0.4) =
4.Transformées de LaplaceDifficile
Un système obéit à y″ + 6y′ + 5y = 15·u(t), avec y(0) = 4 et y′(0) = 0. À t = 0, on applique un échelon d’entrée u(t) = 10. Utilisez la transformée de Laplace pour trouver y(0,75).
Un système a pour fonction de transfert G(s) = 24/(s² + 3s + 12) et reçoit un échelon unitaire. Trouvez le taux d’amortissement ζ.
G(s) = 24/(s² + 3s + 12)
ζ =
6.Réponse indicielleDifficile
Une masse m = 20 kg est suspendue à un ressort k = 5000 N/m avec un amortisseur visqueux c = 126 N·s/m (m·x″ + c·x′ + k·x = F). Au repos à l’équilibre, une force constante F = 5 N est soudain appliquée. Trouvez le taux d’amortissement ζ.
m = 20 kgk = 5000 N/mc = 126 N·s/mF = 5 N
ζ =
7.Diagrammes de Bode et margesMoyen
Trouvez la phase ∠G(jω) de G(s) = 26(s + 4)/(s(s + 9)) à ω = 5 rad/s, en degrés.
G(s) = 26(s + 4)/(s(s + 9))ω = 5 rad/s
∠G(jω) = °
8.Diagrammes de Bode et margesDifficile
Une boucle à retour unitaire a pour fonction de transfert en boucle ouverte L(s) = 6/(s(s + 4)). Trouvez la marge de phase PM en degrés.
L(s) = 6/(s(s + 4))
PM = °
3Solutions
1.Convolution et systèmes LTI
Réponse: y(1.5) = 0,1641
Convolution continue : y(t) = ∫ x(τ)·h(t − τ) dτ. On retourne h, on la fait glisser jusqu’à la position t et on intègre le produit des parties qui se recouvrent. Pour des signaux causaux, l’intégrale va de 0 à t. Résultats utiles : l’aire de y est l’aire de x multipliée par l’aire de h ; si x dure de t₁ à t₂ et h de t₃ à t₄, y dure de t₁ + t₃ à t₂ + t₄. Pour h(t) = e−atu(t) et un échelon de hauteur A, y(t) = (A/a)(1 − e−at).
Réponse impulsionnelle et série de Fourier : la réponse impulsionnelle h(t) décrit entièrement un système LTI : pour y′ + a·y = b·x, c’est h(t) = b·e−atu(t), et pour y[n] = a·y[n − 1] + b·x[n], c’est h[n] = b·aⁿu[n]. Une sinusoïde cos ωt traverse le système avec le gain |H(jω)| et la phase ∠H(jω), où H est la transformée de Fourier de h. Un signal périodique est une somme d’harmoniques en k/T : un signal carré de ±A n’a que des harmoniques impairs d’amplitudes 4A/(kπ), un signal triangulaire de crête A a 8A/(π²k²), et un train d’impulsions de hauteur A et de rapport cyclique D a pour moyenne A·D et pour harmoniques (2A/(kπ))·|sin(kπD)|.
Réponse en fréquence
H(jω) = 2/(10 + jω)
Phase à la fréquence d’entrée
φ = −atan(20/10) = −63,43°
Réponse
φ = −63,43°
3.Transformées de Laplace
Réponse: f(0.4) = 1,526
Transformée de Laplace : F(s) = ∫₀∞ f(t)·e−st dt transforme les dérivées en multiplications par s : ℒ{f′} = sF(s) − f(0), ℒ{f″} = s²F(s) − s·f(0) − f′(0). Couples usuels : 1 ↔ 1/s, e−at ↔ 1/(s + a), t·e−at ↔ 1/(s + a)², sin ωt ↔ ω/(s² + ω²), cos ωt ↔ s/(s² + ω²), et un décalage s → s + σ multiplie la fonction du temps par e−σt. Pour inverser une F(s) rationnelle, décomposez-la en éléments simples : le résidu d’un pôle simple p est A = [(s − p)·F(s)] en s = p. Théorème de la valeur finale : f(∞) = lims→0 s·F(s), valable seulement si tous les pôles de s·F(s) sont dans le demi-plan gauche. Théorème de la valeur initiale : f(0⁺) = lims→∞ s·F(s).
Transformée de Laplace : F(s) = ∫₀∞ f(t)·e−st dt transforme les dérivées en multiplications par s : ℒ{f′} = sF(s) − f(0), ℒ{f″} = s²F(s) − s·f(0) − f′(0). Couples usuels : 1 ↔ 1/s, e−at ↔ 1/(s + a), t·e−at ↔ 1/(s + a)², sin ωt ↔ ω/(s² + ω²), cos ωt ↔ s/(s² + ω²), et un décalage s → s + σ multiplie la fonction du temps par e−σt. Pour inverser une F(s) rationnelle, décomposez-la en éléments simples : le résidu d’un pôle simple p est A = [(s − p)·F(s)] en s = p. Théorème de la valeur finale : f(∞) = lims→0 s·F(s), valable seulement si tous les pôles de s·F(s) sont dans le demi-plan gauche. Théorème de la valeur initiale : f(0⁺) = lims→∞ s·F(s).
Transformer l’équation
(s² + 6s + 5)·Y(s) = 4s + 24 + 150/s
Y(s) = (4s² + 24s + 150)/(s(s + 1)(s + 5))
Décomposition en éléments simples
A = 150/(1·5) = 30
B = 130/((−1)·4) = −32,5
C = 130/((−5)·(−4)) = 6,5
Résultat
y(t) = 30 − 32,5·e−t + 6,5·e−5t
y(0,75) = 14,8
Réponse
y(0,75) = 14,8
Vérification : y(0) = A + B + C = 4.
5.Réponse indicielle
Réponse: ζ = 0,433
Réponse indicielle : un système du premier ordre G(s) = K/(τs + 1) répond à un échelon d’amplitude U par y(t) = K·U·(1 − e−t/τ) : il atteint 63,2 % de la valeur finale K·U après une constante de temps τ et 98 % après 4τ. Un système du second ordre G(s) = K·ωₙ²/(s² + 2ζωₙs + ωₙ²) a pour pôles s = −ζωₙ ± jωd avec ωd = ωₙ√(1 − ζ²). Pour 0 < ζ < 1, la réponse indicielle dépasse sa valeur finale : dépassement en pourcentage %OS = 100·e−ζπ/√(1−ζ²), temps du premier pic tp = π/ωd, et temps de réponse à 2 % ts ≈ 4/(ζωₙ). Comparez le dénominateur avec s² + 2ζωₙs + ωₙ² pour lire ωₙ et ζ. Pour un système masse–ressort–amortisseur, ωₙ = √(k/m) et ζ = c/(2√(km)) ; pour un circuit RLC série, ωₙ = 1/√(LC) et ζ = (R/2)·√(C/L).
Comparez avec la forme standard
Le dénominateur doit correspondre à s² + 2ζωₙs + ωₙ².
ωₙ² = 12 → ωₙ = 3,464 rad/s
2ζωₙ = 3 → ζ = 3/(2·3,464) = 0,433
Résultat
ζ = 0,433
Réponse
ζ = 0,433
6.Réponse indicielle
Réponse: ζ = 0,1992
Réponse indicielle : un système du premier ordre G(s) = K/(τs + 1) répond à un échelon d’amplitude U par y(t) = K·U·(1 − e−t/τ) : il atteint 63,2 % de la valeur finale K·U après une constante de temps τ et 98 % après 4τ. Un système du second ordre G(s) = K·ωₙ²/(s² + 2ζωₙs + ωₙ²) a pour pôles s = −ζωₙ ± jωd avec ωd = ωₙ√(1 − ζ²). Pour 0 < ζ < 1, la réponse indicielle dépasse sa valeur finale : dépassement en pourcentage %OS = 100·e−ζπ/√(1−ζ²), temps du premier pic tp = π/ωd, et temps de réponse à 2 % ts ≈ 4/(ζωₙ). Comparez le dénominateur avec s² + 2ζωₙs + ωₙ² pour lire ωₙ et ζ. Pour un système masse–ressort–amortisseur, ωₙ = √(k/m) et ζ = c/(2√(km)) ; pour un circuit RLC série, ωₙ = 1/√(LC) et ζ = (R/2)·√(C/L).
Pulsation propre et coefficient d’amortissement
ωₙ = √(k/m) = √(5000/20) = 15,81 rad/s
ζ = c/(2√(km)) = 126/(2·√(5000·20)) = 0,1992
Résultat
ζ = 0,1992
Réponse
ζ = 0,1992
7.Diagrammes de Bode et marges
Réponse: ∠G(jω) = −67,71 °
Réponse en fréquence : substituez s = jω dans la fonction de transfert. Le module en décibels est 20·log₁₀|H(jω)| et la phase est la somme des angles des facteurs du numérateur moins ceux des facteurs du dénominateur. Un filtre RC ou RL du premier ordre a sa fréquence de coupure (−3 dB) à ωc = 1/(RC) ou R/L, et fc = ωc/(2π) ; un filtre passe-bas donne |H| = 1/√(1 + (ω/ωc)²) et φ = −atan(ω/ωc). Dans le diagramme de Bode asymptotique (en lignes droites), chaque pôle infléchit le module vers le bas de 20 dB/décade à sa fréquence de coupure, chaque zéro l’infléchit vers le haut de 20 dB/décade, et un intégrateur 1/s donne −20 dB/décade à toutes les fréquences. Pulsation de croisement du gain ωc : |L(jω)| = 1 (0 dB) ; marge de phase PM = 180° + ∠L(jωc). Pulsation de croisement de phase ω180 : ∠L(jω) = −180° ; marge de gain GM = −20·log₁₀|L(jω180)| dB. La boucle fermée est stable quand les deux marges sont positives.
Réponse en fréquence : substituez s = jω dans la fonction de transfert. Le module en décibels est 20·log₁₀|H(jω)| et la phase est la somme des angles des facteurs du numérateur moins ceux des facteurs du dénominateur. Un filtre RC ou RL du premier ordre a sa fréquence de coupure (−3 dB) à ωc = 1/(RC) ou R/L, et fc = ωc/(2π) ; un filtre passe-bas donne |H| = 1/√(1 + (ω/ωc)²) et φ = −atan(ω/ωc). Dans le diagramme de Bode asymptotique (en lignes droites), chaque pôle infléchit le module vers le bas de 20 dB/décade à sa fréquence de coupure, chaque zéro l’infléchit vers le haut de 20 dB/décade, et un intégrateur 1/s donne −20 dB/décade à toutes les fréquences. Pulsation de croisement du gain ωc : |L(jω)| = 1 (0 dB) ; marge de phase PM = 180° + ∠L(jωc). Pulsation de croisement de phase ω180 : ∠L(jω) = −180° ; marge de gain GM = −20·log₁₀|L(jω180)| dB. La boucle fermée est stable quand les deux marges sont positives.