Интеграл Гаусса

редактировать

Интеграл функции Гаусса, равный sqrt (π) График f (x) = e - x 2 {\ displaystyle f (x) = e ^ {- x ^ {2}}}{\ displayst yle f (x) = e ^ {- x ^ {2}}} и область между функцией и x {\ displaystyle x}x -ось, которая равна π {\ displaystyle {\ sqrt {\ pi}}}{\ sqrt {\ pi}} .

интеграл Гаусса, также известный как интеграл Эйлера – Пуассона, является интегралом от функции Гаусса f (x) = e - x 2 {\ displaystyle f (x) = e ^ {- x ^ {2}}}{\ displayst yle f (x) = e ^ {- x ^ {2}}} по всей реальной строке. Названный в честь немецкого математика Карла Фридриха Гаусса, интеграл равен

∫ - ∞ ∞ e - x 2 d x = π. {\ displaystyle \ int _ {- \ infty} ^ {\ infty} e ^ {- x ^ {2}} \, dx = {\ sqrt {\ pi}}.}{\ displaystyle \ int _ {- \ infty} ^ {\ infty} e ^ {- x ^ {2}} \, dx = {\ sqrt {\ pi}}.}

Абрахам де Муавр первоначально открыл этот тип интеграла в 1733 году, а Гаусс опубликовал точный интеграл в 1809 году. У интеграла есть широкий спектр приложений. Например, при небольшом изменении переменных он используется для вычисления нормализующей константы для нормального распределения. Один и тот же интеграл с конечными пределами тесно связан как с функцией ошибок , так и с функцией кумулятивного распределения нормального распределения . В физике этот тип интеграла часто появляется, например, в квантовой механике, чтобы найти плотность вероятности основного состояния гармонического осциллятора. Этот интеграл также используется в формулировке интеграла по путям, чтобы найти пропагатор гармонического осциллятора, и в статистической механике, чтобы найти его статистическую сумму.

, хотя нет элементарной функции существует для функции ошибок, как может быть доказано с помощью алгоритма Риша, интеграл Гаусса может быть решен аналитически с помощью методов многомерного исчисления. То есть не существует элементарного неопределенного интеграла для

∫ e - x 2 dx, {\ displaystyle \ int e ^ {- x ^ {2}} \, dx,}\ int e ^ {- x ^ {2}} \, dx,

, но определенный интеграл

∫ - ∞ ∞ e - x 2 dx {\ displaystyle \ int _ {- \ infty} ^ {\ infty} e ^ {- x ^ {2}} \, dx}\ int _ {- \ infty} ^ {\ infty} e ^ {- x ^ {2}} \, dx

можно оценить. Определенный интеграл произвольной функции Гаусса равен

∫ - ∞ ∞ e - a (x + b) 2 d x = π a. {\ displaystyle \ int _ {- \ infty} ^ {\ infty} e ^ {- a (x + b) ^ {2}} \, dx = {\ sqrt {\ frac {\ pi} {a}}}.}\ int _ {- \ infty} ^ {\ infty} e ^ {- a (x + b) ^ {2}} \, dx = {\ sqrt {\ frac {\ pi} {a}}}.

Содержание

  • 1 Вычисление
    • 1.1 По полярным координатам
      • 1.1.1 Полное доказательство
    • 1.2 По декартовым координатам
  • 2 Связь с гамма-функцией
  • 3 Обобщения
    • 3.1 Интеграл от функции Гаусса
    • 3.2 n-мерное и функциональное обобщение
    • 3.3 n-мерное с линейным членом
    • 3.4 Интегралы аналогичной формы
    • 3.5 Полиномы высшего порядка
  • 4 См. Также
  • 5 Ссылки
    • 5.1 Цитаты
    • 5.2 Источники

Вычисление

По полярным координатам

Стандартный способ вычисления интеграла Гаусса, идея которого восходит к Пуассону., заключается в использовании того свойства, что:

(∫ - ∞ ∞ e - x 2 dx) 2 = ∫ - ∞ ∞ e - x 2 dx ∫ - ∞ ∞ e - y 2 dy = ∫ - ∞ ∞ ∫ - ∞ ∞ е - (х 2 + у 2) dxdy. {\ displaystyle \ left (\ int _ {- \ infty} ^ {\ infty} e ^ {- x ^ {2}} \, dx \ right) ^ {2} = \ int _ {- \ infty} ^ { \ infty} e ^ {- x ^ {2}} \, dx \ int _ {- \ infty} ^ {\ infty} e ^ {- y ^ {2}} \, dy = \ int _ {- \ infty } ^ {\ infty} \ int _ {- \ infty} ^ {\ infty} e ^ {- (x ^ {2} + y ^ {2})} \, dx \, dy.}{\ displaystyle \ left (\ int _ {- \ infty} ^ {\ infty} e ^ {- x ^ {2}} \, dx \ right) ^ {2} = \ int _ {- \ infty} ^ {\ infty} e ^ {- x ^ {2}} \, dx \ int _ {- \ infty} ^ {\ infty} e ^ {- y ^ {2}} \, dy = \ int _ {- \ infty} ^ {\ infty} \ int _ {- \ infty} ^ {\ infty} e ^ {- (x ^ {2} + y ^ {2})} \, dx \, dy.}

Рассмотрим функция e - (x 2 + y 2) = e - r 2 {\ displaystyle e ^ {- (x ^ {2} + y ^ {2})} = e ^ {- r ^ {2}} }{\ displaystyle e ^ {- (x ^ {2} + y ^ {2})} = e ^ {- r ^ {2}}} на плоскости R 2 {\ displaystyle \ mathbb {R} ^ {2}}\ mathbb {R} ^ {2} , и вычислить его интеграл двумя способами:

  1. с одной стороны, при двойном интегрировании в декартовой системе координат его интеграл представляет собой квадрат:
    (∫ e - x 2 dx) 2; {\ displaystyle \ left (\ int e ^ {- x ^ {2}} \, dx \ right) ^ {2};}\ left (\ int e ^ { -x ^ {2}} \, dx \ right) ^ {2};
  2. , с другой стороны, посредством интеграции оболочки (случай двойного интегрирования в полярных координатах ), его интеграл вычисляется как π {\ displaystyle \ pi}\ pi

Сравнение этих двух вычислений дает интеграл, хотя следует позаботиться о задействованы несобственные интегралы.

R 2 e - (x 2 + y 2) dxdy = ∫ 0 2 π ∫ 0 ∞ e - r 2 rdrd θ = 2 π ∫ 0 ∞ re - r 2 dr = 2 π ∫ - ∞ 0 1 2 esdss = - r 2 знак равно π ∫ - ∞ 0 esds = π (e 0 - e - ∞) = π, {\ displaystyle {\ begin {align} \ iint _ {\ mathbf {R} ^ {2}} e ^ {- (x ^ {2} + y ^ {2})} dx \, dy = \ int _ {0} ^ {2 \ pi} \ int _ {0} ^ {\ infty} e ^ {- r ^ {2}} r \, dr \, d \ theta \\ [6pt] = 2 \ pi \ int _ {0} ^ {\ infty} re ^ {- r ^ {2}} \, dr \\ [ 6pt] = 2 \ pi \ int _ {- \ infty} ^ {0} {\ tfrac {1} {2}} e ^ {s} \, ds s = -r ^ {2} \\ [6pt] = \ pi \ int _ {- \ infty} ^ {0} e ^ {s} \, ds \\ [6pt] = \ pi (e ^ {0} -e ^ {- \ infty}) \\ [ 6pt] = \ pi, \ end {align}}}{\ displaystyle {\ begin {align} \ iint _ {\ mathbf {R} ^ {2}} e ^ {- (x ^ {2} + y ^ {2})} dx \, dy = \ int _ {0} ^ {2 \ pi} \ int _ {0} ^ {\ infty} e ^ {- r ^ {2}} r \, dr \, d \ theta \\ [6pt] = 2 \ pi \ int _ {0} ^ {\ infty} re ^ {- r ^ {2}} \, dr \\ [6pt] = 2 \ pi \ int _ {- \ infty} ^ {0} {\ tfrac {1} {2}} e ^ {s} \, ds s = -r ^ {2} \\ [6pt] = \ pi \ int _ {- \ infty} ^ {0} e ^ {s} \, ds \\ [6pt] = \ pi (e ^ {0} -e ^ {- \ infty}) \\ [6pt] = \ pi, \ end {align}}}

где множитель r является определителем Якоби, который появляется из-за преобразования в полярные координаты (r dr dθ - стандартная мера на плоскости, выраженная в полярных координатах Wikibooks: Calculus / Polar Integration # Generalization ), а при замене берется s = −r, поэтому ds = −2r dr.

Объединение этих результатов

(∫ - ∞ ∞ e - x 2 dx) 2 = π, {\ displaystyle \ left (\ int _ {- \ infty} ^ {\ infty} e ^ {- x ^ {2}} \, dx \ right) ^ {2} = \ pi,}\ left (\ int _ {- \ infty} ^ {\ infty} e ^ {- x ^ {2}} \, dx \ right) ^ {2} = \ pi,

, поэтому

∫ - ∞ ∞ e - x 2 dx = π {\ displaystyle \ int _ {- \ infty} ^ {\ infty} e ^ {- x ^ {2}} \, dx = {\ sqrt {\ pi}}}\ int _ {- \ infty} ^ {\ infty} e ^ {- x ^ {2}} \, dx = {\ sqrt {\ pi}} .

Полное доказательство

Чтобы оправдать неправильные двойные интегралы и приравнять два выражения, начнем с аппроксимирующей функции:

I (a) = ∫ - aae - x 2 dx. {\ displaystyle I (a) = \ int _ {- a} ^ {a} e ^ {- x ^ {2}} dx.}I (a) = \ int _ {- a} ^ {a} e ^ {- x ^ {2}} dx.

Если интеграл

∫ - ∞ ∞ e - x 2 dx {\ displaystyle \ int _ {- \ infty} ^ {\ infty} e ^ {- x ^ {2}} \, dx}\ int _ {- \ infty} ^ {\ infty} e ^ {- x ^ {2}} \, dx

были бы абсолютно сходящимися, мы бы получили это Главное значение Коши, то есть предел

lim a → ∞ I (a) {\ displaystyle \ lim _ {a \ to \ infty} I (a)}\ lim _ {a \ to \ infty} I (a)

, будет совпадать с

∫ - ∞ ∞ е - х 2 дх. {\ displaystyle \ int _ {- \ infty} ^ {\ infty} e ^ {- x ^ {2}} \, dx.}\ int _ {- \ infty} ^ {\ infty} e ^ {- x ^ {2}} \, dx.

Чтобы убедиться, что это так, примите во внимание, что

∫ - ∞ ∞ | e - x 2 | dx < ∫ − ∞ − 1 − x e − x 2 d x + ∫ − 1 1 e − x 2 d x + ∫ 1 ∞ x e − x 2 d x < ∞. {\displaystyle \int _{-\infty }^{\infty }|e^{-x^{2}}|\,dx<\int _{-\infty }^{-1}-xe^{-x^{2}}\,dx+\int _{-1}^{1}e^{-x^{2}}\,dx+\int _{1}^{\infty }xe^{-x^{2}}\,dx<\infty.}\ int _ {- \ infty} ^ {\ infty} | e ^ {- x ^ {2}} | \, dx <\ int _ {- \ infty} ^ {- 1} -xe ^ {- x ^ {2}} \, dx + \ int _ {- 1} ^ {1} e ^ {- x ^ {2}} \, dx + \ int _ {1} ^ {\ infty} xe ^ {- x ^ {2}} \, dx <\ infty.

, чтобы мы могли вычислить

∫ - ∞ ∞ e - x 2 dx {\ displaystyle \ int _ {- \ infty} ^ {\ infty} e ^ {- x ^ {2}} \, dx}\ int _ {- \ infty} ^ {\ infty} e ^ {- x ^ {2}} \, dx

просто взяв предел

lim a → ∞ I (a) {\ displaystyle \ lim _ {a \ to \ infty} I (a)}\ lim _ {a \ to \ infty} I (a) .

Взяв квадрат I (a) {\ displaystyle I (a)}I (a) дает

I 2 (a) = (∫ - aae - x 2 dx) (∫ - aae - y 2 dy) = ∫ - aa (∫ - aae - y 2 dy) e - x 2 dx = ∫ - aa ∫ - aae - (x 2 + y 2) dydx. {\ displaystyle {\ begin {align} I ^ {2} (a) = \ left (\ int _ {- a} ^ {a} e ^ {- x ^ {2}} \, dx \ right) \ left (\ int _ {- a} ^ {a} e ^ {- y ^ {2}} \, dy \ right) \\ [6pt] = \ int _ {- a} ^ {a} \ left ( \ int _ {- a} ^ {a} e ^ {- y ^ {2}} \, dy \ right) \, e ^ {- x ^ {2}} \, dx \\ [6pt] = \ int _ {- a} ^ {a} \ int _ {- a} ^ {a} e ^ {- (x ^ {2} + y ^ {2})} \, dy \, dx. \ end {выровнено }}}{\ displaystyle {\ begin { выровнено} I ^ {2} (a) = \ left (\ int _ {- a} ^ {a} e ^ {- x ^ {2}} \, dx \ right) \ left (\ int _ {- a} ^ {a} e ^ {- y ^ {2}} \, dy \ right) \\ [6pt] = \ int _ {- a} ^ {a} \ left (\ int _ {- a} ^ {a} e ^ {- y ^ {2}} \, dy \ right) \, e ^ {- x ^ {2}} \, dx \\ [6pt] = \ int _ {- a} ^ {a} \ int _ {- a} ^ {a} e ^ {- (x ^ {2} + y ^ {2})} \, dy \, dx. \ end {align}}}

Используя теорему Фубини, указанный выше двойной интеграл можно рассматривать как интеграл площадей

∬ [- a, a] × [- a, a] e - (x 2 + y 2) d (x, y), {\ displaystyle \ iint _ {[- a, a] \ times [-a, a]} e ^ {- (x ^ {2} + y ^ {2})} \, d (x, y),}\ iint _ {[- a, a] \ times [-a, a]} e ^ {- (x ^ {2} + y ^ {2})} \, d (x, y),

взято над квадратом с вершинами {(−a, a), (a, a), (a, −a), (−a, −a)} на xy- plane.

Поскольку экспоненциальная функция больше 0 для всех действительных чисел, из этого следует, что интеграл, взятый по вписанной в квадрат окружности, должен быть меньше I (a) 2 {\ displaystyle I (a) ^ {2}}I (a) ^ {2} , и аналогично интеграл, взятый по описанной окружности квадрата, должен быть больше I (a) 2 {\ displaystyle Я (а) ^ {2}}I (a) ^ {2} . Интегралы по двум дискам можно легко вычислить, переключившись с декартовых координат на полярные координаты :

x = r cos ⁡ θ y = r sin ⁡ θ {\ displaystyle {\ begin {align} x = r \ cos \ theta \\ y = r \ sin \ theta \ end {align}}}{\ displaystyle {\ begin {align} x = r \ cos \ theta \\ y = r \ sin \ theta \ end {align}}}
J (r, θ) = [∂ x ∂ r ∂ x ∂ θ ∂ y ∂ r ∂ y ∂ θ] = [cos ⁡ θ - р грех ⁡ θ грех ⁡ θ р соз ⁡ θ] {\ displaystyle \ mathbf {J} (r, \ theta) = {\ begin {bmatrix} {\ dfrac {\ partial x} {\ partial r}} {\ dfrac {\ partial x} {\ partial \ theta}} \\ [1em] {\ dfrac {\ partial y} {\ partial r}} и {\ dfrac {\ partial y} {\ partial \ theta}} \ end {bmatrix}} = {\ begin {bmatrix} \ cos \ theta -r \ sin \ theta \\\ sin \ theta r \ cos \ theta \ end {bmatrix}}}{\ displaystyle \ mathbf {J} (r, \ theta) = {\ begin {bmatrix} {\ dfrac {\ partial x} {\ partial r}} {\ dfrac {\ partial x } {\ partial \ theta}} \\ [1em] {\ dfrac {\ partial y} {\ partial r}} {\ dfrac {\ partial y} {\ partial \ theta}} \ end {bmatrix}} = {\ begin {bmatrix} \ cos \ theta -r \ грех \ theta \\\ sin \ theta r \ cos \ theta \ end {bmatrix}}}
d (x, y) = | J (r, θ) | d (r, θ) = r d (r, θ). {\ displaystyle d (x, y) = | J (r, \ theta) | d (r, \ theta) = r \, d (r, \ theta).}{\ displaystyle d (x, y) = | J (r, \ theta) | d (r, \ theta) = r \, d (r, \ theta).}
∫ 0 2 π ∫ 0 являются - r 2 drd θ < I 2 ( a) < ∫ 0 2 π ∫ 0 a 2 r e − r 2 d r d θ. {\displaystyle \int _{0}^{2\pi }\int _{0}^{a}re^{-r^{2}}\,dr\,d\theta \ int _ {0} ^ {2 \ pi} \ int _ {0} ^ {a} re ^ {- r ^ {2}} \, dr \, d \ theta <I ^ {2} (a) <\ int _ {0} ^ {2 \ pi} \ int _ {0} ^ {a {\ sqrt {2}}} re ^ {- r ^ {2}} \, dr \, d \ the ta.

(См. в полярных координатах от декартовых координат для помощи с полярным преобразованием.)

Интегрирование,

π (1 - e - a 2) < I 2 ( a) < π ( 1 − e − 2 a 2). {\displaystyle \pi (1-e^{-a^{2}})\ pi (1-e ^ {- a ^ {2}}) <I ^ {2} (a) <\ pi (1-e ^ {- 2a ^ {2}}).

По теореме сжатия это дает гауссовский интеграл

∫ - ∞ ∞ e - x 2 dx = π. {\ displaystyle \ int _ {- \ infty} ^ {\ infty} e ^ {- x ^ {2}} \, dx = {\ sqrt {\ pi}}.}\ int _ {- \ infty} ^ {\ infty} e ^ {- х ^ {2}} \, dx = {\ sqrt {\ p i}}.

По декартовым координатам

Другая техника, восходящая к Лапласу (1812 г.), заключается в следующем. Пусть

y = x s d y = x d s. {\ displaystyle {\ begin {align} y = xs \\ dy = x \, ds. \ end {align}}}{ \ begin {align} y = xs \\ dy = x \, ds. \ end {выравнивается}}

Поскольку пределы s при y → ± ∞ зависят от знака x, это упрощает при вычислении используется тот факт, что e является четной функцией , и, следовательно, интеграл по всем действительным числам равен удвоенному интегралу от нуля до бесконечности. То есть

∫ - ∞ ∞ e - x 2 d x = 2 ∫ 0 ∞ e - x 2 d x. {\ displaystyle \ int _ {- \ infty} ^ {\ infty} e ^ {- x ^ {2}} \, dx = 2 \ int _ {0} ^ {\ infty} e ^ {- x ^ {2 }} \, dx.}\ int _ {- \ infty} ^ {\ infty} e ^ {- x ^ {2}} \, dx = 2 \ int _ {0} ^ {\ infty} e ^ {- х ^ {2}} \, dx.

Таким образом, в диапазоне интегрирования x ≥ 0, а переменные y и s имеют одинаковые пределы. Это дает:

I 2 = 4 ∫ 0 ∞ ∫ 0 ∞ e - (x 2 + y 2) dydx = 4 ∫ 0 ∞ (∫ 0 ∞ e - (x 2 + y 2) dy) dx = 4 ∫ 0 ∞ (∫ 0 ∞ e - x 2 (1 + s 2) xds) dx = 4 ∫ 0 ∞ (∫ 0 ∞ e - x 2 (1 + s 2) xdx) ds = 4 ∫ 0 ∞ [1-2 (1 + s 2) e - x 2 (1 + s 2)] x = 0 x = ∞ ds = 4 (1 2 ∫ 0 ∞ ds 1 + s 2) = 2 [arctan ⁡ s] 0 ∞ = π. {\ displaystyle {\ begin {align} I ^ {2} = 4 \ int _ {0} ^ {\ infty} \ int _ {0} ^ {\ infty} e ^ {- (x ^ {2} + y ^ {2})} dy \, dx \\ [6pt] = 4 \ int _ {0} ^ {\ infty} \ left (\ int _ {0} ^ {\ infty} e ^ {- (x ^ {2} + y ^ {2})} \, dy \ right) \, dx \\ [6pt] = 4 \ int _ {0} ^ {\ infty} \ left (\ int _ {0} ^ {\ infty} e ^ {- x ^ {2} (1 + s ^ {2})} x \, ds \ right) \, dx \\ [6pt] = 4 \ int _ {0} ^ {\ infty} \ left (\ int _ {0} ^ {\ infty} e ^ {- x ^ {2} (1 + s ^ {2})} x \, dx \ right) \, ds \\ [6pt] = 4 \ int _ {0} ^ {\ infty} \ left [{\ frac {1} {- 2 (1 + s ^ {2})}} e ^ {- x ^ {2} (1 + s ^ {2})} \ right] _ {x = 0} ^ {x = \ infty} \, ds \\ [6pt] = 4 \ left ({\ frac {1} {2}} \ int _ { 0} ^ {\ infty} {\ frac {ds} {1 + s ^ {2}}} \ right) \\ [6pt] = 2 {\ Big [} \ arctan s {\ Big]} _ {0 } ^ {\ infty} \\ [6pt] = \ pi. \ end {align}}}{\ displaystyle {\ begin {align} I ^ {2} = 4 \ int _ {0} ^ {\ infty} \ int _ {0} ^ {\ infty} e ^ {- (x ^ {2} + y ^ {2})} dy \, dx \\ [6pt ] = 4 \ int _ {0} ^ {\ infty} \ left (\ int _ {0} ^ {\ infty} e ^ {- (x ^ {2} + y ^ {2})} \, dy \ right) \, dx \\ [6pt] = 4 \ int _ {0} ^ {\ infty} \ left (\ int _ {0} ^ {\ infty} e ^ {- x ^ {2} (1 + s ^ {2})} x \, ds \ right) \, dx \\ [6pt] = 4 \ int _ {0} ^ {\ infty} \ left (\ int _ {0} ^ {\ infty} e ^ {- x ^ {2} (1 + s ^ {2})} x \, dx \ right) \, ds \\ [6pt] = 4 \ int _ {0} ^ {\ infty} \ left [{\ frac {1} {- 2 (1 + s ^ {2})}} e ^ {-x ^ {2} (1 + s ^ {2})} \ right] _ {x = 0} ^ {x = \ infty} \, ds \\ [6pt] = 4 \ left ({\ frac {1} {2}} \ int _ {0} ^ {\ infty} {\ frac {ds} {1 + s ^ {2}}} \ right) \\ [6pt] = 2 {\ Big [} \ arctan s {\ Big]} _ {0} ^ {\ infty} \\ [6pt] = \ pi. \ end {align}}}

Следовательно, I = π {\ displaystyle I = {\ sqrt {\ pi}}}I = {\ sqrt {\ pi}} , как и ожидалось.

Связь с гамма-функцией

Подынтегральное выражение - это четная функция,

∫ - ∞ ∞ e - x 2 dx = 2 ∫ 0 ∞ e - x 2 dx {\ displaystyle \ int _ {- \ infty} ^ {\ infty} e ^ {- x ^ {2}} dx = 2 \ int _ {0} ^ {\ infty} e ^ {- x ^ {2}} dx}\ int _ {- \ infty} ^ {\ infty} e ^ {- x ^ {2}} dx = 2 \ int _ {0} ^ {\ infty} e ^ {- x ^ {2}} dx

Таким образом, после замены переменной x = t {\ displaystyle x = {\ sqrt {t}}}x = {\ sqrt {t}} это превращается в интеграл Эйлера

2 ∫ 0 ∞ e - Икс 2 dx знак равно 2 ∫ 0 ∞ 1 2 e - tt - 1 2 dt = Γ (1 2) = π {\ displaystyle 2 \ int _ {0} ^ {\ infty} e ^ {- x ^ {2} } dx = 2 \ int _ {0} ^ {\ infty} {\ frac {1} {2}} \ e ^ {- t} \ t ^ {- {\ frac {1} {2}}} dt = \ Gamma \ left ({\ frac {1} {2}} \ right) = {\ sqrt {\ pi}}}2 \ int _ {0} ^ {\ infty} e ^ {- x ^ {2}} dx = 2 \ int _ {0} ^ {\ infty} {\ frac {1} {2}} \ e ^ {- t} \ t ^ {- {\ frac {1 } {2}}} dt = \ Gamma \ left ({\ frac {1} {2}} \ right) = {\ sqrt {\ pi}}

где Γ (z) = ∫ 0 ∞ tz - 1 e - tdt { \ displaystyle \ Gamma (z) = \ int _ {0} ^ {\ infty} t ^ {z-1} e ^ {- t} dt}{\ displaystyle \ Gamma (г) = \ int _ {0} ^ {\ infty } t ^ {z-1} e ^ {- t} dt} - это гамма-функция. Это показывает, почему факториал полуцелого числа является рациональным кратным π {\ displaystyle {\ sqrt {\ pi}}}{\ sqrt {\ pi}} . В более общем смысле,

∫ 0 ∞ xne - axbdx = Γ ((n + 1) / b) ba (n + 1) / b, {\ displaystyle \ int _ {0} ^ {\ infty} x ^ {n } e ^ {- ax ^ {b}} dx = {\ frac {\ Gamma \ left ((n + 1) / b \ right)} {ba ^ {(n + 1) / b}}},}{\ displaystyle \ int _ {0} ^ {\ infty} x ^ {n} e ^ {- ax ^ {b}} dx = {\ frac {\ Гамма \ left ((n + 1) / b \ right)} {ba ^ {(n + 1) / b}}},}

, который можно получить, подставив t = axb {\ displaystyle t = ax ^ {b}}{\ displaystyle t = ax ^ {b}} в подынтегральное выражение гамма-функции, чтобы получить Γ (z) = azb ∫ 0 ∞ xbz - 1 e - axbdx {\ displaystyle \ Gamma (z) = a ^ {z} b \ int _ {0} ^ {\ infty} x ^ {bz-1} e ^ {- ax ^ {b} } dx}{\ displaystyle \ Gamma (z) = a ^ {z} b \ int _ {0} ^ {\ infty} x ^ {bz-1} e ^ {- ax ^ {b}} dx} .

Обобщения

Интеграл от функции Гаусса

Интеграл от произвольной функции Гаусса равен

∫ - ∞ ∞ e - a ( x + b) 2 dx = π a. {\ displaystyle \ int _ {- \ infty} ^ {\ infty} e ^ {- a (x + b) ^ {2}} \, dx = {\ sqrt {\ frac {\ pi} {a}}}.}\ int _ {- \ infty} ^ {\ infty} e ^ {- a (x + b) ^ {2}} \, dx = {\ sqrt {\ frac {\ pi} {a}}}.

Альтернативная форма:

∫ - ∞ ∞ e - ax 2 + bx + cdx = π aeb 2 4 a + c. {\ displaystyle \ int _ {- \ infty} ^ {\ infty} e ^ {- ax ^ {2} + bx + c} \, dx = {\ sqrt {\ frac {\ pi} {a}}} \, e ^ {{\ frac {b ^ {2}} {4a}} + c}.}{\ displaystyle \ int _ {- \ infty} ^ {\ infty} e ^ {- ax ^ {2} + bx + c} \, dx = {\ sqrt {\ frac {\ pi} {a}}} \, e ^ {{\ frac {b ^ {2}} {4a}} + c}.}

Эта форма полезна для вычисления ожиданий некоторых непрерывных распределений вероятностей, связанных с нормальным распределением, таких как логнормальное распределение, например.

n-мерное и функциональное обобщение

Предположим, что A - симметричная положительно определенная (следовательно, обратимая) матрица точности размера n × n , которая является матрицей, обратной матрице ковариационная матрица. Тогда

∫ - ∞ ∞ exp ⁡ (- 1 2 ∑ i, j = 1 n A ijxixj) dnx = ∫ - ∞ ∞ exp ⁡ (- 1 2 x TA x) dnx = (2 π) n det A Знак равно 1 det (A / 2 π) = det (2 π A - 1) {\ displaystyle \ int _ {- \ infty} ^ {\ infty} \ exp {\ left (- {\ frac {1} {2} } \ sum \ limits _ {i, j = 1} ^ {n} A_ {ij} x_ {i} x_ {j} \ right)} \, d ^ {n} x = \ int _ {- \ infty} ^ {\ infty} \ exp {\ left (- {\ frac {1} {2}} x ^ {T} Ax \ right)} \, d ^ {n} x = {\ sqrt {\ frac {(2 \ pi) ^ {n}} {\ det A}}} = {\ sqrt {\ frac {1} {\ det (A / 2 \ pi)}}} = {\ sqrt {\ det (2 \ pi A ^ {- 1})}}}{\ displaystyle \ int _ {- \ infty} ^ {\ infty} \ exp { \ left (- {\ frac {1} {2}} \ sum \ limits _ {i, j = 1} ^ {n} A_ {ij} x_ {i} x_ {j} \ right)} \, d ^ {n} x = \ int _ {- \ infty} ^ {\ infty} \ exp {\ left (- {\ frac {1} {2}} x ^ {T} Ax \ right)} \, d ^ { n} x = {\ sqrt {\ frac {(2 \ pi) ^ {n}} {\ det A}}} = {\ sqrt {\ frac {1} {\ det (A / 2 \ pi)}} } = {\ sqrt {\ det (2 \ pi A ^ {- 1})}}}

где интеграл считается превышающим R . Этот факт применяется при исследовании многомерного нормального распределения .

Кроме того,

∫ xk 1 ⋯ xk 2 N exp ⁡ (- 1 2 ∑ i, j = 1 n A ijxixj) dnx = (2 π) n det A 1 2 NN! ∑ σ ∈ S 2 N (A - 1) k σ (1) k σ (2) ⋯ (A - 1) k σ (2 N - 1) k σ (2 N) {\ displaystyle \ int x_ {k_ { 1}} \ cdots x_ {k_ {2N}} \, \ exp {\ left (- {\ frac {1} {2}} \ sum \ limits _ {i, j = 1} ^ {n} A_ {ij } x_ {i} x_ {j} \ right)} \, d ^ {n} x = {\ sqrt {\ frac {(2 \ pi) ^ {n}} {\ det A}}} \, {\ frac {1} {2 ^ {N} N!}} \, \ sum _ {\ sigma \ in S_ {2N}} (A ^ {- 1}) _ {k _ {\ sigma (1)} k _ {\ sigma (2)}} \ cdots (A ^ {- 1}) _ {k _ {\ sigma (2N-1)} k _ {\ sigma (2N)}}}{\ displaystyle \ int x_ {k_ {1} } \ cdots x_ {k_ {2N}} \, \ exp {\ left (- {\ frac {1} {2}} \ sum \ limits _ {i, j = 1} ^ {n} A_ {ij} x_ {i} x_ {j} \ right)} \, d ^ {n} x = {\ sqrt {\ frac {(2 \ pi) ^ {n}} {\ det A}}} \, {\ frac { 1} {2 ^ {N} N!}} \, \ Sum _ {\ sigma \ in S_ {2N}} (A ^ {- 1}) _ {k _ {\ sigma (1)} k _ {\ sigma ( 2)}} \ cdots (A ^ {- 1}) _ {k _ {\ sigma (2N-1)} k _ {\ sigma (2N)}}}

где σ - это перестановка из {1,..., 2N}, а дополнительный множитель в правой части - это сумма по всем комбинаторным парам {1,..., 2N} из N копий A.

В качестве альтернативы

∫ f (x →) exp ⁡ (- 1 2 ∑ i, j = 1 n A ijxixj) dnx = (2 π) n det A exp ⁡ (1 2 ∑ i, j = 1 n ( A - 1) ij ∂ ∂ xi ∂ ∂ xj) f (x →) | Икс → знак равно 0 {\ Displaystyle \ int е ({\ vec {x}}) \ exp {\ left (- {\ frac {1} {2}} \ sum \ limits _ {i, j = 1} ^ { n} A_ {ij} x_ {i} x_ {j} \ right)} d ^ {n} x = {\ sqrt {(2 \ pi) ^ {n} \ over \ det A}} \, \ left. \ exp {\ left ({1 \ over 2} \ sum \ limits _ {i, j = 1} ^ {n} (A ^ {- 1}) _ {ij} {\ partial \ over \ partial x_ {i }} {\ partial \ over \ partial x_ {j}} \ right)} f ({\ vec {x}}) \ right | _ {{\ vec {x}} = 0}}{\ displaystyle \ int f ({\ vec {x}}) \ exp {\ left (- {\ frac {1} {2}} \ sum \ limits _ {i, j = 1} ^ {n} A_ {ij} x_ {i} x_ {j} \ right)} d ^ {n} x = {\ sqrt {(2 \ pi) ^ {n} \ over \ det A}} \, \ left. \ exp {\ left ({1 \ over 2} \ sum \ limits _ {i, j = 1} ^ {n} (A ^ {- 1}) _ {ij} {\ partial \ over \ частичный x_ {i}} {\ partial \ over \ partial x_ {j}} \ right)} f ({\ vec {x}}) \ right | _ {{\ vec {x}} = 0}}

для некоторых аналитическая функция f, при условии, что она удовлетворяет некоторым подходящим ограничениям на ее рост и некоторым другим техническим критериям. (Он работает для некоторых функций и не работает для других. Полиномы - это хорошо.) Экспонента над дифференциальным оператором понимается как степенной ряд.

, в то время как функциональные интегралы не имеют строгого определения (или даже в большинстве случаев нестрогий вычислительный), мы можем определить гауссовский функциональный интеграл по аналогии с конечномерным случаем. Однако остается проблема, заключающаяся в том, что (2 π) ∞ {\ displaystyle (2 \ pi) ^ {\ infty}}(2 \ пи) ^ {\ infty} бесконечно, а также функциональный определитель Тоже было бы вообще бесконечно. Об этом можно позаботиться, если мы будем рассматривать только отношения:

∫ f (x 1) ⋯ f (x 2 N) exp ⁡ [- ∬ 1 2 A (x 2 N + 1, x 2 N + 2) f (x 2 N + 1) f (x 2 N + 2) ddx 2 N + 1 ddx 2 N + 2] D f ∫ exp ⁡ [- ∬ 1 2 A (x 2 N + 1, x 2 N + 2) f (x 2 N + 1) f (x 2 N + 2) ddx 2 N + 1 ddx 2 N + 2] D f = 1 2 NN! ∑ σ ∈ S 2 N A - 1 (x σ (1), x σ (2)) ⋯ A - 1 (x σ (2 N - 1), x σ (2 N)). {\ displaystyle {\ frac {\ int f (x_ {1}) \ cdots f (x_ {2N}) \ exp \ left [{- \ iint {\ frac {1} {2}} A (x_ {2N + 1}, x_ {2N + 2}) f (x_ {2N + 1}) f (x_ {2N + 2}) d ^ {d} x_ {2N + 1} d ^ {d} x_ {2N + 2} } \ right] {\ mathcal {D}} f} {\ int \ exp \ left [{- \ iint {\ frac {1} {2}} A (x_ {2N + 1}, x_ {2N + 2})) f (x_ {2N + 1}) f (x_ {2N + 2}) d ^ {d} x_ {2N + 1} d ^ {d} x_ {2N + 2}} \ right] {\ mathcal {D }} f}} = {\ frac {1} {2 ^ {N} N!}} \ sum _ {\ sigma \ in S_ {2N}} A ^ {- 1} (x _ {\ sigma (1)}, x _ {\ sigma (2)}) \ cdots A ^ {- 1} (x _ {\ sigma (2N-1)}, x _ {\ sigma (2N)}).}{\ displaystyle {\ frac {\ int f (x_ {1}) \ cdots f (x_ {2N) }) \ exp \ left [{- \ iint {\ frac {1} {2}} A (x_ {2N + 1}, x_ {2N + 2}) f (x_ {2N + 1}) f (x_ { 2N + 2}) d ^ {d} x_ {2N + 1} d ^ {d} x_ {2N + 2}} \ right] {\ mathcal {D}} f} {\ int \ exp \ left [{- \ iint {\ frac {1} {2}} A (x_ {2N + 1}, x_ {2N + 2}) f (x_ {2N + 1}) f (x_ {2N + 2}) d ^ {d } x_ {2N + 1} d ^ {d} x_ {2N + 2}} \ right] {\ mathcal {D}} f}} = {\ frac {1} {2 ^ {N} N!}} \ сумма _ {\ sigma \ in S_ {2N}} A ^ {- 1} (x _ {\ sigma (1)}, x _ {\ sigma (2)}) \ cdots A ^ {- 1} (x _ {\ sigma (2N-1)}, x _ {\ sigma (2N)}).}

В DeWitt обозначение, уравнение выглядит идентично конечномерному случаю.

n-мерная с линейным членом

Если A снова является симметричной положительно определенной матрицей, то (при условии, что все являются векторами-столбцами)

∫ e - 1 2 ∑ i, j = 1 n A ijxixj + ∑ i = 1 n B ixidnx = ∫ e - 1 2 x → TA x → + B → T x → dnx = (2 π) n det A e 1 2 B → TA - 1 B →. {\ displaystyle \ int e ^ {- {\ frac {1} {2}} \ sum \ limits _ {i, j = 1} ^ {n} A_ {ij} x_ {i} x_ {j} + \ sum \ limits _ {i = 1} ^ {n} B_ {i} x_ {i}} d ^ {n} x = \ int e ^ {- {\ frac {1} {2}} {\ vec {x} } ^ {T} \ mathbf {A} {\ vec {x}} + {\ vec {B}} ^ {T} {\ vec {x}}} d ^ {n} x = {\ sqrt {\ frac {(2 \ pi) ^ {n}} {\ det {A}}}} e ^ {{\ frac {1} {2}} {\ vec {B}} ^ {T} \ mathbf {A} ^ {-1} {\ vec {B}}}.}{\ displaystyle \ int e ^ {- {\ frac {1} {2}} \ sum \ пределы _ {i, j = 1} ^ {n} A_ {ij} x_ {i} x_ {j} + \ sum \ limits _ {i = 1} ^ {n} B_ {i} x_ {i}} d ^ {n} x = \ int e ^ {- {\ frac {1} {2}} {\ vec {x}} ^ {T} \ mathbf {A} {\ vec {x}} + {\ vec { B}} ^ {T} {\ vec {x}}} d ^ {n} x = {\ sqrt {\ frac {(2 \ pi) ^ {n}} {\ det {A}}}} e ^ {{\ frac {1} {2}} {\ vec {B}} ^ {T} \ mathbf {A} ^ {- 1} {\ vec {B}}}.}

Интегралы аналогичной формы

∫ 0 ∞ x 2 ne - x 2 a 2 dx = π a 2 n + 1 (2 n - 1)! ! 2 n + 1 {\ displaystyle \ int _ {0} ^ {\ infty} x ^ {2n} e ^ {- {\ frac {x ^ {2}} {a ^ {2}}}} \, dx = {\ sqrt {\ pi}} {\ frac {a ^ {2n + 1} (2n-1) !!} {2 ^ {n + 1}}}}\ int _ {0} ^ {\ infty} x ^ {2n} e ^ {- {\ frac {x ^ {2}} {a ^ { 2}}}} \, dx = {\ sqrt {\ pi}} {\ frac {a ^ {2n + 1} (2n-1) !!} {2 ^ {n + 1}}}
∫ 0 ∞ x 2 n + 1 e - х 2 а 2 дх = п! 2 a 2 n + 2 {\ displaystyle \ int _ {0} ^ {\ infty} x ^ {2n + 1} e ^ {- {\ frac {x ^ {2}} {a ^ {2}}}} \, dx = {\ frac {n!} {2}} a ^ {2n + 2}}\ int _ {0} ^ {\ infty} x ^ {2n + 1} e ^ {- {\ frac { x ^ {2}} {a ^ {2}}}} \, dx = {\ frac {n!} {2}} a ^ {2n + 2}
∫ 0 ∞ x 2 ne - ax 2 dx = (2 n - 1)! ! an 2 n + 1 π a {\ displaystyle \ int _ {0} ^ {\ infty} x ^ {2n} e ^ {- ax ^ {2}} \, dx = {\ frac {(2n-1)! !} {a ^ {n} 2 ^ {n + 1}}} {\ sqrt {\ frac {\ pi} {a}}}}\ int _ {0} ^ {\ infty} x ^ {2n} e ^ {- ax ^ {2 }} \, dx = {\ frac {(2n-1) !!} {a ^ {n} 2 ^ {n + 1}}} {\ sqrt {\ frac {\ pi} {a}}}
∫ 0 ∞ x 2 n + 1 e - ax 2 dx = п! 2 an + 1 {\ displaystyle \ int _ {0} ^ {\ infty} x ^ {2n + 1} e ^ {- ax ^ {2}} \, dx = {\ frac {n!} {2a ^ { n + 1}}}}{\ displaystyle \ int _ {0} ^ {\ infty} x ^ {2n + 1} e ^ {- ax ^ {2}} \, dx = {\ frac {n!} {2a ^ {n + 1}}}}
∫ 0 ∞ xne - ax 2 dx = Γ (n + 1 2) 2 an + 1 2 {\ displaystyle \ int _ {0} ^ {\ infty} x ^ {n} e ^ {- ax ^ {2}} \, dx = {\ frac {\ Gamma ({\ frac {n + 1} {2}})} {2a ^ {\ frac {n + 1} {2}} }}}{\ displaystyle \ int _ {0} ^ {\ infty} x ^ {n} e ^ {-ax ^ {2}} \, dx = {\ frac {\ Gamma ({\ frac {n + 1} {2}})} {2a ^ {\ frac {n + 1} {2}}}} }

где n {\ displaystyle n}n - положительное целое число, а ! ! {\ displaystyle !!}!! обозначает двойной факториал.

Легкий способ получить их - дифференцировать под знаком интеграла.

∫ - ∞ ∞ x 2 ne - α x 2 dx = (- 1) n ∫ - ∞ ∞ ∂ n ∂ α ne - α x 2 dx = (- 1) n ∂ n ∂ α n ∫ - ∞ ∞ e - α x 2 dx = π (- 1) N ∂ N ∂ α N α - 1 2 знак равно π α (2 N - 1)! ! (2 α) п {\ Displaystyle {\ begin {align} \ int _ {- \ infty} ^ {\ infty} x ^ {2n} e ^ {- \ alpha x ^ {2}} \, dx = \ left (-1 \ справа) ^ {n} \ int _ {- \ infty} ^ {\ infty} {\ frac {\ partial ^ {n}} {\ partial \ alpha ^ {n}}} e ^ {- \ альфа x ^ {2}} \, dx = \ left (-1 \ right) ^ {n} {\ frac {\ partial ^ {n}} {\ partial \ alpha ^ {n}}} \ int _ {- \ infty} ^ {\ infty} e ^ {- \ alpha x ^ {2}} \, dx \\ [6pt] = {\ sqrt {\ pi}} \ left (-1 \ right) ^ {n} {\ frac {\ partial ^ {n}} {\ partial \ alpha ^ {n}}} \ alpha ^ {- {\ frac {1} {2}}} = {\ sqrt {\ frac {\ pi} { \ alpha}}} {\ frac {(2n-1) !!} {\ left (2 \ alpha \ right) ^ {n}}} \ end {align}}}{\ displaystyle {\ begin {align} \ int _ {- \ infty} ^ {\ infty} x ^ {2n} e ^ {- \ alpha x ^ {2}} \, dx = \ left (-1 \ right) ^ {n} \ int _ {- \ infty} ^ {\ infty} {\ frac {\ partial ^ {n}} {\ partial \ alpha ^ {n}}} e ^ {- \ alpha x ^ {2}} \, dx = \ left (- 1 \ right) ^ {n} {\ frac {\ partial ^ {n}} {\ partial \ alpha ^ {n}}} \ int _ {- \ infty} ^ {\ infty} e ^ {- \ alpha x ^ {2}} \, dx \\ [6pt] = {\ sqrt {\ pi}} \ left (-1 \ right) ^ {n} {\ frac {\ partial ^ {n}} {\ partial \ alpha ^ {n}}} \ alpha ^ {- {\ frac {1} {2}}} = {\ sqrt {\ frac {\ pi} {\ alpha}}} {\ frac {(2n-1)! !} {\ left (2 \ alpha \ right) ^ {n}}} \ end {align}}}

Можно также интегрировать по частям и найдите рекуррентное отношение, чтобы решить эту проблему.

Многочлены высшего порядка

Применение линейной замены базиса показывает, что интеграл экспоненты однородного многочлена от n переменных может зависеть только от SL (n) -инварианты полинома. Одним из таких инвариантов является дискриминант , нули которого отмечают особенности интеграла. Однако интеграл может также зависеть от других инвариантов.

Показатели других четных многочленов могут быть решены численно с использованием рядов. Их можно интерпретировать как формальные вычисления, когда нет сходимости. Например, решение интеграла от экспоненты полинома четвертой степени равно

∫ - ∞ ∞ eax 4 + bx 3 + cx 2 + dx + fdx = 1 2 ef ∑ n, m, p = 0 n + p = 0 mod 2 ∞ bnn! с м м! д п п! Г (3 п + 2 м + р + 1 4) (- а) 3 п + 2 м + р + 1 4. {\ displaystyle \ int _ {- \ infty} ^ {\ infty} e ^ {ax ^ {4} + bx ^ {3} + cx ^ {2} + dx + f} \, dx = {\ frac {1 } {2}} e ^ {f} \ sum _ {\ begin {smallmatrix} n, m, p = 0 \\ n + p = 0 \ mod 2 \ end {smallmatrix}} ^ {\ infty} {\ frac {b ^ {n}} {n!}} {\ frac {c ^ {m}} {m!}} {\ frac {d ^ {p}} {p!}} {\ frac {\ Gamma \ left ({\ frac {3n + 2m + p + 1} {4}} \ right)} {(- a) ^ {\ frac {3n + 2m + p + 1} {4}}}}.}{\ displaystyle \ int _ {- \ infty} ^ {\ infty} e ^ {ax ^ {4} + bx ^ {3} + cx ^ {2} + dx + f} \, dx = { \ frac {1} {2}} e ^ {f} \ sum _ {\ begin {smallmatrix} n, m, p = 0 \\ n + p = 0 \ mod 2 \ end {smallmatrix}} ^ {\ infty } {\ frac {b ^ {n}} {n!}} {\ frac {c ^ {m}} {m!}} {\ frac {d ^ {p}} {p!}} {\ frac { \ Gamma \ left ({\ frac {3n + 2m + p + 1} {4}} \ right)} {(- a) ^ {\ frac {3n + 2m + p + 1} {4}}}}. }

Требование n + p = 0 mod 2 состоит в том, что интеграл от −∞ до 0 дает коэффициент (−1) / 2 для каждого члена, а интеграл от 0 до + ∞ дает коэффициент 1/2 для каждого члена. Эти интегралы появляются в таких предметах, как квантовая теория поля.

См. Также

  • icon Математический портал
  • icon Физический портал

Ссылки

Цитаты

Источники

Последняя правка сделана 2021-05-21 13:12:57
Содержание доступно по лицензии CC BY-SA 3.0 (если не указано иное).
Обратная связь: support@alphapedia.ru
Соглашение
О проекте