Outline. Review Small Grid (N = 6, M = 5) Review Finite Differences. Review N = 6, M = 5 Matrix. More numerical elliptic PDEs March 30, 2009

Size: px
Start display at page:

Download "Outline. Review Small Grid (N = 6, M = 5) Review Finite Differences. Review N = 6, M = 5 Matrix. More numerical elliptic PDEs March 30, 2009"

Transcription

1 ore nmercl ellptc Ds rch, 9 5 nneern nlss ddtonl Topcs n mercl ddtonl Topcs n mercl oltons o llptc qtons oltons o llptc qtons Lrr Cretto echncl nneern 5 emnr n nneern nlss rch, 9 tlne Revew lst clss Tretment o bondr condtons econd knd emnn Thrd knd /s b c s s coordnte norml to bondr or Tretment sn derent rdent epressons Compct derence epressons or hher order ccrc Revew Fnte Derences econd-order second dervtves ] [ ε k k & & osson-lplce eqton [, ] k & ltpl b nd dene / Revew mll Grd 6, ondr nodes 5 5 Compttonl olecle 5 Revew 6, 5 tr Zero coecents becse o bondr - 5 6

2 ore nmercl ellptc Ds rch, 9 Revew Itertve oltons Jcob terton ses ll old vles ' ' n ' n ' n b n ' n Gss edel ses most-recent vles ' ' n ' n ' n b n ' n Relton n ' n n ω ω[ ' n ' n ' n ' b ] 7 ecton Tme seconds. ect o Relton Fctor on ecton Tme ξ qre L H 6 b 6 Grd Zero bondr on let, rht nd bottom Top bondr hs,h snπ/l "ther" s derent code wth,h ξ ω opt ξ π cos π cos 66 rd rd 66 rd 88 rd ther code Relton Fctor Revew rrors n n Reltve Chne n Comptble n Resdl ' n n ' n n n Iterton rror ' n n ' n n n ct rror, D ' b 9 rrors ects o Itertons on Lplce qton rrors qre L H 6 b 6 Grd Zero bondr on let, rht nd bottom Top bondr hs,h snπ/l Derence Resdl Iterton error ct rror Itertons ther ondr Condtons Generl condton /s b c, b or Drchlet vle ven, b or emnn rdent ven ed hs both nd b nonzero rte enerl bondr condton sn nte derence epresson or /s Two pproches Usn second order orwrd or bckwrd derence or bondr node dd cttos node otsde bondr nd se centrl derences t bondr otton Lbel the bondres t nd s the est nd st bondres, b, nd c cn be derent t ech otton:, b, c,, b, nd c ondres t nd re the oth nd orth bondres, b, nd c cn be derent t ech otton:, b, c,, b, nd c 5 nneern nlss

3 ore nmercl ellptc Ds rch, 9 Centrl Derence pproch Generl condton /s b c b centrl derences rte enerl eqton or bondr nodes tht contns cttos node otsde reon rte centrl derence eqton or bondr condton olve ths eqton or potentl t cttos node nd se reslt to replce ths vle n enerl eqton t bondr Inclde bondr nodes n tertons orth ondr mple Fnte-derence eqton t Fcttos rom centrl-derence bondr condton eqton t b c c b orth ondr mple II Combne eqtons to elmnte c b b oded coecents t terted bondres: new,, nd ; ; nd mltpled b c 5 CD Generl ondr oded est bondr coecents b c oded north bondr coecents b c 6 CD Generl ondr II oded west bondr coecents b c 7 oded soth bondr coecents b c Forwrd/ckwrd Derence Generl condton, /s b c Use orwrd derences t or nd bckwrd derences t or btn eqton or bondr potentl n terms o two nodes n rom bondr Combne ths eqton wth enerl eqton or rst node n rom the bondr to elmnte nknown bondr potentl o terton on bondr vles, whch re ond t end o tertons Dervton t end o ths presentton 8 5 nneern nlss

4 ore nmercl ellptc Ds rch, 9 5 nneern nlss 9 est nd oth ondr b c b b b c b c b b b c nchned nd nchned nd orth nd st ondr b c b b b c b c b b b c nchned nd nchned nd Compct ppromtons Dene centrl-derence opertor, [ ] Compct second dervtve To ppl ths eqton, we hve to mltpl throh b denomntor nd then ppl the opertor there Fnl Reslt ch vle o s lnked to eht nerest nehbors so ech nte derence eqton hs nne terms Uses wehted vere sorce term Dervton t end o ths presentton [ ] Fnl Reslt or h For, [ ] [ ] h ondr Crosses Grd Dene new nd to dene bondr Get dervtve epressons or neven rd spcn...! '' ' '' h

5 ore nmercl ellptc Ds rch, 9 5 nneern nlss 5 5 ondr Crosses Grd II mpl eqton s ollows '' '' Use or Lplce eqton /, / 6 ondr Crosses Grd III 7 ondr Crosses Grd IV Homn ses wth derent notton Dene nd - - Dene nd ondr Crosses Grd V Use neven nte-derence epressons n derentl eqtons Cn crete problems wth stblt n eplct procedres Cre s reqred n modeln rdent bondr condtons Generll not vored cepton s Flow-D sotwre b C.. Ton Hrt who recommends ths procedre 9 mmr llptc nte-derence ormed n the sme w s prbolc ones qtons reqre tertve solvers mple solvers wll work well or smll problems dvnced solvers reqred or more comple problems eed to tret bondres tht do not ll on rd nodes dterm et ednesd, prl Covers wve eqton, Ds n more thn two ndependent vrbles, clsscton o Ds nd nmercl nlss ntrodcton roblems smlr to prevos mdterm nd homework n qestons?

6 ore nmercl ellptc Ds rch, 9 5 nneern nlss 6 ddtonl terls ldes to present detled dervtons o the enerl bondr condtons sn one-sded derences ee Homn or other dervtons o centrl-derence bondr eqtons The dervton o the hher order derence method sn compct derences s shown on sldes to 7 Generl ondr mple Look t est bondr; other bondres ollow smlr dervton olve one-sded, second-order rst dervtve epresson or c b b c b b c Generl ondr mple II Combne bondr condton wth eqton or rst node rom bondr b c b c b b Generl ondr mple III ondr node no loner n eqton tll solve or nteror nodes onl wth solver Fnd bondr vle ter solton complete Use moded eqton coecents b c b b b c 5 Reslts or ther ondres bsttte revsed coecents or ornl coecents n nte-derence verson o D mlr pproch, wth derent eqtons, t ll bondres Iterte onl on nonbondr nodes nd nd bondr potentls ter tertons complete Dervtons ollow here 6 oth ondr Reslt smlr to est bondr b c b c b b

7 ore nmercl ellptc Ds rch, 9 5 nneern nlss 7 7 orth nd st ondres me pproch wth bckwrd derence epressons Look t st s emple c b b c b b c 8 st ondr Combne bondr condton wth eqton or lst node beore bondr b c b c b b 9 st ondr II Fnl reslts b c b c b b b c b b orth ondr Reslt or bondr smlr to b c b c b b D Compct Derences Compct derences or osson eqton Dene nd drecton derence opertors econd prtl dervtve epressons [ ] [ ] ppl to / /, osson Compct Derences ltpl b / nd se / e know nd, bt wht s

8 ore nmercl ellptc Ds rch, 9 5 nneern nlss 8 ed Derence pertors ppl ech opertor n order Get sme reslt rerdless o order ck to rnl qton Redce bsc nte-derence eqton 5 rnl qton orce Term Contne redcton o bsc eqton Combne Reslts Dvde b nd rerrne Combne Reslts II 6 5 5

Outline. Review Quadrilateral Equation. Review Linear φ i Quadrilateral. Review x and y Derivatives. Review φ Derivatives

Outline. Review Quadrilateral Equation. Review Linear φ i Quadrilateral. Review x and y Derivatives. Review φ Derivatives E 5 Engneerng nlss ore on Fnte Eleents n ore on Fnte Eleents n Two Densons Two Densons Lrr Cretto echncl Engneerng 5 Senr n Engneerng nlss prl 7-9 9 Otlne Revew lst lectre Qrtc ss nctons n two ensons orer

More information

Chapter Newton-Raphson Method of Solving a Nonlinear Equation

Chapter 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 information

ME 501A Seminar in Engineering Analysis Page 1

ME 501A Seminar in Engineering Analysis Page 1 More oundr-vlue Prolems nd genvlue Prolems n Os ovemer 9, 7 More oundr-vlue Prolems nd genvlue Prolems n Os Lrr retto Menl ngneerng 5 Semnr n ngneerng nlss ovemer 9, 7 Outlne Revew oundr-vlue prolems Soot

More information

ORDINARY DIFFERENTIAL EQUATIONS

ORDINARY DIFFERENTIAL EQUATIONS 6 ORDINARY DIFFERENTIAL EQUATIONS Introducton Runge-Kutt Metods Mult-step Metods Sstem o Equtons Boundr Vlue Problems Crcterstc Vlue Problems Cpter 6 Ordnr Derentl Equtons / 6. Introducton In mn engneerng

More information

Chapter Newton-Raphson Method of Solving a Nonlinear Equation

Chapter 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 information

Engineering Tensors. Friday November 16, h30 -Muddy Charles. A BEH430 review session by Thomas Gervais.

Engineering Tensors. Friday November 16, h30 -Muddy Charles. A BEH430 review session by Thomas Gervais. ngneerng Tensors References: BH4 reew sesson b Thoms Gers tgers@mt.ed Long, RR, Mechncs of Solds nd lds, Prentce-Hll, 96, pp - Deen, WD, nlss of trnsport phenomen, Oford, 998, p. 55-56 Goodbod, M, Crtesn

More information

Introduction to Numerical Integration Part II

Introduction 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 information

Outline. Finite Difference Grids. Numerical Analysis. Finite Difference Grids II. Finite Difference Grids III

Outline. Finite Difference Grids. Numerical Analysis. Finite Difference Grids II. Finite Difference Grids III Itrodcto to Nmercl Alyss Mrc, 9 Nmercl Metods or PDEs Lrry Cretto Meccl Egeerg 5B Semr Egeerg Alyss Mrc, 9 Otle Revew mdterm soltos Revew bsc mterl o mercl clcls Expressos or dervtves, error d error order

More information

consider in the case of 1) internal resonance ω 2ω and 2) external resonance Ω ω and small damping

consider in the case of 1) internal resonance ω 2ω and 2) external resonance Ω ω and small damping consder n the cse o nternl resonnce nd externl resonnce Ω nd smll dmpng recll rom "Two_Degs_Frdm_.ppt" tht θ + μ θ + θ = θφ + cos Ω t + τ where = k α α nd φ + μ φ + φ = θ + cos Ω t where = α τ s constnt

More information

Chapter Runge-Kutta 2nd Order Method for Ordinary Differential Equations

Chapter Runge-Kutta 2nd Order Method for Ordinary Differential Equations Cter. Runge-Kutt nd Order Metod or Ordnr Derentl Eutons Ater redng ts cter ou sould be ble to:. understnd te Runge-Kutt nd order metod or ordnr derentl eutons nd ow to use t to solve roblems. Wt s te Runge-Kutt

More information

Dennis Bricker, 2001 Dept of Industrial Engineering The University of Iowa. MDP: Taxi page 1

Dennis 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 information

Principle Component Analysis

Principle 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 information

6 Roots of Equations: Open Methods

6 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 information

8. INVERSE Z-TRANSFORM

8. 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 information

The Number of Rows which Equal Certain Row

The Number of Rows which Equal Certain Row Interntonl Journl of Algebr, Vol 5, 011, no 30, 1481-1488 he Number of Rows whch Equl Certn Row Ahmd Hbl Deprtment of mthemtcs Fcult of Scences Dmscus unverst Dmscus, Sr hblhmd1@gmlcom Abstrct Let be X

More information

International Journal of Pure and Applied Sciences and Technology

International 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 information

8 factors of x. For our second example, let s raise a power to a power:

8 factors of x. For our second example, let s raise a power to a power: CH 5 THE FIVE LAWS OF EXPONENTS EXPONENTS WITH VARIABLES It s no time for chnge in tctics, in order to give us deeper understnding of eponents. For ech of the folloing five emples, e ill stretch nd squish,

More information

Jackson 2.26 Homework Problem Solution Dr. Christopher S. Baird University of Massachusetts Lowell

Jackson 2.26 Homework Problem Solution Dr. Christopher S. Baird University of Massachusetts Lowell Jckson 2.26 Homework Problem Solution Dr. Christopher S. Bird University of Msschusetts Lowell PROBLEM: The two-dimensionl region, ρ, φ β, is bounded by conducting surfces t φ =, ρ =, nd φ = β held t zero

More information

Partially Observable Systems. 1 Partially Observable Markov Decision Process (POMDP) Formalism

Partially 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 information

ME 501A Seminar in Engineering Analysis Page 1

ME 501A Seminar in Engineering Analysis Page 1 umercal Solutons of oundary-value Problems n Os ovember 7, 7 umercal Solutons of oundary- Value Problems n Os Larry aretto Mechancal ngneerng 5 Semnar n ngneerng nalyss ovember 7, 7 Outlne Revew stff equaton

More information

a = Acceleration Linear Motion Acceleration Changing Velocity All these Velocities? Acceleration and Freefall Physics 114

a = Acceleration Linear Motion Acceleration Changing Velocity All these Velocities? Acceleration and Freefall Physics 114 Lner Accelerton nd Freell Phyc 4 Eyre Denton o ccelerton Both de o equton re equl Mgntude Unt Drecton (t ector!) Accelerton Mgntude Mgntude Unt Unt Drecton Drecton 4/3/07 Module 3-Phy4-Eyre 4/3/07 Module

More information

Definition of Tracking

Definition of Tracking Trckng Defnton of Trckng Trckng: Generte some conclusons bout the moton of the scene, objects, or the cmer, gven sequence of mges. Knowng ths moton, predct where thngs re gong to project n the net mge,

More information

M344 - ADVANCED ENGINEERING MATHEMATICS

M344 - ADVANCED ENGINEERING MATHEMATICS M3 - ADVANCED ENGINEERING MATHEMATICS Lecture 18: Lplce s Eqution, Anltic nd Numericl Solution Our emple of n elliptic prtil differentil eqution is Lplce s eqution, lso clled the Diffusion Eqution. If

More information

Review of linear algebra. Nuno Vasconcelos UCSD

Review 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 information

Finite Difference Method

Finite Difference Method 7/0/07 Instructor r. Ramond Rump (9) 747 698 rcrump@utep.edu EE 337 Computatonal Electromagnetcs (CEM) Lecture #0 Fnte erence Method Lecture 0 These notes ma contan coprghted materal obtaned under ar use

More information

Course Review Introduction to Computer Methods

Course Review Introduction to Computer Methods Course Revew Wht you hopefully hve lerned:. How to nvgte nsde MIT computer system: Athen, UNIX, emcs etc. (GCR). Generl des bout progrmmng (GCR): formultng the problem, codng n Englsh trnslton nto computer

More information

FE REVIEW OPERATIONAL AMPLIFIERS (OP-AMPS)

FE REVIEW OPERATIONAL AMPLIFIERS (OP-AMPS) FE EIEW OPEATIONAL AMPLIFIES (OPAMPS) 1 The Opamp An opamp has two nputs and one output. Note the opamp below. The termnal labeled wth the () sgn s the nvertng nput and the nput labeled wth the () sgn

More information

September 13 Homework Solutions

September 13 Homework Solutions College of Engineering nd Computer Science Mechnicl Engineering Deprtment Mechnicl Engineering 5A Seminr in Engineering Anlysis Fll Ticket: 5966 Instructor: Lrry Cretto Septemer Homework Solutions. Are

More information

LOCAL FRACTIONAL LAPLACE SERIES EXPANSION METHOD FOR DIFFUSION EQUATION ARISING IN FRACTAL HEAT TRANSFER

LOCAL 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 information

Censored Models. Li Gan. April, Examples of censored regression: (1) stadium attendance= min(attendance *, capacity)

Censored Models. Li Gan. April, Examples of censored regression: (1) stadium attendance= min(attendance *, capacity) Censore Moels L Gn Aprl 8 mples o ensore regresson: stm ttenne mn(ttenne pt) Here gropng (those oservtons wth ttenne rtes less thn one vs those oservtons wth ttenne rtes t one) re lerl se on oserve tors.

More information

VECTORS VECTORS VECTORS VECTORS. 2. Vector Representation. 1. Definition. 3. Types of Vectors. 5. Vector Operations I. 4. Equal and Opposite Vectors

VECTORS VECTORS VECTORS VECTORS. 2. Vector Representation. 1. Definition. 3. Types of Vectors. 5. Vector Operations I. 4. Equal and Opposite Vectors 1. Defnton A vetor s n entt tht m represent phsl quntt tht hs mgntude nd dreton s opposed to slr tht ls dreton.. Vetor Representton A vetor n e represented grphll n rrow. The length of the rrow s the mgntude

More information

Lecture 4: Piecewise Cubic Interpolation

Lecture 4: Piecewise Cubic Interpolation 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

More information

AE/ME 339. K. M. Isaac. 8/31/2004 topic4: Implicit method, Stability, ADI method. Computational Fluid Dynamics (AE/ME 339) MAEEM Dept.

AE/ME 339. K. M. Isaac. 8/31/2004 topic4: Implicit method, Stability, ADI method. Computational Fluid Dynamics (AE/ME 339) MAEEM Dept. AE/ME 339 Comptatonal Fld Dynamcs (CFD) Comptatonal Fld Dynamcs (AE/ME 339) Implct form of dfference eqaton In the prevos explct method, the solton at tme level n,,n, depended only on the known vales of,

More information

Variable time amplitude amplification and quantum algorithms for linear algebra. Andris Ambainis University of Latvia

Variable 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 information

12 Iterative Methods. Linear Systems: Gauss-Seidel Nonlinear Systems Case Study: Chemical Reactions

12 Iterative Methods. Linear Systems: Gauss-Seidel Nonlinear Systems Case Study: Chemical Reactions HK Km Slghtly moded //9 /8/6 Frstly wrtte t Mrch 5 Itertve Methods er Systems: Guss-Sedel Noler Systems Cse Study: Chemcl Rectos Itertve or ppromte methods or systems o equtos cosst o guessg vlue d the

More information

ESCI 342 Atmospheric Dynamics I Lesson 1 Vectors and Vector Calculus

ESCI 342 Atmospheric Dynamics I Lesson 1 Vectors and Vector Calculus ESI 34 tmospherc Dnmcs I Lesson 1 Vectors nd Vector lculus Reference: Schum s Outlne Seres: Mthemtcl Hndbook of Formuls nd Tbles Suggested Redng: Mrtn Secton 1 OORDINTE SYSTEMS n orthonorml coordnte sstem

More information

UNIVERSITY OF IOANNINA DEPARTMENT OF ECONOMICS. M.Sc. in Economics MICROECONOMIC THEORY I. Problem Set II

UNIVERSITY 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 information

8 Laplace s Method and Local Limit Theorems

8 Laplace s Method and Local Limit Theorems 8 Lplce s Method nd Locl Limit Theorems 8. Fourier Anlysis in Higher DImensions Most of the theorems of Fourier nlysis tht we hve proved hve nturl generliztions to higher dimensions, nd these cn be proved

More information

PART 8. Partial Differential Equations PDEs

PART 8. Partial Differential Equations PDEs he Islamc Unverst of Gaza Facult of Engneerng Cvl Engneerng Department Numercal Analss ECIV 3306 PAR 8 Partal Dfferental Equatons PDEs Chapter 9; Fnte Dfference: Ellptc Equatons Assocate Prof. Mazen Abualtaef

More information

2.12 Pull Back, Push Forward and Lie Time Derivatives

2.12 Pull Back, Push Forward and Lie Time Derivatives Secton 2.2 2.2 Pull Bck Push Forwrd nd e me Dertes hs secton s n the mn concerned wth the follown ssue: n oserer ttched to fxed sy Crtesn coordnte system wll see mterl moe nd deform oer tme nd wll osere

More information

Chapter 3 Single Random Variables and Probability Distributions (Part 2)

Chapter 3 Single Random Variables and Probability Distributions (Part 2) Chpter 3 Single Rndom Vriles nd Proilit Distriutions (Prt ) Contents Wht is Rndom Vrile? Proilit Distriution Functions Cumultive Distriution Function Proilit Densit Function Common Rndom Vriles nd their

More information

Model Fitting and Robust Regression Methods

Model Fitting and Robust Regression Methods Dertment o Comuter Engneerng Unverst o Clorn t Snt Cruz Model Fttng nd Robust Regresson Methods CMPE 64: Imge Anlss nd Comuter Vson H o Fttng lnes nd ellses to mge dt Dertment o Comuter Engneerng Unverst

More information

Matrix Eigenvalues and Eigenvectors September 13, 2017

Matrix Eigenvalues and Eigenvectors September 13, 2017 Mtri Eigenvlues nd Eigenvectors September, 7 Mtri Eigenvlues nd Eigenvectors Lrry Cretto Mechnicl Engineering 5A Seminr in Engineering Anlysis September, 7 Outline Review lst lecture Definition of eigenvlues

More information

ME 501A Seminar in Engineering Analysis Page 1

ME 501A Seminar in Engineering Analysis Page 1 Mtr Trsformtos usg Egevectors September 8, Mtr Trsformtos Usg Egevectors Lrry Cretto Mechcl Egeerg A Semr Egeerg Alyss September 8, Outle Revew lst lecture Trsformtos wth mtr of egevectors: = - A ermt

More information

E-Companion: Mathematical Proofs

E-Companion: Mathematical Proofs E-omnon: Mthemtcl Poo Poo o emm : Pt DS Sytem y denton o t ey to vey tht t ncee n wth d ncee n We dene } ] : [ { M whee / We let the ttegy et o ech etle n DS e ]} [ ] [ : { M w whee M lge otve nume oth

More information

Rank One Update And the Google Matrix by Al Bernstein Signal Science, LLC

Rank 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 information

Name: SID: Discussion Session:

Name: SID: Discussion Session: Nme: SID: Dscusson Sesson: hemcl Engneerng hermodynmcs -- Fll 008 uesdy, Octoer, 008 Merm I - 70 mnutes 00 onts otl losed Book nd Notes (5 ponts). onsder n del gs wth constnt het cpctes. Indcte whether

More information

Introduction to Algebra - Part 2

Introduction to Algebra - Part 2 Alger Module A Introduction to Alger - Prt Copright This puliction The Northern Alert Institute of Technolog 00. All Rights Reserved. LAST REVISED Oct., 008 Introduction to Alger - Prt Sttement of Prerequisite

More information

Kinematics Quantities. Linear Motion. Coordinate System. Kinematics Quantities. Velocity. Position. Don t Forget Units!

Kinematics Quantities. Linear Motion. Coordinate System. Kinematics Quantities. Velocity. Position. Don t Forget Units! Knemtc Quntte Lner Phyc 11 Eyre Tme Intnt t Fundmentl Tme Interl t Dened Poton Fundmentl Dplcement Dened Aerge g Dened Aerge Accelerton g Dened Knemtc Quntte Scler: Mgntude Tme Intnt, Tme Interl nd Speed

More information

APPLIED HEAT-MASS-EXCHANGE PROBLEMS MODELING WITH NONLINEAR COEFFICIENTS IN HEAT AND GAS SUPPLY SYSTEM

APPLIED HEAT-MASS-EXCHANGE PROBLEMS MODELING WITH NONLINEAR COEFFICIENTS IN HEAT AND GAS SUPPLY SYSTEM Interntonl Jornl of Cvl Engneerng nd Technology IJCIET Volme 9 Isse 9 September 8 pp 8 88 Artcle ID: IJCIET_9_9_8 Avlble onlne t http://wwwemecom/jcet/ssessp?jtype=ijciet&vtype=9&itype=9 ISSN Prnt: 97-8

More information

ME 501A Seminar in Engineering Analysis Page 1

ME 501A Seminar in Engineering Analysis Page 1 Phse-plne Anlsis of Ordinr November, 7 Phse-plne Anlsis of Ordinr Lrr Cretto Mechnicl Engineering 5A Seminr in Engineering Anlsis November, 7 Outline Mierm exm two weeks from tonight covering ODEs nd Lplce

More information

ITERATIVE METHODS FOR SOLVING SYSTEMS OF LINEAR ALGEBRAIC EQUATIONS

ITERATIVE METHODS FOR SOLVING SYSTEMS OF LINEAR ALGEBRAIC EQUATIONS Numercl Alyss for Egeers Germ Jord Uversty ITERATIVE METHODS FOR SOLVING SYSTEMS OF LINEAR ALGEBRAIC EQUATIONS Numercl soluto of lrge systems of ler lgerc equtos usg drect methods such s Mtr Iverse, Guss

More information

REGULAR STURM-LIOUVILLE OPERATORS WITH TRANSMISSION CONDITIONS AT FINITE INTERIOR DISCONTINUOUS POINTS

REGULAR STURM-LIOUVILLE OPERATORS WITH TRANSMISSION CONDITIONS AT FINITE INTERIOR DISCONTINUOUS POINTS ournl o Mthetl enes: dvnes nd ppltons Volue Nuer Pes 65-77 REGULR TURM-LIOUVILLE OPERTOR WITH TRNMIION CONDITION T INITE INTERIOR DICONTINUOU POINT XIOLING HO nd IONG UN hool o Mthetl enes Inner Monol

More information

1 nonlinear.mcd Find solution root to nonlinear algebraic equation f(x)=0. Instructor: Nam Sun Wang

1 nonlinear.mcd Find solution root to nonlinear algebraic equation f(x)=0. Instructor: Nam Sun Wang nonlinermc Fin solution root to nonliner lgebric eqution ()= Instructor: Nm Sun Wng Bckgroun In science n engineering, we oten encounter lgebric equtions where we wnt to in root(s) tht stisies given eqution

More information

SOLVING SYSTEMS OF EQUATIONS, ITERATIVE METHODS

SOLVING SYSTEMS OF EQUATIONS, ITERATIVE METHODS ELM Numericl Anlysis Dr Muhrrem Mercimek SOLVING SYSTEMS OF EQUATIONS, ITERATIVE METHODS ELM Numericl Anlysis Some of the contents re dopted from Lurene V. Fusett, Applied Numericl Anlysis using MATLAB.

More information

Remember: Project Proposals are due April 11.

Remember: 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 information

Lecture 36. Finite Element Methods

Lecture 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 information

CISE 301: Numerical Methods Lecture 5, Topic 4 Least Squares, Curve Fitting

CISE 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 information

Answer: First, I ll show how to find the terms analytically then I ll show how to use the TI to find them.

Answer: First, I ll show how to find the terms analytically then I ll show how to use the TI to find them. . CHAPTER 0 SEQUENCE, SERIES, d INDUCTION Secto 0. Seqece A lst of mers specfc order. E / Fd the frst terms : of the gve seqece: Aswer: Frst, I ll show how to fd the terms ltcll the I ll show how to se

More information

We are looking for ways to compute the integral of a function f(x), f(x)dx.

We are looking for ways to compute the integral of a function f(x), f(x)dx. INTEGRATION TECHNIQUES Introdction We re looking for wys to compte the integrl of fnction f(x), f(x)dx. To pt it simply, wht we need to do is find fnction F (x) sch tht F (x) = f(x). Then if the integrl

More information

Announcements. Image Formation: Outline. The course. Image Formation and Cameras (cont.)

Announcements. Image Formation: Outline. The course. Image Formation and Cameras (cont.) nnouncements Imge Formton nd Cmers (cont.) ssgnment : Cmer & Lenses, gd Trnsformtons, nd Homogrph wll be posted lter tod. CSE 5 Lecture 5 CS5, Fll CS5, Fll CS5, Fll The course rt : The phscs of mgng rt

More information

Averaging Level Control of Multiple Tanks: A Passivity Based Approach

Averaging Level Control of Multiple Tanks: A Passivity Based Approach Avergng Level Control of Mltple Tnks: A Pssvty Bse Approch Dnel Sbrbro Deprtment of Electrcl Engneerng, Unvers e Concepcón, Concepcón,, Chle. In collborton wth Romeo Orteg Lbortore es Sgn et Systèmes mes,

More information

4. Eccentric axial loading, cross-section core

4. 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

4.5 JACOBI ITERATION FOR FINDING EIGENVALUES OF A REAL SYMMETRIC MATRIX. be a real symmetric matrix. ; (where we choose θ π for.

4.5 JACOBI ITERATION FOR FINDING EIGENVALUES OF A REAL SYMMETRIC MATRIX. be a real symmetric matrix. ; (where we choose θ π for. 4.5 JACOBI ITERATION FOR FINDING EIGENVALUES OF A REAL SYMMETRIC MATRIX Some reliminries: Let A be rel symmetric mtrix. Let Cos θ ; (where we choose θ π for Cos θ 4 purposes of convergence of the scheme)

More information

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

INTERPOLATION(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 information

First, we will find the components of the force of gravity: Perpendicular Forces (using away from the ramp as positive) ma F

First, we will find the components of the force of gravity: Perpendicular Forces (using away from the ramp as positive) ma F 1.. In Clss or Homework Eercise 1. An 18.0 kg bo is relesed on 33.0 o incline nd ccelertes t 0.300 m/s. Wht is the coeicient o riction? m 18.0kg 33.0? 0 0.300 m / s irst, we will ind the components o the

More information

5.2 Exponent Properties Involving Quotients

5.2 Exponent Properties Involving Quotients 5. Eponent Properties Involving Quotients Lerning Objectives Use the quotient of powers property. Use the power of quotient property. Simplify epressions involving quotient properties of eponents. Use

More information

Problem set 5: Solutions Math 207B, Winter r(x)u(x)v(x) dx.

Problem set 5: Solutions Math 207B, Winter r(x)u(x)v(x) dx. Problem set 5: Soltions Mth 7B, Winter 6. Sppose tht p : [, b] R is continosly differentible fnction sch tht p >, nd q, r : [, b] R re continos fnctions sch tht r >, q. Define weighted inner prodct on

More information

An Ising model on 2-D image

An Ising model on 2-D image School o Coputer Scence Approte Inerence: Loopy Bele Propgton nd vrnts Prolstc Grphcl Models 0-708 Lecture 4, ov 7, 007 Receptor A Knse C Gene G Receptor B Knse D Knse E 3 4 5 TF F 6 Gene H 7 8 Hetunndn

More information

A New Method for Solving Fuzzy Volterra Integro-Differential Equations

A New Method for Solving Fuzzy Volterra Integro-Differential Equations Astrln Jornl of Bsc n Apple Scences, 5(4: 54-64, 2 ISS 99-878 A ew Metho for Solvng Fzzy Volterr Integro-Dfferentl Eqtons T. Allhvrnloo, 2 A. Amrtemoor, M. Khezerloo, S. Khezerloo Deprtment of Mthemtcs,

More information

Lecture 22: Logic Synthesis (1)

Lecture 22: Logic Synthesis (1) Lecture 22: Logc Synthess (1) Sldes courtesy o Demng Chen Some sldes Courtesy o Pro. J. Cong o UCLA Outlne Redng Synthess nd optmzton o dgtl crcuts, G. De Mchel, 1994, Secton 2.5-2.5.1 Overvew Boolen lgebr

More information

I. Theory of Automata II. Theory of Formal Languages III. Theory of Turing Machines

I. Theory of Automata II. Theory of Formal Languages III. Theory of Turing Machines CI 3104 /Winter 2011: Introduction to Forml Lnguges Chpter 16: Non-Context-Free Lnguges Chpter 16: Non-Context-Free Lnguges I. Theory of utomt II. Theory of Forml Lnguges III. Theory of Turing Mchines

More information

Strong Gravity and the BKL Conjecture

Strong Gravity and the BKL Conjecture Introducton Strong Grvty nd the BKL Conecture Dvd Slon Penn Stte October 16, 2007 Dvd Slon Strong Grvty nd the BKL Conecture Introducton Outlne The BKL Conecture Ashtekr Vrbles Ksner Sngulrty 1 Introducton

More information

NAME: MR. WAIN FUNCTIONS

NAME: MR. WAIN FUNCTIONS NAME: M. WAIN FUNCTIONS evision o Solving Polnomil Equtions i one term in Emples Solve: 7 7 7 0 0 7 b.9 c 7 7 7 7 ii more thn one term in Method: Get the right hnd side to equl zero = 0 Eliminte ll denomintors

More information

BAR & TRUSS FINITE ELEMENT. Direct Stiffness Method

BAR & TRUSS FINITE ELEMENT. Direct Stiffness Method BAR & TRUSS FINITE ELEMENT Drect Stness Method FINITE ELEMENT ANALYSIS AND APPLICATIONS INTRODUCTION TO FINITE ELEMENT METHOD What s the nte element method (FEM)? A technqe or obtanng approxmate soltons

More information

Katholieke Universiteit Leuven Department of Computer Science

Katholieke 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 information

Two Coefficients of the Dyson Product

Two 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 information

1 Error Analysis of Simple Rules for Numerical Integration

1 Error Analysis of Simple Rules for Numerical Integration cs41: introduction to numericl nlysis 11/16/10 Lecture 19: Numericl Integrtion II Instructor: Professor Amos Ron Scries: Mrk Cowlishw, Nthnel Fillmore 1 Error Anlysis of Simple Rules for Numericl Integrtion

More information

Chapter 6 Continuous Random Variables and Distributions

Chapter 6 Continuous Random Variables and Distributions Chpter 6 Continuous Rndom Vriles nd Distriutions Mny economic nd usiness mesures such s sles investment consumption nd cost cn hve the continuous numericl vlues so tht they cn not e represented y discrete

More information

1.B Appendix to Chapter 1

1.B Appendix to Chapter 1 Secon.B.B Append o Chper.B. The Ordnr Clcl Here re led ome mporn concep rom he ordnr clcl. The Dervve Conder ncon o one ndependen vrble. The dervve o dened b d d lm lm.b. where he ncremen n de o n ncremen

More information

Chapter 6. Operational Amplifier. inputs can be defined as the average of the sum of the two signals.

Chapter 6. Operational Amplifier.  inputs can be defined as the average of the sum of the two signals. 6 Operatonal mpler Chapter 6 Operatonal mpler CC Symbol: nput nput Output EE () Non-nvertng termnal, () nvertng termnal nput mpedance : Few mega (ery hgh), Output mpedance : Less than (ery low) Derental

More information

ES.182A Topic 32 Notes Jeremy Orloff

ES.182A Topic 32 Notes Jeremy Orloff ES.8A Topic 3 Notes Jerem Orloff 3 Polr coordintes nd double integrls 3. Polr Coordintes (, ) = (r cos(θ), r sin(θ)) r θ Stndrd,, r, θ tringle Polr coordintes re just stndrd trigonometric reltions. In

More information

Lecture 16. Chapter 11. Energy Dissipation Linear Momentum. Physics I. Department of Physics and Applied Physics

Lecture 16. Chapter 11. Energy Dissipation Linear Momentum. Physics I. Department of Physics and Applied Physics Lecture 16 Chapter 11 Physcs I Energy Dsspaton Lnear Momentum Course webste: http://aculty.uml.edu/andry_danylov/teachng/physcsi Department o Physcs and Appled Physcs IN IN THIS CHAPTER, you wll learn

More information

1. Extend QR downwards to meet the x-axis at U(6, 0). y

1. Extend QR downwards to meet the x-axis at U(6, 0). y In the digrm, two stright lines re to be drwn through so tht the lines divide the figure OPQRST into pieces of equl re Find the sum of the slopes of the lines R(6, ) S(, ) T(, 0) Determine ll liner functions

More information

Effects of polarization on the reflected wave

Effects of polarization on the reflected wave Lecture Notes. L Ros PPLIED OPTICS Effects of polrzton on the reflected wve Ref: The Feynmn Lectures on Physcs, Vol-I, Secton 33-6 Plne of ncdence Z Plne of nterfce Fg. 1 Y Y r 1 Glss r 1 Glss Fg. Reflecton

More information

A Modeling System to Combine Optimization and Constraint. Programming. INFORMS, November Carnegie Mellon University.

A Modeling System to Combine Optimization and Constraint. Programming. INFORMS, November Carnegie Mellon University. A Modelng Sstem to Combne Optmzaton and Constrant Programmng John Hooker Carnege Mellon Unverst INFORMS November 000 Based on ont work wth Ignaco Grossmann Hak-Jn Km Mara Axlo Osoro Greger Ottosson Erlendr

More information

Math 61CM - Solutions to homework 9

Math 61CM - Solutions to homework 9 Mth 61CM - Solutions to homework 9 Cédric De Groote November 30 th, 2018 Problem 1: Recll tht the left limit of function f t point c is defined s follows: lim f(x) = l x c if for ny > 0 there exists δ

More information

Let us look at a linear equation for a one-port network, for example some load with a reflection coefficient s, Figure L6.

Let us look at a linear equation for a one-port network, for example some load with a reflection coefficient s, Figure L6. ECEN 5004, prng 08 Actve Mcrowve Crcut Zoy Popovc, Unverty of Colordo, Boulder LECURE 5 IGNAL FLOW GRAPH FOR MICROWAVE CIRCUI ANALYI In mny text on mcrowve mplfer (e.g. the clc one by Gonzlez), gnl flow-grph

More information

Calculus 2: Integration. Differentiation. Integration

Calculus 2: Integration. Differentiation. Integration Clculus 2: Integrtion The reverse process to differentition is known s integrtion. Differentition f() f () Integrtion As it is the opposite of finding the derivtive, the function obtined b integrtion is

More information

Lecture 2: January 27

Lecture 2: January 27 CS 684: Algorithmic Gme Theory Spring 217 Lecturer: Év Trdos Lecture 2: Jnury 27 Scrie: Alert Julius Liu 2.1 Logistics Scrie notes must e sumitted within 24 hours of the corresponding lecture for full

More information

Work and Energy (Work Done by a Varying Force)

Work and Energy (Work Done by a Varying Force) Lecture 1 Chpter 7 Physcs I 3.5.14 ork nd Energy (ork Done y Vryng Force) Course weste: http://fculty.uml.edu/andry_dnylov/techng/physcsi Lecture Cpture: http://echo36.uml.edu/dnylov13/physcs1fll.html

More information

Lecture 26: Angular Momentum II

Lecture 26: Angular Momentum II ecture 6: Angulr Momentum II Ph851 Fll 9 The Angulr Momentum Opertor The ngulr momentum opertor is defined s: r r R r P It is vector opertor: r r e r e r e According to the definition of the crossproduct

More information

Goals: Determine how to calculate the area described by a function. Define the definite integral. Explore the relationship between the definite

Goals: Determine how to calculate the area described by a function. Define the definite integral. Explore the relationship between the definite Unit #8 : The Integrl Gols: Determine how to clculte the re described by function. Define the definite integrl. Eplore the reltionship between the definite integrl nd re. Eplore wys to estimte the definite

More information

Experimental Errors and Error Analysis

Experimental Errors and Error Analysis Expermental Errors and Error Analss Rajee Prabhakar Unerst of Texas at Astn, Astn, TX Freeman Research Grop Meetng March 10, 004 Topcs coered Tpes of expermental errors Redcng errors Descrbng errors qanttatel

More information

Electricity and Magnetism Electric Dipole Continuous Distribution of Charge

Electricity and Magnetism Electric Dipole Continuous Distribution of Charge Electricit nd Mgnetism Electric Dipole Continos Distribtion of Chrge Ln heridn De Anz College Jn 16, 2018 Lst time electric field lines electric field from point chrge net electric field from mn chrges

More information

C PLANE ELASTICITY PROBLEM FORMULATIONS

C PLANE ELASTICITY PROBLEM FORMULATIONS C M.. Tamn, CSMLab, UTM Corse Content: A ITRODUCTIO AD OVERVIEW mercal method and Compter-Aded Engneerng; Phscal problems; Mathematcal models; Fnte element method. B REVIEW OF -D FORMULATIOS Elements and

More information

5.7 Improper Integrals

5.7 Improper Integrals 458 pplictions of definite integrls 5.7 Improper Integrls In Section 5.4, we computed the work required to lift pylod of mss m from the surfce of moon of mss nd rdius R to height H bove the surfce of the

More information

Chapter 1. Chapter 1 1

Chapter 1. Chapter 1 1 Chpter Chpter : Signls nd Systems... 2. Introduction... 2.2 Signls... 3.2. Smpling... 4.2.2 Periodic Signls... 0.2.3 Discrete-Time Sinusoidl Signls... 2.2.4 Rel Exponentil Signls... 5.2.5 Complex Exponentil

More information

Mathematical Derivation of Angular Momenta in Quantum Physics

Mathematical Derivation of Angular Momenta in Quantum Physics ornl of Modern Physcs,,, 9-99 http://d.do.org/.6/jmp..75 Pblshed Onlne ly (http://www.scrp.org/jornl/jmp) Mthemtcl Dervton of Anglr Moment n Qntm Physcs Dnel Grcker Unversté de Strsborg-CNRS, Strsborg,

More information

The Quest for Perfect and Compact Symmetry Breaking for Graph Problems

The Quest for Perfect and Compact Symmetry Breaking for Graph Problems The Quest for Perfect nd Compct Symmetry Breking for Grph Prolems Mrijn J.H. Heule SYNASC Septemer 25, 2016 1/19 Stisfiility (SAT) solving hs mny pplictions... forml verifiction grph theory ioinformtics

More information