5.0

5.0 | 3 ratings Rate this file 21 Downloads (last 30 days) File Size: 136 KB File ID: #13371

Hankel transform

by Marcel Leutenegger

 

13 Dec 2006 (Updated 11 Apr 2007)

Efficient implementations of the Hankel transform and the inverse Hankel transform, respectively.

| Watch this File

File Information
Description

The Hankel transform of order n transforms rotationally symmetric inputs in a computationally efficient manner. In particular, the Hankel transform of order 0 is equivalent to the two-dimensional Fourier transform of a rotationally symmetric input. This package contains four implementations of the Hankel transform and the inverse Hankel transform, respectively.

"hat" and "ihat" perform the Hankel transform of order n with a direct integration using a matrix product. "ht" and "iht" perform the Hankel transform of order 0 by integrating the Bessel kernel a priori. "dht" and "idht" implement the quasi-discrete Hankel transform of integer order n. And, last but not least, "fht" and "ifht" implement the quasi fast Hankel transform of order n.

For more implementation details, please refer to the online documentation at

http://ioalinux1.epfl.ch/~mleutene/MATLABToolbox/HankelTransform.html

MATLAB release MATLAB 6.0 (R12)
Tags for This File  
Everyone's Tags
hankel transform(4), integration, mathematics, quasi fast, quasidiscrete
Tags I've Applied
Add New Tags Please login to tag files.
Please login to add a comment or rating.
Comments and Ratings (5)
12 Apr 2013 Toader

I studied hat.m and ihat.m functions. They work fine if the order it's very small, but the more I increase the order it's not working properly. Why?
Thank you.

14 Jun 2012 Michael

Seems to work nicely, but what is the scaling for (the absolute value of) the output for ht.m?

25 Jan 2012 Nicolas Grisouard  
19 Mar 2009 John Ferguson

Update for documentation:
http://wwwuser.gwdg.de/~mleuten/MATLABToolbox/HankelTransform.html

12 Mar 2007 Eunmi Choi

Thank you!

Updates
11 Apr 2007

Fixed an error in "fht.m" reported by Mark W. Sprague - line 49 should read "if N > 1".

Contact us