Chapter Lagrangian Interpolation

Similar documents
Lagrangian Interpolation

Chapters 2 Kinematics. Position, Distance, Displacement

Notes on the stability of dynamic systems and the use of Eigen Values.

( ) () we define the interaction representation by the unitary transformation () = ()

PHYS 1443 Section 001 Lecture #4

Solution in semi infinite diffusion couples (error function analysis)

Cubic Bezier Homotopy Function for Solving Exponential Equations

Let s treat the problem of the response of a system to an applied external force. Again,

TSS = SST + SSE An orthogonal partition of the total SS

Linear Response Theory: The connection between QFT and experiments

UNIVERSITAT AUTÒNOMA DE BARCELONA MARCH 2017 EXAMINATION

Motion in Two Dimensions

Chapter Floating Point Representation

2/20/2013. EE 101 Midterm 2 Review

CH.3. COMPATIBILITY EQUATIONS. Continuum Mechanics Course (MMC) - ETSECCPB - UPC

The Finite Element Method for the Analysis of Non-Linear and Dynamic Systems

Mechanics Physics 151

Mechanics Physics 151

Volatility Interpolation

GENERATING CERTAIN QUINTIC IRREDUCIBLE POLYNOMIALS OVER FINITE FIELDS. Youngwoo Ahn and Kitae Kim

FI 3103 Quantum Physics

Ordinary Differential Equations in Neuroscience with Matlab examples. Aim 1- Gain understanding of how to set up and solve ODE s

Lecture 18: The Laplace Transform (See Sections and 14.7 in Boas)

CS434a/541a: Pattern Recognition Prof. Olga Veksler. Lecture 4

FTCS Solution to the Heat Equation

Calculus Chapter 1 Introduction to Calculus

[ ] 2. [ ]3 + (Δx i + Δx i 1 ) / 2. Δx i-1 Δx i Δx i+1. TPG4160 Reservoir Simulation 2018 Lecture note 3. page 1 of 5

Should Exact Index Numbers have Standard Errors? Theory and Application to Asian Growth

First-order piecewise-linear dynamic circuits

Variants of Pegasos. December 11, 2009

M. Y. Adamu Mathematical Sciences Programme, AbubakarTafawaBalewa University, Bauchi, Nigeria

Chapter Direct Method of Interpolation

HEAT CONDUCTION PROBLEM IN A TWO-LAYERED HOLLOW CYLINDER BY USING THE GREEN S FUNCTION METHOD

In the complete model, these slopes are ANALYSIS OF VARIANCE FOR THE COMPLETE TWO-WAY MODEL. (! i+1 -! i ) + [(!") i+1,q - [(!

Robust and Accurate Cancer Classification with Gene Expression Profiling

Chapter 5. Circuit Theorems

J i-1 i. J i i+1. Numerical integration of the diffusion equation (I) Finite difference method. Spatial Discretization. Internal nodes.

On computing differential transform of nonlinear non-autonomous functions and its applications

Graduate Macroeconomics 2 Problem set 5. - Solutions

Displacement, Velocity, and Acceleration. (WHERE and WHEN?)

UNIT 1 ONE-DIMENSIONAL MOTION GRAPHING AND MATHEMATICAL MODELING. Objectives

THEORETICAL AUTOCORRELATIONS. ) if often denoted by γ. Note that

Comb Filters. Comb Filters

Response of MDOF systems

Chapter 6: AC Circuits

V.Abramov - FURTHER ANALYSIS OF CONFIDENCE INTERVALS FOR LARGE CLIENT/SERVER COMPUTER NETWORKS

Part II CONTINUOUS TIME STOCHASTIC PROCESSES

. The geometric multiplicity is dim[ker( λi. number of linearly independent eigenvectors associated with this eigenvalue.

Approximate Analytic Solution of (2+1) - Dimensional Zakharov-Kuznetsov(Zk) Equations Using Homotopy

P R = P 0. The system is shown on the next figure:

Homework 8: Rigid Body Dynamics Due Friday April 21, 2017

Econ107 Applied Econometrics Topic 5: Specification: Choosing Independent Variables (Studenmund, Chapter 6)

[Link to MIT-Lab 6P.1 goes here.] After completing the lab, fill in the following blanks: Numerical. Simulation s Calculations

Lecture 11 SVM cont

Performance Analysis for a Network having Standby Redundant Unit with Waiting in Repair

Lecture 6: Learning for Control (Generalised Linear Regression)

TEST CODE O2234OIO FORM TP UNIT2-Paper0l. t hour 30 minutes. 0l JUNE 2016 (a.m.)

On elements with index of the form 2 a 3 b in a parametric family of biquadratic elds

THE PREDICTION OF COMPETITIVE ENVIRONMENT IN BUSINESS

3. OVERVIEW OF NUMERICAL METHODS

Advanced time-series analysis (University of Lund, Economic History Department)

Existence and Uniqueness Results for Random Impulsive Integro-Differential Equation

Bayes rule for a classification problem INF Discriminant functions for the normal density. Euclidean distance. Mahalanobis distance

How about the more general "linear" scalar functions of scalars (i.e., a 1st degree polynomial of the following form with a constant term )?

Energy Storage Devices

Dynamic Team Decision Theory. EECS 558 Project Shrutivandana Sharma and David Shuman December 10, 2005

Online Supplement for Dynamic Multi-Technology. Production-Inventory Problem with Emissions Trading

WebAssign HW Due 11:59PM Tuesday Clicker Information

. The geometric multiplicity is dim[ker( λi. A )], i.e. the number of linearly independent eigenvectors associated with this eigenvalue.

R th is the Thevenin equivalent at the capacitor terminals.

THE method of moment (MOM) is widely used to extract. Fast Green Function Evaluation for Method of Moment. arxiv: v1 [cs.

(1) Cov(, ) E[( E( ))( E( ))]

The ray paths and travel times for multiple layers can be computed using ray-tracing, as demonstrated in Lab 3.

Chapter 2 Linear dynamic analysis of a structural system

John Geweke a and Gianni Amisano b a Departments of Economics and Statistics, University of Iowa, USA b European Central Bank, Frankfurt, Germany

Density Matrix Description of NMR BCMB/CHEM 8190

Example: MOSFET Amplifier Distortion

NATIONAL UNIVERSITY OF SINGAPORE PC5202 ADVANCED STATISTICAL MECHANICS. (Semester II: AY ) Time Allowed: 2 Hours

ISSN MIT Publications

10. A.C CIRCUITS. Theoretically current grows to maximum value after infinite time. But practically it grows to maximum after 5τ. Decay of current :

Density Matrix Description of NMR BCMB/CHEM 8190

Foundations of State Estimation Part II

On One Analytic Method of. Constructing Program Controls

Fuzzy Set Theory in Modeling Uncertainty Data. via Interpolation Rational Bezier Surface Function

Set point control in the state space setting

PHYS 705: Classical Mechanics. Canonical Transformation

Born Oppenheimer Approximation and Beyond

Lecture VI Regression

Research Article Numerical Approximation of Higher-Order Solutions of the Quadratic Nonlinear Stochastic Oscillatory Equation Using WHEP Technique

OMXS30 Balance 20% Index Rules

Including the ordinary differential of distance with time as velocity makes a system of ordinary differential equations.

COMPUTER SCIENCE 349A SAMPLE EXAM QUESTIONS WITH SOLUTIONS PARTS 1, 2

Comparison of Differences between Power Means 1

Physical Simulation Using FEM, Modal Analysis and the Dynamic Equilibrium Equation

(,,, ) (,,, ). In addition, there are three other consumers, -2, -1, and 0. Consumer -2 has the utility function

CHAPTER 10: LINEAR DISCRIMINATION

This document is downloaded from DR-NTU, Nanyang Technological University Library, Singapore.

Normal Random Variable and its discriminant functions

Math 128b Project. Jude Yuen

II. Light is a Ray (Geometrical Optics)

Transcription:

Chaper 5.4 agrangan Inerpolaon Afer readng hs chaper you should be able o:. dere agrangan mehod of nerpolaon. sole problems usng agrangan mehod of nerpolaon and. use agrangan nerpolans o fnd deraes and negrals of dscree funcons. Wha s nerpolaon? Many mes daa s gen only a dscree pons such as x y x y... x n yn x y n n. So how hen does one fnd he alue of y a any oher alue of x? Well a connuous funcon f x may be used o represen he n daa alues wh f x passng hrough he n pons Fgure. Then one can fnd he alue of y a any oher alue of x. Ths s called nerpolaon. Of course f x falls ousde he range of x for whch he daa s gen s no longer nerpolaon bu nsead s called exrapolaon. So wha knd of funcon f x should one choose? A polynomal s a common choce for an nerpolang funcon because polynomals are easy o A ealuae B dfferenae and C negrae relae o oher choces such as a rgonomerc and exponenal seres. Polynomal nerpolaon noles fndng a polynomal of order n ha passes hrough he n daa pons. One of he mehods used o fnd hs polynomal s called he agrangan mehod of nerpolaon. Oher mehods nclude Newon s dded dfference polynomal mehod and he drec mehod. We dscuss he agrangan mehod n hs chaper. 5.5.

5.5. Chaper 5.5 y x y x y x y Fgure Inerpolaon of dscree daa. x y f x x The agrangan nerpolang polynomal s gen by f x n n x f x where n n f n x sands for he n h order polynomal ha approxmaes he funcon y f x gen a x y x y... x y x y and x n n n n daa pons as x x x x x s a weghng funcon ha ncludes a produc of n erms wh erms of omed. The applcaon of agrangan nerpolaon wll be clarfed usng an example. Example The upward elocy of a rocke s gen as a funcon of me n Table. Table Velocy as a funcon of me. s m/s 7.4 5 6.78 57.5.5 6.97 9.67 n n

agrangan Inerpolaon 5.5. Fgure Graph of elocy s. me daa for he rocke example. Deermne he alue of he elocy a 6 seconds usng a frs order agrange polynomal. Soluon For frs order polynomal nerpolaon also called lnear nerpolaon he elocy s gen by

5.5.4 Chaper 5.5 y x y f x x y x Fgure near nerpolaon. Snce we wan o fnd he elocy a 6 and we are usng a frs order polynomal we need o choose he wo daa pons ha are closes o 6 ha also bracke 6 o ealuae. The wo pons are 5 and. Then 5 6. 78 57. 5 ges Hence 5 6.78 57.5 5 5 5 6 6 5 6 6.78 57.5 5 5

agrangan Inerpolaon 5.5.5.86.78.57.5 9.69 m/s You can see ha. 8 and. are lke weghages gen o he eloces a 5 and o calculae he elocy a 6. Quadrac Inerpolaon y x y x y f x x y x Fgure 4 Quadrac nerpolaon. Example The upward elocy of a rocke s gen as a funcon of me n Table. Table Velocy as a funcon of me. s m/s 7.4 5 6.78 57.5.5 6.97 9.67 a Deermne he alue of he elocy a 6 seconds wh second order polynomal nerpolaon usng agrangan polynomal nerpolaon. b Fnd he absolue relae approxmae error for he second order polynomal approxmaon.

5.5.6 Chaper 5.5 Soluon a For second order polynomal nerpolaon also called quadrac nerpolaon he elocy s gen by Snce we wan o fnd he elocy a 6 and we are usng a second order polynomal we need o choose he hree daa pons ha are closes o 6 ha also bracke 6 o ealuae. The hree pons are and 5. Then 4 7. 78 6. 5 5 57. ges Hence 57.5 5 5 6 6 6.78 5 5 6 6 7.4 5 56 6 6.57.5.966.78.87.4 9.9 m/s

agrangan Inerpolaon 5.5.7 b The absolue relae approxmae error a for he second order polynomal s calculaed by consderng he resul of he frs order polynomal Example as he preous approxmaon. 9.9 9.69 a 9.9.84% Example The upward elocy of a rocke s gen as a funcon of me n Table. Table Velocy as a funcon of me s m/s 7.4 5 6.78 57.5.5 6.97 9.67 a Deermne he alue of he elocy a 6 seconds usng hrd order agrangan polynomal nerpolaon. b Fnd he absolue relae approxmae error for he hrd order polynomal approxmaon. c Usng he hrd order polynomal nerpolan for elocy fnd he dsance coered by he rocke from s o 6 s. d Usng he hrd order polynomal nerpolan for elocy fnd he acceleraon of he rocke a 6 s. Soluon a For hrd order polynomal nerpolaon also called cubc nerpolaon he elocy s gen by

5.5.8 Chaper 5.5 Fgure 5 Cubc nerpolaon. Snce we wan o fnd he elocy a 6 and we are usng a hrd order polynomal we need o choose he four daa pons closes o 6 ha also bracke 6 o ealuae. The four pons are 5 and 5.. Then 4 7. 78 6. 5 5 57. 97 6..5 ges x y x y x y x y x f x y

agrangan Inerpolaon 5.5.9 Hence 6.97 5.5.5.5 56 6 6 57.5.5 5.5 56 6 6 6.78.5 5 5 5.5 6 6 6 7.4.5 5.5 6 56 6 6.46.97.57.5.86.78.467.4 9.6 m/s b The absolue percenage relae approxmae error a for he alue obaned for 6 can be obaned by comparng he resul wh ha obaned usng he second order polynomal Example 9.6 9.9 9.6 a.69% c The dsance coered by he rocke beween s o s 6 can be calculaed from he nerpolang polynomal as.5 6.97 5.5.5.5 5 57.5.5 5.5 5 6.78.5 5 5 5.5 7.4.5 5.5 5 6.97.57.5.5 5 5 57.5.5 5.5 5 5 6.78 7.5 5 5.5 7.4.5 5.5 5

5.5. Chaper 5.5 57.5 87.5 675.66 47.5 5.5 7.5 75 4.88 875 45.948 45 65.577 4.45.65.95.544. 5 Noe ha he polynomal s ald beween and. 5 and hence ncludes he lms of and 6. So 6 s 6 s d 6 4.45.65.95.544 d 4 6 4.45.65.95.544 4 65 m d The acceleraon a 6 s gen by 6 d a d 6 Gen ha 4.45.65.95.544. 5 d a d d 4.45.65.95.544 d.65.69.6. 5 a 6.65.696.66 9.665 m/s Noe: There s no need o ge he smplfed hrd order polynomal expresson o conduc he dfferenaon. An expresson of he form ges he derae whou expanson as d d

agrangan Inerpolaon 5.5. INTERPOATION Topc agrange Inerpolaon Summary Texbook noes on he agrangan mehod of nerpolaon Maor General Engneerng Auhors Auar Kaw Mchael Keelas as Resed December 9 Web Se hp://numercalmehods.eng.usf.edu