Методы, которые мы применяем к биномиальным коэффициентам,
оказываются весьма эффективными в тех случаях, когда они работают.
Однако важно помнить, что их успешность часто носит специфический характер — эти методы, как правило,
не универсальны и требуют тщательного подбора под конкретную задачу.
Занимаясь решением разнообразных задач, мы иногда сталкиваемся с тем, что подходы, которые казались перспективными,
неожиданно приводят к тупиковым путям.
Биномиальные коэффициенты ведут себя подобно хамелеонам: их внешний вид легко меняется в зависимости от контекста,
что делает их одновременно гибкими и сложными для систематизации.
Поэтому естественным становится вопрос: существует ли некий унифицирующий объект,
который обобщил бы всё многообразие методов суммирования?
К счастью, такой объект действительно существует — им являются так называемые гипергеометрические ряды.
Начало изучению гипергеометрических рядов было положено в конце 18 – начале 19 века.
Первые систематические исследования появились в работах Леонарда Эйлера в 1760-х годах,
где он рассмотрел ряды с рекуррентным соотношением для коэффициентов.
Окончательную форму гипергеометрической функции и глубокий анализ её свойств дал Карл Фридрих Гаусс
в знаменитой работе 1812 года, посвящённой сходимости и функциональным соотношениям.
Позже, в середине XIX века, Бернард Риман развил теорию с помощью комплексного анализа,
введя геометрический взгляд на гипергеометрическое дифференциальное уравнение.
И до сих пор гипергеометрические ряды остаются предметом активных исследований —
они находят применение в комбинаторике, теории чисел,
математической физике и даже в современных вычислительных алгоритмах.
Гипергеометрическая функция и ряд
Гипергеометрическая функция — функция
от комплексной переменной z с p+q параметрами
Параметры a1,a2,…,ap называются верхними параметрами,
а параметры b1,b2,…,bq называются нижними параметрами.
Для того, чтобы знаменатель не превращался в 0,
каждый нижний параметр не может быть нулём или целым отрицательным числом.
Строчная запись выглядит так: F(a1,a2,…,ap;b1,b2,…,bq∣z).
Обращаю ваше внимание на то, что верхние и нижние параметры разделяются точкой с запятой,
а аргумент, как и в полной нотации, отделяется вертикальной линией.
Конкретное значение гипергеометрической функции в точке z называется гипергеометрическим рядом.
Огромное количество функций являются частными случаями гипергеометрических.
Например, если параметры вообще отсутствуют, мы получаем экспоненту:
F(−−z)=n=0∑∞n!zn=ez
Можно получить обычный геометрический ряд, сделав в числителе множитель n!=1n:
здесь Jν(x) — функция Бесселя первого рода, Iν(x) — модифицированная функция Бесселя, γ(a,x) — неполная гамма-функция, Si(x) — интегральный синус, K(k) и E(k) — эллиптические интегралы первого и второго рода.
Пока гипергеометрические ряды нам не дали ничего, кроме громких имён и красивых определений.
Однако мы увидели, что множество весьма разных функций можно представить именно в виде гипергеометрических рядов.
Это и станет основным объектом нашего внимания.
Какие ряды являются гипергеометрическими?
На этот вопрос легко ответить, если рассмотреть отношение последовательных членов гипергеометрического ряда
У нас получилась рациональная функция по n, умноженная на z.
То есть если отношение соседних членов функционального ряда по z представляет собой рациональную функцию по номеру членов, умноженную на z,
то этот функциональный ряд является гипергеометрической функцией.
Аналогично, если отношение соседних членов обычного числового ряда представляет собой
рациональную функцию по номеру членов, то ряд является гипергеометрическим рядом.
Сразу скажу о способе упрощения выражений.
Поскольку верхние индексы идут в числитель, а нижние в знаменатель, одинаковые числа можно сокращать.
Например,
F(3,41,2,3z)=F(3,41,2,3z)=F(41,2z)
Давайте, чтобы не путаться с обозначениями и определениями, разберём несколько примеров.
Первый пример. Посмотрим на числовой ряд n=0∑∞tn, где
tntn+1=4n2+1n2+7n+10
Это рациональная функция по номеру n.
Значит, разложив числитель и знаменатель на множители, мы сможем найти параметры гипергеометрического ряда,
который выражает исходный числовой ряд.
Тогда наш исходный ряд легко выражается через гипергеометрический ряд
n=0∑∞tn=F(−i/2,i/21,2,541)
Второй пример. Посмотрим на ряд n=0∑∞tn, где
tntn+1=n2+3n+2n2−2nz−2n+6z−3
Это рациональная функция по n, но подходит ли она под наше определение?
Если рассматривать z как аргумент гипергеометрической функции, то нет, она под свойство не подходит.
Здесь z входит линейно в числитель, и его нельзя выделить.
Значит, n=0∑∞tn не является гипергеометрической функцией по z.
Тем не менее, этот ряд является гипергеометрическим рядом:
Теперь, выделив отдельно верхние и нижние параметры, получим, что
n=0∑∞tn=F(21−2z,−31)
Ещё раз подчеркну, что это не гипергеометрическая функция от переменной z,
а гипергеометрический ряд, то есть конкретное значение гипергеометрической функции в точке 1.
Здесь переменная z присутствует в записи одного из параметров.
Третий пример. Рассмотрим разложение тангенса в ряд Тейлора вокруг нуля:
tanz=n=1∑∞(2n)!B2n(−4)n(1−4n)⋅z2n−1
Попробуем вычислить отношение соседних ненулевых коэффициентов.
Обозначим un — коэффициент при z2n−1, то есть
Тогда отношение соседних коэффициентов равно
Асимптотически числа Бернулли растут как ∣B2n∣∼4πn(n/πe)2n,
и их отношение B2n+2/B2n ведёт себя как n2 при больших n.
Но тут даже без асимптотики ясно, что отношение un+1/un не является рациональной функцией от n, поскольку числа Бернулли не задаются рациональной рекурсией.
Следовательно, ряд для tanz не является гипергеометрическим рядом,
поскольку отношение его коэффициентов не является рациональной функцией от номера.
Это особенно примечательно, ведь tanz=sinz/cosz,
а обе функции в правой части гипергеометрические.
Однако, как мы убедились только что, отношение гипергеометрических функций
может и не является гипергеометрической функцией.
Четвёртый пример. А что там с синусом? Там всё гораздо лучше тангенса.
Обратите внимание на то, что я вынес множитель z.
При z=0 значение синуса равно 0,
а значение любого гипергеометрического ряда в точке 0 равно 1.
После вынесения z в разложении синуса я как раз получил ряд,
значение которого в точке 0 равно 1,
а значит приведение этого ряда к гипергеометрическому будет корректно.
Давайте посмотрим, как выглядит правило суммирования по диагонали для биномиальных коэффициентов
в гипергеометрической форме
k=0∑n(kr+k)=(nr+n+1)
Для начала сделаем из конечной суммы бесконечную, заменив k на n−k
Несмотря на то, что сумма формально бесконечная, она на самом деле конечная:
множитель 1/(n−k)! равен 0 при всех k>n.
Подробно этот неестественный факт разобран в статье «Гамма-функция».
Посмотрим теперь на отношение соседних членов.
Будьте внимательны: здесь индексной переменной является k, а не n.
При этом первый (свободный) член ряда равен (nr+n).
Значит, исходная сумма превращается в гипергеометрический ряд
k=0∑n(kr+k)=(nr+n)⋅F(−n−r1,−n1)
Тогда, подставляя значение суммы из правила суммирования по диагонали, получаем
(nr+n)⋅F(−n−r1,−n1)=(nr+n+1)
И можно на поделить на биномиальный коэффициент, и получить чуть более приятную формулу
F(−n−r1,−n1)=r+1r+n+1если(nr+n)=0
Давайте теперь получим явную формулу для одного частного случая гипергеометрической функции —
функции с двумя верхними параметрами a, b и одним нижним параметром c.
Эта функция, на самом деле, исторически была первой ...
Посмотрим на свёртку Вандермонда
k=0∑n(kr)(n−ks)=(nr+s)дляr,s∈Cиn∈Z
Сразу запишем это тождество в гипергеометрической форме.
Отношение соседних членов равно
Здесь предлагаю немного отвлечься и поговорить о сходимости и об обобщениях.
Не знаю, как вам, но мне доставляет моральный дискомфорт полученная формула,
ведь у нас там параметры a и c являются любыми комплексными числами,
а параметр b только отрицательным целым числом.
Хочется какой-то полной симметрии. Давайте её получим.
Давайте узнаем, в каких случаях наш частный гипергеометрический ряд сходится.
Итак, смотрим на гипергеометрический ряд
F(ca,bz)=n=0∑∞ccan⋅bn⋅n!znпри∣z∣=1
Вычислим на асимптотику коэффициента
cc⋅n!an⋅bn∼Γ(a)⋅Γ(b)na+b−c−1приn→∞
Ряд n=0∑∞nα сходится при Reα<−1.
А это значит, что наш гипергеометрический ряд сходится при
Re(a+b−c−1)<−1⟹Re(c−a−b)>0
Далее нам нужно будет использовать частный случай интегрального представления Эйлера,
которое мы в общем виде изучим позднее. А сейчас я просто покажу, как нужно проводить рассуждения.
Заведём переменную t∈[0,1].
Для этой переменной можно разложить выражение (1−t)−a в биномиальный ряд
(1−t)−a=n=0∑∞n!an⋅tn
Этот ряд сходится равномерно при всех значениях t, то есть то, что я дальше буду делать, законно.
Теперь рассмотрим интеграл
0∫1tb−1(1−t)c−b−1(1−t)−adt
Подставим вместо (1−t)−a наш ряд и поменяем знаки интегрирования и суммирования местами
У нас получилась то же самое, что получалось при рассмотрении свёртки Вандермонда,
однако здесь больше нет требования целочисленности параметра b.
Единственное требование на параметры, которое тут есть — гипергеометрический ряд должен сходиться,
то есть Rec>Rea+Reb.
Дифференциальные соотношения
При исследовании функций нам обычно интересны не только их конкретные значения,
но и их рост, поведение при разных аргументах и так далее.
На большинство интересующих нас вопросов могут ответить производные
и дифференциальные уравнения, которые получаются из этих функций.