Lecture 15. Interpolation II. 2 Piecewise polynomial interpolation Hermite splines
|
|
- Patricia Wheeler
- 5 years ago
- Views:
Transcription
1 Lecture 5 Interpolation II Introduction In te previous lecture we focused primarily on polynomial interpolation of a set of n points. A difficulty we observed is tat wen n is large, our polynomial as to be of ig order, namely n. Unfortunately, ig-order polynomials tend to suffer from wiggle, and tis limits teir practical usefulness for interpolation. In tis lecture we will explore ow we can use polynomials of moderate order to acieve smoot interpolations wile avoiding te problems associate wit ig-order polynomials. Piecewise polynomial interpolation Hermite splines Wat we've previously called linear interpolation is more precisely piecewise-linear interpolation. We don't interpolate te entire set of points wit a single line. Instead, we use different line segments over different intervals between sample points Fig.. Te complete interpolation is built by tying togeter tese lines. In fact te sample points were lines join or tie togeter, x i, y i, i=,,, x n, are appropriately called knots. Fig. : Piecewise interpolation: linear left and cubic rigt. Te sample points at wic pieces join or tie togeter are called knots. Te entire interpolation function is described by y= f x =a i x+bi if xi x< xi + and te a i,b i values are determined by te conditions y i=a i xi +bi y i+=ai x i++bi We can express a linear function in different forms, one of wic migt be more convenient for determining coefficients. We could write te it segment in te form of a st order Taylor series expanded about te point x i, EE Numerical Computing
2 Lecture 5: Interpolation II / y= y i+ y i x x i Tis trivially satisfies y= y i wen x= xi. For x= xi + te requirement y i+= y i+ y i x i+ x i provides te value te y i value y i= y i+ y i x i+ x i Alternately, te Lagrange interpolating polynomial can be written by inspection as y= y i x x i+ x x i + yi + x i x i+ x i+ x i We now extend tis idea of piecewise interpolation to polynomials. Wile a piecewise linear interpolation is continuous, te derivative is clearly not continuous at te sample points. Suppose now tat for eac x i we know bot te function value y i and te function slope y i. Let's build a piecewise polynomial interpolation tat as te specified function and slope values at te knots. Tese polynomial pieces are known as splines. Tis term comes from te practice of bending strips of wood or plastic to form smoot curves, a tecnique often used in sip building and precomputer-era drafting. For eac segment we ave four equations to satisfy, te two endpoint function values and te two endpoint slope values. Our interpolation function must terefore ave four unknown coefficients. Since a rd order polynomial as four coefficients we write Fig. f x=a i x +bi x +ci x+d i if x i x<x i+ In eac interval we ave four unknowns a i,b i, c i, d i satisfying four equations y i=a i xi +bi x i +ci x i+d i y i = ai xi + bi x i +c i y i+=ai x i++bi x i++c i x i++d i + y i+= a i x i++ bi x i++c i We could solve tese four equations in four unknowns directly. Instead we'll apply some bookkeeping to obtain a bit cleaner approac. To simplify analysis over an arbitrary interval x i x<x i+ it's a good idea to form te normalized variable u= x x i x x i = xi + x i i As x varies over x i x x i+, u varies over 0 u. Note tat du = dx i so EE Numerical Computing
3 Lecture 5: Interpolation II / y= dy dy du dy = = dx du dx i du or dy =i y du We now write y=a+b u+c u +d u dy =b+ c u+ d u du Our system of equations is now evaluating te above equations at u=0, y i=a i yi =b y i+=a+b+c+d i yi +=b+ c+ d Te simplification from is significant. Tese can easily be solved to give a= yi b= i y i c= y i i y i + y i+ i y i+ d = y i+i y i yi + +i y i+ and te interpolating cubic is y= y i+i y i u+ yi i y i + y i+ i y i+ u + y i+i y i y i++i y i+u It can be convenient to separate out te various y terms to obtain y= y i u + u +i yi u u +u+ y i+ u u +i yi + u +u A bit of factoring puts tis in a more compact form y= u [ y i + u+i yi u]+u [ y i+ u i yi + u] Te result is te so-called Hermite spline interpolation algoritm. Hermite spline interpolation Given x, y, y, x, y, y,, x n, y n, y n, wit increasing x values For a value x Find x i suc tat x i< x<x i+ Set i =x i+ xi and u= x x i/ i Calculate y= u [ y i + u+i yi u]+u [ y i+ u i yi + u] EE Numerical Computing
4 Lecture 5: Interpolation II 4/ x Fig. : y= e sin x. solid circles: sample points, squares: function values, dased line: linear interpolation, solid line: Hermite-spline interpolation. Derivative values were estimated numerically. Considering te entire interpolation algoritm as a function, we write y=s x. A Scilab function to perform Hermite spline interpolation is given in Appendix as interphermite, and an example of Hermite spline interpolation is sown in Fig.. If we ad samples of te form x i, y i, y i, yi we could find 5t order interpolation polynomials for eac interval and so on, in principle, for any number of known derivatives at eac sample point. If we ave function and derivative values up to d m y /dx m, te two endpoints of eac interval will provide m+ equations. A polynomial wit tis many coefficients as order n= m+. Cubic splines If we know function and derivative values at n points, we can interpolate eac interval wit Hermite splines. Often, owever, we only know te function values and not te derivative values. Tis provides only enoug information to uniquely determine a piecewise-linear interpolation. But te smootness of a piecewise-cubic interpolation is igly desirable, and we would like to find a way to keep tat property even wen we lack derivative information. We will refer to piecewise cubic interpolation witout specific derivative values as cubic splines.. Smootest Hermite spline interpolation One approac would be to treat te y i as unknowns and find te values tat optimize some desirable property of te curve. Smootness is an intuitively appealing property to ave. A smoot curve is one in wic te slope does not cange rapidly. A sudden cange in slope produces a kink in te curve, wic is about as unsmoot as you can get. Terefore te second derivative te rate of cange of te slope sould be small. Let's write te integral of te square of te second derivative of S x as a function of te unknown y i EE Numerical Computing
5 Lecture 5: Interpolation II 5/ values: xn ϕ y, y,, y = [ S x ] dx n 4 x Using te notation S i u = u [ yi + u+ i y i u ]+u [ y i+ u i yi + u] for u as given in, we write 4 as n d ϕ= S u du i dt i= i 0 For evenly space samples were i = we sow in Appendix tat minimizing ϕ leads to te equations y + y = y y, y n+ y n = y n y n 5 y i++4 y i + y i = y i+ y i for i=,,, n 6 for te y i values. Te case of nonuniform samples is similar but a bit messier because we ave to keep track of different values.. Continuity of second derivatives Anoter approac is to require tat not only te first, but also te second derivatives of te interpolation be continuous at te knots x, x,, x n. Tere are n knots, since te endpoints x, x n are not knots no oter pieces connect tere. Continuity of te second derivatives at te knots provide n equations. For uniformly spaced samples tese turn out to be 6. Tis tells us tat te smootest Hermite spline interpolation we derived previously also results in continuous second derivatives at all knots. We ten need two more equations to obtain a unique solution. So-called natural end conditions are obtained by setting S x =0 at te endpoints x, x n, in oter words, we let te interpolation go straigt at bot ends. Tis leads to equations 5. Terefore a cubic spline interpolation wit natural end conditions is precisely te optimally smoot Hermite spline interpolation we derived above. Anoter option is to specify te end-point slopes y, y n. Tis is called te fixed-slope end conditions. If we ave a good estimate of tese slopes ten tis makes sense. Oterwise te coice is arbitrary. Finally we can coose te so-called not-a-knot conditions were we require te tird derivative of te interpolation to be continuous at te first and last knots. At tese knots, terefore, te cubic functions and te first, second and tird derivatives are continuous. But cubics tat agree in tis manner are simply te same cubic; tere is no oter possibility. So wat used to be te first and last knots are no longer knots, ence te name not-a-knot. For te uniformly sampled case tese equations read EE Numerical Computing
6 Lecture 5: Interpolation II 6/ Fig. A case were natural conditions produce a more accurate interpolation tan not-aknot conditions. y +4 y + y = y y and y n +4 y n + y n= y n y n 7 Wic end conditions sould we coose? Te natural conditions are attractive because of teir maximally smoot feature. However, in many cases te not-a-knot conditions provide a more accurate interpolation. It depends on te underlying function f x see Fig. and Fig. 4. Actual functions are not necessarily as smoot as possible! Common practice is to use te nota-knot conditions. In practice te two end conditions produce very similar results except, possibly in te first and last intervals.. Optimal smootness of natural cubic spline interpolation We've spent a lot of time working wit cubic splines. We've sown tat natural cubic splines are te smootest-possible piecewise-cubic interpolation of any set of points. If smootness is so desirable, wy not try piecewise interpolation wit even iger-order polynomials? It turns out tat no oter interpolation is smooter tan a natural cubic spline. Tis rater remarkable result tells us tat, as far as smootness is concerned, cubic splines are te best we can do. Suppose tat S x is te natural cubic spline interpolation of te n samples x i, y i. Let f x be any twice-differentiable function tat also interpolates te samples tis could even be te original function from wic te samples were drawn. Ten one can sow tat xn [S x xn x ] dx [ f x ] dx x Tis is te sense in wic we can say tat no function provides a smooter interpolation of a set of data points tan does te natural cubic spline. However, as sown in Fig. and Fig. 4, smooter is not necessarily better. EE Numerical Computing
7 Lecture 5: Interpolation II 7/ Fig. 4 A case were not-a-knot conditions produce a more accurate interpolation tan natural conditions. 4 Te interp function Scilab/Matlab As for most numerical metods we study, Scilab and Matlab ave built-in functions offering state-of-te-art implementations. Te following command yp = interpx,y,xp,str; //str = 'nearest' or 'linear' or 'spline' Implements one-dimensional interpolation of eiter nearest-neigbor, linear or spline type. For spline interpolation not-a-knot end conditions are used. Here te vectors x,y are sample data, xp is te vector of x values were we want interpolations, and yp is te vector of interpolated Fig. 5: Ten randomly selected points dots. Tick green line: nint-order polynomial. Tin blue line: not-a-knot spline. Tin dased red line: natural spline. EE Numerical Computing
8 Lecture 5: Interpolation II 8/ Fig. 6 Sample data represented as impulses or poles of eigt yi. values. Tis is sufficient for almost all one-dimensional interpolation needs. Wit regards to wiggle, te advantage of splines over ig-order polynomials for interpolation is illustrated in Fig. 5. Notice too tat te difference of te two spline end conditions is significant only in te first and last intervals. Tese data points were cosen at random so tere is no actual underlying f x. However, te spline interpolations certainly appear more realistic. 5 Lanczos convolution interpolation Te metods we ave considered so far apply to an arbitrary set of samples x i, y i. If te x values are uniformly spaced, owever, ten some ideas from signal processing can be applied. We turn to tat now. Let's suppose we ave uniformly spaced samples wit x i= x+i. We can visualize our data as sown in Fig. 6. Imagine an impulse or pole of eigt y i erected vertically wit its base on te ground at location x i. For our purposes convolution is te process of replacing eac impulse wit a common impulse response function, centered at te impulse and scaled by te eigt y i. Te triangle function Λ x is sown in Fig. 7. Fig. 7: Te triangle function EE Numerical Computing Λx
9 Lecture 5: Interpolation II 9/ Fig. 8: Convolution of impulses wit triangle function. Tick gray curve is sum of all triangle functions and interpolates te data points. If we replace eac impulse by a stretced version of te triangle function Λ x / = x / x x > 0 { scaled by y i, ten te sum of tese n f x= y i Λ i= x x i produces te interpolation sown in Fig. 8. We recognize tis as te piecewise linear interpolation of te data wit te addition of linear extrapolations at te two ends. Tis naturally leads us to wonder if using a different impulse response function migt produce a better interpolation. Tinking of x as time and y as te amplitude of an audio signal, tere is a remarkable teorem Fig. 9 Te sinc function. EE Numerical Computing
10 Lecture 5: Interpolation II 0/ due to Nyquist wic says tat provided: te signal from wic te audio samples were draw contains frequency components only witin a limited range, and te sample separation is properly cosen, ten a convolution interpolation using te sinc function pronounced sink will exactly recreate te original function f x. Matematically y= f x = y i sinc i= x x i Te sinc function is sincx = sin π x πx and is plotted in Fig. 9. Note tat sinc 0= and sinc n=0 for n a non-zero integer. Unfortunately te sinc function extends to x ±, altoug te amplitude of te bumps drop off as / x. If we are interpolating many points, we'll ave to add a contribution from eac point. A compromise proposed by Lanczos is to window te sinc function by anoter sinc function to produce te Lanczos kernel L x= { sinc x sinc x a 0 x a x >a were a is typically cosen to be a small integer most often or. Tis is plotted in Fig. 0. Te Lanczos interpolation is i+a y= j=i + a yj L x x j Fig. 0 Te Lanczos kernel for EE Numerical Computing a=,,
11 Lecture 5: Interpolation II / Fig. : Seven data points and Lanczos interpolation. were i is te index suc tat x i x<x i+. In tis expression only te a nearest sample points contribute to te interpolation at a given value of x. An example of Lanczos interpolation wit a= is sown in. 6 Appendix Scilab code 6. Hermite spline interpolation ////////////////////////////////////////////////////////////////////// // interphermite.sci // ,, for pedagogic purposes // Given n samples xi,yi,yi in te column vectors x,y,y // were yi=fxi and yi is te derivative of fx at xi, // interpolate at points xp using Hermite splines. // Note: x and xp values must be in ascending order. ////////////////////////////////////////////////////////////////////// function yp=interphermitex, y, y, xp n = lengtx; m = lengtxp; yp = zerosxp; i = ; //start linear searc at first element for j=:m wile xpj>xi+ //find j so tat xj<=ui<=xj+ i = i+; end = xi+-xi; t = xpj-xi/; ypj = t-^*yi**t++yi**t.. +t^*yi+*-*t+yi+**t-; end endfunction EE Numerical Computing
12 Lecture 5: Interpolation II / 7 Appendix Smootest Hermite spline interpolation Assume we ave n samples x i, y i, and te x values are equally spaced x i= x+i. Let n d ϕ= S t dt i i= 0 dt were [ ] [ S i t =t yi t++ y i t +t y i+ t+ y i+ t ] One can sow tat a computer algebra program elps! 0 d S t dt= y i+ y i y i++ yi y i+ y i +4 y i + y i y i++ y i+ i dt [ ] For a knot <i<n, y j appears in S i and S i. Calling d d 4 w= S t dt+ S t dt 0 dt i 0 dt i we ave [ y [ y ] y ] + y + y y w= y i yi y i+ y j y i y i + y i + y i yi + yi + y i+ i i+ + y i y i+ i i i i+ + y i+ Setting w y i = yi ++4 y i + y i + y i y i+=0 we ave y i++4 y i + y i = y i+ y i 8 For te first interval { [ ] } y y y + y y y + y + y y + y =0 y gives us y + y = y y 9 y n+ y n = y n y n 0 wile for te last interval we find EE Numerical Computing
Polynomial Interpolation
Capter 4 Polynomial Interpolation In tis capter, we consider te important problem of approximatinga function fx, wose values at a set of distinct points x, x, x,, x n are known, by a polynomial P x suc
More informationPolynomial Interpolation
Capter 4 Polynomial Interpolation In tis capter, we consider te important problem of approximating a function f(x, wose values at a set of distinct points x, x, x 2,,x n are known, by a polynomial P (x
More information2.8 The Derivative as a Function
.8 Te Derivative as a Function Typically, we can find te derivative of a function f at many points of its domain: Definition. Suppose tat f is a function wic is differentiable at every point of an open
More informationNumerical Differentiation
Numerical Differentiation Finite Difference Formulas for te first derivative (Using Taylor Expansion tecnique) (section 8.3.) Suppose tat f() = g() is a function of te variable, and tat as 0 te function
More informationLecture 21. Numerical differentiation. f ( x+h) f ( x) h h
Lecture Numerical differentiation Introduction We can analytically calculate te derivative of any elementary function, so tere migt seem to be no motivation for calculating derivatives numerically. However
More informationHow to Find the Derivative of a Function: Calculus 1
Introduction How to Find te Derivative of a Function: Calculus 1 Calculus is not an easy matematics course Te fact tat you ave enrolled in suc a difficult subject indicates tat you are interested in te
More informationConsider a function f we ll specify which assumptions we need to make about it in a minute. Let us reformulate the integral. 1 f(x) dx.
Capter 2 Integrals as sums and derivatives as differences We now switc to te simplest metods for integrating or differentiating a function from its function samples. A careful study of Taylor expansions
More informationlecture 26: Richardson extrapolation
43 lecture 26: Ricardson extrapolation 35 Ricardson extrapolation, Romberg integration Trougout numerical analysis, one encounters procedures tat apply some simple approximation (eg, linear interpolation)
More informationAverage Rate of Change
Te Derivative Tis can be tougt of as an attempt to draw a parallel (pysically and metaporically) between a line and a curve, applying te concept of slope to someting tat isn't actually straigt. Te slope
More informationSECTION 3.2: DERIVATIVE FUNCTIONS and DIFFERENTIABILITY
(Section 3.2: Derivative Functions and Differentiability) 3.2.1 SECTION 3.2: DERIVATIVE FUNCTIONS and DIFFERENTIABILITY LEARNING OBJECTIVES Know, understand, and apply te Limit Definition of te Derivative
More informationMVT and Rolle s Theorem
AP Calculus CHAPTER 4 WORKSHEET APPLICATIONS OF DIFFERENTIATION MVT and Rolle s Teorem Name Seat # Date UNLESS INDICATED, DO NOT USE YOUR CALCULATOR FOR ANY OF THESE QUESTIONS In problems 1 and, state
More informationRecall from our discussion of continuity in lecture a function is continuous at a point x = a if and only if
Computational Aspects of its. Keeping te simple simple. Recall by elementary functions we mean :Polynomials (including linear and quadratic equations) Eponentials Logaritms Trig Functions Rational Functions
More informationHOMEWORK HELP 2 FOR MATH 151
HOMEWORK HELP 2 FOR MATH 151 Here we go; te second round of omework elp. If tere are oters you would like to see, let me know! 2.4, 43 and 44 At wat points are te functions f(x) and g(x) = xf(x)continuous,
More information2.11 That s So Derivative
2.11 Tat s So Derivative Introduction to Differential Calculus Just as one defines instantaneous velocity in terms of average velocity, we now define te instantaneous rate of cange of a function at a point
More informationExam 1 Review Solutions
Exam Review Solutions Please also review te old quizzes, and be sure tat you understand te omework problems. General notes: () Always give an algebraic reason for your answer (graps are not sufficient),
More information2.1 THE DEFINITION OF DERIVATIVE
2.1 Te Derivative Contemporary Calculus 2.1 THE DEFINITION OF DERIVATIVE 1 Te grapical idea of a slope of a tangent line is very useful, but for some uses we need a more algebraic definition of te derivative
More informationSin, Cos and All That
Sin, Cos and All Tat James K. Peterson Department of Biological Sciences and Department of Matematical Sciences Clemson University Marc 9, 2017 Outline Sin, Cos and all tat! A New Power Rule Derivatives
More informationOrder of Accuracy. ũ h u Ch p, (1)
Order of Accuracy 1 Terminology We consider a numerical approximation of an exact value u. Te approximation depends on a small parameter, wic can be for instance te grid size or time step in a numerical
More informationLIMITS AND DERIVATIVES CONDITIONS FOR THE EXISTENCE OF A LIMIT
LIMITS AND DERIVATIVES Te limit of a function is defined as te value of y tat te curve approaces, as x approaces a particular value. Te limit of f (x) as x approaces a is written as f (x) approaces, as
More informationThe derivative function
Roberto s Notes on Differential Calculus Capter : Definition of derivative Section Te derivative function Wat you need to know already: f is at a point on its grap and ow to compute it. Wat te derivative
More information1 The concept of limits (p.217 p.229, p.242 p.249, p.255 p.256) 1.1 Limits Consider the function determined by the formula 3. x since at this point
MA00 Capter 6 Calculus and Basic Linear Algebra I Limits, Continuity and Differentiability Te concept of its (p.7 p.9, p.4 p.49, p.55 p.56). Limits Consider te function determined by te formula f Note
More informationCopyright c 2008 Kevin Long
Lecture 4 Numerical solution of initial value problems Te metods you ve learned so far ave obtained closed-form solutions to initial value problems. A closedform solution is an explicit algebriac formula
More information232 Calculus and Structures
3 Calculus and Structures CHAPTER 17 JUSTIFICATION OF THE AREA AND SLOPE METHODS FOR EVALUATING BEAMS Calculus and Structures 33 Copyrigt Capter 17 JUSTIFICATION OF THE AREA AND SLOPE METHODS 17.1 THE
More informationMath 102 TEST CHAPTERS 3 & 4 Solutions & Comments Fall 2006
Mat 102 TEST CHAPTERS 3 & 4 Solutions & Comments Fall 2006 f(x+) f(x) 10 1. For f(x) = x 2 + 2x 5, find ))))))))) and simplify completely. NOTE: **f(x+) is NOT f(x)+! f(x+) f(x) (x+) 2 + 2(x+) 5 ( x 2
More informationTime (hours) Morphine sulfate (mg)
Mat Xa Fall 2002 Review Notes Limits and Definition of Derivative Important Information: 1 According to te most recent information from te Registrar, te Xa final exam will be eld from 9:15 am to 12:15
More informationLecture XVII. Abstract We introduce the concept of directional derivative of a scalar function and discuss its relation with the gradient operator.
Lecture XVII Abstract We introduce te concept of directional derivative of a scalar function and discuss its relation wit te gradient operator. Directional derivative and gradient Te directional derivative
More informationMath 212-Lecture 9. For a single-variable function z = f(x), the derivative is f (x) = lim h 0
3.4: Partial Derivatives Definition Mat 22-Lecture 9 For a single-variable function z = f(x), te derivative is f (x) = lim 0 f(x+) f(x). For a function z = f(x, y) of two variables, to define te derivatives,
More informationIntegral Calculus, dealing with areas and volumes, and approximate areas under and between curves.
Calculus can be divided into two ke areas: Differential Calculus dealing wit its, rates of cange, tangents and normals to curves, curve sketcing, and applications to maima and minima problems Integral
More informationLECTURE 14 NUMERICAL INTEGRATION. Find
LECTURE 14 NUMERCAL NTEGRATON Find b a fxdx or b a vx ux fx ydy dx Often integration is required. However te form of fx may be suc tat analytical integration would be very difficult or impossible. Use
More informationFunction Composition and Chain Rules
Function Composition and s James K. Peterson Department of Biological Sciences and Department of Matematical Sciences Clemson University Marc 8, 2017 Outline 1 Function Composition and Continuity 2 Function
More informationMathematics 5 Worksheet 11 Geometry, Tangency, and the Derivative
Matematics 5 Workseet 11 Geometry, Tangency, and te Derivative Problem 1. Find te equation of a line wit slope m tat intersects te point (3, 9). Solution. Te equation for a line passing troug a point (x
More informationPrecalculus Test 2 Practice Questions Page 1. Note: You can expect other types of questions on the test than the ones presented here!
Precalculus Test 2 Practice Questions Page Note: You can expect oter types of questions on te test tan te ones presented ere! Questions Example. Find te vertex of te quadratic f(x) = 4x 2 x. Example 2.
More information158 Calculus and Structures
58 Calculus and Structures CHAPTER PROPERTIES OF DERIVATIVES AND DIFFERENTIATION BY THE EASY WAY. Calculus and Structures 59 Copyrigt Capter PROPERTIES OF DERIVATIVES. INTRODUCTION In te last capter you
More informationSection 2.7 Derivatives and Rates of Change Part II Section 2.8 The Derivative as a Function. at the point a, to be. = at time t = a is
Mat 180 www.timetodare.com Section.7 Derivatives and Rates of Cange Part II Section.8 Te Derivative as a Function Derivatives ( ) In te previous section we defined te slope of te tangent to a curve wit
More informationDifferentiation in higher dimensions
Capter 2 Differentiation in iger dimensions 2.1 Te Total Derivative Recall tat if f : R R is a 1-variable function, and a R, we say tat f is differentiable at x = a if and only if te ratio f(a+) f(a) tends
More informationHigher Derivatives. Differentiable Functions
Calculus 1 Lia Vas Higer Derivatives. Differentiable Functions Te second derivative. Te derivative itself can be considered as a function. Te instantaneous rate of cange of tis function is te second derivative.
More information= 0 and states ''hence there is a stationary point'' All aspects of the proof dx must be correct (c)
Paper 1: Pure Matematics 1 Mark Sceme 1(a) (i) (ii) d d y 3 1x 4x x M1 A1 d y dx 1.1b 1.1b 36x 48x A1ft 1.1b Substitutes x = into teir dx (3) 3 1 4 Sows d y 0 and states ''ence tere is a stationary point''
More informationMATH745 Fall MATH745 Fall
MATH745 Fall 5 MATH745 Fall 5 INTRODUCTION WELCOME TO MATH 745 TOPICS IN NUMERICAL ANALYSIS Instructor: Dr Bartosz Protas Department of Matematics & Statistics Email: bprotas@mcmasterca Office HH 36, Ext
More informationThe Priestley-Chao Estimator
Te Priestley-Cao Estimator In tis section we will consider te Pristley-Cao estimator of te unknown regression function. It is assumed tat we ave a sample of observations (Y i, x i ), i = 1,..., n wic are
More informationExcursions in Computing Science: Week v Milli-micro-nano-..math Part II
Excursions in Computing Science: Week v Milli-micro-nano-..mat Part II T. H. Merrett McGill University, Montreal, Canada June, 5 I. Prefatory Notes. Cube root of 8. Almost every calculator as a square-root
More informationREVIEW LAB ANSWER KEY
REVIEW LAB ANSWER KEY. Witout using SN, find te derivative of eac of te following (you do not need to simplify your answers): a. f x 3x 3 5x x 6 f x 3 3x 5 x 0 b. g x 4 x x x notice te trick ere! x x g
More informationMAT 145. Type of Calculator Used TI-89 Titanium 100 points Score 100 possible points
MAT 15 Test #2 Name Solution Guide Type of Calculator Used TI-89 Titanium 100 points Score 100 possible points Use te grap of a function sown ere as you respond to questions 1 to 8. 1. lim f (x) 0 2. lim
More informationCHAPTER (A) When x = 2, y = 6, so f( 2) = 6. (B) When y = 4, x can equal 6, 2, or 4.
SECTION 3-1 101 CHAPTER 3 Section 3-1 1. No. A correspondence between two sets is a function only if eactly one element of te second set corresponds to eac element of te first set. 3. Te domain of a function
More information(a) At what number x = a does f have a removable discontinuity? What value f(a) should be assigned to f at x = a in order to make f continuous at a?
Solutions to Test 1 Fall 016 1pt 1. Te grap of a function f(x) is sown at rigt below. Part I. State te value of eac limit. If a limit is infinite, state weter it is or. If a limit does not exist (but is
More informationINTRODUCTION TO CALCULUS LIMITS
Calculus can be divided into two ke areas: INTRODUCTION TO CALCULUS Differential Calculus dealing wit its, rates of cange, tangents and normals to curves, curve sketcing, and applications to maima and
More information5 Ordinary Differential Equations: Finite Difference Methods for Boundary Problems
5 Ordinary Differential Equations: Finite Difference Metods for Boundary Problems Read sections 10.1, 10.2, 10.4 Review questions 10.1 10.4, 10.8 10.9, 10.13 5.1 Introduction In te previous capters we
More informationTaylor Series and the Mean Value Theorem of Derivatives
1 - Taylor Series and te Mean Value Teorem o Derivatives Te numerical solution o engineering and scientiic problems described by matematical models oten requires solving dierential equations. Dierential
More informationCombining functions: algebraic methods
Combining functions: algebraic metods Functions can be added, subtracted, multiplied, divided, and raised to a power, just like numbers or algebra expressions. If f(x) = x 2 and g(x) = x + 2, clearly f(x)
More informationChapter 5 FINITE DIFFERENCE METHOD (FDM)
MEE7 Computer Modeling Tecniques in Engineering Capter 5 FINITE DIFFERENCE METHOD (FDM) 5. Introduction to FDM Te finite difference tecniques are based upon approximations wic permit replacing differential
More informationCubic Functions: Local Analysis
Cubic function cubing coefficient Capter 13 Cubic Functions: Local Analysis Input-Output Pairs, 378 Normalized Input-Output Rule, 380 Local I-O Rule Near, 382 Local Grap Near, 384 Types of Local Graps
More informationNUMERICAL DIFFERENTIATION. James T. Smith San Francisco State University. In calculus classes, you compute derivatives algebraically: for example,
NUMERICAL DIFFERENTIATION James T Smit San Francisco State University In calculus classes, you compute derivatives algebraically: for example, f( x) = x + x f ( x) = x x Tis tecnique requires your knowing
More information1watt=1W=1kg m 2 /s 3
Appendix A Matematics Appendix A.1 Units To measure a pysical quantity, you need a standard. Eac pysical quantity as certain units. A unit is just a standard we use to compare, e.g. a ruler. In tis laboratory
More informationSection 3: The Derivative Definition of the Derivative
Capter 2 Te Derivative Business Calculus 85 Section 3: Te Derivative Definition of te Derivative Returning to te tangent slope problem from te first section, let's look at te problem of finding te slope
More informationDifferential Calculus (The basics) Prepared by Mr. C. Hull
Differential Calculus Te basics) A : Limits In tis work on limits, we will deal only wit functions i.e. tose relationsips in wic an input variable ) defines a unique output variable y). Wen we work wit
More informationChapter 8. Numerical Solution of Ordinary Differential Equations. Module No. 2. Predictor-Corrector Methods
Numerical Analysis by Dr. Anita Pal Assistant Professor Department of Matematics National Institute of Tecnology Durgapur Durgapur-7109 email: anita.buie@gmail.com 1 . Capter 8 Numerical Solution of Ordinary
More information2.3 More Differentiation Patterns
144 te derivative 2.3 More Differentiation Patterns Polynomials are very useful, but tey are not te only functions we need. Tis section uses te ideas of te two previous sections to develop tecniques for
More informationA = h w (1) Error Analysis Physics 141
Introduction In all brances of pysical science and engineering one deals constantly wit numbers wic results more or less directly from experimental observations. Experimental observations always ave inaccuracies.
More informationChapter 4: Numerical Methods for Common Mathematical Problems
1 Capter 4: Numerical Metods for Common Matematical Problems Interpolation Problem: Suppose we ave data defined at a discrete set of points (x i, y i ), i = 0, 1,..., N. Often it is useful to ave a smoot
More information5.1 introduction problem : Given a function f(x), find a polynomial approximation p n (x).
capter 5 : polynomial approximation and interpolation 5 introduction problem : Given a function f(x), find a polynomial approximation p n (x) Z b Z application : f(x)dx b p n(x)dx, a a one solution : Te
More informationERROR BOUNDS FOR THE METHODS OF GLIMM, GODUNOV AND LEVEQUE BRADLEY J. LUCIER*
EO BOUNDS FO THE METHODS OF GLIMM, GODUNOV AND LEVEQUE BADLEY J. LUCIE* Abstract. Te expected error in L ) attimet for Glimm s sceme wen applied to a scalar conservation law is bounded by + 2 ) ) /2 T
More informationTHE IDEA OF DIFFERENTIABILITY FOR FUNCTIONS OF SEVERAL VARIABLES Math 225
THE IDEA OF DIFFERENTIABILITY FOR FUNCTIONS OF SEVERAL VARIABLES Mat 225 As we ave seen, te definition of derivative for a Mat 111 function g : R R and for acurveγ : R E n are te same, except for interpretation:
More informationJournal of Computational and Applied Mathematics
Journal of Computational and Applied Matematics 94 (6) 75 96 Contents lists available at ScienceDirect Journal of Computational and Applied Matematics journal omepage: www.elsevier.com/locate/cam Smootness-Increasing
More informationSolution. Solution. f (x) = (cos x)2 cos(2x) 2 sin(2x) 2 cos x ( sin x) (cos x) 4. f (π/4) = ( 2/2) ( 2/2) ( 2/2) ( 2/2) 4.
December 09, 20 Calculus PracticeTest s Name: (4 points) Find te absolute extrema of f(x) = x 3 0 on te interval [0, 4] Te derivative of f(x) is f (x) = 3x 2, wic is zero only at x = 0 Tus we only need
More informationDerivatives of Exponentials
mat 0 more on derivatives: day 0 Derivatives of Eponentials Recall tat DEFINITION... An eponential function as te form f () =a, were te base is a real number a > 0. Te domain of an eponential function
More informationAnalytic Functions. Differentiable Functions of a Complex Variable
Analytic Functions Differentiable Functions of a Complex Variable In tis capter, we sall generalize te ideas for polynomials power series of a complex variable we developed in te previous capter to general
More informationPractice Problem Solutions: Exam 1
Practice Problem Solutions: Exam 1 1. (a) Algebraic Solution: Te largest term in te numerator is 3x 2, wile te largest term in te denominator is 5x 2 3x 2 + 5. Tus lim x 5x 2 2x 3x 2 x 5x 2 = 3 5 Numerical
More informationClick here to see an animation of the derivative
Differentiation Massoud Malek Derivative Te concept of derivative is at te core of Calculus; It is a very powerful tool for understanding te beavior of matematical functions. It allows us to optimize functions,
More information4. The slope of the line 2x 7y = 8 is (a) 2/7 (b) 7/2 (c) 2 (d) 2/7 (e) None of these.
Mat 11. Test Form N Fall 016 Name. Instructions. Te first eleven problems are wort points eac. Te last six problems are wort 5 points eac. For te last six problems, you must use relevant metods of algebra
More informationChapter 1. Density Estimation
Capter 1 Density Estimation Let X 1, X,..., X n be observations from a density f X x. Te aim is to use only tis data to obtain an estimate ˆf X x of f X x. Properties of f f X x x, Parametric metods f
More informationDerivatives of trigonometric functions
Derivatives of trigonometric functions 2 October 207 Introuction Toay we will ten iscuss te erivates of te si stanar trigonometric functions. Of tese, te most important are sine an cosine; te erivatives
More informationPreface. Here are a couple of warnings to my students who may be here to get a copy of what happened on a day that you missed.
Preface Here are my online notes for my course tat I teac ere at Lamar University. Despite te fact tat tese are my class notes, tey sould be accessible to anyone wanting to learn or needing a refreser
More informationWYSE Academic Challenge 2004 Sectional Mathematics Solution Set
WYSE Academic Callenge 00 Sectional Matematics Solution Set. Answer: B. Since te equation can be written in te form x + y, we ave a major 5 semi-axis of lengt 5 and minor semi-axis of lengt. Tis means
More informationf a h f a h h lim lim
Te Derivative Te derivative of a function f at a (denoted f a) is f a if tis it exists. An alternative way of defining f a is f a x a fa fa fx fa x a Note tat te tangent line to te grap of f at te point
More information2011 Fermat Contest (Grade 11)
Te CENTRE for EDUCATION in MATHEMATICS and COMPUTING 011 Fermat Contest (Grade 11) Tursday, February 4, 011 Solutions 010 Centre for Education in Matematics and Computing 011 Fermat Contest Solutions Page
More information1 Calculus. 1.1 Gradients and the Derivative. Q f(x+h) f(x)
Calculus. Gradients and te Derivative Q f(x+) δy P T δx R f(x) 0 x x+ Let P (x, f(x)) and Q(x+, f(x+)) denote two points on te curve of te function y = f(x) and let R denote te point of intersection of
More informationMATH1151 Calculus Test S1 v2a
MATH5 Calculus Test 8 S va January 8, 5 Tese solutions were written and typed up by Brendan Trin Please be etical wit tis resource It is for te use of MatSOC members, so do not repost it on oter forums
More informationThe Verlet Algorithm for Molecular Dynamics Simulations
Cemistry 380.37 Fall 2015 Dr. Jean M. Standard November 9, 2015 Te Verlet Algoritm for Molecular Dynamics Simulations Equations of motion For a many-body system consisting of N particles, Newton's classical
More informationExercises for numerical differentiation. Øyvind Ryan
Exercises for numerical differentiation Øyvind Ryan February 25, 2013 1. Mark eac of te following statements as true or false. a. Wen we use te approximation f (a) (f (a +) f (a))/ on a computer, we can
More informationName: Answer Key No calculators. Show your work! 1. (21 points) All answers should either be,, a (finite) real number, or DNE ( does not exist ).
Mat - Final Exam August 3 rd, Name: Answer Key No calculators. Sow your work!. points) All answers sould eiter be,, a finite) real number, or DNE does not exist ). a) Use te grap of te function to evaluate
More informationIEOR 165 Lecture 10 Distribution Estimation
IEOR 165 Lecture 10 Distribution Estimation 1 Motivating Problem Consider a situation were we ave iid data x i from some unknown distribution. One problem of interest is estimating te distribution tat
More information1. Consider the trigonometric function f(t) whose graph is shown below. Write down a possible formula for f(t).
. Consider te trigonometric function f(t) wose grap is sown below. Write down a possible formula for f(t). Tis function appears to be an odd, periodic function tat as been sifted upwards, so we will use
More informationBob Brown Math 251 Calculus 1 Chapter 3, Section 1 Completed 1 CCBC Dundalk
Bob Brown Mat 251 Calculus 1 Capter 3, Section 1 Completed 1 Te Tangent Line Problem Te idea of a tangent line first arises in geometry in te context of a circle. But before we jump into a discussion of
More informationlim 1 lim 4 Precalculus Notes: Unit 10 Concepts of Calculus
Syllabus Objectives: 1.1 Te student will understand and apply te concept of te limit of a function at given values of te domain. 1. Te student will find te limit of a function at given values of te domain.
More informationPOLYNOMIAL AND SPLINE ESTIMATORS OF THE DISTRIBUTION FUNCTION WITH PRESCRIBED ACCURACY
APPLICATIONES MATHEMATICAE 36, (29), pp. 2 Zbigniew Ciesielski (Sopot) Ryszard Zieliński (Warszawa) POLYNOMIAL AND SPLINE ESTIMATORS OF THE DISTRIBUTION FUNCTION WITH PRESCRIBED ACCURACY Abstract. Dvoretzky
More informationSection 15.6 Directional Derivatives and the Gradient Vector
Section 15.6 Directional Derivatives and te Gradient Vector Finding rates of cange in different directions Recall tat wen we first started considering derivatives of functions of more tan one variable,
More informationThe total error in numerical differentiation
AMS 147 Computational Metods and Applications Lecture 08 Copyrigt by Hongyun Wang, UCSC Recap: Loss of accuracy due to numerical cancellation A B 3, 3 ~10 16 In calculating te difference between A and
More informationContinuity and Differentiability Worksheet
Continuity and Differentiability Workseet (Be sure tat you can also do te grapical eercises from te tet- Tese were not included below! Typical problems are like problems -3, p. 6; -3, p. 7; 33-34, p. 7;
More informationThe Laplace equation, cylindrically or spherically symmetric case
Numerisce Metoden II, 7 4, und Übungen, 7 5 Course Notes, Summer Term 7 Some material and exercises Te Laplace equation, cylindrically or sperically symmetric case Electric and gravitational potential,
More informationA h u h = f h. 4.1 The CoarseGrid SystemandtheResidual Equation
Capter Grid Transfer Remark. Contents of tis capter. Consider a grid wit grid size and te corresponding linear system of equations A u = f. Te summary given in Section 3. leads to te idea tat tere migt
More informationLesson 6: The Derivative
Lesson 6: Te Derivative Def. A difference quotient for a function as te form f(x + ) f(x) (x + ) x f(x + x) f(x) (x + x) x f(a + ) f(a) (a + ) a Notice tat a difference quotient always as te form of cange
More informationSection 2: The Derivative Definition of the Derivative
Capter 2 Te Derivative Applied Calculus 80 Section 2: Te Derivative Definition of te Derivative Suppose we drop a tomato from te top of a 00 foot building and time its fall. Time (sec) Heigt (ft) 0.0 00
More informationINTRODUCTION AND MATHEMATICAL CONCEPTS
INTODUCTION ND MTHEMTICL CONCEPTS PEVIEW Tis capter introduces you to te basic matematical tools for doing pysics. You will study units and converting between units, te trigonometric relationsips of sine,
More informationOne-Sided Position-Dependent Smoothness-Increasing Accuracy-Conserving (SIAC) Filtering Over Uniform and Non-uniform Meshes
DOI 10.1007/s10915-014-9946-6 One-Sided Position-Dependent Smootness-Increasing Accuracy-Conserving (SIAC) Filtering Over Uniform and Non-uniform Meses JenniferK.Ryan Xiaozou Li Robert M. Kirby Kees Vuik
More information3.1 Extreme Values of a Function
.1 Etreme Values of a Function Section.1 Notes Page 1 One application of te derivative is finding minimum and maimum values off a grap. In precalculus we were only able to do tis wit quadratics by find
More informationPolynomials 3: Powers of x 0 + h
near small binomial Capter 17 Polynomials 3: Powers of + Wile it is easy to compute wit powers of a counting-numerator, it is a lot more difficult to compute wit powers of a decimal-numerator. EXAMPLE
More informationSFU UBC UNBC Uvic Calculus Challenge Examination June 5, 2008, 12:00 15:00
SFU UBC UNBC Uvic Calculus Callenge Eamination June 5, 008, :00 5:00 Host: SIMON FRASER UNIVERSITY First Name: Last Name: Scool: Student signature INSTRUCTIONS Sow all your work Full marks are given only
More informationA MONTE CARLO ANALYSIS OF THE EFFECTS OF COVARIANCE ON PROPAGATED UNCERTAINTIES
A MONTE CARLO ANALYSIS OF THE EFFECTS OF COVARIANCE ON PROPAGATED UNCERTAINTIES Ronald Ainswort Hart Scientific, American Fork UT, USA ABSTRACT Reports of calibration typically provide total combined uncertainties
More informationPoisson Equation in Sobolev Spaces
Poisson Equation in Sobolev Spaces OcMountain Dayligt Time. 6, 011 Today we discuss te Poisson equation in Sobolev spaces. It s existence, uniqueness, and regularity. Weak Solution. u = f in, u = g on
More informationMath 312 Lecture Notes Modeling
Mat 3 Lecture Notes Modeling Warren Weckesser Department of Matematics Colgate University 5 7 January 006 Classifying Matematical Models An Example We consider te following scenario. During a storm, a
More informationIntroduction to Derivatives
Introduction to Derivatives 5-Minute Review: Instantaneous Rates and Tangent Slope Recall te analogy tat we developed earlier First we saw tat te secant slope of te line troug te two points (a, f (a))
More information