Computation of interpolatory splines via triadic subdivision
|
|
- Christiana Perkins
- 5 years ago
- Views:
Transcription
1 Computation of interpolatory splines via triadic subdivision Valery A. Zheludev and Amir Z. Averbuch Abstract We present an algorithm for computation of interpolatory splines of arbitrary order at triadic rational points. The algorithm is based on triadic subdivision of splines. Explicit expressions for the subdivision symbols are established. These are rational functions. The computations are implemented by recursive filtering. 1 Introduction Denote by S p the space of polynomial splines S p x of order p defined on the uniform grid g = {k}, k Z, such that the arrays {S p k}, k Z, belong to l 1. We propose to compute values of the splines via a triadic subdivision. The insertion rule for a spline Sx S p is Triadic Insertion Rule: Let f { = f k = S p k }, k Z. For j =, 1,..., we construct on the grid g j = { k3 j} }, k Z, a spline S p j x, which interpolates the sequence fj = {f j k on the grid g j. Then, = S p j k3 j+1, k Z. f j+1 k Note that the value of a spline at any point can be expressed as a linear combination of its values at grid points. In other words, any value f j+1 k can be derived by some filtering of the sequence f j. We present explicit expressions for these filters for splines of arbitrary order. Their transfer functions are rational functions. Computations are implemented by recursive filtering. Moreover, we prove that for any j N, f j k = S k3 j, k Z. Thus, we obtain a fast algorithm that computes the values of a spline S p x from the space Sp, which interpolates the sequence f on the grid g, at the triadic rational points { k3 j}, k Z. 2 Spline filters 2.1 B-splines The centered B-spline of first order is the characteristic function of the interval [ 1/2, 1/2]. The centered B-spline of order p can be expressed as the convolution M p x = M p 1 x M 1 x, p 2. Note that the B-spline of order p is supported on the interval p/2, p/2. It is positive within its 1
2 support and symmetric about zero. The B-spline M p consists of pieces of polynomials of degree p 1 that are linked to each other at the nodes, such that M p C p 2. Nodes of B-splines of even order are located at points {k} and of odd order at points {k + 1/2}, k Z. The Fourier transform of the B-spline of order p is M p ω = sin ω/2 p e iωx M p x dx =. 2.1 ω/2 The time domain representation of the B-spline is M p x = 1 p 1 p 1 k p 1! k k= x + p p 1 2 k, x + = x + x /2. + Shifts of B-splines form a basis in the space S p. Namely, any spline S p x S p has the following representation: Equation 2.1 implies that for any p 2 M p x k = 1 2π = e 2πiωk m p xω dω, S p x = e iωx k sin ω/2 ω/2 q k M p x k. p dω = e 2πilx l Z where m p xω = e 2πiωx sin πω p l Z e 2πiωx k sin πωp πl + ω p e 2πilx dω πl + ω p. 2.2 The relation 2.2 means that M p x k is a Fourier coefficient of the 1-periodic function m p xω and this function can be represented as the sum: m p xω = e2πikω M p x k. The functions m p xω are 1-periodic with respect to the variable ω and are splines from S p with respect to x. Their shifts along the grid are m p x+l ω = e 2πikω M p x k l = e 2πilω e 2πikω M p x k = e 2πilω m p xω. Following Schoenberg [4], we call the splines m p xω the exponential splines. 2.2 Interpolatory splines Equation 2.2 implies that a spline S p x S p can be represented via the exponential splines: S p x = q k M p x k = q k e 2πiωk m p xω dω = ˆqωm p xω dω, 2.3 where ˆqω = e 2πiωk q k is the discrete Fourier transform DFT of the coefficient sequence q = {q k }, k Z. Integral representation 2.3 was established in [5]. Denote u p ω = m p ω = e 2πikω M p k
3 Due to the compact support and symmetry of the B-splines, the functions u p ω are cosine polynomials. They were extensively studied in [1]. In particular, [1] established that u p ω > for all ω R and p N. In addition, recurrent relations for computation of the polynomials of any order were presented in [1]. However, u p ω can be easily calculated directly from Eq Equation 2.2 provides the following representation of the functions u p ω : u p ω = sin πω p l Z πl + ω p. The values of the spline S p x at the grid points and their DFT are S p n = ˆqωm p nω dω = e 2πiωnˆqωu p ω dω Ŝp ω = ˆqωu p ω. Assume that the spline S p x interpolates a sequence f at the grid points S p n = f n f ω = Ŝp ω ˆqω = f ω u p ω. 2.5 According to Triadic Insertion Rule, which was formulated in Section 1, f 1 n = S p n/3. Then, f 1 ω = n Z e 2πiωn S p n/3 = e 2πiω σ 1 3ω + σ 3ω + e 2πiω σ 1 3ω, 2.6 σ ω = n Z e 2πiωn S p n = ˆqωu p ω = f ω, σ ±1 ω = n Z e 2πiωn S p n ± 1. 3 Due to Eq. 2.2, S p n ± 1 = 3 ˆqωm p n±1/3 ω dω = v p ±1 ω = m p ±1/3 ω = e±2πiω/3 sin πω p l Z e 2πiωnˆqωv p 1 ω dω, e ±2πil/3 πl + ω p. 2.7 Hence, using Eq. 2.5, we obtain σ ±1 ω = vp ±1 ω f u p ω ω Triadic subdivision 3.1 Structure of filters Substituting Eq. 2.8 into Eq. 2.6, we obtain f 1 ω = T p ω f 3ω, T p ω = e 2πiω vp 1 3ω u p 3ω e 2πiω vp 1 3ω u p 3ω. 3.9 Similarly, f j+1 ω = T p ω f j 3ω, j N. Thus, in order to produce the array f j+1, we need to filter the upsampled array f j with the filter T p, whose DFT the frequency response T p ω is given in Eq We establish now the structure of the filter T p. 3
4 Theorem 3.1 The DFT of the filter T p, which originates from the spline of order p, is T p ω = cos πωp u p ω 3 p 1 u p ω Proof: We have from Eqs. 2.7 and 3.9 T p ω = G p ω/u p 3ω, where G p ω = e 2πiω v p 1 3ω + up 3ω + e 2πiω v p 1 3ω = sin 3πωp l Z = cos πω p sin πω p l Z cos 2πl/3. πl + 3ω p 1 + e 2πil/3 + e 2πil/3 πl + 3ω p We split the sum in the last equation into three sub-sums along l = 3m, l = 3m + 1 and l = 3m 1. Respectively, the function G p ω splits as G p ω = cos πω p γ ω + γ 1 ω + γ 1 ω, where γ ω = γ ±1 ω = Hence, Eq. 3.1 follows. sin πωp 3 p m Z sin πωp 3 p m Z 1 mp πm + ω p cos 2πm = 3 p+1 u p ω, 1 m+1p cos 2πm ± 1/3 =. πm + ω ± π/3 p 3.2 Convergence The spline L p x S p, which interpolates the Kronecker delta sequence δk, is called the fundamental spline of the space S p. The spline S p x S p, which interpolates the data sequence f = { f k }, can be represented as S p x = f k Lp x k L p x decays exponentially when x grows. Therefore, Eq represents the spline, which interpolates the sequence { f k }, k Z, even in the case when the sequence has a power growth. In this case, the spline S p x also has a power growth. It was proved in [2] that the spline of power growth, interpolating the data of power growth, is unique. We prove that the values of the spline S p x at any set of triadic rational points S p k3 j can be calculated by successive application of Triadic Insertion Rule, which was formulated in Section 1, to the array f. A similar fact holds for Dyadic Insertion Rule when order p of splines is even. But this is not true for Dyadic Insertion Rule when p is odd. In this case, when p is odd, the spline-based subdivision no longer converges to a spline but rather to a function, which is smoother than the generating spline see [6]. 4
5 Theorem 3.2 Assume the spline S p x belongs to Sp, S p k = f k and f = { fk } l1. Let T p be { } the filter defined by Eq Assume that for j = 1, 2,..., f j = f j k is the array, whose DFT is found from the relation f j+1 ω = T p ω f j 3ω. Then, S p k3 j = f j k, j = 1, 2,..., Proof: We have S p k3 1 = f 1 k from the definition of Triadic Insertion Rule. Next step of subdivision consists in the construction of the spline S p 1 x on the grid g1 = { k3 1}, such that S p 1 k3 1 = f 1 k. Then, fk 2 = Sp 1 k3 2. We prove that fk 2 = Sp k/9. The subsequent relations in Eq are derived by a simple induction. The array f 2 is obtained by repeated application of the filter T p to the array f : f 2 ω = T p ω f 1 3ω = T p ω T p 3ω f 9ω = cos πωp u p ω cos 3πω p u p 3ω 3 p 1 u p 3ω 3 p 1 u p 9ω f 9ω = Hω f u p 9ω 9ω, where Hω = 9 1 p cos πω1 + 2 cos 3πω p u p ω. Denote s k = S p k/9, k Z. Then, the DFT of this array is ŝω = H f 9ωω/u p 9ω, where Hω = 4 e 2νπiω w ν 9ω + e 2νπiω w ν 9ω + u p 9ω, w ±ν ω = m p ±ν/9 ω. ν=1 We can write, using Eq. 2.7, where Hω = sin 9πω p l Z πl + 9ω p e 2νπil/9 + e 2νπil/9 ν=1 = cos 3πω1 + 2 cos πω p sin πω p l Z = cos 3πω1 + 2 cos πω p 4 κ µ ω = µ= 4 4 cos 2πνµ/9 ν=1 κ µ ω, sin πω p 9 p πl + 9ω p m Z 4 cos 2πνl/9 ν=1 1 m+1p πm + ω + µπ/9 p. It is easily verified that for µ = ±1,..., ±4 κ µ ω, while κ ω = 9 1 p u p ω. Hence, it follows that Hω = Hω f 2 k = Sp k3 2, k Z. Repeating our reasoning with the initial data set f 1 instead of f, we prove that f 3 k = Sp 1 k3 3 and so on. 5
6 Remark If the initial data is the delta sequence fk = δ k then we get f j k = Lp k3 l, where L p x is the fundamental spline of the space S p. Therefore, we can extend the assertion of Theorem 3.2 from the splines belonging to S p to the splines that interpolate sequences of power growth. 4 Remarks on the computational aspects of the spline subdivision It was stated in Section 3.1 that in order to produce the array f j+1 we need to filter the upsampled array f j with the filter T p. In the frequency domain, filtering is presented by Eq The frequency response of the filter is given by Eq For computational purposes, we present filtering using the z-transform. We recall that the the z-transform of a sequence a = {a k } l 1 is Az = z k a k. The z-transform of the shifted sequence is z k a k+l = z l Az. Denote z = e 2πiω. Then, the DFT of a sequence a becomes its z-transformâω = Az. The DFT of a filter T becomes its z-transform T ω = T z and it is called the transfer function. The subdivision is represented as F j+1 z = T p zf j z 3, T p z = z z p U p z U p z 3, U p z = z k M p k The function U p z is a Laurent polynomial and the transfer function T p z is a rational function. Such filters are called the infinite impulse response IIR filters as opposite to the finite impulse response FIR filters, whose transfer functions are Laurent polynomials. Proposition 4.1 [2] The roots of the Laurent polynomials U p z are all simple and negative. Each root ζ can be paired with a dual root θ such that ζ θ = 1. Thus, if p = 2r 1, p = 2r then U p z can be represented by: U p z = r n=1 1 γ n 1 + γ n z1 + γ n z 1, < γ 1 < γ 2 <... γ r = e g < 1, g >. The above properties of the polynomials U p z enable to implement filtering via a recursive algorithm, which is widely used in signal processing. There are two ways to implement the described subdivision scheme. Polyphase filtering: One way is to implement the filter T p = {T k }, k Z, using the so-called polyphase representation of the filter: T p z = zt p 1 z3 + T p z3 + z 1 T p 1 z3, Tν p z = z k T p 3k+ν, ν = 1,, 1. Then, the polyphase representation of the array f j+1 is F j+1 z = zf j+1 1 z3 + F j+1 z 3 + z 1 F j+1 1 z 3, Fν j+1 z = z k f j+1 3k+ν = T ν p zf j z, ν = 1,, 1. 6
7 { } { Thus, in order to retrieve the sub-arrays f j+1 3k±1 we have to filter the array { whose transfer functions are T p ±1 z, respectively. Recall that f j+1 3k = f j k } f j k with the filters, }. The polyphase components of the transfer function T p z can be derived from the expression for T p z given in Eq Obviously, T p 1 z = T p 1 z 1. However, it is easily seen from Eq. 3.9 that T p ±1 e2πiω = v p ±1 ω/up ω. Direct filtering: Equation 4.13 suggests that, when several steps of the subdivision must be carried out, a direct application of the filter T p is preferable. It follows from Eq that j 1 F j+1 z = T p z F j z 3 = For example, l=1 T p z 3l F z 3j = U p z j 1 l= z 3l z 3l p U p F z 3j. z 3j F 3 z = U p zz z p z z 3 p z z 9 p U p z 27 F z 27. Thus, the subdivision is implemented via the following steps: 1. The IIR filter with the transfer function 1/U p z is applied to the data array f. 2. The produced array is upsampled 1 and filtered with FIR filter, whose transfer function is z z p repeated j times. 3. The produced array is filtered with FIR filter whose transfer function is U p z. Note that the IIR filtering is applied only once. Examples of filters Linear spline: T 2 z = z z 1 2 /3, T 2 1 z = T 1 2 z 1 = z + 2/3. This is a single FIR filter in the family. All the filters derived from splines of higher orders are IIR. Quadratic spline: T 3 z = z z 1 z z z z 3, T1 3 z = T 3 1z 1 = 25z z 1 9 z z Cubic spline: T 4 z = z z 1 z z z z 3, T1 4 z = T 1z 4 1 = z2 + 6z z 1 27 z z 1. 1 Upsampling means replacing an array {a k } by the array {ã k } such that ã 3k = a k and ã 3k±1 =. 7
8 Spline of fourth degree : z T z z 1 + z 2 z z 1 5 z = 81 z z z 3 + z 6, T1 5 z = T 1z 5 1 = 625z z z 1 + z 2 81 z z z 1 + z 2. Spline of fifth degree : z T z z 1 + z 2 z z 1 6 z = 243 z z z 3 + z 6, T1 6 z = T 1z 6 1 = z z z z z z z z 1 + z 2. We now provide an example of a recursive implementation of an IIR filter. In [6] this is discussed in more details. Implementation of the filter T 3 1. The z-transform of this filter is given in Eq. 4.14: T 3 1 z = 25z z 1 9 z z 1 = α 25z z αz1 + αz 1, where α = Then, application of the filter T 3 1 to an array a = {a k }, whose z-transform is Az, is implemented as a subsequent application of three filters: b = T 3 1a = R l R r F a, which are are defined by their z-transforms: F z = α 9 25z z 1, R r z = Thus, filtering is carried out in three steps: αz 1, R lz = αz. A 1 z = F zaz a 1 k = α 9 25a k a k + a k 1, A 2 z = R r za 1 z 1 + αz 1 A 2 z = A 1 z a 2 k = a1 k αa2 k 1, Bz = R l za 2 z 1 + αz Bz = A 2 z b k = a 2 k αb k+1. The filter F is FIR unlike the filters R l and R r. Application of the filter R r is called causal recursive filtering. Here, for the calculation of the term a 2 k, the previously derived term a2 k 1 is used. Application of the R l is called anti-causal recursive filtering. All these procedures are implemented in a fast way. Computation of splines of higher orders uses filters, which are factorized into longer cascades of the same structure. References [1] I. J. Schoenberg, Contribution to the problem of approximation of equidistant data by analytic functions, Quart.Appl. Math., , 45-99,
9 [2] I. J. Schoenberg, Cardinal interpolation and spline functions, J. Approx. Th., 2, 1969, [3] I. J. Schoenberg, Cardinal interpolation and spline functions II, J. Approx. Theory, , [4] I. J. Schoenberg, Cardinal spline interpolation, CBMS, 12, SIAM, Philadelphia, [5] V. A. Zheludev, Integral representation of slowly growing equidistant splines, Approximation Theory and Applications, 14, no. 4, 1998, [6] V. A. Zheludev, Interpolatory subdivision schemes with infinite masks originated from splines, Advances in Comp. Math. 25,26,
Wavelet transforms generated by splines
Wavelet transforms generated by splines Amir Z. Averbuch Valery A. Zheludev School of Computer Science Tel Aviv University Tel Aviv 69978, Israel Abstract In this paper we design a new family of biorthogonal
More informationA C 2 Four-Point Subdivision Scheme with Fourth Order Accuracy and its Extensions
A C 2 Four-Point Subdivision Scheme with Fourth Order Accuracy and its Extensions Nira Dyn School of Mathematical Sciences Tel Aviv University Michael S. Floater Department of Informatics University of
More informationA C 2 Four-Point Subdivision Scheme with Fourth Order Accuracy and its Extensions
A C 2 Four-Point Subdivision Scheme with Fourth Order Accuracy and its Extensions Nira Dyn Michael S. Floater Kai Hormann Abstract. We present a new four-point subdivision scheme that generates C 2 curves.
More informationA NEW PROOF OF THE SMOOTHNESS OF 4-POINT DESLAURIERS-DUBUC SCHEME
J. Appl. Math. & Computing Vol. 18(2005), No. 1-2, pp. 553-562 A NEW PROOF OF THE SMOOTHNESS OF 4-POINT DESLAURIERS-DUBUC SCHEME YOUCHUN TANG, KWAN PYO KO* AND BYUNG-GOOK LEE Abstract. It is well-known
More informationRECENTLY frames or redundant expansions of signals
2126 IEEE TRANSACTIONS ON SIGNAL PROCESSING, VOL. 54, NO. 6, JUNE 2006 Interpolatory Frames in Signal Space Amir Z. Averbuch, Valery A. Zheludev, Tamir Cohen Abstract We present a new family of frames,
More informationInterpolation with complex B-splines
Interpolation with complex B-splines Brigitte Forster Fakultät für Informatik und Mathematik Universität Passau Brigitte Forster (Passau) Bernried 2016 1 / 32 The first part of a story in three chapters
More informationThe Parks McClellan algorithm: a scalable approach for designing FIR filters
1 / 33 The Parks McClellan algorithm: a scalable approach for designing FIR filters Silviu Filip under the supervision of N. Brisebarre and G. Hanrot (AriC, LIP, ENS Lyon) PEQUAN Seminar, February 26,
More informationApplications of Wavelets and Framelets
Applications of Wavelets and Framelets Bin Han Department of Mathematical and Statistical Sciences University of Alberta, Edmonton, Canada Present at 2017 International Undergraduate Summer Enrichment
More informationInterpolation and Splines
Interpolation and Splines Anna Gryboś October 23, 27 1 Problem setting Many of physical phenomenona are described by the functions that we don t know exactly. Often we can calculate or measure the values
More informationFebruary 2017 (1/20) 2 Piecewise Polynomial Interpolation 2.2 (Natural) Cubic Splines. MA378/531 Numerical Analysis II ( NA2 )
f f f f f (/2).9.8.7.6.5.4.3.2. S Knots.7.6.5.4.3.2. 5 5.2.8.6.4.2 S Knots.2 5 5.9.8.7.6.5.4.3.2..9.8.7.6.5.4.3.2. S Knots 5 5 S Knots 5 5 5 5.35.3.25.2.5..5 5 5.6.5.4.3.2. 5 5 4 x 3 3.5 3 2.5 2.5.5 5
More informationMar. 20 Math 2335 sec 001 Spring 2014
Mar. 20 Math 2335 sec 001 Spring 2014 Chebyshev Polynomials Definition: For an integer n 0 define the function ( ) T n (x) = cos n cos 1 (x), 1 x 1. It can be shown that T n is a polynomial of degree n.
More informationWAVELET transforms have become one of the basic tools
IEEE TRANSACTIONS ON IMAGE PROCESSING, VOL 13, NO 7, JULY 2004 993 A New Family of Spline-Based Biorthogonal Wavelet Transforms and Their Application to Image Compression Amir Z Averbuch and Valery A Zheludev
More informationParameterization of triangular meshes
Parameterization of triangular meshes Michael S. Floater November 10, 2009 Triangular meshes are often used to represent surfaces, at least initially, one reason being that meshes are relatively easy to
More informationFitting to a set of data. Lecture on fitting
Fitting to a set of data Lecture on fitting Linear regression Linear regression Residual is the amount difference between a real data point and a modeled data point Fitting a polynomial to data Could use
More informationRefinable bivariate quartic and quintic C 2 -splines for quadrilateral subdivisions
Refinable bivariate quartic and quintic C 2 -splines for quadrilateral subdivisions Charles K. Chui, Qingtang Jiang Department of Mathematics and Computer Science University of Missouri St. Louis St. Louis,
More informationAn interpolating 4-point C 2 ternary stationary subdivision scheme
Computer Aided Geometric Design 9 (2002) 8 www.elsevier.com/locate/comaid An interpolating 4-point C 2 ternary stationary subdivision scheme M.F Hassan a,, I.P. Ivrissimitzis a, N.A. Dodgson a,m.a.sabin
More informationDURING the past decade, there has been an increasing
IEEE TRANSACTIONS ON SIGNAL PROCESSING, VOL. 53, NO. 4, APRIL 2005 1425 Cardinal Exponential Splines: Part I Theory and Filtering Algorithms Michael Unser, Fellow, IEEE, and Thierry Blu, Member, IEEE Abstract
More informationDeficient Quartic Spline Interpolation
International Journal of Computational Science and Mathematics. ISSN 0974-3189 Volume 3, Number 2 (2011), pp. 227-236 International Research Publication House http://www.irphouse.com Deficient Quartic
More informationHarmonic Spline Series Representation of Scaling Functions
Harmonic Spline Series Representation of Scaling Functions Thierry Blu and Michael Unser Biomedical Imaging Group, STI/BIO-E, BM 4.34 Swiss Federal Institute of Technology, Lausanne CH-5 Lausanne-EPFL,
More informationRepresenting Curves Part II. Foley & Van Dam, Chapter 11
Representing Curves Part II Foley & Van Dam, Chapter 11 Representing Curves Polynomial Splines Bezier Curves Cardinal Splines Uniform, non rational B-Splines Drawing Curves Applications of Bezier splines
More informationPS Geometric Modeling Homework Assignment Sheet I (Due 20-Oct-2017)
Homework Assignment Sheet I (Due 20-Oct-2017) Assignment 1 Let n N and A be a finite set of cardinality n = A. By definition, a permutation of A is a bijective function from A to A. Prove that there exist
More informationRegularity Analysis of Non Uniform Data
Regularity Analysis of Non Uniform Data Christine Potier and Christine Vercken Abstract. A particular class of wavelet, derivatives of B-splines, leads to fast and ecient algorithms for contours detection
More informationTHE GROWTH DEGREE OF VERTEX REPLACEMENT RULES
THE GROWTH DEGREE OF VERTEX REPLACEMENT RULES JOSEPH P. PREVITE AND MICHELLE PREVITE 1. Introduction. 2. Vertex Replacement Rules. 3. Marked Graphs Definition 3.1. A marked graph (G, p) is a graph with
More informationMatrix-valued 4-point Spline and 3-point Non-spline Interpolatory Curve Subdivision Schemes
Matrix-valued 4-point Spline and -point Non-spline Interpolatory Curve Subdivision Schemes Charles K. Chui, Qingtang Jiang Department of Mathematics and Computer Science University of Missouri St. Louis
More informationCS321 Introduction To Numerical Methods
CS3 Introduction To Numerical Methods Fuhua (Frank) Cheng Department of Computer Science University of Kentucky Lexington KY 456-46 - - Table of Contents Errors and Number Representations 3 Error Types
More information4 Generating functions in two variables
4 Generating functions in two variables (Wilf, sections.5.6 and 3.4 3.7) Definition. Let a(n, m) (n, m 0) be a function of two integer variables. The 2-variable generating function of a(n, m) is F (x,
More informationPoint Lattices in Computer Graphics and Visualization how signal processing may help computer graphics
Point Lattices in Computer Graphics and Visualization how signal processing may help computer graphics Dimitri Van De Ville Ecole Polytechnique Fédérale de Lausanne Biomedical Imaging Group dimitri.vandeville@epfl.ch
More informationSpline Solution Of Some Linear Boundary Value Problems
Applied Mathematics E-Notes, 8(2008, 171-178 c ISSN 1607-2510 Available free at mirror sites of http://www.math.nthu.edu.tw/ amen/ Spline Solution Of Some Linear Boundary Value Problems Abdellah Lamnii,
More informationA Study on Subdivision Scheme-Draft. Kwan Pyo Ko Dongseo University Busan, South Korea
A Study on Subdivision Scheme-Draft Kwan Pyo Ko Dongseo University Busan, South Korea April 30, 007 Contents Introduction 9 Subdivision of Univariate Data 3. Definitions and Basic Results.............................
More informationspecified or may be difficult to handle, we often have a tabulated data
Interpolation Introduction In many practical situations, for a function which either may not be explicitly specified or may be difficult to handle, we often have a tabulated data where and for In such
More informationLacunary Interpolation Using Quartic B-Spline
General Letters in Mathematic, Vol. 2, No. 3, June 2017, pp. 129-137 e-issn 2519-9277, p-issn 2519-9269 Available online at http:\\ www.refaad.com Lacunary Interpolation Using Quartic B-Spline 1 Karwan
More informationInterpolation & Polynomial Approximation. Cubic Spline Interpolation II
Interpolation & Polynomial Approximation Cubic Spline Interpolation II Numerical Analysis (9th Edition) R L Burden & J D Faires Beamer Presentation Slides prepared by John Carroll Dublin City University
More informationEdges, interpolation, templates. Nuno Vasconcelos ECE Department, UCSD (with thanks to David Forsyth)
Edges, interpolation, templates Nuno Vasconcelos ECE Department, UCSD (with thanks to David Forsyth) Edge detection edge detection has many applications in image processing an edge detector implements
More informationDiscrete Cubic Interpolatory Splines
Publ RIMS, Kyoto Univ. 28 (1992), 825-832 Discrete Cubic Interpolatory Splines By Manjulata SHRIVASTAVA* Abstract In the present paper, existence, uniqueness and convergence properties of a discrete cubic
More informationHighly Symmetric Bi-frames for Triangle Surface Multiresolution Processing
Highly Symmetric Bi-frames for Triangle Surface Multiresolution Processing Qingtang Jiang and Dale K. Pounds Abstract In this paper we investigate the construction of dyadic affine (wavelet) bi-frames
More information3. Lifting Scheme of Wavelet Transform
3. Lifting Scheme of Wavelet Transform 3. Introduction The Wim Sweldens 76 developed the lifting scheme for the construction of biorthogonal wavelets. The main feature of the lifting scheme is that all
More informationMATH 1A MIDTERM 1 (8 AM VERSION) SOLUTION. (Last edited October 18, 2013 at 5:06pm.) lim
MATH A MIDTERM (8 AM VERSION) SOLUTION (Last edited October 8, 03 at 5:06pm.) Problem. (i) State the Squeeze Theorem. (ii) Prove the Squeeze Theorem. (iii) Using a carefully justified application of the
More informationObjectives and Homework List
MAC 1140 Objectives and Homework List Each objective covered in MAC1140 is listed below. Along with each objective is the homework list used with MyMathLab (MML) and a list to use with the text (if you
More informationInterpolation by Spline Functions
Interpolation by Spline Functions Com S 477/577 Sep 0 007 High-degree polynomials tend to have large oscillations which are not the characteristics of the original data. To yield smooth interpolating curves
More informationGeneralised Mean Averaging Interpolation by Discrete Cubic Splines
Publ. RIMS, Kyoto Univ. 30 (1994), 89-95 Generalised Mean Averaging Interpolation by Discrete Cubic Splines By Manjulata SHRIVASTAVA* Abstract The aim of this work is to introduce for a discrete function,
More informationMath (Spring 2009): Lecture 5 Planes. Parametric equations of curves and lines
Math 18.02 (Spring 2009): Lecture 5 Planes. Parametric equations of curves and lines February 12 Reading Material: From Simmons: 17.1 and 17.2. Last time: Square Systems. Word problem. How many solutions?
More informationGeneralized barycentric coordinates
Generalized barycentric coordinates Michael S. Floater August 20, 2012 In this lecture, we review the definitions and properties of barycentric coordinates on triangles, and study generalizations to convex,
More informationDOWNLOAD PDF BIG IDEAS MATH VERTICAL SHRINK OF A PARABOLA
Chapter 1 : BioMath: Transformation of Graphs Use the results in part (a) to identify the vertex of the parabola. c. Find a vertical line on your graph paper so that when you fold the paper, the left portion
More informationLecture 2.2 Cubic Splines
Lecture. Cubic Splines Cubic Spline The equation for a single parametric cubic spline segment is given by 4 i t Bit t t t i (..) where t and t are the parameter values at the beginning and end of the segment.
More informationQuasi-Quartic Trigonometric Bézier Curves and Surfaces with Shape Parameters
Quasi-Quartic Trigonometric Bézier Curves and Surfaces with Shape Parameters Reenu Sharma Assistant Professor, Department of Mathematics, Mata Gujri Mahila Mahavidyalaya, Jabalpur, Madhya Pradesh, India
More informationFourier Transform of Bernstein-Bézier Polynomials
Fourier Transform of Bernstein-Bézier Polynomials Charles. K. Chui, Tian-Xiao He 2, and Qingtang Jiang 3 Department of Mathematics and Computer Science University of Missouri-St. Louis, St. Louis, MO 632,
More informationCS 450 Numerical Analysis. Chapter 7: Interpolation
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 informationMath 5593 Linear Programming Lecture Notes
Math 5593 Linear Programming Lecture Notes Unit II: Theory & Foundations (Convex Analysis) University of Colorado Denver, Fall 2013 Topics 1 Convex Sets 1 1.1 Basic Properties (Luenberger-Ye Appendix B.1).........................
More informationTHREE LECTURES ON BASIC TOPOLOGY. 1. Basic notions.
THREE LECTURES ON BASIC TOPOLOGY PHILIP FOTH 1. Basic notions. Let X be a set. To make a topological space out of X, one must specify a collection T of subsets of X, which are said to be open subsets of
More informationResearch Article A Family of Even-Point Ternary Approximating Schemes
International Scholarly Research Network ISRN Applied Mathematics Volume, Article ID 97, pages doi:.5//97 Research Article A Family of Even-Point Ternary Approximating Schemes Abdul Ghaffar and Ghulam
More informationGREENWOOD PUBLIC SCHOOL DISTRICT Algebra III Pacing Guide FIRST NINE WEEKS
GREENWOOD PUBLIC SCHOOL DISTRICT Algebra III FIRST NINE WEEKS Framework/ 1 Aug. 6 10 5 1 Sequences Express sequences and series using recursive and explicit formulas. 2 Aug. 13 17 5 1 Sequences Express
More informationGeneralized Lane Riesenfeld algorithms
Generalized Lane Riesenfeld algorithms Thomas J. Cashman a, Kai Hormann a, Ulrich Reif b a Università della Svizzera italiana, Via Giuseppe Buffi 13, CH-6904 Lugano b TU Darmstadt, Schlossgartenstrasse
More informationcorrelated to the Michigan High School Mathematics Content Expectations
correlated to the Michigan High School Mathematics Content Expectations McDougal Littell Algebra 1 Geometry Algebra 2 2007 correlated to the STRAND 1: QUANTITATIVE LITERACY AND LOGIC (L) STANDARD L1: REASONING
More informationThe Euler Equidimensional Equation ( 3.2)
The Euler Equidimensional Equation ( 3.) The Euler Equidimensional Equation ( 3.) The Euler Equidimensional Equation Definition The Euler equidimensional equation for the unknown function y with singular
More informationEdge and local feature detection - 2. Importance of edge detection in computer vision
Edge and local feature detection Gradient based edge detection Edge detection by function fitting Second derivative edge detectors Edge linking and the construction of the chain graph Edge and local feature
More informationCOMPUTER AIDED ENGINEERING DESIGN (BFF2612)
COMPUTER AIDED ENGINEERING DESIGN (BFF2612) BASIC MATHEMATICAL CONCEPTS IN CAED by Dr. Mohd Nizar Mhd Razali Faculty of Manufacturing Engineering mnizar@ump.edu.my COORDINATE SYSTEM y+ y+ z+ z+ x+ RIGHT
More informationConstruction of fractional spline wavelet bases
Construction of fractional spline wavelet bases Michael Unser and Thierry Blu Biomedical Imaging Group, Swiss Federal Institute of Technology Lausanne DMT/IOA, BM 4.7 CH-5 Lausanne EPFL, Switzerland Email:
More informationused to describe all aspects of imaging process: input scene f, imaging system O, and output image g.
SIMG-716 Linear Imaging Mathematics I 03 1 Representations of Functions 1.1 Functions Function: rule that assigns one numerical value (dependent variable) to each of a set of other numerical values (independent
More information7. The Gauss-Bonnet theorem
7. The Gauss-Bonnet theorem 7.1 Hyperbolic polygons In Euclidean geometry, an n-sided polygon is a subset of the Euclidean plane bounded by n straight lines. Thus the edges of a Euclidean polygon are formed
More informationStationary and Non-Stationary Univariate Subdivision Schemes
Punjab University Journal of Mathematics (ISSN 6-56) Vol. 5(3)(8) pp. 5- Stationary and Non-Stationary Univariate Subdivision Schemes Muhammad Asghar Department of Mathematics, The Islamia University of
More informationSampling and Reconstruction
Page 1 Sampling and Reconstruction The sampling and reconstruction process Real world: continuous Digital world: discrete Basic signal processing Fourier transforms The convolution theorem The sampling
More informationTopology 550A Homework 3, Week 3 (Corrections: February 22, 2012)
Topology 550A Homework 3, Week 3 (Corrections: February 22, 2012) Michael Tagare De Guzman January 31, 2012 4A. The Sorgenfrey Line The following material concerns the Sorgenfrey line, E, introduced in
More informationused to describe all aspects of imaging process: input scene f, imaging system O, and output image g.
SIMG-716 Linear Imaging Mathematics I 03 1 Representations of Functions 1.1 Functions Function: rule that assigns one numerical value (dependent variable) to each of a set of other numerical values (independent
More informationMastery. PRECALCULUS Student Learning Targets
PRECALCULUS Student Learning Targets Big Idea: Sequences and Series 1. I can describe a sequence as a function where the domain is the set of natural numbers. Connections (Pictures, Vocabulary, Definitions,
More informationDigital Signal Processing Lecture Notes 22 November 2010
Digital Signal Processing Lecture otes 22 ovember 2 Topics: Discrete Cosine Transform FFT Linear and Circular Convolution Rate Conversion Includes review of Fourier transforms, properties of Fourier transforms,
More informationPoint-Set Topology II
Point-Set Topology II Charles Staats September 14, 2010 1 More on Quotients Universal Property of Quotients. Let X be a topological space with equivalence relation. Suppose that f : X Y is continuous and
More informationConvergence of C 2 Deficient Quartic Spline Interpolation
Advances in Computational Sciences and Technology ISSN 0973-6107 Volume 10, Number 4 (2017) pp. 519-527 Research India Publications http://www.ripublication.com Convergence of C 2 Deficient Quartic Spline
More informationFourier analysis and sampling theory
Reading Required: Shirley, Ch. 9 Recommended: Fourier analysis and sampling theory Ron Bracewell, The Fourier Transform and Its Applications, McGraw-Hill. Don P. Mitchell and Arun N. Netravali, Reconstruction
More information1. How many white tiles will be in Design 5 of the pattern? Explain your reasoning.
Algebra 2 Semester 1 Review Answer the question for each pattern. 1. How many white tiles will be in Design 5 of the pattern Explain your reasoning. 2. What is another way to represent the expression 3.
More informationNumerical Integration
Numerical Integration Numerical Integration is the process of computing the value of a definite integral, when the values of the integrand function, are given at some tabular points. As in the case of
More informationMATH 54 - LECTURE 10
MATH 54 - LECTURE 10 DAN CRYTSER The Universal Mapping Property First we note that each of the projection mappings π i : j X j X i is continuous when i X i is given the product topology (also if the product
More informationLecture 25: Bezier Subdivision. And he took unto him all these, and divided them in the midst, and laid each piece one against another: Genesis 15:10
Lecture 25: Bezier Subdivision And he took unto him all these, and divided them in the midst, and laid each piece one against another: Genesis 15:10 1. Divide and Conquer If we are going to build useful
More information4. Definition: topological space, open set, topology, trivial topology, discrete topology.
Topology Summary Note to the reader. If a statement is marked with [Not proved in the lecture], then the statement was stated but not proved in the lecture. Of course, you don t need to know the proof.
More informationCombinatorics of free product graphs
Combinatorics of free product graphs Gregory Quenell March 8, 994 Abstract We define the return generating function on an abstract graph, and develop tools for computing such functions. The relation between
More informationDigital Signal Processing with Field Programmable Gate Arrays
Uwe Meyer-Baese Digital Signal Processing with Field Programmable Gate Arrays Third Edition With 359 Figures and 98 Tables Book with CD-ROM ei Springer Contents Preface Preface to Second Edition Preface
More informationNon-Holonomic Sequences and Functions
Non-Holonomic Sequences and Functions Stefan Gerhold May 16, 2006 Holonomic Sequences and Functions A function f (z) is holonomic if it satisfies an LODE p 0 (z)f (0) (z) + + p d (z)f (d) (z) = 0 with
More informationTreewidth and graph minors
Treewidth and graph minors Lectures 9 and 10, December 29, 2011, January 5, 2012 We shall touch upon the theory of Graph Minors by Robertson and Seymour. This theory gives a very general condition under
More informationData Analysis 3. Support Vector Machines. Jan Platoš October 30, 2017
Data Analysis 3 Support Vector Machines Jan Platoš October 30, 2017 Department of Computer Science Faculty of Electrical Engineering and Computer Science VŠB - Technical University of Ostrava Table of
More informationThe goal is the definition of points with numbers and primitives with equations or functions. The definition of points with numbers requires a
The goal is the definition of points with numbers and primitives with equations or functions. The definition of points with numbers requires a coordinate system and then the measuring of the point with
More informationTOPOLOGY, DR. BLOCK, FALL 2015, NOTES, PART 3.
TOPOLOGY, DR. BLOCK, FALL 2015, NOTES, PART 3. 301. Definition. Let m be a positive integer, and let X be a set. An m-tuple of elements of X is a function x : {1,..., m} X. We sometimes use x i instead
More informationPantographic polygons
203 Pantographic polygons John Miller and Emanuel Strzelecki Abstract Necessary and sufficient conditions are given for a polygon to be pantographic. The property is held by all regular polygons and by
More informationREGULAR GRAPHS OF GIVEN GIRTH. Contents
REGULAR GRAPHS OF GIVEN GIRTH BROOKE ULLERY Contents 1. Introduction This paper gives an introduction to the area of graph theory dealing with properties of regular graphs of given girth. A large portion
More informationLecture 1. 1 Notation
Lecture 1 (The material on mathematical logic is covered in the textbook starting with Chapter 5; however, for the first few lectures, I will be providing some required background topics and will not be
More informationBernstein-Bezier Splines on the Unit Sphere. Victoria Baramidze. Department of Mathematics. Western Illinois University
Bernstein-Bezier Splines on the Unit Sphere Victoria Baramidze Department of Mathematics Western Illinois University ABSTRACT I will introduce scattered data fitting problems on the sphere and discuss
More informationSubdivision of Curves and Surfaces: An Overview
Subdivision of Curves and Surfaces: An Overview Ben Herbst, Karin M Hunter, Emile Rossouw Applied Mathematics, Department of Mathematical Sciences, University of Stellenbosch, Private Bag X1, Matieland,
More informationCubic spline interpolation
Cubic spline interpolation In the following, we want to derive the collocation matrix for cubic spline interpolation. Let us assume that we have equidistant knots. To fulfill the Schoenberg-Whitney condition
More informationGenerell Topologi. Richard Williamson. May 6, 2013
Generell Topologi Richard Williamson May 6, 2013 1 8 Thursday 7th February 8.1 Using connectedness to distinguish between topological spaces I Proposition 8.1. Let (, O ) and (Y, O Y ) be topological spaces.
More informationFebruary 23 Math 2335 sec 51 Spring 2016
February 23 Math 2335 sec 51 Spring 2016 Section 4.1: Polynomial Interpolation Interpolation is the process of finding a curve or evaluating a function whose curve passes through a known set of points.
More informationOrdinary differential equations solving methods
Radim Hošek, A07237 radhost@students.zcu.cz Ordinary differential equations solving methods Problem: y = y2 (1) y = x y (2) y = sin ( + y 2 ) (3) Where it is possible we try to solve the equations analytically,
More informationNon-Uniform Recursive Doo-Sabin Surfaces (NURDSes)
Non-Uniform Recursive Doo-Sabin Surfaces Zhangjin Huang 1 Guoping Wang 2 1 University of Science and Technology of China 2 Peking University, China SIAM Conference on Geometric and Physical Modeling Doo-Sabin
More informationPLANAR GRAPH BIPARTIZATION IN LINEAR TIME
PLANAR GRAPH BIPARTIZATION IN LINEAR TIME SAMUEL FIORINI, NADIA HARDY, BRUCE REED, AND ADRIAN VETTA Abstract. For each constant k, we present a linear time algorithm that, given a planar graph G, either
More informationFourier transform of Bernstein Bézier polynomials
Illinois Wesleyan University From the SelectedWorks of Tian-Xiao He 27 Fourier transform of Bernstein Bézier polynomials Tian-Xiao He, Illinois Wesleyan University Charles K. Chui Qingtang Jiang Available
More informationC 1 Quintic Spline Interpolation Over Tetrahedral Partitions
C 1 Quintic Spline Interpolation Over Tetrahedral Partitions Gerard Awanou and Ming-Jun Lai Abstract. We discuss the implementation of a C 1 quintic superspline method for interpolating scattered data
More information1.1 Pearson Modeling and Equation Solving
Date:. Pearson Modeling and Equation Solving Syllabus Objective:. The student will solve problems using the algebra of functions. Modeling a Function: Numerical (data table) Algebraic (equation) Graphical
More informationFourier transform of images
Fourier transform of images Stefano Ferrari Università degli Studi di Milano stefano.ferrari@unimi.it Methods for Image Processing academic year 2014 2015 Extension to bidimensional domain The concepts
More information1 Linear programming relaxation
Cornell University, Fall 2010 CS 6820: Algorithms Lecture notes: Primal-dual min-cost bipartite matching August 27 30 1 Linear programming relaxation Recall that in the bipartite minimum-cost perfect matching
More informationCentral issues in modelling
Central issues in modelling Construct families of curves, surfaces and volumes that can represent common objects usefully; are easy to interact with; interaction includes: manual modelling; fitting to
More informationVoluntary State Curriculum Algebra II
Algebra II Goal 1: Integration into Broader Knowledge The student will develop, analyze, communicate, and apply models to real-world situations using the language of mathematics and appropriate technology.
More information1 EquationsofLinesandPlanesin 3-D
1 EquationsofLinesandPlanesin 3-D Recall that given a point P (a, b, c), one can draw a vector from the origin to P. Such a vector is called the position vector of the point P and its coordinates are a,
More informationErasing Haar Coefficients
Recap Haar simple and fast wavelet transform Limitations not smooth enough: blocky How to improve? classical approach: basis functions Lifting: transforms 1 Erasing Haar Coefficients 2 Page 1 Classical
More information