Recherche:Formule explicite des nombres de Bernoulli

Une page de Wikiversité.
Aller à : Navigation, rechercher
Nuvola apps edu mathematics-p.svg
Wikimedia-research-fr.png
Travail de recherche : Formule explicite des nombres de Bernoulli

Cette page dépend de l’espace Travaux de recherche et projets collaboratifs de Wikiversité. En particulier, elle ne fait pas partie de l’espace des cours proposés par Wikiversité. Consultez cette page pour plus de détails.



Introduction

Les nombres de Bernoulli sont parmi les objets le plus fascinants des mathématiques. On les retrouve en arithmétique, en théorie des nombres, en analyse et même en topologie.

De telles expressions sont toujours des polynômes en m, de degré n+1\, et sont appelées polynômes de Bernoulli. Les coefficients des polynômes de Bernoulli sont liés aux nombres de Bernoulli de la façon suivante :

\sum_{k=0}^{m-1} k^n = {1\over{n+1}}\sum_{k=0}^n{n+1\choose{k}} B_k m^{n+1-k}

Par exemple, en donnant à n la valeur 1, on obtient :

0+1+2+\ldots+(m-1) = \frac{1}{2}(B_0 m^2+2 B_1 m^1) = \frac{1}{2}(m^2-m)\,

Les nombres de Bernoulli ont d'abord été étudiés par Jacques Bernoulli, ce qui a conduit Abraham de Moivre à leur donner le nom que nous connaissons aujourd'hui.

Il est possible de calculer les nombres de Bernoulli en utilisant la formule de récurrence suivante :

\sum_{k=0}^n{n+1\choose{k}}B_k = 0\,

avec la condition initiale : B_0 = 1\,.

Rectification

La démonstration suivante a pour but d'expliciter les nombres de Bernoulli et de voir une application utile a travers la série de Taylor de quelques formules trigonométriques.

Suite à un malentendu et a des recherches plus approfondies, il s'est avéré que la démonstration existe à cette page. La démonstration a été établie bien avant par Maurice D'Ocagne: ici

Bonne lecture!

Définitions

Les coefficients K_n^p désignent les nombres de Sterling de seconde espece.

Soit le polynôme x_{{[k]}}\,=x \left( x-1 \right)  \left( x-2 \right) ...\left( x-k+1 \right) avec k≥1, appelé kieme puissance descendante de x.

On définit aussi la kième puissance montante ou factorielle de x en posant: x^{{[k]}}\,=x \left( x+1 \right)  \left( x+2 \right) ...\left( x+k-1 \right).

Les puissances montantes et descendantes sont liées par l'identité: (-x)^{{[k]}}\,=\left( -1 \right) ^{k}x_{{[k]}}.

1er lemme

Pour toute fonction f sur \R, on pose: \Delta f \left( x \right)=f \left( x+1 \right) -f \left( x \right). On a: \Delta\,x_{{[k]}}=k.x_{{[k-1]}}.

En effet: (x+1)_{{[k]}}-x_{{[k]}}=\left( x+1 \right) x_{{[k-1]}}- \left( x-k+1 \right) x_{{[k-1]}}=k.x_{{[k-1]}}.

On a donc, pour k fixé: \sum_{x=0}^{n+1} {(x+1)}^{{[k]}} = \sum_{x=0}^{n+1} {\frac {{\Delta x}^{[k+1]}}{k+1}} . On en déduit par somme télescopique:

1^{{[k]}}+2^{{[k]}}+...+n^{{[k]}} = {\frac {{n}^{[k+1]}}{k+1}}.

Exemple:

Pour k=2: 1*2+2*3+3*4+...+n(n+1) = \frac{n(n+1)(n+2)}{3}

2ème lemme

À partir de là, on arrive à exprimer les puissances ordinaires en fonction des puissances montantes. Exemples:

x^1=x^{[1]}\,

x^2=x^{[2]}-x^{[1]}\,

x^3=x^{[3]}-3x^{[2]}+x^{[1]}\,

x^4=x^{[4]}-6x^{[3]}+7x^{[2]}-x^{[1]}\,

On obtient alors un triangle arithmétique. Il existe donc des coefficients uniques notés K_n^p tel que: x^n=\sum_{p=0}^{n} (-1)^{n-p}.K_n^p.x^{[p]} et (-x)^n=\sum_{p=0}^{n} K_n^p.x_{[p]}

3ème lemme

Soient B_n\, les nombres de Bernoulli: B_n = \sum_{i=1}^{n+1} (-1)^{i-1}.K_n^{i-1}.\frac {(i-1)!}{i}

Demonstration:

x^p = \sum_{i=1}^{p+1} K_p^{i-1}. x_{[i-1]} signifie: \sum_{x=0}^{n-1} x^p = \sum_{x=0}^{n-1} \left ( \sum_{i=1}^{p+1} K_p^{i-1}.x_{[i-1]} \right ) = \sum_{i=1}^{p+1} K_p^{i-1}.\sum_{x=0}^{n-1} x_{[i-1]}

La relation \Delta\,x_{{[i]}}=i*x_{{[i-1]}} entraine par somme telescopique: \sum_{x=0}^{n-1} x_{[i-1]} = \frac{n_{[i]}-0_{[i]}}{i},car i>0. Par suite \sum_{x=0}^{n-1} x^p = \sum_{i=1}^{p+1} K_p^{i-1}.\frac{n_{[i]}}{i}

Revenons maintenant aux polynômes de Bernoulli: Il a été déjà établi que: \frac {x.e^{tx}} {e^x-1} = \sum_{n=0}^{\infty} \frac{B_n(t)}{n!} .x^n Calculons maintenant: B_{{n}} \left( t+1 \right) -B_{{n}} \left( t \right):

\sum_{n=0}^{\infty} \frac {B_n(t+1)-B_n(t)}{n!} .x^n = \frac {x.e^{(t+1)x}}{e^x-1} - \frac {x.e^{tx}} {e^x-1} = x.e^{tx} = \sum_{n=0}^{\infty} \frac{n.t^{n-1}}{n!} .x^n

Par identification: B_{{n}} \left( t+1 \right) -B_{{n}} \left( t \right) = n.t^{n-1}, 	\forall  n\ge0 Donc:

\sum_{k=0}^{n-1} k^p = \frac{1}{p+1}\left( B_{{p+1}} \left( n \right) -B_{{p+1}} \left( 0 \right) \right)

Une identité connue de ces polynômes est: B_{p+1}^' = (p+1).B_p(x) Donc:

\int _{0}^{n} B_{p+1}^'(x) dx = (p+1).\int _{0}^{n} B_{p}(x) dx \Rightarrow \sum_{k=0}^{n-1} k^p = \frac{1}{p+1}\left( B_{{p+1}} \left( n \right) -B_{{p+1}} \left( 0 \right) \right) = \int _{0}^{n} B_{p}(x) dx

\Rightarrow \int _{0}^{n} B_{p}(x) dx = \sum_{i=1}^{p+1} K_p^{i-1}.\frac{n_{[i]}}{i}

En dérivant par rapport a x on obtient: B_p(x) = \sum_{i=1}^{p+1} K_p^{i-1}.\frac{d}{dx}\frac{x_{[i]}}{i} Donc: B_n(0) = \sum_{i=1}^{n+1} K_n^{i-1}.\frac{d}{dx}\frac{x_{[i]}}{i} \Bigg|_{x=0}

Le polynôme x_{[i]} s'annule en x=0, puisque i>0. La valeur de sa dérivée en x=0 est donnée parProxy-Connection: keep-alive Cache-Control: max-age=0

\frac{x_{[i]}}{i} \Bigg|_{x=0} = (-1).(-2)...(-i+1) = (-1)^{(i-1)}.(i-1)! Finalement:

B_n = \sum_{i=1}^{n+1} (-1)^{i-1}.K_n^{i-1}.\frac {(i-1)!}{i}

4ème lemme:

Soit l'opérateur de dérivation: \frac{d}{dx}.x(f) = \frac{d}{dx}(xf) pour f de  C^{\infty}( \R ), alors \forall n: \left( \frac{d}{dx}.x \right)^{(n)}(f) = \sum_{p=0}^{n} K_n^p.x^p.f^{(p)}

Exemples:

\left( \frac{d}{dx}.x \right)^{(2)}(f) = \left( \frac{d}{dx}.x \right) \left( \frac{d}{dx}.x \right) (f) = \left( \frac{d}{dx}.x \right){\left( xf^' \right)} =\ xf^'+x^2f^{''}

En continuant, on obtient:

\left( \frac{d}{dx}.x \right)^{(1)}(f) = xf^{(1)}

\left( \frac{d}{dx}.x \right)^{(2)}(f) = xf^{(1)}+x^2f^{(2)}

\left( \frac{d}{dx}.x \right)^{(3)}(f) = xf^{(1)}+3x^2f^{(2)}+x^3f^{(3)}

\left( \frac{d}{dx}.x \right)^{(4)}(f) = xf^{(1)}+7x^2f^{(2)}+6x^3f^{(3)}+x^4f^{(4)}

On remarque que l'on obtient les mêmes coefficients du triangle arithmétique précédent, d'où:

\left( \frac{d}{dx}.x \right)^{(n)}(f) = \sum_{p=0}^{n} K_n^p.x^p.f^{(p)}

Théorème:

Soient les coefficients K_n^p définis pour n \in \N et p \le n:

K_{n+1}^p = pK_n^p + K_n^{p-1} avec: K_n^n=1 et K_n^0=0, \forall n \in \N

Alors: K_n^p = \frac {1}{p!}.\sum _{i=0}^{p} (-1)^{p-1}.C_p^i.{i}^{n}

En effet, si l'on applique le 4ème lemme à la fonction f(x) = e^x, on obtient: \sum _{p=0}^{\infty }{\frac {{p}^{n}.{x}^{p}}{p!}} = e^x.\sum _{p=0}^{n} K_n^p.x^p

Donc, pour n fixé, on pose a_p = p!.K_n^p, pour tout k \ge 0. Avec cette notation l'égalité devient:

\sum _{k=0}^{\infty }{\frac {a_{{k}}{x}^{k}}{k!}} = {e^{-x}}\sum _{k=0}^{\infty }{\frac {{k}^{n}{x}^{k}}{k!}} \Rightarrow a_p = \sum _{i=0}^{p}{C_{{p}}}^{i} \left( -1 \right) ^{p-i}{i}^{n} \Rightarrow K_n^p = \frac {1}{p!}.\sum _{i=0}^{p} (-1)^{p-1}.C_p^i.{i}^{n}

Corollaire

\forall n, B_n = \sum _{i=0}^{n} \left( \frac{1}{i+1}.\sum_{p=0}^{i} \left( (-1)^p.C_i^p.p^n \right) \right)

Pour l'obtenir, il suffit d'injecter la formule explicite des nombres de Sterling dans celle des nombres de Bernoulli (lemme3).

Applications

Les nombres de Bernoulli sont présents dans le développement limité de la fonction tangente et tangente hyperbolique en série de Taylor(c'est-à-dire sous sa forme générale).

On pose: B(x) = \frac{x}{e^x-1} (génératrice exponentielle). On a: B(x) + \frac{x}{2} = \frac{x}{2}\frac{e^x+1}{e^x-1}.

Commençons par: tanh(x) = {\frac {{e^{x}}-{e^{-x}}}{{e^{x}}+{e^{-x}}}} = {\frac {{e^{2\,x}}-1}{{e^{2\,x}}+1}} = 2{\frac {{e}^{4\,x}+1}{{e}^{4\,x}-1}}-{\frac {{e}^{2\,x}+1}{{e}^{2\,x}-1}} = \left( B(4x)+2x \right) - \left( B(2x)+x \right) = B(4x)-B(2x)+x

Or, on sait que: B(x) = \frac{x}{e^x-1} = \sum_{n=0}^{\infty} B_n.\frac{x^n}{n!}

En remplaçant on obtient: tanh(x) = \sum _{n=1}^{\infty }{\frac {B_{{2\,n}}.{4}^{n} \left( {4}^{n}-1 \right) }{ \left( 2\,n \right) !}{x}^{2\,n-1}}

La tangente s'obtient en remplaçant x par ix dans la nouvelle expression de tanh(x), avec: i2=(-1). On obtient:

tan(x) = \sum _{n=1}^{\infty }{\frac {|B_{{2\,n}}|.{4}^{n} \left( {4}^{n}-1 \right) }{ \left( 2\,n \right) !}{x}^{2\,n-1}}