Русская Википедия:Интеграл Дирихле

Материал из Онлайн справочника
Перейти к навигацииПерейти к поиску

Файл:Dirichlet 3.jpeg
Петер Густав Лежён Дирихле

В математике существует несколько интегралов, известных как интеграл Дирихле, названные в честь немецкого математика Петера Густава Лежена Дирихле, один из которых является несобственным интегралом функции sinc по положительной действительной прямой:

<math>\int_0^\infty \frac{\sin x}{x} \,dx = \frac{\pi}{2}.</math>

Этот интеграл не является абсолютно сходящимся, что означает <math>\Biggl| \frac{\sin x}{x} \Biggl|</math> не интегрируется по Лебегу, и, соответственно, интеграл Дирихле не определен в соответствии с интегрированием Лебега. Однако он определяется в соответствии с несобственным интегралом Римана или обобщенного интеграла Римана или Хенстока — Курцвейла.[1][2] Значение интеграла (в соответствии с интегралом Римана или Хенстока) может быть получено различными способами, включая через преобразование Лапласа, двойное интегрирование, дифференцирование под знаком интеграла, контурное интегрирование и ядро Дирихле.

Определение

Преобразование Лапласа

Пусть <math>f(t)</math> функция, определенная всякий раз, когда <math>t \geq 0</math>. Тогда преобразование Лапласа функции имеет вид

<math>\mathcal{L} \{f(t)\} = F(s) = \int_{0}^{\infty} e^{-st} f(t) \,dt,</math>

если интеграл существует.[3]

Свойство преобразования Лапласа, полезное для вычисления несобственных интегралов:

<math> \mathcal{L} \Biggl [ \frac{f(t)}{t} \Biggl] = \int_{s}^{\infty} F(u) \,du,

</math>

при условии, что <math>\lim_{t \rightarrow 0} \frac{f(t)}{t}</math> существует.

Это свойство можно использовать для вычисления интеграла Дирихле следующим образом:

<math>

\begin{align} \int_{0}^{\infty} \frac{\sin t}{t} \, dt &= \lim_{s \rightarrow 0} \int_{0}^{\infty} e^{-st}\frac{\sin t}{t} \, dt =\lim_{s \rightarrow 0}\mathcal{L} \Biggl [ \frac{\sin t}{t} \Biggl] \\[6pt] &= \lim_{s \rightarrow 0}\int_{s}^{\infty} \frac{du}{u^2 + 1} = \lim_{s \rightarrow 0} \arctan u \Biggl |_{s}^{\infty} \\[6pt] &= \lim_{s \rightarrow 0} \Biggl[ \frac{\pi}{2} - \arctan (s)\Biggl] = \frac{\pi}{2}, \end{align} </math>

так как <math>\mathcal{L}\{\sin t\} = \frac{1}{s^2 + 1}</math> преобразование Лапласа функции <math>\sin t</math>. (См. дифференцирование в разделе «Дифференцирование под знаком интеграла».)

Двойная интеграция

Вычисление интеграла Дирихле с помощью преобразования Лапласа эквивалентно попытке вычислить один и тот же дважды определенный интеграл двумя разными способами, изменив порядок интегрирования на противоположный, а именно:

<math>

\left( I_1 = \int_0^\infty \int _0^\infty e^{-st} \sin t \,dt \,ds \right) = \left( I_2 = \int_0^\infty \int _0^\infty e^{-st} \sin t \,ds \,dt \right),</math>

<math>\left( I_1 = \int_0^\infty \frac{1}{s^2 + 1} \,ds = \frac{\pi}{2} \right) = \left( I_2 = \int_0^\infty \frac{\sin t}{t} \,dt \right), </math> при условии <math> s > 0.

</math>

Дифференцирование под знаком интеграла (трюк Фейнмана)

Сначала перепишем интеграл как функцию дополнительной переменной <math>a</math>. Пусть

<math>f(a)=\int_0^\infty e^{-a\omega} \frac{\sin \omega} \omega \, d\omega.</math>

Чтобы вычислить интеграл Дирихле, нам необходимо определить <math>f(0)</math>.

Продифференцируем по <math>a</math> и применим формулу Лейбница для дифференцирования под знаком интеграла чтобы получить

<math>

\begin{align} \frac{df}{da} & = \frac{d}{da}\int_0^\infty e^{-a\omega} \frac{\sin \omega}{\omega} \, d\omega = \int_0^\infty \frac{\partial}{\partial a}e^{-a\omega}\frac{\sin \omega} \omega \, d\omega \\[6pt] & = -\int_0^\infty e^{-a\omega} \sin \omega \, d\omega. \end{align} </math>

Теперь, используя формулу Эйлера <math>e^{i\omega}=\cos \omega + i\sin \omega,</math> можно выразить синусоиду через комплексные экспоненциальные функции. Таким образом, мы имеем

<math>

\sin(\omega) = \frac{1}{2i} \left( e^{i \omega} - e^{-i \omega}\right). </math>

Следовательно,

<math>

\begin{align} \frac{df}{da} & = -\int_0^\infty e^{-a\omega} \sin \omega \, d\omega = -\int_{0}^{\infty} e^{-a\omega} \frac{e^{i\omega} - e^{-i\omega}}{2i} d\omega \\[6pt] &= -\frac{1}{2i} \int_{0}^{\infty} \left[ e^{-\omega(a-i)} - e^{-\omega(a + i)} \right] d\omega \\[6pt] &= -\frac{1}{2i} \left [ \frac{-1}{a - i} e^{-\omega (a - i)} - \frac{-1}{a + i} e^{-\omega (a + i)}\right] \Biggl|_0^{\infty} \\[6pt] &= -\frac{1}{2i} \left[ 0 - \left( \frac{-1}{a - i} + \frac{1}{a + i} \right) \right] = -\frac{1}{2i} \left( \frac{1}{a - i} - \frac{1}{a + i} \right) \\[6pt] &= -\frac{1}{2i} \left( \frac{a + i - (a -i)}{a^2 + 1} \right) = -\frac{1}{a^2 + 1}. \end{align} </math>

Интегрируя по <math>a</math> дает

<math>f(a) = \int \frac{-da}{a^2 + 1} = A - \arctan a,</math>

Где <math>A</math> постоянная интегрирования, которую необходимо определить. Поскольку <math>\lim_{a \to \infty} f(a) = 0,</math> <math>A = \lim_{a \to \infty} \arctan (a) = \frac{\pi}{2},</math> используя главное значение. Это означает

<math>f(a) = \frac{\pi}{2} - \arctan{a}.</math>

Наконец, для <math>a = 0</math> у нас есть <math>f(0) = \frac{\pi}{2} - \arctan(0) = \frac{\pi}{2}</math>, как прежде.

Комплексное интегрирование

Тот же результат может быть получен путем комплексного интегрирования. Рассмотрим

<math>f(z) = \frac{e^{iz}} z.</math>

Как функция комплексной переменной <math>z</math> оно имеет простой полюс в начале координат, что препятствует применению леммы Жордана, остальные условия которой выполнены.

Определим новую функцию[4]

<math>g(z) = \frac{e^{iz}}{z + i\varepsilon}.</math>

Полюс был перемещен от реальной оси, поэтому <math>g(z)</math> интегрируется по полукругу радиуса <math>R</math> в центре <math>z = 0</math> и замкнута по реальной оси. Затем берем предел <math>\varepsilon \rightarrow 0</math>.

Комплексный интеграл равен нулю по теореме о вычетах, так как внутри пути интегрирования нет полюсов.

<math>0 = \int_\gamma g(z) \,dz = \int_{-R}^R \frac{e^{ix}}{x + i\varepsilon} \, dx + \int_0^\pi \frac{e^{i(Re^{i\theta} + \theta)}}{Re^{i\theta} + i\varepsilon} iR \, d\theta.</math>

Второй член исчезает, когда <math>R</math> стремится к бесконечности. Что касается первого интеграла, то можно использовать одну версию теоремы Сохоцкого — Племеля для интегралов по вещественной прямой: для комплексной функции Шаблон:Mvar, определенной и непрерывно дифференцируемой на вещественной прямой и вещественных константах <math>a</math> и <math>b</math>, зная <math>a < 0 < b</math> можно найти

<math>\lim_{\varepsilon \to 0^+} \int_a^b \frac{f(x)}{x \pm i \varepsilon}\,dx = \mp i \pi f(0) + \mathcal{P} \int_a^b \frac{f(x)}{x} \,dx,</math>

где <math>\mathcal{P}</math> обозначает Шаблон:Нп1. Возвращаясь к приведенному выше исходному расчету, можно написать

<math>0 = \mathcal{P} \int \frac{e^{ix}}{x} \, dx - \pi i.</math>

Взяв мнимую часть с обеих сторон и отметив, что функция <math>\sin(x)/x</math> четная, мы получаем

<math>\int_{-\infty}^{+\infty} \frac{\sin(x)}{x} \,dx = 2 \int_0^{+\infty} \frac{\sin(x)}{x} \,dx.</math>

В заключение,

<math>\lim_{\varepsilon \to 0} \int_\varepsilon^\infty \frac{\sin(x)}{x} \, dx = \int_0^\infty \frac{\sin(x)}{x} \, dx = \frac \pi 2.</math>

В качестве альтернативы можно выбрать в качестве контура интегрирования для <math>f</math> объединение верхних полуплоских полуокружностей радиусов <math>\varepsilon</math> и <math>R</math> вместе с двумя соединяющими их отрезками реальной линии. С одной стороны, контурный интеграл равен нулю независимо от <math>\varepsilon</math> и <math>R</math>; с другой стороны, при <math>\varepsilon \to 0</math> и <math>R \to \infty</math> мнимая часть интеграла сходится к <math>2 I + \Im\big(\ln 0 - \ln(\pi i)\big) = 2I - \pi</math> (<math>\ln z</math> — любая ветвь логарифма на верхней полуплоскости), приводящий к <math>I = \frac{\pi}{2}</math>.

Ядро Дирихле

Пусть

<math>

D_n(x) = 1 + 2\sum_{k=1}^n \cos(2kx) = \frac{\sin[(2n+1)x]}{\sin(x)} </math>

будет ядром Дирихле.[5]

Отсюда следует, что

<math> \int_0^{\frac{\pi}{2}}D_n(x) dx = \frac{\pi}{2}. </math>

Определяем

<math>f(x) =

\begin{cases} \frac{1}{x} - \frac{1}{\sin(x)} & x \neq 0 \\[6pt] 0 & x = 0 \end{cases} </math>

Ясно, что <math>f</math> является непрерывной, когда <math> x \neq 0 </math>, чтобы увидеть её непрерывность при 0, применяется правило Лопиталя.

<math>

\lim_{x\to 0} \frac{\sin(x) - x}{x\sin(x)} = \lim_{x\to 0} \frac{\cos(x) - 1}{\sin(x) + x\cos(x)} = \lim_{x\to 0} \frac{-\sin(x)}{2\cos(x) - x\sin(x)} = 0. </math>

Следовательно, <math>f</math> удовлетворяет требованиям Шаблон:Нп1. Это означает

<math>

\lim_{\lambda \to \infty} \int_a^b f(x)\sin(\lambda x)dx = 0 \Rightarrow \lim_{\lambda \to \infty} \int_a^b \frac{\sin(\lambda x)}{x}dx =

   \lim_{\lambda \to \infty} \int_a^b \frac{\sin(\lambda x)}{\sin(x)}dx.

</math>

(Используемая здесь форма леммы Римана-Лебега доказана в цитируемой статье.)

Выбираем пределы <math>a = 0</math> и <math>b = \pi/2</math>. Мы хотим сказать что

<math>

\begin{align} \int_0^\infty \frac{\sin(t)}{t}dt = & \lim_{\lambda \to \infty} \int_0^{\lambda\frac{\pi}{2}} \frac{\sin(t)}{t}dt \\[6pt] = & \lim_{\lambda \to \infty} \int_0^{\frac{\pi}{2}} \frac{\sin(\lambda x)}{x}dx \\[6pt] = & \lim_{\lambda \to \infty} \int_0^{\frac{\pi}{2}} \frac{\sin(\lambda x)}{\sin(x)}dx \\[6pt] = & \lim_{n\to \infty} \int_0^{\frac{\pi}{2}} \frac{\sin((2n+1)x)}{\sin(x)}dx \\[6pt] = & \lim_{n\to \infty} \int_0^{\frac{\pi}{2}} D_n(x) dx = \frac{\pi}{2} \end{align} </math>

Однако для этого мы должны обосновать переключение реального предела в <math>\lambda</math> на интегральный предел в <math>n</math>. На самом деле это оправдано, если мы можем показать, что предел действительно существует. Докажем это.

Используя интегрирование по частям, мы имеем:

<math>

\int_a^b \frac{\sin(x)}{x}dx = \int_a^b \frac{d(1-\cos(x))}{x}dx = \left. \frac{1-\cos(x)}{x}\right|_a^b + \int_a^b \frac{1-\cos(x)}{x^2}dx </math>

Теперь, так как <math>a \to 0</math> и <math> b \to \infty</math>, член слева сходится без проблем. Смотри cписок пределов тригонометрических функций. Теперь покажем, что <math> \int_{-\infty}^{\infty}\frac{1-\cos(x)}{x^2}dx </math> интегрируем, что означает, что предел существует.[6]

Сначала мы стремимся оценить интеграл вблизи начала координат. Используя разложение косинуса около нуля в ряд Тейлора,

<math>

1 - \cos(x) = 1 - \sum_{k\geq 0}\frac{x^{2k}}{2k!} = -\sum_{k\geq 1}\frac{x^{2k}}{2k!}. </math>

Следовательно,

<math>

\left|\frac{1 - \cos(x)}{x^2}\right|

 =  \left|-\sum_{k\geq 0}\frac{x^{2k}}{2(k+1)!}\right|

\leq \sum_{k\geq 0} \frac{|x|^k}{k!}

 =  e^{|x|}.

</math>

Разбив интеграл на части, получим

<math>
    \int_{-\infty}^{\infty}\left|\frac{1-\cos(x)}{x^2}\right|dx

\leq \int_{-\infty}^{-\varepsilon} \frac{2}{x^2}dx +

    \int_{-\varepsilon}^{\varepsilon} e^{|x|}dx + 
    \int_{\varepsilon}^{\infty} \frac{2}{x^2}dx

\leq K, </math>

для некоторой константы <math>K > 0</math>. Это показывает, что интеграл абсолютно интегрируем, что означает, что исходный интеграл существует, и переход от <math>\lambda</math> к <math>n</math> был фактически оправдан, и доказательство завершено.

См. также

Примечания

Шаблон:Примечания

Шаблон:Интегральное исчисление

  1. Шаблон:Cite journal
  2. Шаблон:Cite book
  3. Шаблон:Cite book
  4. Appel, Walter. Mathematics for Physics and Physicists. Princeton University Press, 2007, p. 226. Шаблон:ISBN.
  5. Шаблон:Cite report
  6. Шаблон:Cite report