QUADRATURES INVOLVING POLYNOMIALS AND DAUBECHIES' WAVELETS *
|
|
- Charles Murphy
- 5 years ago
- Views:
Transcription
1 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 functions with compact supports; in particular, those discovered by Daubechies. It turns out that with a few parameters, which are theoretically exact, these quadratures can be evaluated with algebraic formulae instead of numerical approximations. Those parameters can be obtained with high precision by solving well-conditioned linear systems of equations which involve matrices already seen in the literature of wavelets for other purposes. Key words. scaling equation, wavelet, quadrature, polynomial, Neumann series Subject classications: AMS(MOS): 65A5, 65D3, 65F35, 4C5. Introduction. Because of their similarities with Fourier transforms and nite elements, wavelets are studied for the applicability of being a tool in scientic computation and numerical analysis, such as the numerical solution of dierential equations. Among the recent articles on the later subject, we refer to the works by Bacry, Mallat and Papanicolaou [], Beylkin, Coifman and Rokhlin [4], Glowinski, Lawton, Ravachol and Tenenbaum [8], Jaard [9], Maday, Perrier and Ravel [4], and u and Shann [8]. One of the fundamental computing steps of these applications is the numerical quadratures f(x)(x) dx. As pointed out by Dahmen and Micchelli [5], since the derivatives of also satisfy (other) scaling equations, the techniques derived for the foregoing integrals can be carried to those with derivatives of. Beylkin, Coifman and Rokhlin [] studied numerical methods for the rst two kinds of integrals. There the idea is that (x) can be related to ( j x R? k) by the fast pyramid algorithm [3]; and the integral f(x)( j x? k) dx can be approximated by numerical quadrature rules if j is large enough (that is, if the support of ( j x) is narrow enough). The fundation for the derivation of such quadrature points and parameters is the integrals for which f(x) are polynomials. That is, the moments of the forms R f(x)(x) dx, R f(x)(x)(x? k) dx, or R b a (:) x m (x) dx; and the integrals (:) x m (x)(x? k) dx; (:3) b a x m (x) dx: In [], (:) were approximated by a recursively generated sequence of vectors, based on the -periodic characteristic function of ^(). Sweldens and Piessens [7] went on to show a better * Technical Report 93, Department of Mathematics, National Central University, Taiwan. Department of Mathematics, National Central University, Chung-Li, Taiwan, R.O.C. This work was supported by the National Science Coucil in Taiwan; grant # 8-8-M address: shann@math.ncu.edu.tw. address: yan@math.ncu.edu.tw
2 W. Shann and C. Yen approximation result for the one-point rule in [], and to extend the multi-point rules. They also started with the relation, which states that (:) can be done theoretically exact. Beylkin [3] found the integrals (:) are associated with matrices. In particular, he found the values of R (x) (x? k) dx consist a normalized eigenvector of a certain matrix. This technique turns out to be a special case of the theorems in [5], which guarateed this sort of integrals always associate with eigenproblems and there are always unique solutions with appropriate assumptions. In this article we extend the idea to treat (:3). We also show that the matrices associated with these integrals have uniform bounds in vector `-norm. That is, the upper bound of the norms do not change with the order of wavelets. In particular, the assocated linear systems can be solved with high precision eciently by Neumann series. Wavelets in this article are those discovered by Daubechies [6], which have compact supports and form an orthonormal basis of L (R). In this case, the wavelet (x) is derived from the scaling function (x) which is normalized to have R = and satises the scaling equation (:4) (x) = k c k (x? k): The coecients c k were computed by constructing a certain trigonometric polynomial m () = P ck e?ik, see Daubechies' \Ten Lectures" [7]. However, Strang pointed out [6] that they can also be solved by the following relations: (:5) c k = for k 6 f; ; : : :; p? g; (:6) (:7) and (:8) k k c k? = ; (?) k k m c k = ; for m p? ; ( := ) k c k c k?m? m = : Here p is a given positive integer. We shall call respectively such dened c k, and the scaling coecients, scaling functions, and wavelets of order p. This parameter p plays several roles in the theory of wavelets. What we need now is that supp = [; p? ]. Once (x) is determined, the multiresolution approach asserts that (:9) (x) = k (?) k c?k (x? k) R will be a wavelet. One of the properties of is that x m dx = for m p?. Reducing the computations involving wavelets to the manipulation of c k is a general strategy. It is sought for every kind of application for wavelets. Examples have been seen in numerical algorithms and signal decomposition [, ], among others. Thus the numerical values of c k heavily determine the accuracy of these algorithms. In the next section we discuss the computation of scaling coecients. Then we derive the formulae for R x m. As the immediate applications, polynomials P k m (x? k), for m p?, R are constructed and formulae for x m are derived. Next we evaluate the integrals R on the half-lines [; ) or (?; ]. Finally we present that the similar idea can treat convolutions x m (x)(x?k).
3 Wavelets quadratures 3. Computation of scaling coecients. Note that (:8) contributes exactly p? equations in the system, namely m p? : The situations for m p and m?p are consequences from (:5); those for? p m? are dual to the cases of jmj; and the case of m = can be derived from (:6), the m = case of (:7), and the m 6= part of (:8). Therefore, we can write conditions (:6){(:8) into p equations in x = (x ; : : :; x p? ) T. Let it be a vector-valued function F p (x). Algebraic solutions for p =, and 3 are known; but it is not likely to nd those with higher p. Numerical solutions for F p, p, are listed in a table [7, p. 95] with 6 digits. 3 Let d be the solution of F listed in the table (with (:6) changed to p ). We have kf (d)k 8?6. This is due to the diculty introduced in (:7). For p = and m = 9, the smallest and the largest number in that equation is dierent by a scale of 8. Observe that the set of equations in (:7) can be replaced by (:) k (?) k (k? `) m c k = ; for m p? : for any xed integer `. In particular, we take ` = p to have the most advantage. With (:7) replaced by (:), the foregoing diculty is slightly improved. Let F ~ p (x) be the alternate system of equations. One can write down the Jacobian of F ~ p and perform the Newton's iteration to solve ~F p (x) = (we use the roots of Fp? ~ appended with two zeros as the initial guess). Our FORTRAN code is able to squeeze out more information about the scaling coecients. First, let c be the numerical solution of F ~, we have kf (c)k?8. Secondly, we can calculate the scaling coecients up to p = 4. Some results are listed in Table. All values quoted hereafter are computed in quartic precision on a VA 65. Numerical values of c k of order p 4, together with those values and programs mentioned later in this article, are available from the rst author by request. It makes more sense to query and obtain the data electronically. 3. Integrals on half lines. While solving dierential equations, the domain of interest is usually R a bounded interval, on which boundary conditions are imposed. In this situation, integrals b like a q(x)(j x? k) dx will come in order. In practice one can rescale the problem and only use those wavelets and scaling functions in such ne scales that in the support of each ( j x? k) (or ( j x? k)), there lies at most one boundary point. Since the situations are analogous at either boundary, we now suppose x = is the left boundary point and consider the integrals on [; ). As before, the core formulae for this kind of integrals are N p m;k = x m (x? k) dx: For the translation parameters k? p, the integrals are zero. For k, integration over [; ) is the same as over R. Thus N p m;k = M p m;k. We now concentrate on the cases that? p k?. Again, we start from (:4) with a changing of variable: (3:) x m (x? k) dx = m+ ` c` x m (x? (k + `)) dx: 3 Over there, equation (:6) is replaced by P c k = p. Here we follow the notation in [6].
4 4 W. Shann and C. Yen If? p k?; that is, over half of the support is still inside of [; ). Then the supports of (x? (k + `)), for?k ` p?, actually lies entirely in [; ). That means, for these `, the integrals on the right hand side of (3:) can be evaluated as in Section 3. Over all, we have i c`n p m;k+` + : N p m;k = m+ h `<?k `?k c`m p m;k+` Suppose the part of `?k is known, we have a linear system in N p m;k. Given a p (there is nothing to do with p = ) and x an m, we have unknowns N p m;k for? p k?. Now let 8 >< >: x k = N p m;?k ; b k = p? m+ `=k c`m p m;`?k ; k = ; ; : : :; p? : Let b = (b k ) T and A = (a k`) be the matrix of order p? such that a k` = c k?`. Then x = (x k ) T is the solution of the system (I? A) x = b: m+ Note that A is independent of m and it is a nite submatrix of the low pass lter matrix L dened in [6], in which it is also shown that A has an eigenvalue = with an associated eigenvector ((); (); : : :; (p? )) T. Here we claim another property for A. Proposition 3.. Let A be the matrix dened above. k k be the matrix norm induced by the Euclidean norm. We have kak p : Proof. Since A = S = S where c c... c p?3 c p?4 c A and T = c p? c p? c... c p? c p? C A : By (:8), ST T = T S T =. Thus AA T = SS T T T T : Also by (:8), we have SS T + T T T = I. Thus if is an eigenvalue of SS T with an associated eigenvector v, then and? are eigenvalues of AA T with associated eigenvectors v, respectively. Since SS T and T T T are positive semi-denite, we have. Therefore kaa T k. v and
5 Wavelets quadratures 5 (3:) Thus, for any m and p, the Neumann series n= ( A)n m+ is an ecient method for computing the inverse of I? m+ A. We have prepared values for p 4 and m 4. Those lower order values are listed in Table. (The right most integer in each column is the exponent.) The iterations are computed in quartic precision with the stop criterion being that the `-norm of error <?4. As for the right boundary point, one can translate the problem to the integral on (?; ]. Clearly x m (x? k) dx = M p m;k? N p m;k :? 4. Convolutions. R Now we consider the integrals of the form q(x)( j x? `)( k x? m) dx. If j < k, apply the scaling equation (:4) recursively on ( j x? `) will lift the scale to the same level k. Then a change of variable can go back to level j = k =. Thus in practice we do not need the formula for those 's on dierent scales. Again, after changing of variables, we can concentrate on the core: L p m;k = q(x)(x) (?x) (k) = x m (x)(x? k) dx Since L p ;k = k, we start from m =. The support of (x) of order p is [; p? ], hence for each p and m, we have unknowns L p m;k for? p k p?. Then L p m;k = (4:) = = x m n P n;` cnc` m+ h n;` c n (x? n) ` c`(x? k? `) dx R ( x+n )m (x)(x + n? k? `) dx m c n c`l p m m;k+`?n + c n c` r n;` r= The problem is recursively equivalent to a system of linear equations (4:) (I? A)x = b: m+ n r L m?r;k+`?n i : Let N = p?, then x = (L p m;?n ; : : :; Lp m;n )T, b = (b k ) and A = (a kj ) with b k = p? ` m m m+ c n c` r n= `=` r= n r L p m?r;k+`?n ; where ` = maxfn? k? N; g; ` = minfn? k + N; p? g and a kj = n c n c n+j?k ;?N j; k N: This matrix A was considered by Lawton [] for a sucient and necessary condition on c k (and hence m ()) for constructing orthonormal wavelets basis. It is noted that A has an eigenvalue
6 6 W. Shann and C. Yen with an associated eigenvector v = ( k ). The condition is that is a nondefective eigenvalue of A. Daubechies [7] summerized it as one of the three equivalent conditions. Since (x) is innitely many times dierentiable as p!, dierentiate (:4) to have (s) (x) = s k c k (s) (x? k): R Apply this equation into (x) (s) (x? k) dx and repeat the steps in (4:), we nd that the matrix A = A(p) will asymptotically have eigenvalues,,, 4, : : :, as p!. Actually A has a simple pattern. On the even-numbered columns, by (:8), all entries are except a k;k =. On the odd-numbered columns all entries are except for d j?p e k b j+p c, which is the p-element vector P R k c kc k+(p?) (x)(x? p? ) dx h = P P. P k c kc k+3 P k c kc k+ P k c kc k+ k c kc k+3. k c kc k+(p?) C A = R. (x)(x? 3 ) dx R (x)(x? ) dx R (x)(x? ) dx R (x)(x? 3 ) dx R. (x)(x? p? ) dx Proposition 4.. Let A be the matrix dened in this section. k k be the matrix norm induced by the Euclidean norm. We have kak p : C A : Proof. Dene matrices B = c c c p? c c c p? c p?... c c cp? C A (N+)(N+P ) and D = c p? c p?... c c c3 C A : c c That is, B ij = c j?i and D ij = c (p?)+j?i. Observe that (N+)(N+P ) A = D B T : By (:8), DD T is a diagonal matrix and clearly its spectrum radius is. Hence kdk = p.
7 Wavelets quadratures 7 Then Let x be a (N + P ) dimensional vector. Extend x to (~x k ) ` such that kbxk = ~x k = n xk if k (N + P ) otherwise. = = N+ i= N+ i= N+ i= = ( j;k ( k j j;k j;k 4 kxk : c j?i ~x j c j?i c k?i ~x j ~x k c j c k ~x j+i ~x k+i N+ c j c k ) ( i= N+ c k ) ( i= ~x j+i ~x k+i ) N+ ~x j+i ) = ( i= ~x k+i ) = Therefore kbk. Fortunately the linear system (4:) is not needed for m = (actually b = and I? A is singular). Hence for all cases with p (p = is an empty case) and m, equation (4:) can be solved by the Neumann series (3:). As in the previous section, we list values of L p m;k with smaller p and m on Table 3. For each m, values for k =? p; : : :; p? are listed from top down and from left to right. Except those for m =, only those for k =? p; : : :; are listed. Because L p ;k = (x + k)(x + k) (x) dx = x(x)(x + k) dx = L p ;?k : The integrals x m (x)(x? k) dx and x m (x)(x? k) dx? can be evaluated by the combination of techniques described in Sections 4 and 5. It is slightly tedious but otherwise straight forward. REFERENCES E. Bacry, S. Mallat and G. Papanicolaou, A wavelet based space-time adaptive numerical method for partial dierential equations, RAIRO Math Modelling Numer. Anal., 6: 793{834 (99) G. Beylkin R. Coifman and V. Rokhlin, Fast wavelet transforms and numerical algorithms, Comm. Pure Appl. Math., 44: 4{83 (99) 3 G. Beylkin, On the presentation of operators in bases of compactly supported wavelets, SIAM J. Numer. Anal., 6: 76{74 (99)
8 8 W. Shann and C. Yen 4 G. Beylkin, R. Coifman and V. Rokhlin, Wavelets in numerical analysis, in Wavelets and their applications, Jones and Barlett Publishers, Boston, W. Dahmen and C. A. Micchelli, Using the renement equation for evaluating integrals of wavelets, SIAM J. Num. Anal., 3: 57{537 (993) 6 I. Daubechies, Orthonormal bases of compactly supported wavelets, Comm. Pure Appl. Math., 4: 99{996 (988) 7 I. Daubechies, Ten lectures on wavelets, SIAM, Philadelphia, R. Glowinski, W. M. Lawton, M. Ravachol and E. Tenenbaum, Wavelets solution of linear and nonlinear elliptic, parabolic and hyperbolic problems in one space dimension, in Computing methods in applied sciences and engineering, SIAM, Philadelphia, S. Jaard, Wavelet methods for fast resolution of elliptic problems, SIAM J. Numer. Anal., 9: 965{986 (99) B. Jawerth and W. Sweldens, An overview of wavelet based multiresolution analysis, Math. Rev. to appear. W. Lawton, Necessary and sucient conditions for constructing orthonormal wavelet bases, J. Math. Phys., 3: 57{6 (99) S. Mallat, Multifrequency channel decompositions of images and wavelet models, IEEE Trans. Acoust. Speech Signal Process, 37: 9{ (989) 3 S. Mallat, A theory for multiresolution signal decomposition: The wavelet representation, IEEE Trans. on Patt. Anal. Mach. Intell., : 674{693 (989) 4 Y. Maday, V. Perrier and J.-C. Ravel, Adaptativite dynamique sur bases d'ondelettes pour l'approximation d'equations aux derivees partielles, C. R. Acad. Sci. Paris, Serie I, 3: 45{ 4 (99) 5 V. Perrier, Towards a method for solving partial dierential equations using wavelet bases, in Wavelets, Time-Frequency Methods and Phase Space, Springer-Verlag, Berlin, G. Strang, Wavelets and dilation equations: A brief introduction, SIAM Review, 3: 64{67 (989) 7 W. Sweldens and R. Piessens, Quadrature formulae for the calculation of the wavelet decomposition, SIAM J. Numer. Anal. to appear. 8 J. u and W.-C. Shann, Galerkin-wavelets methods for two-point boundary value problems, Numer. Math., 63: 3{44 (99)
9 Wavelets quadratures 9 p c k mantissa E ? ? ? ?? ? ? ?? ?? ? ? ?3? ? ? ? ? ? ? ? ? ? ?? ?? ? ?? ?? ? ? ?? ?3? ? ?4? ?4? ? ?5? ?6 p c k mantissa E ? ? ? ? ? ? ? ?? ? ? ?? ?3? ? ?3? ? ? ? ? ? ? ? ? ?? ? ? ?? ?? ? ? ?4? ? ? ?4? ? ?6 p c k mantissa E ? ? ? ? ? ?? ? ? ?? ?? ? ?? ?3? ? ? ?4? ? ? ? ? ?? ? ? ?? ? ?? ?? ? ? ?? ? ? ?3? ?3? ? ?4? ?5 Table. Values of c ; : : :; c p? computed by solving ~ F p, for p = 9; : : :; 4.
10 W. Shann and C. Yen m p = ? ?? ? ?? ? ? 3? ? ?3 4? ? ?3 5? ? ?4 m p = ? ? ? ? ?? ? ? ?5? ? ?? ?3? ?5 3? ? ?? ?3? ? ? ? ? ? ? ? ?6 m p = ? ?? ?? ? ? ? ? ?? ?3? ? ?6? ? ?? ? ? ?5? ?8 3? ? ?? ? ? ?5? ?8 4? ? ?? ? ? ?5? ? ? ?? ?? ? ?5? ?9 m p = ? ?? ? ? ? ? ? ? ? ?? ? ? ? ? ? ?? ? ?? ? ? ? ? ? 3? ?? ? ?? ?4? ?5? ?8? ? 4? ? ?? ? ?? ?4? ?5? ?8? ? Table. Values of N p m;k, index k is assumed from top down for k =?;?; : : :;? p.
11 Wavelets quadratures m p = ?? ? ?? ? ? ? ?3? ? ? ? ?3? ? ? m p = 3? ?5? ? ?? ? ?6? ? ?? ? ? ?? ? ? ? ?? ? ? ?? ? ? ? ?? ? ? ?? ? m p = 4? ? ? ?3? ? ?? ? ? ?5? ?5? ? ?? ? ? ? ? ?4? ?7 3? ? ?5? ?3? ? ? ? ? ? ?3? ?6 4? ? ?5? ?3? ? ? ? ? ? ?? ?5 m p = 5? ?? ?8? ?5? ? ?? ? ?? ? ?9? ? ? ?? ? ? ? ? ?? ?3? ?4? ?7? ? Table 3. Values of L p m;k, index k is assumed from top down for k = :p; : : :; p:, except for m =.
2 W. LAWTON, S. L. LEE AND ZUOWEI SHEN is called the fundamental condition, and a sequence which satises the fundamental condition will be called a fu
CONVERGENCE OF MULTIDIMENSIONAL CASCADE ALGORITHM W. LAWTON, S. L. LEE AND ZUOWEI SHEN Abstract. Necessary and sucient conditions on the spectrum of the restricted transition operators are given for the
More informationQuadrature Prefilters for the Discrete Wavelet Transform. Bruce R. Johnson. James L. Kinsey. Abstract
Quadrature Prefilters for the Discrete Wavelet Transform Bruce R. Johnson James L. Kinsey Abstract Discrepancies between the Discrete Wavelet Transform and the coefficients of the Wavelet Series are known
More informationA NOTE ON MATRIX REFINEMENT EQUATIONS. Abstract. Renement equations involving matrix masks are receiving a lot of attention these days.
A NOTE ON MATRI REFINEMENT EQUATIONS THOMAS A. HOGAN y Abstract. Renement equations involving matrix masks are receiving a lot of attention these days. They can play a central role in the study of renable
More informationOctober 7, :8 WSPC/WS-IJWMIP paper. Polynomial functions are renable
International Journal of Wavelets, Multiresolution and Information Processing c World Scientic Publishing Company Polynomial functions are renable Henning Thielemann Institut für Informatik Martin-Luther-Universität
More informationFunctions: A Fourier Approach. Universitat Rostock. Germany. Dedicated to Prof. L. Berg on the occasion of his 65th birthday.
Approximation Properties of Multi{Scaling Functions: A Fourier Approach Gerlind Plona Fachbereich Mathemati Universitat Rostoc 1851 Rostoc Germany Dedicated to Prof. L. Berg on the occasion of his 65th
More 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 information290 J.M. Carnicer, J.M. Pe~na basis (u 1 ; : : : ; u n ) consisting of minimally supported elements, yet also has a basis (v 1 ; : : : ; v n ) which f
Numer. Math. 67: 289{301 (1994) Numerische Mathematik c Springer-Verlag 1994 Electronic Edition Least supported bases and local linear independence J.M. Carnicer, J.M. Pe~na? Departamento de Matematica
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 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 informationMoment Computation in Shift Invariant Spaces. Abstract. An algorithm is given for the computation of moments of f 2 S, where S is either
Moment Computation in Shift Invariant Spaces David A. Eubanks Patrick J.Van Fleet y Jianzhong Wang ẓ Abstract An algorithm is given for the computation of moments of f 2 S, where S is either a principal
More 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 informationVector Space Basics. 1 Abstract Vector Spaces. 1. (commutativity of vector addition) u + v = v + u. 2. (associativity of vector addition)
Vector Space Basics (Remark: these notes are highly formal and may be a useful reference to some students however I am also posting Ray Heitmann's notes to Canvas for students interested in a direct computational
More information. Introduction Let M be an s s dilation matrix, that is, an s s integer matrix whose eigenvalues lie outside the closed unit disk. The compactly suppo
Multivariate Compactly Supported Fundamental Renable Functions, Duals and Biorthogonal Wavelets y Hui Ji, Sherman D. Riemenschneider and Zuowei Shen Abstract: In areas of geometric modeling and wavelets,
More informationOle Christensen 3. October 20, Abstract. We point out some connections between the existing theories for
Frames and pseudo-inverses. Ole Christensen 3 October 20, 1994 Abstract We point out some connections between the existing theories for frames and pseudo-inverses. In particular, using the pseudo-inverse
More informationPositive Denite Matrix. Ya Yan Lu 1. Department of Mathematics. City University of Hong Kong. Kowloon, Hong Kong. Abstract
Computing the Logarithm of a Symmetric Positive Denite Matrix Ya Yan Lu Department of Mathematics City University of Hong Kong Kowloon, Hong Kong Abstract A numerical method for computing the logarithm
More informationConvergence of The Multigrid Method With A Wavelet. Abstract. This new coarse grid operator is constructed using the wavelet
Convergence of The Multigrid Method With A Wavelet Coarse Grid Operator Bjorn Engquist Erding Luo y Abstract The convergence of the two-level multigrid method with a new coarse grid operator is studied.
More information2 W. LAWTON, S. L. LEE AND ZUOWEI SHEN Denote by L 1 (R d ) and L 2 (R d ) the spaces of Lebesque integrable and modulus square integrable functions d
STABILITY AND ORTHONORMALITY OF MULTIVARIATE REFINABLE FUNCTIONS W. LAWTON, S. L. LEE AND ZUOWEI SHEN Abstract. This paper characterizes the stability and orthonormality of the shifts of a multidimensional
More information. Consider the linear system dx= =! = " a b # x y! : (a) For what values of a and b do solutions oscillate (i.e., do both x(t) and y(t) pass through z
Preliminary Exam { 1999 Morning Part Instructions: No calculators or crib sheets are allowed. Do as many problems as you can. Justify your answers as much as you can but very briey. 1. For positive real
More informationThe organization of the paper is as follows. First, in sections 2-4 we give some necessary background on IFS fractals and wavelets. Then in section 5
Chaos Games for Wavelet Analysis and Wavelet Synthesis Franklin Mendivil Jason Silver Applied Mathematics University of Waterloo January 8, 999 Abstract In this paper, we consider Chaos Games to generate
More informationCENTER FOR THE MATHEMATICAL SCIENCES. Characterizations of linear independence and stability. of the shifts of a univariate renable function
UNIVERSITY OF WISCONSIN-MADISON CENTER FOR THE MATHEMATICAL SCIENCES Characterizations of linear independence and stability of the shifts of a univariate renable function in terms of its renement mask
More informationA general theory of discrete ltering. for LES in complex geometry. By Oleg V. Vasilyev AND Thomas S. Lund
Center for Turbulence Research Annual Research Briefs 997 67 A general theory of discrete ltering for ES in complex geometry By Oleg V. Vasilyev AND Thomas S. und. Motivation and objectives In large eddy
More 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 informationThe Erwin Schrodinger International Boltzmanngasse 9. Institute for Mathematical Physics A-1090 Wien, Austria
ESI The Erwin Schrodinger International Boltzmanngasse 9 Institute for Mathematical Physics A-1090 Wien, Austria Space{Frequency Adaptive Approximation for Quantum Hydrodynamic Models Silvia Bertoluzza
More informationWavelet Bases of the Interval: A New Approach
Int. Journal of Math. Analysis, Vol. 1, 2007, no. 21, 1019-1030 Wavelet Bases of the Interval: A New Approach Khaled Melkemi Department of Mathematics University of Biskra, Algeria kmelkemi@yahoo.fr Zouhir
More informationSize properties of wavelet packets generated using finite filters
Rev. Mat. Iberoamericana, 18 (2002, 249 265 Size properties of wavelet packets generated using finite filters Morten Nielsen Abstract We show that asymptotic estimates for the growth in L p (R- norm of
More informationA MULTIGRID ALGORITHM FOR. Richard E. Ewing and Jian Shen. Institute for Scientic Computation. Texas A&M University. College Station, Texas SUMMARY
A MULTIGRID ALGORITHM FOR THE CELL-CENTERED FINITE DIFFERENCE SCHEME Richard E. Ewing and Jian Shen Institute for Scientic Computation Texas A&M University College Station, Texas SUMMARY In this article,
More 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 informationTHE objective of this paper is to introduce a comparative
Proceedings of the World Congress on Engineering 23 Vol I, WCE 23, July 3-5, 23, London, U.K. Numerical Solution of Fractional Integro-differential Equation by Using Cubic B-spline Wavelets Khosrow Maleknejad,
More information2 Section 2 However, in order to apply the above idea, we will need to allow non standard intervals ('; ) in the proof. More precisely, ' and may gene
Introduction 1 A dierential intermediate value theorem by Joris van der Hoeven D pt. de Math matiques (B t. 425) Universit Paris-Sud 91405 Orsay Cedex France June 2000 Abstract Let T be the eld of grid-based
More informationSTABILITY AND LINEAR INDEPENDENCE ASSOCIATED WITH SCALING VECTORS. JIANZHONG WANG y
STABILITY AND LINEAR INDEPENDENCE ASSOCIATED WITH SCALING VECTORS JIANZHONG WANG y Abstract. In this paper, we discuss stability and linear independence of the integer translates of a scaling vector =
More informationAPPROXIMATING THE COMPLEXITY MEASURE OF. Levent Tuncel. November 10, C&O Research Report: 98{51. Abstract
APPROXIMATING THE COMPLEXITY MEASURE OF VAVASIS-YE ALGORITHM IS NP-HARD Levent Tuncel November 0, 998 C&O Research Report: 98{5 Abstract Given an m n integer matrix A of full row rank, we consider the
More informationON NUMERICAL SOLUTIONS OF PARTIAL DIFFERENTIAL EQUATIONS BY THE DECOMPOSITION METHOD. Mustafa Inc
153 Kragujevac J. Math. 26 (2004) 153 164. ON NUMERICAL SOLUTIONS OF PARTIAL DIFFERENTIAL EQUATIONS BY THE DECOMPOSITION METHOD Mustafa Inc Department of Mathematics, Firat University, 23119 Elazig Turkiye
More information1 Ordinary points and singular points
Math 70 honors, Fall, 008 Notes 8 More on series solutions, and an introduction to \orthogonal polynomials" Homework at end Revised, /4. Some changes and additions starting on page 7. Ordinary points and
More informationA theorem on summable families in normed groups. Dedicated to the professors of mathematics. L. Berg, W. Engel, G. Pazderski, and H.- W. Stolle.
Rostock. Math. Kolloq. 49, 51{56 (1995) Subject Classication (AMS) 46B15, 54A20, 54E15 Harry Poppe A theorem on summable families in normed groups Dedicated to the professors of mathematics L. Berg, W.
More informationThe Best Circulant Preconditioners for Hermitian Toeplitz Systems II: The Multiple-Zero Case Raymond H. Chan Michael K. Ng y Andy M. Yip z Abstract In
The Best Circulant Preconditioners for Hermitian Toeplitz Systems II: The Multiple-ero Case Raymond H. Chan Michael K. Ng y Andy M. Yip z Abstract In [0, 4], circulant-type preconditioners have been proposed
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 informationLecture 2: Review of Prerequisites. Table of contents
Math 348 Fall 217 Lecture 2: Review of Prerequisites Disclaimer. As we have a textbook, this lecture note is for guidance and supplement only. It should not be relied on when preparing for exams. In this
More informationThe Fast Wavelet Transform 1. Beyond Fourier transforms. Mac A. Cody 2
The Fast Wavelet Transform 1 Beyond Fourier transforms Mac A. Cody 2 The Fourier transform has been used as a reliable tool in signal analysis for many years. Invented in the early 1800s by its namesake,
More informationThe Mortar Wavelet Method Silvia Bertoluzza Valerie Perrier y October 29, 1999 Abstract This paper deals with the construction of wavelet approximatio
The Mortar Wavelet Method Silvia Bertoluzza Valerie Perrier y October 9, 1999 Abstract This paper deals with the construction of wavelet approximation spaces, in the framework of the Mortar method. We
More informationMax-Planck-Institut fur Mathematik in den Naturwissenschaften Leipzig H 2 -matrix approximation of integral operators by interpolation by Wolfgang Hackbusch and Steen Borm Preprint no.: 04 200 H 2 -Matrix
More information524 Jan-Olov Stromberg practice, this imposes strong restrictions both on N and on d. The current state of approximation theory is essentially useless
Doc. Math. J. DMV 523 Computation with Wavelets in Higher Dimensions Jan-Olov Stromberg 1 Abstract. In dimension d, a lattice grid of size N has N d points. The representation of a function by, for instance,
More informationPROOF OF TWO MATRIX THEOREMS VIA TRIANGULAR FACTORIZATIONS ROY MATHIAS
PROOF OF TWO MATRIX THEOREMS VIA TRIANGULAR FACTORIZATIONS ROY MATHIAS Abstract. We present elementary proofs of the Cauchy-Binet Theorem on determinants and of the fact that the eigenvalues of a matrix
More informationV j+1 W j+1 Synthesis. ψ j. ω j. Ψ j. V j
MORPHOLOGICAL PYRAMIDS AND WAVELETS BASED ON THE QUINCUNX LATTICE HENK J.A.M. HEIJMANS CWI, P.O. Box 94079, 1090 GB Amsterdam, The Netherlands and JOHN GOUTSIAS Center for Imaging Science, Dept. of Electrical
More informationLinear Algebra: Characteristic Value Problem
Linear Algebra: Characteristic Value Problem . The Characteristic Value Problem Let < be the set of real numbers and { be the set of complex numbers. Given an n n real matrix A; does there exist a number
More informationSpurious Chaotic Solutions of Dierential. Equations. Sigitas Keras. September Department of Applied Mathematics and Theoretical Physics
UNIVERSITY OF CAMBRIDGE Numerical Analysis Reports Spurious Chaotic Solutions of Dierential Equations Sigitas Keras DAMTP 994/NA6 September 994 Department of Applied Mathematics and Theoretical Physics
More informationConstrained Leja points and the numerical solution of the constrained energy problem
Journal of Computational and Applied Mathematics 131 (2001) 427 444 www.elsevier.nl/locate/cam Constrained Leja points and the numerical solution of the constrained energy problem Dan I. Coroian, Peter
More informationPeter Deuhard. for Symmetric Indenite Linear Systems
Peter Deuhard A Study of Lanczos{Type Iterations for Symmetric Indenite Linear Systems Preprint SC 93{6 (March 993) Contents 0. Introduction. Basic Recursive Structure 2. Algorithm Design Principles 7
More informationDEN: Linear algebra numerical view (GEM: Gauss elimination method for reducing a full rank matrix to upper-triangular
form) Given: matrix C = (c i,j ) n,m i,j=1 ODE and num math: Linear algebra (N) [lectures] c phabala 2016 DEN: Linear algebra numerical view (GEM: Gauss elimination method for reducing a full rank matrix
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 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 informationNumerical Results on the Transcendence of Constants. David H. Bailey 275{281
Numerical Results on the Transcendence of Constants Involving e, and Euler's Constant David H. Bailey February 27, 1987 Ref: Mathematics of Computation, vol. 50, no. 181 (Jan. 1988), pg. 275{281 Abstract
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 informationHow smooth is the smoothest function in a given refinable space? Albert Cohen, Ingrid Daubechies, Amos Ron
How smooth is the smoothest function in a given refinable space? Albert Cohen, Ingrid Daubechies, Amos Ron A closed subspace V of L 2 := L 2 (IR d ) is called PSI (principal shift-invariant) if it is the
More information1 Outline Part I: Linear Programming (LP) Interior-Point Approach 1. Simplex Approach Comparison Part II: Semidenite Programming (SDP) Concludin
Sensitivity Analysis in LP and SDP Using Interior-Point Methods E. Alper Yldrm School of Operations Research and Industrial Engineering Cornell University Ithaca, NY joint with Michael J. Todd INFORMS
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 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 information1 Vectors. Notes for Bindel, Spring 2017 Numerical Analysis (CS 4220)
Notes for 2017-01-30 Most of mathematics is best learned by doing. Linear algebra is no exception. You have had a previous class in which you learned the basics of linear algebra, and you will have plenty
More informationPerturbation results for nearly uncoupled Markov. chains with applications to iterative methods. Jesse L. Barlow. December 9, 1992.
Perturbation results for nearly uncoupled Markov chains with applications to iterative methods Jesse L. Barlow December 9, 992 Abstract The standard perturbation theory for linear equations states that
More informationthen kaxk 1 = j a ij x j j ja ij jjx j j: Changing the order of summation, we can separate the summands, kaxk 1 ja ij jjx j j: let then c = max 1jn ja
Homework Haimanot Kassa, Jeremy Morris & Isaac Ben Jeppsen October 7, 004 Exercise 1 : We can say that kxk = kx y + yk And likewise So we get kxk kx yk + kyk kxk kyk kx yk kyk = ky x + xk kyk ky xk + kxk
More informationQUASI-UNIFORMLY POSITIVE OPERATORS IN KREIN SPACE. Denitizable operators in Krein spaces have spectral properties similar to those
QUASI-UNIFORMLY POSITIVE OPERATORS IN KREIN SPACE BRANKO CURGUS and BRANKO NAJMAN Denitizable operators in Krein spaces have spectral properties similar to those of selfadjoint operators in Hilbert spaces.
More informationLecture 7 Multiresolution Analysis
David Walnut Department of Mathematical Sciences George Mason University Fairfax, VA USA Chapman Lectures, Chapman University, Orange, CA Outline Definition of MRA in one dimension Finding the wavelet
More 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 information10.34: Numerical Methods Applied to Chemical Engineering. Lecture 2: More basics of linear algebra Matrix norms, Condition number
10.34: Numerical Methods Applied to Chemical Engineering Lecture 2: More basics of linear algebra Matrix norms, Condition number 1 Recap Numerical error Operations Properties Scalars, vectors, and matrices
More information4.1 Eigenvalues, Eigenvectors, and The Characteristic Polynomial
Linear Algebra (part 4): Eigenvalues, Diagonalization, and the Jordan Form (by Evan Dummit, 27, v ) Contents 4 Eigenvalues, Diagonalization, and the Jordan Canonical Form 4 Eigenvalues, Eigenvectors, and
More information2 BJORN JAWERTH AND WIM SWELDENS Our paper diers from these in that it contains some more recent developments and that it focuses on the \multiresolut
AN OVERVIEW OF WAVELET BASED MULTIRESOLUTION ANALYSES BJ ORN JAWERTH yz AND WIM SWELDENS yx Abstract. In this paper we present anoverview of wavelet based multiresolution analyses. First, we briey discuss
More informationLinear Algebra and Eigenproblems
Appendix A A Linear Algebra and Eigenproblems A working knowledge of linear algebra is key to understanding many of the issues raised in this work. In particular, many of the discussions of the details
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 information1 Matrices and Systems of Linear Equations
Linear Algebra (part ) : Matrices and Systems of Linear Equations (by Evan Dummit, 207, v 260) Contents Matrices and Systems of Linear Equations Systems of Linear Equations Elimination, Matrix Formulation
More informationAbsolutely indecomposable symmetric matrices
Journal of Pure and Applied Algebra 174 (2002) 83 93 wwwelseviercom/locate/jpaa Absolutely indecomposable symmetric matrices Hans A Keller a; ;1, A Herminia Ochsenius b;1 a Hochschule Technik+Architektur
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 informationON TRIVIAL GRADIENT YOUNG MEASURES BAISHENG YAN Abstract. We give a condition on a closed set K of real nm matrices which ensures that any W 1 p -grad
ON TRIVIAL GRAIENT YOUNG MEASURES BAISHENG YAN Abstract. We give a condition on a closed set K of real nm matrices which ensures that any W 1 p -gradient Young measure supported on K must be trivial the
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 informationt x 0.25
Journal of ELECTRICAL ENGINEERING, VOL. 52, NO. /s, 2, 48{52 COMPARISON OF BROYDEN AND NEWTON METHODS FOR SOLVING NONLINEAR PARABOLIC EQUATIONS Ivan Cimrak If the time discretization of a nonlinear parabolic
More informationNumerical Methods I: Eigenvalues and eigenvectors
1/25 Numerical Methods I: Eigenvalues and eigenvectors Georg Stadler Courant Institute, NYU stadler@cims.nyu.edu November 2, 2017 Overview 2/25 Conditioning Eigenvalues and eigenvectors How hard are they
More informationReduced Synchronization Overhead on. December 3, Abstract. The standard formulation of the conjugate gradient algorithm involves
Lapack Working Note 56 Conjugate Gradient Algorithms with Reduced Synchronization Overhead on Distributed Memory Multiprocessors E. F. D'Azevedo y, V.L. Eijkhout z, C. H. Romine y December 3, 1999 Abstract
More information/00 $ $.25 per page
Contemporary Mathematics Volume 00, 0000 Domain Decomposition For Linear And Nonlinear Elliptic Problems Via Function Or Space Decomposition UE-CHENG TAI Abstract. In this article, we use a function decomposition
More informationToeplitz-circulant Preconditioners for Toeplitz Systems and Their Applications to Queueing Networks with Batch Arrivals Raymond H. Chan Wai-Ki Ching y
Toeplitz-circulant Preconditioners for Toeplitz Systems and Their Applications to Queueing Networks with Batch Arrivals Raymond H. Chan Wai-Ki Ching y November 4, 994 Abstract The preconditioned conjugate
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 informationExponentials of Symmetric Matrices through Tridiagonal Reductions
Exponentials of Symmetric Matrices through Tridiagonal Reductions Ya Yan Lu Department of Mathematics City University of Hong Kong Kowloon, Hong Kong Abstract A simple and efficient numerical algorithm
More informationNew constructions of domain decomposition methods for systems of PDEs
New constructions of domain decomposition methods for systems of PDEs Nouvelles constructions de méthodes de décomposition de domaine pour des systèmes d équations aux dérivées partielles V. Dolean?? F.
More information= a. a = Let us now study what is a? c ( a A a )
7636S ADVANCED QUANTUM MECHANICS Solutions 1 Spring 010 1 Warm up a Show that the eigenvalues of a Hermitian operator A are real and that the eigenkets of A corresponding to dierent eigenvalues are orthogonal
More informationConjugate Gradient (CG) Method
Conjugate Gradient (CG) Method by K. Ozawa 1 Introduction In the series of this lecture, I will introduce the conjugate gradient method, which solves efficiently large scale sparse linear simultaneous
More information1. Introduction Let the least value of an objective function F (x), x2r n, be required, where F (x) can be calculated for any vector of variables x2r
DAMTP 2002/NA08 Least Frobenius norm updating of quadratic models that satisfy interpolation conditions 1 M.J.D. Powell Abstract: Quadratic models of objective functions are highly useful in many optimization
More informationCOMP 558 lecture 18 Nov. 15, 2010
Least squares We have seen several least squares problems thus far, and we will see more in the upcoming lectures. For this reason it is good to have a more general picture of these problems and how to
More informationContents. 2.1 Vectors in R n. Linear Algebra (part 2) : Vector Spaces (by Evan Dummit, 2017, v. 2.50) 2 Vector Spaces
Linear Algebra (part 2) : Vector Spaces (by Evan Dummit, 2017, v 250) Contents 2 Vector Spaces 1 21 Vectors in R n 1 22 The Formal Denition of a Vector Space 4 23 Subspaces 6 24 Linear Combinations and
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 informationGaussian Elimination without/with Pivoting and Cholesky Decomposition
Gaussian Elimination without/with Pivoting and Cholesky Decomposition Gaussian Elimination WITHOUT pivoting Notation: For a matrix A R n n we define for k {,,n} the leading principal submatrix a a k A
More informationThe average dimension of the hull of cyclic codes
Discrete Applied Mathematics 128 (2003) 275 292 www.elsevier.com/locate/dam The average dimension of the hull of cyclic codes Gintaras Skersys Matematikos ir Informatikos Fakultetas, Vilniaus Universitetas,
More informationS.F. Xu (Department of Mathematics, Peking University, Beijing)
Journal of Computational Mathematics, Vol.14, No.1, 1996, 23 31. A SMALLEST SINGULAR VALUE METHOD FOR SOLVING INVERSE EIGENVALUE PROBLEMS 1) S.F. Xu (Department of Mathematics, Peking University, Beijing)
More informationL. Levaggi A. Tabacco WAVELETS ON THE INTERVAL AND RELATED TOPICS
Rend. Sem. Mat. Univ. Pol. Torino Vol. 57, 1999) L. Levaggi A. Tabacco WAVELETS ON THE INTERVAL AND RELATED TOPICS Abstract. We use an abstract framework to obtain a multilevel decomposition of a variety
More informationSemidefinite and Second Order Cone Programming Seminar Fall 2001 Lecture 2
Semidefinite and Second Order Cone Programming Seminar Fall 2001 Lecture 2 Instructor: Farid Alizadeh Scribe: Xuan Li 9/17/2001 1 Overview We survey the basic notions of cones and cone-lp and give several
More information(z) UN-1(z) PN-1(z) Pn-1(z) - z k. 30(z) 30(z) (a) (b) (c) x y n ...
Minimum Memory Implementations of the Lifting Scheme Christos Chrysas Antonio Ortega HewlettPackard Laboratories Integrated Media Systems Center 151 Page Mill Road, Bldg.3U3 University of Southern California
More informationCompression on the digital unit sphere
16th Conference on Applied Mathematics, Univ. of Central Oklahoma, Electronic Journal of Differential Equations, Conf. 07, 001, pp. 1 4. ISSN: 107-6691. URL: http://ejde.math.swt.edu or http://ejde.math.unt.edu
More informationBilinear Stochastic Elliptic Equations
Bilinear Stochastic Elliptic Equations S. V. Lototsky and B. L. Rozovskii Contents 1. Introduction (209. 2. Weighted Chaos Spaces (210. 3. Abstract Elliptic Equations (212. 4. Elliptic SPDEs of the Full
More information1. The Polar Decomposition
A PERSONAL INTERVIEW WITH THE SINGULAR VALUE DECOMPOSITION MATAN GAVISH Part. Theory. The Polar Decomposition In what follows, F denotes either R or C. The vector space F n is an inner product space with
More informationLECTURE NOTES ELEMENTARY NUMERICAL METHODS. Eusebius Doedel
LECTURE NOTES on ELEMENTARY NUMERICAL METHODS Eusebius Doedel TABLE OF CONTENTS Vector and Matrix Norms 1 Banach Lemma 20 The Numerical Solution of Linear Systems 25 Gauss Elimination 25 Operation Count
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 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 informationRobust linear optimization under general norms
Operations Research Letters 3 (004) 50 56 Operations Research Letters www.elsevier.com/locate/dsw Robust linear optimization under general norms Dimitris Bertsimas a; ;, Dessislava Pachamanova b, Melvyn
More informationlinearly indepedent eigenvectors as the multiplicity of the root, but in general there may be no more than one. For further discussion, assume matrice
3. Eigenvalues and Eigenvectors, Spectral Representation 3.. Eigenvalues and Eigenvectors A vector ' is eigenvector of a matrix K, if K' is parallel to ' and ' 6, i.e., K' k' k is the eigenvalue. If is
More information