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.

近期下载者

相关文件


收藏者