hankel
所属分类:matlab编程
开发工具:matlab
文件大小:114KB
下载次数:120
上传日期:2008-02-09 05:12:03
上 传 者:
DanielWood
说明: Fast Hankel Transform Algorithm
(Fast Hankel Transform Algorithm)
文件列表:
hankel\c1.eps (14863, 2001-07-09)
hankel\c1.m (35, 2001-07-05)
hankel\c2.eps (14957, 2001-07-09)
hankel\c2.m (37, 2001-07-05)
hankel\c3.eps (13944, 2001-07-09)
hankel\c3.m (28, 2001-07-05)
hankel\c4.eps (14317, 2001-07-09)
hankel\c4.m (33, 2001-07-05)
hankel\c5.eps (13955, 2001-07-09)
hankel\c5.m (30, 2001-07-05)
hankel\hankel0.m (45022, 2006-10-11)
hankel\hankel01.m (68407, 2006-10-04)
hankel\hankel1.m (44911, 2006-10-04)
hankel\hankelpts.m (21422, 2001-07-09)
hankel\hankelwts.m (47038, 2001-07-09)
hankel\test.m (2347, 2001-07-05)
hankel (0, 2006-10-11)
The MATLAB routines in this directory implement Hankel transforms of order
0 and 1 by digital filtering. These routines are based on a Fortran program
by Walt Anderson which was published as:
Anderson, W. L., 1979, Computer Program Numerical Integration of Related
Hankel Transforms of Orders 0 and 1 by Adaptive Digital Filtering.
Geophysic, 44(7):1287-1305.
The actual weights used in this code are from a later updated version
of the code that was originally published in 1979.
This digital filtering approach is appropriate for use only with functions
that are smooth. If there are discontinuities in the function or its
derivatives then the code is likely to produce inaccurate results.
For purposes of this code, the Hankel transforms are defined as
H0(B)=int(g(r)*J0(B*r),r=0..infinity)
and
H1(B)=int(g(r)*J1(B*r),r=0..infinity)
The routine hankel0 computes a 0th order transform. The routine hankel1
computes a 1st order transform. The routine hankel01 computes both
transforms. The test script does transforms on five functions with
known transforms and compares the results with known exact values.
The functions c1 through c5 are the five functions that are transformed
by test.
In some cases it can be helpful to compute several related functions
together and then apply the weights to do the transform. The
functions hankelpts and hankelwts can be used to do this.
近期下载者:
相关文件:
收藏者: