У вікні друку виберіть «Зберегти як PDF», щоб зберегти копію на пристрої.
1Формули та методи
Диференціальні рівняння
Експоненціальне зростання/спадання
dy/dt = ky → y = y0ekt
Період напіврозпаду
k = ln 2 / t½
Закон охолодження Ньютона
T(t) = Ta + (T0 − Ta)e−kt
Інтегрувальний множник
y′ + p(x)y = q(x): μ = e∫p dx, (μy)′ = μq
Бак змішування
dm/dt = r·c − (r/V)·m
Матриці та визначники
визначник 2×2
det = ad − bc
Обернена матриця 2×2
A⁻¹ = (1/det)·[d −b; −c a]
Продукт
(AB)ij = Σk AikBkj
Власні значення
det(A − λI) = 0
Правило Крамера
xi = det(Ai)/det(A)
Вектори: скалярний і векторний добутки
Модуль
|a| = √(ax² + ay² + az²)
Одиничний вектор
û = a/|a|
Компоненти за кутом
Fx = F·cos θ, Fy = F·sin θ
Скалярний добуток
a·b = axbx + ayby + azbz = |a||b|cos θ
Скалярна проєкція a на b
a·b/|b|
Робота
W = F·d
Векторний добуток
a × b = (aybz − azby, azbx − axbz, axby − aybx)
Площі
паралелограм |a × b|, трикутник ½|a × b|
Момент відносно O
MO = r × F, r = A − O
Об’єм паралелепіпеда
V = |a·(b × c)|
Метод
Запишіть кожен вектор у компонентах; для сили, заданої модулем і кутом, використовуйте F·cos θ і F·sin θ.
Для довжин і напрямів використовуйте модуль √(x² + y² + z²) та одиничний вектор a/|a|.
Скалярний добуток використовуйте для кутів, проєкцій, перпендикулярності й роботи: a·b = |a||b|cos θ.
Векторний добуток використовуйте для перпендикулярних векторів, площ і моментів: акуратно запишіть три компоненти, зберігаючи порядок a × b.
Інтеграли
Похідна степеня
∫xⁿ dx = xⁿ⁺¹/(n + 1) + C
Визначений інтеграл
∫ab f dx = F(b) − F(a)
Основні інтеграли
∫eˣ dx = eˣ, ∫1/x dx = ln|x|, ∫sin x dx = −cos x
Площа між кривими
∫ab (f − g) dx
Метод
Якщо можливо, перепишіть підінтегральну функцію як суму простих доданків.
Застосуйте степеневе правило: ∫xⁿ dx = xⁿ⁺¹/(n + 1), і табличні інтеграли для ex, 1/x, sin і cos.
Використовуйте підстановку, коли бачите функцію разом з її похідною.
Для визначеного інтеграла обчисліть F(b) − F(a).
Диференціальне числення функцій багатьох змінних
Частинні похідні
fx = ∂f/∂x (y constant), fxy = fyx
Градієнт
∇f = (fx, fy), largest slope |∇f|
Похідна за напрямом
Duf = ∇f·u, u = v/|v|
Дотична площина
z = f(a, b) + fx(x − a) + fy(y − b)
Ланцюгове правило
dz/dt = fx·dx/dt + fy·dy/dt
Тест другої похідної
D = fxxfyy − fxy²: D > 0, fxx > 0 minimum; D > 0, fxx < 0 maximum; D < 0 saddle
Подвійний інтеграл
∬R f dA = ∫ab∫g(x)h(x) f dy dx
Полярні координати
x = r cos θ, y = r sin θ, dA = r dr dθ
Метод
Для частинної похідної трактуйте іншу змінну як сталу й користуйтеся правилами для однієї змінної.
Для задач про градієнт обчисліть ∇f = (fx, fy) у точці; використовуйте |∇f| для найкрутішого нахилу та ∇f·u з одиничним вектором u для похідної за напрямом.
Для екстремумів розв’яжіть fx = 0 і fy = 0 разом, а потім класифікуйте кожну точку за D = fxx·fyy − fxy².
Для подвійного інтеграла накресліть область, оберіть порядок (або полярні координати для кругів), щоб внутрішні межі були обмежувальними кривими, і спершу інтегруйте всередині.
Послідовності та ряди
Арифметична послідовність
aₙ = a₁ + (n − 1)d, Sₙ = n(a₁ + aₙ)/2
Геометрична послідовність
aₙ = a₁rⁿ⁻¹, Sₙ = a₁(1 − rⁿ)/(1 − r)
Нескінченний геометричний ряд
S = a₁/(1 − r), |r| < 1
p-ряд
Σ1/kᵖ converges only for p > 1
Ознака Даламбера
L = lim |ak+1/ak|: L < 1 converges, L > 1 diverges
Радіус збіжності
R = lim |ck/ck+1|
Многочлен Тейлора
Pₙ(x) = Σ f⁽ᵏ⁾(a)(x − a)ᵏ/k!
Ряд Маклорена
eˣ = Σxᵏ/k!, sin x = x − x³/3! + …, cos x = 1 − x²/2! + …, ln(1 + x) = x − x²/2 + …
Оцінка похибки Лагранжа
|Rₙ| ≤ M·|x − a|ⁿ⁺¹/(n + 1)!
Метод
Визначте тип: арифметична (стала різниця), геометрична (сталий знаменник), стандартний ряд, степеневий ряд чи многочлен Тейлора.
Для сум користуйтеся Sₙ = n(a₁ + aₙ)/2, Sₙ = a₁(1 − rⁿ)/(1 − r) або S = a₁/(1 − r) за |r| < 1, де a₁ — перший член суми.
Для збіжності: спершу ознака n-го члена, потім порівняння з Σ1/kᵖ або геометричним рядом, а ознаку Даламбера чи Коші використовуйте для факторіалів і k-х степенів.
Для многочленів Тейлора випишіть f(a), f′(a), f″(a), …, поділіть k-ту похідну на k! і оцініть похибку за M·|x − a|ⁿ⁺¹/(n + 1)!.
Накресліть обидва сигнали й позначте, де кожен з них ненульовий.
Віддзеркаліть h і зсувайте її: для кожного моменту виходу знайдіть діапазон, де обидва сигнали перекриваються.
Перемножуйте й додавайте (дискретно) або перемножуйте й інтегруйте (неперервно) лише на цьому перекритті.
Перевірте результат: вихід починається в суму моментів початку, закінчується в суму моментів завершення, а його площа чи сума дорівнює добутку площ чи сум входів.
2Задачі для практики
1.Диференціальні рівнянняСередня
Розв’яжіть y″ + 3y′ + 2y = 0 при y(0) = 3 і y′(0) = -4. Знайдіть y(2).
y″ + 3y′ + 2y = 0y(0) = 3y′(0) = -4
y(2) =
2.Диференціальні рівнянняСкладна
Система підпорядковується y″ + 4y′ + 29y = 58 для t ≥ 0, починаючи зі стану спокою: y(0) = 0, y′(0) = 0. Знайдіть y(1,5); кути в радіанах.
y″ + 4y′ + 29y = 58y(0) = 0y′(0) = 0
y(1.5) =
3.Матриці та визначникиСередня
Нехай C = AB для поданих матриць. Знайдіть елемент c21 (рядок 2, стовпець 1).
A = [1 −3; 4 4]B = [6 5; 2 4]
c21 =
4.Матриці та визначникиСкладна
Знайдіть більше власне значення λ поданої матриці A.
A = [6 −4; 5 −3]
λ =
5.Вектори: скалярний і векторний добуткиСередня
Знайдіть значення k, за якого вектори a = (0, −5, 5) і b = (−6, 1, k) перпендикулярні.
a = (0, −5, 5)b = (−6, 1, k)
k =
6.Вектори: скалярний і векторний добуткиСкладна
Знайдіть площу A паралелограма, натягнутого на a = (−1, 6, −4) і b = (5, 2, 0).
a = (−1, 6, −4)b = (5, 2, 0)
A =
7.ІнтегралиСередня
Обчисліть I = ∫01 e2x dx.
f(x) = e2x[0, 1]
I =
8.ІнтегралиСкладна
Область між прямою y = 2x і параболою y = 2x² обертається навколо осі x. Знайдіть об’єм V цього тіла.
y = 2xy = 2x²
V =
9.Диференціальне числення функцій багатьох зміннихСередня
Для f(x, y) = −2x³e−y знайдіть похідну за напрямом Duf у (2, −1) у напрямку v = (1, −1).
f(x, y) = −2x³e−yP = (2, −1)v = (1, −1)
Duf =
10.Диференціальне числення функцій багатьох зміннихСкладна
Знайдіть критичні точки f(x, y) = x³ + y³ − 6xy, класифікуйте їх і дайте значення fmin функції f у її локальному мінімумі.
f(x, y) = x³ + y³ − 6xy
fmin =
11.Послідовності та рядиСередня
Рівно один із чотирьох рядів розбігається. Клікніть його.
A: Σ (1/2)ᵏB: Σ (−1)ᵏ/√kC: Σ 1/√kD: Σ 10ᵏ/k!
1A: Σ (1/2)ᵏ
2B: Σ (−1)ᵏ/√k
3C: Σ 1/√k
4D: Σ 10ᵏ/k!
12.Послідовності та рядиСкладна
Скористайтеся многочленом Маклорена степеня 3 для f(x) = ln(1 − x), щоб наблизити f(0,5). Дайте P₃(0,5).
f(x) = ln(1 − x)n = 3x = 0,5
P₃(0.5) =
13.Перетворення ЛапласаСередня
Знайдіть f(t) для F(s) = (5s + 28)/(s² + 10s + 25) і обчисліть f(0,4).
F(s) = (5s + 28)/(s² + 10s + 25)t = 0,4 s
f(0.4) =
14.Перетворення ЛапласаСкладна
Система підпорядковується рівнянню y″ + 7y′ + 10y = 5·u(t). Вона починає зі спокою (y(0) = 0, y′(0) = 0), а вхід — сходинка u(t) = 2 при t ≥ 0. За допомогою перетворення Лапласа знайдіть вихід y(0,2).
Вхід x(t) = 1·e−2tu(t) подано на систему з h(t) = e−1tu(t). Знайдіть вихід y(2) при t = 2 s.
x(t) = 1·e−2tu(t)h(t) = e−1tu(t)t = 2 s
y(2) =
16.Згортка та лінійні стаціонарні системиСкладна
Система описується рівнянням y′(t) + 2·y(t) = 3·x(t). Знайдіть її імпульсну характеристику при t = 0,2 s, h(0,2).
y′ + 2y = 3xt = 0,2 s
h(0.2) =
3Розчини
1.Диференціальні рівняння
Відповідь: y(2) = 0,289
Лінійні рівняння другого порядку: для y″ + b·y′ + c·y = 0 шукайте y = ert: це дає характеристичне рівняння r² + b·r + c = 0. Два дійсні корені r₁ ≠ r₂ дають y = A·er₁t + B·er₂t (надкритичне демпфування); подвійний корінь r дає y = (A + Bt)·ert (критичне демпфування); комплексні корені −α ± jω дають y = e−αt(A cos ωt + B sin ωt) (докритичне демпфування, коливання). Початкові значення y(0) і y′(0) визначають A і B. У вигляді y″ + 2ζωny′ + ωn²y = 0 величина ωn — власна частота, а ζ — коефіцієнт демпфування: при ζ < 1 є коливання. При сталому вході y″ + by′ + cy = c·K розв’язок — усталене значення K плюс вільна складова.
Характеристичні корені
(r + 1)(r + 2) = 0 → r = −1, −2
Сталі за початковими значеннями
A + B = 3, −1A − 2B = -4
A = 2, B = 1
Обчисліть
y(2) = 2·e−2 + 1·e−4 = 0,289
Відповідь
y(2) = 0,289
2.Диференціальні рівняння
Відповідь: y(1.5) = 1,928
Лінійні рівняння другого порядку: для y″ + b·y′ + c·y = 0 шукайте y = ert: це дає характеристичне рівняння r² + b·r + c = 0. Два дійсні корені r₁ ≠ r₂ дають y = A·er₁t + B·er₂t (надкритичне демпфування); подвійний корінь r дає y = (A + Bt)·ert (критичне демпфування); комплексні корені −α ± jω дають y = e−αt(A cos ωt + B sin ωt) (докритичне демпфування, коливання). Початкові значення y(0) і y′(0) визначають A і B. У вигляді y″ + 2ζωny′ + ωn²y = 0 величина ωn — власна частота, а ζ — коефіцієнт демпфування: при ζ < 1 є коливання. При сталому вході y″ + by′ + cy = c·K розв’язок — усталене значення K плюс вільна складова.
Матриці: визначник матриці 2×2 дорівнює ad − bc. Для 3×3 розкладіть за рядком: кожен елемент множиться на визначник матриці 2×2, що залишається після викреслення його рядка і стовпця, зі знаками, що чергуються: + − +. Власні значення λ задовольняють det(A − λI) = 0; для 2×2 це λ² − (слід)·λ + det = 0.
Рядок на стовпець
c21 = 4·6 + 4·2 = 32
Відповідь
c21 = 32
4.Матриці та визначники
Відповідь: λ = 2
Матриці: визначник матриці 2×2 дорівнює ad − bc. Для 3×3 розкладіть за рядком: кожен елемент множиться на визначник матриці 2×2, що залишається після викреслення його рядка і стовпця, зі знаками, що чергуються: + − +. Власні значення λ задовольняють det(A − λI) = 0; для 2×2 це λ² − (слід)·λ + det = 0.
Характеристичне рівняння
λ² − (3)λ + (2) = 0
Розв’яжіть
λ = (3 ± √(1))/2 → λ = 1, 2
λ = 2
Відповідь
λ = 2
5.Вектори: скалярний і векторний добутки
Відповідь: k = 1
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.
Скалярний добуток дорівнює нулю
a·b = 0·(−6) + (−5)·1 + 5·k = 0
−5 + 5·k = 0
Розв’яжіть щодо k
k = −(−5)/5 = 1
Відповідь
k = 1
6.Вектори: скалярний і векторний добутки
Відповідь: A = 38,57
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.
Векторний добуток
(a × b)x = aybz − azby = 6·0 − (−4)·2 = 8
(a × b)y = azbx − axbz = (−4)·5 − (−1)·0 = −20
(a × b)z = axby − aybx = (−1)·2 − 6·5 = −32
Площа
A = |a × b| = √(8² + (−20)² + (−32)²) = 38,57
Відповідь
A = 38,57
7.Інтеграли
Відповідь: I = 3,195
Визначені інтеграли: ∫ₐᵇ f(x) dx — це площа зі знаком між графіком і віссю x від a до b. Знайдіть первісну F (F′ = f) і обчисліть F(b) − F(a). Степеневе правило: ∫xⁿ dx = xⁿ⁺¹/(n + 1). Площа між двома кривими — це ∫(верхня − нижня) dx.
Первісна
F(x) = e2x/2
Підставте межі
I = (e2 − 1)/2 = (7,389 − 1)/2 = 3,195
Відповідь
I = 3,195
8.Інтеграли
Відповідь: V = 1,676
Об’єми тіл обертання: обертання області під y = f(x) від x = a до x = b навколо осі x дає тіло, поперечний переріз якого в точці x — круг радіуса f(x). Сума кругів дає V = π∫ₐᵇ f(x)² dx. Якщо область лежить між зовнішньою кривою R(x) і внутрішньою r(x), кожен зріз — шайба: V = π∫ₐᵇ (R(x)² − r(x)²) dx. Радіуси підносьте до квадрата до віднімання.
Точки перетину
2x = 2x² → x = 0, x = 1
Метод шайб
V = π∫₀1 ((2x)² − (2x²)²) dx
Інтегруйте
V = π(4·1³/3 − 4·1⁵/5) = π·0,5333 = 1,676
Відповідь
V = 1,676
9.Диференціальне числення функцій багатьох змінних
Відповідь: Duf = −76,88
Градієнт: ∇f = (fx, fy) вказує напрямок найшвидшого зростання, а |∇f| — найбільша швидкість зміни. Похідна за напрямом v дорівнює Duf = ∇f·u з одиничним вектором u = v/|v|. Дотична площина в точці (a, b): z = f(a, b) + fx(a, b)(x − a) + fy(a, b)(y − b); вона дає лінійне наближення f поблизу цієї точки. Правило ланцюга: dz/dt = fx·dx/dt + fy·dy/dt.
Градієнт у точці
fx = −6x²e−y → fx(2, −1) = −65,24
fy = 2x³e−y → fy(2, −1) = 43,49
Одиничний вектор
|v| = √(1² + (−1)²) = 1,414
u = (0,7071, −0,7071)
Скалярний добуток
Duf = (−65,24)·0,7071 + 43,49·(−0,7071) = −76,88
Відповідь
Duf = −76,88
10.Диференціальне числення функцій багатьох змінних
Відповідь: fmin = −8
Критичні точки та подвійні інтеграли: критична точка має fx = fy = 0. Тест другої похідної з D = fxxfyy − fxy²: D > 0 і fxx > 0 — локальний мінімум, D > 0 і fxx < 0 — локальний максимум, D < 0 — сідлова точка, а D = 0 нічого не вирішує. Для подвійного інтеграла по довільній області нехай внутрішні межі будуть кривими, що обмежують область; для кругів і секторів використовуйте полярні координати, x = r cos θ, y = r sin θ і dA = r dr dθ (не забудьте множник r).
Критичні точки
Розв’яжіть fx = 0 і fy = 0 разом.
fx = 3x² − 6y = 0
fy = 3y² − 6x = 0
(0, 0), (2, 2)
Тест другої похідної
(0, 0): fxx = 0, fyy = 0, fxy = −6, D = −36 → сідлова точка
Збіжність: якщо члени не прямують до 0, ряд розбігається (необхідна умова збіжності, ознака n-го члена). Σ1/kᵖ збігається лише за p > 1, а Σrᵏ лише за |r| < 1. Порівнюйте ряд з одним із них, коли його члени для великих k поводяться як вони. Ознака Даламбера: L = lim |ak+1/ak|; ознака Коші: L = lim |ak|1/k. L < 1 — збігається, L > 1 — розбігається, L = 1 нічого не вирішує. Степеневий ряд Σck(x − x₀)ᵏ збігається для |x − x₀| < R, де R = lim |ck/ck+1|; знакозмінний ряд із членами, що спадають до 0, збігається.
Перевірте кожен ряд
A: Σ (1/2)ᵏ → Геометричний з r = 1/2, |r| < 1: збігається.
B: Σ (−1)ᵏ/√k → Знакозмінний, із 1/√k, що спадає до 0: збігається (умовно).
C: Σ 1/√k → p-ряд з p = 1/2 ≤ 1: розбігається.
D: Σ 10ᵏ/k! → Ознака Даламбера: L = lim 10/(k + 1) = 0 < 1, збігається.
Відповідь
C: Σ 1/√k
Відповідь
C: Σ 1/√k
12.Послідовності та ряди
Відповідь: P₃(0.5) = −0,6667
Многочлени Тейлора: поблизу x = a гладка функція наближається многочленом Pₙ(x) = f(a) + f′(a)(x − a) + f″(a)(x − a)²/2! + … + f⁽ⁿ⁾(a)(x − a)ⁿ/n!; для a = 0 це многочлен Маклорена. Стандартні ряди: 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ᵏ. Оцінка похибки Лагранжа: |f(x) − Pₙ(x)| ≤ M·|x − a|ⁿ⁺¹/(n + 1)!, де M — найбільше значення |f⁽ⁿ⁺¹⁾| між a і x.
Многочлен Маклорена
P₃(x) = −x − 0,5x² − 0,3333x³
Підставте
P₃(0,5) = −0,5 − 0,125 − 0,0416667 = −0,6667
Відповідь
P₃(0,5) = −0,6667
Точне значення дорівнює f(0,5) = −0,693147; похибка −0,0264805 зменшується, якщо взяти більше членів або точку ближче до 0.
13.Перетворення Лапласа
Відповідь: f(0.4) = 0,8391
Перетворення Лапласа: F(s) = ∫₀∞ f(t)·e−st dt перетворює похідні на множення на s: ℒ{f′} = sF(s) − f(0), ℒ{f″} = s²F(s) − s·f(0) − f′(0). Стандартні пари: 1 ↔ 1/s, e−at ↔ 1/(s + a), t·e−at ↔ 1/(s + a)², sin ωt ↔ ω/(s² + ω²), cos ωt ↔ s/(s² + ω²), а зсув s → s + σ множить оригінал на e−σt. Щоб знайти оригінал раціональної F(s), розкладіть її на простіші дроби: коефіцієнт простого полюса p дорівнює A = [(s − p)·F(s)] при s = p. Теорема про кінцеве значення: f(∞) = lims→0 s·F(s), справедлива лише тоді, коли всі полюси s·F(s) лежать у лівій півплощині. Теорема про початкове значення: f(0⁺) = lims→∞ s·F(s).
Розкладіть знаменник на множники
s² + 10s + 25 = (s + 5)²
Простіші дроби
5s + 28 = 5(s + 5) + 3
F(s) = 5/(s + 5) + 3/(s + 5)²
Обернене перетворення
f(t) = 5·e−5t + 3·t·e−5t
f(0,4) = (5 + 3·0,4)·e−5·0,4 = 0,8391
Відповідь
f(0,4) = 0,8391
14.Перетворення Лапласа
Відповідь: y(0.2) = 0,1281
Перетворення Лапласа: F(s) = ∫₀∞ f(t)·e−st dt перетворює похідні на множення на s: ℒ{f′} = sF(s) − f(0), ℒ{f″} = s²F(s) − s·f(0) − f′(0). Стандартні пари: 1 ↔ 1/s, e−at ↔ 1/(s + a), t·e−at ↔ 1/(s + a)², sin ωt ↔ ω/(s² + ω²), cos ωt ↔ s/(s² + ω²), а зсув s → s + σ множить оригінал на e−σt. Щоб знайти оригінал раціональної F(s), розкладіть її на простіші дроби: коефіцієнт простого полюса p дорівнює A = [(s − p)·F(s)] при s = p. Теорема про кінцеве значення: f(∞) = lims→0 s·F(s), справедлива лише тоді, коли всі полюси s·F(s) лежать у лівій півплощині. Теорема про початкове значення: f(0⁺) = lims→∞ s·F(s).
Перетворіть рівняння
Y(s) = 10/(s(s + 2)(s + 5))
Простіші дроби
A = 10/(2·5) = 1
B = 10/((−2)·3) = −1,667
C = 10/((−5)·(−3)) = 0,6667
Виведення
y(t) = 1 − 1,667·e−2t + 0,6667·e−5t
y(0,2) = 0,1281
Відповідь
y(0,2) = 0,1281
15.Згортка та лінійні стаціонарні системи
Відповідь: y(2) = 0,117
Неперервна згортка: y(t) = ∫ x(τ)·h(t − τ) dτ. Відобразіть h дзеркально, посуньте її в положення t та проінтегруйте добуток частин, що перекриваються. Для причинних сигналів інтеграл береться від 0 до t. Корисні результати: площа y дорівнює площі x, помноженій на площу h; якщо x триває від t₁ до t₂, а h — від t₃ до t₄, то y триває від t₁ + t₃ до t₂ + t₄. Для h(t) = e−atu(t) і східчастого сигналу висотою A маємо y(t) = (A/a)(1 − e−at).
Імпульсна характеристика та ряд Фур’є: імпульсна характеристика h(t) повністю описує лінійну стаціонарну систему: для y′ + a·y = b·x вона має вигляд h(t) = b·e−atu(t), а для y[n] = a·y[n − 1] + b·x[n] — h[n] = b·aⁿu[n]. Синусоїда cos ωt проходить із коефіцієнтом підсилення |H(jω)| і фазою ∠H(jω), де H — перетворення Фур’є функції h. Періодичний сигнал є сумою гармонік на частотах k/T: прямокутна хвиля ±A має лише непарні гармоніки з амплітудами 4A/(kπ), трикутна хвиля з піком A має 8A/(π²k²), а послідовність імпульсів висотою A з коефіцієнтом заповнення D має середнє значення A·D і гармоніки (2A/(kπ))·|sin(kπD)|.