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