В окне печати выберите «Сохранить как 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″ + 2y′ + y = 0 при y(0) = 1 и y′(0) = 4. Найдите y(1,5).
y″ + 2y′ + y = 0y(0) = 1y′(0) = 4
y(1.5) =
2.Дифференциальные уравненияСложная
Решите y″ + 3y′ + 6,25y = 0 при y(0) = 3 и y′(0) = 0 (отпущено из состояния покоя). Найдите y(1,5); углы в радианах.
y″ + 3y′ + 6,25y = 0y(0) = 3y′(0) = 0
y(1.5) =
3.Матрицы и определителиСредняя
Найдите определитель det A показанной матрицы 3×3 A.
A = [4 −5 1; 1 −3 −2; −2 3 2]
det A =
4.Матрицы и определителиСложная
Найдите большее собственное значение λ показанной матрицы A.
A = [3 4; −7 −8]
λ =
5.Векторы: скалярное и векторное произведенияСредняя
Постоянная сила F = (−25, −10, 10) N действует на тело, которое движется по прямой из точки A(5, 0, 4) в точку B(−3, −1, 4). Координаты даны в метрах. Найдите работу W этой силы.
F = (−25, −10, 10) NA(5, 0, 4) mB(−3, −1, 4) m
W = J
6.Векторы: скалярное и векторное произведенияСложная
Сила F = (50, −50, 60) N действует в точке A(0, 0, −4). Координаты даны в метрах. Найдите модуль |M| момента силы относительно точка O(0, −1, 1).
F = (50, −50, 60) NA(0, 0, −4) mO(0, −1, 1) m
|M| = N·m
7.ИнтегралыСредняя
Вычислите I = ∫₀2 x·cos x dx (x в радианах).
∫₀2 x·cos x dx
I =
8.ИнтегралыСложная
Вычислите I = ∫₀2,5 x²·eˣ dx.
∫₀2,5 x²eˣ dx
I =
9.Математический анализ функций нескольких переменныхСредняя
z = −2x³e2y, где x = −t + 2 и y = −t² − 1. Используйте цепное правило, чтобы найти dz/dt при t = 0.
z = −2x³e2yx = −t + 2y = −t² − 1t = 0
dz/dt =
10.Математический анализ функций нескольких переменныхСложная
Функция f(x, y) = 2x² + xy + 2y² − 9x − 6y + 1 имеет критическую точку (2, 1). Вычислите там дискриминант D = fxxfyy − fxy².
f(x, y) = 2x² + xy + 2y² − 9x − 6y + 1P = (2, 1)
D =
11.Последовательности и рядыСредняя
Найдите радиус сходимости R степенного ряда Σk=1∞ k²·2ᵏ·xᵏ.
Σ k²·2ᵏ·xᵏx₀ = 0
R =
12.Последовательности и рядыСложная
Узнайте ряд и найдите его сумму: S = Σk=1∞ (−1)k+1·0,25ᵏ/k.
ak = (−1)k+1·0,25ᵏ/k
S =
13.Преобразования ЛапласаСредняя
F(s) = 12/(s² + 9s + 8). Её оригинал имеет вид f(t) = A·e−t + B·e−8t. Найдите коэффициент A разложением на простые дроби.
F(s) = 12/(s² + 9s + 8)
A =
14.Преобразования ЛапласаСложная
Система подчиняется уравнению y″ + 6y′ + 8y = 32·u(t), где y(0) = 4 и y′(0) = 0. В момент t = 0 подаётся ступенчатый вход u(t) = 1. С помощью преобразования Лапласа найдите y(0,1).
Импульс x(t) = 3 при 0 ≤ t < 1 s (в остальное время нуль) подан на систему с импульсной характеристикой h(t) = 1 при 0 ≤ t < 3 s (в остальное время нуль). Найдите выход y(1) при t = 1 s.
x(t) = 3, 0 ≤ t < 1 sh(t) = 1, 0 ≤ t < 3 st = 1 s
y(1) =
16.Свёртка и линейные стационарные системыСложная
Система подчиняется уравнению y′(t) + 5·y(t) = 10·x(t). Найдите её импульсную характеристику при t = 1 s, h(1).
y′ + 5y = 10xt = 1 s
h(1) =
3Растворы
1.Дифференциальные уравнения
Ответ: y(1.5) = 1,897
Линейные уравнения второго порядка: для 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)² = 0 → r = −1
Константы
A = y(0) = 1
y′(0) = B − 1A → B = 4 + 1 = 5
Вычислите
y(1,5) = (1 + 5·1,5)·e−1,5 = 1,897
Ответ
y(1,5) = 1,897
2.Дифференциальные уравнения
Ответ: y(1.5) = −0,2796
Линейные уравнения второго порядка: для 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,5 ± 2j
Константы
A = 3, B = 1,5·3/2 = 2,25
Вычислите
y(1,5) = e−2,25(3·cos 3 + 2,25·sin 3) = −0,2796
Ответ
y(1,5) = −0,2796
3.Матрицы и определители
Ответ: det A = −13
Матрицы: определитель матрицы 2×2 равен ad − bc. Для 3×3 разложите по строке: каждый элемент умножается на определитель матрицы 2×2, оставшейся после вычёркивания его строки и столбца, с чередующимися знаками + − +. Собственные значения λ удовлетворяют det(A − λI) = 0; для 2×2 это λ² − (след)·λ + det = 0.
Миноры 2×2 первой строки
M₁₁ = 0, M₁₂ = −2, M₁₃ = −3
Разложение по строке
det A = 4·0 − (−5)·(−2) + 1·(−3) = −13
Ответ
det A = −13
4.Матрицы и определители
Ответ: λ = −1
Матрицы: определитель матрицы 2×2 равен ad − bc. Для 3×3 разложите по строке: каждый элемент умножается на определитель матрицы 2×2, оставшейся после вычёркивания его строки и столбца, с чередующимися знаками + − +. Собственные значения λ удовлетворяют det(A − λI) = 0; для 2×2 это λ² − (след)·λ + det = 0.
Характеристическое уравнение
λ² − (−5)λ + (4) = 0
Решите
λ = (−5 ± √(9))/2 → λ = −4, −1
λ = −1
Ответ
λ = −1
5.Векторы: скалярное и векторное произведения
Ответ: W = 210 J
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.
Перемещение
d = B − A = (−3, −1, 4) − (5, 0, 4) = (−8, −1, 0) m
Работа
W = F·d = (−25)·(−8) + (−10)·(−1) + 10·0 = 210 J
Ответ
W = 210 J
6.Векторы: скалярное и векторное произведения
Ответ: |M| = 318 N·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.
Радиус-вектор из O в A
r = A − O = (0, 0, −4) − (0, −1, 1) = (0, 1, −5) m
Интегрирование по частям: ∫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
Пределы
(2·sin 2 + cos 2) − (0 + 1) = 0,4024
Ответ
I = 0,4024
8.Интегралы
Ответ: I = 37,59
Интегрирование по частям: ∫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² требует двух шагов.
Первый шаг
∫x²eˣ dx = x²eˣ − ∫2x eˣ dx
Второй шаг
∫2x eˣ dx = 2x eˣ − 2eˣ
∫x²eˣ dx = eˣ(x² − 2x + 2)
Пределы
12,18·(3,25) − 1·2 = 37,59
Ответ
I = 37,59
9.Математический анализ функций нескольких переменных
Ответ: dz/dt = 3,248
Градиент: ∇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.
Точка и скорости
x(0) = 2, y(0) = −1
dx/dt = −1, dy/dt = −2t = 0
Частные производные
fx = −6x²e2y = −3,248
fy = −4x³e2y = −4,331
Цепное правило
dz/dt = (−3,248)·(−1) + (−4,331)·0 = 3,248
Ответ
dz/dt = 3,248
10.Математический анализ функций нескольких переменных
Ответ: D = 15
Критические точки и двойные интегралы: в критической точке 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).
Вторые производные
fxx = 4, fyy = 4, fxy = 1
Дискриминант
D = 4·4 − 1² = 15
→ локальный минимум
Ответ
D = 15
11.Последовательности и ряды
Ответ: R = 0,5
Сходимость: если члены не стремятся к 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, сходится.
Коэффициенты
ck/ck+1 = (1/2)·(k/(k + 1))² → 1/2
Радиус
R = 0,5
Ответ
R = 0,5
В двух конечных точках признак Даламбера ничего не говорит; их нужно проверять отдельно.
12.Последовательности и ряды
Ответ: S = 0,2231
Многочлены Тейлора: вблизи 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.
Стандартные ряды
Σ (−1)k+1xᵏ/k = ln(1 + x)
Подставьте
x = 0,25
S = ln 1,25 = 0,2231
Ответ
S = 0,2231
13.Преобразования Лапласа
Ответ: A = 1,714
Преобразование Лапласа: 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² + 9s + 8 = (s + 1)(s + 8)
F(s) = A/(s + 1) + B/(s + 8)
Правило закрытия
Умножьте на (s + 1) и положите s = −1: слагаемое с B обращается в нуль.
A = 12/(s + 8) |s = −1 = 12/7 = 1,714
Ответ
A = 1,714
(B = −1,714, поэтому f(t) = 1,714·e−t − 1,714·e−8t)
14.Преобразования Лапласа
Ответ: y(0.1) = 4
Преобразование Лапласа: 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² + 6s + 8)·Y(s) = 4s + 24 + 32/s
Y(s) = (4s² + 24s + 32)/(s(s + 2)(s + 4))
Простые дроби
A = 32/(2·4) = 4
B = 0/((−2)·2) = 0
C = 0/((−4)·(−2)) = 0
Результат
y(t) = 4
y(0,1) = 4
Ответ
y(0,1) = 4
Проверка: y(0) = A + B + C = 4.
15.Свёртка и линейные стационарные системы
Ответ: y(1) = 3
Непрерывная свёртка: 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 − τ) — это импульс h, отражённый и сдвинутый так, что он заканчивается в точке τ = t.
τ ∈ [0, 1] → длина 1 s
Интеграл
y(1) = 3·1·1 = 3
Ответ
y(1) = 3
y(t) имеет форму трапеции: она растёт от 0, держится на уровне 3 и снова падает до 0 при t = 4 s.
16.Свёртка и линейные стационарные системы
Ответ: h(1) = 0,06738
Импульсная характеристика и ряд Фурье: импульсная характеристика 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)|.