Quadrature approximation properties of the spiral-phase quadrature transform

The notion of the 1-D analytic signal is well understood and has found many applications. At the heart of the analytic signal concept is the Hilbert transform. The problem in extending the concept of analytic signal to higher dimensions is that there is no unique multidimensional definition of the Hilbert transform. Also, the notion of analyticity is not so well under stood in higher dimensions. Of the several 2-D extensions of the Hilbert transform, the spiral-phase quadrature transform or the Riesz transform seems to be the natural extension and has attracted a lot of attention mainly due to its isotropic properties. From the Riesz transform, Larkin et al. constructed a vortex operator, which approximates the quadratures based on asymptotic stationary-phase analysis. In this paper, we show an alternative proof for the quadrature approximation property by invoking the quasi-eigenfunction property of linear, shift-invariant systems. We show that the vortex operator comes up as a natural consequence of applying this property. We also characterize the quadrature approximation error in terms of its energy as well as the peak spatial-domain error. Such results are available for 1-D signals, but their counter part for 2-D signals have not been provided. We also provide simulation results to supplement the analytical calculations.

Generalized eigenvalue decomposition of the field autocorrelation in correlation diffusion of photons in turbid media

We propose a novel numerical method based on a generalized eigenvalue decomposition for solving the diffusion equation governing the correlation diffusion of photons in turbid media. Medical imaging modalities such as diffuse correlation tomography and ultrasound-modulated optical tomography have the (elliptic) diffusion equation parameterized by a time variable as the forward model. Hitherto, for the computation of the correlation function, the diffusion equation is solved repeatedly over the time parameter. We show that the use of a certain time-independent generalized eigenfunction basis results in the decoupling of the spatial and time dependence of the correlation function, thus allowing greater computational efficiency in arriving at the forward solution. Besides presenting the mathematical analysis of the generalized eigenvalue problem on the basis of spectral theory, we put forth the numerical results that compare the proposed numerical method with the standard technique for solving the diffusion equation.

Binaural Signal Processing Motivated Generalized Analytic Signal Construction and AM-FM Demodulation

Binaural hearing studies show that the auditory system uses the phase-difference information in the auditory stimuli for localization of a sound source. Motivated by this finding, we present a method for demodulation of amplitude-modulated-frequency-modulated (AM-FM) signals using a ignal and its arbitrary phase-shifted version. The demodulation is achieved using two allpass filters, whose impulse responses are related through the fractional Hilbert transform (FrHT). The allpass filters are obtained by cosine-modulation of a zero-phase flat-top prototype halfband lowpass filter. The outputs of the filters are combined to construct an analytic signal (AS) from which the AM and FM are estimated. We show that, under certain assumptions on the signal and the filter structures, the AM and FM can be obtained exactly. The AM-FM calculations are based on the quasi-eigenfunction approximation. We then extend the concept to the demodulation of multicomponent signals using uniform and non-uniform cosine-modulated filterbank (FB) structures consisting of flat bandpass filters, including the uniform cosine-modulated, equivalent rectangular bandwidth (ERB), and constant-Q filterbanks. We validate the theoretical calculations by considering application on synthesized AM-FM signals and compare the performance in presence of noise with three other multiband demodulation techniques, namely, the Teager-energy-based approach, the Gabor's AS approach, and the linear transduction filter approach. We also show demodulation results for real signals.

Boundary perturbations and the Helmholtz equation in three dimensions

We propose an analytic perturbative scheme in the spirit of Lord Rayleigh's work for determining the eigenvalues of the Helmholtz equation in three dimensions inside an arbitrary boundary where the eigenfunction satisfies either the Dirichlet boundary condition or the Neumann boundary condition. Although numerous works are available in the literature for arbitrary boundaries in two dimensions, to the best of our knowledge the formulation in three dimensions is proposed for the first time. In this novel prescription, we have expanded the arbitrary boundary in terms of spherical harmonics about an equivalent sphere and obtained perturbative closed-form solutions at each order for the problem in terms of corrections to the equivalent spherical boundary for both the boundary conditions. This formulation is in parallel with the standard time-independent Rayleigh-Schrodinger perturbation theory. The efficacy of the method is tested by comparing the perturbative values against the numerically calculated eigenvalues for spheroidal, superegg and superquadric shaped boundaries. It is shown that this perturbation works quite well even for wide departure from spherical shape and for higher excited states too. We believe this formulation would find applications in the field of quantum dots and acoustical cavities.

Unusual eigenvalue spectrum and relaxation in the Levy-Ornstein-Uhlenbeck process

We consider the rates of relaxation of a particle in a harmonic well, subject to Levy noise characterized by its Levy index mu. Using the propagator for this Levy-Ornstein-Uhlenbeck process (LOUP), we show that the eigenvalue spectrum of the associated Fokker-Planck operator has the form (n + m mu)nu where nu is the force constant characterizing the well, and n, m is an element of N. If mu is irrational, the eigenvalues are all nondegenerate, but rational mu can lead to degeneracy. The maximum degeneracy is shown to be 2. The left eigenfunctions of the fractional Fokker-Planck operator are very simple while the right eigenfunctions may be obtained from the lowest eigenfunction by a combination of two different step-up operators. Further, we find that the acceptable eigenfunctions should have the asymptotic behavior vertical bar x vertical bar(-n1-n2 mu) as vertical bar x vertical bar -> infinity, with n(1) and n(2) being positive integers, though this condition alone is not enough to identify them uniquely. We also assert that the rates of relaxation of LOUP are determined by the eigenvalues of the associated fractional Fokker-Planck operator and do not depend on the initial state if the moments of the initial distribution are all finite. If the initial distribution has fat tails, for which the higher moments diverge, one can have nonspectral relaxation, as pointed out by Toenjes et al. Phys. Rev. Lett. 110, 150602 (2013)].

L-infinity norms of holomorphic modular forms in the case of compact quotient

We prove a sub-convex estimate for the sup-norm of L-2-normalized holomorphic modular forms of weight k on the upper half plane, with respect to the unit group of a quaternion division algebra over Q. More precisely we show that when the L-2 norm of an eigenfunction f is one, parallel to f parallel to(infinity) <<(epsilon) k(1/2-1/33+epsilon) for any epsilon > 0 and for all k sufficiently large.

Stress intensity factors calculation in anti-plane fracture problem by orthogonal integral extraction method based on femol

For an anti-plane problem, the differential operator is self-adjoint and the corresponding eigenfunctions belong to the Hilbert space. The orthogonal property between eigenfunctions (or between the derivatives of eigenfunctions) of anti-plane problem is exploited. We developed for the first time two sets of radius-independent orthogonal integrals for extraction of stress intensity factors (SIFs), so any order SIF can be extracted based on a certain known solution of displacement (an analytic result or a numerical result). Many numerical examples based on the finite element method of lines (FEMOL) show that the present method is very powerful and efficient.

Fast numerical methods for mixed, singular Helmholtz boundary value problems and Laplace eigenvalue problems - with applications to antenna design, sloshing, electromagnetic scattering and spectral geometry

This thesis presents a novel class of algorithms for the solution of scattering and eigenvalue problems on general two-dimensional domains under a variety of boundary conditions, including non-smooth domains and certain "Zaremba" boundary conditions - for which Dirichlet and Neumann conditions are specified on various portions of the domain boundary. The theoretical basis of the methods for the Zaremba problems on smooth domains concern detailed information, which is put forth for the first time in this thesis, about the singularity structure of solutions of the Laplace operator under boundary conditions of Zaremba type. The new methods, which are based on use of Green functions and integral equations, incorporate a number of algorithmic innovations, including a fast and robust eigenvalue-search algorithm, use of the Fourier Continuation method for regularization of all smooth-domain Zaremba singularities, and newly derived quadrature rules which give rise to high-order convergence even around singular points for the Zaremba problem. The resulting algorithms enjoy high-order convergence, and they can tackle a variety of elliptic problems under general boundary conditions, including, for example, eigenvalue problems, scattering problems, and, in particular, eigenfunction expansion for time-domain problems in non-separable physical domains with mixed boundary conditions.

The propagation of acoustic waves in a slowly varying duct with radially sheared axial mean flow

We consider the propagation of acoustic waves along a cylindrical duct carrying radially sheared axial mean flow, in which the duct radius is allowed to vary slowly along the axis. In previous work [A.J. Cooper & N. Peake, Journal of Fluid Mechanics 445 (2001) 207-234.] radially sheared axial mean flow with nonzero swirl in a slowly varying duct was considered, but in this paper we set the swirl to zero, thereby allowing simplification of the calculations of both the mean and unsteady flows. In this approach the acoustic wavenumber and corresponding eigenfunction are determined locally, while the wave amplitude is found by solving an evolution equation along the duct. Sample results are presented, including a case in which, perhaps surprisingly, the number of cut-on modes increases as the duct radius decreases. © 2013 Elsevier Ltd. All rights reserved.

On the radiation and diffraction of linear water waves by an infinitely long rectangular structure submerged in oblique seas

The radiation and diffraction of linear water waves by an infinitely long rectangular structure submerged in oblique seas of finite depth is investigated. The analytical expressions for the radiated and diffracted potentials are derived as infinite series by use of the method of separation of variables. The unknown coefficients in the series are determined by the eigenfunction expansion matching method. The expressions for wave forces, hydrodynamic coefficients and reflection and transmission coefficients are given and verified by the boundary element method. Using the present analytical solution, the hydrodynamic influences of the angle of incidence, the submergence, the width and the thickness of the structure on the wave forces, hydrodynamic coefficients, and reflection and transmission coefficients are discussed in detail.

Reflection and transmission properties of double submerged rectangular blocks

An analytic method is used to study the reflection and transmission coefficients of the double submerged rectangular blocks (DSRBs) in oblique waves.. The scattering potentials are obtained by means of the eigenfunction expansion method, and expressions for the reflection and transmission coefficients are determined. The boundary element method is employed to verify the correctness of the present analytical method. The DSRBs have better performance than the single submerged rectangular block (SSRB) in certain cases. The reflection and transmission properties of the DSRBs are investigated for some specific cases, and the influences of the geometric parameters are also presented.

Scattering of linear water waves by a fixed and infinitely long rectangular structure parallel to a vertical wall in oblique seas

The scattering of linear water waves by an infinitely long rectangular structure parallel to a vertical wall in oblique seas is investigated. Analytical expressions for the diffracted potentials are derived using the method of separation of variables. The unknown coefficients in the expressions are determined through the application of the eigenfunction expansion matching method. The expressions for wave forces on the structure are given. The calculated results are compared with those obtained by the boundary element method. In addition, the influences of the wall, the angle of wave incidence, the width of the structure, and the distance between the structure and the wall on wave forces are discussed. The method presented here can be easily extended to the study of the diffraction of obliquely incident waves by multiple rectangular structures.

Optical analysis of AlGaInP laser diodes with real refractive index guided self-aligned structure

Optical modes of AlGaInP laser diodes with real refractive index guided self-aligned (RISA) structure were analyzed theoretically on the basis of two-dimension semivectorial finite-difference methods (SV-FDMs) and the computed simulation results were presented. The eigenvalue and eigenfunction of this two-dimension waveguide were obtained and the dependence of the confinement factor and beam divergence angles in the direction of parallel and perpendicular to the pn junction on the structure parameters such as the number of quantum wells, the Al composition of the cladding layers, the ridge width, the waveguide thickness and the residual thickness of the upper P-cladding layer were investigated. The results can provide optimized structure parameters and help us design and fabricate high performance AlGaInP laser diodes with a low beam aspect ratio required for optical storage applications.

Estimation of Transverse Thermal Conductivity of Doubly-periodic Fiber Reinforced Composites

For steady-state heat conduction a new variational functional for a unit cell of composites with periodic microstructures is constructed by considering the quasi-periodicity of the temperature field and in the periodicity of the heat flux fields. Then by combining with the eigenfunction expansion of complex potential which satisfies the fiber-matrix interface conditions, an eigenfunction expansion-variational method (EEVM) based on a unit cell is developed. The effective transverse thermal conductivities of doubly-periodic fiber reinforced composites are calculated, and the first-order approximation formula for the square and hexagonal arrays is presented,which is convenient for engineering application. The numerical results show a good convergency of the presented method, even through the fiber volume fraction is relatively high. Comparisons with the existing analytical and experimental results are made to demonstrate the accuracy and validity of the first-order approximation formula for the hexagonal array.

指数矩阵近似计算在地震波场延拓中的应用

Seismic wave field numerical modeling and seismic migration imaging based on wave equation have become useful and absolutely necessarily tools for imaging of complex geological objects. An important task for numerical modeling is to deal with the matrix exponential approximation in wave field extrapolation. For small value size matrix exponential, we can approximate the square root operator in exponential using different splitting algorithms. Splitting algorithms are usually used on the order or the dimension of one-way wave equation to reduce the complexity of the question. In this paper, we achieve approximate equation of 2-D Helmholtz operator inversion using multi-way splitting operation. Analysis on Gauss integral and coefficient of optimized partial fraction show that dispersion may accumulate by splitting algorithms for steep dipping imaging. High-order symplectic Pade approximation may deal with this problem, However, approximation of square root operator in exponential using splitting algorithm cannot solve dispersion problem during one-way wave field migration imaging. We try to implement exact approximation through eigenfunction expansion in matrix. Fast Fourier Transformation (FFT) method is selected because of its lowest computation. An 8-order Laplace matrix splitting is performed to achieve a assemblage of small matrixes using FFT method. Along with the introduction of Lie group and symplectic method into seismic wave-field extrapolation, accurate approximation of matrix exponential based on Lie group and symplectic method becomes the hot research field. To solve matrix exponential approximation problem, the Second-kind Coordinates (SKC) method and Generalized Polar Decompositions (GPD) method of Lie group are of choice. SKC method utilizes generalized Strang-splitting algorithm. While GPD method utilizes polar-type splitting and symmetric polar-type splitting algorithm. Comparing to Pade approximation, these two methods are less in computation, but they can both assure the Lie group structure. We think SKC and GPD methods are prospective and attractive in research and practice.