În matematică , o ecuație diferențială obișnuită (prescurtată în EDO , sau ODE din acronimul în engleză Ordinary Differential Equation ) este o ecuație diferențială care implică o funcție a unei variabile și derivatele acesteia de orice ordin: este un obiect matematic larg utilizat în fizică și multe alte domenii ale științei ; de exemplu un sistem dinamic este descris printr-o ecuație diferențială obișnuită.
Ca și în cazul tuturor ecuațiilor diferențiale, de obicei nu este posibil să se rezolve exact un EDO și, în orice caz, nu există metode generale pentru ao face. Prin urmare, diferitele cazuri posibile sunt analizate individual și adesea ne limităm la studierea comportamentului calitativ al soluției fără a fi posibilă obținerea unei expresii analitice. Ecuațiile liniare (de orice ordin) sunt deosebit de simple, deoarece pot fi întotdeauna urmărite înapoi la un sistem de ecuații liniare de primul ordin.
Este {\ displaystyle F \ colon \ Omega \ subseteq \ mathbb {C} ^ {n + 2} \ rightarrow \ mathbb {C}} , cu{\ displaystyle \ Omega \ neq \ varnothing} un întregdeschis și conectat e {\ displaystyle n \ in \ mathbb {N}} .
Se numește ecuație diferențială obișnuită de ordine {\ displaystyle n} o relatie de genul:
{\ displaystyle F \ left (x, u (x), u '(x), \ ldots, u ^ {(n)} (x) \ right) = 0,}
unde cu {\ displaystyle u ^ {(i)} (x)} este indicată derivata{\ displaystyle i} -alea a funcției{\ displaystyle u (x)} .
De sine {\ displaystyle F} este definit într-o regiune {\ displaystyle \ Omega} a spațiului euclidian{\ displaystyle \ mathbb {R} ^ {n + 2}} , atunci considerăm mai corect ecuații diferențiale obișnuite în câmpul real, cu valori reale dacă {\ displaystyle F} este la valori în {\ displaystyle \ mathbb {R}} .
Ordinea unei ecuații este ordinea maximă de derivare care apare în ea, în timp ce adjectivul obișnuit se referă la faptul că necunoscutul este o funcție a unei singure variabile. Dacă necunoscutul depinde de mai multe variabile, avem o ecuație diferențială parțială .
Este {\ displaystyle I} un interval de {\ displaystyle \ mathbb {R}} . O funcție este definită ca soluția sau integrala ecuației diferențiale obișnuite {\ displaystyle \ varphi = \ varphi (x)} astfel încât:
{\ displaystyle \ varphi (x) \ în C ^ {n} (I) \ qquad F \ left (x, \ varphi (x), \ varphi '(x), \ ldots, \ varphi ^ {(n)} (x) \ right) = 0 \ qquad \ forall x \ in I.}
O ecuație diferențială obișnuită se spune că este autonomă dacă {\ displaystyle F} nu depinde în mod explicit de {\ displaystyle x} .
Se spune că o ecuație diferențială obișnuită este scrisă în formă normală dacă poate fi explicitată cu respect {\ displaystyle u ^ {(n)} (x)} :
{\ displaystyle u ^ {(n)} (x) = G \ left (x, u, u ', \ ldots, u ^ {(n-1)} \ right).}
Se mai spune că este liniar dacă {\ displaystyle F} este o combinație liniară de {\ displaystyle u, u ', \ ldots, u ^ {(n)}} , sau:
{\ displaystyle F \ left (x, u, u ', \ ldots, u ^ {(n)} \ right) = s (x) + b_ {0} (x) u + b_ {1} (x) u '+ \ ldots + b_ {n} (x) u ^ {(n)},}
sau, echivalent:
{\ displaystyle u ^ {(n)} = \ sum _ {i = 0} ^ {n-1} a_ {i} (x) u ^ {(i)} + r (x),}
unde este:
{\ displaystyle r (x), a_ {0} (x), a_ {1} (x), \ ldots, a_ {n-1} (x) \ în C ^ {0} (I).}
Termenul {\ displaystyle r (x)} se numește sursă sau forțare și, dacă este zero, se spune că ecuația diferențială liniară este omogenă .
O ecuație obișnuită are soluții liniar independente într-un număr egal cu gradul ecuației și fiecare combinație liniară a acestora este la rândul său o soluție.
Având în vedere o ecuație diferențială obișnuită, dacă este cunoscută o soluție generală a ecuației omogene asociate acesteia, atunci este posibil să se găsească o soluție specială a ecuației „complete”. În acest scop, există mai multe proceduri , inclusiv metoda variațiilor constante și utilizarea transformatei Laplace . Pentru cazurile mai simple există și câteva teorii: de exemplu, pentru ecuațiile de gradul I este posibil să se caute un factor de integrare adecvat, pentru al doilea există teoria Sturm-Liouville . Cu toate acestea, în general, de obicei, singura modalitate posibilă de a studia soluția este utilizarea unei metode de soluție numerică .
Sisteme ODE
Reducerea la un sistem de ecuații de ordinul 1
Un sistem de ecuații ordinare diferențiale ordinare {\ displaystyle 1} în formă normală este o relație vectorială de tipul:
O soluție clasică a unui astfel de sistem este o funcție {\ displaystyle \ mathbf {u} \ colon I \ to \ mathbb {R} ^ {n}} astfel încât:
{\ displaystyle \ mathbf {u} (x) \ în C ^ {1} (I) \ qquad \ mathbf {u} '(x) = \ mathbf {f} \ left (x, \ mathbf {u} (x ) \ right) \ qquad \ forall x \ in I.}
O relevanță deosebită în scopuri practice este reducerea unei ecuații diferențiale obișnuite de ordine {\ displaystyle n} în formă normală la un sistem diferențial de prim ordin. Această tehnică face posibilă simplificarea considerabilă a unor tipuri de probleme, evitând introducerea unor forme complexe de rezoluție. Este:
{\ displaystyle u ^ {(n)} (x) = G \ left (x, u, u ', \ ldots, u ^ {(n-1)} \ right)}
o ecuație diferențială de ordine {\ displaystyle n} tip normal. Ei se definesc:
{\ displaystyle u_ {i} = u ^ {(i-1)} (x) \ qquad \ mathbf {u} = \ left (u_ {i} \ right), \ qquad i \ in \ {1, \ ldots , n \},}
astfel încât {\ displaystyle u_ {i} '= u_ {i + 1}} , și în special {\ displaystyle u_ {n} '= u ^ {(n)} (x)} . Ecuația diferențială este deci echivalentă cu sistemul:
{\ displaystyle {\ begin {cases} u_ {1} '= u_ {2}, \\ u_ {2}' = u_ {3}, \\\ vdots \\ u_ {n-1} '= u_ {n }, \\ u_ {n} '= G \ left (x, u_ {1}, u_ {2}, \ ldots, u_ {n} \ right) = G \ left (x, \ mathbf {u} \ right ). \ end {cases}}}
Prin plasarea:
{\ displaystyle \ mathbf {f} \ left (x, \ mathbf {u} \ right) = {\ begin {pmatrix} u_ {2} \\ u_ {3} \\\ vdots \\ u_ {n} \\ G \ left (x, \ mathbf {u} \ right) \ end {pmatrix}},}
adică, cu alte cuvinte, este întotdeauna posibil să se traducă totul într-o ecuație de ordine 1. Cu o procedură complet similară cu cea urmată, este posibil și invers, adică să se obțină o ecuație de ordine {\ displaystyle n} plecând de la o ecuație de ordinul 1 în care un vector {\ displaystyle \ mathbf {u}} in marime {\ displaystyle n} ambele soluție.
Sistem de ecuații de ordine n
Dacă luăm în considerare un transportator {\ displaystyle \ mathbf {y}} definit ca:
{\ displaystyle F \ left (x, y, y ', \ cdots, y ^ {(n)} \ right) = 0}
o functie {\ displaystyle u \ colon I \ subset \ mathbb {R} \ to \ mathbb {R}} se numește soluția (sau integrala ) ecuației diferențiale obișnuite dacă {\ displaystyle u} este diferențiat{\ displaystyle n} timpul a expirat {\ displaystyle I} și avem:
{\ displaystyle F (x, u, u ', \ \ cdots, \ u ^ {(n)}) = 0, \ quad x \ in I.}
Oferiți două soluții {\ displaystyle u \ colon J \ subset \ mathbb {R} \ to \ mathbb {R}} Și {\ displaystyle v \ colon I \ subset \ mathbb {R} \ to \ mathbb {R}} , {\ displaystyle u} se numește o extensie a {\ displaystyle v} de sine {\ displaystyle I \ subset J} Și:
{\ displaystyle u (x) = v (x), \ quad x \ in I.}
O soluție care nu are extensii se numește o soluție maximă , în timp ce o soluție este definită pe toate {\ displaystyle \ mathbb {R}} se numește o soluție globală .
O soluție generală a unei ecuații de ordine {\ displaystyle n} este o soluție care conține {\ displaystyle n} constante de integrare independente, în timp ce o soluție particulară este obținută din soluția generală prin acordarea unei valori fixe constantelor, de obicei pentru a satisface condițiileinițiale sau condițiile limită . În acest context, o soluție singulară este o soluție care nu poate fi obținută prin atribuirea unei valori definite constantelor de integrare.
O problemă de valoare inițială este o ecuație diferențială obișnuită: [1]
{\ displaystyle y '(t) = f (t, y (t)), \ qquad f: \ Omega \ subset \ mathbb {R} \ times \ mathbb {R} ^ {n},}
asociat cu un punct din domeniul {\ displaystyle f} :
{\ displaystyle (t_ {0}, y_ {0}) \ în \ Omega}
numită condiție inițială . Soluția unei probleme de valoare inițială este deci o funcție {\ displaystyle y} care este soluția ecuației diferențiale și satisface condiția {\ displaystyle y (t_ {0}) = y_ {0}} . Cu alte cuvinte, problema lui Cauchy este de a găsi o curbă {\ displaystyle y} , dintre cele definite de {\ displaystyle y '(t) = f (t, y (t))} , care trece pentru punct {\ displaystyle y (t_ {0}) = y_ {0}} .
Existența locală a unei soluții a fost dovedită de Augustin-Louis Cauchy sub ipoteza continuității și mărginirii lui {\ displaystyle f} într-o regiune a domeniului său ( teorema Peano sau Cauchy-Peano). [2] Ulterior, existența locală și unicitatea au fost arătate de Émile Picard cu ipoteza lipsețitzianității în ceea ce privește {\ displaystyle y} ( teorema existenței și unicității pentru o problemă Cauchy sau Picard-Lindelöf), iar acest rezultat poate fi extins la o formă globală. Daca atunci {\ displaystyle f} este Lipschitz într-o regiune {\ displaystyle D} din domeniu, atunci există cel puțin o soluție-curbă care poate fi diferențiată cu continuitatea care trece prin fiecare punct interior a {\ displaystyle D} . Simplificarea întrebării, dacă {\ displaystyle f} Și {\ displaystyle \ partial f / \ partial y} sunt continue într-un dreptunghi închis în plan {\ displaystyle xy} a formei:
{\ displaystyle R = [x_ {0} -a, x_ {0} + a] \ times [y_ {0} -b, y_ {0} + b],}
unde este {\ displaystyle a, b \ in \ mathbb {R}} Și {\ displaystyle \ times} este produsul cartezian , apoi există un interval:
în care pentru unii {\ displaystyle h \ in \ mathbb {R}} singura soluție poate fi găsită. [3] Acest rezultat se aplică și ecuațiilor neliniare ale formei {\ displaystyle F (x, y)} precum și sisteme de ecuații.
Există mai multe teoreme care permit stabilirea existenței și posibilei unicități locale a soluțiilor pentru problemele inițiale date. Cele două principale sunt teorema existenței și unicității pentru o problemă Cauchy, care presupune Lipschitzianitatea funcției care definește ecuația obișnuită și încheie existența și unicitatea sa locală, și teorema existenței lui Peano, care presupune continuitatea, nu face decât să concluzioneze existența. Prima afirmă că a dat problema valorilor inițiale:
{\ displaystyle y '(t) = f (t, y (t)), \ qquad y (t_ {0}) = y_ {0}, \ qquad t \ in [t_ {0} - \ varepsilon, t_ { 0} + \ varepsilon],}
de sine {\ displaystyle f} este o funcție Lipschitz în {\ displaystyle y} și continuă în {\ displaystyle t} apoi pentru unii {\ displaystyle \ epsilon> 0} există o singură soluție {\ displaystyle y (t)} la problema valorii inițiale pe interval {\ displaystyle [t_ {0} - \ epsilon, t_ {0} + \ epsilon]} . Teorema existenței lui Peano presupune {\ displaystyle F (x, y)} continuă doar, dar nu garantează unicitatea soluției, arătând doar existența sa locală. Dacă există, este fie unic, fie există infinite . O extensie a acestui rezultat seobține cu teorema existenței lui Carathéodory , care se aplică și cazurilor în care ecuația nu este continuă.
Unicitate globală
Dacă ipotezele existenței și teoremei unicității pentru o problemă Cauchy sunt satisfăcute, atunci condiția existenței locale poate fi extinsă la un rezultat global. Pentru fiecare condiție inițială{\ displaystyle (x_ {0}, y_ {0})} există un singur interval maxim deschis (posibil nelimitat):
{\ displaystyle I _ {\ max} = (x _ {-}, x _ {+}), \ qquad x _ {\ pm} \ in \ mathbb {R} \ cup \ {\ pm \ infty \}, \ quad x_ {0} \ în I _ {\ max},}
astfel încât fiecare soluție care îndeplinește condiția inițială este o restricție a soluției care îndeplinește acea condiție inițială care este definită pe domeniu {\ displaystyle I _ {\ max}} . În cazul în care{\ displaystyle x _ {\ pm} \ neq \ pm \ infty} , există doar două posibilități:
{\ displaystyle \ limsup _ {x \ to x _ {\ pm}} \ | y (x) \ | \ rightarrow + \ infty,}
{\ displaystyle \ lim _ {x \ to x _ {\ pm}} y (x) \ in \ partial {\ bar {\ Omega}},}
unde este {\ displaystyle \ Omega} este deschiderea în care {\ displaystyle F} este definit și {\ displaystyle \ partial {\ bar {\ Omega}}} este frontiera sa.
Pentru ca soluția să fie unică, domeniul maxim în care este definită soluția trebuie să fie un interval, care în general depinde de starea inițială. De exemplu, luați în considerare:
{\ displaystyle y '= y ^ {2}.}
De cand {\ displaystyle F (x, y) = y ^ {2}} este Lipschitz, satisface teorema Picard - Lindelöf. Soluția
cu {\ displaystyle f} o funcție continuă definită pe un set deschis de {\ displaystyle \ mathbb {R}} . Problema revine la căutarea primitivelor{\ displaystyle y (t)} din {\ displaystyle f (t)} . În acest caz, dacă {\ displaystyle y (t)} este o soluție și atunci {\ displaystyle y (t) + c} , cu {\ displaystyle c \ in \ mathbb {R}} , este încă soluție. De asemenea, dacă {\ displaystyle y_ {1} (t)} Și {\ displaystyle y_ {2} (t)} sunt soluții pe care le avem{\ displaystyle y_ {2} = y_ {1} + c} pentru unii {\ displaystyle c \ in \ mathbb {R}} .
Având în vedere o condiție inițială {\ displaystyle y (t_ {0}) = y_ {0}} , scriind ecuația în formă {\ displaystyle dy / dt = f (t)} , astfel încât {\ displaystyle dy = f (t) dt} , avem:
{\ displaystyle \ int _ {y_ {0}} ^ {y} d \ xi = \ int _ {t_ {0}} ^ {t} f (\ tau) d \ tau}
{\ displaystyle y (t) = y_ {0} + \ int _ {t_ {0}} ^ {t} f (\ tau) d \ tau.}
Soluția este determinată în mod unic de datele inițiale {\ displaystyle y (t_ {0}) = y_ {0}} : dat orice soluție {\ displaystyle y_ {1} (t)} , da {\ displaystyle y_ {1} (t_ {0}) + c = y_ {0}} prin urmare {\ displaystyle c = y_ {0} -y_ {1} (t_ {0})} .
O ecuație autonomă este o ecuație diferențială obișnuită de tipul:
{\ displaystyle y '(t) = f (y (t)),}
unde este {\ displaystyle f} este o funcție continuă cu prima derivată continuă pe un interval {\ displaystyle I \ subset \ mathbb {R}} , și care nu depinde de variabila independentă {\ displaystyle t} . De sine {\ displaystyle y} este un vector al {\ displaystyle \ mathbb {R} ^ {n}} avem un sistem autonom, adică un sistem de ecuații diferențiale ordinare autonome:
Punctele sunt de o importanță deosebită {\ displaystyle x_ {0}} astfel pentru care {\ displaystyle f (x_ {0}) = 0} , numite puncte de echilibru , cărora le corespunde soluția constantă {\ displaystyle y = x_ {0}} .
Un sistem generic de ecuații diferențiale obișnuite (în care {\ displaystyle f} a depinde de {\ displaystyle t} ):
{\ displaystyle {\ frac {d} {dt}} y (t) = f (y (t), t)}
poate fi autonomizat prin introducerea unei noi necunoscute {\ displaystyle y_ {n + 1} = t} .
Luați în considerare un set care este pur și simplu conectat și deschis{\ displaystyle D \ subset \ mathbb {R} ^ {2}} și două funcții {\ displaystyle I} Și {\ displaystyle J}continuă mai departe {\ displaystyle D} . Ecuația diferențială implicită
{\ displaystyle I (x, y) \, \ mathrm {d} x + J (x, y) \, \ mathrm {d} y = 0}
Termenul „exact” se referă la derivata totală a unei funcții, uneori numită „derivată exactă”, aceea pentru o funcție {\ displaystyle F (x_ {0}, x_ {1}, \ dots, x_ {n-1}, x_ {n})} este dat în {\ displaystyle x_ {0}} din:
În aplicațiile fizice{\ displaystyle I} Și {\ displaystyle J} ele sunt de obicei nu numai continue, ci și continuu diferențiate, iar teorema lui Schwarz oferă apoi o condiție necesară și suficientă pentru existența funcției potențiale {\ displaystyle F} (pentru ecuațiile definite pe un set care nu este pur și simplu conectat, acest criteriu este necesar doar). Există dacă și numai dacă:
Iată câteva cazuri importante de ecuații diferențiale obișnuite care pot fi rezolvate exact.
În tabel funcțiile {\ displaystyle P (x)} , {\ displaystyle Q (x)} , {\ displaystyle P (y)} , {\ displaystyle Q (y)} , {\ displaystyle M (x, y)} Și {\ displaystyle N (x, y)} sunt funcții care pot fi integrate în {\ displaystyle x} , {\ displaystyle y} , in timp ce {\ displaystyle b} Și {\ displaystyle c} li se dau constante reale. În plus, {\ displaystyle C_ {1}, C_ {2}, \ dots} sunt constante arbitrare, în general complexe. Notatia {\ displaystyle \ int ^ {x} F (\ lambda) d \ lambda} indică integrarea {\ displaystyle F (\ lambda)} în comparație cu {\ displaystyle \ lambda} și înlocuirea ulterioară {\ displaystyle \ lambda = x} .
Si moltiplica l'equazione per 2 dy/dx e si sostituisce {\displaystyle 2{\frac {dy}{dx}}{\frac {d^{2}y}{dx^{2}}}={\frac {d}{dx}}\left({\frac {dy}{dx}}\right)^{2}\,\!} quindi si integra due volte.
Funzione complementare y c : assumendo y c = eα x , si sostituisce e si risolve il polinomio in α, ottenendo le funzioni linearmente indipendenti {\displaystyle e^{\alpha _{j}x}} .
Funzione complementare y c : assumendo y c = eα x , si sostituisce e si risolve il polinomio in α, ottenendo le funzioni linearmente indipendenti {\displaystyle 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 ).
{\displaystyle y=y_{c}+y_{p}}
Dato che α j sono soluzioni del polinomio di grado n :
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:
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{\displaystyle m} 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 {\displaystyle 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 {\displaystyle y} :
{\displaystyle my''=-mg.}
Ciò fornisce {\displaystyle y''=-g} . Per determinare la {\displaystyle \varphi =\varphi (x)} , soluzione del problema, occorre conoscere posizione e velocità iniziale del corpo in un certo istante. Integrando:
Come si vede, la soluzione dipende da due parametri{\displaystyle y'(t_{0})} e {\displaystyle y(t_{0})} , rispettivamente velocità e posizione iniziale.
Il teorema di esistenza locale garantisce l'esistenza di almeno una soluzione {\displaystyle \varphi (t)} . Una di queste è banale ed è anche globale, cioè definita su tutto l'insieme {\displaystyle \mathbb {R} } :
Oltre a questa soluzione è tuttavia possibile trovarne un'altra integrando l'equazione, trattandosi di un'equazione differenziale a variabili separabili. Perciò:
che sono anch'esse soluzioni del problema definite sull'insieme {\displaystyle \mathbb {R} } . Inoltre, fissato un {\displaystyle t_{1}>t_{0}} e posto:
definite anch'esse su {\displaystyle \mathbb {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{\displaystyle {df}/{dy}} non è limitata nel punto {\displaystyle y=0} : per far fronte a questo problema vi è il teorema di esistenza globale.
^Elementary Differential Equations and Boundary Value Problems (4th Edition), WE Boyce, RC Diprima, Wiley International, John Wiley & Sons, 1986, ISBN 0-471-83824-1
( 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 ) PL 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 ) AD Polyanin, VF 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.