Formule de Faà di Bruno

La forme la plus connue de cette formule est sans doute

${\displaystyle {d^{n} \over dx^{n}}f(g(x))=\sum {\frac {n!}{m_{1}!\,1!^{m_{1}}\,m_{2}!\,2!^{m_{2}}\,\cdots \,m_{n}!\,n!^{m_{n}}}}f^{(m_{1}+\cdots +m_{n})}(g(x))\prod _{j=1}^{n}\left(g^{(j)}(x)\right)^{m_{j}},}$

où la somme parcourt tous les n-uplets (m₁, ..., m_n) vérifiant la contrainte : ${\displaystyle 1m_{1}+2m_{2}+3m_{3}+\cdots +nm_{n}=n.}$

Pour pouvoir la retenir plus facilement, on l'énonce parfois

${\displaystyle {d^{n} \over dx^{n}}f(g(x))=\sum {\frac {n!}{m_{1}!\,m_{2}!\,\cdots \,m_{n}!}}f^{(m_{1}+\cdots +m_{n})}(g(x))\prod _{j=1}^{n}\left({\frac {g^{(j)}(x)}{j!}}\right)^{m_{j}},}$

forme dans laquelle on reconnaît la composition des séries de Taylor, mais sous cette forme, l'interprétation combinatoire des coefficients discutée plus bas est plus difficile à percevoir.

Combinant les termes correspondants à la même valeur de m₁ + m₂ + ... + m_n = k et remarquant que m_j doit être nul pour j > n − k + 1, on arrive à une autre formule un peu plus simple, exprimée en fonction des polynômes de Bell B_n,k(x₁,...,x_n−k+1) :

${\displaystyle {d^{n} \over dx^{n}}f(g(x))=\sum _{k=1}^{n}f^{(k)}(g(x))B_{n,k}\left(g'(x),g''(x),\dots ,g^{(n-k+1)}(x)\right).}$

La formule peut s'écrire sous la forme « combinatoire » suivante :

${\displaystyle {d^{n} \over dx^{n}}f(g(x))=(f\circ g)^{(n)}(x)=\sum _{\pi \in \Pi }f^{(\left|\pi \right|)}(g(x))\cdot \prod _{B\in \pi }g^{(\left|B\right|)}(x),}$ où

• π parcourt l'ensemble Π de toutes les partitions de l'ensemble { 1, ..., n },

• « B ∈ π » signifie que la variable B parcourt la liste de tous les blocs de la partition π, et

• |A| désigne le cardinal de l'ensemble A (et donc |π| est le nombre de blocs de la partition π, et |B| est la taille du bloc B).

La forme combinatoire peut sembler inutilisable au premier abord, aussi nous allons montrer sur un cas concret à quoi elle ressemble :

${\displaystyle {\begin{aligned}(f\circ g)''''(x)&=f''''(g(x))g'(x)^{4}+6f'''(g(x))g''(x)g'(x)^{2}\\&{}\quad +\;3f''(g(x))g''(x)^{2}+4f''(g(x))g'''(x)g'(x)\\&{}\quad +\;f'(g(x))g''''(x).\end{aligned}}}$

Quelle est la règle ?

${\displaystyle {\begin{aligned}g'(x)^{4}&&\leftrightarrow &&1+1+1+1&&\leftrightarrow &&f''''(g(x))&&\leftrightarrow &&1\\\\g''(x)g'(x)^{2}&&\leftrightarrow &&2+1+1&&\leftrightarrow &&f'''(g(x))&&\leftrightarrow &&6\\\\g''(x)^{2}&&\leftrightarrow &&2+2&&\leftrightarrow &&f''(g(x))&&\leftrightarrow &&3\\\\g'''(x)g'(x)&&\leftrightarrow &&3+1&&\leftrightarrow &&f''(g(x))&&\leftrightarrow &&4\\\\g''''(x)&&\leftrightarrow &&4&&\leftrightarrow &&f'(g(x))&&\leftrightarrow &&1\end{aligned}}}$

Le facteur ${\displaystyle \scriptstyle g''(x)g'(x)^{2}\;}$, par exemple, correspond à la partition 2 + 1 + 1 de l'entier 4 (puisque nous cherchons la dérivée d'ordre 4), et le facteur ${\displaystyle \scriptstyle f'''(g(x))\;}$ qui l'accompagne venant de ce qu'il y a trois termes dans cette partition ; enfin, le coefficient 6 résulte de ce qu'il y a exactement 6 partitions d'un ensemble à 4 éléments de la forme « une paire et deux singletons ».

Les coefficients de Faà di Bruno correspondant à ces comptes de partitions ont une forme explicite ; en effet, le nombre de partitions d'un entier n correspondant à la somme (la partition entière)

${\displaystyle \displaystyle n=\underbrace {1+\cdots +1} _{m_{1}}\,+\,\underbrace {2+\cdots +2} _{m_{2}}\,+\,\underbrace {3+\cdots +3} _{m_{3}}+\cdots }$

est égal à

${\displaystyle {\frac {n!}{m_{1}!\,m_{2}!\,m_{3}!\,\cdots 1!^{m_{1}}\,2!^{m_{2}}\,3!^{m_{3}}\,\cdots }}.}$

Ces coefficients apparaissent aussi dans les polynômes de Bell, utilisés dans l'étude des cumulants.

Si f(x) = e^x, alors toutes les dérivées de f sont les mêmes, et sont facteurs communs de tous les termes de la formule. Quand g(x) est une fonction génératrice de cumulants, f(g(x)) est une fonction génératrice de moments, et le polynôme en les dérivées de g est le polynôme exprimant les moments en fonction des cumulants. Prenant ${\displaystyle g(x)=-x^{2}}$, on retrouve aussi de cette manière les polynômes d'Hermite.

Si f(x) = 1/x, on a ${\displaystyle f^{(n)}(x)=(-1)^{n}n!/x^{n+1}}$ ; on en déduit[5] la formule suivante pour la dérivée n-ième de 1/g :

${\displaystyle (1/g)^{(n)}(x)={\frac {n!}{g^{n+1}(x)}}\sum {\frac {(-1)^{n-m_{0}}(n-m_{0})!}{\displaystyle \prod _{i=1}^{n}(i!)^{m_{i}}\;m_{i}!}}\prod _{i=0}^{n}\left(g^{(i)}(x)\right)^{m_{i}}}$

où ${\displaystyle m_{0}=m_{2}+2m_{3}+\dots +(n-1)m_{n}}$ est tel que : ${\displaystyle n=m_{0}+m_{1}+\dots +m_{n}=m_{1}+2m_{2}+\dots +nm_{n}}$.

Soit y = g(x₁, ..., x_n). Alors l'identité suivante est vraie, que les n variables soient distinctes, ou identiques, ou partitionnées en classes de variables indiscernables (l'exemple concret ci-dessous devrait rendre cela clair) :

${\displaystyle {\partial ^{n} \over \partial x_{1}\cdots \partial x_{n}}f(y)=\sum _{\pi \in \Pi }f^{(\left|\pi \right|)}(y)\cdot \prod _{B\in \pi }{\partial ^{\left|B\right|}y \over \prod _{j\in B}\partial x_{j}}}$[6],

où, comme précédemment, π parcourt l'ensemble Π de toutes les partitions de { 1, ..., n }, B ∈ π signifie que B parcourt les blocs de la partition π, et |A| est le cardinal de l'ensemble A. Une généralisation supplémentaire, correspondant au cas où Y est une variable vectorielle, est due à Tsoy-Wo Ma[7].

Si ${\displaystyle g(X)=\sum _{\ell =1}^{\infty }{b_{\ell } \over \ell$ !}X^{\ell }{\text{ et }}f(X)=\sum _{n=1}^{\infty }{a_{n} \over n!}X^{n}{\text{ et }}g(f(X))=h(X)=\sum _{n=1}^{\infty }{c_{n} \over n!}X^{n},} alors le coefficient c_n (qui est la dérivée formelle n-ième en 0 de la série formelle h) est donné par : ${\displaystyle c_{n}=\sum _{\pi =\left\{\,B_{1},\,\dots ,\,B_{k}\,\right\}}a_{\left|B_{1}\right|}\cdots a_{\left|B_{k}\right|}b_{k}}$ (avec les mêmes conventions que précédemment : π parcourt l'ensemble des partitions de l'ensemble { 1, ..., n }, B₁, ..., B_k sont les blocs de la partition π, et | B_j | est le nombre d'éléments du j-ème bloc, pour j = 1, ..., k).

Cette version de la formule est particulièrement bien adaptée aux méthodes de l'analyse combinatoire[8].

Elle se réécrit : ${\displaystyle c_{n}=\sum _{k=1}^{n}B_{n,k}(a_{1},\dots ,a_{n-k+1})b_{k},}$ où les B_n,k(a₁,...,a_n−k+1) sont des polynômes de Bell.

Par exemple, pour tout entier ℓ ≥ 1 : ${\displaystyle \left(\sum _{n=1}^{\infty }{a_{n} \over n!}X^{n}\right)^{\ell }=\sum _{n=\ell }^{\infty }B_{n,\ell }(a_{1},\dots ,a_{n-\ell +1}){\ell$ ! \over n!}X^{n},} ou encore : ${\displaystyle \left(\sum _{m=0}^{\infty }d_{m}X^{m}\right)^{\ell }=\sum _{k=0}^{\infty }B_{k+\ell ,\ell }\left({\tfrac {d_{0}}{1!}},\dots ,{\tfrac {d_{k}}{(k+1)!}}\right){\ell$ ! \over (k+\ell )!}X^{k}.}