Closed Form for Some Gaussian Convolutions

Total Page:16

File Type:pdf, Size:1020Kb

Closed Form for Some Gaussian Convolutions Closed Form for Some Gaussian Convolutions Hossein Mobahi Computer Science & Artificial Intelligence Lab. Massachusetts Institute of Technology Cambridge, MA, USA [email protected] Abstract The convolution of a function with an isotropic Gaussian appears in many contexts such as differential equations, computer vision, signal pro- cessing, and numerical optimization. Although this convolution does not always have a closed form expression, there are important family of func- tions for which closed form exists. This article investigates some of such cases. 1 Introduction Motivation. The convolution of a function with an isotropic Gaussian, aka Weierstrass transform, appears in many contexts. In differential equations, it appears as the solution of diffusion equations [Widder, 1975]. In computer vision and image processing, it is used to blur images for generating multi- scale or scale-space representations [Koenderink, 1984]. In addition, deriva- tives of the blurred images are used for defining visual operations and fea- tures such as in SIFT [Lowe, 1999]. In signal processing, it defines filters for noise reduction [Blinchikoff and Zverev, 1976] and preventing aliasing in sam- pled data [Foley et al., 1990]. Finally, in optimization theory, it is used to simplify a nonconvex objective function by suppressing spurious critical points [Mobahi and Fisher, 2015]. arXiv:1602.05610v2 [math.CA] 7 Mar 2016 Numerical Approximation. In general, the Weierstrass transform does not have a closed form expression. When that is the case, the transformation needs to be computed or approximated numerically. Recent efforts, mainly moti- vated from optimization applications, have resulted in efficient algorithms for such numerical approximations, e.g. by adaptive importance sampling [Maggiar et al., 2015]. Exact Closed Form. Putting the general arguments aside, there are impor- tant family of functions for which the Weierstrass transform does have a closed form. This article investigates two such families, namely, multivariate polyno- mials and Gaussian radial basis functions. These two classes are rich, as they 1 are able to approximate any arbitrary function (under some mild conditions), up to a desired accuracy. In addition, we present a result which provides the closed form Weierstrass transform of a multivariate function, given the closed form transform of its univariate cases. Notation. We use x for scalars, x for vectors, . for the ℓ2 norm, and i for √ 1. We also use , for equality by definition. Wek define the isotropic Gaussian − kernel kσ(x) as below. x 2 dim(x) k k kσ(x) , (√2πσ)− e− 2σ2 . 2 Closed Form Expressions 2.1 Polynomials Given a multivariate polynomial function f : Rn R. We show that the → convolution of f with kσ has a closed form expression. Suppose f consists of m monomials, m f(x) , aj fj (x) , Xj=1 where fj is the j’th term and aj is its coefficient. Each monomial has the , n pd,j form fj (x) Πd=1xd with xd being the d’th component of the vector x, and pd,j being a nonnegative integer exponent of xd in fj . Then we have, m m m n pd,j n [f⋆kσ](x)= aj [fj⋆kσ](x) = ajΠd=1 [ d ⋆kσ](xd) = aj Πd=1u(xd,pd,j, σ) . Xj=1 Xj=1 Xj=1 The function u has the following form (see the Appendix for the proof), x u(x, p, σ) , [p ⋆ k ](x) = (iσ)p h ( ) , σ p iσ where hp is the p’th Hermite polynomial, i.e., 2 p 2 p x d x hp(x) , ( 1) e 2 e− 2 . − dxp The following table lists u(x, p, σ) for 0 p 10. ≤ ≤ 2 p u(x, p, σ) 0 1 1 x 2 σ2 + x2 3 3σ2x + x3 4 3σ4 +6σ2x2 + x4 5 15σ4x + 10σ2x3 + x5 6 15σ6 + 45σ4x2 + 15σ2x4 + x6 7 105σ6x + 105σ4x3 + 21σ2x5 + x7 8 105σ8 + 420σ6x2 + 210σ4x4 + 28σ2x6 + x8 9 945σ8x + 1260σ6x3 + 378σ4x5 + 36σ2x7 + x9 10 945σ10 + 4725σ8x2 + 3150σ6x4 + 630σ4x6 + 45σ2x8 + x10 2 3 4 4 Example Let f(x, y) , x y +λ(x +y ). Then, using the table, [f ⋆kσ](x, y) has the following form, 2 2 2 3 4 2 2 4 4 2 2 4 [f ⋆ kσ] = (σ + x )(3σ y + y )+ λ(3σ +6σ x + x +3σ +6σ y + y ) . 2.2 Gaussian Radial Basis Functions Isotropic Gaussian Radial-Basis-Functions (RBFs) are general function approx- imators [Hartman et al., 1990]. That is, any function can be approximated to a desired accuracy by the following form, 2 m kx−xj k − 2δ2 f(x) , aj e j , Xj=1 n where xj R is the center and δj R+ is the width of the j’th basis function, for j∈ = 1,...,m. We have shown∈ in [Mobahi et al., 2012] that the convolution of f with kσ has the following form, 2 m kx−xj k δj n − 2(δ2 +σ2) [f ⋆ kσ](x)= aj ( ) e j . 2 2 Xj=1 δj + σ q In order to be self-contained, we include our proof as presented in [Mobahi et al., 2012] in the appendix of this article. 2.3 Trigonometric Functions σ2 sin(x) e− 2 sin(x) σ2 cos(x) e− 2 cos(x) 2 1 2σ2 sin (x) 2 1 e− cos(2x) 2 1 − 2σ2 cos (x) 2 1+ e− cos(2x) 2σ2 sin(x)cos(x) e− sin(x)cos(x) 3 Example Let f(x, y) , 2 cos(x) 3 sin(y)+4cos(x) sin(y) 8 sin(x)cos(y)+ 2 − − 4 . Compute [f ⋆ kσ](x, y). We expand the quadratic form and then use the table, as shown below. 2 2 cos(x) 3 sin(y)+4cos(x) sin(y) 8 sin(x)cos(y)+4 (1) − − = 4cos 2(x) + 9 sin2(y)+16cos2(x) sin2(y) + 64 sin2(x)cos2(y)+16 (2) 12 cos(x) sin(y)+16cos2(x) sin(y) 32 cos(x) sin(x)cos(y) (3) − − +16 cos(x) 24 sin(y) cos(x) + 48 sin(y) sin(x)cos(y) 24 sin(y) (4) − − 64 cos(x) sin(y) sin(x) + 32cos(x) sin(y) 64 sin(x)cos(y) (5) − − 2.4 Functions with Linear Argument T This result tells us about the form of [f x ⋆ kσ()](w) given [f ⋆ kσ](y). Lemma 1 Let f : R R have well-defined Gaussian convolution f˜(y ; σ) , → T [f ⋆kσ](y). Let x and w be any real vectors of the same length. Then [f x ⋆ T kσ()](w)= f˜(w x ; σ x ). We now use this resultk ink the following examples. These examples are related to the earlier Table about the smoothed response functions. Example Suppose f(y) , sign(y). It is easy to check that f˜(y; σ) = erf( y ). √2σ T Applying the lemma to this implies that [sign(T x)⋆k ()](w) = erf( w x ). σ √2σ x k k Example Suppose f(y) , max(0,y). It is easy to check that f˜(y; σ) = y2 σ 1 y T e− 2σ2 + y 1+erf( ) . Applying the lemma to this implies that [max(0, x)⋆ √2π 2 √2σ T 2 (w x) T σ x 2 2 1 T w x k ()](w)= k k e− 2σ kxk + w x 1+erf( ) . σ √2π 2 √2σ x k k Example Suppose f(y) , sin(y). From trigonometric table in 2.3 we know 2 σ T f˜(y; σ) = e− 2 sin(y). Applying the lemma to this implies that [sin( x) ⋆ 2 2 σ kxk T kσ()](w)= e− 2 sin(w x). In particular, when x , (x, 1), then smoothed σ2(1+x2) version of sin(w1x + w2) w.r.t. w is e− 2 sin(w1x + w2). References [Blinchikoff and Zverev, 1976] Blinchikoff, H. and Zverev, A. (1976). Filtering in the time and frequency domains. R.E. Krieger Pub. Co. [Foley et al., 1990] Foley, J. D., van Dam, A., Feiner, S. K., and Hughes, J. F. (1990). Computer Graphics: Principles and Practice (2Nd Ed.). Addison- Wesley Longman Publishing Co., Inc., Boston, MA, USA. 4 [Hartman et al., 1990] Hartman, E., Keeler, J. D., and Kowalski, J. M. (1990). Layered neural networks with gaussian hidden units as universal approxima- tions. Neural Comput., 2(2):210–215. [Koenderink, 1984] Koenderink, J. J. (1984). The structure of images. Biological Cybernetics, 50(5):363–370. [Lowe, 1999] Lowe, D. G. (1999). Object recognition from local scale-invariant features. In Proceedings of the International Conference on Computer Vision, ICCV ’99. [Maggiar et al., 2015] Maggiar, A., Waechter, A., Dolinskaya, I. S., and Staum, J. (2015). A derivative-free trust-region algorithm for the optimization of functions smoothed via gaussian convolution using adaptive multiple impor- tance sampling. Technical Report IEMS Department, Northwestern Univer- sity. [Mobahi and Fisher, 2015] Mobahi, H. and Fisher, III, J. W. (2015). On the link between gaussian homotopy continuation and convex envelopes. In Lecture Notes in Computer Science (EMMCVPR 2015). Springer. [Mobahi et al., 2012] Mobahi, H., Zitnick, C. L., and Ma, Y. (2012). Seeing through the blur. In 2012 IEEE Computer Society Conference on Computer Vision and Pattern Recognition (CVPR 2012), Rhode Island, USA. IEEE Computer Society. [Widder, 1975] Widder, D. V. (1975). The Heat Equation. Academic Press. 5 6 Appendices A Polynomials p p x , Here we prove that [ ⋆ kσ](x) = (iσ) hp( iσ ), where i √ 1 and hp is the 2 2 − x dp x , p 2 2 p’th Hermite polynomial hp(x) ( 1) e dxp e− . Using Fourier Transform F we proceed as below, − p [ ⋆ kσ](x) (6) 1 p = F − F [ ⋆ kσ](x) (7) 1{ { p }} = F − F x F kσ(x) (8) { { } { }} 1 p (p) 1 1 σ2ω2 = F − √2π( i) δ (ω) e− 2 (9) { − √2π } p 1 (p) 1 1 σ2ω2 = √2π( i) F − δ (ω) e− 2 (10) − { √2π } p (p) 1 1 σ2ω2 iωx = √2π( i) δ (ω) e− 2 e dω (11) − ZR √2π p p d 1 1 σ2ω2+iωx √ 2 = 2π( i) p e− (12) − dω √2π ω=0 p p d 1 1 σ2ω2+iωx √ 2 = 2π( i) p e− (13) − dω √2π ω=0 p d 1 2 2 p 2 σ ω +iωx = ( i) p e− (14) − dω ω=0 (ω− 1 ix)2 p x2 σ2 p d 2σ2 2 − − σ2 = ( i) p e (15) − dω ω=0 ω 1 ix 2 ( − 2 ) 2 p σ p x d − 2 = ( i) e− 2σ2 e σ2 (16) − dωp ω=0 2 1 2 2 σ (ω− ix) p x 1 p σ2 = ( i) e− 2σ2 hp(σ(ω ix))( σ) e− 2 (17) − − σ2 − ω=0 x = (iσ)ph ( ) .
Recommended publications
  • Ballistic Electron Transport in Graphene Nanodevices and Billiards
    Ballistic Electron Transport in Graphene Nanodevices and Billiards Dissertation (Cumulative Thesis) for the award of the degree Doctor rerum naturalium of the Georg-August-Universität Göttingen within the doctoral program Physics of Biological and Complex Systems of the Georg-August University School of Science submitted by George Datseris from Athens, Greece Göttingen 2019 Thesis advisory committee Prof. Dr. Theo Geisel Department of Nonlinear Dynamics Max Planck Institute for Dynamics and Self-Organization Prof. Dr. Stephan Herminghaus Department of Dynamics of Complex Fluids Max Planck Institute for Dynamics and Self-Organization Dr. Michael Wilczek Turbulence, Complex Flows & Active Matter Max Planck Institute for Dynamics and Self-Organization Examination board Prof. Dr. Theo Geisel (Reviewer) Department of Nonlinear Dynamics Max Planck Institute for Dynamics and Self-Organization Prof. Dr. Stephan Herminghaus (Second Reviewer) Department of Dynamics of Complex Fluids Max Planck Institute for Dynamics and Self-Organization Dr. Michael Wilczek Turbulence, Complex Flows & Active Matter Max Planck Institute for Dynamics and Self-Organization Prof. Dr. Ulrich Parlitz Biomedical Physics Max Planck Institute for Dynamics and Self-Organization Prof. Dr. Stefan Kehrein Condensed Matter Theory, Physics Department Georg-August-Universität Göttingen Prof. Dr. Jörg Enderlein Biophysics / Complex Systems, Physics Department Georg-August-Universität Göttingen Date of the oral examination: September 13th, 2019 Contents Abstract 3 1 Introduction 5 1.1 Thesis synopsis and outline.............................. 10 2 Fundamental Concepts 13 2.1 Nonlinear dynamics of antidot superlattices..................... 13 2.2 Lyapunov exponents in billiards............................ 20 2.3 Fundamental properties of graphene......................... 22 2.4 Why quantum?..................................... 31 2.5 Transport simulations and scattering wavefunctions................
    [Show full text]
  • Weierstrass Transform on Howell's Space-A New Theory
    IOSR Journal of Mathematics (IOSR-JM) e-ISSN: 2278-5728, p-ISSN: 2319-765X. Volume 16, Issue 6 Ser. IV (Nov. – Dec. 2020), PP 57-62 www.iosrjournals.org Weierstrass Transform on Howell’s Space-A New Theory 1 2* 3 V. N. Mahalle , S.S. Mathurkar , R. D. Taywade 1(Associate Professor, Dept. of Mathematics, Bar. R.D.I.K.N.K.D. College, Badnera Railway (M.S), India.) 2*(Assistant Professor, Dept. of Mathematics, Govt. College of Engineering, Amravati, (M.S), India.) 3(Associate Professor, Dept. of Mathematics, Prof. Ram Meghe Institute of Technology & Research, Badnera, Amravati, (M.S), India.) Abstract: In this paper we have explained a new theory of Weierstrass transform defined on Howell’s space. The definition of the Weierstrass transform on this space of test functions is discussed. Some Fundamental theorems are proved and some results are developed. Weierstrass transform defined on Howell’s space is linear, continuous and one-one mapping from G to G c is also discussed. Key Word: Weierstrass transform, Howell’s Space, Fundamental theorem. --------------------------------------------------------------------------------------------------------------------------------------- Date of Submission: 10-12-2020 Date of Acceptance: 25-12-2020 --------------------------------------------------------------------------------------------------------------------------------------- I. Introduction A new theory of Fourier analysis developed by Howell, presented in an elegant series of papers [2-6] n supporting the Fourier transformation of the functions on R . This theory is similar to the theory for tempered distributions and to some extent, can be viewed as an extension of the Fourier theory. Bhosale [1] had extended Fractional Fourier transform on Howell’s space. Mahalle et. al. [7, 8] described a new theory of Mellin and Laplace transform defined on Howell’s space.
    [Show full text]
  • Memristor-Based Approximation of Gaussian Filter
    Memristor-based Approximation of Gaussian Filter Alex Pappachen James, Aidyn Zhambyl, Anju Nandakumar Electrical and Computer Engineering department Nazarbayev University School of Engineering Astana, Kazakhstan Email: [email protected], [email protected], [email protected] Abstract—Gaussian filter is a filter with impulse response of fourth section some mathematical analysis will be provided to Gaussian function. These filters are useful in image processing obtained results and finally the paper will be concluded by of 2D signals, as it removes unnecessary noise. Also, they could stating the result of the done work. be helpful for data transmission (e.g. GMSK modulation). In practice, the Gaussian filters could be approximately designed II. GAUSSIAN FILTER AND MEMRISTOR by several methods. One of these methods are to construct Gaussian-like filter with the help of memristors and RLC circuits. A. Gaussian Filter Therefore, the objective of this project is to find and design Gaussian filter is such filter which convolves the input signal appropriate model of Gaussian-like filter, by using mentioned devices. Finally, one possible model of Gaussian-like filter based with the impulse response of Gaussian function. In some on memristor designed and analysed in this paper. sources this process is also known as the Weierstrass transform [2]. The Gaussian function is given as in equation 1, where µ Index Terms—Gaussian filters, analog filter, memristor, func- is the time shift and σ is the scale. For input signal of x(t) tion approximation the output is the convolution of x(t) with Gaussian, as shown in equation 2.
    [Show full text]
  • Eigenfunctions and Eigenoperators of Cyclic Integral Transforms with Application to Gaussian Beam Propagation
    NORTHWESTERN UNIVERSITY Eigenfunctions and Eigenoperators of Cyclic Integral Transforms with Application to Gaussian Beam Propagation A DISSERTATION SUBMITTED TO THE GRADUATE SCHOOL IN PARTIAL FULFILLMENT OF THE REQUIREMENTS for the degree DOCTOR OF PHILOSOPHY Field of Applied Mathematics By Matthew S. McCallum EVANSTON, ILLINOIS December 2008 2 c Copyright by Matthew S. McCallum 2008 All Rights Reserved 3 ABSTRACT Eigenfunctions and Eigenoperators of Cyclic Integral Transforms with Application to Gaussian Beam Propagation Matthew S. McCallum An integral transform which reproduces a transformable input function after a finite number N of successive applications is known as a cyclic transform. Of course, such a transform will reproduce an arbitrary transformable input after N applications, but it also admits eigenfunction inputs which will be reproduced after a single application of the transform. These transforms and their eigenfunctions appear in various applications, and the systematic determination of eigenfunctions of cyclic integral transforms has been a problem of interest to mathematicians since at least the early twentieth century. In this work we review the various approaches to this problem, providing generalizations of published expressions from previous approaches. We then develop a new formalism, differential eigenoperators, that reduces the eigenfunction problem for a cyclic transform to an eigenfunction problem for a corresponding ordinary differential equation. In this way we are able to relate eigenfunctions of integral equations to boundary-value problems, 4 which are typically easier to analyze. We give extensive examples and discussion via the specific case of the Fourier transform. We also relate this approach to two formalisms that have been of interest to the math- ematical physics community – hyperdifferential operators and linear canonical transforms.
    [Show full text]
  • A Unified Framework for the Fourier, Laplace, Mellin and Z Transforms
    The Generalized Fourier Transform: A Unified Framework for the Fourier, Laplace, Mellin and Z Transforms Pushpendra Singha,∗, Anubha Guptab, Shiv Dutt Joshic aDepartment of ECE, National Institute of Technology Hamirpur, Hamirpur, India. bSBILab, Department of ECE, Indraprastha Institute of Information Technology Delhi, Delhi, India cDepartment of EE, Indian Institute of Technology Delhi, Delhi, India Abstract This paper introduces Generalized Fourier transform (GFT) that is an extension or the generalization of the Fourier transform (FT). The Unilateral Laplace transform (LT) is observed to be the special case of GFT. GFT, as proposed in this work, contributes significantly to the scholarly literature. There are many salient contribution of this work. Firstly, GFT is applicable to a much larger class of signals, some of which cannot be analyzed with FT and LT. For example, we have shown the applicability of GFT on the polynomially decaying functions and super exponentials. Secondly, we demonstrate the efficacy of GFT in solving the initial value problems (IVPs). Thirdly, the generalization presented for FT is extended for other integral transforms with examples shown for wavelet transform and cosine transform. Likewise, generalized Gamma function is also presented. One interesting application of GFT is the computation of generalized moments, for the otherwise non-finite moments, of any random variable such as the Cauchy random variable. Fourthly, we introduce Fourier scale transform (FST) that utilizes GFT with the topological isomorphism of an exponential map. Lastly, we propose Generalized Discrete-Time Fourier transform (GDTFT). The DTFT and unilateral z-transform are shown to be the special cases of the proposed GDTFT. The properties of GFT and GDTFT have also been discussed.
    [Show full text]
  • (1.1) «-V (1.2) «"" «"II 7:-7Z—R-77
    A DISCRETE ANALOGUE OF THE WEIERSTRASS TRANSFORM1 EDWARD NORMAN 1. Introduction and definitions. In this note we consider a discrete analogue of the transform ,TÍ-g}*ttg(y)£yi /(*)=/_; This transform may be considered as a convolution; accordingly our analogue will take the form 00 f(n) = £ k(n — m)g(m) —00 where the kernel sequence k(n) is to be appropriately chosen. The choice of k(n) is motivated as follows. One way in which the kernel e~y '* arises in the theory of the con- volution transform is by consideration of the class of kernels whose bilateral Laplace transform is of the form 6*ñ(i--)x- (1.1) «-V -tlak i \ aj When c = 0 we obtain the class of kernels whose analogue was con- sidered by Pollard and Standish [4] where the kernels formed a sub- class of the totally positive sequences. The totally positive sequences can be characterized as sequences having a generating function of the form (1.2) «""■«"II 1 7:-7Z—r-77(1 - a,z)(i - bjZ-1) (see [l]). The factor eaz+h'~ in (1.2) corresponds to the factor e" in (1.1). Letting a = b = v/2 we have the familiar expansion eC/2) (H-*"1) = Y,Im(v)zm Received by the editors February 10, 1959 and, in revised form, October 22, 1959. 1 This work is part of a doctoral thesis done under the direction of Professor Harry Pollard at Cornell University. * See Hirschman and Widder [3] for the theory of the Weierstrass transform. License or copyright restrictions may apply to redistribution; see https://www.ams.org/journal-terms-of-use596 DISCRETE ANALOGUEOF WEIERSTRASS TRANSFORM 597 where Im(v) is the modified Bessel coefficient of the first kind.
    [Show full text]
  • Introduction to Hyperfunctions and Their Integral Transforms
    [chapter] [chapter] [chapter] 1 2 Introduction to Hyperfunctions and Their Integral Transforms Urs E. Graf 2009 ii Contents Preface ix 1 Introduction to Hyperfunctions 1 1.1 Generalized Functions . 1 1.2 The Concept of a Hyperfunction . 3 1.3 Properties of Hyperfunctions . 13 1.3.1 Linear Substitution . 13 1.3.2 Hyperfunctions of the Type f(φ(x)) . 15 1.3.3 Differentiation . 18 1.3.4 The Shift Operator as a Differential Operator . 25 1.3.5 Parity, Complex Conjugate and Realness . 25 1.3.6 The Equation φ(x)f(x) = h(x)............... 28 1.4 Finite Part Hyperfunctions . 33 1.5 Integrals . 37 1.5.1 Integrals with respect to the Independent Variable . 37 1.5.2 Integrals with respect to a Parameter . 43 1.6 More Familiar Hyperfunctions . 44 1.6.1 Unit-Step, Delta Impulses, Sign, Characteristic Hyper- functions . 44 1.6.2 Integral Powers . 45 1.6.3 Non-integral Powers . 49 1.6.4 Logarithms . 51 1.6.5 Upper and Lower Hyperfunctions . 55 α 1.6.6 The Normalized Power x+=Γ(α + 1) . 58 1.6.7 Hyperfunctions Concentrated at One Point . 61 2 Analytic Properties 63 2.1 Sequences, Series, Limits . 63 2.2 Cauchy-type Integrals . 71 2.3 Projections of Functions . 75 2.3.1 Functions satisfying the H¨olderCondition . 77 2.3.2 Projection Theorems . 78 2.3.3 Convergence Factors . 87 2.3.4 Homologous and Standard Hyperfunctions . 88 2.4 Projections of Hyperfunctions . 91 2.4.1 Holomorphic and Meromorphic Hyperfunctions . 91 2.4.2 Standard Defining Functions .
    [Show full text]
  • Hermite Polynomials 1 Hermite Polynomials
    Hermite polynomials 1 Hermite polynomials In mathematics, the Hermite polynomials are a classical orthogonal polynomial sequence that arise in probability, such as the Edgeworth series; in combinatorics, as an example of an Appell sequence, obeying the umbral calculus; in numerical analysis as Gaussian quadrature; in finite element methods as shape functions for beams; and in physics, where they give rise to the eigenstates of the quantum harmonic oscillator. They are also used in systems theory in connection with nonlinear operations on Gaussian noise. They were defined by Laplace (1810) [1] though in scarcely recognizable form, and studied in detail by Chebyshev (1859).[2] Chebyshev's work was overlooked and they were named later after Charles Hermite who wrote on the polynomials in 1864 describing them as new.[3] They were consequently not new although in later 1865 papers Hermite was the first to define the multidimensional polynomials. Definition There are two different ways of standardizing the Hermite polynomials: • The "probabilists' Hermite polynomials" are given by , • while the "physicists' Hermite polynomials" are given by . These two definitions are not exactly identical; each one is a rescaling of the other, These are Hermite polynomial sequences of different variances; see the material on variances below. The notation He and H is that used in the standard references Tom H. Koornwinder, Roderick S. C. Wong, and Roelof Koekoek et al. (2010) and Abramowitz & Stegun. The polynomials He are sometimes denoted by H , n n especially in probability theory, because is the probability density function for the normal distribution with expected value 0 and standard deviation 1.
    [Show full text]
  • The Weierstrass Transform on Some Spaces of Type S
    P: ISSN No. 0976-8602 RNI No.UPENG/2012/42622 VOL.-V, ISSUE-III, July-2016 E: ISSN No. 2349 - 9443 Asian Resonance The Weierstrass Transform on Some Spaces of Type S Abstract The main objective of this paper is to study the classical theory of the Weierstrass transform is extended to the Dirac delta function and some of their basic properties. The characterizations of the Weierstrass transform on the Gel'fand-Shilov spaces of type S is shown. Keywords: Weierstrass transform, Fourier transform, Green's function, Dirac delta function. Introduction This present paper, my contribution is to define the Weierstrass transform in the form of Dirac delta function and study some of their basic properties. The characterizations of the Weierstrass transform on the Gel'fand-Shilov spaces of type S is shown. The conventional Weierstrass transform of a suitably restricted function fx() on the real axis R is defined by [5, p. 267] as: 1 (tx )2 (1) W f( x ) f ( x )exp dx . 2 R 4 The Green's function G t x, k of the heat equation for the infinite interval is 1 (tx )2 , where k is positive constant. G t x, k exp 2 k 4k Ashutosh Mahato The limit of as 0 represent the Dirac delta Assistant Professor, 2 function 1 (tx ) Deptt.of Mathematics, (tx ) lim exp . C.V.Raman College of 0 2 k 4k Engineering, Bhubaneswar, So, in terms of Dirac delta function ()tx , the Weierstrass Odisha transform (1) can be rewritten as Wfx ()(,)()(). fxtˆ fx txdx (2) R Aim of the Study/ The problem / Objective of the Study The main objective of this paper is to develop the classical theory of the Weierstrass transform is extended to the Dirac delta function and some of their basic properties.
    [Show full text]
  • Weierstrass Transforms
    JAC : A JOURNAL OF COMPOSITION THEORY ISSN : 0731-6755 GENERALIZED PROPERTIES ON MELLIN - WEIERSTRASS TRANSFORMS Amarnath Kumar Thakur Department of Mathematics Dr.C.V.Raman University, Bilaspur, State, Chhattisgarh Email- [email protected] Gopi Sao Department of Mathematics Dr.C.V.Raman University, Bilaspur, State, Chhattisgarh Email- [email protected] Hetram Suryavanshi Department of Mathematics Dr.C.V.Raman University, Bilaspur, State, Chhattisgarh Email- [email protected] Abstract- In the present paper, we introduce definite integral transforms in operational calculus for the generalized properties on Mellin- Weierstrass associated transforms. These results are derived from the MW Properties These formulas may be considered as promising approaches in expressing in fractional calculus. Keywords – Laplace Transform , Convolution Transform Mellin Transforms , Weierstrass Transforms I. Introduction In mathematical terms the Fourier and Laplace transformations were introduced to solve physical problems. In fact, the first event of change is found by Reiman in which he used to study the famous zeta function. The Finnish mathematician Hjalmar Mellin (1854-1933), who was the first to give an orderly appearance change and its reversal,working in the theory of special functions, he developed applications towards a solution of hypergeometric differential equations and derivation of asymptotic expansions. Mellin's contribution gives a major place to the theory of analytic functions and essentially depends on the Cauchy theorem and the method of residuals. Actually, the Mellin transformation can also be placed in another framework, in which some cases is more closely related to Volume XIII, Issue IX, SEPTEMBER 2020 Page No: 98 JAC : A JOURNAL OF COMPOSITION THEORY ISSN : 0731-6755 original ideas of Reimann's original.In addition to its use in mathematics, Mellin's transformation has been applied in many different ways and areas of Physics and Engineering.
    [Show full text]
  • Asymptotic Expansions of Some Integral Transforms by Using Generalized Functions by Ahmed I, Zayed
    TRANSACTIONS OF THE AMERICAN MATHEMATICAL SOCIElY Volume 272, Number 2, August 1982 ASYMPTOTIC EXPANSIONS OF SOME INTEGRAL TRANSFORMS BY USING GENERALIZED FUNCTIONS BY AHMED I, ZAYED ABSTRACT. The technique devised by Wong to derive the asymptotic expansions of multiple Fourier transforms by using the theory of Schwartz distributions is ex- tended to a large class of integral transforms. The extension requires establishing a general procedure to extend these integral transforms to generalized functions. Wong's technique is then applied to some of these integral transforms to obtain their asymptotic expansions. This class of integral transforms encompasses, among others, the Laplace, the Airy, the K and the Hankel transforms. 1. Introduction. Recently, the theory of generalized functions has been introduced into the field of asymptotic expansions of integral transforms; see [18, 19,6]. One of the advantages of this approach is to be able to interpret and to assign values to some divergent integrals. This approach goes back to Lighthill [10] who was the first to use generalized functions to study the asymptotic expansions of Fourier transforms. Borrowing Lighthill's idea, several people were able to derive asymptotic expan- sions for other types of integral transforms such as Laplace [2], Stieltjes [11], Riemann-Liouville fractional integral transforms [12] and Mellin convolution of two algebraically dominated functions [17]. The use of generalized functions in the theory of asymptotic expansions has another advantage; that is, in some cases it can provide explicit expressions for the error terms in shorter and more direct ways than the classical methods [9,11]. However, one should not overemphasize the importance of the generalized functions approach since, for the most part, it is an alternative but not an exclusive technique.
    [Show full text]
  • On the Summability of Series of Hermite Polynomials
    JOURNAL OF MATHEMATICAL ANALYSIS AND APPLICATIONS 8,406-422 (1964) On the Summability of Series of Hermite Polynomials G. G. BILODEAU Boston College, Chestnut Hill, Massachusetts Submitted by R. P. Boas I. INTRODUCTION Let d(x) be in L for every finite interval and assume that exists for n = 0, 1, 2, mm*where HJx) is the Hermite polynomial defined by H,(x) = (- 1)” exeDn e& . (1.2) Then r+(x) - 3 a,(2%!)-1 H,(x) . (1.3) n=O It is known (see [l, p. 3811 or [2, p. 2401) that the series in (1.3) converges to q%(x)for suitable restrictions on 4 near x and f 00. The series is also sum- mable to 4(.x) by various methods and under less restrictive conditions on 4. Abel summability was considered in 13, p. 4581, Bore1 summabilities in [4], Cesaro and others in [5]. (See footnote in [2, p. 2341 on ref. 5.) The purpose of this paper is to find a summability method which will sum the series to $(x) under conditions on 4 near x = f 03 which do little more than assume the existence of each of the coefficients a,. In Section II, we discuss a previously used summability and show that the result already known cannot be greatly improved. In Section III another summability method is introduced and the main result of the paper is proved. An example is then given to show that this result is very nearly “best possible” and an application of the result to the inversion of the Weierstrass transform is indicated.
    [Show full text]