28.06.2014 Views

General form of the finite Hankel transform Forma general ... - SciELO

General form of the finite Hankel transform Forma general ... - SciELO

General form of the finite Hankel transform Forma general ... - SciELO

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

Rev. Téc. Ing. Univ. Zulia. Vol. 31, Edición Especial, 70 - 78, 2008<br />

<strong>General</strong> <strong>form</strong> <strong>of</strong> <strong>the</strong> <strong>finite</strong> <strong>Hankel</strong> trans<strong>form</strong><br />

Ovidio Pirela Cristalino 1 y Alfredo Villalobos Mena 2<br />

1 Facultad de Arquitectura y Diseño. 2 Centro de Investigación de Matemática Aplicada (CIMA).<br />

Facultad de Ingeniería, Universidad del Zulia, Apartado 10482. Maracaibo, Venezuela.<br />

Telefax: (58) 0261-7598490/7598806. opirela@hotmail.com,<br />

alfredovillalobosmena@yahoo.es.<br />

Abstract<br />

In this work <strong>the</strong> <strong>the</strong>ory <strong>of</strong> Sturm-Liouville Finite Trans<strong>form</strong>s is used to obtain a <strong>General</strong> Form <strong>of</strong> <strong>the</strong><br />

Finite <strong>Hankel</strong> Trans<strong>form</strong> in <strong>the</strong> interval [a, b]. Differents trans<strong>form</strong>s are obtained, as particular cases. Finally,<br />

one <strong>of</strong> this trans<strong>form</strong>s is applied to solve one boundary value problem.<br />

Key words: Sturm-Liouville <strong>finite</strong> trans<strong>form</strong>s, <strong>finite</strong> <strong>Hankel</strong> trans<strong>form</strong>s.<br />

<strong>Forma</strong> <strong>general</strong> de la trans<strong>form</strong>ada finita de <strong>Hankel</strong><br />

Resumen<br />

En este trabajo de investigación se obtiene una <strong>form</strong>a <strong>general</strong> de la Trans<strong>form</strong>ada Finita de <strong>Hankel</strong><br />

en el intervalo [a, b], usando la teoría de las trans<strong>form</strong>adas finitas de Sturm-Liouville. A partir de los resultados<br />

<strong>general</strong>es se obtienen, como casos particulares, distintas trans<strong>form</strong>adas finitas de <strong>Hankel</strong>. Finalmente<br />

se aplica una de estas trans<strong>form</strong>adas para resolver un problema de contorno.<br />

Palabras clave: Trans<strong>form</strong>adas finitas de Sturm-Liouville, trans<strong>form</strong>adas finitas de <strong>Hankel</strong>.<br />

Introducción<br />

La técnica de las trans<strong>form</strong>adas integrales<br />

es uno de los métodos más poderosos para resolver<br />

problemas en Ingeniería y otras ramas de la<br />

Ciencia, ya que las ecuaciones diferenciales ordinarias<br />

y en derivadas parciales con condiciones<br />

iniciales y de borde, pueden ser resueltas de una<br />

manera directa y sistemática, mediante este método.<br />

Dentro del amplio campo de las trans<strong>form</strong>adas<br />

integrales existen las denominadas trans<strong>form</strong>adas<br />

finitas de Sturm-Liouville [17], las cuales<br />

tienen múltiples aplicaciones. En particular,<br />

las trans<strong>form</strong>adas finitas de <strong>Hankel</strong> y algunas<br />

<strong>general</strong>izaciones de éstas, han sido ampliamente<br />

usadas [2, 3, 7-10, 13, 14, 16, 18] en la resolución<br />

de problemas de contorno. Recientemente,<br />

algunos autores han publicado trabajos relacionados<br />

con el tema aquí tratado [1, 5, 6, 11].<br />

En este trabajo se obtiene una <strong>form</strong>a <strong>general</strong><br />

para la Trans<strong>form</strong>ada Finita de <strong>Hankel</strong> en el<br />

intervalo [ ab, , ] la cual involucra las trans<strong>form</strong>adas<br />

comúnmente usadas y otras trans<strong>form</strong>adas<br />

adicionales de gran utilidad en algunas aplicaciones.<br />

por<br />

Trans<strong>form</strong>adas finitas<br />

de Sturm-Liouville<br />

Sea L el operador diferencial lineal definido<br />

1 d df ( x ) <br />

Lf( x)<br />

qx ( ) r( x) f( x)<br />

px ( ) dx <br />

dx <br />

<br />

(1)<br />

donde qx ( ) C 2 [ ab , ], r( x) C [ a, b]<br />

y qx ( )<br />

px ( ) 0<br />

para x [ a, b].<br />

Rev. Téc. Ing. Univ. Zulia. Vol. 31, Edición Especial, 2008


<strong>Forma</strong> <strong>general</strong> de la trans<strong>form</strong>ada finita de <strong>Hankel</strong> 71<br />

Sean M y N los operadores de contorno definidos<br />

por<br />

Mf a f( a) a f( a), Nf b f( b) b f( b) (2)<br />

1 2<br />

1 2<br />

donde a 1<br />

, a 2<br />

, b 1<br />

y b 2<br />

son constantes conocidas tales<br />

que al menos una de a 1<br />

y a 2<br />

y al menos una de<br />

b 1<br />

y b 2<br />

es distinta de cero.<br />

Se demuestra que si n y Kn ( x)son los valores<br />

propios y las funciones propias del problema<br />

de Sturm-Liouville [4]<br />

( L )<br />

K 0 , MK<br />

NK<br />

0 (3)<br />

entonces la trans<strong>form</strong>ada finita definida por<br />

<br />

<br />

b<br />

T f( x); n f( n) p( x) K ( x) f( x)<br />

dx (4)<br />

a<br />

tiene la propiedad que si f( x) C 1 [ a, b]<br />

, entonces<br />

<br />

<br />

T Lf( x); n f( n)<br />

Mf Nf<br />

(5)<br />

donde<br />

n n n<br />

<br />

1<br />

a2<br />

q( a) Kn( a), ( a2<br />

<br />

0)<br />

n 1<br />

a q( a) K<br />

( a), ( a 0)<br />

y<br />

n<br />

T 1<br />

1<br />

n<br />

1<br />

b2<br />

q( b) K ( b), ( b2<br />

0)<br />

<br />

1<br />

b q( b) K ( b), ( b 0)<br />

<br />

donde<br />

1<br />

n<br />

n<br />

2<br />

n<br />

2<br />

(6)<br />

(7)<br />

La fórmula de inversión correspondiente es<br />

f( n) Kn( x)<br />

f( n); x<br />

fx)<br />

<br />

2<br />

n1<br />

K ( x)<br />

<br />

n<br />

(8)<br />

K ( x) p( x) K ( x)<br />

dx<br />

(9)<br />

n<br />

2<br />

b<br />

2<br />

n<br />

a<br />

Trans<strong>form</strong>ada finita de <strong>Hankel</strong><br />

en el intervalo [a, b]<br />

Las trans<strong>form</strong>adas de <strong>Hankel</strong> [17] corresponden<br />

al operador lineal<br />

2<br />

d 1 d <br />

L <br />

<br />

2<br />

<br />

dx x dx x<br />

2<br />

2<br />

(10)<br />

el cual es de la <strong>form</strong>a (1) con<br />

px ( ) x, qx ( ) x, r( x) 2 2<br />

(11)<br />

x<br />

Tomando <br />

2 , el problema de Sturm-<br />

Liouville (3) queda con<strong>form</strong>ado por la ecuación de<br />

Bessel [12, 19]<br />

2<br />

d K( x) 1 dK( x ) <br />

2<br />

2<br />

<br />

<br />

dx x dx x<br />

2<br />

2<br />

<br />

Kx<br />

( ) 0<br />

<br />

(12)<br />

en el intervalo [ ab, , ] donde 0 a b, junto con<br />

las condiciones de contorno<br />

aKa ( ) aKa ( ) , bKb ( ) bK( b)<br />

(13)<br />

1 2<br />

0<br />

1 2<br />

0<br />

Se deduce [14] que los valores propios n<br />

( n 123 , , , ) de este problema son las raíces positivas<br />

de la ecuación trascendente<br />

<br />

<br />

<br />

2<br />

a b J ( b) Y ( a) J ( a) Y ( b)<br />

1 1 <br />

<br />

<br />

<br />

<br />

a1b2 J( b) Y( a) J( a) Y ( b)<br />

<br />

a2b1 J( b) Y ( a) J( a) Y( b)<br />

<br />

a b J( b) Y( a) J( a) Y( b) 0 (14)<br />

2 2<br />

<br />

y las funciones propias correspondientes vienen<br />

dadas por<br />

<br />

<br />

K ( x) a Y ( a) a Y( a) J ( x)<br />

<br />

<br />

<br />

n 1 n 2 n n n<br />

aJ ( a) a J ( a) Y( x)<br />

(15)<br />

1 n 2 n n n<br />

De acuerdo con (4), se define la <strong>form</strong>a <strong>general</strong><br />

de la trans<strong>form</strong>ada finita de <strong>Hankel</strong> en el intervalo<br />

[a, b] por<br />

<br />

<br />

h f( x); n f ( n)<br />

<br />

b<br />

<br />

a<br />

<br />

<br />

<br />

x a Y ( a) a Y( a) J ( x)<br />

<br />

1 n 2 n n <br />

<br />

aJ ( a) a J ( a) Y( x) f( xdx ) (16)<br />

1 n 2 n n n<br />

Usando (5), se tiene que<br />

2<br />

<br />

( ); <br />

n ( ) n 1<br />

( )<br />

2<br />

( )<br />

<br />

<br />

h f x n f n a f a a f a <br />

<br />

donde<br />

n b f( b) b f ( b)<br />

(17)<br />

1 2<br />

n<br />

Rev. Téc. Ing. Univ. Zulia. Vol. 31, Edición Especial, 2008


72 Pirela y Villalobos<br />

1<br />

a2<br />

aa1<br />

Y( na)<br />

<br />

2<br />

<br />

<br />

2 2 2<br />

1<br />

2 2<br />

2<br />

n<br />

<br />

1<br />

a <br />

2nY( na) J( na)<br />

na <br />

a a <br />

(21)<br />

<br />

( a<br />

2<br />

0)<br />

aJ<br />

1 ( na)<br />

<br />

Casos particulares<br />

a <br />

n<br />

<br />

<br />

2 nJ( na) Y( na) ,<br />

<br />

(18)<br />

1<br />

a1<br />

an<br />

a1Y<br />

( na)<br />

<br />

Tomando valores particulares para las<br />

<br />

a2nY<br />

( <br />

constantes a<br />

<br />

na) J<br />

( na)<br />

1<br />

, a 2<br />

, b 1<br />

y b 2<br />

en las fórmulas <strong>general</strong>es<br />

(14) a (21), se obtienen distintos casos de la<br />

<br />

( a<br />

2<br />

0)<br />

<br />

<br />

aJ<br />

1 ( na)<br />

<br />

trans<strong>form</strong>ada finita de <strong>Hankel</strong> en el intervalo [a,<br />

<br />

a <br />

2 nJ( na) Y( na) ,<br />

b], los cuales se resumen a continuación.<br />

y<br />

Caso 1<br />

1<br />

a<br />

b2<br />

b<br />

a1Y<br />

( na)<br />

<br />

1<br />

b1 1 , a2 b2<br />

0<br />

<br />

a <br />

2 nY<br />

( na) J( nb)<br />

Operadores de contorno:<br />

<br />

( b2<br />

0)<br />

aJ<br />

1 ( na)<br />

<br />

Mf f( a) , Nf f( b) (22)<br />

a <br />

2nJ( na) Y( nb) ,<br />

n<br />

<br />

(19)<br />

1<br />

b1<br />

bn<br />

a1Y<br />

( na)<br />

<br />

Trans<strong>form</strong>ada:<br />

<br />

a2nY<br />

( <br />

na) J<br />

( nb)<br />

b<br />

<br />

( b2 0)<br />

h <br />

1, <br />

f( x); n f1,<br />

( n) <br />

x Y( na) J( nx)<br />

<br />

<br />

<br />

aJ<br />

1 ( na)<br />

<br />

a<br />

<br />

a <br />

2 nJ( na) Y( nb) ,<br />

J <br />

( na) Y( nx) f( x) dx (23)<br />

Según (8), la fórmula de inversión viene dada<br />

por<br />

donde n ( n 12 , , ) son las raíces positivas de la<br />

ecuación<br />

<br />

1<br />

f ( n) Kn( x)<br />

J( nb) Y( na) J( na) Y( nb)<br />

0 (24)<br />

h<br />

f ( n); x<br />

f( x)<br />

<br />

2<br />

(20)<br />

n1<br />

Kn( x)<br />

Trans<strong>form</strong>ada de f( x)<br />

Para la determinación de<br />

2<br />

2<br />

h <br />

1, <br />

f( x); n n<br />

f1,<br />

( n) f( a)<br />

<br />

<br />

b<br />

2<br />

2 J( na)<br />

K x xa Y a a Y a <br />

n( ) 1 ( n ) <br />

2n ( n<br />

) J( n<br />

x)<br />

<br />

a<br />

J<br />

b f ( b ) (25)<br />

( n )<br />

2<br />

aJ a a J a 1 ( n ) <br />

2n ( n ) Y( nx)<br />

dx<br />

Fórmula de inversión:<br />

se utilizan algunas integrales y fórmulas de recurrencia<br />

de las funciones de Bessel [12, 15, 19] y, h1, f1, ( n); x<br />

f( x)<br />

<br />

2<br />

1<br />

<br />

2<br />

luego de un laborioso trabajo [14], se obtiene finalmente<br />

2 2<br />

nJ ( nb) f1, ( n) Y( na) J( nx) J( na) Y( nx)<br />

<br />

<br />

2<br />

2<br />

n1<br />

J ( na) J<br />

( nb)<br />

<br />

2<br />

2 2 aJ ( a a J a<br />

1 n ) <br />

2n ( n<br />

)<br />

(26)<br />

Kn( x)<br />

<br />

2 2<br />

n<br />

bJ 1 ( nb)<br />

b2<br />

<br />

nJ<br />

( n b) 2<br />

donde la suma es tomada sobre todas las raíces<br />

<br />

2 <br />

positivas de la ecuación (24)<br />

<br />

2 2 2<br />

<br />

1<br />

2 2 2 1<br />

<br />

n<br />

<br />

nb b b <br />

Caso 2<br />

a b 0 , a b<br />

1.<br />

1 2 2 1<br />

Rev. Téc. Ing. Univ. Zulia. Vol. 31, Edición Especial, 2008


<strong>Forma</strong> <strong>general</strong> de la trans<strong>form</strong>ada finita de <strong>Hankel</strong> 73<br />

Operadores de contorno:<br />

Mf f ( a) , Nf f ( b) (27)<br />

<br />

Trans<strong>form</strong>ada:<br />

<br />

h f( x); n f ( n)<br />

<br />

b<br />

2, <br />

2,<br />

<br />

<br />

xY( aJ ) ( x) J<br />

( aY ) ( x) f( xdx ) (28)<br />

a<br />

n n n n<br />

donde n ( n 12 , , ) son las raíces positivas de la<br />

ecuación<br />

J ( b) Y( a) J ( a) Y ( b)<br />

0 (29)<br />

n n n n<br />

Trans<strong>form</strong>ada de f( x)<br />

2<br />

2<br />

h f x n<br />

2, <br />

( ); n<br />

f<br />

2,<br />

( n) f( a)<br />

<br />

<br />

2n<br />

J ( na)<br />

J b f (<br />

( )<br />

b ) (30)<br />

Fórmula de inversión:<br />

1<br />

<br />

h2, f 2, ( n); x<br />

f( x)<br />

<br />

2<br />

2<br />

J ( nb) f2<br />

, ( n) Y ( na) J( nx) J( na) Y( nx)<br />

<br />

<br />

2<br />

2<br />

<br />

n1<br />

J ( na<br />

) <br />

<br />

2<br />

1 J<br />

( nb)<br />

2 2<br />

a <br />

<br />

n<br />

2<br />

n<br />

<br />

(31)<br />

donde la suma es tomada sobre todas las raíces<br />

positivas de la ecuación (29).<br />

Caso 3<br />

a b 1 , a b<br />

0,<br />

1 2 2 1<br />

Operadores de contorno:<br />

Mf f ( a) , Nf f ( b) (32)<br />

<br />

Trans<strong>form</strong>ada:<br />

<br />

h f( x); n f ( n)<br />

<br />

b<br />

<br />

a<br />

3, <br />

3,<br />

<br />

<br />

xY ( aJ ) ( x) J ( aY ) ( x) f( xdx ) (33)<br />

n n n n<br />

donde n ( n 12 , , ) son las raíces positivas de la<br />

ecuación<br />

<br />

J ( b) Y ( a) J ( a) Y( b) 0 (34)<br />

n n n n<br />

Trans<strong>form</strong>ada de f( x)<br />

2<br />

2<br />

h f x n<br />

3 , <br />

( ); n<br />

f<br />

3 , ( n) f( a)<br />

<br />

<br />

2 J( na)<br />

J b f (<br />

( )<br />

b )<br />

(35)<br />

<br />

Fórmula de inversión:<br />

1<br />

<br />

h3<br />

, f 3 , ( n); x<br />

f( x)<br />

<br />

2<br />

n<br />

<br />

2 2<br />

nJ ( nb) f3,<br />

( n) Y( na) J( nx) J( na) Y( nx)<br />

<br />

<br />

2<br />

1<br />

2<br />

<br />

n<br />

2<br />

J ( na<br />

) <br />

1<br />

J ( )<br />

2 2 <br />

nb<br />

b <br />

n<br />

n<br />

2<br />

(36)<br />

donde la suma es tomada sobre todas las raíces<br />

positivas de la ecuación (34)<br />

Caso 4<br />

a b 0 , a b<br />

1<br />

1 1 2 2<br />

Operadores de contorno:<br />

Mf f( a) , Nf f( b) (37)<br />

<br />

Trans<strong>form</strong>ada:<br />

<br />

h f( x); n f ( n)<br />

<br />

b<br />

4, <br />

4,<br />

<br />

<br />

xY( aJ ) ( x) J<br />

( aY ) ( x) f( xdx ) (38)<br />

a<br />

n n n n<br />

donde n ( n 12 , , ) son las raíces positivas de la<br />

ecuación<br />

J( b) Y( a) J( a) Y( b) 0 (39)<br />

n n n n<br />

Trans<strong>form</strong>ada de f( x)<br />

2<br />

2<br />

h f x n<br />

4 , <br />

( ); n<br />

f<br />

4 , ( n) f( a)<br />

<br />

<br />

2 J<br />

( na)<br />

<br />

J<br />

b f (<br />

( )<br />

b ) (40)<br />

Fórmula de inversión:<br />

1<br />

<br />

h4<br />

, f 4 , ( n); x f( x)<br />

<br />

2<br />

<br />

n<br />

2<br />

<br />

Rev. Téc. Ing. Univ. Zulia. Vol. 31, Edición Especial, 2008


74 Pirela y Villalobos<br />

2<br />

J ( nb) f4,<br />

( n) Y( na) J( nx) J( na) Y( nx)<br />

<br />

<br />

2<br />

2<br />

2<br />

<br />

2<br />

<br />

n1<br />

J<br />

( na<br />

) 1<br />

J ( ) 1<br />

2 2 <br />

nb<br />

<br />

2 2 <br />

b a <br />

n<br />

n<br />

(41)<br />

donde la suma es tomada sobre todas las raíces<br />

positivas de la ecuación (39)<br />

Caso 5<br />

a h , a b 1 , b 0<br />

1 2 1 2<br />

Operadores de contorno:<br />

Mf hf ( a) f ( a) , Nf f ( b) (42)<br />

<br />

Trans<strong>form</strong>ada:<br />

<br />

h f( x); n f ( n)<br />

<br />

b<br />

5, <br />

5,<br />

<br />

<br />

x hY ( a) Y( a) J ( x)<br />

<br />

a<br />

n n n n<br />

<br />

hJ ( a) J ( a) Y ( x) f ( x)<br />

dx (43)<br />

n n n n<br />

donde n ( n 12 , , ) son las raíces positivas de la<br />

ecuación<br />

<br />

<br />

<br />

<br />

hJ( bY ) ( a) J( aY ) ( b)<br />

<br />

n n n n<br />

J ( b) Y( a) J ( a) Y ( b)<br />

0 (44)<br />

n n n n n<br />

Trans<strong>form</strong>ada de f( x)<br />

2<br />

2<br />

h f x n f n hf a f a <br />

5 , <br />

( ); n<br />

5 , ( ) ( ) ( ) <br />

<br />

4 f( b)<br />

(45)<br />

2<br />

aJ ( aY ) ( b) J( bY ) ( a)<br />

<br />

n n n n<br />

Fórmula de inversión (ver ec. 46 abajo),<br />

donde la suma es tomada sobre todas las raíces<br />

positivas de la ecuación (44).<br />

Caso 6<br />

a h , a b 1 , b 0<br />

1 2 2 1<br />

Operadores de contorno:<br />

Mf hf ( a) f ( a) , Nf f ( b) (47)<br />

<br />

Trans<strong>form</strong>ada:<br />

<br />

h f( x); n f ( n)<br />

<br />

b<br />

6, <br />

6,<br />

<br />

<br />

x hY ( a) Y( a) J ( x)<br />

<br />

a<br />

n n n n<br />

<br />

hJ ( a) J ( a) Y ( x) f ( x)<br />

dx (48)<br />

n n n n<br />

donde n ( n 12 , , ) son las raíces positivas de la<br />

ecuación<br />

<br />

<br />

<br />

<br />

hJ ( bY ) ( a) J( aY ) ( b)<br />

<br />

n n n n<br />

J( b) Y( a) J( a) Y( b) 0 (49)<br />

n n n n n<br />

Trans<strong>form</strong>ada de f( x)<br />

2<br />

2<br />

h f x n f n hf a f a <br />

6 , <br />

( ); n<br />

6 , ( ) ( ) ( ) <br />

<br />

4 f( b)<br />

(50)<br />

2<br />

aJ ( aY ) ( b) J<br />

( bY ) ( a)<br />

<br />

n n n n n<br />

Fórmula de inversión (ver ec. 51 abajo),<br />

donde la suma es tomada sobre todas las raíces<br />

positivas de la ecuación (49)<br />

Caso 7<br />

a b 1 , a 0<br />

, b h<br />

1 2 2 1<br />

Operadores de contorno:<br />

Mf f ( a) , Nf hf ( b) f ( b) (52)<br />

1<br />

<br />

h5<br />

, <br />

f<br />

5 , ( n); x f( x)<br />

<br />

<br />

<br />

<br />

nJ ( nb) f<br />

5 , ( n) hY( na) nY<br />

( na) J(<br />

nx) hJ( na) nJ( na) Y( nx)<br />

f( x)<br />

(46)<br />

2<br />

2<br />

2<br />

n<br />

<br />

1<br />

2 2<br />

2<br />

hJ ( a J a <br />

n ) n ( n<br />

) J ( nb) n<br />

<br />

2<br />

h<br />

<br />

a <br />

1<br />

<br />

2 2 2<br />

h6<br />

, <br />

f<br />

6 , ( n); x f( x)<br />

<br />

<br />

<br />

<br />

nJ ( nb) f<br />

6 , ( n) hY( na) nY<br />

( na) J(<br />

nx) hJ( na) nJ( na) Y( nx)<br />

f( x)<br />

(51)<br />

2<br />

2<br />

2 2 <br />

n1<br />

hJ ( <br />

na ) nJ ( na<br />

)<br />

1<br />

2 2<br />

2<br />

<br />

2 2 <br />

J <br />

2<br />

<br />

( nb)<br />

n<br />

h <br />

b a <br />

2 2 2<br />

Rev. Téc. Ing. Univ. Zulia. Vol. 31, Edición Especial, 2008<br />

n


<strong>Forma</strong> <strong>general</strong> de la trans<strong>form</strong>ada finita de <strong>Hankel</strong> 75<br />

<br />

Trans<strong>form</strong>ada:<br />

<br />

h f( x); n f ( n)<br />

<br />

b<br />

7, <br />

7,<br />

<br />

<br />

a<br />

<br />

xY ( aJ ) ( x) J ( aY ) ( x) f( xdx ) (53)<br />

n n n n<br />

donde n ( n 12 , , ) son las raíces positivas de la<br />

ecuación<br />

<br />

<br />

<br />

<br />

hJ( bY ) ( a) J( aY ) ( b)<br />

<br />

n n n n<br />

J ( b) Y ( a) J ( a) Y( b) 0 (54)<br />

n n n n n<br />

Trans<strong>form</strong>ada de f( x)<br />

2<br />

2<br />

h f x n<br />

7 , <br />

( ); n<br />

f7<br />

, ( n) f( a)<br />

b<br />

<br />

Y ( a) J ( b) J ( a) Y ( b)<br />

<br />

n n n n<br />

<br />

hf ( b) f ( b)<br />

(55)<br />

Fórmula de inversión (ver ec. 56 abajo),<br />

donde la suma es tomada sobre todas las raíces<br />

positivas de la ecuación (54).<br />

Caso 8<br />

a 0 , a b 1<br />

, b h<br />

1 2 2 1<br />

Operadores de contorno:<br />

Mf f ( a) , Nf hf ( b) f ( b) (57)<br />

<br />

Trans<strong>form</strong>ada:<br />

<br />

h f( x); n f ( n)<br />

<br />

b<br />

8, <br />

8,<br />

<br />

<br />

xY( aJ ) ( x) J<br />

( aY ) ( x) f( xdx ) (58)<br />

a<br />

n n n n<br />

donde n ( n 12 , , ) son las raíces positivas de la<br />

ecuación<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

hJ( bY ) ( a) J ( aY ) ( b)<br />

<br />

n n n n<br />

J( b) Y( a) J( a) Y( b) 0 (59)<br />

n n n n n<br />

Trans<strong>form</strong>ada de f( x)<br />

2<br />

2<br />

h f x n<br />

8 , <br />

( ); n<br />

f<br />

8 , ( n) f( a)<br />

bn<br />

<br />

J ( a) Y ( b) J ( b) Y( a)<br />

<br />

n n n n<br />

<br />

hf ( b) f ( b)<br />

(60)<br />

Fórmula de inversión (ver ec. 61 abajo),<br />

donde la suma es tomada sobre todas las raíces<br />

positivas de la ecuación (59)<br />

Caso 9<br />

a h , a b 1<br />

, b h<br />

1 1 2 2 1 2<br />

Operadores de contorno:<br />

Mf h f ( a) f ( a) , Nf h f ( b) f ( b) (62)<br />

<br />

1 2<br />

Trans<strong>form</strong>ada:<br />

<br />

h f( x); n f ( n)<br />

<br />

b<br />

9, <br />

9,<br />

<br />

<br />

x h Y ( a) Y( a) J ( x)<br />

<br />

a<br />

1 n n n n<br />

<br />

hJ ( a) J ( a) Y( x) f( xdx ) (63)<br />

1 n n n n<br />

donde n ( n 12 , , ) son las raíces positivas de la<br />

ecuación<br />

<br />

<br />

<br />

<br />

h h J ( b) Y ( a) J ( a) Y ( b)<br />

<br />

1 2 n n n n<br />

<br />

<br />

<br />

<br />

h1 n J( nb) Y( na) J( na) Y ( nb)<br />

<br />

h2 n J( nb) Y ( na) J( na) Y( nb)<br />

<br />

2<br />

n n n n n<br />

J( b) Y( a) J( a) Y( b) 0 (64)<br />

Trans<strong>form</strong>ada de f( x)<br />

<br />

1<br />

<br />

<br />

<br />

h7<br />

, <br />

f7<br />

, ( n); x f( x)<br />

<br />

<br />

n hJ ( nb) nJ ( nb) f7<br />

, ( n) Y ( na ) J (<br />

nx) J( na) Y( nx)<br />

f( x)<br />

<br />

2<br />

2 <br />

n1<br />

2<br />

J ( na) 2<br />

<br />

2 2<br />

n h<br />

<br />

<br />

2 hJ(<br />

<br />

<br />

nb) nJ<br />

( nb)<br />

b <br />

2 2 2<br />

<br />

<br />

(56)<br />

<br />

2<br />

<br />

2 <br />

hJ b J b f n<br />

1<br />

( n ) n ( n<br />

)<br />

8 , (<br />

) Y<br />

<br />

( na) J( nx) J( na) Y( nx)<br />

h8<br />

, f 8 , ( n); x<br />

f( x)<br />

(61)<br />

2<br />

2<br />

<br />

n 2 2<br />

2 2<br />

2 <br />

1<br />

J<br />

( na)<br />

n<br />

<br />

2<br />

h <br />

<br />

hJ b J b<br />

n n <br />

n 1<br />

<br />

( ) ( )<br />

2 2<br />

b <br />

a <br />

n<br />

Rev. Téc. Ing. Univ. Zulia. Vol. 31, Edición Especial, 2008


76 Pirela y Villalobos<br />

<br />

<br />

h9, <br />

<br />

f( x);<br />

n <br />

2<br />

2<br />

<br />

<br />

n f<br />

9, ( n) h1f( a) f ( a)<br />

<br />

<br />

b h J ( b) Y ( a) J ( a) Y ( b)<br />

<br />

<br />

<br />

<br />

1 n n n n<br />

J ( b) Y( a) J<br />

( a) Y ( b)<br />

n n n n n<br />

<br />

h2 f ( b ) f ( b )<br />

(65)<br />

Fórmula de inversión (ver ec. 66 abajo),<br />

donde la suma es tomada sobre todas las raíces<br />

positivas de la ecuación (64).<br />

Los cuatro primeros casos dados son tratados<br />

separadamente por Sneddon [17], mientras<br />

que aquí resultan como casos especiales de la<br />

<strong>form</strong>a <strong>general</strong> de la trans<strong>form</strong>ada finita de <strong>Hankel</strong>,<br />

al igual que los otros cinco casos restantes.<br />

<br />

Un problema de contorno<br />

Considérese la ecuación diferencial<br />

ur ( , t) 2<br />

ur ( , t) 1 ur ( , t)<br />

<br />

k<br />

2<br />

( a r b,<br />

t<br />

0)<br />

t<br />

r<br />

r r<br />

<br />

(67)<br />

con las condiciones<br />

hu( a, t) ur<br />

( a, t) f ( t)<br />

; ubt ( , ) gt ( );<br />

ur ( , 0) wr ( )<br />

(68)<br />

De acuerdo con el Caso 5 de la sección anterior,<br />

se define la trans<strong>form</strong>ada<br />

h u r t r n<br />

5, 0<br />

( , ); u( n, t)<br />

<br />

b<br />

ru( r, t) hY ( a) Y ( a) J ( r ) <br />

a<br />

0 n n 0 n 0 n<br />

( ) ( ) ( ) <br />

hJ a J<br />

a Y r dr<br />

0 n n 0 n 0 n<br />

h u r t r n<br />

5, 0<br />

( , ); u( n, t)<br />

<br />

b<br />

ru( r, t) hY a Y a <br />

0( n ) n 1( n ) J<br />

0( nr<br />

) <br />

a<br />

hJ ( a) J ( a) Y ( r ) dr<br />

(69)<br />

0 n n 1 n 0 n<br />

donde n ( n 12 , , ) son las raíces positivas de la<br />

ecuación<br />

<br />

<br />

<br />

<br />

hJ ( bY ) ( a) J ( aY ) ( b)<br />

<br />

0 n 0 n 0 n 0 n<br />

J ( b) Y( a) J ( a) Y ( b)<br />

<br />

n 0 n 0 n 0 n 0 n 0<br />

esto es<br />

<br />

<br />

<br />

<br />

hJ ( bY ) ( a) J ( aY ) ( b)<br />

<br />

0 n 0 n 0 n 0 n<br />

J ( a) Y ( b) J ( b) Y ( a)<br />

(70)<br />

n 1 n 0 n 0 n 1 n 0<br />

Aplicando la trans<strong>form</strong>ada (69) en la ecuación<br />

(67), se tiene<br />

ur ( , t) <br />

h<br />

5,<br />

0 ; r n<br />

<br />

t<br />

<br />

2<br />

ur ( , t) 1 ur ( , t) <br />

kh<br />

5,<br />

0<br />

2<br />

; r n<br />

r<br />

r r<br />

<br />

unt<br />

( , ) <br />

2<br />

2<br />

k<br />

unt huat u at<br />

n ( , ) ( , ) r( , ) <br />

t<br />

<br />

<br />

4<br />

<br />

2<br />

aJ ( aY ) ( b) J ( bY ) ( a)<br />

<br />

ubt ( , ) <br />

<br />

<br />

0 n 0 n 0 n 0 n<br />

unt<br />

( , ) <br />

2<br />

2<br />

k nu( n, t) f( t)<br />

<br />

t<br />

<br />

<br />

4<br />

2<br />

<br />

aJ ( aY ) ( b) J ( bY ) ( a)<br />

<br />

gt<br />

0 n 0 n 0 n 0 n<br />

<br />

() <br />

<br />

obteniéndose la ecuación diferencial lineal de<br />

primer orden<br />

unt<br />

( , )<br />

k<br />

knu( n, t) <br />

2<br />

2<br />

f ( t ) <br />

t<br />

<br />

2<br />

aJ ( bY ) ( a) J ( aY ) ( b)<br />

<br />

gt<br />

<br />

0 n 0 n 0 n 0 n<br />

<br />

() <br />

<br />

(71)<br />

La solución <strong>general</strong> de la ecuación (71) viene<br />

dada por [20]<br />

1<br />

<br />

<br />

h 9, f9, ( n); x f( x)<br />

2 2<br />

2<br />

nh 2 J( nb) nJ( nb) f9,<br />

( n) h1Y( na)<br />

nY ( na) J( nx) h1J( na) nJ( na) Y( nx)<br />

<br />

f( x)<br />

<br />

2<br />

2 <br />

2<br />

h1J(<br />

na) nJ( na) <br />

2 2<br />

n h <br />

<br />

h J( nb)<br />

nJ<br />

2 2<br />

<br />

2<br />

<br />

b <br />

2<br />

n<br />

<br />

1<br />

2 2<br />

2<br />

( nb<br />

)<br />

n h<br />

2 1 <br />

a <br />

Rev. Téc. Ing. Univ. Zulia. Vol. 31, Edición Especial, 2008


<strong>Forma</strong> <strong>general</strong> de la trans<strong>form</strong>ada finita de <strong>Hankel</strong> 77<br />

t<br />

2 <br />

k<br />

t<br />

2k<br />

<br />

n<br />

unt ( , ) e<br />

f( x)<br />

<br />

<br />

0<br />

2<br />

aJ ( bY ) ( a) J ( aY ) ( b)<br />

<br />

gx<br />

<br />

0 n 0 n 0 n 0 n<br />

<br />

( ) <br />

<br />

kn<br />

x<br />

e 2 dx c<br />

(72)<br />

Al tomar trans<strong>form</strong>ada en la condición<br />

ur ( , 0) wr ( ), se tiene que<br />

un ( , 0) wn ( )<br />

y, de (72), se deduce que c w( n). Luego:<br />

t<br />

k <br />

2 <br />

kn<br />

t<br />

2<br />

unt ( , ) e<br />

f( x)<br />

<br />

0 <br />

2<br />

<br />

aJ <br />

<br />

0<br />

bY0 a J0 aY0<br />

b gx<br />

<br />

( ) <br />

( n ) ( n ) ( n ) ( n ) <br />

k x<br />

<br />

n<br />

e 2<br />

dx<br />

w( n) <br />

(73)<br />

<br />

Usando la correspondiente fórmula de inversión,<br />

se tiene (ver ec. 74 abajo), donde unt ( , )<br />

está dado por la ecuación (73) y la suma es tomada<br />

sobre todas las raíces positivas de la ecuación<br />

(70).<br />

En el caso particular de condiciones de contorno<br />

constantes<br />

f() t A , g()<br />

t B<br />

la ecuación (73) se reduce a<br />

unt e k n<br />

( , ) 2<br />

t 2<br />

<br />

2<br />

<br />

n<br />

<br />

2B<br />

<br />

A<br />

<br />

<br />

aJ bY a J aY b<br />

0( n )<br />

0( n ) <br />

0( n )<br />

0( n<br />

) <br />

k t<br />

<br />

n<br />

( e<br />

2 1) w( n)<br />

(75)<br />

<br />

Si además, wr ( ) W(constante), se deduce<br />

que [14]<br />

W <br />

wn ( ) hJ bY a J aY b <br />

2 1( n )<br />

0( n ) <br />

0( n )<br />

1( n<br />

)<br />

<br />

n<br />

2h<br />

<br />

<br />

<br />

n J na Y nb J nb Y na<br />

<br />

1( )<br />

1( )<br />

1( )<br />

1( ) (76)<br />

a<br />

Sustituyendo (76) en (75), se tiene<br />

1 2<br />

kn<br />

t<br />

unt ( , ) <br />

2<br />

e<br />

n<br />

<br />

<br />

2 2B<br />

<br />

<br />

A<br />

<br />

<br />

aJ <br />

<br />

<br />

0( nbY )<br />

0( na) J0( naY )<br />

0( nb)<br />

<br />

2<br />

k ( e n t <br />

1) WhJ1( nb) Y0( na) J0( na) Y1( nb)<br />

<br />

<br />

2h<br />

<br />

<br />

<br />

<br />

n J1( na) Y1( nb) J1( nb) Y1( na) a<br />

(77)<br />

<br />

Para A<br />

B<br />

0, la ecuación (77) se reduce a<br />

W 2<br />

unt e<br />

kn<br />

t<br />

( , ) <br />

2<br />

n<br />

<br />

hJ bY a J aY b<br />

1( n )<br />

0( n ) <br />

0( n )<br />

1( n<br />

)<br />

<br />

2h<br />

<br />

<br />

<br />

n J na Y nb J nb Y na<br />

<br />

1( )<br />

1( )<br />

1( )<br />

1( ) (78)<br />

a<br />

n<br />

n<br />

n<br />

y la solución (74) queda<br />

1<br />

ur ( , t) h unt ( , ); n<br />

r<br />

5,<br />

0<br />

<br />

<br />

<br />

nJ0<br />

( nb) u( n, t) hY0( na) nY1( na) J0(<br />

nr) hJ0( na) nJ1( na) Y0( nr)<br />

ur ( , t)<br />

<br />

2<br />

2<br />

hJ ( a J a <br />

n 1 0<br />

n ) n 0( n<br />

) J<br />

2 2 2<br />

0( n<br />

b )( n<br />

h )<br />

2 2<br />

y, finalmente, resulta<br />

2<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

nJ0<br />

( nb) u( n, t) hY0( na) nY1( na) J0(<br />

nr) hJ0( na) nJ1( na) Y0( nr)<br />

ur ( , t)<br />

<br />

2<br />

2<br />

2 2 2<br />

n1 hJ ( a) J ( a) (<br />

n h ) J ( n b)<br />

2 2<br />

2<br />

0 n n 1 n<br />

<br />

<br />

<br />

0<br />

<br />

<br />

(74)<br />

Rev. Téc. Ing. Univ. Zulia. Vol. 31, Edición Especial, 2008


78 Pirela y Villalobos<br />

<br />

<br />

<br />

nJ0<br />

( nb) u( n, t) hY0( na) nY1( na) J0(<br />

nr) hJ0( na) nJ1( na) Y0( nr)<br />

ur ( , t)<br />

<br />

2<br />

2<br />

2 2 2<br />

n1 hJ ( a) J ( a) (<br />

n h ) J ( n b)<br />

2 2<br />

2<br />

0 n n 1 n<br />

<br />

<br />

<br />

0<br />

<br />

<br />

(79)<br />

Referencias Bibliográficas<br />

1. Al-Hajri, M. and Kalla, S. L.: On an integral<br />

trans<strong>form</strong> involving Bessel functions. Bull.<br />

Soc. Math. Macedonie 28 (2004), 5-18.<br />

2. Bermúdez, M.: Sobre una Trans<strong>form</strong>ada Finita<br />

<strong>General</strong>izada de <strong>Hankel</strong>; Rev. Téc. Ing.<br />

Univ. Zulia, Vol. 17, Nº 1, (1994).<br />

3. Bermúdez, M.: Una nota sobre unas Trans<strong>form</strong>adas<br />

Finitas <strong>General</strong>izadas de <strong>Hankel</strong><br />

de segunda y tercera clase; Rev. Téc. Ing.<br />

Univ. Zulia, Vol. 18, Nº 2, (1995).<br />

4. Churchill, R. V.: Series de Fourier y Problemas<br />

de Contorno, Ediciones del Castillo, S.<br />

A., Madrid, (1966).<br />

5. Eldabe N. T., El-Shahed M. and Shawkey M.:<br />

An extension <strong>of</strong> <strong>the</strong> <strong>finite</strong> <strong>Hankel</strong> trans<strong>form</strong>.<br />

Appl. Math. Comput. 151 (2004), 713-717.<br />

6. Garg Mridula; Rao Alka and Kalla S. L.: On a<br />

<strong>general</strong>ized <strong>finite</strong> <strong>Hankel</strong> trans<strong>form</strong>. Appl.<br />

Math. Comput. (2007). Doi. 10.1016/j.amc.<br />

2007..01.076.<br />

7. Kalla, S. L. and Kushwaha, R. S.: Production<br />

<strong>of</strong> Heat in an In<strong>finite</strong> Cylinder, I y II, Acta<br />

Mexicana Cient. Téc., 4, (1970), 89-93.<br />

8. Kalla, S. L.: Conduction <strong>of</strong> Heat in a Finite<br />

Circular Cylinder, Rev. Math. Univ. Lourenco<br />

Marques, 3, (1972), 11-18.<br />

9. Kalla S. L., Battig A. and Luccioni R.: Conduction<br />

<strong>of</strong> Heat in a Semi-in<strong>finite</strong> Hollow Cylinder,<br />

Anal. Sti. Univ. Al. I. Cuza (Física), 19,<br />

(1973), 9-16.<br />

10. Kalla S. L., Battig A. and Luccioni R.: Motion<br />

<strong>of</strong> Viscous Fluid Contained between Two Coaxial<br />

Cylinders, Ranchi Math. Jour., (1973).<br />

11. Khajah, H.G.: A modified <strong>Hankel</strong> trans<strong>form</strong>.<br />

Integral Trans<strong>form</strong>s and Special Functions,<br />

14 (2003), 403-412.<br />

12. Lebedev, N. N.: Special Functions and Their<br />

Applications, Dover Publications, Inc., New<br />

York, (1972).<br />

13. Manstretta, R.: Algunas Aplicaciones de las<br />

Trans<strong>form</strong>adas Finitas <strong>General</strong>izadas de<br />

<strong>Hankel</strong>; Tesis de Maestría en Matemática<br />

Aplicada, Universidad del Zulia, (1999).<br />

14. Pirela Cristalino, Ovidio José. <strong>Forma</strong> <strong>General</strong><br />

de la Trans<strong>form</strong>ada Finita de <strong>Hankel</strong> y<br />

sus Aplicaciones a Problemas de Contorno.<br />

(2006) Trabajo de Grado. Universidad del<br />

Zulia. Facultad de Ingeniería. División de<br />

Postgrado. Maracaibo, Tutor: Pr<strong>of</strong>. Alfredo<br />

Villalobos Mena.<br />

15. Prudnikov A. P., Brichkov Yu. and Marichev<br />

O. I.: Integrals and Series, Vol. 2, Gordon and<br />

Breach Science Publishers, New York,<br />

(1986).<br />

16. Ramírez, I.: Trans<strong>form</strong>adas Finitas <strong>General</strong>izadas<br />

de <strong>Hankel</strong> de Cuarta, Quinta y Sexta<br />

clase; Tesis de Maestría en Matemática Aplicada,<br />

Universidad del Zulia, (1996).<br />

17. Sneddon I. N.: The Use <strong>of</strong> Integral Trans<strong>form</strong>s,<br />

Tata Mc Graw-Hill Publishing Company<br />

Ltd., New Delhi, (1974).<br />

18. Villalobos, A.: Temperaturas Estacionarias<br />

en un Cilindro Hueco con Conductividad<br />

Térmica Variable usando una Trans<strong>form</strong>ada<br />

de <strong>Hankel</strong>; Rev. Téc. Ing. Univ. Zulia, Vol. 25,<br />

Nº 3, (2002).<br />

19. Watson G. N.: A Treatise on <strong>the</strong> Theory <strong>of</strong><br />

Bessel Functions, Cambridge University<br />

Press, London, (1980).<br />

20. Zill D. G.: Ecuaciones Diferenciales con Aplicaciones<br />

de Modelado, International Thomson<br />

Editores, México, (1997).<br />

Recibido el 30 de Junio de 2007<br />

En <strong>form</strong>a revisada el 31 de Julio de 2008<br />

Rev. Téc. Ing. Univ. Zulia. Vol. 31, Edición Especial, 2008

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!