Facultad de Matemática, Astronomı́a y Fı́sica

Universidad Nacional de Córdoba

Métodos Matemáticos de la Fı́sica

Guı́as de Trabajos Prácticos

Compilador:
Pedro A. Pury
Fa.M.A.F. c 1996

1
Presentación:

El presente trabajo de edición de las guı́as de Trabajos Prácticos de Métodos Matemáticos de la Fı́sica, lo
llevé a cabo durante mi participación en los prácticos de dicha materia durante los años 1993 y 1994, en los
cuales fue dictada por el Profesor C. B. Briozzo. Por lo tanto, el orden y el contenido de las sucesivas guı́as
es el correspondiente al programa del curso implementado por el Profesor Briozzo en dichas opotunidades.
Los problemas que componen la presente recopilación fueron extraı́dos, en su mayor parte, del libro de
Jon Mathews and R. L. Walker: Mathematical Methods of Physics (Addison Wesley, Second Edition, 1969).
Considero que esta obra es una muy buena referencia como libro (“de cocina”) para la parte práctica de la
materia; si bien, no cubre las espectativas de rigor, ni de profundidad, que se pretenden para el dictado teórico
en un curso de tercer año de la Licenciatura. Es ası́ que que mi crédito sólo cubre la categorı́a de compilador del
material. Sin embargo, los problemas que integran la última guı́a, corresponden a una selección de enunciados
que elaboré para diferentes fechas de exámenes correspondientes a los años académicos 1993 y 1994. La
exhaustiva discusión de estos problemas con Carlos Briozzo, está en mi recuerdo como una experiencia grata
y estimulante. El motivo de la inclusión de estos problemas adicionales, es el de proporcionar al estudiante
una ejercitación adicional al prepararse para rendir.
Hasta la presente recopilación, no existı́a una edición de la totalidad de los problemas con un formato
uniforme y con cierta unidad de estilo en cuanto a la redacción y a la notación. Creo, sinceramente, que esto
constituye mi pequeña contribución como compilador a mejorar la calidad de este material didáctico. Creo
tambien que la mayor utilidad de este trabajo radica en que se encuentra en un formato digital. El material
fue procesado en su totalidad en LATEX incorporando inclusive las figuras, quedando a disposición de los
docentes de la Facultad los archivos .tex. Considero que de esta manera se simplificará en gran medida el
trabajo de revisión anual de las guı́as por parte de las sucesivas cátedras que estén a cargo de la materia.
Es mi deseo proponer la creación de una base de datos digital para docencia en el ámbito de la Facultad.
Esta se integrarı́a con el aporte voluntario por parte de los docentes de las guı́as de trabajos prácticos y
los exámenes de las diferentes materias. La disponibilidad del material, en una base digital previamente
organizada, simplificarı́a la búsqueda de problemas y contribuirı́a a mejorar sucesivamente el nivel de los
prácticos de todas las materias.
Por otro lado considero que este material es patrimonio intelectual de nuestra Facultad, y como tal deben
ser registrados sus derechos en esta Institución. De esta forma su uso total o parcial en otras facultades de
esta Universidad o en otras instituciones deberá realizarse con el explı́cito reconocimiento de los créditos
correspondientes a la Fa.M.A.F..

Pedro A. Pury

2
 Facultad de Matemática, Astronomı́a y Fı́sica
Universidad Nacional de Córdoba

Métodos Matemáticos de la Fı́sica

Guı́a N◦ 1: Álgebra Lineal

Problema 1: Sean ~a y ~b dos vectores arbitrarios en un espacio vectorial con producto interno. Demuestre
la desigualdad de Cauchy-Schwartz:
(~a · ~a) (~b · ~b) ≥ |~a · ~b|2

Determine cuando es válida la igualdad.

Ayuda: Utilice el vector: ~c = ~a + λ~b y exija que ~c · ~c ≥ 0, ∀λ.
Considere una base ortonormal {~ei } en un espacio vectorial de dimensión infinita. Por cuestiones de
convergencia, el desarrollo de un vector arbitrario en términos de la base:
∞
X
~x = xi ~ei
i=1

puede carecer de significado. No obstante, utilizando la desigualdad de Cauchy-Schwartz, demuestre la desi-

gualdad de Bessel:
n
X
|xi |2 ≤ ~x · ~x
i=1

resultado que es válido para todo n, por grande que este sea.

Problema 2: Siendo A, B y C matrices cuadradas, probar que:

a) (AB)−1 = B−1 A−1

b) ˜
(AB) = B̃Ã
c) (AB)† = B † A†
d) Tr(ABC) = Tr(BCA)
e) det(AB) = det(A) det(B).

Problema 3: Supongamos que las matrices A y B son hermiteanas y las matrices C y D son unitarias.
Demostrar que:
a) C−1 AC es hermiteana
b) C−1 DC es unitaria
c) i(AB − BA) es hermiteana.

3
Problema 4: Sea T : V → V la transformación lineal dada por T (~x) = c~x, donde c es un escalar fijo.
Demostrar que T es unitaria si y sólo si |c| = 1.
Si V es un espacio unidimensional, demostrar que las únicas transformaciones unitarias en V son las des-
criptas anteriormente. En particular, si V es un espacio real unidimensional, existen sólo dos transformaciones
ortogonales, T (x) = x y T (x) = −x.

Problema 5: Las siguientes matrices 2 × 2 en C:

     
0 1 0 −1 1 0
σx =   σy = i   σz =  
1 0 1 0 0 −1

son las matrices de Pauli. Probar que:

P
a) lm klm σl σm = 2iσk . Es decir que puede utilizarse la notación: ~σ × ~σ = 2i~σ , con ~σ = (σx , σy , σz ).
b) σj2 = I, ∀j; donde I es la matriz identidad.
c) El conjunto de matrices de Pauli mas I forman una base para el espacio de matrices 2 × 2.

Problema 6: Dada la transformación de semejanza A0 = S−1 AS, probar que:

a) Tr(A0 ) = Tr(A)
b) det(A0 ) = det(A)
c) Si A es hermiteana y S unitaria entonces A’ es hermiteana
d) Si A es ortogonal y S ortogonal entonces A’ es ortogonal
e) Si A es unitaria y S unitaria entonces A’ es unitaria

Dada una función analı́tica f , es usual utilizar la notación f (A), donde A es una matriz cuadrada, para
representar al desarrollo en serie de Taylor de f “evaluado” en A:
∞
X f (n) (0) n
f (A) = A
n=0
n!

Utilizando esta notación pruebe que:

i) f (A0 ) = S−1 f (A)S

ii) det(eA ) = eTr A .

Problema 7: Si los autovectores ~vi , correspondientes a los autovalores λi , de la matriz A forman una base,
probar que:
Q
a) det(A) = i λi
P
b) Tr(A) = i λi .

Problema 8: Sea U una matriz unitaria y ~x1 , ~x2 dos vectores propios de U pertenecientes a los autovalores
λ1 , λ2 , respectivamente. Demostrar que:

i) |λ1 | = |λ2 | = 1
ii) Si λ1 6= λ2 , x†1 x2 = 0

4
Problema 9: Encontrar los autovalores y autovectores normalizados de la matriz:
 
1 2 3
 
 4 5 6 
 
 
7 8 9

Problema 10: En una dada base {~ei } de un espacio vectorial abstracto, una transformación lineal y un
dado vector de dicho espacio, quedan respectivamente determinados por:
   
2 1 0 1
   
 1 2 0 ,  2 
   

   
0 0 5 3

Encontrar las representaciones matriciales de la transformación y del vector en una nueva base tal que
la antigua queda representada por:
     
1 1 0
     
~e1 =  1  ~e2 =  −1  ~e3 =  0 
     
     
0 0 1

Problema 11: Escribir la matriz A y el vector ~x en la base en la cual A es diagonal.

   
0 −i 0 0 0 1
   
 i 0 0 0 0  a 
   

   
A= 0 0 3 0 0  ~x =  i 
   
   
 0 0 0 1 −i  b 
   

   
0 0 0 i 1 −1

Problema 12: Demostrar cada una de las siguientes proposiciones relativas a una matriz A, n × n, real y
ortogonal:

a) Si λ es un autovalor real de A, entonces λ = 1 o λ = −1.

b) Si λ es un autovalor complejo de A, el complejo conjugado λ∗ también es autovalor de A.

Es decir, los autovalores de A no reales son conjugados a pares.

c) Si n es impar, A tiene por lo menos un autovalor real.

Problema 13: Sea V un espacio euclı́deo real de dimensión n. Una transformación ortogonal T : V → V
con determinante igual a 1 se llama rotación. Si n es impar demostrar que 1 es autovalor para T . Esto prueba
que toda rotación de un espacio de dimensión impar tiene un eje fijo.

5
Problema 14: Las cónicas con centro en el origen de un sistema coordenado obedecen la ecuación:

ax2 + 2hxy + by 2 = k

˜ M ~x = k.
a) Escribir la ecuación en forma matricial: ~x

b) Encuentre la transformación que diagonaliza M e interprete geométricamente.

c) Ejemplifique con 5x2 − 4xy + 2y 2 = 30.

Problema 15: El tensor de inercia de un cuerpo rı́gido en un particular sistema de coordenadas está expre-
sado por la matriz:
 q q 
5 3 3
 q 2 2 4 
q
 3 7 1 
 q2 3 18
 
q 
3 1 13
4 18 6

Encuentre los momentos principales de inercia y las direcciones de los ejes principales. Recuerde que en un
sistema coordenado orientado en la dirección de los ejes principales la matriz asociada al tensor es diagonal.

Problema 16: Demostrar que las oscilaciones unidimensionales del sistema mecánico de la figura están
¨ = A~x; donde:
regidas por la ecuación vectorial: m ~x
 
−(k1 + k2 ) k2 0
 
A= k2 −(k2 + k3 ) k3
 

 
0 k3 −(k3 + k4 )

a) Resolver la ecuación mediante la substitución ~x(t) = eiωt ~y, demostrando que esto conduce a un problema
de autovalores.

b) Suponiendo que k1 = k2 = k3 = k4 = k encontrar las frecuencias de los modos normales de oscilación.

c) Escribir la solución de la ecuación bajo las condiciones iniciales: x1 (0) = x2 (0) = x3 (0) = x0 y ẋ1 (0) =
ẋ2 (0) = ẋ3 (0) = 0.

k1 k2 k3 k4

 
 B B B TT z B B B TT z B B B TT z B B B TT 
B  B  B  B B B B B B B B B
 B B B B B B B B B B B B 
 
m m m

Problema 17: Si T es un operador lineal sobre un espacio vectorial de dimensión finita, diagonalizable,
probar el teorema de Cayley-Hamilton: Si f es el polinomio caracterı́stico de T entonces f (T ) = 0.
Ayuda: Probar que dado un autovector ~x de T de autovalor λ, para todo polinomio q se cumple: q(T ) ~x =
q(λ) ~x.

6
Problema 18: Una transformación lineal T se llama nilpotente si para algún número natural k se cumple
que T k = 0. Probar:

a) Si λ es un autovalor de T entonces λ = 0.

b) λ = 0 es un autovalor de T .

Se concluye entonces que el conjunto de autovalores de un operador nilpotente es exactamente {0}.

Problema 19: Sea V = {p(x), polinomios de grado menor o igual que(n − 1)}. Considere sobre este espacio
el operador lineal derivación D.

a) Probar que D es nilpotente.

b) Construya la forma de Jordan correspondiente a D.

Problema 20: Si A es una matriz n × n cuyo polinomio caracterı́stico resulta: f (x) = (x − c 1 )d1 . . . (x − ck )dk
Pk
con i=1 di = n. ¿Cual es la traza de A?

Problema 21: Sea A una matriz triangular n × n tal que aii 6= ajj , ∀i 6= j. Escribir la forma de Jordan
asociada con A.

Problema 22: Escribir las formas de Jordan de las siguientes matrices:

 
1 1 1 1 1 1 1 1
   
2 1 −1 0 1 0  0 0 0 0 0 0 0 1 
 
   
 0 2 0 1 1 0   0 0 0 0 0 0 0 −1 
   
   
 0 0 2 0 1 0   0 1 1 0 0 0 0 1 
   
   
 0 0 0 2 1 0   0 0 0 1 1 0 0 0 
   
   
 0 0 0 0 2 1   0 1 1 1 1 1 0 1 
   
   
0 0 0 0 0 −1  0 −1 −1 −1 −1 0 1 −1 
 
 
0 0 0 0 0 0 0 0

Fa.M.A.F c 1995

7
 Facultad de Matemática, Astronomı́a y Fı́sica
Universidad Nacional de Córdoba

Métodos Matemáticos de la Fı́sica

Guı́a N◦2: Análisis Tensorial

Problema 1: Si ijk es el sı́mbolo antisimétrico de Levi-Civita:

P
a) Demostrar que i ijk ilm = δjl δkm − δjm δkl .
P
b) Calcular ij ijk ijl
P P
c) Dada la matriz M, calcular: ijk lmn ijk lmn Mil Mjm Mkn

Problema 2: Considere un tensor covariante simétrico arbitrario gij . Mediante gij , podemos definir una
operación de diferenciación covariante. Demostrar que la derivada covariante de g ij definida de esta forma
se anula.
 
 i 
Problema 3: Consideremos el sı́mbolo de Christoffel formado con la ayuda de un tensor covariante
 jk 
simétrico arbitrario glm . Demostrar que:
 
 i  ∂ p
= ln det g
 ik  ∂xk

Problema 4: Se definen las coordenadas α, β, γ mediante las relaciones:

x = βγ, y = αγ, z = αβ

donde x, y, z son las coordenadas cartesianas. Dada una función arbitraria u, calcular ∇ 2 u en función de las
derivadas de u respecto de α, β, γ.

Problema 5: Se definen las coordenadas r, θ, φ mediante las relaciones:

x = r senθ cosφ
y = r senθ senφ
z = r cosθ

donde x, y, z son las coordenadas cartesianas. Dada una función arbitraria u, calcular ∇ 2 u en función de las
derivadas de u respecto de r, θ, φ.

Fa.M.A.F c 1995

8
 Facultad de Matemática, Astronomı́a y Fı́sica
Universidad Nacional de Córdoba

Métodos Matemáticos de la Fı́sica

Guı́a N◦ 3: Ecuaciones Diferenciales Ordinarias

Parte I: Soluciones por Cuadraturas

Problema 1: Verifique que las siguientes ecuaciones diferenciales se reducen a ecuaciones diferenciales de
variables separables:
dy
a) = f (ax + by + c) donde a, b, y c son constantes.
dx
dy y
b) =f
dx x
 
dy a1 x + b 1 y + c 1
c) =f
dx a2 x + b 2 y + c 2
Problema 2: Calcular las soluciones particulares de las siguientes ecuaciones diferenciales de variables
separables, que satisfacen las condiciones especificadas:
dy
a) (1 + ex ) y = ex ; y(x = 0) = 1.
dx
dy  π  π
b) sen x = y ln y; (i) y x = = e, (ii) y x = = 1.
dx 2 2
dy π
c) x3 sen y = 2; y(x → ∞) = .
dx 3
Problema 3: Calcular la solución general de las siguientes ecuaciones diferenciales reducibles a ecuaciones
diferenciales homogéneas:
dy
a) x + y − 2 + (x − y + 4) = 0.
dx
dy
b) x + y + 1 + (2 x + 2 y − 1) = 0.
dx
dy
c) (x2 y 2 − 1) + 2 x y 3 = 0, (Ayuda: Utilizar el cambio de variable y = z α ).
dx
Problema 4: Resolver las siguientes ecuaciones diferenciales utilizando la sustitución: y(x) = u(x) v(x),
dy
a) + 2 y = x2 + 2 x (Ecuación lineal inhomogénea).
dx
dy
b) x2 + 2 x3 y = y 2 (1 + 2 x2 ) (Ecuación de Bernoulli).
dx

9
Problema 5: Resolver las siguientes ecuaciones diferenciales teniendo en cuenta si son diferenciales exactas
o bien si existe un factor integrante:
dy
a) x (2 x2 + y 2 ) + y (x2 + 2 y 2 ) = 0.
dx
dy
b) (x2 + y) − x = 0.
dx
dy
c) (x + y 2 ) − 2 x y = 0.
dx
Problema 6: Determine todas las curvas del plano tales que toda recta tangente a cada curva forma siempre
con los ejes coordenados un triángulo de area A = 2 a2 , siendo a una constante.

Problema 7: Calcular la solución general de las siguientes ecuaciones diferenciales ordinarias:

dy dy
a) x2 + y2 = x y
dx dx
p
dy x 1 + y2
b) = √
dx y 1 + x2
dy a2
c) =
dx (x + y)2
dy 1
d) + y cos x = sen 2x
dx 2
dy
e) (1 − x2 ) − x y = x y2
dx
dy p
f ) 2 x3 = 1 + 1 + 4 x2 y
dx
2
 2
d y dy
g) 2
+ +1=0
dx dx

d2 y
h) = ey
dx2
d2 y dy
i) x (1 − x) +4 + 2y = 0
dx2 dx
dy
j) (1 − x) y 2 − x3 =0
dx
dy
k) x + y + x 4 y 4 ex = 0
dx
dy
l) (1 + x2 ) + y = arctg x
dx
 2
2 dy dy
m) x − 2 (x y − 4) + y 2 = 0 (soluciones general y singular)
dx dx
 2
d2 y dy
n) y 2 − − 6 x y2 = 0
dx dx

10
 2
d2 y

dy dy
o) x4 y + x4 + 3 x3 y −1=0
dx2 dx dx

d2 y
p) x2 − 2y = x
dx2
d3 y d2 y dy
q) − 2 − + 2 y = sen x
dx3 dx2 dx
d4 y d2 y
r) 4
+ 2 2 + y = cos x
dx dx
d2 y dy
s) +3 + 2 y = exp (ex )
dx2 dx
 2 2   2 !3
2 d y dy
t) a = 1+
dx2 dx

d2 y 2 dy 1
u) + +y =
dx2 x dx x

Problema 8: En la activación de una lámina de Indio por un flujo constante de neutrones, el número N de
átomos radiactivos obedece la ecuación:
dN
= λ Ns − λ N
dt
donde Ns es el número constante que se alcanza luego de la “saturación”. Calcular N (t) si N (t = 0) = 0.

Problema 9: Encontrar la solución general de:

d2 y dA dy y(x)
A(x) + + =0
dx2 dx dx A(x)
donde A(x) es una fiunción conocida e y(x) es la incógnita.

Problema 10: Calcular la solución general de la ecuación:

d2 y dy
x +2 + n2 x y = sen(ωx)
dx2 dx
Ayuda: Eliminar el término en derivada primera.

Problema 11: Notar que: y = x es una solución de la ecuación:

d2 y dy
(1 − x) +x − y = (1 − x)2
dx2 dx
si el segundo miembro se iguala a cero. Use este hecho para obtener la solución general de la ecuación dada.

Problema 12: Considere la ecuación diferencial:

d2 y dy
2
+ p(x) + q(x) y = 0
dx dx
en el intervalo a ≤ x ≤ b. Suponga que se conocen dos soluciones, y1 (x) e y2 (x), tales que:

y1 (a) = 0 y2 (a) 6= 0
y1 (b) 6= 0 y2 (b) = 0

11
Dar la solución de la ecuación:
d2 y dy
+ p(x) + q(x) y = f (x)
dx2 dx
que obedece las condiciones y(a) = y(b) = 0, en la forma:
Z b
y(x) = G(x, x0 ) f (x0 ) dx0
a

donde G(x, x0 ), la llamada función de Green, se construye sólo en términos de las soluciones y 1 e y2 y asume
diferentes formas funcionales para x0 < x y x0 > x.
Ilustrar este problema resolviendo:

d2 y
+ k 2 y = f (x); y(a) = y(b) = 0
dx2

Parte II: Soluciones en Series de Potencias

Problema 13: Considere la ecuación diferencial:

d2 y
 
2 dy 2 l(l + 1)
+ + K + − y=0
dx2 x dx x x2

0 ≤ x < ∞, donde l es un entero no negativo.

Encontrar todos los valores de la constante K los cuales generan soluciones finitas en todo el rango de x.
Una ecuación como esta se obtiene al resolver la ecuación de Schrödinger para el átomo de Hidrógeno.

Problema 14: ¿Para que valores de la constante K la ecuación diferencial:

d2 y
 
1 K
− + y=0
dx2 4 x

(0 < x < ∞) tiene una solución no trivial que se anula para x → 0 y x → ∞?

Problema 15: ¿Para que valores de la constante K la ecuación diferencial:

d2 y dy
x − 2x + (K − 3 x) y = 0
dx2 dx

posee una solución acotada en el rango 0 ≤ x < ∞?

Problema 16: Se desea una solución de la ecuación diferencial:

d2 y dy
x +2 + (E − x) y = 0
dx2 dx

tal que y(0) = 1, y(∞) = 0. ¿Para que valores de E es esto posible?

12
Problema 17: La ecuación de Bessel para m = 0 es:
d2 y dy
x2 +x + x2 y = 0
dx2 dx
Una solución de esta ecuación es:
x2 x4
J0 (x) = 1 − + − ···
4 64
Mostrar que existe una segunda solución de la forma:

J0 (x) ln x + A x2 + B x4 + C x6 + · · ·

y encuentre los tres primeros coeficientes A, B, C.

Problema 18: Considere la ecuación diferencial:

d2 y dy
x 2
+ (2 − x) − 2y = 0
dx dx
Dar dos soluciones, una regular y de valor 1 en el origen y la otra de la forma:
1
+ A(x) ln x + B(x)
x
donde A(x) y B(x) son regulares en el origen. Dar los primeros tres términos de los desarrollos en serie de
A(x) y B(x).

Parte III: Sistemas de Ecuaciones Diferenciales

Problema 19: Encontrar los todos los puntos de equilibrio de los siguientes sistemas de ecuaciones diferen-
ciales: 
 
 ẋ = ax− bxy
 ẋ = x − x2 − 2 x y 

ẏ = −c y + d x y
 ẏ = 2 y − 2 y 2 − 3 x y 

 ż = z + x2 + y 2


Problema 20: Considere el sistema de ecuaciones diferenciales:

dx dy
= a x + b y, = cx + dy
dt dt
a) Mostrar que x = 0, y = 0 es el único punto de equilibrio si: ad − bc 6= 0.
b) Mostrar que el sistema posee una lı́nea de puntos de equilibrio si: ad − bc = 0.

Problema 21: Encuentre las soluciones de equilibrio de los siguientes sistemas de ecuaciones diferenciales
y determine si son estables o inestables:
  
 ẋ = x − x3 − x y 2  ẋ = x2 + y 2 − 1  ẋ = tg(x + y)
 ẏ = 2 y − y 5 − x4 y  ẏ = 2 x y  ẏ = x + x3

13
Problema 22: Considere un sistema presa-predador donde el predador posee un medio alternativo de
sustento. Este sistema puede ser modelado por el sistema de ecuaciones diferenciales:

 x˙ = α x (β − x ) + γ x x
1 1 1 1 1 1 1 2
 x˙ = α x (β − x ) − γ x x
2 2 2 2 2 2 1 2

donde x(t) e y(t) son las poblaciones al tiempo t de predadores y presas respectivamente y α i , βi y γi son
constantes.
a) Mostrar que el cambio de coordenadas
 
t
βi yi (t) = xi
α i βi

reduce el sistema de ecuaciones a:


 y˙ = y1 (1 − y1 ) + a1 y1 y2
1
 y˙ = y2 (1 − y2 ) − a2 y1 y2
2

donde
γ 1 β2 γ 2 β1
a1 = ; a2 = .
α 1 β1 α 2 β2

b) ¿Cuales son las poblaciones de equilibrio estables cuando: (i) 0 < a2 < 1, (ii) a2 > 1?

c) Se observa que a1 = 3a2 (a2 es una medida de la agresividad del predador). ¿Cuál es el valor de a2 si el
instinto del predador consiste en maximizar su población de equilibrio estable?

Fa.M.A.F c 1995

14
 Facultad de Matemática, Astronomı́a y Fı́sica
Universidad Nacional de Córdoba

Métodos Matemáticos de la Fı́sica

Guı́a N◦ 4: Distribuciones y Transformadas Integrales

Problema 1: Sea la distribución Tf : C0∞ (<) → <, es decir una funcional lineal continua tal que:
Z ∞
Tf (g) = f (x) g(x) dx; ∀g ∈ C0∞ (<)
−∞

Probar que si f es continua entonces: Soporte(f ) = Soporte(Tf ).

Problema 2: ¿Cómo extenderı́a al concepto de distribución la noción de función par (f (x) = f (−x)) y el
de función impar (f (x) = −f (−x))?. ¿Qué propiedades conserva la extensión?

Problema 3: Sea: 
 x, x ≥ 0
g(x) =
 0, x ≤ 0

Claramente g(x) es continua pero no diferenciable (en el sentido clásico). Encuentre las tres primeras deri-
vadas de g en el sentido distribucional.

Problema 4: La parte principal, en el sentido de Cauchy, de una función:

 
1 1
Z
P (f ) = lı́m f (x) dx
x →0 |x|≥ x

es una distribución. ¿Cómo debe interpretarse la expresión:

 
1 1
lı́m =P − iπδ(x − x0 )?
→0 x − x0 + i x − x0

Problema 5: Desarrolle la función dibujada a continuación en una serie de Fourier.

f (x)
6
1
QQ
 Q f (x + L) = f (x)
 Q
 Q
 Q
b  Q "
b  QQ""-
b
x
− L2 L
2

15
Problema 6: Considere la serie de Fourier para la función:

 +1 0<θ<π
f (θ) =
 −1 π < θ < 2π

f (θ + 2π) = f (θ).
a) A la derecha de θ = 0, la suma de los primeros n términos de la serie lucen como se detalla en la figura a
continuación. Encontrar δn , el “overshoot” del primer máximo.

δn 6
?

-
0 θ

b) Mostrar que lı́mn→∞ δn ≈ 0,18. Este comportamiento es conocido como fenómeno de Gibbs.

Problema 7: Una función f (x) es igual a e−x para 0 < x < 1.

P
a) Desarrollle f (x) como una serie de Fourier de la forma n Bn sen(n π x).
b) Desarrollle f (x) como una serie de Fourier de perı́odo 1.

Problema 8: Un sistema lineal es exitado por una señal periódica f (t), tal que f (t + T ) = f (t). La respuesta
del sistema es tal que una entrada sinusoidal de frecuencia angular ω es multiplicada por (ω 0 /ω)2 , a menos
que ω = 0, en cuyo caso no hay salida. La salida puede ser escrita como:
Z T
1
g(t) = G(t − t0 ) f (t0 ) dt0
T 0

Calcule la función G(t).

Problema 9: ¿De qué función es el desarrollo:

cos 3θ cos 5θ
cos θ + + + ...
9 25

una serie de Fourier?

Problema 10: Encontrar la derivada generalizada de la siguiente función de perı́odo 2π:


1
 (π − x) si 0<x≤π
 2


f (x) = 0 si x=0


 1
 − 2 (π + x) si −π ≤ x < 0

16
 1
P
Verificar que esta función posee el desarrollo en serie de Fourier: n n sen(nx).
P
Comentario: Derivando término a término esta serie obtenemos la serie divergente: n cos(nx). Entonces,
el concepto de función generalizada nos permite adscribir un significado a series divergentes en el sentido
tradicional.

Problema 11: Mostrar que:

b
f (x0 )
Z
f (x) δ(g(x)) dx =
a |g 0 (x0 )|
asumiendo que la ecuación g(x) = 0 posee una sola raiz en el intervalo a < x < b.

Problema 12: Evaluar la siguiente expresión:

Z π Z 2
dx dy δ(sen(x)) δ(x2 − y 2 ).
0 1

Problema 13: Calcular la transformada de Fourier de la función de onda para un electrón 2p en el átomo
de hidrógeno:  
1 r
φ(~x) = p z exp −
32πa5o 2ao
donde ao es el radio de Bohr y z la coordenada rectangular.

Problema 14: Un sistema lineal tiene como respuesta G(ω)e−iωt a una señal de entrada e−iωt , donde ω es
arbitrario. Si la entrada tiene la forma particular:

 0 (t < 0)
f (t) =
 e−λt (t > 0)

donde λ es una constante fija, la salida resultante es:


 0 (t < 0)
F (t) =
 (1 − e−αt )e−λt (t > 0)

donde α es otra constante fija.

a) Calcular G(ω).
b) Calcular la respuesta del sistema a la entrada f (t) = Aδ(t).

Problema 15: Calcular la transformada de Laplace L[f (x)] de la función dibujada a continuación:
f (x)
6
1 ···

-
1 2 3 4 5 6 x

17
Problema 16: Una función f (x) tiene el siguiente desarrollo en serie:
∞
X c n xn
f (x) =
n=0
n!
P∞
Escribir la función g(y) = n=0 cn xn en forma cerrada en términos de f (x).

Problema 17: Utilizando la representación integral:

1 π
Z
Jo (x) = cos(xsenθ) dθ
π 0

calcular la transformada de Laplace de Jo (x).

Problema 18: ¿De qué función es:

1
(s2 + 1)(s − 1)
la transformada de Laplace?

Problema 19: Tres nucleos radioactivos decaen sucesivamente en serie, de forma tal que los números N i (t)
de cada tipo obedecen las ecuaciones:

dN1 dN2 dN3

= −λ1 N1 ; = λ1 N1 − λ2 N2 ; = λ2 N2 − λ3 N3
dt dt dt

Si inicialmente N1 = N , N2 = 0, N3 = n, calcular N3 (t) utilizando transformadas de Laplace.

Fa.M.A.F c 1995

18
 Facultad de Matemática, Astronomı́a y Fı́sica
Universidad Nacional de Córdoba

Métodos Matemáticos de la Fı́sica

Guı́a N◦5: Funciones Especiales

Problema 1: Sea: 
 +1 0<x<1
f (x) =
 −1 −1 < x < 0

Desarrolle f (x) como serie de Polinomios de Legendre Pl (x).

Problema 2: Utilizando una función generatriz, o bien por otro método, evalue la serie:
∞
X xn+1
Pn (x)
n=0
n+1

donde Pn (x) son los polinomios de Legendre.

Problema 3: Evalue Pn0 (x),

a) directamente de la fórmula de Rodrigues:
 n
1 d
n
Pn (x) = x2 − 1 ;
2n n! dx

b) a partir de la función generatriz:

∞
1 X
F (h, z) = √ = hn Pn (z) .
1 − 2hz + h2 n=0

Problema 4: Verifique la siguiente representación para los armónicos esféricos:

 21
dl+m

1 2l + 1 (l − m)!
l
Ylm (Ω) = eimφ (−senθ)m cos2 θ − 1
2l l! 4π (l + m)! d(cosθ) l+m

Problema 5: Los polinomios de Hermite Hn (x) pueden ser definidos por la función generatriz:
∞
2hx−h2
X hn
e = Hn (x)
n=0
n!

a) Escribir la relación de recursión que conecta: Hn−1 (x), Hn (x) y Hn+1 (x).
b) Evaluar: Z ∞
2
e−x /2
Hn (x) dx
−∞

19
Problema 6: Considere las funciones fn (x) definidas por:
∞
X xn
(a) f0 (x) =
n=0
(n!)2

(b) (n + 1)fn+1 = xfn − fn+2

(c) fn0 = fn−1

Encontrar una función generatriz G(x, t) tal que:

∞
X
G(x, t) = fn (x)tn
n=−∞

Problema 7: Sean f (x) y g(x) dos soluciones de la ecuación de Bessel (con el mismo m). Mostrar que su
Wronskiano tiene la forma:
constante
f g 0 − gf 0 =
x

Problema 8: Calcular el Wronskiano de Jm (x) y Ym (x) (m arbitrario).

Problema 9: Mostrar que la definición:

∞
X (−1)r  x m+2r
Jm (x) =
r=0
r! Γ(m + r + 1) 2

de los Jm (x) implica que J−m (x) = (−1)m Jm (x) para todo m entero.

Problema 10: Verificar directamente que:

 1/2
2
J1/2 (x) = senx
πx
 1/2
2
J−1/2 (x) = cosx
πx

Fa.M.A.F c 1995

20
 Facultad de Matemática, Astronomı́a y Fı́sica
Universidad Nacional de Córdoba

Métodos Matemáticos de la Fı́sica

Guı́a N◦ 6: Ecuaciones Diferenciales en Derivadas Parciales

Problema 1: Resuelva las siguientes ecuaciones diferenciales en derivadas parciales de primer orden, consi-
derando el problema de valores iniciales especificado:
∂u ∂u
a) +c =0; u(x, 0) = h(x)
∂t ∂x
∂u ∂u
b) +u =0; u(x, 0) = h(x)
∂t ∂x
En el caso b), ¿Hasta que valor del tiempo t se pueden extender las soluciones?
Ayudas:
i) Especifique el campo de direcciones caracterı́sticas asociado con cada ecuación.
ii) Considere la curva γ(s) = (s, 0, h(s)) en <3 . Obtenga la curva caracterı́stica γ(s, t) que pasa por cada γ(s)
en t = 0, resolviendo la ecuación diferencial ordinaria que determina el campo de direcciones caracterı́sticas.

Problema 2: Considere un cı́rculo de radio a con centro en el origen de coordenadas. Sean (r, φ) las
coordenadas polares y (x, y) las correspondientes coordenadas rectangulares del plano. Calcular la solución
del problema de Dirichlet (interior) para la ecuación de Laplace (∇2 u = 0) con las siguientes condiciones de
contorno:

a) u(r = a) = A

b) u(r = a) = A cosφ

c) u(r = a) = A + B y

d) u(r = a) = A x y

e) u(r = a) = A + B senφ

f ) u(r = a) = A sen2 φ + B cos2 φ; A, B: ctes.

∂u
Problema 3: Considere en el problema anterior condiciones de contorno de Neumann: = f (φ) (interior)
∂n c
para los siguientes casos particulares:

a) f (φ) = A

21
b) f (φ) = A x

c) f (φ) = A (x2 − y 2 )

d) f (φ) = A cosφ + B

e) f (φ) = A senφ + B sen3 φ

Identifique los problemas incorrectamente formulados.

Problema 4: Determinar la distribución estacionaria de la temperatura dentro de la capa esférica a < r < b
si la esfera r = a se mantiene a la temperatura T1 y la esfera r = b a la temperatura T2 .

Problema 5: Calcular la solución de la ecuación de Poisson en el plano:

∇2 u(r, φ) = 1

dentro del cı́rculo de radio r = a si u(r = a) = 0.

Problema 6: Resolver el problema de las oscilaciones transversales propias de una cuerda homogénea de
longitud L si:

a) Los extremos de la cuerda están fijos rı́gidamente.

du
b) Los estremos de la cuerda están libres. Es decir = 0 en los extremos de la cuerda. Esto tiene lugar
dx
cuando los extremos de la cuerda están sujetos mediante anillos (de masas despreciables), que deslizan sin
rozamiento sobre barras paralelas.

c) Un extremo de la cuerda está fijo rı́gidamente y el otro libre.

du
d) Los extremos de la cuerda están fijos elásticamente, es decir que: − h u|x=0,L = 0.
dx x=0,L

e) Un extremo de la cuerda está fijo rı́gidamente y el otro elásticamente.

f ) Un extremo de la cuerda está fijo elásticamente y el otro está libre.

Problema 7: En los instrumentos musicales de cuerda percutida (por ejemplo, el piano) las oscilaciones
transversales de la cuerda tensa se generan mediante un golpe, que le da a la cuerda una velocidad inicial sin
desviación inicial. Si el martillo percutor es plano, rı́gido y tiene un ancho 2δ, el perfil de velocidad inicial
de la cuerda estará dado por:


 vo |x − xo | < δ


∂u
(x, t = 0) =
∂t 
 0 |x − xo | > δ



asumiendo que el martillito golpea a una distancia xo del extremo de la cuerda y que esta tiene un largo L.

a) Determine las oscilaciones de la cuerda u(x, t) para todo t > 0, asumiendo que utt = c2 uxx .

22
b) Para la energı́a de la cuerda oscilando se tiene la expresión:
L
1
Z
ρ c2 u2x + u2t dx


E=
2 0

donde la densidad ρ = M/L, si M es la masa de la cuerda. Calcule la energı́a asociada con cada modo de
oscilación o armónico.

c) Determine en que puntos de la cuerda xo la percusión no excita el m-ésimo modo de oscilación.

d) Bajo la razonable aproximación δ << L determine cuál es el armónico que más resuena (es decir que
posee mayor energı́a) si se percute la cuerda a una distancia xo de uno de sus extremos.

Problema 8: Calcular las oscilaciones propias de una membrana circular,

a) con la frontera fija rı́gidamente.
b) con la frontera libre.

Problema 9: Calcular la frecuencia más baja de oscilación de una onda acústica en una esfera hueca de
∂φ
radio R. La condición de contorno es: = 0 y φ obedece la ecuación de ondas homogénea.
∂r r=R
Problema 10: La cantidad u satisface la ecuación de ondas homogénea dentro de un cilindro de radio a, con
la condición de contorno u = 0 en las paredes del cilindro. En el extremo z = 0 tiene que u = u 0 e−iω0 t ; es
decir, se están enviando ondas a lo largo del tubo con diferentes distribuciones espaciales (modos). Encuentre
la velocidad de fase del modo fundamental como función de la frecuencia ω0 e interprete el resultado para
pequeños valores de ω0 .

Problema 11: Encontrar las frecuencias de oscilación más bajas de un parche de tambor en forma de
√
triángulo isósceles recto de lados a, a, a 2.

Problema 12: Considere un parche de tambor con la forma de un sector circular de radio R y ángulo β.
a) ¿Cuál modo es el ilustrado en la figura?. Dé su respuesta para β = π/2, π, 3/2π.

e
e
e
e  nodo
e β
e
R

b) Esquematice en un gráfico las frecuencias de la primer media docena de modos como funciones de β.

23
Problema 13: La ecuación que describe las ondas elásticas en un medio isotrópico es:

∂ 2~a
(λ + 2µ)∇∇ · ~a − µ∇ ∧ (∇ ∧ ~a) = ρ
∂t2

donde ~a mide los desplazamientos desde el equilibrio, ρ es la densidad, y λ y µ son las constantes elásticas
del medio. Encuentre la frecuencia más baja de oscilación de una esfera elástica isotrópica de radio R, dado
que:
a) ~a(~x) es de la forma f (r)êr , y
b) la condición de contorno en r = R es:

∂ar
λ∇ · ~a + 2µ =0
∂r

Con el fin de calcular un número concreto, considere al final λ = µ.

Problema 14: En el exterior de un cilindro infinitamente largo de radio a la función potencial u(~x, t)
satisface la ecuación de ondas. El cilindro se encuentra dividido longitudinalmente, y sobre su superficie


 u0 e−iω0 t (0 < φ < π)


u=

 −u0 e−iω0 t (π < φ < 2π)



Calcular u(~x, t) fuera del cilindro, si sólo existen ondas salientes para r >> a.

Problema 15: Asuma que la densidad de neutrones n dentro del U235 obedece la ecuación diferencial:

1 ∂n
∇2 n + λn =
κ ∂t

a) Encontrar el radio crı́tico R0 para el cual una esfera de U235 con un radio igual o mayor a dicho valor es
inestable, es decir que densidad de neutrones en su interior crece exponencialmente con el tiempo.
b) Suponga dos hemisferios de U235 tales que se encuentren en el lı́mite de estabilidad. Si se juntan para
formar una esfera, esta resulta inestable, y se cumple que:

n ∼ et/τ

Encontrar la constante de tiempo τ de la explosión resultante.

Problema 16: Una esfera de radio R está uniformemente a temperatura T = 0. En el instante t = 0, es su-
mergida en un baño lı́quido a temperatura constante T0 . Calcular la subsecuente distribución de temperatura
T (~x, t) dentro de la esfera. (Sea κ = conductividad térmica / (densidad × calor especı́fico)).

Problema 17: La temperatura en una esfera homogénea de radio a obedece la ecuación de difusión con
constante κ. Por acciones externas, la temperatura de la superficie de la esfera es forzada a comportarse
según lo ilustrado en la figura. Caltular la temperatura T (t) en el centro de la esfera.

24
 T (r = a)
6

To

-
t
−3 τ −2 τ −τ τ 2τ 3τ 4τ 5τ

Las alteraciones se extienden para t = ±∞.

Problema 18: Encontrar los tres autovalores más pequeños de la ecuación de Schrödinger para una partı́cula
confinada en una caja cilı́ndrica de radio a y altura h.
h̄2 2
− ∇ ψ = Eψ (ψ = 0 sobre las paredes) a ' h
2m

Problema 19: Utilice la transformada coseno con respecto a y para calcular la distribución de temperatura
estacionaria en un sólido semi-infinito (x > 0) cuando la temperatura en la superficie (x = 0) se mantiene
constante e igual a la unidad para −a < y < a, y cero fuera de esa franja.

Problema 20: Un alambre recto de radio a es sumergido en un volumen infinito de lı́quido. Inicialmente el
alambre y el lı́quido tienen temperatura T = 0. En t = 0, el alambre es súbitamente llevado a la temperatura
T0 y mantenido en ella. Encuentre la transformada de Laplace F (r, s) de la distribución de temperatura
T (r, t) resultante en el lı́quido.

Problema 21: Utilice transformadas de Fourier para calcular el movimiento de una cuerda estirada, infinita-
mente larga, con desplazamientos iniciales dados y(x, 0) = φ(x) y velocidad inicial nula. Los desplazamientos
satisfacen la ecuación de ondas homogénea.

Problema 22: Las funciones de Green pueden calcularse para ecuaciones diferenciales homogéneas con
condiciones de contorno inhomogéneas. Para ilustrar esto, considere la ecuación de Helmholtz:

∇2 u(r, θ) + k 2 u(r, θ) = 0

dentro del cı́rculo de radio R, con la condición de contorno dada: u(R, θ) = f (θ). La solución puede escribirse
como: Z 2π
u(~x) = G(~x, θ 0 ) f (θ 0 ) dθ 0
0
Calcule G(~x, θ 0 ).

Fa.M.A.F c 1995

25
 Facultad de Matemática, Astronomı́a y Fı́sica
Universidad Nacional de Córdoba

Métodos Matemáticos de la Fı́sica

Guı́a N◦ 7: Ejemplos de Problemas de Examen

p
Problema 1: Considere una cuerda estirada XXXX
XXX
 XXX
 h X

con sus extremos fijos, la cual en el instante XXX
inicial está sujeta del punto p (ver figura). d L

a) Indique que condición debe satisfacer la distancia d a un extremo fijo de manera tal que, si la cuerda en
reposo se libera de p en un dado instante, en la evolución temporal resultante no se encuentre presente el
enésimo modo normal de oscilación. Demuestre que no es posible eliminar el modo fundamental.

b) Grafique la forma de la cuerda en el instante T̄ = L/c, donde c es la velocidad de propagación de las

perturbaciones en la cuerda. Interprete fı́sicamente el tiempo T̄ .

Problema 2: Considere una cuerda estirada

con un extremo fijo y otro libre (desliza sin
roce sobre un eje). Si la cuerda es sujeta en p
(ver figura), indique que condición debe satis- p

libre
facer la distancia d al extremo fijo de manera
  h
tal que; si la cuerda en reposo se libera de p  d
fijo
en un dado instante, en la evolución temporal L

resultante no se encuentre presente el enésimo

modo normal de oscilación. Demuestre que no
es posible eliminar el modo fundamental.

Problema 3: Considere un cilindro hueco conductor de radio interior a y exterior b el cual se encuentra a
una temperatura T = 0 homogénea. En el instante t = 0 por el interior del cilindro comienza a circular un
fluido con temperatura constante T = T0 .
Asumiendo que la superficie exterior del cilindro se mantiene a temperatura constante T = 0, que el cilindro
es muy largo y que la velocidad del fluido es lo suficientemente grande como para despreciar el tiempo que
emplea en llenarlo, calcule la temperatura en todo punto interior de la pared del cilindro para t 6= 0.

26
Asuma ahora que a es muy grande pero b − a se mantiene finito. Calcule cuál es el tiempo caracterı́stico en
la relajación al estado estacionario del cilindro. Explicite sus expresiones con el mayor detalle posible.

Problema 4: Considere un péndulo ideal que cuelga bajo la acción de la gravedad y es perturbado por una
fuerza tangente a la trayectoria dada por:

F (t) = Fo cos(ω t)

a) Determine en la aproximación de oscilaciones pequeñas el ángulo de desviación del equilibrio utilizando

la transformada de Laplace.

b) Indique para qué valores de ω no se mantendrá válida para todo t la aproximación de pequeñas oscilaciones
y explique porqué.

Problema 5: Considere una membrana elástica circular de radio R sujeta en su periferia, e inmersa en un
medio viscoso que ofrece una resistencia al movimiento proporcional a la velocidad. El movimiento de la
membrana resulta ası́ regido por la ecuación:

∂u 1 ∂2u
∇2 u − 2 h − 2 2 =0.
∂t c ∂t

Asumiendo que R h c = 5, determine la solución u(~r, t) de la ecuación para todo t > 0 dado que en el instante
inicial la membrana se encuentra en reposo y su forma está dada por:
 r
u(r, t = 0) = A 1 − .
R

Determine explı́citamente los coeficientes de los distintos modos presentes en términos de las constantes del
problema; e interprete que clase de evolución temporal satisface cada uno de dichos modos.

Problema 6: Considere una barra semi-infinita homogénea de sección transversal uniforme, cuya superficie
se encuentra térmicamente aislada, y su extremo se mantiene siempre a temperatura cero. Si se aplica un
pulso de calor a una distancia d del extremo, de manera tal que la distribución de temperatura inicial queda
expresada por T (x, t = 0) = δ(x − d) (asumiendo el origen de coordenadas en el extremo de la barra)
determine:

a) La distribución de temperatura en la barra para todo t > 0.

b) El flujo de calor en el extremo de la barra y en x = d para todo t > 0.

Ayudas: La solución del problema con condición de contorno es igual, para x ≥ 0, a la solución del problema
de la barra infinita (sin contorno) pero con una distribución inicial, la cual contempla un adecuado pulso
imagen a la izquierda del extremo. La solución del problema infinito puede tratarse mediante una adecuada
transformada integral. Verifique que la solución obtenida por este método satisface la condición de contorno
especificada.

27
Problema 7: En los instrumentos musicales de cuerda percutida (por ejemplo, el piano) las oscilaciones
transversales de la cuerda tensa se generan mediante un golpe; el cual le da a la cuerda una velocidad inicial
sin desviación inicial. Si el martillo percutor es plano, rı́gido y tiene un ancho 2δ, el perfil de velocidad inicial
de la cuerda estará dado por: 

 vo |x − xo | < δ


∂u
(x, t = 0) =
∂t 
 0 |x − xo | > δ



asumiendo que el martillito golpea a una distancia xo del extremo de la cuerda y que esta tiene un largo L.

a) Determine las oscilaciones de la cuerda u(x, t) para todo t > 0, asumiendo que utt = c2 uxx .

b) Para la energı́a de la cuerda oscilando se tiene la expresión:

1 L
Z
ρ c2 u2x + u2t dx


E=
2 0
donde la densidad ρ = M/L, si M es la masa de la cuerda. Calcule la energı́a asociada con cada modo de
oscilación o armónico.

c) Determine en que puntos de la cuerda xo la percusión no excita el m-ésimo modo de oscilación.

d) Bajo la razonable aproximación δ << L determine cuál es el armónico que mas resuena (es decir que
posee mayor energı́a) si se percute la cuerda a una distancia xo de uno de sus extremos.

Problema 8: La aceleración tangencial at de un péndulo de masa m y longitud l que cuelga bajo la acción
de la gravedad g, viene dada por:
mat = −m g senα − κ v + F (t)

donde α el ángulo de desviación desde la vertical. El término −κ v, donde κ es constante y v es la velocidad del
péndulo (con su signo), representa la fuerza de fricción del péndulo con el aire. F (t) es una fuerza periódica
externa que perturba el péndulo:
F (t) = Fo cos(ω0 t)

A los fines de no tener oscilaciones sobreamortiguadas asuma que el parámetro adimensional:

 κ 2 l
Λ=
m g
es menor que 1.

a) Escriba en la aproximación de oscilaciones pequeñas la ecuación diferencial que satisface el ángulo de

desviación α como función del tiempo.

b) Utilizando la transformada integral de Laplace, resuelva la ecuación diferencial bajo las condiciones
iniciales: α(t = 0) = 0 y α̇(t = 0) = 0
p
c) Explicite su resultado para el caso en que la frecuencia ω0 = g/l. Simplifique sus expresiones asumiendo
que ∆ << 1 y retenga sólo las contribuciones de primer orden en ∆.

28
i) Determine el tiempo caracterı́stico asociado con la relajación desde la condición inicial a un regimen
estacionario.
ii) Calcule la amplitud y frecuencia del movimiento asintótico resultante.

Problema 9: Un cilindro de radio R, infinitamente largo,

se encuentra semi-sumergido en un lı́quido a temperatura '$
T1
constante T0 , mientras que su parte superior se encuentra R

expuesto a una atmósfera cuya temperatura es constante &%

T0
igual a T1 . El cilindro a su vez, rota sobre su eje con velo-
cidad angular constante Ω.
La interfase entre el aire y el lı́quido es perfectamente plana y el eje del cilindro se encuentra sobre ella.
El lı́quido no moja la superficie del cilindro, y en consecuencia no es arrastrado por este en su rotación. El
cilindro está fabricado con un material homogéneo cuya conductividad térmica es κ, y el contacto térmico
que presenta la superficie del cilindro con ambos medios es perfecto.

a) Determine la distribución de temperatura en todo punto del espacio en el regimen estacionario desde
el sistema de laboratorio (donde está descripto el problema).

b) En particular, calcule explı́citamente la temperatura sobre el eje del cilindro.

Ayudas: En primer lugar considere el problema respecto de un referencial solidario al cilindro. Desde dicho
referencial el cilindro se halla en reposo pero la condición de contorno es ahora una función del tiempo que
debe determinar. Resuelva la ecuación de conducción del calor mediante una transformada de Fourier en el
tiempo, con el fin de evitar la especificación de condiciones iniciales. Para el cálculo de la transformada de
Fourier de la condición de contorno conviene previamente calcular la serie de Fourier correspondiente a la
dependencia temporal. Finalmente transforme las expresiones al sistema de laboratorio. Tenga presente que
3
las funciones de Bessel de argumento con fase 4π pueden expresarse en término de las funciones reales ber
y bei según:
 √ 
Jn i i x = bern (x) + i bein (x)

Problema 10: Considere un parche circular de tambor de radio R inicialmente en reposo al cual se lo
percute con una aguja (el area de contacto es despreciable frente a la del parche).

a) Determine en que puntos del parche la percusión genera oscilaciones en las cuales no está presente el
segundo modo normal de vibración.

b) Determine en que puntos del parche la percusión genera oscilaciones en las cuales no está presente el
cuarto modo normal de vibración.

Considere como primer modo normal al fundamental. Para dar sus respuestas utilice todos los argumentos
geométricos y fı́sicos que considere pertinentes para reducir u omitir cálculos matemáticos.

29
Problema 11: Considere un cilindro de radio R, semi-infinito, que se encuentra rodeado por un lı́quido
a temperatura constante T = 0, mientras que su tapa (en el plano z = 0) está expuesta a una variación
periódica de temperatura dada por la función:

n T ≤ t ≤ n + 21 T


T1 ,







H(t) = (∀n entero).





 T0 , n + 1 T ≤ t ≤ (n + 1) T



2

El cilindro está fabricado con un material homogéneo cuya conductividad térmica es κ, y el contacto térmico
que presenta toda la superficie del cilindro es perfecto.

a) Determine la distribución de temperatura en el interior del cilindro en el regimen estacionario.

b) Calcule el flujo total de calor, promediado en un perı́odo, a través de la tapa plana en z = 0.

c) Calcule el flujo total de calor, promediado en un perı́odo, a través de la superficie lateral del cilindro.

d) Determine la relación entre los resultados de los puntos b) y c). Analice el caso particular T 1 = −T0 .
Explique todas las afirmaciones que realice en este item.

Ayuda: La densidad de flujo de calor (o flujo por unidad de area) a través de una supercie S viene dado
por:
∂T
j=−
∂n
∂
donde: ∂n es la derivada normal a la superfie S. El signo resultante en la densidad de flujo ası́ calculada,
determina el sentido del flujo de calor a través de la superficie con respecto al sentido del versor n̂ nornal a
la superficie.

Problema 12: A requisito de la barra brava de su club de futbol favorito, un famoso luthier italiano debe
construir un corneta cónica cuya frecuencia fundamental sea exáctamente 500 Hz.

a) Determine el espectro de frecuencias correspondiente a los modos normales de vibración de una cavidad
cónica de base circular abierta.

b) Calcule las dimensiones de la corneta que debe fabricar el luthier.

Ayudas: A todo fin práctico la corneta es un cono de base circular, fabricado en un material rı́gido, abierto
en el extremo opuesto al vértice. Soplando por un pequeño orificio practicado en el vértice se inducen ondas
de presión en el aire. Dichas ondas de presión obedecen la ecuación de ondas usual, en la cual la velocidad
m
de propagación es la velocidad del sonido en el aire: 330 . Las condiciones de contorno bajo las cuales debe
s
resolverse la ecuación son tales que:

i) El gradiente de presión normal a las paredes rı́gidas es nulo, y

30
ii) Sobre el extremo abierto la presión debe mantenerse constante igual a la atmosférica. Esta condición
puede simplificarse imponiéndose sobre un casquete esférico con centro en el vértice del cono y radio igual a
la generatriz del mismo. Además reescribiendo la ecuación para la variable adecuada, esta es una condición
dde contorno homogénea.

Problema 13: En muchas ocasiones, son los miembros adultos de una población de presas quienes son prin-
cipalmente atacados por los predadores, mientras que los miembros jóvenes se encuentran mejor protegidos
de los ataques; o bien por su tamaño más pequeño, o bien porque viven en una región inaccesible a los
predadores. Sea x1 el número de presas adultas y x2 el número de presas jóvenes en la población de presas,
e y el número de predadores. Se cumple que:

dx1
= −a1 x1 + a2 x2 − b x1 y
dt

dx2
= n x1 − (a1 + a2 ) x2
dt

dy
= −c y + d x1 y
dt

donde a1 y c son las tasas de mortalidad por unidad de tiempo de la población de presas y de predadores
respectivamente; n es la tasa de natalidad de la población de presas y a2 es la tasa por unidad de tiempo
de promoción de individuos jóvenes a adultos en la población de presas (un individuo se considera adulto
cuando está en condiciones de procrear). Por último, b y d son los coeficientes usuales de acople entre las
poblaciones de presas y predadores en el modelo de Volterra.

a) Determine todos los puntos de equilibrio del sistema.

b) Establezca bajo que condiciones sobre los coeficientes del problema, es posible que dichos puntos de
equilibrio sean estables.

Problema 14: La temperatura de la superficie de la tierra tiene una periodicidad diaria y anual muy
marcada. Con el fin de estudiar el problema de la propagación del calor en un terreno, suponga que la tierra
es un semi–espacio infinito, homogéneo y uniforme caracterizado por la conductividad térmica κ; y que su
superficie es un plano perpendicular al eje coordenado x sobre el cual se tiene una condición de contorno
periódica:
T (x = 0, t) = A cos(ω0 t)

a) Resuelva la ecuación de difusión del calor correspondiente a este problema en forma completa (i.e. espe-
cificando todas las constantes de la solución).
Ayuda: Dado que se tiene un problema sin condiciones iniciales, es conveniente efectuar una transformada
de Fourier en el tiempo de la ecuación y su condición de contorno, antes de proceder a resolverla.

31
b) Se verifica experimentalmente que la amplitud de las oscilaciones resultantes en el terreno disminuye de
manera exponencial con la profundidad; y que existe un retardo entre el tiempo en que ocurre un máximo
(o mı́nimo) en la superficie y el tiempo correspondiente en el que este aparece a una dada profundidad en el
terreno. Especifique la constante de atenuación de la amplitud y la diferencia de fase resultante en la solución
del problema.

c) En la construcción de una bodega de vino es importante lograr la mayor estabilidad posible en la tempe-
ratura durante perı́odos muy prolongados. Suponiendo que la amplitud térmica anual en la superficie de la
tierra es de 19,5o C, calcule la profundidad a partir de la cual debe construirse la bodega para que la ampli-
tud térmica no supere los 2,6o C durante el año. Calcule además el tiempo de retardo entre la temperatura
máxima registrada en la superficie y la correspondiente registrada en la bodega. Considere κ = 0, 004 cm 2 /s.

Problema 15: Dada una ecuación diferencial de segundo orden, lineal, de coeficientes constantes y no
homogénea:
y 00 + p1 y 0 + p2 y = f (x)

donde f (x) es una función periódica (perı́odo 2π) dada, desarrollable en serie de Fourier:
∞
a0 X
f (x) = + (an cos(nx) + bn sen(nx))
2 n=1

es posible que se admitan soluciones periódicas.

a) A fin de estudiar esta posibilidad, proponga como solución un desarrollo en serie de Fourier:
∞
A0 X
y(x) = + (An cos(nx) + Bn sen(nx))
2 n=1

y determine bajo que condiciones de los parámetros de la ecuación, existen tales soluciones.

b) Determine si las siguientes ecuaciones admiten soluciones periódicas, y en caso afirmativo construya dichas
soluciones.
i) y 00 + y = cosx ii) y 00 − y = |senx|

Problema 16: La versión discreta (en la coordenada espacial, pero de tiempo continuo) de la ecuación de
difusión, es la correspondiente a las “caminatas aleatorias simétricas”:
dPn (t)
= Pn+1 (t) + Pn−1 (t) − 2 Pn (t)
dt
con n entero (−∞ < n < ∞).

Construya la solución Pn (t) de esta ecuación con la condición inicial Pn (0) = δn,o .

Ayuda: La solución puede escribirse de la forma: Pn (t) = f (t) Bn (z); con z = z(t) y Bn (z) es una apropiada
función de Bessel.

Fa.M.A.F c 1995

32