Computing the Laplace transform and the convolution for more functions adjoined

Similar documents
Computing the inverse Laplace transform for rational functions vanishing at infinity

f(t)e st dt. (4.1) Note that the integral defining the Laplace transform converges for s s 0 provided f(t) Ke s 0t for some constant K.

Generalized sources (Sect. 6.5). The Dirac delta generalized function. Definition Consider the sequence of functions for n 1, Remarks:

Laplace Transform. Chapter 4

Laplace Transform Introduction

Chapter 6: The Laplace Transform 6.3 Step Functions and

The Laplace transform

Reading assignment: In this chapter we will cover Sections Definition and the Laplace transform of simple functions

(f g)(t) = Example 4.5.1: Find f g the convolution of the functions f(t) = e t and g(t) = sin(t). Solution: The definition of convolution is,

e st f (t) dt = e st tf(t) dt = L {t f(t)} s

Ordinary differential equations

2.161 Signal Processing: Continuous and Discrete Fall 2008

Computing inverse Laplace Transforms.

Name: Solutions Final Exam

ADVANCED ENGINEERING MATHEMATICS MATLAB

Math Shifting theorems

MULTI-RESTRAINED STIRLING NUMBERS

ENGIN 211, Engineering Math. Laplace Transforms

LTI Systems (Continuous & Discrete) - Basics

Unit 2: Modeling in the Frequency Domain Part 2: The Laplace Transform. The Laplace Transform. The need for Laplace

Finite-dimensional spaces. C n is the space of n-tuples x = (x 1,..., x n ) of complex numbers. It is a Hilbert space with the inner product

HIGHER-ORDER LINEAR ORDINARY DIFFERENTIAL EQUATIONS IV: Laplace Transform Method David Levermore Department of Mathematics University of Maryland

Math 2C03 - Differential Equations. Slides shown in class - Winter Laplace Transforms. March 4, 5, 9, 11, 12, 16,

HIGHER-ORDER LINEAR ORDINARY DIFFERENTIAL EQUATIONS IV: Laplace Transform Method. David Levermore Department of Mathematics University of Maryland

be the set of complex valued 2π-periodic functions f on R such that

Special Mathematics Laplace Transform

(an improper integral)

Lecture 1 January 5, 2016

Lecture 4: Numerical solution of ordinary differential equations

Differential equations, comprehensive exam topics and sample questions

Reading assignment: In this chapter we will cover Sections Definition and the Laplace transform of simple functions

We denote the space of distributions on Ω by D ( Ω) 2.

FORMAL GROUPS OF CERTAIN Q-CURVES OVER QUADRATIC FIELDS

11.8 Power Series. Recall the geometric series. (1) x n = 1+x+x 2 + +x n +

HW 9, Math 319, Fall 2016

Practice Problems For Test 3

MA Ordinary Differential Equations and Control Semester 1 Lecture Notes on ODEs and Laplace Transform (Parts I and II)

2.161 Signal Processing: Continuous and Discrete Fall 2008

Vectors in Function Spaces

Notes on uniform convergence

Math 308 Exam II Practice Problems

ECEN 420 LINEAR CONTROL SYSTEMS. Lecture 2 Laplace Transform I 1/52

Practice Problems For Test 3

Math 201 Lecture 17: Discontinuous and Periodic Functions

Laplace Transforms Chapter 3

MATRIX INTEGRALS AND MAP ENUMERATION 2

Complex Inversion Formula for Stieltjes and Widder Transforms with Applications

Math 3313: Differential Equations Laplace transforms

Exam Basics. midterm. 1 There will be 9 questions. 2 The first 3 are on pre-midterm material. 3 The next 1 is a mix of old and new material.

MATH 1A, Complete Lecture Notes. Fedor Duzhin

arxiv: v1 [math.ca] 6 May 2016

JUHA KINNUNEN. Partial Differential Equations

SOME FACTORIZATIONS IN THE TWISTED GROUP ALGEBRA OF SYMMETRIC GROUPS. University of Rijeka, Croatia

MATH 205B NOTES 2010 COMMUTATIVE ALGEBRA 53

CITY UNIVERSITY LONDON

1 Math 241A-B Homework Problem List for F2015 and W2016

IMPROVEMENTS OF COMPOSITION RULE FOR THE CANAVATI FRACTIONAL DERIVATIVES AND APPLICATIONS TO OPIAL-TYPE INEQUALITIES

Lecture 7: Laplace Transform and Its Applications Dr.-Ing. Sudchai Boonto

Taylor and Maclaurin Series

9.2 The Input-Output Description of a System

Part V. 17 Introduction: What are measures and why measurable sets. Lebesgue Integration Theory

1 Review of di erential calculus

Topics in Fourier analysis - Lecture 2.

L2 gains and system approximation quality 1

Recall that any inner product space V has an associated norm defined by

Mathematical Foundations of Signal Processing

1 Signals and systems

Considering our result for the sum and product of analytic functions, this means that for (a 0, a 1,..., a N ) C N+1, the polynomial.

2 Classification of Continuous-Time Systems

Lecture 6 January 21, 2016

2tdt 1 y = t2 + C y = which implies C = 1 and the solution is y = 1

Notions such as convergent sequence and Cauchy sequence make sense for any metric space. Convergent Sequences are Cauchy

Name: Solutions Exam 3

ALEXEY V. GAVRILOV. I p (v) : T p M T Expp vm the operator of parallel transport along the geodesic γ v. Let E p be the map

Math 341 Fall 2008 Friday December 12

e (x y)2 /4kt φ(y) dy, for t > 0. (4)

On Tornheim's double series

Hilbert Spaces. Contents

Fourier Series and Transforms. Revision Lecture

Convolution. Define a mathematical operation on discrete-time signals called convolution, represented by *. Given two discrete-time signals x 1, x 2,

Math 20D: Form B Final Exam Dec.11 (3:00pm-5:50pm), Show all of your work. No credit will be given for unsupported answers.

X b n sin nπx L. n=1 Fourier Sine Series Expansion. a n cos nπx L 2 + X. n=1 Fourier Cosine Series Expansion ³ L. n=1 Fourier Series Expansion

4. Noether normalisation

THE S 1 -EQUIVARIANT COHOMOLOGY RINGS OF (n k, k) SPRINGER VARIETIES

7: FOURIER SERIES STEVEN HEILMAN

Hecke Operators, Zeta Functions and the Satake map

MATH MEASURE THEORY AND FOURIER ANALYSIS. Contents

A sufficient condition for the existence of the Fourier transform of f : R C is. f(t) dt <. f(t) = 0 otherwise. dt =

Chapter 5. Modular arithmetic. 5.1 The modular ring

MA 201, Mathematics III, July-November 2016, Laplace Transform

Systems Analysis and Control

Laplace Transform Theory - 1

COMPOSITIONAL INVERSION OF TRIANGULAR SETS OF SERIES. Marilena BARNABEI

Harmonic Analysis on the Cube and Parseval s Identity

Conformal maps. Lent 2019 COMPLEX METHODS G. Taylor. A star means optional and not necessarily harder.

Existence Theory: Green s Functions

An approach to the Biharmonic pseudo process by a Random walk

Math 210B. Artin Rees and completions

Arithmetic Funtions Over Rings with Zero Divisors

A CONSTRUCTION FOR ABSOLUTE VALUES IN POLYNOMIAL RINGS. than define a second approximation V 0

Transcription:

CUBO A Mathematical Journal Vol.17, N ō 3, (71 9). October 15 Computing the Laplace transform and the convolution for more functions adjoined Takahiro Sudo Department of Mathematical Sciences, Faculty of Science, University of the Ryukyus, Senbaru 1, Nishihara, Okinawa 93-13, Japan. sudo@math.u-ryukyu.ac.jp ABSTRACT We compute explicitly the Laplace transform and the convolution for more functions on the real than the continuous functions that obtained as the inverse Laplace transforms of rational functions on the complex plane vanishing at infinity. We also consider the algebraic structure with more functions adjoined. RESUMEN Calculamos explícitamente la transformada de Laplace y la convolución para más funciones en los reales que las funciones continuas obtenidas como transformadas de Laplace inversas de funciones racionales en el plano complejo que se anulan en el infinito. También consideramos la estructura algebraica con más funciones adjuntadas. Keywords and Phrases: Laplace tranform, Gaussian function, Dirac delta function, convolution. 1 AMS Mathematics Subject Classification: 44A1, 44A35, 6C15, 6A6, 6A9.

7 Takahiro Sudo CUBO 1 Introduction 17, 3 (15) This paper is a continuation from the paper [3] of the author, in which we compute explicitly the inverse Laplace transform for rational functions on the complex plane vanishing at infinity and consider the algebraic structure for the corresponding algebra of continuous elementary functions on the real R for which the Laplace transform on the infinite interval [, ) is defined. In this paper, we compute explicitly the Laplace transform for more functions on the real than the continuous functions that obtained as the inverse Laplace transforms of rational functions on the complex plane vanishing at infinity. For this, we compute the convolution product with more functions adjoined, such as Gaussian functions, the Dirac delta functions, and the shifted delta functions. We also consider the algebraic structure with more functions adjoined. This paper after this introduction is organized as follows: Adjoining Gaussian functions; 3 Algebraic structure with Gaussian functions; 4 Adjoining the delta function; 5 Algebraic structure with the delta function; 6 Adjoining the shifted delta functions; 7 Algebraic structure with the shifted delta functions; 8 More with the weak derivative. Our elementary explicit computation results obtained and algebraic consideration results obtained should be useful as a reference. Refer to [4] or [] (or [1]) for some basics about the Laplace transform. Notation. We denote by A(R) the real algebra with convolution, generated by the three sets of elementary continuous functions on the real line R with t real variable: {t n n N = {, 1,, }} of monomials and {e µt µ R} of exponential functions based to e and {sin t, cos t R} of trigonometric functions and their point-wise multiples such as t n e µt sin t. The (extended in our sense) convolution product for A(R) (but the same as the usual one by restriction to [, )) is defined as the (Riemann) integral: f g =(f g)(t) = f(t τ)g(τ)dτ, (t R) for f, g A(R), which makes sense and is a commutative and associative operation as well known as the usual case on the interval [, ). Denote by C the complex plane. For f = f(t) a suitable function on R (as in A(R)), the (extended in our sense) Laplace transform of f (but the same as the usual one of f restricted to [, )) is defined by the (Riemann) integral: L(f(t)) = L(f(t))(s) = e st f(t)dt for s C (formally or in the domain of convergence).

CUBO 17, 3 (15) Computing the Laplace transform and the convolution for more... 73 Adjoining Gaussian functions Quite well known, but first of all we consider Gaussian functions as follows: Lemma.1. For t R, s C and R with >,wehave π L(e t )= e s 4, Proof. Check that e st e t = e (t+ s ) + s 4 and then L(e t )=e s 4 e (t+ s ) dt. Using change of variables and Cauchy s theorem, we obtain that the integral is computed to be equal to s + e u s + = 1 du with u = (t + s )) π e u du =. Lemma.. For t R, s C and, µ R positive, we have L(e t e µt )= π 4 µ e +µ 4µ s. We need to check the convolution explicitly as follows. Proposition 1. For t R and positive, ν R, 1 (e t e µt )(t) = e (1 + µ = e (1 +µ )t n= µt +µ +µ )t t +µ e u du n+1 + µ n+1 ( 1) n 1 ( + µ) n+1 n! n + 1 tn+1 Proof. We compute (e t e µt )(t) = e t t (+µ)(τ e +µ ) + t e (t τ) e µτ dτ = +µ dτ = e (1 +µ )t t (+µ)(τ e +µ ) dτ

74 Takahiro Sudo CUBO 17, 3 (15) with 1 +µ >, and the integral is computed by change of variables as: g(t) = t (+µ)(τ e +µ ) dτ µt 1 +µ + µ t +µ e u du with u = + µ(τ t + µ ). Using Taylor s expansion of e y we compute the integral above as: µt +µ t +µ ( 1) n n= µt e u +µ du = t +µ n! µt +µ ( 1) n u n du = n! n= u n du = t +µ ( 1) n µ n+1 + n+1 1 n! ( + µ) n+ 1 n + 1 tn+1. n= Corollary 1. The convolutions e t e µt for, µ positive reals are point-wise limits of elements of A(R). Proof. Note that any Gaussian functions are point-wise limits of elements of A(R) by Taylor series. Proposition. For R positive, n N, andt R, (1) (t n e t )= with I k τk e τ dτ given by and I k+1 = e t k l=1 I k = e t + n k= ( ) n t t n k ( 1) k τ k e τ dτ k (k)!! () l ((k l + 1))!! t(k l+1) + (k)!! () k+1 (1 ), e t k l=1 (k 1)!! () k (k 1)!! () l ((k l + 1) 1)!! t(k l+1) 1 m= Also, for µ R non-zero and R positive, () (e µt e t )=e µ 4 e µt n= ( ) m m!(m + 1) tm+1. ( 1) n n+ 1 n!(n + 1) {(t + µ )n+1 ( µ )n+1 }.

CUBO 17, 3 (15) Computing the Laplace transform and the convolution for more... 75 Also, for µ R non-zero and R positive, with (3) (sin µt e t )=sin µt cos µτe τ dτ = sin µτe τ dτ = cos µτe τ dτ cos µt ( ) n n= n= n! ( ) n n! τ n cos µτdτ, τ n sin µτdτ sin µτe τ dτ and each integral term is computed inductively as: and similarly, I c,n I s,n τ n cos µτdτ = 1 µ tn sin µt + n µ tn 1 cos µτ τ n sin µτdτ = 1 µ tn cos µt + n µ tn 1 sin µτ (4) (cos µt e t )=cos µt with the integrals the same as above. cos µτe τ dτ sin µt (n)(n 1) µ I c,n, (n)(n 1) µ I s,n, sin µτe τ dτ Proof. For (1) we compute (t n e t )(t) = (t τ) n e τ dτ = n k= ( ) n t t n k ( 1) k τ k e τ dτ k and the integral is computed by integration by parts as: { } 1 I k τ k e τ dτ = τ k 1 ) dτ (e τ = 1 tk 1 e t + k 1 I k

76 Takahiro Sudo CUBO 17, 3 (15) so that, inductively, I k+1 = 1 tk e t k(k ) 4 () k t e t + k k(k ) () tk e t () 3 t k 4 e t k(k ) () k I 1 with and with I 1 τe τ dτ = 1 (1 e t ), I k = 1 tk 1 e t k 1 () tk 3 e t (k 1)(k 3) 3 () k te t + I = m= e τ dτ = ( ) m m! m= (k 1)(k 3) () 3 (k 1)(k 3) 1 () k I ( ) m τ m dτ m! τ m dτ = m= ( ) m m!(m + 1) tm+1. t k 5 e t For () we compute e µt e t = e µ(t τ) e τ dτ = e µ 4 e µt e (τ+ µ ) dτ. By change of variables as u = (τ + µ ), the integral above is computed as = e (τ+ µ ) dτ = ( 1) n n= n! (t+ µ ) ( µ ) (t+ µ ) ( µ ) u n du = e u (t+ µ ) du = n= ( µ ) n= 1 n! ( u ) n du ( 1) n n+ 1 n!(n + 1) {(t + µ )n+1 ( µ )n+1 }. For (3) we next compute (sin µt e t )(t) = = sin µt sin µ(t τ)e t dτ cos µτe τ dτ cos µt sin µτe τ dτ,

CUBO 17, 3 (15) Computing the Laplace transform and the convolution for more... 77 and the first integral is computed as: cos µτe τ dτ = cos µτ n= and each integral term is computed inductively as: I c,n τ n cos µτdτ 1 n! ( τ ) n dτ = = 1 µ tn sin µt + n µ tn 1 cos µτ by using integration by parts. Similarly, we obtain and I s,n sin µτe τ dτ = τ n sin µτdτ ( ) n n= n! = 1 µ tn cos µt + n µ tn 1 sin µτ ( ) n n= n! (n)(n 1) µ I c,n τ n sin µτdτ τ n cos µτdτ (n)(n 1) µ I s,n. Remark. We omit the case of convolutions with general monomials such as multiples t n e µt sin µt, but which are certainly computable by the similar computation technique as above, with some more integration. Refer to [3]. 3 Algebraic structure with Gaussian functions We denote by A w (R) the set of all functions on R that are written as point-wise limits in R of elements of A(R). Lemma 3.1. A w (R) is an algebra over R with point-wise multiplication. algebra over R with convolution. Also, A w (R) is an Proof. Suppose that f, g A w (R) with f = lim f n, g = lim g n (n ) point-wise limits of f n,g n A w (R). Since the limits at any t R are finite, then and hence, f g A w (R). (f g)(t) =lim f n (t) lim g n (t) =lim(f n g n )(t) R,

78 Takahiro Sudo CUBO 17, 3 (15) Note also that we define the convolution of f, g A w (R) as f g = lim(f n g n ) A w (R), where it is shown in general ([3]) that A(R) is closed under convolution. This is well defined, because if f = lim h n and g = lim k n another limits with h n,k n A(R), and since the limits at any t R are finite, then (f n g n )(t) (h n k n )(t) =((f n f) g n )(t)+(f (g n g))(t) +((f h n ) g)(t)+(h n (g k n ))(t), and the right hand side goes to zero as n by applying the dominated convergence theorem to each convolution, and hence, lim(f n g n ) lim(h n k n ) lim(f n g n ) f n g n + f n g n h n k n + h n k n lim(h n k n ), which go to zero as n. We denote by G(R) the algebra generated by A(R) and the set {e t R, > }, with either point-wise multiplication or convolution as a product. Corollary. The algebra G(R) is contained in A w (R) and is a subalgebra of A w (R) with either point-wise multiplication or convolution. We denote by R (C) the algebra of rational functions on C vanishing at infinity with point-wise multiplication and by E(C) the algebra generated by both R (C) and the set {e s R, > } of inverse Gaussian functions on C with s complex variable. Corollary 3. The algebra G(R) with convolution is isomorphic to the algebra E(C) by the Laplace transform and by the inverse Laplace transform. Proof. Note that L(f g) =L(f) L(g) and L is injective for continuous functions on R. 4 Adjoining the delta function We consider the Dirac delta function. Define the delta function as δ(t) = lim g (t) lim π e t = { t = t with R positive. Note that the delta function has the meaning in integration only in the weak sense, i.e. as the weak limit of g (t) by taking the inner product or the convolution with test

CUBO 17, 3 (15) Computing the Laplace transform and the convolution for more... 79 functions, and the delta function is also viewed as a distribution, i.e., a functional (or a scalarvalued linear map) on the linear space of test functions. Also, one can take other functions of different types to define the delta function. Note that g (t)dt = but we need to have the next: Lemma 4.1. The functions g are in G(R), andfors C, which goes to 1 as. L(g )(s) =e s 4, Proof. Use Lemma.1. Lemma 4.. For s C, wehave which goes to 1 as and µ. L(g g µ )(s) =e (+µ)s 4µ, Proof. Use Lemma 4.1 and L(g g µ )=L(g ) L(g µ ). Proposition 3. For positive µ, R, (1) lim (e µt g )(t) = e µt if t>, if t =, e µt if t<; and for n N and >, and for µ R and >, () lim (tn g )(t) = (3) lim (eµt g )(t) = t n if t>, if t =, t n if t<; e µt if t>, if t =, if t<; e µt and for µ R and >, (4) lim (sin µt g )(t) = sin µt if t>, if t =, sin µt if t<;

8 Takahiro Sudo CUBO 17, 3 (15) and (5) lim (cos µt g )(t) = cos µt if t>, if t =, cos µt if t<. Proof. The first equation (1) follows from the first equation in Proposition 1 by taking the limit with respect to. For () we next use Proposition. It follows that lim (tn g )(t) = lim tn g (τ)dτ, namely, other terms in the convolution go to zero as,andthen which goes to g (τ)dτ = t 1 if t>, if t =, 1 if t<, π e u du as. (u = τ), For (3), it also follows as in parts of the proof of Proposition that (e µt g )(t) =e µ 4 e µt µ ) dτ = e µ 4 e µt π e (τ+ (t+ µ ) µ π e u du (u = (τ + µ )), which goes to e µt if t>, if t =, if t<, e µt as. For (4) we next compute by allowing complex coefficients for this case and using the Euler formula as: (sin µt g )(t) =( eµti e µti g )(t) i = 1 e µ(t τ)i g (τ)dτ 1 e µ(t τ)i g (τ)dτ i i = eµti i e µτi g (τ)dτ e µti i e µτi g (τ)dτ

CUBO 17, 3 (15) Computing the Laplace transform and the convolution for more... 81 and the first integral is converted to which goes to = e µ 4 µi π e (τ+ ) µ 4 dτ (t+ µi ) iµ e u du (u = (τ + µi π )) 1 if t>, if t =, 1 if t< as via Cauchy integral theorem. Similarly, the second integral above is converted to 4 dτ which goes to = e µ 4 µi π e (τ ) µ (t µi ) iµ e u du (u = (τ µi π )) 1 if t>, if t =, 1 if t< as. Therefore, we obtain one of the last two equations in the statement. The other case (5) is obtained similarly, by cos µt = eµti +e µti. 5 Algebraic structure with the delta function We denote by G δ (R) the algebra generated by G(R) and the delta function δ(t) with convolution, where we assume that f δ = f on [, ) for any f G(R) and δ δ = δ on [, ), which follows from our definition of the delta function in the weak sense, so that G δ (R) restricted to [, ) is isomorphic to the unitization of G(R) by R, restricted to [, ), as a real algebra. We define the Laplace transform for δ as the unit function on C. L(δ(t))(s) =1(s) Note that G δ (R) is also viewed as the generated by G(R) and the point-wise limit function of Gaussian functions on R. We denote by E 1 (C) the algebra generated by E(C) and the unit function 1 on C with pointwise multiplication.

8 Takahiro Sudo CUBO 17, 3 (15) Corollary 4. The algebra G δ (R) with convolution, restricted to [, ), is isomorphic to the algebra E 1 (C) by the Laplace transform and by the inverse Laplace transform. Remark. Our Proposition 3 says that the convolution f δ on R for f G δ (R) is defined as f(t) if t>, (f δ)(t) = if t =, f(t) if t<. 6 Adjoining the shifted delta functions It is also defined as that for t R and s C, L(δ(t µ)) = e µs for µ>. We extend this definition for µ R. Indeed,checkthat Lemma 6.1. We have for any µ R. lim L(g (t µ))(s) =e µs Proof. We compute L(g (t µ))(s) = e st π e (t µ) dt = e µs e s 4 e (t+ s µ) π dt and changing variables as u = (t + s u) and using Cauchy integral theorem and taking the limit as we obtain the result in the statement. Lemma 6.. We have lim L(g (t µ) g (t ρ))(s) =e (µ+ρ)s for any µ, ρ R. It follows by defining the left hand side to be equal to Proof. Note that L(δ(t µ) δ(t ρ))(s) =e (µ+ρ)s. L(g (t µ) g (t ρ))(s) =L(g (t µ))(s) L(g (t ρ))(s) and the both factors are computed in the proof of Lemma 6.1.

CUBO 17, 3 (15) Computing the Laplace transform and the convolution for more... 83 Proposition 4. For positive µ, R and ρ R, if ρ>,then e µ(t ρ) t>ρ, (1) + lim (e µt g (t ρ))(t) = 1 t = ρ, t < ρ, and if ρ<,then (1) lim (e µt g (t ρ))(t) = t > ρ, 1 t = ρ, t<ρ; e µ(t ρ) and for n N and ρ R, if ρ>,then and if ρ<,then and for µ R, if ρ>,then () + lim (tn g (t ρ))(t) = () lim (tn g (t ρ))(t) = (3) + lim (eµt g (t ρ))(t) = t n t>ρ, ρ n t = ρ, t < ρ, t > ρ, ρ n t = ρ, t n t<ρ; e µ(t ρ) if t>ρ, 1 if t = ρ, if t<ρ, and if ρ<,then (3) lim (eµt g (t ρ))(t) = if t>ρ, 1 if t = ρ, if t<ρ; e µ(t ρ) and for µ R, if ρ>,then and if ρ<,then (4) + lim (sin µt g (t ρ)(t) = (4) lim (sin µt g (t ρ)(t) = sin µ(t ρ) if t>ρ, if t = ρ, if t<ρ, if t>ρ, if t = ρ, sin µ(t ρ) if t<ρ;

84 Takahiro Sudo CUBO 17, 3 (15) and if ρ>,then and if ρ>,then (5) + lim (cos µt g (t ρ))(t) = (5) lim (cos µt g (t ρ))(t) = cos µ(t ρ) if t>ρ, 1 if t = ρ, if t<ρ, if t>ρ, 1 if t = ρ, cos µ(t ρ) if t<ρ. Proof. For (1) ± we compute (e µt g (t ρ))(t) = e π e (t τ ρ) µτ dτ = e µ +µ (t ρ) π e (+µ)(τ +µ(t (t ρ) = e µ +µ (t ρ) +µ ) (t ρ) +µ (t ρ) +µ ) dτ π( + µ) e u du where u = + µ(τ (t ρ) µt+ρ +µ ) and if τ = t, thenu =, which goes to as for ρ>. +µ If t>ρ,then (t ρ) +µ goes to and if t<ρ, then it goes to. Therefore, if ρ>,then lim g (e µt (t ρ))(t) = e µ(t ρ) t>ρ, 1 t = ρ, t < ρ, and if ρ<,then lim (e µt g (t ρ))(t) = t > ρ, 1 t = ρ, t<ρ, e µ(t ρ) For () ± we next compute (t n g (t ρ))(t) = n k= and the following integrals are converted as: τ k e (τ ρ) dτ = ( ) n t t n k ( 1) k τ k dρ k π e (τ ρ) (t ρ) ρ (ρ + u ) k e u du

CUBO 17, 3 (15) Computing the Laplace transform and the convolution for more... 85 where u = (τ ρ). We then use Proposition and it follows that if ρ>,then t n t>ρ, lim (tn g (t ρ))(t) = ρ n t = ρ, t < ρ, and if ρ<,then lim (tn g (t ρ))(t) = t > ρ, ρ n t = ρ, t n t<ρ. For (3) ± we next compute (e µt g (t ρ))(t) = e µ(t τ) dτ π e (τ ρ) = e µ 4ρµ 4 e µt π e (τ+ = e µ 4ρµ 4 e µt (t ρ+ µ ) µ ρ µ ρ ) dτ π e u du (u = (τ + µ ρ )), which goes to, if ρ>,then as, and if ρ<,then e µ(t ρ) if t>ρ, 1 if t = ρ, if t<ρ if t>ρ, 1 if t = ρ, if t<ρ. e µ(t ρ) For (4) ± we finally compute by allowing complex coefficients for this case and using the Euler formula as: (sin µt g (t ρ))(t) =( eµti e µti g (t ρ))(t) i = 1 e µ(t τ)i g (τ ρ)dτ 1 e µ(t τ)i g (τ ρ)dτ i i = eµti i e µτi g (τ ρ)dτ e µti i e µτi g (τ ρ)dτ

86 Takahiro Sudo CUBO 17, 3 (15) and the first integral is converted to µ ρµi e 4 µ ρµi = e 4 which goes to, if ρ>,then, µi π e (τ ρ+ (t ρ+ µi ) ( ρ+ µi ) ) dτ π e u du (u = (τ ρ + µi )) e ρµi if t>ρ, e ρµi if t = ρ, if t<ρ as via Cauchy integral theorem, and if ρ<,then, if t>ρ, e ρµi if t = ρ, e ρµi if t<ρ. Similarly, the second integral above is converted to have the similar limits as above, where e ρµi is replaced with e ρµi since u is replaced with u = (τ ρ µi ). Therefore, we obtain one of the last two equations in the statement. The other case (5) ± is obtained similarly, by cos µt = e µti +e µti. Remark. The Heaviside unit function u(t) on R is defined by u(t) =1 if t>and u(t) = if t<aswell as u() = 1. For ρ R positive, the shifted Heaviside unit function u(t ρ) is defined similarly. It is known that L(u(t ρ))(s) = e ρs. s We may extend the definition of u(t ρ) for ρ<. But we always have Note also that L 1 ( 1 s e ρs )= L(u(t ρ))(s) = 1 s, ρ. 1 δ(τ ρ)dτ = { 1 t > ρ, t < ρ, where note that this shifted delta function comes from the usual Gaussian functions, not equal to ours multiplied by.

CUBO 17, 3 (15) Computing the Laplace transform and the convolution for more... 87 7 Algebraic structure with the shifted delta functions We denote by G sδ (R) the algebra generated by G(R) and the shifted delta functions δ(t µ) for µ R, with convolution. For convenience, we define G sδ (R + ) to be the algebra generated by G δ (R) and the positive shifted delta functions δ(t µ) for µ>, restricted to [, ). We define the following convolution as (f δ(t µ))(t) =u(t µ)f(t µ) for f G sδ (R + ) and µ>,whereu(t µ) is the shifted Heaviside unit function on [, ). Note that this definition is not perfectly compatible with Proposition 4 but either changing g (t µ) with it multiplied by 1 on [, ) or cutting off functions to be zero on (,) make sense almost everywhere. In fact, such a scalar multiplication is allowed to define the convolution with the shifted delta functions in the weak sense. We denote by E s1 (C) the algebra generated by E 1 (C) and the set {e ρs ρ>}, with pointwise mulitiplication. We define the Laplace transform L(f g) to be L(f) L(g) for f, g G sδ (R + ). Note that it is known as the first shift law that and L(δ(t µ)) = e µs. L(u(t µ)f(t µ)) = e µs L(f)(s) Corollary 5. The algebra G sδ (R + ) with convolution is isomorphic to E s1 (C) with pointwise multiplication, as a real algebra, via the Laplace transform and the inverse Laplace transform. Proof. The Laplace transform L is injective on G δ (R). The injectiveness of L extends to G sδ (R) by definition. Also, G sδ (R + ) is mapped onto E s1 (C) by L as an algebra homomorphism, by definition. Note that, by definition, G sδ (R + ) contains discontinuous functions such as the shifted Heaviside unit functions. Lemma 7.1. For f, g G sδ (R + ) and for µ, s >, wehave (u(t µ)f(t µ) u(t s)g(t s))(t) { if t<µor t<s, = u f(t τ µ)g(τ s)dτ if t>µand t>s. s Note that in the second non-trivial case, the integral is just the sub-integral of (f(t µ) g(t s))(t) = f(t τ µ)g(τ s)dτ restricted to the sub-interval [s, t µ] of [, t], sothatwemay call the convolution in the statement the sub-convolution. Therefore,

88 Takahiro Sudo CUBO 17, 3 (15) Corollary 6. The convolution in G sδ (R + ) generates such sub-convolutions, which generate the algebra. Remark. However, the convolution product in G sδ (R + ) is somewhat difficult to see what it is, but the corresponding point-wise multiplication in E s1 (C) is relatively easy to see what it is. As well known in applications to solving ordinary differential equations with initial values, a problem related to G sδ (R + ) is relatively easily solved in E s1 (C) via the Laplace transform, and then the last task is to compute the corresponding convolution via the inverse Laplace transform. 8 More with the weak derivative Recall that Lemma 8.1. We have L(δ (t)) = s in the weak sense, i.e., as a functional or a distribution, so that we may definel 1 (s) =δ (t) for t R and s C. Proof. We compute L(δ (t)) = = e st δ (t)dt ( 1) t e st δ(t)dt = sl(δ(t)) = s. by definition Note that the algebra G(R) is closed under the usual derivative and that G(R) is contained in C (R) the algebra of all infinitely many times differentiable (or smooth) functions on R. Lemma 8.. For f G(R), wehave in the weak sense. Proof. We compute (f δ )(t) = = (f δ )(t) =f (t) f(t τ)δ (τ)dτ ( 1) τ f(t τ) δ(τ)dτ =(f δ)(t) =f (t).

CUBO 17, 3 (15) Computing the Laplace transform and the convolution for more... 89 We define the n-th derivative δ (n) to be the n-times convolution with δ. For f G(R), we define f δ (n) to be f (n) the n-th derivative of f. We also define L(δ (n) ) to be s n for s C. We denote by D A (R) the algebra generated by A(R), the delta function δ, and the set {δ (n) n Z,n 1}, with convolution. We also denote by R(C) the algebra of all rational functions on C with point-wise multiplication. Proposition 5. There is an algebra isomorphism by the Laplace transform from D A (R) to R(C). We denote by D G (R) the algebra defined by replacing A(R) with G(R) in the definition of D A (R). We also denote by R E (C) the algebra generated by R(C) and E(C). Proposition 6. There is an algebra isomorphism by the Laplace transform from D G (R) to R E (C). Lemma 8.3. For f G sδ (R + ) and µ>,wehave almost everywhere in the weak sense. (f δ (t µ))(t) =u(t µ)f (t) Proof. We compute (f δ (t µ))(t) = = f(t τ)δ (τ µ)dτ ( 1) τ f(t τ) δ(τ µ)dτ =(f δ(t µ))(t) = u(t µ)f (t) for f G sδ (R + ),wherenotethatf is defined almost everywhere (except some finite points of R + ). We denote by D G sδ(r + ) the algebra generated by G sδ (R + ) and the weak derivatives δ (n) (t µ) for n positive integers and µ>reals. We also denote by R E s1(c) the algebra generated by R(C) and E s1 (C). Proposition 7. There is an algebra isomorphism by the Laplace transform from D G sδ(r + ) to R E s1(c). Received: May 13. Accepted: May 13.

9 Takahiro Sudo CUBO 17, 3 (15) References [1] Hiroshi Fukawa, Laplace transformation and ordinary differential equations (in Japanese), Sho-Ko-dou (1995). [] M. S. J., Mathematics Dictionary, (Sugaku Jiten, in Japanese), Math. Soc. Japan, 4th edition, Iwanami (7). [3] T. Sudo, Computing the inverse Laplace transform for rational functions vanishing at infinity, Cubo A Math. J., 16, No 3 (14), 97-117. [4] A. Vretblad, Fourier Analysis and Its Applications, GTM 3, Springer (3).