ON A NUMERICAL SOLUTION OF THE LAPLACE EQUATION
|
|
- Marlene Ray
- 6 years ago
- Views:
Transcription
1 ON A NUMERICAL SOLUTION OF THE LAPLACE EQUATION Jasmina Veta Buralieva (joint work with E. Hadzieva and K. Hadzi-Velkova Saneva ) GFTA 2015 Ohrid, August 2015 (Faculty of Informatics,UGD-Stip) On a numerical solution of the Laplace equation GFTA, / 34
2 Outline - Classical Galerkin method for ODE - Wavelets and MRA - Wavelet-Galerkin method for ODE - Transformation of the Laplace DE - Application of W-G method on the Laplace DE (Faculty of Informatics,UGD-Stip) On a numerical solution of the Laplace equation GFTA, / 34
3 Classical Galerkin metod for ODE Classical Galerkin metod for ODE Sturm-Liouville equation with BC Lu(t) d dt ( g(t) du dt ) + h(t)u(t) = f (t), a t b (1) u(a) = c, u(b) = d. (2) 1) {v j } - complete orthonormal system for L 2 ([a, b]) 2) every v j C 2 ([a, b]) 3) v j (a) = c, v j (b) = d. Approximation u s of the exact solution u u s = k Λ x k v k. (3) (Faculty of Informatics,UGD-Stip) On a numerical solution of the Laplace equation GFTA, / 34
4 Classical Galerkin metod for ODE Criterion for coefficients x k < Lu s, v j >=< f, v j >, j Λ. (4) If we substitute the equation (3) in (4) we obtain Lv k, v j x k = f, v j, j Λ. (5) k Λ A = [a j,k ] j,k Λ, a j,k = Lv k, v j ; X = (x k ) k Λ ; Y = (y k ) k Λ, y k = f, v k AX = Y. (6) Wavelet-Galerkin method: functions v j are wavelets (Faculty of Informatics,UGD-Stip) On a numerical solution of the Laplace equation GFTA, / 34
5 Wavelets Wavelets wavelet ψ: L 2 function which satisfy the admissibility condition C ψ = The condition (7) implies that ˆψ(0) = ˆψ(ω) 2 dω <. (7) ω ψ(t)dt = 0. wavelets ψ a,b (t) = 1 a ψ ( t b a ), a > 0, b R. (Faculty of Informatics,UGD-Stip) On a numerical solution of the Laplace equation GFTA, / 34
6 Multiresolution analysis Multiresolution analysis (MRA) Multiresolution analysis of the space L 2 (R) consists of a sequence of closed subspace {V j } j= with the following properties: 1. V j V j+1 2. j Z V j = L 2 (R) 3. j Z V j = {0} 4. f (t) V j f (2t) V j+1 5. f (t) V j f (t k) V j, k Z 6. there exists a function φ (called scaling function or father wavelet) such that φ j,k (t) = 2 j/2 φ(2 j t k), k Z constitute orthonormal basis for corresponding subspace V j. (Faculty of Informatics,UGD-Stip) On a numerical solution of the Laplace equation GFTA, / 34
7 Multiresolution analysis Let φ L 2 (R) be compactly supported scaling function of MRA. Then 1) φ(t)dt 0 2)φ(t) = k Z a kφ(2t k), where a k are real coefficients and a k 0 for only finitely many k Z (the number of nonzero coefficients a k is denoted by L). 3) φ j,k (t) = 2 j/2 φ(2 j t k), j, k Z are orthonormal in L 2 (R) i.e. φ(t n)φ(t k)dt = δ k,n, (8) where δ n,k = { 0, n k 1, n = k. (9) (Faculty of Informatics,UGD-Stip) On a numerical solution of the Laplace equation GFTA, / 34
8 Multiresolution analysis One can construct wavelet ψ such that ψ j,k (t) = 2 j/2 φ(2 j t k), j, k Z constitute an orthonormal basis for L 2 (R). Daubechies scaling function Daubechies wavelet function ψ(t) = L 1 φ(t) = a k φ(2t k) (10) 1 k=2 L k=0 ( 1) k a 1 k φ(2t k) (11) where L is a positive even integer and denotes the genus of the Daubechies wavelet. (Faculty of Informatics,UGD-Stip) On a numerical solution of the Laplace equation GFTA, / 34
9 Wavelet-Galerkin method for ordinary differential equations Wavelet-Galerkin method for ODE g(t)u (t) + g (t)u (t) + h(t)u(t) = 0, t [a, b], (12) with BC u(a) = c, u(b) = d. (13) Approximate solution u j (t) = 2 j k=1 L where φ is the scaling function of MRA. c k φ j,k (t), k Z, (14) (Faculty of Informatics,UGD-Stip) On a numerical solution of the Laplace equation GFTA, / 34
10 Wavelet-Galerkin method for ordinary differential equations Remark There are no closed-form formulas for the Daubechies wavelets and scaling functions. W-G method with Daubechies scaling functions: homogeneous differential equations. (Faculty of Informatics,UGD-Stip) On a numerical solution of the Laplace equation GFTA, / 34
11 Wavelet-Galerkin method for ordinary differential equations For j = 0 and L = 4 g(t) d 2 dt 2 u 0 (t) = 1 k= 3 1 k= 3 c k φ(t k), t [a, b]. (15) c k φ(t k) + g (t) d dt +h(t) 1 k= 3 1 k= 3 c k φ(t k)+ c k φ(t k) = 0. (16) (Faculty of Informatics,UGD-Stip) On a numerical solution of the Laplace equation GFTA, / 34
12 Wavelet-Galerkin method for ordinary differential equations Taking inner product with φ(t n), n { 3, 2, 1, 0, 1}, we obtain c k Ω n k + c k a n,k + c k s n,k = 0, (17) k= 3 k= 3 k= 3 where Ω n k = a n,k = g(t)φ (t k)φ(t n)dt, (18) g (t)φ (t k)φ(t n)dt, (19) s n,k = 4 3 h(t)φ(t k)φ(t n)dt. (20) (Faculty of Informatics,UGD-Stip) On a numerical solution of the Laplace equation GFTA, / 34
13 Wavelet-Galerkin method for ordinary differential equations By using BC (13) we obtain and u 0 (a) = u 0 (b) = 1 k= 3 1 k= 3 c k φ(a k) = c (21) c k φ(b k) = d (22) We replace the first and the last equation of system (17) by (21) and (22), respectively and obtain the matrix equation TC = B (23) (Faculty of Informatics,UGD-Stip) On a numerical solution of the Laplace equation GFTA, / 34
14 Wavelet-Galerkin method for ordinary differential equations T = φ(a + 3) φ(a + 2) Ω a 2, 3 + s 2, 3 Ω a 2, 2 + s 2, 2 Ω a 1, 3 + s 1, 3 Ω a 1, 2 + s 1, 2 Ω a 0, 3 + s 0, 3 Ω a 0, 2 + s 0, 2 φ(b + 3) φ(b + 2) φ(a + 1) φ(a) φ(a 1) Ω a 2, 1 + s 2, 1 Ω a 2,0 + s 2,0 Ω a 2,1 + s 2,1 Ω a 1, 1 + s 1, 1 Ω a 1,0 + s 1,0 Ω a 1,1 + s 1,1 Ω a 0, 1 + s 0, 1 Ω 0 + a 0,0 + s 0,0 Ω a 0,1 + a 0,1 φ(b + 1) φ(b) φ(b 1) (Faculty of Informatics,UGD-Stip) On a numerical solution of the Laplace equation GFTA, / 34
15 Wavelet-Galerkin method for ordinary differential equations C = c 3 c 2 c 1 c 0 c 1, B = a b. (Faculty of Informatics,UGD-Stip) On a numerical solution of the Laplace equation GFTA, / 34
16 Transformation of the Laplace DE Transformation of the Laplace DE The famous Laplace equation with substitutions u 2 u x u y u z 2 = 0 (24) has the form x = r sin θ cos ϕ, y = r sin θ sin ϕ, z = r cos θ, u r ( r 2 u ) + 1 r sin θ ( sin θ u ) + 1 θ θ sin 2 θ 2 u = 0. (25) ϕ2 (Faculty of Informatics,UGD-Stip) On a numerical solution of the Laplace equation GFTA, / 34
17 Transformation of the Laplace DE Fourier method subsumes u = u(r, θ, ϕ) can be represented u(r, ϕ, θ) = R(r)Φ(ϕ)Θ(θ) ΦΘ d dr (r 2 R ) + RΦ 1 d sin θ dθ (sin θθ ) + 1 sin 2 RΘΦ = 0, θ 1 R d [ 1 dr (r 2 R ) = Θ 1 sin θ d dθ (sin θ Θ ) + 1 ] sin 2 θ Φ, (26) Φ 1 R d dr (r 2 R ) = λ, 1 Θ 1 sin θ d dθ (sin θ Θ ) + 1 sin 2 θ Φ Φ = λ. (27) (Faculty of Informatics,UGD-Stip) On a numerical solution of the Laplace equation GFTA, / 34
18 Transformation of the Laplace DE System of three ODEs 1 Θ sin θ d dθ dθ (sin θ 1 R d dr (r 2 R ) = λ Φ Φ = µ dθ ) + λ sin2 θ = µ (28) The first ODE is Cauchy-Euler equation r 2 R + 2rR λr = 0. (29) Its exact solution is for R(r) = C 1 r n + C 2 1 r n+1 λ = n(n + 1). (Faculty of Informatics,UGD-Stip) On a numerical solution of the Laplace equation GFTA, / 34
19 Transformation of the Laplace DE The second ODE is Sturm-Liouville equation general solution sin θ d dθ Its exact solution is Φ + µφ = 0, (30) Φ(ϕ) = A cos mϕ + B sin mϕ, for µ = m 2, where m = 1, 2,.... The third ODE is Sturm-Liouville equation ( sin θ dθ dθ ) ( ) + Θ λ sin 2 θ µ = 0. (31) P n,m = (1 x 2 ) m d m P n(x) 2 dx m = (1 x 2 ) d n+m n!2 n dx n+m [(x 2 1) n ], (32) for λ = n(n + 1) and µ = m 2 where P n (x) = 1 n!2 n d n dx n [(x 2 1) n ], x = cos θ. (Faculty of Informatics,UGD-Stip) On a numerical solution of the Laplace equation GFTA, / 34 m 2
20 Application of the W-G method on the Laplace equation Application of the W-G method on the Laplace equation satisfies 1 2 t2, t [0, 1] t φ(t) = 2 3t 3 2, t [1, 2] 1 2 t2 3t + 9 2, t [2, 3] 0, t / [0, 3] (33) so L = 4. φ(t) = 1 4 φ(2t) φ(2t 1) φ(2t 2) + 1 φ(2t 3), 4 (Faculty of Informatics,UGD-Stip) On a numerical solution of the Laplace equation GFTA, / 34
21 Application of the W-G method on the Laplace equation 1. λ = 0 and µ = 4 The first ODE is with the BC: R(1) = 1, R(3) = 0. Approximate solution is r 2 R + 2rR = 0 (34) R 0 (r) = { c 1 φ(r + 1) + c 0 φ(r) + c 1 φ(r 1), r [1, 2] c 0 φ(r) + c 1 φ(r 1), r [2, 3] where c 1 = , c 0 = , c 1 = 0. (Faculty of Informatics,UGD-Stip) On a numerical solution of the Laplace equation GFTA, / 34
22 Application of the W-G method on the Laplace equation Table: Comparison of results Case t numerical solution R 0 exact solution R absolute error (Faculty of Informatics,UGD-Stip) On a numerical solution of the Laplace equation GFTA, / 34
23 Application of the W-G method on the Laplace equation (Faculty of Informatics,UGD-Stip) On a numerical solution of the Laplace equation GFTA, / 34
24 Application of the W-G method on the Laplace equation The second ODE is with BC: Φ(0) = 1, Φ( π 4 ) = 1. Approximate solution is Φ + 4Φ = 0 (35) Φ 0 (ϕ) = c 2 φ(ϕ + 2) + c 1 φ(ϕ + 1) + c 0 φ(ϕ), ϕ [0, π 4 ] where c 2 = , c 1 = , c 0 = (Faculty of Informatics,UGD-Stip) On a numerical solution of the Laplace equation GFTA, / 34
25 Application of the W-G method on the Laplace equation Table: Comparison of results Case t numerical solution Φ 0 exact solution Φ absolute error π/ (Faculty of Informatics,UGD-Stip) On a numerical solution of the Laplace equation GFTA, / 34
26 Application of the W-G method on the Laplace equation (Faculty of Informatics,UGD-Stip) On a numerical solution of the Laplace equation GFTA, / 34
27 Application of the W-G method on the Laplace equation The third ODE is with BC: Θ(1) = 1, Θ(2) = 2. sin 2 (θ)θ + cos(θ) sin(θ)θ 4Θ = 0 (36) (Faculty of Informatics,UGD-Stip) On a numerical solution of the Laplace equation GFTA, / 34
28 Application of the W-G method on the Laplace equation 2. For λ = 0 and µ = 2 The second ODE is with BC: Φ(0) = 1, Φ( π 4 ) = 1. Approximate solution is Φ + 2Φ = 0 (37) where Φ 0 (ϕ) = c 2 φ(ϕ + 2) + c 1 φ(ϕ + 1) + c 0 φ(ϕ), ϕ [0, π 4 ] c 2 = , c 1 = , c 0 = (Faculty of Informatics,UGD-Stip) On a numerical solution of the Laplace equation GFTA, / 34
29 Application of the W-G method on the Laplace equation The third ODE is with BC: Θ(1) = 1, Θ(2) = 2. Approximate solution is sin 2 (θ)θ + cos(θ) sin(θ)θ 2Θ = 0 (38) Θ 0 (θ) = c 1 φ(θ + 1) + c 0 φ(θ) + c 1 φ(θ 1), θ [1, 2] where c 1 = 2.427, c 0 = , c 1 = (Faculty of Informatics,UGD-Stip) On a numerical solution of the Laplace equation GFTA, / 34
30 Application of the W-G method on the Laplace equation Table: Numerical results Case t numerical solution Φ 0 Case t numerical solution Θ π/ // // // // 2 2 (Faculty of Informatics,UGD-Stip) On a numerical solution of the Laplace equation GFTA, / 34
31 Application of the W-G method on the Laplace equation 3. λ = 1 and µ = 4 The first ODE is with BC: R(1) = 1, R(3) = 0. Approximate solution is r 2 R + 2rR R = 0 (39) where R 0 (r) = { c 1 φ(r + 1) + c 0 φ(r) + c 1 φ(r 1), r [1, 2] c 0 φ(r) + c 1 φ(r 1), r [2, 3] c 1 = , c 0 = , c 1 = 0. (Faculty of Informatics,UGD-Stip) On a numerical solution of the Laplace equation GFTA, / 34
32 Application of the W-G method on the Laplace equation The third ODE is sin 2 (θ)θ + cos(θ) sin(θ)θ + Θ with BC: Θ(1) = 1, Θ(2) = 2. Approximate solution is ( ) sin 2 (θ) 4 = 0 (40) Θ 0 (θ) = c 1 φ(θ + 1) + c 0 φ(θ) + c 1 φ(θ 1), θ [1, 2] where c 1 = , c 0 = , c 1 = (Faculty of Informatics,UGD-Stip) On a numerical solution of the Laplace equation GFTA, / 34
33 Application of the W-G method on the Laplace equation Table: Numerical results Case t numerical solution Φ 0 Case t numerical solution Θ (Faculty of Informatics,UGD-Stip) On a numerical solution of the Laplace equation GFTA, / 34
34 References A. H. Siddiqi, Applied Functional Analysis: Numerical Methods, Wavelet Methods and Image Processing, Markel Dekker, New York, Anandita D., A wavelet-galerkin method for the solution of partial differential equation, master thesis, 2011 C. Qian, J. Weiss, Wavelets and the numerical solution of partial differential equations, Journal of Computational Physics, Vol. 106, Issue 1, 1993, pp D. S. Mitrinovič, J. D. Kečki: Jednačine matematičke fizike, Beograd, I. Daubeshies, Ten lectures on Wavelets, Philadelphia: SIAM, I. Daubechies, Orthonormal bases of compactly suppoted wavelets, Commun. Pure Appl. Math., 41, 1988, pp G. G. Walter, X. Shen, Wavelets and Other Orthogonal Systems With Application, CRS Press, Secon Editiotn,2000 M. W. Frazier, An Introduction to Wavelets Through Linear Algebra, Springer-Verlag, New York, S. Kostadinova, J. Veta Buralieva, K. Hadzi-Velkova Saneva (2013), Wavelet-Galerkin solution of some ordinary differential equation,proceedings of XI Internation Conference ETAI 2013, September 2013, Ohrid, Macedonia. S. Mallat, Multiresolution approximation and wavelets, Trans. Amer. Math. Soc., 315, 1989, pp T. Lofti, K. Mahdiani, Numerical solution of boundary value problem by using wavelet-galerkin methods, Mathematical Sciences, Vol. 1, No. 3, 2007, pp Vladimirov V. S., Uravneninija matematicheskori fiziki, Nauka, Moskva,1967. V. Mishra, Sabina, Wavelet Galerkin solutions of ordinary differential equations, Int. Joirnal of Math.Analysis, Vol. 5(9), 2001, (Faculty of Informatics,UGD-Stip) On a numerical solution of the Laplace equation GFTA, / 34
Ninth 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 informationBiorthogonal 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 informationConstruction of Biorthogonal B-spline Type Wavelet Sequences with Certain Regularities
Illinois Wesleyan University From the SelectedWorks of Tian-Xiao He 007 Construction of Biorthogonal B-spline Type Wavelet Sequences with Certain Regularities Tian-Xiao He, Illinois Wesleyan University
More informationLecture 7 Multiresolution Analysis
David Walnut Department of Mathematical Sciences George Mason University Fairfax, VA USA Chapman Lectures, Chapman University, Orange, CA Outline Definition of MRA in one dimension Finding the wavelet
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 informationWAVELET EXPANSIONS OF DISTRIBUTIONS
WAVELET EXPANSIONS OF DISTRIBUTIONS JASSON VINDAS Abstract. These are lecture notes of a talk at the School of Mathematics of the National University of Costa Rica. The aim is to present a wavelet expansion
More information1. Fourier Transform (Continuous time) A finite energy signal is a signal f(t) for which. f(t) 2 dt < Scalar product: f(t)g(t)dt
1. Fourier Transform (Continuous time) 1.1. Signals with finite energy A finite energy signal is a signal f(t) for which Scalar product: f(t) 2 dt < f(t), g(t) = 1 2π f(t)g(t)dt The Hilbert space of all
More informationCourse and Wavelets and Filter Banks
Course 8.327 and.30 Wavelets and Filter Banks Multiresolution Analysis (MRA): Requirements for MRA; Nested Spaces and Complementary Spaces; Scaling Functions and Wavelets Scaling Functions and Wavelets
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 informationInner products. Theorem (basic properties): Given vectors u, v, w in an inner product space V, and a scalar k, the following properties hold:
Inner products Definition: An inner product on a real vector space V is an operation (function) that assigns to each pair of vectors ( u, v) in V a scalar u, v satisfying the following axioms: 1. u, v
More informationFrom Fourier to Wavelets in 60 Slides
From Fourier to Wavelets in 60 Slides Bernhard G. Bodmann Math Department, UH September 20, 2008 B. G. Bodmann (UH Math) From Fourier to Wavelets in 60 Slides September 20, 2008 1 / 62 Outline 1 From Fourier
More informationLEGENDRE POLYNOMIALS AND APPLICATIONS. We construct Legendre polynomials and apply them to solve Dirichlet problems in spherical coordinates.
LEGENDRE POLYNOMIALS AND APPLICATIONS We construct Legendre polynomials and apply them to solve Dirichlet problems in spherical coordinates.. Legendre equation: series solutions The Legendre equation is
More informationBoundary Value Problems in Cylindrical Coordinates
Boundary Value Problems in Cylindrical Coordinates 29 Outline Differential Operators in Various Coordinate Systems Laplace Equation in Cylindrical Coordinates Systems Bessel Functions Wave Equation the
More informationConnection to Laplacian in spherical coordinates (Chapter 13)
Connection to Laplacian in spherical coordinates (Chapter 13) We might often encounter the Laplace equation and spherical coordinates might be the most convenient 2 u(r, θ, φ) = 0 We already saw in Chapter
More informationMath 337, Summer 2010 Assignment 5
Math 337, Summer Assignment 5 Dr. T Hillen, University of Alberta Exercise.. Consider Laplace s equation r r r u + u r r θ = in a semi-circular disk of radius a centered at the origin with boundary conditions
More informationSpherical Coordinates and Legendre Functions
Spherical Coordinates and Legendre Functions Spherical coordinates Let s adopt the notation for spherical coordinates that is standard in physics: φ = longitude or azimuth, θ = colatitude ( π 2 latitude)
More informationSufficient conditions for functions to form Riesz bases in L 2 and applications to nonlinear boundary-value problems
Electronic Journal of Differential Equations, Vol. 200(200), No. 74, pp. 0. ISSN: 072-669. URL: http://ejde.math.swt.edu or http://ejde.math.unt.edu ftp ejde.math.swt.edu (login: ftp) Sufficient conditions
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 informationExpansion of 1/r potential in Legendre polynomials
Expansion of 1/r potential in Legendre polynomials In electrostatics and gravitation, we see scalar potentials of the form V = K d Take d = R r = R 2 2Rr cos θ + r 2 = R 1 2 r R cos θ + r R )2 Use h =
More informationMATH 590: Meshfree Methods
MATH 590: Meshfree Methods Chapter 9: Conditionally Positive Definite Radial Functions Greg Fasshauer Department of Applied Mathematics Illinois Institute of Technology Fall 2010 fasshauer@iit.edu MATH
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 informationMath 300 Winter 2011 Advanced Boundary Value Problems I. Bessel s Equation and Bessel Functions
Math 3 Winter 2 Avance Bounary Value Problems I Bessel s Equation an Bessel Functions Department of Mathematical an Statistical Sciences University of Alberta Bessel s Equation an Bessel Functions We use
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 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 informationBAND-LIMITED REFINABLE FUNCTIONS FOR WAVELETS AND FRAMELETS
BAND-LIMITED REFINABLE FUNCTIONS FOR WAVELETS AND FRAMELETS WEIQIANG CHEN AND SAY SONG GOH DEPARTMENT OF MATHEMATICS NATIONAL UNIVERSITY OF SINGAPORE 10 KENT RIDGE CRESCENT, SINGAPORE 119260 REPUBLIC OF
More informationPHYS 404 Lecture 1: Legendre Functions
PHYS 404 Lecture 1: Legendre Functions Dr. Vasileios Lempesis PHYS 404 - LECTURE 1 DR. V. LEMPESIS 1 Legendre Functions physical justification Legendre functions or Legendre polynomials are the solutions
More informationEvolution equations with spectral methods: the case of the wave equation
Evolution equations with spectral methods: the case of the wave equation Jerome.Novak@obspm.fr Laboratoire de l Univers et de ses Théories (LUTH) CNRS / Observatoire de Paris, France in collaboration with
More informationMath 4263 Homework Set 1
Homework Set 1 1. Solve the following PDE/BVP 2. Solve the following PDE/BVP 2u t + 3u x = 0 u (x, 0) = sin (x) u x + e x u y = 0 u (0, y) = y 2 3. (a) Find the curves γ : t (x (t), y (t)) such that that
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 informationFourier-like Transforms
L 2 (R) Solutions of Dilation Equations and Fourier-like Transforms David Malone December 6, 2000 Abstract We state a novel construction of the Fourier transform on L 2 (R) based on translation and dilation
More informationQuantum Mechanics in 3-Dimensions
Quantum Mechanics in 3-Dimensions Pavithran S Iyer, 2nd yr BSc Physics, Chennai Mathematical Institute Email: pavithra@cmi.ac.in August 28 th, 2009 1 Schrodinger equation in Spherical Coordinates 1.1 Transforming
More informationACM 126a Solutions for Homework Set 4
ACM 26a Solutions for Homewor Set 4 Laurent Demanet March 2, 25 Problem. Problem 7.7 page 36 We need to recall a few standard facts about Fourier series. Convolution: Subsampling (see p. 26): Zero insertion
More informationCONSTRUCTION OF ORTHONORMAL WAVELETS USING KAMPÉ DE FÉRIET FUNCTIONS
PROCEEDINGS OF THE AMERICAN MATHEMATICAL SOCIETY Volume 1 Number 1 Pages 89 94 S -999669-X Article electronically published on May 1 CONSTRUCTION OF ORTHONORMAL WAVELETS USING KAMPÉ DE FÉRIET FUNCTIONS
More informationMRA Frame Wavelets with Certain Regularities Associated with the Refinable Generators of Shift Invariant Spaces
Chapter 6 MRA Frame Wavelets with Certain Regularities Associated with the Refinable Generators of Shift Invariant Spaces Tian-Xiao He Department of Mathematics and Computer Science Illinois Wesleyan University,
More informationJournal of Computational and Applied Mathematics. Relations among eigenvalues of left-definite Sturm Liouville problems
Journal of Computational and Applied Mathematics 236 (2012) 3426 3433 Contents lists available at SciVerse ScienceDirect Journal of Computational and Applied Mathematics journal homepage: www.elsevier.com/locate/cam
More informationBessel s Equation. MATH 365 Ordinary Differential Equations. J. Robert Buchanan. Fall Department of Mathematics
Bessel s Equation MATH 365 Ordinary Differential Equations J. Robert Buchanan Department of Mathematics Fall 2018 Background Bessel s equation of order ν has the form where ν is a constant. x 2 y + xy
More informationWavelet-Based Numerical Homogenization for Scaled Solutions of Linear Matrix Equations
International Journal of Discrete Mathematics 2017; 2(1: 10-16 http://www.sciencepublishinggroup.com/j/dmath doi: 10.11648/j.dmath.20170201.13 Wavelet-Based Numerical Homogenization for Scaled Solutions
More informationA Tilt at TILFs. Rod Nillsen University of Wollongong. This talk is dedicated to Gary H. Meisters
A Tilt at TILFs Rod Nillsen University of Wollongong This talk is dedicated to Gary H. Meisters Abstract In this talk I will endeavour to give an overview of some aspects of the theory of Translation Invariant
More informationTwo special equations: Bessel s and Legendre s equations. p Fourier-Bessel and Fourier-Legendre series. p
LECTURE 1 Table of Contents Two special equations: Bessel s and Legendre s equations. p. 259-268. Fourier-Bessel and Fourier-Legendre series. p. 453-460. Boundary value problems in other coordinate system.
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 informationAMATH 353 Lecture 9. Weston Barger. How to classify PDEs as linear/nonlinear, order, homogeneous or non-homogeneous.
AMATH 353 ecture 9 Weston Barger 1 Exam What you need to know: How to classify PDEs as linear/nonlinear, order, homogeneous or non-homogeneous. The definitions for traveling wave, standing wave, wave train
More informationIntroduction to Sturm-Liouville Theory and the Theory of Generalized Fourier Series
CHAPTER 5 Introduction to Sturm-Liouville Theory and the Theory of Generalized Fourier Series We start with some introductory examples. 5.. Cauchy s equation The homogeneous Euler-Cauchy equation (Leonhard
More informationUNIVERSITY of LIMERICK OLLSCOIL LUIMNIGH
UNIVERSITY of LIMERICK OLLSCOIL LUIMNIGH College of Informatics and Electronics END OF SEMESTER ASSESSMENT PAPER MODULE CODE: MS425 SEMESTER: Autumn 25/6 MODULE TITLE: Applied Analysis DURATION OF EXAMINATION:
More informationFinal Exam May 4, 2016
1 Math 425 / AMCS 525 Dr. DeTurck Final Exam May 4, 2016 You may use your book and notes on this exam. Show your work in the exam book. Work only the problems that correspond to the section that you prepared.
More informationIn this worksheet we will use the eigenfunction expansion to solve nonhomogeneous equation.
Eigen Function Expansion and Applications. In this worksheet we will use the eigenfunction expansion to solve nonhomogeneous equation. a/ The theory. b/ Example: Solving the Euler equation in two ways.
More informationFINAL EXAM, MATH 353 SUMMER I 2015
FINAL EXAM, MATH 353 SUMMER I 25 9:am-2:pm, Thursday, June 25 I have neither given nor received any unauthorized help on this exam and I have conducted myself within the guidelines of the Duke Community
More informationAssignment #09 - Solution Manual
Assignment #09 - Solution Manual 1. Choose the correct statements about representation of a continuous signal using Haar wavelets. 1.5 points The signal is approximated using sin and cos functions. The
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 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 informationHamburger Beiträge zur Angewandten Mathematik
Hamburger Beiträge zur Angewandten Mathematik Error Estimates for Filtered Back Projection Matthias Beckmann and Armin Iske Nr. 2015-03 January 2015 Error Estimates for Filtered Back Projection Matthias
More informationx 3y 2z = 6 1.2) 2x 4y 3z = 8 3x + 6y + 8z = 5 x + 3y 2z + 5t = 4 1.5) 2x + 8y z + 9t = 9 3x + 5y 12z + 17t = 7
Linear Algebra and its Applications-Lab 1 1) Use Gaussian elimination to solve the following systems x 1 + x 2 2x 3 + 4x 4 = 5 1.1) 2x 1 + 2x 2 3x 3 + x 4 = 3 3x 1 + 3x 2 4x 3 2x 4 = 1 x + y + 2z = 4 1.4)
More informationMathematical Modeling using Partial Differential Equations (PDE s)
Mathematical Modeling using Partial Differential Equations (PDE s) 145. Physical Models: heat conduction, vibration. 146. Mathematical Models: why build them. The solution to the mathematical model will
More informationQUALIFYING EXAMINATION Harvard University Department of Mathematics Tuesday August 31, 2010 (Day 1)
QUALIFYING EXAMINATION Harvard University Department of Mathematics Tuesday August 31, 21 (Day 1) 1. (CA) Evaluate sin 2 x x 2 dx Solution. Let C be the curve on the complex plane from to +, which is along
More informationSurvey on Wavelet Transform and Application in ODE and Wavelet Networks
Advances in Dynamical Systems and Applications. ISSN 0973-5321 Volume 1 Number 2 (2006), pp. 129 162 c Research India Publications http://www.ripublication.com/adsa.htm Survey on Wavelet Transform and
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 informationA Present Position-Dependent Conditional Fourier-Feynman Transform and Convolution Product over Continuous Paths
International Journal of Mathematical Analysis Vol. 9, 05, no. 48, 387-406 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/0.988/ijma.05.589 A Present Position-Dependent Conditional Fourier-Feynman Transform
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 informationQuantum Mechanics in Three Dimensions
Physics 342 Lecture 21 Quantum Mechanics in Three Dimensions Lecture 21 Physics 342 Quantum Mechanics I Monday, March 22nd, 21 We are used to the temporal separation that gives, for example, the timeindependent
More informationHilbert Space Problems
Hilbert Space Problems Prescribed books for problems. ) Hilbert Spaces, Wavelets, Generalized Functions and Modern Quantum Mechanics by Willi-Hans Steeb Kluwer Academic Publishers, 998 ISBN -7923-523-9
More informationPhysics 116C Solutions to the Practice Final Exam Fall The method of Frobenius fails because it works only for Fuchsian equations, of the type
Physics 6C Solutions to the Practice Final Exam Fall 2 Consider the differential equation x y +x(x+)y y = () (a) Explain why the method of Frobenius method fails for this problem. The method of Frobenius
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 informationMath 241 Final Exam Spring 2013
Name: Math 241 Final Exam Spring 213 1 Instructor (circle one): Epstein Hynd Wong Please turn off and put away all electronic devices. You may use both sides of a 3 5 card for handwritten notes while you
More informationNotes on Special Functions
Spring 25 1 Notes on Special Functions Francis J. Narcowich Department of Mathematics Texas A&M University College Station, TX 77843-3368 Introduction These notes are for our classes on special functions.
More informationDISCRETE CDF 9/7 WAVELET TRANSFORM FOR FINITE-LENGTH SIGNALS
DISCRETE CDF 9/7 WAVELET TRANSFORM FOR FINITE-LENGTH SIGNALS D. Černá, V. Finěk Department of Mathematics and Didactics of Mathematics, Technical University in Liberec Abstract Wavelets and a discrete
More informationWavelet analysis as a p adic spectral analysis
arxiv:math-ph/0012019v3 23 Feb 2001 Wavelet analysis as a p adic spectral analysis S.V.Kozyrev February 3, 2008 Institute of Chemical Physics, Russian Academy of Science Abstract New orthonormal basis
More informationFrame Wavelet Sets in R d
Frame Wavelet Sets in R d X. DAI, Y. DIAO Department of Mathematics University of North Carolina at Charlotte Charlotte, NC 28223 xdai@uncc.edu Q. GU Department of Mathematics Each China Normal University
More informationThe Legendre Wavelet Method for Solving Singular Integro-differential Equations
Computational Methods for Differential Equations http://cmde.tabrizu.ac.ir Vol. 2, No. 2, 2014, pp. 62-68 The Legendre Wavelet Method for Solving Singular Integro-differential Equations Naser Aghazadeh,
More informationThe purpose of this lecture is to present a few applications of conformal mappings in problems which arise in physics and engineering.
Lecture 16 Applications of Conformal Mapping MATH-GA 451.001 Complex Variables The purpose of this lecture is to present a few applications of conformal mappings in problems which arise in physics and
More informationImage Compression by Using Haar Wavelet Transform and Singular Value Decomposition
Master Thesis Image Compression by Using Haar Wavelet Transform and Singular Value Decomposition Zunera Idrees 9--5 Eliza Hashemiaghjekandi 979-- Subject: Mathematics Level: Advance Course code: 4MAE Abstract
More informationMATH 31BH Homework 5 Solutions
MATH 3BH Homework 5 Solutions February 4, 204 Problem.8.2 (a) Let x t f y = x 2 + y 2 + 2z 2 and g(t) = t 2. z t 3 Then by the chain rule a a a D(g f) b = Dg f b Df b c c c = [Dg(a 2 + b 2 + 2c 2 )] [
More informationwhere u 0 is given. In all our applications it is possible to subtract a known function v(x, t) from u(x, t) so that w(x, t) = u(x, t) v(x, t)
CHAPTER PDE Partial Differential Equations in Two Independent Variables D. An Overview Drawing on the Sturm-Liouville eigenvalue theory and the approximation of functions we are now ready to develop the
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 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 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 informationDevelopment and Applications of Wavelets in Signal Processing
Development and Applications of Wavelets in Signal Processing Mathematics 097: Senior Conference Paper Published May 014 David Nahmias dnahmias1@gmailcom Abstract Wavelets have many powerful applications
More informationVibrating-string problem
EE-2020, Spring 2009 p. 1/30 Vibrating-string problem Newton s equation of motion, m u tt = applied forces to the segment (x, x, + x), Net force due to the tension of the string, T Sinθ 2 T Sinθ 1 T[u
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 informationCh 4: The Continuous-Time Fourier Transform
Ch 4: The Continuous-Time Fourier Transform Fourier Transform of x(t) Inverse Fourier Transform jt X ( j) x ( t ) e dt jt x ( t ) X ( j) e d 2 Ghulam Muhammad, King Saud University Continuous-time aperiodic
More informationLecture 15: Time and Frequency Joint Perspective
WAVELETS AND MULTIRATE DIGITAL SIGNAL PROCESSING Lecture 15: Time and Frequency Joint Perspective Prof.V.M.Gadre, EE, IIT Bombay Introduction In lecture 14, we studied steps required to design conjugate
More informationDiscontinuous Galerkin methods for fractional diffusion problems
Discontinuous Galerkin methods for fractional diffusion problems Bill McLean Kassem Mustapha School of Maths and Stats, University of NSW KFUPM, Dhahran Leipzig, 7 October, 2010 Outline Sub-diffusion Equation
More informationSemi-orthogonal wavelet frames on positive half-line using the Walsh Fourier transform
NTMSCI 6, No., 175-183 018) 175 New Trends in Mathematical Sciences http://dx.doi.org/10.085/ntmsci.018.83 Semi-orthogonal wavelet frames on positive half-line using the Walsh Fourier transform Abdullah
More informationSpotlight on Laplace s Equation
16 Spotlight on Laplace s Equation Reference: Sections 1.1,1.2, and 1.5. Laplace s equation is the undriven, linear, second-order PDE 2 u = (1) We defined diffusivity on page 587. where 2 is the Laplacian
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 informationHow smooth is the smoothest function in a given refinable space? Albert Cohen, Ingrid Daubechies, Amos Ron
How smooth is the smoothest function in a given refinable space? Albert Cohen, Ingrid Daubechies, Amos Ron A closed subspace V of L 2 := L 2 (IR d ) is called PSI (principal shift-invariant) if it is the
More informationChapter 7 Wavelets and Multiresolution Processing
Chapter 7 Wavelets and Multiresolution Processing Background Multiresolution Expansions Wavelet Transforms in One Dimension Wavelet Transforms in Two Dimensions Image Pyramids Subband Coding The Haar
More informationIntroduction to Signal Spaces
Introduction to Signal Spaces Selin Aviyente Department of Electrical and Computer Engineering Michigan State University January 12, 2010 Motivation Outline 1 Motivation 2 Vector Space 3 Inner Product
More informationWavelets and multiresolution representations. Time meets frequency
Wavelets and multiresolution representations Time meets frequency Time-Frequency resolution Depends on the time-frequency spread of the wavelet atoms Assuming that ψ is centred in t=0 Signal domain + t
More informationINVARIANCE OF A SHIFT-INVARIANT SPACE
INVARIANCE OF A SHIFT-INVARIANT SPACE AKRAM ALDROUBI, CARLOS CABRELLI, CHRISTOPHER HEIL, KERI KORNELSON, AND URSULA MOLTER Abstract. A shift-invariant space is a space of functions that is invariant under
More informationMATH FALL 2014 HOMEWORK 10 SOLUTIONS
Problem 1. MATH 241-2 FA 214 HOMEWORK 1 SOUTIONS Note that u E (x) 1 ( x) is an equilibrium distribution for the homogeneous pde that satisfies the given boundary conditions. We therefore want to find
More informationWavelets on Z N. Milos Savic
Student-Faculty Seminar Wavelets on Z N Milos Savic Milos Savic graduated in May 004 with a major in Mathematics and minors in Computer Science and Business. He has played soccer and water polo during
More informationMath 322. Spring 2015 Review Problems for Midterm 2
Linear Algebra: Topic: Linear Independence of vectors. Question. Math 3. Spring Review Problems for Midterm Explain why if A is not square, then either the row vectors or the column vectors of A are linearly
More informationMultiscale Frame-based Kernels for Image Registration
Multiscale Frame-based Kernels for Image Registration Ming Zhen, Tan National University of Singapore 22 July, 16 Ming Zhen, Tan (National University of Singapore) Multiscale Frame-based Kernels for Image
More information(3) Let Y be a totally bounded subset of a metric space X. Then the closure Y of Y
() Consider A = { q Q : q 2 2} as a subset of the metric space (Q, d), where d(x, y) = x y. Then A is A) closed but not open in Q B) open but not closed in Q C) neither open nor closed in Q D) both open
More informationIntroduction to Orthogonal Transforms. with Applications in Data Processing and Analysis
i Introduction to Orthogonal Transforms with Applications in Data Processing and Analysis ii Introduction to Orthogonal Transforms with Applications in Data Processing and Analysis October 14, 009 i ii
More informationMATH 423 Linear Algebra II Lecture 33: Diagonalization of normal operators.
MATH 423 Linear Algebra II Lecture 33: Diagonalization of normal operators. Adjoint operator and adjoint matrix Given a linear operator L on an inner product space V, the adjoint of L is a transformation
More informationAnalytical and Numerical Aspects of Wavelets
An- Najah National University Faculty of Graduate Studies Analytical and Numerical Aspects of Wavelets By Noora Hazem Abdel-Hamid Janem Supervisor Prof. Naji Qatanani This Thesis is Submitted in Partial
More informationPHZ 6607 Fall 2004 Homework #4, Due Friday, October 22, 2004
Read Chapters 9, 10 and 20. PHZ 6607 Fall 2004 Homework #4, Due Friday, October 22, 2004 1. The usual metric of four-dimensional flat Minkowski-space in spherical-polar coordinates is ds 2 = dt 2 + dr
More information# Points Score Total 100
Name: PennID: Math 241 Make-Up Final Exam January 19, 2016 Instructions: Turn off and put away your cell phone. Please write your Name and PennID on the top of this page. Please sign and date the pledge
More informationTwo-channel sampling in wavelet subspaces
DOI: 10.1515/auom-2015-0009 An. Şt. Univ. Ovidius Constanţa Vol. 23(1),2015, 115 125 Two-channel sampling in wavelet subspaces J.M. Kim and K.H. Kwon Abstract We develop two-channel sampling theory in
More informationLecture 2. Classification of Differential Equations and Method of Integrating Factors
Math 245 - Mathematics of Physics and Engineering I Lecture 2. Classification of Differential Equations and Method of Integrating Factors January 11, 2012 Konstantin Zuev (USC) Math 245, Lecture 2 January
More information