Dans la fenêtre d’impression, choisissez « Enregistrer au format PDF » pour garder une copie sur votre appareil.
1Formules et méthodes
Équations différentielles
Croissance/décroissance exponentielle
dy/dt = ky → y = y0ekt
Demi-vie
k = ln 2 / t½
Loi de refroidissement de Newton
T(t) = Ta + (T0 − Ta)e−kt
Facteur intégrant
y′ + p(x)y = q(x): μ = e∫p dx, (μy)′ = μq
Cuve de mélange
dm/dt = r·c − (r/V)·m
Matrices et déterminants
déterminant 2×2
det = ad − bc
Inverse d’une matrice 2×2
A⁻¹ = (1/det)·[d −b; −c a]
Produit
(AB)ij = Σk AikBkj
Valeurs propres
det(A − λI) = 0
Règle de Cramer
xi = det(Ai)/det(A)
Vecteurs : produit scalaire et produit vectoriel
Norme
|a| = √(ax² + ay² + az²)
Vecteur unitaire
û = a/|a|
Composantes à partir d’un angle
Fx = F·cos θ, Fy = F·sin θ
Produit scalaire
a·b = axbx + ayby + azbz = |a||b|cos θ
Projection scalaire de a sur b
a·b/|b|
Travail
W = F·d
Produit vectoriel
a × b = (aybz − azby, azbx − axbz, axby − aybx)
Aires
parallélogramme |a × b|, triangle ½|a × b|
Moment par rapport à O
MO = r × F, r = A − O
Volume d’un parallélépipède
V = |a·(b × c)|
Méthode
Écrivez chaque vecteur en composantes ; pour une force donnée par une norme et un angle, utilisez F·cos θ et F·sin θ.
Pour les longueurs et les directions, utilisez la norme √(x² + y² + z²) et le vecteur unitaire a/|a|.
Utilisez le produit scalaire pour les angles, les projections, la perpendicularité et le travail : a·b = |a||b|cos θ.
Utilisez le produit vectoriel pour les vecteurs perpendiculaires, les aires et les moments : écrivez soigneusement les trois composantes, en gardant l’ordre a × b.
Intégrales
Règle des puissances
∫xⁿ dx = xⁿ⁺¹/(n + 1) + C
Intégrale définie
∫ab f dx = F(b) − F(a)
Primitives usuelles
∫eˣ dx = eˣ, ∫1/x dx = ln|x|, ∫sin x dx = −cos x
Aire entre deux courbes
∫ab (f − g) dx
Méthode
Réécrivez l’intégrande comme une somme de termes simples si possible.
Appliquez la règle des puissances : ∫xⁿ dx = xⁿ⁺¹/(n + 1), et les primitives usuelles de ex, 1/x, sin et cos.
Utilisez un changement de variable quand une fonction et sa dérivée apparaissent ensemble.
Pour une intégrale définie, calculez F(b) − F(a).
Analyse à plusieurs variables
Dérivées partielles
fx = ∂f/∂x (y constant), fxy = fyx
Gradient
∇f = (fx, fy), largest slope |∇f|
Dérivée directionnelle
Duf = ∇f·u, u = v/|v|
Plan tangent
z = f(a, b) + fx(x − a) + fy(y − b)
Dérivation en chaîne
dz/dt = fx·dx/dt + fy·dy/dt
Test des dérivées secondes
D = fxxfyy − fxy²: D > 0, fxx > 0 minimum; D > 0, fxx < 0 maximum; D < 0 saddle
Intégrale double
∬R f dA = ∫ab∫g(x)h(x) f dy dx
Coordonnées polaires
x = r cos θ, y = r sin θ, dA = r dr dθ
Méthode
Pour une dérivée partielle, traitez l’autre variable comme une constante et utilisez les règles à une variable.
Pour les questions de gradient, calculez ∇f = (fx, fy) au point ; utilisez |∇f| pour la pente maximale et ∇f·u avec un vecteur unitaire u pour une dérivée directionnelle.
Pour les extremums, résolvez fx = 0 et fy = 0 ensemble, puis classez chaque point avec D = fxx·fyy − fxy².
Pour une intégrale double, esquissez la région, choisissez l’ordre (ou les coordonnées polaires pour les disques) de sorte que les bornes intérieures soient les courbes qui la délimitent, et intégrez d’abord l’intérieur.
Suites et séries
Suite arithmétique
aₙ = a₁ + (n − 1)d, Sₙ = n(a₁ + aₙ)/2
Suite géométrique
aₙ = a₁rⁿ⁻¹, Sₙ = a₁(1 − rⁿ)/(1 − r)
Série géométrique infinie
S = a₁/(1 − r), |r| < 1
Série de Riemann
Σ1/kᵖ converges only for p > 1
Test du rapport
L = lim |ak+1/ak|: L < 1 converges, L > 1 diverges
Rayon de convergence
R = lim |ck/ck+1|
Polynôme de Taylor
Pₙ(x) = Σ f⁽ᵏ⁾(a)(x − a)ᵏ/k!
Série de Maclaurin
eˣ = Σxᵏ/k!, sin x = x − x³/3! + …, cos x = 1 − x²/2! + …, ln(1 + x) = x − x²/2 + …
Majoration de l’erreur de Lagrange
|Rₙ| ≤ M·|x − a|ⁿ⁺¹/(n + 1)!
Méthode
Identifiez le type : arithmétique (différence constante), géométrique (rapport constant), une série usuelle, une série entière ou un polynôme de Taylor.
Pour les sommes, utilisez Sₙ = n(a₁ + aₙ)/2, Sₙ = a₁(1 − rⁿ)/(1 − r), ou S = a₁/(1 − r) quand |r| < 1, avec a₁ le premier terme de la somme.
Pour la convergence : d’abord le test du terme général, puis comparez avec Σ1/kᵖ ou une série géométrique, et utilisez le test du rapport ou de la racine pour les factorielles et les puissances k-ièmes.
Pour les polynômes de Taylor, listez f(a), f′(a), f″(a), …, divisez la dérivée k-ième par k!, et majore l’erreur par M·|x − a|ⁿ⁺¹/(n + 1)!.
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.Équations différentiellesMoyen
Une pièce d’acier forgé à 300 °C refroidit dans un atelier à 20 °C. Elle suit la loi de refroidissement de Newton, dT/dt = −k·(T − Tₐ). Avec k = 0,02 min⁻¹, trouvez sa température T au bout de 10 min.
T0 = 300 °CTa = 20 °Ck = 0,02 min⁻¹t = 10 min
T = °C
2.Équations différentiellesDifficile
Résolvez y″ + y′ + 25,25y = 0 avec y(0) = 3 et y′(0) = 0 (lâché sans vitesse initiale). Trouvez y(0,75) ; angles en radians.
y″ + y′ + 25,25y = 0y(0) = 3y′(0) = 0
y(0.75) =
3.Matrices et déterminantsMoyen
Soit C = AB pour les matrices représentées. Trouve le coefficient c12 (ligne 1, colonne 2).
A = [5 −2; 2 −1]B = [−5 1; −3 5]
c12 =
4.Matrices et déterminantsDifficile
Trouve la plus grande valeur propre λ de la matrice A représentée.
A = [2 −4; −3 3]
λ =
5.Vecteurs : produit scalaire et produit vectorielMoyen
Trouvez l’angle θ entre les vecteurs a = (3, −5, 2) et b = (0, 3, −3), en degrés.
a = (3, −5, 2)b = (0, 3, −3)
θ = °
6.Vecteurs : produit scalaire et produit vectorielDifficile
Pour a = (−4, 1, 4) et b = (2, 6, 2), trouvez la composante z (a × b)z du produit vectoriel a × b.
a = (−4, 1, 4)b = (2, 6, 2)
(a × b)z =
7.IntégralesMoyen
Calcule I = ∫02 2x·(x² + 1)1 dx.
f(x) = 2x·(x² + 1)1[0, 2]
I =
8.IntégralesDifficile
Calculez I = ∫₀1 x²·eˣ dx.
∫₀1 x²eˣ dx
I =
9.Analyse à plusieurs variablesMoyen
Pour f(x, y) = ln(x² + 2y²), trouvez la dérivée directionnelle Duf en P(−1, 0,5) dans la direction de Q(−4, 0,5).
f(x, y) = ln(x² + 2y²)P = (−1, 0,5)Q = (−4, 0,5)
Duf =
10.Analyse à plusieurs variablesDifficile
Utilisez les coordonnées polaires pour calculer I = ∬R (x² + y²) dA, où R est la moitié supérieure du disque x² + y² ≤ 2,25 (y ≥ 0).
f = x² + y²0 ≤ r ≤ 1,50 ≤ θ ≤ π
I =
11.Suites et sériesMoyen
Appliquez le test du rapport à Σk=1∞ k³·5ᵏ/3ᵏ. Trouvez L = lim |ak+1/ak|.
ak = k³·5ᵏ/3ᵏ
L =
12.Suites et sériesDifficile
Utilisez le polynôme de Maclaurin de degré 4 de f(x) = e−x pour approcher f(1,5). Donnez P₄(1,5).
f(x) = e−xn = 4x = 1,5
P₄(1.5) =
13.Transformées de LaplaceMoyen
Trouvez f(t) pour F(s) = 45/(s(s + 2)(s + 6)) par décomposition en éléments simples, et calculez f(0,75).
F(s) = 45/(s(s + 2)(s + 6))t = 0,75 s
f(0.75) =
14.Transformées de LaplaceDifficile
Résolvez y″ + 11y′ + 30y = 0 avec y(0) = 6 et y′(0) = −15 à l’aide de la transformée de Laplace, et trouvez y(0,2).
y″ + 11y′ + 30y = 0y(0) = 6y′(0) = −15t = 0,2 s
y(0.2) =
15.Convolution et systèmes LTIMoyen
L’impulsion x(t) = 10 pour 0 ≤ t < 2 s (nulle ailleurs) est appliquée à un système de réponse h(t) = e−1tu(t). Trouvez la sortie y(2,5) à t = 2,5 s, une fois l’impulsion terminée.
x(t) = 10, 0 ≤ t < 2 sh(t) = e−1tu(t)t = 2,5 s
y(2.5) =
16.Convolution et systèmes LTIDifficile
Un système obéit à y[n] = 0,8·y[n − 1] + 2·x[n] et part du repos. Trouvez l’échantillon h[5] de sa réponse impulsionnelle.
y[n] = 0,8y[n − 1] + 2x[n]x[n] = δ[n]
h[5] =
3Solutions
1.Équations différentielles
Réponse: T = 249,2 °C
Équations différentielles du premier ordre : dy/dt = k·y a pour solution y = y₀·ekt (croissance pour k > 0, décroissance pour k < 0) ; une demi-vie T donne k = ln 2/T. La loi de refroidissement de Newton, dT/dt = −k·(T − Tₐ), donne T(t) = Tₐ + (T₀ − Tₐ)·e−kt. Une équation linéaire y′ + p(x)·y = q(x) se résout avec le facteur intégrant μ = e∫p dx : (μy)′ = μq, donc μy = ∫μq dx + C.
Solution de la loi de refroidissement
T(t) = 20 + 280·e−0,02t
Remplacer t
T(10) = 20 + 280·0,8187 = 249,2 °C
Réponse
T = 249,2 °C
2.Équations différentielles
Réponse: y(0.75) = −1,81
Équations linéaires du second ordre : pour y″ + b·y′ + c·y = 0, essayez y = ert : on obtient l’équation caractéristique r² + b·r + c = 0. Deux racines réelles r₁ ≠ r₂ donnent y = A·er₁t + B·er₂t (suramorti) ; une racine double r donne y = (A + Bt)·ert (amortissement critique) ; des racines complexes −α ± jω donnent y = e−αt(A cos ωt + B sin ωt) (sous-amorti, oscillant). Les valeurs initiales y(0) et y′(0) fixent A et B. Écrite sous la forme y″ + 2ζωny′ + ωn²y = 0, ωn est la pulsation propre et ζ le coefficient d’amortissement : ζ < 1 donne des oscillations. Avec une entrée constante, y″ + by′ + cy = c·K, la solution est la valeur stationnaire K plus la réponse libre.
Matrices : le déterminant d'une matrice 2×2 est ad − bc. Pour une 3×3, développe selon une ligne : chaque coefficient multiplié par le déterminant de la matrice 2×2 qui reste après avoir rayé sa ligne et sa colonne, avec des signes alternés + − +. Les valeurs propres λ vérifient det(A − λI) = 0 ; en 2×2 c'est λ² − (trace)·λ + det = 0.
Ligne fois colonne
c12 = 5·1 + (−2)·5 = −5
Réponse
c12 = −5
4.Matrices et déterminants
Réponse: λ = 6
Matrices : le déterminant d'une matrice 2×2 est ad − bc. Pour une 3×3, développe selon une ligne : chaque coefficient multiplié par le déterminant de la matrice 2×2 qui reste après avoir rayé sa ligne et sa colonne, avec des signes alternés + − +. Les valeurs propres λ vérifient det(A − λI) = 0 ; en 2×2 c'est λ² − (trace)·λ + det = 0.
Équation caractéristique
λ² − (5)λ + (−6) = 0
Résoudre
λ = (5 ± √(49))/2 → λ = −1, 6
λ = 6
Réponse
λ = 6
5.Vecteurs : produit scalaire et produit vectoriel
Réponse: θ = 143,4 °
Dot product:a·b = axbx + ayby + azbz = |a|·|b|·cos θ, a scalar. It gives the angle between two vectors, cos θ = a·b/(|a|·|b|), and it is zero when they are perpendicular. The scalar projection of a onto b is a·b/|b|, and a constant force F moving through a displacement d does work W = F·d.
Produit scalaire
a·b = 3·0 + (−5)·3 + 2·(−3) = −21
Normes
|a| = √(3² + (−5)² + 2²) = 6,164
|b| = √(0² + 3² + (−3)²) = 4,243
Angle
cos θ = −21/(6,164·4,243) = −0,803
θ = arccos(−0,803) = 143,4°
Réponse
θ = 143,4°
6.Vecteurs : produit scalaire et produit vectoriel
Réponse: (a × b)z = −26
Cross product:a × b = (aybz − azby, azbx − axbz, axby − aybx) is perpendicular to both vectors, and |a × b| = |a|·|b|·sin θ is the area of the parallelogram they span (the triangle has half of it). Order matters: b × a = −(a × b). The moment of a force F applied at point A about point O is MO = r × F, with r = A − O.
Produit vectoriel, composante par composante
(a × b)x = aybz − azby = 1·2 − 4·6 = −22
(a × b)y = azbx − axbz = 4·2 − (−4)·2 = 16
(a × b)z = axby − aybx = (−4)·6 − 1·2 = −26
Composante demandée
(a × b)z = −26
Réponse
(a × b)z = −26
Vérification : a × b est perpendiculaire aux deux vecteurs, donc son produit scalaire avec a et avec b est nul.
7.Intégrales
Réponse: I = 12
Intégrales définies : ∫ₐᵇ f(x) dx est l'aire algébrique entre le graphe et l'axe des x de a à b. Trouve une primitive F (avec F′ = f) et calcule F(b) − F(a). Règle des puissances : ∫xⁿ dx = xⁿ⁺¹/(n + 1). L'aire entre deux courbes est ∫(haute − basse) dx.
Changement de variable
u = x² + 1, du = 2x dx
I = ∫15 u1 du
Intégrer et évaluer
I = [u2/2]15 = (25 − 1)/2 = 12
Réponse
I = 12
8.Intégrales
Réponse: I = 0,7183
Intégration par parties : ∫u·dv = u·v − ∫v·du, issue de la règle du produit. Choisissez pour u le facteur qui se simplifie quand on le dérive (une puissance de x, ou ln x), et pour dv le facteur que vous savez intégrer (eˣ, sin x, cos x, ou x lorsque u = ln x). Pour une intégrale définie, évaluez u·v aux deux bornes : ∫ₐᵇ u dv = [u·v]ₐᵇ − ∫ₐᵇ v du. Un facteur x² demande deux étapes.
Première étape
∫x²eˣ dx = x²eˣ − ∫2x eˣ dx
Deuxième étape
∫2x eˣ dx = 2x eˣ − 2eˣ
∫x²eˣ dx = eˣ(x² − 2x + 2)
Limites
2,718·(1) − 1·2 = 0,7183
Réponse
I = 0,7183
9.Analyse à plusieurs variables
Réponse: Duf = 1,333
Gradient : ∇f = (fx, fy) pointe dans la direction de la plus forte montée, et |∇f| est le taux de variation maximal. La dérivée directionnelle suivant une direction v est Duf = ∇f·u avec le vecteur unitaire u = v/|v|. Le plan tangent en (a, b) est z = f(a, b) + fx(a, b)(x − a) + fy(a, b)(y − b) ; il donne l’approximation linéaire de f près de ce point. Dérivation en chaîne : dz/dt = fx·dx/dt + fy·dy/dt.
Gradient au point
fx = 2x/(x² + 2y²) → fx(−1, 0,5) = −1,333
fy = 4y/(x² + 2y²) → fy(−1, 0,5) = 1,333
Vecteur unitaire
|v| = √((−3)² + 0²) = 3
u = (−1, 0)
Produit scalaire
Duf = (−1,333)·(−1) + 1,333·0 = 1,333
Réponse
Duf = 1,333
10.Analyse à plusieurs variables
Réponse: I = 3,976
Points critiques et intégrales doubles : un point critique vérifie fx = fy = 0. Test des dérivées secondes avec D = fxxfyy − fxy² : D > 0 et fxx > 0 est un minimum local, D > 0 et fxx < 0 un maximum local, D < 0 un point selle, et D = 0 ne permet pas de conclure. Pour une intégrale double sur une région générale, prenez comme bornes intérieures les courbes qui délimitent la région ; pour des disques et des secteurs, utilisez les coordonnées polaires, x = r cos θ, y = r sin θ et dA = r dr dθ (n’oubliez pas le facteur r).
Forme polaire
I = ∫0π ∫01,5 r²·r dr dθ
Intégrer
∫ r²·r dr = [r⁴/4]
I = 3,976
Réponse
I = 3,976
11.Suites et séries
Réponse: L = 1,667
Convergence : si les termes ne tendent pas vers 0, la série diverge (test du terme général). Σ1/kᵖ ne converge que pour p > 1, et Σrᵏ que pour |r| < 1. Comparez une série à l’une d’elles quand ses termes se comportent comme celles-ci pour k grand. Test du rapport : L = lim |ak+1/ak| ; test de la racine : L = lim |ak|1/k. L < 1 converge, L > 1 diverge, L = 1 ne décide rien. Une série entière Σck(x − x₀)ᵏ converge pour |x − x₀| < R avec R = lim |ck/ck+1| ; une série alternée dont les termes décroissent vers 0 converge.
Rapport de termes consécutifs
|ak+1/ak| = ((k + 1)/k)³·5/3 → 5/3
Conclusion
L = 1,667 > 1
Réponse
L = 1,667
L > 1, donc la série diverge.
12.Suites et séries
Réponse: P₄(1.5) = 0,2734
Polynômes de Taylor : près de x = a, une fonction régulière est approchée par Pₙ(x) = f(a) + f′(a)(x − a) + f″(a)(x − a)²/2! + … + f⁽ⁿ⁾(a)(x − a)ⁿ/n! ; pour a = 0, c’est le polynôme de Maclaurin. Séries usuelles : eˣ = Σxᵏ/k!, sin x = x − x³/3! + x⁵/5! − …, cos x = 1 − x²/2! + x⁴/4! − …, ln(1 + x) = x − x²/2 + x³/3 − …, 1/(1 − x) = Σxᵏ. Majoration de l’erreur de Lagrange : |f(x) − Pₙ(x)| ≤ M·|x − a|ⁿ⁺¹/(n + 1)!, où M est la plus grande valeur de |f⁽ⁿ⁺¹⁾| entre a et x.
La valeur exacte est f(1,5) = 0,22313 ; l’erreur −0,0503073 diminue quand on prend plus de termes ou un point plus proche de 0.
13.Transformées de Laplace
Réponse: f(0.75) = 2,516
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²Y − 6s − (−15) + 11(sY − 6) + 30Y = 0
Y(s) = (6s + 51)/((s + 5)(s + 6))
Décomposition en éléments simples
A = 21/1 = 21
B = 15/(−1) = −15
Résultat
y(t) = 21·e−5t − 15·e−6t
y(0,2) = 3,208
Réponse
y(0,2) = 3,208
Vérification : y(0) = A + B = 6.
15.Convolution et systèmes LTI
Réponse: y(2.5) = 5,244
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).
Sortie à la fin de l’impulsion
y(2) = (10/1)(1 − e−2) = 8,647
Décroissance après l’impulsion
y(2,5) = 8,647·e−1·0,5 = 5,244
Réponse
y(2,5) = 5,244
16.Convolution et systèmes LTI
Réponse: h[5] = 0,6554
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)|.