Теория вычетов

Интеграл $\int_{-\infty}^{+\infty} \frac{\sin x}{x}\, dx$ не берётся стандартными методами. Через теорему вычетов - две строки: расширить до комплексной плоскости, замкнуть контур, подобрать полюс. Ответ $\pi$. Тот же трюк стоит за обратным преобразованием Лапласа в каждом MATLAB, за расчётом устойчивости IIR-фильтров в scipy.signal, за фейнмановскими диаграммами в квантовой электродинамике. Вычет в полюсе - это "вес" особенности. Теорема Коши превращает сложный контурный интеграл в конечную сумму.

• **Обратное преобразование Лапласа:** $f(t) = \sum_k \mathrm{Res}(F(s)e^{st}, s_k)$ - переходные процессы в управлении и электронике
• **Устойчивость IIR-фильтров:** полюсы Z-преобразования внутри единичного круга - критерий для audio DSP, scipy.signal, MATLAB
• **Разложение на простые дроби:** коэффициенты = вычеты - основа проектирования фильтров (Баттерворт, Чебышёв)
• **Несобственные интегралы:** $\int_{-\infty}^{+\infty} f(x)dx$ через контур Бромвича - физика, финансовая математика

Предварительные знания

• Ряды Тейлора и Лорана

Полюсы и вычеты

**1825 год. Огюстен-Луи Коши.** Париж после Наполеона, Академия наук, кофе и страницы формул. Коши замечает: интегралы по замкнутым контурам чувствительны только к точкам, где функция "взрывается". Всё остальное - тихий фон. Через 100 лет Дирак использует это для вычисления квантовомеханических амплитуд. Через 200 - в каждом IIR-фильтре в наушниках.

**Изолированная особая точка** $z_0$ - место, где $f(z)$ перестаёт быть аналитической, но вокруг неё есть проколотый круг, где всё хорошо. Три типа:

**Классификация изолированных особенностей** (через коэффициенты ряда Лорана $\sum_{n=-\infty}^{\infty} c_n (z-z_0)^n$): - **Устранимая**: главная часть отсутствует ($c_n = 0$ при $n < 0$). Функция достраивается до аналитической. Вычет = 0. - **Полюс порядка $m$**: конечное число отрицательных степеней, $c_{-m} \neq 0$. Функция ведёт себя как $1/(z-z_0)^m$. - **Существенная особенность**: бесконечно много отрицательных степеней. Например, $e^{1/z}$ в нуле.

**Вычет** - это коэффициент $c_{-1}$ при $(z-z_0)^{-1}$ в ряде Лорана. Именно это число контролирует значение контурного интеграла. Для полюса порядка $m$ есть формула, не требующая раскладывать весь ряд:

**Формулы вычисления вычета в полюсе $z_0$:** **Простой полюс** ($m=1$): $$\mathrm{Res}(f, z_0) = \lim_{z \to z_0} (z - z_0)\, f(z)$$ Если $f = p/q$ и $q(z_0)=0$, $q'(z_0) \neq 0$, то $\mathrm{Res}(f, z_0) = p(z_0)/q'(z_0)$. **Полюс порядка $m$:** $$\mathrm{Res}(f, z_0) = \frac{1}{(m-1)!} \lim_{z \to z_0} \frac{d^{m-1}}{dz^{m-1}} \left[(z - z_0)^m f(z)\right]$$

Теорема Коши о вычетах

Сложный контурный интеграл по закрытой кривой вокруг нескольких полюсов - и всё сводится к сложению нескольких чисел. Интеграл зависит только от особых точек внутри контура, а не от самой формы пути. Это и есть ключевой удар:

**Теорема Коши о вычетах:** Если $f$ аналитична внутри и на простом замкнутом контуре $C$ (обходимом против часовой стрелки) кроме конечного числа изолированных особых точек $z_1, \ldots, z_n$ внутри $C$, то: $$\oint_C f(z)\, dz = 2\pi i \sum_{k=1}^{n} \mathrm{Res}(f, z_k)$$ **Частный случай (теорема Коши):** если особых точек нет - интеграл равен нулю.

**Интуиция:** каждый полюс - это "воронка" в потоке. Интеграл по контуру суммирует силу всех воронок внутри. Форма контура не важна - важно, что именно оказалось внутри.

**Знак имеет значение.** Если контур обходится по часовой стрелке - результат умножается на $-1$. В инженерии (Найквист, Брович) часто используют специфические контуры - всегда отслеживать направление обхода.

Вычисление вещественных интегралов

**Парадокс**: $\int_{-\infty}^{+\infty} \frac{\sin x}{x}\, dx$ - вещественный интеграл. Элементарными методами не берётся: нет первообразной в замкнутой форме. Метод вычетов решает это за две строки - уходя в комплексную плоскость и возвращаясь обратно.

Стандартная схема: расширить $\mathbb{R} \to \mathbb{C}$, выбрать контур (обычно полуокружность), применить теорему вычетов, взять предел при $R \to \infty$, выделить вещественную часть:

**Три класса вещественных интегралов через вычеты:** **1. Тригонометрические** $\int_0^{2\pi} R(\cos\theta, \sin\theta)\, d\theta$: Замена $z = e^{i\theta}$, $\cos\theta = (z+z^{-1})/2$, $\sin\theta = (z-z^{-1})/(2i)$, $d\theta = dz/(iz)$. Контур - единичная окружность. **2. Несобственные** $\int_{-\infty}^{+\infty} f(x)\, dx$: Замкнуть верхней полуокружностью (если $f \to 0$ достаточно быстро). Лемма Жордана: при $|z| = R \to \infty$ вклад дуги $\to 0$. **3. С ветвью логарифма или $x^\alpha$**: Контур-"ключ" (keyhole contour) вокруг разреза.

**Лемма Жордана** - ключ к обоснованию: если $|f(z)| \to 0$ при $|z| \to \infty$ (равномерно по направлению в верхней полуплоскости), то интеграл по дуге $\to 0$. Для $e^{iz} f(z)$ достаточно $f \to 0$ без ограничения на скорость. Именно поэтому интегралы с $\sin$ и $\cos$ берутся через $e^{iz}$.

Приложения: Лаплас, Z-преобразование, фильтры

Разложение рациональной функции на простые дроби - это нахождение вычетов. IIR-фильтр в вашем плеере - это передаточная функция $H(z)$, обратное Z-преобразование через вычеты. Контур Бромвича - это обратное преобразование Лапласа через вычеты. Везде один инструмент.

**Обратное преобразование Лапласа через вычеты (интеграл Бромвича):** $$f(t) = \frac{1}{2\pi i} \int_{c-i\infty}^{c+i\infty} F(s)\, e^{st}\, ds = \sum_k \mathrm{Res}\!\left(F(s) e^{st},\, s_k\right)$$ где $s_k$ - полюсы $F(s)$. **Устойчивость:** - **Непрерывное время** (Лаплас): устойчиво $\iff$ все $\mathrm{Re}(s_k) < 0$ - **Дискретное время** (Z): устойчиво $\iff$ все $|z_k| < 1$

**Разложение на простые дроби через вычеты**: дробная функция $F(s) = P(s)/Q(s)$ при $\deg P < \deg Q$ с простыми полюсами $s_k$ раскладывается в $\sum_k c_k/(s-s_k)$, где $c_k = \mathrm{Res}(F, s_k)$. Это основа проектирования IIR-фильтров в scipy.signal: функция `residue` делает именно это.

Ключевые идеи

• **Вычет** = коэффициент $c_{-1}$ в ряде Лорана. Для простого полюса: $\mathrm{Res}(f, z_0) = \lim_{z \to z_0}(z-z_0)f(z)$
• **Теорема Коши:** $\oint_C f\, dz = 2\pi i \sum \mathrm{Res}(f, z_k)$ - контурный интеграл через особые точки внутри
• **Вещественные интегралы:** расширение до $\mathbb{C}$, замкнутый контур, лемма Жордана, выделение вещественной части
• **Устойчивость:** $\mathrm{Re}(s_k) < 0$ для непрерывного времени; $|z_k| < 1$ для дискретного - оба через вычеты

Что дальше

Теория вычетов - центральный инструмент комплексного анализа:

• Вычисление интегралов вычетами — Конкретные техники: тригонометрические, несобственные, с логарифмом - детальный разбор
• Z-преобразование и DSP — Полюсы Z-преобразования = вычеты, устойчивость цифровых фильтров

Вопросы для размышления

• Почему теорема о вычетах является следствием теоремы Коши? Что меняет присутствие полюса внутри контура?
• Интеграл $\int_{-\infty}^{+\infty} e^{-x^2}dx = \sqrt{\pi}$ не берётся контурным методом - почему? Что делает функцию $e^{-x^2}$ особенной?
• IIR-фильтр оказался неустойчивым в продакшне. Как связать это с расположением полюсов Z-преобразования?

Связанные уроки

• ca-06 — Конкретные техники вычисления вещественных интегралов через контур
• dsp-03 — Полюсы Z-преобразования = вычеты - устойчивость цифровых фильтров
• ca-13 — Приложения комплексного анализа в CS - прямое развитие теории вычетов
• calc-14-improper — Несобственные интегралы вычисляются контурным методом через вычеты
• ml-13-svm — SVM через KKT - та же идея: сложное глобальное через локальные особые точки
• calc-13-techniques