Dynamics of four families of methods with the same weight function to solve nonlinear equations

Titolo del Lavoro

Dinamica di quattro famiglie di metodi con la stessa funzione di peso per la risoluzione di equazioni non lineari.
Autrici: Livia J. Quiñonez T. e Carlos E. Cadenas R.

1. Il Problema

Lo studio affronta la necessità di analizzare il comportamento dinamico di metodi iterativi utilizzati per risolvere equazioni non lineari della forma $f(x) = 0$. Sebbene esistano metodi classici ben noti (come Newton, Halley e Chebyshev), la letteratura richiede un'analisi più approfondita delle famiglie di metodi generati tramite funzioni di peso.

Il problema specifico trattato è la dinamica complessa di quattro diverse famiglie di metodi ottenuti combinando convessamente il metodo di Newton e un metodo di tipo Newton-Halley. Queste famiglie, sebbene definite da equazioni iterative diverse per casi generali (radici multiple o metodi a due punti), mostrano un comportamento dinamico identico quando applicate a polinomi quadratici con due radici distinte. L'obiettivo è identificare i parametri che garantiscono un comportamento stabile e convergente, evitando comportamenti caotici o la convergenza verso punti fissi indesiderati.

2. Metodologia

Gli autori hanno adottato un approccio basato sulla dinamica complessa e sulla teoria dei sistemi dinamici razionali:

• Definizione delle Famiglie: Sono state considerate quattro famiglie di metodi (etichettate come (1), (3), (5) e (7) nel testo), tutte dipendenti da un parametro libero $A$ e basate sulla stessa funzione di peso $H(t)$.
• Riduzione tramite Teorema di Scala: Utilizzando il Scaling Theorem, gli autori hanno dimostrato che lo studio della dinamica per polinomi generici di grado 2 può essere ridotto allo studio di un'unica mappa razionale universale $S(z)$, indipendente dalle radici specifiche del polinomio. La mappa coniugata è:
  $$S(z) = \frac{z^3((A - 1)z + 2A - 1)}{(2A - 1)z + A - 1}$$
• Analisi dei Punti Fissi e Critici:
  • Sono state derivate le espressioni analitiche per i punti fissi (soluzioni di $S(z)=z$) e i punti critici (soluzioni di $S'(z)=0$).
  • Sono stati classificati i punti fissi in base alla loro stabilità (attrattori, repulsivi, parabolici, superattrattori) calcolando il modulo della derivata prima $|S'(z)|$.
• Spazi dei Parametri e Piani Dinamici:
  • Sono stati generati gli spazi dei parametri iterando il punto critico libero $z_{c1}$ su una griglia complessa di $1000 \times 1000$ punti per il parametro $A$.
  • Sono stati costruiti piani dinamici (basins of attraction) per valori specifici di $A$ selezionati dallo spazio dei parametri, visualizzando le regioni di convergenza verso le radici (0 e $\infty$) e verso i punti fissi "strani" (non associati alle radici).
  • È stata condotta un'analisi numerica per identificare orbite periodiche di periodo due.

3. Contributi Chiave

1. Unificazione Dinamica: Dimostrazione che quattro diverse famiglie di metodi (metodo a un punto, combinazione per radici multiple, metodi a due punti) condividono la stessa classe di coniugazione analitica quando applicati a polinomi quadratici.
2. Classificazione della Stabilità:
   • Identificazione analitica delle condizioni di stabilità per i punti fissi strani $z=1$ e $z_{2,3}$.
   • Dimostrazione che $z=1$ è un attrattore superattrattivo per $A = 3/4$ e che i punti $z_{2,3}$ sono superattrattivi per $A = 5/6$.
   • Definizione delle regioni di stabilità nel piano complesso del parametro $A$ (dischi di raggio $2/55$ e $1/35$ centrati rispettivamente su $42/55$ e $29/35$).
3. Mappatura dei Comportamenti: Creazione di mappe dettagliate che distinguono tra:
   • Convergenza alle radici (comportamento desiderato).
   • Convergenza a punti fissi strani (comportamento indesiderato).
   • Comportamento caotico o non convergente (punti blu nelle figure).
4. Scoperta di Orbite Periodiche: Identificazione e calcolo di orbite periodiche di periodo due per specifici valori del parametro, evidenziando la complessità del sistema in certe regioni.

4. Risultati Principali

• Casi Speciali:
  • $A=1$ (Metodo di Newton): $S(z)=z^2$, punti fissi strani repulsivi.
  • $A=0$ (Metodo di Halley): $S(z)=z^3$, punti fissi strani repulsivi.
  • $A=1/2$ (Metodo Super-Halley): $S(z)=z^4$, punti fissi strani repulsivi.
  • $A=2/3$: $S(z)=-z^3$, punti fissi strani repulsivi.
• Stabilità dei Punti Fissi Strani:
  • Il punto fisso $z=1$ è stabile (attrattore) solo se $|A - 42/55| < 2/55$. Al di fuori di questo intervallo, diventa repulsivo.
  • I punti fissi $z_{2,3}$ sono stabili solo se $|A - 29/35| < 1/35$.
• Spazi dei Parametri: Le figure (4 e 5) mostrano che esistono ampie regioni del piano complesso di $A$ dove il metodo non converge alle radici, ma verso punti fissi strani o non converge affatto. Le regioni rosse indicano la convergenza alle radici (buon comportamento), mentre le regioni blu indicano non-convergenza.
• Piani Dinamici: Le analisi visive (Figure 6-10) confermano che per valori reali di $A$ lontani dalle zone di stabilità dei punti fissi strani, i bacini di attrazione delle radici sono ben definiti e separati da una circonferenza unitaria. Tuttavia, per valori complessi o specifici reali, si osservano distorsioni significative e la comparsa di orbite periodiche.

5. Significato e Implicazioni

Questo lavoro è significativo per la comunità della matematica computazionale e dell'analisi numerica perché:

• Guida alla Selezione dei Parametri: Fornisce criteri rigorosi per scegliere il parametro $A$ in modo da massimizzare la probabilità di convergenza alle radici e minimizzare il rischio di convergenza a punti fissi strani o cicli periodici.
• Comprensione Teorica: Chiarisce la relazione tra diverse famiglie di metodi iterativi, mostrando che la loro essenza dinamica è la stessa in contesti specifici, semplificando l'analisi futura.
• Visualizzazione del Rischio: Le mappe dei parametri e i piani dinamici offrono uno strumento visivo immediato per identificare "zone pericolose" dove l'uso del metodo potrebbe fallire, permettendo agli ingegneri e ai ricercatori di evitare tali configurazioni parametriche.
• Scoperta di Comportamenti Complessi: L'identificazione di orbite periodiche di periodo due evidenzia la ricchezza dinamica di queste famiglie, suggerendo che anche metodi di alto ordine di convergenza possono esibire comportamenti caotici se i parametri non sono scelti con cura.

In sintesi, il paper offre una mappa completa della dinamica di queste famiglie di metodi, trasformando una scelta parametrica spesso empirica in una decisione basata su solide fondamenta teoriche e visive.