Метод вычисления некоторых интегралов по траекториям в комплексной плоскости
В математической области комплексного анализа контурное интегрирование представляет собой метод оценки определенных интегралов вдоль путей в комплексной плоскости . [1] [2] [3]
Одним из применений контурных интегралов является оценка интегралов вдоль действительной прямой, которые нелегко найти, используя только методы действительной переменной. [5]
Методы контурной интеграции включают в себя:
прямое интегрирование комплекснозначной функции вдоль кривой в комплексной плоскости;
Для нахождения этих интегралов или сумм можно использовать один метод или комбинацию этих методов, или различные предельные процессы.
Кривые в комплексной плоскости
В комплексном анализе контур — это тип кривой в комплексной плоскости . В контурном интегрировании контуры дают точное определение кривых, на которых интеграл может быть соответствующим образом определен. Кривая в комплексной плоскости определяется как непрерывная функция из замкнутого интервала действительной прямой в комплексную плоскость: .
Это определение кривой совпадает с интуитивным понятием кривой, но включает параметризацию непрерывной функцией из замкнутого интервала. Это более точное определение позволяет нам рассмотреть, какими свойствами должна обладать кривая, чтобы быть полезной для интегрирования. В следующих подразделах мы сужаем набор кривых, которые мы можем интегрировать, чтобы включить только те, которые могут быть построены из конечного числа непрерывных кривых, которым можно задать направление. Более того, мы ограничим «части» от пересечения самих себя и потребуем, чтобы каждая часть имела конечную (неисчезающую) непрерывную производную. Эти требования соответствуют требованию, чтобы мы рассматривали только кривые, которые можно начертить, например, пером, в последовательности ровных, устойчивых штрихов, которые останавливаются только для начала нового фрагмента кривой, и все это без взятия пера. [6]
Направленные плавные кривые
Контуры часто определяются в терминах направленных гладких кривых. [6] Они дают точное определение «части» гладкой кривой, из которой состоит контур.
Гладкая кривая — это кривая с неисчезающей, непрерывной производной, такая, что каждая точка пересекается только один раз ( z является взаимно-однозначной), за возможным исключением кривой, у которой конечные точки совпадают ( ). В случае, когда конечные точки совпадают, кривая называется замкнутой, и функция должна быть взаимно-однозначной везде в остальном, а производная должна быть непрерывной в идентифицированной точке ( ). Гладкая кривая, которая не является замкнутой, часто называется гладкой дугой. [6]
Параметризация кривой обеспечивает естественный порядок точек на кривой: предшествует , если . Это приводит к понятию направленной гладкой кривой . Наиболее полезно рассматривать кривые независимо от конкретной параметризации. Это можно сделать, рассмотрев классы эквивалентности гладких кривых с тем же направлением. Направленная гладкая кривая затем может быть определена как упорядоченный набор точек в комплексной плоскости, который является образом некоторой гладкой кривой в их естественном порядке (согласно параметризации). Обратите внимание, что не все упорядочения точек являются естественным порядком гладкой кривой. Фактически, данная гладкая кривая имеет только два таких упорядочения. Кроме того, одна замкнутая кривая может иметь любую точку в качестве своей конечной точки, в то время как гладкая дуга имеет только два выбора для своих конечных точек.
Контуры
Контуры — это класс кривых, на которых мы определяем контурную интеграцию. Контур — это направленная кривая, которая состоит из конечной последовательности направленных гладких кривых, конечные точки которых сопоставлены для получения одного направления. Это требует, чтобы последовательность кривых была такой, чтобы конечная точка совпадала с начальной точкой для всех таких, что . Это включает в себя все направленные гладкие кривые. Кроме того, одна точка в комплексной плоскости считается контуром. Символ часто используется для обозначения соединения кривых вместе для формирования новой кривой. Таким образом, мы могли бы записать контур , который состоит из кривых, как
Контурные интегралы
Контурный интеграл сложной функции является обобщением интеграла для вещественных функций. Для непрерывных функций в комплексной плоскости контурный интеграл можно определить по аналогии с криволинейным интегралом , сначала определив интеграл по направленной гладкой кривой через интеграл по вещественному параметру. Более общее определение можно дать в терминах разбиений контура по аналогии с разбиением интервала и интегралом Римана . В обоих случаях интеграл по контуру определяется как сумма интегралов по направленным гладким кривым, составляющим контур.
Для непрерывных функций
Чтобы определить контурный интеграл таким образом, необходимо сначала рассмотреть интеграл по действительной переменной от комплексной функции. Пусть будет комплексной функцией действительной переменной, . Действительная и мнимая части часто обозначаются как и , соответственно, так что
Тогда интеграл от комплексной функции по интервалу определяется как
Теперь, чтобы определить контурный интеграл, пусть будет непрерывной функцией на направленной гладкой кривой . Пусть будет любой параметризацией , согласованной с ее порядком (направлением). Тогда интеграл вдоль обозначается
и задается как [6]
Это определение хорошо определено. То есть результат не зависит от выбранной параметризации. [6] В случае, когда действительный интеграл в правой части не существует, говорят, что интеграл вдоль не существует.
Как обобщение интеграла Римана
Обобщение интеграла Римана на функции комплексной переменной выполняется в полной аналогии с его определением для функций от действительных чисел. Разбиение направленной гладкой кривой определяется как конечный упорядоченный набор точек на . Интеграл по кривой является пределом конечных сумм значений функции, взятых в точках разбиения, в пределе, когда максимальное расстояние между любыми двумя последовательными точками разбиения (в двумерной комплексной плоскости), также называемое сеткой, стремится к нулю.
Прямые методы
Прямые методы предполагают вычисление интеграла с помощью методов, аналогичных тем, которые используются при вычислении линейных интегралов в многомерном исчислении. Это означает, что мы используем следующий метод:
параметризация контура
Контур параметризуется дифференцируемой комплекснозначной функцией действительных переменных, либо контур разбивается на части и параметризуется по отдельности.
подстановка параметризации в подынтегральное выражение
Подстановка параметризации в подынтегральное выражение преобразует интеграл в интеграл одной действительной переменной.
прямая оценка
Интеграл вычисляется методом, аналогичным методу вычисления интеграла действительной переменной.
Пример
Фундаментальным результатом комплексного анализа является то, что контурный интеграл 1/з равно 2π i , где траектория контура принимается за единичную окружность, пройденную против часовой стрелки (или любую положительно ориентированную жорданову кривую около 0). В случае единичной окружности существует прямой метод вычисления интеграла
При оценке этого интеграла используйте единичную окружность | z | = 1 в качестве контура, параметризованного z ( t ) = e it , где t ∈ [0, 2π] , тогда дз/дт = т.е. это и
что является значением интеграла. Этот результат применим только к случаю, когда z возводится в степень -1. Если степень не равна -1, то результат всегда будет равен нулю.
Приложения интегральных теорем
Приложения интегральных теорем также часто используются для вычисления контурного интеграла по контуру, что означает, что действительный интеграл вычисляется одновременно с вычислением контурного интеграла.
Контур выбирается таким образом, чтобы он следовал той части комплексной плоскости, которая описывает действительный интеграл, а также охватывал особенности подынтегральной функции, чтобы было возможно применение интегральной формулы Коши или теоремы о вычетах.
Применение этих интегральных формул дает нам значение интеграла по всему контуру.
разделение контура на контур по действительной части и мнимой части
Весь контур можно разделить на контур, который следует за частью комплексной плоскости, описывающей действительный интеграл, выбранный ранее (назовем его R ), и интеграл, который пересекает комплексную плоскость (назовем его I ). Интеграл по всему контуру равен сумме интегралов по каждому из этих контуров.
демонстрация того, что интеграл, пересекающий комплексную плоскость, не играет никакой роли в сумме
Если можно показать, что интеграл I равен нулю или если искомый действительный интеграл является несобственным, то если мы покажем, что интеграл I, как описано выше, стремится к 0 , то интеграл по R будет стремиться к интегралу по контуру R + I.
заключение
Если мы сможем продемонстрировать вышеописанный шаг, то мы сможем напрямую вычислить R , действительный интеграл.
Пример 1
Рассмотрим интеграл
Чтобы оценить этот интеграл, рассмотрим комплекснозначную функцию
которая имеет особенности в i и − i . Выберем контур, который будет охватывать действительный интеграл, здесь будет удобна полуокружность с граничным диаметром на действительной прямой (идущей, скажем, от − a до a ). Назовем этот контур C .
Поскольку единственная особенность контура — это точка i , то мы можем записать
что приводит функцию к виду для непосредственного применения формулы. Затем, используя интегральную формулу Коши,
Мы берем первую производную в приведенных выше шагах, потому что полюс является полюсом второго порядка. То есть, ( z − i ) возводится во вторую степень, поэтому мы используем первую производную от f ( z ) . Если бы это была ( z − i ), возведенная в третью степень, мы бы использовали вторую производную и делили на 2! и т. д. Случай ( z − i ) в первой степени соответствует производной нулевого порядка — просто самой f ( z ) .
Нам нужно показать, что интеграл по дуге полуокружности стремится к нулю при a → ∞ , используя оценочную лемму
где M — верхняя граница | f ( z ) | вдоль дуги, а L — длина дуги. Теперь,
Итак
Использование метода остатков
Рассмотрим ряд Лорана f ( z ) относительно i , единственной особенности , которую нам нужно рассмотреть. Тогда мы имеем
(См. пример расчета Лорана из ряда Лорана для вывода этого ряда.)
При осмотре становится ясно, что остаток — − я/4 , поэтому по теореме о вычетах имеем
Таким образом, мы получаем тот же результат, что и раньше.
Контурная заметка
В качестве отступления может возникнуть вопрос, не берем ли мы полуокружность, чтобы включить другую особенность, охватывающую − i . Чтобы интеграл вдоль действительной оси двигался в правильном направлении, контур должен двигаться по часовой стрелке, т. е. в отрицательном направлении, меняя знак интеграла в целом.
Это не влияет на использование метода остатков по сериям.
Пример 2 – Распределение Коши
Интеграл
(которая возникает в теории вероятностей как скалярное множитель характеристической функции распределения Коши ) сопротивляется методам элементарного исчисления . Мы оценим ее, выразив ее как предел контурных интегралов вдоль контура C , который идет вдоль действительной линии от − a до a , а затем против часовой стрелки вдоль полуокружности с центром в 0 от a до − a . Возьмем a больше 1, так что мнимая единица i будет заключена внутри кривой. Контурный интеграл равен
Поскольку e itz является целой функцией (не имеющей особенностей ни в одной точке комплексной плоскости), эта функция имеет особенности только там, где знаменатель z 2 + 1 равен нулю. Поскольку z 2 + 1 = ( z + i )( z − i ) , это происходит только там, где z = i или z = − i . Только одна из этих точек находится в области, ограниченной этим контуром. Вычет f ( z ) при z = i равен
Аналогичное рассуждение с дугой, которая обвивается вокруг −i , а не i, показывает, что если t < 0 , то
, наконец, мы имеем это:
(Если t = 0 , то интеграл немедленно поддается методам действительного исчисления и его значение равно π .)
Пример 3 – тригонометрические интегралы
В интегралы, включающие тригонометрические функции , можно внести определенные замены , так что интеграл преобразуется в рациональную функцию комплексной переменной, а затем для вычисления интеграла можно использовать вышеприведенные методы.
В качестве примера рассмотрим
Мы пытаемся сделать замену z = e it . Теперь вспомним
и
Принимая C за единичную окружность, подставляем и получаем:
Рассматриваемые особенности находятся при Пусть C 1 — малый круг около , а C 2 — малый круг около Тогда мы приходим к следующему:
Пример 3а – тригонометрические интегралы, общая процедура
Вышеуказанный метод может быть применен ко всем интегралам типа
где P и Q — многочлены, т.е. интегрируется рациональная функция в тригонометрических терминах. Обратите внимание, что пределы интегрирования могут быть также π и − π , как в предыдущем примере, или любой другой парой конечных точек, отстоящих друг от друга на 2 π .
Хитрость заключается в использовании подстановки z = e it , где dz = ie it dt и, следовательно,
Эта подстановка отображает интервал [0, 2π] в единичную окружность. Более того,
и
так, что рациональная функция f ( z ) по z получается из подстановки, и интеграл становится
который, в свою очередь, вычисляется путем суммирования остатков f ( z ) 1/из внутри единичной окружности.
Изображение справа иллюстрирует это, для
чего мы сейчас вычисляем. Первый шаг — признать, что
Замена дает
Полюса этой функции находятся в точках 1 ± √ 2 и −1 ± √ 2 . Из них 1 + √ 2 и −1 − √ 2 находятся вне единичной окружности (показаны красным цветом, не в масштабе), тогда как 1 − √ 2 и −1 + √ 2 находятся внутри единичной окружности (показаны синим цветом). Соответствующие остатки оба равны − я √ 2/16 , так что значение интеграла равно
Пример 4 – обрезка веток
Рассмотрим действительный интеграл
Мы можем начать с формулировки комплексного интеграла
Мы можем снова использовать интегральную формулу Коши или теорему о вычетах, чтобы получить соответствующие вычеты. Однако важно отметить, что z 1/2 = e (Log z )/2 , поэтому z 1/2 имеет ветвь среза . Это влияет на наш выбор контура C . Обычно логарифмическая ветвь среза определяется как отрицательная вещественная ось, однако это немного усложняет вычисление интеграла, поэтому мы определяем ее как положительную вещественную ось.
Затем мы используем так называемый контур замочной скважины , который состоит из небольшой окружности вокруг начала координат радиуса ε , простирающейся до отрезка прямой, параллельного и близкого к положительной действительной оси, но не касающегося ее, до почти полной окружности, возвращающейся к отрезку прямой, параллельному, близкому и расположенному ниже положительной действительной оси в отрицательном смысле, возвращаясь к малой окружности в середине.
Обратите внимание, что z = −2 и z = −4 находятся внутри большого круга. Это два оставшихся полюса, выводимых путем разложения знаменателя подынтегральной функции. Точка ветвления при z = 0 была исключена путем обхода вокруг начала координат.
Пусть γ — малый круг радиуса ε , Γ — больший, радиуса R , тогда
Можно показать, что интегралы по Γ и γ стремятся к нулю при ε → 0 и R → ∞ , по оценочному аргументу выше, который оставляет два члена. Теперь, поскольку z 1/2 = e (Log z )/2 , на контуре вне разреза ветви мы получили 2 π в аргументе вдоль γ . (По тождеству Эйлера , e i π представляет единичный вектор, который, следовательно, имеет π в качестве своего логарифма. Это π и есть то, что подразумевается под аргументом z . Коэффициент 1/2 заставляет нас использовать 2π . ) Итак
Поэтому:
Используя теорему о вычетах или интегральную формулу Коши (сначала применяя метод простейших дробей для вывода суммы двух простых контурных интегралов), получаем
Пример 5 – квадрат логарифма
В этом разделе рассматривается тип интеграла,
примером которого является.
Для вычисления этого интеграла используется функция
и ветвь логарифма, соответствующая −π < arg z ≤ π .
Мы вычислим интеграл f ( z ) вдоль контура замочной скважины, показанного справа. Как оказалось, этот интеграл является кратным исходному интегралу, который мы хотим вычислить, и по теореме Коши о вычетах мы имеем
Пусть R — радиус большого круга, а r — радиус малого. Обозначим верхнюю линию через M , а нижнюю — через N. Как и прежде, возьмем предел при R → ∞ и r → 0. Вклады от двух кругов исчезают. Например, с леммой ML имеем следующую верхнюю границу :
Чтобы вычислить вклады M и N, мы устанавливаем z = − x + iε на M и z = − x − iε на N , причем 0 < x < ∞ :
что дает
Пример 6 – логарифмы и остаток на бесконечности
Мы стремимся оценить
Это требует тщательного изучения
Мы построим f ( z ) так, чтобы она имела ветвь, срезанную на [0, 3] , показанную на диаграмме красным цветом. Для этого мы выбираем две ветви логарифма, устанавливая
и
Следовательно, разрез z 3 ⁄ 4 равен (−∞, 0] , а разрез (3 − z ) 1/4 равен (−∞, 3] . Легко видеть, что разрез произведения этих двух функций, т. е. f ( z ) , равен [0, 3] , поскольку f ( z ) на самом деле непрерывен на (−∞, 0) . Это происходит потому, что когда z = − r < 0 и мы приближаемся к разрезу сверху, f ( z ) имеет значение
При подходе снизу f ( z ) имеет значение
Но
так что мы имеем непрерывность через разрез. Это проиллюстрировано на диаграмме, где два черных ориентированных круга помечены соответствующим значением аргумента логарифма, используемого в z 3 ⁄ 4 и (3 − z ) 1/4 .
Мы будем использовать контур, показанный на диаграмме зеленым цветом. Для этого мы должны вычислить значение f ( z ) вдоль отрезков линии чуть выше и чуть ниже разреза.
Пусть z = r (в пределе, т.е. когда два зеленых круга сжимаются до радиуса ноль), где 0 ≤ r ≤ 3. Вдоль верхнего сегмента мы находим, что f ( z ) имеет значение
, а вдоль нижнего сегмента,
Отсюда следует, что интеграл ф ( з )/5 − з вдоль верхнего отрезка в пределе − iI , а вдоль нижнего отрезка I .
Если мы можем показать, что интегралы по двум зеленым кругам исчезают в пределе, то мы также имеем значение I , по теореме о вычетах Коши . Пусть радиус зеленых кругов равен ρ , где ρ < 0,001 и ρ → 0 , и применим неравенство ML . Для круга C L слева находим
Аналогично для окружности C R справа имеем
Теперь, используя теорему о вычетах Коши , мы имеем
где знак минус обусловлен направлением по часовой стрелке вокруг вычетов. Используя ветвь логарифма из предыдущего, ясно,
Полюс показан синим цветом на диаграмме. Значение упрощается до
Для остатка на бесконечности мы используем следующую формулу:
Подставляя, находим
и
где мы использовали тот факт, что −1 = e π i для второй ветви логарифма. Далее применяем биномиальное разложение, получая
Вывод таков:
Наконец, следует, что значение I равно
что дает
Оценка с теоремой о вычетах
Используя теорему о вычетах , мы можем оценить замкнутые контурные интегралы. Ниже приведены примеры оценки контурных интегралов с помощью теоремы о вычетах.
Используя теорему о вычетах, оценим этот контурный интеграл.
Напомним, что теорема о вычетах гласит:
где — остаток , а — особенности , лежащие внутри контура (при этом ни одна из них не лежит непосредственно на ).
имеет только один полюс, . Из этого мы определяем, что остаток должен быть
Соответствующий двойной контурный интеграл будет иметь следующий вид:
Теперь оценим . Между тем, составим соответствующий тройной интеграл:
Пример 2
Пусть векторное поле , и заметим, что в этом случае имеется 4 параметра. Пусть это векторное поле ограничено следующим:
Чтобы оценить это, мы должны использовать теорему о расходимости, как указано ранее, и мы должны оценить . Пусть
Таким образом, мы можем вычислить контурный интеграл с . Мы можем использовать тот же метод для вычисления контурных интегралов для любого векторного поля с .
Интегральное представление
Интегральное представление функции — это выражение функции, включающее контурный интеграл. Известны различные интегральные представления для многих специальных функций . Интегральные представления могут быть важны по теоретическим причинам, например, для аналитического продолжения или функциональных уравнений , или иногда для численных оценок .
^ Сталкер, Джон (1998). Комплексный анализ: основы классической теории функций. Springer. стр. 77. ISBN 0-8176-4038-X.
^ Бак, Джозеф; Ньюман, Дональд Дж. (1997). "Главы 11 и 12". Комплексный анализ . Springer. стр. 130–156. ISBN0-387-94756-6.
^ Кранц, Стивен Джордж (1999). "Глава 2". Справочник по комплексным переменным . Springer. ISBN0-8176-4011-8.
^ Митринович, Драгослав С.; Кечкич, Йован Д. (1984). "Глава 2". Метод вычетов Коши: теория и приложения. Springer. ISBN90-277-1623-4.
^ Митринович, Драгослав С.; Кечкич, Йован Д. (1984). "Глава 5". Метод вычетов Коши: теория и приложения. Springer. ISBN90-277-1623-4.
^ abcde Сафф, Эдвард Б.; Снайдер, Артур Дэвид (2003). "Глава 4". Основы комплексного анализа с приложениями к инженерии, науке и математике (3-е изд.). Prentice Hall. ISBN0-1390-7874-6.