DISCRETE DIFFERENTIAL OPERATORS IN MULTIDIMENSIONAL HAAR WAVELET SPACES
|
|
- Darren Oliver Thornton
- 5 years ago
- Views:
Transcription
1 IJMMS 2004: PII. S Hindawi Publishing Corp. DISCRETE DIFFERENTIAL OPERATORS IN MULTIDIMENSIONAL HAAR WAVELET SPACES CARLO CATTANI and LUIS M. SÁNCHEZ RUIZ Received 30 July 2003 We consider a class of discrete differential operators acting on multidimensional Haar wavelet basis with the aim of finding wavelet approximate solutions of partial differential problems. Although these operators depend on the interpolating method used for the Haar wavelets regularization and the scale dimension space they can be easily used to define the space of approximate wavelet solutions Mathematics Subject Classification: 35A35 42C40 41A Introduction. Advantage has extensively taken upon wavelets in order to analyze a number of different applications such as image processing signal detection geophysics medicine or turbulent flows [ ]. More mathematically focussed differential equations and even nonlinear problems have also been studied with this cargo of possibilities that wavelets have brought out into scene [ ]. In this paper we study a Cauchy problem with the PDE u t = Lu where the unknown function utx L 2 [01] for each t [0T] x [01] d d N and L is a nonnecessarily linear partial differential operator whose derivatives act on the space variable. Situations like this happen quite commonly in evolution processes such as fluid dynamics heat transfer elasticity traffic flow shock propagation and so on [5 12]. According to the nature and properties of L different approaches have been developed and even simple difference approximations may provide good results [3 11]. However in the presence of discontinuities some difficulties arise. The development of Haar wavelets with its extreme suitability for dealing with experimental problems having piecewise constant functions as initial conditions or functions with sharp discontinuities such as those that appear in the Riemann problem or within smooth environments seems to provide a natural scheme to study these problems and wavelet approximate solutions have been obtained by using some regular mostly Daubechies bases of wavelets [11]. In the following we propose an algorithm based upon Haar wavelet basis in order to define a discrete operator that maps piecewise constant functions into piecewise constant functions which projects the continuous operator L therein and enables to get wavelet approximate solutions in L 2 [01] just working with Haar wavelets. 2. Haar wavelet basis. Haar wavelets are piecewise constant functions with compact support and finite jumps at their extremal points. The drawback of not being smooth enough has been overcome by a smoothing process based on spline interpolation [3].
2 2348 C. CATTANI AND L. M. SÁNCHEZ RUIZ Haar functions ϕ n k x := 2n/2 ϕ2 n x k are generated from the scaling function ϕx which is given by the characteristic function χ [01]. Then the Haar family of wavelets ψ n k x :={2n/2 ψ2 n x k} kn Z fulfills [ k 2 n/2 x 2 k+1/2 n 2 n ψ n k x = [ k+1/2 2 n/2 x k+1 2 n 2 n 0 elsewhere ψ n k x 2 = 1 and the multiresolution axioms [7] as well as the recursive conditions n+1/2 ϕ n k x = ϕn+1 k x+ϕ n+1 k+1 x 2n+1/2 ψ n k x = ϕn+1 k x ϕ n+1 k+1 x kn Z. 2.2 Haar wavelets form a complete orthonormal system for the finite-energy L 2 R- functions [7] so that L 2 R = W n = V q W j q Z n Z j q V n+1 = V n W n. 2.3 Here for each fixed n Z V n is the subspace { V n := f L 2 R : f = nk Z c n k χ D n k cn k = constant } 2.4 where D n k = [k/2n k+1/2 n and W n are orthogonal subspaces in L 2 R the so-called wavelet spaces. As usual we consider the scalar product + fg := fxgxdx fg L 2 R. 2.5 According to the first equation of 2.3 we have that fx= nk Z β n k ψn k x. Ifwe fix a resolution value N = 2 M M Nin2.3 then an approximation of the L 2 R-space is obtained for example by L 2 N R V 0 n=0 W n thatis M 1 fx π N fx:= α0 0 + n=0 2 n 1 k=0 β n k ψn k x 2.6 π N being a projection operator into V N so that π N : L 2 R V N. In general the coefficients αn k β n k are given by α n k = fϕ n k β n k = fψ n k. 2.7 The dyadic discretization is the operator N : L 2 R K N l 2 where K stands for the complex or real field and N fx = f N = f 0 f 1...f N 1 T K N with f k := fk/n 1 0 k N 1.
3 DIFFERENTIAL OPERATORS IN MULTIDIMENSIONAL HAAR WAVELET The above notation may be simplified if we denote β N := α0 0 β0 0 β1 0...βM and Ψ N x := ϕxψ0 0 xψ1 0 x...ψm 1 2 x so that 1 π N fx= β N ΨN T x fn = β N N ΨN T x. 2.8 If P 2 i denotes the permutation or shuffle matrix cf. [3 4 7] and I is the identity matrix then the discrete Haar wavelet transform is the one-to-one linear operator of K N onto K N such that [ M f N = β T ] N = P2 k I 2 M 2 k H2 k I 2 M 2 k f N. 2.9 k=1 The matrices H 2 i have size 2 i 2 i and are the direct sum of the lattice coefficients of the system /2 2 1/2 2 1/2 2 1/2 H 2 i := 2 1/2 2 1/2 2 1/2 2 1/2 where for 1 j 2 i 1 and O pq standing for the p q null matrix 1 2 i O 2j 22j 2 O 2j 22 O 2j /2 2 1/2 i 2j 2 1/2 2 1/2 = 2 1/2 2 1/2 O 22j 2 j 2 1/2 2 1/2 O 22 i 2j O 2 i 2j2j 2 O 2 i 2j2 O 2 i 2j2 i 2j Thus having in mind that V N+1 K N by identifying α0 0 + M 1 2 n 1 n=0 k=0 β n k ψn k x with {α0 0 βn k }n=0...m 1 k= the projection operator π N : L 2 R V N+1 may be factorized by means of the following diagram: L 2 R N K N π N 2.12 V N+1 K N The operator π N maps each f L 2 R into a piecewise constant function π N f V N+1 being easy to check lim N π N f = f pointwise cf. [7]. Note that the inverse discrete Haar transform 1 of 2.9 is univocally defined [7] but the inverse N := N 1 of the operator N is not because there are infinitely many functions interpolating a finite set of values. 3. Haar series regularization. In what follows we will consider the interpolating operator p : K N C p 1 Ω which maps the vector f N = f i N 1 i=0 based on the dyadic nodes of Ω R into the p-differentiable function sx := p f N where C 1 Ω stands for the piecewise continuous functions defined on Ω. Hencesx is an interpolating function of order p through i/n 1f i i = 0...N 1.
4 2350 C. CATTANI AND L. M. SÁNCHEZ RUIZ Definition 3.1 interpolating p-space. For a given resolution N the interpolating p-space σ p is the space σ p L 2 Ω of interpolating functions with nonredundant fixed conditions at the boundary { σ p := sx = p f N with p := p 1 = N f N K N}. 3.1 There follows that fixing the resolution N and the interpolating method each function in σ p corresponds one-to-one to a discrete set f N making the following diagrams commutative: L 2 Ω σ p N K N L 2 Ω σ p p K N π N π N V N+1 V N+1 Theorem 3.2. For a fixed resolution N and a given interpolating method the inverse π N of the projection operator π N is univocally defined when acting on the σ p space functions by means of π N = p 1 π N = N. 3.3 Proof. Note that N p = p N = I and 1 is univocally defined [7]. Definition 3.3. The q-derivative of a p-interpolating function 0 q p atagiven dimension N is the linear operator δ pq : K N V N+1 such that dq δ pq := π N dx q p. 3.4 This definition provides an algorithm that maps piecewise constant functions into piecewise constant functions. If we look at it we may note that firstly it smoothens piecewise constant functions with a suitable interpolation; secondly it finds out the derivative of the interpolating function and then this derivative is transformed into a piecewise constant function by the discrete Haar wavelet transform. From [7 page 13] and since d q /dx q p N f belongs to L 2 Ω for each f L 2 Ω it follows that d q /dx q p N f δ pq N f 2 L 2 2 2M 1 2 µω K + ε where µ stands for the Lebesgue measure K = M 1 k=0 m 0kN 2 + m 1kN 2 1/2 withm 0kN and m 1kN being the average value of d q /dx q p N f over the first and the second half of D N k respectively and ε may be taken as small as desired since any function in L 2 Ω may be approximated by functions with compact support which are piecewise constant. It is easy to see that if f = M 1 k=0 f k χ Dk g = M 1 k=0 g k χ Dk V N where f k g k are constants then there are unique c k such that f k = g k +c k. Hence there exists a matrix A pq depending on N that makes δ pq N f = A pq N f δ pq f N = A pq β N N Ψ T N x. 3.5
5 DIFFERENTIAL OPERATORS IN MULTIDIMENSIONAL HAAR WAVELET At the end of this paper we provide an appendix with some of these matrices A for different spline interpolations. 4. Multidimensional wavelet solution. In a d-dimensional variable space x = x 1...x d we restrict ourselves to considering the same dyadic discretization on each coordinate of x so that the natural generalization of the dyadic discretization operator N : L 2 R d d k=1 K N N = 2 M is given by N i1 fx := f N 1 i 2 N 1... i d N 1 0 i j N 1 1 j d. 4.1 By using tensor products we define V N x := d i=1 V N x i where the orthonormal wavelet basis is {ψ n k x 1 ψ n k x d} 0 k N 1cf.[17]. Thus V N x is the space of the piecewise constant functions on the d-dimensional product D N k DN h 0 h k 2 N 1 which can be represented as π N fx := + 1 i d α i ϕ x i + 0 k i 2 n i 1 0 n i M 1 i d 1 ij d i j 0 k i 2 n i 1 0 n i M 1 i d 1 i d β n i k i ψ n i k i xi ω n i k i j ϕ x j ψ n i k i xi 4.2 in short π N fx := 0 i j M 1 1 j d β i 1i 2 i d N N ψ N x i1 ψ N x id. Then the discrete d-dimensional Haar wavelet transform is the one-to-one linear operator d = d of K N d onto K N d and the discrete derivative of N fx is the tensor product of the corresponding 1-dimensional discrete derivative. Hence in the 2-dimensional case this leads to a matrix A pq which in fact is the same as the one of the 1-dimensional case such that δ pq i1 x f δ pq y f N 1 i 2 N 1 i1 N 1 i 2 N 1 = A pq f = A pq i1 N 1 i 2 0 i 2 N 1 N 1 0 i 1 N proving that if x 1 = i 1 /N 1 x 2 = i 2 /N 1 δ pq x N f x 1 x 2 = A pq β N N Ψ T N x1 N Ψ T N x2 4.4 and analogously for δ pq y. Taking this into account we obtain the following result.
6 2352 C. CATTANI AND L. M. SÁNCHEZ RUIZ Proposition 4.1. problem An approximate wavelet solution ũtx=π N utx of the Cauchy u t = Lu u0x = u 0 x x [01] d t [0T] 4.5 is obtained by taking ũtx := 0 i j M 1 1 j d β i 1i 2 i d N t N ψ N xi1 ψn xid 4.6 where the wavelet coefficients are explicit functions of t. Proof. From 4.5 it follows that π N u t = π N Lu. HenceifL δ denotes the operator acting on piecewise functions obtained by replacing the partial derivatives that appear on L with the corresponding δ pq i j N for example if L = L / x 2 / x y then L δ = Lδ p1 x N δ p1 y δ p1 x N it happens that the projection in V N x must verify ũ t tx = L δ ũtx and 0 i j M 1 1 j d β i 1i 2 i d N 0 N ψ N xi1 ψn xid = N fx. 4.7 Orthonormality of the scaling functions ϕ and ψ n k gives an easily solved first-order ordinary system in the unknown wavelet coefficients. Note that due to the definition of Haar wavelets the boundary conditions are automatically fulfilled. 5. Spline Haar solution of a parabolic equation. In this section we will compare the Haar regularized solution of the following 1-dimensional parabolic problem with the exact solution and the one obtained with Crank-Nicholson s method. Example 5.1. Compareattime0.04 and 0.08 the above-mentioned methods in the heat equation without exchanges at the boundary and a periodic function as initial function given by u t = u xx 0 x 1 t 0 ux0 = sinπx u0t= u1t= In Figure 5.1 the exact solution uxt = sinπxe π2t at t = 0.01 and the one obtained by the Haar regularized method based on 64 nodes are shown as functions of the variable distance x. Table 5.1 compares the set of values that we obtain HR at time t = 0.04 t = 0.08 with the exact ones E and those obtained by Crank-Nicholson s method CN.
7 DIFFERENTIAL OPERATORS IN MULTIDIMENSIONAL HAAR WAVELET u heat x distance 1 Figure 5.1. Exact and Haar regularized solution with 64 nodes. Table 5.1 x t = 0.04 t = 0.08 E HR CN E HR CN Annex. Finally we enclose a list of matrices A pq announced at the end of Section 3 which satisfy δ pq N f = A pq N f. We fix our attention on the case p = 3 q = 1. For this and different resolutions N = 2 M M N the matrices A pq are computed according to the following steps: 1 discretize the interval [01] at the dyadic points k/n 1 0 k N 1 2 consider the points whose first coordinates are the dyadic points and whose second coordinates are arbitrary values a 0 a 1...a N 1 that is 0a 0 1/N 1a 1...1a N 1 3 build up the cubic spline that goes through the aforementioned set of points 4 compute the first derivative of this spline and evaluate it at the dyadic points 5 assuming the above evaluations generate the points 0a 0 1/N 1a a N 1 thena31 is the N N matrix that satisfies a 0 a 1.. a N 1 = A 31 a 0 a 1. a N
8 2354 C. CATTANI AND L. M. SÁNCHEZ RUIZ With a bound error of 10 5 routine calculations make A 31 become N = 2 1 = N = 2 2 = N = 2 3 = Acknowledgment. The second author is supported by Generalitat Valenciana. References [1] E. Bacry S. Mallat and G. Papanicolaou A wavelet based space-time adaptive numerical method for partial differential equations RAIRO Modél. Math. Anal. Numér no [2] W. Cai and J. Wang Adaptive multiresolution collocation methods for initial-boundary value problems of nonlinear PDEs SIAM J. Numer. Anal no [3] C. Cattani Haar wavelet splines J. Interdiscip. Math no [4] Haar wavelet-based technique for sharp jumps classification Math. Comput. Modelling no [5] R. Courant and K. O. Friedrichs Supersonic Flow and Shock Waves Interscience Publishers New York [6] W. Dahmen Wavelet methods for PDEs some recent developments J.Comput.Appl.Math no [7] I. Daubechies Ten Lectures on Wavelets CBMS-NSF Regional Conference Series in Applied Mathematics vol. 61 Society for Industrial and Applied Mathematics Pennsylvania [8] L. Debnath Wavelet transforms and their applications Proc. Indian Nat. Sci. Acad. Part A no [9] R. A. DeVore B. Jawerth and B. J. Lucier Image compression through wavelet transform coding IEEE Trans. Inform. Theory no [10] D. L. Donoho Nonlinear wavelet methods for recovery of signals densities and spectra from indirect and noisy data Different Perspectives on Wavelets San Antonio Tex
9 DIFFERENTIAL OPERATORS IN MULTIDIMENSIONAL HAAR WAVELET Proc. Sympos. Appl. Math. vol. 47 American Mathematical Society Rhode Island 1993 pp [11] S. Lazaar P. Ponenti J. Liandrat and P. Tchamitchian Wavelet algorithms for numerical resolution of partial differential equations Comput. Methods Appl. Mech. Engrg no [12] R. J. LeVeque Numerical Methods for Conservation Laws 2nd ed. Lectures in Mathematics ETH Zürich Birkhäuser Basel [13] P. J. Oonincx A wavelet method for detecting S-waves in seismic data Comput. Geosci no [14] W. Schempp Wavelet modelling of high resolution radar imaging and clinical magnetic resonance tomography Results Math no [15] J. Segman and Y. Y. Zeevi Image analysis by wavelet-type transforms: group-theoretic approach J. Math. Imaging Vision no [16] O. V. Vasilyev and S. Paolucci A fast adaptive wavelet collocation algorithm for multidimensional PDEs J. Comput. Phys no [17] P. Wojtaszczyk A Mathematical Introduction to Wavelets London Mathematical Society Student Texts vol. 37 Cambridge University Press Cambridge [18] V. Zimin and F. Hussain Wavelet based model for small-scale turbulence Phys. Fluids no Carlo Cattani: Dipartimento di Scienze Farmaceutiche DiFARMA Università degli Studi di Salerno Via Ponte Don Melillo Invariante 11/C Fisciano Salerno Italy address: ccattani@unisa.it Luis M. Sánchez Ruiz: Departamento de Matemática Aplicada Escuela Técnica Superior de Ingeniería de Diseño ETSID Universidad Politécnica de Valencia Valencia Spain address: lmsr@mat.upv.es
Biorthogonal Spline Type Wavelets
PERGAMON Computers and Mathematics with Applications 0 (00 1 0 www.elsevier.com/locate/camwa Biorthogonal Spline Type Wavelets Tian-Xiao He Department of Mathematics and Computer Science Illinois Wesleyan
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 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 informationA Very Brief Introduction to Conservation Laws
A Very Brief Introduction to Wen Shen Department of Mathematics, Penn State University Summer REU Tutorial, May 2013 Summer REU Tutorial, May 2013 1 / The derivation of conservation laws A conservation
More informationPOINT VALUES AND NORMALIZATION OF TWO-DIRECTION MULTIWAVELETS AND THEIR DERIVATIVES
November 1, 1 POINT VALUES AND NORMALIZATION OF TWO-DIRECTION MULTIWAVELETS AND THEIR DERIVATIVES FRITZ KEINERT AND SOON-GEOL KWON,1 Abstract Two-direction multiscaling functions φ and two-direction multiwavelets
More informationWe have to prove now that (3.38) defines an orthonormal wavelet. It belongs to W 0 by Lemma and (3.55) with j = 1. We can write any f W 1 as
88 CHAPTER 3. WAVELETS AND APPLICATIONS We have to prove now that (3.38) defines an orthonormal wavelet. It belongs to W 0 by Lemma 3..7 and (3.55) with j =. We can write any f W as (3.58) f(ξ) = p(2ξ)ν(2ξ)
More informationLittlewood Paley Spline Wavelets
Proceedings of the 6th WSEAS International Conference on Wavelet Analysis & Multirate Systems, Bucharest, Romania, October 6-8, 26 5 Littlewood Paley Spline Wavelets E. SERRANO and C.E. D ATTELLIS Escuela
More informationHow much should we rely on Besov spaces as a framework for the mathematical study of images?
How much should we rely on Besov spaces as a framework for the mathematical study of images? C. Sinan Güntürk Princeton University, Program in Applied and Computational Mathematics Abstract Relations between
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 informationFinite Volume Schemes: an introduction
Finite Volume Schemes: an introduction First lecture Annamaria Mazzia Dipartimento di Metodi e Modelli Matematici per le Scienze Applicate Università di Padova mazzia@dmsa.unipd.it Scuola di dottorato
More informationON SAMPLING RELATED PROPERTIES OF B-SPLINE RIESZ SEQUENCES
ON SAMPLING RELATED PROPERTIES OF B-SPLINE RIESZ SEQUENCES SHIDONG LI, ZHENGQING TONG AND DUNYAN YAN Abstract. For B-splineRiesz sequencesubspacesx span{β k ( n) : n Z}, thereis an exact sampling formula
More informationWavelets and Image Compression. Bradley J. Lucier
Wavelets and Image Compression Bradley J. Lucier Abstract. In this paper we present certain results about the compression of images using wavelets. We concentrate on the simplest case of the Haar decomposition
More informationrecent developments of approximation theory and greedy algorithms
recent developments of approximation theory and greedy algorithms Peter Binev Department of Mathematics and Interdisciplinary Mathematics Institute University of South Carolina Reduced Order Modeling in
More informationHigh-resolution finite volume methods for hyperbolic PDEs on manifolds
High-resolution finite volume methods for hyperbolic PDEs on manifolds Randall J. LeVeque Department of Applied Mathematics University of Washington Supported in part by NSF, DOE Overview High-resolution
More informationAPPLICATIONS OF THE THICK DISTRIBUTIONAL CALCULUS
Novi Sad J. Math. Vol. 44, No. 2, 2014, 121-135 APPLICATIONS OF THE THICK DISTRIBUTIONAL CALCULUS Ricardo Estrada 1 and Yunyun Yang 2 Abstract. We give several applications of the thick distributional
More informationWavelets and regularization of the Cauchy problem for the Laplace equation
J. Math. Anal. Appl. 338 008440 1447 www.elsevier.com/locate/jmaa Wavelets and regularization of the Cauchy problem for the Laplace equation Chun-Yu Qiu, Chu-Li Fu School of Mathematics and Statistics,
More informationApplied Math Qualifying Exam 11 October Instructions: Work 2 out of 3 problems in each of the 3 parts for a total of 6 problems.
Printed Name: Signature: Applied Math Qualifying Exam 11 October 2014 Instructions: Work 2 out of 3 problems in each of the 3 parts for a total of 6 problems. 2 Part 1 (1) Let Ω be an open subset of R
More informationA NOTE ON WAVELET EXPANSIONS FOR DYADIC BMO FUNCTIONS IN SPACES OF HOMOGENEOUS TYPE
EVISTA DE LA UNIÓN MATEMÁTICA AGENTINA Vol. 59, No., 208, Pages 57 72 Published online: July 3, 207 A NOTE ON WAVELET EXPANSIONS FO DYADIC BMO FUNCTIONS IN SPACES OF HOMOGENEOUS TYPE AQUEL CESCIMBENI AND
More informationWavelets and modular inequalities in variable L p spaces
Wavelets and modular inequalities in variable L p spaces Mitsuo Izuki July 14, 2007 Abstract The aim of this paper is to characterize variable L p spaces L p( ) (R n ) using wavelets with proper smoothness
More informationResearch Article Fractional Calculus and Shannon Wavelet
Mathematical Problems in Engineering Volume 2012, Article ID 502812, 26 pages doi:10.1155/2012/502812 Research Article Fractional Calculus and Shannon Wavelet Carlo Cattani Department of Mathematics, University
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 informationMLISP: Machine Learning in Signal Processing Spring Lecture 8-9 May 4-7
MLISP: Machine Learning in Signal Processing Spring 2018 Prof. Veniamin Morgenshtern Lecture 8-9 May 4-7 Scribe: Mohamed Solomon Agenda 1. Wavelets: beyond smoothness 2. A problem with Fourier transform
More information446 SCIENCE IN CHINA (Series F) Vol. 46 introduced in refs. [6, ]. Based on this inequality, we add normalization condition, symmetric conditions and
Vol. 46 No. 6 SCIENCE IN CHINA (Series F) December 003 Construction for a class of smooth wavelet tight frames PENG Lizhong (Λ Π) & WANG Haihui (Ξ ) LMAM, School of Mathematical Sciences, Peking University,
More informationSparse Multidimensional Representation using Shearlets
Sparse Multidimensional Representation using Shearlets Demetrio Labate a, Wang-Q Lim b, Gitta Kutyniok c and Guido Weiss b, a Department of Mathematics, North Carolina State University, Campus Box 8205,
More informationNon-linear Wave Propagation and Non-Equilibrium Thermodynamics - Part 3
Non-linear Wave Propagation and Non-Equilibrium Thermodynamics - Part 3 Tommaso Ruggeri Department of Mathematics and Research Center of Applied Mathematics University of Bologna January 21, 2017 ommaso
More informationGalerkin method for the numerical solution of the RLW equation using quintic B-splines
Journal of Computational and Applied Mathematics 19 (26) 532 547 www.elsevier.com/locate/cam Galerkin method for the numerical solution of the RLW equation using quintic B-splines İdris Dağ a,, Bülent
More informationCHAPTER III HAAR WAVELET METHOD FOR SOLVING FISHER S EQUATION
CHAPTER III HAAR WAVELET METHOD FOR SOLVING FISHER S EQUATION A version of this chapter has been published as Haar Wavelet Method for solving Fisher s equation, Appl. Math.Comput.,(ELSEVIER) 211 (2009)
More informationA First Course in Wavelets with Fourier Analysis
* A First Course in Wavelets with Fourier Analysis Albert Boggess Francis J. Narcowich Texas A& M University, Texas PRENTICE HALL, Upper Saddle River, NJ 07458 Contents Preface Acknowledgments xi xix 0
More informationA Riesz basis of wavelets and its dual with quintic deficient splines
Note di Matematica 25, n. 1, 2005/2006, 55 62. A Riesz basis of wavelets and its dual with quintic deficient splines F. Bastin Department of Mathematics B37, University of Liège, B-4000 Liège, Belgium
More informationLet p 2 ( t), (2 t k), we have the scaling relation,
Multiresolution Analysis and Daubechies N Wavelet We have discussed decomposing a signal into its Haar wavelet components of varying frequencies. The Haar wavelet scheme relied on two functions: the Haar
More informationDeep Learning: Approximation of Functions by Composition
Deep Learning: Approximation of Functions by Composition Zuowei Shen Department of Mathematics National University of Singapore Outline 1 A brief introduction of approximation theory 2 Deep learning: approximation
More informationNontechnical introduction to wavelets Continuous wavelet transforms Fourier versus wavelets - examples
Nontechnical introduction to wavelets Continuous wavelet transforms Fourier versus wavelets - examples A major part of economic time series analysis is done in the time or frequency domain separately.
More informationUNIVERSITY OF WISCONSIN-MADISON CENTER FOR THE MATHEMATICAL SCIENCES. Tight compactly supported wavelet frames of arbitrarily high smoothness
UNIVERSITY OF WISCONSIN-MADISON CENTER FOR THE MATHEMATICAL SCIENCES Tight compactly supported wavelet frames of arbitrarily high smoothness Karlheinz Gröchenig Amos Ron Department of Mathematics U-9 University
More informationWavelet Bases of the Interval: A New Approach
Int. Journal of Math. Analysis, Vol. 1, 2007, no. 21, 1019-1030 Wavelet Bases of the Interval: A New Approach Khaled Melkemi Department of Mathematics University of Biskra, Algeria kmelkemi@yahoo.fr Zouhir
More informationApplications of Polyspline Wavelets to Astronomical Image Analysis
VIRTUAL OBSERVATORY: Plate Content Digitization, Archive Mining & Image Sequence Processing edited by M. Tsvetkov, V. Golev, F. Murtagh, and R. Molina, Heron Press, Sofia, 25 Applications of Polyspline
More informationANALYSIS AND NUMERICAL METHODS FOR SOME CRACK PROBLEMS
INTERNATIONAL JOURNAL OF NUMERICAL ANALYSIS AND MODELING, SERIES B Volume 2, Number 2-3, Pages 155 166 c 2011 Institute for Scientific Computing and Information ANALYSIS AND NUMERICAL METHODS FOR SOME
More informationLecture No 1 Introduction to Diffusion equations The heat equat
Lecture No 1 Introduction to Diffusion equations The heat equation Columbia University IAS summer program June, 2009 Outline of the lectures We will discuss some basic models of diffusion equations and
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 informationBIHARMONIC WAVE MAPS INTO SPHERES
BIHARMONIC WAVE MAPS INTO SPHERES SEBASTIAN HERR, TOBIAS LAMM, AND ROLAND SCHNAUBELT Abstract. A global weak solution of the biharmonic wave map equation in the energy space for spherical targets is constructed.
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 informationVectors in Function Spaces
Jim Lambers MAT 66 Spring Semester 15-16 Lecture 18 Notes These notes correspond to Section 6.3 in the text. Vectors in Function Spaces We begin with some necessary terminology. A vector space V, also
More informationResearch Article A Shannon-Runge-Kutta-Gill Method for Convection-Diffusion Equations
Mathematical Problems in Engineering Volume 203, Article ID 63734, 5 pages http://dx.doi.org/0.55/203/63734 Research Article A Shannon-Runge-Kutta-Gill Method for Convection-Diffusion Equations Xiaoming
More informationApproximation of Integrable Functions by Wavelet Expansions
Results Math 72 27, 23 2 c 26 The Authors. This article is published with open access at Springerlink.com 422-6383/7/323-9 published online October 25, 26 DOI.7/s25-6-64-z Results in Mathematics Approximation
More informationRafał Kapica, Janusz Morawiec CONTINUOUS SOLUTIONS OF ITERATIVE EQUATIONS OF INFINITE ORDER
Opuscula Mathematica Vol. 29 No. 2 2009 Rafał Kapica, Janusz Morawiec CONTINUOUS SOLUTIONS OF ITERATIVE EQUATIONS OF INFINITE ORDER Abstract. Given a probability space (, A, P ) and a complete separable
More informationQuintic B-Spline Galerkin Method for Numerical Solutions of the Burgers Equation
5 1 July 24, Antalya, Turkey Dynamical Systems and Applications, Proceedings, pp. 295 39 Quintic B-Spline Galerkin Method for Numerical Solutions of the Burgers Equation İdris Dağ 1,BülentSaka 2 and Ahmet
More informationModule 4 MULTI- RESOLUTION ANALYSIS. Version 2 ECE IIT, Kharagpur
Module 4 MULTI- RESOLUTION ANALYSIS Lesson Theory of Wavelets Instructional Objectives At the end of this lesson, the students should be able to:. Explain the space-frequency localization problem in sinusoidal
More informationEXAMPLES OF REFINABLE COMPONENTWISE POLYNOMIALS
EXAMPLES OF REFINABLE COMPONENTWISE POLYNOMIALS NING BI, BIN HAN, AND ZUOWEI SHEN Abstract. This short note presents four examples of compactly supported symmetric refinable componentwise polynomial functions:
More informationIntroduction. J.M. Burgers Center Graduate Course CFD I January Least-Squares Spectral Element Methods
Introduction In this workshop we will introduce you to the least-squares spectral element method. As you can see from the lecture notes, this method is a combination of the weak formulation derived from
More information( nonlinear constraints)
Wavelet Design & Applications Basic requirements: Admissibility (single constraint) Orthogonality ( nonlinear constraints) Sparse Representation Smooth functions well approx. by Fourier High-frequency
More informationQUADRATURES INVOLVING POLYNOMIALS AND DAUBECHIES' WAVELETS *
QUADRATURES INVOLVING POLYNOMIALS AND DAUBECHIES' WAVELETS * WEI-CHANG SHANN AND JANN-CHANG YAN Abstract. Scaling equations are used to derive formulae of quadratures involving polynomials and scaling/wavelet
More information2 (Bonus). Let A X consist of points (x, y) such that either x or y is a rational number. Is A measurable? What is its Lebesgue measure?
MA 645-4A (Real Analysis), Dr. Chernov Homework assignment 1 (Due 9/5). Prove that every countable set A is measurable and µ(a) = 0. 2 (Bonus). Let A consist of points (x, y) such that either x or y is
More informationA PROJECTED HESSIAN GAUSS-NEWTON ALGORITHM FOR SOLVING SYSTEMS OF NONLINEAR EQUATIONS AND INEQUALITIES
IJMMS 25:6 2001) 397 409 PII. S0161171201002290 http://ijmms.hindawi.com Hindawi Publishing Corp. A PROJECTED HESSIAN GAUSS-NEWTON ALGORITHM FOR SOLVING SYSTEMS OF NONLINEAR EQUATIONS AND INEQUALITIES
More informationLucio Demeio Dipartimento di Ingegneria Industriale e delle Scienze Matematiche
Scuola di Dottorato THE WAVE EQUATION Lucio Demeio Dipartimento di Ingegneria Industriale e delle Scienze Matematiche Lucio Demeio - DIISM wave equation 1 / 44 1 The Vibrating String Equation 2 Second
More informationWavelet transforms and compression of seismic data
Mathematical Geophysics Summer School, Stanford, August 1998 Wavelet transforms and compression of seismic data G. Beylkin 1 and A. Vassiliou 2 I Introduction Seismic data compression (as it exists today)
More information3 (Due ). Let A X consist of points (x, y) such that either x or y is a rational number. Is A measurable? What is its Lebesgue measure?
MA 645-4A (Real Analysis), Dr. Chernov Homework assignment 1 (Due ). Show that the open disk x 2 + y 2 < 1 is a countable union of planar elementary sets. Show that the closed disk x 2 + y 2 1 is a countable
More informationLECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES)
LECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES) RAYTCHO LAZAROV 1 Notations and Basic Functional Spaces Scalar function in R d, d 1 will be denoted by u,
More informationOn the Hilbert Transform of Wavelets
On the Hilbert Transform of Wavelets Kunal Narayan Chaudhury and Michael Unser Abstract A wavelet is a localized function having a prescribed number of vanishing moments. In this correspondence, we provide
More informationApplied Analysis (APPM 5440): Final exam 1:30pm 4:00pm, Dec. 14, Closed books.
Applied Analysis APPM 44: Final exam 1:3pm 4:pm, Dec. 14, 29. Closed books. Problem 1: 2p Set I = [, 1]. Prove that there is a continuous function u on I such that 1 ux 1 x sin ut 2 dt = cosx, x I. Define
More informationECE 901 Lecture 16: Wavelet Approximation Theory
ECE 91 Lecture 16: Wavelet Approximation Theory R. Nowak 5/17/29 1 Introduction In Lecture 4 and 15, we investigated the problem of denoising a smooth signal in additive white noise. In Lecture 4, we considered
More informationInitial Boundary Value Problems for Scalar and Vector Burgers Equations
Initial Boundary Value Problems for Scalar and Vector Burgers Equations By K. T. Joseph and P. L. Sachdev In this article we stu Burgers equation and vector Burgers equation with initial and boundary conditions.
More informationNonstationary Subdivision Schemes and Totally Positive Refinable Functions
Nonstationary Subdivision Schemes and Totally Positive Refinable Functions Laura Gori and Francesca Pitolli December, 2007 Abstract In this paper we construct a class of totally positive refinable functions,
More informationWAVELETS WITH COMPOSITE DILATIONS
ELECTRONIC RESEARCH ANNOUNCEMENTS OF THE AMERICAN MATHEMATICAL SOCIETY Volume 00, Pages 000 000 (Xxxx XX, XXXX S 1079-6762(XX0000-0 WAVELETS WITH COMPOSITE DILATIONS KANGHUI GUO, DEMETRIO LABATE, WANG-Q
More informationSome Aspects of Solutions of Partial Differential Equations
Some Aspects of Solutions of Partial Differential Equations K. Sakthivel Department of Mathematics Indian Institute of Space Science & Technology(IIST) Trivandrum - 695 547, Kerala Sakthivel@iist.ac.in
More informationRay equations of a weak shock in a hyperbolic system of conservation laws in multi-dimensions
Proc. Indian Acad. Sci. (Math. Sci.) Vol. 126, No. 2, May 2016, pp. 199 206. c Indian Academy of Sciences Ray equations of a weak in a hyperbolic system of conservation laws in multi-dimensions PHOOLAN
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 informationDeterminantal point processes and random matrix theory in a nutshell
Determinantal point processes and random matrix theory in a nutshell part I Manuela Girotti based on M. Girotti s PhD thesis and A. Kuijlaars notes from Les Houches Winter School 202 Contents Point Processes
More informationConstruction of Multivariate Compactly Supported Orthonormal Wavelets
Construction of Multivariate Compactly Supported Orthonormal Wavelets Ming-Jun Lai Department of Mathematics The University of Georgia Athens, GA 30602 April 30, 2004 Dedicated to Professor Charles A.
More informationORTHONORMAL SAMPLING FUNCTIONS
ORTHONORMAL SAMPLING FUNCTIONS N. KAIBLINGER AND W. R. MADYCH Abstract. We investigate functions φ(x) whose translates {φ(x k)}, where k runs through the integer lattice Z, provide a system of orthonormal
More informationAffine and Quasi-Affine Frames on Positive Half Line
Journal of Mathematical Extension Vol. 10, No. 3, (2016), 47-61 ISSN: 1735-8299 URL: http://www.ijmex.com Affine and Quasi-Affine Frames on Positive Half Line Abdullah Zakir Husain Delhi College-Delhi
More informationTHE objective of this paper is to introduce a comparative
Proceedings of the World Congress on Engineering 23 Vol I, WCE 23, July 3-5, 23, London, U.K. Numerical Solution of Fractional Integro-differential Equation by Using Cubic B-spline Wavelets Khosrow Maleknejad,
More informationOn the fast algorithm for multiplication of functions in the wavelet bases
Published in Proceedings of the International Conference Wavelets and Applications, Toulouse, 1992; Y. Meyer and S. Roques, edt., Editions Frontieres, 1993 On the fast algorithm for multiplication of functions
More informationRadon measure solutions for scalar. conservation laws. A. Terracina. A.Terracina La Sapienza, Università di Roma 06/09/2017
Radon measure A.Terracina La Sapienza, Università di Roma 06/09/2017 Collaboration Michiel Bertsch Flavia Smarrazzo Alberto Tesei Introduction Consider the following Cauchy problem { ut + ϕ(u) x = 0 in
More informationEntropic Schemes for Conservation Laws
CONSTRUCTVE FUNCTON THEORY, Varna 2002 (B. Bojanov, Ed.), DARBA, Sofia, 2002, pp. 1-6. Entropic Schemes for Conservation Laws Bojan Popov A new class of Godunov-type numerical methods (called here entropic)
More informationHigh order geometric smoothness for conservation laws
High order geometric smoothness for conservation laws Martin Campos-Pinto, Albert Cohen, and Pencho Petrushev Abstract The smoothness of the solutions of 1D scalar conservation laws is investigated and
More informationfor all subintervals I J. If the same is true for the dyadic subintervals I D J only, we will write ϕ BMO d (J). In fact, the following is true
3 ohn Nirenberg inequality, Part I A function ϕ L () belongs to the space BMO() if sup ϕ(s) ϕ I I I < for all subintervals I If the same is true for the dyadic subintervals I D only, we will write ϕ BMO
More informationFUNCTION BASES FOR TOPOLOGICAL VECTOR SPACES. Yılmaz Yılmaz
Topological Methods in Nonlinear Analysis Journal of the Juliusz Schauder Center Volume 33, 2009, 335 353 FUNCTION BASES FOR TOPOLOGICAL VECTOR SPACES Yılmaz Yılmaz Abstract. Our main interest in this
More informationConservation Laws and Finite Volume Methods
Conservation Laws and Finite Volume Methods AMath 574 Winter Quarter, 2011 Randall J. LeVeque Applied Mathematics University of Washington January 3, 2011 R.J. LeVeque, University of Washington AMath 574,
More informationCONVERGENCE THEORY. G. ALLAIRE CMAP, Ecole Polytechnique. 1. Maximum principle. 2. Oscillating test function. 3. Two-scale convergence
1 CONVERGENCE THEOR G. ALLAIRE CMAP, Ecole Polytechnique 1. Maximum principle 2. Oscillating test function 3. Two-scale convergence 4. Application to homogenization 5. General theory H-convergence) 6.
More informationRing-like structures of frequency domains of wavelets
Ring-like structures of frequency domains of wavelets Zhihua Zhang and Naoki aito Dept. of Math., Univ. of California, Davis, California, 95616, UA. E-mail: zzh@ucdavis.edu saito@math.ucdavis.edu Abstract.
More informationQuintic Functional Equations in Non-Archimedean Normed Spaces
Journal of Mathematical Extension Vol. 9, No., (205), 5-63 ISSN: 735-8299 URL: http://www.ijmex.com Quintic Functional Equations in Non-Archimedean Normed Spaces A. Bodaghi Garmsar Branch, Islamic Azad
More informationSchemes. Philipp Keding AT15, San Antonio, TX. Quarklet Frames in Adaptive Numerical. Schemes. Philipp Keding. Philipps-University Marburg
AT15, San Antonio, TX 5-23-2016 Joint work with Stephan Dahlke and Thorsten Raasch Let H be a Hilbert space with its dual H. We are interested in the solution u H of the operator equation Lu = f, with
More informationPractice Exercises on Differential Equations
Practice Exercises on Differential Equations What follows are some exerices to help with your studying for the part of the final exam on differential equations. In this regard, keep in mind that the exercises
More informationSolutions of Burgers Equation
ISSN 749-3889 (print, 749-3897 (online International Journal of Nonlinear Science Vol.9( No.3,pp.9-95 Solutions of Burgers Equation Ch. Srinivasa Rao, Engu Satyanarayana Department of Mathematics, Indian
More informationChapter Two: Numerical Methods for Elliptic PDEs. 1 Finite Difference Methods for Elliptic PDEs
Chapter Two: Numerical Methods for Elliptic PDEs Finite Difference Methods for Elliptic PDEs.. Finite difference scheme. We consider a simple example u := subject to Dirichlet boundary conditions ( ) u
More informationStrengthened Sobolev inequalities for a random subspace of functions
Strengthened Sobolev inequalities for a random subspace of functions Rachel Ward University of Texas at Austin April 2013 2 Discrete Sobolev inequalities Proposition (Sobolev inequality for discrete images)
More informationA New Fourth-Order Non-Oscillatory Central Scheme For Hyperbolic Conservation Laws
A New Fourth-Order Non-Oscillatory Central Scheme For Hyperbolic Conservation Laws A. A. I. Peer a,, A. Gopaul a, M. Z. Dauhoo a, M. Bhuruth a, a Department of Mathematics, University of Mauritius, Reduit,
More informationFast Marching Methods for Front Propagation Lecture 1/3
Outline for Front Propagation Lecture 1/3 M. Falcone Dipartimento di Matematica SAPIENZA, Università di Roma Autumn School Introduction to Numerical Methods for Moving Boundaries ENSTA, Paris, November,
More informationAdaptive Control with Multiresolution Bases
Adaptive Control with Multiresolution Bases Christophe P. Bernard Centre de Mathématiques Appliquées Ecole Polytechnique 91128 Palaiseau cedex, FRANCE bernard@cmapx.polytechnique.fr Jean-Jacques E. Slotine
More informationNoncommutative Uncertainty Principle
Noncommutative Uncertainty Principle Zhengwei Liu (joint with Chunlan Jiang and Jinsong Wu) Vanderbilt University The 12th East Coast Operator Algebras Symposium, Oct 12, 2014 Z. Liu (Vanderbilt) Noncommutative
More informationNONCLASSICAL SHOCK WAVES OF CONSERVATION LAWS: FLUX FUNCTION HAVING TWO INFLECTION POINTS
Electronic Journal of Differential Equations, Vol. 2006(2006), No. 149, pp. 1 18. ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu (login: ftp) NONCLASSICAL
More informationRational Chebyshev pseudospectral method for long-short wave equations
Journal of Physics: Conference Series PAPER OPE ACCESS Rational Chebyshev pseudospectral method for long-short wave equations To cite this article: Zeting Liu and Shujuan Lv 07 J. Phys.: Conf. Ser. 84
More informationEnergy method for wave equations
Energy method for wave equations Willie Wong Based on commit 5dfb7e5 of 2017-11-06 13:29 Abstract We give an elementary discussion of the energy method (and particularly the vector field method) in the
More informationREGULARITY AND CONSTRUCTION OF BOUNDARY MULTIWAVELETS
REGULARITY AND CONSTRUCTION OF BOUNDARY MULTIWAVELETS FRITZ KEINERT Abstract. The conventional way of constructing boundary functions for wavelets on a finite interval is to form linear combinations of
More informationSOME SMOOTH COMPACTLY SUPPORTED TIGHT WAVELET FRAMES WITH VANISHING MOMENTS
SOME SMOOTH COMPACTLY SUPPORTED TIGHT WAVELET FRAMES WITH VANISHING MOMENTS A. SAN ANTOLÍN AND R. A. ZALIK Abstract. Let A R d d, d 1 be a dilation matrix with integer entries and det A = 2. We construct
More informationA CLASS OF EVEN DEGREE SPLINES OBTAINED THROUGH A MINIMUM CONDITION
STUDIA UNIV. BABEŞ BOLYAI, MATHEMATICA, Volume XLVIII, Number 3, September 2003 A CLASS OF EVEN DEGREE SPLINES OBTAINED THROUGH A MINIMUM CONDITION GH. MICULA, E. SANTI, AND M. G. CIMORONI Dedicated to
More informationError Analysis for H 1 Based Wavelet Interpolations
Error Analysis for H 1 Based Wavelet Interpolations Tony F. Chan Hao-Min Zhou Tie Zhou Abstract We rigorously study the error bound for the H 1 wavelet interpolation problem, which aims to recover missing
More informationJim Lambers ENERGY 281 Spring Quarter Lecture 5 Notes
Jim ambers ENERGY 28 Spring Quarter 27-8 ecture 5 Notes These notes are based on Rosalind Archer s PE28 lecture notes, with some revisions by Jim ambers. Fourier Series Recall that in ecture 2, when we
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 informationMath 660-Lecture 23: Gudonov s method and some theories for FVM schemes
Math 660-Lecture 3: Gudonov s method and some theories for FVM schemes 1 The idea of FVM (You can refer to Chapter 4 in the book Finite volume methods for hyperbolic problems ) Consider the box [x 1/,
More informationWeighted Essentially Non-Oscillatory limiters for Runge-Kutta Discontinuous Galerkin Methods
Weighted Essentially Non-Oscillatory limiters for Runge-Kutta Discontinuous Galerkin Methods Jianxian Qiu School of Mathematical Science Xiamen University jxqiu@xmu.edu.cn http://ccam.xmu.edu.cn/teacher/jxqiu
More information