Risolutore di Sistemi di EDO

Risolutore di Sistemi di EDO

Il Risolutore di Sistemi di EDO è un kit completo per equazioni differenziali per sistemi accoppiati lineari e non lineari. Inserisci una matrice dei coefficienti 2×2 o 3×3 e lo strumento eseguirà l'analisi completa di autovalori / autovettori, scriverà la soluzione generale e particolare in LaTeX, classificherà l'equilibrio all'origine come sella, nodo, spirale o centro e disegnerà un ritratto di fase interattivo con una traiettoria animata. Per i sistemi planari non lineari è possibile digitare i membri destri arbitrari \(f(x,y)\) e \(g(x,y)\) e lo strumento produrrà un ritratto di fase RK4 ad alta precisione.

Cos'è un Sistema di EDO?

Un sistema di equazioni differenziali ordinarie accoppia diverse funzioni incognite di una singola variabile — solitamente il tempo \(t\) — attraverso le loro derivate. Nella forma più compatta,

Quando \(\mathbf{F}(t, \mathbf{x}) = A\mathbf{x}\) per una matrice costante \(A\), il sistema è lineare e autonomo — ed è qui che la teoria è più affascinante: l'intero comportamento a lungo termine è determinato dagli autovalori di \(A\).

La Ricetta degli Autovalori per i Sistemi Lineari

Per \(\mathbf{x}' = A\mathbf{x}\) il metodo standard è:

1. Calcolare il polinomio caratteristico \(\det(\lambda I - A) = 0\).
2. Risolvere per gli autovalori \(\lambda_1, \lambda_2, \dots\).
3. Per ogni autovalore, trovare un autovettore \(v\) risolvendo \((A - \lambda I) v = 0\).
4. Assemblare la soluzione generale come combinazione lineare: \(\mathbf{x}(t) = c_1 v_1 e^{\lambda_1 t} + c_2 v_2 e^{\lambda_2 t} + \cdots\).
5. Determinare le costanti \(c_i\) inserendo la condizione iniziale \(\mathbf{x}(0)\) nella soluzione generale.

Tre Casi per i Sistemi 2×2

Il Piano Traccia-Determinante

Per una matrice 2×2 con traccia \(T = a_{11} + a_{22}\) e determinante \(D = a_{11} a_{22} - a_{12} a_{21}\), l'intera classificazione rientra in un unico schema:

Ecco perché il pannello dei risultati mostra in evidenza \(T\), \(D\) e \(\Delta = T^2 - 4D\) — tre numeri sono sufficienti per nominare l'equilibrio.

Sistemi Non Lineari e Ritratto di Fase

La maggior parte delle EDO del mondo reale sono non lineari e non hanno una soluzione in forma chiusa. Lo strumento le gestisce integrando numericamente le equazioni con un metodo Runge–Kutta del 4° ordine (RK4), che ha un errore di troncamento locale \(O(h^5)\) ed è lo standard per i campi vettoriali regolari.

Il ritratto di fase sovrappone:

• Un campo vettoriale campionato su una griglia 13×13, che mostra la direzione del flusso in ogni punto.
• La traiettoria dalla tua condizione iniziale, disegnata in rosso con un corridore arancione animato che mostra la direzione del tempo.
• Diverse linee di flusso da un anello di punti di partenza, fornendo un quadro globale della dinamica.
• Per i sistemi lineari 2×2, gli assi degli autovettori (ciano tratteggiato) — queste sono le direzioni invarianti lungo le quali le soluzioni scorrono esponenzialmente.

Come Usare Questo Risolutore

1. Scegli una modalità — Lineare 2×2, Lineare 3×3 o Non lineare 2D — tramite le schede nella parte superiore del modulo.
2. Inserisci i coefficienti o le equazioni. Clicca su qualsiasi Esempio Rapido per precompilare un sistema canonico (nodo stabile, centro, sella, pendolo, Van der Pol, ecc.).
3. Inserisci una condizione iniziale \((x_0, y_0)\) e un intervallo di tempo \(T\). Valori di \(T\) tipici sono 6–20 per gli oscillatori e 3–6 per i sistemi stabili a decadimento rapido.
4. Clicca su Risolvi. Apparirà la pagina completa dei risultati con la classificazione, gli autovalori, gli autovettori, la soluzione in forma chiusa (modalità lineari), il ritratto di fase animato e il grafico delle serie temporali.
5. Riavvia la traiettoria usando il pulsante sotto il ritratto di fase se desideri guardare di nuovo il corridore attraversare la curva CI.

Esempio Svolto — Oscillatore Armonico Smorzato

L'oscillatore smorzato \(\ddot{x} + 2\zeta \omega \dot{x} + \omega^2 x = 0\) può essere riscritto come un sistema 2D ponendo \(y = \dot{x}\):

Per \(\omega = 1\) e \(\zeta = 0.2\) (sottosmorzato), la matrice è \(A = \begin{pmatrix} 0 & 1 \\ -1 & -0.4 \end{pmatrix}\). Traccia \(T = -0.4\), determinante \(D = 1\), discriminante \(\Delta = 0.16 - 4 = -3.84 < 0\), quindi otteniamo una spirale stabile con autovalori \(-0.2 \pm 0.9798\,i\). La traiettoria spiraleggia verso l'origine e la serie temporale mostra sinusoidi a decadimento esponenziale.

Applicazioni

• Sistemi meccanici — sistemi molla-massa accoppiati, pendoli, giroscopi.
• Circuiti elettrici — reti RLC, filtri con amplificatori operazionali, controllo in spazio di stato.
• Dinamica delle popolazioni — predatore-preda di Lotka–Volterra, specie in competizione, epidemiologia (SIR, SIS).
• Cinetica chimica — reti di reazione, oscillatori di Belousov–Zhabotinsky.
• Neuroscienze — modello di neurone di FitzHugh–Nagumo, riduzioni di Hodgkin–Huxley.
• Teoria del controllo — modelli di impianto linearizzati, progettazione di osservatori, margini di stabilità.

Suggerimenti e Accorgimenti

• Se la traiettoria diverge rapidamente, riduci l'intervallo di tempo T — un sistema instabile può superare i limiti di qualsiasi visualizzazione in poche unità di tempo.
• Per un autovalore ripetuto, il risolutore trova automaticamente l'autovettore generalizzato \(w\) risolvendo \((A - \lambda I)w = v\), così ottieni il termine \(tv\) senza calcoli manuali.
• Per i sistemi non lineari, le frecce del campo vettoriale rivelano anche equilibri non nell'origine come punti ciano — osserva il ritratto per individuare le regioni a magnitudo zero.
• Per i sistemi 3×3 non è presente il ritratto di fase (il 3D è difficile da visualizzare in una pagina 2D), ma la soluzione in forma chiusa e il verdetto di stabilità sono comunque validi.
• Le condizioni iniziali e gli intervalli di tempo sono separati dalla classificazione: cambiarli sposta solo la traiettoria rossa, non il verdetto sugli autovalori.

Domande Frequenti

Cos'è un sistema di equazioni differenziali ordinarie?

Un sistema di equazioni differenziali ordinarie (EDO) è un insieme di equazioni accoppiate che mettono in relazione le derivate di diverse funzioni incognite di una singola variabile indipendente, tipicamente il tempo. La forma classica è \( \mathbf{x}'(t) = F(t, \mathbf{x}(t)) \), dove \( \mathbf{x} \) è un vettore di stati e \(F\) è il campo vettoriale. I sistemi lineari possono essere scritti in modo compatto come \( \mathbf{x}' = A\mathbf{x} + \mathbf{b} \), e il loro comportamento è determinato quasi interamente dagli autovalori della matrice dei coefficienti \(A\).

In che modo gli autovalori classificano l'equilibrio di un sistema lineare 2×2?

Per un sistema 2×2 \( \mathbf{x}' = A\mathbf{x} \) l'origine è classificata dalla traccia \(T\) e dal determinante \(D\) di \(A\): \(D < 0\) fornisce un punto di sella (instabile); \(D > 0\) con \(T^2 > 4D\) fornisce un nodo (stabile se \(T < 0\), instabile se \(T > 0\)); \(D > 0\) con \(T^2 < 4D\) fornisce una spirale (stabile se \(T < 0\), instabile se \(T > 0\), un centro puro se \(T = 0\)). Il caso limite \(T^2 = 4D\) produce un nodo degenere.

Com'è fatta la soluzione in forma chiusa quando gli autovalori sono complessi?

Se \(A\) ha autovalori complessi coniugati \( \alpha \pm i\beta \) con autovettore complesso \( v = p + iq \), la soluzione generale reale è \( \mathbf{x}(t) = e^{\alpha t} \left[ c_1 (p \cos\beta t - q \sin\beta t) + c_2 (p \sin\beta t + q \cos\beta t) \right] \). L'esponenziale \(e^{\alpha t}\) controlla l'ampiezza (crescente, decrescente o costante), mentre il seno e il coseno gestiscono la rotazione.

Cosa succede quando la matrice ha un autovalore ripetuto?

Se la matrice ha un autovalore ripetuto \(\lambda\) ma un solo autovettore linearmente indipendente \(v\), è necessario anche un autovettore generalizzato \(w\) che risolva \( (A - \lambda I) w = v \). La soluzione generale assume quindi la forma \( \mathbf{x}(t) = c_1 v e^{\lambda t} + c_2 (tv + w) e^{\lambda t} \). Se lo spazio proprio risulta essere bidimensionale, la matrice è un multiplo scalare dell'identità su quel sottospazio invariante e la soluzione si riduce a \( \mathbf{x}(t) = (c_1 v_1 + c_2 v_2) e^{\lambda t} \).

Questo strumento può risolvere sistemi non lineari simbolicamente?

La modalità non lineare risolve il sistema numericamente utilizzando un integratore Runge–Kutta del quarto ordine (RK4) e traccia il ritratto di fase. La maggior parte dei sistemi non lineari non ha una soluzione in forma chiusa, quindi questo è l'approccio standard. È comunque possibile leggere il comportamento locale vicino agli equilibri linearizzando, operazione gestita dalla modalità lineare 2×2 — calcola lo Jacobiano nel punto fisso e inseriscilo come \(A\).

Cos'è un ritratto di fase?

Un ritratto di fase è una rappresentazione geometrica delle soluzioni di un sistema 2D nel piano \(x\)–\(y\). Ogni soluzione traccia una curva chiamata traiettoria, e l'insieme delle traiettorie insieme alle frecce del campo vettoriale rivela il comportamento qualitativo: se le soluzioni spiraleggiano verso l'interno, si allontanano a sella, oscillano o si stabilizzano su equilibri. I ritratti di fase rendono visibile a colpo d'occhio la struttura globale di un sistema.

Ulteriori letture

• Sistema di equazioni differenziali — Wikipedia
• Ritratto di fase — Wikipedia
• Autovettori e autovalori — Wikipedia
• Metodi di Runge–Kutta — Wikipedia
• Oscillatore di Van der Pol — Wikipedia