Métodos Matemáticos Continuación
Métodos Matemáticos Continuación
Métodos Matemáticos Continuación
ya que
∞
f ( z + 2πn) = ∑A e
m = −∞
m
im ( z + 2πn )
∞
= ∑A e
m = −∞
m
imz
∫e
−im ' z
f ( z )dz = ∑ A ∫e
m = −∞
m
− im ' z imz
e dz (2)
0 0
2π 2π 2π
e i ( m− m ') z
∫e ∫ dz = 2π , por lo
i ( m − m ') z
pero dz = = 0 si m ≠ m’ y si m = m’, entonces
0
i ( m − m' ) 0 0
2π
tanto ∫ e i ( m− m ') z dz = 2πδ m.m ' quedando, entonces, la Ec.(2) como
0
2π ∞
∫e
− im ' z
f ( z ) dz = ∑A
m = −∞
m' 2πδ m ,m ' = 2πAm o sea que,
0
2π
1
∫e
−im ' z
Am ' = f ( z )dz (3)
2π 0
En resumen, las Ecs. (1) y (3) constituyen las expresiones básicas para la expansión de
una función periódica en serie de Fourier.
2π
Si ahora hacemos z = z ' donde el periodo de z’ es a, entonces la Ec. (1) llega a ser
a
∞ 2π
im z'
F ( z') = ∑A e
m = −∞
n
a
(1’)
1
EJEMPLO: Si g (k ) = evalúe f(z):
k + a2
2
∞
1
f(z)= ∫k
−∞
2
+ a2
e ikz dk z ∈ R , lo cual se puede evaluar por el teorema del residuo en el
e ikz e iz ( ai )
= 2πi lim (k − ai ) = 2πi
k → ai (k − ai )(k + ai ) 2ai
π
∴ f ( z) = e − az , z > 0
a
π
Si z<0, entonces el contorno en el semiplano inferior y f ( z ) = − e az
a
I.2 Transformada de Laplace
Si F (t ) = f (t )e − ct H (t ) (6)
1 para t > 0
donde H (t ) =
0 para t < 0
entonces, su transformada de Fourier, por la Ec. (5), es (haciendo, por conveniencia,
G(ω) = 2πg(ω)
∞ ∞
G ( w) = ∫ F (t )e −iwt
dt = ∫ f (t )e
− ct
H (t )e −iwt dt
−∞ −∞
∞
= ∫ f (t )e −t ( c+iw) dt
0
1
Ejemplo: Si G ( s) = evalúe f(t):
( s + a) 2 + b 2
c + i∞
1 1
H (t ) f (t ) = ∫
2π .i c−i∞ ( s + a ) + b
2 2
e ts ds a, b, t > 0
c
2
+b 2
e tz e tz
lim ( z + a + ib) + lim ( z + a − ib)
z → − a −ib ( z + a + ib)( z + a − ib) z →− a+ib ( z + a + ib)( z + a − ib)
− at −ibt
e e e −at e ibt e − at sin bt
= + =
− 2ib 2ib b
c + i∞
1 1 e − at sin bt
2π .i c−∫i∞ ( s + a ) 2 + b 2
Por lo tanto e ds =
ts
H (t )
b
df
I.3 Transformada de Laplace de la derivada: f’ (t)=
dt
d2 f
I.4 La transformada de Laplace de la segunda derivada : = f ' ' (t ) es
dt 2
(12)
t
I.5 Transformada de Laplace de la integral ∫ f (t ' )dt '
0
Donde en el cambio de los límites de integración hemos usado la identidad en las áreas
de integración como se muestra en la figura:
dI 1
t
∴ RI (t ) + L + Q0 + ∫ I (t ' )dt ' = ξ 0 H (t )
dt C 0
1 Q0 1 I ( s ) ξ 0
RI ( s ) + L[SI ( s ) − I (0)] + + = (13)
C s C s s
∞
1
Nótese que £(1)= ∫ e − st dt =
0
s
R 1 R2
Haciendo a ≡ ; b≡ − 2 , entonces, usando (10):
2L LC 4 L
c + i∞
1
= ∫ I (s)e
ts
I (t ) H (t ) ds
2π .i c−i∞
Q
ξ 0 − 0 − at
= c e sin bt
L b
Q0
ξ0 − 1 R 2
R
c 1 − t
∴ I (t ) H (t ) = e 2L
sin − t
L 1 R2 LC 4 L2
− 2
LC 4 L
donde hemos usado el resultado del ejemplo que sigue la Ec. (10).
Otros ejemplos: Ver problema 36 del capitulo 7 del texto de Spiegel o, en este mismo
texto, el problema 117.
II. DESARROLLO DE LAGRANGE
Usando el teorema del argumento: Si f(z) es analítica dentro y sobre una curva simple
cerrada C excepto para un numero finito de polos dentro de C, entonces
1 f '( z)
2π .i ∫
dz = N − P (1)
C
f ( z)
g(z) es analítica dentro y sobre C, f(z) tiene ceros en ar (r=1…l) y polos en br(r=1…k)
Usando la Ec. (2) para g(z)=z y f(z)=z-a-ξΦ(z) que solo tiene un cero simple (n=1) en
z=a+ ξΦ(z) dentro de C, entonces
1 1 − ξφ ' ( z ' )
2π .i ∫ z '− a − ξφ ( z ' )
z' dz ' = z (3)
C
ξ n d n−1
[
]
∞
z =a+∑ n−1 Φ ( z )
n
(4)
n =1 n! dz z =a
Una generalización a este desarrollo es
ξ n d n−1
[
]
∞
F ( z ) = F (a ) + ∑ n−1 F ' ( z )Φ ( z )
n
(5)
n =1 n! dz z =a
Probar Ec. (5):
Usando g(z)=F(z) y f(z)=z-a-ξΦ(z) en (2):
1 1 − ξφ ' ( z ' ) 1F (z') 1 − ξφ ' ( z ' )
2π .i ∫ 2π .i ∫
F ( z) = F (z') dz ' = ⋅ dz '
C
z '− a − ξφ ( z ' ) C
z '−a 1 − ξφ ( z ' ) ( z '− a)
1 F ( z' ) ∞ ξ nφ n ( z ' )
( )∑
2π .i C∫ z '− a
= 1 − ξφ ' ( z ' ) n
dz '
n=0 ( z '− a)
1 F ( z' ) ∞
ξ n F ( z ' )Φ n ( z ' ) F ( z ' )φ n−1 ( z' )φ ' ( z ' )
dz '+ ∑
2π .i C∫ z '− a ∫
= − dz '
n =1 2π .i C ( z '− a )
n +1
( z − a) n
∞
ξ n F ( z' ) d φ n ( z ' )
= F (a) − ∑ ∫ n dz ' ( z'−a) n dz'
n =1 2π .i C
∞
ξ n dF ( z' ) φ n ( z' )
= F ( a) + ∑ ∫ dz ⋅ ( z '−a) n dz'
n =1 2π .in C
ξ n d n−1
[
]
∞
= F (a) + ∑ n−1 ⋅ F ' ( z )φ n ( z )
n =1 n! dz z =a
(5)
Ejemplo 1: Sea z = 1 + ζz p donde φ ( z ) = z p , a = 1 → φ (1) = 1,
d 2 d 2p
dz φ ( z ) = dz z = 2 p z
2 p −1
[ ]z =1 = 2p
z =1 z =1
d2 3 d 2 3p
2 φ ( z ) = 2 z = 3 p(3 p − 1) z
3 p −2
[ ]z =1 = 3 p3 p − 1) , etc….
dz z =1 dz z =1
reemplazando en la ec. (4), obtenemos:
ζ ζ2 ζ3 ζ4
z = 1+ + 2p + 3 p(3 p − 1) + 4 p (4 p − 1)(4 p − 2) + ......
1! 2! 3! 4!
para el caso p = 12 , ζ = 1 → z = 1 + 1 + 1 / 2 +
3 3
2 2 ( − 1)
+ ... ≈ 2.62
6
Por otro lado, el cálculo exacto puede obtenerse de z = 1 + z1 / 2 , o sea,
( z − 1) 2 = z → z 2 − 3 z + 1 = 0 → z = 32 ± 9
4 − 1 = 32 ± 2
5
≈ 2.62 (tomando la raíz mayor)
z − a − 12 ζ ( z 2 − 1) = 0 → z 2 − 2 z / ζ − 1 + 2a / ζ = 0 → z = ζ1 ± 1
ζ2
− 2ζa + 1 = ζ1 ± ζ1 1 − 2aζ + ζ 2
1 − zζ = 1 − 2aζ + ζ 2 .
1 1
Consideremos ahora F ( z ) = =
1 − 2 aζ + ζ 2 1 − zζ
1 1 1 1 ∞
( 1
ζ ( z '− 1) n )
=
2 π .i ∫ . ( z '− a ) (1 − 1 / 2 ζ ( z '2 − 1)
) dz ' = 2π .i ∫ ∑
( z '− a ) n = 0
. 2
( z '− a ) n
dz '
C ( z '− a ) C
1 dz ' ∞
1 2 ζ ( z ' −1)
1 2 n
( ∞
)
ζ n 1 dn 2 n
= ∫ + ∑ ∫
2π .i C ( z '− a) n=1 2πi C ( z '− a) n +1
= 1 + ∑ n n (a − 1) ,
n =1 2 n! da
1 dn 2
Pn (t ) = n n
(t − 1) n , siendo P0 (t ) = 1 . Por consiguiente,
2 n! dt
∞
1
= ∑ ζ n Pn (a ) (7)
1 − 2aζ + ζ 2 n =0
1
Aplicar este resultado para la expansión de r r cuando r>>r’. Note que
r − r'
r r
r − r ' = r 2 + r '2 −2rr ' cosθ = r 1 − 2( rr' )cosθ + ( rr' ) .
2
Por lo tanto, identificando
1 1 1 ∞
= ∑ ( rr' ) Pn (cosθ )
n
r r = (8)
r − r ' r 1 − 2( r ' ) cosθ + ( r ' )2 r n=0
r r
1 1 ( z 2 − 1) n 1 dn 2
Pn (t ) = ∫
2 2πi C ( z − t )
n n +1
dz = n
2 n! dt n
(t − 1) n (9)
1 1
∫
2π [(t 2
− 1) + 2t t 2 − 1l iθ + (t 2 − 1)l 2lθ l −inθ dθ ]
n
( t −1)
= n
2 2π 0 2
n
1 1
∫
2π [(t 2
− 1)l −lθ + 2t t 2 − 1 + (t 2 − 1)l lθ dθ ]
n
( t − 1)
= n
2 2π 0 2
n
1 1
∫
2π [2(t 2
− 1) cosθ + 2t t 2 − 1 dθ ] n
( t − 1)
= n
2 2π 0 2
n
1 2π
=
2π ∫
0
(t + t 2 − 1 cosθ ) n dθ (10)
Consideremos ahora
∞
∑ζ
n =0
n
Pn (t ) ≡ F (ζ , t ) (11)
1
lo que es idéntica a la Ec.(7). La función F (ζ , t ) = es llamada la función
1 − 2tζ − ζ 2
Agrupando
∞ ∞ ∞
Cambiando en el segundo término n+1 por n, y en el tercer término n-1 por, e igualando
los coeficientes hn a cero por ser linealmente independientes, obtenemos:
Una segunda relación de recurrencia se obtiene derivando la función generatriz (14) con
respecto a z:
∂F
(1 − 2hz − h 2 ) = hF
∂z
lo cual da
∞ ∞
(1 − 2hz − h 2 )∑ h n Pn' ( z ) = h ∑ h n Pn ( z )
0 0
ó
∞ ∞ ∞ ∞
∑h 0
n
Pn' ( z ) − 2 z ∑ h n+1 Pn' ( z ) + ∑ h n+ 2 Pn' ( z ) = ∑ h n+1 Pn ( z )
0 0 0