Lecture 4: Piecewise Cubic Interpolation
|
|
- Abigayle Daniels
- 5 years ago
- Views:
Transcription
1 Lecture notes on Vrtonl nd Approxmte Methods n Appled Mthemtcs - A Perce UBC Lecture 4: Pecewse Cubc Interpolton Compled 6 August 7 In ths lecture we consder pecewse cubc nterpolton n whch cubc polynoml pproxmton s ssumed over ech subntervl In order to obtn suffcent nformton to determne these coeffcents, we requre contnuty of the nterpoltng polynomls n neghborng ntervls s well s the contnuty of number of dervtves Dependng on the number of dervtves tht we requre to be contnuous t the sttchng ponts we obtn dfferent pproxmton schemes n whch ddtonl dt bout the functon fx beng pproxmted mght hve to be furnshed n order tht there s suffcent dt to determne the coeffcents of the cubc polynomls Key Concepts: Pecewse cubc nterpolton, Cubc Splnes, Cubc Hermte Interpolton 4 Pecewse Cubc Interpolton 4 Degree of freedom nlyss of pecewse cubc nterpolnts Consder the domn [, b] tht s prttoned nto N ntervls hvng N + nodes nd N nternl nodes In ech of these subntervls ssume tht dfferent cubc polynoml s to be constructed Let p 3,N x be the combnton of ll these cubc polynomls We now dscuss the lterntve pproxmton schemes tht one obtns dependng on the dfferent schemes of constrnt Schemes of constrnt: Pecewse Cubc Hermte polynomls: p 3,N x nd p 3,N x re contnuous t nteror nodes There re four unknown coeffcents for ech of the N ntervls nd N constrnts due to the contnuty requrements on p 3,N x nd p 3,N x DOF : 4N N = N + specfy functon vlue nd ts dervtve t ll N + nodes Hve to specfy f & f t ll N + nodes! Cubc Splne: p 3,N x, p 3,N x, nd p 3,N x re ll contnuous t nteror ponts DOF : 4N 3N = N + 3 = N + }{{} + specfy f t ll N + nodes nd mpose EXTRA CONDITIONS The extr condtons re up to the user to prescrbe dependng on the pplcton, eg p 3x = = p 3x N whcs clled the nturl splne
2 4 Pecewse Cubc Hermte polynomls x = x x N = L DOF: 4N N = N + prescrbe fx nd f x t =,, N Representton of f n terms of Hermte bss functons h x nd h x: N N fx hx = fx h x + f x h x where h d x j = δ j nd dx h x j = δ j d dx h x j = h x j = Constructng bss functons on the cnoncl ntervl [, ]: = = } We use the lner Lgrnge bss functons N ξ = +ξ ξ : N ξ b = δ b, ξ,b = ± h h h Let h ξ = α ξ + β [N ξ] ; h ξ = α ξ + β[n ξ] h ξ = γ ξ + δ [N ξ] ; h ξ = γ ξ + δ [N ξ] - h ξ To fnd α, β, γ nd δ we mpose the condtons : = h = β α = h = α + β α = α β α = nd β = h ξ = 4 + ξ ξ Smlrly h ξ = 4 ξ + ξ For the dervtve bss functons = h = δ γ = h = γ δ } Smlrly h ξ = 4 ξ + ξ δ = γ = h ξ = 4 + ξ ξ
3 Expressons for bss functons on n rbtrry ntervl [x, x + ] Interpolton 3 Use the lner trnsformton xξ = x N ξ + x +N ξ = x ξ + x + + ξ x + x + = + x + x ξ The nverse trnsformton s: ξx = x x + x + x + x = x x + x + x + ξ = x + x + x x x + = x x x + x x + x + ξ = x + x x x ξ = x + x ξ = x + x + x x + x x Here x = x + x nd observe tht h x = [ x + x x ] x + x x 3 4 h + x = [ x + x + x] x x x 3 4 h x = x x x + x x 43 h + x = x + xx x x 44 d dξ ξ = d dx ξx dx dξ = d dx ξx x Error nvolved: For functon vlues the error s gven by: whle for dervtves the error s : fx hx f 4 x f x h x f 4 3 x 3 6
4 4 43 Pecewse Cubc Splne nterpolton NDOF: 4N 3N = N + + specfy fx t x,, x N + extr condtons 43 Dervton usng Cubc Herme nterpolton Snce we hve smlr pecewse cubc polynomls to the Pecewse Cubc Hermte polynomls on ech subntervl but wth ddtonl contnuty requred t the N nteror nodes, our strtng pont s the Hermte cubc bss expnson We then mpose ddtonl condtons to mke up for the dervtves f x whch re not known or requred n the cse of splnes On [x k, x k+ ] x k x k x k+ [ x k + x x k ] x k+ x [ x k + x k+ x] x x k sx = f k x k 3 + f k+ x k 3 +s x x k x k+ x k x k + s x x k+ x x k k+ x k where the f k nd f k+ from the Hermte expnson hve been replced by the unknown qunttes s k nd s k+ whch re to be determned by system of equtons whch ensure tht x s contnuous t nternl nodes The frst nd second dervtves of sx must be contnuous t the ntereor ponts x j, j = N Contnuty of the frst dervtve s lredy obtned by our choce of bss functons So we pply contnuty of x to get equtons for the coeffcents s k Tht s, we clculte s x on the two ntervls [x k, x k ] nd [x k, x k+ ] nd requre contnuty t x k After some lgebr x k s k + x k + x k s k + x k s k+ = 3f[x k, x k+ ] x k + f[x k, x k ] x k We hve N equtons nd N + unknowns s, s,, s N so we need more condtons Sy we specfy I f = s nd s N = f N then on unform mesh x k = x: 4 s 4 s 4 4 s N = 3 f[x, x ] + f[x, x ] f[x, x 3 ] + f[x, x ] f[x N, x N ] + f[x N, x N ] s s N A trdgonl system of equtons esy to solve!
5 Interpolton 5 43 Alterntve cubc splne dervton for dervtve prmry vrbles y y + Y x x x + Let Y x = x x 3 + b x x + c x x + d Y x = 3 x x + b x x + c Y x = 6 x x + b y = Y x = d y + = Y x + = h 3 + b h + c + y where = x + x Let s = Y x = c s + = Y x + = 3 h + b + s Y x = b = Y x = 6 + b [ ] [ ] [ h 3 h y+ y 3h = s h ] b s + s / / = y + y s h 4 s + s h h 4 6 y { b = h 3 s + s 3h y + y s } / h 4 y+ y = h 3 + s + + s h y+ y s b = s h s + + s = y + 6 s + s + 6 y s + s 6 y + y = s s + 4 s + s s s + s = 3f[x, x ] + f[x, x + ]
6 6 433 Another splne dervton wth y y y + s prmry vrbles Y x Y x Y + x x x x + Let Y x = x x 3 + b x x + c x x + d Y x = 3 x x + b x x + c Y x = 6 x x + b t x : y = Y x = d nterpoltes the functon t x 45 y + = Y x + = h 3 + b h + c + y enforce contnuty 46 Introduce prmry vrble : = Y x = b b = / 47 Buld n contnuty of Y + = Y x + = 6 + = s Impose contnuty of frst dervtves: c = y + s s + Y x = s x x 3 + s 6 x x + c x x + y s + y + = s h 3 + s 6 + c + y y+ y s + c = s s 6 h = y s + + s 6 s Y + + s 6 Y x = 3 x x + b x x + c = c x = 3 h + b + c s = 3 s 6 h + + { y 6 48 s + s 6 } or = 6 y y = 6 y y =, N EXTRA CONDITIONS N EQ
7 434 Extr Condtons y y y n x h x x h n x N X N Nturl Splne: = = N two extr condtons Interpolton 7 h + h h h h + h h h N h N h N + h N N = 6 y y h y y h y 3 y h y y h y N y N h N y Ny N h N Clmped Splne: Specfed frst dervtves Gven f x = A nd f x n = B A = Y x = c = y y h h s + h 6 = A h + h = 6 y y h A B = Y Nx n = 3 n h n + b n h n + c n sn s n = 3 h n + s n h n + y n y n hn n + h n n 6h n h n 6 6 B y n y n = 3h n n 3h n n + 6h n n h n n h n n h n 6 B y ny n h n = h n n + h n n But h + h + h + h = 6 {[f [x, x ] f [x, x ]} h h h h + h h h n h n + h n h n h n h n N N = 6 y y /h A y y /h y y h y n y n /h n y n y n /h n B ynyn y h n
8 8 3 Qudrtc boundry ntervl representton = s s 6h qudrtc = N = N = n = s N s N 6h N qudrtc 3h + h h h h + h h 4 Lner extrpolton h N3 + h N h N h N h N + 3h N N = 6 f[x, x ] f[x, x ] f[x N, x N f[x N, x N ] S x S S S x x x = h + h h h = s = h h h h + h + h + h = 6f[x, x ] f[x, x ] h h + h h + h h + h + h = RHS h h h h + h h h + h + h = RHS h h +h h h h + h h h h h h + h h h h
9 Interpolton 9 5 Not--knot condton Y x Y x Y N x Y N x Unknowns = 4N Constrnts = 3N 3 4N 3N 3 = N + unknowns the vlues t x,, x N = 6f[x, x ] f[x, x ] Y x = x x 3 + b x x + c x x + d Y x = 3 x x + b x x + c Y x = 6 x x + b, x x = h s = Y x = 6 x x + 6h / = h + h = h + h h h h But h + h + h + h = 6f[x, x ] f[x, x ] + h h + h + h h h Also + h N h N + h N N + h N h N h N h + h = 6f[x, x ] f[x, x ] h N + h N N = 6f[x N, x N ] f[x N, x N ] + h h h + h h h h + h h h + h h 6 A perodc splne: h N h N h N + h N x + = x N x + = x N h N + h N h N h N + h N N N
10 III Smoothness property of nturl splne x = = x N In order to mpose ths condton t s convenent to consder the k s unknowns n whch cse the equtons become: x k k + x k + x k+ k + x k+ k+ = 6f[x k, x k+ ] f[x k, x k ] k =,, N NOTE: Importnt dentty: y [y + ] = y y + = y Let yx be ny other nterpolnt of fx t x,, x N then y dx Now dx = y dx = y s y dx + = y s = y s b b b N y dx x k y s dx = const on ech subntervl k= x k N x k c k y s botnterpolnts k= If we choose s : = = b for exmple or f y nd s botnterpolte f t nd b then y dx = y dx + x k dx dx Thus s s the nterpolnt wth the MINIMUM CURVATURE e smoothest Error nvolved n splne nterpolton: fx = sx f 4 5 x4 384 f x s x f 4 x 3 4 f x s x f 5 x 4 6 y dx = x = mx x : for nonunform ponts nd ll x [x, x N ] for nonunform ponts nd ll x [x, x N ] for unform smple ponts x cubc splne provdes n excellent technque for NUMERICAL DIFFERENTIATION
11 Interpolton How do we solve for the { k }? d c d c 3 d 3 d n c n n d n x x x n = b b n d d x + c x = b 4 d x + d x + c x 3 = b 4 d c x + c x 3 = b b d d d x + c x 3 = b d = d d c b = b d b Smlrly d k+x k+ + c k+ x k+ = b k+ 43 where d k+ = d k+ k+ d k c k b k+ = b k+ k+ d b k 44 k d c d c d n c n d n x x n = b b b n Bck Substtuton x n = b n/d n x n = b n c n x/d n x n = b n/d n 45 x k = b k c k x k+ /d k k = n,, 46
12 Note, f we let m k = k d k then A = m m n d c c n d n LU = LU Forwrd Ly = b y = L b Bck Ux = y Substtuton
Lecture 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 informationDCDM BUSINESS SCHOOL NUMERICAL METHODS (COS 233-8) Solutions to Assignment 3. x f(x)
DCDM BUSINESS SCHOOL NUMEICAL METHODS (COS -8) Solutons to Assgnment Queston Consder the followng dt: 5 f() 8 7 5 () Set up dfference tble through fourth dfferences. (b) Wht s the mnmum degree tht n nterpoltng
More informationINTERPOLATION(1) ELM1222 Numerical Analysis. ELM1222 Numerical Analysis Dr Muharrem Mercimek
ELM Numercl Anlss Dr Muhrrem Mercmek INTEPOLATION ELM Numercl Anlss Some of the contents re dopted from Lurene V. Fusett, Appled Numercl Anlss usng MATLAB. Prentce Hll Inc., 999 ELM Numercl Anlss Dr Muhrrem
More informationLecture 36. Finite Element Methods
CE 60: Numercl Methods Lecture 36 Fnte Element Methods Course Coordntor: Dr. Suresh A. Krth, Assocte Professor, Deprtment of Cvl Engneerng, IIT Guwht. In the lst clss, we dscussed on the ppromte methods
More informationPolynomial Regression Models
LINEAR REGRESSION ANALYSIS MODULE XII Lecture - 6 Polynomal Regresson Models Dr. Shalabh Department of Mathematcs and Statstcs Indan Insttute of Technology Kanpur Test of sgnfcance To test the sgnfcance
More informationTHE COMBINED SHEPARD ABEL GONCHAROV UNIVARIATE OPERATOR
REVUE D ANALYSE NUMÉRIQUE ET DE THÉORIE DE L APPROXIMATION Tome 32, N o 1, 2003, pp 11 20 THE COMBINED SHEPARD ABEL GONCHAROV UNIVARIATE OPERATOR TEODORA CĂTINAŞ Abstrct We extend the Sheprd opertor by
More informationScatterplot Smoothing
1 Sttstcs 540, Smoothng Sctterplot Smoothng Overvew Problem... The usul settng for sctterplot smoothng s the delzed regresson model y = f(x )+σɛ, ɛ N(0, 1), =1,...,n, where the observtons (x,y ) re ndependent.
More informationPrinciple Component Analysis
Prncple Component Anlyss Jng Go SUNY Bufflo Why Dmensonlty Reducton? We hve too mny dmensons o reson bout or obtn nsghts from o vsulze oo much nose n the dt Need to reduce them to smller set of fctors
More informationPyramid Algorithms for Barycentric Rational Interpolation
Pyrmd Algorthms for Brycentrc Rtonl Interpolton K Hormnn Scott Schefer Astrct We present new perspectve on the Floter Hormnn nterpolnt. Ths nterpolnt s rtonl of degree (n, d), reproduces polynomls of degree
More informationInternational Journal of Pure and Applied Sciences and Technology
Int. J. Pure Appl. Sc. Technol., () (), pp. 44-49 Interntonl Journl of Pure nd Appled Scences nd Technolog ISSN 9-67 Avlle onlne t www.jopst.n Reserch Pper Numercl Soluton for Non-Lner Fredholm Integrl
More informationUNIVERSITY OF IOANNINA DEPARTMENT OF ECONOMICS. M.Sc. in Economics MICROECONOMIC THEORY I. Problem Set II
Mcroeconomc Theory I UNIVERSITY OF IOANNINA DEPARTMENT OF ECONOMICS MSc n Economcs MICROECONOMIC THEORY I Techng: A Lptns (Note: The number of ndctes exercse s dffculty level) ()True or flse? If V( y )
More information6 Roots of Equations: Open Methods
HK Km Slghtly modfed 3//9, /8/6 Frstly wrtten t Mrch 5 6 Roots of Equtons: Open Methods Smple Fed-Pont Iterton Newton-Rphson Secnt Methods MATLAB Functon: fzero Polynomls Cse Study: Ppe Frcton Brcketng
More informationRank One Update And the Google Matrix by Al Bernstein Signal Science, LLC
Introducton Rnk One Updte And the Google Mtrx y Al Bernsten Sgnl Scence, LLC www.sgnlscence.net here re two dfferent wys to perform mtrx multplctons. he frst uses dot product formulton nd the second uses
More informationChapter Newton-Raphson Method of Solving a Nonlinear Equation
Chpter.4 Newton-Rphson Method of Solvng Nonlner Equton After redng ths chpter, you should be ble to:. derve the Newton-Rphson method formul,. develop the lgorthm of the Newton-Rphson method,. use the Newton-Rphson
More informationThe Schur-Cohn Algorithm
Modelng, Estmton nd Otml Flterng n Sgnl Processng Mohmed Njm Coyrght 8, ISTE Ltd. Aendx F The Schur-Cohn Algorthm In ths endx, our m s to resent the Schur-Cohn lgorthm [] whch s often used s crteron for
More informationNice plotting of proteins II
Nce plottng of protens II Fnal remark regardng effcency: It s possble to wrte the Newton representaton n a way that can be computed effcently, usng smlar bracketng that we made for the frst representaton
More information7.2 Volume. A cross section is the shape we get when cutting straight through an object.
7. Volume Let s revew the volume of smple sold, cylnder frst. Cylnder s volume=se re heght. As llustrted n Fgure (). Fgure ( nd (c) re specl cylnders. Fgure () s rght crculr cylnder. Fgure (c) s ox. A
More informationNUMERICAL DIFFERENTIATION
NUMERICAL DIFFERENTIATION 1 Introducton Dfferentaton s a method to compute the rate at whch a dependent output y changes wth respect to the change n the ndependent nput x. Ths rate of change s called the
More information4. Eccentric axial loading, cross-section core
. Eccentrc xl lodng, cross-secton core Introducton We re strtng to consder more generl cse when the xl force nd bxl bendng ct smultneousl n the cross-secton of the br. B vrtue of Snt-Vennt s prncple we
More information( ) [ ( k) ( k) ( x) ( ) ( ) ( ) [ ] ξ [ ] [ ] [ ] ( )( ) i ( ) ( )( ) 2! ( ) = ( ) 3 Interpolation. Polynomial Approximation.
3 Interpolaton {( y } Gven:,,,,,, [ ] Fnd: y for some Mn, Ma Polynomal Appromaton Theorem (Weerstrass Appromaton Theorem --- estence ε [ ab] f( P( , then there ests a polynomal
More informationConvexity preserving interpolation by splines of arbitrary degree
Computer Scence Journal of Moldova, vol.18, no.1(52), 2010 Convexty preservng nterpolaton by splnes of arbtrary degree Igor Verlan Abstract In the present paper an algorthm of C 2 nterpolaton of dscrete
More informationJens Siebel (University of Applied Sciences Kaiserslautern) An Interactive Introduction to Complex Numbers
Jens Sebel (Unversty of Appled Scences Kserslutern) An Interctve Introducton to Complex Numbers 1. Introducton We know tht some polynoml equtons do not hve ny solutons on R/. Exmple 1.1: Solve x + 1= for
More informationChapter Newton-Raphson Method of Solving a Nonlinear Equation
Chpter 0.04 Newton-Rphson Method o Solvng Nonlner Equton Ater redng ths chpter, you should be ble to:. derve the Newton-Rphson method ormul,. develop the lgorthm o the Newton-Rphson method,. use the Newton-Rphson
More informationMany-Body Calculations of the Isotope Shift
Mny-Body Clcultons of the Isotope Shft W. R. Johnson Mrch 11, 1 1 Introducton Atomc energy levels re commonly evluted ssumng tht the nucler mss s nfnte. In ths report, we consder correctons to tomc levels
More informationFitting a Polynomial to Heat Capacity as a Function of Temperature for Ag. Mathematical Background Document
Fttng Polynol to Het Cpcty s Functon of Teperture for Ag. thetcl Bckground Docuent by Theres Jul Zelnsk Deprtent of Chestry, edcl Technology, nd Physcs onouth Unversty West ong Brnch, J 7764-898 tzelns@onouth.edu
More informationCISE 301: Numerical Methods Lecture 5, Topic 4 Least Squares, Curve Fitting
CISE 3: umercl Methods Lecture 5 Topc 4 Lest Squres Curve Fttng Dr. Amr Khouh Term Red Chpter 7 of the tetoo c Khouh CISE3_Topc4_Lest Squre Motvton Gven set of epermentl dt 3 5. 5.9 6.3 The reltonshp etween
More informationReview of linear algebra. Nuno Vasconcelos UCSD
Revew of lner lgebr Nuno Vsconcelos UCSD Vector spces Defnton: vector spce s set H where ddton nd sclr multplcton re defned nd stsf: ) +( + ) (+ )+ 5) λ H 2) + + H 6) 3) H, + 7) λ(λ ) (λλ ) 4) H, - + 8)
More informationRemember: Project Proposals are due April 11.
Bonformtcs ecture Notes Announcements Remember: Project Proposls re due Aprl. Clss 22 Aprl 4, 2002 A. Hdden Mrov Models. Defntons Emple - Consder the emple we tled bout n clss lst tme wth the cons. However,
More informationIntroduction to Vapor/Liquid Equilibrium, part 2. Raoult s Law:
CE304, Sprng 2004 Lecture 4 Introducton to Vapor/Lqud Equlbrum, part 2 Raoult s Law: The smplest model that allows us do VLE calculatons s obtaned when we assume that the vapor phase s an deal gas, and
More informationModule 2. Random Processes. Version 2 ECE IIT, Kharagpur
Module Random Processes Lesson 6 Functons of Random Varables After readng ths lesson, ou wll learn about cdf of functon of a random varable. Formula for determnng the pdf of a random varable. Let, X be
More informationApplied Statistics Qualifier Examination
Appled Sttstcs Qulfer Exmnton Qul_june_8 Fll 8 Instructons: () The exmnton contns 4 Questons. You re to nswer 3 out of 4 of them. () You my use ny books nd clss notes tht you mght fnd helpful n solvng
More informationKatholieke Universiteit Leuven Department of Computer Science
Updte Rules for Weghted Non-negtve FH*G Fctorzton Peter Peers Phlp Dutré Report CW 440, Aprl 006 Ktholeke Unverstet Leuven Deprtment of Computer Scence Celestjnenln 00A B-3001 Heverlee (Belgum) Updte Rules
More information8. INVERSE Z-TRANSFORM
8. INVERSE Z-TRANSFORM The proce by whch Z-trnform of tme ere, nmely X(), returned to the tme domn clled the nvere Z-trnform. The nvere Z-trnform defned by: Computer tudy Z X M-fle trn.m ued to fnd nvere
More informationIntroduction to Numerical Integration Part II
Introducton to umercl Integrton Prt II CS 75/Mth 75 Brn T. Smth, UM, CS Dept. Sprng, 998 4/9/998 qud_ Intro to Gussn Qudrture s eore, the generl tretment chnges the ntegrton prolem to ndng the ntegrl w
More informationLecture 2: Numerical Methods for Differentiations and Integrations
Numercal Smulaton of Space Plasmas (I [AP-4036] Lecture 2 by Lng-Hsao Lyu March, 2018 Lecture 2: Numercal Methods for Dfferentatons and Integratons As we have dscussed n Lecture 1 that numercal smulaton
More informationOnline Appendix to. Mandating Behavioral Conformity in Social Groups with Conformist Members
Onlne Appendx to Mndtng Behvorl Conformty n Socl Groups wth Conformst Members Peter Grzl Andrze Bnk (Correspondng uthor) Deprtment of Economcs, The Wllms School, Wshngton nd Lee Unversty, Lexngton, 4450
More informationTwo Coefficients of the Dyson Product
Two Coeffcents of the Dyson Product rxv:07.460v mth.co 7 Nov 007 Lun Lv, Guoce Xn, nd Yue Zhou 3,,3 Center for Combntorcs, LPMC TJKLC Nnk Unversty, Tnjn 30007, P.R. Chn lvlun@cfc.nnk.edu.cn gn@nnk.edu.cn
More informationMore metrics on cartesian products
More metrcs on cartesan products If (X, d ) are metrc spaces for 1 n, then n Secton II4 of the lecture notes we defned three metrcs on X whose underlyng topologes are the product topology The purpose of
More informationModelli Clamfim Equazioni differenziali 7 ottobre 2013
CLAMFIM Bologna Modell 1 @ Clamfm Equazon dfferenzal 7 ottobre 2013 professor Danele Rtell danele.rtell@unbo.t 1/18? Ordnary Dfferental Equatons A dfferental equaton s an equaton that defnes a relatonshp
More informationChemical Reaction Engineering
Lecture 20 hemcl Recton Engneerng (RE) s the feld tht studes the rtes nd mechnsms of chemcl rectons nd the desgn of the rectors n whch they tke plce. Lst Lecture Energy Blnce Fundmentls F 0 E 0 F E Q W
More informationPartially Observable Systems. 1 Partially Observable Markov Decision Process (POMDP) Formalism
CS294-40 Lernng for Rootcs nd Control Lecture 10-9/30/2008 Lecturer: Peter Aeel Prtlly Oservle Systems Scre: Dvd Nchum Lecture outlne POMDP formlsm Pont-sed vlue terton Glol methods: polytree, enumerton,
More informationMultiple view geometry
EECS 442 Computer vson Multple vew geometry Perspectve Structure from Moton - Perspectve structure from moton prolem - mgutes - lgerc methods - Fctorzton methods - Bundle djustment - Self-clrton Redng:
More informationModelli Clamfim Equazioni differenziali 22 settembre 2016
CLAMFIM Bologna Modell 1 @ Clamfm Equazon dfferenzal 22 settembre 2016 professor Danele Rtell danele.rtell@unbo.t 1/22? Ordnary Dfferental Equatons A dfferental equaton s an equaton that defnes a relatonshp
More informationALGORITHM FOR THE CALCULATION OF THE TWO VARIABLES CUBIC SPLINE FUNCTION
ANALELE ŞTIINŢIFICE ALE UNIVERSITĂŢII AL.I. CUZA DIN IAŞI (S.N.) MATEMATICĂ, Tomul LIX, 013, f.1 DOI: 10.478/v10157-01-00-y ALGORITHM FOR THE CALCULATION OF THE TWO VARIABLES CUBIC SPLINE FUNCTION BY ION
More informationChemical Reaction Engineering
Lecture 20 hemcl Recton Engneerng (RE) s the feld tht studes the rtes nd mechnsms of chemcl rectons nd the desgn of the rectors n whch they tke plce. Lst Lecture Energy Blnce Fundmentls F E F E + Q! 0
More informationLOCAL FRACTIONAL LAPLACE SERIES EXPANSION METHOD FOR DIFFUSION EQUATION ARISING IN FRACTAL HEAT TRANSFER
Yn, S.-P.: Locl Frctonl Lplce Seres Expnson Method for Dffuson THERMAL SCIENCE, Yer 25, Vol. 9, Suppl., pp. S3-S35 S3 LOCAL FRACTIONAL LAPLACE SERIES EXPANSION METHOD FOR DIFFUSION EQUATION ARISING IN
More informationMath 497C Sep 17, Curves and Surfaces Fall 2004, PSU
Mth 497C Sep 17, 004 1 Curves nd Surfces Fll 004, PSU Lecture Notes 3 1.8 The generl defnton of curvture; Fox-Mlnor s Theorem Let α: [, b] R n be curve nd P = {t 0,...,t n } be prtton of [, b], then the
More information1 Matrix representations of canonical matrices
1 Matrx representatons of canoncal matrces 2-d rotaton around the orgn: ( ) cos θ sn θ R 0 = sn θ cos θ 3-d rotaton around the x-axs: R x = 1 0 0 0 cos θ sn θ 0 sn θ cos θ 3-d rotaton around the y-axs:
More informationPhysics 5153 Classical Mechanics. Principle of Virtual Work-1
P. Guterrez 1 Introducton Physcs 5153 Classcal Mechancs Prncple of Vrtual Work The frst varatonal prncple we encounter n mechancs s the prncple of vrtual work. It establshes the equlbrum condton of a mechancal
More informationModeling curves. Graphs: y = ax+b, y = sin(x) Implicit ax + by + c = 0, x 2 +y 2 =r 2 Parametric:
Modelng curves Types of Curves Graphs: y = ax+b, y = sn(x) Implct ax + by + c = 0, x 2 +y 2 =r 2 Parametrc: x = ax + bxt x = cos t y = ay + byt y = snt Parametrc are the most common mplct are also used,
More informationAn Introduction to Support Vector Machines
An Introducton to Support Vector Mchnes Wht s good Decson Boundry? Consder two-clss, lnerly seprble clssfcton problem Clss How to fnd the lne (or hyperplne n n-dmensons, n>)? Any de? Clss Per Lug Mrtell
More informationModelli Clamfim Equazione del Calore Lezione ottobre 2014
CLAMFIM Bologna Modell 1 @ Clamfm Equazone del Calore Lezone 17 15 ottobre 2014 professor Danele Rtell danele.rtell@unbo.t 1/24? Convoluton The convoluton of two functons g(t) and f(t) s the functon (g
More informationLecture 13 APPROXIMATION OF SECOMD ORDER DERIVATIVES
COMPUTATIONAL FLUID DYNAMICS: FDM: Appromaton of Second Order Dervatves Lecture APPROXIMATION OF SECOMD ORDER DERIVATIVES. APPROXIMATION OF SECOND ORDER DERIVATIVES Second order dervatves appear n dffusve
More informationReproducing Kernel Hilbert Space for. Penalized Regression Multi-Predictors: Case in Longitudinal Data
Interntonl Journl of Mthemtcl Anlyss Vol. 8, 04, no. 40, 95-96 HIKARI Ltd, www.m-hkr.com http://dx.do.org/0.988/jm.04.47 Reproducng Kernel Hlbert Spce for Penlzed Regresson Mult-Predctors: Cse n Longudnl
More informationDemand. Demand and Comparative Statics. Graphically. Marshallian Demand. ECON 370: Microeconomic Theory Summer 2004 Rice University Stanley Gilbert
Demnd Demnd nd Comrtve Sttcs ECON 370: Mcroeconomc Theory Summer 004 Rce Unversty Stnley Glbert Usng the tools we hve develoed u to ths ont, we cn now determne demnd for n ndvdul consumer We seek demnd
More informationEconomics 130. Lecture 4 Simple Linear Regression Continued
Economcs 130 Lecture 4 Contnued Readngs for Week 4 Text, Chapter and 3. We contnue wth addressng our second ssue + add n how we evaluate these relatonshps: Where do we get data to do ths analyss? How do
More informationPoint-based methods for estimating the length of a parametric curve
Journl of Computtonl nd Appled Mthemtcs 196 (006) 51 5 www.elsever.com/locte/cm Pont-bsed methods for estmtng the length of prmetrc curve Mchel S. Floter, Atgerr F. Rsmussen CMA, Unversty of Oslo, P.O.
More informationGAUSS ELIMINATION. Consider the following system of algebraic linear equations
Numercl Anlyss for Engneers Germn Jordnn Unversty GAUSS ELIMINATION Consder the followng system of lgebrc lner equtons To solve the bove system usng clsscl methods, equton () s subtrcted from equton ()
More informationLAPLACE TRANSFORM SOLUTION OF THE PROBLEM OF TIME-FRACTIONAL HEAT CONDUCTION IN A TWO-LAYERED SLAB
Journl of Appled Mthemtcs nd Computtonl Mechncs 5, 4(4), 5-3 www.mcm.pcz.pl p-issn 99-9965 DOI:.75/jmcm.5.4. e-issn 353-588 LAPLACE TRANSFORM SOLUTION OF THE PROBLEM OF TIME-FRACTIONAL HEAT CONDUCTION
More informationHaddow s Experiment:
schemtc drwng of Hddow's expermentl set-up movng pston non-contctng moton sensor bems of sprng steel poston vres to djust frequences blocks of sold steel shker Hddow s Experment: terr frm Theoretcl nd
More informationCMDA 4604: Intermediate Topics in Mathematical Modeling Lecture 19: Interpolation and Quadrature
CMDA 4604: Intermedite Topics in Mthemticl Modeling Lecture 19: Interpoltion nd Qudrture In this lecture we mke brief diversion into the res of interpoltion nd qudrture. Given function f C[, b], we sy
More informationModule 3: Element Properties Lecture 1: Natural Coordinates
Module 3: Element Propertes Lecture : Natural Coordnates Natural coordnate system s bascally a local coordnate system whch allows the specfcaton of a pont wthn the element by a set of dmensonless numbers
More information18.7 Artificial Neural Networks
310 18.7 Artfcl Neurl Networks Neuroscence hs hypotheszed tht mentl ctvty conssts prmrly of electrochemcl ctvty n networks of brn cells clled neurons Ths led McCulloch nd Ptts to devse ther mthemtcl model
More informationPART 1: VECTOR & TENSOR ANALYSIS
PART : VECTOR & TENSOR ANALYSIS wth LINEAR ALGEBRA Obectves Introduce the concepts, theores, nd opertonl mplementton of vectors, nd more generlly tensors, n dvnced engneerng nlyss. The emphss s on geometrc
More informationON SIMPSON S INEQUALITY AND APPLICATIONS. 1. Introduction The following inequality is well known in the literature as Simpson s inequality : 2 1 f (4)
ON SIMPSON S INEQUALITY AND APPLICATIONS SS DRAGOMIR, RP AGARWAL, AND P CERONE Abstrct New neultes of Smpson type nd ther pplcton to udrture formule n Numercl Anlyss re gven Introducton The followng neulty
More information8.1 Arc Length. What is the length of a curve? How can we approximate it? We could do it following the pattern we ve used before
.1 Arc Length hat s the length of a curve? How can we approxmate t? e could do t followng the pattern we ve used before Use a sequence of ncreasngly short segments to approxmate the curve: As the segments
More informationPHYS 705: Classical Mechanics. Calculus of Variations II
1 PHYS 705: Classcal Mechancs Calculus of Varatons II 2 Calculus of Varatons: Generalzaton (no constrant yet) Suppose now that F depends on several dependent varables : We need to fnd such that has a statonary
More informationINTRODUCTION TO COMPLEX NUMBERS
INTRODUCTION TO COMPLEX NUMBERS The numers -4, -3, -, -1, 0, 1,, 3, 4 represent the negtve nd postve rel numers termed ntegers. As one frst lerns n mddle school they cn e thought of s unt dstnce spced
More informationLeast squares cubic splines without B-splines S.K. Lucas
Least squares cubc splnes wthout B-splnes S.K. Lucas School of Mathematcs and Statstcs, Unversty of South Australa, Mawson Lakes SA 595 e-mal: stephen.lucas@unsa.edu.au Submtted to the Gazette of the Australan
More informationStatistics 423 Midterm Examination Winter 2009
Sttstcs 43 Mdterm Exmnton Wnter 009 Nme: e-ml: 1. Plese prnt your nme nd e-ml ddress n the bove spces.. Do not turn ths pge untl nstructed to do so. 3. Ths s closed book exmnton. You my hve your hnd clcultor
More informationON A DETERMINATION OF THE INITIAL FUNCTIONS FROM THE OBSERVED VALUES OF THE BOUNDARY FUNCTIONS FOR THE SECOND-ORDER HYPERBOLIC EQUATION
Advanced Mathematcal Models & Applcatons Vol.3, No.3, 2018, pp.215-222 ON A DETERMINATION OF THE INITIAL FUNCTIONS FROM THE OBSERVED VALUES OF THE BOUNDARY FUNCTIONS FOR THE SECOND-ORDER HYPERBOLIC EUATION
More informationInner Product. Euclidean Space. Orthonormal Basis. Orthogonal
Inner Product Defnton 1 () A Eucldean space s a fnte-dmensonal vector space over the reals R, wth an nner product,. Defnton 2 (Inner Product) An nner product, on a real vector space X s a symmetrc, blnear,
More informationCHI-SQUARE DIVERGENCE AND MINIMIZATION PROBLEM
CHI-SQUARE DIVERGENCE AND MINIMIZATION PROBLEM PRANESH KUMAR AND INDER JEET TANEJA Abstrct The mnmum dcrmnton nformton prncple for the Kullbck-Lebler cross-entropy well known n the lterture In th pper
More informationn f(x i ) x. i=1 In section 4.2, we defined the definite integral of f from x = a to x = b as n f(x i ) x; f(x) dx = lim i=1
The Fundmentl Theorem of Clculus As we continue to study the re problem, let s think bck to wht we know bout computing res of regions enclosed by curves. If we wnt to find the re of the region below the
More informationProfessor Terje Haukaas University of British Columbia, Vancouver The Q4 Element
Professor Terje Haukaas Unversty of Brtsh Columba, ancouver www.nrsk.ubc.ca The Q Element Ths document consders fnte elements that carry load only n ther plane. These elements are sometmes referred to
More informationQuiz: Experimental Physics Lab-I
Mxmum Mrks: 18 Totl tme llowed: 35 mn Quz: Expermentl Physcs Lb-I Nme: Roll no: Attempt ll questons. 1. In n experment, bll of mss 100 g s dropped from heght of 65 cm nto the snd contner, the mpct s clled
More informationHAMILTON-JACOBI TREATMENT OF LAGRANGIAN WITH FERMIONIC AND SCALAR FIELD
AMION-JACOBI REAMEN OF AGRANGIAN WI FERMIONIC AND SCAAR FIED W. I. ESRAIM 1, N. I. FARAA Dertment of Physcs, Islmc Unversty of Gz, P.O. Box 18, Gz, Plestne 1 wbrhm 7@hotml.com nfrht@ugz.edu.s Receved November,
More information4. More general extremum principles and thermodynamic potentials
4. More generl etremum prncples nd thermodynmc potentls We hve seen tht mn{u(s, X )} nd m{s(u, X)} mply one nother. Under certn condtons, these prncples re very convenent. For emple, ds = 1 T du T dv +
More information10-801: Advanced Optimization and Randomized Methods Lecture 2: Convex functions (Jan 15, 2014)
0-80: Advanced Optmzaton and Randomzed Methods Lecture : Convex functons (Jan 5, 04) Lecturer: Suvrt Sra Addr: Carnege Mellon Unversty, Sprng 04 Scrbes: Avnava Dubey, Ahmed Hefny Dsclamer: These notes
More informationThe Quadratic Trigonometric Bézier Curve with Single Shape Parameter
J. Basc. Appl. Sc. Res., (3541-546, 01 01, TextRoad Publcaton ISSN 090-4304 Journal of Basc and Appled Scentfc Research www.textroad.com The Quadratc Trgonometrc Bézer Curve wth Sngle Shape Parameter Uzma
More informationINSTITUTE OF AERONAUTICAL ENGINEERING Dundigal, Hyderabad
INSTITUTE OF AERONAUTICAL ENGINEERING Dundgl, Hyderbd - 5 3 FRESHMAN ENGINEERING TUTORIAL QUESTION BANK Nme : MATHEMATICS II Code : A6 Clss : II B. Te II Semester Brn : FRESHMAN ENGINEERING Yer : 5 Fulty
More informationVariable time amplitude amplification and quantum algorithms for linear algebra. Andris Ambainis University of Latvia
Vrble tme mpltude mplfcton nd quntum lgorthms for lner lgebr Andrs Ambns Unversty of Ltv Tlk outlne. ew verson of mpltude mplfcton;. Quntum lgorthm for testng f A s sngulr; 3. Quntum lgorthm for solvng
More informationStudy of Trapezoidal Fuzzy Linear System of Equations S. M. Bargir 1, *, M. S. Bapat 2, J. D. Yadav 3 1
mercn Interntonl Journl of Reserch n cence Technology Engneerng & Mthemtcs vlble onlne t http://wwwsrnet IN (Prnt: 38-349 IN (Onlne: 38-3580 IN (CD-ROM: 38-369 IJRTEM s refereed ndexed peer-revewed multdscplnry
More information90 S.S. Drgomr nd (t b)du(t) =u()(b ) u(t)dt: If we dd the bove two equltes, we get (.) u()(b ) u(t)dt = p(; t)du(t) where p(; t) := for ll ; t [; b]:
RGMIA Reserch Report Collecton, Vol., No. 1, 1999 http://sc.vu.edu.u/οrgm ON THE OSTROWSKI INTEGRAL INEQUALITY FOR LIPSCHITZIAN MAPPINGS AND APPLICATIONS S.S. Drgomr Abstrct. A generlzton of Ostrowsk's
More informationZbus 1.0 Introduction The Zbus is the inverse of the Ybus, i.e., (1) Since we know that
us. Introducton he us s the nverse of the us,.e., () Snce we now tht nd therefore then I V () V I () V I (4) So us reltes the nodl current njectons to the nodl voltges, s seen n (4). In developng the power
More informationModule 3 LOSSY IMAGE COMPRESSION SYSTEMS. Version 2 ECE IIT, Kharagpur
Module 3 LOSSY IMAGE COMPRESSION SYSTEMS Verson ECE IIT, Kharagpur Lesson 6 Theory of Quantzaton Verson ECE IIT, Kharagpur Instructonal Objectves At the end of ths lesson, the students should be able to:
More informationCanonical transformations
Canoncal transformatons November 23, 2014 Recall that we have defned a symplectc transformaton to be any lnear transformaton M A B leavng the symplectc form nvarant, Ω AB M A CM B DΩ CD Coordnate transformatons,
More information6.6 The Marquardt Algorithm
6.6 The Mqudt Algothm lmttons of the gdent nd Tylo expnson methods ecstng the Tylo expnson n tems of ch-sque devtves ecstng the gdent sech nto n tetve mtx fomlsm Mqudt's lgothm utomtclly combnes the gdent
More informationDennis Bricker, 2001 Dept of Industrial Engineering The University of Iowa. MDP: Taxi page 1
Denns Brcker, 2001 Dept of Industrl Engneerng The Unversty of Iow MDP: Tx pge 1 A tx serves three djcent towns: A, B, nd C. Ech tme the tx dschrges pssenger, the drver must choose from three possble ctons:
More informationChapter 5. Solution of System of Linear Equations. Module No. 6. Solution of Inconsistent and Ill Conditioned Systems
Numercal Analyss by Dr. Anta Pal Assstant Professor Department of Mathematcs Natonal Insttute of Technology Durgapur Durgapur-713209 emal: anta.bue@gmal.com 1 . Chapter 5 Soluton of System of Lnear Equatons
More informationAPPENDIX A Some Linear Algebra
APPENDIX A Some Lnear Algebra The collecton of m, n matrces A.1 Matrces a 1,1,..., a 1,n A = a m,1,..., a m,n wth real elements a,j s denoted by R m,n. If n = 1 then A s called a column vector. Smlarly,
More informationJean Fernand Nguema LAMETA UFR Sciences Economiques Montpellier. Abstract
Stochstc domnnce on optml portfolo wth one rsk less nd two rsky ssets Jen Fernnd Nguem LAMETA UFR Scences Economques Montpeller Abstrct The pper provdes restrctons on the nvestor's utlty functon whch re
More informationPhysics 5153 Classical Mechanics. D Alembert s Principle and The Lagrangian-1
P. Guterrez Physcs 5153 Classcal Mechancs D Alembert s Prncple and The Lagrangan 1 Introducton The prncple of vrtual work provdes a method of solvng problems of statc equlbrum wthout havng to consder the
More informationNumerical Heat and Mass Transfer
Master degree n Mechancal Engneerng Numercal Heat and Mass Transfer 06-Fnte-Dfference Method (One-dmensonal, steady state heat conducton) Fausto Arpno f.arpno@uncas.t Introducton Why we use models and
More informationLecture notes. Fundamental inequalities: techniques and applications
Lecture notes Fundmentl nequltes: technques nd pplctons Mnh Hong Duong Mthemtcs Insttute, Unversty of Wrwck Eml: m.h.duong@wrwck.c.uk Jnury 4, 07 Abstrct Inequltes re ubqutous n Mthemtcs (nd n rel lfe.
More informationFTCS Solution to the Heat Equation
FTCS Soluton to the Heat Equaton ME 448/548 Notes Gerald Recktenwald Portland State Unversty Department of Mechancal Engneerng gerry@pdx.edu ME 448/548: FTCS Soluton to the Heat Equaton Overvew 1. Use
More informationCalculus of Variations Basics
Chapter 1 Calculus of Varatons Bascs 1.1 Varaton of a General Functonal In ths chapter, we derve the general formula for the varaton of a functonal of the form J [y 1,y 2,,y n ] F x,y 1,y 2,,y n,y 1,y
More informationNumerical Analysis: Trapezoidal and Simpson s Rule
nd Simpson s Mthemticl question we re interested in numericlly nswering How to we evlute I = f (x) dx? Clculus tells us tht if F(x) is the ntiderivtive of function f (x) on the intervl [, b], then I =
More informationGeorgia Tech PHYS 6124 Mathematical Methods of Physics I
Georga Tech PHYS 624 Mathematcal Methods of Physcs I Instructor: Predrag Cvtanovć Fall semester 202 Homework Set #7 due October 30 202 == show all your work for maxmum credt == put labels ttle legends
More informationComputation of Fifth Degree of Spline Function Model by Using C++ Programming
www.ijci.org 89 Computton o Ft Degree o plne Functon Model b Usng C Progrmmng Frdun K. Hml, Aln A. Abdull nd Knd M. Qdr Mtemtcs Dept, Unverst o ulmn, ulmn, IRAQ Mtemtcs Dept, Unverst o ulmn, ulmn, IRAQ
More information