Local approximation on surfaces with discontinuities,given limited order Fourier coefficients

We present a spline approximation method for a piece of a surface where jump discontinuities occur along curves. The data for the surface is assumed to be Fourier coefficients which are limited in order and possibly contaminated with noise. The support of the approximation is bounded by three sides of a rectangle with a fourth boundary possibly curved. Discontinuities of the surface may occur across the curved side and linear sides adjacent to it. The approximation uses a small number of lines through the support and parallel to the straight boundary lines that are adjacent to the curve. Along each line a one-dimensional spline approximation is done for a section of the surface over the line. This approximation uses two-dimensional Fourier coefficient data, localizing spline functions, and a technique which we developed earlier for one-dimensional analogues of the problem. We use a spline quasi-interpolation scheme to create a surface approximation from the section approximations. The result is accurate even when the surface is discontinuous across the curved boundary and adjacent side boundaries.     

Regular rapidly decreasing nonlinear generalized functions. Application to microlocal regularity

We present new types of regularity for nonlinear generalized functions, based on the notion of regular growth with respect to the regularizing parameter of the Colombeau simplified model. This generalizes the notion of G^∞-regularity introduced by M. Oberguggenberger. A key point is that these regularities can be characterized, for compactly supported generalized functions, by a property of their Fourier transform. This opens the door to microanalysis of singularities of generalized functions, with respect to these regularities. We present a complete study of this topic, including properties of the Fourier transform (exchange and regularity theorems) and relationship with classical theory, via suitable results of embeddings.     

On a convergence of the rational-trigonometric-polynomial approximations realized by the roots of the Laguerre polynomials

The paper considers a problem of approximation of functions by means of their finite number of Fourier coefficients. Convergence acceleration of approximations by the truncated Fourier series is achieved by application of polynomial and rational correction functions. Rational corrections include unknown parameters whose determination is a crucial problem. We consider an approach connected with the roots of the Laguerre polynomials and study the rates of convergence of such approximations.     

Effective numerical evaluation of the double Hilbert transform

In this paper, we propose two methods to compute the double Hilbert transform of periodic functions. First, we establish the quadratic formula of trigonometric interpolation type for double Hilbert transform and obtain an estimation of the remainder. We call this method 2D mechanical quadrature method (2D-MQM). Numerical experiments show that 2D-MQM outperforms the library function “hilbert” in Matlab when the values of the functions being handled are very large or approach to infinity. Second, we propose a complex analytic method to calculate the double Hilbert transform, which is based on the 2D adaptive Fourier decomposition, and the method is called as 2D-HAFD. In contrast to the pointwise approximation, 2D-HAFD provides explicit rational functional approximations and is valid for all signals of finite energy.     

A note on the discrete Cauchy‐Kovalevskaya extension

In this paper, we exploit the umbral calculus framework to reformulate the so‐called discrete Cauchy‐Kovalevskaya extension in the scope of hypercomplex variables. The key idea is to consider not only formal power series representation for the underlying solution, but also integral representations for the Chebyshev polynomials of first and second kind by means of its Cauchy principal values. It turns out that the resulting integral representation associated to our toy problem is a space‐time Fourier type inversion formula. Moreover, with the aid of some Laplace transform identities involving the generalized Mittag‐Leffler function, we are able to establish a link with a Cauchy problem of differential‐difference type.     

Fast computation of the field of a vertical magnetic dipole above a layered medium

We consider the direct problem of aeroelectrosounding for a layered medium. The specific feature of this problem is that its solution is representable as a Bessel transform of a function ƒ of a certain type. We have developed and implemented a method for fast computation of the Bessel transform of such functions. The method relies on the approximation of ƒ by a linear combination of exponentials with prespecified exponents. The Bessel transform of the original function ƒ is computed analytically as the Bessel transform of the approximating linear combination. Numerical computations have been carried out. __________ Translated from Prikladnaya Matematika i Informatika, No. 26, pp. 39–48, 2007.     

Windowed Fourier transform of two-dimensional quaternionic signals

In this paper, we generalize the classical windowed Fourier transform (WFT) to quaternion-valued signals, called the quaternionic windowed Fourier transform (QWFT). Using the spectral representation of the quaternionic Fourier transform (QFT), we derive several important properties such as reconstruction formula, reproducing kernel, isometry, and orthogonality relation. Taking the Gaussian function as window function we obtain quaternionic Gabor filters which play the role of coefficient functions when decomposing the signal in the quaternionic Gabor basis. We apply the QWFT properties and the (right-sided) QFT to establish a Heisenberg type uncertainty principle for the QWFT. Finally, we briefly introduce an application of the QWFT to a linear time-varying system.     

Generalized fourier transform and its applications

We consider the generalized Fourier transform treated as an operator on the dual of an arbitrary locally convex space. We give a definition of this operator and establish its basic properties. Special attention is paid to cases in which the range of the generalized Fourier transform coincides with a weighted space of entire functions. The results are applied to finding the orders and types of operators in various spaces.     

Decomposition theorems and approximation by a ``floating" system of exponentials

The main problem considered in this paper is the approximation of a trigonometric polynomial by a trigonometric polynomial with a prescribed number of harmonics. The method proposed here gives an opportunity to consider approximation in different spaces, among them the space of continuous functions, the space of functions with uniformly convergent Fourier series, and the space of continuous analytic functions. Applications are given to approximation of the Sobolev classes by trigonometric polynomials with prescribed number of harmonics, and to the widths of the Sobolev classes. This work supplements investigations by Maiorov, Makovoz and the author where similar results were given in the integral metric.

     

A wavelet method for solving singular integral equation of MHD

This paper presents Haar wavelet approximation to solve a singular integral equation which has singularities on a diagonal of the domain R=[-1,1]×[-1,1]. The singularities arise basically due to modified Bessel function K₀ which appears as a part of the kernel. Thus the integral equation is weakly (logarithmically) singular only. The problem is solved considering all the singularities of the kernel and results are examined for approximations of different orders. Our interest to solve the problem using Haar wavelet basis is due to its simplicity and efficiency in numerical approximation. The results show convergence trend as mesh is refined. Comparison is made with some available results obtained earlier by partial consideration of the singularities.     

Computation of the fractional Fourier transform

In this paper we make a critical comparison of some programs for the digital computation of the fractional Fourier transform that are freely available and we describe our own implementation that filters the best out of the existing ones. Two types of transforms are considered: first, the fast approximate fractional Fourier transform algorithm for which two algorithms are available. The method is described in [H.M. Ozaktas, M.A. Kutay, G. Bozda i, IEEE Trans. Signal Process. 44 (1996) 2141–2150]. There are two implementations: one is written by A.M. Kutay, the other is part of package written by J. O'Neill. Second, the discrete fractional Fourier transform algorithm described in the master thesis by Ç. Candan [Bilkent University, 1998] and an algorithm described by S.C. Pei, M.H. Yeh, and C.C. Tseng [IEEE Trans. Signal Process. 47 (1999) 1335–1348].     

Integral representation of slowly growing equidistant splines

In this paper we consider polynomial splines S(x) with equidistant nodes which may grow as O (|x|^s). We present an integral representation of such splines with a distribution kernel. This representation is related to the Fourier integral of slowly growing functions. The part of the Fourier exponentials herewith play the so called exponential splines by Schoenberg. The integral representation provides a flexible tool for dealing with the growing equidistant splines. First, it allows us to construct a rich library of splines possessing the property that translations of any such spline form a basis of corresponding spline space. It is shown that any such spline is associated with a dual spline whose translations form a biorthogonal basis. As examples we present solutions of the problems of projection of a growing function onto spline spaces and of spline interpolation of a growing function. We derive formulas for approximate evaluation of splines projecting a function onto the spline space and establish therewith exact estimations of the approximation errors.     

The Mehler Formula for the Generalized Clifford-Hermite Polynomials

The Mehler formula for the Hermite polynomials allows for an integral representation of the one-dimensional Fractional Fourier transform. In this paper, we introduce a multi-dimensional Fractional Fourier transform in the framework of Clifford analysis. By showing that it coincides with the classical tensorial approach we are able to prove Mehler＇s formula for the generalized Clifford-Hermite polynomials of Clifford analysis.     

Growth properties of Nevanlinna matrices for rational moment problems

We consider rational moment problems on the real line with their associated orthogonal rational functions. There exists a Nevanlinna-type parameterization relating to the problem, with associated Nevanlinna matrices of functions having singularities in the closure of the set of poles of the rational functions belonging to the problem. We prove results related to the growth at the singularities of the functions in a Nevanlinna matrix, and in particular provide bounds on the growth analogous to the corresponding result in the classical polynomial case, when the number of singularities is finite.     

A generalization of the Wiener rational basis functions on infinite intervals,Part II — Numerical investigation

In Part I we introduced the generalized Wiener rational basis functions, and here in Part II we continue our investigation with numerical experiments. Wiener’s generalized basis can utilize the fast Fourier transform for integer values of the decay parameter s$s$; we outline two algorithms for doing so. In addition, the issue of Galerkin representations for polynomial nonlinearities of expansions is addressed.     

Fourier transform and distributional representation of the generalized gamma function with some applications

We present a Fourier transform representation of the generalized gamma functions, which leads to a distributional representation for them as a series of Dirac-delta functions. Applications of these representations are shown in evaluation of the integrals of products of the generalized gamma function with other functions. The results for Euler’s gamma function are deduced as special cases. The relation of the generalized gamma function with the Macdonald function is exploited to deduce new identities for it.     

Finite element method for elliptic problems with edge singularities

We consider tangentially regular solution of the Dirichlet problem for an homogeneous strongly elliptic operator with constant coefficients, on an infinite vertical polyhedral cylinder based on a bounded polygonal domain in the horizontal xy-plane. The usual complex blocks of singularities in the non-tensor product singular decomposition of the solution are made more explicit by a suitable choice of the regularizing kernel. This permits to design a well-posed semi-discrete singular function method (SFM), which differs from the usual SFM in that the dimension of the space of trial and test functions is infinite. Partial Fourier transform in the z-direction (of edges) enables us to overcome the difficulty of an infinite dimension and to obtain optimal orders of convergence in various norms for the semi-discrete solution. Asymptotic error estimates are also proved for the coefficients of singularities. For practical computations, an optimally convergent full-disc! ! ! retization approach, which consists in coupling truncated Fourier series in the z-direction with the SFM in the xy-plane, is implemented. Other good (though not optimal) schemes, which are based on a tensor product form of singularities are investigated. As an illustration of the results, we consider the Laplace operator.     

An exact far-field inversion for the Born approximation and plane wave decompositions for Hankel functions

We consider the Born approximation (representative for first-order approximations) of the scattering problem for the scalar Helroholtz equation with a fixed real-valued free-space wavenumber and a complex-valued compactly supported potential. The boundary condition is the Sommerfeld radiation condition. We derive an exact series-integral representation of the potential from the Fourier coefficients of its far-field pattern, suitable for discussion of the connected stability problem. Furthermore we stress the connection between this representation and some plane wave decompositions for Hankel functions. Without loss of generality we restrict ourselves to the case of two space dimensions.     

Fourier transform over finite field and identities between Gauss sums

This paper is a sequel to [7]. Here we study identities for the Fourier transform of "elementary functions" over finite fields containing "exponentials" of rational monomial functions. It turns out that these identities are governed by monomial identities between Gauss sums. We show that similar to the case of complex numbers such identities correspond to linear relations between certain divisors on the space of multiplicative characters.