Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                

EDO Exactas

Descargar como pdf o txt
Descargar como pdf o txt
Está en la página 1de 12

Ecuaciones Diferenciales Exactas

Coordinación de Ecuaciones Diferenciales y Métodos


Numéricos, DMCC

Ecuaciones Diferenciales Exactas.


Ecuaciones que se reducen a EDO exactas: Factor integrante.

DMCC, Facultad de Ciencia, USACH


Diferencial Total

Definición: Si f = f (x, y ) es una función con primeras derivadas continuas en una


región del plano XY , se define su Diferencial Total o Exacto por:

∂f ∂f
df = dx + dy . (1)
∂x ∂y
Observación:
1 Se deduce que si:
df = 0 ⇒ f (x, y ) = C

2 Dada una familia de curvas f (x, y ) = C , la ecuación diferencial asociada a esta


familia es una ecuación de la forma diferencial total igual cero, es decir: df = 0

2/12
Ecuaciones Diferenciales Exactas

Consideremos ecuaciones diferenciales ordinarias de la forma

M(x, y )dx + N(x, y )dy = 0 . (2)

Definición: Diremos que la ecuación (2) es una ecuación diferencial exacta en D ⊂ R2 ,


si existe una función f (x, y ) tal que

∂f ∂f
(x, y ) = M(x, y ) y (x, y ) = N(x, y ) ∀(x, y ) ∈ D .
∂x ∂y

En tal caso

∂f ∂f
df = dx + dy = M(x, y )dx + N(x, y )dy = 0 ,
∂x ∂y

y (2) es equivalente a
df = 0 .
Luego

y = y (x) , x ∈I es solución de (2) ⇐⇒ f (x, y (x)) = c , ∀x ∈ I .

3/12
Criterio para una ecuación diferencial exacta:
Sean M, N funciones con derivadas parciales continuas de primer orden en B ⊂ R2 . Entonces

∂M ∂N
M(x, y )dx + N(x, y )dy = 0 es exacta en B ⇐⇒ ≡ en B .
∂y ∂x

∂f = M y ∂f = N en B. Entonces
Demostración: =⇒ Supongamos existe f tal que ∂x ∂y

∂M ∂ 2f ∂ 2f ∂ ∂f ∂N
 
= = = = .
∂y ∂y ∂x ∂x ∂y ∂x ∂y ∂x

⇐= Sabemos que ∂M ∂y
= ∂N
∂x
en B. Debemos encontrar f tal que ∂x∂f = M y ∂f = N.
∂y
Integremos, por ejemplo, M(x, y ) con respecto a x, manteniendo y constante:
Z
f (x, y ) = M(x, y )dx + g (y ) , (3)

∂f = N. Entonces
donde g (y ) es la ”constante” de integración a determinar, tal que se cumpla ∂y

∂f ∂
Z
0
(x, y ) = M(x, y )dx + g (y ).
∂y ∂y

Por lo tanto debemos tener



Z
0
g (y ) = N(x, y ) − M(x, y )dx.
∂y

Por último, se integra con respecto a y y se sustituye en (3) para obtener la función buscada. 4/12
Ejemplo 1:
Encontrar la solución de
2 2 2 3
(3x + 6xy )dx + (6x y + 4y )dy = 0 .

Tenemos M(x, y ) = 3x 2 + 6xy 2 y N(x, y ) = 6x 2 y + 4y 3 . Comprobemos si es exacta:

∂M ∂N 2
(x, y ) = 12xy = (x, y ) , ∀(x, y ) ∈ R Es Exacta!!
∂y ∂x

∂f = M y ∂f = N. Comencemos por integrar ∂f = M con respecto a x


Luego existe f tal que ∂x ∂y ∂x
Z
2 2 3 2 2
f (x, y ) = (3x + 6xy )dx + g (y ) = x + 3x y + g (y ) .

∂f = N, para eso derivemos con respecto a y la expresión anterior e igualamos a N:


Usemos la otra condición ∂y

∂f 2 0
(x, y ) = 6x y + g (y ) , y como
∂y
2 3
N(x, y ) = 6x y + 4y , debemos tener
0 3
g (y ) = 4y
4
=⇒ g (y ) = y .

Por lo tanto
3 2 2 4
f (x, y ) = x + 3x y + y .
Ası́, las curvas soluciones (x, y (x)) de la ecuación satisfacen (f (x, y ) = c)

3 2 2 4
x + 3x y + y = c, c ∈ R. 5/12
Ecuaciones que se reducen a EDO exactas:
Factor integrante.

6/12
Factor integrante

Definición: Si la ecuación

M(x, y )dx + N(x, y )dy = 0

no es exacta en D, se llama factor integrante a toda función u = u(x, y ) definida en


D tal que
u(x, y )M(x, y )dx + u(x, y )N(x, y )dy = 0
es una ecuación exacta en D.
Ejemplo: La ecuación

y (1 + xy )dx − xdy = 0 No es exacta!! (4)

En efecto, si M(x, y ) = y (1 + xy ) y N(x, y ) = −x, se tiene

∂M ∂N
= 1 + 2xy 6= = −1.
∂y ∂x

7/12
Factor integrante
Sin embargo para y 6= 0 la ecuación multiplicada por µ(x, y ) = 12
y

y (1 + xy ) x
dx − dy = 0 ,
y2 y2

es decir
1 x
( + x)dx − dy = 0 , Es exacta!! (5)
y y2

Ası́ u(x, y ) = 12 es factor integrante de (4) en el conjunto D = {(x, y )/y 6= 0}.


y
Para resolver (5) ponemos

1 x x2
Z  
f (x, y ) = + x dx + g (y ) = + + g (y ) .
y y 2

Pero como
∂f x 0 x
(x, y ) = − + g (y ) = − ,
∂y y2 y2
tenemos
0
g (y ) = 0 =⇒ g (y ) = c .
Por lo tanto, toda solución y = y (x) de (5) verifica la ecuación implı́cita

x x2
+ = c, para algún c ∈ R .
y 2

Despejando y obtenemos

2x 2
y (x) = , c ∈ R , para x tal que 2c − x 6= 0 .
2c − x 2 8/12
Factor integrante

Observación: Para que u = u(x, y ) sea factor integrante de Mdx + Ndy = 0 , debemos
tener
∂ ∂
(uM) = (uN). (6)
∂y ∂x

Casos en que es fácil encontrar el factor integrante

1) Si  
1 ∂M ∂N
− = f (x) (depende sólo de x) ,
N ∂y ∂x
entonces nuestra ecuación admite factor integrante que depende solo de x.
En efecto si u = u(x) la ecuación (6) queda de la forma

u 0 (x)
 
∂M ∂N 1 ∂M ∂N
u = u 0 (x)N + u =⇒ = − = f (x) ,
∂y ∂x u N ∂y ∂x

lo que implica
Z R
f (x)dx
ln(u(x)) = f (x)dx o bien u(x) = e .

9/12
Factor integrante

2) Si  
1 ∂M ∂N
− − = g (y ) (depende sólo de y) ,
M ∂y ∂x
entonces nuestra ecuación admite factor integrante que depende solo de y .

En efecto si u = u(y ) la ecuación (6) queda de la forma

u 0 (y )
 
∂M ∂N 1 ∂M ∂N
u 0 (y )M + u =u =⇒ = − − = g (y ) ,
∂y ∂x u M ∂y ∂x

lo que implica
Z R
g (y )dy
ln(u(y )) = g (y )dy o bien u(y ) = e .

10/12
Ejemplo 2
Encuentre la solución del problema de valores iniciales:

y 3
dx + (y − ln(x))dy = 0 y (1) = 1.
x

Solución: Sea
y 3
M(x, y ) = , y N(x, y ) = y − ln(x) ,
x
tenemos
∂M 1 ∂N −1
= , y = , No es exacta!!
∂y x ∂x x

Pero
1 ∂M ∂N x 2 2
 
− − = − · = − ,
M ∂y ∂x y x y

implica que el factor integrante es

ln( 12 )
R 1
−2 dy −2 ln(y ) 1
u(y ) = e y = e = e y = ·
y2

Nuestra ecuación multiplicada por el factor integrante es

1 ln(x)
 
dx + y− dy = 0 Es exacta!!
xy y2

En efecto,
∂ 1 ∂ ln(x) 1
 
= y− =− .
∂y xy ∂x y2 xy 2 11/12
Ejemplo 2
Luego, existe existe f tal que
∂f 1 ∂f ln(x)
= y =y− ,
∂x xy ∂y y2
integrando con respecto a x la primera relación se tiene:

1 1
Z
f (x, y ) = dx + g (y ) = ln(x) + g (y ) .
xy y

Usando la segunda relación


∂f ln(x) 0 ln(x)
(x, y ) = − + g (y ) = y − ,
∂y y2 y2
lo que implica
0 y2
g (y ) = y =⇒ g (y ) = .
2
Por lo tanto
1 y2
f (x, y ) = ln(x) + ,
y 2
y la solución general satisface la ecuación implı́cita

1 y2
ln(x) + = c, c ∈ R.
y 2

Para encontrar la solución que verifica la condición inicial y (1) = 1, se sustituye esta en la expresión anterior para
obtener c = 1/2, luego la solución del PVI es

1 y2 1
ln(x) + = .
y 2 2 12/12

También podría gustarte