Calcolare La Funzione Di Stabilita Di Runge Kutta Matlav

Guida Completa alla Funzione di Stabilità di Runge-Kutta in MATLAB

La funzione di stabilità è un concetto fondamentale nell’analisi dei metodi numerici per la risoluzione di equazioni differenziali ordinarie (ODE). Quando si utilizzano metodi come Runge-Kutta in MATLAB, comprendere la funzione di stabilità R(hλ) è essenziale per valutare la stabilità numerica del metodo per diversi valori del passo h e degli autovalori λ del problema.

1. Cos’è la Funzione di Stabilità?

La funzione di stabilità R(z) (dove z = hλ) descrive come il metodo numerico risponde all’equazione test:

y’ = λy, y(0) = y₀

Per un metodo numerico, la soluzione approssimata dopo un passo è data da:

yₙ₊₁ = R(hλ) yₙ

La stabilità richiede che |R(hλ)| ≤ 1 per garantire che la soluzione numerica non diverga.

2. Funzioni di Stabilità per Diverse Varianti di Runge-Kutta

Ogni variante del metodo Runge-Kutta ha una funzione di stabilità diversa:

Metodo	Funzione di Stabilità R(z)	Regione di Stabilità Assoluta
Eulero Esplicito	R(z) = 1 + z	|1 + z| ≤ 1 (cerchio di raggio 1 centrato in -1)
RK2 (Midpoint)	R(z) = 1 + z + z²/2	Regione limitata nell’emisfero sinistro
RK4 Classico	R(z) = 1 + z + z²/2 + z³/6 + z⁴/24	Regione estesa ma non A-stabile
RKF45 (Adattivo)	Complessa (dipende dall’implementazione)	Variabile, spesso meno stabile di RK4

3. Come Calcolare la Funzione di Stabilità in MATLAB

In MATLAB, è possibile calcolare la funzione di stabilità per un metodo Runge-Kutta specifico utilizzando le seguenti tecniche:

1. Definizione della Funzione R(z): Implementare manualmente la funzione di stabilità per il metodo desiderato.
2. Uso di ode45 e Analisi Posteriore: Eseguire la simulazione e analizzare la risposta per diversi valori di hλ.
3. Toolbox Specializzati: Utilizzare toolbox come Chebfun per analisi avanzate.

% Esempio: Funzione di stabilità per RK4 R = @(z) 1 + z + z.^2/2 + z.^3/6 + z.^4/24; % Valutazione per z = h*lambda z_values = linspace(-5, 1, 1000) + 1i*linspace(-3, 3, 1000)’; z_values = z_values(:); stability = abs(R(z_values)) <= 1; % Plot della regione di stabilità figure; hold on; plot(z_values(stability), '.b', 'MarkerSize', 1); plot(z_values(~stability), '.r', 'MarkerSize', 1); xlabel('Re(z)'); ylabel('Im(z)'); title('Regione di Stabilità Assoluta per RK4'); axis equal; grid on;

4. Analisi della Stabilità Assoluta

La stabilità assoluta è cruciale per problemi stiff, dove gli autovalori λ hanno parti reali molto negative. Un metodo è A-stabile se la sua regione di stabilità include tutto il semipiano sinistro (Re(z) ≤ 0).

• Metodi Espliciti (RK4, Eulero): Non sono A-stabili. La regione di stabilità è limitata.
• Metodi Impliciti: Possono essere A-stabili (es. Euler Implicito, Trapezoidale).
• Metodi Adattivi (RKF45): La stabilità dipende dal controllo del passo.

Per problemi stiff, MATLAB offre solutori dedicati come ode15s o ode23s, che sono A-stabili o L-stabili.

5. Confronto tra Metodi Runge-Kutta

La seguente tabella confronta le proprietà di stabilità dei metodi più comuni:

Metodo	Ordine	A-Stabile	L-Stabile	Regione di Stabilità (Re(z) ≤ x)
Eulero Esplicito	1	No	No	-2 ≤ Re(z) ≤ 0
Eulero Implicito	1	Sì	Sì	Tutto Re(z) ≤ 0
RK2 (Midpoint)	2	No	No	-2 ≤ Re(z) ≤ 0
RK4 Classico	4	No	No	-2.78 ≤ Re(z) ≤ 0
Trapezoidale (RK2 Implicito)	2	Sì	No	Tutto Re(z) ≤ 0

6. Applicazioni Pratiche in MATLAB

Per analizzare la stabilità di un sistema usando Runge-Kutta in MATLAB:

1. Definisci il Problema: Specificare la funzione f(t, y) e le condizioni iniziali.
2. Scegli il Solutore: Selezionare ode45 (RK4) per problemi non stiff o ode15s per problemi stiff.
3. Analizza la Stabilità: Variare il passo h e osservare la risposta.
4. Visualizza i Risultati: Usare plot per confrontare soluzioni con diversi passi.

% Esempio: Analisi di stabilità per y’ = -20y + 20cos(t) tspan = [0 10]; y0 = 1; f = @(t, y) -20*(y – cos(t)); % Soluzione con ode45 (RK4) [ t1, y1 ] = ode45(f, tspan, y0); % Soluzione con passo più piccolo options = odeset(‘MaxStep’, 0.01); [ t2, y2 ] = ode45(f, tspan, y0, options); % Plot dei risultati figure; plot(t1, y1, ‘b-‘, t2, y2, ‘r–‘); xlabel(‘Tempo t’); ylabel(‘Soluzione y(t)’); legend(‘Passo Default’, ‘Passo Ridotto’); title(‘Influenza del Passo sulla Stabilità’);

7. Errori Comuni e Soluzioni

• Divergenza della Soluzione: Causata da un passo h troppo grande. Ridurre h o usare un metodo implicito.
• Oscillazioni Numeriche: Tipiche di metodi espliciti con autovalori immaginari. Usare metodi con migliore stabilità.
• Errore di Arrotondamento: Può dominare per passi molto piccoli. Usare precisione doppia (ode45 la usa di default).

8. Risorse Esterne Autorevoli

Per approfondire:

• MIT – Analisi di Stabilità per Runge-Kutta (PDF)
• UCLA – Stabilità dei Metodi Numerici (PDF)
• Documentazione Ufficiale MATLAB su ODE

Conclusione

La funzione di stabilità di Runge-Kutta è uno strumento potente per valutare l’affidabilità dei metodi numerici in MATLAB. Mentre i metodi espliciti come RK4 sono efficienti per problemi non stiff, i metodi impliciti o i solutori dedicati come ode15s sono essenziali per problemi stiff. Comprendere R(hλ) permette di scegliere il metodo e il passo ottimali, evitando instabilità e garantendo risultati accurati.