Иллюстрированный самоучитель по MatLab

Гамма-функция и ее варианты

Гамма-функция определяется выражением:

Иллюстрированный самоучитель по MatLab › Специальные математические функции › Гамма-функция и ее варианты

Неполная гамма-функция определяется как:

Иллюстрированный самоучитель по MatLab › Специальные математические функции › Гамма-функция и ее варианты

  • gamma (А) – возвращает гамма-функцию элементов А. Аргумент А должен быть вещественным.
  • gamma inc(X,А) – возвращает неполную гамма-функцию соответствующих элементов X и А. Аргументы X и А должны быть вещественными и иметь одинаковый размер (или любой из них может быть скалярным).
  • gammaln(A) – возвращает логарифмическую гамма-функцию, gammaln(A) = log(gamma(A)). Команда gammaln позволяет избежать переполнения, которое может происходить, если вычислять логарифмическую гамма-функцию непосредственно, используя log(gamma(A)).

Примеры:

>> f=[5.3];d=gamma(f)
d =
24 2
>> h=gammaln(f)
h =
3.1781 0.6931

Гамма-функция имеет довольно "сложный", график, заслуживающий построения (рис. 9.2).

Иллюстрированный самоучитель по MatLab › Специальные математические функции › Гамма-функция и ее варианты
Рис. 9.2. График гамма-функции

Это можно осуществить с помощью следующего файла-сценария:

%Gamma function graphicclear syms x
ezplot(gamma(x).[-4 4]) grid on

Гамма-функция вычисляется по известному алгоритму W. J. Kody (1989 г.). Для вычисления неполной гамма-функции используются рекуррентные формулы.

Если Вы заметили ошибку, выделите, пожалуйста, необходимый текст и нажмите CTRL + Enter, чтобы сообщить об этом редактору.