7. Piecewise Polynomial (Spline) Interpolation
|
|
- Hillary Lang
- 6 years ago
- Views:
Transcription
1 Piecewise Polynomial (Spline Interpolation Single polynomial interpolation has two major disadvantages First, it is not computationally efficient when the number of data points is large When the number of data points is large, degree of the interpolating polynomial becomes large too and, therefore, it takes much more time to compute the value of the interpolating polynomial at an intermediate point Second, it generates relatively large numerical error when the number of data points is large because the number of arithmetic operations involved in the computation of a high degree polynomial is large too We will present an interpolation technique that overcomes the problems of simple polynomial interpolation in this chapter Given a set of m points (x 1, f 1, (x, f,, (x m, f m, and (x m, f m with a x 1 < x < < x m < x m b, the idea is to construct a function g (x on [a,b ] such that g ( f i, i 1,,, m (71 and g (x (x, x [ ]; i 1,,, m (7 where (x are polynomials (of low degree In general, (x are different and they are usually chosen to be of degree three g (x constructed this way is called a piecewise polynomial and is a continuous function on [a, b ] (since ( g ( (, i 1,,, m 1 Advantages of piecewise polynomial interpolation are disadvantages of single polynomial interpolation Two types of piecewise polynomial interpolation techniques will be discussed 71 Hermite Interpolation For each data point (, f i, if the derivative of the unknown function f at, f (, is also given, and the piecewise interpolant is required to interpolate the derivatives as well, then the resulting piecewise interpolant is called a piecewise Hermite polynomial This means that for each interval [ ], (x needs to satisfy four conditions: If (x is a cubic polynomial, ie, ( f i ; ( f i ; ( f i ; ( f i (73 (a i,, c i, d i to be determined then we have (x a i + x + c i x + d i x 3
2 a i + + c i + d i 3 f i a i + + c i + c i + c i 3 + d i + 3d i + 3d i f i f i f i (74 or a i c i d i f i f i f i f i (75 Therefore by solving this system, say using Gaussian elimination with partial pivoting, for a i,, c i, and d i, we get the polynomial (x An algorithm can be given as follows Algorithm: Piecewise Cubic Hermite Interpolation (PCHI Input: (, f i, f i, i 1,,,m, a x 1 < x < < x m < x m b, m : number of segments x : a number between a and b Output: the value of the piecewise cubic Hermite polynomial at x var: C : array[1m,14] of real; {to store coefficients of all } k, i : integer; 1 [Construct the piecewise cubic Hermite polynomial] for i :1 to m do begin Gauss(, f i, f i, f i, f i, a i,, c i, d i ; D [i, 1] : a i ; D [i, ] : ; D [i, 3] : c i ; D [i, 4] : d i ;
3 end; [Compute the value of the piecewise polynomial at x] k : 1; while (x > x k do k : k ; {Find the interval I i [ ] that contains x } b : D [k,4]; {b d k } for i : 3 downto 1 do b : D [k,x ] + x*b ; 3 Return b p k (x In step 1 of the above algorithm, Gauss(, f i, f i, f i, f i, a i,, c i, d i is a procedure that solves the equation in (75 for a i,, c i, and d i using Gaussian elimination with partial pivoting (x, however, can be computed in a more efficient way Note that if h 1 (t, h (t, h 3 (t and h 4 (t are defined as follows: h 1 (t 1 3t + t 3 ; h (t 3t t 3 ; h 3 (t t t + t 3 ; t [0, 1] (76 h 4 (t t + t 3 ; h 1 (t h (t 0 1 Figure 71 Graphs of h 1 (t and h (t then it is easy to see that the value of h 1 (t at t 0 is 1 and its value at t 1 and its derivatives at t 0 and t 1 are all equal to 0, while the value of h (t at t 1 is 1 and its value at t 0 and its derivatives at t 0 and t 1 are all equal to 0 (see Figure 71 for the graphs of h 1 (t and h (t On the other hand, the derivative of h 3 (t at t 0 is 1 and its derivative at t 1 and its values at
4 h 3 (t h 4 (t 0 1 Figure 7 Graphs of h 3 (t and h 4 (t t 0 and t 1 are all equal to 0, while the derivative of h 4 (t at t 1 is 1 and its derivative at t 0 and its values at t 0 and t 1 are all equal to 0 (see Figure 7 for the graphs of h 3 (t and h 4 (t Namely, h 1 (0 1; h 1 (1 0; h 1 (0 0; h 1 (1 0; h (0 0; h (1 1; h (0 0; h (1 0; h 3 (0 0; h 3 (1 0; h 3 (0 1; h 3 (1 0; h 4 (0 0; h 4 (1 0; h 4 (0 0; h 4 (1 1; Hence, if the values and derivatives of a function g (t are given at t 0 and t 1, and denoted by A, B, C, and D, respectively, ie, A g (0, B g (1, C g (0 and D g (1, then a polynomial of degree three which interpolates g (t and g (t at t 0 and t 1 can be defined as follows: p (t A h 1 (t + B h (t + C h 3 (t + D h 4 (t (77 It is easy to see that p (0 A, p (1 B, p (0 C and p (1 D h 1 (t and h (t are called Hermite basis functions of the first kind and h 3 (t and h 4 (t are called Hermite basis functions of the second kind To construct the polynomial (x which satisfies the conditions in (73, one first defines H 1 (x, H (x, H 3 (x and H 4 (x as follows: x H 1 (x h 1 ( x 1 3( x + ( 3 ; x x x 3 ; H (x h ( x H 3 (x h 3 ( 3( x 3( x ( x + ( 3 ; (78 x H 4 (x h 4 ( x ( x + ( 3 ;
5 where x [ ] H 1 (x, H (x, H 3 (x and H 4 (x satisfy the following conditions: H 1 ( 1; H 1 ( 0; H 1 ( 0; H 1 ( 0; H ( 0; H ( 1; H ( 0; H ( 0; H 3 ( 0; H 3 ( 0; H 3 ( 1; H 3 ( 0; H 4 ( 0; H 4 ( 0; H 4 ( 0; H 4 ( 1; Namely, H 1 (x, H (x, H 3 (x and H 4 (x are Hermite basis functions defined on the interval [ ] Therefore, following the p (t defined in (77, one can define the (x which satisfies the conditions in (73 as follows: (x f i H 1 (x + f i H (x + f i H 3 (x + f i H 4 (x (79 It is easy to see that ( f i, ( f i, ( f i and ( f i 7 Spline Interpolation A piecewise cubic interpolating polynomial g (x is called a cubic spline if it has continuous first and second derivatives The process of determining such a piecewise polynomial is called cubic spline interpolation Given a set of n points (x 1, f 1, (x, f,, (x m, f m, and (x m, f m with a x 1 < x < < x m < x m b, a cubic spline g (x that interpolates these points must satisfy the following conditions: (1 g ( f i, i 1,,, m ; ( for each subinterval [ ], i 1,,, m, g (x is equal to some cubic polynomial (x for x [ ]; (3 for each, i, 3,, m, the left derivative of g at equals the right derivative of g at, ie, g ( g ( +, and the left second derivative of g at equals the right second derivative of g at, ie, g ( g ( + Condition (3 is equivalent to the following conditions: 1 ( (, 1 ( ( ; i, 3,, m (710 For each i, 3,, m, by using ( 1, f i 1, ( 1, f i 1, (, f i and (, f i as interpolation points for the interval [ 1 ], we can put 1 (x in Newton s form as follows: 1 (x f [ 1 ] + f [ 1 1 ](x 1 + f [ 1 1 ](x 1 (x 1 + f [ 1 1 ](x 1 (x 1 (x (711 By (418 we know that if g (x is a (continuous function on the interval [a, b ] then the second
6 divided difference of g (x, g [a, b ], can be expressed as follows If g (x is differentiable at x a then we have This follows from the fact that g [b ] g [a ] g (b g (a g [a, b ] b a b a g [a, a ] g (a (71 g [a, a ] lim g [a, b ] lim b a b a g (b g (a b a g (a Therefore by using (418, (71 and the fact that f [ 1 ] f ( 1 f i 1, (711 can be expressed as 1 (x f i 1 + f ( 1 (x f [ 1 ] f [ 1 1 ] 1 f [ 1 ] f ( 1 1 (x 1 (x 1 [(x 1 ( 1 ] f i 1 + f ( 1 (x 1 + f [ 1 ] f ( 1 1 (x 1 + (f [ 1 ] f [ 1 1 ] (x 1 + f [ 1 ] f [ 1 1 ] 1 (x 1 3 By combining like terms and using the notation, 1 (x can be expressed as 1 (x a i (x 1 + c i 1 (x 1 + d i 1 (x 1 3 ; i, 3,,m (71 for where a i 1 f i 1 (713 1 f ( 1 (714 c i 1 f [ 1 ] 1 d i 1 (715
7 d i 1 f [ 1 ] f [ 1 1 ] f [ 1 ] + 1 ( (716 Note that in (71, 1 are the only unknowns; c i 1 and d i 1 are computable once we have for all i 1,,, m, m To find a solution for (71, recall that g (x must have continuous second derivative at, i, 3,, m, which is equivalent to the the second part of (710, ie, Since 1 ( ( ; i, 3,, m (717 and 1 ( c i 1 + 6d i 1 ( c i (717 is equivalent to c i 1 + 3d i 1 c i ; i, 3,, m (718 Using the definition of c k in (715, (718 can be written as f [ 1 ] 1 d i 1 + 3d i 1 f [ ] d i By combining like terms and using the definition of d k in (716, we have f [ 1 ] 1 f [ 1 ] ( f [ ] f [ ] + ( or Therefore, we have f [ 1 ] f [ 1 ] f [ ] + f [ ] 1 + ( + + 3(f [ 1 ] + f [ ] ; i, 3,, m (719
8 where, i 1,,, m, are to be determined (Note that f ( But we don t know the values of f ( (719 can also be put in in matrix form, as follows x ( x 1 + x x 3 x 1 ( x + x 3 x m 1 x ( x m 1 + x m x m x m ( x m + x m 1 x m 1 b 1 b b ṁ b m e e 3 e m 1 e m (70 where e i 3(f [ 1 ] + f [ ] ; i, 3,, m (71 With m unknows and m 1 equations, this is an underdetermined system We need two extra conditions (equations to solve this system We will consider two cases here Case 1: f 1 and f m are given In this case, the system of equations in (70 becomes of the following form: ( x 1 + x x 3 x 1 ( x + x 3 x m 1 x ( x m 1 + x m x m x m ( x m + x m 1 b b 3 b m 1 b m e x b 1 e 3 e m 1 e m x m 1 b m (7 where e i, i, 3,, m, are defined in (71 This is a system of m 1 equations in m 1 unknowns and it is strictly diagonally dominant Therefore, (7 is solvable The process of constructing the spline function g (x can be done in four steps (1 Construct e i, i, 3,, m, using the following foumula: e i 3(f [ 1 ] + f [ ] f [ ] f [ 1 ] f [x 3( x i ] f [ ] i + x i 1 1 f i f 3( x i 1 i f i f + x i i 1
9 ( Solve (7 using a special form of Gaussian elimination without partial pivoting for b, b 3,, b m (3 Use (716 to compute d i for i 1,,, m d i f [ 1 ] + 1 ( f i f i ( 1 (4 Use (715 to compute c i for i 1,,, m, and them form (x for i 1,,, m An algorithm can ge given as follows For m given points (x 1, f 1, (x, f,, (x m, f m, and the tangents b 1 and b m at x 1 and x m, this algorithm constructs a cubic spline that interpolates (, f i, i 1,,, m, b 1, and b m, and computes the value of the cubic spline at x z Algorithm: CUBIC_SPLINE_1 Input: m : number of interpolation points (, f i, i 1,,, m, a x 1 < x < < x m < x m b b 1, b m : tangents at x 1 and x m z : a point different from all Output: the value of the cubic spline at z var: X, {array to store all } F, {array to store all f i } B : array[1m] of real; {array to store all } C, {array to store all c i } D, {array to store all d i } : array[1m] of real; {array to store all } 1 [Construct the cubic spline] for i :1 to m do begin X [i ]: ; F [i ]: f i ; end; B [1]: b 1 ; B [m ]: b m ; SPECIAL_GAUSS(X, F, B ;
10 {This procedure solves (7 for, i, 3,, m, and returns these values in B [i ], i, 3,, m } for i :1 to m do [i ]: X [i ] X [i ]; for i :1 to m do D [i ]: (B [i ] * (F [i ] F [i ]/ [i ] + B [i ]/( [i ] ; for i :1 to m do C [i ]: ((F [i ] F [i ]/ [i ] B [i ]/ [i ] [i ]*D [i ] ; [Compute the value of the cubic spline at x z ] k : 1; while (z > X [k ] do k : k ; {Find the interval I k [x k, x k ] that contains z } g : F [k ] + B [k ]* (z X [k ] + C [k ]* (z X [k ] + D [k ]* (z X [k ] 3 ; 3 Reture g To design a special form of Gaussian elimination without partial pivoting, first observe that, since most of the entries in the coefficient matrix of (7 are zero, we only need a (m 1* 3 array to store all the non-zero entries of the coefficient matrix (see Figure and an one dimensional array to store the constants Besides, each equation, except the last one, will be used to eliminate the equation below it only For instance, equation 1 will be used to eliminate equation only Therefore, procedure SPECIAL_GAUSS can be designed as follows ( x + x 1 x 1 x 3 ( x 3 + x x x 4 ( x 4 + x 3 x 3 * * * * * * * * * x m 1 ( x m 1 + x m x m x m ( x m + x m 1 Figure 73 An (m 1 3 array for the coefficient matrix
11 Procedure SPECIAL_GAUSS (X, F : array[1m] of real; var B :array[1m] of real; var: M : array[1m-1,13] of real; H : array[1m-1] of real; a : real; 1 [Forward Elimination] for i :1 to m 1 do begin M [i,1]: X [i +] X [i ]; M [i,3]: X [i ] X [i ]; M [i,]: * (M [i,1] + M [i,3]; H [i ]: 3* (M [i,1]* (F [i ] F [i ]/M [i,3] + M [i,3]* (F [i +] F [i ]/M [i,1]; end; H [1]:H [1] M [1,1]*B [1]; H [m 1]:H [m 1] M [m 1,3]*B [m ]; for i :1 to m do begin a : M [i,1]/m [i,]; {Multiplier} M [i,]: M [i,] a*m [i,3]; H [i ]: H [i ] a*h [i ]; end; [Back Substitution] B [m ]: H [m 1]/M [m 1,]; for i :m downto 1 do B [i ]: (H [i ] M [i,3]*b [i +]/M [i,];
12 References 1 Ward Cheney and David Kincaid, Numerical Mathematics and Computing, 3rd Edition, Brooks/Cole, Pacific Grove, California, 1994 Samuel D Conte and Carl de Boor, Elementary Numerical Analysis: An Algorithmic Approach, 3rd Edition, McGraw-Hill, New York, Carl-Erik Froberg, Introduction to Numerical Analysis, nd edition, Svenska Bokforlaget Bonniers, Sweden, Shoichiro Nakamura, Applied Numerical Methods in C, Prentice Hall, Englewood Cliffs, New Jersey, Pat H Sterbenz, Floating-Point Computation, Prentice Hall, Englewood Cliffs, New Jersey
NUMERICAL MATHEMATICS & COMPUTING 7th Edition
NUMERICAL MATHEMATICS & COMPUTING 7th Edition Ward Cheney/David Kincaid c UT Austin Engage Learning: Thomson-Brooks/Cole wwwengagecom wwwmautexasedu/cna/nmc7 October 16, 2011 Ward Cheney/David Kincaid
More informationHani Mehrpouyan, California State University, Bakersfield. Signals and Systems
Hani Mehrpouyan, Department of Electrical and Computer Engineering, Lecture 26 (LU Factorization) May 30 th, 2013 The material in these lectures is partly taken from the books: Elementary Numerical Analysis,
More informationDM559 Linear and Integer Programming. Lecture 2 Systems of Linear Equations. Marco Chiarandini
DM559 Linear and Integer Programming Lecture Marco Chiarandini Department of Mathematics & Computer Science University of Southern Denmark Outline 1. Outline 1. 3 A Motivating Example You are organizing
More informationNUMERICAL MATHEMATICS AND COMPUTING
NUMERICAL MATHEMATICS AND COMPUTING Fourth Edition Ward Cheney David Kincaid The University of Texas at Austin 9 Brooks/Cole Publishing Company I(T)P An International Thomson Publishing Company Pacific
More informationNumerical Integration Schemes for Unequal Data Spacing
American Journal of Applied Mathematics 017; () 48- http//www.sciencepublishinggroup.com/j/ajam doi 10.1148/j.ajam.01700.1 ISSN 330-0043 (Print); ISSN 330-00X (Online) Numerical Integration Schemes for
More informationElectrical Engg. Department Semester :4
Semester :4 Teaching Scheme for Electrical Engineering Semester I Course Code Course Title Lecture Hours Tutorial Hours Practical Hours Credit MA 2002 Mathematics I: Introduction to Numerical Analysis
More information. =. a i1 x 1 + a i2 x 2 + a in x n = b i. a 11 a 12 a 1n a 21 a 22 a 1n. i1 a i2 a in
Vectors and Matrices Continued Remember that our goal is to write a system of algebraic equations as a matrix equation. Suppose we have the n linear algebraic equations a x + a 2 x 2 + a n x n = b a 2
More informationx x2 2 + x3 3 x4 3. Use the divided-difference method to find a polynomial of least degree that fits the values shown: (b)
Numerical Methods - PROBLEMS. The Taylor series, about the origin, for log( + x) is x x2 2 + x3 3 x4 4 + Find an upper bound on the magnitude of the truncation error on the interval x.5 when log( + x)
More informationEmpirical Models Interpolation Polynomial Models
Mathematical Modeling Lia Vas Empirical Models Interpolation Polynomial Models Lagrange Polynomial. Recall that two points (x 1, y 1 ) and (x 2, y 2 ) determine a unique line y = ax + b passing them (obtained
More informationInterpolation Theory
Numerical Analysis Massoud Malek Interpolation Theory The concept of interpolation is to select a function P (x) from a given class of functions in such a way that the graph of y P (x) passes through the
More informationRoots of Polynomials
Roots of Polynomials (Com S 477/577 Notes) Yan-Bin Jia Sep 26, 2017 A direct corollary of the fundamental theorem of algebra is that p(x) can be factorized over the complex domain into a product a n (x
More informationMATRICES. a m,1 a m,n A =
MATRICES Matrices are rectangular arrays of real or complex numbers With them, we define arithmetic operations that are generalizations of those for real and complex numbers The general form a matrix of
More informationTABLE OF CONTENTS INTRODUCTION, APPROXIMATION & ERRORS 1. Chapter Introduction to numerical methods 1 Multiple-choice test 7 Problem set 9
TABLE OF CONTENTS INTRODUCTION, APPROXIMATION & ERRORS 1 Chapter 01.01 Introduction to numerical methods 1 Multiple-choice test 7 Problem set 9 Chapter 01.02 Measuring errors 11 True error 11 Relative
More informationLectures 9-10: Polynomial and piecewise polynomial interpolation
Lectures 9-1: Polynomial and piecewise polynomial interpolation Let f be a function, which is only known at the nodes x 1, x,, x n, ie, all we know about the function f are its values y j = f(x j ), j
More informationSimultaneous Linear Equations
Simultaneous Linear Equations PHYSICAL PROBLEMS Truss Problem Pressure vessel problem a a b c b Polynomial Regression We are to fit the data to the polynomial regression model Simultaneous Linear Equations
More informationSection Gaussian Elimination
Section. - Gaussian Elimination A matrix is said to be in row echelon form (REF) if it has the following properties:. The first nonzero entry in any row is a. We call this a leading one or pivot one..
More informationarxiv: v2 [math.na] 9 Jun 2009
ON THE CHOLESKY METHOD arxiv:0906065v [mathna] 9 Jun 009 Christian Rakotonirina Institut Supérieur de Technologie d Antananarivo, IST-T, BP 8, Madagascar e-mail: rakotopierre@refermg June 9, 009 Abstract
More informationMA2501 Numerical Methods Spring 2015
Norwegian University of Science and Technology Department of Mathematics MA5 Numerical Methods Spring 5 Solutions to exercise set 9 Find approximate values of the following integrals using the adaptive
More informationInterpolation & Polynomial Approximation. Hermite Interpolation I
Interpolation & Polynomial Approximation Hermite Interpolation I Numerical Analysis (th Edition) R L Burden & J D Faires Beamer Presentation Slides prepared by John Carroll Dublin City University c 2011
More informationMAT 460: Numerical Analysis I. James V. Lambers
MAT 460: Numerical Analysis I James V. Lambers January 31, 2013 2 Contents 1 Mathematical Preliminaries and Error Analysis 7 1.1 Introduction............................ 7 1.1.1 Error Analysis......................
More informationCHAPTER 4. Interpolation
CHAPTER 4 Interpolation 4.1. Introduction We will cover sections 4.1 through 4.12 in the book. Read section 4.1 in the book on your own. The basic problem of one-dimensional interpolation is this: Given
More informationDirect Methods for Solving Linear Systems. Matrix Factorization
Direct Methods for Solving Linear Systems Matrix Factorization Numerical Analysis (9th Edition) R L Burden & J D Faires Beamer Presentation Slides prepared by John Carroll Dublin City University c 2011
More information5. Hand in the entire exam booklet and your computer score sheet.
WINTER 2016 MATH*2130 Final Exam Last name: (PRINT) First name: Student #: Instructor: M. R. Garvie 19 April, 2016 INSTRUCTIONS: 1. This is a closed book examination, but a calculator is allowed. The test
More informationJim Lambers MAT 460/560 Fall Semester Practice Final Exam
Jim Lambers MAT 460/560 Fall Semester 2009-10 Practice Final Exam 1. Let f(x) = sin 2x + cos 2x. (a) Write down the 2nd Taylor polynomial P 2 (x) of f(x) centered around x 0 = 0. (b) Write down the corresponding
More informationIntro Polynomial Piecewise Cubic Spline Software Summary. Interpolation. Sanzheng Qiao. Department of Computing and Software McMaster University
Interpolation Sanzheng Qiao Department of Computing and Software McMaster University July, 2012 Outline 1 Introduction 2 Polynomial Interpolation 3 Piecewise Polynomial Interpolation 4 Natural Cubic Spline
More informationMathematics for Engineers. Numerical mathematics
Mathematics for Engineers Numerical mathematics Integers Determine the largest representable integer with the intmax command. intmax ans = int32 2147483647 2147483647+1 ans = 2.1475e+09 Remark The set
More informationNUMERICAL MATHEMATICS & COMPUTING 7th Edition
NUMERICAL MATHEMATICS & COMPUTING 7th Edition Ward Cheney/David Kincaid c UT Austin Engage Learning: Thomson-Brooks/Cole wwwengagecom wwwmautexasedu/cna/nmc6 October 16, 2011 Ward Cheney/David Kincaid
More informationCAM Ph.D. Qualifying Exam in Numerical Analysis CONTENTS
CAM Ph.D. Qualifying Exam in Numerical Analysis CONTENTS Preliminaries Round-off errors and computer arithmetic, algorithms and convergence Solutions of Equations in One Variable Bisection method, fixed-point
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 informationMATH 3511 Lecture 1. Solving Linear Systems 1
MATH 3511 Lecture 1 Solving Linear Systems 1 Dmitriy Leykekhman Spring 2012 Goals Review of basic linear algebra Solution of simple linear systems Gaussian elimination D Leykekhman - MATH 3511 Introduction
More informationCurve Fitting. 1 Interpolation. 2 Composite Fitting. 1.1 Fitting f(x) 1.2 Hermite interpolation. 2.1 Parabolic and Cubic Splines
Curve Fitting Why do we want to curve fit? In general, we fit data points to produce a smooth representation of the system whose response generated the data points We do this for a variety of reasons 1
More information3.1 Interpolation and the Lagrange Polynomial
MATH 4073 Chapter 3 Interpolation and Polynomial Approximation Fall 2003 1 Consider a sample x x 0 x 1 x n y y 0 y 1 y n. Can we get a function out of discrete data above that gives a reasonable estimate
More informationLecture Note 3: Polynomial Interpolation. Xiaoqun Zhang Shanghai Jiao Tong University
Lecture Note 3: Polynomial Interpolation Xiaoqun Zhang Shanghai Jiao Tong University Last updated: October 24, 2013 1.1 Introduction We first look at some examples. Lookup table for f(x) = 2 π x 0 e x2
More information1 Review of Interpolation using Cubic Splines
cs412: introduction to numerical analysis 10/10/06 Lecture 12: Instructor: Professor Amos Ron Cubic Hermite Spline Interpolation Scribes: Yunpeng Li, Mark Cowlishaw 1 Review of Interpolation using Cubic
More informationEBG # 3 Using Gaussian Elimination (Echelon Form) Gaussian Elimination: 0s below the main diagonal
EBG # 3 Using Gaussian Elimination (Echelon Form) Gaussian Elimination: 0s below the main diagonal [ x y Augmented matrix: 1 1 17 4 2 48 (Replacement) Replace a row by the sum of itself and a multiple
More informationChapter 1: Systems of linear equations and matrices. Section 1.1: Introduction to systems of linear equations
Chapter 1: Systems of linear equations and matrices Section 1.1: Introduction to systems of linear equations Definition: A linear equation in n variables can be expressed in the form a 1 x 1 + a 2 x 2
More information7.6 The Inverse of a Square Matrix
7.6 The Inverse of a Square Matrix Copyright Cengage Learning. All rights reserved. What You Should Learn Verify that two matrices are inverses of each other. Use Gauss-Jordan elimination to find inverses
More informationGeneralization Of The Secant Method For Nonlinear Equations
Applied Mathematics E-Notes, 8(2008), 115-123 c ISSN 1607-2510 Available free at mirror sites of http://www.math.nthu.edu.tw/ amen/ Generalization Of The Secant Method For Nonlinear Equations Avram Sidi
More information1 - Systems of Linear Equations
1 - Systems of Linear Equations 1.1 Introduction to Systems of Linear Equations Almost every problem in linear algebra will involve solving a system of equations. ü LINEAR EQUATIONS IN n VARIABLES We are
More informationLECTURE NOTES ELEMENTARY NUMERICAL METHODS. Eusebius Doedel
LECTURE NOTES on ELEMENTARY NUMERICAL METHODS Eusebius Doedel TABLE OF CONTENTS Vector and Matrix Norms 1 Banach Lemma 20 The Numerical Solution of Linear Systems 25 Gauss Elimination 25 Operation Count
More informationReview. Numerical Methods Lecture 22. Prof. Jinbo Bi CSE, UConn
Review Taylor Series and Error Analysis Roots of Equations Linear Algebraic Equations Optimization Numerical Differentiation and Integration Ordinary Differential Equations Partial Differential Equations
More information2 Systems of Linear Equations
2 Systems of Linear Equations A system of equations of the form or is called a system of linear equations. x + 2y = 7 2x y = 4 5p 6q + r = 4 2p + 3q 5r = 7 6p q + 4r = 2 Definition. An equation involving
More informationMATHEMATICAL METHODS INTERPOLATION
MATHEMATICAL METHODS INTERPOLATION I YEAR BTech By Mr Y Prabhaker Reddy Asst Professor of Mathematics Guru Nanak Engineering College Ibrahimpatnam, Hyderabad SYLLABUS OF MATHEMATICAL METHODS (as per JNTU
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 informationLecture Note 3: Interpolation and Polynomial Approximation. Xiaoqun Zhang Shanghai Jiao Tong University
Lecture Note 3: Interpolation and Polynomial Approximation Xiaoqun Zhang Shanghai Jiao Tong University Last updated: October 10, 2015 2 Contents 1.1 Introduction................................ 3 1.1.1
More informationNote on Newton Interpolation Formula
Int. Journal of Math. Analysis, Vol. 6, 2012, no. 50, 2459-2465 Note on Newton Interpolation Formula Ramesh Kumar Muthumalai Department of Mathematics D.G. Vaishnav College, Arumbaam Chennai-600106, Tamilnadu,
More informationNUMERICAL METHODS C. Carl Gustav Jacob Jacobi 10.1 GAUSSIAN ELIMINATION WITH PARTIAL PIVOTING
. Gaussian Elimination with Partial Pivoting. Iterative Methods for Solving Linear Systems.3 Power Method for Approximating Eigenvalues.4 Applications of Numerical Methods Carl Gustav Jacob Jacobi 8 4
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 informationMath 365 Homework 5 Due April 6, 2018 Grady Wright
Math 365 Homework 5 Due April 6, 018 Grady Wright 1. (Polynomial interpolation, 10pts) Consider the following three samples of some function f(x): x f(x) -1/ -1 1 1 1/ (a) Determine the unique second degree
More informationMath 1314 Week #14 Notes
Math 3 Week # Notes Section 5.: A system of equations consists of two or more equations. A solution to a system of equations is a point that satisfies all the equations in the system. In this chapter,
More informationCS 257: Numerical Methods
CS 57: Numerical Methods Final Exam Study Guide Version 1.00 Created by Charles Feng http://www.fenguin.net CS 57: Numerical Methods Final Exam Study Guide 1 Contents 1 Introductory Matter 3 1.1 Calculus
More informationPARTIAL DIFFERENTIAL EQUATIONS
MATHEMATICAL METHODS PARTIAL DIFFERENTIAL EQUATIONS I YEAR B.Tech By Mr. Y. Prabhaker Reddy Asst. Professor of Mathematics Guru Nanak Engineering College Ibrahimpatnam, Hyderabad. SYLLABUS OF MATHEMATICAL
More informationPolynomial Interpolation
Polynomial Interpolation (Com S 477/577 Notes) Yan-Bin Jia Sep 1, 017 1 Interpolation Problem In practice, often we can measure a physical process or quantity (e.g., temperature) at a number of points
More informationChapter 9: Gaussian Elimination
Uchechukwu Ofoegbu Temple University Chapter 9: Gaussian Elimination Graphical Method The solution of a small set of simultaneous equations, can be obtained by graphing them and determining the location
More informationONE - DIMENSIONAL INTEGRATION
Chapter 4 ONE - DIMENSIONA INTEGRATION 4. Introduction Since the finite element method is based on integral relations it is logical to expect that one should strive to carry out the integrations as efficiently
More informationIntroduction to Numerical Analysis
Introduction to Numerical Analysis S. Baskar and S. Sivaji Ganesh Department of Mathematics Indian Institute of Technology Bombay Powai, Mumbai 400 076. Introduction to Numerical Analysis Lecture Notes
More informationReferences. Microcomputer, Macmillan. Angell, I. O. and Jones, B. J. (1983b). Advanced Graphics with the Sinclair ZX
References Adams, A. G. (1969).'Areas under the normal curve', Algorithm 39, Computer J., 12. American National Standards Institute (1978). American National Standard for Minimal BASIC, ANSI X3.60. Angell,
More informationSPLINE INTERPOLATION
Spline Background SPLINE INTERPOLATION Problem: high degree interpolating polynomials often have extra oscillations. Example: Runge function f(x = 1 1+4x 2, x [ 1, 1]. 1 1/(1+4x 2 and P 8 (x and P 16 (x
More informationLectures on Linear Algebra for IT
Lectures on Linear Algebra for IT by Mgr. Tereza Kovářová, Ph.D. following content of lectures by Ing. Petr Beremlijski, Ph.D. Department of Applied Mathematics, VSB - TU Ostrava Czech Republic 2. Systems
More informationChapter 7. Tridiagonal linear systems. Solving tridiagonal systems of equations. and subdiagonal. E.g. a 21 a 22 a A =
Chapter 7 Tridiagonal linear systems The solution of linear systems of equations is one of the most important areas of computational mathematics. A complete treatment is impossible here but we will discuss
More informationNumerical Mathematics & Computing, 7 Ed. 4.1 Interpolation
Numerical Mathematics & Computing, 7 Ed. 4.1 Interpolation Ward Cheney/David Kincaid c UT Austin Engage Learning: Thomson-Brooks/Cole www.engage.com www.ma.utexas.edu/cna/nmc6 November 7, 2011 2011 1 /
More informationMatrix notation. A nm : n m : size of the matrix. m : no of columns, n: no of rows. Row matrix n=1 [b 1, b 2, b 3,. b m ] Column matrix m=1
Matrix notation A nm : n m : size of the matrix m : no of columns, n: no of rows Row matrix n=1 [b 1, b 2, b 3,. b m ] Column matrix m=1 n = m square matrix Symmetric matrix Upper triangular matrix: matrix
More informationCS 450 Numerical Analysis. Chapter 8: Numerical Integration and Differentiation
Lecture slides based on the textbook Scientific Computing: An Introductory Survey by Michael T. Heath, copyright c 2018 by the Society for Industrial and Applied Mathematics. http://www.siam.org/books/cl80
More informationComputer Aided Design of Thermal Systems (ME648)
Computer Aided Design of Thermal Systems (ME648) PG/Open Elective Credits: 3-0-0-9 Updated Syallabus: Introduction. Basic Considerations in Design. Modelling of Thermal Systems. Numerical Modelling and
More informationIllustration of Gaussian elimination to find LU factorization. A = a 11 a 12 a 13 a 14 a 21 a 22 a 23 a 24 a 31 a 32 a 33 a 34 a 41 a 42 a 43 a 44
Illustration of Gaussian elimination to find LU factorization. A = a 21 a a a a 31 a 32 a a a 41 a 42 a 43 a 1 Compute multipliers : Eliminate entries in first column: m i1 = a i1 a 11, i = 2, 3, 4 ith
More informationNumerical Analysis: Solving Systems of Linear Equations
Numerical Analysis: Solving Systems of Linear Equations Mirko Navara http://cmpfelkcvutcz/ navara/ Center for Machine Perception, Department of Cybernetics, FEE, CTU Karlovo náměstí, building G, office
More informationGlobal polynomial interpolants suffer from the Runge Phenomenon if the data sites (nodes) are not chosen correctly.
Piecewise polynomial interpolation Global polynomial interpolants suffer from the Runge Phenomenon if the data sites (nodes) are not chosen correctly. In many applications, one does not have the freedom
More information9.1 - Systems of Linear Equations: Two Variables
9.1 - Systems of Linear Equations: Two Variables Recall that a system of equations consists of two or more equations each with two or more variables. A solution to a system in two variables is an ordered
More informationNUMERICAL METHODS C. Carl Gustav Jacob Jacobi 10.1 GAUSSIAN ELIMINATION WITH PARTIAL PIVOTING
0. Gaussian Elimination with Partial Pivoting 0.2 Iterative Methods for Solving Linear Systems 0.3 Power Method for Approximating Eigenvalues 0.4 Applications of Numerical Methods Carl Gustav Jacob Jacobi
More informationSystem of Linear Equations
Chapter 7 - S&B Gaussian and Gauss-Jordan Elimination We will study systems of linear equations by describing techniques for solving such systems. The preferred solution technique- Gaussian elimination-
More informationScientific Computing
2301678 Scientific Computing Chapter 2 Interpolation and Approximation Paisan Nakmahachalasint Paisan.N@chula.ac.th Chapter 2 Interpolation and Approximation p. 1/66 Contents 1. Polynomial interpolation
More informationSKMM 3023 Applied Numerical Methods
UNIVERSITI TEKNOLOGI MALAYSIA SKMM 3023 Applied Numerical Methods Numerical Differentiation ibn Abdullah Faculty of Mechanical Engineering Òº ÙÐÐ ÚºÒÙÐÐ ¾¼½ SKMM 3023 Applied Numerical Methods Numerical
More informationY = ax + b. Numerical Applications Least-squares. Start with Self-test 10-1/459. Linear equation. Error function: E = D 2 = (Y - (ax+b)) 2
Ch.10 Numerical Applications 10-1 Least-squares Start with Self-test 10-1/459. Linear equation Y = ax + b Error function: E = D 2 = (Y - (ax+b)) 2 Regression Formula: Slope a = (N ΣXY - (ΣX)(ΣY)) / (N
More informationMA 323 Geometric Modelling Course Notes: Day 07 Parabolic Arcs
MA 323 Geometric Modelling Course Notes: Day 07 Parabolic Arcs David L. Finn December 9th, 2004 We now start considering the basic curve elements to be used throughout this course; polynomial curves and
More informationBACHELOR OF COMPUTER APPLICATIONS (BCA) (Revised) Term-End Examination December, 2015 BCS-054 : COMPUTER ORIENTED NUMERICAL TECHNIQUES
No. of Printed Pages : 5 BCS-054 BACHELOR OF COMPUTER APPLICATIONS (BCA) (Revised) Term-End Examination December, 2015 058b9 BCS-054 : COMPUTER ORIENTED NUMERICAL TECHNIQUES Time : 3 hours Maximum Marks
More informationCurve Fitting and Interpolation
Chapter 5 Curve Fitting and Interpolation 5.1 Basic Concepts Consider a set of (x, y) data pairs (points) collected during an experiment, Curve fitting: is a procedure to develop or evaluate mathematical
More informationON NUMERICAL METHODS FOR ELLIPTIC TRANSMISSION/BOUNDARY VALUE PROBLEMS
ON NUMERICAL METHODS FOR ELLIPTIC TRANSMISSION/BOUNDARY VALUE PROBLEMS ANIRBAN ROY 349 LAFAYETTE DR, BETHLEHEM, PA, 182 Project Supervisors: Anna L. Mazzucato and Victor Nistor Department of Mathematics,
More informationLU Factorization. Marco Chiarandini. DM559 Linear and Integer Programming. Department of Mathematics & Computer Science University of Southern Denmark
DM559 Linear and Integer Programming LU Factorization Marco Chiarandini Department of Mathematics & Computer Science University of Southern Denmark [Based on slides by Lieven Vandenberghe, UCLA] Outline
More informationThe Exponential of a Matrix
The Exponential of a Matrix 5-8- The solution to the exponential growth equation dx dt kx is given by x c e kt It is natural to ask whether you can solve a constant coefficient linear system x A x in a
More informationPivoting. Reading: GV96 Section 3.4, Stew98 Chapter 3: 1.3
Pivoting Reading: GV96 Section 3.4, Stew98 Chapter 3: 1.3 In the previous discussions we have assumed that the LU factorization of A existed and the various versions could compute it in a stable manner.
More informationA fourth order method for finding a simple root of univariate function
Bol. Soc. Paran. Mat. (3s.) v. 34 2 (2016): 197 211. c SPM ISSN-2175-1188 on line ISSN-00378712 in press SPM: www.spm.uem.br/bspm doi:10.5269/bspm.v34i1.24763 A fourth order method for finding a simple
More information2.1 Gaussian Elimination
2. Gaussian Elimination A common problem encountered in numerical models is the one in which there are n equations and n unknowns. The following is a description of the Gaussian elimination method for
More information2.29 Numerical Fluid Mechanics Spring 2015 Lecture 4
2.29 Spring 2015 Lecture 4 Review Lecture 3 Truncation Errors, Taylor Series and Error Analysis Taylor series: 2 3 n n i1 i i i i i n f( ) f( ) f '( ) f ''( ) f '''( )... f ( ) R 2! 3! n! n1 ( n1) Rn f
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 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 informationMatrix decompositions
Matrix decompositions How can we solve Ax = b? 1 Linear algebra Typical linear system of equations : x 1 x +x = x 1 +x +9x = 0 x 1 +x x = The variables x 1, x, and x only appear as linear terms (no powers
More informationNumerical Analysis Fall. Gauss Elimination
Numerical Analysis 2015 Fall Gauss Elimination Solving systems m g g m m g x x x k k k k k k k k k 3 2 1 3 2 1 3 3 3 2 3 2 2 2 1 0 0 Graphical Method For small sets of simultaneous equations, graphing
More informationAIMS Exercise Set # 1
AIMS Exercise Set #. Determine the form of the single precision floating point arithmetic used in the computers at AIMS. What is the largest number that can be accurately represented? What is the smallest
More informationFinite Mathematics Chapter 2. where a, b, c, d, h, and k are real numbers and neither a and b nor c and d are both zero.
Finite Mathematics Chapter 2 Section 2.1 Systems of Linear Equations: An Introduction Systems of Equations Recall that a system of two linear equations in two variables may be written in the general form
More informationFitting a Natural Spline to Samples of the Form (t, f(t))
Fitting a Natural Spline to Samples of the Form (t, f(t)) David Eberly, Geometric Tools, Redmond WA 9852 https://wwwgeometrictoolscom/ This work is licensed under the Creative Commons Attribution 4 International
More informationHermite normal form: Computation and applications
Integer Points in Polyhedra Gennady Shmonin Hermite normal form: Computation and applications February 24, 2009 1 Uniqueness of Hermite normal form In the last lecture, we showed that if B is a rational
More information4.1 Classical Formulas for Equally Spaced Abscissas. 124 Chapter 4. Integration of Functions
24 Chapter 4. Integration of Functions of various orders, with higher order sometimes, but not always, giving higher accuracy. Romberg integration, which is discussed in 4.3, is a general formalism for
More informationInverses. Stephen Boyd. EE103 Stanford University. October 28, 2017
Inverses Stephen Boyd EE103 Stanford University October 28, 2017 Outline Left and right inverses Inverse Solving linear equations Examples Pseudo-inverse Left and right inverses 2 Left inverses a number
More informationIntroduction to Decision Sciences Lecture 6
Introduction to Decision Sciences Lecture 6 Andrew Nobel September 21, 2017 Functions Functions Given: Sets A and B, possibly different Definition: A function f : A B is a rule that assigns every element
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 informationMAC Module 2 Systems of Linear Equations and Matrices II. Learning Objectives. Upon completing this module, you should be able to :
MAC 0 Module Systems of Linear Equations and Matrices II Learning Objectives Upon completing this module, you should be able to :. Find the inverse of a square matrix.. Determine whether a matrix is invertible..
More informationCS412: Introduction to Numerical Methods
CS412: Introduction to Numerical Methods MIDTERM #1 2:30PM - 3:45PM, Tuesday, 03/10/2015 Instructions: This exam is a closed book and closed notes exam, i.e., you are not allowed to consult any textbook,
More informationThere is a unique function s(x) that has the required properties. It turns out to also satisfy
Numerical Analysis Grinshpan Natural Cubic Spline Let,, n be given nodes (strictly increasing) and let y,, y n be given values (arbitrary) Our goal is to produce a function s() with the following properties:
More informationINTERPOLATION. and y i = cos x i, i = 0, 1, 2 This gives us the three points. Now find a quadratic polynomial. p(x) = a 0 + a 1 x + a 2 x 2.
INTERPOLATION Interpolation is a process of finding a formula (often a polynomial) whose graph will pass through a given set of points (x, y). As an example, consider defining and x 0 = 0, x 1 = π/4, x
More informationStudent Study Session. Theorems
Students should be able to apply and have a geometric understanding of the following: Intermediate Value Theorem Mean Value Theorem for derivatives Extreme Value Theorem Name Formal Statement Restatement
More information