Особые свойства Гамма-функции Эйлера (85973)

Посмотреть архив целиком

Реферат


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

Пояснительная записка к курсовой работе выполнена в объёме 36 листов. Она содержит таблицу значений гамма-функции при некоторых значениях переменных и тексты программ для вычисления значений Гамма-функции и для построения графика, а также 2 рисунка.

Для написания курсовой работы было использовано 7 источников.


Введение


Выделяют особый класс функций, представимых в виде собственого либо несобственого интеграла, который зависит не только от формальной переменной, а и от параметра.

Такие функции называются интегралами зависящими от параметра. К их числу относятся гамма и бета функции Эйлера.

Бета функции представимы интегралом Эйлера первого рода:

Гамма функция представляется интегралом Эйлера второго рода:

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

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

Эйлеровы интегралы представляют собой хорошо изученные неэлементарные функции. Задача считается решённой, если она приводится к вычислению эйлеровых интегралов.

  1. Бэта-функция Эйлера


Бэта – функции определяются интегралом Эйлера первого рода:

=(1.1)

Он представляет функцию от двух переменных параметров и : функцию B. Если эти параметры удовлетворяют условиям и ,то интеграл (1.1) будет несобственным интегралом, зависящим от параметров и ,причём особыми точками этого интеграла будут точки и


Интеграл (1.1) сходятся при .Полагая получим:


= - =


т.e. аргумент и входят в симметрично. Принимая во внимание тождество


по формуле интегрирования почестям имеем



Откуда получаем


=

(1.2)

При целом b = n последовательно применяя (1.2)

Получим


(1.3)



при целых = m,= n, имеем



но B(1,1) = 1,следовательно:


Положим в (1.1) .Так как график функции симметрична относительно прямой ,то




и в результате подстановки , получаем


полагая в(1.1) ,откуда , получим


(1.4)


разделяя интеграл на два в пределах от 0 до 1 и от 1 до и применение ко второму интегралу подстановки ,получим


2. Гамма-функция


2.1 Определение


Восклицательный знак в математических трудах обычно означает взятие факториала какого-либо целого неотрицательного числа:

n! = 1·2·3·...·n.

Функцию факториал можно еще записать в виде рекурсионного соотношения:

(n+1)! = (n+1)·n!.

Это соотношение можно рассматривать не только при целых значениях n.

Рассмотрим разностное уравнение

G(z+1)=zG(z).

(2.1)

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


    1. Интегральное представление


Перейдем к решению этого уравнения. Будем искать решение в виде интеграла Лапласа:

В этом случае правая часть уравнения (2.1) может быть записана в виде:


Эта формула справедлива, если существуют пределы для внеинтегрального члена. Заранее нам не известно поведение образа [(G)\tilde](p) при p® ±¥. Предположим, что образ гамма-функции таков, что внеинтегральное слагаемое равно нулю. После того, как будет найдено решение, надо будет проверить, верно ли предположение о внеинтегральном слагаемом, иначе придется искать G(z) как-нибудь по-другому.

Левая часть равенства (2.1) записывается следующим образом:

Тогда уравнение (2.1) для образа гамма-функции имеет вид:


Это уравнение легко решить:


(2.2)


Нетрудно заметить, что найденная функция [(Г)\tilde](p) на самом деле такова, что внеинтегральный член в формуле (2.2) равен нулю.

Зная образ гамма-функции, легко получить и выражение для прообраза:

Это неканоническая формула, для того, чтобы привести ее к виду, полученному Эйлером, надо сделать замену переменной интегрирования: t = exp(-p), тогда интеграл примет вид:

Постоянная C выбирается так, чтобы при целых значениях z гамма-функция совпадала с функцией факториал: Г(n+1) = n!, тогда:

следовательно C = 1. Окончательно, получаем формулу Эйлера для гамма-функции:

(2.3)

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

Проверить, что функция, определенная формулой (2.3), действительно удовлетворяет уравнению (2.1), можно, проинтегрировав интеграл в правой части этой формулы по частям:


2.3 Область определения и полюсы

В подынтегральной функции интеграла (2.3) при экспонента exp(-tz) при R(z) > 0 убывает гораздо быстрее, чем растет алгебраическая функция t(z-1). Особенность в нуле - интегрируемая, поэтому несобственный интеграл в (2.3) сходится абсолютно и равномерно при R (z) > 0. Более того, последовательным дифференцированием по параметру z легко убедиться, что Г(z) - голоморфная функция при R (z) > 0. Однако, непригодность интегрального представления (2.3) при R (z) 0 не означает, что там не определена сама гамма-функция - решение уравнения (2.1).

Рассмотрим поведение Г(z) в окрестности нуля. Для этого представим:

где - голоморфная функция в окрестности z = 0. Из формулы (2.1) следует:


Тогда

то есть Г(z) имеет полюс первого порядка при z = 0.

Также легко получить:

то есть в окрестности точки функция Г(z) также имеет полюс первого порядка.

Таким же образом можно получить формулу:

(2.4)

Из этой формулы следует, что точки z = 0,-1,-2,... - простые полюсы гамма-функции и других полюсов на вещественной оси эта функция не имеет. Нетрудно вычислить вычет в точке z = -n, n = 0,1,2,...:


    1. Представление Ганкеля через интеграл по петле


Выясним, имеет ли гамма-функция нули. Для этого рассмотрим функцию

Полюсы этой функции и есть нули функции Г(z).

Разностное уравнение для I(z) легко получить, воспользовавшись выражением для Г(z):

Выражение для решения этого уравнения в виде интеграла можно получить так же, как было получено интегральное выражение для гамма-функции - через преобразование Лапласа. Ниже приведены вычисления.ни такие же, как и в п.1).ии теграла будут точки

или


После разделения переменных получим:


Проинтегрировав получаем:


или


Переход к прообразу Лапласа дает:

В полученном интеграле сделаем замену переменной интегрирования:

тогда

Здесь важно заметить, что подынтегральная функция при нецелых значениях z имеет точку ветвления t = 0. На комплексной плоскости переменной t проведем разрез по отрицательной вещественной полуоси. Интеграл по этой полуоси представим как сумму интеграла по верхнему берегу этого разреза от до 0 и интеграла от 0 до по нижнему берегу разреза. Чтобы интеграл не проходил через точку ветвления, устроим вокруг нее петлю.




Рис1: Петля в интегральном представлении Ганкеля.


В результате получим:

Чтобы выяснить значение постоянной, вспомним, что I(1) = 1, с другой стороны:

Интегральное представление

(2.5)

называется представлением Ганкеля по петле.

Легко видеть, что функция 1/Г(z) не имеет полюсов в комплексной плоскости, следовательно, гамма-функция не имеет нулей.

С помощью этого интегрального представления можно получить формулу для произведения гамма-функций. Для этого в интеграле сделаем замену переменной , тогда:


то есть


    1. Предельная форма Эйлера


Гамма-функцию можно представить в виде бесконечного произведения. Это можно заметить, если в интеграле (2.3) представить

Тогда интегральное представление гамма-функции:

В этой формуле мы можем поменять пределы - предел интегрирования в несобственном интеграле и предел при внутри интеграла. Приведем результат:

Возьмем по частям этот интеграл:


Если провести эту процедуру n раз, получим:



Переходя к пределу, получим предельную форму Эйлера для гамма-функции:



Случайные файлы

Файл
150259.rtf
31590.rtf
29909.rtf
181796.rtf
22311.rtf




Чтобы не видеть здесь видео-рекламу достаточно стать зарегистрированным пользователем.
Чтобы не видеть никакую рекламу на сайте, нужно стать VIP-пользователем.
Это можно сделать совершенно бесплатно. Читайте подробности тут.