Cubic Splines. Antony Jameson. Department of Aeronautics and Astronautics, Stanford University, Stanford, California, 94305
|
|
- Blaze Baker
- 6 years ago
- Views:
Transcription
1 Cubic Splines Antony Jameson Department of Aeronautics and Astronautics, Stanford University, Stanford, California, References on splines 1. J. H. Ahlberg, E. N. Nilson, J. H. Walsh. Theory of splines and their applications, Academic Press, I. J. Schoenberg (ed.). Spline functions, (Symposium at U. Wisconsin, 1969), Academic Press, I. J. Schoenberg. Quart. Appl. Math. 4, 1946, pp , pp C. H. Reinsch. Numer. Math. 10, 1967, pp Smoothing by spline functions. 5. Schulz, Spline Analysis. 1
2 2 Definition of spline A spline is a piecewise polynomial in which the coefficients of each polynomial are fixed between knots or joints. Figure 1: Typically cubics are used. Then the coefficients are chosen to match the function and its first and second derivatives at eacoint. There remain one free condition at each end, or two conditions at one end. However, using only starting conditions the spline is unstable. In general with n th degree polynomials one can obtain continuity up to the n 1 derivative. The most common spline is a cubic spline. Then the spline function y(x) satisfies y (4) (x) = 0, y(3)(x) = const, y (x) = a(x) + h. But for a beam between simple supports y (x) = M(x) EI where M(x) varies linearly. Thus a spline is the curve obtained from a draughtsman s spline. 2
3 Figure 2: Draughtsman s spline 3 Equations of cubic spline Let data be given at x 0, x 1, x n with values y 0, y 1, y n. Let S(x) bet the spline. Let M j = S (x j ), = x j x j 1 be moment at j th point. Then between x j 1 and x j S (x) = M j 1 x j x + M j x x j 1 But S (x) = M j 1 (x j x) M j (x x j 1 ) A (x j x) 3 (x x j 1 ) 3 S(x) = M j 1 + M j + Ax + B 6 6 S(x j 1 ) = y j 1, S(x j ) = y j. We obtain ( (x j x) 3 (x x j 1 ) 3 h 2 ) ( j xj x h 2 ) j x xj S(x) = M j 1 +M j + y j 1 M j 1 + y j + M j
4 S (x) = M j 1 (x j x) M j (x x j 1 ) y j y j 1 M j M j 1 6 We now have continuity of S (x) and S(x) at x j. We also require continuity of S (x). Now S (x j ) = 6 M j M j + y j y j 1 S (x + j ) = +1 M j M j+1 + y j+1 y j +1 Equating these gives equations for the M j at every interior node: 6 M j M j M j+1 = y j+1 y j +1 y j y j 1 To complete the system we need 2 additional equations. Setting S (x 0 ) = M 0 = 0, S (x n ) = M n = 0, corresponds to simple supports at the ends of the beam. Better accuracy is obtained if S (x 0 ) and S (x n ) are known. Then from the equation for S (x n ) we have h 1 3 M 0 + h 1 6 M 1 = y 1 y 0 h 1 y 0 h n 6 M n 1 + h n 3 M n = y n y n y n 1 h n If neither S (x) nor S (x) are known at the ends, one may set S (x) = 0, or M 1 M 0 = 0 M n 1 M n = 0 Alternatively one may use some linear combination of these end conditions. In any 4
5 case one obtains a tridiagonal set of equations which may be written 2 λ 0 µ 1 2 λ 1 µ n 1 2 λ n 1 µ n 2 M 0 M 1 M n 1 M n = d 0 d 1 d n 1 d n where at interior points λ j = , µ j = 1 λ j d j = 6 (y j+1 y j )/+1 (y j y j 1 )/ + +1 and λ 0, µ n, d 0, d n depend on the end conditions. These equations are diagonally dominant, leading to stable behaviour. Spline equations for first derivatives Let m j = S (x j ) Then S(x) = m j 1 (x j x) 2 (x x j 1 ) h 2 j +y j 1 (x j x) 2 [2(x x j 1 ) + ] h 3 j 5 m j (x x j 1 ) 2 (x j x) h 2 j + y j (x x j 1 ) 2 [2(x j x) + ] h 3 j
6 S (x j x)(2x j 1 + x j 3x) (x x j 1 )(2x j x j 1 3x) (x) = m j 1 m h 2 j j h 2 j +6 y j y j 1 (x h 3 j x)(x x j 1 ) j S (x) = 2m j 1 2x j + x j 1 3x h 2 j 2m j 2x j 1 + x j 3x) h 2 j +6 y j y j 1 (x h 3 j + x j 1 2x) j S (x j ) = 2m j 1 + 4m j 6 y j y j 1 h 2 j S (x + j ) = 4m j 2m j y j+1 y j h 2 j+1 Equating these gives ( 1 1 m j ) m j + 1 m j+1 = 3 y j y j y j+1 y j h 2 j h 2 j+1 Then 2 µ 0 λ 1 2 µ 1 λ N 1 2 µ N 1 λ N 2 m 0 m 1 m N 1 m N = c 0 c 1 c N 1 c N where at interior points c j = 3λ j y j y j 1 + 3µ j y j+1 y j +1 6
7 The end condition S = 0 gives µ 0 = 1, c 0 = 3 y j y 0 λ N = 1, c N = 3 y N y N 1 h N 1 4 Minimum curvature property of splines Consider all functions f(x) such that f(x i ) = y i Then the spline S(x) with S (x 0 ) = S (x n ) = 0 is that function f(x) which minimize: xn x 0 f (x) 2 dx (Ahlberg, p. 76) Proof = xn xn x 0 f (x) 2 dx 2 x 0 (f (x) S (x)) 2 dx (f (x) S (x)) S (x)dx S (x) 2 dx The middle term is n i=1 xi Integrating by parts it becomes x i 1 (f (x) S (x)) S (x)dx n n [(f (x) S (x)) S (x)] i i 1 i=1 i=1 xi x i 1 (f (x) S (x)) S (x)dx 7
8 S (x) is constant in each segment, and can be taken out of each integral in the second term. But xi which vanishes because x i 1 (f (x) S (x)) dx = [f(x) S(x)] i i 1 f(x i ) = S(x i ) = y i, i = 0, 1,... n Also the contributions to the sum in the first term vanish at every interior node, leaving (f (x n ) S (x n )) S (x n ) (f (x 0 ) S (x 0 )) S (x 0 ) Thus xn x 0 (f (x) S (x)) S (x)dx = 0 if either S (x) = 0 at each end or f (x) = S (x). Consequently xn x 0 f (x) 2 dx xn xn S (x) 2 dx = (f (x) S (x)) 2 dx > 0 x 0 x 0 unless f(x) = S(x). 8
9 5 Errors in spline approximation If the maximum interval is h then f(x) S(x) h 4 f (x) S (x) h 3 f (x) S (x) h 2 f (x) S (x) h This is proved by Ahlberg, Nilson and Walsh, p. 29, provided that f (4) is continuous. To study convergence properties of splines we need an estimate of the norm of the inverse of the matrix B of the equations for the spline. Taking x = max x i the induced norm of B is Bx B = sup x 0 x = max i n b ij j=1 9
10 With a diagonally dominant matrix we have n y = Bx = max b ij x j i j=1 n b ij x j j=1 where x = x k, or y { b kk j k { min i b kj b ii j i } b ij x } x Thus B 1 B 1 y = sup y 0 y [ min i { b i i j i b ij }] 1 or for the spline matrix B 1 max [ (2 λ 0 ) 1, (2 λ n ) 1, 1 ] Now let δ = max and let max δ β <. Then the spline has the following convergence propety as the maximum mesh spacing 10
11 δ is reduced: Theorem A: (Ahlberg, p. 22) Let f(x) be continuous on [a, b] and let max [ λ 0, µ n ] < 2. Then if S(x) is a spline fit to f(x) with x 0 = a, x n = b, then f(x) S(x) 0 as δ 0. Also if f(x) f(x ) κ x x α, 0 α 1 then f(x) S(x) = O (δ α ). Proof: On [x j 1, x j ] the spline equations can be arranged as S(x) f(x) = M j 1(x j x) [ ] (x j x) 2 h 2 j + M [ j(x x j 1 ) (x xj 1 ) 2 hj] f j + f j 1 2 f(x) + (f j f j 1 ) x j + x j 1 2x 2 By maximizing the coefficients of M j and M j 1 we find that they do not exceed h 2 j/3 5/2. Let A i j be the elements of B 1 where B is the coefficient matrix. Then n 1 (f i+1 f i )/h i+1 (f i f i 1 )/h i M j = 6 A ji + A j0 d 0 + A jn d n. h i + h i+1 i=1 11
12 Also let µ(f, δ) be the modulus of continuity of f(x) on [a, b] defined as µ(f, δ) = sup f(x) f(x ) x x <δ Then remembering the condition that δ β, we find that h 2 j (f i+1 f i )/h i+1 (f i f i 1 )/h i h i + h i+1 β2 µ(f, δ) Thus h 2 j ( M j 1 + M j ) 2 B 1 { 6β 2 µ(f, δ) + δ 2 ( d 0 + d n ) } and finally S(x) f(x) B 1 { 6β 2 µ(f, δ) + δ 2 ( d 0 + d n ) } ( + µ(f, δ) ) 2 If the function satisfies f(x) f(x ) κ(x x ) α, 0 < α 1 then µ(f, δ) κδ α establishing the required result, in light of the inequality satisfied by B 1, provided that λ < 2, µ n < 2. If f is in C 2 we can prove Theorem B: (Ahlberg, p. 27) 12
13 f (p) x S (p) (x) = O(δ 2 p ) where S(x) satisfies end conditions 2M 0 + M 1 = 6 ( ) f1 f 0 f 0 h 1 h 1 or M 0 = f 0. Proof: If B is the matrix of the spline we have BM = d whence B ( M d ) ( = I B ) d 3 3 The right side is 1 3 d 0 d 1 µ 1 (d 1 d 0 ) λ 1 (d 2 d 1 ) d N d N 1 Now at interior points d j 6 = (f j+1 f j )/+1 (f j f j 1 )/ + +1 = f[x j 1, x j, x j+1 ] = 1 2 f (ξ j ), x j 1 ξ j x j+1 13
14 Also by the Taylor theorem d 0 6 = (f 1 f 0 )/h 1 f 0 h 1 = 1 2 f (ξ 0 ), x 0 ξ 0 x 1 Thus if µ is the modulus of continuity of f, f (x) f (x ) µ(f, δ), x x δ then here since λ j + µ j = 1, ( I B ) d 3 3µ(f, δ) Also with the given end conditions B 1 1 so M d ( 3 B 1 I B ) d 3 3µ(f, δ) Now since d j 3 = f (ξ j ) it also follows that f d 3 µ(f, δ) whence M f 4µ(f, δ) and since S is linear on each interval S f 5µ(f, δ) 14
15 Now if f (x) f (x ) κ x x α, 0 < α 1 we have S f 5κδ α Since S(x j ) = f(x j ), by Rolle s theorem, there is a point ξ j in the interval [x j 1, x j ] at which S (ξ j ) = f (ξ j ). Then on this interval f (x) S (x) = x ξ j (f (x) S (x)) 5δµ(f, δ) and integrating a second time f(x) S(x) 5 2 δ2 µ(f, δ) Even stronger convergence theorems can be proved. Ahlberg (p. 29) gives: Theorem C: Let f (4) (x) be continuous then f (p) (x) S (p) (x) = O ( δ 3 p), p = 0, 1, 2, 3 B spline 15
16 6 Smoothing splines Figure 3: B spline Instead of requiring a spline to coincide with the function at specidied points one may simply require it to pass closer and permit some smoothing. To do this let S(x) be chosen to minimize R 1 xn x 0 f (x) 2 dx + R 2 n (f(x i ) y i ) 2 Q i=0 i Then by the calculus of variations the solution can be shown to be a spline witoints at x 1, x 2, x n 1 which does not necessarily pass through thte points. Reinsch (Numer. Math 10, 1967, pp ) shows that the equations now reduce to a five diagonal system, which can still be solved with a number of operations proportional to n. 16
17 7 Change of independent variable With any type of polynomial fitting or spline fitting difficulties occur if the function y(x) has too large a slope. In extreme cases the given points may represent a closed curve or other curve that bends back on itself. Figure 4: In such a case the x i do not form a monotone sequence, and also dy is inifinite at dx the nose, while dx is infinite at the top and bottom. To overcome this difficulty we dy need to parameterize the curve and set x = x(t) y = y(t) where t is a new independent variable that increases monotonely. Then form separate splines for x(t) and y(t). 17
18 8 Integration Splines are useful for integration when function is defined by a table of coordinates at unequal intervals. Conclusions Acknowledgements References 18
Cubic 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 informationG-SPLINE INTERPOLATION FOR APPROXIMATING THE SOLUTION OF FRACTIONAL DIFFERENTIAL EQUATIONS USING LINEAR MULTI- STEP METHODS
Journal of Al-Nahrain University Vol.0(), December, 00, pp.8- Science G-SPLINE INTERPOLATION FOR APPROXIMATING THE SOLUTION OF FRACTIONAL DIFFERENTIAL EQUATIONS USING LINEAR MULTI- STEP METHODS Osama H.
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 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 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 informationCubic Splines; Bézier Curves
Cubic Splines; Bézier Curves 1 Cubic Splines piecewise approximation with cubic polynomials conditions on the coefficients of the splines 2 Bézier Curves computer-aided design and manufacturing MCS 471
More 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 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 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 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 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 informationLinear Algebra Practice Problems
Linear Algebra Practice Problems Math 24 Calculus III Summer 25, Session II. Determine whether the given set is a vector space. If not, give at least one axiom that is not satisfied. Unless otherwise stated,
More informationReview for the Final Exam
Math 171 Review for the Final Exam 1 Find the limits (4 points each) (a) lim 4x 2 3; x x (b) lim ( x 2 x x 1 )x ; (c) lim( 1 1 ); x 1 ln x x 1 sin (x 2) (d) lim x 2 x 2 4 Solutions (a) The limit lim 4x
More informationMATH 215/255 Solutions to Additional Practice Problems April dy dt
. For the nonlinear system MATH 5/55 Solutions to Additional Practice Problems April 08 dx dt = x( x y, dy dt = y(.5 y x, x 0, y 0, (a Show that if x(0 > 0 and y(0 = 0, then the solution (x(t, y(t of the
More informationChapter 2 Interpolation
Chapter 2 Interpolation Experiments usually produce a discrete set of data points (x i, f i ) which represent the value of a function f (x) for a finite set of arguments {x 0...x n }. If additional data
More informationIntegration - Past Edexcel Exam Questions
Integration - Past Edexcel Exam Questions 1. (a) Given that y = 5x 2 + 7x + 3, find i. - ii. - (b) ( 1 + 3 ) x 1 x dx. [4] 2. Question 2b - January 2005 2. The gradient of the curve C is given by The point
More informationMaxima and Minima. (a, b) of R if
Maxima and Minima Definition Let R be any region on the xy-plane, a function f (x, y) attains its absolute or global, maximum value M on R at the point (a, b) of R if (i) f (x, y) M for all points (x,
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 informationThe Closed Form Reproducing Polynomial Particle Shape Functions for Meshfree Particle Methods
The Closed Form Reproducing Polynomial Particle Shape Functions for Meshfree Particle Methods by Hae-Soo Oh Department of Mathematics, University of North Carolina at Charlotte, Charlotte, NC 28223 June
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 informationGradient Based Optimization Methods
Gradient Based Optimization Methods Antony Jameson, Department of Aeronautics and Astronautics Stanford University, Stanford, CA 94305-4035 1 Introduction Consider the minimization of a function J(x) where
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 informationThere is a unique function s(x) that has the required properties. It turns out to also satisfy
Numerical Analysis Grinshpan Natural Cubic Spline Let,, n be given nodes (strictly increasing) and let y,, y n be given values (arbitrary) Our goal is to produce a function s() with the following properties:
More 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 informationOn a max norm bound for the least squares spline approximant. Carl de Boor University of Wisconsin-Madison, MRC, Madison, USA. 0.
in Approximation and Function Spaces Z Ciesielski (ed) North Holland (Amsterdam), 1981, pp 163 175 On a max norm bound for the least squares spline approximant Carl de Boor University of Wisconsin-Madison,
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 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 5630: Iterative Methods for Systems of Equations Hung Phan, UMass Lowell March 22, 2018
1 Linear Systems Math 5630: Iterative Methods for Systems of Equations Hung Phan, UMass Lowell March, 018 Consider the system 4x y + z = 7 4x 8y + z = 1 x + y + 5z = 15. We then obtain x = 1 4 (7 + y z)
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 informationFinite difference method for elliptic problems: I
Finite difference method for elliptic problems: I Praveen. C praveen@math.tifrbng.res.in Tata Institute of Fundamental Research Center for Applicable Mathematics Bangalore 560065 http://math.tifrbng.res.in/~praveen
More informationChordal cubic spline interpolation is fourth order accurate
Chordal cubic spline interpolation is fourth order accurate Michael S. Floater Abstract: It is well known that complete cubic spline interpolation of functions with four continuous derivatives is fourth
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 informationLinear Algebra Massoud Malek
CSUEB Linear Algebra Massoud Malek Inner Product and Normed Space In all that follows, the n n identity matrix is denoted by I n, the n n zero matrix by Z n, and the zero vector by θ n An inner product
More informationx x2 2 + x3 3 x4 3. Use the divided-difference method to find a polynomial of least degree that fits the values shown: (b)
Numerical Methods - PROBLEMS. The Taylor series, about the origin, for log( + x) is x x2 2 + x3 3 x4 4 + Find an upper bound on the magnitude of the truncation error on the interval x.5 when log( + x)
More 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 informationOn orthogonal polynomials for certain non-definite linear functionals
On orthogonal polynomials for certain non-definite linear functionals Sven Ehrich a a GSF Research Center, Institute of Biomathematics and Biometry, Ingolstädter Landstr. 1, 85764 Neuherberg, Germany Abstract
More informationMaria Cameron Theoretical foundations. Let. be a partition of the interval [a, b].
Maria Cameron 1 Interpolation by spline functions Spline functions yield smooth interpolation curves that are less likely to exhibit the large oscillations characteristic for high degree polynomials Splines
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 informationCIV-E1060 Engineering Computation and Simulation Examination, December 12, 2017 / Niiranen
CIV-E16 Engineering Computation and Simulation Examination, December 12, 217 / Niiranen This examination consists of 3 problems rated by the standard scale 1...6. Problem 1 Let us consider a long and tall
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 informationA Note on Multiplicity Weight of Nodes of Two Point Taylor Expansion
Applied Mathematical Sciences, Vol, 207, no 6, 307-3032 HIKARI Ltd, wwwm-hikaricom https://doiorg/02988/ams2077302 A Note on Multiplicity Weight of Nodes of Two Point Taylor Expansion Koichiro Shimada
More informationMATH 1014 Tutorial Notes 8
MATH4 Calculus II (8 Spring) Topics covered in tutorial 8:. Numerical integration. Approximation integration What you need to know: Midpoint rule & its error Trapezoid rule & its error Simpson s rule &
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 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 informationLecture 9 Approximations of Laplace s Equation, Finite Element Method. Mathématiques appliquées (MATH0504-1) B. Dewals, C.
Lecture 9 Approximations of Laplace s Equation, Finite Element Method Mathématiques appliquées (MATH54-1) B. Dewals, C. Geuzaine V1.2 23/11/218 1 Learning objectives of this lecture Apply the finite difference
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 informationApplied/Numerical Analysis Qualifying Exam
Applied/Numerical Analysis Qualifying Exam August 9, 212 Cover Sheet Applied Analysis Part Policy on misprints: The qualifying exam committee tries to proofread exams as carefully as possible. Nevertheless,
More informationChapter Two: Numerical Methods for Elliptic PDEs. 1 Finite Difference Methods for Elliptic PDEs
Chapter Two: Numerical Methods for Elliptic PDEs Finite Difference Methods for Elliptic PDEs.. Finite difference scheme. We consider a simple example u := subject to Dirichlet boundary conditions ( ) u
More informationGeometric Stiffness Effects in 2D and 3D Frames
Geometric Stiffness Effects in D and 3D Frames CEE 41. Matrix Structural Analsis Department of Civil and Environmental Engineering Duke Universit Henri Gavin Fall, 1 In situations in which deformations
More informationarxiv: v1 [math.na] 1 May 2013
arxiv:3050089v [mathna] May 03 Approximation Properties of a Gradient Recovery Operator Using a Biorthogonal System Bishnu P Lamichhane and Adam McNeilly May, 03 Abstract A gradient recovery operator based
More information11 a 12 a 21 a 11 a 22 a 12 a 21. (C.11) A = The determinant of a product of two matrices is given by AB = A B 1 1 = (C.13) and similarly.
C PROPERTIES OF MATRICES 697 to whether the permutation i 1 i 2 i N is even or odd, respectively Note that I =1 Thus, for a 2 2 matrix, the determinant takes the form A = a 11 a 12 = a a 21 a 11 a 22 a
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 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 informationOR MSc Maths Revision Course
OR MSc Maths Revision Course Tom Byrne School of Mathematics University of Edinburgh t.m.byrne@sms.ed.ac.uk 15 September 2017 General Information Today JCMB Lecture Theatre A, 09:30-12:30 Mathematics revision
More informationMath 4377/6308 Advanced Linear Algebra I Dr. Vaughn Climenhaga, PGH 651A HOMEWORK 3
Math 4377/6308 Advanced Linear Algebra I Dr. Vaughn Climenhaga, PGH 651A Fall 2013 HOMEWORK 3 Due 4pm Wednesday, September 11. You will be graded not only on the correctness of your answers but also on
More informationMath221: HW# 7 solutions
Math22: HW# 7 solutions Andy Royston November 7, 25.3.3 let x = e u. Then ln x = u, x2 = e 2u, and dx = e 2u du. Furthermore, when x =, u, and when x =, u =. Hence x 2 ln x) 3 dx = e 2u u 3 e u du) = e
More informationPolynomial Functions
Polynomial Functions Equations and Graphs Characteristics The Factor Theorem The Remainder Theorem http://www.purplemath.com/modules/polyends5.htm 1 A cross-section of a honeycomb has a pattern with one
More informationREVIEW OF DIFFERENTIAL CALCULUS
REVIEW OF DIFFERENTIAL CALCULUS DONU ARAPURA 1. Limits and continuity To simplify the statements, we will often stick to two variables, but everything holds with any number of variables. Let f(x, y) be
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 informationMultivariable Calculus
2 Multivariable Calculus 2.1 Limits and Continuity Problem 2.1.1 (Fa94) Let the function f : R n R n satisfy the following two conditions: (i) f (K ) is compact whenever K is a compact subset of R n. (ii)
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 informationOptimization using Calculus. Optimization of Functions of Multiple Variables subject to Equality Constraints
Optimization using Calculus Optimization of Functions of Multiple Variables subject to Equality Constraints 1 Objectives Optimization of functions of multiple variables subjected to equality constraints
More informationMATH 205C: STATIONARY PHASE LEMMA
MATH 205C: STATIONARY PHASE LEMMA For ω, consider an integral of the form I(ω) = e iωf(x) u(x) dx, where u Cc (R n ) complex valued, with support in a compact set K, and f C (R n ) real valued. Thus, I(ω)
More informationMecE 390 Final examination, Winter 2014
MecE 390 Final examination, Winter 2014 Directions: (i) a double-sided 8.5 11 formula sheet is permitted, (ii) no calculators are permitted, (iii) the exam is 80 minutes in duration; please turn your paper
More informationMath 0230 Calculus 2 Lectures
Math 00 Calculus Lectures Chapter 8 Series Numeration of sections corresponds to the text James Stewart, Essential Calculus, Early Transcendentals, Second edition. Section 8. Sequences A sequence is a
More informationSECTION A. f(x) = ln(x). Sketch the graph of y = f(x), indicating the coordinates of any points where the graph crosses the axes.
SECTION A 1. State the maximal domain and range of the function f(x) = ln(x). Sketch the graph of y = f(x), indicating the coordinates of any points where the graph crosses the axes. 2. By evaluating f(0),
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 informationLecture 4: Piecewise Cubic Interpolation
Lecture notes on Vritionl nd Approximte Methods in Applied Mthemtics - A Peirce UBC Lecture 4: Piecewise Cubic Interpoltion Compiled 5 September In this lecture we consider piecewise cubic interpoltion
More informationTime Series and Forecasting Lecture 4 NonLinear Time Series
Time Series and Forecasting Lecture 4 NonLinear Time Series Bruce E. Hansen Summer School in Economics and Econometrics University of Crete July 23-27, 2012 Bruce Hansen (University of Wisconsin) Foundations
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 informationReview: Power series define functions. Functions define power series. Taylor series of a function. Taylor polynomials of a function.
Taylor Series (Sect. 10.8) Review: Power series define functions. Functions define power series. Taylor series of a function. Taylor polynomials of a function. Review: Power series define functions Remarks:
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 informationMATH3283W LECTURE NOTES: WEEK 6 = 5 13, = 2 5, 1 13
MATH383W LECTURE NOTES: WEEK 6 //00 Recursive sequences (cont.) Examples: () a =, a n+ = 3 a n. The first few terms are,,, 5 = 5, 3 5 = 5 3, Since 5
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 information3.2 Bases of C[0, 1] and L p [0, 1]
70 CHAPTER 3. BASES IN BANACH SPACES 3.2 Bases of C[0, ] and L p [0, ] In the previous section we introduced the unit vector bases of `p and c 0. Less obvious is it to find bases of function spaces like
More information3.4 Introduction to power series
3.4 Introduction to power series Definition 3.4.. A polynomial in the variable x is an expression of the form n a i x i = a 0 + a x + a 2 x 2 + + a n x n + a n x n i=0 or a n x n + a n x n + + a 2 x 2
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 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 informationMath 116 Final Exam December 15, 2011
Math 6 Final Exam December 5, 2 Name: EXAM SOLUTIONS Instructor: Section:. Do not open this exam until you are told to do so. 2. This exam has 4 pages including this cover. There are problems. Note that
More informationOn the Stability of the Best Reply Map for Noncooperative Differential Games
On the Stability of the Best Reply Map for Noncooperative Differential Games Alberto Bressan and Zipeng Wang Department of Mathematics, Penn State University, University Park, PA, 68, USA DPMMS, University
More informationFinding discontinuities of piecewise-smooth functions
Finding discontinuities of piecewise-smooth functions A.G. Ramm Mathematics Department, Kansas State University, Manhattan, KS 66506-2602, USA ramm@math.ksu.edu Abstract Formulas for stable differentiation
More informationThe Infinity Norm of a Certain Type of Symmetric Circulant Matrix
MATHEMATICS OF COMPUTATION, VOLUME 31, NUMBER 139 JULY 1977, PAGES 733-737 The Infinity Norm of a Certain Type of Symmetric Circulant Matrix By W. D. Hoskins and D. S. Meek Abstract. An attainable bound
More information1 Discretizing BVP with Finite Element Methods.
1 Discretizing BVP with Finite Element Methods In this section, we will discuss a process for solving boundary value problems numerically, the Finite Element Method (FEM) We note that such method is a
More informationMATH 124B: HOMEWORK 2
MATH 24B: HOMEWORK 2 Suggested due date: August 5th, 26 () Consider the geometric series ( ) n x 2n. (a) Does it converge pointwise in the interval < x
More informationChapter 8 ~ Quadratic Functions and Equations In this chapter you will study... You can use these skills...
Chapter 8 ~ Quadratic Functions and Equations In this chapter you will study... identifying and graphing quadratic functions transforming quadratic equations solving quadratic equations using factoring
More information8.7 Taylor s Inequality Math 2300 Section 005 Calculus II. f(x) = ln(1 + x) f(0) = 0
8.7 Taylor s Inequality Math 00 Section 005 Calculus II Name: ANSWER KEY Taylor s Inequality: If f (n+) is continuous and f (n+) < M between the center a and some point x, then f(x) T n (x) M x a n+ (n
More informationChapter 1 Vector Spaces
Chapter 1 Vector Spaces Per-Olof Persson persson@berkeley.edu Department of Mathematics University of California, Berkeley Math 110 Linear Algebra Vector Spaces Definition A vector space V over a field
More informationCubic B-spline Collocation Method for Fourth Order Boundary Value Problems. 1 Introduction
ISSN 1749-3889 print, 1749-3897 online International Journal of Nonlinear Science Vol.142012 No.3,pp.336-344 Cubic B-spline Collocation Method for Fourth Order Boundary Value Problems K.N.S. Kasi Viswanadham,
More informationLecture 1 INF-MAT3350/ : Some Tridiagonal Matrix Problems
Lecture 1 INF-MAT3350/4350 2007: Some Tridiagonal Matrix Problems Tom Lyche University of Oslo Norway Lecture 1 INF-MAT3350/4350 2007: Some Tridiagonal Matrix Problems p.1/33 Plan for the day 1. Notation
More informationCurve Fitting. 1 Interpolation. 2 Composite Fitting. 1.1 Fitting f(x) 1.2 Hermite interpolation. 2.1 Parabolic and Cubic Splines
Curve Fitting Why do we want to curve fit? In general, we fit data points to produce a smooth representation of the system whose response generated the data points We do this for a variety of reasons 1
More informationPH.D. PRELIMINARY EXAMINATION MATHEMATICS
UNIVERSITY OF CALIFORNIA, BERKELEY SPRING SEMESTER 207 Dept. of Civil and Environmental Engineering Structural Engineering, Mechanics and Materials NAME PH.D. PRELIMINARY EXAMINATION MATHEMATICS Problem
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 informationLECTURE 5: THE METHOD OF STATIONARY PHASE
LECTURE 5: THE METHOD OF STATIONARY PHASE Some notions.. A crash course on Fourier transform For j =,, n, j = x j. D j = i j. For any multi-index α = (α,, α n ) N n. α = α + + α n. α! = α! α n!. x α =
More informationExercises for Multivariable Differential Calculus XM521
This document lists all the exercises for XM521. The Type I (True/False) exercises will be given, and should be answered, online immediately following each lecture. The Type III exercises are to be done
More informationCHAPTER 3 Further properties of splines and B-splines
CHAPTER 3 Further properties of splines and B-splines In Chapter 2 we established some of the most elementary properties of B-splines. In this chapter our focus is on the question What kind of functions
More information2 Two-Point Boundary Value Problems
2 Two-Point Boundary Value Problems Another fundamental equation, in addition to the heat eq. and the wave eq., is Poisson s equation: n j=1 2 u x 2 j The unknown is the function u = u(x 1, x 2,..., x
More informationManprit Kaur and Arun Kumar
CUBIC X-SPLINE INTERPOLATORY FUNCTIONS Manprit Kaur and Arun Kumar manpreet2410@gmail.com, arun04@rediffmail.com Department of Mathematic and Computer Science, R. D. Univerity, Jabalpur, INDIA. Abtract:
More informationInterpolation and extrapolation
Interpolation and extrapolation Alexander Khanov PHYS6260: Experimental Methods is HEP Oklahoma State University October 30, 207 Interpolation/extrapolation vs fitting Formulation of the problem: there
More informationA collocation method for solving some integral equations in distributions
A collocation method for solving some integral equations in distributions Sapto W. Indratno Department of Mathematics Kansas State University, Manhattan, KS 66506-2602, USA sapto@math.ksu.edu A G Ramm
More information