Nsryjdtyk

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.

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.

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.

Indice

1 Definizione
- 1.1 Sistemi di ODE
  - 1.1.1 Riduzione a sistema di equazioni di ordine 1
  - 1.1.2 Sistema di equazioni di ordine n
- 1.2 Soluzioni

2 Esistenza della soluzione e problema di Cauchy
- 2.1 Esistenza e unicità locale
- 2.2 Unicità globale

3 Equazioni del primo ordine
- 3.1 Sistemi autonomi
- 3.2 Equazione differenziale esatta

4 Soluzioni esatte

5 Applicazioni

6 Esempio

7 Note

8 Bibliografia

9 Voci correlate

10 Collegamenti esterni

Definizione |

Sia $F:Omega subseteq {mathbb {C}}^{{n+2}}rightarrow {mathbb {C}}$ , con $Omega neq varnothing$ un insieme aperto e connesso e $nin mathbb{N}$ .

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

Fleft(x,u(x),u'(x),ldots ,u^{{(n)}}(x)right)=0

dove con $u^{{(i)}}(x)$ si indica la derivata i-esima della funzione $u(x)$ .

Se $F!$ è definita in una regione $Omega !$ dello spazio euclideo $mathbb{R} ^{{n+2}}$ , allora si considerano più propriamente equazioni differenziali ordinarie nel campo reale, a valori reali se $F!$ è a valori in $mathbb {R}$ .

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 $I$ un intervallo di $mathbb {R}$ . Si definisce soluzione o integrale dell'equazione differenziale ordinaria una funzione $varphi =varphi (x)$ tale che:

varphi (x)in C^{n}(I)qquad Fleft(x,varphi (x),varphi '(x),ldots ,varphi ^{{(n)}}(x)right)=0qquad forall xin I

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

Un'equazione differenziale ordinaria si dice scritta in forma normale se può essere esplicitata rispetto $u^{{(n)}}(x)$ :

u^{{(n)}}(x)=Gleft(x,u,u',ldots ,u^{{(n-1)}}right)

Si dice inoltre lineare se $F$ è combinazione lineare di $u,u',ldots ,u^{{(n)}}$ , ovvero:

Fleft(x,u,u',ldots ,u^{{(n)}}right)=s(x)+b_{0}(x)u+b_{1}(x)u'+ldots +b_{{n}}(x)u^{{(n)}}

o, in modo equivalente:

u^{{(n)}}=sum _{{i=0}}^{{n-1}}a_{i}(x)u^{{(i)}}+r(x)

dove:

r(x),a_{0}(x),a_{1}(x),ldots ,a_{{n-1}}(x)in C^{0}(I)

Il termine $r(x)$ è detto sorgente o forzante, e se è nullo l'equazione differenziale lineare si dice omogenea.

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

Data un'equazione differenziale ordinaria, 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, tra cui il metodo delle variazioni delle costanti e l'utilizzo della trasformata di Laplace. 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 |

Riduzione a sistema di equazioni di ordine 1 |

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

{mathbf {u}}'={mathbf {f}}left(x,{mathbf {u}}right)

Una soluzione classica di un tale sistema è una funzione ${mathbf {u}}:Ito mathbb{R} ^{n}$ tale che:

{mathbf {u}}(x)in C^{1}(I)qquad {mathbf {u}}'(x)={mathbf {f}}left(x,{mathbf {u}}(x)right)qquad forall xin I

Di particolare rilevanza ai fini pratici è la riduzione di un'equazione differenziale ordinaria di ordine $n$ 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:

u^{{(n)}}(x)=Gleft(x,u,u',ldots ,u^{{(n-1)}}right)

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

u_{i}=u^{{(i-1)}}(x)qquad {mathbf {u}}=left(u_{i}right)qquad iin {1,ldots ,n}

in modo che $u_{{i}}'=u_{{i+1}}$ , ed in particolare $u_{{n}}'=u^{{(n)}}(x)$ . L'equazione differenziale è dunque equivalente al sistema:

{begin{cases}u_{1}'=u_{2}\u_{2}'=u_{3}\vdots \u_{{n-1}}'=u_{n}\u_{n}'=Gleft(x,u_{1},u_{2},ldots ,u_{n}right)=Gleft(x,{mathbf {u}}right)!end{cases}}

Ponendo:

{mathbf {f}}left(x,{mathbf {u}}right)={begin{pmatrix}u_{2}\u_{3}\vdots \u_{n}\Gleft(x,{mathbf {u}}right)end{pmatrix}}

si ottiene:

{mathbf {u}}'={mathbf {f}}left(x,{mathbf {u}}right)

ovvero, 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, ovvero ottenere un'equazione di ordine n partendo da una equazione di ordine 1 in cui un vettore $mathbf{u}$ di dimensione $n$ sia soluzione.

Sistema di equazioni di ordine n |

Se si considera un vettore ${mathbf y}$ definito come:

{mathbf y}(x)=(y_{1}(x),y_{2}(x),dots ,y_{m}(x))

e una funzione $mathbf F$ che agisce su ${mathbf y}$ e le sue derivate, allora la scrittura:

{mathbf {y}}^{{(n)}}={mathbf {F}}left(x,{mathbf {y}},{mathbf {y}}',{mathbf {y}}'',cdots {mathbf {y}}^{{(n-1)}}right)

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

{displaystyle {begin{pmatrix}y_{1}^{(n)}\y_{2}^{(n)}\vdots \y_{m}^{(n)}end{pmatrix}}={begin{pmatrix}F_{1}left(x,mathbf {y} ,mathbf {y} ',mathbf {y} '',cdots mathbf {y} ^{(n-1)}right)\F_{2}left(x,mathbf {y} ,mathbf {y} ',mathbf {y} '',cdots mathbf {y} ^{(n-1)}right)\vdots \F_{m}left(x,mathbf {y} ,mathbf {y} ',mathbf {y} '',cdots mathbf {y} ^{(n-1)}right)\end{pmatrix}}}

L'analogo sistema implicito è:

{mathbf {F}}left(x,{mathbf {y}},{mathbf {y}}',{mathbf {y}}'',cdots {mathbf {y}}^{{(n)}}right)={boldsymbol {0}}

dove ${mathbf 0}=(0,0,dots )$ . In forma matriciale:

{begin{pmatrix}F_{1}(x,{mathbf {y}},{mathbf {y}}',{mathbf {y}}'',cdots {mathbf {y}}^{{(n)}})\F_{2}(x,{mathbf {y}},{mathbf {y}}',{mathbf {y}}'',cdots {mathbf {y}}^{{(n)}})\vdots \F_{m}(x,{mathbf {y}},{mathbf {y}}',{mathbf {y}}'',cdots {mathbf {y}}^{{(n)}})\end{pmatrix}}={begin{pmatrix}0\0\vdots \0\end{pmatrix}}

Soluzioni |

Data un'equazione:

Fleft(x,y,y',cdots ,y^{{(n)}}right)=0

una funzione $u:Isubset mathbb{R} to mathbb{R}$ è detta soluzione (o integrale) dell'equazione differenziale ordinaria se $u$ è differenziabile n volte su $I$ e si ha:

F(x,u,u', cdots , u^{{(n)}})=0quad xin I

Date due soluzioni $u:Jsubset mathbb{R} to mathbb{R}$ e $v:Isubset mathbb{R} to mathbb{R}$ , $u$ è detta estensione di $v$ se $Isubset J$ e:

u(x)=v(x)quad xin I

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

Una soluzione generale di un'equazione di ordine n è una soluzione contenente n 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 |

.mw-parser-output .vedi-anche{border:1px solid #CCC;font-size:95%;margin-bottom:.5em}.mw-parser-output .vedi-anche td:first-child{padding:0 .5em}.mw-parser-output .vedi-anche td:last-child{width:100%}

Lo stesso argomento in dettaglio: Problema di Cauchy.

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

y'(t)=f(t,y(t))qquad f:Omega subset {mathbb {R}}times {mathbb {R}}^{n}

cui è associato un punto nel dominio di $f$ :

(t_{0},y_{0})in Omega

chiamato condizione iniziale. La soluzione di un problema ai valori iniziali è quindi una funzione $y$ che è soluzione dell'equazione differenziale e soddisfa la condizione $y(t_{0})=y_{0}$ . In altri termini, il problema di Cauchy consiste nel trovare una curva $y$ , tra quelle definite da $y'(t)=f(t,y(t))$ , che passi per il punto $y(t_{0})=y_{0}$ .

L'esistenza locale di una soluzione fu provata da Augustin-Louis Cauchy sotto l'ipotesi di continuità e limitatezza di $f$ 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 $y$ (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 $f$ è lipschitziana in una regione $D$ del dominio allora esiste almeno una curva-soluzione differenziabile con continuità passante per ogni punto interno a $D$ . Semplificando la questione, se $f$ e $partial f/partial y$ sono continue in un rettangolo chiuso nel piano x-y della forma:

R=[x_{0}-a,x_{0}+a]times [y_{0}-b,y_{0}+b]

dove $a,bin mathbb{R}$ e $times$ è il prodotto cartesiano, allora vi è un intervallo:

I=[x_{0}-h,x_{0}+h]subset [x_{0}-a,x_{0}+a]

dove per qualche $hin mathbb{R}$ l'unica soluzione può essere trovata.^[3] Questo risultato si applica anche ad equazioni non lineari della forma $F(x,y)$ 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 $y$ può assumere valori su infiniti spazi dimensionali, come gli spazi di Banach o spazi di distribuzioni.

Esistenza e unicità locale |

Lo stesso argomento in dettaglio: Teorema di esistenza e unicità per un problema di Cauchy e Teorema di esistenza di Peano.

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:

y'(t)=f(t,y(t))qquad y(t_{0})=y_{0}qquad tin [t_{0}-varepsilon ,t_{0}+varepsilon ]

se $f$ è una funzione lipschitziana in $y$ e continua in $t$ allora per qualche $epsilon >0$ esiste un'unica soluzione $y(t)$ al problema ai valori iniziali sull'intervallo $[t_{0}-epsilon ,t_{0}+epsilon ]$ . Il teorema di esistenza di Peano assume $F(x,y)$ 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 |

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 $(x_0,y_0)$ esiste un unico massimo intervallo aperto (che può avere estensione infinita):

I_{{max}}=(x_{-},x_{+})qquad x_{pm }in {mathbb {R}}quad x_{0}in I_{{max}}

tale per cui ogni soluzione che soddisfa la condizione iniziale è una restrizione della soluzione soddisfacente tale condizione iniziale che è definita sul dominio $I_{{max}}$ . Nel caso in cui $x_{pm }nrightarrow pm infty$ , vi sono solo due possibilità:

lim _{{xto x_{pm }}}|y(x)|rightarrow infty qquad lim _{{xto x_{pm }}}in partial {bar {Omega }}

dove $Omega$ è l'aperto in cui $F$ è definita e $partial {bar {Omega }}$ è 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:

y'=y^{2}

Dato che $F(x,y)=y^{2}$ è lipschitziana, soddisfa il teorema di Picard–Lindelöf. La soluzione:

y(x)={frac {y_{0}}{(x_{0}-x)y_{0}+1}}

ha come massimo intervallo per il suo dominio:

{begin{cases}{mathbb {R}}&y_{0}=0\(-infty ,x_{0}+{frac {1}{y_{0}}})&y_{0}>0\(x_{0}+{frac {1}{y_{0}}},+infty )&y_{0}<0end{cases}}

Equazioni del primo ordine |

Il caso più semplice è quello in cui:

y'(t)=f(t)

con $f$ una funzione continua definita su un aperto di $R$ . Il problema si riconduce alla ricerca delle primitive $y(t)$ di $f(t)$ . In tal caso se $y(t)$ è una soluzione allora anche $y(t)+c$ , con $cin mathbb{R}$ , è ancora soluzione. Inoltre, se $y_{1}(t)$ e $y_{2}(t)$ sono soluzioni si ha che $y_{2}=y_{1}+c$ per qualche $cin mathbb{R}$ .

Data una condizione iniziale $y(t_{0})=y_{0}$ , scrivendo l'equazione nella forma $dy/dt=f(t)$ , in modo che $dy=f(t)dt$ , si ha:

int _{{y_{0}}}^{y}dxi =int _{{t_{0}}}^{t}f(tau )dtau

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

y(t)=y_{0}+int _{{t_{0}}}^{t}f(tau )dtau

La soluzione è univocamente determinata dal dato iniziale $y(t_{0})=y_{0}$ : data una qualsiasi soluzione $y_{1}(t)$ , si ha $y_{1}(t_{0})+c=y_{0}$ e pertanto $c=y_{0}-y_{1}(t_{0})$ .

Sistemi autonomi |

Lo stesso argomento in dettaglio: Sistema autonomo (matematica).

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

y'(t)=f(y(t))

dove $f$ è una funzione continua con derivata prima continua in tutto un intervallo $Isubset mathbb{R}$ , e che non dipende dalla variabile indipendente $t$ . Se $y$ è un vettore di $mathbb {R} ^{n}$ si ha un sistema autonomo, ovvero un sistema di equazioni differenziali ordinarie autonome:

y'_{i}=f_{i}(y_{1},y_{2},dots ,y_{n})qquad i=1,dots ,n

Di particolare importanza sono i punti $x_0$ tali per cui $f(x_{0})=0$ , detti punti di equilibrio, ai quali corrisponde la soluzione costante $y=x_{0}$ .

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

{frac {d}{dt}}y(t)=f(y(t),t)

può essere reso autonomo introducendo una nuova incognita $y_{{n+1}}=t$ .

Equazione differenziale esatta |

Lo stesso argomento in dettaglio: Equazione differenziale esatta.

Si considerino un insieme semplicemente connesso e aperto $Dsubset mathbb{R} ^{2}$ e due funzioni $I$ e $J$ continue su $D$ . L'equazione differenziale implicita:

I(x,y),{mathrm {d}}x+J(x,y),{mathrm {d}}y=0

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

{frac {partial F}{partial x}}(x,y)=Iqquad {frac {partial F}{partial y}}(x,y)=J

Il termine "esatta" si riferisce alla derivata totale di una funzione, detta talvolta "derivata esatta", che per una funzione $F(x_{0},x_{1},dots ,x_{{n-1}},x_{n})$ è data in $x_0$ da:

{frac {{mathrm {d}}F}{{mathrm {d}}x_{0}}}={frac {partial F}{partial x_{0}}}+sum _{{i=1}}^{{n}}{frac {partial F}{partial x_{i}}}{frac {{mathrm {d}}x_{i}}{{mathrm {d}}x_{0}}}

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

{frac {partial I}{partial y}}(x,y)={frac {partial J}{partial x}}(x,y)

Soluzioni esatte |

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

Nella tabella le funzioni $P(x)$ , $Q(x)$ , $P(y)$ , $Q(y)$ , $M(x,y)$ e $N(x,y)$ sono funzioni integrabili in $x$ , $y$ , mentre $b$ e $c$ sono costanti reali date. Inoltre, $C_{1},C_{2},dots$ sono costanti arbitrarie, in generale complesse. La notazione $int ^{x}F(lambda )dlambda$ indica l'integrazione di $F(lambda )$ rispetto a $lambda$ e la successiva sostituzione $lambda =x$ .

Equazione differenziale	Metodo di soluzione	Soluzione generale
Equazioni separabili
Primo ordine, separabile in x e y, caso generale. $P_{1}(x)Q_{1}(y)+P_{2}(x)Q_{2}(y),{frac {dy}{dx}}=0,!$ $P_{1}(x)Q_{1}(y),dx+P_{2}(x)Q_{2}(y),dy=0,!$	Separazione delle variabili (divisione per P₂Q₁).	$int ^{x}{frac {P_{1}(lambda )}{P_{2}(lambda )}},dlambda +int ^{y}{frac {Q_{2}(lambda )}{Q_{1}(lambda )}},dlambda =C,!$
Primo ordine, separabile in x. ${frac {dy}{dx}}=F(x),!$ $dy=F(x),dx,!$	Integrazione diretta.	$y=int ^{x}F(lambda ),dlambda +C,!$
Primo ordine, autonoma, separabile in y. ${frac {dy}{dx}}=F(y),!$ $dy=F(y),dx,!$	Separazione delle variabili (divisione per F).	$x=int ^{y}{frac {dlambda }{F(lambda )}}+C,!$
Primo ordine, separabile in x e y. $P(y){frac {dy}{dx}}+Q(x)=0,!$ $P(y),dy+Q(x),dx=0,!$	Integrazione su tutto il dominio.	$int ^{y}P(lambda ),{dlambda }+int ^{x}Q(lambda ),dlambda =C,!$
Equazioni del primo ordine
Primo ordine, omogenea. ${frac {dy}{dx}}=Fleft({frac {y}{x}}right),!$	Si pone y = ux e si procede per separazione delle variabili in u e x.	$ln(Cx)=int ^{{y/x}}{frac {dlambda }{F(lambda )-lambda }},!$
Primo ordine, separabile. $yM(xy)+xN(xy),{frac {dy}{dx}}=0,!$ $yM(xy),dx+xN(xy),dy=0,!$	Separazione delle variabili (divisione per xy).	$ln(Cx)=int ^{{xy}}{frac {N(lambda ),dlambda }{lambda [N(lambda )-M(lambda )]}},!$ Se N = M la soluzione è xy = C.
Differenziale esatto, primo ordine. $M(x,y){frac {dy}{dx}}+N(x,y)=0,!$ $M(x,y),dy+N(x,y),dx=0,!$ dove ${frac {partial M}{partial x}}={frac {partial N}{partial y}},!$	Integrazione su tutto il dominio.	${begin{aligned}F(x,y)&=int ^{y}M(x,lambda ),dlambda +int ^{x}N(lambda ,y),dlambda \&+Y(y)+X(x)=Cend{aligned}},!$
Differenziale inesatto, primo ordine. $M(x,y){frac {dy}{dx}}+N(x,y)=0,!$ $M(x,y),dy+N(x,y),dx=0,!$ dove ${frac {partial M}{partial x}}neq {frac {partial N}{partial y}},!$	Il fattore di integrazione μ(x, y) soddisfa: ${frac {partial (mu M)}{partial x}}={frac {partial (mu N)}{partial y}},!$	Se si può trovare μ(x, y): ${begin{aligned}F(x,y)&=int ^{y}mu (x,lambda )M(x,lambda ),dlambda +int ^{x}mu (lambda ,y)N(lambda ,y),dlambda \&+Y(y)+X(x)=C\end{aligned}},!$
Equazioni del secondo ordine
Secondo ordine, autonoma. ${frac {d^{2}y}{dx^{2}}}=F(y),!$	Si moltiplica l'equazione per 2dy/dx e si sostituisce $2{frac {dy}{dx}}{frac {d^{2}y}{dx^{2}}}={frac {d}{dx}}left({frac {dy}{dx}}right)^{2},!$ quindi si integra due volte.	$x=pm int ^{y}{frac {dlambda }{{sqrt {2int ^{lambda }F(varepsilon ),dvarepsilon +C_{1}}}}}+C_{2},!$
Equazioni lineari
Primo ordine, lineare, non omogenea. ${frac {dy}{dx}}+P(x)y=Q(x),!$	Fattore di integrazione: $e^{{int ^{x}P(lambda ),dlambda }}$ .	$y=e^{{-int ^{x}P(lambda ),dlambda }}left[int ^{x}e^{{int ^{lambda }P(varepsilon ),dvarepsilon }}Q(lambda ),{dlambda }+Cright]$
Secondo ordine, lineare, non omogenea, a coefficienti costanti. ${frac {d^{2}y}{dx^{2}}}+b{frac {dy}{dx}}+cy=r(x),!$	Funzione complementare y_c: assumendo y_c = e^αx, si sostituisce e si risolve il polinomio in α, ottenendo le funzioni linearmente indipendenti $e^{{alpha _{j}x}}$ . Integrale particolare y_p: in generale si applica il metodo delle variazioni delle costanti, sebbene talvolta è sufficiente studiare r(x).	$y=y_{c}+y_{p}$ Se b² > 4c, allora: $y_{c}=C_{1}e^{{left(-b+{sqrt {b^{2}-4c}}right){frac {x}{2}}}}+C_{2}e^{{-left(b+{sqrt {b^{2}-4c}}right){frac {x}{2}}}},!$ Se b² = 4c, allora: $y_{c}=(C_{1}x+C_{2})e^{{-bx/2}},!$ Se b² < 4c, allora: $y_{c}=e^{{-b{frac {x}{2}}}}left[C_{1}sin {left({sqrt {left\|b^{2}-4cright\|}}{frac {x}{2}}right)}+C_{2}cos {left({sqrt {left\|b^{2}-4cright\|}}{frac {x}{2}}right)}right],!$
Ordine n, lineare, non omogenea, a coefficienti costanti. $sum _{{j=0}}^{n}b_{j}{frac {d^{j}y}{dx^{j}}}=r(x),!$	Funzione complementare y_c: assumendo y_c = e^αx, si sostituisce e si risolve il polinomio in α, ottenendo le funzioni linearmente indipendenti $e^{{alpha _{j}x}}$ . Integrale particolare y_p: in generale si applica il metodo delle variazioni delle costanti, sebbene talvolta sia sufficiente studiare r(x).	$y=y_{c}+y_{p}$ Dato che α_j sono soluzioni del polinomio di grado n: $prod _{{j=1}}^{n}left(alpha -alpha _{j}right)=0,!$ allora per α_j tutti diversi: $y_{c}=sum _{{j=1}}^{n}C_{j}e^{{alpha _{j}x}},!$ per ogni radice α_j ripetuta k_j volte: $y_{c}=sum _{{j=1}}^{n}left(sum _{{ell =1}}^{{k_{j}}}C_{ell }x^{{ell -1}}right)e^{{alpha _{j}x}},!$ 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: $C_{j}e^{{alpha _{j}x}}=C_{j}e^{{chi _{j}x}}cos(gamma _{j}x+phi _{j}),!$ dove ϕ_j è una costante arbitraria.

Applicazioni |

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 $m$ in caduta libera, sotto l'azione della forza di gravità. Attraverso le leggi descritte da Newton sulla dinamica dei corpi, si ha che:

{mathbf {F}}=m{mathbf {a}}=m{begin{bmatrix}x''\y''\z''end{bmatrix}}

Riferendosi ad un sistema di coordinate cartesiane, con l'asse $y$ 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 $y$ :

my''=-mg

Ciò fornisce $y''=-g$ . Per determinare la $varphi =varphi (x)$ , soluzione del problema, occorre conoscere posizione e velocità iniziale del corpo in un certo istante. Integrando:

int _{{t_{0}}}^{{t}}y''(t),dt=-int _{{t_{0}}}^{{t}}g,dt

da cui:

y'(t)=y'(t_{0})-g(t-t_{0})

e integrando nuovamente:

int _{{t_{0}}}^{{t}}y'(t),dt=int _{{t_{0}}}^{{t}}left(y'(t_{0})-g(t-t_{0})right),dt

da cui si ottiene:

y(t)=y(t_{0})+y'(t_{0})(t-t_{0})-g{frac {(t-t_{0})^{2}}{2}}

Come si vede, la soluzione dipende da due parametri $y'(t_{0})$ e $y(t_{0})$ , rispettivamente velocità e posizione iniziale.

Esempio |

Si consideri il problema di Cauchy:

{begin{cases}y'=y^{{{frac {1}{3}}}}\y(t_{0})=0end{cases}}

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

varphi (t)=0quad forall tin mathbb{R}

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

{frac {dy}{dt}}=y^{{{frac {1}{3}}}}

da cui:

y^{{-{frac {1}{3}}}},dy=dt

e quindi:

int _{0}^{y}y^{{-{frac {1}{3}}}},dy=int _{{t_{0}}}^{t},dtqquad yneq 0

dalla quale si ottiene:

varphi _{pm }^{0}(t)=pm left({frac {2}{3}}(t-t_{0})right)^{{{frac {3}{2}}}}qquad t-t_{0}>0

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

psi _{+}^{0}(t)={begin{cases}varphi _{+}^{0}(t)&t>t_{0}\0&tleq t_{o}end{cases}}qquad psi _{-}^{0}(t)={begin{cases}varphi _{-}^{0}(t)&t>t_{0}\0&tleq t_{o}end{cases}}

che sono anch'esse soluzioni del problema definite sull'insieme $R$ . Inoltre, fissato un $t_{1}>t_{0}$ e posto:

varphi _{pm }^{1}(t)=pm left({frac {2}{3}}(t-t_{1})right)^{{{frac {3}{2}}}}qquad t-t_{1}>0

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

psi _{pm }^{1}(t)={begin{cases}varphi _{pm }^{1}(t)&t>t_{1}\0&tleq t_{1}end{cases}}

definite anch'esse su $R$ . 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 ${df}/{dy}$ non è limitata nel punto $y=0$ : per far fronte a questo problema vi è il teorema di esistenza globale.

Note |

^ Encyclopedia of Mathematics - Cauchy problem, su encyclopediaofmath.org. URL consultato il 06-01-2013.

^ Encyclopedia of Mathematics - Peano theorem, su encyclopediaofmath.org. URL consultato il 06-01-2013.

^ 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

^ Encyclopedia of Mathematics - Kovalevskaya theorem, su encyclopediaofmath.org. URL consultato il 06-01-2013.

Bibliografia |

(EN) Vladimir Igorevich Arnold (1988): Geometrical Methods in the Theory of Ordinary Differential Equations, 2nd ed., Springer, ISBN 0-387-96649-8

(EN) Vladimir Igorevich Arnold (1992): Ordinary Differential Equations, Springer, ISBN 3-540-54813-0

(EN) Martin Braun (1993): Differential Equations and their Applications. An Introduction to Applied Mathematics, 4th ed., Springer, ISBN 0-387-97894-1

(EN) Coddington, Earl A.; Levinson, Norman (1955) Theory of ordinary differential equations, New York, Toronto, London: McGraw-Hill Book Company, Inc. XII, 429 p. (1955).

(EN) D. Zwillinger (1997): Handbook of Differential Equations 3rd ed., Academic Press

(EN) P. L. Sachdev (1997): A Compendium on Nonlinear Ordinary Differential Equations, John Wiley, ISBN 0-471-53134-0

(EN) Po-Fang Hsieh, Yasutaka Sibuya (1999): Basic Theory of Ordinary Differential Equations, Springer, ISBN 0-387-98699-5

(EN) Philip Hartman (2002): Ordinary Differential Equations, 2nd ed., SIAM, ISBN 0-89871-510-5

(EN) A. D. Polyanin, V. F. Zaitsev (2003): Handbook of Exact Solutions for Ordinary Differential Equations, 2nd ed., Chapman & Hall / CRC Press, ISBN 1-58488-297-2

(EN) Refaat El Attar (2005): Ordinary Differential Equations, Lulu Press, ISBN 1-4116-3920-0

Voci correlate |

34-XX sigla della sezione della MSC dedicata alle equazioni differenziali ordinarie.

Condizioni al contorno

Equazione differenziale esatta

Equazione differenziale lineare

Equazione differenziale lineare del secondo ordine

Equazione differenziale lineare di ordine superiore al primo

Funzione lipschitziana

Metodi di soluzione analitica per equazioni differenziali ordinarie

Metodi di soluzione numerica per equazioni differenziali ordinarie

Metodo delle differenze finite

Metodo delle variazioni delle costanti

Problema ai valori iniziali

Problema di Cauchy

Teorema di esistenza di Peano

Teorema di esistenza e unicità per un problema di Cauchy

Wronskiano

Collegamenti esterni |

Equazione differenziale ordinaria, su thes.bncf.firenze.sbn.it, Biblioteca Nazionale Centrale di Firenze.

(EN) Equazione differenziale ordinaria, su Enciclopedia Britannica, Encyclopædia Britannica, Inc.

(EN) E.F. Mishchenko, Differential equation, ordinary, in Encyclopaedia of Mathematics, Springer e European Mathematical Society, 2002.

$Matematica$ Portale Matematica: accedi alle voci di Wikipedia che trattano di matematica

[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.

[EDEBVP-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.

搜尋此網誌