PowerPoints organized by Dr. Michael R. Gustafson II, Duke University Revised by Prof. Jang, CAU

Size: px
Start display at page:

Download "PowerPoints organized by Dr. Michael R. Gustafson II, Duke University Revised by Prof. Jang, CAU"

Transcription

1 Part 4 Capter 6 Sple ad Peewe Iterpolato PowerPot orgazed y Dr. Mael R. Gutao II Duke Uverty Reved y Pro. Jag CAU All mage opyrgt Te MGraw-Hll Compae I. Permo requred or reproduto or dplay.

2 Capter Ojetve Udertadg d tat t ple mmze ollato y ttg lower-order polyomal to data a peewe ao. Kowg ow to develop ode to perorm tale lookup. Reogzg wy u polyomal are preerale to quadrat ad ger-order ple. Udertadg d te odto tat t uderle a u t. Udertadg te deree etwee atural lamped ad ot-a-kot ed odto. Kowg ow to t a ple to data wt MATLAB ult- uto. Udertadg d ow multdmeoal l terpolato t mplemeted wt MATLAB.

3 Itroduto to Sple - t order polyomal l a lead to erroeou reult eaue o roud-o error ad ollato. A alteratve approa to apply lower- order polyomal a peewe ao to uet o data pot. Tee oetg polyomal are alled ple uto. Trd order urve are alled u ple. Sple mmze ollato ad redue roud-o error due to ter lower-order ature.

4 Hger Order v. Sple Sple elmate ollato y ug mall uet o pot or ea terval rater ta every pot. T epeally ueul we tere are jump te data: a rd order polyomal 5 t order polyomal l 7 t order polyomal d Lear ple eve t order polyomal geerated y ug par o pot at a tme A mu more aeptale auray

5 Lear Sple Lear Sple For data pot tere are - terval ad ea terval a t ow l t F l l t l t ple uto. For lear ple ea uto merely te tragt le oetg te two pot et ea ed o te terval. Lear ple uto a e ormulated a a a Iterept ad te lope are Kot Lear ple uto amout to Newto rt-order polyomal Lear ple uto amout to Newto rt-order polyomal.

6 Eample 6. Q. Ft te data elow wt rt-order ple ad Evaluate te uto at = Soluto or te eod terval rom 4.5 to 7. te uto

7 Sple Developmet a Frt-order ple d tragt-le equato etwee ea par o pot tat Go troug te pot Prmary dadvatage: tey are ot moot. - > at te pot kot were te two ple meet te lope age aruptly. rt dervatve are ot otuou at te kot Seod-order ple d quadrat equato etwee ea par o pot tat Go troug te pot Cotuou rt dervatve at te teror pot Trd-order d ple d u equato etwee ea par o pot tat Go troug te pot Cotuou rt ad eod dervatve at te teror pot Note tat te reult o u ple terpolato are deret rom te reult o a terpolatg u.

8 Quadrat Sple To eure tat te t dervatve are otuou at te kot a ple o at leat + order mut e ued. Te quadrat ple ould ave otuou rt dervatve at te kot. For data pot tere are - terval ad oequetly - ukow oeet tat eed to e oud. a

9 Quadrat Sple Quadrat Sple. Cotuty odto: Te uto mut pa troug all te t T t d d t pot. Te umer o oeet redued to -. a or l dj l l l a a. Te uto value o adjaet polyomal mut e equal at te kot. or equato o tere are - remag odto. o tere are remag odto.

10 Quadrat Sple. Te rt dervatve at te teror ode mut e otuou. -> - equato -> > = -> oly oe more equato eed e determed. or 4. Aume tat t te eod dervatve zero at te rt pot. -> te rt two pot are oeted wt a tragt le. -> Lear equato A = rgt ad de vetor were ot o ad.

11 Eample 6. / Q. Ft te quadrat ple to te ame data ued E. 6.. Ue te reult to etmate te value at =

12 Eample 6. / Soluto we ave our pot ad = terval. Ater applyg te otuty equato ad te zero eod-dervatve odto -=5 oeet e eed to e olved

13 Eample 6. / Tere are two ortomg rom quadrat ple:. Stragt le oetg te rt two pot.. Te ple or te lat terval eem to wg too g. -> Cu ple do ot et tee prolem.

14 Cu Sple Wle data o a partular ze preet may opto or te order o ple uto u ple are preerred eaue tey provde te mplet repreetato tat et te dered appearae o moote. Lear ple ave dotuou rt dervatve Quadrat ple ave dotuou eod dervatve ad requre ettg te eod dervatve at ome pot to a pre-determed value *ut* Quadrat or ger-order ple ted to et te talte eret ger order polyomal ll-odtog or ollato

15 Cu Sple ot I geeral te t ple uto or a u ple a e wrtte a: a d For data pot tere are - terval ad tu 4- ukow to evaluate to olve all te ple uto oeet. Cotuty odto eed to e ated ad rt dervatve ad eod dervatve eed to e otuou at teror pot. -> greatly eae te moote.

16 Solvg Sple Coeet g p Oe odto requre tat te ple uto goe troug te rt ad lat pot o te terval yeldg - equato o te orm: equato o te orm: or or a Aoter requre tat te rt dervatve otuou at d a Aoter requre tat te rt dervatve otuou at ea teror pot yeldg - equato o te orm: or d A trd requre tat te eod dervatve otuou at ea teror pot yeldg - equato o te orm: or d ea teror pot yeldg equato o te orm: 6 or d Tee gve 4-6 total equato ad two addtoal equato are eeded!

17 Solvg Sple Coeet Solvg Sple Coeet Ug te two equato Ug te two equato d Te mpled equato are gve y -> - multaeou trdagoal equato or g q wt - ukow oeet. I we ave two addtoal odto we a ave two addtoal odto we a olve or te.

18 Two Addtoal Equato Tere are everal opto or te al two equato: Natural ed odto - aume te eod ddervatve at tte ed kot are zero. Clamped ed odto - aume te rt dervatve at te rt ad lat kot are kow. Not-a-kot ed odto - ore otuty o te trd dervatve at te eod ad te et-to-lat to kot. reult te rt two terval avg te ame ple uto ad te lat two terval avg te ame ple uto

19 Commoly Ued Ed Codto Codto Frt ad Lat Equato Natural ] [ ] [ Clamped were ' ad ' are te peed rt dervatve at te rt ad lat ode repetvely ] [ p y Not-a-kot

20 Natural Cu Sple Natural Cu Sple Te eod dervatve at te rt ode equal to zero a 6 d 6 d d q Te eod dervatve at te lat ode equal to zero a 6 d d Te al equato a ow e wrtte matr orm a ] [ ] [ ] [ ] [

21 Eample 6. / p / Q. Ft u ple to te ame data ued E.5.. Ue te reult to etmate te value at = 5. Sol. ] [ ] [ ] [ ] [

22 d d Eample 6. /.5 d

23 Peewe Iterpolato MATLAB opto MATLAB a everal ult- uto to mplemet peewe terpolato. Te rt ple: yy=ple y T perorm u ple terpolato geerally ug ot-a-kot odto. I y ota two more value ta a etre te te rt ad lat value y are ued a te dervatve at te ed pot.e. lamped

24 Not-a-kot Sple Eample opto Geerate data: = lpae- 9; y =./+5*.^; Calulate model pot ad determe ot-a-kot terpolato = lpae- ; yy = ple y ; Calulate atual uto value at model pot ad data pot te 9-pot ot-a-kot terpolato old ad te atual uto daed d yr =./+5*.^ plot y o yy - yr --

25 Clamped Sple Eample opto Geerate data w/ rt dervatve ormato: = lpae- 9; y =./+5*. ^; y = [ y -4] Calulate model pot ad determe ot-a-kot terpolato = lpae- ; yy = ple y ; Calulate atual uto value at model pot ad data pot te 9-pot lamped terpolato old ad te atual uto daed yr =./+5*.^ plot y o yy - yr --

26 MATLAB terp Futo opto Wle ple a oly perorm u ple MATLAB terp uto a perorm everal deret kd o terpolato: y = terp y metod & y ota te orgal data ota te pot at w to terpolate metod a trg otag te dered metod: earet - earet egor terpolato lear - oet te pot wt tragt le ple - ot-a-kot u ple terpolato pp or u - peewe u Hermte terpolato

27 Peewe Polyomal Comparo opto

Chapter 5. Curve fitting

Chapter 5. Curve fitting Chapter 5 Curve ttg Assgmet please use ecell Gve the data elow use least squares regresso to t a a straght le a power equato c a saturato-growthrate equato ad d a paraola. Fd the r value ad justy whch

More information

PGE 310: Formulation and Solution in Geosystems Engineering. Dr. Balhoff. Interpolation

PGE 310: Formulation and Solution in Geosystems Engineering. Dr. Balhoff. Interpolation PGE 30: Formulato ad Soluto Geosystems Egeerg Dr. Balhoff Iterpolato Numercal Methods wth MATLAB, Recktewald, Chapter 0 ad Numercal Methods for Egeers, Chapra ad Caale, 5 th Ed., Part Fve, Chapter 8 ad

More information

Given a table of data poins of an unknown or complicated function f : we want to find a (simpler) function p s.t. px (

Given a table of data poins of an unknown or complicated function f : we want to find a (simpler) function p s.t. px ( Iterpolato 1 Iterpolato Gve a table of data pos of a ukow or complcated fucto f : y 0 1 2 y y y y 0 1 2 we wat to fd a (smpler) fucto p s.t. p ( ) = y for = 0... p s sad to terpolate the table or terpolate

More information

Basic Concepts in Numerical Analysis November 6, 2017

Basic Concepts in Numerical Analysis November 6, 2017 Basc Cocepts Nuercal Aalyss Noveber 6, 7 Basc Cocepts Nuercal Aalyss Larry Caretto Mecacal Egeerg 5AB Sear Egeerg Aalyss Noveber 6, 7 Outle Revew last class Mdter Exa Noveber 5 covers ateral o deretal

More information

3/20/2013. Splines There are cases where polynomial interpolation is bad overshoot oscillations. Examplef x. Interpolation at -4,-3,-2,-1,0,1,2,3,4

3/20/2013. Splines There are cases where polynomial interpolation is bad overshoot oscillations. Examplef x. Interpolation at -4,-3,-2,-1,0,1,2,3,4 // Sples There re ses where polyoml terpolto s d overshoot oslltos Emple l s Iterpolto t -,-,-,-,,,,,.... - - - Ide ehd sples use lower order polyomls to oet susets o dt pots mke oetos etwee djet sples

More information

Mathematically, integration is just finding the area under a curve from one point to another. It is b

Mathematically, integration is just finding the area under a curve from one point to another. It is b Numerl Metods or Eg [ENGR 9] [Lyes KADEM 7] CHAPTER VI Numerl Itegrto Tops - Rem sums - Trpezodl rule - Smpso s rule - Rrdso s etrpolto - Guss qudrture rule Mtemtlly, tegrto s just dg te re uder urve rom

More information

Data Processing Techniques

Data Processing Techniques Uverstas Gadjah Mada Departmet o Cvl ad Evrometal Egeerg Master o Egeerg Natural Dsaster Maagemet Data Processg Techques Curve Fttg: Regresso ad Iterpolato 3Oct7 Curve Fttg Reerece Chapra, S.C., Caale

More information

A class of efficient Ratio type estimators for the Estimation of Population Mean Using the auxilliary information in survey sampling

A class of efficient Ratio type estimators for the Estimation of Population Mean Using the auxilliary information in survey sampling 0 IJE olume Iue IN: - A la o eet ato tpe etmator or te Etmato o Populato Mea Ug te aullar ormato urve amplg Mr uzar. Maqool ad T. A. aja vo o Agrultural tatt KUAT-Kamr 00 Ida. Atrat: - Etmato o te populato

More information

CS5620 Intro to Computer Graphics

CS5620 Intro to Computer Graphics CS56 Itro to Computer Graphcs Geometrc Modelg art II Geometrc Modelg II hyscal Sples Curve desg pre-computers Cubc Sples Stadard sple put set of pots { } =, No dervatves specfed as put Iterpolate by cubc

More information

Numerical Analysis Formulae Booklet

Numerical Analysis Formulae Booklet Numercal Aalyss Formulae Booklet. Iteratve Scemes for Systems of Lear Algebrac Equatos:.... Taylor Seres... 3. Fte Dfferece Approxmatos... 3 4. Egevalues ad Egevectors of Matrces.... 3 5. Vector ad Matrx

More information

SYSTEMS OF NON-LINEAR EQUATIONS. Introduction Graphical Methods Close Methods Open Methods Polynomial Roots System of Multivariable Equations

SYSTEMS OF NON-LINEAR EQUATIONS. Introduction Graphical Methods Close Methods Open Methods Polynomial Roots System of Multivariable Equations SYSTEMS OF NON-LINEAR EQUATIONS Itoduto Gaphal Method Cloe Method Ope Method Polomal Root Stem o Multvaale Equato Chapte Stem o No-Lea Equato /. Itoduto Polem volvg o-lea equato egeeg lude optmato olvg

More information

Introduction to Numerical Differentiation and Interpolation March 10, !=1 1!=1 2!=2 3!=6 4!=24 5!= 120

Introduction to Numerical Differentiation and Interpolation March 10, !=1 1!=1 2!=2 3!=6 4!=24 5!= 120 Itroducto to Numercal Deretato ad Iterpolato Marc, Itroducto to Numercal Deretato ad Iterpolato Larr Caretto Mecacal Egeerg 9 Numercal Aalss o Egeerg stems Marc, Itroducto Iterpolato s te use o a dscrete

More information

L5 Polynomial / Spline Curves

L5 Polynomial / Spline Curves L5 Polyomal / Sple Curves Cotets Coc sectos Polyomal Curves Hermte Curves Bezer Curves B-Sples No-Uform Ratoal B-Sples (NURBS) Mapulato ad Represetato of Curves Types of Curve Equatos Implct: Descrbe a

More information

Taylor s Series and Interpolation. Interpolation & Curve-fitting. CIS Interpolation. Basic Scenario. Taylor Series interpolates at a specific

Taylor s Series and Interpolation. Interpolation & Curve-fitting. CIS Interpolation. Basic Scenario. Taylor Series interpolates at a specific CIS 54 - Iterpolato Roger Crawfs Basc Scearo We are able to prod some fucto, but do ot kow what t really s. Ths gves us a lst of data pots: [x,f ] f(x) f f + x x + August 2, 25 OSU/CIS 54 3 Taylor s Seres

More information

ENGI 4430 Numerical Integration Page 5-01

ENGI 4430 Numerical Integration Page 5-01 ENGI 443 Numercal Itegrato Page 5-5. Numercal Itegrato I some o our prevous work, (most otaly the evaluato o arc legth), t has ee dcult or mpossle to d the dete tegral. Varous symolc algera ad calculus

More information

Chapter 12-b Integral Calculus - Extra

Chapter 12-b Integral Calculus - Extra C - Itegrl Clulus Cpter - Itegrl Clulus - Etr Is Newto Toms Smpso BONUS Itroduto to Numerl Itegrto C - Itegrl Clulus Numerl Itegrto Ide s to do tegrl smll prts, lke te wy we preseted tegrto: summto. Numerl

More information

INTERPOLATION(2) ELM1222 Numerical Analysis. ELM1222 Numerical Analysis Dr Muharrem Mercimek

INTERPOLATION(2) ELM1222 Numerical Analysis. ELM1222 Numerical Analysis Dr Muharrem Mercimek ELM Numerl Alss Dr Murrem Merme INTEROLATION ELM Numerl Alss Some of te otets re dopted from Luree V. Fusett Appled Numerl Alss usg MATLAB. rete Hll I. 999 ELM Numerl Alss Dr Murrem Merme Tod s leture

More information

Lecture 5: Interpolation. Polynomial interpolation Rational approximation

Lecture 5: Interpolation. Polynomial interpolation Rational approximation Lecture 5: Iterpolato olyomal terpolato Ratoal appromato Coeffcets of the polyomal Iterpolato: Sometme we kow the values of a fucto f for a fte set of pots. Yet we wat to evaluate f for other values perhaps

More information

ESS Line Fitting

ESS Line Fitting ESS 5 014 17. Le Fttg A very commo problem data aalyss s lookg for relatoshpetwee dfferet parameters ad fttg les or surfaces to data. The smplest example s fttg a straght le ad we wll dscuss that here

More information

Chapter Simpson s 1/3 Rule of Integration. ( x)

Chapter Simpson s 1/3 Rule of Integration. ( x) Cpter 7. Smpso s / Rule o Itegrto Ater redg ts pter, you sould e le to. derve te ormul or Smpso s / rule o tegrto,. use Smpso s / rule t to solve tegrls,. develop te ormul or multple-segmet Smpso s / rule

More information

Numerical Analysis Topic 4: Least Squares Curve Fitting

Numerical Analysis Topic 4: Least Squares Curve Fitting Numerl Alss Top 4: Lest Squres Curve Fttg Red Chpter 7 of the tetook Alss_Numerk Motvto Gve set of epermetl dt: 3 5. 5.9 6.3 The reltoshp etwee d m ot e ler. Fd futo f tht est ft the dt 3 Alss_Numerk Motvto

More information

Analysis of Lagrange Interpolation Formula

Analysis of Lagrange Interpolation Formula P IJISET - Iteratoal Joural of Iovatve Scece, Egeerg & Techology, Vol. Issue, December 4. www.jset.com ISS 348 7968 Aalyss of Lagrage Iterpolato Formula Vjay Dahya PDepartmet of MathematcsMaharaja Surajmal

More information

4 Round-Off and Truncation Errors

4 Round-Off and Truncation Errors HK Km Slgtly moded 3//9, /8/6 Frstly wrtte at Marc 5 4 Roud-O ad Trucato Errors Errors Roud-o Errors Trucato Errors Total Numercal Errors Bluders, Model Errors, ad Data Ucertaty Recallg, dv dt Δv v t Δt

More information

Analyzing Control Structures

Analyzing Control Structures Aalyzg Cotrol Strutures sequeg P, P : two fragmets of a algo. t, t : the tme they tae the tme requred to ompute P ;P s t t Θmaxt,t For loops for to m do P t: the tme requred to ompute P total tme requred

More information

Johns Hopkins University Department of Biostatistics Math Review for Introductory Courses

Johns Hopkins University Department of Biostatistics Math Review for Introductory Courses Johs Hopks Uverst Departmet of Bostatstcs Math Revew for Itroductor Courses Ratoale Bostatstcs courses wll rel o some fudametal mathematcal relatoshps, fuctos ad otato. The purpose of ths Math Revew s

More information

General Method for Calculating Chemical Equilibrium Composition

General Method for Calculating Chemical Equilibrium Composition AE 6766/Setzma Sprg 004 Geeral Metod for Calculatg Cemcal Equlbrum Composto For gve tal codtos (e.g., for gve reactats, coose te speces to be cluded te products. As a example, for combusto of ydroge wt

More information

Chapter Newton-Raphson Method of Solving Simultaneous Nonlinear Equations

Chapter Newton-Raphson Method of Solving Simultaneous Nonlinear Equations Chapter 7 Newto-Rapho Method o Solg Smltaeo Nolear Eqato Ater readg th chapter o hold be able to: dere the Newto-Rapho method ormla or mltaeo olear eqato deelop the algorthm o the Newto-Rapho method or

More information

Johns Hopkins University Department of Biostatistics Math Review for Introductory Courses

Johns Hopkins University Department of Biostatistics Math Review for Introductory Courses Johs Hopks Uverst Departmet of Bostatstcs Math Revew for Itroductor Courses Ratoale Bostatstcs courses wll rel o some fudametal mathematcal relatoshps, fuctos ad otato. The purpose of ths Math Revew s

More information

Numerical Differentiation

Numerical Differentiation College o Egeerg ad Computer Scece Mecacal Egeerg Departmet Numercal Aalyss Notes November 4, 7 Istructor: Larry Caretto Numercal Deretato Itroducto Tese otes provde a basc troducto to umercal deretato

More information

Supervised learning: Linear regression Logistic regression

Supervised learning: Linear regression Logistic regression CS 57 Itroducto to AI Lecture 4 Supervsed learg: Lear regresso Logstc regresso Mlos Hauskrecht mlos@cs.ptt.edu 539 Seott Square CS 57 Itro to AI Data: D { D D.. D D Supervsed learg d a set of eamples s

More information

Numerical Differentiation and Integration

Numerical Differentiation and Integration Numerl Deretto d Itegrto Overvew Numerl Deretto Newto-Cotes Itegrto Formuls Trpezodl rule Smpso s Rules Guss Qudrture Cheyshev s ormul Numerl Deretto Forwrd te dvded deree Bkwrd te dvded deree Ceter te

More information

Computer Graphics. Geometric Modeling. Geometric Modeling. Page. Copyright Gotsman, Elber, Barequet, Karni, Sheffer Computer Science - Technion

Computer Graphics. Geometric Modeling. Geometric Modeling. Page. Copyright Gotsman, Elber, Barequet, Karni, Sheffer Computer Science - Technion Computer Graphcs Geometrc Modelg Geometrc Modelg A Example 4 Outle Objectve: Develop methods ad algorthms to mathematcally model shape of real world objects Categores: Wre-Frame Represetato Object s represeted

More information

Numerical Differentiation

Numerical Differentiation College o Egeerg ad Computer Scece Mecacal Egeerg Departmet ME 9 Numercal Aalyss Marc 4, 4 Istructor: Larry Caretto Numercal Deretato Itroducto Tese otes provde a basc troducto to umercal deretato usg

More information

Homework Assignment Number Eight Solutions

Homework Assignment Number Eight Solutions D Keer MSE 0 Dept o Materals Scece & Egeerg Uversty o Teessee Kovlle Homework Assgmet Number Eght Solutos Problem Fd the soluto to the ollowg system o olear algebrac equatos ear () Soluto: s Sce ths s

More information

Computational Methods CMSC/AMSC/MAPL 460. Quadrature: Integration

Computational Methods CMSC/AMSC/MAPL 460. Quadrature: Integration Computatioal Metods CMSC/AMSC/MAPL 6 Quadrature: Itegratio Ramai Duraiswami, Dept. o Computer Siee Some material adapted rom te olie slides o Eri Sadt ad Diae O Leary Numerial Itegratio Idea is to do itegral

More information

A Conventional Approach for the Solution of the Fifth Order Boundary Value Problems Using Sixth Degree Spline Functions

A Conventional Approach for the Solution of the Fifth Order Boundary Value Problems Using Sixth Degree Spline Functions Appled Matheatcs, 1, 4, 8-88 http://d.do.org/1.4/a.1.448 Publshed Ole Aprl 1 (http://www.scrp.org/joural/a) A Covetoal Approach for the Soluto of the Ffth Order Boudary Value Probles Usg Sth Degree Sple

More information

Lecture 12 APPROXIMATION OF FIRST ORDER DERIVATIVES

Lecture 12 APPROXIMATION OF FIRST ORDER DERIVATIVES FDM: Appromato of Frst Order Dervatves Lecture APPROXIMATION OF FIRST ORDER DERIVATIVES. INTRODUCTION Covectve term coservato equatos volve frst order dervatves. The smplest possble approach for dscretzato

More information

Simple Linear Regression

Simple Linear Regression Statstcal Methods I (EST 75) Page 139 Smple Lear Regresso Smple regresso applcatos are used to ft a model descrbg a lear relatoshp betwee two varables. The aspects of least squares regresso ad correlato

More information

Reaction Time VS. Drug Percentage Subject Amount of Drug Times % Reaction Time in Seconds 1 Mary John Carl Sara William 5 4

Reaction Time VS. Drug Percentage Subject Amount of Drug Times % Reaction Time in Seconds 1 Mary John Carl Sara William 5 4 CHAPTER Smple Lear Regreo EXAMPLE A expermet volvg fve ubject coducted to determe the relatohp betwee the percetage of a certa drug the bloodtream ad the legth of tme t take the ubject to react to a tmulu.

More information

MMJ 1113 FINITE ELEMENT METHOD Introduction to PART I

MMJ 1113 FINITE ELEMENT METHOD Introduction to PART I MMJ FINITE EEMENT METHOD Cotut requremets Assume that the fuctos appearg uder the tegral the elemet equatos cota up to (r) th order To esure covergece N must satsf Compatblt requremet the fuctos must have

More information

Lecture 26 Finite Differences and Boundary Value Problems

Lecture 26 Finite Differences and Boundary Value Problems 4//3 Leture 6 Fnte erenes and Boundar Value Problems Numeral derentaton A nte derene s an appromaton o a dervatve - eample erved rom Talor seres 3 O! Negletng all terms ger tan rst order O O Tat s te orward

More information

Chapter 3. Differentiation 3.3 Differentiation Rules

Chapter 3. Differentiation 3.3 Differentiation Rules 3.3 Dfferetato Rules 1 Capter 3. Dfferetato 3.3 Dfferetato Rules Dervatve of a Costat Fucto. If f as te costat value f(x) = c, te f x = [c] = 0. x Proof. From te efto: f (x) f(x + ) f(x) o c c 0 = 0. QED

More information

7.0 Equality Contraints: Lagrange Multipliers

7.0 Equality Contraints: Lagrange Multipliers Systes Optzato 7.0 Equalty Cotrats: Lagrage Multplers Cosder the zato of a o-lear fucto subject to equalty costrats: g f() R ( ) 0 ( ) (7.) where the g ( ) are possbly also olear fuctos, ad < otherwse

More information

CS321. Introduction to Numerical Methods

CS321. Introduction to Numerical Methods CS Itroducto to Numercl Metods Lecture Revew Proessor Ju Zg Deprtmet o Computer Scece Uversty o Ketucky Legto, KY 6 6 Mrc 7, Number Coverso A geerl umber sould be coverted teger prt d rctol prt seprtely

More information

1 0, x? x x. 1 Root finding. 1.1 Introduction. Solve[x^2-1 0,x] {{x -1},{x 1}} Plot[x^2-1,{x,-2,2}] 3

1 0, x? x x. 1 Root finding. 1.1 Introduction. Solve[x^2-1 0,x] {{x -1},{x 1}} Plot[x^2-1,{x,-2,2}] 3 Adrew Powuk - http://www.powuk.com- Math 49 (Numercal Aalyss) Root fdg. Itroducto f ( ),?,? Solve[^-,] {{-},{}} Plot[^-,{,-,}] Cubc equato https://e.wkpeda.org/wk/cubc_fucto Quartc equato https://e.wkpeda.org/wk/quartc_fucto

More information

On the convergence of derivatives of Bernstein approximation

On the convergence of derivatives of Bernstein approximation O the covergece of dervatves of Berste approxmato Mchael S. Floater Abstract: By dfferetatg a remader formula of Stacu, we derve both a error boud ad a asymptotc formula for the dervatves of Berste approxmato.

More information

r 1,2 = a 1 (a 2 2 4a 2 a 0 )

r 1,2 = a 1 (a 2 2 4a 2 a 0 ) Fall emester- 6 Numeral Methods Mehaal Egeerg - ME Dr. aeed J. Almalow, smalow@taahu.edu.sa CHAPTER. Roots of Equatos, Lear ystems ad Algera ystems.. Roots of Polyomal The polyomal equato s: f( ) = a +

More information

Singular Value Decomposition. Linear Algebra (3) Singular Value Decomposition. SVD and Eigenvectors. Solving LEs with SVD

Singular Value Decomposition. Linear Algebra (3) Singular Value Decomposition. SVD and Eigenvectors. Solving LEs with SVD Sgular Value Decomosto Lear Algera (3) m Cootes Ay m x matrx wth m ca e decomosed as follows Dagoal matrx A UWV m x x Orthogoal colums U U I w1 0 0 w W M M 0 0 x Orthoormal (Pure rotato) VV V V L 0 L 0

More information

Ordinary Least Squares Regression. Simple Regression. Algebra and Assumptions.

Ordinary Least Squares Regression. Simple Regression. Algebra and Assumptions. Ordary Least Squares egresso. Smple egresso. Algebra ad Assumptos. I ths part of the course we are gog to study a techque for aalysg the lear relatoshp betwee two varables Y ad X. We have pars of observatos

More information

Line Fitting and Regression

Line Fitting and Regression Marquette Uverst MSCS6 Le Fttg ad Regresso Dael B. Rowe, Ph.D. Professor Departmet of Mathematcs, Statstcs, ad Computer Scece Coprght 8 b Marquette Uverst Least Squares Regresso MSCS6 For LSR we have pots

More information

LECTURE 8: Topics in Chaos Ricker Equation. Period doubling bifurcation. Period doubling cascade. A Quadratic Equation Ricker Equation 1.0. x x 4 0.

LECTURE 8: Topics in Chaos Ricker Equation. Period doubling bifurcation. Period doubling cascade. A Quadratic Equation Ricker Equation 1.0. x x 4 0. LECTURE 8: Topcs Chaos Rcker Equato (t ) = (t ) ep( (t )) Perod doulg urcato Perod doulg cascade 9....... A Quadratc Equato Rcker Equato (t ) = (t ) ( (t ) ). (t ) = (t ) ep( (t )) 6. 9 9. The perod doulg

More information

Multivariate Transformation of Variables and Maximum Likelihood Estimation

Multivariate Transformation of Variables and Maximum Likelihood Estimation Marquette Uversty Multvarate Trasformato of Varables ad Maxmum Lkelhood Estmato Dael B. Rowe, Ph.D. Assocate Professor Departmet of Mathematcs, Statstcs, ad Computer Scece Copyrght 03 by Marquette Uversty

More information

Chapter 3. Differentiation 3.2 Differentiation Rules for Polynomials, Exponentials, Products and Quotients

Chapter 3. Differentiation 3.2 Differentiation Rules for Polynomials, Exponentials, Products and Quotients 3.2 Dfferetato Rules 1 Capter 3. Dfferetato 3.2 Dfferetato Rules for Polyomals, Expoetals, Proucts a Quotets Rule 1. Dervatve of a Costat Fucto. If f as te costat value f(x) = c, te f x = [c] = 0. x Proof.

More information

Linear Approximating to Integer Addition

Linear Approximating to Integer Addition Lear Approxmatg to Iteger Addto L A-Pg Bejg 00085, P.R. Cha apl000@a.com Abtract The teger addto ofte appled cpher a a cryptographc mea. I th paper we wll preet ome reult about the lear approxmatg for

More information

REVIEW OF SIMPLE LINEAR REGRESSION SIMPLE LINEAR REGRESSION

REVIEW OF SIMPLE LINEAR REGRESSION SIMPLE LINEAR REGRESSION REVIEW OF SIMPLE LINEAR REGRESSION SIMPLE LINEAR REGRESSION I lear regreo, we coder the frequecy dtrbuto of oe varable (Y) at each of everal level of a ecod varable (X). Y kow a the depedet varable. The

More information

. The set of these sums. be a partition of [ ab, ]. Consider the sum f( x) f( x 1)

. The set of these sums. be a partition of [ ab, ]. Consider the sum f( x) f( x 1) Chapter 7 Fuctos o Bouded Varato. Subject: Real Aalyss Level: M.Sc. Source: Syed Gul Shah (Charma, Departmet o Mathematcs, US Sargodha Collected & Composed by: Atq ur Rehma (atq@mathcty.org, http://www.mathcty.org

More information

Simple Linear Regression Analysis

Simple Linear Regression Analysis LINEAR REGREION ANALYSIS MODULE II Lecture - 5 Smple Lear Regreo Aaly Dr Shalabh Departmet of Mathematc Stattc Ida Ittute of Techology Kapur Jot cofdece rego for A jot cofdece rego for ca alo be foud Such

More information

Outline. Numerical Heat Transfer. Review All Black Surfaces. Review View Factor, F i j or F ij. Review Gray Diffuse Opaque II

Outline. Numerical Heat Transfer. Review All Black Surfaces. Review View Factor, F i j or F ij. Review Gray Diffuse Opaque II umercao Heat raser ay 9 ad, 7 umercal Heat raser arry Caretto ecacal geerg 75 Heat raser ay 9 ad, 7 Outle Wat s umercal aalyss Cosderatos o coducto, covecto ad radato evew umercal aalyss bascs ervatve

More information

Numerical Differentiation

Numerical Differentiation Part 5 Capter 19 Numercal Derentaton PowerPonts organzed by Dr. Mcael R. Gustason II, Duke Unversty Revsed by Pro. Jang, CAU All mages copyrgt Te McGraw-Hll Companes, Inc. Permsson requred or reproducton

More information

Numerical Interpolation with Polynomials Approximation and Curve Fitting, Focus the MATLAB

Numerical Interpolation with Polynomials Approximation and Curve Fitting, Focus the MATLAB IOSR Joural of Mathematcs (IOSR-JM) e-issn: 78-578, p-issn: 319-765X. Volume 1, Issue 1 Ver. IV (Ja. - Feb. 016), PP 10-0 www.osrjourals.org Numercal Iterpolato wth Polyomals Appromato ad Curve Fttg, Focus

More information

QR Factorization and Singular Value Decomposition COS 323

QR Factorization and Singular Value Decomposition COS 323 QR Factorzato ad Sgular Value Decomposto COS 33 Why Yet Aother Method? How do we solve least-squares wthout currg codto-squarg effect of ormal equatos (A T A A T b) whe A s sgular, fat, or otherwse poorly-specfed?

More information

( ) 2 2. Multi-Layer Refraction Problem Rafael Espericueta, Bakersfield College, November, 2006

( ) 2 2. Multi-Layer Refraction Problem Rafael Espericueta, Bakersfield College, November, 2006 Mult-Layer Refracto Problem Rafael Espercueta, Bakersfeld College, November, 006 Lght travels at dfferet speeds through dfferet meda, but refracts at layer boudares order to traverse the least-tme path.

More information

2/20/2013. Topics. Power Flow Part 1 Text: Power Transmission. Power Transmission. Power Transmission. Power Transmission

2/20/2013. Topics. Power Flow Part 1 Text: Power Transmission. Power Transmission. Power Transmission. Power Transmission /0/0 Topcs Power Flow Part Text: 0-0. Power Trassso Revsted Power Flow Equatos Power Flow Proble Stateet ECEGR 45 Power Systes Power Trassso Power Trassso Recall that for a short trassso le, the power

More information

Chapter 5: Take Home Test

Chapter 5: Take Home Test Chapter : Take Home Test AB Calulus - Hardtke Name Date: Tuesday, / MAY USE YOUR CALCULATOR FOR THIS PAGE. Roud aswers to three plaes. Sore: / Show diagrams ad work to justify eah aswer.. Approimate the

More information

8-4 P 2. = 12 kw. AIR T = const. Therefore, Q &

8-4 P 2. = 12 kw. AIR T = const. Therefore, Q & 8-4 8-4 Air i compreed teadily by a compreor. e air temperature i mataed contant by eat rejection to te urroundg. e rate o entropy cange o air i to be determed. Aumption i i a teady-low proce ce tere i

More information

Problem Set 4 Solutions

Problem Set 4 Solutions 4 Eoom Altos of Gme Theory TA: Youg wg /08/0 - Ato se: A A { B, } S Prolem Set 4 Solutos - Tye Se: T { α }, T { β, β} Se Plyer hs o rte formto, we model ths so tht her tye tke oly oe lue Plyer kows tht

More information

Log1 Contest Round 2 Theta Complex Numbers. 4 points each. 5 points each

Log1 Contest Round 2 Theta Complex Numbers. 4 points each. 5 points each 01 Log1 Cotest Roud Theta Complex Numbers 1 Wrte a b Wrte a b form: 1 5 form: 1 5 4 pots each Wrte a b form: 65 4 4 Evaluate: 65 5 Determe f the followg statemet s always, sometmes, or ever true (you may

More information

Initial-Value Problems for ODEs. numerical errors (round-off and truncation errors) Consider a perturbed system: dz dt

Initial-Value Problems for ODEs. numerical errors (round-off and truncation errors) Consider a perturbed system: dz dt Ital-Value Problems or ODEs d GIVEN: t t,, a FIND: t or atb umercal errors (roud-o ad trucato errors) Cosder a perturbed sstem: dz t, z t, at b z a a Does z(t) (t)? () (uqueess) a uque soluto (t) exsts

More information

UNIT 2 SOLUTION OF ALGEBRAIC AND TRANSCENDENTAL EQUATIONS

UNIT 2 SOLUTION OF ALGEBRAIC AND TRANSCENDENTAL EQUATIONS Numercal Computg -I UNIT SOLUTION OF ALGEBRAIC AND TRANSCENDENTAL EQUATIONS Structure Page Nos..0 Itroducto 6. Objectves 7. Ital Approxmato to a Root 7. Bsecto Method 8.. Error Aalyss 9.4 Regula Fals Method

More information

8 The independence problem

8 The independence problem Noparam Stat 46/55 Jame Kwo 8 The depedece problem 8.. Example (Tua qualty) ## Hollader & Wolfe (973), p. 87f. ## Aemet of tua qualty. We compare the Huter L meaure of ## lghte to the average of coumer

More information

Assignment 7/MATH 247/Winter, 2010 Due: Friday, March 19. Powers of a square matrix

Assignment 7/MATH 247/Winter, 2010 Due: Friday, March 19. Powers of a square matrix Assgmet 7/MATH 47/Wter, 00 Due: Frday, March 9 Powers o a square matrx Gve a square matrx A, ts powers A or large, or eve arbtrary, teger expoets ca be calculated by dagoalzg A -- that s possble (!) Namely,

More information

Current Programmed Control (i.e. Peak Current-Mode Control) Lecture slides part 2 More Accurate Models

Current Programmed Control (i.e. Peak Current-Mode Control) Lecture slides part 2 More Accurate Models Curret Progred Cotrol.e. Pek Curret-Mode Cotrol eture lde prt More Aurte Model ECEN 5807 Drg Mkovć Sple Frt-Order CPM Model: Sury Aupto: CPM otroller operte delly, Ueful reult t low frequee, well uted

More information

Chapter 5 Properties of a Random Sample

Chapter 5 Properties of a Random Sample Lecture 6 o BST 63: Statstcal Theory I Ku Zhag, /0/008 Revew for the prevous lecture Cocepts: t-dstrbuto, F-dstrbuto Theorems: Dstrbutos of sample mea ad sample varace, relatoshp betwee sample mea ad sample

More information

Objectives of Multiple Regression

Objectives of Multiple Regression Obectves of Multple Regresso Establsh the lear equato that best predcts values of a depedet varable Y usg more tha oe eplaator varable from a large set of potetal predctors {,,... k }. Fd that subset of

More information

CH E 374 Computational Methods in Engineering Fall 2007

CH E 374 Computational Methods in Engineering Fall 2007 CH E 7 Computatoal Methods Egeerg Fall 007 Sample Soluto 5. The data o the varato of the rato of stagato pressure to statc pressure (r ) wth Mach umber ( M ) for the flow through a duct are as follows:

More information

r y Simple Linear Regression How To Study Relation Between Two Quantitative Variables? Scatter Plot Pearson s Sample Correlation Correlation

r y Simple Linear Regression How To Study Relation Between Two Quantitative Variables? Scatter Plot Pearson s Sample Correlation Correlation Maatee Klled Correlato & Regreo How To Study Relato Betwee Two Quattatve Varable? Smple Lear Regreo 6.11 A Smple Regreo Problem 1 I there relato betwee umber of power boat the area ad umber of maatee klled?

More information

15-381: Artificial Intelligence. Regression and neural networks (NN)

15-381: Artificial Intelligence. Regression and neural networks (NN) 5-38: Artfcal Itellece Reresso ad eural etorks NN) Mmck the bra I the earl das of AI there as a lot of terest develop models that ca mmc huma thk. Whle o oe ke eactl ho the bra orks ad, eve thouh there

More information

Topic 9 - Taylor and MacLaurin Series

Topic 9 - Taylor and MacLaurin Series Topic 9 - Taylor ad MacLauri Series A. Taylors Theorem. The use o power series is very commo i uctioal aalysis i act may useul ad commoly used uctios ca be writte as a power series ad this remarkable result

More information

Binary classification: Support Vector Machines

Binary classification: Support Vector Machines CS 57 Itroducto to AI Lecture 6 Bar classfcato: Support Vector Maches Mlos Hauskrecht mlos@cs.ptt.edu 539 Seott Square CS 57 Itro to AI Supervsed learg Data: D { D, D,.., D} a set of eamples D, (,,,,,

More information

Optimal MinMax Problems Regarding Capitalization of Compound Interest

Optimal MinMax Problems Regarding Capitalization of Compound Interest Proeedg of the 5th WEA Iteratoal Coferee o Eoomy ad Maagemet Traformato (Volume II Optmal MMax Problem Regardg Captalzato of Compoud Iteret ILIE MITRAN Departmet of Mathemat ILIE RĂCOLEAN, IMOLA DRIGĂ

More information

ENGI 4421 Propagation of Error Page 8-01

ENGI 4421 Propagation of Error Page 8-01 ENGI 441 Propagato of Error Page 8-01 Propagato of Error [Navd Chapter 3; ot Devore] Ay realstc measuremet procedure cotas error. Ay calculatos based o that measuremet wll therefore also cota a error.

More information

CHAPTER 5 Vectors and Vector Space

CHAPTER 5 Vectors and Vector Space HAPTE 5 Vetors d Vetor Spe 5. Alger d eometry of Vetors. Vetor A ordered trple,,, where,, re rel umers. Symol:, B,, A mgtude d dreto.. Norm of vetor,, Norm =,, = = mgtude. Slr multplto Produt of slr d

More information

12.2 Estimating Model parameters Assumptions: ox and y are related according to the simple linear regression model

12.2 Estimating Model parameters Assumptions: ox and y are related according to the simple linear regression model 1. Estmatg Model parameters Assumptos: ox ad y are related accordg to the smple lear regresso model (The lear regresso model s the model that says that x ad y are related a lear fasho, but the observed

More information

Mu Sequences/Series Solutions National Convention 2014

Mu Sequences/Series Solutions National Convention 2014 Mu Sequeces/Seres Solutos Natoal Coveto 04 C 6 E A 6C A 6 B B 7 A D 7 D C 7 A B 8 A B 8 A C 8 E 4 B 9 B 4 E 9 B 4 C 9 E C 0 A A 0 D B 0 C C Usg basc propertes of arthmetc sequeces, we fd a ad bm m We eed

More information

Lecture 07: Poles and Zeros

Lecture 07: Poles and Zeros Lecture 07: Poles ad Zeros Defto of poles ad zeros The trasfer fucto provdes a bass for determg mportat system respose characterstcs wthout solvg the complete dfferetal equato. As defed, the trasfer fucto

More information

Block-Based Compact Thermal Modeling of Semiconductor Integrated Circuits

Block-Based Compact Thermal Modeling of Semiconductor Integrated Circuits Block-Based Compact hermal Modelg of Semcoductor Itegrated Crcuts Master s hess Defese Caddate: Jg Ba Commttee Members: Dr. Mg-Cheg Cheg Dr. Daqg Hou Dr. Robert Schllg July 27, 2009 Outle Itroducto Backgroud

More information

Model Fitting, RANSAC. Jana Kosecka

Model Fitting, RANSAC. Jana Kosecka Model Fttg, RANSAC Jaa Kosecka Fttg: Issues Prevous strateges Le detecto Hough trasform Smple parametrc model, two parameters m, b m + b Votg strateg Hard to geeralze to hgher dmesos a o + a + a 2 2 +

More information

y y ˆ i i Difference between

y y ˆ i i Difference between Fttg Equatos Idea: Te varale of terest (depedet varale, y ) s ard to measure. Tere are easy to measure varales (predctor/ depedet) tat are related to te varale of terest, laeled,,... m We measure te y

More information

Chapter Gauss-Seidel Method

Chapter Gauss-Seidel Method Chpter 04.08 Guss-Sedel Method After redg ths hpter, you should be ble to:. solve set of equtos usg the Guss-Sedel method,. reogze the dvtges d ptflls of the Guss-Sedel method, d. determe uder wht odtos

More information

Section 3. Measurement Errors

Section 3. Measurement Errors eto 3 Measuremet Errors Egeerg Measuremets 3 Types of Errors Itrs errors develops durg the data aqusto proess. Extrs errors foud durg data trasfer ad storage ad are due to the orrupto of the sgal y ose.

More information

1 Lyapunov Stability Theory

1 Lyapunov Stability Theory Lyapuov Stablty heory I ths secto we cosder proofs of stablty of equlbra of autoomous systems. hs s stadard theory for olear systems, ad oe of the most mportat tools the aalyss of olear systems. It may

More information

NUMERICAL DIFFERENTIAL 1

NUMERICAL DIFFERENTIAL 1 NUMERICAL DIFFERENTIAL Ruge-Kutta Metods Ruge-Kutta metods are ver popular ecause o teir good eiciec; ad are used i most computer programs or dieretial equatios. Te are sigle-step metods as te Euler metods.

More information

Multiple Regression. More than 2 variables! Grade on Final. Multiple Regression 11/21/2012. Exam 2 Grades. Exam 2 Re-grades

Multiple Regression. More than 2 variables! Grade on Final. Multiple Regression 11/21/2012. Exam 2 Grades. Exam 2 Re-grades STAT 101 Dr. Kar Lock Morga 11/20/12 Exam 2 Grades Multple Regresso SECTIONS 9.2, 10.1, 10.2 Multple explaatory varables (10.1) Parttog varablty R 2, ANOVA (9.2) Codtos resdual plot (10.2) Trasformatos

More information

Comparison of Four Methods for Estimating. the Weibull Distribution Parameters

Comparison of Four Methods for Estimating. the Weibull Distribution Parameters Appled Mathematal Sees, Vol. 8, 14, o. 83, 4137-4149 HIKARI Ltd, www.m-hkar.om http://dx.do.org/1.1988/ams.14.45389 Comparso of Four Methods for Estmatg the Webull Dstrbuto Parameters Ivaa Pobočíková ad

More information

Phys 2310 Fri. Oct. 23, 2017 Today s Topics. Begin Chapter 6: More on Geometric Optics Reading for Next Time

Phys 2310 Fri. Oct. 23, 2017 Today s Topics. Begin Chapter 6: More on Geometric Optics Reading for Next Time Py F. Oct., 7 Today Topc Beg Capte 6: Moe o Geometc Optc eadg fo Next Tme Homewok t Week HW # Homewok t week due Mo., Oct. : Capte 4: #47, 57, 59, 6, 6, 6, 6, 67, 7 Supplemetal: Tck ee ad e Sytem Pcple

More information

Maximum Likelihood Estimation

Maximum Likelihood Estimation Marquette Uverst Maxmum Lkelhood Estmato Dael B. Rowe, Ph.D. Professor Departmet of Mathematcs, Statstcs, ad Computer Scece Coprght 08 b Marquette Uverst Maxmum Lkelhood Estmato We have bee sag that ~

More information

Applied Fitting Theory VII. Building Virtual Particles

Applied Fitting Theory VII. Building Virtual Particles Appled Fttg heory II Paul Avery CBX 98 38 Jue 8, 998 Apr. 7, 999 (rev.) Buldg rtual Partcles I Statemet of the problem I may physcs aalyses we ecouter the problem of mergg a set of partcles to a sgle partcle

More information

Signal Recovery - Prof. S. Cova - Exam 2016/02/16 - P1 pag.1

Signal Recovery - Prof. S. Cova - Exam 2016/02/16 - P1 pag.1 gal Recovery - Pro.. Cova - Exam 06/0/6 - P pag. PROBEM Data ad Note Appled orce F rt cae: tep ple ecod cae: rectaglar ple wth drato p = 5m Pezoelectrc orce eor A q =0pC/N orce-to-charge covero C = 500pF

More information

MONOPOLISTIC COMPETITION MODEL

MONOPOLISTIC COMPETITION MODEL MONOPOLISTIC COMPETITION MODEL Key gredets Cosumer utlty: log (/ ) log (taste for varety of dfferetated goods) Produto of dfferetated produts: y (/ b) max[ f, ] (reasg returs/fxed osts) Assume that good,

More information

Interval Valued Fuzzy Neutrosophic Soft Structure Spaces

Interval Valued Fuzzy Neutrosophic Soft Structure Spaces Neutrosop Sets ad Systems Vol 5 0 6 Iterval Valued Fuzzy Neutrosop Sot Struture Spaes Iroara & IRSumat Nrmala College or Wome Combatore- 608 Tamladu Ida E-mal: sumat_rama005@yaooo bstrat I ts paper we

More information