Русская Википедия:Производящая функция последовательности

Материал из Онлайн справочника
Перейти к навигацииПерейти к поиску

Шаблон:Другие значения Производя́щая фу́нкция после́довательности — алгебраическое понятие, которое позволяет работать с разными комбинаторными объектами аналитическими методами. Они дают гибкий способ описывать соотношения в комбинаторике, а иногда помогают вывести явные формулы для числа комбинаторных объектов определённого типа.

Если дана последовательность <math>\{ a_n \}</math> чисел <math>a_0, a_1, a_2, a_3, \ldots</math>, то из них можно построить формальный степенной ряд

<math>A(t)=\sum_{n=0}^\infty a_n t^n=a_0+a_1t+a_2t^2+a_3t^3+\dotsb</math>,

который называется производящей функцией <math>A(t)</math> этой последовательности.

Близким понятием является экспоненциальная производящая функция последовательности <math>\{a_n\}</math> — степенной ряд

<math>A(t)=\sum_{n=0}^\infty \frac{a_n t^n}{n!}=a_0+a_1t+\frac{a_2}{2}t^2+\frac{a_3}{6}t^3+\dotsb</math>,

у которого коэффициент перед <math>a_n</math> поделён на факториал числа <math>n</math>.

Замечания

Зачастую производящая функция последовательности чисел является рядом Тейлора некоторой аналитической функции, что может использоваться для изучения свойств самой последовательности. Однако, в общем случае производящая функция не обязана быть аналитической. Например, оба ряда

<math>\sum_{n=0}^\infty (3^n)^n x^n</math> и <math>\sum_{n=0}^\infty (2^n)^n x^n</math>

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

Свойства

  • Производящая функция суммы (или разности) двух последовательностей равна сумме (или разности) соответствующих производящих функций.
  • Произведение производящих функций <math>A(x)=\sum_{n=0}^\infty a_n x^n</math> и <math>B(x)=\sum_{n=0}^\infty b_n x^n</math> последовательностей <math>\{a_n\}</math> и <math>\{b_n\}</math> является производящей функцией свёртки <math>c_n = \sum_{k=0}^n a_k b_{n-k}</math> этих последовательностей:
<math>A(x)B(x)=\sum_{n=0}^\infty c_n x^n.</math>
  • Если <math>A(x)=\sum_{n=0}^\infty a_n \frac{x^n}{n!}</math> и <math>B(x)=\sum_{n=0}^\infty b_n \frac{x^n}{n!}</math> — экспоненциальные производящие функции последовательностей <math>\{a_n\}</math> и <math>\{b_n\}</math>, то их произведение <math>A(x)\cdot B(x)=\sum_{n=0}^\infty c_n \frac{x^n}{n!}</math> является экспоненциальной производящей функцией последовательности <math>c_n = \sum_{k=0}^n {n\choose k} a_k b_{n-k}</math>.

Примеры использования

В комбинаторике

Число композиций

Пусть <math>B_m^n</math> — это количество композиций неотрицательного целого числа n длины m, то есть, представлений n в виде <math>n=k_1+k_2+\cdots+k_m</math>, где <math>\{k_i\}</math> — неотрицательные целые числа. Число <math>B_m^n</math> также является числом сочетаний с повторениями из m по n, то есть, количество выборок n возможно повторяющих элементов из множества <math>\{ 1, 2, \dots, m\}</math> (при этом каждый член <math>k_i</math> в композиции можно трактовать как количество элементов i в выборке).

При фиксированном m производящей функцией последовательности <math>B_m^0, B_m^1, \dots</math> является:

<math>\sum_{n=0}^\infty B_m^n x^n=(1+x+x^2+\cdots)^m=(1-x)^{-m}.</math>

Поэтому число <math>B_m^n</math> может быть найдено как коэффициент при <math>x^n</math> в разложении <math>(1-x)^{-m}</math> по степеням x. Для этого можно воспользоваться определением биномиальных коэффициентов или же непосредственно взять n раз производную в нуле:

<math>B_m^n = (-1)^n {-m\choose n} = \frac{m(m+1)\cdots(m+n-1)}{n!} = {m+n-1 \choose n}.</math>
Число связных графов

Обозначим через <math>a_n</math> число всех графов с вершинами <math>\{1,\dots,n\}</math> и через <math>c_n</math> число всех связных графов с этими вершинами.

Заметим, что <math>a_n=2^{\binom n 2}</math>. В частности легко посчитать первые члены этой последовательности

<math>1,\ 2,\ 8,\ 64,\ 1024,\ 32768,\ 2097152,\ \dots</math>

Рассмотрим экспоненциальные производящие функции этих последовательностей:

<math>a(x)=a_1\cdot x+\tfrac12\cdot a_2\cdot x^2+\dots+\tfrac1{n!}\cdot a_n\cdot x^n+\dots</math>
<math>c(x)=c_1\cdot x+\tfrac12\cdot c_2\cdot x^2+\dots+\tfrac1{n!}\cdot c_n\cdot x^n+\dots</math>

Оба ряда расходятся при <math>x\ne 0</math>, тем не менее их можно рассматривать как формальные степенные ряды и для этих рядов выполняется следующее соотношение:

<math>1+a(x)=e^{c(x)},</math>

из которого следует простое рекуррентное соотношение для <math>c_n</math>, позволяющее быстро найти первые члены этой последовательности[1]

<math>1,\ 1,\ 4,\ 38,\ 728,\ 26704,\ 1866256,\ 251548592,\ 66296291072,\ \dots</math>

В теории вероятностей

<math>\mathbb{P}(X=j) = p_j,\; j=0,1,...;\quad \sum\limits_{j=0}^{\infty} p_j = 1</math>

то её математическое ожидание может быть выражено через производящую функцию последовательности <math>\{p_i\}</math>

<math>P(s)=\sum_{k=0}^\infty\;p_k s^k, </math>

как значение первой производной в единице: <math>M[X] = P'(1)</math> (стоит отметить, что ряд для P(s) сходится, по крайней мере, при <math>|s|\le 1</math>). Действительно,

<math>P'(s)=\sum_{k=1}^\infty{k p_k s^{k-1}}.</math>

При подстановке <math>s=1</math> получим величину <math>P'(1)=\sum_{k=1}^\infty{k p_k}</math>, которая по определению является математическим ожиданием дискретной случайной величины. Если этот ряд расходится, то <math>\lim_{s\to 1}P'(s)=\infty</math> — а <math>X</math> имеет бесконечное математическое ожидание, <math>P'(1)=M[X]=\infty</math>

  • Теперь возьмём производящую функцию <math>Q(s)</math> последовательности «хвостов» распределения <math>\{q_k\}</math>
<math>q_k=\mathbb{P}(X>k)=\sum_{j=k+1}^\infty{p_j};\quad Q(s)=\sum_{k=0}^\infty\;q_k s^k.</math>

Эта производящая функция связана с определённой ранее функцией <math>P(s)</math> свойством: <math>Q(s)=\frac{1-P(s)}{1-s}</math> при <math>|s|<1</math>. Из этого по теореме о среднем следует, что математическое ожидание равно просто значению этой функции в единице:

<math>M[X]=P'(1)=Q(1)</math>
  • Дифференцируя <math>P'(s)=\sum_{k=1}^\infty{k p_k s^{k-1}}</math> и используя соотношение <math>P'(s)=Q(s)-(1-s)Q'(s)</math>, получим:
<math>M[X(X-1)]=\sum{k(k-1)p_k}=P(1)=2Q'(1)</math>

Чтобы получить дисперсию <math>D[X]</math>, к этому выражению надо прибавить <math>M[X]-M^2[X]</math>, что приводит к следующим формулам для вычисления дисперсии:

<math>D[X]=P(1)+P'(1)-P'^2(1)=2Q'(1)+Q(1)-Q^2(1)</math>.

В случае бесконечной дисперсии <math>\lim_{s\to 1}P(s)=\infty</math>.

Вариации и обобщения

Производящая функция Дирихле

Производящая функция Дирихле последовательности <math>\{a_n\}</math> — это формальный ряд Дирихле

<math>\sum_{n=1}^\infty \frac{a_n}{n^s}</math>.
  • Производящей функцией Дирихле последовательности единиц 1,1,… является дзета-функция Римана:
    <math>\zeta(s) = \sum_{n=1}^\infty \frac{1}{n^s}.</math>
  • Если <math>\alpha(s)=\sum_{n=1}^\infty \frac{a_n}{n^s}</math> и <math>\beta(s)=\sum_{n=1}^\infty \frac{b_n}{n^s}</math> — производящие функции Дирихле последовательностей <math>\{a_n\}</math> и <math>\{b_n\}</math>, то их произведение <math>\alpha(s)\cdot \beta(s)=\sum_{n=0}^\infty \frac{c_n}{n^s}</math> является производящей функцией Дирихле свёртки Дирихле — последовательности <math>c_n = \sum_{d\mid n} a_d b_{n/d}</math>.

История

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

Примечания

Шаблон:Примечания

Литература

Ссылки