Interpolation and Deformations A short cookbook
|
|
- Avice Anthony
- 6 years ago
- Views:
Transcription
1 Interpolation and Deformations A short cookbook Fall 2000; Updated: 29 September 205
2 Linear Interpolation p ρ 2 2 = [ ] = 20 T p ρ = [ ] = 5 p 3 = ρ3 =? 0?? T [ ] T Fall 2000; Updated: 29 September 205
3 Linear Interpolation p3 = q3 = a??? + 3 b a p = [ ] T q = a p q [ ] ( ) 2 2 T = = [ ] T b Fall 2000; Updated: 29 September 205
4 Linear Interpolation p 2, A2 λ p, A p = p + λ p p ( ) ( ) 3 2 A 3 =??? A + λ A A λ = 2 ( λ) A λ A 2 ( p3 p) ( p2 p) ( p p ) ( p p ) = Fall 2000; Updated: 29 September 205
5 Bilinear Interpolation u ( λ) u λ i+, j+ i, j+ ( µ ) u i, j λ µ ( λ) u i+, j Fall 2000; Updated: 29 September 205
6 Bilinear Interpolation u u i+, j+ i, j+ ( µ ) µ u i, j λ Fall 2000; Updated: 29 September 205 ( λ) ( µ ) u µ i+, j u( λ, µ ) = λ ( µ u ( µ ) ) u + + ( λ ) µ u + + ( µ ) u = u + λ u u + µ u u + λµ u u ( ) i, j i, j i, j i, j ( ) ( ) ( ) i, j i, j i, j i, j i, j i, j i, j
7 i, j + i, j + Bilinear Interpolation u, A u i+, j +, Ai +, j + u, A i, j i, j u ; A i+, j i+, j u u u u u ( λ, µ ) = interpolate( { λ, µ }{, } i, j, i+, j, i+, j +, i, j + ) ( λ, µ ) = interpolate( { λ, µ }{, } i, j, i+, j, i+, j +, i, j + ) A A A A A Fall 2000; Updated: 29 September 205
8 Let { } { A,, A } 2 -linear Interpolation Λ = λ,, λ, with 0 λ k be a set of interpolation parameters, and let A = be a set of constants. Then we define: linearinterpolate( Λ, A) = ( λ ) A( Λ ) = A( λ,, λ ) linearinterpolate( Λ, A,, A ) { } 2 { A } A + λ linearinterpolate( Λ,,, ) = linearinterpolate( Λ, A) OTE: Sometimes in this situation we will use notation Fall 2000; Updated: 29 September 205
9 Barycentric Interpolation p ( λ) p 2 p λ p p 2 p 2 p p p = λ p + µ p ( λ ) = ( λ ) 2 + λ λ + µ = Fall 2000; Updated: 29 September 205
10 Barycentric Interpolation p 3 A 3 p p 2 A p ( ( λ,, µ ) ) =???? p 3 + λ( ( p p ) 3 ) + µ ( ( p 2 p ) 3 ) = λ p + µ p2 2 + ( ( λ µ ) ) p3 3 p λ, µ, ν = λ p + µ p + ν p where λ + µ + ν = ( ) (,, ) A A λ µ ν = λ A + µ A + νa Fall 2000; Updated: 29 September 205
11 Barycentric Interpolation p 3 A 3 p p 2 A 2 A p p p2 p3 λ = µ ν Fall 2000; Updated: 29 September 205
12 Let Barycentric Interpolation Λ = λ,, λ, with 0 λ and λ = { } { A,, A } 2 k k k = be a set of interpolation parameters, and let A = be a set of constants. Then we define: BarycentricInterpolate( Λ, A) = Λ i A = λ k = k A k OTE: Sometimes in this situation we will use notation A( Λ ) = A( λ,, λ ) = BarycentricInterpolate( Λ, A) Fall 2000; Updated: 29 September 205 OTE: This is a special case of barycentric Bezier st polynomial interpolations (here, degree)
13 Interpolation of functions y( v) 0 v Fall 2000; Updated: 29 September 205
14 Fitting of interpolation curves The discussion below follows (in part) G. Farin, Curves and surfaces for computer-aided geometric design, a practical guide, Academic Press, Boston, 990, chapter 0 and pp Fall 2000; Updated: 29 September 205
15 ote that many forms of polynomial may be used for the P, k,k is the power basis: P ( v ) v a moment 0 ( v ). One common (not very good) choice = k -D Interpolation 0 k, k k = 0 { y v y v } Given set of known values ( ),..., ( ), 0 0 find an approximating polynomial y P( c,..., c ; v ) P( c,..., c ; v ) = c P ( v ) Better choices are the Bernstein plynomials and the b-spline basis functions, which we will discuss in m 0 m v Fall 2000; Updated: 29 September 205
16 -D Interpolation 0 k, k k = 0 { y v y v } Given set of known values ( ),..., ( ), find an approximating polynomial y P( c,..., c ; v ) P( c,..., c ; v ) = c P ( v ) 0 0 m 0 m To do this, solve: P,0( v0 ) P, ( v0 ) c0 y0 " # " " " P,0 ( vm ) P, ( vm ) c ym Fall 2000; Updated: 29 September 205
17 Bezier and Bernstein Polynomials ( ) ( ) = k P c0,, c; v ck v v k = 0 k = k = 0, k ( v ) ( ) k where B, k ( v ) = v v k Excellent numerical stability for 0<v< There exist good ways to convert to more conventional power basis c B k k Fall 2000; Updated: 29 September 205
18 Barycentric Bezier Polynomials ( ) = c k P c 0,,c ;u,v where B,k (u,v) = k =0 k =0 Excellent numerical stability for 0<u,v< There exist good ways to convert to more conventional power basis k = c k B,k (u,v) k u k v k u k v k u+v= Fall 2000; Updated: 29 September 205
19 Bezier Curves " Suppose that the coefficients c are multi-dimensional vectors (e.g., 2D or 3D points). Then the polynomial " " " P c,, c ; v = c B ( v ) ( ) 0 k, k k = 0 computed over the range 0 v curve with control vertices c ". j j generates a Bezier c c 2 c Fall 2000; Updated: 29 September 205 c 3
20 Bezier Curves: de Casteljau Algorithm " Given coefficients c, Bezier curves can be generated recursively by repeated linear interpolation: " " P c c v = b where b (,, ; ) = 0 0 " c 0 j j k k ( ) k j = j + j + b v b vb j c0 v c ( v) v ( v) c2 v ( v) c Fall 2000; Updated: 29 September 205
21 Iterative Form of decasteljau Algorithm Step : b c for 0 j Step 2 : for k step until k = do for j 0 step until j = k do Step 3: return j j b ( v ) b + vb b j j j Fall 2000; Updated: 29 September 205
22 Advantages of Bezier Curves umerically very robust Many nice mathematical properties Smooth Global (may be viewed as a disadvantage) Fall 2000; Updated: 29 September 205
23 B-splines Given coefficient values C = { c 0,", c L+D } "knot points" u = {u 0,",u L+2D 2 } with u i u i+ D = "degree" of desired B-spline Can define an interpolated curve P(C,u; u) on u D u < u L+D Then P(C;u) = L+D c j j =0 D (u) j where j D (u) are B-spline basis polynomials (discussed later) Fall 2000; Updated: 29 September 205
24 B-Spline Polynomials Some useful references include L-09_BSplines_URBS.pdf Fall 2000; Updated: 29 September 205
25 B-spline polynomials & B-spline basis functions Given C, u,d as before where P(C,u;u) = j 0 (u) = L+D c j j =0 u j u u j 0 Otherwise j D (u) j k (u) = u u j k u j +k u j (u) + u u j +k k j u j +k u j + (u) for k > 0 j Fall 2000; Updated: 29 September 205
26 For a B-spline polynomial P(C,u;t) = L+D c j j=0 B-Spline Polynomials j D (u,t) the basis functions j D (u,t) are a function of the degree of the polynomial and the vector u = u 0,",u n of "knot points". The polynomial is "uniform" if the distance between knot points is evenly spaced and "non-uniform" otherwise Fall 2000; Updated: 29 September 205
27 deboor Algorithm (Farin) Given u, c, D as before, can evaluate P(c,u;u) recursively as follows: Step : Determine index i such that u i u < u i+ Step 2: Determine multplicity r such that u i r = u i r + = = u i Step 3: Set " d j 0 = c j for i D + j i + Step 4: Compute P(c,u;u) = d D r i+ recursively, where " d j k = u j+d k u " k dj + u j+d k u j u u j " k dj = α k " k j d j u j+d k u k j γ j + β j k " d j k γ j k Fall 2000; Updated: 29 September 205
28 deboor Algorithm: Example D=3, r=0 d j k = u j +D k u k + u j +D k u j dj 3 0 d i+ = p( { ", ci,"},3;u) 3 2 = α i+ d i + 3 βi+ 2 d i+ ( 2 ) 3 d i+ γ i+ ( ) 2 d i+ γ i+ 2 = α i+ d i + 2 βi+ k u u j α k k j d j dj = + β k k j d j u j +D k u k k j γ j γ j (( ) d 2 i + (u ui ) d 2 ) i+ (u i+ u i ) ( ) d i + (u ui ) d i+ = u i+ u ( ) = u i+2 u (u i+2 u i ) ( 2 + β ) 2 i d i γ i = ( (u i+ u) d i + (u u i ) d ) i (u i+ u i ) ( d 0 i + βi+ d 0 ) ( i+ γ i+ = ( u i+3 u) d 0 i + (u ui ) d 0 ) i+ (u i+3 u i ) ( 0 + β 0 ) i d i γ i = ( (u i+2 u) d 0 i + (u u i ) d 0 ) (u i i+2 u i ) ( d 0 i 2 + β d 0 ) ( i γ i i = (u i+ u) d 0 i 2 + (u u i 2 ) d 0 ) (u i i+ u i 2 ) d i 2 = αi 2 d i d i+ d i = α i+ d i = αi d i = α i Fall 2000; Updated: 29 September 205
29 deboor Algorithm (alternative formula) An alternative formulation from Wikipedia is given as: Given u, c, D as before, can evaluate P(c,u;u) recursively as follows: Step : Determine index i such that u i u < u i+ Step 2: Determine multplicity r such that u i r = u i r + = = u i Step 3: Set " d j 0 = c j for i D + j i + Step 4: Compute P(c,u;u) = d D r i+ recursively, where " d k j = ( α k,j ) d " k " j + α k k,j dj where α k,i = u u j u j+d+ k u j Source: Fall 2000; Updated: 29 September 205
30 Uniform B-Spline Polynomials Third degree uniform B-spline P(C,u;t) = c j 2 j (u,t) with t j = j j 3 j (u,t) = 6 (t j)2 if j t < j ( t j ) 3 + 3( t j ) + 3( t j ) + if j + t < j ( t j ) 3 6( t j ) + 4 if j + 2 t < j ( t j 3 ) if j + 3 t < j otherwise Fall 2000; Updated: 29 September 205
31 Some advantages of B-splines Efficient umerically stable Smooth Local Fall 2000; Updated: 29 September 205
32 2D Interpolation (tensor form) Consider the 2D polynomial m n P( u, v ) = c A ( u) B ( v ) i = 0 j = 0 ij i j where A ( u) and B ( v ) can be arbitrary i c00 c0n B0( v ) = [ A0 ( u),, Am ( u)] " # " " c c B ( v ) j m0 mn n functions (good choices Bernstein polynomials or B-Spline basis functions. Suppose that we have samples ys = y( u, v ) for s = 0,..., s s s We want to find an approximating polynomial P Fall 2000; Updated: 29 September 205
33 2D Interpolation: Finding the best fit Given a set of sample values y s (u s,v s ) corresponding to 2D coordinates (u s,v s ), left hand side basis functions A 0 (u),,a m (u) and right hand side basis functions B 0 (v),,b n (v), the goal is to find the matrix C of coefficients c ij. To do this, solve the least squares problem " y s (u s,v s ) " " A 0 (u s )B 0 (v s ) " " A 0 (u s )B (v s ) " " A i (u s )B j (v s ) " " " A m (u s )B n (v s ) " " c 00 c 0 c ij c mn Fall 2000; Updated: 29 September 205
34 2D Interpolation: Sampling on a regular grid A common special case arises when the ( u, v ) form a regular grid. In this { 0 } { 0 } { 0 } case we have u u,, u and v v,, v. For each value v v,, v solve the row least squares problem j s v s s u " " " X j 0 ys( us, v j ) A0 ( us ) Am ( us ) " " " " X jm for the unknown m-vector X. Then solve m n-variable least squares problems X X X X X X " " " " " # " " " " X X X m m 0 m v v v j s s v B ( v ) B( v0) Bn( v0) c00 c0 cm0 B0( v) B( v) Bn( v) c0 c c m " " " " " # " " " # " c0n c n cmn " " " B0( v ) B( v ) ( ) v B v n v v for the vectors c j 0,, c jn. ote that this latter step requires only SVD or similar matrix computation Fall 2000; Updated: 29 September 205
35 2D Interpolation: Sampling on a regular grid There are a number of caveats to the grid method on the previous slide. (E.g., you need enough data for each of the least squares problems). But where applicable the method can save computation time since it replaces a number of m and n variable least squares problems for one big m x n problem ote that there is a similar trick that you can play by grouping all the common u i elements together. ote that the y s and the c s do not have to be scalar numbers. They can be Vectors, Matrices, or other objects that have appropriate algebraic properties Fall 2000; Updated: 29 September 205
36 -dimensional interpolation The image cannot be displayed. Your computer may not have enough memory to open the image, or the image may have been corrupted. Restart your computer, and then open the file again. If the red x still appears, you may have to delete the image and then insert it again. Define F " ( u) = A ( u ) A ( u ) i i i i Then solve the least squares problem c00 0 # # " " " c 0 0 " F00 0( us ) F0 0( us ) Fm m ( u ) n s y # s # # cm mn Fall 2000; Updated: 29 September 205
37 -dimensional interpolation The methods described earlier generalize naturally to dimensions. " P( u) = P( u,, u ) = c A ( u ) A ( u ) where A K i m m i in i i i = 0 i = 0 ( u) can be arbitrary functions (good choices are Bernstein polynomials or B-Spline basis functions). Suppose that we have samples " y = y ( u ) for s = 0,..., s s s We want to find coefficients of c approximating polynomial P. i in Fall 2000; Updated: 29 September 205
38 Example: 3D Calibration of Distortion Fall 2000; Updated: 29 September 205
39 Example: 3D Calibration of Distortion x y z c 000 c000 c000 " " x y z F 000( us ) # F555 ( us ) p s ps ps x y z c555 c555 c Fall 2000; Updated: 29 September 205
40 Example: 3D Calibration of Distortion The correction function will then look like this: p = CorrectDistortion( q) min max { u = ScaleToBox( q, q, q ) } return c B ( u ) B ( u ) B ( u ) i=0 j=0 k=0 i, j, k 5, i x 5, j y 5, k z Fall 2000; Updated: 29 September 205
Interpolation and Deformations A short cookbook
Interpolation and Deformations A short cookbook 600.445 Fall 2000; Updated: 0 October 2002 Linear Interpolation p ρ 2 2 = [ 40 30 20] = 20 T p ρ = [ 0 5 20] = 5 p 3 = ρ3 =? 0?? T [ 20 20 20] T 2 600.445
More informationLagrange Interpolation and Neville s Algorithm. Ron Goldman Department of Computer Science Rice University
Lagrange Interpolation and Neville s Algorithm Ron Goldman Department of Computer Science Rice University Tension between Mathematics and Engineering 1. How do Mathematicians actually represent curves
More informationThe Essentials of CAGD
The Essentials of CAGD Chapter 4: Bézier Curves: Cubic and Beyond Gerald Farin & Dianne Hansford CRC Press, Taylor & Francis Group, An A K Peters Book www.farinhansford.com/books/essentials-cagd c 2000
More informationInterpolation and polynomial approximation Interpolation
Outline Interpolation and polynomial approximation Interpolation Lagrange Cubic Splines Approximation B-Splines 1 Outline Approximation B-Splines We still focus on curves for the moment. 2 3 Pierre Bézier
More informationCubic Splines; Bézier Curves
Cubic Splines; Bézier Curves 1 Cubic Splines piecewise approximation with cubic polynomials conditions on the coefficients of the splines 2 Bézier Curves computer-aided design and manufacturing MCS 471
More informationArsène Pérard-Gayot (Slides by Piotr Danilewski)
Computer Graphics - Splines - Arsène Pérard-Gayot (Slides by Piotr Danilewski) CURVES Curves Explicit y = f x f: R R γ = x, f x y = 1 x 2 Implicit F x, y = 0 F: R 2 R γ = x, y : F x, y = 0 x 2 + y 2 =
More informationIntroduction to Computer Graphics. Modeling (1) April 13, 2017 Kenshi Takayama
Introduction to Computer Graphics Modeling (1) April 13, 2017 Kenshi Takayama Parametric curves X & Y coordinates defined by parameter t ( time) Example: Cycloid x t = t sin t y t = 1 cos t Tangent (aka.
More informationKatholieke Universiteit Leuven Department of Computer Science
Interpolation with quintic Powell-Sabin splines Hendrik Speleers Report TW 583, January 2011 Katholieke Universiteit Leuven Department of Computer Science Celestijnenlaan 200A B-3001 Heverlee (Belgium)
More informationCurves. Hakan Bilen University of Edinburgh. Computer Graphics Fall Some slides are courtesy of Steve Marschner and Taku Komura
Curves Hakan Bilen University of Edinburgh Computer Graphics Fall 2017 Some slides are courtesy of Steve Marschner and Taku Komura How to create a virtual world? To compose scenes We need to define objects
More informationCurves, Surfaces and Segments, Patches
Curves, Surfaces and Segments, atches The University of Texas at Austin Conics: Curves and Quadrics: Surfaces Implicit form arametric form Rational Bézier Forms and Join Continuity Recursive Subdivision
More informationCS 598: Communication Cost Analysis of Algorithms Lecture 9: The Ideal Cache Model and the Discrete Fourier Transform
CS 598: Communication Cost Analysis of Algorithms Lecture 9: The Ideal Cache Model and the Discrete Fourier Transform Edgar Solomonik University of Illinois at Urbana-Champaign September 21, 2016 Fast
More informationGEOMETRIC MODELLING WITH BETA-FUNCTION B-SPLINES, I: PARAMETRIC CURVES
International Journal of Pure and Applied Mathematics Volume 65 No. 3 2010, 339-360 GEOMETRIC MODELLING WITH BETA-FUNCTION B-SPLINES, I: PARAMETRIC CURVES Arne Lakså 1, Børre Bang 2, Lubomir T. Dechevsky
More informationLinear Algebra (Review) Volker Tresp 2017
Linear Algebra (Review) Volker Tresp 2017 1 Vectors k is a scalar (a number) c is a column vector. Thus in two dimensions, c = ( c1 c 2 ) (Advanced: More precisely, a vector is defined in a vector space.
More informationMAT300/500 Programming Project Spring 2019
MAT300/500 Programming Project Spring 2019 Please submit all project parts on the Moodle page for MAT300 or MAT500. Due dates are listed on the syllabus and the Moodle site. You should include all neccessary
More informationTetrahedral C m Interpolation by Rational Functions
Tetrahedral C m Interpolation by Rational Functions Guoliang Xu State Key Laboratory of Scientific and Engineering Computing, ICMSEC, Chinese Academy of Sciences Chuan I Chu Weimin Xue Department of Mathematics,
More informationBézier Curves and Splines
CS-C3100 Computer Graphics Bézier Curves and Splines Majority of slides from Frédo Durand vectorportal.com CS-C3100 Fall 2017 Lehtinen Before We Begin Anything on your mind concerning Assignment 1? CS-C3100
More informationApproximation of Circular Arcs by Parametric Polynomial Curves
Approximation of Circular Arcs by Parametric Polynomial Curves Gašper Jaklič Jernej Kozak Marjeta Krajnc Emil Žagar September 19, 005 Abstract In this paper the approximation of circular arcs by parametric
More informationNovel polynomial Bernstein bases and Bézier curves based on a general notion of polynomial blossoming
See discussions, stats, and author profiles for this publication at: https://www.researchgate.net/publication/283943635 Novel polynomial Bernstein bases and Bézier curves based on a general notion of polynomial
More informationSample Exam 1 KEY NAME: 1. CS 557 Sample Exam 1 KEY. These are some sample problems taken from exams in previous years. roughly ten questions.
Sample Exam 1 KEY NAME: 1 CS 557 Sample Exam 1 KEY These are some sample problems taken from exams in previous years. roughly ten questions. Your exam will have 1. (0 points) Circle T or T T Any curve
More informationKasetsart University Workshop. Multigrid methods: An introduction
Kasetsart University Workshop Multigrid methods: An introduction Dr. Anand Pardhanani Mathematics Department Earlham College Richmond, Indiana USA pardhan@earlham.edu A copy of these slides is available
More informationSolutions of the Harmonic Oscillator Equation in a B-Polynomial Basis
Adv. Studies Theor. Phys., Vol. 3, 2009, no. 12, 451-460 Solutions of the Harmonic Oscillator Equation in a B-Polynomial Basis Muhammad I. Bhatti Department of Physics and Geology University of Texas-Pan
More information1 Computing the Invariant Distribution
LECTURE NOTES ECO 613/614 FALL 2007 KAREN A. KOPECKY 1 Computing the Invariant Distribution Suppose an individual s state consists of his current assets holdings, a [a, ā] and his current productivity
More informationRational Bézier Patch Differentiation using the. Rational Forward Difference Operator
Rational Bézier Patch Differentiation using the Rational Forward Difference Operator Xianming Chen Richard Riesenfeld Elaine Cohen School of Computing University of Utah CGI 05 p.1/37 Introduction For
More informationElementary maths for GMT
Elementary maths for GMT Linear Algebra Part 2: Matrices, Elimination and Determinant m n matrices The system of m linear equations in n variables x 1, x 2,, x n a 11 x 1 + a 12 x 2 + + a 1n x n = b 1
More informationCHAPTER 3 Further properties of splines and B-splines
CHAPTER 3 Further properties of splines and B-splines In Chapter 2 we established some of the most elementary properties of B-splines. In this chapter our focus is on the question What kind of functions
More information1 f(u)b n i (u)du, (1.1)
AMRX Applied Mathematics Research express 25, No. 4 On a Modified Durrmeyer-Bernstein Operator and Applications Germain E. Randriambelosoa 1 Introduction Durrmeyer [6] has introduced a Bernstein-type operator
More informationNotes on Cellwise Data Interpolation for Visualization Xavier Tricoche
Notes on Cellwise Data Interpolation for Visualization Xavier Tricoche urdue University While the data (computed or measured) used in visualization is only available in discrete form, it typically corresponds
More informationApproximation of Circular Arcs by Parametric Polynomials
Approximation of Circular Arcs by Parametric Polynomials Emil Žagar Lecture on Geometric Modelling at Charles University in Prague December 6th 2017 1 / 44 Outline Introduction Standard Reprezentations
More information6. Iterative Methods for Linear Systems. The stepwise approach to the solution...
6 Iterative Methods for Linear Systems The stepwise approach to the solution Miriam Mehl: 6 Iterative Methods for Linear Systems The stepwise approach to the solution, January 18, 2013 1 61 Large Sparse
More informationJournal of Computational and Applied Mathematics
Journal of Computational and Applied Mathematics 233 2010) 1554 1561 Contents lists available at ScienceDirect Journal of Computational and Applied Mathematics journal homepage: wwwelseviercom/locate/cam
More informationPractical Linear Algebra: A Geometry Toolbox
Practical Linear Algebra: A Geometry Toolbox Third edition Chapter 12: Gauss for Linear Systems Gerald Farin & Dianne Hansford CRC Press, Taylor & Francis Group, An A K Peters Book www.farinhansford.com/books/pla
More informationM2R IVR, October 12th Mathematical tools 1 - Session 2
Mathematical tools 1 Session 2 Franck HÉTROY M2R IVR, October 12th 2006 First session reminder Basic definitions Motivation: interpolate or approximate an ordered list of 2D points P i n Definition: spline
More informationComputer Aided Design. B-Splines
1 Three useful references : R. Bartels, J.C. Beatty, B. A. Barsky, An introduction to Splines for use in Computer Graphics and Geometric Modeling, Morgan Kaufmann Publications,1987 JC.Léon, Modélisation
More informationLecture 20: Bezier Curves & Splines
Lecture 20: Bezier Curves & Splines December 6, 2016 12/6/16 CSU CS410 Bruce Draper & J. Ross Beveridge 1 Review: The Pen Metaphore Think of putting a pen to paper Pen position described by time t Seeing
More informationCHAPTER 10 Shape Preserving Properties of B-splines
CHAPTER 10 Shape Preserving Properties of B-splines In earlier chapters we have seen a number of examples of the close relationship between a spline function and its B-spline coefficients This is especially
More informationLeast Squares Approximation
Chapter 6 Least Squares Approximation As we saw in Chapter 5 we can interpret radial basis function interpolation as a constrained optimization problem. We now take this point of view again, but start
More information2 The De Casteljau algorithm revisited
A new geometric algorithm to generate spline curves Rui C. Rodrigues Departamento de Física e Matemática Instituto Superior de Engenharia 3030-199 Coimbra, Portugal ruicr@isec.pt F. Silva Leite Departamento
More informationDeterminants. Chia-Ping Chen. Linear Algebra. Professor Department of Computer Science and Engineering National Sun Yat-sen University 1/40
1/40 Determinants Chia-Ping Chen Professor Department of Computer Science and Engineering National Sun Yat-sen University Linear Algebra About Determinant A scalar function on the set of square matrices
More informationFDTD for 1D wave equation. Equation: 2 H Notations: o o. discretization. ( t) ( x) i i i i i
FDTD for 1D wave equation Equation: 2 H = t 2 c2 2 H x 2 Notations: o t = nδδ, x = iδx o n H nδδ, iδx = H i o n E nδδ, iδx = E i discretization H 2H + H H 2H + H n+ 1 n n 1 n n n i i i 2 i+ 1 i i 1 = c
More informationJim Lambers MAT 610 Summer Session Lecture 1 Notes
Jim Lambers MAT 60 Summer Session 2009-0 Lecture Notes Introduction This course is about numerical linear algebra, which is the study of the approximate solution of fundamental problems from linear algebra
More informationKevin James. MTHSC 3110 Section 2.1 Matrix Operations
MTHSC 3110 Section 2.1 Matrix Operations Notation Let A be an m n matrix, that is, m rows and n columns. We ll refer to the entries of A by their row and column indices. The entry in the i th row and j
More informationNonlinear Means in Geometric Modeling
Nonlinear Means in Geometric Modeling Michael S. Floater SINTEF P. O. Box 124 Blindern, 0314 Oslo, Norway E-mail: Michael.Floater@math.sintef.no Charles A. Micchelli IBM Corporation T.J. Watson Research
More informationCONTROL POLYGONS FOR CUBIC CURVES
On-Line Geometric Modeling Notes CONTROL POLYGONS FOR CUBIC CURVES Kenneth I. Joy Visualization and Graphics Research Group Department of Computer Science University of California, Davis Overview B-Spline
More informationMatrices: 2.1 Operations with Matrices
Goals In this chapter and section we study matrix operations: Define matrix addition Define multiplication of matrix by a scalar, to be called scalar multiplication. Define multiplication of two matrices,
More informationApplied Math 205. Full office hour schedule:
Applied Math 205 Full office hour schedule: Rui: Monday 3pm 4:30pm in the IACS lounge Martin: Monday 4:30pm 6pm in the IACS lounge Chris: Tuesday 1pm 3pm in Pierce Hall, Room 305 Nao: Tuesday 3pm 4:30pm
More informationLinear Algebra for Machine Learning. Sargur N. Srihari
Linear Algebra for Machine Learning Sargur N. srihari@cedar.buffalo.edu 1 Overview Linear Algebra is based on continuous math rather than discrete math Computer scientists have little experience with it
More informationn n matrices The system of m linear equations in n variables x 1, x 2,..., x n can be written as a matrix equation by Ax = b, or in full
n n matrices Matrices Definitions Diagonal, Identity, and zero matrices Addition Multiplication Transpose and inverse The system of m linear equations in n variables x 1, x 2,..., x n a 11 x 1 + a 12 x
More informationQuaternion Cubic Spline
Quaternion Cubic Spline James McEnnan jmcennan@mailaps.org May 28, 23 1. INTRODUCTION A quaternion spline is an interpolation which matches quaternion values at specified times such that the quaternion
More informationLecture 1 INF-MAT : Chapter 2. Examples of Linear Systems
Lecture 1 INF-MAT 4350 2010: Chapter 2. Examples of Linear Systems Tom Lyche Centre of Mathematics for Applications, Department of Informatics, University of Oslo August 26, 2010 Notation The set of natural
More informationFunctional SVD for Big Data
Functional SVD for Big Data Pan Chao April 23, 2014 Pan Chao Functional SVD for Big Data April 23, 2014 1 / 24 Outline 1 One-Way Functional SVD a) Interpretation b) Robustness c) CV/GCV 2 Two-Way Problem
More informationCMSC427 Parametric curves: Hermite, Catmull-Rom, Bezier
CMSC427 Parametric curves: Hermite, Catmull-Rom, Bezier Modeling Creating 3D objects How to construct complicated surfaces? Goal Specify objects with few control points Resulting object should be visually
More informationINF4130: Dynamic Programming September 2, 2014 DRAFT version
INF4130: Dynamic Programming September 2, 2014 DRAFT version In the textbook: Ch. 9, and Section 20.5 Chapter 9 can also be found at the home page for INF4130 These slides were originally made by Petter
More informationLinear Algebra Review. Vectors
Linear Algebra Review 9/4/7 Linear Algebra Review By Tim K. Marks UCSD Borrows heavily from: Jana Kosecka http://cs.gmu.edu/~kosecka/cs682.html Virginia de Sa (UCSD) Cogsci 8F Linear Algebra review Vectors
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 informationInterpolation and polynomial approximation Interpolation
Outline Interpolation and polynomial approximation Interpolation Lagrange Cubic Approximation Bézier curves B- 1 Some vocabulary (again ;) Control point : Geometric point that serves as support to the
More informationIn practice one often meets a situation where the function of interest, f(x), is only represented by a discrete set of tabulated points,
1 Interpolation 11 Introduction In practice one often meets a situation where the function of interest, f(x), is only represented by a discrete set of tabulated points, {x i, y i = f(x i ) i = 1 n, obtained,
More informationWilliam Stallings Copyright 2010
A PPENDIX E B ASIC C ONCEPTS FROM L INEAR A LGEBRA William Stallings Copyright 2010 E.1 OPERATIONS ON VECTORS AND MATRICES...2 Arithmetic...2 Determinants...4 Inverse of a Matrix...5 E.2 LINEAR ALGEBRA
More informationUMCS. Fast multidimensional Bernstein-Lagrange algorithms. The John Paul II Catholic University of Lublin, ul. Konstantynow 1H, Lublin, Poland
Annales Informatica AI XII, 2 2012) 7 18 DOI: 10.2478/v10065-012-0002-6 Fast multidimensional Bernstein-Lagrange algorithms Joanna Kapusta 1, Ryszard Smarzewski 1 1 Institute of Mathematics and Computer
More informationComputergrafik. Matthias Zwicker Universität Bern Herbst 2016
Computergrafik Matthias Zwicker Universität Bern Herbst 2016 2 Today Curves Introduction Polynomial curves Bézier curves Drawing Bézier curves Piecewise curves Modeling Creating 3D objects How to construct
More informationCHANGE OF BASIS AND ALL OF THAT
CHANGE OF BASIS AND ALL OF THAT LANCE D DRAGER Introduction The goal of these notes is to provide an apparatus for dealing with change of basis in vector spaces, matrices of linear transformations, and
More informationConvergence under Subdivision and Complexity of Polynomial Minimization in the Simplicial Bernstein Basis
Convergence under Subdivision and Complexity of Polynomial Minimization in the Simplicial Bernstein Basis Richard Leroy IRMAR, University of Rennes 1 richard.leroy@ens-cachan.org Abstract In this article,
More informationMATLAB for Engineers
MATLAB for Engineers Adrian Biran Moshe Breiner ADDISON-WESLEY PUBLISHING COMPANY Wokingham, England Reading, Massachusetts Menlo Park, California New York Don Mills, Ontario Amsterdam Bonn Sydney Singapore
More informationGeometric meanings of the parameters on rational conic segments
Science in China Ser. A Mathematics 005 Vol.48 No.9 09 09 Geometric meanings of the parameters on rational conic segments HU Qianqian & WANG Guojin Department of Mathematics, Zhejiang University, Hangzhou
More informationReading. w Foley, Section 11.2 Optional
Parametric Curves w Foley, Section.2 Optional Reading w Bartels, Beatty, and Barsky. An Introduction to Splines for use in Computer Graphics and Geometric Modeling, 987. w Farin. Curves and Surfaces for
More information4. Multi-linear algebra (MLA) with Kronecker-product data.
ect. 3. Tensor-product interpolation. Introduction to MLA. B. Khoromskij, Leipzig 2007(L3) 1 Contents of Lecture 3 1. Best polynomial approximation. 2. Error bound for tensor-product interpolants. - Polynomial
More informationMTH5112 Linear Algebra I MTH5212 Applied Linear Algebra (2017/2018)
MTH5112 Linear Algebra I MTH5212 Applied Linear Algebra (2017/2018) COURSEWORK 3 SOLUTIONS Exercise ( ) 1. (a) Write A = (a ij ) n n and B = (b ij ) n n. Since A and B are diagonal, we have a ij = 0 and
More informationMA 323 Geometric Modelling Course Notes: Day 12 de Casteljau s Algorithm and Subdivision
MA 323 Geometric Modelling Course Notes: Day 12 de Casteljau s Algorithm and Subdivision David L. Finn Yesterday, we introduced barycentric coordinates and de Casteljau s algorithm. Today, we want to go
More informationAn O(h 2n ) Hermite approximation for conic sections
An O(h 2n ) Hermite approximation for conic sections Michael Floater SINTEF P.O. Box 124, Blindern 0314 Oslo, NORWAY November 1994, Revised March 1996 Abstract. Given a segment of a conic section in the
More informationChapter 2 Interpolation
Chapter 2 Interpolation Experiments usually produce a discrete set of data points (x i, f i ) which represent the value of a function f (x) for a finite set of arguments {x 0...x n }. If additional data
More informationCSE 167: Lecture 11: Bézier Curves. Jürgen P. Schulze, Ph.D. University of California, San Diego Fall Quarter 2012
CSE 167: Introduction to Computer Graphics Lecture 11: Bézier Curves Jürgen P. Schulze, Ph.D. University of California, San Diego Fall Quarter 2012 Announcements Homework project #5 due Nov. 9 th at 1:30pm
More informationApril 26, Applied mathematics PhD candidate, physics MA UC Berkeley. Lecture 4/26/2013. Jed Duersch. Spd matrices. Cholesky decomposition
Applied mathematics PhD candidate, physics MA UC Berkeley April 26, 2013 UCB 1/19 Symmetric positive-definite I Definition A symmetric matrix A R n n is positive definite iff x T Ax > 0 holds x 0 R n.
More informationSignatures of GL n Multiplicity Spaces
Signatures of GL n Multiplicity Spaces UROP+ Final Paper, Summer 2016 Mrudul Thatte Mentor: Siddharth Venkatesh Project suggested by Pavel Etingof September 1, 2016 Abstract A stable sequence of GL n representations
More informationUniversity of Houston, Department of Mathematics Numerical Analysis, Fall 2005
4 Interpolation 4.1 Polynomial interpolation Problem: LetP n (I), n ln, I := [a,b] lr, be the linear space of polynomials of degree n on I, P n (I) := { p n : I lr p n (x) = n i=0 a i x i, a i lr, 0 i
More informationLock-Free Approaches to Parallelizing Stochastic Gradient Descent
Lock-Free Approaches to Parallelizing Stochastic Gradient Descent Benjamin Recht Department of Computer Sciences University of Wisconsin-Madison with Feng iu Christopher Ré Stephen Wright minimize x f(x)
More informationComputing the Hausdorff Distance between Two B-Spline Curves. Zachi Shtain
Computing the Hausdorff Distance between Two B-Spline Curves Zachi Shtain Based on the work of: Chen et al. 2010 Definition Given two curves C 1, C 2, their Hausdorff distance is defined as: H Where:,
More informationMatrices Gaussian elimination Determinants. Graphics 2009/2010, period 1. Lecture 4: matrices
Graphics 2009/2010, period 1 Lecture 4 Matrices m n matrices Matrices Definitions Diagonal, Identity, and zero matrices Addition Multiplication Transpose and inverse The system of m linear equations in
More informationRobust Stability. Robust stability against time-invariant and time-varying uncertainties. Parameter dependent Lyapunov functions
Robust Stability Robust stability against time-invariant and time-varying uncertainties Parameter dependent Lyapunov functions Semi-infinite LMI problems From nominal to robust performance 1/24 Time-Invariant
More informationMultigrid Methods and their application in CFD
Multigrid Methods and their application in CFD Michael Wurst TU München 16.06.2009 1 Multigrid Methods Definition Multigrid (MG) methods in numerical analysis are a group of algorithms for solving differential
More informationReview of linear algebra
Review of linear algebra 1 Vectors and matrices We will just touch very briefly on certain aspects of linear algebra, most of which should be familiar. Recall that we deal with vectors, i.e. elements of
More informationResearch Article On the Derivatives of Bernstein Polynomials: An Application for the Solution of High Even-Order Differential Equations
Hindawi Publishing Corporation Boundary Value Problems Volume 211, Article ID 829543, 16 pages doi:1.1155/211/829543 Research Article On the Derivatives of Bernstein Polynomials: An Application for the
More informationThe Kernel Trick, Gram Matrices, and Feature Extraction. CS6787 Lecture 4 Fall 2017
The Kernel Trick, Gram Matrices, and Feature Extraction CS6787 Lecture 4 Fall 2017 Momentum for Principle Component Analysis CS6787 Lecture 3.1 Fall 2017 Principle Component Analysis Setting: find the
More informationREFRESHER. William Stallings
BASIC MATH REFRESHER William Stallings Trigonometric Identities...2 Logarithms and Exponentials...4 Log Scales...5 Vectors, Matrices, and Determinants...7 Arithmetic...7 Determinants...8 Inverse of a Matrix...9
More informationA B-SPLINE-LIKE BASIS FOR THE POWELL-SABIN 12-SPLIT BASED ON SIMPLEX SPLINES
MATHEMATICS OF COMPUTATION Volume, Number, Pages S 5-578(XX)- A B-SPLINE-LIKE BASIS FOR THE POWELL-SABIN -SPLIT BASED ON SIMPLEX SPLINES ELAINE COHEN, TOM LYCHE, AND RICHARD F. RIESENFELD Abstract. We
More informationM3: Multiple View Geometry
M3: Multiple View Geometry L18: Projective Structure from Motion: Iterative Algorithm based on Factorization Based on Sections 13.4 C. V. Jawahar jawahar-at-iiit.net Mar 2005: 1 Review: Reconstruction
More informationElements of Continuum Elasticity. David M. Parks Mechanics and Materials II February 25, 2004
Elements of Continuum Elasticity David M. Parks Mechanics and Materials II 2.002 February 25, 2004 Solid Mechanics in 3 Dimensions: stress/equilibrium, strain/displacement, and intro to linear elastic
More informationMatrices. Chapter Definitions and Notations
Chapter 3 Matrices 3. Definitions and Notations Matrices are yet another mathematical object. Learning about matrices means learning what they are, how they are represented, the types of operations which
More informationConstruction of Smooth Fractal Surfaces Using Hermite Fractal Interpolation Functions. P. Bouboulis, L. Dalla and M. Kostaki-Kosta
BULLETIN OF THE GREEK MATHEMATICAL SOCIETY Volume 54 27 (79 95) Construction of Smooth Fractal Surfaces Using Hermite Fractal Interpolation Functions P. Bouboulis L. Dalla and M. Kostaki-Kosta Received
More informationMath 4377/6308 Advanced Linear Algebra
2.3 Composition Math 4377/6308 Advanced Linear Algebra 2.3 Composition of Linear Transformations Jiwen He Department of Mathematics, University of Houston jiwenhe@math.uh.edu math.uh.edu/ jiwenhe/math4377
More informationHessenberg eigenvalue eigenvector relations and their application to the error analysis of finite precision Krylov subspace methods
Hessenberg eigenvalue eigenvector relations and their application to the error analysis of finite precision Krylov subspace methods Jens Peter M. Zemke Minisymposium on Numerical Linear Algebra Technical
More informationNumerical Methods I: Interpolation (cont ed)
1/20 Numerical Methods I: Interpolation (cont ed) Georg Stadler Courant Institute, NYU stadler@cims.nyu.edu November 30, 2017 Interpolation Things you should know 2/20 I Lagrange vs. Hermite interpolation
More informationB-Spline Interpolation on Lattices
B-Spline Interpolation on Lattices David Eberly, Geometric Tools, Redmond WA 98052 https://www.geometrictools.com/ This work is licensed under the Creative Commons Attribution 4.0 International License.
More informationMath 671: Tensor Train decomposition methods
Math 671: Eduardo Corona 1 1 University of Michigan at Ann Arbor December 8, 2016 Table of Contents 1 Preliminaries and goal 2 Unfolding matrices for tensorized arrays The Tensor Train decomposition 3
More informationFast Matrix Product Algorithms: From Theory To Practice
Introduction and Definitions The τ-theorem Pan s aggregation tables and the τ-theorem Software Implementation Conclusion Fast Matrix Product Algorithms: From Theory To Practice Thomas Sibut-Pinote Inria,
More informationComputer Graphics Keyframing and Interpola8on
Computer Graphics Keyframing and Interpola8on This Lecture Keyframing and Interpola2on two topics you are already familiar with from your Blender modeling and anima2on of a robot arm Interpola2on linear
More informationRecurrences COMP 215
Recurrences COMP 215 Analysis of Iterative Algorithms //return the location of the item matching x, or 0 if //no such item is found. index SequentialSearch(keytype[] S, in, keytype x) { index location
More informationMATH 320, WEEK 7: Matrices, Matrix Operations
MATH 320, WEEK 7: Matrices, Matrix Operations 1 Matrices We have introduced ourselves to the notion of the grid-like coefficient matrix as a short-hand coefficient place-keeper for performing Gaussian
More informationGeometric Modeling Summer Semester 2010 Mathematical Tools (1)
Geometric Modeling Summer Semester 2010 Mathematical Tools (1) Recap: Linear Algebra Today... Topics: Mathematical Background Linear algebra Analysis & differential geometry Numerical techniques Geometric
More informationChapter 4: Interpolation and Approximation. October 28, 2005
Chapter 4: Interpolation and Approximation October 28, 2005 Outline 1 2.4 Linear Interpolation 2 4.1 Lagrange Interpolation 3 4.2 Newton Interpolation and Divided Differences 4 4.3 Interpolation Error
More informationScientific Computing WS 2017/2018. Lecture 18. Jürgen Fuhrmann Lecture 18 Slide 1
Scientific Computing WS 2017/2018 Lecture 18 Jürgen Fuhrmann juergen.fuhrmann@wias-berlin.de Lecture 18 Slide 1 Lecture 18 Slide 2 Weak formulation of homogeneous Dirichlet problem Search u H0 1 (Ω) (here,
More informationBernstein polynomials of degree N are defined by
SEC. 5.5 BÉZIER CURVES 309 5.5 Bézier Curves Pierre Bézier at Renault and Paul de Casteljau at Citroën independently developed the Bézier curve for CAD/CAM operations, in the 1970s. These parametrically
More information