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

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

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!