20.07.2013 Views

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

SHOW MORE
SHOW LESS

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

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!