Slepian functions on the sphere, generalized Gaussian quadrature rule
|
|
- Whitney Higgins
- 5 years ago
- Views:
Transcription
1 INSTITUTE OF PHYSICS PUBLISHING Inverse Problems ( INVERSE PROBLEMS PII: S66-56( Slepian functions on the sphere, generalized Gaussian quadrature rule L Miranian Department of Mathematics, University of California, Berkeley CA, 947, USA luiza@math.berkeley.edu Received October 3, in final form March 4 Published April 4 Online at stacks.iop.org/ip//877 (DOI:.88/66-56//3/4 Abstract Denote by K the operator of time band time limiting on the surface of the sphere and consider the problem of computing singular vectors of K. This problem can be reduced to a simpler task of computing eigenfunctions of a differential operator, if a differential operator, which commutes with K and has a simple spectrum, can be exhibited. In Grünbaum et al (98 SIAM J. Appl. Math such a second-order differential operator commuting with K on the appropriate subspaces was constructed. In this paper, this algebraic property of commutativity is used to produce an efficient numerical scheme for computing a convenient basis for the space of singular vectors of K. The basis forms an extended Chebyshev system, and a generalized Gaussian quadrature rule for such a basis is presented.. Introduction The fundamental problem of recovering a time-limited function from the knowledge of its Fourier transform on a certain band of frequencies is a central chapter in signal processing. This problem plays an important role in many aspects of image processing since it underlies the question of how to make optimal use of the available information that is always limited and corrupted by noise. The remarkable series of papers [ 8], by Slepian, Landau and Pollak in connection with the issue of time band-limited signals has had a tremendous influence on many areas of engineering, science and mathematics. This work puts some of the pioneering work of Shannon [9] on firmer ground. Their starting points were fairly applied aspects of communication theory, optics, lasers, etc, but it became apparent that the ideas were applicable to many other situations. The work presented in this paper deals with the case when the real line is replaced by the surface of the sphere. Here, the mathematical and computational issue is to get good approximations to the appropriate Slepian functions. In this instance Slepian functions refer to a basis for the space of eigenfunctions of the operator K which is obtained by the 66-56/4/3877+6$3. 4 IOP Publishing Ltd Printed in the UK 877
2 878 L Miranian successive application of the operations of time, band and time limiting. In the classical case of the real line, the computation of Slepian functions is done using the fact that the classical differential operator, resulting from the Laplacian by separation of variables in prolate spheroidal coordinates, happens to commute with the time band-limiting integral operator. This commuting differential operator has a simple spectrum, hence its eigenfunctions form a basis for the space of eigenfunctions of the integral operator. In the case of the surface of the sphere, it was shown in [], that for a polar cap as well as for two symmetrically placed caps (one at each pole, a certain second-order differential operator D = d [ ( x (b x d ] L(L +x m (b x dx dx x defined on the interval [b, ] commutes with K on the spaces of functions whose dependence on φ is of the form e imφ. The operator D has a simple spectrum, hence its eigenfunctions, which are also eigenfunctions of K, can serve as a basis for the space of eigenfunctions of K. The same applies to the complement, in the sphere, of one or two polar caps. If the region in question has less symmetry, then one can always consider the integral operator, but the search for a commuting local operator has proved elusive. This algebraic property of existence of a commutative differential operator holds the key to a good algorithm. Since the object to be produced is a basis for the time band time-limited functions on a certain region of the sphere, one needs to expand easily in this basis, which requires an efficient quadrature rule for evaluating inner products. It is well known that for a system of functions that forms an extended Chebyshev system, a generalized Gaussian quadrature rule always exists. In [] a method of obtaining such a quadrature rule using the appropriate continuation scheme and well-chosen starting points for Newton s method is described. The operator K happens to be a finite rank Fredholm operator, which implies that it has only a finite number of non-zero eigenvalues. There are many ways to compute the null space of the operator K, but the algebraic property discussed above makes it possible to replace the computation of the eigenfunctions of K by the computation of the eigenfunctions of D. This not only simplifies the task from the numerical point of view, but also produces orthogonal functions which form an extended Chebyshev system, assuring the existence of an efficient quadrature rule. In order to numerically compute the eigenfunctions (of an appropriate self-adjoint extension of D we expand them in the basis of the shifted Legendre polynomials, and reduce the problem to the computation of generalized eigenvalues and eigenvectors of certain sparse matrices. We then use the method that has been advocated recently in []. The object produced is a good basis for the space of the eigenfunctions of K in the case when the regions of interest are either a polar cap or a spherical belt bounded by two parallels. Among the applications envisaged is the problem discussed in [] involving gravity field missions. Due to launch conditions or engineering reasons, the sampling is done not on the whole surface of the earth but on one of the type of regions mentioned above. The computation of the so-called Slepian functions, in this case, has been attempted in the geodesy community without taking recourse to the mathematical and computational advantages that are derived from exploiting the results in [] and the very recent note []. For other geodesy applications, see [3, 4]. In the first of these papers the region where the function is known is the surface of the oceans. This paper is organized as follows. In section. the problem of computing eigenfunctions of the integral operator K is discussed. In section. some properties of shifted Legendre polynomials are recalled. Sections.3 and.4 describe a method for computing Slepian
3 Slepian functions on the sphere, generalized Gaussian quadrature rule 879 functions using the differential operator. In section.5 the generalized Gaussian quadrature rule for the Slepian functions on the spherical cap is presented and concluding remarks are in section 3.. Numerical computation of the eigenproblem.. Direct computation of eigenfunctions of the integral operator K In this section, an attempt to compute eigenvalues/eigenvectors of the integral operator directly is described, and a more efficient alternative is suggested. As discussed in [], denote by A the polar cap θ arccos(b, φ π. Then the operator K is the finite convolution integral operator. Denote then (Kf (u = u = (sin θ cos φ,sin θ sin φ,cos θ, L P l ( u, u f (u du = A l= A ( L l= m= l where Y lm (u is the usual spherical harmonic with l Y lm (uy lm (u f(u du, and or π Y lm (x, φy lm (x, φ dx dφ =, l +(l m! Y lm (u = Y lm (cos θ,φ = 4π (l + m! P l m (cos θe imφ, l +(l m! Y lm (x, φ = 4π (l + m! P l m (x e imφ, x = cos θ. In the formulae above Pl m (x denotes the associated Legendre polynomial. The operator K is a singular Fredholm operator with rank (L +, hence it has only (L + non-zero eigenvalues. In the following proposition, some properties of the eigenvalues and eigenvectors of K are summarized. Proposition. Consider the finite (L + rank symmetric Fredholm operator as defined above and let H m be the subspaces of functions on the polar cap whose φ dependence is given by e imφ. Then (i There are L +linearly independent orthogonal eigenfunctions of K that belong to the space H ; consequently K has only L +distinct non-zero eigenvalues that correspond to eigenfunctions in H. (ii The L(L +/ non-zero eigenvalues of K have multiplicity : in both subspaces H m and H m K has a simple spectrum, i.e. L m +non-zero distinct eigenvalues, where m =,...,L. (iii L m +eigenfunctions of K belong to H m for all m =,...,L and are orthogonal.
4 88 L Miranian Proof. Let us fix L and m, and see what the integral operator looks like on H m, i.e. take g(u = f(cos θe i mφ = f(xe i mφ and apply the operator K to it: ( L l Kf(u= Y lm (uy lm (u g(u du A l= m= l π ( L l = Y lm (x, φy lm (x,φ f(x e i mφ dx dφ b l= m= l π ( L l l +(l m! = 4π (l + m! P l m (x e imφ Pl m (x e imφ f(x e i mφ dx dφ = b b L l= m= l l= m= l l l +(l m! 4π (l + m! P l m L = e i mφ l +(l m! b π l= m To simplify the notation denote (l + m! P l m L K(m,x,x l +(l m! = π l= m as the kernel of K in subspace H m. The following observations can be made: π (xpl m (x f (x (e imφ e iφ ( m+ m dφ dx (xp m l (x f (x dx. (l + m! P l m (xp m l (x (i The associated Legendre polynomials Pl m are linearly independent and K(m,x,x = K(m,x,xdefines a singular symmetric Fredholm operator which has only L m + distinct non-zero eigenvalues for all m =,...,L. (ii If m =, then K has only L + distinct non-zero eigenvalues in H, and because of the symmetry of the kernel K(,x,x the corresponding eigenvectors are orthogonal. (iii K(m,x,x = K( m, x, x. The symmetry of the kernel K(m,x,x implies that the eigenvectors of K in H m corresponding to different eigenvalues are orthogonal. (iv If m, then in both H m and H m the operator K has the same kernel, hence L m + eigenvalues of K will be duplicated for every m =,...,L. Attempts to compute the eigenfunctions of the integral operator KF n (u = µ n F n (u directly have not been fruitful. In the experiments discussed below the integral operator with L = 3 was discretized using the Gaussian quadrature rule with N x = -point grid in x = cos(θ variable and N φ = -point grid in φ variable. Let T be the matrix obtained after discretization of the integral operator K. The disadvantages of this method are: (i The size of the matrix T is N x N φ = 44; it depends quadratically on the grid size, which makes computing its eigenvalues/eigenvectors an intensive task. In the alternative approach a grid of 4 points in the x variable is used, whereas in the direct approach it is computationally infeasible because of the reason just mentioned. (ii Only eigenfunctions corresponding to non-zero eigenvalues could be computed directly, hence a procedure that would produce an orthogonal basis for the null space of K is needed. Below is the summary of some numerical experiments. Only the eigenfunctions of T that correspond to non-zero eigenvalues are meaningful, so only these are considered in the text below.
5 Slepian functions on the sphere, generalized Gaussian quadrature rule Figure. Results of the direct discretization (four left-hand figures and eigenfunctions obtained by using the differential operator D (four right-hand figures for m =,,, 3. Functions presented are numbered according to the number of roots, i.e. G 3,G,G, G ; L = 3. (i Exactly four eigenfunctions of T are in subspace H ; three eigenfunctions in H ± ; two eigenfunctions in H ± and one eigenfunction in H ±3. (ii Eigenfunctions F j (x e ±imφ with j =,...,L m +oft in each H ±m are very close to these produced by the alternative procedure below (see figure.
6 88 L Miranian (iii Eigenfunctions of T corresponding to the eigenvalue cannot be computed by direct discretization. There is an efficient alternative to the procedure described above. It produces an orthonormal basis for the space of eigenfunctions of K. Consider the following second-order differential operator on the interval [b, ]: D = d dx [ ( x (b x d dx ] L(L +x m (b x x. A look at [] will show that this is the appropriate operator D that commutes with the K built there when acting on H m. The operator D defines a Sturm Liouville problem, hence it has a simple spectrum and orthogonal eigenfunctions. Because it commutes with K we can say that the eigenfunctions of D are also eigenfunctions of K when acting on H m. Since we search for an orthonormal basis in the space of eigenfunctions of K, the eigenfunctions of D provide us with such a basis. In order to compute numerically the eigenfunctions (of an appropriate self-adjoint extension of D we expand them in the basis of the shifted Legendre polynomials, and reduce the problem to (generalized eigenproblem for some well-structured matrices... Shifted Legendre polynomials In this section, an overview of some facts about shifted Legendre polynomials is presented. These facts will be used to design a procedure for computing eigenvalues/eigenvectors of the differential operator D. Define shifted Legendre polynomials to be the solutions of the following second-order differential equation: (b x( xs n +(x b S n n(n +S n =. Denote b := (+b/,b := ( b/. The following properties of S n will be useful later: (i recursion relation S n = b S n + b (n + n + S n+ + b n n + S n ; (ii derivative ( x(b xs n = b n(n + n + (S n+ S n ; (iii normalized shifted Legendre polynomials S k S k (k +/b..3. Computation of the eigenproblem DF n = λ n F n : case m = In this section, the problem of computing the eigenfunctions of the differential operator D with m = is reduced to the problem of computing eigenvectors of a certain symmetric tridiagonal matrix. Consider the following eigenproblem: ( d dx [ ( x (b x d dx ] L(L +x F n = λ n F n. (
7 Slepian functions on the sphere, generalized Gaussian quadrature rule 883 Let ( a n,an,an,... be the coefficients of the expansion of F n (x in the basis of shifted Legendre polynomials F n = ak n S k. k= After substituting F n (x into (, using properties (i, (ii and the linear independence of S k one obtains a recursion relation C k a k + B k+ a k+ + A k a k λ n a k =, ( where A k = k(k +(+b L(L +b B k = k k + b [(k (k + L(L +] C k = k + k + b [k(k + L(L +]. After rewriting ( for the normalized polynomials using property (iii, the recursion relation can be written in the matrix form as Mā n = λ n ā n, (3 where ā n = ( ā n, ān,... T are the coefficients of the expansion in the basis of the normalized shifted Legendre polynomials; M k,k = k(k +(+b L(L +b, b (k + M k,k+ = [k(k + L(L +], (k +3(k + M k+,k = b (k + (k +3(k + [k(k + L(L +], k =,,... and the remainder of the entries of the matrix being zero. Using ā n obtained from (3, we express F n = k= ān S k k. The numerical evidence strongly suggests that the coefficients āk n decay very fast. In practice, this allows us to compute F n = N k= ān S k k,for certain large values of N..4. Computation of the eigenproblem DG n = µ n G n : case m> In this section, the problem of computing eigenfunctions of the differential operator D with m>isreduced to a generalized matrix eigenproblem. Consider the following eigenproblem: [( x (b xg n ] L(L +xg n m (b x G x n = µ n G n. (4 Introduce the following simplifying notation: (i Matrix S = ( S, S, S 3,..., where S k are normalized shifted Legendre polynomials. (ii Matrix A with columns A k = ( ā k, āk, āk,... T, where ā k j are the coefficients of the expansion of F k (x in the basis of normalized shifted Legendre polynomials. (iii F k = S j āj k = SAk. j=
8 884 L Miranian (iv Also, xf k = = xs j aj k = j= j= ā k j = SQA k, j= ( b S j + a k j ( b S j + b (j + j + S j+ + b j j + S j b (j + (j +(j +3 S j+ + where the symmetric tridiagonal matrix Q has entries Q k,k = b, Q k,k+ = Q k+,k = b (k + (k +3(k +, with =,,,... Similarly x F k = SQ A k. (v Denote c n = ( c n,cn,cn,... T. (vi Matrix = diag(λ,λ,...; A = MA. (vii Denote I to be the identity matrix. b j S j (j +(j Since the functions F k form an orthonormal basis for L [b, ], one could expand G n in the basis of F k, i.e. G n = ck n F k(x = ck n SAk = SAc n = γk n S k (x, (5 k= k= where γ n = (γ n,γn,...t = Ac n. After substituting (5 into(4 and using the fact that F k satisfies ( one obtains c n k (λ k F k m (b x F x k µ n F k =, or ck n (λ kf k ( x m (b xf k µ n F k ( x =. (6 k= Using the notation above, (6 can be written as ck n (λ ks(i Q A k m S(b QA k µ n S(I Q A k k= = S ck n (λ k(i Q A k m (b QA k µ n (I Q A k =. k= Since the functions S k are linearly independent, one obtains = ck n (λ k(i Q A k m (b QA k µ n (I Q A k k= k= = (I Q A c n m (b QAc n µ n (I Q Ac n = (I Q MAc n m (b QAc n µ n (I Q Ac n = ((I Q M m (b Q µ n (I Q Ac n,
9 Slepian functions on the sphere, generalized Gaussian quadrature rule 885 which implies the generalized eigenproblem (I Q (M µ n Iγ n = m (b Qγ n. (7 Although the coefficients γ n can be computed from the elegant matrix eigenproblem (7, numerical experiments have shown that the coefficients of the expansion G n (x = k= γ k n S k (x decay slowly. However, the scheme above can be significantly improved, as we explain now. Taking into account boundary conditions, write G n = ( x m/ g n (x, for a certain function g n (x. After rewriting (4 in terms of g n (x we arrive at the following eigenproblem: [( x (b xg n ] +mx(x bg n + ((m +m L Lx mb(m +g n = µ n g n. (8 Now, the scheme described above can be applied to the function g n (x = k= cn k F k for some coefficients ck n using the differential equation (8. In this case the derivation of the generalized matrix eigenproblem is similar to that performed at the beginning of this section. After elaborate calculations one arrives at (I Q(M µ n Iα n = (mb(m +(I Q +mq 3 m(m +(Q Q α n, (9 where b Q 3 (k, k = (k +3(k, b Q 3 (k +,k= Q 3 (k, k + = (k +3(k, b (k + Q 3 (k, k + = (k +3 (k +(k +5, b (k + Q 3 (k +,k= (k +3 (k +(k +5, E(k, k = b k(k +, Q 3 = Q 3 E for k =,,, 3,...; matrices Q and M were defined before, and Ac n = α n. Observe that g n (x = ck n F k(x = ck n SAk = S ck n Ak = SAc n = Sα n, k= k= which means that while c n are coefficients of the expansion of g n (x in the basis of F k (x s, ᾱ n s are coefficients of the expansion of g n (x in the basis of normalized shifted Legendre polynomials S k s. From (9 we can compute g n = k= αn S k k, and obtain G n (x = ( x m/ g n (x. Experiments suggest that coefficients of the expansions of g n = k= αn S k k decay very rapidly. Moreover, it is not hard to see that if G n (x = ( x m/ g n (x = ( x m/ k= αn S k k is an eigenfunction of K that corresponds to a non-zero eigenvalue, then ᾱp n = for all p>l m. In order to observe this recall that on the subspace H m kernel of K is K(x,x = ( x m/ ( x m/ Z L m (x, x, where Z L m (x, x is a symmetric polynomial in x and x of degree L m. Then, α n S p (xg n (x ( p = dx b ( x m/ = S p (x K(m,x,x G b ( x m/ n (x dx dx ξ n b = ( S ξ n G n (x p (x ( x m/ K(m,x,x dx dx b b k=
10 886 L Miranian ( K(m,x,x ξ n b b ( x S m/ p (x dx dx = ( ( x ξ n m/ Z L m (x, x S p (x dx dx = b b for all p>l m, since the Legendre polynomial S p (x of degree p is orthogonal to x k for all k<p. Because of the very rapid decay of the coefficients αk n, eigenfunctions G n can be computed efficiently. It is very important to sort eigenvalues (along with corresponding eigenvectors of (9 in the ascending order before computing the sum g n = N k= ān S k k for some appropriate finite N. In figure one can see the decay of the coefficients of the expansion in the case m =,L = and b = for the eigenfunctions G 5,G, G 3,G 8 ; m = 4,L = 7 and b = /for the eigenfunction G 5,G, G 3,G 8. In figure 3 eigenfunctions G k of D are presented, for k =, 5,,, 35, 5,L = 5, m =,b = /..5. Construction of generalized Gaussian quadratures for the eigenfunctions of D In this section, an algorithm for constructing the generalized Gaussian quadrature rule for the eigenfunctions of D is described and results of some numerical experiments are presented. Functions G n (x form a complete orthonormal basis for the space L [b, ], and a quadrature rule for computing integrals of the form b f(xg n(x dx efficiently for various functions f(xis needed. In this instance, an efficient quadrature rule refers to a quadrature such that the error decreases at least exponentially as a function of the number of nodes used in the integration. The eigenfunctions G n form an extended Chebyshev system, and according to the principal result of [6], there exists a unique n-point generalized Gaussian quadrature rule with weight W :[b, ] R +. Nodes (x,...,x n and weights (w,...,w n of the quadrature satisfy a system of nonlinear equations n w i G j (x i = G j (xw(x dx; j =,...,n. ( i= b Newton s method for this system of equations is always quadratically convergent, since the fact that eigenfunctions form an extended Chebyshev system implies the Jacobian of the system being always nonsingular ([5, lemma.6]. In order to provide a good starting point for Newton s method we use the continuation scheme suggested in [5]. A numerical algorithm for obtaining weights and nodes of n-point generalized Gaussian quadrature is described below. (i As a starting point of the algorithm take the n roots of G n, denote them as r = (r,...,r n. (ii Use a continuation scheme, i.e. let the weight functions ω: [, ] [b, ] R + be defined by the formula n ω(α,x = αw(x + ( α δ(x r j, j= where δ denotes the Dirac delta function. Observe that when α =, the weight function is equal to the desired weight function W(x, and when α = the Gaussian weights and nodes are w i =,x i = r i for i =,...,n.
11 Slepian functions on the sphere, generalized Gaussian quadrature rule Figure. Top four plots: L =,m =,b = ; bottom four plots: L = 7,m = 4,b = /; graphs correspond to absolute values of the first 35 coefficients of the expansions for eigenfunctions G 5,G,G 3,G 8 versus their index. (iii Damped Newton s method (at every iteration of Newton s method search for an appropriate step size along the direction prescribed by regular Newton s method is used to solve the system on every step of the continuation scheme.
12 888 L Miranian Figure 3. Eigenfunctions G k of D for k =, 5,,, 35, 5; b = /,L= 5,m=. In the numerical experiments conducted using MATLAB, the quadrature rule ( for various weight functions was obtained. In particular, quadratures with W(x =, x b were constructed and the resulting weights and nodes are presented in tables and. To see the accuracy of the quadratures, the test function f(x= sin 5 (x cos(7x( x, for instance, is integrated with relative and absolute errors of 7.8 and.4 6 with a 5-point quadrature, where W(x =. The test function f(x= x (x b 3/ ( x 5/ is integrated with relative and absolute errors of 5.65 and.4 with a -point quadrature, where W(x = x b.
13 Slepian functions on the sphere, generalized Gaussian quadrature rule 889 Table. Nodes and weights of a n = 5-, -, -point generalized Gaussian quadrature; L = 4; m = ; b = ; W(x =. N Nodes x i Weights w i On average, the algorithm described above performs only one step of the continuation scheme. On every step of the continuation scheme it does about seven steps of Newton s iteration with around three step-size adjustments per each iteration. The functions G n e imφ form a complete, orthonormal basis for functions on the spherical cap, hence we need to compute double integrals in colatitudinal and longitudinal variables. The procedure described above produces the nodes and weights for the colatitudinal variable. Integration with respect to the longitudinal variable can be done by a Gaussian quadrature as well, so the final quadrature rule is π n φ ( n f(x,φdx dφ = f(x k,φ i wk x, b i= w φ i k=
14 89 L Miranian Table. Nodes and weights of a n = 5-, -, -point generalized Gaussian quadrature; L = 3; m = ; b = ; W(x = x b. N Nodes x i Weights w i where nodes and weights x k,wk x are products of the scheme presented before, and nodes and weights φ i,w φ i are those that correspond to the quadrature rule for longitudinal variable with n φ being the desired number of nodes. Figures 4 and 5 show the nodes of the 4-point and -point quadratures on the sphere, where L =,m =,b = /and L = 5,m =,b = / correspondingly. The nodes tend to concentrate more at the bottom of the spherical cap rather than around the pole. A consequence of this fact is that the quadrature scheme described in this section does not produce the familiar north pole oversampling problem.
15 Slepian functions on the sphere, generalized Gaussian quadrature rule Figure 4. Nodes for 4-point quadrature, b = /,L=,m= Figure 5. Nodes for -point quadrature, b = /,L= 5,m=. 3. Conclusions An efficient numerical scheme for evaluating a basis for the set of singular vectors of a time band time-limiting operator for certain regions on the surface has been presented. The basis functions form an orthonormal set, as well as an extended Chebyshev system, which guarantees the existence of a generalized Gaussian quadrature rule. An algorithm for computing weights and nodes for such a quadrature rule is presented. The nodes produced tend to concentrate not at the pole of the sphere, but at the bottom of the spherical cap, not exhibiting a common problematic effect called north pole oversampling.
16 89 L Miranian Acknowledgments The author is very grateful to Professor F A Grünbaum for many invaluable discussions. The author would also like to thank Professor W Kahan for a discussion on the Gaussian quadrature rule and for providing very useful references on the subject. References [] Grünbaum F A, Longhi L and Perlstadt M 98 Differential operators commuting with finite convolution integral operators: some non-abelian examples SIAM J. Appl. Math [] Slepian D and Pollak H O 96 Prolate spheroidal wave functions, Fourier analysis and uncertainty, I Bell Syst. Tech. J [3] Landau H J and Pollak H O 96 Prolate spheroidal wave functions, Fourier analysis and uncertainty, II Bell Syst. Tech. J [4] Landau H J and Pollak H O 96 Prolate spheroidal wave functions, Fourier analysis and uncertainty, III Bell Syst. Tech. J [5] Slepian D 964 Prolate spheroidal wave functions, Fourier analysis and uncertainty, IV Bell Syst. Tech. J [6] Slepian D 983 Some comments on Fourier analysis, uncertainty and modeling SIAM Rev [7] Slepian D 976 On bandwidth Proc. IEEE [8] Slepian D 978 Prolate spheroidal wave functions, Fourier analysis and uncertainty, V Bell Syst. Tech. J [9] Shannon C E 948 A mathematical theory of communication Bell Syst. Tech. J Shannon C E 948 Bell Syst. Tech. J [] Xiao H, Rokhlin V and Yarvin N Prolate spheroidal wave functions, quadrature, and interpolation Inverse Problems [] Albertella A, Sanso F and Sneeuw N 999 Band-limited functions on a bounded spherical domain: the Slepian problem on the sphere J. Geodesy [] Grünbaum F A and Miranian L The magic of the prolate spheroidal functions in various setups Proc. SPIE [3] Hwang Ch 99 Spectral analysis using orthonormal functions with a case study on the sea surface topography Geophys. J. Int [4] Albertella A and Sneeuw N The analysis of radiometric data with Slepian functions Phys. Chem. Earth A [5] Cheng H, Rokhlin V and Yarvin N 999 Nonlinear optimization, quadrature, and interpolation SIAM. J. Optim [6] Karlin S and Studden W 966 Tchebyshev Systems with Applications in Analysis and Statistics (New York: Wiley-Interscience
The magic of the prolate spheroidal functions in various setups
The magic of the prolate spheroidal functions in various setups F. A. Grflnbaum and L. Miranian Department of Mathematics, University of California, Berkeley CA, 94720 ABSTRACT The singular functions for
More informationSlepian Functions Why, What, How?
Slepian Functions Why, What, How? Lecture 1: Basics of Fourier and Fourier-Legendre Transforms Lecture : Spatiospectral Concentration Problem Cartesian & Spherical Cases Lecture 3: Slepian Functions and
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 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 informationThe Slepian approach revisited:
he Slepian approach revisited: dealing with the polar gap in satellite based geopotential recovery Oliver Baur, Nico Sneeuw Geodätisches Institut Universität Stuttgart 3 rd GOCE user workshop 6-8 Noveber
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 informationFast Algorithms for Spherical Harmonic Expansions
Fast Algorithms for Spherical Harmonic Expansions Vladimir Rokhlin and Mark Tygert Research Report YALEU/DCS/RR-1309 December 17, 004 Abstract An algorithm is introduced for the rapid evaluation at appropriately
More informationSpherical Harmonics on S 2
7 August 00 1 Spherical Harmonics on 1 The Laplace-Beltrami Operator In what follows, we describe points on using the parametrization x = cos ϕ sin θ, y = sin ϕ sin θ, z = cos θ, where θ is the colatitude
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 information256 Summary. D n f(x j ) = f j+n f j n 2n x. j n=1. α m n = 2( 1) n (m!) 2 (m n)!(m + n)!. PPW = 2π k x 2 N + 1. i=0?d i,j. N/2} N + 1-dim.
56 Summary High order FD Finite-order finite differences: Points per Wavelength: Number of passes: D n f(x j ) = f j+n f j n n x df xj = m α m dx n D n f j j n= α m n = ( ) n (m!) (m n)!(m + n)!. PPW =
More informationc 2006 Society for Industrial and Applied Mathematics
SIAM J. SCI. COMPUT. Vol. 7, No. 6, pp. 1903 198 c 006 Society for Industrial and Applied Mathematics FAST ALGORITHMS FOR SPHERICAL HARMONIC EXPANSIONS VLADIMIR ROKHLIN AND MARK TYGERT Abstract. An algorithm
More informationAngular Momentum. Classically the orbital angular momentum with respect to a fixed origin is. L = r p. = yp z. L x. zp y L y. = zp x. xpz L z.
Angular momentum is an important concept in quantum theory, necessary for analyzing motion in 3D as well as intrinsic properties such as spin Classically the orbital angular momentum with respect to a
More informationApplied Linear Algebra
Applied Linear Algebra Peter J. Olver School of Mathematics University of Minnesota Minneapolis, MN 55455 olver@math.umn.edu http://www.math.umn.edu/ olver Chehrzad Shakiban Department of Mathematics University
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 information1 Distributions (due January 22, 2009)
Distributions (due January 22, 29). The distribution derivative of the locally integrable function ln( x ) is the principal value distribution /x. We know that, φ = lim φ(x) dx. x ɛ x Show that x, φ =
More informationd 1 µ 2 Θ = 0. (4.1) consider first the case of m = 0 where there is no azimuthal dependence on the angle φ.
4 Legendre Functions In order to investigate the solutions of Legendre s differential equation d ( µ ) dθ ] ] + l(l + ) m dµ dµ µ Θ = 0. (4.) consider first the case of m = 0 where there is no azimuthal
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 informationPreliminary/Qualifying Exam in Numerical Analysis (Math 502a) Spring 2012
Instructions Preliminary/Qualifying Exam in Numerical Analysis (Math 502a) Spring 2012 The exam consists of four problems, each having multiple parts. You should attempt to solve all four problems. 1.
More informationIndex. C 2 ( ), 447 C k [a,b], 37 C0 ( ), 618 ( ), 447 CD 2 CN 2
Index advection equation, 29 in three dimensions, 446 advection-diffusion equation, 31 aluminum, 200 angle between two vectors, 58 area integral, 439 automatic step control, 119 back substitution, 604
More informationArchiv der Mathematik Holomorphic approximation of L2-functions on the unit sphere in R^3
Manuscript Number: Archiv der Mathematik Holomorphic approximation of L-functions on the unit sphere in R^ --Manuscript Draft-- Full Title: Article Type: Corresponding Author: Holomorphic approximation
More informationFast Multipole Methods for The Laplace Equation. Outline
Fast Multipole Methods for The Laplace Equation Ramani Duraiswami Nail Gumerov Outline 3D Laplace equation and Coulomb potentials Multipole and local expansions Special functions Legendre polynomials Associated
More informationarxiv: v1 [math.ca] 21 Mar 2014
SLEPIAN SPATIAL-SPECTRAL CONCENTRATION ON THE BALL ZUBAIR KHALID, RODNEY A. KENNEDY, AND JASON D. MCEWEN arxiv:1403.5553v1 [math.ca] 21 Mar 2014 Abstract. We formulate and solve the Slepian spatial-spectral
More informationTotal least squares. Gérard MEURANT. October, 2008
Total least squares Gérard MEURANT October, 2008 1 Introduction to total least squares 2 Approximation of the TLS secular equation 3 Numerical experiments Introduction to total least squares In least squares
More informationClassical Field Theory: Electrostatics-Magnetostatics
Classical Field Theory: Electrostatics-Magnetostatics April 27, 2010 1 1 J.D.Jackson, Classical Electrodynamics, 2nd Edition, Section 1-5 Electrostatics The behavior of an electrostatic field can be described
More informationAPPENDIX A. Background Mathematics. A.1 Linear Algebra. Vector algebra. Let x denote the n-dimensional column vector with components x 1 x 2.
APPENDIX A Background Mathematics A. Linear Algebra A.. Vector algebra Let x denote the n-dimensional column vector with components 0 x x 2 B C @. A x n Definition 6 (scalar product). The scalar product
More informationON INTERPOLATION AND INTEGRATION IN FINITE-DIMENSIONAL SPACES OF BOUNDED FUNCTIONS
COMM. APP. MATH. AND COMP. SCI. Vol., No., ON INTERPOLATION AND INTEGRATION IN FINITE-DIMENSIONAL SPACES OF BOUNDED FUNCTIONS PER-GUNNAR MARTINSSON, VLADIMIR ROKHLIN AND MARK TYGERT We observe that, under
More informationBefore you begin read these instructions carefully:
NATURAL SCIENCES TRIPOS Part IB & II (General) Friday, 31 May, 2013 9:00 am to 12:00 pm MATHEMATICS (2) Before you begin read these instructions carefully: You may submit answers to no more than six questions.
More informationSPECTRAL METHODS: ORTHOGONAL POLYNOMIALS
SPECTRAL METHODS: ORTHOGONAL POLYNOMIALS 31 October, 2007 1 INTRODUCTION 2 ORTHOGONAL POLYNOMIALS Properties of Orthogonal Polynomials 3 GAUSS INTEGRATION Gauss- Radau Integration Gauss -Lobatto Integration
More informationTheory and Computation for Bilinear Quadratures
Theory and Computation for Bilinear Quadratures Christopher A. Wong University of California, Berkeley 15 March 2015 Research supported by NSF, AFOSR, and SIAM C. Wong (UC Berkeley) Bilinear Quadratures
More informationRecall that any inner product space V has an associated norm defined by
Hilbert Spaces Recall that any inner product space V has an associated norm defined by v = v v. Thus an inner product space can be viewed as a special kind of normed vector space. In particular every inner
More informationOn the Analytical and Numerical Properties of the Truncated Laplace Transform II.
The Laplace Transform is frequently encountered in mathematics, physics, engineering and other fields. However, the spectral properties of the Laplace Transform tend to complicate its numerical treatment;
More informationMath 307 Learning Goals. March 23, 2010
Math 307 Learning Goals March 23, 2010 Course Description The course presents core concepts of linear algebra by focusing on applications in Science and Engineering. Examples of applications from recent
More informationWAVELET MULTIPLIERS AND SIGNALS
J. Austral. Math. Soc. Sen B 40(1999), 437-446 WAVELET MULTIPLIERS AND SIGNALS Z. HE 1 and M. W. WONG 1 (Received 26 March 1996; revised 6 August 1996) Abstract The Schatten-von Neumann property of a pseudo-differential
More informationEXAM MATHEMATICAL METHODS OF PHYSICS. TRACK ANALYSIS (Chapters I-V). Thursday, June 7th,
EXAM MATHEMATICAL METHODS OF PHYSICS TRACK ANALYSIS (Chapters I-V) Thursday, June 7th, 1-13 Students who are entitled to a lighter version of the exam may skip problems 1, 8-11 and 16 Consider the differential
More informationSolving the 3D Laplace Equation by Meshless Collocation via Harmonic Kernels
Solving the 3D Laplace Equation by Meshless Collocation via Harmonic Kernels Y.C. Hon and R. Schaback April 9, Abstract This paper solves the Laplace equation u = on domains Ω R 3 by meshless collocation
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 informationLinear Algebra. and
Instructions Please answer the six problems on your own paper. These are essay questions: you should write in complete sentences. 1. Are the two matrices 1 2 2 1 3 5 2 7 and 1 1 1 4 4 2 5 5 2 row equivalent?
More informationACM/CMS 107 Linear Analysis & Applications Fall 2017 Assignment 2: PDEs and Finite Element Methods Due: 7th November 2017
ACM/CMS 17 Linear Analysis & Applications Fall 217 Assignment 2: PDEs and Finite Element Methods Due: 7th November 217 For this assignment the following MATLAB code will be required: Introduction http://wwwmdunloporg/cms17/assignment2zip
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 informationSolutions Preliminary Examination in Numerical Analysis January, 2017
Solutions Preliminary Examination in Numerical Analysis January, 07 Root Finding The roots are -,0, a) First consider x 0 > Let x n+ = + ε and x n = + δ with δ > 0 The iteration gives 0 < ε δ < 3, which
More informationA fast randomized algorithm for overdetermined linear least-squares regression
A fast randomized algorithm for overdetermined linear least-squares regression Vladimir Rokhlin and Mark Tygert Technical Report YALEU/DCS/TR-1403 April 28, 2008 Abstract We introduce a randomized algorithm
More informationOn classical orthogonal polynomials and differential operators
INSTITUTE OF PHYSICS PUBLISHING JOURNAL OF PHYSICS A: MATHEMATICAL AND GENERAL J. Phys. A: Math. Gen. 38 2005) 6379 6383 oi:10.1088/0305-4470/38/28/010 On classical orthogonal polynomials an ifferential
More informationLecture IX. Definition 1 A non-singular Sturm 1 -Liouville 2 problem consists of a second order linear differential equation of the form.
Lecture IX Abstract When solving PDEs it is often necessary to represent the solution in terms of a series of orthogonal functions. One way to obtain an orthogonal family of functions is by solving a particular
More informationADVANCED ENGINEERING MATHEMATICS MATLAB
ADVANCED ENGINEERING MATHEMATICS WITH MATLAB THIRD EDITION Dean G. Duffy Contents Dedication Contents Acknowledgments Author Introduction List of Definitions Chapter 1: Complex Variables 1.1 Complex Numbers
More informationSpherical Slepian functions and the polar gap in geodesy
Geophys. J. Int. (6) 166, 139 161 doi: 1.1111/j.1365-46X.6.365.x Spherical Slepian functions and the polar gap in geodesy Frederik J. Simons 1 and F. A. Dahlen 1 Department of Earth Sciences, University
More informationyou expect to encounter difficulties when trying to solve A x = b? 4. A composite quadrature rule has error associated with it in the following form
Qualifying exam for numerical analysis (Spring 2017) Show your work for full credit. If you are unable to solve some part, attempt the subsequent parts. 1. Consider the following finite difference: f (0)
More informationFixed point iteration and root finding
Fixed point iteration and root finding The sign function is defined as x > 0 sign(x) = 0 x = 0 x < 0. It can be evaluated via an iteration which is useful for some problems. One such iteration is given
More informationBASIC EXAM ADVANCED CALCULUS/LINEAR ALGEBRA
1 BASIC EXAM ADVANCED CALCULUS/LINEAR ALGEBRA This part of the Basic Exam covers topics at the undergraduate level, most of which might be encountered in courses here such as Math 233, 235, 425, 523, 545.
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 informationLecture 11 Spin, orbital, and total angular momentum Mechanics. 1 Very brief background. 2 General properties of angular momentum operators
Lecture Spin, orbital, and total angular momentum 70.00 Mechanics Very brief background MATH-GA In 9, a famous experiment conducted by Otto Stern and Walther Gerlach, involving particles subject to a nonuniform
More informationKasetsart University Workshop. Multigrid methods: An introduction
Kasetsart University Workshop Multigrid methods: An introduction Dr. Anand Pardhanani Mathematics Department Earlham College Richmond, Indiana USA pardhan@earlham.edu A copy of these slides is available
More informationQuantum Mechanics Solutions
Quantum Mechanics Solutions (a (i f A and B are Hermitian, since (AB = B A = BA, operator AB is Hermitian if and only if A and B commute So, we know that [A,B] = 0, which means that the Hilbert space H
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 informationCompleteness of the Dirac oscillator eigenfunctions
INSTITUTE OF PHYSICS PUBLISHING JOURNAL OF PHYSICS A: MATHEMATICAL AND GENERAL J. Phys. A: Math. Gen. 34 (001) 4991 4997 www.iop.org/journals/ja PII: S0305-4470(01)011-X Completeness of the Dirac oscillator
More informationLecture 10: October 27, 2016
Mathematical Toolkit Autumn 206 Lecturer: Madhur Tulsiani Lecture 0: October 27, 206 The conjugate gradient method In the last lecture we saw the steepest descent or gradient descent method for finding
More informationLectures 15: Cayley Graphs of Abelian Groups
U.C. Berkeley CS294: Spectral Methods and Expanders Handout 15 Luca Trevisan March 14, 2016 Lectures 15: Cayley Graphs of Abelian Groups In which we show how to find the eigenvalues and eigenvectors of
More informationMATH 31 - ADDITIONAL PRACTICE PROBLEMS FOR FINAL
MATH 3 - ADDITIONAL PRACTICE PROBLEMS FOR FINAL MAIN TOPICS FOR THE FINAL EXAM:. Vectors. Dot product. Cross product. Geometric applications. 2. Row reduction. Null space, column space, row space, left
More informationM.A. Botchev. September 5, 2014
Rome-Moscow school of Matrix Methods and Applied Linear Algebra 2014 A short introduction to Krylov subspaces for linear systems, matrix functions and inexact Newton methods. Plan and exercises. M.A. Botchev
More informationExtrapolation of Bandlimited Signals Using Slepian Functions
, 23-25 October, 2013, San Francisco, USA Extrapolation of Bandlimited Signals Using Slepian Functions Amal Devasia 1, Member, IAENG and Michael Cada 2 Abstract An efficient and reliable method to extrapolate
More informationElectromagnetism HW 1 math review
Electromagnetism HW math review Problems -5 due Mon 7th Sep, 6- due Mon 4th Sep Exercise. The Levi-Civita symbol, ɛ ijk, also known as the completely antisymmetric rank-3 tensor, has the following properties:
More informationParameter Selection Techniques and Surrogate Models
Parameter Selection Techniques and Surrogate Models Model Reduction: Will discuss two forms Parameter space reduction Surrogate models to reduce model complexity Input Representation Local Sensitivity
More informationDetailed analysis of prolate quadratures and interpolation formulas
As demonstrated by Slepian et. al. in a sequence of classical papers see [33], [34], [7], [35], [36]), prolate spheroidal wave functions PSWFs) provide a natural and efficient tool for computing with bandlimited
More informationClass Averaging in Cryo-Electron Microscopy
Class Averaging in Cryo-Electron Microscopy Zhizhen Jane Zhao Courant Institute of Mathematical Sciences, NYU Mathematics in Data Sciences ICERM, Brown University July 30 2015 Single Particle Reconstruction
More informationIntroduction to Mathematical Physics
Introduction to Mathematical Physics Methods and Concepts Second Edition Chun Wa Wong Department of Physics and Astronomy University of California Los Angeles OXFORD UNIVERSITY PRESS Contents 1 Vectors
More informationSummary: angular momentum derivation
Summary: angular momentum derivation L = r p L x = yp z zp y, etc. [x, p y ] = 0, etc. (-) (-) (-3) Angular momentum commutation relations [L x, L y ] = i hl z (-4) [L i, L j ] = i hɛ ijk L k (-5) Levi-Civita
More informationSupporting Information
Supporting Information Harig and Simons 1.173/pnas.178519 SI Text Determination of Noise. Gravity Recovery and Climate Experiment (GRACE) data are released as spherical harmonic coefficients along with
More informationMath 102, Winter Final Exam Review. Chapter 1. Matrices and Gaussian Elimination
Math 0, Winter 07 Final Exam Review Chapter. Matrices and Gaussian Elimination { x + x =,. Different forms of a system of linear equations. Example: The x + 4x = 4. [ ] [ ] [ ] vector form (or the column
More informationClass notes: Approximation
Class notes: Approximation Introduction Vector spaces, linear independence, subspace The goal of Numerical Analysis is to compute approximations We want to approximate eg numbers in R or C vectors in R
More informationEigenfunctions on the surface of a sphere. In spherical coordinates, the Laplacian is. u = u rr + 2 r u r + 1 r 2. sin 2 (θ) + 1
Eigenfunctions on the surface of a sphere In spherical coordinates, the Laplacian is u = u rr + 2 r u r + 1 r 2 [ uφφ sin 2 (θ) + 1 sin θ (sin θ u θ) θ ]. Eigenfunctions on the surface of a sphere In spherical
More informationIntroduction. Chapter One
Chapter One Introduction The aim of this book is to describe and explain the beautiful mathematical relationships between matrices, moments, orthogonal polynomials, quadrature rules and the Lanczos and
More informationarxiv:math/ v1 [math.st] 11 Mar 2006
Geophys. J. Int. (XXXX) XXX, 1 6 Spherical Slepian functions and the polar gap in geodesy arxiv:math/6371v1 [math.st] 11 Mar 6 Frederik J. Simons 1 and F. A. Dahlen 1 Department of Earth Sciences, University
More informationFoundations of Matrix Analysis
1 Foundations of Matrix Analysis In this chapter we recall the basic elements of linear algebra which will be employed in the remainder of the text For most of the proofs as well as for the details, the
More informationLinear Algebra using Dirac Notation: Pt. 2
Linear Algebra using Dirac Notation: Pt. 2 PHYS 476Q - Southern Illinois University February 6, 2018 PHYS 476Q - Southern Illinois University Linear Algebra using Dirac Notation: Pt. 2 February 6, 2018
More informationDiagonalization of the Coupled-Mode System.
Diagonalization of the Coupled-Mode System. Marina Chugunova joint work with Dmitry Pelinovsky Department of Mathematics, McMaster University, Canada Collaborators: Mason A. Porter, California Institute
More informationNumerical Methods I Solving Nonlinear Equations
Numerical Methods I Solving Nonlinear Equations Aleksandar Donev Courant Institute, NYU 1 donev@courant.nyu.edu 1 MATH-GA 2011.003 / CSCI-GA 2945.003, Fall 2014 October 16th, 2014 A. Donev (Courant Institute)
More informationMATH 205C: STATIONARY PHASE LEMMA
MATH 205C: STATIONARY PHASE LEMMA For ω, consider an integral of the form I(ω) = e iωf(x) u(x) dx, where u Cc (R n ) complex valued, with support in a compact set K, and f C (R n ) real valued. Thus, I(ω)
More informationContents. Preface to the Third Edition (2007) Preface to the Second Edition (1992) Preface to the First Edition (1985) License and Legal Information
Contents Preface to the Third Edition (2007) Preface to the Second Edition (1992) Preface to the First Edition (1985) License and Legal Information xi xiv xvii xix 1 Preliminaries 1 1.0 Introduction.............................
More informationComputing Phase Noise Eigenfunctions Directly from Steady-State Jacobian Matrices
Computing Phase Noise Eigenfunctions Directly from Steady-State Jacobian Matrices Alper Demir David Long Jaijeet Roychowdhury Bell Laboratories Murray Hill New Jersey USA Abstract The main effort in oscillator
More informationLegendre s Equation. PHYS Southern Illinois University. October 18, 2016
Legendre s Equation PHYS 500 - Southern Illinois University October 18, 2016 PHYS 500 - Southern Illinois University Legendre s Equation October 18, 2016 1 / 11 Legendre s Equation Recall We are trying
More informationThe Prolate Spheroidal Phenomenon as a Consequence of Bispectrality
Centre de Recherches Mathématiques CRM Proceedings and Lecture Notes Volume 37, 2004 The Prolate Spheroidal Phenomenon as a Consequence of Bispectrality F. Alberto Grünbaum and Milen Yakimov Abstract.
More information(a) Determine the general solution for φ(ρ) near ρ = 0 for arbitary values E. (b) Show that the regular solution at ρ = 0 has the series expansion
Problem 1. Curious Wave Functions The eigenfunctions of a D7 brane in a curved geometry lead to the following eigenvalue equation of the Sturm Liouville type ρ ρ 3 ρ φ n (ρ) = E n w(ρ)φ n (ρ) w(ρ) = where
More informationEquations (3) and (6) form a complete solution as long as the set of functions fy n (x)g exist that satisfy the Properties One and Two given by equati
PHYS/GEOL/APS 661 Earth and Planetary Physics I Eigenfunction Expansions, Sturm-Liouville Problems, and Green's Functions 1. Eigenfunction Expansions In seismology in general, as in the simple oscillating
More informationarxiv: v1 [cs.lg] 26 Jul 2017
Updating Singular Value Decomposition for Rank One Matrix Perturbation Ratnik Gandhi, Amoli Rajgor School of Engineering & Applied Science, Ahmedabad University, Ahmedabad-380009, India arxiv:70708369v
More informationExamination paper for TMA4215 Numerical Mathematics
Department of Mathematical Sciences Examination paper for TMA425 Numerical Mathematics Academic contact during examination: Trond Kvamsdal Phone: 93058702 Examination date: 6th of December 207 Examination
More informationSome recent results on controllability of coupled parabolic systems: Towards a Kalman condition
Some recent results on controllability of coupled parabolic systems: Towards a Kalman condition F. Ammar Khodja Clermont-Ferrand, June 2011 GOAL: 1 Show the important differences between scalar and non
More informationSeparation of Variables in Linear PDE: One-Dimensional Problems
Separation of Variables in Linear PDE: One-Dimensional Problems Now we apply the theory of Hilbert spaces to linear differential equations with partial derivatives (PDE). We start with a particular example,
More informationMidterm for Introduction to Numerical Analysis I, AMSC/CMSC 466, on 10/29/2015
Midterm for Introduction to Numerical Analysis I, AMSC/CMSC 466, on 10/29/2015 The test lasts 1 hour and 15 minutes. No documents are allowed. The use of a calculator, cell phone or other equivalent electronic
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 informationAngular momentum. Quantum mechanics. Orbital angular momentum
Angular momentum 1 Orbital angular momentum Consider a particle described by the Cartesian coordinates (x, y, z r and their conjugate momenta (p x, p y, p z p. The classical definition of the orbital angular
More informationExercise Sheet 1.
Exercise Sheet 1 You can download my lecture and exercise sheets at the address http://sami.hust.edu.vn/giang-vien/?name=huynt 1) Let A, B be sets. What does the statement "A is not a subset of B " mean?
More informationChapter 4: Interpolation and Approximation. October 28, 2005
Chapter 4: Interpolation and Approximation October 28, 2005 Outline 1 2.4 Linear Interpolation 2 4.1 Lagrange Interpolation 3 4.2 Newton Interpolation and Divided Differences 4 4.3 Interpolation Error
More informationAlgebra C Numerical Linear Algebra Sample Exam Problems
Algebra C Numerical Linear Algebra Sample Exam Problems Notation. Denote by V a finite-dimensional Hilbert space with inner product (, ) and corresponding norm. The abbreviation SPD is used for symmetric
More informationOrthogonal Polynomials and Gaussian Quadrature
Orthogonal Polynomials and Gaussian Quadrature 1. Orthogonal polynomials Given a bounded, nonnegative, nondecreasing function w(x) on an interval, I of the real line, we consider the Hilbert space L 2
More informationFun With Carbon Monoxide. p. 1/2
Fun With Carbon Monoxide p. 1/2 p. 1/2 Fun With Carbon Monoxide E = 0.25 ± 0.05 ev Electron beam results p. 1/2 Fun With Carbon Monoxide E = 0.25 ± 0.05 ev Electron beam results C V (J/K-mole) 35 30 25
More informationMATH 590: Meshfree Methods
MATH 590: Meshfree Methods The Connection to Green s Kernels Greg Fasshauer Department of Applied Mathematics Illinois Institute of Technology Fall 2014 fasshauer@iit.edu MATH 590 1 Outline 1 Introduction
More informationLinear Algebra Practice Problems
Linear Algebra Practice Problems Math 24 Calculus III Summer 25, Session II. Determine whether the given set is a vector space. If not, give at least one axiom that is not satisfied. Unless otherwise stated,
More informationNumerical Methods I Solving Square Linear Systems: GEM and LU factorization
Numerical Methods I Solving Square Linear Systems: GEM and LU factorization Aleksandar Donev Courant Institute, NYU 1 donev@courant.nyu.edu 1 MATH-GA 2011.003 / CSCI-GA 2945.003, Fall 2014 September 18th,
More informationCayley-Hamilton Theorem
Cayley-Hamilton Theorem Massoud Malek In all that follows, the n n identity matrix is denoted by I n, the n n zero matrix by Z n, and the zero vector by θ n Let A be an n n matrix Although det (λ I n A
More informationComputation of the scattering amplitude in the spheroidal coordinates
Computation of the scattering amplitude in the spheroidal coordinates Takuya MINE Kyoto Institute of Technology 12 October 2015 Lab Seminar at Kochi University of Technology Takuya MINE (KIT) Spheroidal
More informationMA3025 Course Prerequisites
MA3025 Course Prerequisites MA 3025 (4-1) MA3025 (4-1) Logic and Discrete Mathematics: Provides a rigorous foundation in logic and elementary discrete mathematics. Topics from logic include modeling English
More information