Fast Hankel Transform And Its Application To Beam

Wei-ci Liu,Fa-qiang Wang,Rui-sheng Liang

DOI: https://doi.org/10.12783/dtcse/cmsam2018/26573

2018-01-01

Abstract:In this paper, the numerical technique that is Fast Hankel Transform is presented for the beam propagation in two transverse dimensions with cylindrical symmetry. The Fast Hankel Transform is based on Fastest Fourier Transform in the West, and its application is discussed.

What problem does this paper attempt to address?

Hexagonal Array Fast Fourier Transform Beamforming Algorithm

PAN Xiang,ZHOU Bin,LI Jian-long,GONG Xian-yi

DOI: https://doi.org/10.3969/j.issn.1004-1699.2007.05.031

2007-01-01

Abstract:While the amount of hexagonal array elements is high,the conventional phase shift-sum beamforming approach requires a heavy computational load.To solve this problem,a beamforming algorithm utilizing hexagonal fast Fourier transform(HFFT)is presented.Using hexagonal discrete Fourier transform(HDFT)to accomplish beamforming will significantly reduce the computational complexity because of the existence of HFFT algorithm to compute HDFT.We firstly use fast Fourier transform(FFT)on each array element to transform the input into frequency domain,then,rearrange the FFT results and obtain conventional frequency domain beamforming output using HFFT on all array elements for each frequency bin.Theoretical analysis suggests the presented algorithm reduce the computational load by an extent of 95% compared with the conventional approach for narrow-band signal beamforming.The results of computer simulation and experimental data analysis illustrate that the algorithm can dramatically reduce the computational complexity and is applicable for engineering implementation while maintaining the performance of the array processing.
A FOURIER TRANSFORM SOLUTION FOR AN INFINITE BEAM ON A SEMI-INFINITE SPACE LOADED BY A CONCENTRATED FORCE

Tong Genshu,Xuan Zejun

DOI: https://doi.org/10.3969/j.issn.1007-9963.2013.08.005

2013-01-01

Steel Construction

Abstract:This paper made a study on the classic problem of an infinite beam on a semi-infinite space acted by a concentrated load on the beam.One of the basic equation to determine the settlement of the semi-space below the beam was based on the well-known Boussinesq's solution of a normal load on a semi-infinite space,the other was the beam differential equation.Fourier transform technique was used to solve the equations.Because the Fourier transform was involved with an integration of function with singularity,it was found that a constant could not be determined analytically.So finite element analysis was used.Different values of the constant were tried to get a best fit to the FE results,and the undetermined coefficient was found to be C=γ+e≈3.295 5.The characteristic length of the problem was found to be a fourth root expression rather than a third root one in the current available literature.Finally an equivalent bearing length was proposed to simplify the computation of local maximum bearing compressive stess of the beam under the load.The formulae of bending moment and shear force of the beam could be used strength calculation of track beam.
Two-dimensional Solutions for Orthotropic Materials by the State Space Method

Xu Rongqiao,Ding Haojiang

DOI: https://doi.org/10.1016/j.compstruct.2005.10.006

IF: 6.3

2005-01-01

Composite Structures

Abstract:A state space method for two-dimensional problems of orthotropic materials is proposed which provides an alternative method to investigate the mechanical behavior of homogeneous and laminated beams. The method is illustrated for a homogenous/laminated beam subjected to time-dependent transverse loads. A particular solution for a straight-crested wave propagating along an infinite beam is obtained. The bending, buckling and free vibrations of composite beams are also studied. Numerical results are presented and compared with those available in the literature.
Method Based on Fast Fourier Transform for Calculating Conical Refraction of Beams with Noncircular Symmetry

Pengcheng Nie,Dagong Jia,Chuang Du,Hongxia Zhang,Tianhua Xu,Tiegeng Liu

DOI: https://doi.org/10.1109/jphot.2017.2669520

IF: 2.4

2017-01-01

IEEE Photonics Journal

Abstract:Conical refraction of optical beams with circular symmetry is often analyzed using Belsky-Khapalyuk-Berry (BKB) theory; however, for beams with noncircular symmetry, it is difficult to obtain analytical expressions for far-field diffraction patterns. We propose a method, based on fast Fourier transform (FFT), for calculating conical refraction of beams with noncircular symmetry and verify it experimentally using a quasi-plane wave passing through a square aperture and focusing lens. Excellent agreement between theoretical and experimental results has been achieved.
Analytical solution for fixed-end orthotropic beams subjected to uniform load

De-jin HUANG,Hao-jiang DING,Hui-ming WANG

DOI: https://doi.org/10.3785/j.issn.1008-973X.2006.03.031

2006-01-01

Abstract:For the plane stress problem of fixed end orthotropic beams subjected to uniform load, the analytical solutions of stress and displacement were resolved. A stress function involving several unknown coefficients was constructed in the solution process, and the stress and displacement expressions were obtained through traditional Airy stress function methodology. After the fixed end boundary condition was treated with two methods, the unknown coefficients in stress function were determined with stress and displacement boundary conditions, and the analytical solutions were achieved. The results show that compared with finite element method (FEM) numerical solutions by Nastran codes, the analytical solutions are close to and agree well with the FEM solutions. The analytical solutions can reach high precision for the beams with large span height ratio, and can be degenerated into isotropic results.
Fast computation of radar cross-section by fast multipole method in conjunction with lifting wavelet-like transform

Mingsheng Chen,Wei E. I. Sha,Q. Wu,Zhixiang Huang,Xianliang Wu

DOI: https://doi.org/10.1049/iet-map.2009.0417

2010-01-01

Abstract:The fast multipole method (FMM) in conjunction with the lifting wavelet-like transform scheme is proposed for the scattering analysis of differently shaped three-dimensional perfectly electrical conducting objects. As a flexible and efficient matrix compression technique, the proposed method can sparsify the aggregation matrix and disaggregation matrix in real time with compression ratio about 30%...
Fast Hankel tensor–vector product and its application to exponential data fitting

Weiyang Ding,Liqun Qi,Yimin Wei

DOI: https://doi.org/10.1002/nla.1970

2015-01-01

Numerical Linear Algebra with Applications

Abstract:This paper is contributed to a fast algorithm for Hankel tensor-vector products. First, we explain the necessity of fast algorithms for Hankel and block Hankel tensor-vector products by sketching the algorithm for both one-dimensional and multi-dimensional exponential data fitting. For proposing the fast algorithm, we define and investigate a special class of Hankel tensors that can be diagonalized by the Fourier matrices, which is called anti-circulant tensors. Then, we obtain a fast algorithm for Hankel tensor-vector products by embedding a Hankel tensor into a larger anti-circulant tensor. The computational complexity is about O(m2nlogmn) for a square Hankel tensor of order m and dimension n, and the numerical examples also show the efficiency of this scheme. Moreover, the block version for multi-level block Hankel tensors is discussed. Copyright (c) 2015 John Wiley & Sons, Ltd.
Fast solution of Foldy-Lax equations for two-dimensional radiation and scattering problems

Yao-Jiang Zhang,Er-Ping Li

DOI: https://doi.org/10.1109/EMCZUR.2006.214906

2006-01-01

Abstract:The general fast multipole expressions of arbitrary order Hankel functions are derived by using lowering and raising operators of cylindrical harmonics. These expressions are then used to transform the dense matrix in Foldy-Lax equations into a combination of sparse matrices (aggregation, translation and disaggregation matrices). Thus, the computational complexity of such an algorithm is found to be of O(N 1.5) instead of O(N2) of the traditional scattering matrix method, where N denotes the total harmonics number used to expand scattered of all the cylinders. The details of the implementation issues are investigated and, the accuracy and efficiency of this novel fast algorithm are verified by several numerical examples
A Nonuniform Fast Hankel Transform

Paul G. Beckman,Michael O'Neil

2024-11-15

Abstract:We describe a fast algorithm for computing discrete Hankel transforms of moderate orders from $n$ nonuniform points to $m$ nonuniform frequencies in $O((m+n)\log\min(n,m))$ operations. Our approach combines local and asymptotic Bessel function expansions with nonuniform fast Fourier transforms. The order of each expansion is adjusted automatically according to error analysis to obtain any desired precision $\varepsilon$. Several numerical examples are provided which demonstrate the speed and accuracy of the algorithm in multiple regimes and applications.

Numerical Analysis
Fast multiscale Gaussian beam methods for wave equations in bounded convex domains

Gang Bao,Jun Lai,Jianliang Qian

DOI: https://doi.org/10.1016/j.jcp.2013.12.034

IF: 4.645

2014-01-01

Journal of Computational Physics

Abstract:Motivated by fast multiscale Gaussian wavepacket transforms and multiscale Gaussian beam methods which were originally designed for pure initial-value problems of wave equations, we develop fast multiscale Gaussian beam methods for initial boundary value problems of wave equations in bounded convex domains in the high frequency regime. To compute the wave propagation in bounded convex domains, we have to take into account reflecting multiscale Gaussian beams, which are accomplished by enforcing reflecting boundary conditions during beam propagation and carrying out suitable reflecting beam summation. To propagate multiscale beams efficiently, we prove that the ratio of the squared magnitude of beam amplitude and the beam width is roughly conserved, and accordingly we propose an effective indicator to identify significant beams. We also prove that the resulting multiscale Gaussian beam methods converge asymptotically. Numerical examples demonstrate the accuracy and efficiency of the method.
Fast Multipole Accelerated Scattering Matrix Method for Multiple Scattering of A Large Number of Cylinders

Y. J. Zhang,E. P. Li

DOI: https://doi.org/10.2528/pier07030503

2007-01-01

Abstract:The lowering and raising operators of cylindrical harmonics are used to derive the general fast multipole expressions of arbitrary order Hankel functions. These expressions are then employed to transform the dense matrix in the scattering matrix method (SMM) into a combination of sparse matrices (aggregation, translation and disaggregation matrices). The novel method is referred to as fast multipole accelerated scattering matrix method (FMA-SMM). Theoretical study shows FMA-SMM has lower complexity O(N-1.5) instead of SMM's O(N-2), where N stands for total harmonics number used. An empirical formula is derived to relate the minimum group size in FMA-SMM to the highest order Hankel functions involved. The various implementation parameters are carefully investigated to guarantee the algorithm's accuracy and efficiency. The impact of the cylinders density on convergence rate of iterative solvers (BiCGStab(2) here), memory cost as well as CPU time is also investigated. Up to thousands of cylinders can be easily simulated and potential applications in photonic crystal devices are illustrated.
Fast Tensor-Vector Products

Weiyang Ding,Yimin Wei

DOI: https://doi.org/10.1016/B978-0-12-803953-3.50003-4

2016-01-01

Abstract:Hankel structures are widely employed in data analysis and signal processing. Not only Hankel matrices but also higher-order Hankel tensors arise frequently in disciplines such as exponential data fitting [13, 39, 93, 94, 89], frequency domain subspace identification [111], multidimensional seismic trace interpolation [118], and so on. As far as we know, the term "Hankel tensor" was first introduced by Luque and Thibon [80]. Boyer et al. [5, 11] discussed the higher-order singular value decompositions (HOSVD) of structured tensors, including symmetric tensors, Hankel tensors, and Toeplitz tensors in more detail. Moreover, Papy et al. employed Hankel-type tensors in exponential data fitting [13, 93, 94]. De Lathauwer [34] also concerned the "separation" of signals that can be modeled as sums of exponentials (or more generally, as exponential polynomials) by Hankel tensor approaches. As for the properties of Hankel tensors, Qi [102] recently investigated the spectral properties of Hankel tensor largely via the generating function. Song and Qi [112] investigated the spectral properties of Hilbert tensors, which are special Hankel tensors.
Computation of beam propagation in turbulent field

Yao Yang,Zhaofeng Cen,Xiaotong Li

DOI: https://doi.org/10.1117/12.2246998

2016-01-01

Abstract:The split-step Fourier method (SSFM) is introduced to analyze the beam propagation in a relatively large-sized turbulent filed, whose refractive-index profile is already detected. The numerical method is achieved by fast Fourier transform (FFT). To obtain the optimal sampling number, we propose an adaptive spread-spectrum method as an optimization. The SSFM is widely used for solving the nonlinear Schrodinger equation ([1]). The advantage of the SSFM is apparently its simple formalism and suitability to our situation. The direct numerical solution of the Helmholtz equation, derived from this method, yields detailed information of the spatial and angular properties of the propagation beam. On the other hand, a set of approximations restrict its applicability, the requirements for the accurate application of the method are summarized and a set of formulas is generalized in this paper. The efficiency of the SSFM depends on the sampling number, the adaptive spread-spectrum method yields optimal sampling number to increase the computational efficiency. To testify the accuracy of our algorithm, we use graded-index medium as the turbulent filed, for the reason that the beam propagation in turbulent field with random refractive-index profile is ruleless and has no unified reference. The simulation result testifies our algorithm is tremendously accurate, capable of selecting the optimal N automatically and much more computationally efficient than the original algorithm.
Propagation of Laguerre-Gaussian and Elegant Laguerre-Gaussian Beams in Apertured Fractional Hankel Transform Systems.

ZR Mei,DM Zhao

DOI: https://doi.org/10.1364/josaa.21.002375

2004-01-01

Journal of the Optical Society of America A

Abstract:On the basis of the fact that a hard-edged-aperture function can be expanded into a finite sum of complex Gaussian functions, approximate analytical expressions for the output field distribution of a Laguerre-Gaussian beam and an elegant Laguerre-Gaussian beam passing through apertured fractional Hankel transform systems are derived. Some numerical simulation comparisons are done, by using the approximate analytical formulas and diffraction integral formulas, and it is shown that our method can significantly improve the numerical calculation efficiency.
A unified framework for asymptotic analysis and computation of finite Hankel transform

Haiyong Wang

DOI: https://doi.org/10.48550/arXiv.1801.06950

2018-04-12

Abstract:In this paper we present a unified framework for asymptotic analysis and computation of the finite Hankel transform. This framework enables us to derive asymptotic expansions of the transform, including the cases where the oscillator has zeros and stationary points. As a consequence, two efficient and affordable methods for computing the transform numerically are developed and a detailed analysis of their asymptotic error estimate is carried out. Numerical examples are provided to confirm our analysis.

Numerical Analysis
Optical implement of Hartley transform on elliptical Gaussian beams

Xinyue Du,Daomu Zhao

DOI: https://doi.org/10.1016/j.optcom.2009.11.064

IF: 2.4

2010-01-01

Optics Communications

Abstract:The optical Hartley transform (HT) is expressed in matrix form. The analytical formula of an off-axial elliptical Gaussian beam passing through the HT system is obtained. Numerical examples show that the information about the displacement magnitude and direction of the source can be represented in the form of fringes at the output plane of the HT system.
On Quantum Mechanical Hankel Transform and Its Applications

Hong-Yi Fan,Ji-Suo Wang

DOI: https://doi.org/10.1088/0253-6102/48/2/020

2007-01-01

Communications in Theoretical Physics

Abstract:Based on our preceding works of how to relate the mathematical Hankel transform to quantum mechanicalrepresentation transform and how to express the Bessel equation by an operator identity in some appropriate represen-tations we propose the concept of quantum mechanical Hankel transform with regard to quantum state vectors.Thenwe discuss its new applications.
Fractional Hankel transform and the diffraction of misaligned optical systems

Fan Ge,Daomu Zhao,Shaomin Wang

DOI: https://doi.org/10.1080/09500340410001674394

IF: 1.3

2006-01-01

Journal of Modern Optics

Abstract:Under the condition of rotational symmetry fractional Fourier transforms become fractional Hankel transforms. The relationships between the fractional Hankel transform and the diffraction of a misaligned optical system are established in this paper. By settling and discussing the diffraction integral formulae of a misaligned ABCD optical system, the fractional Hankel transform of an input function can be obtained by selecting a suitable reference plane under given conditions.
Fast Finite Hilbert Transform Via DE Quadrature Scheme

Changguo Ji,Yidong Cui,Wentian Cao,Shanglian Bao

DOI: https://doi.org/10.3233/xst-2010-0237

2010-01-01

Abstract:Purpose: The finite Hilbert transform (FHT) or inverse finite Hilbert transform (IFHT) is recently found to have some important applications in computerized tomography (CT) arena [1-6], where they are used to filter the derivatives of back-projected data in the chord-line based CT reconstruction algorithms. In this paper, we implemented, improved and validated a fast numerical solution to the FHT via a double exponential (DE) integration scheme. A same strategy can be used to compute IFHT.Methods: To overcome the Underflow of floating-point numbers, we first determined the range of variable transformation from the minimum positive value of single or double precision floating point number, the integration step can be further determined by the range of variable transformation and the integration level. Two functions with their known analytical FHTs are used to validate the implementation of the FHT via DE scheme. The surface map and 2D contour of the FHT transformation error with respect to integration level and the range of the variable transformation are used to numerically determine the optimal numbers for a fast FHT.Results: Given a specific precision, the lowest integration level and the optimal range of variable transformation, which are used to transform a signal with a certain degree Of fluctuation, can be numerically determined by the surface map and 2D contour of the standard deviation of transformation error. These two numbers can then be taken to efficiently compute the FHT for other signals with the same or less degree of fluctuation.Conclusions: The FHT via DE scheme and the numerical method to determine the integration level and the range of transformation can be used for fast FHT in certain applications, such as data filtering in chord-line based CT reconstruction algorithms.
New method to calculate Sommerfeld integral - fast Hankel transform

Jun Hu,Zaiping Nie

1998-01-01

Tien Tzu Hsueh Pao/Acta Electronica Sinica

Abstract:Fast Hankel transform is applied to calculate Sommerfeld integral. It attains the value of integral rapidly and accurately by the convolution of input function with filter coefficients under the finite length of the filter and avoids the tedious numerical integration and computation of Bessel function, so the computational efficiency and accuracy are improved greatly.

Fast Hankel Transform And Its Application To Beam

Hexagonal Array Fast Fourier Transform Beamforming Algorithm

A FOURIER TRANSFORM SOLUTION FOR AN INFINITE BEAM ON A SEMI-INFINITE SPACE LOADED BY A CONCENTRATED FORCE

Two-dimensional Solutions for Orthotropic Materials by the State Space Method

Method Based on Fast Fourier Transform for Calculating Conical Refraction of Beams with Noncircular Symmetry

Analytical solution for fixed-end orthotropic beams subjected to uniform load

Fast computation of radar cross-section by fast multipole method in conjunction with lifting wavelet-like transform

Fast Hankel tensor–vector product and its application to exponential data fitting

Fast solution of Foldy-Lax equations for two-dimensional radiation and scattering problems

A Nonuniform Fast Hankel Transform

Fast multiscale Gaussian beam methods for wave equations in bounded convex domains

Fast Multipole Accelerated Scattering Matrix Method for Multiple Scattering of A Large Number of Cylinders

Fast Tensor-Vector Products

Computation of beam propagation in turbulent field

Propagation of Laguerre-Gaussian and Elegant Laguerre-Gaussian Beams in Apertured Fractional Hankel Transform Systems.

A unified framework for asymptotic analysis and computation of finite Hankel transform

Optical implement of Hartley transform on elliptical Gaussian beams

On Quantum Mechanical Hankel Transform and Its Applications

Fractional Hankel transform and the diffraction of misaligned optical systems

Fast Finite Hilbert Transform Via DE Quadrature Scheme

New method to calculate Sommerfeld integral - fast Hankel transform