Métodos Matemáticos Continuación

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

I.

TRANSFORMADAS DE FOURIER Y DE LAPLACE


Antes se tenía, que si una función f(z) es periódica
f ( z + 2πn) = f ( z ) , n entero
entonces,

f ( z) = ∑A e
m = −∞
m
imz
(1)

ya que

f ( z + 2πn) = ∑A e
m = −∞
m
im ( z + 2πn )


= ∑A e
m = −∞
m
imz

por ser e imn 2π = e iN 2π = 1


De la Ec. (1) podemos calcular Am usando
2π ∞ 2π

∫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


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


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.

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’)

y la Ec. (3) es ahora


a 2π
1 −im z '
Am = ∫ e a F ( z ' )dz ' (3’)
a0
y por ser f(z) periódica podemos usar también
a

1 2 −im z
Am = ∫ e a F ( z )dz (3’’)
a −a
2

I. 1 Transformada de Fourier. En el límite en que a→∞, podemos hacer el cambio de



2π ∞
variables en las Ecs. (1) y (3) por
a
n → k; ∑
n = −∞
An → ∫ g (k )dk ; k ∈ R, de tal forma
−∞

que la Ec. (1) se transforma en



f ( z ) = ∫ g (k )e ikz dk , z ∈ R (4)
−∞

y la Ec. (3) en:



1
∫ f ( z )e
−ikz
g (k ) = dz (5)
2π −∞

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

plano complejo k, y usando el lema de Jordan: si z>0 usamos el contorno (semicírculo


de radio R→∞) en el semiplano superior, obteniéndose,

1
∫k
−∞
2
+a 2
e ikz dk = 2πi Re sF (k = ai )

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

haciendo s = c+iw (7)



entonces G ( s ) = ∫ f (t )e −st dt (8)
0

Además, por la Ec. (4) y (7)



1
F (t ) = f (t )e −ct H (t ) = ∫ G(w)e
iwt
dw
2π −∞
c + i∞ (9)
1 ds

t ( s −c )
= G ( s )e
2π c − i∞
i

ó, simplificando el factor común e-ct a ambos lados de la ecuación (9):


c + i∞
1
2π .i c−∫i∞
f (t ) H (t ) = G ( s )e ts ds (10)

Las Ecs. (8) y (10) definen la transformación entre .

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

Use el contorno en el plano z indicado en la figura


1
∫ ( z + a) e t dz = 2π .i ∑ Re sF ( z )
2

c
2
+b 2

Polos simples en z1+a+ib=0 y z2+a-ib=0 cuyos


residuos son:

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

∴ £[f ’(t)]= - f(0)+s£[f(t)] (11)


donde hemos integrado por partes y usando la Ec. (8).
Análogamente, integrando dos veces por parte, demostrar que

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:

i) 0 <t’< ∞ y luego t’→∞ área


superior para la primera integral
ii) 0<t<t y luego t’ →∞ área inferior
para la segunda integral.
Ejemplo: Calcule la corriente I(t) en el circuito de la figura si en t=0 se cierra el
interruptor:
dI Q 1, t > 0
RI + L + = ξ 0 H (t ) , H (t ) = 
dt C 0, t < 0
t
dQ
I= → Q(t ) = Q0 + ∫ I (t ' )dt '
dt 0

dI 1  
t
∴ RI (t ) + L + Q0 + ∫ I (t ' )dt ' = ξ 0 H (t )
dt C  0 

Condición inicial I(0)=0:


Tomando la transformada de Laplace de todos los términos de la ecuación:

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

Despejando I(s) de la ecuación (13):


Q0 Q0
ξ0 − ξ0 −
I ( s) = = C C
 R 1   2 2 
L s 2 + s +  L  s + R  + 1 − R 
 L LC  2 L  LC 4 L2 


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)

Donde N y P son, respectivamente, el numero de ceros y polos dentro de C, contando


l k
multiplicidades: N = ∑ nr P = ∑ p r (ver problema 16, cap. 5 texto de Spiegel)
r =1 r =1

Para la generalización de este teorema, ver problema 90, cap. 5, Spiegel:


l k
1 f '( z)
∫ g ( z) dz = ∑ nr g (ar ) − ∑ p r g (br ) (2)
2π .i C f ( z) r =1 r =1

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

1 z' 1 − ξφ ' ( z ' )


2π .i ∫ z '− a 1 − ξΦ ( z ' ) ( z '−a )
ó z = ⋅ dz '
C
1 z'  ∞ n Φ n ( z' ) 
(1 − ξφ ' ( z' ) )∑ ξ
2π .i C∫ z '− a
= n 
dz '
 n = 0 ( z '− a ) 
1 z' ∞
ξ n  z ' Φ n ( z ' ) z ' Φ n−1 ( z' )Φ ' ( z' ) 
=
2π .i C∫ z '− a
dz ' + ∑ ∫
n =1 2π .i C  ( z '− a )
n +1

( z '− a ) n
dz


ξ n
z' d  Φ ( z')  n
= a−∑ ∫  n 
dz ' , Integrando por partes
n =1 2 π .i C
n dz '  ( z '− a ) 

ξ n
Φ ( z' )
n ∞
ξ n  d n−1 n 
=a+∑ ∫ dw = a + ∑  n Φ ( z )
n =1 2π .in C ( z − a )
n
n =1 n!  dz  z =a
Finalmente,

ξ 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)

Ejemplo 2: Sea z = a + 12 ζ ( z 2 − 1) , ó f ( z ) = z − a − 12 ζ ( z 2 − 1) , cuya derivada es

f ' ( z ) = 1 − ζz , y además tiene que tiene un cero en f ( z ) = 0, esto es

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

de lo cual se obtiene: zζ − 1 = ± 1 − 2aζ + ζ 2 , y tomando la raíz negativa, identificar

1 − zζ = 1 − 2aζ + ζ 2 .

1 1
Consideremos ahora F ( z ) = =
1 − 2 aζ + ζ 2 1 − zζ

Reemplazando en la primera línea de la ecuación (5), se obtiene:


1 1 − ξφ ' ( z ' ) 1 1 1 − z 'ξ
2π .i ∫ 2π .i ∫ (1 − z 'ζ ) z '− a − ξ
F ( z) = F ( z') dz ' = . dz '
C
z '− a − ξφ ( z ' ) C 2 ( z ' −1)
1 2

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 

pero, los polinomios de Legendre están definidos por

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

( rr' ) ≡ ζ , cosθ ≡ a, obtenemos la expansión:

1 1 1 ∞
= ∑ ( rr' ) Pn (cosθ )
n
r r = (8)
r − r ' r 1 − 2( r ' ) cosθ + ( r ' )2 r n=0
r r

OTRA FORMA DE DEMOSTRAR LA ECUACIÓN (7):

Use la representación de los polinomios de Legendre

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)

Sea C un círculo de radio t 2 − 1 centrado en z0 = t , esto es z = t + t 2 − 1l iθ , donde

0 ≤ θ ≤ 2π , dz = t 2 − 1l iθ idθ = i ( z − t )dθ , por consiguiente, la integral en la Ec. (9)


puede escribirse como:
( )
n
 t + t 2 − 1l iθ 2
− 1 t 2 − 1l iθ idθ
1 1 2π  
2 n 2πi ∫0
Pn (t ) =
(t 2
− 1l iθ ) n +1

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)

Usando la Ec. (10), obtenemos


1 2π ∞ 1 2π dθ
F (ζ , t ) = ∫ ∑ζ (t + t 2 − 1 cosθ ) n dθ = ∫
n
(12)
2π 0
n =0 2π 0
1 − ζ (t − t 2 − 1 cosθ )
2π dθ 2π
Para evaluar la última integral, use: ∫ 0 a − b cosθ
=
a 2 − b2
, con lo cual la Ec.(12)

(usando también la Ec. (11)) llega a ser:



1 1
∑ζ n
Pn (t ) = = | (13)
n =0 (1 − tζ ) − ζ (t − 1) 2 2 2
1 − 2t ζ − ζ 2

1
lo que es idéntica a la Ec.(7). La función F (ζ , t ) = es llamada la función
1 − 2tζ − ζ 2

generatriz de los polinomios de Legendre.


Fórmulas de recurrencia:
Usando las Ecs. (11) y (12), con los cambios de variables z=t y ζ=h, obtenemos

1
∑h n
Pn ( z ) = ≡ F (h, z ) (14)
n =0 1 − 2 zh − h 2
Derivando con respecto a h:
∂F (h, z ) z−h
= F ( h, z )
∂h 1 − 2hz + h 2
Por consiguiente,
∂F (h, z )
(1 − 2hz − h 2 ) = ( z − h) F
∂h
ó
∞ ∞
( z − h)∑ h n Pn ( z) = (1 − 2hz − h 2 )∑ nh n −1 Pn ( z )
0 0

Agrupando
∞ ∞ ∞

∑ h n Pn ( z )( z + 2nz) − ∑ h n +1 Pn ( z )(1 + n) − ∑ h n−1 Pn ( z )n = 0


0 0 0

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:

(n + 1) Pn+1 (z) − (2n + 1) zPn ( z ) + nPn−1 (z) = 0 (15)

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

Igualando coeficientes de hn en ambos lados, obtenemos

Pn' ( z ) − 2 zPn'−1 ( z ) + Pn'−2 ( z ) = Pn−1 ( z ) (16)

También podría gustarte