В окне печати выберите «Сохранить как 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.Дифференциальные уравненияСредняя
Банку газировки температурой 6 °C достали из холодильника в комнату при 23 °C. Процесс подчиняется закону охлаждения Ньютона dT/dt = −k·(T − Tₐ). При k = 0,09 min⁻¹ найдите температуру T через 18 min.
T0 = 6 °CTa = 23 °Ck = 0,09 min⁻¹t = 18 min
T = °C
2.Дифференциальные уравненияСложная
Система подчиняется уравнению y″ + 2y′ + 5y = 50 при t ≥ 0, начиная из состояния покоя: y(0) = 0, y′(0) = 0. Найдите y(1,5); углы в радианах.
y″ + 2y′ + 5y = 50y(0) = 0y′(0) = 0
y(1.5) =
3.Матрицы и определителиСредняя
Пусть C = AB для показанных матриц. Найдите элемент c11 (строка 1, столбец 1).
A = [−5 6; 5 1]B = [−5 −5; 0 −5]
c11 =
4.Матрицы и определителиСложная
Решите систему −x + 2y = 4 и 6x − 6y = −6. Найдите x.
−x + 2y = 46x − 6y = −6
x =
5.Векторы: скалярное и векторное произведенияСредняя
Найдите скалярную проекцию compb a вектора a = (2, 5, −4) на b = (0, −5, 1), то есть компоненту a в направлении b.
a = (2, 5, −4)b = (0, −5, 1)
compb a =
6.Векторы: скалярное и векторное произведенияСложная
Найдите площадь A треугольника с вершинами A(−4, −1, −3), B(2, 1, 1) и C(−2, −1, 4). Координаты даны в метрах.
A(−4, −1, −3)B(2, 1, 1)C(−2, −1, 4)
A = m²
7.ИнтегралыСредняя
Вычислите I = ∫₀π x·cos x dx (x в радианах).
∫₀π x·cos x dx
I =
8.ИнтегралыСложная
Найдите площадь A области между параболой y = 4 − x² и осью x.
y = 4 − x²
A =
9.Математический анализ функций нескольких переменныхСредняя
Для f(x, y) = sin(3x − 3y) найдите наибольшую скорость изменения f в точке (2, −1), то есть |∇f|.
f(x, y) = sin(3x − 3y)P = (2, −1)
|∇f| =
10.Математический анализ функций нескольких переменныхСложная
Интеграл I = ∫01 ∫y1 1/(1 + x²)² dx dy нельзя вычислить в таком порядке. Измените порядок интегрирования и вычислите его.
f = 1/(1 + x²)²0 ≤ y ≤ x ≤ 1
I =
11.Последовательности и рядыСредняя
Ровно один из четырёх рядов расходится. Нажмите на него.
A: Σ 3/k²B: Σ 1/(3ᵏ + k)C: Σ (−1)ᵏ/kD: Σ (7/4)ᵏ
1A: Σ 3/k²
2B: Σ 1/(3ᵏ + k)
3C: Σ (−1)ᵏ/k
4D: Σ (7/4)ᵏ
12.Последовательности и рядыСложная
Узнайте ряд и найдите его сумму: S = Σk=0∞ 2ᵏ/k!.
ak = 2ᵏ/k!
S =
13.Преобразования ЛапласаСредняя
F(s) = (2s + 27)/(s² + 3s + 2). Её оригинал имеет вид f(t) = A·e−t + B·e−2t. Найдите коэффициент A разложением на простые дроби.
F(s) = (2s + 27)/(s² + 3s + 2)
A =
14.Преобразования ЛапласаСложная
Система описывается дифференциальным уравнением y″ + 6y′ + 8y = 2u′ + 7u, где u(t) — вход, y(t) — выход. Найдите её передаточную функцию G(s) = Y(s)/U(s) (при нулевых начальных условиях) и вычислите импульсную характеристику g(t) при t = 0,25 s. Найдите g(0,25).
y″ + 6y′ + 8y = 2u′ + 7ut = 0,25 s
g(0.25) =
15.Свёртка и линейные стационарные системыСредняя
Система имеет импульсную характеристику h(t) = e−2tu(t). На входе ступенька x(t) = 1·u(t). Найдите выход y(0,2) при t = 0,2 s.
h(t) = e−2tu(t)x(t) = 1·u(t)t = 0,2 s
y(0.2) =
16.Свёртка и линейные стационарные системыСложная
Треугольная волна линейно нарастает и спадает между −3 и +3 s периодом T = 5 ms. Найдите амплитуду c5 её фурье-компоненты на частоте 1000 Hz.
±3T = 5 msf = 1000 Hz
c5 =
3Растворы
1.Дифференциальные уравнения
Ответ: T = 19,64 °C
Дифференциальные уравнения первого порядка: dy/dt = k·y имеет решение y = y₀·ekt (рост при k > 0, убывание при k < 0); период полураспада T даёт k = ln 2/T. Закон охлаждения Ньютона dT/dt = −k·(T − Tₐ) даёт T(t) = Tₐ + (T₀ − Tₐ)·e−kt. Линейное уравнение y′ + p(x)·y = q(x) решается с интегрирующим множителем μ = e∫p dx: (μy)′ = μq, поэтому μy = ∫μq dx + C.
Решение закона охлаждения
T(t) = 23 + (−17)·e−0,09t
Подставьте t
T(18) = 23 + (−17)·0,1979 = 19,64 °C
Ответ
T = 19,64 °C
2.Дифференциальные уравнения
Ответ: y(1.5) = 12,05
Линейные уравнения второго порядка: для 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.
Строка на столбец
c11 = (−5)·(−5) + 6·0 = 25
Ответ
c11 = 25
4.Матрицы и определители
Ответ: x = 2
Матрицы: определитель матрицы 2×2 равен ad − bc. Для 3×3 разложите по строке: каждый элемент умножается на определитель матрицы 2×2, оставшейся после вычёркивания его строки и столбца, с чередующимися знаками + − +. Собственные значения λ удовлетворяют det(A − λI) = 0; для 2×2 это λ² − (след)·λ + det = 0.
Определитель матрицы коэффициентов
det A = (−1)·(−6) − 2·6 = −6
Замените столбец x
det Aₓ = 4·(−6) − 2·(−6) = −12
x = −12/(−6) = 2
Ответ
x = 2
y = 3
5.Векторы: скалярное и векторное произведения
Ответ: compb a = −5,687
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 = 2·0 + 5·(−5) + (−4)·1 = −29
Длина b
|b| = √(0² + (−5)² + 1²) = 5,099
Скалярная проекция
compb a = a·b/|b| = −29/5,099 = −5,687
Ответ
compb a = −5,687
Отрицательная: у a есть компонента, направленная против b.
6.Векторы: скалярное и векторное произведения
Ответ: A = 18,49 m²
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.
Векторы сторон
AB = B − A = (6, 2, 4)
AC = C − A = (2, 0, 7)
Векторное произведение
(AB × AC)x = AByACz − ABzACy = 2·7 − 4·0 = 14
(AB × AC)y = ABzACx − ABxACz = 4·2 − 6·7 = −34
(AB × AC)z = ABxACy − AByACx = 6·0 − 2·2 = −4
|AB × AC| = √(14² + (−34)² + (−4)²) = 36,99
Площадь треугольника
A = ½·36,99 = 18,49 m²
Ответ
A = 18,49 m²
7.Интегралы
Ответ: I = −2
Интегрирование по частям: ∫u·dv = u·v − ∫v·du, из правила произведения. Выбирайте u как множитель, который упрощается при дифференцировании (степень x или ln x), а dv — как множитель, который можно проинтегрировать (eˣ, sin x, cos x или x, когда u = ln x). Для определённого интеграла вычисляйте u·v на обоих пределах: ∫ₐᵇ u dv = [u·v]ₐᵇ − ∫ₐᵇ v du. Множитель x² требует двух шагов.
Выберите u и dv
u = x, du = dx; dv = cos x dx, v = sin x
Первообразная
∫x cos x dx = x sin x − ∫sin x dx = x sin x + cos x
Пределы
(π·sin π + cos π) − (0 + 1) = −2
Ответ
I = −2
8.Интегралы
Ответ: A = 10,67
Определённые интегралы: ∫ₐᵇ f(x) dx — это площадь со знаком между графиком и осью x от a до b. Найдите первообразную F (F′ = f) и вычислите F(b) − F(a). Степенное правило: ∫xⁿ dx = xⁿ⁺¹/(n + 1). Площадь между двумя кривыми — ∫(верхняя − нижняя) dx.
Нули
4 − x² = 0 → x = ±2
Интегрируйте
A = [4x − x³/3]−22 = 5,333 − (−5,333) = 10,67
Ответ
A = 10,67
9.Математический анализ функций нескольких переменных
Ответ: |∇f| = 3,866
Градиент: ∇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 = 3 cos(3x − 3y) → fx(2, −1) = −2,733
fy = −3 cos(3x − 3y) → fy(2, −1) = 2,733
Длина градиента
|∇f| = √((−2,733)² + 2,733²) = 3,866
Ответ
|∇f| = 3,866
10.Математический анализ функций нескольких переменных
Ответ: I = 0,25
Критические точки и двойные интегралы: в критической точке 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).
Область
0 ≤ y ≤ x, 0 ≤ x ≤ 1
Изменённый порядок
I = ∫01 ∫0x 1/(1 + x²)² dy dx = ∫01 x·1/(1 + x²)² dx
Замена u = x²
I = [−1/(2(1 + x²))]01 = 0,25
Ответ
I = 0,25
11.Последовательности и ряды
Ответ: D: Σ (7/4)ᵏ
Сходимость: если члены не стремятся к 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: Σ 3/k² → p-ряд с p = 2 > 1: сходится.
B: Σ 1/(3ᵏ + k) → Меньше (1/3)ᵏ, сходящегося геометрического ряда: сходится по признаку сравнения.
D: Σ (7/4)ᵏ → Геометрический ряд с r = 7/4 > 1: члены растут, поэтому ряд расходится.
Ответ
D: Σ (7/4)ᵏ
Ответ
D: Σ (7/4)ᵏ
12.Последовательности и ряды
Ответ: S = 7,389
Многочлены Тейлора: вблизи 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.
Стандартные ряды
Σ xᵏ/k! = eˣ
Подставьте
x = 2
S = e2 = 7,389
Ответ
S = 7,389
13.Преобразования Лапласа
Ответ: A = 25
Преобразование Лапласа: 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² + 3s + 2 = (s + 1)(s + 2)
F(s) = A/(s + 1) + B/(s + 2)
Правило закрытия
Умножьте на (s + 1) и положите s = −1: слагаемое с B обращается в нуль.
A = (2s + 27)/(s + 2) |s = −1 = 25/1 = 25
Ответ
A = 25
(B = −23, поэтому f(t) = 25·e−t − 23·e−2t)
14.Преобразования Лапласа
Ответ: g(0.25) = 1,094
Преобразование Лапласа: 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).
Передаточная функция
G(s) = Y(s)/U(s) = (2s + 7)/(s² + 6s + 8)
Простые дроби
s² + 6s + 8 = (s + 2)(s + 4)
A = 1,5, B = 0,5
Импульсная характеристика
g(t) = 1,5·e−2t + 0,5·e−4t
g(0,25) = 1,094
Ответ
g(0,25) = 1,094
15.Свёртка и линейные стационарные системы
Ответ: y(0.2) = 0,1648
Непрерывная свёртка: 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).
Свёртка со ступенькой
y(t) = 1·∫₀ᵗ e−2τ dτ = (1/2)·(1 − e−2t)
Вычислите
y(0,2) = 0,5·(1 − e−0,4) = 0,1648
Ответ
y(0,2) = 0,1648
16.Свёртка и линейные стационарные системы
Ответ: c5 = 0,09727
Импульсная характеристика и ряд Фурье: импульсная характеристика 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)|.