Approximation numérique de racines isolées multiples de systèmes analytiques
Annales Henri Lebesgue, Volume 3 (2020), pp. 901-957.

Metadata

Keywords systèmes d’équations, racines singulières, déflation, rang numérique, évaluation

Abstract

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 de 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 ε.


References

[AH11] Argyros, Ionnis K.; Hilout, Said Semilocal convergence of Newton’s method for singular systems with constant rank derivatives, J. Korean Soc. Math. Educ., Ser. B, Pure Appl. Math., Volume 18 (2011) no. 2, pp. 97-111 | MR | Zbl

[BCSS97] Blum, Leonore; Cucker, Felipe; Shub, Michael; Smale, Steve Complexity and Real Computation, Springer, 1997 (foreword by Richard M. Karp) | Zbl

[CLO05] Cox, David A.; Little, John; O’shea, Donald B. Using algebraic geometry, Graduate Texts in Mathematics, Volume 185, Springer, 2005 | MR | Zbl

[Ded06] Dedieu, Jean-Pierre Points fixes, zéros et la méthode de Newton, Mathématiques & Applications, Volume 54, Springer, 2006 | Zbl

[DK80a] Decker, Dwight W.; Kelley, Carl T. Newton’s method at singular points. I., SIAM J. Numer. Anal., Volume 17 (1980) no. 1, pp. 66-70 | DOI | MR | Zbl

[DK80b] Decker, Dwight W.; Kelley, Carl T. Newton’s method at singular points. II, SIAM J. Numer. Anal., Volume 17 (1980) no. 3, pp. 465-471 | DOI | MR | Zbl

[DKK83] Decker, Dwight W.; Keller, Herbert B.; Kelley, Carl T. Convergence rates for Newton’s method at singular points, SIAM J. Numer. Anal., Volume 20 (1983) no. 2, pp. 296-314 | DOI | MR | Zbl

[DLZ11] Dayton, Barry H.; Li, Tien-Yien; Zeng, Zhonggang Multiple zeros of nonlinear systems, Math. Comput., Volume 80 (2011) no. 276, pp. 2143-2168 | DOI | MR | Zbl

[DS01] Dedieu, Jean-Pierre; Shub, Mike On simple double zeros and badly conditioned zeros of analytic functions of n variables, Math. Comput., Volume 70 (2001) no. 233, pp. 319-327 | DOI | MR | Zbl

[DZ05] Dayton, Barry H.; Zeng, Zhonggang Computing the multiplicity structure in solving polynomial systems, Proceedings of the 2005 international symposium on Symbolic and algebraic computation (2005), pp. 116-123 | Zbl

[EI95] Eisenstat, Stanley C.; Ipsen, Ilse C. F. Relative perturbation techniques for singular value problems, SIAM J. Numer. Anal., Volume 32 (1995) no. 6, pp. 1972-1988 | DOI | MR | Zbl

[Ems78] Emsalem, Jacques Géométrie des points épais, Bull. Soc. Math. Fr., Volume 106 (1978), pp. 399-416 | DOI | Numdam | MR | Zbl

[EY36] Eckart, Carl; Young, Gale The approximation of one matrix by another of lower rank, Psychometrika, Volume 1 (1936), pp. 211-218 | DOI | Zbl

[FH05] Fierro, Ricardo D.; Hansen, Per Christian UTV expansion pack : Special-purpose rank-revealing algorithms, Numer. Algorithms, Volume 40 (2005) no. 1, pp. 47-66 | DOI | MR | Zbl

[GLSY05] Giusti, Marc; Lecerf, Grégoire; Salvy, Bruno; Yakoubsohn, Jean-Claude On location and approximation of clusters of zeros of analytic functions, Found. Comput. Math., Volume 5 (2005) no. 3, pp. 257-311 | DOI | MR

[GLSY07] Giusti, Marc; Lecerf, Grégoire; Salvy, Bruno; Yakoubsohn, Jean-Claude On location and approximation of clusters of zeros : Case of embedding dimension one, Found. Comput. Math., Volume 7 (2007) no. 1, pp. 1-58 | DOI | MR | Zbl

[GO81] Griewank, Andreas; Osborne, Michael R. Newton’s method for singular problems when the dimension of the null space is >1, SIAM J. Numer. Anal., Volume 18 (1981) no. 1, pp. 145-149 | DOI | MR | Zbl

[GO83] Griewank, Andreas; Osborne, Michael R. Analysis of Newton’s method at irregular singularities, SIAM J. Numer. Anal., Volume 20 (1983) no. 4, pp. 747-773 | DOI | MR | Zbl

[Gri85] Griewank, Andreas n solving nonlinear equations with simple singularities or nearly singular solutions, SIAM Rev., Volume 27 (1985) no. 4, pp. 537-563 | DOI | MR | Zbl

[GY14] Giusti, Marc; Yakoubsohn, Jean-Claude Multiplicity hunting and approximating multiple roots of polynomial systems, Recent Advances in Real Complexity and Computation (Contemporary Mathematics) Volume 604, American Mathematical Society, 2014, pp. 105-128 | DOI

[HMS17] Hauenstein, Jonathan D.; Mourrain, Bernard; Szanto, Agnes On deflation and multiplicity structure, J. Symb. Comput., Volume 83 (2017), pp. 228-253 | DOI | MR | Zbl

[Kat95] Kato, Tosio Perturbation theory for linear operators, Classics in Mathematics, Volume 132, Springer, 1995 | Zbl

[Kra13] Krantz, Steven G. Geometric analysis of the Bergman kernel and metric, Graduate Texts in Mathematics, Volume 268, Springer, 2013 | MR | Zbl

[KS83] Kelley, Carl T.; Suresh, Ram A new acceleration method for Newton’s method at singular points, SIAM J. Numer. Anal., Volume 20 (1983) no. 5, pp. 1001-1009 | DOI | MR | Zbl

[Lec02] Lecerf, Grégoire Quadratic Newton iteration for systems with multiplicity, Found. Comput. Math., Volume 2 (2002) no. 3, pp. 247-293 | DOI | MR | Zbl

[LVZ06] Leykin, Anton; Verschelde, Jan; Zhao, Ailing Newton’s method with deflation for isolated singularities of polynomial systems, Theor. Comput. Sci., Volume 359 (2006) no. 1-3, pp. 111-122 | DOI | MR | Zbl

[LZ12a] Li, Nan; Zhi, Lihong Computing isolated singular solutions of polynomial systems : case of breadth one, SIAM J. Numer. Anal., Volume 50 (2012) no. 1, pp. 354-372 | MR | Zbl

[LZ12b] Li, Nan; Zhi, Lihong Computing the multiplicity structure of an isolated singular solution : Case of breadth one, J. Symbolic Comput., Volume 47 (2012) no. 6, pp. 700-710 | MR | Zbl

[LZ14] Li, Nan; Zhi, Lihong Verified error bounds for isolated singular solutions of polynomial systems, SIAM J. Numer. Anal., Volume 52 (2014) no. 4, pp. 1623-1640 | MR | Zbl

[Mar11] Markovsky, Ivan Low rank approximation : algorithms, implementation, applications, Communications and Control Engineering Series, Springer, 2011

[MG03] Miranian, L.; Gu, Ming Strong rank revealing LU factorizations, Linear Algebra Appl., Volume 367 (2003) no. 1, pp. 1-16 | DOI | MR | Zbl

[Mir60] Mirsky, Leon Symmetric gauge functions and unitarily invariant norms, Q. J. Math., Oxf. II. Ser., Volume 11 (1960), pp. 50-59 | DOI | MR | Zbl

[MM11] Mantzaflaris, Angelos; Mourrain, Bernard Deflation and certified isolation of singular zeros of polynomial systems, ISSAC 2011—Proceedings of the 36th international symposium on Symbolic and Algebraic Computation (2011), pp. 249-256 | Zbl

[Mou97] Mourrain, Bernard Isolated points, duality and residues, J. Pure Appl. Algebra, Volume 117/118 (1997), pp. 469-493 | DOI | MR | Zbl

[Oji87] Ojika, Takeo Modified deflation algorithm for the solution of singular problems. I. A system of nonlinear algebraic equations, J. Math. Anal. Appl., Volume 123 (1987) no. 1, pp. 199-221 | DOI | MR | Zbl

[OWM83] Ojika, Takeo; Watanabe, Satoshi; Mitsui, Taketommo Deflation algorithm for the multiple roots of a system of nonlinear equations, J. Math. Anal. Appl., Volume 96 (1983) no. 2, pp. 463-479 | DOI | MR | Zbl

[Pan00] Pan, Chin-Tsuan On the existence and computation of rank-revealing LU factorizations, Conference Celebrating the 60th Birthday of Robert J. Plemmons (Winston–Salem, NC, 1999) (Linear Algebra and its Applications) Volume 316 (2000) no. 1-3, pp. 199-222 | MR | Zbl

[Ral66] Rall, Louis B. Convergence of the Newton process to multiple solutions, Numer. Math., Volume 9 (1966) no. 1, pp. 23-37 | DOI | MR | Zbl

[Red78] Reddien, Georges W. On Newton’s method for singular problems, SIAM J. Numer. Anal., Volume 15 (1978) no. 5, pp. 993-996 | DOI | MR | Zbl

[Red79] Reddien, Georges W. Newton’s method and high order singularities, Comput. Math. Appl., Volume 5 (1979) no. 2, pp. 79-86 | DOI | MR | Zbl

[Rud08] Rudin, Walter Function theory in the unit ball of C n , Classics in Mathematics, Springer, 2008 (Reprint of the 1980 original) | Zbl

[Sch70] Schröder, Ernst Ueber unendliche viele Algorithmen zur Auflösung der Gleichungen, Math. Ann., Volume 2 (1870), pp. 317-365 | DOI

[Sha67] Shamanskiĭ, V. On the application of Newton’s method in a singular case, USSR Computational Mathematics and Mathematical Physics, Volume 7 (1967) no. 4, pp. 72-85 | DOI | MR

[SS90] Stewart, Gilbert W.; Sun, Ji-Guang Matrix Perturbation Theory, Computer Science and Scientific Computing, Academic press, 1990 | Zbl

[SY05] Shen, Yun-Qiu; Ypma, Tjalling J. Newton’s method for singular nonlinear equations using approximate left and right nullspaces of the jacobian, Appl. Numer. Math., Volume 54 (2005) no. 2, pp. 256-265 | DOI | MR | Zbl

[Wey12] Weyl, Hermann Das asymptotische Verteilungsgesetz der Eigenwerte linearer partieller Differentialgleichungen (mit einer Anwendung auf die Theorie der Hohlraumstrahlung), Math. Ann., Volume 71 (1912) no. 4, pp. 441-479 | DOI | MR | Zbl

[XL08] Xu, Xiubin; Li, Chong Convergence criterion of Newton’s method for singular systems with constant rank derivatives, J. Math. Anal. Appl., Volume 345 (2008) no. 2, pp. 689-701 | MR | Zbl

[Yak90] Yakoubsohn, Jean-Claude On Newton’s rule and Sylvester’s theorems, J. Pure Appl. Algebra, Volume 65 (1990) no. 3, pp. 293-309 | DOI | MR | Zbl

[Yam83] Yamamoto, Norio Newton’s method for singular problems and its application to boundary value problems, J. Math., Tokushima Univ, Volume 17 (1983), pp. 27-88 | MR | Zbl