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
Le signal x(t) n’est non nul qu’entre t = −2 s et t = 1 s, et h(t) qu’entre t = 2 s et t = 3 s. Trouvez l’instant tend où y(t) = x(t) ∗ h(t) se termine.
x: −2 … 1 sh: 2 … 3 s
tend = s
2.Convolution et systèmes LTIDifficile
Un signal carré commute entre 1 et 0 avec la période T = 4 ms (une demi-période à chaque niveau). Trouvez l’amplitude c3 de sa composante de Fourier à 750 Hz.
0 … 1T = 4 msf = 750 Hz
c3 =
3.Transformées de LaplaceMoyen
Trouvez f(t) pour F(s) = (4s + 36)/(s² + 2s + 1) et calculez f(1).
F(s) = (4s + 36)/(s² + 2s + 1)t = 1 s
f(1) =
4.Transformées de LaplaceDifficile
Un système obéit à y″ + 6y′ + 8y = 24·u(t), avec y(0) = 6 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,1).
Un système a pour fonction de transfert G(s) = 35,5/(s² + 5s + 71) et reçoit un échelon unitaire. Trouvez le temps de réponse à 2 % ts de la réponse indicielle, avec ts ≈ 4/(ζωₙ).
G(s) = 35,5/(s² + 5s + 71)
ts = s
6.Réponse indicielleDifficile
Une masse m = 5 kg est suspendue à un ressort k = 400 N/m avec un amortisseur visqueux c = 27 N·s/m (m·x″ + c·x′ + k·x = F). Au repos à l’équilibre, une force constante F = 100 N est soudain appliquée. Trouvez le dépassement %OS en pourcentage du déplacement.
m = 5 kgk = 400 N/mc = 27 N·s/mF = 100 N
%OS = %
7.Diagrammes de Bode et margesMoyen
Trouvez la phase ∠G(jω) de G(s) = 37(s + 6)/(s(s + 3)) à ω = 1 rad/s, en degrés.
G(s) = 37(s + 6)/(s(s + 3))ω = 1 rad/s
∠G(jω) = °
8.Diagrammes de Bode et margesDifficile
Une boucle à retour unitaire a pour fonction de transfert en boucle ouverte L(s) = 34/(s(s + 1)(s + 6)). Trouvez la marge de gain GM en dB.
L(s) = 34/(s(s + 1)(s + 6))
GM = dB
3Solutions
1.Convolution et systèmes LTI
Réponse: tend = 4 s
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).
Début et fin
tstart = −2 + 2 = 0 s
tend = 1 + 3 = 4 s
Réponse
tend = 4 s
Réponse
tend = 4 s
2.Convolution et systèmes LTI
Réponse: c3 = 0,2122
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)|.
Rang de l’harmonique
f₀ = 1/T = 250 Hz → k = 750/250 = 3
Amplitude
c3 = 4·0,5/(3π) = 0,2122
Réponse
c3 = 0,2122
3.Transformées de Laplace
Réponse: f(1) = 13,24
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).
Factoriser le dénominateur
s² + 2s + 1 = (s + 1)²
Décomposition en éléments simples
4s + 36 = 4(s + 1) + 32
F(s) = 4/(s + 1) + 32/(s + 1)²
Transformée inverse
f(t) = 4·e−t + 32·t·e−t
f(1) = (4 + 32·1)·e−1 = 13,24
Réponse
f(1) = 13,24
4.Transformées de Laplace
Réponse: y(0.1) = 6,789
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 + 8)·Y(s) = 6s + 36 + 240/s
Y(s) = (6s² + 36s + 240)/(s(s + 2)(s + 4))
Décomposition en éléments simples
A = 240/(2·4) = 30
B = 192/((−2)·2) = −48
C = 192/((−4)·(−2)) = 24
Résultat
y(t) = 30 − 48·e−2t + 24·e−4t
y(0,1) = 6,789
Réponse
y(0,1) = 6,789
Vérification : y(0) = A + B + C = 6.
5.Réponse indicielle
Réponse: ts = 1,6 s
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 + ωₙ².
ωₙ² = 71 → ωₙ = 8,426 rad/s
2ζωₙ = 5 → ζ = 5/(2·8,426) = 0,2967
Résultat
ζωₙ = 2,5 → ts ≈ 4/2,5 = 1,6 s
Réponse
ts = 1,6 s
6.Réponse indicielle
Réponse: %OS = 36,98 %
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) = √(400/5) = 8,944 rad/s
ζ = c/(2√(km)) = 27/(2·√(400·5)) = 0,3019
Résultat
%OS = 100·e−ζπ/√(1−ζ²) = 36,98 %
Réponse
%OS = 36,98 %
7.Diagrammes de Bode et marges
Réponse: ∠G(jω) = −98,97 °
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.