Solucion Del Examn Parcial

1.
Universidad
Nacional de Callao
ECUACIONE
S
DIFERENCIA
LES
Ciencia y Tecnologa rumbo al tercer milenio
GRUPO:
7
ECUACIONES DIFERENCIALES PARCIALES
Profesor:
CASTRO VIDAL RAUL PEDRO
Integrantes:
..
ABANTO LPEZ MIJAIL ANDRS

ADVINCULA-EVARISTO-IVAN RICHARD
RICALDI CRISTOBAL JUNIOR ARTURO
SUCASAIRE-MERMA-CHRISTIAN LUCAS
WEDEMEYER DE LAS CASAS JOS ANDRES
YABAR-PRINCIPE-CRISTIAN ENRIQUE
Pgina 0
1323220337
1323220587
1323220061
1323220195
1323220212
1323210197
ECUACIONES DIFERENCIALES
PARCIALES
I.
OBJETIVOS
Reconocer y clasificar una ecuacin diferencial parcial.
Resolver una ecuacin diferencial parcial con el mtodo de
separacin de variables.
Resolver las aplicaciones de una ecuacin diferencial parcial con el
mtodo de separacin de variables.
II.
INTRODUCCION
Muchas veces nos presentamos ante problemas matemticos, fsicos,

qumicos, etc. que parecen no tener solucin, analizamos nuestros
datos, nuestras variables y tratamos de resolver el problema de forma
convencional, as que debemos buscar otro tipo de soluciones un tanto
ms complejas. La solucin de ecuaciones diferenciales parciales
muchas veces se puede simplificar a resolver una ecuacin diferencial
parcial de segundo orden. Problemas tpicos son la propagacin del
sonido o del calor, la electrosttica, la electrodinmica, la dinmica de
fluidos, la elasticidad, la mecnica cuntica y muchos otros. Se las
conoce tambin como ecuaciones diferenciales parciales.
III.
DEFINICION DE ECUACION DIFERENCIAL PARCIAL,

ORDEN Y LINEALIDAD
Una ecuacin diferencial parcial (EDP) es aquella cuyas incgnitas son

funciones de diversas variables, con la peculiaridad de que en dicha
ecuacin figuran no solo las propias funciones sino tambin sus
derivadas. Una variable independiente pueden ser tiempo y posicin, o
bien dos coordenadas de posicin o una combinacin de varias variables
con el tiempo.
1. Entonces, qu es una ecuacin diferencial parcial?
UNIVERSIDAD NACIONAL DEL CALLAOPgina 1
Es una expresin matemtica que contiene una o ms variables

dependientes y dos o ms variables dependientes. Una ecuacin
diferencial en derivadas parciales por su semejanza con las EDO, es una
ecuacin donde una cierta funcin incgnita u viene definida por una
relacin entre sus derivadas parciales con respecto a las variables
independientes.
u=u (x , y , z ),
Si
F x, y ,z ,
una ecuacin diferencial en derivadas parciales sera:
u u u 2 u 2 u 2 u 2 u
,
,
,
,
,
,
, =0
x y z x2 y 2 z 2 x y
Ejemplos:
a)
2 u 2 u
+
=0
x2 y2
b)
2 u 2 u
u
= 2 2
2
t
x
t
2. Orden una Ecuacin Diferencial Parcial

Las ecuaciones diferenciales parciales se pueden clasificar de acuerdo a
su orden, linealidad y tipo de condiciones de frontera.
El orden de una ecuacin diferencial parcial est determinado por la
derivada parcial de mayor orden presente en la expresin como se
muestra en los siguientes ejemplos:
Ejemplo:
a)
2 u 2 u
+
=0
x2 y2
Es una EDP de segundo orden
b) La siguiente EDP
u
u u
+ =0
x y
Es una PDE de primer orden.
3. Existencia y unicidad
Aunque el asunto de la existencia y unicidad de las soluciones de
las ecuaciones diferenciales ordinarias tiene una respuesta muy
satisfactoria resumida en el teorema de Picard-Lindelf, el mismo asunto
para las ecuaciones en derivadas parciales est lejos de estar
satisfactoriamente resuelto. Aunque existe un teorema general,
el teorema de Cauchy-Kovalevskaya, que afirma que para una EDP, que
es analtica en la funcin incgnita y sus derivadas, tiene una nica
solucin analtica.
Aunque este resultado que parece establecer la existencia y unicidad de
la soluciones, aparecen ejemplos de EDP de primer orden cuyos
coeficientes tienen derivadas de cualquier orden (aunque sin ser
analticas) pero que no tienen solucin. Incluso si la solucin de una EDP
existe y es nica, sta puede tener propiedades indeseables.
Un ejemplo de comportamiento patolgico es la secuencia de problemas
de Cauchy dependientes del parmetro n para la ecuacin de Laplace:
2
u u
+
=0
x2 y2
Con condiciones iniciales
u ( x , 0 )=0,
u
sin nx
( x ,0 )=
y
n
es un entero. La derivada de u con respecto a y se

aproxima a 0 uniformemente en x a medida que n se incrementa,
Donde
pero la solucin es:
u ( x , y )=
( sinh ny ) ( sinh nx )
2
n
Esta solucin se aproxima a infinito si
nx
no es un entero mltiplo de
para cualquier valor de y. El problema de Cauchy para la ecuacin de

Laplace se denomina mal propuesto o mal definido, puesto que la
solucin no depende continuamente de los datos del problema. Estos
problemas mal definidos no son usualmente satisfactorios para las
aplicaciones fsicas.
IV.
CLASIFICACIN
ORDEN.
DE
LAS
EDP
DE
SEGUNDO
Las EDP de segundo orden se clasifican habitualmente dentro de cuatro

tipos de EDP que son de inters fundamental, a continuacin se dan
ejemplos de estos cuatro tipos:
Ecuacin
Nombre
Tipo
2 u=0
Laplace
Elptica
2 u 2 2
=c u
t2
Onda
Hiperblica
u
2
=k u
t
Difusin
Parablicas
2 u=ku
Helmholtz
Elptica
Con mayor generalidad, si se tiene una ecuacin de segundo orden del

tipo:
A uxx +2 Bu xy +C u yy + D ux + E u y + F=0
se dice que es elptica si la matriz
[ ]
Z= A B
B C
tiene un
[ ]
tiene un
[ ]
tiene un
determinante mayor a 0.
Se dice que es parablica si la matriz
Z= A B
B C
determinante igual a 0.
Se dice que es hiperblica si la matriz
Z= A B
B C
determinante menor a 0.
V.
LA TRANSFORMADA DE FOURIER
1. DEFINICIN Y PROPIEDADES FUNDAMENTALES

Una de las herramientas ms tiles en el estudio de las Ecuaciones en
Derivadas Parciales (EDP) es la transformada de Fourier. Dada una
funcin f L1 (Rn) definimos su transformacin de Fourier del modo
siguiente:
f ( )=1 /
Conviene sealar que f est bien definida para todo Rn puesto que
e ix f(x) L1(Rn). Esto es efectivamente cierto ya que:
|ei x f (x)|=|f ( x)| L1 ( Rn )

De esta propiedad se deduce que:
^f L (Rn )f L1 ( Rn )
Por otra parte, el Teorema de la Convergencia Dominada (TCD) permite
ver que f es tambin una funcin continua. Por tanto, f BC(Rn).
La transformada de Fourier tiene varias propiedades fundamentales que
la hacen no solamente una herramienta muy til en el contexto de las
EDP sino tambin en muchos otros mbitos de las Matemticas y de las
otras Ciencias y Tecnologas. Una de ellas est relacionada con el
comportamiento de la transformada de Fourier en relacin a los
operadores de derivacin. Tenemos:
^
f
( )=
x
eix
n
( 2 )2
f
i
( x ) dx=
eix f ( x ) dx=i f^ ( )
n
x
(2 )2 R
n
Gracias a esta propiedad, las EDP lineales con coeficientes constantes se

reducen a ecuaciones algebraicas lineales en el espacio de Fourier. En
efecto, consideremos por ejemplo la ecuacin de Laplace:
u=f en R
Tomando la transformada de Fourier obtenemos:

2
|| u^ ( ) =f^ ( ) en Rn (1)
En la obtencin de la ecuacin anterior

esencial que
hemos utilizado de manera
n
2
^
u
=
2
xi
i=1
n
n
2
^
u
^
u= 2 =
i=1 xi
i=1
De (1) Podemos obtener de manera inmediata una expresin para la

transformada de Fourier u^ de la solucin u:
2
u^ ( )=f^ ( ) /||
Pero, para obtener la solucin necesitamos invertir de la transformada

de Fourier. La transformada inversa de Fourier se define del modo
siguiente:
f^ ( x )=
eix f ( ) d
(2 )
n
n
2 R
Nuevamente, si f L1 (Rn ),
est bien definida y es una funcin de
BC(Rn ).
La Gaussiana es un ejemplo particularmente importante en el que la
transformada de Fourier es fcil de calcular explcitamente. Tenemos
^
1
et x ( )=
2t
2
n /2
e| |/4 t
2
En efecto,
^
et|x| ( ) =
2
e
n
(2 )2
ix t|x|
dx=
j=1
t x
ei x e | | d x j
1
(2 )2
Vamos a calcular la siguiente integral:

i
2t
x +
2
t
e
dx=
ix t |x|
Rn
Rn
Donde es el contorno Im (z) = /2t del plano complejo. Deformando los

contornos la integral, por el teorema de los residuos, puede reducirse al
clculo de la misma sobre el eje real:
t z
dz= et x dx= /t
De las dos ltimas ecuaciones obtenemos:
ix t x
2
4t
dx=e
En efecto, dada f Rn, si definimos la funcin

ft ( x )=f (tx)
Tenemos
1
f^ ( )= n
t
n /2
f^ ( /t )
Para comprobar hacemos lo siguiente:

f^t ( ) =
eix f t ( x ) dx=
n
(2 ) 2
eix f ( tx ) dx=
n
(2 ) 2
eiy f ( y ) dy
t
n
n
2 n R
(2 ) t
La identidad sobre la transformada de Fourier de la Gaussiana es muy

til. En particular proporciona la solucin fundamental de la ecuacin
del calor.
En efecto, recordemos que G es la solucin del sistema
2
|x| G/4 t
n
G ( x , t ) =( 4 t ) 2 exp
(*)
eltaG =0 en Rn x ( 0, )
G ( x ,0 )= ( x ) en Rn
Tomando la transformada de Fourier en la primera ecuacin

(**)
^ t +||2 ^ =0, R n , t >0
G
Por lo tanto
2
^ ( , t )=e|| t G
^ ( , 0 ) ()
G
Conviene observar que en (**) slo se toma la transformada de Fourier
en la variable espacial, lo que transforma la ecuacin en derivadas
parciales en una ecuacin diferencial ordinaria con parmetro .
Por otra parte, en virtud del dato inicial de (*) tenemos
^ ( , 0 )= ^ ()
G
Por otra parte
1
^
n
2
(2 )
Son varias las maneras de justificar la identidad de la ecuacin anterior.
Una primera es simplemente aplicar la definicin de la transformada de
Fourier, teniendo en cuenta que es una medida. Se obtiene:
1
1
^0 ( )=
eix 0 ( x ) dx=
n
n
(2 )2 R
( 2 )2
n
Otra posibilidad es observar que la familia

n/ 2
t
( ) e|x| t
Es una aproximacin de la identidad cuando t de modo que

.
n /2
( t )
R
e|x| t ( x ) dx (0)
Para todo BC (Rn). Entonces

t
^0 ( )=lim
t
( )
n
2
lim 1
^
1
t|x|
e
( )= t n e|| /4 t =
n
(2 )2
(2 ) 2
2
Combinando (***) y (****) obtenemos

^ ( , t ) =
1
(2 )
n
2
e|| t
Basta por ultimo aplicar la transformada inversa de Fourier para

obtener G(x, t) a partir de su transformada de Fourier de la ecuacin
anterior. Obtenemos entonces
G ( x , t ) =4 t
n
2
( )
|x|
exp
4t
Esta ltima formula se obtiene de manera fcil gracias a esta otra

^
1 n /2 | |/ 4 t
t x
(
)
e
=
e
2t
2
ecuacin:
En virtud de ^y^ se obtiene:

f^ ( )=f ()
Sigamos ahora estudiando las propiedades fundamentales de la

transformada de Fourier. La identidad de Plancherel garantiza que
f L ( R )=f^L ( R )
2
Para todo f L 1 (R n) L 2 (R n).

Antes de probar (la ltima ecuacin) observamos que esta identidad
garantiza que la transformada de Fourier define una isometra de L 2 (R
n ) en L 2 (R n ). Adems, la transformada de Fourier puede extenderse
por densidad a una isometra definida en todo el espacio L2 (Rn).
Para probar se necesita:
v ( x ) w^ ( x )dx= v^ ( x ) w (x) dx
R
Para todo v, w L1 (Rn). Conviene en primer lugar observar que ambas

integrales estn bien definidas puesto que como v, w L1 (Rn ),
entonces v, w L (Rn ), por lo que vw L1 (Rn).
La identidad anterior es fcil de comprobar. En efecto
v ( x ) w^ ( x )dx= v ( x )
R
eixy w ( y ) dy =
n
(2 ) 2
(2 )
eixy v ( x ) w ( y ) dxdy
n
n
2 R
Que coincide con la expresin de
^v ( x ) w( x )dx
R
Aplicando esta identidad con v(x) = e |x| 2 y usando la expresin de

v antes calculada
^
w()e|| d=
2
y
1
(2 )
w (x) e|x| / 4 dx
n
n
2 R
Obtenemos que
Dada u L 1 (R n ) L 2 (R n ), definimos v(x) := u(x) y w = u v

(que pertenece a L 1 (R n ) C(R n )). Tenemos
n /2
n
^ =u^
w
v=(2 ) u^ v^ L ( R ) (V )
Esta es otra de las propiedades importantes de la transformada de

Fourier: la transformada de Fourier de la convolucion es, mdulo una
constante multiplicativa, igual al producto de las transformadas de
Fourier.
En este caso, como v(x) = u(x), tenemos
v^ ( )=
1
(2 )
eix u (x ) dx=u ()
n
n
2 R
y por tanto
u ()2 .
w
^ =
(2 ) n/ 2
Como w es continua tenemos tambin

2
n/2
lim
Adems
n/2
|u^| d=
R
Lo cual concluye la prueba de la identidad de Plancherel.

Comprobamos ahora la formula (V) que pone en relacin la convolucion
y la transformada de Fourier. Tenemos
ix
( uv)( x ) dx
w
^ ( )= ^
uv ( )=
Despus de desarrollar, integrar vamos a obtener la siguiente ecuacin:
^ ( )=
w
Probemos algunas propiedades ms de la transformada de Fourier. La
identidad de Parseval garantiza que
( u^ v^ d)
( u v dx)=
R
Rn
Para comprobarlo, dados u, v y C, aplicamos la identidad de

Plancherel a la funcin u + v de modo que
2
u+v L ( R )=u^ +^v L ( R )

2
Expandiendo esta identidad tenemos
[|u|2 +||2 +u ( v ) +u( v )] dx= [|u^|2+|^ v|2+u^ ( ^ v ) +u^ ( v^ )] d

R
Tomando sucesivamente = 1 y = i en esta identidad obtenemos

.Comprobamos por ultimo la frmula de inversin de la transformada de
Fourier, es decir,
( f^ )=f .
Dado z R n y > 0 consideramos la funcin

v^ ( ) =
(2 )
n
n
2 R
ixz x2
v ( x)=e
ix ( y z) |x|
1
dx=
(2)
n
2
e|yz| /4
En este punto hemos usado la frmula:

Que identifica la transformada de Fourier de la Gaussiana
Tenemos
ix tx
dx=e
/ 4 t
1 /2
()
Deducimos que
^f ( ) e
iz||
d=
(2 )
f (x)e
n
n
2 R
|xz|
4
dx .
Ahora bien, la expresin de la derecha de esta identidad converge a f(z)

para casi todo z R n , mientras que la de la izquierda converge a ( f)
(z). Esto concluye la prueba de la frmula de inversin. Otra de las
propiedades ms clsicas de la Transformada de Fourier es la que se
obtiene en relacin al operador de traslacin. En efecto, dada una
funcin f L 1 (R) definimos una traslacin de la misma de amplitud a
R n:
f a ( x )=f (x +a)
La transformada de Fourier de la trasladada fa viene dada por:
ia.
f^ a ( )=e f ( )
En efecto, mediante un simple cambio de variables se obtiene
n/2
1
f^ a ( )=
f ( y ) dy=e ia. f^ ( )
^f a ( )= 1
De manera anloga, se comprueba que

^
e ia. f ( )=f a ( )
Para la Transformada
semejantes:
Inversa
de
Fourier
fa (x)=eia . x f (x)
Y
e
ia.
se
obtienen
frmulas
VI.
LA ECUACIN DE TRANSPORTE LINEAL
Las ecuaciones que modelan fenmenos de propagacin de ondas y

vibraciones tpicamente Ecuaciones en Derivadas parciales (EDP) de
orden dos. Sin embrago en todas ellas subyacen las ecuaciones de
Transporte de orden en uno que analizamos en esta seccin.
El modelo ms sencillo es:
(1)
ut
ux
=0
Como se vio anteriormente es fcil comprobar que
u=u(x , t)
es solucin
de esta ecuacin si y solo si es constante a lo largo de las lneas

caractersticas.
(2)
x + t = cte
De este modo deducimos que las ecuaciones de (1) son de la forma:

(3)
u=f (x +t)
Donde f es el perfil de la solucin en el instante t=0 :

u ( x , 0 )=f (x )
(4)
Como se dijo anteriormente obtuvimos esta frmula de representacin

de la solucin mediante la transformada de Fourier.
La solucin (3) de (1) es entonces una simple onda pura en el que el
perfil f se transporta (avanza) en el eje real a velocidad constante uno.
Al invertir el sentido del tiempo (haciendo el cambio de variable t a t) la
ecuacin se transforma en:
(5)
ut +u x =0
Cuyas soluciones son de la forma:
u=g ( x , t)
(6)
Tratndose de ondas viajeras que se propagan en direccin opuesta con

velocidad uno. Vemos por tanto que las soluciones de la ecuacin de
transporte pueden calcularse de manera explcita y que en ellas se
observa un sencillo fenmeno de transporte lineal sin deformacin.
Como decamos esta ecuacin de transporte de orden uno subyace en
muchas de las ecuaciones d onda que intervienen en la teora de las
vibraciones.
As
por
ejemplo
el
operador
de
ondas
2t 2x = ( t x )( t + x )=(t + x )(t x )
En el caso en el que el transporte se produzca en un medio no
homogneo, la ecuacin correspondiente adopta la forma:
(7)
ut +C ( x ) u x =0
En este caso, nuevamente la solucin es constante a lo largo de las

caractersticas que son soluciones de:
(8)
x ( t )=c ( x ( t ) ) , t> 0 ; x ( 0 )=x 0
En efecto es fcil comprobar que si u es solucin de (7) se tiene que

u(x(t), t) es independiente de t para cualquier curva caracterstica.
En realidad, tal y como ocurran en los coeficientes constantes, esta
propiedad permite definir la solucin. En efecto suponiendo que C = C(x)
x0
es una funcin Lipshchitz, para cada
existe una nica solucin de
(8). Adems las caractersticas estn globalmente definidas y no pueden
cortarse por unicidad de la solucin de la EDO (8).Por ltimo por cada
1
punto del espacio tiempo ( x 1 , t
pasa una nica caracterstica. Para
comprobarlo basta resolver la ecuacin:

x ( t )=c (x ( t ))
Con dato inicial
x ( t 1 )=x 1
1
(8) pasa por ( x 1 , t
y observar que la caracterstica solucin de
si el dato inicial
x0
coincide con el valor de esta
solucin en t=0, x(0).

De este modo vemos que la solucin de la ecuacin de transporte (7)
esta globalmente definida por su dato inicial
u ( x , 0 )=f ( x ) .
VII. LA ECUACIN DEL CALOR
Como hemos mencionado en la introduccin la ecuacin del calor es el

prototipo de ecuacin de evolucin de tipo parablico cuyas variantes
estn presentes de manera sistemtica en todos los modelos
matemticos de la difusin y de la Mecnica de Fluidos. Se trata de un
modelo fuertemente irreversible en tiempo en el que la informacin se
propaga a velocidad infinita. Estas propiedades quedaran claramente de
manifiesto a lo largo de esta seccin en la que recordamos sus
principales propiedades analticas.
El problema de valores iniciales en R

Considere el problema de valores iniciales para la ecuacin del calor.
ut u x =0 xR ,t >0
u ( x , 0 )=f ( x ) xR
(13)
Se trata evidentemente de una ecuacin con orden dos lineal con

coeficientes constantes y por tanto susceptibles de que se aplique el
teorema de C-K.
A pesar de que el problemas de los valores iniciales (1) no entra en el
marco del TEOREMA DE LA C-K cabe abrdalo a partir del mtodo de
Cauchy basado en el desarrollo de las soluciones en serie de potencias.
Supongamos que el dato inicial f = f(x) es una funcin analtica de modo
que:
f ( x )= f k x k
k ..0
Y busquemos una solucin analtica de la forma:
u ( x , t )= uk ( t )x
k0
Para que la ecuacin del calor se satisfaga es necesario.
u ( x , t )= u k (t) x k
k0
k (k 1)uk ( t ) x k2=0
k0
Esta identidad se verifica si y solo si:
u k ( t )( k+ 2 )( k +1 ) u k+2 (t )=0,k 0
Por otra parte el lado inicial u(x,0)=f(x) de la EDP permite identificar el

uk
dato inicial de cada uno de los coeficientes
.
uk ( 0 )=f k
, k 0
Se trata de un sistema de ecuaciones diferenciales con un nmero

infinito de incgnitas. Su resolucin excede obviamente de las tcnicas
que la Teora EDO proporciona.
Queda pendiente entender si los coeficientes que hemos identificado de
este modo garantizan la convergencia de las series de potencias
correspondientes.
La solucin puede ser construida de manera explcita mediante la
solucin fundamental o ncleo de GAUS.
G ( x , t ) =4 t
1 /2
2
exp ( x /4 t
Es fcil comprobar que el ncleo G verifica las siguientes propiedades.
AUTOSEMEJANZA: G es de la forma
1
x
G( x , t)= F
t t
( )
Dnde:
F( z ) = 4 t 1/ 2 exp ( z 2 / 4 t
CONSERVACIN DE LA MASA.
G ( x , t ) dx =1,t> 0
C
La funcin G es de clase
>0.
La amplitud mxima de la funcin G decrece con el factor
multiplicativo del orden
(de hecho analtica real para todo t
t 1 /2 .
G es una funcin par con respecto a la variable espacial x.
Es adems fcil comprobar que el ncleo de Gauss es una solucin de

la ecuacin del calor. En efecto en vista de la estructura auto
semejante de G tenemos.
t=
1
2t
3
2
( xt ) 2tx F ( xt )
2
G
Y
x / t
1
Gzz = 3/ 2 F )
t
Nos queda identificar el dato inicial que el ncleo de Gauss toma en t=0.
En vista de su estructura auto semejante es fcil comprobar que cuando
+
t-- 0 , G(t) converge a la masa o delta de DIRAC en x=0 en el sentido
que:
x
t
ty
+
F () ( x ) dx= F ( y )
G ( x , t ) ( x ) dx= 1t
Para toda funcin continua y acotada
Tomando a
x , y , t
como infinitsimos, la ecuacin
escribirse as:
{ }
T
2 t 2 k
=k 2
k = a
t
c
x2
puede
Esta es la Ecuacin de conduccin del Calor , el factor

coeficiente de trasmisin del calor .
es el
En tres dimensiones la ecuacin (19) puede generalizarse as:
2
=k
+ 2+ 2
2
t
x y z
Una vez obtenida la ecuacin del calor se procede a su solucin

analtica; recordemos que lo primero que hacemos es comprobar, que
realmente, es una ecuacin parablica, tomando la forma general de
una ecuacin diferencial parcial de segundo orden lineal.
2 u
2 u
2 u
u
u
A 2 +B
+C
+D
+E
+ Fu=G
2
X Y
x
y
x
y
Verificamos que s de qu
B 4 AC =0 , esta se puede llevar bajo la
transformacin a fin. A su forma cannica
U t =U xx +F ( U x , U , X ,t )
Donde F es una relacin de estas cuatro variables.

Nota: Recordemos que se dice parablica ya que su forma cannica
U t =U xx
tiene forma de parbola, es decir,
, va la transformada de
Fourier esto es S=k2, se escapa del alcance de estas notas, pero la
transformada de Fourier cambia (en algn sentido cambia las derivadas
por polinomios) ver que

solucin de
F ( n u ) =( ik ) F( u)
donde i es el imaginario
i =1
As basta hacer un estudio de este caso para las ecuaciones parablicas.

EL PROBLEMA DE DIRICHLET
El problema de Dirichlet es un problema que consiste en hallar una
funcin que es la solucin de una ecuacin en derivadas parciales (EDP)
en el interior de un dominio de
Rn
(o ms generalmente una variedad
diferenciable) que tome valores prescritos sobre el contorno de dicho

dominio.
Para un dominio D teniendo la suavidad suficiente sobre el contorno
D , la solucin general al problema de Dirichlet es:
u ( x ) = vs
Dnde:
G( x , s) ds
n
G( x , s)
es la ecuacin de Green para la ecuacin den
derivadas parciales, y :
G( X , S)
G(x , s )
=n . s G ( x , s )=
n
si
i
, es la derivada de la funcin de
Green a lo largo del vector unitario normal apuntando a su interior

La integracin se realiza sobre el contorno con la medida
n .
s . La
funcin v(s) est dada por la solucin nica de la ecuacin de Fredholm

de segunda clase.
F( x ) =-
v (x)
G(x , s)ds
+ v (s)
2
n
, la funcin de Green a utilizar en la
integral de arriba desaparece en el contorno.

G(x, s)=0
Para s D Y x D . Tal funcin de Green usualmente es una suma de

las funciones de Green del campo libre una solucin armnica a la
ecuacin diferencial.
El problema de Dirichlet para funciones armnicas siempre tiene
solucin, y esa solucin es nica cuando el contorno es suficientemente
suave y f(s) es continua. Ms precisamente, tiene solucin cuando:
D C(1,) , Para 0< donde
C(1,)
es la condicin de Holder.
Las condiciones anteriores pueden relajarse pudindose probar que el

problema de Dirichlet admite solucin continua para un dominio convexo
(sin requerir condiciones de suavidad) o cuando cada punto del contorno
pertenece a una bola cerrada ntegramente contenida en el conjunto
complementario del interior del dominio (condicin de la bola cerrada
exterior).
Supngase la existencia de dos funciones armnicas en una regin R
simplemente conexa que cumplen:
u1= u2
x R
x R
u1 R =u2 R =F ( x)
Se construye la siguiente funcin.

u0=u 1u2
x R
Por construccin la funcin dad cumple.

u
( 1u 2)= u1u 2=0 y u 0 R =u1 R u2 R =0
u0=
Integrando la norma del gradiente de
u0
en todo el volumen
delimitado por el contorno y aplicando el teorema de Green se tiene

que:
u0
( u0 )dv u0 u0 dv= u0 u0 ds u0 u0 dv
u 0 u0 dv=
Pero:
u 0=0 entonces : u0 u0 dv=0
u0 R =0 entonces : u0 u0 ds=0
De modo que:
2
u0 u 0 dv =||u0|| dv=0 entonc es u0=0 x R entonces

u0=cte x R
Dado que
u0=0
u0 R =0 se tiene finalmente:
x R
Entonces:
u1=u2
x R
Los problemas de Dirichlet son tpicos de las ecuaciones en derivadas

parciales elpticas, la teora del potencial, y la ecuacin de Laplace en
particular. Otros ejemplos son la ecuacin bi-armnica y las ecuaciones
relacionadas con la teora de la elasticidad.
Este es uno de los problemas de varios tipos de clases problemas de
EDP definidos por la informacin dada en el contorno, entre los cuales
estn tambin el problema de Neumann y el problema de Cauchy.
VIII. LA ECUACIN DE BURGES

Este ejercicio tiene una finalidad puramente didctica. Veremos algunas
de las particularidades de los mtodos pseudoespectrales para la
resolucin de ecuaciones en derivadas parciales en dominios peridicos.
El uso de un mtodo espectral, de formulacin mucho ms complicada,
se justifica por ser mucho ms rpido. Aunque parezca sorprendente la
transformada rpida de Fourier o FFT (por sus siglas en ingls Fast
Fourier Transform) es una operacin especialmente eficiente en trminos

computacionales. En este primer ejemplo, en el que se resuelve una
ecuacin con trminos de conveccin y difusin con un mtodo explcito,
no se hacen evidentes las ventajas. Se hacen obvias cuando, resolviendo
las ecuaciones de Navier Stokes en dominios bidimensionales o
tridimensionales, se usan mtodos de integracin temporal un poco ms
sofisticados. He aqu la ecuacin de Burgers:
t u+u . u=v 2 u
En una dimensin:
ut =v u xxu u x
La ecuacin se esquematiza de la siguiente manera:
t u=L ( u ) D (u )
Como el dominio es peridico podemos suponer que la forma de la
solucin es u=u^ exp ( ikx ) , que es equivalente a plantear la ecuacin para
el espectro de velocidad du transformada de Fourier.
u^
t u^ =L(^u )D )
En este caso la forma del operador lineal es:

2
L ( u^ )=k v u^
Es obvio que sin el trmino no lineal la resolucin de la EDO sera trivial.
Bastara con transformar la condicin inicial al espacio de frecuencias e
integrar sobre l.
No existe ninguna manera exacta de convertir la operacin
u ut
en un
operador D( u^ ), este clculo debe realizarse en el espacio real. Para

cada paso temporal ser necesario anti-transformar la solucin, calcular
el operador D(u), transformarlo y sumarlo en el espacio de frecuencias.
La ecuacin de Burgers nos reserva un fenmeno no esperado: el

aliasing. Cuando un operador diferencial es lineal no se mezcla el
espectro de ninguna manera. Por ejemplo, el operador derivada segunda
xx u^ =k 2x u^ no introduce ninguna operacin entre cada uno de las
componentes de
u^ ; lo mismo sucede con los filtros proporcionales
como el pasa-altos o el pasa-bajos. El operador

introducir aliasing en
IX.
D (u^ ) no es lineal e
u^.
TRANSFORMADAS INTEGRALES DE FUNCIONES
TRANSFORMADA DE LAPLACE
DEFINICION 1: Dada una funcin f(t) definida (al menos) para cada t >
0, (si se prefiere se puede suponer que f(t) = 0 para cada t < 0) la
transformada de Laplace de f es otra funcin (que designamos por
L(f) o F) y que viene dada por
+
L ( f ) ( s )=F ( s ) = f ( t ) est dt ,
0
DEFINICION 2 Se dice que una funcin f(t) es de orden exponencial

cuando t +, si existen M > 0 y T > 0 tales que se verifica
|f(t)| Met para cada t > T.
EJEMPLO 2.1:
1) Toda funcin acotada es de orden exponencial cuando t +,
para todo 0, dado que
|f(t)| M Met para cada t 0.
En particular, las funciones constantes, seno y coseno son de orden
exponencial cuando t +, para todo 0.
2) Si f(t) et 0 cuando t +, entonces f(t) de orden exponencial
cuando t +. En particular, cualquier polinomio es de orden
exponencial cuando t +, para todo > 0.
3) Sin embargo, la funcin f(t) = et2 no es de orden exponencial

cuando t +, para ningn 0.
Supuesto que lo fuera, existiran IR, M > 0 y T > 0 tales que et2
Met para cada t > T. Entonces, et(t) M para cada t > T y tomando
limites cuando t +, llegamos a la contradiccin de que + M.
EJEMPLO 2.2:
1) Probar que no existe la transformada de Laplace L (f)(s) de la
funcin f(t) = et2 para ningn valor de s. Indicacin: Probar que para
cada s IR se verifica
et2st es2/4 t > 0.
2) Probar que no existe ninguna funcin f(t) o funcin continua a trozos
en [0,T] para todo T > 0 y de orden exponencial cuando t + tal
que
L ( f ) ( s )=
s
3 s1 .
Las principales propiedades de la transformada de Laplace que vamos a

utilizar son las siguientes:
i) La transformada de Laplace es lineal: dadas dos funciones f1 y f2,
L(f1 + f2) = L(f1) + L(f2), , IR.
ii) Supongamos que f es continua en [0,+), de orden exponencial
cuando t + y que f(t) es continua a trozos en [0,T] para cada T > 0.
Entonces, si s > , L(f)(s) = sL(f)(s) f(0).
iii) Supongamos que f y f son continuas en [0,+), de orden exponencial
cuando t + y que f(t) es continua a trozos en [0,T] para cada T > 0.
Entonces, si s > ,
L(f)(s) = s2L(f)(s) sf(0) f(0).
TEOREMA Supongamos que f(t) y g(t) son funciones continuas a trozos
en [0,T] para todo T > 0 y de orden exponencial cuando t +.
Entonces
L(f * g)(s) = L(f)(s) L(g)(s),
si s > .
APLICACIN DE LA TRANSFORMADA DE LAPLACE A LAS EDP

Consideremos el problema de la vibraciones de una cuerda infinita
colocada en el eje x > 0, supuesta conocida la posicin del extremo x =
0 en cada momento, junto a las condiciones iniciales en cada punto
utt ( x , t )=u xx ( x ,t ) , x >0, t> 0 Ec . de Ondas
u ( x , 0 )=f ( t ) t >0 Condicinfrontera
u ( x , 0 )=0, x >0 Posicininicial ( en reposo )
u ( x , 0 )=0, x >0 Velocidad inicial (nula)
Desde un punto de vista fsico, nos interesan las soluciones u(x,t) que
estn acotadas cuando x +. En lo que sigue, volvemos a suponer
que todas las funciones que aparecen, verifican las hiptesis necesarias
para que tengan sentido los clculos que vamos a llevar a cabo.
Aplicando al problema la transformada de Laplace respecto de la
variable t, llamando
+
U ( x , s )=L ( u ( x , t ) ) ( s )= u ( x , t ) est dt
0
y utilizando la propiedad iii) sobre la transformada de la derivada

segunda respecto de t, se obtiene
s 2 U ( x , s )su ( x ,0 )u1 ( x , 0 ) =U xx ( x , s ) , x> 0, s> y

U ( 0, s ) =L ( f ) ( s ) , s> y
(3.34)
Teniendo en cuenta las condiciones iniciales del problema (3.33) nos
queda
s 2 U ( x , s )=U xx ( x , s ) , x >0, s> y

U ( 0, s )=L ( f ) ( s ) , s> y
(3.35)
Fijado s > 0, la EDP anterior podemos resolverla como una EDO lineal
de segundo orden respecto de x, por lo que integrando resulta
U(x, s) = c1(s)esx + c2(s)esx,
donde c1 y c2 son funciones arbitrarias. Si u(x,t) est acotada cuando x
+, tambin debe estar
U(x, s) acotada cuando x +, por lo que obligatoriamente c1(s) 0.
Por ello,
U(x, s) = c2(s)esx,
y utilizando la condicin frontera, concluimos que c2(s) = U(0,s) = L(f)(s),
de donde resulta
U(x, s) = L(f)(s) esx.
Queremos ahora aplicar transformada inversa y obtener la solucin del
problema inicial (3.33). Para ello, volvemos a usar la propiedad indicada
en el apartado anterior, pero ahora con a = x y f(t) = f(t) , resultando
u( x ,t )=
0 , sit <0
f ( ts ) , si t>0
Esta solucin se puede interpretar del siguiente modo: la posicin del

extremo x = 0 en cada instante de tiempo t se transmite al punto x0 > 0
de la cuerda en el instante t + x0, lo cual significa en particular que la
onda se mueve hacia la derecha con velocidad 1.
PROBLEMAS
ENUNCIADO
Ecuacin de Laplace (en polares)
Distribucin de temperaturas en la corona circular:
2 u=0
u ( 1, )=f ( ) =sen ( )
u ( 2, )=g ( )= ( 2 )
SOL.
Observando el enunciado, podemos ver que en el contorno hay dos
condiciones no nulas. Esto significa que este problema no puede ser
resuelto en forma directa aplicando el mtodo de separacin de
variables. Por lo tanto, lo resolveremos aplicando el principio de

superposicin. En este caso, el mismo consiste en descomponer el
problema en dos ms simples, resolver cada uno de ellos por separado y
luego obtener la solucin del problema original haciendo la suma de las
soluciones obtenidas para cada uno de esos dos problemas.
Entonces proponemos los siguientes dos problemas:
Problema 1
2 u 1=0 u 1 ( 1, )=f ( )=sen ( )

u1 ( 2, )=0
Problema 2
2
u 2=0
u2 ( 1, )=0
u 2 ( 2, )=g ( )= (2 )
Una vez resueltos estos dos problemas y obtenidas
u1 ( r , )
u2 ( r , ) ,
la solucin del problema del enunciado ser la siguiente:

u ( r , )=u1 ( r , ) +u2 ( r , )
SOL. (Problema 1)
Nota: vamos a llamar
u1=u
para no arrastrar el subindice 1 durante la
resolucin del problema.

Para resolver este problema vamos a utilizar el mtodo de separacin de
variables, que parte de suponer que podemos expresar la funcin
u ( r , ) como el producto de una funcin que depende slo de la
variable
y otra funcin que depende slo de la variable
u ( r , )=R ( r ) ( ) (1)
El Laplaciano expresado en coordenadas polares es el siguiente:

1
1
2 u ( r , ) =urr + ur + 2 u (2)
r
r
El problema busca hallar una funcin que cumpla con
2 u=0 ,
entonces:
1
1
urr + ur + 2 u =0(3)
r
r
A continuacin evaluaremos en la ecuacin (3) las distintas derivadas
parciales de u ( r , ) , a partir de la expresin (1):
ur =R' ( r ) ( )
urr =R ' ' ( r ) ( )
u =R ( r ) ' ' ( )
Reemplazamos en (3):
1
1
R' ' (r) ( ) + R' ( r ) ( ) + 2 R ( r ) ' ' ( )=0
r
r
Dividimos ambos miembros por
R (r) ()
donde
R (r )
ser distintas de la funcin nula ya que de lo contrario

funcin nula.
()
u (r , )
deben
sera la
R ' ' (r ) ( ) 1 R ' ( r ) ( ) 1 R ( r ) ' ' ( )

+
+
=0
R ( r ) ( ) r R ( r ) () r 2 R ( r ) ()
''
'
''
R (r ) 1 R ( r ) 1 ( )
+
+ 2
=0
R ( r ) r R ( r ) r ( )
Separamos de cada lado de la igualdad las funciones de distinta

variable:
''
'
''
R (r ) 1 R ( r ) 1 ( )
+
=
R ( r ) r R ( r ) r 2 ( )
R' ' (r) R ' ( r ) '' ( )

r
+r
=
(4)
()
R (r )
R (r )
2
En la expresin (4) el miembro izquierdo de la igualdad depende slo de

la variable r y el miembro derecho depende slo de la variable .
Para que se cumpla la igualdad para todo
r ( 1,2 )
y para todo
( 0,2 ) , entonces ambos miembros deben ser iguales a una misma

constante :
r2
R' ' ( r )
R' (r ) ' ' ( )
+r
=
= (5)
( )
R (r )
R (r )
De aqu se desprenden las siguientes ecuaciones diferenciales:
}{
' ' ( )
=
()
' ' ( )+ ()=0
(6)
''
'
2 ''
'
R (r )
R (r )
r R ( r )+ r R ( r ) R ( r )=0 (7)
r2
+r
=
R (r )
R (r )
Por las caractersticas del recinto del problema, un punto con

y con un ngulo
=0
y otro con el mismo
r 0 ( 1,2 )
y con
r 0 ( 1,2 )
=2
representan al mismo punto., Por lo tanto la funcin y sus derivadas

deben tomar el mismo valor. Bajo esta condicin, utilizando la expresin
(1), planteamos:
u ( r 0 , 0 )=u ( r 0 , 2 ) R ( r 0 ) (0)=R ( r 0 ) (2 )
u ( r 0 , 0 )=u ( r 0 ,2 ) R ( r 0 ) ' ( 0 )=R ( r 0 ) ' (2 )
Estas condiciones se deben cumplir para todo

R (r)
r ( 1,2 ) . Sabemos que
no puede ser la funcin nula, ya que en ese caso
funcin trivial y sta no cumple con la condicin
u (r , )
sera la
u ( 1, )=f ( )=sen ( ) .
Entonces deducimos que:

( 0 ) = ( 2 ) ( I )
'
'
(0)= (2 )(II )
Adems, las condiciones de contorno
u ( 1, )=f ()
u ( 2, )=0
con la
expresin (1) resultan:

u ( 2, )=R ( 2 )
( )=0 R ( 2 )=0( III )
u ( 1, )=f ( )=sen ( ) (IV )
La ecuacin (6) y las condiciones
(I)
( II)
forman un problema de
Sturm Liouville Simplificado (PSLS). Por lo tanto los posibles valores de

son 0
Resolvemos el caso
=0
para
( )
La solucin de la ecuacin diferencial (6) con
=0
es:
'
( )=A 0 + B0 ( )= A0
Aplicando la condicin de contorno
(I)
tenemos:
( 0 ) = ( 2 ) B0 =2 A0 + B0 A 0=0
Aplicando la condicin de contorno
( II)
obtenemos:
'
'
( 0 )= ( 2 ) A0 A 0
No tenemos restricciones sobre
B0
y por lo tanto:
0 ( ) =B 0 (8)
=0
Concluimos entonces que

autofuncin
es autovalor y est asociado a la
0 ( ) =1 .
Resolvemos el caso
>0
para
( )
La solucin de la ecuacin diferencial (6) con

( )=A cos ( ) + B sen ( ) (9)
>0
es:
' ( )=A sen ( ) + B cos ( )
Utilizando la condicin de contorno
(I ) :
( 0 ) = ( 2 ) A= A cos ( 2 ) +B sen ( 2 )
A [ cos ( 2 ) 1 ] + B sen ( 2 ) =0(10)
Utilizando la condicin de contorno
(II) :
'
'
( 0 )= ( 2 ) B =A sen ( 2 )+ B cos ( 2 )
>0 B=A sen ( 2 ) + B cos ( 2 )
A sen ( 2 ) + B [ cos ( 2 )1 ]=0(11)
Las ecuaciones (10) y (11) forman un sistema de dos ecuaciones

homogneo con incgnitas A y B . Expresado en forma matricial
queda:
[ cos ( 2 )1 ]
sen ( 2 )
)(
sen ( 2 )
A =0
[ cos ( 2 ) 1 ] B 0
)()
Dado que se trata de un sistema homogneo, si hay solucin nica, sta

es la solucin trivial. Para que haya una solucin diferente de la trivial el
determinante de la matriz debe ser cero. Entonces buscamos que:
det
[ cos ( 2 ) 1 ]
sen ( 2 )
sen ( 2 )
=0
[ cos ( 2 )1 ]
[ cos ( 2 )1 ] [ cos ( 2 )1 ]+ sen ( 2 ) sen ( 2 )=0
cos2 ( 2 )2cos ( 2 ) +1+sen 2 ( 2 )=0

2
cos ( 2 ) + sen ( 2 ) 2cos ( 2 ) +1=0

1
12 cos ( 2 ) +1=0 22 cos ( 2 )=0

cos ( 2 ) =1 2 =2 n
=n
=n2 (12)
Utilizando la expresin (12) en la ecuacin (9) obtenemos:

( )=A cos ( n ) + B sen ( n )
Para cada valor de
n=1,2,3 .
que da lugar a la solucin
, existe un par de valores
n ( )
correspondiente. Por lo tanto tenemos:
n ( ) =A n cos ( n ) + Bn sen ( n ) n N (13)
Ahora buscaremos la solucin
Resolvemos el caso
=0
para
R (r)
de la ecuacin diferencial (7).
R (r)
En este caso la ecuacin diferencial (7) se reduce a:

r 2 R ' ' ( r )+ r R ' ( r )=0
La solucin de esta ecuacin es la siguiente:

R0 (r )=C 0 ln (r )+ D0
Aplicamos la condicin de contorno
( III ) :
R0 ( 2 )=C 0 ln (2)+ D0 =0 D 0=C 0 ln (2) R0 ( r )=C 0 ln ( r )C 0 ln (2)
[ ( )]
R0 ( r )=C 0 ln
r
(14 )
2
Resolvemos el caso
>0
R (r)
para
La solucin general de la ecuacin diferencial (7) para

R ( r )=C r + D r
Reemplazamos
a partir de (12):
Rn ( r )=C r n+ D r n
Aplicamos la condicin de contorno

n
(III ) :
2n
Rn ( 2 )=C 2 + D 2 =0 D 2 =C 2 D=C 2
Por lo tanto:
Rn ( r )=C r nC 22 n rn
Sacando factor comn
C 2n :
>0
es:
Rn ( r )=2n C ( 2n r n2n r n ) =E
Para cada
Rn ( r )=En
r n r
2
2
(( ) ( ) )
existe una solucin
r n r
2
2
(( ) ( ) )
(15)
u (r , )
La solucin
R ( r ) , entonces tenemos:
ser una combinacin lineal de las
un ( r , )
obtenidas.
u ( r , )=R0 ( r ) 0 ( ) + Rn ( r ) n ( )
n=1
Utilizando las soluciones (8), (13), (14) y (15) en la expresin anterior:
r
r n r
u ( r , )=B0 C0 ln
+ [ A cos ( n ) + Bn sen ( n ) ] En
2 n=1 n
2
2
()
(( ) ( ) )
Agrupando los coeficientes queda:
u ( r , )=B
0 C 0 ln
F0
u ( r , )=F 0 ln
r
+ An En cos ( n )+ B
n En sen ( n )
2 n=1
() [
Fn
Gn
r n r
2
2
] (( ) ( ) )
Para hallar los coeficientes
(( ) ( ) )
r
r
r
+ [ F n cos ( n )+ G n sen ( n ) ]
2 n=1
2
2
()
F0
Fn
(I V ) :
Gn
(16)
utilizamos la condicin
u ( 1, )=f ( )=sen( )
1
1 n 1 n
1 n 1 n
u ( 1, )=F 0 ln
+ F n
cos ( n ) +Gn
sen ( n ) =sen( )
2 n=1
2
2
2
2
()
H0
(( ) ( ) )
(( ) ( ) )
Hn
In
H 0 + H n cos ( n ) + I n sen ( n ) =sen ()(17)

n=1
H0 ,
Hn , e
In
pueden interpretarse como los coeficientes de un
desarrollo en Serie Trigonomtrica de Fourier (DSTF) de la funcin

sen .
Recordemos que
{ 1 ; cos ( n ) ; sen ( n ) } n N , [ 0 , 2 ) y p ( )=1 forman una
base ortogonal para el DSTF. Como la funcin que buscamos desarrollar (

sen () ) es uno de los elementos de esa base, el desarrollo en serie
consiste en que todos los coeficientes
correspondiente a la autofuncin. Por lo tanto:
H 0=0 F0 =
H n=0 Fn =
H0
=0
1
ln
2
Hn
n
(( ) ( ) )
1
1
2
2
si n=1
I 1=1
G1 =
I n=
si n>1
I n=0
=0
I1
1
(( 12 ) ( 12 ) )
Gn =
1 2
=
3 3
2
In
n
(( ) ( ) )
1
1
2
2
=0
son
cero,
excepto
el
Remplazamos
2
u ( r , )=
3
F0
r 1 r
2
2
Fn
Gn
en la ecuacin (16), obteniendo:
(( ) ( ) )
sen(1 )
Recordamos que esta es la solucin del problema 1 y entonces

retomamos el subndice correspondiente de la solucin:
( 43 1r 13 r ) sen( )(18)
u1 ( r , )=
SOL. (Problema 2)
Aplicando el mtodo de separacin de variables proponemos que:
u2 ( r , )=R ( r ) ( )
Resolviendo en forma anloga como para
u1
llegamos a las mismas
ecuaciones (6) y (7), y las condiciones de contorno
(I )
' ' ( ) + ()=0

(6)
2 ''
'
r R ( r ) +r R ( r ) R ( r )=0 (7)
( 0 ) = ( 2 ) ( I )
' (0)= ' (2 )(II )
Las condiciones de contorno restantes para
u2
son:
( II) .
u2 (1, )=R ( 1 )
( )=0
R ( 1 ) =0(V )
u2 (2, )=g ( )= ( 2 ) (VI )
Nuevamente la ecuacin (6) y las condiciones
(I)
(II)
representan
un PSLS y la resolucin, igual a la desarrollada en el problema anterior,

nos da:
=0
0 ( ) =B 0
> 0 n ( )= An cos ( n )+ B n sen ( n )
Con
=n2 .
La solucin a la ecuacin (7) para
=0
es:
R0 (r )=C 0 ln(r )+ D0
Aplicamos
(V ) :
R0 ( 1 )=C 0 ln ( 1 ) + D0=0
D0=0
R0 (r )=C 0 ln(r )
La solucin de la ecuacin (7) para
>0
es:
Rn ( r )=C n r n+ D n r n
Aplicamos
(V ) :
Rn ( 1 )=C n 1n + Dn 1n=0
Dn=Cn
Rn ( r )=C n ( r nr n )
Escribimos la solucin general
u2 ( r , )
en forma anloga a lo realizado
en el problema 1:
u2 ( r , )=R0 ( r ) 0 ( )+ Rn ( r ) n ( )
n=1
Reemplazando nos queda:
n
n
u2 ( r , )=B0 C 0 ln (r )+ [ An cos ( n )+ B n sen ( n ) ] Cn ( r r )
n=1
( n n )
u2 ( r , )=B
A n Cn cos ( n )+ B
0 C0 ln (r)+
n C n sen ( n ) r r
n=1
E0
En
Fn
u2 ( r , )=E0 ln(r)+ [ E n cos ( n ) + Fn sen ( n ) ] ( r nr n ) (19)

n=1
Usamos la condicin
Fn
(VI )
para encontrar los coeficientes
E0
En
u2 (2, )= ( 2 ) , 0 < 2
u2 (2, )=E
E n ( 2n 2n ) cos ( n ) +
F n ( 2n2n ) sen (n )
0 ln ( 2 ) +
n=1
G0
Gn
Hn
u2 (2, )=G 0 + G n cos ( n ) + H n sen(n )= ( 2 ) (20)

n=1
La ltima ecuacin puede interpretarse como un DSTF de la funcin

g ( )= ( 2 ) . Por lo tanto a continuacin calculamos los coeficientes
G0
Gn
Hn
. Comenzamos con
G0
1
2 3
( 2 ) d=
2
2
2 3
)|
( 2 )2 ( 2 )3 ( 2 )3 1 1 4 2 2 2
1
=
2
=
=
=
2
2
3
2 2 3
6
3
1
1
G0= g ( ) d=
T 0
2 0
G0=E 0 ln ( 2 )=
E0=
2
3
2 2
(21)
3 ln ( 2 )
A continuacin calculamos
T
Gn :
G n=
2
g ( ) cos ( n ) d= 1 ( 2 ) cos ( n ) d
T 0
0
G n=
1
[ 2 cos ( n ) 2 cos ( n ) ] d
0
Gn=2 cos ( n ) d
0
1
2 cos ( n ) d
0
Resolviendo por tabla:

cos ( n ) sen (n)
Gn=2
+
n
n2
Gn=2
)| (
0
1 2 cos ( n ) 2 2
+ 3 sen( n)
n n
n2
cos ( 2 n ) 2
sen (2 n )
cos ( 0 )
+
n2
n2
) (
)|
cos ( 2 n ) 4 2 2
1 4
+
3 sen
(2 n )
n
n2
n
G n=
1 4 4
= 2
n2
n
( )
n
Gn=E n ( 2 2 )=
En=
4
n2
4
(22)
n ( 2n2n)
2
Hn
Finalmente calculamos
T
2
1
H n= g ( ) sen ( n ) d= ( 2 ) sen ( n ) d
T 0
0
1
H n= [ 2 sen ( n )2 sen ( n ) ] d
0
1
H n=2 sen ( n ) d 2 sen ( n ) d
0
0
Resolvemos por tabla:

sen ( n ) cos(n)
H n=2
n
n2
H n=2
)| (
0
1 2 sen ( n ) 2 2
cos(n)
n n3
n2
sen ( 2 n ) 2
cos(2 n )
) (
)|
sen ( 2 n ) 4 2 2
1 4
2
3 cos
( 2 n ) + 3
cos(0)
2
n
n
n
n
( )
H n=
((
) ( ))
4 1
4 2 2
2
4 4

3 + 3 =
+
=0
n
n
n
n
n
n
H n=F n ( 2n2n )=0
Fn =0(23)
Reemplazando (21), (22) y (23) en la ecuacin (19) tenemos:
(4 )
2 2
u2 ( r , )=
ln(r )+ 2 n n cos ( n ) ( r nrn ) (24)
3 ln (2 )
n=1 n ( 2 2 )
Solucin general del problema

Finalmente utilizando las expresiones (18) y (24) obtenemos la solucin
del ejercicio:
u ( r , )=u1 ( r , ) +u2 ( r , )
u ( r , )=
2
(4 )
4 r
2
n
n
sen()+
ln ( r)+ 2 n n cos ( n ) ( r r )
3r 3
3 ln ( 2 )
n=1 n ( 2 2 )
Problema 3
Cuando una cuerda vibrante se somete a una fuerza vertical estndar
que vara con la distancia horizontal desde el extremo izquierdo, la
ecuacin de onda de la forma:
a
d2 U
d2U
+
Ax=
dx2
dt 2
Donde A es una constante

sujeta a:
U (0,t)=0
U(x, 0)=0
resuelva la ecuacin diferencial parcial

;
;
U(1,t)=0
dU
=0
dt t=0
f(x)=0
t>0
t<0
Al no ser homognea se propone un cambio de variable:

a
d2 U
d2U
+
Ax=
2
2
dx
dt
U(x, t)= w(x, t) +(x)

La segunda derivada con respecto de x es:
2
d U d w
= 2 + ' ' (x)
dx2
dx
d2U d2 w
= 2
2
dt
dt
w t ( x , t ) =ut ( x , t )
w t ( x , 0 )=u t ( x , 0 )=0
w ( x ,0 )=U ( x , 0 ) ( x )
Y sabemos que
U ( x , 0 )=0
por lo que:
w ( x ,0 )=0 ( x )
Sustituyendo los valores en la ecuacin original:
d2 w
d2 w
=a
+ ' ' (x ) + Ax
dt 2
dx2
As se obtiene la ecuacin homognea.

d2w
d2w
=a
dt 2
dx 2
Por lo que
2
''
a ( x )+ Ax=0
Encontrando el valor de
''
( x )=
' ( x )=
( x) =
''
(x)
Ax
a2
Ax
+
2a 2
Ax
+x+
2
2 3 a
De:
U(x,t)= w(x,t) +(x)
condicin:
se despeja
w(x,t) y se aplica la primer
w(x,t) = U(x,t) (x)

w (0,t)= 0 = U(0,t) - (0)
Sustituyendo el valor de cero en las ecuaciones obtenidas para U y
respectivamente se obtiene:
W(0,t)= 0 = U(0,t) - (0) = 0 + =0
Por lo tanto =0
Aplicando la segunda condicin:

W(1,t)= 0 = U(1,t) - (1)
Sustituyendo el valor
en las ecuaciones obtenidas para U y
respectivamente se obtiene:
W(1,t)= 0 = U(1,t) - (1) =
A
2
6a
=0
por lo que =
Sustituyendo en la ecuacin encontrada para (x) obtenemos:

( x) =
Ax Ax
+
x
6 a2 6 a2
( x) =
A
( xx 3 )
2
6a
d2w
d2w
=a
dt 2
dx 2
w ( x ,t )=( C1 cos ( x )+ C2 sen ( x ) ) (C 3 cos ( at )+ C4 sen ( at ))
Aplicando la condicin: w(0,t)=0
C1 cos ( 0 )+C 2 sen ( 0 ) =0 ; C1 cos ( 0 ) =0 ; C1=0

Por lo que la ecuacin se reduce:
w ( x ,t )=( C2 sen ( x ) ) (C 3 cos ( at )+C 4 sen ( at ))
Aplicando la condicin:
W t ( x ,0 )=0
w t ( x , t ) =( C 2 sen ( x )) (C 3 sen ( at ) +C 4 cos ( at ))
A
2
6a
C 3 sen ( 0 )+C 4 cos ( 0 )=0 ; C 4 cos ( 0 )=0 ;C 4 =0

w t ( x , 0 ) =
Por lo que la ecuacin se reduce:
w ( x ,t )=C 2 sen ( x ) C 3 cos ( at )
Agrupando las constantes:
w ( x ,t )=bn sen ( x ) cos ( at )
Aplicando la condicin w(L,t)=0
sen ( L) =0
Entonces:
=n
y sabemos que
sen ( n )=0
L= n por lo que
=n/L
; como L=1
entonces:
Sustituyndola en la ecuacin nos queda:

w ( x ,t )=bn sen ( nx ) cos ( nat )
Pero como n puede tomar cualquier valor entonces se convierte en una
sumatoria:
w ( x ,t )= bn sen ( nx ) cos ( nat )

n=1
Aplicando w(x,0)
w ( x ,0 )= b n sen ( nx )=0
n=1
w ( x ,0 )=U ( x , 0 ) ( x ) =0
w ( x ,0 )=
A
( xx 3 )=f ( x )=0
2
6a
Se resuelve la seria de Fourier seno
2
bn = f ( x ) sen ( nx ) dx
L0
Se aplica la ltima condicin de frontera f(x)=0 y se sustituye en valor
de para la serie de Fourier:
L=1
1
bn =
A
( xx ) s en ( nx ) dx
3a 0
A
A
bn =
xsen ( nx ) dx+
x sen ( nx ) dx
3a 0
3a 0
PRIMERA PARTE DE LA INTEGRAL:
1
A
xsen ( nx ) dx
3a 0
Por partes:
u=x
du=dx
dv=sen(nx)
v=-cos(nx)(1/n)
1 1 1
x
cos ( n x ) + cos ( n x)dx
n
n 0
0
1 (1 )n
1
1 2
cos ( n ) +
sen ( n x ) =
n
n
n
0
( )
A (1 ) A (1 )
=
resultado de la primer integral
3a
n
3 an
SEGUNDA PARTE DE LA INTEGRAL:

1
A
x3 sen ( nx ) dx
2
3a 0
Por partes:
u=x
du= 3xdx
dv= sen(nx)
v=-cos(nx)(1/n)
1 3 1
x
cos ( n x ) + x cos ( n x)dx
n
n 0
0
Por partes:
u=x
du=2xdx
dv=cos(nx)
v=sen((nx)(1/n)
n
1 2 1
(1 ) 3 x 2
+
sen ( n x ) x sen ( n x ) dx
n
n n
n 0
0
Por partes:
u=x
du=dx
dv=sen(nx)
((
v=-cos(nx)(1/n)
1 1 1
(1 )n 3 2 x
+
cos ( n x ) cos ( n x ) dx
n
n n n
n 0
0
))
n
(1 )n
6 (1 ) 1
+
n
n
n
(n )
((
1
n
1
sen ( n x )
0
(1 ) 6 (1 )
+
n
(n )
n
A (1 ) 6 (1 )
2
+
n
(n )
3a
A (1 )n 2 A (1 )n
+ 2
resultado de la segunda parte de la integral
3 a2 n
a (n )
A (1 )n A (1 )n 2 A (1 )n
bn =
+
3 a n 3 a2 n a 2( n)
2 A (1 )n
bn = 2
a ( n)
Sustituyendo en la serie de Fourier y aplican la condicin w(x,0):
w ( x ,0 )=
n=1
2 A (1 )
sen ( n )
a2 (n )
w ( x ,0 )=
(1 )
2A
sen ( n )
2
a n=1 n
(1 )
2A
w ( x ,t )= 2
sen ( nx ) cos ( nat )
a n =1 n
w(x,t) = U(x,t) (x)
U(x,t)= w(x,t) +(x) sustituyendo para encontrar el resultado final:

(1 )n
2A
A
U ( x , t )= 2
sen ( nx ) cos ( nat ) + 2 ( x x3 )
a n=1 n
6a
Problema 4
a) En el problema 1, suponga que
a=b=
f ( x )=100 x ( x) . Sin
usar la solucin u(x,y), bosqueje la apariencia que tendra la

superficie en la regin rectangular definida mediante 0 x ,
0 y .
b) Cul es el valor mximo de la temperatura u para

0 y ?
a)
Condiciones de Frontera:
u ( 0, y ) =0
u ( a , y )=0
u ( x , 0 )=0
u ( x , b ) =f (x )
0x ,
2 u 2 u
+
=0
x2 y2
La solucin general de la ecuacin de calor de este caso es:
u ( x , y )=(c 1 cos ( x ) + c2 sen ( x ) )( c 3 e y + c 4 ey )
Aplicando las condiciones

u ( 0, y ) =0, c 1 ( c3 e y + c 4 ey )=0 c 1=0
u ( x , 0 )=0, c2 sen
u ( a , y )=0, c 2 sen ( a )=0 a=n , =
( na x )( c +c )=0 c =c
3
n
x
a
n
y
a
( ( ))( ( ))
u ( x , y )= c 2 sen
u ( x , b ) =f ( x )= b n sen
n=1
senh
u ( x , y )= bn sen
n=1
( na x) senh ( na b)
b)
Suponiendo que
a=b= y
u ( x , b ) =f ( x )= b n sen
n=1
f ( x )=100 x ( x)
( na x) senh ( na b)
Se integra para obtener el valor de

a
bn
2
n
bn = f ( x ) sen
x dx
a0
a
( )
( na x) senh ( na y )
n
a
bn =
asenh
bn =
( na b)
200
n
asenh
b
a
100 x (x ) sen ( na x ) dx
0
[)
200
I1 =
n
asenh
b
a
( )
200
I1 =
n
asenh
b
a
( )
xsen
0
n
n
x dx x 2 sen
x dx
a
a
0
( )
x dx =
x sen ( n
a )
0
x dx=
xsen n
a
0
( )
( )
[
( )
( )
[
( )
( ) ( ) ( )]
200
xa
n
a
n
cos
x + cos
x dx
n
a
n
a
n
asenh
b
a
Para
I2 =
( )
200
xa
n
a
cos
x +
n
a
n
n
asenh
b
a
Evaluando y aplicando la condicin

I1 =
n
sen
x
a
a=
200
cos ( n )
senh( n ) n
I2
200
n
asenh
b
a
( )
200
I2 =
n
asenh
b
a
( )
[
( )
( ) xcos nxa dx ]
n
x dx=
a
200
x 2 a
nx
a
co s
+2
n
a
n
n
asenh
b
a
x dx=
x 2 sen ( n
a )
200
x 2 a
nx
a
cos
+2
n
a
n
n
asenh
b
a
x 2 sen
0
( )
( )
( )[ nxa sen nxa na sen nxa dx ]]
200
I2 =
n
asenh
b
a
n
x sen a x dx=
0
2
( )
( )
Evaluando y aplicando la condicin

I2 =
[
( )
200
x 2 a
nx
a 2
nx 2 a3
nx
cos
+2 x
sen
+ 3 3 cos
n
a
n
a
a
n
n
asenh
b
a
( )
a=
200
400
400
cos ( n ) 3
cos ( n ) + 3
nsenh (n )
n senh(n )
n senh(n)
Por lo tanto:
bn =I 1 + I 2
bn =
200
200
400
400
cos ( n ) +
cos ( n ) 3
cos ( n )+ 3
senh( n )n
nsenh(n )
n senh(n )
n senh (n )
bn =
400
[ 1(1)n ]
n senh (n )
3
Por lo tanto u(x,y) es:
u ( x , y )=
n=1
400
[ 1(1)n ] sen ( nx ) senh ( ny )
n senh(n)
3
Se realiza la derivada parcial con respecto a x y y, para comprobar

en qu punto se obtiene una temperatura mxima de acuerdo a la
grfica de las derivadas parciales.
u x ( x , y )=
n=1
u y ( x , y )=
n=1
400
[ 1(1)n ] cos ( nx ) senh ( ny )
n senh (n )
3
400
[ 1(1)n ] sen ( nx ) cosh ( ny )
n senh( n )
3
X.
APLICACIONES DE LAS EDP LINEALES DE PRIMER

ORDEN:
ECUACIN DE TRANSPORTE
Consideremos un fluido que se mueve con una velocidad constante V en
un tubo recto, fino y con seccin transversal A. Supongamos que el
fluido contiene un contaminante cuya concentracin en el punto x y en
el instante t denotaremos por u(x, t). Para simplificar supondremos
que no hay otras fuentes de contaminantes en el tubo y que el
contaminante no puede escapar a travs de sus paredes. Entonces, en
el instante t, la cantidad de contaminante en la seccin del tubo entre
las posiciones x1 y x2 viene dada por
x2
u(x ,t ) Adx
x1
Por otro lado, podemos expresar la cantidad de contaminante que fluye

a travs de un plano situado en el punto x, durante el intervalo de
tiempo [t1, t2], como
t2
u(x ,t ) AVdt
t1
Estamos ahora en condiciones de hacer el siguiente balance: la cantidad

de contaminante en la seccin [x1, x2] en el instante t2 es igual a la
cantidad en dicha seccin en el instante anterior t1 ms la cantidad que
entro a travs del plano en la posicin x1 durante el intervalo de tiempo
[t1, t2] tenemos la cantidad que sali a travs del plano en la posicin x2
durante el mismo intervalo de tiempo. Esto es,
t2
t2
u(x ,t ) Adx+ u(x

1
t1
,t )
AVdt u(x , t) AVdt

t1
x2
x2
u(x ,t ) Adx=
2
x1
x1
Por el Teorema Fundamental del Clculo sabemos que
x2
x2
t 2 ; x2
u(x ,t ) Adx u(x ,t ) Adx= u t (x ,t ) Adxdt

2
x1
t2
x1
t 1 ; x1
t2
x 2 ;t 2
u(x ,t ) AVdt u(x , t) AVdt= ux ( x ,t ) AVdtdx

t1
t1
x 1 ;t 1
Combinando estas identidades llegamos a que
u
(
)
( t x ,t +V ux ( x , t )) Adxdt =0
t 2 ;x 2
t 1 ;x 1
Suponiendo que la igualdad anterior se verifica en cada segmento del

tubo y en cada intervalo de tiempo y que la funcin u(x, t) y sus
derivadas parciales de primer orden son continuas, obtenemos
ut (x , t )+V ux (x ,t )=0.
LA ECUACIN DE TRANSPORTE SIMPLE
Toma la forma,
ut + cu x =0
(a)
Admite como soluciones en R2 a los llamados frentes de onda (travelling

waves"),
u(x ,t )=h(xct )
(b)
Donde se conoce a c como velocidad de propagacin.

Como en el caso de las EDOS, un problema de Cauchy permite
determinar todas las soluciones de la primera ecuacin (a). A tal efecto
es ms sugestivo escribir (a) en la forma,
ut =cu x
(c)
E imaginarse que el valor inicial es toda una funcin de x mientras que

el lugar de los datos inciales es, en vez de un punto, todo el eje x.
TEOREMA (a)
Para cada C1(R) el problema:
u t +cu x =0
u( x ,0 )=( x )
Solo admite
u=(xct )
como solucin.
Ejemplo 1. Aplicacin del mtodo de las caractersticas para resolver el

problema de Cauchy para la ecuacin de transporte
u t +cu x =0,
u( x ,0 )=( x ) ,
1
x R ( C ( R ))
Solucin u(x, t) = (x ct): propagacin de a lo largo de las

caractersticas, onda fija que se mueve a velocidad constante.
Ejemplo 2. Problema de Cauchy para la ecuacin de transporte con

condicin inicial sobre una proyeccin caracterstica:
u t +cu x =0,
u( ct, t )=h(t ) ,
t R (h C 1 ( R ) )
Infinitas o ninguna solucin dependiendo de si h es constante o no,

respectivamente.
Ejemplo 3. El problema del trfico en una autopista (Haberman [1]).

Deduccin del modelo. La simplificacin del modelo con la hiptesis de
que la velocidad depende exclusivamente de la densidad del trfico (x,
t) lleva a la ecuacin:
t + c ( ) x =0
Modelo para cualquier fenmeno fsico en el que se conserva la masa y

el flujo de materia es funcin de la densidad. Aplicacin del mtodo de
las caractersticas en el problema de Cauchy no lineal dado por la
ecuacin de Burgers,
u t +uu x =0,
u( x ,0 )=h( x ) ,
xR
Primer caso, con h una funcin creciente: ondas expansivas. Segundo

caso, con h una funcin decreciente en algn intervalo, por ejemplo
si x 0
1
h( x ) 1x si 0 x 1
0
si x 1
Cortes en las proyecciones caractersticas. Ondas compresivas y onda de

choque: singularidad, no existe solucin clsica globalmente definida.
Discordancia entre la ecuacin y la realidad (el flujo continua, la
ecuacin es demasiado exigente). Concepto de solucin dbil como
solucin de la ley de conservacin:
x2
1
1
u(x ,t ) dx+ u(x ,t ) u(x
t
2
2
x
2
,t )
=0
Desarrollo terico para una ecuacin general tipo
A (u)t + B(u) x =0
Con ley de conservacin asociada
u
A ( ( x , t))dx+ B(u(x ,t ) )B(u(x ,t ))=0
2
x2
t x
1
Si u(x, t) es una solucin regular salvo en una curva regular del plano x
= (t) a travs de la cual presenta una discontinuidad de salto finito, se
ha de cumplir la condicin de choque:
+ ( t )
u
( t )
+ ( t )
( t )
u
u B
B
( t ) =
Donde
x (t ) u (x , t)
( t )=lim
x (t )+ u (x , t)
+ ( t )=lim
En el ejemplo anterior, la velocidad de la curva de choque ha de ser

entonces 1/2, lo que permite determinar la curva (t) = (t+1)/2 y definir
globalmente la solucin dbil
si x < ( t )
u ( x , t )= 1
0 si x (t )
XI.
DISTRIBUCIONES Y ESPACIOS DE SOBOLEV
Problema 1
Demostrar que la desigualdad de Poincar no se cumple en toda la recta
real.
SOL.
Tenemos que probar que no existe ninguna constante positiva C > 0 tal
que
2 dx C || dx , H 1 ( R)
Para ello vamos a construir una sucesin de funciones test n C

tales que
2 dx /| | dx
Dada una funcin test no trivial C
n ( x )=
(R) definimos
( xn )
Entonces
(R)
n2 ( x ) dx= ( nx ) dx= 2( y)dy
Mientras que
| n (x)| dx= 12 ( ) ( xn ) dx= 1n ( ) ( y )dy

2
Por lo tanto
n2 ( x ) dx
R
2 ( y)dy
= R
n2 , n
( n ) (x )dx ( ) ( y) dy
R
Problema 2
Comprobar que la desigualdad de Poincar tampoco se cumple en un
intervalo semi-infinito.
SOL: Por traslacin y simetra basta considerar el caso del intervalo (0,
) donde la misma prueba del problema anterior se aplica.
Problema 3
Sea = (0, 1) (0, 1) R2 . Demostrar que
2 dxdy | x| dxdy , C 1c ()
SOL. Para cualquier punto (x, y) , como (0, y) tenemos que

(0, y) = 0. Integrando en la variable x en el segmento horizontal que
une (0, y) con (x, y) tenemos entonces
(x , y)= x (s , y) ds
0
Por lo tanto
x
|(x , y)||x (s , y)|dx

0
(|
x (s , y)| dx
| x (s , y)| dx
0
1
2
)
)
1
2
Entonces
|(x , y)| | x (s , y)| dx , (x , y)

0
Integrando en (x, y) obtenemos que
|( x, y)|2 dxdy |x (s , y)|2 dxdy
Problema 4. Consideramos la banda infinita; demostrar que se cumple

la desigualdad del ejercicio anterior.
= (0, 1) Ry R2
SOL. La prueba del ejercicio anterior se aplica sin ninguna alteracin
pues y no juega ms que el papel de un parmetro de integracin.
Problema 5
Consideramos el problema de Neumann
u xx=f , 0< x< 1

u x ( 0 )=u x ( 1 )=0
(1)
Probar que la condicin necesaria para que(1) tenga solucin es

1
que
f 0
0
Probar que la solucin de (1) no puede ser nica.

Probar, usando la desigualdad de Poincar -Wirtinger, que la
solucin de (1) de media cero es nica.
SOL
Integrando la ecuacin de (1) y usando las condiciones de

contorno tenemos
fdx=uxx dx=u x (1)+u x (0)=0

0
En efecto, si u es una solucin de (1) cualquier otra funcin de la

forma u = u + c con c R constante es tambin solucin de(1)
puesto que
uc =(u+ c) =u =f
y que
uc (x)=( u+c ) (x)=u (x)=0, x=( 0, 1 )
Supongamos que u1 y u2 son dos soluciones de (1). Entonces v =

u1 u2 es solucin de
v 2x dx=0
0
Ahora bien, si u1 y u2 son funciones de media nula tambin v lo es por lo

1
que, por la desigualdad de Poincare-Wirtinger, si
v 2x dx=0
0
tambin se
tiene
v 2 dx=0
0
y, por consiguiente, v 0. Es decir, u1 u2 lo cual
garantiza el resultado de unicidad enunciado.
BIBLIOGRAFA
[1] R. Haberman, Ecuaciones en Derivadas Parciales con Series
de Fourier y Problemas de Contorno, Prentice Hall, Madrid,
2003.
[2] I. Peral, Primer Curso de Ecuaciones en Derivadas Parciales,
Addison-Wesley, 1995. (Disponible en la pgina web del profesor
Ireneo Peral)
[3] Y. Pinchover, J. Rubinstein, An Introduction to Partial
Differential Equations, Cambridge University Press, 2005.

Solucion Del Examn Parcial

Cargado por

Copyright:

Formatos disponibles

Solucion Del Examn Parcial

Cargado por

Información del documento

Descripción original:

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Solucion Del Examn Parcial

Cargado por

Copyright:

Formatos disponibles

1.

CASTRO VIDAL RAUL PEDRO

ABANTO LPEZ MIJAIL ANDRS

Muchas veces nos presentamos ante problemas matemticos, fsicos,

DEFINICION DE ECUACION DIFERENCIAL PARCIAL,

Una ecuacin diferencial parcial (EDP) es aquella cuyas incgnitas son

UNIVERSIDAD NACIONAL DEL CALLAOPgina 1

Es una expresin matemtica que contiene una o ms variables

una ecuacin diferencial en derivadas parciales sera:

2. Orden una Ecuacin Diferencial Parcial

UNIVERSIDAD NACIONAL DEL CALLAOPgina 2

Es una PDE de primer orden.

UNIVERSIDAD NACIONAL DEL CALLAOPgina 3

es un entero. La derivada de u con respecto a y se

pero la solucin es:

Esta solucin se aproxima a infinito si

para cualquier valor de y. El problema de Cauchy para la ecuacin de

Las EDP de segundo orden se clasifican habitualmente dentro de cuatro

UNIVERSIDAD NACIONAL DEL CALLAOPgina 4

Con mayor generalidad, si se tiene una ecuacin de segundo orden del

se dice que es elptica si la matriz

Se dice que es parablica si la matriz

Se dice que es hiperblica si la matriz

1. DEFINICIN Y PROPIEDADES FUNDAMENTALES

UNIVERSIDAD NACIONAL DEL CALLAOPgina 5

|ei x f (x)|=|f ( x)| L1 ( Rn )

Gracias a esta propiedad, las EDP lineales con coeficientes constantes se

Tomando la transformada de Fourier obtenemos:

UNIVERSIDAD NACIONAL DEL CALLAOPgina 6

En la obtencin de la ecuacin anterior

hemos utilizado de manera

De (1) Podemos obtener de manera inmediata una expresin para la

Pero, para obtener la solucin necesitamos invertir de la transformada

est bien definida y es una funcin de

UNIVERSIDAD NACIONAL DEL CALLAOPgina 7

Vamos a calcular la siguiente integral:

Donde es el contorno Im (z) = /2t del plano complejo. Deformando los

De las dos ltimas ecuaciones obtenemos:

En efecto, dada f Rn, si definimos la funcin

UNIVERSIDAD NACIONAL DEL CALLAOPgina 8

Para comprobar hacemos lo siguiente:

La identidad sobre la transformada de Fourier de la Gaussiana es muy

Tomando la transformada de Fourier en la primera ecuacin

UNIVERSIDAD NACIONAL DEL CALLAOPgina 9

Otra posibilidad es observar que la familia

Es una aproximacin de la identidad cuando t de modo que

Para todo BC (Rn). Entonces

Combinando (***) y (****) obtenemos

UNIVERSIDAD NACIONAL DEL CALLAOPgina 10

Basta por ultimo aplicar la transformada inversa de Fourier para

Esta ltima formula se obtiene de manera fcil gracias a esta otra

En virtud de ^y^ se obtiene:

Sigamos ahora estudiando las propiedades fundamentales de la

Para todo f L 1 (R n) L 2 (R n).

Para todo v, w L1 (Rn). Conviene en primer lugar observar que ambas

UNIVERSIDAD NACIONAL DEL CALLAOPgina 11

Que coincide con la expresin de

Aplicando esta identidad con v(x) = e |x| 2 y usando la expresin de

Dada u L 1 (R n ) L 2 (R n ), definimos v(x) := u(x) y w = u v

Combinando (*) y (**) obtenemos