5 Fourier and Laplace Transforms

Total Page:16

File Type:pdf, Size:1020Kb

5 Fourier and Laplace Transforms 5 Fourier and Laplace Transforms “There is no branch of mathematics, however abstract, which may not some day be applied to phenomena of the real world.”, Nikolai Lobatchevsky (1792-1856) 5.1 Introduction In this chapter we turn to the study of Fourier transforms, which provide integral representations of functions defined on the entire real line. Such functions can represent analog signals. Recall that analog signals are continuous signals which are sums over a continuous set of fre- quencies. Our starting point will be to rewrite Fourier trigonometric series as Fourier exponential series. The sums over discrete frequencies will lead to a sum (integral) over continuous frequencies. The resulting integrals will be complex integrals, which can be evaluated using contour methods. We will investigate the properties of these Fourier transforms and get prepared to ask how the analog signal representations are related to the Fourier se- ries expansions over discrete frequencies which we had seen in Chapter 2. Fourier series represented functions which were defined over finite do- mains such as x 2 [0, L]. Our explorations will lead us into a discussion of the sampling of signals in the next chapter. We will also discuss a related integral transform, the Laplace transform. In this chapter we will explore the use Laplace transforms are useful in solving initial value problems in differen- of integral transforms. Given a function f (x), we define an integral transform to tial equations and can be used to relate the input to the output of a linear a new function F(k) as system. Both transforms provide an introduction to a more general theory Z b F(k) = f (x)K(x, k) dx. of transforms, which are used to transform specific problems to simpler a ones. Here K(x, k) is called the kernel of the In Figure 5.1 we summarize the transform scheme for solving an initial transform. We will concentrate specifi- cally on Fourier transforms, value problem. One can solve the differential equation directly, evolving the Z ¥ initial condition y(0) into the solution y(t) at a later time. fˆ(k) = f (x)eikx dx, −¥ However, the transform method can be used to solve the problem indi- and Laplace transforms rectly. Starting with the differential equation and an initial condition, one Z ¥ F(s) = f (t)e−st dt. computes its Transform (T) using 0 Z ¥ Y(s) = y(t)e−st dt. 0 170 fourier and complex analysis Figure 5.1: Schematic of using trans- ODE, y(0) T Alg Eqn forms to solve a linear ordinary differ- ential equation. solve solve y(t) IT Y(s) Applying the transform to the differential equation, one obtains a simpler (algebraic) equation satisfied by Y(s), which is simpler to solve than the original differential equation. Once Y(s) has been found, then one applies the Inverse Transform (IT) to Y(s) in order to get the desired solution, y(t). We will see how all of this plays out by the end of the chapter. We will begin by introducing the Fourier transform. First, we need to see how one can rewrite a trigonometric Fourier series as complex exponential series. Then we can extend the new representation of such series to ana- log signals, which typically have infinite periods. In later chapters we will highlight the connection between these analog signals and their associated digital signals. 5.2 Complex Exponential Fourier Series Before deriving the Fourier transform, we will need to rewrite the trigonometric Fourier series representation as a complex exponential Fourier series. We first recall from Chapter 2 the trigonometric Fourier se- ries representation of a function defined on [−p, p] with period 2p. The Fourier series is given by ¥ a0 f (x) ∼ + ∑ (an cos nx + bn sin nx) ,( 5.1) 2 n=1 where the Fourier coefficients were found as 1 Z p an = f (x) cos nx dx, n = 0, 1, . , p −p 1 Z p bn = f (x) sin nx dx, n = 1, 2, . .( 5.2) p −p In order to derive the exponential Fourier series, we replace the trigono- metric functions with exponential functions and collect like exponential terms. This gives ¥ inx −inx inx −inx a0 e + e e − e f (x) ∼ + ∑ an + bn 2 n=1 2 2i a ¥ a − ib ¥ a + ib = 0 + ∑ n n einx + ∑ n n e−inx.( 5.3) 2 n=1 2 n=1 2 fourier and laplace transforms 171 The coefficients of the complex exponentials can be rewritten by defining 1 c = (a + ib ), n = 1, 2, . .( 5.4) n 2 n n This implies that 1 c¯ = (a − ib ), n = 1, 2, . .( 5.5) n 2 n n So far, the representation is rewritten as ¥ ¥ a0 inx −inx f (x) ∼ + ∑ c¯ne + ∑ cne . 2 n=1 n=1 Re-indexing the first sum, by introducing k = −n, we can write −¥ ¥ a0 −ikx −inx f (x) ∼ + ∑ c¯−ke + ∑ cne . 2 k=−1 n=1 Since k is a dummy index, we replace it with a new n as −¥ ¥ a0 −inx −inx f (x) ∼ + ∑ c¯−ne + ∑ cne . 2 n=−1 n=1 We can now combine all the terms into a simple sum. We first define cn for negative n’s by cn = c¯−n, n = −1, −2, . a0 Letting c0 = 2 , we can write the complex exponential Fourier series repre- sentation as ¥ −inx f (x) ∼ ∑ cne ,( 5.6) n=−¥ where 1 c = (a + ib ), n = 1, 2, . , n 2 n n 1 c = (a− − ib− ), n = −1, −2, . , n 2 n n a c = 0 .( 5.7) 0 2 Given such a representation, we would like to write out the integral forms of the coefficients, cn. So, we replace the an’s and bn’s with their integral representations and replace the trigonometric functions with complex expo- nential functions. Doing this, we have for n = 1, 2, . , 1 c = (a + ib ) n 2 n n 1 1 Z p i Z p = f (x) cos nx dx + f (x) sin nx dx 2 p −p p −p 1 Z p einx + e−inx i Z p einx − e−inx = f (x) dx + f (x) dx 2p −p 2 2p −p 2i 1 Z p = f (x)einx dx.( 5.8) 2p −p 172 fourier and complex analysis It is a simple matter to determine the cn’s for other values of n. For n = 0, we have that Z p a0 1 c0 = = f (x) dx. 2 2p −p For n = −1, −2, . ., we find that Z p Z p 1 −inx 1 inx cn = c¯n = f (x)e dx = f (x)e dx. 2p −p 2p −p Therefore, we have obtained the complex exponential Fourier series coeffi- cients for all n. Now we can define the complex exponential Fourier series for the function f (x) defined on [−p, p] as shown below. Complex Exponential Series for f (x) Defined on [−p, p] ¥ −inx f (x) ∼ ∑ cne ,( 5.9) n=−¥ Z p 1 inx cn = f (x)e dx.( 5.10) 2p −p We can easily extend the above analysis to other intervals. For example, for x 2 [−L, L] the Fourier trigonometric series is ¥ a0 npx npx f (x) ∼ + ∑ an cos + bn sin 2 n=1 L L with Fourier coefficients 1 Z L npx an = f (x) cos dx, n = 0, 1, . , L −L L 1 Z L npx bn = f (x) sin dx, n = 1, 2, . L −L L This can be rewritten as an exponential Fourier series of the form Complex Exponential Series for f (x) Defined on [−L, L] ¥ −inpx/L f (x) ∼ ∑ cne ,( 5.11) n=−¥ Z L 1 inpx/L cn = f (x)e dx.( 5.12) 2L −L We can now use this complex exponential Fourier series for function de- fined on [−L, L] to derive the Fourier transform by letting L get large. This will lead to a sum over a continuous set of frequencies, as opposed to the sum over discrete frequencies, which Fourier series represent. 5.3 Exponential Fourier Transform Both the trigonometric and complex exponential Fourier series provide us with representations of a class of functions of finite period in fourier and laplace transforms 173 terms of sums over a discrete set of frequencies. In particular, for functions defined on x 2 [−L, L], the period of the Fourier series representation is 2L. We can write the arguments in the exponentials, e−inpx/L, in terms of −iw x the angular frequency, wn = np/L, as e n . We note that the frequencies, np nn, are then defined through wn = 2pnn = L . Therefore, the complex exponential series is seen to be a sum over a discrete, or countable, set of frequencies. We would now like to extend the finite interval to an infinite interval, x 2 (−¥, ¥), and to extend the discrete set of (angular) frequencies to a continuous range of frequencies, w 2 (−¥, ¥). One can do this rigorously. n It amounts to letting L and n get large and keeping L fixed. p We first define Dw = L , so that wn = nDw. Inserting the Fourier coeffi- cients (5.12) into Equation (5.11), we have ¥ −inpx/L f (x) ∼ ∑ cne n=−¥ ¥ 1 Z L = ∑ f (x)einpx/L dx e−inpx/L n=−¥ 2L −L ¥ Dw Z L = ∑ f (x)eiwnx dx e−iwn x.( 5.13) n=−¥ 2p −L Now, we let L get large, so that Dw becomes small and wn approaches the angular frequency w. Then, 1 ¥ Z L f (x) ∼ lim ∑ f (x)eiwnx dx e−iwn xDw Dw!0,L!¥ 2p n=−¥ −L 1 Z ¥ Z ¥ = f (x)eiwx dx e−iwx dw.( 5.14) 2p −¥ −¥ Looking at this last result, we formally arrive at the definition of the Definitions of the Fourier transform and Fourier transform.
Recommended publications
  • The Laplace Transform of the Psi Function
    PROCEEDINGS OF THE AMERICAN MATHEMATICAL SOCIETY Volume 138, Number 2, February 2010, Pages 593–603 S 0002-9939(09)10157-0 Article electronically published on September 25, 2009 THE LAPLACE TRANSFORM OF THE PSI FUNCTION ATUL DIXIT (Communicated by Peter A. Clarkson) Abstract. An expression for the Laplace transform of the psi function ∞ L(a):= e−atψ(t +1)dt 0 is derived using two different methods. It is then applied to evaluate the definite integral 4 ∞ x2 dx M(a)= , 2 2 −a π 0 x +ln (2e cos x) for a>ln 2 and to resolve a conjecture posed by Olivier Oloa. 1. Introduction Let ψ(x) denote the logarithmic derivative of the gamma function Γ(x), i.e., Γ (x) (1.1) ψ(x)= . Γ(x) The psi function has been studied extensively and still continues to receive attention from many mathematicians. Many of its properties are listed in [6, pp. 952–955]. Surprisingly, an explicit formula for the Laplace transform of the psi function, i.e., ∞ (1.2) L(a):= e−atψ(t +1)dt, 0 is absent from the literature. Recently in [5], the nature of the Laplace trans- form was studied by demonstrating the relationship between L(a) and the Glasser- Manna-Oloa integral 4 ∞ x2 dx (1.3) M(a):= , 2 2 −a π 0 x +ln (2e cos x) namely, that, for a>ln 2, γ (1.4) M(a)=L(a)+ , a where γ is the Euler’s constant. In [1], T. Amdeberhan, O. Espinosa and V. H. Moll obtained certain analytic expressions for M(a) in the complementary range 0 <a≤ ln 2.
    [Show full text]
  • 7 Laplace Transform
    7 Laplace transform The Laplace transform is a generalised Fourier transform that can handle a larger class of signals. Instead of a real-valued frequency variable ω indexing the exponential component ejωt it uses a complex-valued variable s and the generalised exponential est. If all signals of interest are right-sided (zero for negative t) then a unilateral variant can be defined that is simple to use in practice. 7.1 Development The Fourier transform of a signal x(t) exists if it is absolutely integrable: ∞ x(t) dt < . | | ∞ Z−∞ While it’s possible that the transform might exist even if this condition isn’t satisfied, there are a whole class of signals of interest that do not have a Fourier transform. We still need to be able to work with them. Consider the signal x(t)= e2tu(t). For positive values of t this signal grows exponentially without bound, and the Fourier integral does not converge. However, we observe that the modified signal σt φ(t)= x(t)e− does have a Fourier transform if we choose σ > 2. Thus φ(t) can be expressed in terms of frequency components ejωt for <ω< . −∞ ∞ The bilateral Laplace transform of a signal x(t) is defined to be ∞ st X(s)= x(t)e− dt, Z−∞ where s is a complex variable. The set of values of s for which this transform exists is called the region of convergence, or ROC. Suppose the imaginary axis s = jω lies in the ROC. The values of the Laplace transform along this line are ∞ jωt X(jω)= x(t)e− dt, Z−∞ which are precisely the values of the Fourier transform.
    [Show full text]
  • FOURIER TRANSFORM Very Broadly Speaking, the Fourier Transform Is a Systematic Way to Decompose “Generic” Functions Into
    FOURIER TRANSFORM TERENCE TAO Very broadly speaking, the Fourier transform is a systematic way to decompose “generic” functions into a superposition of “symmetric” functions. These symmetric functions are usually quite explicit (such as a trigonometric function sin(nx) or cos(nx)), and are often associated with physical concepts such as frequency or energy. What “symmetric” means here will be left vague, but it will usually be associated with some sort of group G, which is usually (though not always) abelian. Indeed, the Fourier transform is a fundamental tool in the study of groups (and more precisely in the representation theory of groups, which roughly speaking describes how a group can define a notion of symmetry). The Fourier transform is also related to topics in linear algebra, such as the representation of a vector as linear combinations of an orthonormal basis, or as linear combinations of eigenvectors of a matrix (or a linear operator). To give a very simple prototype of the Fourier transform, consider a real-valued function f : R → R. Recall that such a function f(x) is even if f(−x) = f(x) for all x ∈ R, and is odd if f(−x) = −f(x) for all x ∈ R. A typical function f, such as f(x) = x3 + 3x2 + 3x + 1, will be neither even nor odd. However, one can always write f as the superposition f = fe + fo of an even function fe and an odd function fo by the formulae f(x) + f(−x) f(x) − f(−x) f (x) := ; f (x) := . e 2 o 2 3 2 2 3 For instance, if f(x) = x + 3x + 3x + 1, then fe(x) = 3x + 1 and fo(x) = x + 3x.
    [Show full text]
  • Subtleties of the Laplace Transform
    Troubles at the Origin: Consistent Usage and Properties of the Unilateral Laplace Transform Kent H. Lundberg, Haynes R. Miller, and David L. Trumper Massachusetts Institute of Technology The Laplace transform is a standard tool associated with the analysis of signals, models, and control systems, and is consequently taught in some form to almost all engineering students. The bilateral and unilateral forms of the Laplace transform are closely related, but have somewhat different domains of application. The bilateral transform is most frequently seen in the context of signal processing, whereas the unilateral transform is most often associated with the study of dynamic system response where the role of initial conditions takes on greater significance. In our teaching we have found some significant pitfalls associated with teaching our students to understand and apply the Laplace transform. These confusions extend to the presentation of this material in many of the available mathematics and engineering textbooks as well. The most significant confusion in much of the textbook literature is how to deal with the origin in the application of the unilateral Laplace transform. That is, many texts present the transform of a time function f(t) as Z ∞ L{f(t)} = f(t)e−st dt (1) 0 without properly specifiying the meaning of the lower limit of integration. Said informally, does the integral include the origin fully, partially, or not at all? This issue becomes significant as soon as singularity functions such as the unit impulse are introduced. While it is not possible to devote full attention to this issue within the context of a typical undergraduate course, this “skeleton in the closet” as Kailath [8] called it needs to be brought out fully into the light.
    [Show full text]
  • A Note on Laplace Transforms of Some Particular Function Types
    A Note on Laplace Transforms of Some Particular Function Types Henrik Stenlund∗ Visilab Signal Technologies Oy, Finland 9th February, 2014 Abstract This article handles in a short manner a few Laplace transform pairs and some extensions to the basic equations are developed. They can be applied to a wide variety of functions in order to find the Laplace transform or its inverse when there is a specific type of an implicit function involved.1 0.1 Keywords Laplace transform, inverse Laplace transform 0.2 Mathematical Classification MSC: 44A10 1 Introduction 1.1 General The following two Laplace transforms have appeared in numerous editions and prints of handbooks and textbooks, for decades. arXiv:1402.2876v1 [math.GM] 9 Feb 2014 1 ∞ 3 s2 2 2 4u Lt[F (t )],s = u− e− f(u)du (F ALSE) (1) 2√π Z0 u 1 f(ln(s)) ∞ t f(u)du L− [ ],t = (F ALSE) (2) s s ln(s) Z Γ(u + 1) · 0 They have mostly been removed from the latest editions and omitted from other new handbooks. However, one fresh edition still carries them [1]. Obviously, many people have tried to apply them. No wonder that they are not accepted ∗The author is obliged to Visilab Signal Technologies for supporting this work. 1Visilab Report #2014-02 1 anymore since they are false. The actual reasons for the errors are not known to the author; possibly it is a misprint inherited from one print to another and then transported to other books, believed to be true. The author tried to apply these transforms, stumbling to a serious conflict.
    [Show full text]
  • An Introduction to Fourier Analysis Fourier Series, Partial Differential Equations and Fourier Transforms
    An Introduction to Fourier Analysis Fourier Series, Partial Differential Equations and Fourier Transforms Notes prepared for MA3139 Arthur L. Schoenstadt Department of Applied Mathematics Naval Postgraduate School Code MA/Zh Monterey, California 93943 August 18, 2005 c 1992 - Professor Arthur L. Schoenstadt 1 Contents 1 Infinite Sequences, Infinite Series and Improper Integrals 1 1.1Introduction.................................... 1 1.2FunctionsandSequences............................. 2 1.3Limits....................................... 5 1.4TheOrderNotation................................ 8 1.5 Infinite Series . ................................ 11 1.6ConvergenceTests................................ 13 1.7ErrorEstimates.................................. 15 1.8SequencesofFunctions.............................. 18 2 Fourier Series 25 2.1Introduction.................................... 25 2.2DerivationoftheFourierSeriesCoefficients.................. 26 2.3OddandEvenFunctions............................. 35 2.4ConvergencePropertiesofFourierSeries.................... 40 2.5InterpretationoftheFourierCoefficients.................... 48 2.6TheComplexFormoftheFourierSeries.................... 53 2.7FourierSeriesandOrdinaryDifferentialEquations............... 56 2.8FourierSeriesandDigitalDataTransmission.................. 60 3 The One-Dimensional Wave Equation 70 3.1Introduction.................................... 70 3.2TheOne-DimensionalWaveEquation...................... 70 3.3 Boundary Conditions ............................... 76 3.4InitialConditions................................
    [Show full text]
  • Laplace Transform
    Chapter 7 Laplace Transform The Laplace transform can be used to solve differential equations. Be- sides being a different and efficient alternative to variation of parame- ters and undetermined coefficients, the Laplace method is particularly advantageous for input terms that are piecewise-defined, periodic or im- pulsive. The direct Laplace transform or the Laplace integral of a function f(t) defined for 0 ≤ t< 1 is the ordinary calculus integration problem 1 f(t)e−stdt; Z0 succinctly denoted L(f(t)) in science and engineering literature. The L{notation recognizes that integration always proceeds over t = 0 to t = 1 and that the integral involves an integrator e−stdt instead of the usual dt. These minor differences distinguish Laplace integrals from the ordinary integrals found on the inside covers of calculus texts. 7.1 Introduction to the Laplace Method The foundation of Laplace theory is Lerch's cancellation law 1 −st 1 −st 0 y(t)e dt = 0 f(t)e dt implies y(t)= f(t); (1) R R or L(y(t)= L(f(t)) implies y(t)= f(t): In differential equation applications, y(t) is the sought-after unknown while f(t) is an explicit expression taken from integral tables. Below, we illustrate Laplace's method by solving the initial value prob- lem y0 = −1; y(0) = 0: The method obtains a relation L(y(t)) = L(−t), whence Lerch's cancel- lation law implies the solution is y(t)= −t. The Laplace method is advertised as a table lookup method, in which the solution y(t) to a differential equation is found by looking up the answer in a special integral table.
    [Show full text]
  • Laplace Transforms: Theory, Problems, and Solutions
    Laplace Transforms: Theory, Problems, and Solutions Marcel B. Finan Arkansas Tech University c All Rights Reserved 1 Contents 43 The Laplace Transform: Basic Definitions and Results 3 44 Further Studies of Laplace Transform 15 45 The Laplace Transform and the Method of Partial Fractions 28 46 Laplace Transforms of Periodic Functions 35 47 Convolution Integrals 45 48 The Dirac Delta Function and Impulse Response 53 49 Solving Systems of Differential Equations Using Laplace Trans- form 61 50 Solutions to Problems 68 2 43 The Laplace Transform: Basic Definitions and Results Laplace transform is yet another operational tool for solving constant coeffi- cients linear differential equations. The process of solution consists of three main steps: • The given \hard" problem is transformed into a \simple" equation. • This simple equation is solved by purely algebraic manipulations. • The solution of the simple equation is transformed back to obtain the so- lution of the given problem. In this way the Laplace transformation reduces the problem of solving a dif- ferential equation to an algebraic problem. The third step is made easier by tables, whose role is similar to that of integral tables in integration. The above procedure can be summarized by Figure 43.1 Figure 43.1 In this section we introduce the concept of Laplace transform and discuss some of its properties. The Laplace transform is defined in the following way. Let f(t) be defined for t ≥ 0: Then the Laplace transform of f; which is denoted by L[f(t)] or by F (s), is defined by the following equation Z T Z 1 L[f(t)] = F (s) = lim f(t)e−stdt = f(t)e−stdt T !1 0 0 The integral which defined a Laplace transform is an improper integral.
    [Show full text]
  • Fourier Transforms & the Convolution Theorem
    Convolution, Correlation, & Fourier Transforms James R. Graham 11/25/2009 Introduction • A large class of signal processing techniques fall under the category of Fourier transform methods – These methods fall into two broad categories • Efficient method for accomplishing common data manipulations • Problems related to the Fourier transform or the power spectrum Time & Frequency Domains • A physical process can be described in two ways – In the time domain, by h as a function of time t, that is h(t), -∞ < t < ∞ – In the frequency domain, by H that gives its amplitude and phase as a function of frequency f, that is H(f), with -∞ < f < ∞ • In general h and H are complex numbers • It is useful to think of h(t) and H(f) as two different representations of the same function – One goes back and forth between these two representations by Fourier transforms Fourier Transforms ∞ H( f )= ∫ h(t)e−2πift dt −∞ ∞ h(t)= ∫ H ( f )e2πift df −∞ • If t is measured in seconds, then f is in cycles per second or Hz • Other units – E.g, if h=h(x) and x is in meters, then H is a function of spatial frequency measured in cycles per meter Fourier Transforms • The Fourier transform is a linear operator – The transform of the sum of two functions is the sum of the transforms h12 = h1 + h2 ∞ H ( f ) h e−2πift dt 12 = ∫ 12 −∞ ∞ ∞ ∞ h h e−2πift dt h e−2πift dt h e−2πift dt = ∫ ( 1 + 2 ) = ∫ 1 + ∫ 2 −∞ −∞ −∞ = H1 + H 2 Fourier Transforms • h(t) may have some special properties – Real, imaginary – Even: h(t) = h(-t) – Odd: h(t) = -h(-t) • In the frequency domain these
    [Show full text]
  • STATISTICAL FOURIER ANALYSIS: CLARIFICATIONS and INTERPRETATIONS by DSG Pollock
    STATISTICAL FOURIER ANALYSIS: CLARIFICATIONS AND INTERPRETATIONS by D.S.G. Pollock (University of Leicester) Email: stephen [email protected] This paper expounds some of the results of Fourier theory that are es- sential to the statistical analysis of time series. It employs the algebra of circulant matrices to expose the structure of the discrete Fourier transform and to elucidate the filtering operations that may be applied to finite data sequences. An ideal filter with a gain of unity throughout the pass band and a gain of zero throughout the stop band is commonly regarded as incapable of being realised in finite samples. It is shown here that, to the contrary, such a filter can be realised both in the time domain and in the frequency domain. The algebra of circulant matrices is also helpful in revealing the nature of statistical processes that are band limited in the frequency domain. In order to apply the conventional techniques of autoregressive moving-average modelling, the data generated by such processes must be subjected to anti- aliasing filtering and sub sampling. These techniques are also described. It is argued that band-limited processes are more prevalent in statis- tical and econometric time series than is commonly recognised. 1 D.S.G. POLLOCK: Statistical Fourier Analysis 1. Introduction Statistical Fourier analysis is an important part of modern time-series analysis, yet it frequently poses an impediment that prevents a full understanding of temporal stochastic processes and of the manipulations to which their data are amenable. This paper provides a survey of the theory that is not overburdened by inessential complications, and it addresses some enduring misapprehensions.
    [Show full text]
  • Fourier Transform, Convolution Theorem, and Linear Dynamical Systems April 28, 2016
    Mathematical Tools for Neuroscience (NEU 314) Princeton University, Spring 2016 Jonathan Pillow Lecture 23: Fourier Transform, Convolution Theorem, and Linear Dynamical Systems April 28, 2016. Discrete Fourier Transform (DFT) We will focus on the discrete Fourier transform, which applies to discretely sampled signals (i.e., vectors). Linear algebra provides a simple way to think about the Fourier transform: it is simply a change of basis, specifically a mapping from the time domain to a representation in terms of a weighted combination of sinusoids of different frequencies. The discrete Fourier transform is therefore equiv- alent to multiplying by an orthogonal (or \unitary", which is the same concept when the entries are complex-valued) matrix1. For a vector of length N, the matrix that performs the DFT (i.e., that maps it to a basis of sinusoids) is an N × N matrix. The k'th row of this matrix is given by exp(−2πikt), for k 2 [0; :::; N − 1] (where we assume indexing starts at 0 instead of 1), and t is a row vector t=0:N-1;. Recall that exp(iθ) = cos(θ) + i sin(θ), so this gives us a compact way to represent the signal with a linear superposition of sines and cosines. The first row of the DFT matrix is all ones (since exp(0) = 1), and so the first element of the DFT corresponds to the sum of the elements of the signal. It is often known as the \DC component". The next row is a complex sinusoid that completes one cycle over the length of the signal, and each subsequent row has a frequency that is an integer multiple of this \fundamental" frequency.
    [Show full text]
  • M. A. Pathan, O. A. Daman LAPLACE TRANSFORMS of the LOGARITHMIC FUNCTIONS and THEIR APPLICATIONS 1. Introduction the Aim of This
    DEMONSTRATIO MATHEMATICA Vol. XLVI No 3 2013 M. A. Pathan, O. A. Daman LAPLACE TRANSFORMS OF THE LOGARITHMIC FUNCTIONS AND THEIR APPLICATIONS Abstract. This paper deals with theorems and formulas using the technique of Laplace and Steiltjes transforms expressed in terms of interesting alternative logarith- mic and related integral representations. The advantage of the proposed technique is illustrated by logarithms of integrals of importance in certain physical and statistical problems. 1. Introduction The aim of this paper is to obtain some theorems and formulas for the evaluation of finite and infinite integrals for logarithmic and related func- tions using technique of Laplace transform. Basic properties of Laplace and Steiltjes transforms and Parseval type relations are explicitly used in combi- nation with rules and theorems of operational calculus. Some of the integrals obtained here are related to stochastic calculus [6] and common mathemati- cal objects, such as the logarithmic potential [3], logarithmic growth [2] and Whittaker functions [2, 3, 4, 6, 7] which are of importance in certain physical and statistical applications, in particular in energies, entropies [3, 5, 7 (22)], intermediate moment problem [2] and quantum electrodynamics. The ad- vantage of the proposed technique is illustrated by the explicit computation of a number of different types of logarithmic integrals. We recall here the definition of the Laplace transform 1 −st (1) L ff (t)g = L[f(t); s] = \ e f (t) dt: 0 Closely related to the Laplace transform is the generalized Stieltjes trans- form 2010 Mathematics Subject Classification: 33C05, 33C90, 44A10. Key words and phrases: Laplace and Steiltjes transforms, logarithmic functions and integrals, stochastic calculus and Whittaker functions.
    [Show full text]