Piecewise polynomial interpolation
|
|
- Shawn Sparks
- 5 years ago
- Views:
Transcription
1 Chapter 2 Piecewise polynomial interpolation In ection.6., and in Lab, we learned that it is not a good idea to interpolate unctions by a highorder polynomials at equally spaced points. However, it transpires that it is possible to obtain a very good approximations using a very simple method. The trick is to use a spline: a piecewise polynomial interpolating unction. We ll consider three important example o splines:. linear splines 2. (natural) cubic splines. 3. Hermite piecewise cubics. In this section, we always have N equally spaced points: let h = (b a)/n, then a = x, b = x N and x i = x + ih. Oten these are reerred to as knots points (or simply as knots), and denote the set o knot points by ω N := {x i } N i=. For more details about splines, have a look at [M3, Chap. ], and [98, Lectures and ]. 2.. Construction Deinition 2... Let be a unction that is continuous on [a, b]. The linear spline interpolant to is the continuous unction l such that (i) l(x i ) = (x i ) or each i =,,..., N, (ii) l is a linear unction l i on each interval [x i, x i ]. It is easy to write down a ormula or the l i, based on Lagrange polynomials. et h = (b a)/n. Then, or x [x i, x i ] l i (x) = (x i ) x i x h + (x i ) x x i. (2.) h Example Write down the linear spline interpolant to (x) = e x at the knot points {,, }. 2. Linear Interpolating plines We irst study the piecewise linear interpolant, also called a linear spline. We will see that they have important properties, including (a) they are easy to construct and analyse; (b) the bound on the error decreases as the number o interpolation points increases; (c) the error we get using a linear spline is no more than twice the error using the best possible (piecewise linear) approximation; amd (d) o all the interpolants to at a given set o points, the linear spline is the one with the smallest st derivative Analysis We know that i p N is the polynomial o degree N that interpolates at N equally spaced points, it does not ollow that p N as n. But as we will see, the piecewise linear interpolant to converges to, albeit slowly. This is veriied in the ollowing theorem, which is a direct consequence o Theorem.3.3.
2 2 Theorem uppose that, and are all continuous and deined on the interval [a, b]. Let l be the linear spline interpolant to on the N + equally spaced points a = x < x < x N = b with h = x i x i = (b a)/n. Then l h2 8, (Here, as usual, g is deined as max a x b g(x).) (x) l 4 (x) 5 5 (x).9 l 8 (x) (x).9 l 6 (x) It now ollows directly rom Theorem 2..3 that lim l =. n Example Figure 2. shows linear spline interpolations o Runge s example: (x) = on [ 5, 5]. + x2 These diagrams appear to support our assertion that the error tends to zero as n. Example uppose you are interpolating (x) = e x on N equally spaced intervals between x = and x N =. What value o N would you have to take to ensure that the maximum error is less than 2? 5 5 Fig. 2.: Linear spline interpolants to /( + x 2 ) on [ 5, 5]. Contrast with Figure Best approximation For the next part o the analysis it will help to think o piecewise linear interpolation as an operator. Then we can compare the linear spline to all the other piecewise linear approximations. First, observe that one can deine an ininite number o piecewise linear unctions on a given set o knot points ω N. We ll call the set o these unctions L. Deinition For a ixed set o knot points ω N, let L be the operator that maps the continuous unction to its linear spline interpolant l L. Now suppose that g L. Then L(g) = g. That is L is a projection: L(L()) = L(). It is not hard to see that one could ind a dierent unction ˆl L that is a better approximation o in sense that max (x) ˆl(x) < max (x) l(x). x x x n x x x n However, l is very easy to ind, and the associated error is no worse than twice (x) ˆl(x). Theorem 2..7 (tewart s Aternotes goes to grad school,
3 2 Lecture ). Let l = L(). For all ˆl L, l 2 ˆl. (That L is a projection is key to the proo.) 2..5 Exercises Exercise 2.. Page 28 o the Department o Education s old Mathematics Tables ( The Log Tables ) reports that ln() =, ln(.5) =.455 and ln(2) =.693. (i) Write down the linear spline l that interpolates (x) = ln(x) at the points x =, x =.5 and x 2 = 2. (ii) Use this to estimate ln(x) at x =.2. How does this compare to the value in the tables? (.823) (iii) Give an estimate or the maximum error: max (x) l(x). x 2 (iv) What value o n would you choose to ensure that (x) l(x). or all x [, 2] Minimum Energy The inal interesting property o l that we will study is called the minimum energy property. Deinition Let u be a unction that is continuous and deined on the interval [a, b] except, maybe, at the (countable set) ω N o knot points Then the 2-norm o u is u 2,[a,b] := ( b u 2 (x)dx ) /2. Usually we just write this as u 2. Let H be the set o all unctions u that are continuous on [a, b] and have u 2 <. Note that l H, even though we have not properly deined l at the mesh points ω N. Theorem 2..9 (üli and Mayers, Thm..2). Let w be any unction in H that interpolates the unction at the points in ω N. Let l we the linear spline interpolant o. Then l 2 w 2. a Exercise 2.2. As an alternative to (2.), one can deine the linear spline interpolant to a unction is as a linear combination o a set o piecewise linear basis unctions {ψ i } N i= : { i = j ψ i (x j ) = i j (i) Write down a ormula or the ψ i (x); (ii) derive a ormula or l(x) in terms o the ψ i. This exercise is useul: we ll be using these basis unctions (called hat unctions) in the inal section o the course. ψ i (x) ψ i (x) ψ i+ (x) x i x i x i+ Fig. 2.2: ome hat unctions More precisely, we should say everywhere, except on a set o measure zero. What that means is very important to the areas o measure theory and unctional analysis. However, since some o you are not taking those courses until next year, we ll air-brush over the exact deinition.
4 Cubic plines The Cubic pline 2 is perhaps the most useul and popular interpolating unction used in numerical analysis, automotive and aeronautical engineering, computer and ilm animation, digital photography, inancial modelling, and as many other areas as there are human endeavours where continuous processes are modelled by discrete ones. One could argue that the most undamental shortcoming o the linear spline interpolant to (see, or example, Figure 2.) is that they are not smooth : that is, although l i (x i ) = l i+ (x i ), it is not generally the case that l i (x i) = l i+ (x i). Also, the interpolating unction cannot capture the curvature o. (I you think about this last statement, you ll see that it can be expressed as l (x) = or all x [a, b].) Cubic splines try to balance the simplicity o the linear spline approach by using a low-order polynomial on each interval with the desire or curvature and more smoothness by using cubics, and by orcing adjacent ones, and their irst and second derivatives, to agree at knot points. Deinition Let be a unction that is continuous on [a, b]. The cubic spline interpolant to is the continuous unction such that (i) or i =,..., N, on each interval [x i, x i ] let (x) = s i (x), where each o the s i is a cubic polynomial. (ii) s i (x i ) = (x i ) or i =,..., N, (iii) s i (x i ) = (x i ) or i =,..., N, (iv) s i (x i) = s i+ (x i) or i =,..., N, (v) s i (x i) = s i+ (x i) or i =,..., N. o, we have deined the cubic spline as a unction that interpolates at N + points, has continuous irst and second derivatives on [x, x N ] and is a cubic polynomial on each o the n intervals [x i, x i ]. That is, it is piecewise cubic. We know that one can write a cubic as p 3 (x) = a + a x + a 2 x 2 + a 3 x 3. Thus it takes 4 terms to uniquely deine a single cubic. To uniquely determine N cubics requires 4N terms. They can be ound by solving 4N (linearly independent) equations. But an inspection o (ii) (iv) above reveals only 4N 2 equations. 2 It is generally accepted that (mathematical) splines were irst described by Isaac choenberg (93-99) during WW2. However their physical realisation had been used in the ship-building and aircrat industries prior to that. The missing equations can be chosen in a number o ways: (i) by setting (x ) = and (x N ) =. This is called a natural spline, and is the approach we ll take. (ii) by setting (x ) = and (x N ) =. called a clamped spline. This is (iii) set (x ) = (x N ) and (x ) = (x N ). This is the periodic spline and is used or interpolating, say, trigonometric unctions. (iv) only use N 2 components o the spline: s 2,..., s N. But extend to two end ones so that s 2 (x ) = (x ) and s N = (x N ). This is called the not-aknot condition. In Figure 2.3 we show the natural cubic spline interpolants to (x) = /( + x 2 ) (with a = 5 and b = 5 as usual) on N = 4, 8 and 6 intervals. Compared with Figure 2., we seem to get signiicantly better approximation. Moreover, the rate at which the error decreases with respect to h is much aster..5 Cubic pline on 4 intervals Cubic pline on 8 intervals Cubic pline on 6 intervals Fig. 2.3: Cubic spline interpolants to /( + x 2 ). Compare with Figure 2.
5 Constructing Cubic plines ome notation: ˆ h = x i x i = (x N x )/N or all i. ˆ i := (x i ). To construct the spline, irst observe that i is piecewise cubic, then is piecewise quadratic, and is piecewise linear. (i) Let σ i = (x i ) or i =,..., N. and so, or i =, 2,..., N α i = i h h 6 σ i, Notice that this gives β i = α i. β i = i h h 6 σ i. (2.3) (iv) o, now we just need the equations or σ i. Two o these come rom the act that this is a natural cubic spline, with (x ) = and (x N ) =, thereore σ = and σ N =. The remaining N equations come rom the act that is continuous at x, x 2,..., x N. o we set that s i (x i ) = s i+ (x i ). Ater some work (see Exercise 2.3), we can show this means that the system is: σ =, (2.4a) 6 (σ i + 4σ i + σ i+ ) = h 2 ( i 2 i + i+ ) (2.4b) or i =,..., N, (ii) Integrate twice: σ N = A cubic spline example (2.4c) Example Find the natural cubic spline interpolant to at the points x =, x =, x 2 = 2 and x 3 = 3 where =, = 2, 2 = and 3 =. [We won t do all the details in class, so here they are] olution: From (2.2) we see we are looking or (x)= We get s i (x) = α i (x x i ) + β i (x i x)+ σ i 6h (x i x) 3 + σ i 6h (x x i ) 3, (2.2) or x [x i, x i ]. (iii) The α i and β i arose as the constants o integration. To ind them: α x + β ( x) + σ 6h ( x)3 + σ 6h x3,x [, ], α 2 (x ) + β 2 (2 x) + σ 6h (2 x)3 + σ 2 6h (x )3,x [, 2], α 3 (x x 2 ) + β 3 (3 x) + σ 2 6h (3 x)3 + σ 3 6h (x 2)3,x [2, 3]. We irst solve or the σ i using (2.4): σ 4 σ 8 = 6 This gives 4 σ 2 σ 3 σ =, σ = 24 5, σ 2 = 6 5, σ 3 =. Now use (2.3) to get α = 4 5, α 2 = 4 5, α 3 =. and β = β 2 = 4 5, β 3 = 4 5.
6 24 The answer is 4 x x3, x [, ], 4 4 (x ) + (2 x) 5 5 (x) = 4 5 (2 x)3 + (x 5 )3, x [, 2], 4 (3 x) + (3 5 5 x)3. x [2, 3]. This is shown in Figure 2.4. Example What is the smallest value o N that you must take to ensure that, i interpolating (x) = e x on N equally sized intervals, the error is less that 8? How does this compare with a linear spline interpolation (see Example 2..5)? (x i,y i ) Fig. 2.4: pline interpolant rom Example Error Estimates We state the ollowing error estimates without proo. You don t need to know these, or be able to prove them. But you should be able to use them i required. Theorem I C 4 [a, b] and is its cubic spline interpolant on N + equally spaced points, then where M 4 = (iv) M 4h 4, 24 M 4h 3, 3 8 M 4h 2, Example Give an upper bound on the error or the cubic spline interpolant to = e x on the interval [, ] with N = mesh points. Recall Theorem 2..9 or linear splines. There is an analogous result or natural cubic splines. Theorem Let u be any unction that interpolates at ω N, and is such that u H 2 (x, x N ). Then 2 u 2. The proo is not hard - it is analogous to the proo o Theorem Exercises Exercise 2.3. When deducing the system o equations or the natural cubic spline, we showed how to construct the ormulation given in (2.2) and the relationship between σ i, α i and β i in (2.3). Now careully show to deduce the system (2.4). Exercise 2.4. (For students who did MA385). Write the equations in (2.4) as a linear system Aσ = b, where A is an n n matrix. how that A is nonsingular, and hence that the system as a unique solution. Exercise 2.5. Find the natural cubic spline interpolant to (x) = sin(πx/2) at the nodes {x i } 3 i= = {,, 2, 3}. Calculate value o the interpolant at x = 2.5. What is the error at this point? Exercise 2.6. Take (x) = ln(x), x =, x N = 2. What value o N would you have to take to ensure that ln(x) (x) 4 or all x [, 2]?
7 Piecewise Hermite Interpolation In this section we introduce another type o cubic spline. It is not as smooth as the natural spline o the precious section only it and its irst derivative are continuous on [x, x N ] and it requires that we know (x i ). But its easier to construct (don t have to solve a linear system) and analyse than that natural spline. As beore, or short-hand, we write (x i ) as i and (x i ) as i. And, as ever, we ll simpliy the analysis by taking the points to be equally spaced: x i x i = h or each i) PCHIP Deinition Given a set o interpolation points x < x < < x N, the Piecewise Cubic Hermite pline Interpolant (PCHIP),, to the unction, satisies (i) C [x, x N ], (ii) (x i ) = (x i ) and (x i ) = (x i ) or i =,,..., N. (iii) On each interval [x i, x i ], is a cubic polynomial. We can construct these splines as ollows. On each interval [x i, x i ], let be the cubic polynomial i given by i (x) = c + c (x x i )+ c 2 (x x i ) 2 + c 3 (x x i ) 3, (2.5) Then we can show... Figure 2.5 shows some PCHIP interpolants to (x) = /( + x 2 ) on the interval [ 5, 5]. Compare with Figures 2. and PCHIP pline on 4 intervals PCHIP pline on 8 intervals PCHIP pline on 6 intervals Fig. 2.5: PCHIP interpolants to /( + x 2 ). Compare with Figure 2.3 We now want to prove an error estimate. The norm we use is := max x x x N (x) (x) Theorem Let C 4 [x, x N ] and let be the Hermite Cubic pline interpolant to it at the N equally spaced points a = x < x < < x N = b. Then h4 384 (iv). This gives that c = i, c = i, c 2 = 3 h 2 ( i i ) h ( i + 2 i ), c 3 = h 2 ( i + i ) 2 h 3 ( i i ).
8 Estimating derivatives You can ignore this section. It is included or completeness. It should make sense to anyone who took MA385, given that it relies only on Taylor series. As we deined the PCHIP interpolant one needs to know in order to be able to construct it. However, you might notice that the MATLAB unction >> pchip only requires, not. How does it do this, one might wonder? It turns out that there are important variants on the PHCIP scheme that don t involve knowing,,..., N, but instead uses approximations or. There are several approaches or deriving these estimates, including (i) Geometrically. Recall, all we are really looking or is the slope o the tangent to at x = x i. (ii) By inding a polynomial that interpolates and dierentiating that. (iii) Taylor s Theorem. For more details see [96, Lecture 24]. Or look back at your notes or MA385, Lecture 4. i lope is (x i ) (x i ) x i x i i lope is (x = lim i ) (x i σ) σ σ lope is i+ x i x i x i+ Fig. 2.6: Estimating (x i )dx As suggested by Figure 2.6 one could take (x i ) i i x i x i = h ( i i ) (x i+ ) (x i ) x i+ x i Or (x i ) i+ i = x i+ x i h ( i+ i ). ince one seems to over estimate (x i ), and the other seems to under-estimate, we could take the average: (x i ) i + i+ x i+ x i = 2h ( i + i ) These approximations can also be easily derived rom Taylor s Theorem: or some η (a, b). Using this we can get i+ = i + h (x i ) + h2 2 (x i ) + h3 6 (τ ), i = i h (x i ) + h2 2 (x i ) h3 6 (τ 2 ), or τ (x i, x i+ ) and τ 2 (x i, x i ). Now subtract the st equation rom the second to get the ormula and error estimate Exercises Exercise 2.7. Recall Exercise 2.5. Calculate the value to the PCHIP interpolant to (x) = sin(πx/2) at the nodes {x i } 3 i= = {,, 2, 3} at the point x = 2.5. What is the error at this point? Exercise 2.8. Let (x) = ln(x). Let l and be the piecewise linear and Hermite cubic spline interpolants (respectively) to on n+ equally spaced points = x < x < < x N = 2. What value o n would you have to take to ensure that (i) max x 2 (x) l(x) 4? (ii) max x 2 (x) (x) 4? Also, using the code you developed in Lab 2, compare these theoretical results with the value needed on practice. Exercise 2.9. There are ways o constructing the PCHIP, other than (2.5). For example, let s = x x k, then (x) = h3 3hs 2 + 2s 3 h 3 k + 3hs2 2s 3 h 3 k + s(s h) 2 h 2 k + s2 (s h) h 2 k how that this is the same as the PCHIP. (x) = (a) + (x a) (x a)2 (a) + (a)+ 2 (x a) 3 (x a)n (a) + + (n) (a)+ 3! n! (x a) n+ (n+) (η), (n + )!
February 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 informationES 240: Scientific and Engineering Computation. a function f(x) that can be written as a finite series of power functions like
Polynomial Deinition a unction () that can be written as a inite series o power unctions like n is a polynomial o order n n ( ) = A polynomial is represented by coeicient vector rom highest power. p=[3-5
More information8 Piecewise Polynomial Interpolation
Applied Math Notes by R. J. LeVeque 8 Piecewise Polynomial Interpolation 8. Pitfalls of high order interpolation Suppose we know the value of a function at several points on an interval and we wish to
More informationPolynomials tend to oscillate (wiggle) a lot, even when our true function does not.
AMSC/CMSC 460 Computational Methods, Fall 2007 UNIT 2: Spline Approximations Dianne P O Leary c 2001, 2002, 2007 Piecewise polynomial interpolation Piecewise polynomial interpolation Read: Chapter 3 Skip:
More informationGenerell Topologi. Richard Williamson. May 6, 2013
Generell Topologi Richard Williamson May 6, Thursday 7th January. Basis o a topological space generating a topology with a speciied basis standard topology on R examples Deinition.. Let (, O) be a topological
More informationSome Inequalities Involving Fuzzy Complex Numbers
Theory Applications o Mathematics & Computer Science 4 1 014 106 113 Some Inequalities Involving Fuzzy Complex Numbers Sanjib Kumar Datta a,, Tanmay Biswas b, Samten Tamang a a Department o Mathematics,
More informationConcavity. Notice the location of the tangents to each type of curve.
Concavity We ve seen how knowing where a unction is increasing and decreasing gives a us a good sense o the shape o its graph We can reine that sense o shape by determining which way the unction bends
More informationLab 9 - GEOMETRICAL OPTICS
161 Name Date Partners Lab 9 - GEOMETRICAL OPTICS OBJECTIVES Optics, developed in us through study, teaches us to see - Paul Cezanne Image rom www.weidemyr.com To examine Snell s Law To observe total internal
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 informationGEOMETRICAL OPTICS OBJECTIVES
Geometrical Optics 207 Name Date Partners OBJECTIVES OVERVIEW GEOMETRICAL OPTICS To examine Snell s Law and observe total internal relection. To understand and use the lens equations. To ind the ocal length
More information9.8 Graphing Rational Functions
9. Graphing Rational Functions Lets begin with a deinition. Deinition: Rational Function A rational unction is a unction o the orm P where P and Q are polynomials. Q An eample o a simple rational unction
More informationAPPM/MATH Problem Set 4 Solutions
APPM/MATH 465 Problem Set 4 Solutions This assignment is due by 4pm on Wednesday, October 16th. You may either turn it in to me in class on Monday or in the box outside my office door (ECOT 35). Minimal
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 informationSkill Sets Chapter 5 Functions
Skill Sets Chapter 5 Functions No. Skills Examples o questions involving the skills. Sketch the graph o the (Lecture Notes Example (b)) unction according to the g : x x x, domain. x, x - Students tend
More informationScientific Computing: Interpolation
Scientific Computing: Interpolation Aleksandar Donev Courant Institute, NYU donev@courant.nyu.edu Course MATH-GA.243 or CSCI-GA.22, Fall 25 October 22nd, 25 A. Donev (Courant Institute) Lecture VIII /22/25
More informationDerivative. Bernstein polynomials: Jacobs University Visualization and Computer Graphics Lab : ESM4A - Numerical Methods 313
Derivative Bernstein polynomials: 120202: ESM4A - Numerical Methods 313 Derivative Bézier curve (over [0,1]): with differences. being the first forward 120202: ESM4A - Numerical Methods 314 Derivative
More informationMultiple-Choice Test Spline Method Interpolation COMPLETE SOLUTION SET
Multiple-Choice Test Spline Method Interpolation COMPLETE SOLUTION SET 1. The ollowing n data points, ( x ), ( x ),.. ( x, ) 1, y 1, y n y n quadratic spline interpolation the x-data needs to be (A) equally
More informationL ENSES. Lenses Spherical refracting surfaces. n 1 n 2
Lenses 2 L ENSES 2. Sherical reracting suraces In order to start discussing lenses uantitatively, it is useul to consider a simle sherical surace, as shown in Fig. 2.. Our lens is a semi-ininte rod with
More informationMAPI Computer Vision. Multiple View Geometry
MAPI Computer Vision Multiple View Geometry Geometry o Multiple Views 2- and 3- view geometry p p Kpˆ [ K R t]p Geometry o Multiple Views 2- and 3- view geometry Epipolar Geometry The epipolar geometry
More informationConsider functions such that then satisfies these properties: So is represented by the cubic polynomials on on and on.
1 of 9 3/1/2006 2:28 PM ne previo Next: Trigonometric Interpolation Up: Spline Interpolation Previous: Piecewise Linear Case Cubic Splines A piece-wise technique which is very popular. Recall the philosophy
More informationNeighbourhood Operations
Neighbourhood Operations Neighbourhood operations simply operate on a larger neighbourhood o piels than point operations Origin Neighbourhoods are mostly a rectangle around a central piel Any size rectangle
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 informationChapter 3 Image Enhancement in the Spatial Domain
Chapter 3 Image Enhancement in the Spatial Domain Yinghua He School o Computer Science and Technology Tianjin University Image enhancement approaches Spatial domain image plane itsel Spatial domain methods
More informationLecture 8. Divided Differences,Least-Squares Approximations. Ceng375 Numerical Computations at December 9, 2010
Lecture 8, Ceng375 Numerical Computations at December 9, 2010 Computer Engineering Department Çankaya University 8.1 Contents 1 2 3 8.2 : These provide a more efficient way to construct an interpolating
More informationA Proposed Approach for Solving Rough Bi-Level. Programming Problems by Genetic Algorithm
Int J Contemp Math Sciences, Vol 6, 0, no 0, 45 465 A Proposed Approach or Solving Rough Bi-Level Programming Problems by Genetic Algorithm M S Osman Department o Basic Science, Higher Technological Institute
More informationUsing a Projected Subgradient Method to Solve a Constrained Optimization Problem for Separating an Arbitrary Set of Points into Uniform Segments
Using a Projected Subgradient Method to Solve a Constrained Optimization Problem or Separating an Arbitrary Set o Points into Uniorm Segments Michael Johnson May 31, 2011 1 Background Inormation The Airborne
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 informationAsymptotic Error Analysis
Asymptotic Error Analysis Brian Wetton Mathematics Department, UBC www.math.ubc.ca/ wetton PIMS YRC, June 3, 2014 Outline Overview Some History Romberg Integration Cubic Splines - Periodic Case More History:
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 informationLarger K-maps. So far we have only discussed 2 and 3-variable K-maps. We can now create a 4-variable map in the
EET 3 Chapter 3 7/3/2 PAGE - 23 Larger K-maps The -variable K-map So ar we have only discussed 2 and 3-variable K-maps. We can now create a -variable map in the same way that we created the 3-variable
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 informationCurve Representation ME761A Instructor in Charge Prof. J. Ramkumar Department of Mechanical Engineering, IIT Kanpur
Curve Representation ME761A Instructor in Charge Prof. J. Ramkumar Department of Mechanical Engineering, IIT Kanpur Email: jrkumar@iitk.ac.in Curve representation 1. Wireframe models There are three types
More informationCurve fitting using linear models
Curve fitting using linear models Rasmus Waagepetersen Department of Mathematics Aalborg University Denmark September 28, 2012 1 / 12 Outline for today linear models and basis functions polynomial regression
More informationComputer Data Analysis and Plotting
Phys 122 February 6, 2006 quark%//~bland/docs/manuals/ph122/pcintro/pcintro.doc Computer Data Analysis and Plotting In this lab we will use Microsot EXCEL to do our calculations. This program has been
More informationParameterization. Michael S. Floater. November 10, 2011
Parameterization Michael S. Floater November 10, 2011 Triangular meshes are often used to represent surfaces, at least initially, one reason being that meshes are relatively easy to generate from point
More informationFall CSCI 420: Computer Graphics. 4.2 Splines. Hao Li.
Fall 2014 CSCI 420: Computer Graphics 4.2 Splines Hao Li http://cs420.hao-li.com 1 Roller coaster Next programming assignment involves creating a 3D roller coaster animation We must model the 3D curve
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 informationlecture 10: B-Splines
9 lecture : -Splines -Splines: a basis for splines Throughout our discussion of standard polynomial interpolation, we viewed P n as a linear space of dimension n +, and then expressed the unique interpolating
More informationUNIT #2 TRANSFORMATIONS OF FUNCTIONS
Name: Date: UNIT # TRANSFORMATIONS OF FUNCTIONS Part I Questions. The quadratic unction ollowing does,, () has a turning point at have a turning point? 7, 3, 5 5, 8. I g 7 3, then at which o the The structure
More informationThree-Dimensional Object Representations Chapter 8
Three-Dimensional Object Representations Chapter 8 3D Object Representation A surace can be analticall generated using its unction involving the coordinates. An object can be represented in terms o its
More informationIn some applications it may be important that the extrema of the interpolating function are within the extrema of the given data.
Shape-preserving piecewise poly. interpolation In some applications it may be important that the extrema of the interpolating function are within the extrema of the given data. For example: If you the
More informationThe Rational Zero Theorem
The Rational Zero Theorem Our goal in this section is to learn how we can ind the rational zeros o the polynomials. For example: x = x 4 + x x x + ( ) We could randomly try some actors and use synthetic
More informationRemark. Jacobs University Visualization and Computer Graphics Lab : ESM4A - Numerical Methods 331
Remark Reconsidering the motivating example, we observe that the derivatives are typically not given by the problem specification. However, they can be estimated in a pre-processing step. A good estimate
More informationTHE FINANCIAL CALCULATOR
Starter Kit CHAPTER 3 Stalla Seminars THE FINANCIAL CALCULATOR In accordance with the AIMR calculator policy in eect at the time o this writing, CFA candidates are permitted to use one o two approved calculators
More informationReflection and Refraction
Relection and Reraction Object To determine ocal lengths o lenses and mirrors and to determine the index o reraction o glass. Apparatus Lenses, optical bench, mirrors, light source, screen, plastic or
More informationKnow it. Control points. B Spline surfaces. Implicit surfaces
Know it 15 B Spline Cur 14 13 12 11 Parametric curves Catmull clark subdivision Parametric surfaces Interpolating curves 10 9 8 7 6 5 4 3 2 Control points B Spline surfaces Implicit surfaces Bezier surfaces
More information/633 Introduction to Algorithms Lecturer: Michael Dinitz Topic: Approximation algorithms Date: 11/27/18
601.433/633 Introduction to Algorithms Lecturer: Michael Dinitz Topic: Approximation algorithms Date: 11/27/18 22.1 Introduction We spent the last two lectures proving that for certain problems, we can
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 informationGlobal Constraints. Combinatorial Problem Solving (CPS) Enric Rodríguez-Carbonell (based on materials by Javier Larrosa) February 22, 2019
Global Constraints Combinatorial Problem Solving (CPS) Enric Rodríguez-Carbonell (based on materials by Javier Larrosa) February 22, 2019 Global Constraints Global constraints are classes o constraints
More informationRough Connected Topologized. Approximation Spaces
International Journal o Mathematical Analysis Vol. 8 04 no. 53 69-68 HIARI Ltd www.m-hikari.com http://dx.doi.org/0.988/ijma.04.4038 Rough Connected Topologized Approximation Spaces M. J. Iqelan Department
More informationx n x n stepnumber k order r error constant C r+1 1/2 5/12 3/8 251/720 abs. stab. interval (α,0) /11-3/10
MATH 573 LECTURE NOTES 77 13.8. Predictor-corrector methods. We consider the Adams methods, obtained from the formula xn+1 xn+1 y(x n+1 y(x n ) = y (x)dx = f(x,y(x))dx x n x n by replacing f by an interpolating
More informationLecture 9: Introduction to Spline Curves
Lecture 9: Introduction to Spline Curves Splines are used in graphics to represent smooth curves and surfaces. They use a small set of control points (knots) and a function that generates a curve through
More informationMath 2A Vector Calculus Chapter 11 Test Fall 07 Name Show your work. Don t use a calculator. Write responses on separate paper.
Math A Vector Calculus Chapter Test Fall 7 Name Show our work. Don t use a calculator. Write responses on separate paper.. Consider the nice, smooth unction z, whose contour map is shown at right. a. Estimate
More informationAutomated Theorem Proving and Proof Checking. Proof Checking. Motivation. Cunning Plan. Standard Architecture. History
Automated Theorem Proving and Proo Checking Engler: Automatically Generating Malicious Disks using Symex IEEE Security and Privacy 2006 Use CIL and Symbolic Execution on Linux FS code Special model o memory,
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 informationThe Bisection Method versus Newton s Method in Maple (Classic Version for Windows)
The Bisection Method versus (Classic Version for Windows) Author: Barbara Forrest Contact: baforres@uwaterloo.ca Copyrighted/NOT FOR RESALE version 1.1 Contents 1 Objectives for this Lab i 2 Approximate
More information(Refer Slide Time: 00:02:24 min)
CAD / CAM Prof. Dr. P. V. Madhusudhan Rao Department of Mechanical Engineering Indian Institute of Technology, Delhi Lecture No. # 9 Parametric Surfaces II So these days, we are discussing the subject
More informationTutte Embeddings of Planar Graphs
Spectral Graph Theory and its Applications Lectre 21 Ttte Embeddings o Planar Graphs Lectrer: Daniel A. Spielman November 30, 2004 21.1 Ttte s Theorem We sally think o graphs as being speciied by vertices
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 informationMath-3 Lesson 1-7 Analyzing the Graphs of Functions
Math- Lesson -7 Analyzing the Graphs o Functions Which unctions are symmetric about the y-axis? cosx sin x x We call unctions that are symmetric about the y -axis, even unctions. Which transormation is
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 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 informationAn introduction to interpolation and splines
An introduction to interpolation and splines Kenneth H. Carpenter, EECE KSU November 22, 1999 revised November 20, 2001, April 24, 2002, April 14, 2004 1 Introduction Suppose one wishes to draw a curve
More informationFractal Art: Fractal Image and Music Generator
Proceedings o the 7th WSEAS Int. Con. on Signal Processing, Computational Geometry & Artiicial Vision, Athens, Greece, August 4-6, 007 159 Fractal Art: Fractal Image and Music Generator RAZVAN TANASIE
More informationPentagon contact representations
Pentagon contact representations Stean Felsner Institut ür Mathematik Technische Universität Berlin Hendrik Schrezenmaier Institut ür Mathematik Technische Universität Berlin August, 017 Raphael Steiner
More informationAM205: lecture 2. 1 These have been shifted to MD 323 for the rest of the semester.
AM205: lecture 2 Luna and Gary will hold a Python tutorial on Wednesday in 60 Oxford Street, Room 330 Assignment 1 will be posted this week Chris will hold office hours on Thursday (1:30pm 3:30pm, Pierce
More informationHandout 4 - Interpolation Examples
Handout 4 - Interpolation Examples Middle East Technical University Example 1: Obtaining the n th Degree Newton s Interpolating Polynomial Passing through (n+1) Data Points Obtain the 4 th degree Newton
More informationCOMP 250 Fall graph traversal Nov. 15/16, 2017
Graph traversal One problem we oten need to solve when working with graphs is to decide i there is a sequence o edges (a path) rom one vertex to another, or to ind such a sequence. There are various versions
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 informationInterpolation. TANA09 Lecture 7. Error analysis for linear interpolation. Linear Interpolation. Suppose we have a table x x 1 x 2...
TANA9 Lecture 7 Interpolation Suppose we have a table x x x... x n+ Interpolation Introduction. Polynomials. Error estimates. Runge s phenomena. Application - Equation solving. Spline functions and interpolation.
More informationInterpolation. Chapter Information from Data
Chapter 5 Interpolation 5. Information from Data The topic of this chapter is interpolation, which relates to passing a curve through a set of data points. To put this in context, extracting information
More informationMaximizing an interpolating quadratic
Week 11: Monday, Apr 9 Maximizing an interpolating quadratic Suppose that a function f is evaluated on a reasonably fine, uniform mesh {x i } n i=0 with spacing h = x i+1 x i. How can we find any local
More informationSet 5, Total points: 100 Issued: week of
Prof. P. Koumoutsakos Prof. Dr. Jens Walther ETH Zentrum, CLT F 1, E 11 CH-809 Zürich Models, Algorithms and Data (MAD): Introduction to Computing Spring semester 018 Set 5, Total points: 100 Issued: week
More informationMATRIX ALGORITHM OF SOLVING GRAPH CUTTING PROBLEM
UDC 681.3.06 MATRIX ALGORITHM OF SOLVING GRAPH CUTTING PROBLEM V.K. Pogrebnoy TPU Institute «Cybernetic centre» E-mail: vk@ad.cctpu.edu.ru Matrix algorithm o solving graph cutting problem has been suggested.
More informationLecture VIII. Global Approximation Methods: I
Lecture VIII Global Approximation Methods: I Gianluca Violante New York University Quantitative Macroeconomics G. Violante, Global Methods p. 1 /29 Global function approximation Global methods: function
More informationSplines. Parameterization of a Curve. Curve Representations. Roller coaster. What Do We Need From Curves in Computer Graphics? Modeling Complex Shapes
CSCI 420 Computer Graphics Lecture 8 Splines Jernej Barbic University of Southern California Hermite Splines Bezier Splines Catmull-Rom Splines Other Cubic Splines [Angel Ch 12.4-12.12] Roller coaster
More informationCalculus I Review Handout 1.3 Introduction to Calculus - Limits. by Kevin M. Chevalier
Calculus I Review Handout 1.3 Introduction to Calculus - Limits by Kevin M. Chevalier We are now going to dive into Calculus I as we take a look at the it process. While precalculus covered more static
More informationEngineering Analysis ENG 3420 Fall Dan C. Marinescu Office: HEC 439 B Office hours: Tu-Th 11:00-12:00
Engineering Analysis ENG 3420 Fall 2009 Dan C. Marinescu Office: HEC 439 B Office hours: Tu-Th 11:00-12:00 1 Lecture 24 Attention: The last homework HW5 and the last project are due on Tuesday November
More informationChapter 6. Curves and Surfaces. 6.1 Graphs as Surfaces
Chapter 6 Curves and Surfaces In Chapter 2 a plane is defined as the zero set of a linear function in R 3. It is expected a surface is the zero set of a differentiable function in R n. To motivate, graphs
More informationTheoretical calculation of the natural flock restoration time of a cotton comber
Indian Journal o Fibre & Textile esearch Vol. 40, September 05, pp. 88-9 Theoretical calculation o the natural lock restoration time o a cotton comber Li Xin-rong a, Jiang Xiu-ming, Yang Jian-cheng & Wang
More informationParametric curves. Reading. Curves before computers. Mathematical curve representation. CSE 457 Winter Required:
Reading Required: Angel 10.1-10.3, 10.5.2, 10.6-10.7, 10.9 Parametric curves CSE 457 Winter 2014 Optional Bartels, Beatty, and Barsky. An Introduction to Splines for use in Computer Graphics and Geometric
More informationMA 323 Geometric Modelling Course Notes: Day 10 Higher Order Polynomial Curves
MA 323 Geometric Modelling Course Notes: Day 10 Higher Order Polynomial Curves David L. Finn December 14th, 2004 Yesterday, we introduced quintic Hermite curves as a higher order variant of cubic Hermite
More informationInterpolation - 2D mapping Tutorial 1: triangulation
Tutorial 1: triangulation Measurements (Zk) at irregular points (xk, yk) Ex: CTD stations, mooring, etc... The known Data How to compute some values on the regular spaced grid points (+)? The unknown data
More informationMTH 122 Calculus II Essex County College Division of Mathematics and Physics 1 Lecture Notes #11 Sakai Web Project Material
MTH Calculus II Essex County College Division of Mathematics and Physics Lecture Notes # Sakai Web Project Material Introduction - - 0 - Figure : Graph of y sin ( x y ) = x cos (x + y) with red tangent
More informationEXERCISES SHORTEST PATHS: APPLICATIONS, OPTIMIZATION, VARIATIONS, AND SOLVING THE CONSTRAINED SHORTEST PATH PROBLEM. 1 Applications and Modelling
SHORTEST PATHS: APPLICATIONS, OPTIMIZATION, VARIATIONS, AND SOLVING THE CONSTRAINED SHORTEST PATH PROBLEM EXERCISES Prepared by Natashia Boland 1 and Irina Dumitrescu 2 1 Applications and Modelling 1.1
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 informationLecture 9. Curve fitting. Interpolation. Lecture in Numerical Methods from 28. April 2015 UVT. Lecture 9. Numerical. Interpolation his o
Curve fitting. Lecture in Methods from 28. April 2015 to ity Interpolation FIGURE A S Splines Piecewise relat UVT Agenda of today s lecture 1 Interpolation Idea 2 3 4 5 6 Splines Piecewise Interpolation
More informationMaster Thesis. Comparison and Evaluation of Didactic Methods in Numerical Analysis for the Teaching of Cubic Spline Interpolation
Master Thesis Comparison and Evaluation of Didactic Methods in Numerical Analysis for the Teaching of Cubic Spline Interpolation Abtihal Jaber Chitheer supervised by Prof. Dr. Carmen Arévalo May 17, 2017
More informationSYSTEMS OF NONLINEAR EQUATIONS
SYSTEMS OF NONLINEAR EQUATIONS Widely used in the mathematical modeling of real world phenomena. We introduce some numerical methods for their solution. For better intuition, we examine systems of two
More informationSung-Eui Yoon ( 윤성의 )
CS480: Computer Graphics Curves and Surfaces Sung-Eui Yoon ( 윤성의 ) Course URL: http://jupiter.kaist.ac.kr/~sungeui/cg Today s Topics Surface representations Smooth curves Subdivision 2 Smooth Curves and
More informationConnected Minimal Acceleration Trigonometric Curves
Connected Minimal Acceleration Trigonometric Curves Tony Barrera Barrera Kristiansen AB Anders Hast Creative Media Lab, University of Gävle Ewert Bengtsson Centre for Image Analysis Uppsala University
More informationThe λ-calculus. 1 Background on Computability. 2 Some Intuition for the λ-calculus. 1.1 Alan Turing. 1.2 Alonzo Church
The λ-calculus 1 Background on Computability The history o computability stretches back a long ways, but we ll start here with German mathematician David Hilbert in the 1920s. Hilbert proposed a grand
More informationImage Reconstruction. Prof. George Wolberg Dept. of Computer Science City College of New York
Image Reconstruction Pro. George Wolberg Dept. o Computer Science Cit College o New Yor Objectives In this lecture we describe image reconstruction: - Interpolation as convolution - Interpolation ernels
More informationParameterization for curve interpolation
Working title: Topics in Multivariate Approximation and Interpolation 101 K. Jetter et al., Editors c 2005 Elsevier B.V. All rights reserved Parameterization for curve interpolation Michael S. Floater
More informationDifferentiability and Tangent Planes October 2013
Differentiability and Tangent Planes 14.4 04 October 2013 Differentiability in one variable. Recall for a function of one variable, f is differentiable at a f (a + h) f (a) lim exists and = f (a) h 0 h
More informationObjects 2: Curves & Splines Christian Miller CS Fall 2011
Objects 2: Curves & Splines Christian Miller CS 354 - Fall 2011 Parametric curves Curves that are defined by an equation and a parameter t Usually t [0, 1], and curve is finite Can be discretized at arbitrary
More informationThe Tangent Line as The Best Linear Approximation
The Tangent Line as The Best Linear Approximation Mark Howell Gonzaga High School Arlington, VA A tangent line to a curve, y f (x), at a point where x a has two important properties: it contains the point
More information4 Visualization and. Approximation
4 Visualization and Approximation b A slope field for the differential equation y tan(x + y) tan(x) tan(y). It is not always possible to write down an explicit formula for the solution to a differential
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 informationLinear Interpolating Splines
Jim Lambers MAT 772 Fall Semester 2010-11 Lecture 17 Notes Tese notes correspond to Sections 112, 11, and 114 in te text Linear Interpolating Splines We ave seen tat ig-degree polynomial interpolation
More information