APPROXIMATION NUMÉRIQUE DE RACINES ISOLÉES MULTIPLES DE SYSTÈMES ANALYTIQUES
Résumé
Abstract. The approximation of a multiple isolated root is a difficult problem. In fact the root can even be a repulsive root for a fixed point method like the Newton method. However there exists a huge literature on this topic but the answers given are not satisfactory. Numerical methods allowing a local convergence analysis work often under specific hypotheses. This viewpoint favouring numerical analysis forgets the geometry and the structure of the local algebra. Thus appeared so-called symbolic-numeric methods, yet full of lessons, but their precise numerical analysis is still missing. We propose in this paper a method of symbolic-numeric kind, whose numerical treatment is certified. The general idea is to construct a finite sequence of systems, admitting the same root, and called the deflation sequence, so that the multiplicity of the root drops strictly between two successive systems. So the root becomes regular. Then we can extract a regular square system we call deflated system. We described already the construction of this deflated sequence when the singular root is known. The originality of this paper consists on one hand to construct a deflation sequence from a point close to the root and on the other hand to give a numerical analysis of this method. Analytic square integrable functions build the functional frame. Using the Bergman kernel, reproducing kernel of this functional frame, we are able to give a α-theory à la Smale. Furthermore we present new results on the determination of the numerical rank of a matrix and the closeness to zero of the evaluation map. As an important consequence we give an algorithm computing a deflation sequence free of ε, threshold quantity measuring the numerical approximation, meaning that the entry of this algorithm does not involve the variable ε.
L'approximation d'une racine isolée multiple est un problème difficile. En effet la racine peut même être répulsive pour une méthode de point fixe comme la méthode de Newton. La littérature sur le sujet est vaste mais les réponses proposées pour résoudre ce problème ne sont pas satisfaisantes. Des méthodes numériques qui permettent de faire une analyse locale de convergence sont souvent élaborées sous des hypothèses particulières. Ce point de vue privilégiant l'analyse numérique néglige la géométrie et la structure de l'algèbre locale. C'est ainsi qu'ont émergé des méthodes qualifiés de symboliques-numériques. Mais l'analyse numérique précise de ces méthodes pourtant riches d'enseignement n'a pas été faite. Nous proposons dans cet article une méthode de type symbolique-numérique dont le traitement numérique est certifié. L'idée générale est de construire une suite finie de systèmes admettant la même racine, appelée suite de déflation, telle que la multiplicité de la racine chute strictement entre deux systèmes successifs. La racine devient ainsi régulière lors du dernier système. Il suffit alors d'en extraire un système carré régulier pour obtenir ce que nous appelons système déflaté. Nous avions déjà décrit la construction de cette suite de déflation quand la racine est connue. L'originalité de cette étude consiste d'une part à définir une suite de déflation à partir d'un point proche de la racine et d'autre part à donner une analyse numérique de cette méthode. Le cadre fonctionnel de cette analyse est celui des systèmes analytiques constitués de fonctions de carré intégrable. En utilisant le noyau Bergman, noyau reproduisant de cet espace fonctionnel, nous donnons une α-théorie à la Smale de cette suite de déflation. De plus nous présentons des résultats nouveaux relatifs à la détermination du rang numérique d'une matrice et à celle de la proximité à zéro de l'application évaluation. Comme conséquence importante nous donnons un algorithme de calcul d'une suite de déflation qui est libre de ε, quantité-seuil qui mesure l'approximation numérique, dans le sens que les entrées de cet algorithme ne comportent pas la variable ε.
Classification mathématique par sujets (2010). 65F30, 65H10, 65Y20, 68Q25, 68W30.
Origine : Fichiers produits par l'(les) auteur(s)
Loading...