Equazione differenziale ordinaria

Da Wikipedia, l'enciclopedia libera.
Vai alla navigazione Vai alla ricerca

In matematica, un'equazione differenziale ordinaria (abbreviata in EDO, oppure ODE dall'acronimo inglese Ordinary Differential Equation) è un'equazione differenziale che coinvolge una funzione di una variabile e le sue derivate di ordine qualsiasi: si tratta di un oggetto matematico estensivamente utilizzato in fisica e in molti altri ambiti della scienza; ad esempio un sistema dinamico viene descritto da un'equazione differenziale ordinaria.

Come succede per tutte le equazioni differenziali, solitamente non è possibile risolvere esattamente una EDO e comunque non esistono metodi generali per farlo. I diversi casi possibili sono pertanto analizzati singolarmente, e spesso ci si limita a studiare il comportamento qualitativo della soluzione senza che sia possibile ottenerne un'espressione analitica. Particolarmente semplici risultano le equazioni lineari (di qualunque ordine) poiché si possono sempre ricondurre ad un sistema di equazioni lineari del primo ordine.

Sia , con un insieme aperto e connesso e .

Si definisce equazione differenziale ordinaria di ordine una relazione del tipo:

dove con si indica la derivata -esima della funzione .

Se è definita in una regione dello spazio euclideo , allora si considerano più propriamente equazioni differenziali ordinarie nel campo reale, a valori reali se è a valori in .

L'ordine di un'equazione è l'ordine massimo di derivazione che vi compare, mentre l'aggettivo ordinario si riferisce al fatto che l'incognita è una funzione di una sola variabile. Nel caso l'incognita dipenda da più variabili si ha un'equazione differenziale alle derivate parziali.

Sia un intervallo di . Si definisce soluzione o integrale dell'equazione differenziale ordinaria una funzione tale che:

Un'equazione differenziale ordinaria si dice autonoma se non dipende esplicitamente da .

Un'equazione differenziale ordinaria si dice scritta in forma normale se può essere esplicitata rispetto :

Si dice inoltre lineare se è combinazione lineare di , ossia:

o, in modo equivalente:

dove:

Il termine è detto sorgente o forzante, e se è nullo l'equazione differenziale lineare si dice omogenea, se non è nullo l'equazione differenziale lineare si dice completa.

Un'equazione differenziale ordinaria lineare possiede soluzioni linearmente indipendenti in numero pari al grado dell'equazione, e ogni loro combinazione lineare è a sua volta soluzione.

Data un'equazione differenziale ordinaria lineare, nel caso si conosca una soluzione generale dell'equazione omogenea ad essa associata allora è possibile trovare una soluzione particolare dell'equazione "completa". A tal fine esistono diverse procedure (che riguardano anche le soluzioni di equazioni ordinarie non lineari), tra cui il metodo delle variazioni delle costanti e l'utilizzo della trasformata di Laplace per le equazioni ordinarie lineari. Per i casi più semplici vi sono inoltre alcune teorie: ad esempio, per le equazioni di primo grado è possibile cercare un opportuno fattore di integrazione, per quelle di secondo vi è la teoria di Sturm-Liouville. In generale, tuttavia, di solito l'unico modo possibile per studiare la soluzione è l'utilizzo di un metodo di soluzione numerica.

Sistemi di ODE

[modifica | modifica wikitesto]

Riduzione a sistema di equazioni di ordine 1

[modifica | modifica wikitesto]

Un sistema di equazioni differenziali ordinarie di ordine in forma normale è una relazione vettoriale del tipo:

Una soluzione classica di un tale sistema è una funzione tale che:

Di particolare rilevanza ai fini pratici è la riduzione di un'equazione differenziale ordinaria di ordine in forma normale ad un sistema differenziale del primo ordine. Questa tecnica permette di semplificare notevolmente alcuni tipi di problemi, evitando l'introduzione di complesse forme di risoluzione. Sia:

un'equazione differenziale di ordine di tipo normale. Si definiscono:

in modo che , ed in particolare . L'equazione differenziale è dunque equivalente al sistema:

Ponendo:

si ottiene:

ossia, detto in altri termini, si può sempre tradurre tutto in un'equazione di ordine 1. Con un procedimento del tutto analogo a quello seguito è possibile anche il viceversa, ossia ottenere un'equazione di ordine partendo da una equazione di ordine 1 in cui un vettore di dimensione sia soluzione.

Sistema di equazioni di ordine n

[modifica | modifica wikitesto]

Se si considera un vettore definito come:

e una funzione che agisce su e le sue derivate, allora la scrittura:

denota un sistema esplicito di equazioni differenziali ordinarie di ordine . In forma di vettori colonna si ha:

L'analogo sistema implicito è:

dove . In forma matriciale:

Data un'equazione:

una funzione è detta soluzione (o integrale) dell'equazione differenziale ordinaria se è differenziabile volte su e si ha:

Date due soluzioni e , è detta estensione di se e:

Una soluzione che non possiede estensioni è detta soluzione massimale, mentre una soluzione definita su tutto è detta soluzione globale.

Una soluzione generale di un'equazione di ordine è una soluzione contenente costanti di integrazione indipendenti, mentre una soluzione particolare è ottenuta dalla soluzione generale conferendo un valore fissato alle costanti, solitamente in modo da soddisfare le condizioni iniziali o condizioni al contorno. In tale contesto, una soluzione singolare è una soluzione che non può essere ottenuta assegnando un valore definito alle costanti di integrazione.

Esistenza della soluzione e problema di Cauchy

[modifica | modifica wikitesto]
Lo stesso argomento in dettaglio: Problema di Cauchy.

Un problema ai valori iniziali è un'equazione differenziale ordinaria:[1]

cui è associato un punto nel dominio di :

chiamato condizione iniziale. La soluzione di un problema ai valori iniziali è quindi una funzione che è soluzione dell'equazione differenziale e soddisfa la condizione . In altri termini, il problema di Cauchy consiste nel trovare una curva , tra quelle definite da , che passi per il punto .

L'esistenza locale di una soluzione fu provata da Augustin-Louis Cauchy sotto l'ipotesi di continuità e limitatezza di in una regione del suo dominio (teorema di Peano o Cauchy-Peano).[2] Successivamente l'esistenza e l'unicità locali furono mostrate da Émile Picard con l'ipotesi di lipschitzianità rispetto a (teorema di esistenza e unicità per un problema di Cauchy o Picard–Lindelöf), e tale risultato può essere esteso ad una forma globale. Se quindi è lipschitziana in una regione del dominio allora esiste almeno una curva-soluzione differenziabile con continuità passante per ogni punto interno a . Semplificando la questione, se e sono continue in un rettangolo chiuso nel piano della forma:

dove e è il prodotto cartesiano, allora vi è un intervallo:

nel quale per qualche l'unica soluzione può essere trovata.[3] Questo risultato si applica anche ad equazioni non lineari della forma così come a sistemi di equazioni.

Il teorema di Cauchy-Kovalevskaya, che si applica anche per le equazioni alle derivate parziali, mostra in modo più generale che se l'incognita e le condizioni iniziali di un'equazione differenziale sono localmente funzioni analitiche allora una soluzione analitica esiste ed è unica.[4] La funzione incognita può assumere valori su infiniti spazi dimensionali, come gli spazi di Banach o spazi di distribuzioni.

Esistenza e unicità locale

[modifica | modifica wikitesto]

Vi sono diversi teoremi che consentono di stabilire l'esistenza e l'eventuale unicità locale di soluzioni per dati problemi iniziali. I due principali sono il teorema di esistenza e unicità per un problema di Cauchy, il quale assume la lipschitzianità della funzione che definisce l'equazione ordinaria e ne conclude l'esistenza e unicità locale, e il teorema di esistenza di Peano, che assume la continuità conclude soltanto l'esistenza. Il primo afferma che dato il problema ai valori iniziali:

se è una funzione lipschitziana in e continua in allora per qualche esiste un'unica soluzione al problema ai valori iniziali sull'intervallo . Il teorema di esistenza di Peano assume soltanto continua, ma non garantisce l'unicità della soluzione, mostrandone solo la sua esistenza locale. Se essa esiste, o è unica o ne esistono infinite. Un'estensione di questo risultato si ha con il teorema di esistenza di Carathéodory, che si applica anche ai casi in cui l'equazione non è continua.

Unicità globale

[modifica | modifica wikitesto]

Se le ipotesi del teorema di esistenza e unicità per un problema di Cauchy sono soddisfatte, allora la condizione di esistenza locale può essere estesa ad un risultato globale. Per ogni condizione iniziale esiste un unico massimo intervallo aperto (eventualmente illimitato):

tale per cui ogni soluzione che soddisfa la condizione iniziale è una restrizione della soluzione soddisfacente tale condizione iniziale che è definita sul dominio . Nel caso in cui , vi sono solo due possibilità:

dove è l'aperto in cui è definita e è la sua frontiera.

Affinché la soluzione sia unica, il massimo dominio in cui è definita la soluzione deve essere un intervallo, che in generale dipende dalla condizione iniziale. Ad esempio, si consideri:

Dato che è lipschitziana, soddisfa il teorema di Picard–Lindelöf. La soluzione

ha come massimo intervallo per il suo dominio:

Equazioni del primo ordine

[modifica | modifica wikitesto]

Il caso più semplice è quello in cui:

con una funzione continua definita su un aperto di . Il problema si riconduce alla ricerca delle primitive di . In tal caso se è una soluzione allora anche , con , è ancora soluzione. Inoltre, se e sono soluzioni si ha che per qualche .

Data una condizione iniziale , scrivendo l'equazione nella forma , in modo che , si ha:

e la soluzione è fornita dal teorema fondamentale del calcolo integrale:

La soluzione è univocamente determinata dal dato iniziale : data una qualsiasi soluzione , si ha e pertanto .

Sistemi autonomi

[modifica | modifica wikitesto]
Lo stesso argomento in dettaglio: Sistema autonomo (matematica).

Un'equazione autonoma è un'equazione differenziale ordinaria del tipo:

dove è una funzione continua con derivata prima continua in tutto un intervallo , e che non dipende dalla variabile indipendente . Se è un vettore di si ha un sistema autonomo, ovvero un sistema di equazioni differenziali ordinarie autonome:

Di particolare importanza sono i punti tali per cui , detti punti di equilibrio, ai quali corrisponde la soluzione costante .

Un generico sistema di equazioni differenziali ordinarie (in cui dipende da ):

può essere reso autonomo introducendo una nuova incognita .

Equazione differenziale esatta

[modifica | modifica wikitesto]
Lo stesso argomento in dettaglio: Equazione differenziale esatta.

Si considerino un insieme semplicemente connesso e aperto e due funzioni e continue su . L'equazione differenziale implicita

è un'equazione differenziale esatta se esiste una funzione differenziabile con continuità , detta potenziale, tale che:

Il termine "esatta" si riferisce alla derivata totale di una funzione, detta talvolta "derivata esatta", che per una funzione è data in da:

Nelle applicazioni fisiche e non sono solitamente solo continue, ma anche differenziabili con continuità, e il teorema di Schwarz fornisce allora una condizione necessaria e sufficiente per l'esistenza della funzione potenziale (per equazioni definite su un insieme non semplicemente connesso tale criterio è solo necessario). Esso esiste se e solo se:

Soluzioni esatte

[modifica | modifica wikitesto]

Nel seguito si riportano alcuni importanti casi di equazioni differenziali ordinarie risolvibili esattamente.

Nella tabella le funzioni , , , , e sono funzioni integrabili in , , mentre e sono costanti reali date. Inoltre, sono costanti arbitrarie, in generale complesse. La notazione indica l'integrazione di rispetto a e la successiva sostituzione .

Equazione differenziale Metodo di soluzione Soluzione generale
Equazioni separabili
Primo ordine, separabile in x e y, caso generale.

Separazione delle variabili (divisione per P2Q1).
Primo ordine, separabile in x.

Integrazione diretta.
Primo ordine, autonoma, separabile in y.

Separazione delle variabili (divisione per F).
Primo ordine, separabile in x e y.

Integrazione su tutto il dominio.
Equazioni del primo ordine
Primo ordine, omogenea.

Si pone y = ux e si procede per separazione delle variabili in u e x.
Primo ordine, separabile.

Separazione delle variabili (divisione per xy).

Se N = M la soluzione è xy = C.

Differenziale esatto, primo ordine.

dove

Integrazione su tutto il dominio.
Differenziale inesatto, primo ordine.

dove

Il fattore di integrazione μ(x, y) soddisfa:

Se si può trovare μ(x, y):

Equazioni del secondo ordine
Secondo ordine, autonoma.

Si moltiplica l'equazione per 2dy/dx e si sostituisce quindi si integra due volte.
Equazioni lineari
Primo ordine, lineare, non omogenea.

Fattore di integrazione: .
Secondo ordine, lineare, non omogenea, a coefficienti costanti.

Funzione complementare yc: assumendo yc = eαx, si sostituisce e si risolve il polinomio in α, ottenendo le funzioni linearmente indipendenti .

Integrale particolare yp: in generale si applica il metodo delle variazioni delle costanti, sebbene talvolta è sufficiente studiare r(x).

Se b2 > 4c, allora:

Se b2 = 4c, allora:

Se b2 < 4c, allora:

Ordine n, lineare, non omogenea, a coefficienti costanti.

Funzione complementare yc: assumendo yc = eαx, si sostituisce e si risolve il polinomio in α, ottenendo le funzioni linearmente indipendenti .

Integrale particolare yp: in generale si applica il metodo delle variazioni delle costanti, sebbene talvolta sia sufficiente studiare r(x).

Dato che αj sono soluzioni del polinomio di grado n:

allora per αj tutti diversi:

per ogni radice αj ripetuta kj volte:

per qualche αj complesso, si pone α = χj + iγj e si usa la formula di Eulero per ottenere che alcuni termini del precedente risultato si possono scrivere come:

dove ϕj è una costante arbitraria.

In fisica, la modellizzazione di un sistema richiede spesso la risoluzione di un problema ai valori iniziali. In questo contesto, ad esempio, l'equazione differenziale può descrivere l'evoluzione di un sistema dinamico nel tempo in funzione delle condizioni iniziali: si consideri un punto materiale di massa in caduta libera, sotto l'azione della forza di gravità. Attraverso le leggi descritte da Newton sulla dinamica dei corpi, si ha che:

Riferendosi ad un sistema di coordinate cartesiane, con l'asse parallelo e discorde al verso dell'accelerazione di gravità e proiettando la relazione vettoriale precedente sugli assi coordinati, si ha che l'unica equazione significativa è quella rispetto all'asse :

Ciò fornisce . Per determinare la , soluzione del problema, occorre conoscere posizione e velocità iniziale del corpo in un certo istante. Integrando:

da cui:

e integrando nuovamente:

da cui si ottiene:

Come si vede, la soluzione dipende da due parametri e , rispettivamente velocità e posizione iniziale.

Si consideri il problema di Cauchy:

Il teorema di esistenza locale garantisce l'esistenza di almeno una soluzione . Una di queste è banale ed è anche globale, cioè definita su tutto l'insieme :

Oltre a questa soluzione è tuttavia possibile trovarne un'altra integrando l'equazione, trattandosi di un'equazione differenziale a variabili separabili. Perciò:

da cui:

e quindi:

dalla quale si ottiene:

Si possono così costruire delle funzioni che conservano la continuità a partire dalle due soluzioni trovate:

che sono anch'esse soluzioni del problema definite sull'insieme . Inoltre, fissato un e posto:

si possono costruire nuove soluzioni dello stesso problema a partire a questa:

definite anch'esse su . Tali soluzioni sono infinite, e se si rappresentano graficamente si ottiene una figura detta pennello di Peano. Il fatto che un problema di Cauchy può possedere infinite soluzioni è detto talvolta fenomeno di Peano, ed è dovuto al fatto che la derivata non è limitata nel punto : per far fronte a questo problema vi è il teorema di esistenza globale.

  1. ^ Encyclopedia of Mathematics - Cauchy problem, su encyclopediaofmath.org. URL consultato il 06-01-2013.
  2. ^ Encyclopedia of Mathematics - Peano theorem, su encyclopediaofmath.org. URL consultato il 06-01-2013.
  3. ^ Elementary Differential Equations and Boundary Value Problems (4th Edition), W.E. Boyce, R.C. Diprima, Wiley International, John Wiley & Sons, 1986, ISBN 0-471-83824-1
  4. ^ Encyclopedia of Mathematics - Kovalevskaya theorem, su encyclopediaofmath.org. URL consultato il 06-01-2013.

Voci correlate

[modifica | modifica wikitesto]

Altri progetti

[modifica | modifica wikitesto]

Collegamenti esterni

[modifica | modifica wikitesto]
Controllo di autoritàThesaurus BNCF 45336 · NDL (ENJA00574993
  Portale Matematica: accedi alle voci di Wikipedia che trattano di matematica