A Wavelet-based Method for Overcoming the Gibbs Phenomenon
|
|
- Ophelia Richard
- 5 years ago
- Views:
Transcription
1 AMERICAN CONFERENCE ON APPLIED MATHEMATICS (MATH '08), Harvard, Massachusetts, USA, March 4-6, 008 A Wavelet-based Method for Overcoming the Gibbs Phenomenon NATANIEL GREENE Department of Mathematics Computer Science Kingsborough Community College, CUNY 00 Oriental Boulevard, Brooklyn, NY 35 UNITED STATES ngreene.math@gmail.com Abstract: The Gibbs phenomenon refers to the lack of uniform convergence which occurs in many orthogonal basis approximations to piecewise smooth functions. This lack of uniform convergence manifests itself in spurious oscillations near the points of discontinuity a low order of convergence away from the discontinuities. Here we describe a numerical procedure for overcoming the Gibbs phenomenon called the inverse wavelet reconstruction method. The method takes the Fourier coefcients of an oscillatory partial sum uses them to construct the wavelet coefcients of a non-oscillatory wavelet series. Key Words: Gibbs phenomenon, wavelets, Gegenbauer reconstruction, inverse polynomial reconstruction. Introduction Fourier orthogonal polynomial series are known for their highly accurate expansions for smooth functions. In fact it is known that the more derivatives a function has, the faster the approximation will converge. However, when a function possesses jumpdiscontinuities the approximation will fail to converge uniformly. In addition, spurious oscillations will cause a loss of accuracy throughout the entire domain. This lack of uniform convergence is known as the Gibbs phenomenon. Methods for post-processing approximations which suffer from the Gibbs phenomenon include the Gegenbauer reconstruction method of Gottlieb Shu [7,9,0], the method of Pade approximants due to Driscoll Fornberg [], the method of spectral molliers due to Gottlieb Tadmor [8] Tadmor Tanner [7, 8], the inverse polynomial reconstruction method of Shizgal Jung [3,4,5,6], the Freund polynomial reconstruction method of Gelb Tanner [6]. These reconstruction methods can be combined with an effective method for edge-detection developed by Gelb Tadmor [3,4,5,6], to yield an exponentially accurate reconstruction of the original function. In this paper we describe a new numerical method for overcoming the Gibbs phenomenon following the work of Shizgal Jung, called the inverse wavelet reconstruction method. We begin with a brief review of the essential denitions of wavelets which we will need. Recall that a wavelet is a function L (R) satisfying: Z Z (x) dx 0 () b () d < ; () jj b here is the Fourier transform of. The function is also known as an analyzing wavelet or a mother wavelet since any function f L (R) can be expressed as a continuous sum of translations dilations involving according to the continuous wavelet transform. The continuous wavelet transform is given by Z t b (W f) (b; a) jaj f (t) dt a the inverse continuous wavelet transform is given by f (x) C Z Z b;a (t) jaj C Z [(W f) (b; a)] t b a b () d: jj b;a (t) a dadb (3) ISSN: ISBN:
2 A discrete wavelet series is given by f (x) X X j k c j;k j;k (x) (4) the discrete wavelet series coefcients are given by k c j;k (W f) j ; j AMERICAN CONFERENCE ON APPLIED MATHEMATICS (MATH '08), Harvard, Massachusetts, USA, March 4-6, 008 j;k (x) j j x k. For more information see, for example, Chui []. We will make use of wavelet series for our reconstruction technique further on. The Fourier coefcients b f (n) can be expressed in terms of Gegenbauer- coefcients as follows. bf (n) Z Z X X f (x) e inx dx X bf (m) bf (m) C m (x) e inx dx Z Cm (x) e inx dx bf (m) c m;n (7) The connection coefcients c m;n are given by The Inverse Wavelet Reconstruction Method. Inverse Polynomial Reconstruction One numerical technique for dealing with the problem of the Gibbs phenomenon is the inverse polynomial reconstruction method of Shizgal Jung [3, 4, 5, 6]. Their inverse method is itself an alternative approach to the original direct Gegenbauer reconstruction method of Gottlieb Shu [7, 9, 0]. In the inverse method, one solves a system of linear equations for the Gegenbauer polynomial reconstruction coefcients in terms original Fourier coefcients, whose expansion suffers from the Gibbs phenomenon. The inverse polynomial reconstruction method begins with a Fourier partial sum S N f (x) NX n N bf (n) e inx : (5) The function also has a Gegenbauer- expansion given by bf (m) h m f (x) Z X bf (m) C m (x) h m (m + ) m! () f (x) C m (x) x dx (6) c m;n Z C m (x) e inx dx (8) can be computed numerically. An explicit formula for c m;n was derived by this author, Greene [, ], is given for the cases n 6 0 n 0 respectively by: c m;n c m;0 [m] X j0 8 < : ( i) m j p J n m j+ (n) (9) j () 3 m j m j j! 3 m j m j ( ) m () m j!( 3 ) m ; m even 0; m odd. (0) The inverse polynomial reconstruction procedure is obtained by truncating the innite series above solving the system of equations bf (n) NX bf (m) c m;n, n N:::N: () One then computes the Gegenbauer reconstruction approximation S Nf (x) NX bf (m) C m (x) : () + () (m + ) :. Inverse Wavelet Reconstruction We explore here an analogous approach which seeks to reconstruct the original function in terms of ISSN: ISBN:
3 wavelets. We begin with a Fourier partial sum as before assume that the function f (x) can also be expressed as a discrete wavelet series: f (x) X X m l Z m;l (x) f (x) m;l (x)dx m;l (x) m ( m x l) : We now derive a formula expressing the Fourier coef- cients in terms of wavelet coefcients. bf (n) Z Z X f (x) e inx dx X X m l X m l X X m l m;l (x) e inx dx Z m;l (x) e inx dx b m;l (n) (3) The inverse wavelet reconstruction method is obtained by truncating the doubly innite sum described above solving for the wavelet coefcients. We suggest solving the following system of equations. AMERICAN CONFERENCE ON APPLIED MATHEMATICS (MATH '08), Harvard, Massachusetts, USA, March 4-6, 008 bf (n) m M l L b m;l (n) (4) n N:::N ; N ML (5) for the wavelet coefcients : One then computes the wavelet reconstruction approximation S M;L f (x) m M l L m;l (x) : (6) The terms b m;l (n) are the nth Fourier coefcients of m;l (x) : Since we are solving a system of N equations in (M) (L) 4ML unknown wavelet coefcients, in order for the system to be invertible we must have N 4ML or N ML: It is this convention which we will use for our numerical implementation below, owing to the simplicity of the relation N ML: Alternatively, one may solve the system bf (n) m M l L b m;l (n) (7) n N:::N; N ML + M + L: (8) then compute S M;L f (x) m M l L m;l (x) : (9) The condition in (8) is due to the fact that we are solving a system of N + equations in (M + ) (L + ) unknown wavelet coefcients. In order for the system to be invertible we must have N + (M + ) (L + ) or N ML+M +L: A detailed proof of convergence is the subject of concurrent work. However, we illustrate the apparent convergence for certain wavelets with some numerical examples. 3 Numerical Results Numerical experiments show that the inverse wavelet reconstruction approach does yield accurate uniform approximations for a variety of wavelet families. We illustrate this with three wavelet families: Poisson wavelets given by (x) Mexican hat wavelets given by (x) p 3 4 ( + x ) ; (0) x e x ; () Morlet wavelets given by (x) cos x e x : () For each case we begin with a Fourier series of N 6; 36; 64; 00 expansion coefcients of an analytic non-periodic test function f(x) 4 tan (x) whose partial sum suffers from the Gibbs phenomenon. We then compute the corresponding wavelet series reconstructions. We chose the number of Fourier coefcients such that N ML M L: The absolute error of these reconstructions, jf(x) S M;L f (x)j, are displayed in the gures provided. For the three wavelets shown the Morlets perform the best for a xed number of Fourier coefcients, the Mexican hat wavelets perform comparably, though slightly less well, the Poisson wavelets perform the least well. ISSN: ISBN:
4 AMERICAN CONFERENCE ON APPLIED MATHEMATICS (MATH '08), Harvard, Massachusetts, USA, March 4-6, 008 Figure : The graphs illustrate the Gibbs phenomenon for a 6 term Fourier partial sum. The 6 Fourier coefcients are used to reconstruct a 6 term Morlet wavelet series accurate to four decimal places. Figure 3: The absolute error of a Poisson wavelet reconstruction based on 6, 36, 64, 00 Fourier coefcients. Figure : The absolute error of a Morlet wavelet reconstruction based on 6, 36, 64, 00 Fourier coefcients. Figure 4: The absolute error of a Mexican hat wavelet reconstruction based on 6, 36, 64, 00 Fourier coefcients. ISSN: ISBN:
5 4 Conclusions The numerical results indicate that the inverse wavelet reconstruction method yields an accurate uniformly converging reconstruction approximation for a variety of wavelets. Current work underway includes a study of the technique for a broader spectrum of wavelets, reconstruction from series other than Fourier series, such as orthogonal polynomials or wavelets, comparison with other methods, an analytic estimation of error proof of convergence. Acknowledgements: This research was supported in part by a PSC-CUNY Research Award (grant No ). References AMERICAN CONFERENCE ON APPLIED MATHEMATICS (MATH '08), Harvard, Massachusetts, USA, March 4-6, 008 [] C. K. Chui, An Introduction to Wavelets, Academic Press, San Diego, 99. [] T. A. Driscoll, B. Fornberg, A Pade-based Algorithm for Overcoming the Gibbs Phenomenon, Numerical Algorithms, 6, 00, pp [3] A. Gelb, E. Tadmor, Detection of Edges in Spectral Data, Appl. Comp. Harmonic Anal., 7, 999, pp [4] A. Gelb, E. Tadmor, Detection of Edges in Spectral Data II. Nonlinear Enhancement, SIAM J. Numer. Anal., Vol. 38, No. 4, 000, pp [5] A. Gelb, E. Tadmor, Spectral Reconstruction of Piecewise Smooth Functions from their Discrete Data, Mathematical Modeling Numerical Analysis, 36:, 00, pp [6] A. Gelb, J. Tanner, Robust Reprojection Methods for the Resolution of the Gibbs Phenomenon., Applied Computational Harmonic Analysis, Vol. 0,, 006, pp [7] D. Gottlieb, C.-W. Shu, A. Solomonoff, H. Veven, On the Gibbs's Phenomenon I: Recovering Exponential Accuracy from the Fourier Partial Sum of a Nonperiodic Analytic Function, J. Comput. Appl. Math., 43, 99, pp [8] D. Gottlieb, E. Tadmor, Recovering Pointwise Values of Discontinuous Data within Spectral Accuracy, in: Progress Supercomputing in Computational Fluid Dynamics, Proceedings of a 984 U.S.-Israel Workshop, Progress in Scientic Computing, Vol. 6 (E. M. Murman S. S. Abarbanel eds.), Birkhauer, Boston, 985, pp [9] D. Gottlieb C.-W. Shu, On the Gibbs Phenomenon its Resolution, SIAM Rev., Vol. 39, 997, pp [0] D. Gottlieb C.-W. Shu, A General Theory for the Resolution of the Gibbs Phenomenon, in: Tricomi's Ideas Contemporary Applied Mathematics, National Italian Academy of Science, 997. [] N. Greene, On the Recovery of Piecewise Smooth Functions from their Integral Transforms Spectral Data, Ph.D. Dissertation, SUNY Stony Brook, 004. [] N. Greene, Fourier Series of Orthogonal Polynomials, in: Proceedings of the American Conference on Applied Mathematics, Cambridge, Massachusetts, USA, March 4-6, 008. [3] J.-H. Jung, B. D. Shizgal, Generalization of the Inverse Polynomial Reconstruction Method in the Resolution of the Gibbs Phenomenon, J. Comp. Appl. Math., v. 7, n., 004, pp.3-5. [4] J.-H. Jung, B. D. Shizgal, Inverse polynomial reconstruction of two dimensional Fourier images, J. Scientic Computing, v.5, n.3, 005, pp [5] B. D. Shizgal J.-H. Jung, Towards the resolution of the Gibbs Phenomena, J. Comput. Appl. Math., 6, 003, pp [6] R. Pasquetti, On Inverse methods for the Resolution of the Gibbs Phenomenon, Journal of Computational Applied Mathematics, Vol. 70, no., 004, pp [7] E. Tadmor J. Tanner, Adaptive Molliers - High Resolution Recovery of Piecewise Smooth Data from its Spectral Information, J. Foundations of Comp. Math., 00, [8] E. Tadmor J. Tanner, Adaptive Filters for Piecewise Smooth Spectral Data, IMA J. Numerical Anal., Vol. 5, 4, 005, pp ISSN: ISBN:
Recovery of high order accuracy in radial basis function approximation for discontinuous problems
Recovery of high order accuracy in radial basis function approximation for discontinuous problems Chris L. Bresten, Sigal Gottlieb 1, Daniel Higgs, Jae-Hun Jung* 2 Abstract Radial basis function(rbf) methods
More information(4) k = 2N + I. A Numerical Study on the Accuracy of Fourier Spectral Methods Applied to the Nonlinear Burgers Equation.
A Numerical Study on the Accuracy of Fourier Spectral Methods Applied to the Nonlinear Burgers Equation Chi-Wang Shu* Peter S. Wong* Abstract It is well known that spectral method can achieve exponential
More informationAPPROXIMATING THE INVERSE FRAME OPERATOR FROM LOCALIZED FRAMES
APPROXIMATING THE INVERSE FRAME OPERATOR FROM LOCALIZED FRAMES GUOHUI SONG AND ANNE GELB Abstract. This investigation seeks to establish the practicality of numerical frame approximations. Specifically,
More informationA Singular Integral Transform for the Gibbs-Wilbraham Effect in Inverse Fourier Transforms
Universal Journal of Integral Equations 4 (2016), 54-62 www.papersciences.com A Singular Integral Transform for the Gibbs-Wilbraham Effect in Inverse Fourier Transforms N. H. S. Haidar CRAMS: Center for
More informationAn Accurate Fourier-Spectral Solver for Variable Coefficient Elliptic Equations
An Accurate Fourier-Spectral Solver for Variable Coefficient Elliptic Equations Moshe Israeli Computer Science Department, Technion-Israel Institute of Technology, Technion city, Haifa 32000, ISRAEL Alexander
More informationApplied and Computational Harmonic Analysis 7, (1999) Article ID acha , available online at
Applied and Computational Harmonic Analysis 7, 0 35 (999) Article ID acha.999.06, available online at http://www.idealibrary.com on Detection of Edges in Spectral Data Anne Gelb Department of Mathematics,
More informationPAPER Resolution of the Gibbs Phenomenon for Fractional Fourier Series
57 PAPER Resolution of the Gibbs Phenomenon for Fractional Fourier Series Hongqing ZHU a), Member, Meiyu DING, and Daqi GAO, Nonmembers SUMMARY The nth partial sums of a classical Fourier series have large
More informationA FAST SOLVER FOR ELLIPTIC EQUATIONS WITH HARMONIC COEFFICIENT APPROXIMATIONS
Proceedings of ALGORITMY 2005 pp. 222 229 A FAST SOLVER FOR ELLIPTIC EQUATIONS WITH HARMONIC COEFFICIENT APPROXIMATIONS ELENA BRAVERMAN, MOSHE ISRAELI, AND ALEXANDER SHERMAN Abstract. Based on a fast subtractional
More informationScientific Computing: An Introductory Survey
Scientific Computing: An Introductory Survey Chapter 12 Prof. Michael T. Heath Department of Computer Science University of Illinois at Urbana-Champaign Copyright c 2002. Reproduction permitted for noncommercial,
More informationSingular Padé-Chebyshev Reconstruction: The Case When Jump Locations are Unknown
Singular Padé-Chebyshev Reconstruction: The Case When Jump Locations are Unknown Arnel L. Tampos Jose Ernie C. Lope Jan S. Hesthaven Abstract We present a singularity-based approach to resolve the Gibbs
More informationFOURIER PADÉ APPROXIMATIONS AND FILTERING FOR SPECTRAL SIMULATIONS OF AN INCOMPRESSIBLE BOUSSINESQ CONVECTION PROBLEM
MATHEMATICS OF COMPUTATION Volume 7, Number 9, July 7, Pages 7 9 S -78(7)8- Article electronically published on February, 7 FOURIER PADÉ APPROXIMATIONS AND FILTERING FOR SPECTRAL SIMULATIONS OF AN INCOMPRESSIBLE
More informationEnhanced spectral viscosity approximations for conservation laws
Applied Numerical Mathematics 33 (000) 3 1 Enhanced spectral viscosity approximations for conservation laws Anne Gelb a,1, Eitan Tadmor b, a Department of Mathematics, Arizona State University, Tempe,
More informationSolving Poisson s Equations Using Buffered Fourier Spectral Method
Solving Poisson s Equations Using Buffered Fourier Spectral Method Yinlin Dong Hassan Abd Salman Al-Dujaly Chaoqun Liu Technical Report 2015-12 http://www.uta.edu/math/preprint/ Solving Poisson s Equations
More informationAdaptive Edge Detectors for Piecewise Smooth Data Based on the minmod Limiter
Journal of Scientific Computing, Vol. 8, os. /3, September 6 ( 6) DOI:.7/s95-988 Adaptive Edge Detectors for Piecewise Smooth Data Based on the minmod Limiter A. Gelb and E. Tadmor Received April 9, 5;
More information3 rd class Mech. Eng. Dept. hamdiahmed.weebly.com Fourier Series
Definition 1 Fourier Series A function f is said to be piecewise continuous on [a, b] if there exists finitely many points a = x 1 < x 2
More informationPADÉ-BASED INTERPRETATION AND CORRECTION OF THE GIBBS PHENOMENON
PADÉ-BASED INTERPRETATION AND CORRECTION OF THE GIBBS PHENOMENON TOBIN A. DRISCOLL AND BENGT FORNBERG Abstract. The convergence of a Fourier series on an interval can be interpreted naturally as the convergence
More informationSpace-Frequency Atoms
Space-Frequency Atoms FREQUENCY FREQUENCY SPACE SPACE FREQUENCY FREQUENCY SPACE SPACE Figure 1: Space-frequency atoms. Windowed Fourier Transform 1 line 1 0.8 0.6 0.4 0.2 0-0.2-0.4-0.6-0.8-1 0 100 200
More informationSpectral collocation and waveform relaxation methods with Gegenbauer reconstruction for nonlinear conservation laws
Spectral collocation and waveform relaxation methods with Gegenbauer reconstruction for nonlinear conservation laws Z. Jackiewicz and B. Zubik Kowal November 2, 2004 Abstract. We investigate Chebyshev
More information* ICASE G 11t w. t )t Operated by Universities Space Research Association
[ NASA Contractor Report 194959 0 ICASE Report No. 94-67 AD-A284 063 * ICASE A NOTE ON THE ACCURACY OF SPECTRAL METHOD APPLIED TO NONLINEAR CONSERVATION LAWS Chi-Wang Shu Peter S. Wong ( ELECTED biic SEP
More informationThe Resolution of the Gibbs Phenomenon for Fourier Spectral Methods
The Resolution of the Gibbs Phenomenon for Fourier Spectral Methods Anne Gelb Sigal Gottlieb December 5, 6 Introduction Fourier spectral methods have emerged as powerful computational techniques for the
More informationHighly Oscillatory Problems: Computation, Theory and Applications
Highly Oscillatory Problems: Computation, Theory and Applications Isaac Newton Institute for Mathematical Sciences Edited by 1 Rapid function approximation by modified Fourier series Daan Huybrechs and
More informationSpace-Frequency Atoms
Space-Frequency Atoms FREQUENCY FREQUENCY SPACE SPACE FREQUENCY FREQUENCY SPACE SPACE Figure 1: Space-frequency atoms. Windowed Fourier Transform 1 line 1 0.8 0.6 0.4 0.2 0-0.2-0.4-0.6-0.8-1 0 100 200
More information1 Introduction to Wavelet Analysis
Jim Lambers ENERGY 281 Spring Quarter 2007-08 Lecture 9 Notes 1 Introduction to Wavelet Analysis Wavelets were developed in the 80 s and 90 s as an alternative to Fourier analysis of signals. Some of the
More informationLinear Algebra: Characteristic Value Problem
Linear Algebra: Characteristic Value Problem . The Characteristic Value Problem Let < be the set of real numbers and { be the set of complex numbers. Given an n n real matrix A; does there exist a number
More informationApplied and Computational Harmonic Analysis
Appl. Comput. Harmon. Anal. 27 (2009) 351 366 Contents lists available at ScienceDirect Applied and Computational Harmonic Analysis www.elsevier.com/locate/acha Letter to the Editor Nonlinear inversion
More informationMath 489AB A Very Brief Intro to Fourier Series Fall 2008
Math 489AB A Very Brief Intro to Fourier Series Fall 8 Contents Fourier Series. The coefficients........................................ Convergence......................................... 4.3 Convergence
More informationCONTROL SYSTEMS, ROBOTICS, AND AUTOMATION Vol. VI - System Identification Using Wavelets - Daniel Coca and Stephen A. Billings
SYSTEM IDENTIFICATION USING WAVELETS Daniel Coca Department of Electrical Engineering and Electronics, University of Liverpool, UK Department of Automatic Control and Systems Engineering, University of
More informationTime-Frequency Analysis
Time-Frequency Analysis Basics of Fourier Series Philippe B. aval KSU Fall 015 Philippe B. aval (KSU) Fourier Series Fall 015 1 / 0 Introduction We first review how to derive the Fourier series of a function.
More informationWavelets and applications
Chapter 3 Wavelets and applications 3. Multiresolution analysis 3.. The limits of Fourier analysis Although along this chapter the underlying Hilbert space will be L 2 (R), we start with a completely explicit
More informationEdge Detection Using Fourier Coefficients
Edge Detection Using Fourier Coefficients Shlomo Engelberg. INTRODUCTION. Edge detection and the detection of discontinuities are important in many fields. In image processing, for example, one often needs
More informationExponential Approximations Using Fourier Series Partial Sums
NASA/CR-201751 ICASE Report No. 97-56 Exponential Approximations Using Fourier Series Partial Sums Nana S. Baneree Advanta Corporation and James F. Geer State University of New York Institute for Computer
More informationKey Words: cardinal B-spline, coefficients, moments, rectangular rule, interpolating quadratic spline, hat function, cubic B-spline.
M a t h e m a t i c a B a l a n i c a New Series Vol. 24, 2, Fasc.3-4 Splines in Numerical Integration Zlato Udovičić We gave a short review of several results which are related to the role of splines
More informationStrong stability preserving high order time discretizations.
SIGAL GOTTLIEB Mathematics Department 285 Old Westport Road North Dartmouth, MA 02747 sgottlieb@umassd.edu 322 Cole Avenue Providence, RI 02906 Phone: (401) 751-9416 sigalgottlieb@yahoo.com Current Research
More informationRecovering Exponential Accuracy From Non-Harmonic Fourier Data Through Spectral Reprojection
Recovering Eponential Accuracy From Non-Harmonic Fourier Data Through Spectral Reprojection Anne Gelb Taylor Hines November 5, 2 Abstract Spectral reprojection techniques make possible the recovery of
More informationWeighted ENO Schemes
Xiaolei Chen Advisor: Prof. Xiaolin Li Department of Applied Mathematics and Statistics Stony Brook University, The State University of New York February 7, 014 1 3 Mapped WENO-Z Scheme 1D Scalar Hyperbolic
More informationarxiv: v1 [math.cv] 18 Aug 2015
arxiv:508.04376v [math.cv] 8 Aug 205 Saddle-point integration of C bump functions Steven G. Johnson, MIT Applied Mathematics Created November 23, 2006; updated August 9, 205 Abstract This technical note
More informationNumerical Approximation Methods for Non-Uniform Fourier Data
Numerical Approximation Methods for Non-Uniform Fourier Data Aditya Viswanathan aditya@math.msu.edu 2014 Joint Mathematics Meetings January 18 2014 0 / 16 Joint work with Anne Gelb (Arizona State) Guohui
More informationMore on Fourier Series
More on Fourier Series R. C. Trinity University Partial Differential Equations Lecture 6.1 New Fourier series from old Recall: Given a function f (x, we can dilate/translate its graph via multiplication/addition,
More informationHierarchical Reconstruction with up to Second Degree Remainder for Solving Nonlinear Conservation Laws
Hierarchical Reconstruction with up to Second Degree Remainder for Solving Nonlinear Conservation Laws Dedicated to Todd F. Dupont on the occasion of his 65th birthday Yingjie Liu, Chi-Wang Shu and Zhiliang
More informationSpectral Methods for Hyperbolic Problems
Spectral Methods for Hyperbolic Problems J. S. Hesthaven EPFL-SB-MATHICSE-MCSS, Ecole Polytechnique Fédérale de Lausanne CH-1015 Lausanne, Switzerland E-mail: Jan.Hesthaven@epfl.ch We review spectral methods
More informationConstructing Approximation Kernels for Non-Harmonic Fourier Data
Constructing Approximation Kernels for Non-Harmonic Fourier Data Aditya Viswanathan aditya.v@caltech.edu California Institute of Technology SIAM Annual Meeting 2013 July 10 2013 0 / 19 Joint work with
More informationEmily Jennings. Georgia Institute of Technology. Nebraska Conference for Undergraduate Women in Mathematics, 2012
δ 2 Transform and Fourier Series of Functions with Multiple Jumps Georgia Institute of Technology Nebraska Conference for Undergraduate Women in Mathematics, 2012 Work performed at Kansas State University
More informationComputing Derivatives of Noisy Signals Using Orthogonal Functions Expansions.
Computing Derivatives of Noisy Signals Using Orthogonal Functions Expansions. Adi Ditkowski, Abhinav Bhandari, Brian W. Sheldon 16 Nov. 8 Abstract In many applications noisy signals are measured. These
More informationAn Introduction to Wavelets and some Applications
An Introduction to Wavelets and some Applications Milan, May 2003 Anestis Antoniadis Laboratoire IMAG-LMC University Joseph Fourier Grenoble, France An Introduction to Wavelets and some Applications p.1/54
More informationNEW BOUNDS FOR TRUNCATION-TYPE ERRORS ON REGULAR SAMPLING EXPANSIONS
NEW BOUNDS FOR TRUNCATION-TYPE ERRORS ON REGULAR SAMPLING EXPANSIONS Nikolaos D. Atreas Department of Mathematics, Aristotle University of Thessaloniki, 54006, Greece, e-mail:natreas@auth.gr Abstract We
More informationFourier Transform Chapter 10 Sampling and Series
Fourier Transform Chapter 0 Sampling and Series Sampling Theorem Sampling Theorem states that, under a certain condition, it is in fact possible to recover with full accuracy the values intervening between
More information14 Fourier analysis. Read: Boas Ch. 7.
14 Fourier analysis Read: Boas Ch. 7. 14.1 Function spaces A function can be thought of as an element of a kind of vector space. After all, a function f(x) is merely a set of numbers, one for each point
More informationWavelets: Theory and Applications. Somdatt Sharma
Wavelets: Theory and Applications Somdatt Sharma Department of Mathematics, Central University of Jammu, Jammu and Kashmir, India Email:somdattjammu@gmail.com Contents I 1 Representation of Functions 2
More informationFourier series of sums of products of Bernoulli functions and their applications
Available online at www.isr-publications.com/jnsa J. Nonlinear Sci. Appl., 7, 798 85 Research Article Journal Homepage: www.tjnsa.com - www.isr-publications.com/jnsa Fourier series of sums of products
More informationIMAGE RECONSTRUCTION FROM UNDERSAMPLED FOURIER DATA USING THE POLYNOMIAL ANNIHILATION TRANSFORM
IMAGE RECONSTRUCTION FROM UNDERSAMPLED FOURIER DATA USING THE POLYNOMIAL ANNIHILATION TRANSFORM Anne Gelb, Rodrigo B. Platte, Rosie Renaut Development and Analysis of Non-Classical Numerical Approximation
More informationMath 115 ( ) Yum-Tong Siu 1. Derivation of the Poisson Kernel by Fourier Series and Convolution
Math 5 (006-007 Yum-Tong Siu. Derivation of the Poisson Kernel by Fourier Series and Convolution We are going to give a second derivation of the Poisson kernel by using Fourier series and convolution.
More informationMathematical Methods for Computer Science
Mathematical Methods for Computer Science Computer Laboratory Computer Science Tripos, Part IB Michaelmas Term 2016/17 Professor J. Daugman Exercise problems Fourier and related methods 15 JJ Thomson Avenue
More informationarxiv: v1 [math.na] 21 Nov 2017
High Order Finite Difference Schemes for the Heat Equation Whose Convergence Rates are Higher Than Their Truncation Errors, A. Ditkowski arxiv:7.0796v [math.na] Nov 07 Abstract Typically when a semi-discrete
More informationGeneralized Least-Squares Regressions III: Further Theory and Classication
Generalized Least-Squares Regressions III: Further Theor Classication NATANIEL GREENE Department of Mathematics Computer Science Kingsborough Communit College, CUNY Oriental Boulevard, Brookln, NY UNITED
More informationMoment Computation in Shift Invariant Spaces. Abstract. An algorithm is given for the computation of moments of f 2 S, where S is either
Moment Computation in Shift Invariant Spaces David A. Eubanks Patrick J.Van Fleet y Jianzhong Wang ẓ Abstract An algorithm is given for the computation of moments of f 2 S, where S is either a principal
More informationPHYS 502 Lecture 3: Fourier Series
PHYS 52 Lecture 3: Fourier Series Fourier Series Introduction In mathematics, a Fourier series decomposes periodic functions or periodic signals into the sum of a (possibly infinite) set of simple oscillating
More informationMath Numerical Analysis
Math 541 - Numerical Analysis Joseph M. Mahaffy, jmahaffy@mail.sdsu.edu Department of Mathematics and Statistics Dynamical Systems Group Computational Sciences Research Center San Diego State University
More informationChapter 6. Legendre and Bessel Functions
Chapter 6 Legendre and Bessel Functions Legendre's Equation Legendre's Equation (order n): legendre d K y''k y'c n n C y = : is an important ode in applied mathematics When n is a non-negative integer,
More informationApproximation theory
Approximation theory Xiaojing Ye, Math & Stat, Georgia State University Spring 2019 Numerical Analysis II Xiaojing Ye, Math & Stat, Georgia State University 1 1 1.3 6 8.8 2 3.5 7 10.1 Least 3squares 4.2
More informationOPTIMAL FILTER AND MOLLIFIER FOR PIECEWISE SMOOTH SPECTRAL DATA
MATHEMATICS OF COMPUTATIO Volume, umber, Pages S 5-578(XX)- OPTIMAL FILTER AD MOLLIFIER FOR PIECEWISE SMOOTH SPECTRAL DATA JARED TAER This paper is dedicated to Eitan Tadmor for his direction. Abstract.
More informationPeriodic functions: simple harmonic oscillator
Periodic functions: simple harmonic oscillator Recall the simple harmonic oscillator (e.g. mass-spring system) d 2 y dt 2 + ω2 0y = 0 Solution can be written in various ways: y(t) = Ae iω 0t y(t) = A cos
More informationHarmonic Wavelet Transform and Image Approximation
Harmonic Wavelet Transform and Image Approximation Zhihua Zhang and Naoki Saito Dept of Math, Univ of California, Davis, California, 95616, USA Email: zhangzh@mathucdavisedu saito@mathucdavisedu Abstract
More informationApproximation by Conditionally Positive Definite Functions with Finitely Many Centers
Approximation by Conditionally Positive Definite Functions with Finitely Many Centers Jungho Yoon Abstract. The theory of interpolation by using conditionally positive definite function provides optimal
More informationModified Fourier Spectral Method for Non-periodic CFD
Modified Fourier Spectral Method for Non-periodic CFD Huankun Fu Chaoqun Liu Technical Report 2011-14 http://www.uta.edu/math/preprint/ Modified Fourier Spectral Method for Non-periodic CFD Huankun Fu
More informationSuper Viscosity and Spectral Approximations of Nonlinear Conservation Laws
umerical Methods for Fluid Dynamics. IV M. J. Baines and K. W. Morton, eds. Clarendon Press, Oxford, 1993, pp. 69-82 Super Viscosity and Spectral Approximations of onlinear Conservation Laws EITA TADMOR
More informationNinth International Water Technology Conference, IWTC9 2005, Sharm El-Sheikh, Egypt 673
Ninth International Water Technology Conference, IWTC9 2005, Sharm El-Sheikh, Egypt 673 A NEW NUMERICAL APPROACH FOR THE SOLUTION OF CONTAMINANT TRANSPORT EQUATION Mohamed El-Gamel Department of Mathematical
More informationRecurrence Relations and Fast Algorithms
Recurrence Relations and Fast Algorithms Mark Tygert Research Report YALEU/DCS/RR-343 December 29, 2005 Abstract We construct fast algorithms for decomposing into and reconstructing from linear combinations
More informationDiscrete Wavelet Transform: A Technique for Image Compression & Decompression
Discrete Wavelet Transform: A Technique for Image Compression & Decompression Sumit Kumar Singh M.Tech Scholar, Deptt. of Computer Science & Engineering Al-Falah School of Engineering & Technology, Faridabad,
More informationWavelets in Scattering Calculations
Wavelets in Scattering Calculations W. P., Brian M. Kessler, Gerald L. Payne polyzou@uiowa.edu The University of Iowa Wavelets in Scattering Calculations p.1/43 What are Wavelets? Orthonormal basis functions.
More informationCompletion Date: Monday February 11, 2008
MATH 4 (R) Winter 8 Intermediate Calculus I Solutions to Problem Set #4 Completion Date: Monday February, 8 Department of Mathematical and Statistical Sciences University of Alberta Question. [Sec..9,
More informationPolynomial Approximation: The Fourier System
Polynomial Approximation: The Fourier System Charles B. I. Chilaka CASA Seminar 17th October, 2007 Outline 1 Introduction and problem formulation 2 The continuous Fourier expansion 3 The discrete Fourier
More informationStability of Kernel Based Interpolation
Stability of Kernel Based Interpolation Stefano De Marchi Department of Computer Science, University of Verona (Italy) Robert Schaback Institut für Numerische und Angewandte Mathematik, University of Göttingen
More informationDecomposition of Riesz frames and wavelets into a finite union of linearly independent sets
Decomposition of Riesz frames and wavelets into a finite union of linearly independent sets Ole Christensen, Alexander M. Lindner Abstract We characterize Riesz frames and prove that every Riesz frame
More informationResearch Article A Note on the Solutions of the Van der Pol and Duffing Equations Using a Linearisation Method
Mathematical Problems in Engineering Volume 1, Article ID 693453, 1 pages doi:11155/1/693453 Research Article A Note on the Solutions of the Van der Pol and Duffing Equations Using a Linearisation Method
More informationThe Gibbs Phenomenon
The Gibbs Phenomenon Eli Dean & Britton Girard Math 572, Fall 2013 December 4, 2013 1. Introduction 1 S N f(x) := f(n)e inx := 1 2π n N n N denotes the partial Fourier sum for f at the point x f(x)e inx
More informationSpectral methods for compressible reactive flows
C. R. Mecanique 333 (2005) 3 16 http://france.elsevier.com/direct/cras2b/ High-Order Methods for the Numerical Simulation of Vortical and Turbulent Flows Spectral methods for compressible reactive flows
More informationInfinite series, improper integrals, and Taylor series
Chapter Infinite series, improper integrals, and Taylor series. Determine which of the following sequences converge or diverge (a) {e n } (b) {2 n } (c) {ne 2n } (d) { 2 n } (e) {n } (f) {ln(n)} 2.2 Which
More informationFourier and Partial Differential Equations
Chapter 5 Fourier and Partial Differential Equations 5.1 Fourier MATH 294 SPRING 1982 FINAL # 5 5.1.1 Consider the function 2x, 0 x 1. a) Sketch the odd extension of this function on 1 x 1. b) Expand the
More informationWavelets Marialuce Graziadei
Wavelets Marialuce Graziadei 1. A brief summary 2. Vanishing moments 3. 2D-wavelets 4. Compression 5. De-noising 1 1. A brief summary φ(t): scaling function. For φ the 2-scale relation hold φ(t) = p k
More informationExamples of the Fourier Theorem (Sect. 10.3). The Fourier Theorem: Continuous case.
s of the Fourier Theorem (Sect. 1.3. The Fourier Theorem: Continuous case. : Using the Fourier Theorem. The Fourier Theorem: Piecewise continuous case. : Using the Fourier Theorem. The Fourier Theorem:
More informationLectures notes. Rheology and Fluid Dynamics
ÉC O L E P O L Y T E C H N IQ U E FÉ DÉR A L E D E L A U S A N N E Christophe Ancey Laboratoire hydraulique environnementale (LHE) École Polytechnique Fédérale de Lausanne Écublens CH-05 Lausanne Lectures
More informationPAijpam.eu NUMERICAL SOLUTION OF WAVE EQUATION USING HAAR WAVELET Inderdeep Singh 1, Sangeeta Arora 2, Sheo Kumar 3
International Journal of Pure and Applied Mathematics Volume 98 No. 4 25, 457-469 ISSN: 3-88 (printed version); ISSN: 34-3395 (on-line version) url: http://www.ijpam.eu doi: http://dx.doi.org/.2732/ijpam.v98i4.4
More informationMath 121A: Homework 6 solutions
Math A: Homework 6 solutions. (a) The coefficients of the Fourier sine series are given by b n = π f (x) sin nx dx = x(π x) sin nx dx π = (π x) cos nx dx nπ nπ [x(π x) cos nx]π = n ( )(sin nx) dx + π n
More informationIntroduction to Discrete-Time Wavelet Transform
Introduction to Discrete-Time Wavelet Transform Selin Aviyente Department of Electrical and Computer Engineering Michigan State University February 9, 2010 Definition of a Wavelet A wave is usually defined
More informationFunctions: A Fourier Approach. Universitat Rostock. Germany. Dedicated to Prof. L. Berg on the occasion of his 65th birthday.
Approximation Properties of Multi{Scaling Functions: A Fourier Approach Gerlind Plona Fachbereich Mathemati Universitat Rostoc 1851 Rostoc Germany Dedicated to Prof. L. Berg on the occasion of his 65th
More informationHierarchical Reconstruction with up to Second Degree Remainder for Solving Nonlinear Conservation Laws
Hierarchical Reconstruction with up to Second Degree Remainder for Solving Nonlinear Conservation Laws Dedicated to Todd F. Dupont on the occasion of his 65th birthday Yingjie Liu, Chi-Wang Shu and Zhiliang
More information1. Numerical Analysis of Spectral Methods/Theory and Applications, with S. Orszag, CBMS-SIAM No. 26, 1977, 170 pages.
DAVID GOTTLIEB Ford Foundation Professor Professor of Applied Mathematics Brown University Publications Updated February 2007 Books 1. Numerical Analysis of Spectral Methods/Theory and Applications, with
More informationFourier Series. (Com S 477/577 Notes) Yan-Bin Jia. Nov 29, 2016
Fourier Series (Com S 477/577 otes) Yan-Bin Jia ov 9, 016 1 Introduction Many functions in nature are periodic, that is, f(x+τ) = f(x), for some fixed τ, which is called the period of f. Though function
More informationA MULTIGRID ALGORITHM FOR. Richard E. Ewing and Jian Shen. Institute for Scientic Computation. Texas A&M University. College Station, Texas SUMMARY
A MULTIGRID ALGORITHM FOR THE CELL-CENTERED FINITE DIFFERENCE SCHEME Richard E. Ewing and Jian Shen Institute for Scientic Computation Texas A&M University College Station, Texas SUMMARY In this article,
More informationAP Calculus Testbank (Chapter 9) (Mr. Surowski)
AP Calculus Testbank (Chapter 9) (Mr. Surowski) Part I. Multiple-Choice Questions n 1 1. The series will converge, provided that n 1+p + n + 1 (A) p > 1 (B) p > 2 (C) p >.5 (D) p 0 2. The series
More informationA Chebyshev Polynomial Rate-of-Convergence Theorem for Stieltjes Functions
mathematics of computation volume 39, number 159 july 1982, pages 201-206 A Chebyshev Polynomial Rate-of-Convergence Theorem for Stieltjes Functions By John P. Boyd Abstract. The theorem proved here extends
More informationA gradient recovery method based on an oblique projection and boundary modification
ANZIAM J. 58 (CTAC2016) pp.c34 C45, 2017 C34 A gradient recovery method based on an oblique projection and boundary modification M. Ilyas 1 B. P. Lamichhane 2 M. H. Meylan 3 (Received 24 January 2017;
More informationMA 201: Differentiation and Integration of Fourier Series Applications of Fourier Series Lecture - 10
MA 201: Differentiation and Integration of Fourier Series Applications of Fourier Series ecture - 10 Fourier Series: Orthogonal Sets We begin our treatment with some observations: For m,n = 1,2,3,... cos
More informationTwo Channel Subband Coding
Two Channel Subband Coding H1 H1 H0 H0 Figure 1: Two channel subband coding. In two channel subband coding A signal is convolved with a highpass filter h 1 and a lowpass filter h 0. The two halfband signals
More informationDETECTION OF EDGES IN SPECTRAL DATA II. NONLINEAR ENHANCEMENT
SIAM J. UMER. AAL. Vol. 38, o. 4, pp. 389 48 c 2 Society for Industrial and Applied Mathematics DETECTIO OF EDGES I SPECTRAL DATA II. OLIEAR EHACEMET AE GELB AD EITA TADMOR Abstract. We discuss a general
More informationA general theory of discrete ltering. for LES in complex geometry. By Oleg V. Vasilyev AND Thomas S. Lund
Center for Turbulence Research Annual Research Briefs 997 67 A general theory of discrete ltering for ES in complex geometry By Oleg V. Vasilyev AND Thomas S. und. Motivation and objectives In large eddy
More informationCOMPACTLY SUPPORTED ORTHONORMAL COMPLEX WAVELETS WITH DILATION 4 AND SYMMETRY
COMPACTLY SUPPORTED ORTHONORMAL COMPLEX WAVELETS WITH DILATION 4 AND SYMMETRY BIN HAN AND HUI JI Abstract. In this paper, we provide a family of compactly supported orthonormal complex wavelets with dilation
More informationQuadrature Prefilters for the Discrete Wavelet Transform. Bruce R. Johnson. James L. Kinsey. Abstract
Quadrature Prefilters for the Discrete Wavelet Transform Bruce R. Johnson James L. Kinsey Abstract Discrepancies between the Discrete Wavelet Transform and the coefficients of the Wavelet Series are known
More informationA Novel Fast Computing Method for Framelet Coefficients
American Journal of Applied Sciences 5 (11): 15-157, 008 ISSN 1546-939 008 Science Publications A Novel Fast Computing Method for Framelet Coefficients Hadeel N. Al-Taai Department of Electrical and Electronic
More informationTHE S-ELEMENTARY WAVELETS ARE PATH-CONNECTED. D. M. Speegle
THE S-ELEMENTARY WAVELETS ARE PATH-CONNECTED D. M. Speegle Abstract. A construction of wavelet sets containing certain subsets of R is given. The construction is then modified to yield a continuous dependence
More information