Quadrature integration for orthogonal wavelet systems
|
|
- George Russell
- 5 years ago
- Views:
Transcription
1 Quadrature integration for orthogonal wavelet systems Bruce R. Johnson, a Jason P. Modisette, b Peter J. Nordlander, b and James L. Kinsey a a Department of Chemistry and Rice Quantum Institute Rice University Houston, TX b Department of Physics and Rice Quantum Institute Rice University Houston, TX
2 Abstract Wavelet systems can be used as bases in quantum mechanical applications where localization and scale are both important. General quadrature formulae are developed for accurate evaluation of integrals involving compact support wavelet families, and their use is demonstrated in examples of spectral analysis and integrals over anharmonic potentials. In contrast to usual expectations for these uniformly-spaced basis functions, it is shown that nonuniform spacings of sample points are readily allowed. Adaptive wavelet quadrature schemes are also presented for the purpose of meeting specific accuracy criteria without excessive oversampling. 2
3 I. Introduction There is now an extensive wavelet literature within a variety of applications, e.g., digital signal processing, image processing, fingerprint classification, remote sensing, target recognition, etc. 1-7 There has been some exploration of the use of wavelets in chemical applications, 8-20 but of only limited extent to date. Nevertheless, wavelet methods promise a distinct advantage over Fourier methods for problems requiring different levels of resolution in different locations, e.g., spectral analysis and compression, molecular calculations of electronic or nuclear motion, etc. Regarded as basis functions for calculation or analysis of quantum wave functions, wavelet systems offer certain unusual properties. In fact, only with the introduction of the orthonormal compact support wavelets by Daubechies 1,3 did it become clear that there exist basis sets that are capable of being simultaneously localized, multiscale, multicenter, and orthonormal. The purpose of this paper is to provide efficient and accurate quadrature methods for required integrals over compact support wavelet families and their generalizations. Such functions do not have simple closed forms, but are instead defined by recursion between different scales. The original Daubechies bases are generated by two related functions, the scaling function φ(x) and the wavelet ψ(x). Both are nonzero only for 0 < x < L 1, where L is an even integer. Figure 1 shows these functions for L = 8. From these functions one obtains others that are ust copies widened by factors 2 and translated by steps k 2, where and k are integers, φ k ( x) = 2 /2 φ ( 2 x k ), (1) ψ k ( x) = 2 /2 ψ ( 2 x k ), (2) 3
4 with φ 0 0 = φ, ψ 0 0 = ψ. By construction, the set of scaling functions at a particular resolution level, { φ k }, form a basis which is orthogonal though not complete. A more 1 complete set { φ k } is also orthogonal but is twice as dense. The difference in detail between these two bases is spanned by the set of wavelets at level, so that 1 { φ k } = { φ k } { ψ k }. This is expressed in terms of the recursion relations, φ k L 1 ( x) = h k φ 1 2 k+ k ( x ), (3) k = 0 ψ k L 1 ( x) = g k φ 1 2 k + k ( x ), (4) k =0 where the constants are related by g k = ( 1) k +1 h L 1 k. 2 Ultimately, every quantity that must be calculated boils down to some function of the h k, so that the latter take on a central role in calculations. An orthonormal multiresolution basis can be comprised of J scaling functions { φ k } on a coarsest level = J and wavelet functions { ψ k } from this coarsest level down to a finest scale = 0 which provide the various additional octaves of detail. The approximate expansion of a function f in this basis then takes the form ( ) ψ k f x J J J f ψ k ( x) + φ k f φ k ( x). (5) = 0 k Detail functions { ψ k } at each scale which are unimportant (have small expansion k coefficients ψ k f ) may simply be dropped, leaving a finite basis with customized resolution. This immediately generalizes to multidimensional Cartesian cases by using products of such basis functions. For quantum applications, the needs arise to calculate the proection integrals φ k f, ψ k f for the vector basis and the matrix elements φ k f φ k, φ k f ψ k, 4
5 ψ k f ψ k, where f may be a general function of x. Use of the recursion relations in Eqs. (3) and (4) allows all of these integrals to be expressed in terms of integrals with a scaling function (or functions) on a single finest scale, but accurate methods are still required for calculation of the remaining set of integrals. One method mentioned by Daubechies 3 (p. 166) for calculation of the vector proection integrals uses Fourier convolution. While this method is not widely used, Wei and Chou 17 have adapted it to their density functional calculations. Alternatively, a quadrature method with uniformlyspaced sample points was derived for the vector elements by Sweldens and Piessens. 21,22 The latter method was used by Modisette, et al., 15 for wavelet expansion of a semisingular double-well potential; by expanding the wave functions in the same basis, matrix elements could then be evaluated in terms of exactly-calculable "connection coefficients," or integrals of products of three functions from the wavelet basis. 7,23,24 The portability of this approach is complicated by the large number of connection coefficients generally required, however. A recent adaptive wavelet collocation method developed for O(N) solution of PDE's 25 directly addresses calculation of matrix elements, though a linear system of equations depending on the sample grid must generally be solved. In the present paper, parallel treatments are given of polynomial-based quadrature formulae for direct calculation of φ k f and φ k f φ k. An unusual aspect is that there is no requirement that sample points lie on a uniformly spaced grid. For a given uniform or irregular grid, the quadrature weights are explicitly expressed in terms of Lagrange interpolating polynomials and polynomial moments of the scaling functions. The latter are independent of the grid choice and need only be calculated once. Furthermore, use of the scaling function recursion in Eq. (3) leads naturally to adaptive evaluation whereby a 5
6 chosen accuracy goal can be reached without excessive over-sampling. Straightforward generalizations of these developments also hold for biorthogonal compact support wavelet systems. 26 II. Proection Integral Quadrature For a given and k, a change of variables transforms the general vector integral for scaling functions into a standardized form, φ k f = 2 /2 dx φ ( 2 x k ) f x ( ) = 2 / 2 dx φ( x ) f [ 2 ( x + k )]. (6) Thus it is only necessary to consider the = k = 0 case, for which the integral is approximated by quadrature, ( ) dx φ( x ) f ( x ) ω q f x q r. (7) q =1 The only restriction imposed on the sample points is that the x q are assumed to be ordered, x 1 < x 2 < x r. One can specify theω q in terms of the x q by requiring that Eq. (6) be exact for f equal to a polynomial of order r 1 or less: r p ω q x q = m p dx φ( x ) x p, 0 p r 1. (8) q= 1 The general solution to these equations (see Appendix A) is ω q = r 1 p= 0 m p p! ( p ) L r,q ( 0 ), (9) using the Lagrange interpolating polynomials L r, q ( x ) = r x x q (10) q q x q x q 6
7 ( p and their derivatives L ) r, q at x = 0. The Lagrange polynomials are the coefficients for interpolation of smooth functions, 27 r L r,q q= 1 ( ) ( x ) f x q f ( x), (11) with equality holding for f(x) a polynomial in x of order less than r. Their connection with numerical quadrature is well known, 28,29 although the simple result Eq. (9) (which is not really restricted to wavelet systems) appears to be new. Lagrange polynomials are also familiar in wavelet contexts from their use in interpolatory refinement schemes, e.g., the symmetric iterative interpolation process of Deslauriers and Dubuq. 30 The latter scheme uses a "fundamental solution" that was subsequently shown 31 to coincide with the autocorrelation function of the compact support scaling functions (see Appendix B). The Lagrange polynomials 29,32 and their derivatives are easily evaluated once the quadrature nodes are specified. As for the moments m p, their calculation may be accomplished by the simple recursive equation derived by Gopinath and Burrus, 33 m p = 1 2 p ( ) 1 p 1 p m p p µ p p, (12) p = 0 where the µ p are discrete moments of the scaling coefficients of Eq. (3), µ p = 1 2 L 1 h k k p. (13) k =0 The leading moment m 0 = 1 as a matter of normalization, providing the starting point for the recursion. (Not all of the resulting moments m p are independent of each other, as is discussed in Appendix B.) The moments are thus determined completely by the scaling coefficients h k, whereas the quadrature coefficients are determined in terms of both the h k 7
8 and the nodes x q. In calculations for a given compact support basis, the moments need only be calculated once. Sweldens and Piessens 21 considered the important case of uniformly-spaced x q. Such an arrangement allows the sample points for scaling functions of adacent k to overlap, requires only a single set of quadrature coefficients for all k, minimizes the number of function evaluations needed, and yields a higher order error of O( x r ), where x is the spacing. This is a Newton-Cotes type of quadrature 28 for which is allowed a single shift of the entire set of sample points with respect to the origin; this degree of freedom can be used to make the quadrature exact for polynomials up to order r instead of r 1. One may achieve even greater accuracy by defining a Gauss quadrature for which all r of the x q are chosen (non-uniformly) so as to exactly include polynomials up to order 2r 1; in this case there is no overlap between the sample points for adacent scaling functions. In the present work, the disposition of the x q is left open for purposes of adaptability. For compact support wavelets, it is known 21,33 that m 2 = m 1 2, although higher moments are not simple powers of m 1 (Appendix B). Introducing the differences m p = m p m 1 p, (14) Eq. (8) can be rearranged, ω q = r 1 p= 0 m 1 p p! ( p) L r, q ( 0)+ r 1 p= 3 m p p! ( p ) L r,q ( 0) ( ) + = L r,q m 1 r 1 p= 3 m p p! ( p ) L r, q ( 0), (15) 8
9 where the first series has been recognized as the Taylor expansion of the Lagrange interpolation polynomial evaluated at m 1. To the extent that the second series can be ignored, Eqs. (7) and (11) show that ( ) dx φ( x ) f ( x ) L r,q ( m 1 ) f x q r f m 1 q= 1 ( ) (16) This directly reflects the fact determined by Gopinath and Burrus that a one-point quadrature using m 1 as the sample point has O( x 3 ) error. Higher order error generally requires a multi-point quadrature. An exception is found in the Coifman wavelet family, 34,35 where m 1 = m 2 = = m M -1 =0 for chosen integer M, leading to a one-point quadrature error of O( x M ). For many applications, it is worthwhile to use a more flexible framework than orthogonal wavelet families. In biorthogonal families, 26 the functions φ and ψ are orthogonal to dual functions φ and ψ rather than to themselves. The scaled and translated versions of these functions, φ k Eqs. (3) and (4), and ψ k, satisfy their own scaling relations, cf., φ k ( x) = h 1 k k ( x), (17) k φ 2k+ ψ k ( x) = g 1 k k ( x), (18) k φ 2k+ The extra freedom within biorthogonal families allows, for instance, φ and ψ to possess definite symmetry, an impossibility for the orthogonal compact support families. 1,3 It also provides the framework for the Lifting scheme introduced by Sweldens 36,37 for 9
10 purposes of adaptable multiresolution analyses. The general expansion of a function in a biorthogonal basis may be expressed as ( ) f x J ψ k f ψ k ( x) + J φ k = 0 k k rather than Eq. (5). Similarly the analogs of Eqs. (7) and (8) are r dx φ ( x ) ( ) f x r ( ) J f φ k ( x) (19) ω q f x q, (20) q =1 ω x p = m q q p dx φ ( x ) x p, 0 p r 1. (21) q =1 The structure of the equations is otherwise unchanged so that the quadrature coefficent expressions immediately adapt to the biorthogonal wavelet case. III. Wavelet Transform of UV Ozone Absorption Spectrum One of the potential uses of wavelet analysis is to molecular spectra with complex features corresponding to action on different scales. The ozone molecule provides a particular example in its Hartley UV absorption band (see Johnson, et al., and references therein), as shown in Fig. 2. Photoexcitation to a steeply repulsive dissociative region of the potential energy surface for the D 1 B 2 state leads to a single intense band which is continuous across several thousand cm 1. In addition, smaller features appear across much of the band as minor irregular structure with average spacing 250 cm 1. Under Fourier transformation, the time domain autocorrelation function exhibits completely resolved features at shorter (< 10 fs) and longer ( fs) times corresponding, respectively, to direct and indirect dissociation contributions. Dynamical analysis of the indirect contribution, corresponding to Franck-Condon-region exit and 10
11 return, has progressed in recent years with the aid of both classical 40 and quantum propagations. More definitive information is currently being obtained by Le, et al., 44, in the form of detailed Raman excitation profiles for the different accessible final vibrational states; these exhibit structure that reflects upon the same excited-state dynamics, but in a vibration-specific manner. One limitation of the Fourier transform (FT) is that it completely defocuses information about the old variable ν. The tremendous breadth of the total band implies that there are many different energetic regimes mixed together in the autocorrelation function, and that a oint time-frequency analysis with partial resolution in both variables is to be preferred. The first such application to ozone used a sliding window on the spectral data prior to the FT operation, 39 thereby allowing the following of time features between partially-localized regions of frequency. Gabor 45 and Vibrogram/Husimi 46 transform techniques have also proven to be valuable alternatives for time-frequency analysis of classical and quantum dynamics in other systems, though they have not been applied to ozone. In the present work, an alternative separation of the distinct dynamical regimes is obtained through wavelet multiresolution analysis of the absorption spectrum. Using one of the orthonormal Daubechies wavelet families (L = 8), it is possible to calculate proection integrals of translated and scaled scaling functions on a level = 0, dν 1 ν φ ν ν 0 ν k σ( ν) ( ) = ν dx φ x r σ ( ν 0 + k ν + x ν ) ( ) ν ω q σ ν 0 + k ν + x q ν. (22) q=1 11
12 An absolutely typical initial complication arises from the fact that the absorption crosssection σ was measured 47 in uniform wavelength steps of 0.1 Å, i.e., the frequency variable ν has non-equidistant spacings. One may perform interpolation of different types in order to obtain uniform steps in ν, but this is unnecessary if the quadrature described above is used. A total of r irregularly spaced experimental points (with average spacing close to ν) were selected from within the support of each φ k 0 and the associated proections then calculated by the rules described above. With the average spacing chosen to be 3 cm 1, comparison of integrals for r = 7 and r = 8 generally showed convergence to 5 significant figures. This reflects, of course, a combination of both the inherent quadrature errors and any point-to-point experimental uncertainties. Using the recursions in Eqs. (3) and (4), the proections onto an orthonormal basis J of wavelets { ψ k }, J, and scaling functions { φ k } on level J may then be calculated in the normal fashion. In this way, the spectrum is decomposed simultaneously into different scales and positions. The smoothed version of the absorption cross-section shown in Fig. 2 was obtained by summing over all k the proections onto the φ k J for J = 8. This partial reconstruction process requires an algorithm that accurately calculates the values of the scaling functions at particular points (see Burrus, et al. 6 ). Figure 3 shows the corresponding summed ψ k proections, representing the differences in detail between different scales. It is seen that the spectral structure is separated from the smooth component and is concentrated primarily in three octaves of scale change spanned by the wavelets for = 5 7. One obtains localized frequency information even from these summed quantities, e.g., it is clear from = 5 that higher-frequency components of the structure are more prominent on the low-frequency side of band center. 12
13 The corresponding wavelet coefficients shown in Fig. 4 of course contain much more detailed local frequency information. All of the information so obtained is preserved under Fourier transformation of the wavelets and scaling functions to the time domain, so that one can determine how the longer-time features in the autocorrelation function change with both scale and general frequency region, subect only to the constraints of the Uncertainty Principle. (Visual analysis, e.g., for comparison with the windowed FT analysis of Johnson and Kinsey, 39 would be simpler using the Continuous Wavelet Transform 3 although the latter is redundant and less efficient to implement.) A more important advantage of this type of decomposition should be realized in the analysis of the Raman experiments in ozone. Transformation of a Raman excitation profile to the time domain is not a simple linear operation since the phase is not measured and must be recovered using Hilbert transform techniques The wavelet decomposition should prove useful here for separation of the early and later time domain features; from these results the shapes of the excited state potential and transition dipole in the Franck-Condon geometry region and beyond can be determined 53,54 (see also Shapiro 55 ). It merits emphasis that the irregular sampling is accommodated easily in this application. This type of irregularity stands in contrast to that in the Lifting scheme, 36,37 where the biorthogonal basis itself is tailored. In the present case, the underlying wavelet basis is regularly-spaced as usual, but the sample points are not tied to that arrangement. The numerical quadrature performed for the expansion coefficients on the finest scale does not change the O( N) nature of the wavelet recursion, but does increase the work needed for initialization by a factor proportional to the quadrature order r. In the language of filter banks used in engineering contexts, this corresponds to a "prefiltering" 13
14 operation on the data, 5,56,57 and the quadrature scheme used here represents a polynomial prefilter for irregularly-spaced samples. IV. Morse Potential Expansion and Adaptive Quadrature In wavelet applications, easy compression is one of the usual goals. If a particular wavelet coefficient is calculated to be small, that element of the basis may be eliminated. On the other hand, efficiency and/or storage in numerical applications may require that these coefficients should not even be calculated if they are known to be small, for instance, if overlapping elements from coarser scales have already been excluded. This argues for a refinement approach 3,30 where finer scales are used only as needed by the problem under investigation. In this section, the results of Section II are embedded in an adaptive wavelet quadrature. A prototypical anharmonic molecular potential function is given by the Morse potential, ( ) = D e exp 2a( x x e ) V x [ ] 2D e exp[ a( x x e )], (23) characterized by a steeply repulsive inner wall, a more slowly-varying neighborhood of the minimum, and a flat outer region. Parameters appropriate to molecular hydrogen are used, a = Å 1, x e = Å, D e = ev. Proections are calculated of this function on the scaling functions from the L = 10 Daubechies system, φ k J ( x) = 2 J λ ( ) 1/2 φ x /2 J λ k ( ), (24) with spacing 2 J λ = 0.4 Å (λ and J are not independent parameters here). 14
15 For the sake of simplicity (but not necessity), the uniformly-spaced quadrature associated with the basis functions on each level is used, x q = ( k + q 1) 2 λ. The same quadrature weights are therefore used for each functionφ k, and r 1 of the points for neighboring values of k overlap. In Fig. 5 the absolute values of the differences between integrals for orders r and r + 1 are displayed; it is very clear that significant lack of convergence occurs for low r when only the coarsest quadrature is used. The highly inaccurate case r = 1 is the usual basis of the Fast Wavelet Transform. Given the O( x r )error, proceeding to finer scales unfortunately provides only slow linear convergence with respect to spacing (or cubic in the special case that the sampling is moved to the position of the first moment). Errors are exponentially reduced for higher r in the range covered, r 10. (Even higher orders r are allowed but, as is generally the case with Newton-Cotes quadratures, can eventually lead to divergence.) The asymptotic rates of convergence (slopes on the right-hand side of the semi-log plot) all conform. The absolute errors, however, differ systematically for different k. As is evident from the top of Fig. 6, the source of the differences is the strongly varying behavior of the anharmonic Morse potential over the region spanned by the different φ k 's. Greater accuracy for level J can be obtained by evaluation of the scaling function integrals on levels J 1, J 2,, until the desired level of accuracy is achieved, followed by use of Eqs. (3) and (4) to recur back to level J. The best procedure is clearly to use a value of r large enough to obtain satisfactory convergence rates and thus to avoid needlessly high numbers of recursions, but coupled with selectivity over which regions are subected to refinement. Convergence may be measured either by (i) comparison of the numerical results using quadrature at level with those of quadrature at 15
16 level 1 followed by recursion back to, or (ii) direct quadrature at level using orders r and r + 1 as done for the highest level in Fig. 5. The latter (simpler) choice was found to be adequate in the calculations here, though this need not always be the case. In the bottom section of Fig. 6 are shown the results for a demanding coarsestlevel convergence tolerance of 10 8 ev. For each level of refinement, the tolerance is correspondingly reduced by a factor of 2. The markers in each row correspond to the left-hand sides of included scaling functions; these functions are shown explicitly for the coarsest level = 4. On the finest level required, = 0, the adaptive procedure has reduced the 199 potentially-required integrals to only 20. One may obtain greater compression for a smaller choice of L, although this comes at the price of reduced smoothness in the basis functions. Efficiency in the calculations can be gained by utilizing the fact that only approximately half of the sample points required on each finer level are new. Furthermore, each φ k for which the integral is recalculated using refinement requires the integrals for φ 1 2k, φ 2k ,, φ 2k +L 1, but there is a multiplexing advantage sinceφ k ±1 require L 2 of the same integrals, etc. Finally, when performing calculations with a basis consisting of functions from different levels of resolution, many of the intermediate integrals calculated will be of use for basis function integrals at more than one scale. Thus, the intermediate results can be stored and re-used to speed execution time. If storage becomes an issue (e.g., in multidimensional wavelet calculations), however, it may be more desirable to allow for a certain amount of recalculation of integrals to trade speed for storage savings. 16
17 V. Matrix Element Quadrature The considerations above have only considered integrals linear in the vector basis. Bilinear integrals are standard in quantum mechanics, where one typically evaluates matrix elements of kinetic and potential energy operators between two basis functions. While kinetic energy operators in a compact support basis may be calculated by straightforward methods, 23,52 the same is not true for matrix elements of a potential which is some general function of the coordinate (or, in multiple dimensions, coordinates). Nevertheless, the quadrature methods described above translate straightforwardly to the case of such bilinear integrals. As before, the problem is first simplified by virtue of the recursions in Eqs. (3) and (4). For a particular finite multiresolution basis of wavelets and scaling functions, all necessary matrix elements can be reduced to integrals between scaling functions on the same scale. It is therefore sufficient to consider the = 0 case, 0 φ k f φ k 0 ( ) f x = dxφ x k ( )φ( x k ). (25) This matrix is automatically banded. Since the scaling functions vanish for arguments outside (0, L 1), functions of index k and k' spatially overlap only for k k L 2. Due to symmetry, there are L 1 distinct cases, exemplified by the specific choices k' = 0, k = 0, 1, 2,, L 2. A quadrature strategy is adopted which allows for distinct sets of quadrature points x q,k and weights ω q,k for each of these L 1 values of k, dx φ( x ) f ( x )φ( x k ) ω q,k f x q, k r k. (26) q= 1 ( ) 17
18 The x q,k may be chosen to overlap for different values of k, reducing the number of samples required. Insisting that the quadrature be exact for powers 0 through r k 1 leads to r k ω q, k x p q,k = m p,k dx φ x q =1 ( ) x p φ( x k), 0 p r k 1. (27) While the moments are different from before, the structure of the equations is identical. Therefore the weights are given by ω q, k = r k 1 m p,k L ( p ) p! rk,q,k i =0 ( 0). (28) The particular moments in Eq. (27) may be obtained in the manner used by Beylkin 52 to obtain matrix representations of operators. This method uses change of scales for both scaling functions in Eq. (27) and change of integration variables to obtain linear equations that the moments must satisfy, p m p,k = 2 p p L 1 m p p,2 k + k a p p, k. (29) p = 0 k =1 L The a p,k are the easily calculated discrete sums a p, k = h n +k h n n p, (30) n where the scaling function coefficients h n are zero unless 0 n L 1. Also, from Eq. (27) the moments for k < 0 can be shown to be related to those for k > 0 by m p, k = p p =0 p p ( k) p p m p,k. (31) Thus Eqs. (29) need only be solved for the moments with non-negative k (using standard methods), and this is required only once. 18
19 The polynomial-based quadrature above can be adapted straightforwardly to matrix elements between the dual sets of functions in a biorthogonal basis. The latter type of integral can arise, for instance, in solution of differential equations using the biorthogonal systems in the Wavelet Galerkin method In quantum contexts, biorthogonal matrices may have application to problems in which non-hermitian Hamiltonian-type operators arise, e.g., in the use of complex rotation methods for L 2 basis calculations of resonance eigenvalues. 62,63 One loses the advantages of symmetry of the matrices in these cases, but gains in flexibility. Following the Morse potential example of the last section, an adaptive quadrature may be implemented to evaluate in tandem the potential matrix elements for the L = 10 scaling function basis. The same uniformly-spaced quadrature is used, but this time the starting point of the grid for each matrix element must be chosen for two scaling functions which partially overlap. A reasonable choice is to start at the left-hand side of their region of intersection, i.e., at 2 λ max(k', k ). The multiplexing advantage stemming from overlapping grids still applies despite the fact that a matrix is evaluated in this case. The number of samples required at the first stage is only linearly proportional to the distance spanned by the scaling function basis at level. Convergence may again be monitored by comparing results for adacent orders r k and r k + 1. All matrix elements with the same starting samples, i.e., the same max(k', k ), are regarded ointly; if any of them fail to converge, then the entire subset is earmarked for more accurate evaluation by recursion of the relevant integrals from the next finer level. This starts a recursive refinement effectively focusing on the local behavior of the potential sample points. The procedure is continued to finer scales, tightening the 19
20 convergence threshold by a factor of 2 each time, until the desired convergence is reached and the results may be backtracked to the original scale. The alternative convergence criterion mentioned in the vector case applies to the matrix case as well: one may monitor convergence by agreement of quadrature at level and of quadrature at level 1 with recursion to level. This would require one extra scale as well as a greater amount of computation at intermediate steps. While the matrix integrals are thus capable of being calculated in close analogy to the vector integrals, there is one chief difference. The recursion between scales must be implemented for both the bra and ket scaling functions at the same time, thereby increasing the computational effort. In this regard, it is in principle possible to streamline things by combining the two separate matrix operations into a single compound operation. For the level basis, 0 k 4, the required matrix elements were sought with a convergence threshold of 10 5 ev. Quadrature orders 6 and 7 were used for checking of convergence for all matrix elements. Table 2 shows the rapid convergence as a function of scale for all of the matrix elements needed to evaluate the top level integrals. The number of distinct function evaluations is kept minimal by the current procedure. The required matrix elements φ k V φ k are indicated explicitly in Fig. 7, with the starting points of the two scaling functions mapped onto separate axes. Thus one can see immediately the spatial distribution of the matrix elements required and, in particular, the adaptive increase in density for the region near the steep inner wall of the potential. 20
21 VI. Conclusion Integration by polynomial-based quadrature for orthonormal (or bi-orthonormal) compact support wavelet families has been investigated for the purpose of use in chemical and physical applications. In the above, a quadrature formula has been developed in terms of easily calculated derivatives of Lagrange interpolation polynomials for the chosen coordinate grid and particular moments which are grid independent. Easy application to irregular grids was demonstrated in the wavelet decomposition of the Hartley absorption spectrum of ozone. Adaptive quadrature was developed for integrals of products of reasonably smooth functions and one or two scaling functions, and applied to integrals of the anharmonic Morse potential of the types that will arise in numerical quantum mechanical calculations using wavelet bases. Acknowledgments This work was supported by the National Science Foundation and the Robert A. Welch Foundation. The authors also wish to thank R. O. Wells, C. S. Burrus, J. Tian and W. Sweldens for conversations. 21
22 Appendix A Derivation of Eq. (9) for the quadrature weights proceeds directly from the interpolation formula in Eq. (11). The equality sign holds for the latter if f ( x) is a power of x less than r, r n L r,q ( x ) x q = x n. (A1) q=1 Differentiating this with respect to x p times and setting x equal to zero yields r ( L p ) r,q q= 1 n ( 0) x q = p!δ p, n, (A2) or, multiplying both sides by m p / p!, m p r p! L ( p ) r,q ( 0 ) x n q = m p δ p,n. (A3) q =1 If this is summed over p up to r 1, the Kronecker delta is unity for p = n and zero otherwise. Thus, r r 1 m p p! L ( p ) r,q ( 0) p = 0 x n q = m n. (A4) q=1 Since this is true for all n < r, the quantity in brackets is the solution of Eq. (8). Appendix B Nonlinear relations exist between various moments of the scaling function in orthogonal compact support wavelet systems. This is most easily seen by considering the autocorrelation function (the fundamental function of the Deslauriers-Dubuq scheme 30 ) ( ) = dx Φ y φ( x )φ ( x + y ). (B1) 22
23 The moments of Φ( y ) may be determined after a little algebra in terms of the moments m p of the scaling function defined in Eq. (8), 64 M n = dy y n Φ y n n ( ) = ( 1) n p p =0 m p n p m p (B2) On the other hand, Beylkin 52 and Tian and Wells 35 have shown that M n = 0 for n = 0, 1,, L 1. For odd n this yields no new information; for even n, however, one finds m 2 = m 2 1, L 4, (B3) m 4 = 4 m 3 m 1 3m 4 1, L 6, (B4) m 6 = 6m 5 m 1 +10m m 3 m m 6 1, L 8, (B5) m 8 = 8m 7 m 1 +56m 5 m 3 168m 5 m m 3 m m 2 3 m m 8 1, L 10, (B6) etc. The first of these relations was derived by Gopinath and Burrus 33 for the Daubechies scaling functions and later shown to be more general by Sweldens and Piessens. 21 To the knowledge of the authors, the higher-order relations do not appear to have been noted before. They are valid for all of the compact support wavelet families since the latter all share the same autocorrelation function. 23
24 Bibliography 1 I. Daubechies, Comm. Pure Appl. Math. 41, 909 (1988). 2 S. G. Mallat, IEEE Trans. Pattern Anal. Machine Intell. 11, 674 (1989). 3 I. Daubechies, Ten Lectures on Wavelets (SIAM Publications, Philadelphia, 1992). 4 Wavelets: Mathematics and Applications, Vol., edited by J. J. Benedetto and M. W. Frazier (CRC Press, Boca Raton, 1994). 5 G. Strang and T. Nguyen, Wavelets and Filter Banks (Wellesley-Cambridge Press, Wellesley, 1996). 6 C. S. Burrus, R. A. Gopinath, and H. Guo, Introduction to Wavelets and Wavelet Transforms (Prentice-Hall, 1996). 7 H. L. Resnikoff and R. O. Wells, Jr., Wavelet Analysis: The Scalable Structure of Information (Springer Verlag, 1998). 8 T. Paul, in Wavelets: Time-Frequency Methods and Phase Space, edited by J. M. Combes, A. Grossman, and P. Tchamitchian (Springer Verlag, Berlin, 1989), pp D. Permann, Phys. Rev. Lett. 69, 2607 (1992). 10 Z. Li, A. Borrmann, and C. C. Martens, Chem. Phys. Lett. 214, 362 (1993). 11 K. Cho, T. A. Arias, J. D. Joannopolous, and P. K. Lam, Phys. Rev. Lett. 71, 1808 (1993). 12 P. Fischer and M. Defranceschi, Int. J. Quant. Chem. 45, 619 (1993). 13 D. Permann and I. Hamilton, J. Chem. Phys. 100, 379 (1994). 14 L. Gagnon and J. M. Lina, J. Phys. A 27, 8207 (1994). 24
25 15 J. P. Modisette, P. Nordlander, J. L. Kinsey, and B. R. Johnson, Chem. Phys. Lett. 250, 485 (1996). 16 A. Askar, A. E. Cetin, and H. Rabitz, J. Phys. Chem. 100, (1996). 17 S. Wei and M. Y. Chou, Phys. Rev. Lett. 76, 2650 (1996). 18 P. Fischer and M. Defranceschi, SIAM J. Numer. Anal. 35, 1 (1996). 19 S. Goedecker and O. V. Ivanov, Sol. State Comm. 105, 665 (1998). 20 D. J. Kouri, (to be published). 21 W. Sweldens and R. Piessens, SIAM J. Numer. Anal. 31, 1240 (1994). 22 W. Sweldens and R. Piessens, Numerische Mathematik 68, 377 (1994). 23 A. Latto, H. L. Resnikoff, and E. Tenenbaum, Compactly Supported Wavelets, AD91078, Aware, Inc. (1991). 24 W. Dahmen and C. A. Micchelli, SIAM J. Numer. Anal. 30, 507 (1993). 25 O. V. Vasilyev and S. Paolucci, J. Comp. Phys. 125, 498 (1996). 26 A. Cohen, I. Daubechies, and J. C. Feauveau, Comm. Pure Appl. Math. 55, 485 (1992). 27 M. Abramowitz and I. A. Stegun, Handbook of Mathematical Functions with Formulas, Graphs, and Mathematical Tables (U. S. Government Printing Office, Washington, D.C., 1972). 28 A. Ralston and P. Rabinowitz, A First Course in Numerical Analysis, Second ed. (McGraw-Hill, New York, 1978). 29 H. R. Schwarz and J. Waldvogel, Numerical Analysis: A Comprehensive Introduction (Wiley, Chichester, 1989). 25
26 30 G. Deslauriers and S. Dubuc, Constr. Approx. 5, 49 (1989). 31 M. J. Shensa, IEEE Trans. Signal Processing 40, 2464 (1992). 32 H. Rutishauser, Lectures on Numerical Mathematics (Birkhaüser, Boston, 1990). 33 R. A. Gopinath and C. S. Burris, in Proceedings of the ISCAS-92 (San Diego, 1992). 34 G. Beylkin, R. Coifman, and V. Rokhlin, Comm. Pure and Appl. Math 44, 141 (1991). 35 J. Tian and R. O. Wells, Jr., Vanishing Moments and Wavelet Approximation, Computer Mathematics Laboratory Report CML9501, Rice University (1995). 36 W. Sweldens, in Wavelet Applications in Signal and Image Processing III, Proc. SPIE 2569, edited by A. F. Laine and M. Unser (1995), pp W. Sweldens, SIAM J. Math. Anal. 29, 511 (1998). 38 B. R. Johnson and J. L. Kinsey, Phys. Rev. Lett. 62, 1607 (1989). 39 B. R. Johnson and J. L. Kinsey, J. Chem. Phys. 91, 7638 (1989). 40 B. R. Johnson, B.-Y. Chang, C.-W. Hsiao, L. Le, and J. L. Kinsey, J. Chem. Phys. 108, 7670 (1998). 41 C. Leforestier, F. LeQuéré, K. Yamashita, and K. Morokuma, J. Chem. Phys. 101, 3806 (1994). 42 N. Balakrishnan and G. D. Billing, J. Chem. Phys. 101, 2968 (1994). 43 C. Leforestier, (to be published). 44 L. D. Le, C.-W. Hsiao, E. S. Lotfi, C. Kittrell, B. R. Johnson, and J. L. Kinsey, (in progress). 26
27 45 D. C. Rouben and G. S. Ezra, J. Chem. Phys. 103, 1375 (1995). 46 K. Hirai, E. J. Heller, and P. Gaspard, J. Chem. Phys. 103, 5970 (1995). 47 J. Brion, A. Chakir, D. Daumont, J. Malicet, and C. Parisse, Chem. Phys. Lett. 213, 610 (1993). 48 F. Remacle, R. D. Levine, and J. L. Kinsey, Chem. Phys. Lett. 205, 267 (1993). 49 F. Remacle and R. D. Levine, J. Chem. Phys. 99, 4908 (1993). 50 F. Remacle and R. D. Levine, J. Phys. Chem. 97, (1993). 51 S.-Y. Lee, Z. W. Feng, and R. C. K. Yeo, J. Raman Spec. 28, 411 (1997). 52 G. Beylkin, SIAM J. Numer. Anal. 6, 1716 (1992). 53 B. R. Johnson and J. L. Kinsey, J. Chem. Phys. 99, 7267 (1993). 54 B. R. Johnson and J. L. Kinsey, in Femtosecond Chemistry, edited by J. Manz and L. Wöste (Verlag Chemie, New York, 1995), pp M. Shapiro, Chem. Phys. Lett. 242, 548 (1995). 56 X.-P. Peng, L.-S. Tian, and Y.-N. Peng, IEEE Trans. Signal Processing 44, 129 (1996). 57 P. Abry and P. Flandrin, IEEE Sig. Proc. Lett. 1, 32 (1994). 58 R. Glowinski, W. M. Lawton, M. Ravachol, and E. Tenenbaum, Wavelet solutions of linear and nonlinear elliptic, parabolic and hyperbolic problems in one space dimension, AD890527, Aware, Inc. (1989). 59 J. Ko, A. J. Kurdila, R. O. Wells, and X. Zhou, Comm. Num. Methods Eng. 12, 281 (1996). 60 J. R. Williams and K. Amaratunga, J. Comp. Phys. 122, 30 (1995). 27
28 61 W. Dahmen, A. Kunoth, and K. Urban, Computing 56, 259 (1996). 62 W. P. Reinhardt, Int. J. Quant. Chem. 21, 133 (1982). 63 K. McCormick and P. J. Nordlander, (unpublished). 64 M. Unser, IEEE Trans. Signal Proc. 44, 519 (1996). 28
29 Tables Table 1. Scaling coefficients and lowest moments for the L = 10 Daubechies scaling functions. k h k m k e e e e e e e3 Table 2. Convergence of Morse potential matrix element calculations in adaptive quadrature. For each scale, Max, the maximum discrepancy between r = 6 and r = 7 evaluations of the matrix elements, occurs for k' = k = 0. Only approximately half of the samples needed at each scale refinement require new evaluations. No. of No. of New Scale Max Tolerance Samples Samples e e e e e e e e e e
30 Figures (x) (x) x Figure 1. Scaling function φ (x ) and wavelet ψ(x ) for Daubechies family with L = x10 3 Frequency (cm 1 ) Figure 2. Hartley absorption band of ozone measured at 218 K by Brion, et al., 47 with smoothed reconstruction obtained from expansion in scaling functions at scale cm 1. 30
31 = x10 3 Frequency (cm 1 ) Figure 3. The significant detail components of the Hartley absorption cross-section at scales 3. 2 obtained by expansion in wavelets ψ k and summation over all k contributions for each. Vertical offsets are given to traces other than that for = 6. = Frequency (cm 1 ) 44x10 3 Figure 4. The wavelet coefficients corresponding to the components of the Hartley absorption cross-section shown in Fig. 3. For each and k, the position on the frequency scale marks the left hand side of the corresponding wavelet. These positions are marked by dots for the less congested scales. Vertical offsets are given to coefficients other than those for = 6. 31
32 k Quadrature order r Figure 5. Absolute values of the differences for order r and order r + 1 quadrature evaluation of proection integrals φ k V for the Morse potential. The L = 10 Daubechies scaling functions of spacing x = 0.4 Å are used. The sample spacing here is the same as the basis function spacing. V(x) φ k x (Å) 2 3 Figure 6. Morse potential and L = 10 Daubechies scaling functions for which proection integrals are calculated. For each scale, the + signs represent the left hand endpoints of scaling functions included in the adaptive quadrature. The greatest detail is required in the vicinity of the steep repulsive wall. 32
33 k 2 λ (Å) 3 4 Figure 7. Required integrals φ k V φ k from adaptive calculation of the Morse potential matrix elements. Each symbol represents a matrix element involving the product φ k φ k ; the positions along the left and bottom axes represent the initial points of the k' and k scaling functions, respectively. The circles correspond to the target integrals on level = 4, while the black dots correspond to the matrix elements required on finer levels. Symmetry reduces the number of independent integrals required by approximately a factor of 2. 33
Quadrature Prefilters for the Discrete Wavelet Transform. Bruce R. Johnson. James L. Kinsey. Abstract
Quadrature Prefilters for the Discrete Wavelet Transform Bruce R. Johnson James L. Kinsey Abstract Discrepancies between the Discrete Wavelet Transform and the coefficients of the Wavelet Series are known
More informationSolution of Cartesian and Curvilinear Quantum. Equations via Multiwavelets on the Interval
Solution of Cartesian and Curvilinear Quantum Equations via Multiwavelets on the Interval Bruce R. Johnson, Jeffrey. Mackey, a and James. Kinsey Department of Chemistry and Rice Quantum Institute Rice
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 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 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 informationMoments conservation in adaptive Vlasov solver
12 Moments conservation in adaptive Vlasov solver M. Gutnic a,c, M. Haefele b,c and E. Sonnendrücker a,c a IRMA, Université Louis Pasteur, Strasbourg, France. b LSIIT, Université Louis Pasteur, Strasbourg,
More informationQUADRATURES INVOLVING POLYNOMIALS AND DAUBECHIES' WAVELETS *
QUADRATURES INVOLVING POLYNOMIALS AND DAUBECHIES' WAVELETS * WEI-CHANG SHANN AND JANN-CHANG YAN Abstract. Scaling equations are used to derive formulae of quadratures involving polynomials and scaling/wavelet
More 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 informationNinth International Water Technology Conference, IWTC9 2005, Sharm El-Sheikh, Egypt 673
Ninth International Water Technology Conference, IWTC9 2005, Sharm El-Sheikh, Egypt 673 A NEW NUMERICAL APPROACH FOR THE SOLUTION OF CONTAMINANT TRANSPORT EQUATION Mohamed El-Gamel Department of Mathematical
More informationThe Application of Legendre Multiwavelet Functions in Image Compression
Journal of Modern Applied Statistical Methods Volume 5 Issue 2 Article 3 --206 The Application of Legendre Multiwavelet Functions in Image Compression Elham Hashemizadeh Department of Mathematics, Karaj
More informationOn the fast algorithm for multiplication of functions in the wavelet bases
Published in Proceedings of the International Conference Wavelets and Applications, Toulouse, 1992; Y. Meyer and S. Roques, edt., Editions Frontieres, 1993 On the fast algorithm for multiplication of functions
More 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 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 informationIntroduction to Wavelets and Wavelet Transforms
Introduction to Wavelets and Wavelet Transforms A Primer C. Sidney Burrus, Ramesh A. Gopinath, and Haitao Guo with additional material and programs by Jan E. Odegard and Ivan W. Selesnick Electrical and
More informationWavelet Transform. Figure 1: Non stationary signal f(t) = sin(100 t 2 ).
Wavelet Transform Andreas Wichert Department of Informatics INESC-ID / IST - University of Lisboa Portugal andreas.wichert@tecnico.ulisboa.pt September 3, 0 Short Term Fourier Transform Signals whose frequency
More informationSymmetric Wavelet Tight Frames with Two Generators
Symmetric Wavelet Tight Frames with Two Generators Ivan W. Selesnick Electrical and Computer Engineering Polytechnic University 6 Metrotech Center, Brooklyn, NY 11201, USA tel: 718 260-3416, fax: 718 260-3906
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 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 in curvilinear coordinate quantum
Wavelets in curvilinear coordinate quantum calculations: H + electronic states A. Maloney, a James L. Kinsey and Bruce R. Johnson Department of Chemistry and Rice Quantum Institute MS 600, Rice University,
More informationWavelet Bi-frames with Uniform Symmetry for Curve Multiresolution Processing
Wavelet Bi-frames with Uniform Symmetry for Curve Multiresolution Processing Qingtang Jiang Abstract This paper is about the construction of univariate wavelet bi-frames with each framelet being symmetric.
More informationA Tutorial on Wavelets and their Applications. Martin J. Mohlenkamp
A Tutorial on Wavelets and their Applications Martin J. Mohlenkamp University of Colorado at Boulder Department of Applied Mathematics mjm@colorado.edu This tutorial is designed for people with little
More informationarxiv:quant-ph/ v1 29 Mar 2003
Finite-Dimensional PT -Symmetric Hamiltonians arxiv:quant-ph/0303174v1 29 Mar 2003 Carl M. Bender, Peter N. Meisinger, and Qinghai Wang Department of Physics, Washington University, St. Louis, MO 63130,
More informationSOLUTION OF 1D AND 2D POISSON'S EQUATION BY USING WAVELET SCALING FUNCTIONS
SOLUION OF D AND D POISSON'S EQUAION BY USING WAVELE SCALING FUNCIONS R. B. Burgos a, and H. F. C. Peixoto b a Universidade do Estado do Rio de Janeiro Departamento de Estruturas e Fundações Rua S. Francisco
More informationFrame Diagonalization of Matrices
Frame Diagonalization of Matrices Fumiko Futamura Mathematics and Computer Science Department Southwestern University 00 E University Ave Georgetown, Texas 78626 U.S.A. Phone: + (52) 863-98 Fax: + (52)
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 informationINFINITUDE OF MINIMALLY SUPPORTED TOTALLY INTERPOLATING BIORTHOGONAL MULTIWAVELET SYSTEMS WITH LOW APPROXIMATION ORDERS. Youngwoo Choi and Jaewon Jung
Korean J. Math. (0) No. pp. 7 6 http://dx.doi.org/0.68/kjm.0...7 INFINITUDE OF MINIMALLY SUPPORTED TOTALLY INTERPOLATING BIORTHOGONAL MULTIWAVELET SYSTEMS WITH LOW APPROXIMATION ORDERS Youngwoo Choi and
More informationDiscrete Simulation of Power Law Noise
Discrete Simulation of Power Law Noise Neil Ashby 1,2 1 University of Colorado, Boulder, CO 80309-0390 USA 2 National Institute of Standards and Technology, Boulder, CO 80305 USA ashby@boulder.nist.gov
More informationCS 450 Numerical Analysis. Chapter 8: Numerical Integration and Differentiation
Lecture slides based on the textbook Scientific Computing: An Introductory Survey by Michael T. Heath, copyright c 2018 by the Society for Industrial and Applied Mathematics. http://www.siam.org/books/cl80
More information1 Introduction to Wavelet Analysis
Jim Lambers ENERGY 281 Spring Quarter 2007-08 Lecture 9 Notes 1 Introduction to Wavelet Analysis Wavelets were developed in the 80 s and 90 s as an alternative to Fourier analysis of signals. Some of the
More informationNonstationary Subdivision Schemes and Totally Positive Refinable Functions
Nonstationary Subdivision Schemes and Totally Positive Refinable Functions Laura Gori and Francesca Pitolli December, 2007 Abstract In this paper we construct a class of totally positive refinable functions,
More informationOn the Hilbert Transform of Wavelets
On the Hilbert Transform of Wavelets Kunal Narayan Chaudhury and Michael Unser Abstract A wavelet is a localized function having a prescribed number of vanishing moments. In this correspondence, we provide
More informationSchemes. Philipp Keding AT15, San Antonio, TX. Quarklet Frames in Adaptive Numerical. Schemes. Philipp Keding. Philipps-University Marburg
AT15, San Antonio, TX 5-23-2016 Joint work with Stephan Dahlke and Thorsten Raasch Let H be a Hilbert space with its dual H. We are interested in the solution u H of the operator equation Lu = f, with
More information1 Mathematical preliminaries
1 Mathematical preliminaries The mathematical language of quantum mechanics is that of vector spaces and linear algebra. In this preliminary section, we will collect the various definitions and mathematical
More information( nonlinear constraints)
Wavelet Design & Applications Basic requirements: Admissibility (single constraint) Orthogonality ( nonlinear constraints) Sparse Representation Smooth functions well approx. by Fourier High-frequency
More informationWavelets in Scattering Calculations
Wavelets in Scattering Calculations W. P., Brian M. Kessler, Gerald L. Payne polyzou@uiowa.edu The University of Iowa Wavelets in Scattering Calculations p.1/43 What are Wavelets? Orthonormal basis functions.
More informationLifting Parameterisation of the 9/7 Wavelet Filter Bank and its Application in Lossless Image Compression
Lifting Parameterisation of the 9/7 Wavelet Filter Bank and its Application in Lossless Image Compression TILO STRUTZ Deutsche Telekom AG, Hochschule für Telekommunikation Institute of Communications Engineering
More informationToward a Realization of Marr s Theory of Primal Sketches via Autocorrelation Wavelets: Image Representation using Multiscale Edge Information
Toward a Realization of Marr s Theory of Primal Sketches via Autocorrelation Wavelets: Image Representation using Multiscale Edge Information Naoki Saito 1 Department of Mathematics University of California,
More informationConstruction of wavelets. Rob Stevenson Korteweg-de Vries Institute for Mathematics University of Amsterdam
Construction of wavelets Rob Stevenson Korteweg-de Vries Institute for Mathematics University of Amsterdam Contents Stability of biorthogonal wavelets. Examples on IR, (0, 1), and (0, 1) n. General domains
More informationWavelets and Multiresolution Processing
Wavelets and Multiresolution Processing Wavelets Fourier transform has it basis functions in sinusoids Wavelets based on small waves of varying frequency and limited duration In addition to frequency,
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 informationParametrizing orthonormal wavelets by moments
Parametrizing orthonormal wavelets by moments 2 1.5 1 0.5 0 0.5 1 1.5 0 0.5 1 1.5 2 2.5 3 1.4 1.2 1 0.8 0.6 0.4 0.2 0 0.2 0.4 0 1 2 3 4 5 Georg Regensburger Johann Radon Institute for Computational and
More informationCOMPACTLY SUPPORTED ORTHONORMAL COMPLEX WAVELETS WITH DILATION 4 AND SYMMETRY
COMPACTLY SUPPORTED ORTHONORMAL COMPLEX WAVELETS WITH DILATION 4 AND SYMMETRY BIN HAN AND HUI JI Abstract. In this paper, we provide a family of compactly supported orthonormal complex wavelets with dilation
More 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 informationORTHONORMAL SAMPLING FUNCTIONS
ORTHONORMAL SAMPLING FUNCTIONS N. KAIBLINGER AND W. R. MADYCH Abstract. We investigate functions φ(x) whose translates {φ(x k)}, where k runs through the integer lattice Z, provide a system of orthonormal
More informationThe quantum state as a vector
The quantum state as a vector February 6, 27 Wave mechanics In our review of the development of wave mechanics, we have established several basic properties of the quantum description of nature:. A particle
More informationOn Riesz-Fischer sequences and lower frame bounds
On Riesz-Fischer sequences and lower frame bounds P. Casazza, O. Christensen, S. Li, A. Lindner Abstract We investigate the consequences of the lower frame condition and the lower Riesz basis condition
More informationDigital Image Processing
Digital Image Processing Wavelets and Multiresolution Processing (Wavelet Transforms) Christophoros Nikou cnikou@cs.uoi.gr University of Ioannina - Department of Computer Science 2 Contents Image pyramids
More informationA Reverse Technique for Lumping High Dimensional Model Representation Method
A Reverse Technique for Lumping High Dimensional Model Representation Method M. ALPER TUNGA Bahçeşehir University Department of Software Engineering Beşiktaş, 34349, İstanbul TURKEY TÜRKİYE) alper.tunga@bahcesehir.edu.tr
More informationRecent Developments in Numerical Methods for 4d-Var
Recent Developments in Numerical Methods for 4d-Var Mike Fisher Slide 1 Recent Developments Numerical Methods 4d-Var Slide 2 Outline Non-orthogonal wavelets on the sphere: - Motivation: Covariance Modelling
More informationWavelets in Image Compression
Wavelets in Image Compression M. Victor WICKERHAUSER Washington University in St. Louis, Missouri victor@math.wustl.edu http://www.math.wustl.edu/~victor THEORY AND APPLICATIONS OF WAVELETS A Workshop
More informationMathematical Introduction
Chapter 1 Mathematical Introduction HW #1: 164, 165, 166, 181, 182, 183, 1811, 1812, 114 11 Linear Vector Spaces: Basics 111 Field A collection F of elements a,b etc (also called numbers or scalars) with
More informationElliptic Problems / Multigrid. PHY 604: Computational Methods for Physics and Astrophysics II
Elliptic Problems / Multigrid Summary of Hyperbolic PDEs We looked at a simple linear and a nonlinear scalar hyperbolic PDE There is a speed associated with the change of the solution Explicit methods
More informationDigital Image Processing
Digital Image Processing Wavelets and Multiresolution Processing () Christophoros Nikou cnikou@cs.uoi.gr University of Ioannina - Department of Computer Science 2 Contents Image pyramids Subband coding
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 informationMATH 590: Meshfree Methods
MATH 590: Meshfree Methods Chapter 33: Adaptive Iteration Greg Fasshauer Department of Applied Mathematics Illinois Institute of Technology Fall 2010 fasshauer@iit.edu MATH 590 Chapter 33 1 Outline 1 A
More informationMATH 590: Meshfree Methods
MATH 590: Meshfree Methods Chapter 33: Adaptive Iteration Greg Fasshauer Department of Applied Mathematics Illinois Institute of Technology Fall 2010 fasshauer@iit.edu MATH 590 Chapter 33 1 Outline 1 A
More informationFast matrix algebra for dense matrices with rank-deficient off-diagonal blocks
CHAPTER 2 Fast matrix algebra for dense matrices with rank-deficient off-diagonal blocks Chapter summary: The chapter describes techniques for rapidly performing algebraic operations on dense matrices
More informationChem 3502/4502 Physical Chemistry II (Quantum Mechanics) 3 Credits Spring Semester 2006 Christopher J. Cramer. Lecture 9, February 8, 2006
Chem 3502/4502 Physical Chemistry II (Quantum Mechanics) 3 Credits Spring Semester 2006 Christopher J. Cramer Lecture 9, February 8, 2006 The Harmonic Oscillator Consider a diatomic molecule. Such a molecule
More informationPhysics 221A Fall 2018 Notes 22 Bound-State Perturbation Theory
Copyright c 2018 by Robert G. Littlejohn Physics 221A Fall 2018 Notes 22 Bound-State Perturbation Theory 1. Introduction Bound state perturbation theory applies to the bound states of perturbed systems,
More informationAnalysis of dynamic light scattering signals with discrete wavelet transformation
EUROPHYSICS LETTERS 10 September 1996 Europhys. Lett., 35 (8), pp. 621-626 (1996) Analysis of dynamic light scattering signals with discrete wavelet transformation M. Kroon( ),G.H.Wegdamand R. Sprik Van
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 informationMGA Tutorial, September 08, 2004 Construction of Wavelets. Jan-Olov Strömberg
MGA Tutorial, September 08, 2004 Construction of Wavelets Jan-Olov Strömberg Department of Mathematics Royal Institute of Technology (KTH) Stockholm, Sweden Department of Numerical Analysis and Computer
More informationApplications of Polyspline Wavelets to Astronomical Image Analysis
VIRTUAL OBSERVATORY: Plate Content Digitization, Archive Mining & Image Sequence Processing edited by M. Tsvetkov, V. Golev, F. Murtagh, and R. Molina, Heron Press, Sofia, 25 Applications of Polyspline
More informationTime series denoising with wavelet transform
Paper Time series denoising with wavelet transform Bartosz Kozłowski Abstract This paper concerns the possibilities of applying wavelet analysis to discovering and reducing distortions occurring in time
More informationThe discrete and fast Fourier transforms
The discrete and fast Fourier transforms Marcel Oliver Revised April 7, 1 1 Fourier series We begin by recalling the familiar definition of the Fourier series. For a periodic function u: [, π] C, we define
More informationA Novel Fast Computing Method for Framelet Coefficients
American Journal of Applied Sciences 5 (11): 15-157, 008 ISSN 1546-939 008 Science Publications A Novel Fast Computing Method for Framelet Coefficients Hadeel N. Al-Taai Department of Electrical and Electronic
More informationPage 404. Lecture 22: Simple Harmonic Oscillator: Energy Basis Date Given: 2008/11/19 Date Revised: 2008/11/19
Page 404 Lecture : Simple Harmonic Oscillator: Energy Basis Date Given: 008/11/19 Date Revised: 008/11/19 Coordinate Basis Section 6. The One-Dimensional Simple Harmonic Oscillator: Coordinate Basis Page
More informationWe have to prove now that (3.38) defines an orthonormal wavelet. It belongs to W 0 by Lemma and (3.55) with j = 1. We can write any f W 1 as
88 CHAPTER 3. WAVELETS AND APPLICATIONS We have to prove now that (3.38) defines an orthonormal wavelet. It belongs to W 0 by Lemma 3..7 and (3.55) with j =. We can write any f W as (3.58) f(ξ) = p(2ξ)ν(2ξ)
More informationMultidimensional partitions of unity and Gaussian terrains
and Gaussian terrains Richard A. Bale, Jeff P. Grossman, Gary F. Margrave, and Michael P. Lamoureu ABSTRACT Partitions of unity play an important rôle as amplitude-preserving windows for nonstationary
More informationQuantum Mechanics- I Prof. Dr. S. Lakshmi Bala Department of Physics Indian Institute of Technology, Madras
Quantum Mechanics- I Prof. Dr. S. Lakshmi Bala Department of Physics Indian Institute of Technology, Madras Lecture - 4 Postulates of Quantum Mechanics I In today s lecture I will essentially be talking
More informationIntroduction to Electronic Structure Theory
Introduction to Electronic Structure Theory C. David Sherrill School of Chemistry and Biochemistry Georgia Institute of Technology June 2002 Last Revised: June 2003 1 Introduction The purpose of these
More informationNonlinear Stationary Subdivision
Nonlinear Stationary Subdivision Michael S. Floater SINTEF P. O. Box 4 Blindern, 034 Oslo, Norway E-mail: michael.floater@math.sintef.no Charles A. Micchelli IBM Corporation T.J. Watson Research Center
More information1 Infinite-Dimensional Vector Spaces
Theoretical Physics Notes 4: Linear Operators In this installment of the notes, we move from linear operators in a finitedimensional vector space (which can be represented as matrices) to linear operators
More informationLinear Algebra (part 1) : Vector Spaces (by Evan Dummit, 2017, v. 1.07) 1.1 The Formal Denition of a Vector Space
Linear Algebra (part 1) : Vector Spaces (by Evan Dummit, 2017, v. 1.07) Contents 1 Vector Spaces 1 1.1 The Formal Denition of a Vector Space.................................. 1 1.2 Subspaces...................................................
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 informationPerformance Evaluation of Generalized Polynomial Chaos
Performance Evaluation of Generalized Polynomial Chaos Dongbin Xiu, Didier Lucor, C.-H. Su, and George Em Karniadakis 1 Division of Applied Mathematics, Brown University, Providence, RI 02912, USA, gk@dam.brown.edu
More informationNUMERICAL METHODS FOR ENGINEERING APPLICATION
NUMERICAL METHODS FOR ENGINEERING APPLICATION Second Edition JOEL H. FERZIGER A Wiley-Interscience Publication JOHN WILEY & SONS, INC. New York / Chichester / Weinheim / Brisbane / Singapore / Toronto
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 informationCOMPUTATION OF BESSEL AND AIRY FUNCTIONS AND OF RELATED GAUSSIAN QUADRATURE FORMULAE
BIT 6-85//41-11 $16., Vol. 4, No. 1, pp. 11 118 c Swets & Zeitlinger COMPUTATION OF BESSEL AND AIRY FUNCTIONS AND OF RELATED GAUSSIAN QUADRATURE FORMULAE WALTER GAUTSCHI Department of Computer Sciences,
More informationWavelet transforms and compression of seismic data
Mathematical Geophysics Summer School, Stanford, August 1998 Wavelet transforms and compression of seismic data G. Beylkin 1 and A. Vassiliou 2 I Introduction Seismic data compression (as it exists today)
More informationaxioms Construction of Multiwavelets on an Interval Axioms 2013, 2, ; doi: /axioms ISSN
Axioms 2013, 2, 122-141; doi:10.3390/axioms2020122 Article OPEN ACCESS axioms ISSN 2075-1680 www.mdpi.com/journal/axioms Construction of Multiwavelets on an Interval Ahmet Altürk 1 and Fritz Keinert 2,
More informationQuantum Mechanics- I Prof. Dr. S. Lakshmi Bala Department of Physics Indian Institute of Technology, Madras
Quantum Mechanics- I Prof. Dr. S. Lakshmi Bala Department of Physics Indian Institute of Technology, Madras Lecture - 6 Postulates of Quantum Mechanics II (Refer Slide Time: 00:07) In my last lecture,
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 informationBesov regularity for operator equations on patchwise smooth manifolds
on patchwise smooth manifolds Markus Weimar Philipps-University Marburg Joint work with Stephan Dahlke (PU Marburg) Mecklenburger Workshop Approximationsmethoden und schnelle Algorithmen Hasenwinkel, March
More informationFINITE-DIMENSIONAL LINEAR ALGEBRA
DISCRETE MATHEMATICS AND ITS APPLICATIONS Series Editor KENNETH H ROSEN FINITE-DIMENSIONAL LINEAR ALGEBRA Mark S Gockenbach Michigan Technological University Houghton, USA CRC Press Taylor & Francis Croup
More informationChem 3502/4502 Physical Chemistry II (Quantum Mechanics) 3 Credits Fall Semester 2006 Christopher J. Cramer. Lecture 5, January 27, 2006
Chem 3502/4502 Physical Chemistry II (Quantum Mechanics) 3 Credits Fall Semester 2006 Christopher J. Cramer Lecture 5, January 27, 2006 Solved Homework (Homework for grading is also due today) We are told
More informationDigital Image Processing
Digital Image Processing, 2nd ed. Digital Image Processing Chapter 7 Wavelets and Multiresolution Processing Dr. Kai Shuang Department of Electronic Engineering China University of Petroleum shuangkai@cup.edu.cn
More informationVector Spaces in Quantum Mechanics
Chapter 8 Vector Spaces in Quantum Mechanics We have seen in the previous Chapter that there is a sense in which the state of a quantum system can be thought of as being made up of other possible states.
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 informationPhysics 221A Fall 2017 Notes 27 The Variational Method
Copyright c 2018 by Robert G. Littlejohn Physics 221A Fall 2017 Notes 27 The Variational Method 1. Introduction Very few realistic problems in quantum mechanics are exactly solvable, so approximation methods
More informationScientific Computing: An Introductory Survey
Scientific Computing: An Introductory Survey Chapter 12 Prof. Michael T. Heath Department of Computer Science University of Illinois at Urbana-Champaign Copyright c 2002. Reproduction permitted for noncommercial,
More informationAPPROXIMATING THE INVERSE FRAME OPERATOR FROM LOCALIZED FRAMES
APPROXIMATING THE INVERSE FRAME OPERATOR FROM LOCALIZED FRAMES GUOHUI SONG AND ANNE GELB Abstract. This investigation seeks to establish the practicality of numerical frame approximations. Specifically,
More informationExploring the energy landscape
Exploring the energy landscape ChE210D Today's lecture: what are general features of the potential energy surface and how can we locate and characterize minima on it Derivatives of the potential energy
More informationChem 3502/4502 Physical Chemistry II (Quantum Mechanics) 3 Credits Spring Semester 2006 Christopher J. Cramer. Lecture 7, February 1, 2006
Chem 350/450 Physical Chemistry II (Quantum Mechanics) 3 Credits Spring Semester 006 Christopher J. Cramer ecture 7, February 1, 006 Solved Homework We are given that A is a Hermitian operator such that
More informationThe Discrete Kalman Filtering of a Class of Dynamic Multiscale Systems
668 IEEE TRANSACTIONS ON CIRCUITS AND SYSTEMS II: ANALOG AND DIGITAL SIGNAL PROCESSING, VOL 49, NO 10, OCTOBER 2002 The Discrete Kalman Filtering of a Class of Dynamic Multiscale Systems Lei Zhang, Quan
More informationPhysics 202 Laboratory 5. Linear Algebra 1. Laboratory 5. Physics 202 Laboratory
Physics 202 Laboratory 5 Linear Algebra Laboratory 5 Physics 202 Laboratory We close our whirlwind tour of numerical methods by advertising some elements of (numerical) linear algebra. There are three
More informationPage 52. Lecture 3: Inner Product Spaces Dual Spaces, Dirac Notation, and Adjoints Date Revised: 2008/10/03 Date Given: 2008/10/03
Page 5 Lecture : Inner Product Spaces Dual Spaces, Dirac Notation, and Adjoints Date Revised: 008/10/0 Date Given: 008/10/0 Inner Product Spaces: Definitions Section. Mathematical Preliminaries: Inner
More informationA Wavelet-Based Technique for Identifying, Labeling, and Tracking of Ocean Eddies
MARCH 2002 LUO AND JAMESON 381 A Wavelet-Based Technique for Identifying, Labeling, and Tracking of Ocean Eddies JINGJIA LUO Department of Earth and Planetary Physics, Graduate School of Science, University
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 informationWavelet Filter Transforms in Detail
Wavelet Filter Transforms in Detail Wei ZHU and M. Victor WICKERHAUSER Washington University in St. Louis, Missouri victor@math.wustl.edu http://www.math.wustl.edu/~victor FFT 2008 The Norbert Wiener Center
More information