Interpolation Theory
|
|
- Ralph Lucas
- 6 years ago
- Views:
Transcription
1 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 given data points (x i, y i ), i 1,,, n The points may arise as measurements in a physical problem, or they may be obtained from a known function The interpolating function is usually chosen from a restricted class of functions, with polynomials being the most commonly used class Theorem 1 (Weierstrass Approximation Theorem) If f(x) is a continuous function on [a, b] and ɛ > is given, then there exists a polynomial P (x), defined on [a, b], with the property that f(x) P (x) < ɛ for all x [a, b] Theorem Given n + 1 distinct points x, x 1,, x n and n + 1 ordinates y, y 1,, y n, there is a polynomial P (x) of degree n that interpolates to y i at x i, i, 1,,, n This polynomial P (x) is unique among the set of all polynomials of degree at most n Proof Consider the following linear system of n + 1 equations and n + 1 unknowns a + a 1 x + a x + a n x n y a + a 1 x n + a x n + a n x n n y n The unknowns are (a, a 1,, a n ) The proof follows from the fact that this system has a unique solution Taylor Polynomials Let f(x) be n + 1 times differentiable near a point x Then the n-degree Taylor polynomial of f(x) at x is The function P n (x) f(x ) + f (x )(x x ) + f (x ) (x x ) + f (n) (x ) (x x ) n! n! R n (x) P n (x) f(x) f (n+1) (ξ(x)) (x x ) (n+1) (n + 1)! is called the Remainder function or the Error function Example Approximate ln 11 using a Taylor polynomial of degree 3 Note that if f(x) 1 1+x, then 1 f(x)dx ln 11 Thus we need to integrate the polynomial P 3 (x) 1 x + x x 3 from zero to 1 The error is 1 (1 x + x x 3 )dx 1 1 [ f(x) P 3 (x) x 4 dx x 5 5 ] 1 [x x + x3 3 x ] 1 Bernstein Polynomials Given a function defined on [, 1], the Bernstein polynomial of degree n for f(x) is given by n ( ) n B n (x) f( k k n )xk (1 x) n k k
2 Massoud Malek Interpolation Theory Page where ( ) n k denotes n! k!(n k)! It can be shown that, if f(x) is continuous on [, 1] and x [, 1], then lim B n(x ) f(x ) n The polynomial can be used in a constructive proof of Weierstrass Theorem Lagrange Polynomial Given two points (x, y ) and (x 1, y 1 ) with x x 1, we can draw a straight line through them The line is the graph of the linear polynomial P 1 (x) (x 1 x) x 1 x y + (x x ) x 1 x y 1 Given n + 1 distinct points x, x 1,, x n and n + 1 ordinates y, y 1, ldots, y n, according to Theorem, there is a unique polynomial P n (x) that interpolates to y k at x k, k, 1,,, n For k 1,,, n, define the polynomial L n,k (x) (x x n ) (x x k 1 )(x x k+1 ) (x x n ) (x k x ) (x k x k 1 )(x k x k+1 ) (x k x n ) Then the polynomial L n (x) n y k L n,k (x) k is called the Lagrange Interpolation Polynomial Note that L n,k (x j ) δ k,j, where δ k,j is called Kronecker delta function, { 1, j δ k,j k;, j k k, j n j k n (x x j ) (x k x j ) Theorem 3 If x, x 1, x n are distinct numbers in the interval [a, b] and if f(x) C n+1 [a, b], then, for each x [a, b], there exist a number ξ(x) (a, b) such that f(x) L n (x) + f (n+1) (ξ(x)) (x x )(x x 1 ) (x x n ) (n + 1)! Example Let f(x) 1 1+x, the find the Lagrange polynomial of degree two passing through the points (, f()), (1, f(1)) and (, f()) Hence We have f() 1, f(1) 1, f() 1 3 and L (x) k (x 1)(x ) L, (x) ( 1)( ) L,1 (x) L, (x) (x 1)(x ) x 3x + (x )(x ) (1 )(1 ) (x )(x ) x + x (x )(x 1) ( )( 1) L,k (x)f(k) x 3x + (x )(x 1) x x (1) + ( x + x)( 1 ) + x x ( 1 3 ) x 6 x 3 + 1
3 Massoud Malek Interpolation Theory Page 3 Inverse Interpolation Suppose f(x) C 1 [a, b], f (x) on [a, b] and f(x) has one zero x in [a, b] Then clearly f(x) is a one to one function on [a, b] By choosing few numbers x k s in [a, b], we create the points (y k, x k ) s and using these points, we define the Lagrange polynomial n L n (y) L n,k (y)x k This way L() becomes an approximation to x k Newton Divided Differences The Lagrange form of the interpolation polynomial can be used for interpolating a function But there are other forms that are much more convenient Suppose P n 1 (x) is a polynomial of degree n that interpolate f(x) We would like to pass from P n 1 (x) to another interpolation polynomial of degree one greater; so we need to find P n (x) P n 1 (x) + C(x) C(x) correction term Then in general C(x) is a polynomial of degree n, since usually degree(p n 1 ) n 1 and degree(p n ) n Also we have C(x k ) P n (x k ) P n 1 (x k ) f(x k ) f(x k ) k, 1,, n 1 Since P n (x) f(x n ), we have f(x n ) P n 1 (x n ) a n (x n x ) (x n x n 1 ) This coefficient a n is called the n-th order Newton Divided Difference of f(x), and it is denoted by a n f[x, x 1,, x n ] Thus our interpolation formula becomes P n (x) P n 1 (x) + (x x ) (x x n 1 )f[x,, x n ] Using the Lagrange formula, we obtain f[x, x 1,, x n ] f[x 1, x,, x n ] f[x, x 1,, x n 1 ] x n x We have P (x) f(x ) P 1 (x) f(x ) + (x x )f[x, x 1 ] P n (x) f(x ) + (x x )f[x, x 1 ] + (x x )(x x 1 )f[x, x 1, x ] + (x x ) (x x n 1 )f[x, x 1,, x n ] This is called Newton s divided difference formula for the interpolation polynomial It is much better for computation that the Lagrange formula To construct the divided differences, use the format shown in the following table x i f(x i ) f[x i, x i+1 ] D f[x i ] D 3 f[x i ] D 4 f[x i ] x y f[x, x 1 ] x 1 y 1 f[x, x 1, x ] f[x 1, x ] f[x, x 1, x, x 3 ] x y f[x 1, x, x 3 ] f[x, x 4 ] f[x, x 3 ] f[x 1, x, x 3, x 4 ] x 3 y 3 f[x, x 3, x 4 ] f[x 1,, x 5 ]
4 Massoud Malek Interpolation Theory Page 4 Example We construct a divided difference table for f(x) x, shown in the next table x i f(x i ) f[x i, x i+1 ] D f[x i ] D 3 f[x i ] D 4 f[x i ] The polynomial P 4 (x) (x )[(3494) + (x 1)[( 411) + (x )[(9167) + (x 3)[( 84)]]]] is known as the Newton Forward Divided-Difference form The following polynomial is known as the Newton Backward Divided-Difference form: Q 4 (x) (x 4)[(3618)+(x 3)[( 3585)+(x )[(8333)+(x 1)[( 84)]]]] Algorithm DIVDIF(N,X,D) Remark: On entrance, X (X(), X(N)) and D (F (X 1 ) F (X n )) On exit, D(J) will contain F [X, X 1, X J ] STEP1 INPUT N a positive integer; STEP INPUT X a vector of points x(), X(N); STEP3 INPUT D a vector of values of some function evaluated at the nodes in X; STEP4 FOR I : 1, N; FOR J : N, I, 1; D(J) : D(J) D(J 1) X(J) X(J I) RETURN; STEP5 STOP Algorithm FORWDINTERP(N,X,F) Remark: On entrance, X (X(), X(N)) and D (D(), D(1), D(n) On exit, P will contain the value P N (T ) of the N-th degree polynomial interpolation of F (X) STEP1 P : D N ; STEP FOR K : N 1,, 1; P : D(K) + (T X(K) )P RETURN; STEP3 STOP Note By using Divided-difference and Inverse Interpolation, we may approximate the zero of certain functions Hermite Interpolation For a variety of applications, it is convenient to consider polynomials P (x) that interpolate to a function f(x) at certain points x, x 1, x n and in addition to have the derivative polynomial P (x) interpolate to the derivative function f (x) Since the first derivative of f(x) and P (x) agree at the interpolation points, they have basically the same shape To construct the Hermite interpolation polynomial we use the functions L n,k (x) defined in the Lagrange formula Consider the function f(x) and n + 1 points x, x 1, x n We use the Newton s Divided Differences by choosing the following points: (x, y ), (x, y ), (x 1, y 1 ), (x 1, y 1 ), (x, y ), (x, y ),, (x n, y n ), (x n, y n )
5 Massoud Malek Interpolation Theory Page 5 and replacing f[x i, x i ] with f (x i ) The following table illustrates the interpolation method: x i f(x i ) f[x i, x i+1 ] D f[x i ] D 3 f[x i ] D 4 f[x i ] x y f (x ) x y f[x, x, x 1 ] f[x, x 1 ] f[x, x, x 1, x 1 ] x 1 y 1 f[x, x 1, x 1 ] f[x, x, x 1, x 1, x ] f (x 1 ) f[x, x 1, x 1, x ] x 1 y 1 f[x 1, x 1, x ] f[x, x 1, x 1, x, x ] f[x 1, x ] f[x, x 1, x 1, x ] x y f[x 1, x, x ] f[x 1, x 1, x, x, x 3 ] f (x ) f[x 1, x 1, x, x ] x y f[x, x, x 3 ] f[x 1, x, x, x 3, x 3 ] Interpolation Using Spline Functions The simplest method of interpolation is to connect the points points by straight line segments; this is called Piecewise Linear Interpolation, and the associated interpolation function is denoted by S(x) Although S(x) agrees with the data, but it has the disadvantage of not having a smooth graph Another choice is to connect the data points by using succession of quadratic interpolation polynomials Although the resulting curve is smoother than the first method, but the graph can still have some corners For many applications, this interpolating function may be completely adequate there, but we would like to obtain a smooth curve We shall see that a cubic spline produces the desired curve Cubic Spline Suppose n + 1 data points (x, y ), (x 1, y 1 ),, (x n, y n ) are given For simplicity, assume x < x 1 < x n and let a x, b x n We seek a function S(x) defined on [a, b] that interpolates the data S(x i ) y i, i, 1,, n For smoothness of S(x), we need to require that S (x) and S (x) be continuous In addition, we would like to have the curve follow the general shape given by the piecewise linear function connecting the data points (x i, y i ) The standard way in which this has been done is to ask that the derivative S (x) not change too rapidly between note points This has been carried out by requiring the second derivative S (x) to be as small as possible and, more precisely, by requiring that b a [S (x)] dx be also as small as possiblethis may not be a perfect mathematical realization of the idea of a smooth shape-preserving interpolation function for the data {(x i, y i )}, but it usually gives a very good interpolating function There is a unique solution S(x) to this problem, and it satisfies the following: S1 S(x) is a cubic polynomial, denoted S j (x) on each [x j, x j+1 ], for j, 1, n 1; S S(x j ) f(x j ) for each j, 1,, n; S3 S j+1 (x j+1 ) S j (x j+1 ) for each j, 1,, n ;
6 Massoud Malek Interpolation Theory Page 6 S4 S j+1 (x j+1) S j (x j+1) for each j, 1,, n ; S5 S j+1 (x j+1) S j (x j+1) for each j, 1,, n ; S6 one of the following set of boundary conditions is satisfied: (i) S (x ) S (x n ) (Free boundary) (ii) S (x ) f (x ) and S (x n ) f (x n ) (Clamped boundary) To construct the cubic spline interpolant for a given function f(x), the conditions (S 1 through S 6 ) can be applied to the cubic polynomials S j (x) A j + B j (x x j ) + C j (x x j ) + D j (x x j ) 3 for each j, 1,, n 1 For j, 1,, n 1, define h j x j+1 x j, then we have A j S j (x j ) f(x j ) and A j+1 A j + B j h j + C j h j + D j h 3 j B f (x ), B j S j(x j ) (j > ) and B j+1 B j + C j h j + 3D j h j B j S (x j ), C n S (x n ) C j+1 C j + 3D j h j Since h j s and A j s are given, once the values of C j s are known, then it is a simple matter to find B j s and D j s Free Boundary When S (a) S (b), then the vector C (C, C 1,, C n ) t is the unique solution to the equation 1 h (h + h 1 ) h 1 h 1 (h 1 + h ) h h n (h n + h n 1 ) h n 1 1 C C 1 C n 1 C n 3(A A 1) h 1 3(A n A n 1) h n 1 3(A1 A) h 3(An 1 An ) h n Clamped Boundary For S (a) f (a) and S (b) f (b), the vector C (C, C 1,, C n ) t is the unique solution to the equation h h h (h + h 1 ) h 1 h 1 (h 1 + h ) h h n (h n + h n 1 ) h n 1 h n 1 h n 1 C C 1 C n 1 C n 3(A 1 A ) h 3(A A 1) h 1 3(A n A n 1) h n 1 3f (a) 3(A1 A) h 3(An 1 An ) h n 3f (b) 3(An An 1) h n 1
3.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 informationChapter 3 Interpolation and Polynomial Approximation
Chapter 3 Interpolation and Polynomial Approximation Per-Olof Persson persson@berkeley.edu Department of Mathematics University of California, Berkeley Math 128A Numerical Analysis Polynomial Interpolation
More informationCubic Splines MATH 375. J. Robert Buchanan. Fall Department of Mathematics. J. Robert Buchanan Cubic Splines
Cubic Splines MATH 375 J. Robert Buchanan Department of Mathematics Fall 2006 Introduction Given data {(x 0, f(x 0 )), (x 1, f(x 1 )),...,(x n, f(x n ))} which we wish to interpolate using a polynomial...
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 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 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 informationExam 2. Average: 85.6 Median: 87.0 Maximum: Minimum: 55.0 Standard Deviation: Numerical Methods Fall 2011 Lecture 20
Exam 2 Average: 85.6 Median: 87.0 Maximum: 100.0 Minimum: 55.0 Standard Deviation: 10.42 Fall 2011 1 Today s class Multiple Variable Linear Regression Polynomial Interpolation Lagrange Interpolation Newton
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 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 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 informationMath 651 Introduction to Numerical Analysis I Fall SOLUTIONS: Homework Set 1
ath 651 Introduction to Numerical Analysis I Fall 2010 SOLUTIONS: Homework Set 1 1. Consider the polynomial f(x) = x 2 x 2. (a) Find P 1 (x), P 2 (x) and P 3 (x) for f(x) about x 0 = 0. What is the relation
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 informationSome notes on Chapter 8: Polynomial and Piecewise-polynomial Interpolation
Some notes on Chapter 8: Polynomial and Piecewise-polynomial Interpolation See your notes. 1. Lagrange Interpolation (8.2) 1 2. Newton Interpolation (8.3) different form of the same polynomial as Lagrange
More informationInput: A set (x i -yy i ) data. Output: Function value at arbitrary point x. What for x = 1.2?
Applied Numerical Analysis Interpolation Lecturer: Emad Fatemizadeh Interpolation Input: A set (x i -yy i ) data. Output: Function value at arbitrary point x. 0 1 4 1-3 3 9 What for x = 1.? Interpolation
More informationNumerical Methods. King Saud University
Numerical Methods King Saud University Aims In this lecture, we will... find the approximate solutions of derivative (first- and second-order) and antiderivative (definite integral only). Numerical Differentiation
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 information(0, 0), (1, ), (2, ), (3, ), (4, ), (5, ), (6, ).
1 Interpolation: The method of constructing new data points within the range of a finite set of known data points That is if (x i, y i ), i = 1, N are known, with y i the dependent variable and x i [x
More informationData Analysis-I. Interpolation. Soon-Hyung Yook. December 4, Soon-Hyung Yook Data Analysis-I December 4, / 1
Data Analysis-I Interpolation Soon-Hyung Yook December 4, 2015 Soon-Hyung Yook Data Analysis-I December 4, 2015 1 / 1 Table of Contents Soon-Hyung Yook Data Analysis-I December 4, 2015 2 / 1 Introduction
More informationChapter 1 Numerical approximation of data : interpolation, least squares method
Chapter 1 Numerical approximation of data : interpolation, least squares method I. Motivation 1 Approximation of functions Evaluation of a function Which functions (f : R R) can be effectively evaluated
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 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 informationLecture 10 Polynomial interpolation
Lecture 10 Polynomial interpolation Weinan E 1,2 and Tiejun Li 2 1 Department of Mathematics, Princeton University, weinan@princeton.edu 2 School of Mathematical Sciences, Peking University, tieli@pku.edu.cn
More informationNovember 20, Interpolation, Extrapolation & Polynomial Approximation
Interpolation, Extrapolation & Polynomial Approximation November 20, 2016 Introduction In many cases we know the values of a function f (x) at a set of points x 1, x 2,..., x N, but we don t have the analytic
More informationQ 0 x if x 0 x x 1. S 1 x if x 1 x x 2. i 0,1,...,n 1, and L x L n 1 x if x n 1 x x n
. - Piecewise Linear-Quadratic Interpolation Piecewise-polynomial Approximation: Problem: Givenn pairs of data points x i, y i, i,,...,n, find a piecewise-polynomial Sx S x if x x x Sx S x if x x x 2 :
More informationHermite Interpolation
Jim Lambers MAT 77 Fall Semester 010-11 Lecture Notes These notes correspond to Sections 4 and 5 in the text Hermite Interpolation Suppose that the interpolation points are perturbed so that two neighboring
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 informationNumerical Analysis: Interpolation Part 1
Numerical Analysis: Interpolation Part 1 Computer Science, Ben-Gurion University (slides based mostly on Prof. Ben-Shahar s notes) 2018/2019, Fall Semester BGU CS Interpolation (ver. 1.00) AY 2018/2019,
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 information6x 2 8x + 5 ) = 12x 8
Example. If f(x) = x 3 4x + 5x + 1, then f (x) = 6x 8x + 5 Observation: f (x) is also a differentiable function... d dx ( f (x) ) = d dx ( 6x 8x + 5 ) = 1x 8 The derivative of f (x) is called the second
More informationCurve Fitting. Objectives
Curve Fitting Objectives Understanding the difference between regression and interpolation. Knowing how to fit curve of discrete with least-squares regression. Knowing how to compute and understand the
More informationApplied Numerical Analysis Quiz #2
Applied Numerical Analysis Quiz #2 Modules 3 and 4 Name: Student number: DO NOT OPEN UNTIL ASKED Instructions: Make sure you have a machine-readable answer form. Write your name and student number in the
More informationIntroduction Linear system Nonlinear equation Interpolation
Interpolation Interpolation is the process of estimating an intermediate value from a set of discrete or tabulated values. Suppose we have the following tabulated values: y y 0 y 1 y 2?? y 3 y 4 y 5 x
More informationNewton's forward interpolation formula
Newton's Interpolation Formulae Interpolation is the process of approximating a given function, whose values are known at N + 1 tabular points, by a suitable polynomial, P N (x) of degree N which takes
More informationUNIT-II INTERPOLATION & APPROXIMATION
UNIT-II INTERPOLATION & APPROXIMATION LAGRANGE POLYNAMIAL 1. Find the polynomial by using Lagrange s formula and hence find for : 0 1 2 5 : 2 3 12 147 : 0 1 2 5 : 0 3 12 147 Lagrange s interpolation formula,
More informationIntroductory Numerical Analysis
Introductory Numerical Analysis Lecture Notes December 16, 017 Contents 1 Introduction to 1 11 Floating Point Numbers 1 1 Computational Errors 13 Algorithm 3 14 Calculus Review 3 Root Finding 5 1 Bisection
More informationMath Numerical Analysis
Math 541 - Numerical Analysis Joseph M. Mahaffy, jmahaffy@mail.sdsu.edu Department of Mathematics and Statistics Dynamical Systems Group Computational Sciences Research Center San Diego State University
More informationInterpolation and the Lagrange Polynomial
Interpolation and the Lagrange Polynomial MATH 375 J. Robert Buchanan Department of Mathematics Fall 2013 Introduction We often choose polynomials to approximate other classes of functions. Theorem (Weierstrass
More informationNUMERICAL METHODS. x n+1 = 2x n x 2 n. In particular: which of them gives faster convergence, and why? [Work to four decimal places.
NUMERICAL METHODS 1. Rearranging the equation x 3 =.5 gives the iterative formula x n+1 = g(x n ), where g(x) = (2x 2 ) 1. (a) Starting with x = 1, compute the x n up to n = 6, and describe what is happening.
More informationInterpolation APPLIED PROBLEMS. Reading Between the Lines FLY ROCKET FLY, FLY ROCKET FLY WHAT IS INTERPOLATION? Figure Interpolation of discrete data.
WHAT IS INTERPOLATION? Given (x 0,y 0 ), (x,y ), (x n,y n ), find the value of y at a value of x that is not given. Interpolation Reading Between the Lines Figure Interpolation of discrete data. FLY ROCKET
More information1 Piecewise Cubic Interpolation
Piecewise Cubic Interpolation Typically the problem with piecewise linear interpolation is the interpolant is not differentiable as the interpolation points (it has a kinks at every interpolation point)
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 informationInterpolation and Polynomial Approximation I
Interpolation and Polynomial Approximation I If f (n) (x), n are available, Taylor polynomial is an approximation: f (x) = f (x 0 )+f (x 0 )(x x 0 )+ 1 2! f (x 0 )(x x 0 ) 2 + Example: e x = 1 + x 1! +
More informationMath 473: Practice Problems for Test 1, Fall 2011, SOLUTIONS
Math 473: Practice Problems for Test 1, Fall 011, SOLUTIONS Show your work: 1. (a) Compute the Taylor polynomials P n (x) for f(x) = sin x and x 0 = 0. Solution: Compute f(x) = sin x, f (x) = cos x, f
More informationDepartment of Applied Mathematics and Theoretical Physics. AMA 204 Numerical analysis. Exam Winter 2004
Department of Applied Mathematics and Theoretical Physics AMA 204 Numerical analysis Exam Winter 2004 The best six answers will be credited All questions carry equal marks Answer all parts of each question
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 informationAn idea how to solve some of the problems. diverges the same must hold for the original series. T 1 p T 1 p + 1 p 1 = 1. dt = lim
An idea how to solve some of the problems 5.2-2. (a) Does not converge: By multiplying across we get Hence 2k 2k 2 /2 k 2k2 k 2 /2 k 2 /2 2k 2k 2 /2 k. As the series diverges the same must hold for the
More informationComputational Physics
Interpolation, Extrapolation & Polynomial Approximation Lectures based on course notes by Pablo Laguna and Kostas Kokkotas revamped by Deirdre Shoemaker Spring 2014 Introduction In many cases, a function
More informationChapter 5: Numerical Integration and Differentiation
Chapter 5: Numerical Integration and Differentiation PART I: Numerical Integration Newton-Cotes Integration Formulas The idea of Newton-Cotes formulas is to replace a complicated function or tabulated
More informationInterpolation. P. Sam Johnson. January 30, P. Sam Johnson (NITK) Interpolation January 30, / 75
Interpolation P. Sam Johnson January 30, 2015 P. Sam Johnson (NITK) Interpolation January 30, 2015 1 / 75 Overview One of the basic ideas in Mathematics is that of a function and most useful tool of numerical
More informationSection Taylor and Maclaurin Series
Section.0 Taylor and Maclaurin Series Ruipeng Shen Feb 5 Taylor and Maclaurin Series Main Goal: How to find a power series representation for a smooth function us assume that a smooth function has a power
More informationMATH 1372, SECTION 33, MIDTERM 3 REVIEW ANSWERS
MATH 1372, SECTION 33, MIDTERM 3 REVIEW ANSWERS 1. We have one theorem whose conclusion says an alternating series converges. We have another theorem whose conclusion says an alternating series diverges.
More informationCOURSE Numerical integration of functions
COURSE 6 3. Numerical integration of functions The need: for evaluating definite integrals of functions that has no explicit antiderivatives or whose antiderivatives are not easy to obtain. Let f : [a,
More information6x 2 8x + 5 ) = 12x 8. f (x) ) = d (12x 8) = 12
AMS/ECON 11A Class Notes 11/6/17 UCSC *) Higher order derivatives Example. If f = x 3 x + 5x + 1, then f = 6x 8x + 5 Observation: f is also a differentiable function... d f ) = d 6x 8x + 5 ) = 1x 8 dx
More informationMADHAVA MATHEMATICS COMPETITION, December 2015 Solutions and Scheme of Marking
MADHAVA MATHEMATICS COMPETITION, December 05 Solutions and Scheme of Marking NB: Part I carries 0 marks, Part II carries 30 marks and Part III carries 50 marks Part I NB Each question in Part I carries
More informationMATH ASSIGNMENT 07 SOLUTIONS. 8.1 Following is census data showing the population of the US between 1900 and 2000:
MATH4414.01 ASSIGNMENT 07 SOLUTIONS 8.1 Following is census data showing the population of the US between 1900 and 2000: Years after 1900 Population in millions 0 76.0 20 105.7 40 131.7 60 179.3 80 226.5
More informationAdditional exercises with Numerieke Analyse
Additional exercises with Numerieke Analyse March 10, 017 1. (a) Given different points x 0, x 1, x [a, b] and scalars y 0, y 1, y, z 1, show that there exists at most one polynomial p P 3 with p(x i )
More informationLECTURE 16 GAUSS QUADRATURE In general for Newton-Cotes (equispaced interpolation points/ data points/ integration points/ nodes).
CE 025 - Lecture 6 LECTURE 6 GAUSS QUADRATURE In general for ewton-cotes (equispaced interpolation points/ data points/ integration points/ nodes). x E x S fx dx hw' o f o + w' f + + w' f + E 84 f 0 f
More informationNumerical Methods of Approximation
Contents 31 Numerical Methods of Approximation 31.1 Polynomial Approximations 2 31.2 Numerical Integration 28 31.3 Numerical Differentiation 58 31.4 Nonlinear Equations 67 Learning outcomes In this Workbook
More informationExample 1 Which of these functions are polynomials in x? In the case(s) where f is a polynomial,
1. Polynomials A polynomial in x is a function of the form p(x) = a 0 + a 1 x + a 2 x 2 +... a n x n (a n 0, n a non-negative integer) where a 0, a 1, a 2,..., a n are constants. We say that this polynomial
More informationEngg. Math. II (Unit-IV) Numerical Analysis
Dr. Satish Shukla of 33 Engg. Math. II (Unit-IV) Numerical Analysis Syllabus. Interpolation and Curve Fitting: Introduction to Interpolation; Calculus of Finite Differences; Finite Difference and Divided
More informationSection 3.1. Best Affine Approximations. Difference Equations to Differential Equations
Difference Equations to Differential Equations Section 3.1 Best Affine Approximations We are now in a position to discuss the two central problems of calculus as mentioned in Section 1.1. In this chapter
More informationDRAFT - Math 101 Lecture Note - Dr. Said Algarni
2 Limits 2.1 The Tangent Problems The word tangent is derived from the Latin word tangens, which means touching. A tangent line to a curve is a line that touches the curve and a secant line is a line that
More informationAP Calculus (BC) Chapter 9 Test No Calculator Section Name: Date: Period:
WORKSHEET: Series, Taylor Series AP Calculus (BC) Chapter 9 Test No Calculator Section Name: Date: Period: 1 Part I. Multiple-Choice Questions (5 points each; please circle the correct answer.) 1. The
More informationLecture 32: Taylor Series and McLaurin series We saw last day that some functions are equal to a power series on part of their domain.
Lecture 32: Taylor Series and McLaurin series We saw last day that some functions are equal to a power series on part of their domain. For example f(x) = 1 1 x = 1 + x + x2 + x 3 + = ln(1 + x) = x x2 2
More information(f(x) P 3 (x)) dx. (a) The Lagrange formula for the error is given by
1. QUESTION (a) Given a nth degree Taylor polynomial P n (x) of a function f(x), expanded about x = x 0, write down the Lagrange formula for the truncation error, carefully defining all its elements. How
More informationApplied Numerical Analysis Homework #3
Applied Numerical Analysis Homework #3 Interpolation: Splines, Multiple dimensions, Radial Bases, Least-Squares Splines Question Consider a cubic spline interpolation of a set of data points, and derivatives
More informationNewton-Raphson Type Methods
Int. J. Open Problems Compt. Math., Vol. 5, No. 2, June 2012 ISSN 1998-6262; Copyright c ICSRS Publication, 2012 www.i-csrs.org Newton-Raphson Type Methods Mircea I. Cîrnu Department of Mathematics, Faculty
More information1 Lecture 8: Interpolating polynomials.
1 Lecture 8: Interpolating polynomials. 1.1 Horner s method Before turning to the main idea of this part of the course, we consider how to evaluate a polynomial. Recall that a polynomial is an expression
More informationFunction Approximation
1 Function Approximation This is page i Printer: Opaque this 1.1 Introduction In this chapter we discuss approximating functional forms. Both in econometric and in numerical problems, the need for an approximating
More informationCubic Splines. Antony Jameson. Department of Aeronautics and Astronautics, Stanford University, Stanford, California, 94305
Cubic Splines Antony Jameson Department of Aeronautics and Astronautics, Stanford University, Stanford, California, 94305 1 References on splines 1. J. H. Ahlberg, E. N. Nilson, J. H. Walsh. Theory of
More informationCALCULUS JIA-MING (FRANK) LIOU
CALCULUS JIA-MING (FRANK) LIOU Abstract. Contents. Power Series.. Polynomials and Formal Power Series.2. Radius of Convergence 2.3. Derivative and Antiderivative of Power Series 4.4. Power Series Expansion
More informationGENG2140, S2, 2012 Week 7: Curve fitting
GENG2140, S2, 2012 Week 7: Curve fitting Curve fitting is the process of constructing a curve, or mathematical function, f(x) that has the best fit to a series of data points Involves fitting lines and
More informationPreliminary Examination in Numerical Analysis
Department of Applied Mathematics Preliminary Examination in Numerical Analysis August 7, 06, 0 am pm. Submit solutions to four (and no more) of the following six problems. Show all your work, and justify
More information100 CHAPTER 4. SYSTEMS AND ADAPTIVE STEP SIZE METHODS APPENDIX
100 CHAPTER 4. SYSTEMS AND ADAPTIVE STEP SIZE METHODS APPENDIX.1 Norms If we have an approximate solution at a given point and we want to calculate the absolute error, then we simply take the magnitude
More informationNumerical integration and differentiation. Unit IV. Numerical Integration and Differentiation. Plan of attack. Numerical integration.
Unit IV Numerical Integration and Differentiation Numerical integration and differentiation quadrature classical formulas for equally spaced nodes improper integrals Gaussian quadrature and orthogonal
More informationMath Numerical Analysis Mid-Term Test Solutions
Math 400 - Numerical Analysis Mid-Term Test Solutions. Short Answers (a) A sufficient and necessary condition for the bisection method to find a root of f(x) on the interval [a,b] is f(a)f(b) < 0 or f(a)
More informationChapter 3: Root Finding. September 26, 2005
Chapter 3: Root Finding September 26, 2005 Outline 1 Root Finding 2 3.1 The Bisection Method 3 3.2 Newton s Method: Derivation and Examples 4 3.3 How To Stop Newton s Method 5 3.4 Application: Division
More informationTaylor and Maclaurin Series. Approximating functions using Polynomials.
Taylor and Maclaurin Series Approximating functions using Polynomials. Approximating f x = e x near x = 0 In order to approximate the function f x = e x near x = 0, we can use the tangent line (The Linear
More informationHomework and Computer Problems for Math*2130 (W17).
Homework and Computer Problems for Math*2130 (W17). MARCUS R. GARVIE 1 December 21, 2016 1 Department of Mathematics & Statistics, University of Guelph NOTES: These questions are a bare minimum. You should
More informationScientific Computing: An Introductory Survey
Scientific Computing: An Introductory Survey Chapter 7 Interpolation Prof. Michael T. Heath Department of Computer Science University of Illinois at Urbana-Champaign Copyright c 2002. Reproduction permitted
More informationLecture 34: Recall Defn: The n-th Taylor polynomial for a function f at a is: n f j (a) j! + f n (a)
Lecture 34: Recall Defn: The n-th Taylor polynomial for a function f at a is: n f j (a) P n (x) = (x a) j. j! j=0 = f(a)+(f (a))(x a)+(1/2)(f (a))(x a) 2 +(1/3!)(f (a))(x a) 3 +... + f n (a) (x a) n n!
More informationLECTURE 10: REVIEW OF POWER SERIES. 1. Motivation
LECTURE 10: REVIEW OF POWER SERIES By definition, a power series centered at x 0 is a series of the form where a 0, a 1,... and x 0 are constants. For convenience, we shall mostly be concerned with the
More informationPolynomial Approximations and Power Series
Polynomial Approximations and Power Series June 24, 206 Tangent Lines One of the first uses of the derivatives is the determination of the tangent as a linear approximation of a differentiable function
More information8.5 Taylor Polynomials and Taylor Series
8.5. TAYLOR POLYNOMIALS AND TAYLOR SERIES 50 8.5 Taylor Polynomials and Taylor Series Motivating Questions In this section, we strive to understand the ideas generated by the following important questions:
More informationTaylor and Maclaurin Series. Approximating functions using Polynomials.
Taylor and Maclaurin Series Approximating functions using Polynomials. Approximating f x = e x near x = 0 In order to approximate the function f x = e x near x = 0, we can use the tangent line (The Linear
More informationROOT FINDING REVIEW MICHELLE FENG
ROOT FINDING REVIEW MICHELLE FENG 1.1. Bisection Method. 1. Root Finding Methods (1) Very naive approach based on the Intermediate Value Theorem (2) You need to be looking in an interval with only one
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 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 informationSection 5.8. Taylor Series
Difference Equations to Differential Equations Section 5.8 Taylor Series In this section we will put together much of the work of Sections 5.-5.7 in the context of a discussion of Taylor series. We begin
More informationHow might we evaluate this? Suppose that, by some good luck, we knew that. x 2 5. x 2 dx 5
8.4 1 8.4 Partial Fractions Consider the following integral. 13 2x (1) x 2 x 2 dx How might we evaluate this? Suppose that, by some good luck, we knew that 13 2x (2) x 2 x 2 = 3 x 2 5 x + 1 We could then
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 information1. Let A R be a nonempty set that is bounded from above, and let a be the least upper bound of A. Show that there exists a sequence {a n } n N
Applied Analysis prelim July 15, 216, with solutions Solve 4 of the problems 1-5 and 2 of the problems 6-8. We will only grade the first 4 problems attempted from1-5 and the first 2 attempted from problems
More informationLecture 4: Numerical solution of ordinary differential equations
Lecture 4: Numerical solution of ordinary differential equations Department of Mathematics, ETH Zürich General explicit one-step method: Consistency; Stability; Convergence. High-order methods: Taylor
More informationInterpolation. Escuela de Ingeniería Informática de Oviedo. (Dpto. de Matemáticas-UniOvi) Numerical Computation Interpolation 1 / 34
Interpolation Escuela de Ingeniería Informática de Oviedo (Dpto. de Matemáticas-UniOvi) Numerical Computation Interpolation 1 / 34 Outline 1 Introduction 2 Lagrange interpolation 3 Piecewise polynomial
More informationChapter 11. Taylor Series. Josef Leydold Mathematical Methods WS 2018/19 11 Taylor Series 1 / 27
Chapter 11 Taylor Series Josef Leydold Mathematical Methods WS 2018/19 11 Taylor Series 1 / 27 First-Order Approximation We want to approximate function f by some simple function. Best possible approximation
More informationReview I: Interpolation
Review I: Interpolation Varun Shankar January, 206 Introduction In this document, we review interpolation by polynomials. Unlike many reviews, we will not stop there: we will discuss how to differentiate
More informationMath 1310 Section 4.1: Polynomial Functions and Their Graphs. A polynomial function is a function of the form ...
Math 1310 Section 4.1: Polynomial Functions and Their Graphs A polynomial function is a function of the form... where 0,,,, are real numbers and n is a whole number. The degree of the polynomial function
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 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 informationNumerical Analysis Preliminary Exam 10.00am 1.00pm, January 19, 2018
Numerical Analysis Preliminary Exam 0.00am.00pm, January 9, 208 Instructions. You have three hours to complete this exam. Submit solutions to four (and no more) of the following six problems. Please start
More information