(0, 0), (1, ), (2, ), (3, ), (4, ), (5, ), (6, ).
|
|
- Lenard Knight
- 5 years ago
- Views:
Transcription
1 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 1, x N ], then in order to find the value of y at x = x 0, and x 0 [x 1, x N ] then we can use interpolation If x 0 / [x 1, x N ] then we have extrapolation 2 Consider the points (x i, y i ) given by (0, 0), (1, 08415), (2, 09093), (3, 01411), (4, 07568), (5, 09589), (6, 02794) Interpolation means finding a value for y given any x (0, 6) 3 A function y = f(x) is said to interpolate the data points (x 1, y 1 ),,(x n, y n ) if f(x i ) = y i for each 1 i n Since f is a function, this means that the x i s must all be distinct - there is no such restriction on the y i s The main theorem of polynomial interpolation shows that: Ifx 0, x 1, x n are distinct real numbers, then for arbitrary values, y 1,, y n, there is a unique polynomial f n (x), of degree n such that f n (x i ) = y i, (0 i n) a) Example: Find an interpolation ploynomial for the set of points (5, 1), ( 7, 23) b) Find an interpolation polynomial for (x 0, y 0 ), (x 0, y 1 ) We repeat the above arguments to find: Letting f 1 (x) = y 1 x x 0 x x 0 + y 0 x x 1 x 0 x 1 l 0 (x) = x x 1 x 0 x 1, l 1 (x) = x x 0 x 1 x 0, we have f 1 (x) = y 0 l 0 (x) + y 1 l 1 (x) 4 Piecewise constant interpolation - assign the value to be the same as the nearest known datapoint eg at x = 025 the piecewise constant interpolated value would be Obviously not that accurate but can be used in multidimensional problems 5 Linear Interpolation - straight line between two data points If you have (x a, y a ), (x b, y b ), then the interpolant at x [x a, x b ] is y = y a + (y b y a ) x x a x b x a Not that precise though its fast and easy The interpolant is not differentiable at the point x b, x a The error is f(x) g(x) C(x b x a ) 2, 1
2 where C = 1 8 max y [x a,x b ] g (y) Problem - show this The error in linear interpolation is proportional to the square of the distance between the data points 6 Fit N +1 points with a N th degree polynomial f(x) = exact function of which only N + 1 discrete values are known and used to establish an interpolating or approximating function p(x) p(x) is the approximating or interpolating function and will pass through all the N + 1 specified interpolation points or nodes 7 Polynomial Interpolation - replace linear function with a polynomial function If we have N data points, then there is one polynomial of degree N 1 that goes through all the data points So for the 7 data points here, there is just one polynomial of degree 6 that goes through all the data points a) (x i, y i ), i (0, 1,,,, n) b) p(x) = a n x n + a n 1 x n 1 + a n 2 x n a 1 x + a 0, with p(x i ) = y i, i (0, 1,, n) c) Substituting this into the above we get the matrix equation Ax = d, where the matrix A is the Vandermonde matrix Problem Find the form of this matrix equation We can solve this matrix equation for the coefficients a n (contained in the vector x, using Gaussian elimination, but if the matrix A is large, this can lead to large errors in computing the coefficients a i d) For Lagrange Interpolation, note that p(x) = N f i V i (x), i=0 where V i (x) is a polynomil of degree N associated with each node i such that V i (x j ) = 0, i j, V i (x j ) = 1, i = j, and f i = p(x i ) With 5 interpolation nodes, we have p(x 3 ) = f 0 V 0 (x 3 ) + f 1 V 1 (x 3 ) + f 2 V 2 (x 3 ) + f 3 V 3 (x 3 ) + f 4 V 4 (x 3 ), 2
3 which means V 0 (x 3 ) = 0, V 1 (x 3 ) = 0, V 2 (x 3 ) = 0, V 3 (x 3 ) = 1, V 4 (x 3 ) = 0 and so p(x 3 ) = f 3 In order to construct the V i (x) we need a degree N polynomial with roots at x 0, x 1, x i 1, x i+1,, x N and V i (x i ) = 1 The function W i (x) = (x x 0 )(x x 1 )(x x 2 )(x x i 1 )(x x i+1 )(x x N ) has the required roots - it equals zero at the nodes x 0, x 1, x 2,, x N except at node x i The degree of W i (x) is N After the following normalization, we get the Lagrange basis functions V i (x), V i (x) = (x x 0)(x x 1 )(x x 2 )(x x i 1 )(x x i+1 )(x x N ) (x i x 0 )(x i x 1 )(x i x 2 )(x i x i 1 )(x i x i+1 )(x i x N ) and V i (x i ) = 1 Also V i (x j ) = 0, i j e) For Linear Lagrange (N=1), we have where p(x) = f 0 V 0 (x) + f 1 V 1 (x), V 0 (x) = (x x 0) (x 0 x 1 ), V 1(x) = (x x 0) (x 1 x 0 ) For example, given the following data: x 0 = 2, x 1 = 5, f 0 = 15, f 1 = 40, find the linear interpolating function p(x) Lagrange basis functions are V 0 (x) = (5 x), V 3 1 (x) = (x 2) The interpolating polynomial is p(x) = 3 15V V 1 f) For the following data, x 0 = 3, x 1 = 4, x 2 = 5, f 0 = 1, f 1 = 2, f 2 = 4 we find the quadratic interpolating function: Lagrange basis functions are V 0 (x) = (x 4)(x 5) (3 4)(3 5), V 1 (x) = V 2 (x) = (x 3)(x 5) (4 3)(4 5), (x 3)(x 4) (5 3)(5 4), and the interpolating function is p(x) = 10V 0 (x) + 20V 1 (x) + 40V 2 (x) e) For cubic Lagrange interpolation (N = 3), use the following data to construct the interpolating polynomial (040, ), (05, ), (07, ), (080, ) 3
4 Whats the interpolated value at x = 06? g) Using Taylor series type arguments, the error can be shown to be e(x) = f(x) p(x), where p(x) is the interpolating polynomial and f(x) is the exact function you want to interpolate e(x) = L(x)f (N+1) (ξ), x 0 ξ x N, where f N+1 (ξ) is the N + 1 th derivative of f wrt x evaluated at ξ, and L(x) = (x x 0)(x x 1 )(x x N ) (N + 1)! L(x) is minimum at the center of [x 0, x N ] and maximum near the edges, L(x) = 0 at all data points As the size of the interpolating domain [x 0, x N ] increases, so does the maximum error in the interval As N increases, L max decreases But for N > N crit, L max increases Properties of f (N+1) (ξ) will also influence the situation h) Estimate the error made in the previous example knowing that f(x) = ln(x) e(x) L(x)f N+1 (x m ), where N = 3, x m = 06, then e(06) = and E(x) = ln(060) g(060) = d) In terms of Lagrange polynomials, we can write the interpolating polynomial as p(x) = n i=0 Π 0 j n,j i x x j x i x j This interpolating polynomial is unique e) Let f be a continuous function in [a, b] and let P be a ploynomial of degree smaller than or equal to n that interpolates the function f at n+1 distinct points x 0, x 1,, x n [a, n] Then to each x [a, b], there exists a ξ x [a, b] such that f(x) p(x) = 1 (n + 1)! f(n+1) (ξ x )Π n i=0(x x i ) Put another way, if the function f satisfies the bound f (n+1) (x) M for x [a, b], then f(x) p n (x) 4 1 4(n + 1) M(b a n )n+1
5 f) Let f(x) = sinx and approximate it with an interpolation polynomial with 10 equally spaced nodes in [0, 16875]n = 9, a = 0, b = 16875, f (10) (x) = sinx, f (10) (x) 1, forallx [a, b] Then sinx p n (x) 1 40 (1)(16875 ) If all the derivatives of f can be uniformly bounded for all n, then we can force the error to zero by increasing the number of interpolation nodes But if we cannot uniformly bound all the derivatives, then increasing the number of nodes does not necessarily drive the error to zero even if f and all its derivatives are continuous g) Computationally expensive, can exhibit oscillatory artifacts (Runge s phenomena) Can restrict attention to Chebyshev polynomials h) Newton Polynomials: a) Let q 0 (x) = 1, q 1 (x) = x x 0, q 2 = (x x 0 )(x x 1 ),, q n (x) = (x x 0 )(x x 1 )(x x n 1 ), then the q 0 (x), q 1 (x),, q n (x) is a basis span of 1, x, x 2,, x n This leads to the Newton form of the interpolating polynomial: f n (x) = c 0 q 0 (x) + c 1 q 1 (x) + + c n q n (x) Given the data (x 0, f(x 0 )), (x 1, f(x 1 )),,(x n, f(x n )), we have to specify the following coefficients: c 1, i = 0, 1,, n b) Divided difference: The zeroth order divided difference is f[x i ] = f(x i ) Then given x i, x j and f(x i ), f(x j ), the first order divided difference is f[x i, x j ] = f(x i) f(x j ) x i x j Note that f[x i, x j ] is an approximation of the derivatives f (x i ) and f (x j ) if they exist Higher order divided differences are defined recursively: f[x 0, x 1, x 2, x n ] = f[x 0, x 1,,,,, x n 1 ] f[x 1, x 2,, x n ] x 0 x n c) Whats f[x 0, x 1 ], f[x 0, x 1, x 2 ], f[x 0, x 1, x 2, x 3 ]? If we have (3, 1), (1, 3), (5, 2), (6, 4) what are the numerical values of these divided differences? d) For any x, f[x, x 0 ] = f(x) f(x 0) x x 0, and f(x) = f(x 0 ) + f[x, x 0 ](x x 0 ), 5
6 and and So then f[x, x 0, x 1 ] = f[x, x 0] f[x 0, x 1 ] x x 1, f[x, x 0 ] = f[x 0, x 1 ] + f[x, x 0, x 1 ](x x 1 ) f(x) = f(x 0 ) + f[x 0, x 1 ](x x 0 ) + f[x, x 0, x 1 ](x x 0 )(x x 1 ) Carrying on the with the third order divided difference, we have f(x) = f(x 0 )+f[x 0, x 1 ](x x 0 )+f[x 0, x 1, x 2 ](x x 0 )(x x 1 )+f[x, x 0, x 1, x 2, x 3 ](x x 0 )(x Hence f(x) = f(x 0 ) + f[x 0, x 1 ](x x 0 ) + + f[x 0,, x n ](x x 0 )(x x n 1 ) Take the Newton form of the interpolation polynomial as N n (x) = f(x 0 ) + f[x 0, x 1 ](x x 0 ) + + f[x 0,, x n ](x x 0 )(x x n 1 ) e) Check that N n (x i ) = f(x i ), 1 i n f) So c 0 = f(x 0 ), c 1 = f[x 0, x 1 ],, c n = f[x 0, x 1,, x n ] g) Whats the difference between Lagrange and Newton form? h) Example: f(x) = sinx(0, f(0)) = (0, 0), (π/2, f(π/2)) given Then L 1 (x) = 0 x π/2) 0 π/2 + 1 x 0 π/2 0, N 1 (x) = f(x 0 ) + f[x 0, x 1 ](x x 0 ) = (x 0), 0 π/2 This can be used to compute the approximation at x = π/8 Now if we add the point (π/4, f(π/4), we can see that for the Newton form, we only have to add one term but for the Lagrange form we have to start again a) ) Forward Newton-Gregory: Forward differences are defined as follows: 0 f i = f i = f(x i ), f i = f i+1 f i, 2 f i = f i+1 f i, 6
7 2 f i = (f i+2 f i+1 ) (f i+1 f i ) = f i+2 2f i+1 f i, 3 f i = 2 f i+1 2 f i, 3 f i = f i+3 3f i+2 + 3f i+1 f i k f i = k 1 f i+1 k 1 f i N + 1 data points develop upto N th order forward differences i) Develop a forward difference table for (2, 7), (4, 3), (6, 6), (8, 25), (10, 62), (12, 129) b) The Taylor series expansion for f(x) about x 0 is f(x) = f(x 0 )+(x x 0 ) df dx x=x ! (x x 0) 2 d2 f dx 2 x = x ! (x x 0) d3 f dx 3 x=x 0 +O(x x 4 ) 4 f(x) = f 0 + (x x 0 )f (1) ! (x x 0) 2 f (2) ! (x x 0) 4 f (3) 0 + O(x x 0 ) 4, Since f 0 = f 1 f 0 and because f 1 = f(x 1 ), we can write, f 1 = f 0 +(x 1 x 0 )f (1) ! (x 1 x 0 ) 2 f (2) ! (x 1 x 0 ) 3 f (3) 0 +O(x 1 x 0 ) 4, If h = x 1 x 0, then so that f 1 = f 0 + hf (1) ! h2 f (2) ! h3 f (3) 0 + O(h 4 ), f (1) 0 = f 0 h 1 2! f(2) 0 1 3! h2 f (3) 0 O(h 3 ) Now express the second order forward difference in terms of f 0, f (1) 0,, 2 f 0 = f 2 2f 1 + f 0, Noting that f 1 = f(x 1 ) and we already have an expression for it and f 2 = f(x 2 ), f 2 = f 0 +(x 2 x 0 )f (1) ! (x 2 x 0 ) 2 f (2) ! (x 2 x 0 ) 3 f (3) 0 +O(x 2 x 0 4 ) Note that x 2 x 0 = 2h, f 2 = f 0 + 2hf (1) ! h2 f (2) ! f(3) 0 + O(h 4 ), 7
8 and hence f (2) 0 = 2 f 0 h 2 hf (3) 0 + O(h) 2, Now we express the third order difference in terms of f 0, f (1) 0, so that f (3) 0 = 3 f 0 h 3 Noting that f(x) = g(x) + e(x), we have + O(h) g(x) = f 0 + (x x 0 ) f 0 h + 1 2! (x x 0)(x x 1 ) 2 f 0 h ! 3 f 0 h N! (x x 0)(x x 1 )(x x 2 )(x x N 1 ) N f 0 h N, the N th degree polynomial approximation to N + 1 data points and is identical to that derived for Lagrange interpolation l) The error function is e(x) = f(x) g(x) = (x x 0)(x x 1 )(x x N ) f (N+1) (ξ), x 0 ξ x N (N + 1)! This is identical to the error for Lagrange Interpolation But note that f (N+1) (x 0 ) (N+1 )f 0 h N+1 8 Newton Backward Interpolation: This essentially the same as Newton forward interpolation except that backward differences are used a) These are defined as: 0 f i = f i, f i = f i f i 1, 2 f i = f i f i 1, k f i = k 1 f i k 1 f i 1, b) For N + 1 data points fitted with an Nth degree polynomial, we have g(x) = f N +(x x N ) f N h + 1 2! (x x N)(x x N 1 ) 2 f N h ! (x x N)(x x N 1 (x x N 2 ) 8
9 9 These Newtonian formulae are easily developed for non-unifrmly spaced nodes 10 Fourier Interpolation: Interpolation with trigonometric polynomials a) A trigonometric polynomial is of the form p(x) = a 0 + K a k cos(kx) + k=1 K b k sin(kx) There are 2K + 1 coefficients, a 0, a 1, a 2, a K, b 1, b K, and we wish to compute these coefficients so that the function passes through N points, (x i, y i ), 0 i N 1 Since the trigonometric polynomial is periodic with period 2π, the N points can be distributed and ordered in one period as 0 < x 0 < x 1 < x 2 < < x N 1 < 2π, the problem is to find the coefficients such that the trigonometric polynomial p satisfies these conditions We assume N 2K + 1 a) Similar to the formula for Lagrange interpolation, the solution is: p(x) = 2K Π 2K m=0,m k k=0 k=1 sin(05(x x m ) sin(05(x k x m )) b) When the points x n are equally spaced, we have x n = 2π n N, 0 < n < N And Y k = N 1 n=0 y n = p(x n ) = 1 N nk i2π y n e N, N 1 k=0 nk i2π Y k e N This is the discrete Fourier transform 6 Can sometimes reduce errors by Chebyshev Interpolation - where Chebyshev polynomials are used instead of ordinary polynomials 7 Rational Function interpolation - this is interpolation using a rational function - a function which is just a ratio of two polynomials: f(x) = p 0 + p 1 x + + p L x L q 0 + q 1 x + + q N x N Polynomial interpolation methods can sometimes be unstable on equidistant grids If we can choose the grids, we can sometimes get around this with either Chebysehv polynomials or rational function interpolation However though we can always find an interpolating polynomial, not all sets of points have an interpolating rational function and poles can sometimes also be a big problem 9
10 5 Extrapolation: find p(x M ) when x M is outside of [x 0, x N ] 6 Linear extrapolation: y(x M ) y k 1 + (x k x k 1 ) x k x k 1 (y k y k 1 ), where (x k, y k ), (x k 1, y k 1 ) are the two data points nearest x M 7 A polynomial curve using Lagrange or Newton methods can be fitted through the entire data or just near the end and the resulting polynomial used to extrapolate the data 10
November 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 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 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 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 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 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 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 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 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 informationApplied Numerical Analysis (AE2220-I) R. Klees and R.P. Dwight
Applied Numerical Analysis (AE0-I) R. Klees and R.P. Dwight February 018 Contents 1 Preliminaries: Motivation, Computer arithmetic, Taylor series 1 1.1 Numerical Analysis Motivation..........................
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(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 informationWe consider the problem of finding a polynomial that interpolates a given set of values:
Chapter 5 Interpolation 5. Polynomial Interpolation We consider the problem of finding a polynomial that interpolates a given set of values: x x 0 x... x n y y 0 y... y n where the x i are all distinct.
More informationMidterm Review. Igor Yanovsky (Math 151A TA)
Midterm Review Igor Yanovsky (Math 5A TA) Root-Finding Methods Rootfinding methods are designed to find a zero of a function f, that is, to find a value of x such that f(x) =0 Bisection Method To apply
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 informationNotation Nodes are data points at which functional values are available or at which you wish to compute functional values At the nodes fx i
LECTURE 6 NUMERICAL DIFFERENTIATION To find discrete approximations to differentiation (since computers can only deal with functional values at discrete points) Uses of numerical differentiation To represent
More informationIntegration, differentiation, and root finding. Phys 420/580 Lecture 7
Integration, differentiation, and root finding Phys 420/580 Lecture 7 Numerical integration Compute an approximation to the definite integral I = b Find area under the curve in the interval Trapezoid Rule:
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 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 informationInterpolation. Chapter Interpolation. 7.2 Existence, Uniqueness and conditioning
76 Chapter 7 Interpolation 7.1 Interpolation Definition 7.1.1. Interpolation of a given function f defined on an interval [a,b] by a polynomial p: Given a set of specified points {(t i,y i } n with {t
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 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 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 informationMath 578: Assignment 2
Math 578: Assignment 2 13. Determine whether the natural cubic spline that interpolates the table is or is not the x 0 1 2 3 y 1 1 0 10 function 1 + x x 3 x [0, 1] f(x) = 1 2(x 1) 3(x 1) 2 + 4(x 1) 3 x
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 informationFixed point iteration and root finding
Fixed point iteration and root finding The sign function is defined as x > 0 sign(x) = 0 x = 0 x < 0. It can be evaluated via an iteration which is useful for some problems. One such iteration is given
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 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 informationMA2501 Numerical Methods Spring 2015
Norwegian University of Science and Technology Department of Mathematics MA2501 Numerical Methods Spring 2015 Solutions to exercise set 7 1 Cf. Cheney and Kincaid, Exercise 4.1.9 Consider the data points
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 informationTaylor Series and Numerical Approximations
Taylor Series and Numerical Approximations Hilary Weller h.weller@reading.ac.uk August 7, 05 An introduction to the concept of a Taylor series and how these are used in numerical analysis to find numerical
More informationNumerical Analysis Solution of Algebraic Equation (non-linear equation) 1- Trial and Error. 2- Fixed point
Numerical Analysis Solution of Algebraic Equation (non-linear equation) 1- Trial and Error In this method we assume initial value of x, and substitute in the equation. Then modify x and continue till we
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 informationInterpolation and Approximation
Interpolation and Approximation The Basic Problem: Approximate a continuous function f(x), by a polynomial p(x), over [a, b]. f(x) may only be known in tabular form. f(x) may be expensive to compute. Definition:
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 informationCS 323: Numerical Analysis and Computing
CS 323: Numerical Analysis and Computing MIDTERM #2 Instructions: This is an open notes exam, i.e., you are allowed to consult any textbook, your class notes, homeworks, or any of the handouts from us.
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 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 informationA first order divided difference
A first order divided difference For a given function f (x) and two distinct points x 0 and x 1, define f [x 0, x 1 ] = f (x 1) f (x 0 ) x 1 x 0 This is called the first order divided difference of f (x).
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 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 informationNumerical Methods in Physics and Astrophysics
Kostas Kokkotas 2 October 17, 2017 2 http://www.tat.physik.uni-tuebingen.de/ kokkotas Kostas Kokkotas 3 TOPICS 1. Solving nonlinear equations 2. Solving linear systems of equations 3. Interpolation, approximation
More informationFunction approximation
Week 9: Monday, Mar 26 Function approximation A common task in scientific computing is to approximate a function. The approximated function might be available only through tabulated data, or it may be
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 informationNumerical interpolation, extrapolation and fi tting of data
Chapter 6 Numerical interpolation, extrapolation and fi tting of data 6.1 Introduction Numerical interpolation and extrapolation is perhaps one of the most used tools in numerical applications to physics.
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 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 informationMarch Algebra 2 Question 1. March Algebra 2 Question 1
March Algebra 2 Question 1 If the statement is always true for the domain, assign that part a 3. If it is sometimes true, assign it a 2. If it is never true, assign it a 1. Your answer for this question
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 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 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 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 informationOutline. 1 Interpolation. 2 Polynomial Interpolation. 3 Piecewise Polynomial Interpolation
Outline Interpolation 1 Interpolation 2 3 Michael T. Heath Scientific Computing 2 / 56 Interpolation Motivation Choosing Interpolant Existence and Uniqueness Basic interpolation problem: for given data
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 informationBSM510 Numerical Analysis
BSM510 Numerical Analysis Polynomial Interpolation Prof. Manar Mohaisen Department of EEC Engineering Review of Precedent Lecture Polynomial Regression Multiple Linear Regression Nonlinear Regression Lecture
More informationNumerical Methods in Physics and Astrophysics
Kostas Kokkotas 2 October 20, 2014 2 http://www.tat.physik.uni-tuebingen.de/ kokkotas Kostas Kokkotas 3 TOPICS 1. Solving nonlinear equations 2. Solving linear systems of equations 3. Interpolation, approximation
More informationyou expect to encounter difficulties when trying to solve A x = b? 4. A composite quadrature rule has error associated with it in the following form
Qualifying exam for numerical analysis (Spring 2017) Show your work for full credit. If you are unable to solve some part, attempt the subsequent parts. 1. Consider the following finite difference: f (0)
More informationNumerical Analysis Preliminary Exam 10 am to 1 pm, August 20, 2018
Numerical Analysis Preliminary Exam 1 am to 1 pm, August 2, 218 Instructions. You have three hours to complete this exam. Submit solutions to four (and no more) of the following six problems. Please start
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 information7x 5 x 2 x + 2. = 7x 5. (x + 1)(x 2). 4 x
Advanced Integration Techniques: Partial Fractions The method of partial fractions can occasionally make it possible to find the integral of a quotient of rational functions. Partial fractions gives us
More informationTwo hours. To be provided by Examinations Office: Mathematical Formula Tables. THE UNIVERSITY OF MANCHESTER. 29 May :45 11:45
Two hours MATH20602 To be provided by Examinations Office: Mathematical Formula Tables. THE UNIVERSITY OF MANCHESTER NUMERICAL ANALYSIS 1 29 May 2015 9:45 11:45 Answer THREE of the FOUR questions. If more
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 informationInterpolation. 1. Judd, K. Numerical Methods in Economics, Cambridge: MIT Press. Chapter
Key References: Interpolation 1. Judd, K. Numerical Methods in Economics, Cambridge: MIT Press. Chapter 6. 2. Press, W. et. al. Numerical Recipes in C, Cambridge: Cambridge University Press. Chapter 3
More information1 Solutions to selected problems
Solutions to selected problems Section., #a,c,d. a. p x = n for i = n : 0 p x = xp x + i end b. z = x, y = x for i = : n y = y + x i z = zy end c. y = (t x ), p t = a for i = : n y = y(t x i ) p t = p
More information, applyingl Hospital s Rule again x 0 2 cos(x) xsinx
Lecture 3 We give a couple examples of using L Hospital s Rule: Example 3.. [ (a) Compute x 0 sin(x) x. To put this into a form for L Hospital s Rule we first put it over a common denominator [ x 0 sin(x)
More informationPartial Fractions. Calculus 2 Lia Vas
Calculus Lia Vas Partial Fractions rational function is a quotient of two polynomial functions The method of partial fractions is a general method for evaluating integrals of rational function The idea
More informationFourth Order RK-Method
Fourth Order RK-Method The most commonly used method is Runge-Kutta fourth order method. The fourth order RK-method is y i+1 = y i + 1 6 (k 1 + 2k 2 + 2k 3 + k 4 ), Ordinary Differential Equations (ODE)
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 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 informationExamination paper for TMA4130 Matematikk 4N: SOLUTION
Department of Mathematical Sciences Examination paper for TMA4 Matematikk 4N: SOLUTION Academic contact during examination: Morten Nome Phone: 9 84 97 8 Examination date: December 7 Examination time (from
More informationConvergence of sequences and series
Convergence of sequences and series A sequence f is a map from N the positive integers to a set. We often write the map outputs as f n rather than f(n). Often we just list the outputs in order and leave
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 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 informationLEAST SQUARES APPROXIMATION
LEAST SQUARES APPROXIMATION One more approach to approximating a function f (x) on an interval a x b is to seek an approximation p(x) with a small average error over the interval of approximation. A convenient
More informationSYSTEMS OF NONLINEAR EQUATIONS
SYSTEMS OF NONLINEAR EQUATIONS Widely used in the mathematical modeling of real world phenomena. We introduce some numerical methods for their solution. For better intuition, we examine systems of two
More informationPolynomial Interpolation
Chapter Polynomial Interpolation. Introduction Suppose that we have a two sets of n + real numbers {x i } n+ i= and {y i} n+ i=, and that the x i are strictly increasing: x < x < x 2 < < x n. Interpolation
More informationPower series solutions for 2nd order linear ODE s (not necessarily with constant coefficients) a n z n. n=0
Lecture 22 Power series solutions for 2nd order linear ODE s (not necessarily with constant coefficients) Recall a few facts about power series: a n z n This series in z is centered at z 0. Here z can
More informationCOURSE Numerical integration of functions (continuation) 3.3. The Romberg s iterative generation method
COURSE 7 3. Numerical integration of functions (continuation) 3.3. The Romberg s iterative generation method The presence of derivatives in the remainder difficulties in applicability to practical problems
More informationAy190 Computational Astrophysics
CALIFORNIA INSTITUTE OF TECHNOLOGY Division of Physics, Mathematics, and Astronomy Ay190 Computational Astrophysics Christian D. Ott and Andrew Benson cott@tapir.caltech.edu, abenson@tapir.caltech.edu
More informationNewton s Method and Linear Approximations
Newton s Method and Linear Approximations Curves are tricky. Lines aren t. Newton s Method and Linear Approximations Newton s Method for finding roots Goal: Where is f (x) = 0? f (x) = x 7 + 3x 3 + 7x
More informationFINITE DIFFERENCES. Lecture 1: (a) Operators (b) Forward Differences and their calculations. (c) Backward Differences and their calculations.
FINITE DIFFERENCES Lecture 1: (a) Operators (b) Forward Differences and their calculations. (c) Backward Differences and their calculations. 1. Introduction When a function is known explicitly, it is easy
More information(x + 1)(x 2) = 4. x
dvanced Integration Techniques: Partial Fractions The method of partial fractions can occasionally make it possible to find the integral of a quotient of rational functions. Partial fractions gives us
More informationApproximation theory
Approximation theory Xiaojing Ye, Math & Stat, Georgia State University Spring 2019 Numerical Analysis II Xiaojing Ye, Math & Stat, Georgia State University 1 1 1.3 6 8.8 2 3.5 7 10.1 Least 3squares 4.2
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 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 informationExamples 2: Composite Functions, Piecewise Functions, Partial Fractions
Examples 2: Composite Functions, Piecewise Functions, Partial Fractions September 26, 206 The following are a set of examples to designed to complement a first-year calculus course. objectives are listed
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 informationCalculus II. Monday, March 13th. WebAssign 7 due Friday March 17 Problem Set 6 due Wednesday March 15 Midterm 2 is Monday March 20
Announcements Calculus II Monday, March 13th WebAssign 7 due Friday March 17 Problem Set 6 due Wednesday March 15 Midterm 2 is Monday March 20 Today: Sec. 8.5: Partial Fractions Use partial fractions to
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 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 informationERROR IN LINEAR INTERPOLATION
ERROR IN LINEAR INTERPOLATION Let P 1 (x) be the linear polynomial interpolating f (x) at x 0 and x 1. Assume f (x) is twice continuously differentiable on an interval [a, b] which contains the points
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 information2
1 Notes for Numerical Analysis Math 54 by S. Adjerid Virginia Polytechnic Institute and State University (A Rough Draft) 2 Contents 1 Polynomial Interpolation 5 1.1 Review...............................
More informationMTH4101 CALCULUS II REVISION NOTES. 1. COMPLEX NUMBERS (Thomas Appendix 7 + lecture notes) ax 2 + bx + c = 0. x = b ± b 2 4ac 2a. i = 1.
MTH4101 CALCULUS II REVISION NOTES 1. COMPLEX NUMBERS (Thomas Appendix 7 + lecture notes) 1.1 Introduction Types of numbers (natural, integers, rationals, reals) The need to solve quadratic equations:
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 informationApplied Math for Engineers
Applied Math for Engineers Ming Zhong Lecture 15 March 28, 2018 Ming Zhong (JHU) AMS Spring 2018 1 / 28 Recap Table of Contents 1 Recap 2 Numerical ODEs: Single Step Methods 3 Multistep Methods 4 Method
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 informationn 1 f n 1 c 1 n+1 = c 1 n $ c 1 n 1. After taking logs, this becomes
Root finding: 1 a The points {x n+1, }, {x n, f n }, {x n 1, f n 1 } should be co-linear Say they lie on the line x + y = This gives the relations x n+1 + = x n +f n = x n 1 +f n 1 = Eliminating α and
More informationi x i y i
Department of Mathematics MTL107: Numerical Methods and Computations Exercise Set 8: Approximation-Linear Least Squares Polynomial approximation, Chebyshev Polynomial approximation. 1. Compute the linear
More informationOutline. 1 Numerical Integration. 2 Numerical Differentiation. 3 Richardson Extrapolation
Outline Numerical Integration Numerical Differentiation Numerical Integration Numerical Differentiation 3 Michael T. Heath Scientific Computing / 6 Main Ideas Quadrature based on polynomial interpolation:
More information