Méthode de Frobenius

En analyse, la méthode de Frobenius, du nom du mathématicien allemand Ferdinand Georg Frobenius, est une technique d'obtention du développement en série entière des solutions d'une équation différentielle linéaire de la forme :
 

la variable  étant en général complexe, au voisinage du point , sous réserve que  et  soient analytiques, ou possèdent un point singulier dit régulier en ce point. Si ces conditions sont respectées, la méthode de Frobenius permet alors de déterminer au moins une solution de la forme :

 

Cette méthode se généralise à une équation différentielle linéaire d'ordre p quelconque, sous réserve des conditions de régularité suffisantes sur les fonctions  apparaissant devant chacune des dérivées .

Principes généraux 
Il est fréquent qu'il ne soit pas possible d'intégrer directement les équations différentielles, même linéaires, et d'en exprimer les solutions à partir de polynômes ou de fonctions transcendantales "ordinaires" (par exemple exponentielle, logarithme, fonctions trigonométriques…). Par ailleurs, même si une solution analytique peut être obtenue, celle-ci peut avoir une forme très complexe, et peu utilisable en pratique.

Pour toutes ces raisons, il est utile de disposer de méthodes permettant d'obtenir des formes approchées des solutions d'une équation différentielle. Ces méthodes peuvent être classées en deux grandes catégories :
 des méthodes locales, où l'on cherche à exprimer la solution approchée au voisinage d'un point particulier ;
 des méthodes globales, où l'on cherche à exprimer cette solution approchée sur un certain intervalle : il s'agit essentiellement des méthodes de perturbation, où l'on simplifie l'équation initiale en négligeant certains termes, ce qui rend l'équation facilement intégrable, avant de tenir compte des termes négligés en utilisant la solution « d'ordre 0 » ainsi obtenue. Ce type de méthode est fréquemment utilisée en astronomie ou en physique quantique.

Principe de la méthode 
La méthode de Frobenius appartient à la première catégorie de méthodes. Comme indiqué en introduction, elle s'applique aux équations différentielles linéaires du second ordre de la forme générale :
 

et permet d'obtenir au voisinage d'un point  un développement d'au moins une solution en série entière de la forme 

Le principe sera de substituer cette expression dans l'équation différentielle, et par identification terme à terme de déterminer  puis une expression des coefficients , qui fait généralement intervenir une relation de récurrence liant entre ceux-ci.

Existence et convergence de la série 
Si du fait de la linéarité de l'équation différentielle il est toujours possible, formellement, d'effectuer cette substitution, la question de la convergence de la série, et du rayon associé, se pose. Il est utile de se placer pour cette question dans le plan complexe. 

En effet, une fonction  de variable complexe (et uniforme) dérivable au voisinage d'un point  est en fait indéfiniment dérivable dans ce voisinage. Elle est alors dite analytique autour de ce point, et pour tout point contenu dans le domaine autour de  où  est analytique, elle est développable en série entière, avec une convergence uniforme dans tout le domaine considéré.

Points réguliers et singuliers 
De façon générale, pour un point  une fonction de variable complexe, analytique dans un domaine  autour de ce point, est développable en série de Laurent : 

ce développement étant unique, et la convergence uniforme dans tout le domaine considéré. La série  est la partie régulière du développement, la série  étant la partie principale de celui-ci.

Si  est analytique en , tous les coefficients de la partie principale de son développement de Laurent sont nuls : ce type de point est dit régulier. Dans le cas contraire, le point est dit singulier, et il est possible de distinguer : 
 les points singuliers essentiels, pour lesquels la partie principale du développement comprend une infinité de termes : dans ce cas,  n'existe pas ;
 les pôles d'ordre p, pour lesquels la partie principale comprend un nombre fini de termes et est de la forme : .

Forme générale des solutions - Théorème de Fuchs 
Les notions précédentes permettent de préciser les conditions d'existence d'au moins une solution de l'équation différentielle admettant, dans un domaine de rayon non nul autour du point considéré, un développement en série entière de la forme proposée plus haut. 

De fait, il est possible de démontrer alors le théorème suivant relatif à l'existence d'une telle solution :
Si, dans l'équation différentielle de la forme , les fonctions  et  sont analytiques dans le domaine , il existe alors une solution de l'équation de la forme générale , avec  constant, et  analytique sur tout le domaine considéré. 

Si alors  et  sont en plus analytiques sur le domaine  incluant , ce point sera dit ordinaire. Si ce n'est pas le cas, le point sera qualifié de singulier. Si cependant il existe, pour toute solution de l'équation différentielle, un réel  tel que , le point singulier est dit régulier. Dans le cas contraire, il sera irrégulier.

Le  précise alors que  sera un point singulier régulier si  et  sont analytiques dans tout le domaine . 

Il est aussi possible d'étudier le caractère régulier ou singulier à l'infini, en effectuant dans l'équation différentielle le changement de variable (transformation de Möbius) , et d'étudier le comportement en . 

  soit l'équation différentielle . Clairement  et  sont analytiques sur tout le plan complexe sauf en l'origine, dont il est évident qu'elle constitue un pôle d'ordre 1. Comme en revanche  et  sont analytiques à l'origine, le point  constitue un point singulier régulier. 
 Si l'on effectue alors le changement de variable  l'équation différentielle se met sous la forme , donc dans ce cas il est clair que  et  ne sont pas analytiques en  : on en déduit donc que  est un point singulier irrégulier pour cette équation différentielle.
Par suite, si en  les fonctions  et  sont analytiques, ou sont telles que ce point soit un point singulier régulier, alors il existera au moins une solution de la forme générale , avec  réel et  analytique dans un voisinage . Il est clair d'après ce qui précède que ce rayon de convergence  sera égal au minimum à la distance au point singulier le plus proche de . Tout ceci justifie la forme du développement proposé dans la méthode de Frobenius.

Exposé de la méthode 
La méthode de Frobenius comporte les étapes suivantes :
 Soit l'équation différentielle  de la forme générale  et un point ordinaire, ou singulier régulier en , on pose  et , lesquelles fonctions sont analytiques au point considéré ;
 On substitue dans l'équation différentielle l'expression de la solution sous la forme de la série entière (dite de Frobenius):  avec , appelé exposant indiciel. 
 Il vient alors successivement les expressions suivantes pour les différents termes de l'équation différentielle :
 
 
 
 Par suite l'équation différentielle se met sous la forme :
 
 or par identification cette égalité n'est vérifiée que si chacun des termes de même degré en  de la somme est nul.
 Par convention, on suppose que  et l'on considère le terme d'ordre le plus bas, soit celui en , évalué en  où par définition  et  sont analytiques, il vient une équation du second degré en  appelée équation indicielle :
 
 dont les racines sont notées de façon générale  et  avec .
 La substitution dans les termes suivants d'une des valeurs de  ainsi obtenues permet d'obtenir une expression des coefficients ,  restant arbitraire (le plus souvent on prend ). Généralement cette expression prend la forme d'une relation de récurrence entre ,  et .

Exemples d'application 
Deux exemples permettent d'illustrer la méthode :
 Soit l'équation différentielle , qui se met aussitôt sous la forme 

 Par suite  et , ces fonctions admettant bien un point singulier régulier en , et , .
 D'après ce qui précède, la substitution dans l'équation du développement en série de Frobenius pris autour de  donne :
 ,
 soit encore, après réarrangement des termes de même degré  :
,
 l'équation indicielle étant, elle, donnée par , laquelle admet la racine double .
 La substitution de cette valeur de l'exposant indiciel dans l'expression précédente donne alors :
 
 ce qui implique la relation de récurrence , soit en fonction de  l'expression .
 Soit l'équation différentielle : , pour laquelle  est un point singulier régulier de façon évidente. Dans ce cas, ,  et l'équation indicielle est donnée par , laquelle admet deux racines distinctes  et .
 La substitution du développement en série entière de  donne, elle, l'expression :
  ce qui implique que  pour  pour toutes les valeurs de .
Il n'y a donc pas de relation de récurrence, toutefois en prenant  on obtient cependant deux solutions possibles  et . Il est facile de vérifier qu'il s'agit bien de solutions de l'équation, toutefois la seconde n'est pas analytique en l'origine (pôle d'ordre 2).

Notes et références

Notes

Références 

Série (mathématiques)
Analyse complexe
Équation différentielle
Physique mathématique
Méthode mathématique de la physique