Méthodes multipas linéaires

Les méthodes multipas linéaires calculent chaque nouvelle valeur de solution à partir d'une combinaison linéaire de plusieurs valeurs de solution et de dérivées précédentes, réutilisant le travail passé pour atteindre un ordre élevé à faible coût par pas.

Definition

Une méthode multipas linéaire est une méthode pour les équations différentielles ordinaires qui détermine la prochaine valeur de solution par une relation linéaire fixe entre un certain nombre de valeurs de solution précédentes et d'évaluations du second membre.

Scope

Ce sujet couvre les familles d'Adams-Bashforth (explicites) et d'Adams-Moulton (implicites), les formules de différenciation rétrograde pour les problèmes raides, l'implémentation des prédicteur-correcteur, les polynômes caractéristiques et la condition de racine qui définissent la zéro-stabilité, ainsi que les barrières d'ordre de Dahlquist qui limitent ce que de telles méthodes peuvent accomplir.

Core questions

• Comment les méthodes multipas réutilisent-elles les valeurs passées pour atteindre un ordre élevé avec une seule nouvelle évaluation de fonction par pas ?
• Qu'est-ce que la zéro-stabilité, et comment la condition de racine sur le polynôme caractéristique l'exprime-t-elle ?
• Comment les paires prédicteur-correcteur combinent-elles les formules explicites et implicites en pratique ?
• Que disent les barrières d'ordre de Dahlquist sur les limites de la précision et de la stabilité des méthodes multipas ?

Key theories

Zéro-stabilité et la condition de racine

Une méthode multipas est zéro-stable, et donc convergente lorsqu'elle est consistante, précisément lorsque les racines de son premier polynôme caractéristique se trouvent dans le disque unité fermé avec seulement des racines simples sur la frontière ; cette condition de racine est l'analogue multipas de la stabilité.

Barrières de Dahlquist

La première barrière de Dahlquist limite l'ordre d'une méthode à k pas zéro-stable, et sa seconde barrière montre qu'aucune méthode multipas linéaire A-stable ne peut avoir un ordre supérieur à deux, c'est pourquoi les solveurs de problèmes raides d'ordre élevé s'appuient sur le compromis BDF de stabilité relative plutôt qu'absolue.

Mechanisms

Les méthodes d'Adams intègrent un polynôme d'interpolation à travers les valeurs de dérivées passées : Adams-Bashforth utilise uniquement des valeurs connues (explicites), Adams-Moulton inclut la nouvelle valeur inconnue (implicite) pour une plus grande précision et stabilité. En pratique, les deux sont associées comme un prédicteur-correcteur : la formule explicite prédit, l'implicite corrige, généralement en une ou deux itérations. Les formules de différenciation rétrograde, quant à elles, différencient les valeurs de solution passées pour approximer la dérivée au nouveau point, donnant les méthodes stables pour les problèmes raides qui sont au cœur des codes d'EDO raides. Étant donné que les méthodes multipas nécessitent plusieurs valeurs de départ, elles sont amorcées par une méthode à un pas.

Clinical relevance

Les méthodes multipas linéaires, en particulier les formules de différenciation rétrograde, sont à la base des solveurs d'EDO raides de production utilisés en cinétique chimique, en simulation de circuits électroniques et dans les grands systèmes différentiels-algébriques, où l'évaluation du second membre est coûteuse et la réutilisation des évaluations passées via des formules multipas génère des gains d'efficacité majeurs.

History

Adams et Bashforth ont introduit les formules multipas au XIXe siècle, Moulton ajoutant des variantes implicites ; l'analyse de Dahlquist dans les années 1950-1960 a établi la théorie de la stabilité et les barrières d'ordre qui régissent le domaine, et les travaux de C. William Gear dans les années 1970 ont fait des codes basés sur les formules de différenciation rétrograde la norme pour les problèmes raides.

Key figures

• John Couch Adams
• Francis Bashforth
• Forest Ray Moulton
• Germund Dahlquist
• C. William Gear

Seminal works

• hairer1993
• iserles2008

Frequently asked questions

En quoi les méthodes multipas diffèrent-elles des méthodes de Runge-Kutta ?

Les méthodes de Runge-Kutta effectuent plusieurs nouvelles évaluations de dérivées à l'intérieur de chaque pas mais les écartent ensuite, tandis que les méthodes multipas réutilisent les valeurs de dérivées des pas précédents. Les méthodes multipas sont donc moins coûteuses par pas mais nécessitent des valeurs de départ supplémentaires et une gestion spéciale des changements de taille de pas.

Qu'est-ce que la condition de racine ?

C'est l'exigence selon laquelle les racines du premier polynôme caractéristique de la méthode se trouvent à l'intérieur ou sur le cercle unité, avec des racines simples sur la frontière. Elle garantit que les petites erreurs ne sont pas amplifiées à mesure que les pas s'accumulent, assurant ainsi que la méthode est zéro-stable et donc convergente.