Ecuaciones en
devivadas parciales
Ecuación general
(X2 n4f(x
-
, y , z) = skx , y ,
z)
Problema de Poissan
-2f(xy ,
z) =
s(x , y ,
z)
e me
potencial frente
Problema de Laplace
-21 =
-I
byz
-> O
cond. de
f(0r) =
frontera
cte
3 diferencias
finitas
Caso unidimensional
Consideranos el oscilader armónico forzado
ji(t) =
-
wk() -
zbxt) +
#
Cambio de variables
to escalas
q = so
y san para
dimensiones
s =
E quitar
to
= 2 = s
= I dt
2 I =
=- w- +
q"(s) = -
witg-2bto q'+ ,
to
q"(s) = -
witg-2bto q'+ ,
to
Se definen to =
+/zb
w2
= -
4b2
E
xt
=
q"(s) = -
wq -
q(s + Fls)
q"(s) +
q((s) +
wq =
F(s)
Caso particular de
f"(x) +
y(x)f((x) +
q(x)fx) = s(x)
Df =
s - b
operador matricial
Para resolver la ec diferencial
nxx N
3
x xn = n
=
0,1 - -
, ,
f(x) fr = f(xn)
Aproximaciones para las derivadas
x)
f fa) -
= fibin) + fkin)x 0(x2)
2Ax
+
fbl+ f(x 1x)
-
= f"(xn) + f kin) se e
+ 0 (1)
f"(x) +
y(x)f((x) +
q(x)fx) = s(x)
Evaluando en xn -9(n) =
pr
q(xn) =
qn
(xn) Sn
futt S
=
2f(n) f( (x)
(() - + -
Ax2
f(xn fx)
9n() anfr
- -
+
+ =
Su
2Ax
frt/sic Ex) + +
ful-iz+qn)+ fr / -Es) = Sa
1 Esa 9 - 2
Ests(l=
-
San su
frt/sic Ex) + +
ful-iz+qn)+ fr / -Es) = Sa
Para n = 1 ...,
N
,
+z(5 si 2) +
+ +- (q1 -
z) +
to(Ents) =
S1
(Estsiz) +z (92 z) + f
(- +se) Se
=
Es
-
+
:
. +sl+ fre/qu-)+ frt-ts)
=
Sa
trt (25+ siz)
+
fr (qu 52) - +
fr -) -
Est 2) =
Se
Solución {fu fo ,
Art condiciones de frontera
En forma del arreglo
DI =
*
DI =
a ⑪(i
dQue hacemos ca to y frutt
+z(5 si 2) +
+ +- (q1 -
z) = s1-
tof Ents)
fr(qu 552)- +
fr =
( Es2)
-
=
Sr- tr (+ siz)
La ecuación matrical se cawerte en
/Ents)
I I
to
DE =
I -
B can
=
Irt (+ siz)
Operador matricial
a
Ecuación diferencial
f"(x) +
y(x)f((x) +
q(x)fx) = s(x)
Se convierte
(
en
proyección de f
=
D
DI =
-
sobre la malla de Exi]
-
=
fuentes
contra anes
·
de fratera
Solución
=
D" -5)
Ejem! (3H12f(x) =
f =
foq
x
=
xoS =
fo(q" +
Gr a) = 0
+
7CO
Biq
+ 3 r
En 1D
f"(x) +
y(x)f((x) +
q(x)fx) = s(x)
Solución I =
D -
5) I =
{ f(x)3
to y fut San las cad defrontera
En ID
↑2f(xy) =
s(x , y)
Solución I =
D-5") I {f(xi
=
, 9; 3i ,j = 0
3
A
Condiciones de frontera
f(0 y)
,
=
c f(x ,
0) =
c3
f(a y)
,
=
cf(x a) ,
=
Dx
Caso bidimensional
Tenemos la matriz
-
fum =
+(xn Bu) ,
i, =
1 ,2 3
,
Para simplificar el cálalo
3
*
.
trabajamos con el vector
-
- - -
Ny 3 f=
-----
=
fz 1
f32
-D
f
I
+
z
-
fez
732
+1 3
La posición de fum en es
fee
-33
g(n m)
,
=
(m 1)Nx-
+
n
El problema so ruelve undimensional
↑2f(xy) =
s(x , y)
Aproximaciones yora
y en x =
xn , y
=
ym
fla- ya) fania) fant e
Del
- +
2
Ax
fla -- fania) f+ y
Del
+
2
Ay
Usando :
g(n m) ,
=
(m 1)Nx
-
+
n
#genima-fginins
fgintin) +
& x2 Ax2
-
orfginzfgunstfglnintel
& z
-
y
syz
Entonces g(n , m) =
(n 1)Net n
-
+- te
a 2
dx
1x
2
(-z .
1) inte
en -
z 7fs
#
-(- .
-
/ /
, / /
/
--
/
Matriz
(NxNg(X (NxNg)
"
Para construir Doc
Di =
0 - Llenomos las submatrices NaXNy
Condiciones de frontera
- =
552(- 2 , 1 ,
0)()f3
+
r
g(n m)
,
=
(n 1)Net n
-
= 0x2
s=
est en
-((0 1
2)(i)
-
, , +
z
fo f fu
: fe ts
Las derivadas en g(n m) (n 1)Net n
y
=
-
gfgs+
fac
Oy2 Ayz
74
=Enet1 +
syz
-
(,0.
142
(i) +
-
114
2fg(a)
Afgiz
Ay
=t -
c
= 2f z + ts
-00
Ayz
, ,
1) (i) ür
t o ene
0
·
I
y
a t
1
0
1
0 0
0
Finalmente D =
H x
+
DY