Специальные функции

Мотивация: функции, выходящие за рамки элементарных

Решение многих задач математической физики — дифференциальных уравнений с переменными коэффициентами — приводит к функциям, не выражаемым через многочлены, тригонометрические функции и экспоненту. Их называют специальными: гамма-функция, функции Бесселя, гипергеометрические функции. Комплексный анализ даёт единую теорию этих функций.

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

Определение (при Re s > 0): Γ(s) = ∫₀^∞ t^{s−1} e^{−t} dt.

Аналитическое продолжение: Γ(s) продолжается на ℂ {0, −1, −2, ...} с простыми полюсами в s = 0, −1, −2, ... и вычетами Res_{s=−n} Γ = (−1)^n/n!

Функциональное уравнение: Γ(s+1) = s·Γ(s). Следовательно: Γ(n) = (n−1)! для натуральных n.

Формула отражения (Эйлер): Γ(s)·Γ(1−s) = π/sin(πs).

Формула дублирования (Лежандр): Γ(s)·Γ(s+1/2) = √π · Γ(2s) / 2^{2s−1}.

Бета-функция

Определение: B(a,b) = ∫₀¹ t^{a−1}(1−t)^{b−1} dt = Γ(a)Γ(b)/Γ(a+b).

Используется в статистике (бета-распределение), комбинаторике.

Функции Бесселя

Уравнение Бесселя: z²w'' + zw' + (z²−ν²)w = 0.

Возникает при разделении переменных уравнения Гельмгольца в полярных координатах. Решения — функции Бесселя первого рода Jν(z):

Jν(z) = Σₘ₌₀^∞ (−1)^m (z/2)^{2m+ν} / (m! Γ(m+ν+1)).

Нули Jν(z) — дискретные (определяют собственные частоты барабана, трубы с круглым сечением).

Численный пример

Задача: Вычислить Γ(1/2) используя формулу отражения.

Шаг 1. По формуле отражения при s = 1/2: Γ(1/2) · Γ(1 − 1/2) = π / sin(π/2) = π/1 = π.

Шаг 2. Γ(1/2) · Γ(1/2) = π → Γ(1/2)² = π.

Шаг 3. Γ(1/2) = √π ≈ 1.7725.

Проверка через интеграл: Γ(1/2) = ∫₀^∞ t^{-1/2} e^{-t} dt. Замена t = u²: = 2∫₀^∞ e^{-u²} du = 2 · (√π/2) = √π ✓.

Вычисление Γ(3/2): Γ(3/2) = (1/2)·Γ(1/2) = √π/2 ≈ 0.886.

Бета-функция: B(1/2, 1/2) = Γ(1/2)²/Γ(1) = π/1 = π. Это также равно ∫₀¹ t^{-1/2}(1−t)^{-1/2}dt = ∫₀¹ 1/√(t(1−t))dt — замена t = sin²θ → 2∫₀^{π/2}dθ = π ✓.

Приложение: Площадь n-мерного единичного шара: Vₙ = π^{n/2}/Γ(n/2+1). При n=2: V₂ = π^1/Γ(2) = π (площадь единичного круга) ✓. При n=3: V₃ = π^{3/2}/Γ(5/2) = π^{3/2}/(3√π/4) = 4π/3 ✓.

Реальное приложение

Физика: квантовый гармонический осциллятор описывается функциями Эрмита Hₙ(x) = (−1)ⁿ e^{x²} dⁿ/dxⁿ e^{−x²}, которые ортогональны с весом e^{−x²}. Нормировочный интеграл ∫_{-∞}^∞ Hₙ² e^{-x²} dx = 2ⁿ n! √π — вычисляется через Γ-функцию. Дискретные энергетические уровни Eₙ = ℏω(n+1/2) соответствуют этим функциям.

Дополнительные аспекты

Гамма-функция Γ(z) = ∫₀^∞ t^{z−1} e^{−t} dt продолжается через функциональное уравнение Γ(z+1) = z·Γ(z) на всю плоскость, кроме z = 0, −1, −2, …, где имеет простые полюса с вычетами (−1)^n/n!. Формула отражения Γ(z)·Γ(1−z) = π/sin πz связывает Γ с тригонометрическими функциями. Бета-функция B(p,q) = Γ(p)Γ(q)/Γ(p+q) выражает интегралы вида ∫₀¹ t^{p−1}(1−t)^{q−1} dt и активно используется в статистике (бета-распределение, байесовский анализ). Дзета-функция ζ(s) удовлетворяет функциональному уравнению, симметричному относительно прямой Re s = 1/2; знаменитая гипотеза Римана утверждает, что все нетривиальные нули лежат на этой прямой.

Связь с другими разделами математики

Теория специальных функций тесно переплетена с линейными дифференциальными уравнениями. Классическая схема Фукса описывает уравнения с регулярными особенными точками; работа Эмиля Пикара и Эрнста Шарля развивает классификацию таких уравнений по монодромии. Гипергеометрическая функция Гаусса появляется как универсальное решение второго порядка с тремя регулярными особенностями; многие специальные функции (Лежандра, Бесселя, Якоби) реализуются как частные случаи, что делает теорему Римана о гипергеометрическом уравнении ключевым связующим результатом.

Алгебраические структуры возникают через теорию представлений. Полиномы Лежандра, Чебышёва, полиномы Якоби появляются как матричные элементы неприводимых представлений групп SO(3), SU(2), SU(1,1). Работы Германа Вейля и Гарольда Бумби связывают сферические функции с гармоническим анализом на компактных группах. В более современной постановке, развитой Жаком Диксомье и Николаем Виленкиным, специальные функции описываются как собственные функции касимирова оператора в универсальной обёртывающей алгебре.

Связь с топологией проявляется через интегралы по контурам, циклы в гомологии и теорию периодов. Исследования Картье и Делиня показывают, что многие специальные функции можно рассматривать как периоды мотивов; гипергеометрические интегралы дают примеры периодов, описывающих монодромию над многообразиями, задаваемыми алгебраическими уравнениями.

В теории вероятностей гамма- и бета-функции лежат в основе гамма- и бета-распределений, а полиномы Эрмита, Лагерра и Чебышёва реализуются как ортогональные полиномы для законов Гаусса, распределения Пуассона, арксин-распределения. Теорема Мигдала–Крамера о предельных распределениях в ортогональных разложениях использует эти системы в изучении линейных функционалов от случайных величин.

Численные методы опираются на асимптотические разложения, разработанные Джеффрисом и Ольвером, и на рекуррентные соотношения. В современных библиотеках (например, спецификация Boost.Math) реализация гамма-функции и функций Бесселя базируется на комбинации рядов, асимптотик и рациональных аппроксимаций Паде. Теорема Стирлинга служит фундаментом для устойчивых алгоритмов вычисления логарифма гамма-функции.

Историческая справка и развитие идеи

Начало систематическому изучению специальных функций положили работы Леонарда Эйлера в 1720–1740-х годах: он вводит гамма-функцию, исследует её продолжение и формулу отражения. Жозеф Фурье, изучая теплопроводность (Mémoire sur la propagation de la chaleur, 1822), выводит тригонометрические ряды и сталкивается с уравнениями, приводящими к функциям Бесселя, хотя терминология ещё не устоялась. Фридрих Бессель в 1820–1830-х годах исследует орбиты планет и кометы Галлея и публикует статьи в Astronomische Nachrichten, где появляются функции, носящие его имя. Карл Густав Якоб Якоби и Адольф Гурвиц разрабатывают теорию эллиптических функций и интегралов, которые можно рассматривать как ранний прототип современных специальных функций, связанных с комплексными многообразиями. Карл Гаусс в труде Disquisitiones generales circa seriem infinitam (1813) вводит гипергеометрический ряд и закладывает основу унифицированной теории. В середине XIX века Эрмит, Лежандр, Лагерр вводят системы ортогональных полиномов как решения собственных задач для дифференциальных операторов, мотивируемые задачами по небесной механике и теории потенциальных полей. В конце XIX – начале XX века Лебег и Хилберта формализуют понятие ортогонального разложения в пространстве L2, что придаёт специальным функциям статус базисов в функциональных пространствах.