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é à 460 °C refroidit dans un atelier à 23 °C. Elle suit la loi de refroidissement de Newton, dT/dt = −k·(T − Tₐ). Avec k = 0,04 min⁻¹, trouvez sa température T au bout de 39 min.
T0 = 460 °CTa = 23 °Ck = 0,04 min⁻¹t = 39 min
T = °C
2.Équations différentiellesDifficile
Un système obéit à y″ + 4y′ + 13y = 26 pour t ≥ 0, à partir du repos : y(0) = 0, y′(0) = 0. Trouvez la valeur de pic ymax de sa réponse indicielle.
y″ + 4y′ + 13y = 26y(0) = 0y′(0) = 0
ymax =
3.Matrices et déterminantsMoyen
Calcule le déterminant det A de la matrice 3×3 A représentée.
A = [−3 1 5; 2 3 1; 5 4 −1]
det A =
4.Matrices et déterminantsDifficile
Trouve la plus grande valeur propre λ de la matrice A représentée.
A = [−2 −1; 7 6]
λ =
5.Vecteurs : produit scalaire et produit vectorielMoyen
Trouvez l’angle θ entre les vecteurs a = (6, 0, 3) et b = (1, 4, 3), en degrés.
a = (6, 0, 3)b = (1, 4, 3)
θ = °
6.Vecteurs : produit scalaire et produit vectorielDifficile
Le vecteur unitaire n̂ est perpendiculaire à a = (−3, −3, −2) et à b = (0, −5, −1) et a une composante z positive. Trouvez sa composante y n̂y.
a = (−3, −3, −2)b = (0, −5, −1)n̂z > 0
n̂y =
7.IntégralesMoyen
Calcule I = ∫03 ex dx.
f(x) = ex[0, 3]
I =
8.IntégralesDifficile
Calcule l'aire A du domaine entre la parabole y = 1 − x² et l'axe des x.
y = 1 − x²
A =
9.Analyse à plusieurs variablesMoyen
z = 2x³e3y avec x = t − 1 et y = −t². Utilisez la dérivation en chaîne pour trouver dz/dt en t = −1.
z = 2x³e3yx = t − 1y = −t²t = −1
dz/dt =
10.Analyse à plusieurs variablesDifficile
f(x, y) = x⁴ + y⁴ − 16xy a un point critique en (2, 2). Calculez-y le discriminant D = fxxfyy − fxy².
f(x, y) = x⁴ + y⁴ − 16xyP = (2, 2)
D =
11.Suites et sériesMoyen
Appliquez le test du rapport à Σk=1∞ k·3ᵏ/5ᵏ. Trouvez L = lim |ak+1/ak|.
ak = k·3ᵏ/5ᵏ
L =
12.Suites et sériesDifficile
Utilisez le polynôme de Maclaurin de degré 5 de f(x) = sin(3x) pour approcher f(0,8). Donnez P₅(0,8).
f(x) = sin(3x)n = 5x = 0,8
P₅(0.8) =
13.Transformées de LaplaceMoyen
F(s) = (5s + 28)/(s² + 14s + 48). Sa transformée inverse a la forme f(t) = A·e−6t + B·e−8t. Trouvez le coefficient A par décomposition en éléments simples.
F(s) = (5s + 28)/(s² + 14s + 48)
A =
14.Transformées de LaplaceDifficile
Un système est décrit par l’équation différentielle y″ + 4y′ + 5y = 2u′ + 17u, d’entrée u(t) et de sortie y(t). Trouvez sa fonction de transfert G(s) = Y(s)/U(s) (conditions initiales nulles) et calculez la réponse impulsionnelle g(t) à t = 1 s. Trouvez g(1).
y″ + 4y′ + 5y = 2u′ + 17ut = 1 s
g(1) =
15.Convolution et systèmes LTIMoyen
L’entrée x(t) = 2·e−4tu(t) est appliquée à un système de réponse h(t) = e−2tu(t). Trouvez la sortie y(1) à t = 1 s.
x(t) = 2·e−4tu(t)h(t) = e−2tu(t)t = 1 s
y(1) =
16.Convolution et systèmes LTIDifficile
Un système obéit à y[n] = 0,6·y[n − 1] + 2·x[n] et part du repos. L’entrée est l’échelon unité x[n] = u[n]. Trouvez y[4].
y[n] = 0,6y[n − 1] + 2x[n]x[n] = u[n]
y[4] =
3Solutions
1.Équations différentielles
Réponse: T = 114,8 °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) = 23 + 437·e−0,04t
Remplacer t
T(39) = 23 + 437·0,2101 = 114,8 °C
Réponse
T = 114,8 °C
2.Équations différentielles
Réponse: ymax = 2,246
É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.
Valeur stationnaire et racines
yss = 2, r = −2 ± 3j
Constantes à partir des valeurs initiales
y(0) = 2 + A = 0 → A = −2
y′(0) = −2A + 3B = 0 → B = −1,333
Pic
t = π/3 = 1,047
ymax = 2·(1 + e−2·π/3) = 2,246
Réponse
ymax = 2,246
3.Matrices et déterminants
Réponse: det A = −7
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.
Mineurs 2×2 de la première ligne
M₁₁ = −7, M₁₂ = −7, M₁₃ = −7
Développement par cofacteurs
det A = (−3)·(−7) − 1·(−7) + 5·(−7) = −7
Réponse
det A = −7
4.Matrices et déterminants
Réponse: λ = 5
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
λ² − (4)λ + (−5) = 0
Résoudre
λ = (4 ± √(36))/2 → λ = −1, 5
λ = 5
Réponse
λ = 5
5.Vecteurs : produit scalaire et produit vectoriel
Réponse: θ = 63,99 °
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 = 6·1 + 0·4 + 3·3 = 15
Normes
|a| = √(6² + 0² + 3²) = 6,708
|b| = √(1² + 4² + 3²) = 5,099
Angle
cos θ = 15/(6,708·5,099) = 0,4385
θ = arccos(0,4385) = 63,99°
Réponse
θ = 63,99°
6.Vecteurs : produit scalaire et produit vectoriel
Réponse: n̂y = −0,1783
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.
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.
Primitive
F(x) = ex/1
Évaluer aux bornes
I = (e3 − 1)/1 = (20,09 − 1)/1 = 19,09
Réponse
I = 19,09
8.Intégrales
Réponse: A = 1,333
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.
Zéros
1 − x² = 0 → x = ±1
Intégrer
A = [1x − x³/3]−11 = 0,6667 − (−0,6667) = 1,333
Réponse
A = 1,333
9.Analyse à plusieurs variables
Réponse: dz/dt = −3,585
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.
Point et taux de variation
x(−1) = −2, y(−1) = −1
dx/dt = 1, dy/dt = −2t = 2
Dérivées partielles
fx = 6x²e3y = 1,195
fy = 6x³e3y = −2,39
Dérivation en chaîne
dz/dt = 1,195·1 + (−2,39)·2 = −3,585
Réponse
dz/dt = −3,585
10.Analyse à plusieurs variables
Réponse: D = 2048
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).
Dérivées secondes
fxx = 48, fyy = 48, fxy = −16
Discriminant
D = 48·48 − (−16)² = 2048
→ minimum local
Réponse
D = 2048
11.Suites et séries
Réponse: L = 0,6
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)·3/5 → 3/5
Conclusion
L = 0,6 < 1
Réponse
L = 0,6
L < 1, donc la série converge.
12.Suites et séries
Réponse: P₅(0.8) = 0,7596
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.
Polynôme de Maclaurin
P₅(x) = 3x − 4,5x³ + 2,025x⁵
Remplacer
P₅(0,8) = 2,4 − 2,304 + 0,663552 = 0,7596
Réponse
P₅(0,8) = 0,7596
La valeur exacte est f(0,8) = 0,675463 ; l’erreur −0,0840888 diminue quand on prend plus de termes ou un point plus proche de 0.
13.Transformées de Laplace
Réponse: A = −1
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² + 14s + 48 = (s + 6)(s + 8)
F(s) = A/(s + 6) + B/(s + 8)
Méthode des résidus
Multipliez par (s + 6) et posez s = −6 : le terme en B disparaît.
A = (5s + 28)/(s + 8) |s = −6 = −2/2 = −1
Réponse
A = −1
(B = 6, donc f(t) = −e−6t + 6·e−8t)
14.Transformées de Laplace
Réponse: g(1) = 1,627
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).
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)|.