Ecuatii diferentiale ˆın form˘a simetric˘a. Integrale prime
Ecuatii diferentiale ˆın form˘a simetric˘a. Integrale prime
Ecuatii diferentiale ˆın form˘a simetric˘a. Integrale prime
Create successful ePaper yourself
Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.
Ecuat¸ii diferent¸iale în formă simetrică. <strong>Integrale</strong> <strong>prime</strong><br />
Considerăm o relat¸ie de forma<br />
(1) P (x, y)dx + Q(x, y)dy = 0<br />
unde P, Q : D → R sunt funct¸ii continue pe o submult¸ime nevidă D ⊂ R 2 . Aici<br />
notat¸iile dx, respectiv dy reprezintă diferent¸ialele funct¸iilor de două variable (x, y) ↦→<br />
x, respectiv (x, y) ↦→ y. Dacă notăm<br />
ω = P (x, y)dx + Q(x, y)dy<br />
spunem că ω este o formă diferent¸ială de ordinul întâi sau, mai simplu, este o 1-<br />
formă.<br />
Fără a preciza deocamdată ,,necunoscuta”, numim (1) ecuat¸ie Pfaff sau ecuat¸ie<br />
diferent¸ială în formă simetrică. In legătură cu (1) dăm următoarele definit¸ii.<br />
Definit¸ia 1 Fie U ⊂ R 2 o mult¸ime nevidă, deschisă ¸si conexă. Dacă există o funct¸ie<br />
continuă µ : U → R ¸si o funct¸ie de clasă C 1 neconstantă H : U → R astfel încât<br />
dH = µω în U<br />
atunci spunem că (1) este o ecuat¸ie integrabilă în U.<br />
Funct¸ia H se nume¸ste integrală primă în U a ecuat¸iei (1), iar µ se nume¸ste<br />
factor integrant corespunzător integralei <strong>prime</strong> H.<br />
In cazul în care există o funct¸ie de clasă C 1 neconstantă H : U → R astfel încât<br />
dH = ω în U<br />
atunci spunem în plus că ω este o 1-formă exactă în U, ¸si că (1) este o ecuat¸ie<br />
exactă în U.<br />
Mult¸imile de nivel ale unei integrale <strong>prime</strong> H în U pentru (1) se numesc curbe<br />
integrale în U ale ecuat¸iei diferent¸iale în formă simetrică (1).<br />
1
Observat¸ii. 1. Relat¸ia<br />
este echivalentă cu<br />
∂H<br />
(x, y) = µ(x, y)P (x, y) ¸si<br />
∂x<br />
dH = µω în U<br />
∂H<br />
(x, y) = µ(x, y)Q(x, y) pentru orice (x, y) ∈ U.<br />
∂y<br />
2. Deoarece curbele de nivel ale unei integrale <strong>prime</strong> H în U a ecuat¸iei (1) se numesc<br />
curbe integrale, mai spunem că solut¸ia generală în U a ecuat¸iei (1) este dată de<br />
formula<br />
Exemplul 1. Considerăm ecuat¸ia Pfaff<br />
H(x, y) = c, c ∈ R constantă arbitrară.<br />
xdx + ydy = 0.<br />
Este u¸sor de văzut că funct¸ia H : R 2 → R dată prin H(x, y) = (x 2 + y 2 )/2 are<br />
∂H/∂x = x ¸si ∂H/∂y = y, adică dH = xdx + ydy. Prin urmare ecuat¸ia dată este<br />
integrabilă în R 2 , chiar exactă de fapt. Funct¸ia H este o integrală primă în R 2 a<br />
ecuat¸iei date, iar solut¸ia ei generală este dată de formula<br />
x 2 + y 2 = c, c ∈ R constantă arbitrară.<br />
Familia curbelor integrale ale ecuat¸iei Pfaff xdx + ydy = 0 este familia de cercuri cu<br />
centrul in origine de rază arbitrară √ c cu c > 0.<br />
Se poate arăta că ¸si funct¸iile ˜ H, ˆ H, ˇ H : R 2 → R date prin ˜ H(x, y) = x 2 + y 2 ,<br />
ˆH(x, y) = x 2 + y 2 + 10 respectiv ˇ H(x, y) = e x2 +y 2<br />
sunt integrale <strong>prime</strong> în R 2 ale<br />
acestei ecuat¸ii. Dar să observăm că acestea se pot scrie ,,în funct¸ie” de H, adică<br />
notând ˜ F , ˆ F , ˇ F : R → R funct¸iile date prin ˜ F (h) = 2h, ˆ F (h) = 2h + 10 respectiv<br />
ˇF (h) = e 2h avem ˜ H = ˜ F (H), ˆ H = ˆ F (H) respectiv ˇ H = ˇ F (H). Observăm de<br />
asemenea că toate aceste integrale <strong>prime</strong> au acelea¸si curbe de nivel ca ¸si H.<br />
Exercit¸iul 1. Arătat¸i că ecuat¸ia Pfaff<br />
ydx + xdy = 0<br />
este exactă ¸si reprezentat¸i curbele ei integrale.<br />
2
Dăm în continuare un rezultat important ¸si foarte util în a stabili dacă o ecuat¸ie<br />
este sau nu exactă.<br />
Teorema 1 Fie I, J ⊂ R intervale nevide deschise ¸si notăm U = I × J. Fie P, Q ∈<br />
C 1 (U). Ecuat¸ia (1) is exactă în U dacă ¸si numai dacă<br />
(2)<br />
∂P<br />
∂y<br />
∂Q<br />
(x, y) = (x, y) pentru orice (x, y) ∈ U.<br />
∂x<br />
Exemplul 2. Fie I, J ⊂ R intervale nevide deschise, (x0, y0) ∈ I × J fixat ¸si<br />
P ∈ C 1 (I), Q ∈ C 1 (J). Ecuat¸ia Pfaff<br />
este exactă în I × J cu integrala primă<br />
H(x, y) =<br />
P (x)dx + Q(y)dy = 0<br />
x<br />
x0<br />
P (s)ds +<br />
Exemplul 3. Ne propunem să arătăm că ecuat¸ia<br />
y<br />
y0<br />
Q(s)ds.<br />
(x 2 + y)dx + (x − 3y)dy = 0<br />
este exactă în R 2 ¸si să determinăm o integrală primă.<br />
Notând P (x, y) = x 2 +y ¸si Q(x, y) = x−3y avem ∂P/∂y = 1 ¸si ∂Q/∂x = 1. Pe baza<br />
Teoremei 1 deducem că ecuat¸ia este exactă în R 2 . Prin urmare, există H : R 2 → R<br />
astfel încât dH = (x 2 +y)dx+(x−3y)dy, adică ∂H/∂x = x 2 +y ¸si ∂H/∂y = x−3y.<br />
Integrând relat¸ia ∂H/∂x = x 2 +y în raport cu variabila x avem că există ϕ ∈ C 1 (R)<br />
astfel încât H(x, y) = x 3 /3+xy+ϕ(y). Inlocuind această expresie în ∂H/∂y = x−3y<br />
obt¸inem x + ϕ ′ (y) = x − 3y, adică ϕ ′ (y) = −3y. Alegem ϕ(y) = −3y 2 /2 ¸si obt¸inem<br />
că<br />
H(x, y) = x 3 /3 + xy − 3y 2 /2<br />
este o integrală primă în R 2 a ecuat¸iei date.<br />
3
Dăm în continuare o caracterizare a not¸iunii de integrală primă.<br />
Teorema 2 Fie U ⊂ R 2 o mult¸ime nevidă, deschisă ¸si conexă (de exemplu U = R 2 ),<br />
¸si P, Q ∈ C(U) astfel încât (P (x, y), Q(x, y)) = (0, 0) pentru orice (x, y) ∈ U.<br />
Fie H ∈ C 1 (U). Avem că H este o integrală primă în U a ecuat¸iei Pfaff<br />
(1) P (x, y)dx + Q(x, y)dy = 0<br />
dacă ¸si numai dacă<br />
(3) Q(x, y) ∂H<br />
∂x<br />
(x, y) − P (x, y)∂H (x, y) = 0 pentru orice (x, y) ∈ U.<br />
∂y<br />
Demonstrat¸ie. Dacă H e o integrală primă în U a ecuat¸iei Pfaff (1) atunci (a se<br />
vedea ¸si Observat¸ia 1) există o funct¸ie continuă µ : U → R astfel încât ∂H/∂x = µP<br />
¸si ∂H/∂y = µQ. De aici imediat se poate observa că (3) are loc.<br />
Dacă H satisface (3) atunci avem că pentru fiecare (x, y) ∈ U vectorii (Q, −P )<br />
¸si ∇H = (∂H/∂x, ∂H/∂y) sunt ortogonali. Prin urmare vectorii ∇H ¸si (P, Q)<br />
sunt paraleli. Atunci există µ astfel încât ∇H = µ(P, Q). Să ne reamintim<br />
că acest rat¸ionament l-am făcut pentru fiecare (x, y) ∈ U. Prin urmare am<br />
obt¸inut funct¸ia µ : U → R despre care ne propunem să demonstrăm că este<br />
continuă. Pentru aceasta definim mult¸imile U1 = {(x, y) ∈ U : P (x, y) = 0}<br />
¸si U2 = {(x, y) ∈ U : Q(x, y) = 0}. Acestea sunt mult¸imi deschise (deoarece<br />
P ¸si Q sunt continue) care satisfac relat¸ia U1 ∪ U2 = U, însă există posibili-<br />
tatea pentru fiecare dintre ele să coincidă cu U sau chiar să fie vide. Definim<br />
funct¸iile µ1 : U1 → R ¸si µ2 : U2 → R prin µ1(x, y) = ∂H<br />
(x, y)/P (x, y), respectiv<br />
∂x<br />
µ2(x, y) = ∂H<br />
∂y (x, y)/Q(x, y). Este u¸sor de văzut că µ1 ¸si µ2 sunt continue. Deoarece<br />
∇H = µ(P, Q) în U, avem că µ1 = µ|U1 ¸si µ2 = µ|U2, ceea ce, împreună cu faptul că<br />
U1 ∪ U2 = U, ne asigură continuitatea lui µ pe U. Deci, într-adevăr, există o funct¸ie<br />
continuă µ : U → R astfel încât dH = µP dx + µQdy, adică H e o integrală primă<br />
în U a ecuat¸iei (1). <br />
Observat¸ia 3. Un punct (x0, y0) ∈ U astfel încât (P (x0, y0), Q(x0, y0)) = (0, 0) se<br />
nume¸ste punct singular pentru ecuat¸ia Pfaff (1). Vă invit să descoperit¸i unde s-a<br />
folosit în demonstrat¸ia Teoremei anterioare că (1) nu are puncte singulare în U.<br />
4
Revenim la ecuat¸iile diferent¸iale de forma<br />
(4) y ′ = f(x, y)<br />
unde f : D → R este o funct¸ie continuă pe o submult¸ime nevidă D ⊂ R 2 . Ne<br />
propunem mai întâi să dăm not¸iunea de integrală primă pentru ecuat¸ia (4).<br />
Definit¸ia 2 Fie U ⊂ R 2 o mult¸ime nevidă, deschisă ¸si conexă ¸si H : U → R o<br />
funct¸ie de clasă C 1 neconstantă. Dacă oricare ar fi ϕ : I → R o solut¸ie în U pentru<br />
(4) există o constantă c ∈ R astfel încât<br />
H(x, ϕ(x)) = c pentru orice x ∈ I<br />
atunci spunem că H este o integrală primă în U a ecuat¸iei (4).<br />
Exemplul 4. H(x, y) = ye −x este o integrală primă în R 2 a ecuat¸iei y ′ = y. Stim<br />
că solut¸ia generală a acestei ecuat¸ii este y = ce x . Atunci oricare ar fi ϕ : R → R o<br />
solut¸ie există c ∈ R astfel încât ϕ(x) = ce x pentru orice x ∈ R. Mai mult, evident<br />
avem H(x, ϕ(x)) = c pentru orice x ∈ R.<br />
Considerăm acum ecuat¸ia diferent¸ială în formă simetrică<br />
(5) f(x, y)dx − dy = 0<br />
care se obt¸ine formal din (4) înlocuind y ′ = dy/dx.<br />
Următorul rezultat arată că ecuat¸iile diferent¸iale (4) ¸si (5) sunt geometric echiva-<br />
lente, adică au acelea¸si curbe integrale.<br />
Teorema 3 Fie U ⊂ R 2 o mult¸ime nevidă, deschisă ¸si conexă ce cont¸ine doar<br />
puncte de existent¸ă ¸si unicitate pentru ecuat¸ia (4). Fie H ∈ C 1 (U). Avem că H este<br />
o integrală primă în U pentru (4) dacă ¸si numai dacă H este o integrală primă în<br />
U pentru (5).<br />
Demonstrat¸ie. Mai întâi presupunem că H este o integrală primă în U pentru (4).<br />
Fie (x0, y0) ∈ U arbitrar fixat ¸si ϕ : I → R unica solut¸ie (în U) a Problemei Cauchy<br />
y ′ = f(x, y), y(x0) = y0.<br />
5
Atunci există c ∈ R astfel încât H(x, ϕ(x)) = c pentru orice x ∈ I. Derivând<br />
această relat¸ie obt¸inem<br />
∂H<br />
∂H<br />
(x, ϕ(x)) +<br />
∂x ∂y (x, ϕ(x))ϕ′ (x) = 0 pentru orice x ∈ I.<br />
In particular pentru x = x0 ¸si t¸inând cont că ϕ(x0) = y0 ¸si ϕ ′ (x0) = f(x0, y0) avem<br />
∂H<br />
∂x (x0, y0) + ∂H<br />
∂y (x0, y0)f(x0, y0) = 0.<br />
Având în vedere că (x0, y0) ∈ U este un punct arbitrar fixat deducem că<br />
(6)<br />
∂H<br />
∂x<br />
∂H<br />
(x, y) + (x, y)f(x, y) = 0 pentru orice (x, y) ∈ U.<br />
∂y<br />
De aici, pe baza Teoremei 2, avem că H este integrală primă în U a ecuat¸iei (5).<br />
Acum presupunem că H este o integrală primă în U pentru (5). Tot pe baza<br />
Teoremei 2 avem că (6) are loc. Fie ϕ : I → R o solut¸ie în U a ecuat¸iei (4). Avem<br />
d<br />
∂H<br />
∂H<br />
H(x, ϕ(x)) = (x, ϕ(x)) +<br />
dx ∂x ∂y (x, ϕ(x))ϕ′ (x) pentru orice x ∈ I.<br />
Deoarece ϕ ′ (x) = f(x, ϕ(x)) pentru orice x ∈ I, folosind (6) obt¸inem că<br />
d<br />
H(x, ϕ(x)) = 0 pentru orice x ∈ I,<br />
dx<br />
adică H(x, ϕ(x)) este o funct¸ie constantă pe I. Cu aceasta demonstrat¸ia este<br />
încheiată. <br />
6