36.1 Why is it important to be able to find roots to systems of equations? Up to this point, we have discussed how to find the solution to
|
|
- Bruce Potter
- 6 years ago
- Views:
Transcription
1 ChE Lecture Notes - D. Keer, 5/9/98 Lecture 6,7,8 - Rootndng n systems o equatons (A) Theory (B) Problems (C) MATLAB Applcatons Tet: Supplementary notes rom Instructor 6. Why s t mportant to be able to nd roots to systems o equatons? Up to ths pont, we have dscussed how to nd the soluton to sngle non-lnear equaton (Newton Raphson Method) systems o lnear equatons (Lnear Algebra) However, requently n chemcal engneerng, we need to nd the soluton to a system o non-lnear equatons. Ths orces us to combne the Newton-Raphson method wth Lnear Algebra. The technque we wll dscuss n ths secton s Multvarate Newton Raphson Method. 6. Multvarate Newton-Raphson - Theory Recall that when we wanted to nd the soluton to a sngle non-lnear equaton o the orm () (.) The bass o the Newton-Raphson method lay n the act that we can appromate the dervatve o () numercally. ( ) d() d () ( ) (.) leadng to the teratve ormula () ( ) (.5) Now consder the system o n non-lnear equatons and n unknowns.
2 ChE Lecture Notes - D. Keer, 5/9/98... n n (,,,...n, n) (,,,...n, n) (,,,..., ) (,,,...n, n) (,,,..., ) n n n n (6.) Ths s the most general orm o the problems that ace chemcal engneers and encompass lnear equatons (when all the unctons,, are lnear) and sngle equatons (when n ). One technque used to solve ths problem s called the multvarate Newton Raphson Method (MNRM). The dea ollows rom the sngle-varable case. The basc dea agan stems rom the act that the total dervatve o a uncton,, s d d d (6.) or the case o two varables or d n d (6.) or the n varable case. We can dscretze ths epresson and wrte: () () ({ } ) { } n () () ( ) ( ) (6.4) where the s the nde over unctons, the s the nde over varables and the superscrpt n parentheses stands or the teraton. Note that we have only one varable (n), ths reduces to the Newton-Raphson method that we have already learned. Now consder that we have n equatons so that to n. We have a system o equatons whch we can wrte as: () ({ } ) { } n () () () ( ) ( ) (6.5)
3 ChE Lecture Notes - D. Keer, 5/9/98 Just as n the sngle varable case, we want our net teraton to take us to the root so we assume (). We can then wrte ths system n matr notaton as: ( ) that { } J δ R (6.6) wherer s called the resdual vector at the k th teraton and s dened as ( k) R ({ } ) (6.7) wherej s called the Jacoban matr at the k th teratonand s dened as and where ( J ), (6.8) δ (k ) (6.9) so that the new guess or the s (k ) δ (6.) The algorthm or solvng the multvarate Newton-Raphson ollows analogously rom the sngle varable NRM. The steps are as ollows:. Make an ntal guess or. calculate the Jacoban and the Resdual.. Solve equaton Calculate new rom equaton I the soluton has not converged, loop back to step. The multvarate Newton-Raphson Method suers rom the same short-comngs as the sngle-varable Newton-Raphson Method. () You need a good ntal guess. () You don t get quadratc convergence untl you are close to the soluton. () I the partal dervatves are zero, the method blows up. I the partal dervatves are close to zero, the method may not converge.
4 ChE Lecture Notes - D. Keer, 5/9/98 6. Multvarate Newton-Raphson - Problems As wth all the numercal technques we have been usng, t s necessary to practce wth them beore we become procent wth them. Eample One: ( (,, ) ) ( ) ( ) 4 ( ) ( ) (6.) The vs plot o the ollowng unctons s gven below. The rst uncton s that o a crcle, centered at the orgn wth radus. The second uncton s a parabola. We see there are two and only two solutons to ths systems o equatons
5 ChE Lecture Notes - D. Keer, 5/9/98 In order to use MNRM, we must rst determne the unctonal orm o the partal dervatves ( J) ( J), ( J) ( J),,, Then ollowng the algorthm outlned above: Step One. Make an ntal guess. One o the solutons looks to be at. and.. Step Two. Usng that ntal guess, calculate the resdual and the Jacoban. ( 4 and R () J ) Step Three. Solve J δ R (Usng Lnear Algebra) () δ.. Step Four. Calculate new values or va equaton (6.) ().9.8 Step Fve. Loop back to Step. and repeat untl converged. Here are what urther teratons yeld teraton J R δ
6 ChE Lecture Notes - D. Keer, 5/9/ e.79e.559e e e e e e - 8 The other root n to the problem s located at by symmetry. Eample Two. Lnear systems are a subset o non-lnear systems. The multvarate Newton-Raphson solve lnear systems eactly n one teraton, ust as was the case n the sngle-varable problem. Consder the system o lnear equatons: ( (,, ) 5 ) ( ) ( ) 4 ( ) ( ) (6.)
7 ChE Lecture Notes - D. Keer, 5/9/98 In order to use MNRM, we must rst determne the unctonal orm o the partal dervatves ( J) 5 ( J), ( J) ( J),,, Then ollowng the algorthm outlned above: Step One. Make an ntal guess. One o the solutons looks to be at. and.. Step Two. Usng that ntal guess, calculate the resdual and the Jacoban. ( 5 8 and R () J ) Step Three. Solve J δ R (Usng Lnear Algebra) () δ Step Four. Calculate new values or va equaton (6.) () A second teraton wll show that ths s the eact soluton. 6. Multvarate Newton-Raphson - MATLAB Very quckly problems lke ths become too dcult to solve by hand. We rely on sotware lke MATLAB to solve these problems or us. MATLAB does not necessarly use a multvarate Newton-Raphson method to solve a system o non-lnear equatons but t uses some smlar numercal technque. Ths alternate technque appromates the partal dervatves, so that the user only has to enter the unctonal orm o the unctons and not that o the partal dervatves. 7
8 ChE Lecture Notes - D. Keer, 5/9/98 On the webste, you can download a routne called syseqn.m. Ths routne wll allow you to solve a system o non-lnear algebrac equatons. The descrpton or how to use the le can be obtaned by openng MATLAB, movng to the drectory where you have downloaded the syseqn.m le, and typng help syseqn Ths yelds: syseqn(o) syseqn solves a system o nonlnear algebrac equatons. The system o equatons are stored n the le syseqnnput.m The ntal guesses are gven as a vector n o. The soluton s wrtten to the screen and to the le "syseqn.out" Author: Davd Keer Date: October, 998 The routne uses the MATLAB nstrnsc uncton zero to solve the system o equatons. At the MATLAB command lne, the routne s nvoke by typng, or eample, syseqn() you had one equaton and you wanted your ntal guess to be o. Or, you had a system o three non-lnear equaton, you could nvoke syseqn.m by typng syseqn([;4;6]) you had one equaton and you wanted your ntal guess to be, o, 6. o,, 4 o, The algebrac equatons are entered n the le syseqnnput.m Ths le can be as smple as uncton [] syseqnnput() ^-; to solve or the roots o -. The syseqnnput.m can also be much more complcated. Two more nvolved eamples o the syseqnnput.m or (a) solvng a system o mass balances n an etractor, and (b) determnng the chemcal equlbra compostons and temperature n a vessel where multple reactons are occurng n a non-adabatc and non-sothermal envronment. In the latter case, you wll see that the syseqnnput.m le can be very complcated. It doesn t matter how complcated t appears, so long as, at the end o the le, you have evaluated the uncton o nterest. Eample: Consder the ollowng problem rom chemcal engneerng. You want to use an etracton process to remove a contamnant rom a eed stream. The process dagram looks lke ths: 8
9 ChE Lecture Notes - D. Keer, 5/9/98 Etract, E Solvent, S mass transer Feed, F Ranate, R The data you are gven s F mol / hr F,b F,c F,..9. S S,b S, 5mol / hr S,c R? R,??. E? E,c E,?.? You have s unknown varables. You need s equatons to solve these. You have three mass balances: F S R E F S F,c R E S, E, You have two constrants on the compostons E,c You have one separaton rato: K.68 9
10 ChE Lecture Notes - D. Keer, 5/9/98 These are your s equatons. The second, thrd, and sth equatons are non-lnear. You need to use a technque lke multvarate Newton-Raphson to solve ths problem. In order to use MNRM to solve ths system o non-lnear equatons, rearrange the equatons so that the rght hand sde s zero F F S E,c S F,c S,.68 R R E E E, MATLAB can solve ths problem easly, usng the syseqn.m routne provded on the web-ste. You don t alter the syseqn.m le. You smply type the equatons nto the syseqnnput.m le and then nvoke the syseqn.m routne wth reasonable ntal guesses. Ths routne can then be used to do parameter studes.. Look at derent compostons n the ranate as a uncton o derent eed rates or solvent rates or eed or solvent compostons. (In order to change parameters, lke F, you only have to change the value o F n etractnput.m.). Look at derent cases. Perhaps you want to specy the Ranate composton and determne the eedrate. (In order to change whch unctons are varables, you smply make sure that you assgn your ntal guesses to the approprate varables at the begnnng o the routne.. Look at any mass balance system. Ths code gves the skeleton or solvng any system o nonlnear equatons. All that s requred s changng the nput le. Wth a lttle study o the program, you wll be able to mody t and use t n all your classes. Usng ths code solves the above eample n a racton o a second: F mol / hr F,b F,c F,..9. S S,b S, 5mol / hr S,c R 94.74mol / hr R, E 55.6mol / hr E, E,c
OPTIMISATION. Introduction Single Variable Unconstrained Optimisation Multivariable Unconstrained Optimisation Linear Programming
OPTIMIATION Introducton ngle Varable Unconstraned Optmsaton Multvarable Unconstraned Optmsaton Lnear Programmng Chapter Optmsaton /. Introducton In an engneerng analss, sometmes etremtes, ether mnmum or
More informationChapter 3 Differentiation and Integration
MEE07 Computer Modelng Technques n Engneerng Chapter Derentaton and Integraton Reerence: An Introducton to Numercal Computatons, nd edton, S. yakowtz and F. zdarovsky, Mawell/Macmllan, 990. Derentaton
More informationSingle Variable Optimization
8/4/07 Course Instructor Dr. Raymond C. Rump Oce: A 337 Phone: (95) 747 6958 E Mal: rcrump@utep.edu Topc 8b Sngle Varable Optmzaton EE 4386/530 Computatonal Methods n EE Outlne Mathematcal Prelmnares Sngle
More informationLecture 2 Solution of Nonlinear Equations ( Root Finding Problems )
Lecture Soluton o Nonlnear Equatons Root Fndng Problems Dentons Classcaton o Methods Analytcal Solutons Graphcal Methods Numercal Methods Bracketng Methods Open Methods Convergence Notatons Root Fndng
More information: Numerical Analysis Topic 2: Solution of Nonlinear Equations Lectures 5-11:
764: Numercal Analyss Topc : Soluton o Nonlnear Equatons Lectures 5-: UIN Malang Read Chapters 5 and 6 o the tetbook 764_Topc Lecture 5 Soluton o Nonlnear Equatons Root Fndng Problems Dentons Classcaton
More informationSummary with Examples for Root finding Methods -Bisection -Newton Raphson -Secant
Summary wth Eamples or Root ndng Methods -Bsecton -Newton Raphson -Secant Nonlnear Equaton Solvers Bracketng Graphcal Open Methods Bsecton False Poston (Regula-Fals) Newton Raphson Secant All Iteratve
More informationCHAPTER 4d. ROOTS OF EQUATIONS
CHAPTER 4d. ROOTS OF EQUATIONS A. J. Clark School o Engneerng Department o Cvl and Envronmental Engneerng by Dr. Ibrahm A. Assakka Sprng 00 ENCE 03 - Computaton Methods n Cvl Engneerng II Department o
More informationEEE 241: Linear Systems
EEE : Lnear Systems Summary #: Backpropagaton BACKPROPAGATION The perceptron rule as well as the Wdrow Hoff learnng were desgned to tran sngle layer networks. They suffer from the same dsadvantage: they
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 informationCISE301: Numerical Methods Topic 2: Solution of Nonlinear Equations
CISE3: Numercal Methods Topc : Soluton o Nonlnear Equatons Dr. Amar Khoukh Term Read Chapters 5 and 6 o the tetbook CISE3_Topc c Khoukh_ Lecture 5 Soluton o Nonlnear Equatons Root ndng Problems Dentons
More informationChapter Newton s Method
Chapter 9. Newton s Method After readng ths chapter, you should be able to:. Understand how Newton s method s dfferent from the Golden Secton Search method. Understand how Newton s method works 3. Solve
More informationFinite 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 informationSolution of Linear System of Equations and Matrix Inversion Gauss Seidel Iteration Method
Soluton of Lnear System of Equatons and Matr Inverson Gauss Sedel Iteraton Method It s another well-known teratve method for solvng a system of lnear equatons of the form a + a22 + + ann = b a2 + a222
More informationLinear Feature Engineering 11
Lnear Feature Engneerng 11 2 Least-Squares 2.1 Smple least-squares Consder the followng dataset. We have a bunch of nputs x and correspondng outputs y. The partcular values n ths dataset are x y 0.23 0.19
More informationNumerical Methods Solution of Nonlinear Equations
umercal Methods Soluton o onlnear Equatons Lecture Soluton o onlnear Equatons Root Fndng Prolems Dentons Classcaton o Methods Analytcal Solutons Graphcal Methods umercal Methods Bracketng Methods Open
More informationSection 8.3 Polar Form of Complex Numbers
80 Chapter 8 Secton 8 Polar Form of Complex Numbers From prevous classes, you may have encountered magnary numbers the square roots of negatve numbers and, more generally, complex numbers whch are the
More informationLecture 21: Numerical methods for pricing American type derivatives
Lecture 21: Numercal methods for prcng Amercan type dervatves Xaoguang Wang STAT 598W Aprl 10th, 2014 (STAT 598W) Lecture 21 1 / 26 Outlne 1 Fnte Dfference Method Explct Method Penalty Method (STAT 598W)
More informationME 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 informationCinChE Problem-Solving Strategy Chapter 4 Development of a Mathematical Model. formulation. procedure
nhe roblem-solvng Strategy hapter 4 Transformaton rocess onceptual Model formulaton procedure Mathematcal Model The mathematcal model s an abstracton that represents the engneerng phenomena occurrng n
More information1 GSW Iterative Techniques for y = Ax
1 for y = A I m gong to cheat here. here are a lot of teratve technques that can be used to solve the general case of a set of smultaneous equatons (wrtten n the matr form as y = A), but ths chapter sn
More informationA MODIFIED METHOD FOR SOLVING SYSTEM OF NONLINEAR EQUATIONS
Journal of Mathematcs and Statstcs 9 (1): 4-8, 1 ISSN 1549-644 1 Scence Publcatons do:1.844/jmssp.1.4.8 Publshed Onlne 9 (1) 1 (http://www.thescpub.com/jmss.toc) A MODIFIED METHOD FOR SOLVING SYSTEM OF
More informationChapter 12. Ordinary Differential Equation Boundary Value (BV) Problems
Chapter. Ordnar Dfferental Equaton Boundar Value (BV) Problems In ths chapter we wll learn how to solve ODE boundar value problem. BV ODE s usuall gven wth x beng the ndependent space varable. p( x) q(
More informationUniversity of Washington Department of Chemistry Chemistry 452/456 Summer Quarter 2014
Lecture 12 7/25/14 ERD: 7.1-7.5 Devoe: 8.1.1-8.1.2, 8.2.1-8.2.3, 8.4.1-8.4.3 Unversty o Washngton Department o Chemstry Chemstry 452/456 Summer Quarter 2014 A. Free Energy and Changes n Composton: The
More informationA Tale of Friction Basic Rollercoaster Physics. Fahrenheit Rollercoaster, Hershey, PA max height = 121 ft max speed = 58 mph
A Tale o Frcton Basc Rollercoaster Physcs Fahrenhet Rollercoaster, Hershey, PA max heght = 11 t max speed = 58 mph PLAY PLAY PLAY PLAY Rotatonal Movement Knematcs Smlar to how lnear velocty s dened, angular
More informationCSci 6974 and ECSE 6966 Math. Tech. for Vision, Graphics and Robotics Lecture 21, April 17, 2006 Estimating A Plane Homography
CSc 6974 and ECSE 6966 Math. Tech. for Vson, Graphcs and Robotcs Lecture 21, Aprl 17, 2006 Estmatng A Plane Homography Overvew We contnue wth a dscusson of the major ssues, usng estmaton of plane projectve
More informationChapter 5 rd Law of Thermodynamics
Entropy and the nd and 3 rd Chapter 5 rd Law o hermodynamcs homas Engel, hlp Red Objectves Introduce entropy. Derve the condtons or spontanety. Show how S vares wth the macroscopc varables,, and. Chapter
More information2 Finite difference basics
Numersche Methoden 1, WS 11/12 B.J.P. Kaus 2 Fnte dfference bascs Consder the one- The bascs of the fnte dfference method are best understood wth an example. dmensonal transent heat conducton equaton T
More informationReview of Taylor Series. Read Section 1.2
Revew of Taylor Seres Read Secton 1.2 1 Power Seres A power seres about c s an nfnte seres of the form k = 0 k a ( x c) = a + a ( x c) + a ( x c) + a ( x c) k 2 3 0 1 2 3 + In many cases, c = 0, and the
More informationWork is the change in energy of a system (neglecting heat transfer). To examine what could
Work Work s the change n energy o a system (neglectng heat transer). To eamne what could cause work, let s look at the dmensons o energy: L ML E M L F L so T T dmensonally energy s equal to a orce tmes
More informationIntegrals and Invariants of Euler-Lagrange Equations
Lecture 16 Integrals and Invarants of Euler-Lagrange Equatons ME 256 at the Indan Insttute of Scence, Bengaluru Varatonal Methods and Structural Optmzaton G. K. Ananthasuresh Professor, Mechancal Engneerng,
More informationEE 330 Lecture 24. Small Signal Analysis Small Signal Analysis of BJT Amplifier
EE 0 Lecture 4 Small Sgnal Analss Small Sgnal Analss o BJT Ampler Eam Frda March 9 Eam Frda Aprl Revew Sesson or Eam : 6:00 p.m. on Thursda March 8 n Room Sweene 6 Revew rom Last Lecture Comparson o Gans
More informationRoot Finding
Root Fndng 886307 What s Computer Scence? Computer scence s a dscplne that spans theory and practce. It requres thnkng both n abstract terms and n concrete terms. The practcal sde o computng can be seen
More informationUniversity of Washington Department of Chemistry Chemistry 452/456 Summer Quarter 2014
Lecture 16 8/4/14 Unversty o Washngton Department o Chemstry Chemstry 452/456 Summer Quarter 214. Real Vapors and Fugacty Henry s Law accounts or the propertes o extremely dlute soluton. s shown n Fgure
More informationTEST 5 (phy 240) 2. Show that the volume coefficient of thermal expansion for an ideal gas at constant pressure is temperature dependent and given by
ES 5 (phy 40). a) Wrte the zeroth law o thermodynamcs. b) What s thermal conductvty? c) Identyng all es, draw schematcally a P dagram o the arnot cycle. d) What s the ecency o an engne and what s the coecent
More informationSome Comments on Accelerating Convergence of Iterative Sequences Using Direct Inversion of the Iterative Subspace (DIIS)
Some Comments on Acceleratng Convergence of Iteratve Sequences Usng Drect Inverson of the Iteratve Subspace (DIIS) C. Davd Sherrll School of Chemstry and Bochemstry Georga Insttute of Technology May 1998
More informationGeneral Tips on How to Do Well in Physics Exams. 1. Establish a good habit in keeping track of your steps. For example, when you use the equation
General Tps on How to Do Well n Physcs Exams 1. Establsh a good habt n keepng track o your steps. For example when you use the equaton 1 1 1 + = d d to solve or d o you should rst rewrte t as 1 1 1 = d
More informationInexact Newton Methods for Inverse Eigenvalue Problems
Inexact Newton Methods for Inverse Egenvalue Problems Zheng-jan Ba Abstract In ths paper, we survey some of the latest development n usng nexact Newton-lke methods for solvng nverse egenvalue problems.
More informationCollege of Computer & Information Science Fall 2009 Northeastern University 20 October 2009
College of Computer & Informaton Scence Fall 2009 Northeastern Unversty 20 October 2009 CS7880: Algorthmc Power Tools Scrbe: Jan Wen and Laura Poplawsk Lecture Outlne: Prmal-dual schema Network Desgn:
More informationChapter 5 Function-based Monte Carlo
Chapter 5 Functon-based Monte Carlo 5.1 Four technques or estmatng ntegrals Our net set o mathematcal tools that we wll develop nvolve Monte Carlo ntegraton. In the grand scheme o thngs, our study so ar
More informationGrid Generation around a Cylinder by Complex Potential Functions
Research Journal of Appled Scences, Engneerng and Technolog 4(): 53-535, 0 ISSN: 040-7467 Mawell Scentfc Organzaton, 0 Submtted: December 0, 0 Accepted: Januar, 0 Publshed: June 0, 0 Grd Generaton around
More informationFE 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 informationSection 3.6 Complex Zeros
04 Chapter Secton 6 Comple Zeros When fndng the zeros of polynomals, at some pont you're faced wth the problem Whle there are clearly no real numbers that are solutons to ths equaton, leavng thngs there
More informationComplex Variables. Chapter 18 Integration in the Complex Plane. March 12, 2013 Lecturer: Shih-Yuan Chen
omplex Varables hapter 8 Integraton n the omplex Plane March, Lecturer: Shh-Yuan hen Except where otherwse noted, content s lcensed under a BY-N-SA. TW Lcense. ontents ontour ntegrals auchy-goursat theorem
More information10.34 Numerical Methods Applied to Chemical Engineering Fall Homework #3: Systems of Nonlinear Equations and Optimization
10.34 Numercal Methods Appled to Chemcal Engneerng Fall 2015 Homework #3: Systems of Nonlnear Equatons and Optmzaton Problem 1 (30 ponts). A (homogeneous) azeotrope s a composton of a multcomponent mxture
More information15-381: Artificial Intelligence. Regression and cross validation
15-381: Artfcal Intellgence Regresson and cross valdaton Where e are Inputs Densty Estmator Probablty Inputs Classfer Predct category Inputs Regressor Predct real no. Today Lnear regresson Gven an nput
More informationSpring Force and Power
Lecture 13 Chapter 9 Sprng Force and Power Yeah, energy s better than orces. What s net? Course webste: http://aculty.uml.edu/andry_danylov/teachng/physcsi IN THIS CHAPTER, you wll learn how to solve problems
More informationForce = F Piston area = A
CHAPTER III Ths chapter s an mportant transton between the propertes o pure substances and the most mportant chapter whch s: the rst law o thermodynamcs In ths chapter, we wll ntroduce the notons o heat,
More information(1) The saturation vapor pressure as a function of temperature, often given by the Antoine equation:
CE304, Sprng 2004 Lecture 22 Lecture 22: Topcs n Phase Equlbra, part : For the remander of the course, we wll return to the subject of vapor/lqud equlbrum and ntroduce other phase equlbrum calculatons
More informationUnit 5: Quadratic Equations & Functions
Date Perod Unt 5: Quadratc Equatons & Functons DAY TOPIC 1 Modelng Data wth Quadratc Functons Factorng Quadratc Epressons 3 Solvng Quadratc Equatons 4 Comple Numbers Smplfcaton, Addton/Subtracton & Multplcaton
More informationMathematical Economics MEMF e ME. Filomena Garcia. Topic 2 Calculus
Mathematcal Economcs MEMF e ME Flomena Garca Topc 2 Calculus Mathematcal Economcs - www.seg.utl.pt/~garca/economa_matematca . Unvarate Calculus Calculus Functons : X Y y ( gves or each element X one element
More informationECEN 5005 Crystals, Nanocrystals and Device Applications Class 19 Group Theory For Crystals
ECEN 5005 Crystals, Nanocrystals and Devce Applcatons Class 9 Group Theory For Crystals Dee Dagram Radatve Transton Probablty Wgner-Ecart Theorem Selecton Rule Dee Dagram Expermentally determned energy
More informationDifference Equations
Dfference Equatons c Jan Vrbk 1 Bascs Suppose a sequence of numbers, say a 0,a 1,a,a 3,... s defned by a certan general relatonshp between, say, three consecutve values of the sequence, e.g. a + +3a +1
More informationLecture 12: Discrete Laplacian
Lecture 12: Dscrete Laplacan Scrbe: Tanye Lu Our goal s to come up wth a dscrete verson of Laplacan operator for trangulated surfaces, so that we can use t n practce to solve related problems We are mostly
More informationORDINARY DIFFERENTIAL EQUATIONS EULER S METHOD
Numercal Analss or Engneers German Jordanan Unverst ORDINARY DIFFERENTIAL EQUATIONS We wll eplore several metods o solvng rst order ordnar derental equatons (ODEs and we wll sow ow tese metods can be appled
More informationRelaxation Methods for Iterative Solution to Linear Systems of Equations
Relaxaton Methods for Iteratve Soluton to Lnear Systems of Equatons Gerald Recktenwald Portland State Unversty Mechancal Engneerng Department gerry@pdx.edu Overvew Techncal topcs Basc Concepts Statonary
More informationSome modelling aspects for the Matlab implementation of MMA
Some modellng aspects for the Matlab mplementaton of MMA Krster Svanberg krlle@math.kth.se Optmzaton and Systems Theory Department of Mathematcs KTH, SE 10044 Stockholm September 2004 1. Consdered optmzaton
More informationSOLVING NON-LINEAR SYSTEMS BY NEWTON s METHOD USING SPREADSHEET EXCEL Tay Kim Gaik Universiti Tun Hussein Onn Malaysia
SOLVING NON-LINEAR SYSTEMS BY NEWTON s METHOD USING SPREADSHEET EXCEL Tay Km Gak Unverst Tun Hussen Onn Malaysa Kek Se Long Unverst Tun Hussen Onn Malaysa Rosmla Abdul-Kahar
More informationMMA and GCMMA two methods for nonlinear optimization
MMA and GCMMA two methods for nonlnear optmzaton Krster Svanberg Optmzaton and Systems Theory, KTH, Stockholm, Sweden. krlle@math.kth.se Ths note descrbes the algorthms used n the author s 2007 mplementatons
More informationCHARACTERISTICS OF COMPLEX SEPARATION SCHEMES AND AN ERROR OF SEPARATION PRODUCTS OUTPUT DETERMINATION
Górnctwo Geonżynera Rok 0 Zeszyt / 006 Igor Konstantnovch Mladetskj * Petr Ivanovch Plov * Ekaterna Nkolaevna Kobets * Tasya Igorevna Markova * CHARACTERISTICS OF COMPLEX SEPARATION SCHEMES AND AN ERROR
More informationCHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE
CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE Analytcal soluton s usually not possble when exctaton vares arbtrarly wth tme or f the system s nonlnear. Such problems can be solved by numercal tmesteppng
More informationAdvanced Circuits Topics - Part 1 by Dr. Colton (Fall 2017)
Advanced rcuts Topcs - Part by Dr. olton (Fall 07) Part : Some thngs you should already know from Physcs 0 and 45 These are all thngs that you should have learned n Physcs 0 and/or 45. Ths secton s organzed
More informationU.C. Berkeley CS294: Beyond Worst-Case Analysis Luca Trevisan September 5, 2017
U.C. Berkeley CS94: Beyond Worst-Case Analyss Handout 4s Luca Trevsan September 5, 07 Summary of Lecture 4 In whch we ntroduce semdefnte programmng and apply t to Max Cut. Semdefnte Programmng Recall that
More informationThe Fundamental Theorem of Algebra. Objective To use the Fundamental Theorem of Algebra to solve polynomial equations with complex solutions
5-6 The Fundamental Theorem of Algebra Content Standards N.CN.7 Solve quadratc equatons wth real coeffcents that have comple solutons. N.CN.8 Etend polnomal denttes to the comple numbers. Also N.CN.9,
More informationPHYS 1441 Section 002 Lecture #15
PHYS 1441 Secton 00 Lecture #15 Monday, March 18, 013 Work wth rcton Potental Energy Gravtatonal Potental Energy Elastc Potental Energy Mechancal Energy Conservaton Announcements Mdterm comprehensve exam
More informationLecture 10 Support Vector Machines II
Lecture 10 Support Vector Machnes II 22 February 2016 Taylor B. Arnold Yale Statstcs STAT 365/665 1/28 Notes: Problem 3 s posted and due ths upcomng Frday There was an early bug n the fake-test data; fxed
More informationLeast Squares Fitting of Data
Least Squares Fttng of Data Davd Eberly Geoetrc Tools, LLC http://www.geoetrctools.co/ Copyrght c 1998-2014. All Rghts Reserved. Created: July 15, 1999 Last Modfed: February 9, 2008 Contents 1 Lnear Fttng
More informationGeneralized Linear Methods
Generalzed Lnear Methods 1 Introducton In the Ensemble Methods the general dea s that usng a combnaton of several weak learner one could make a better learner. More formally, assume that we have a set
More informationONE-DIMENSIONAL COLLISIONS
Purpose Theory ONE-DIMENSIONAL COLLISIONS a. To very the law o conservaton o lnear momentum n one-dmensonal collsons. b. To study conservaton o energy and lnear momentum n both elastc and nelastc onedmensonal
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 informationExercises of Fundamentals of Chemical Processes
Department of Energ Poltecnco d Mlano a Lambruschn 4 2056 MILANO Exercses of undamentals of Chemcal Processes Prof. Ganpero Gropp Exercse 7 ) Estmaton of the composton of the streams at the ext of an sothermal
More informationDEMO #8 - GAUSSIAN ELIMINATION USING MATHEMATICA. 1. Matrices in Mathematica
demo8.nb 1 DEMO #8 - GAUSSIAN ELIMINATION USING MATHEMATICA Obectves: - defne matrces n Mathematca - format the output of matrces - appl lnear algebra to solve a real problem - Use Mathematca to perform
More informationMA 323 Geometric Modelling Course Notes: Day 13 Bezier Curves & Bernstein Polynomials
MA 323 Geometrc Modellng Course Notes: Day 13 Bezer Curves & Bernsten Polynomals Davd L. Fnn Over the past few days, we have looked at de Casteljau s algorthm for generatng a polynomal curve, and we have
More informationDUE: WEDS FEB 21ST 2018
HOMEWORK # 1: FINITE DIFFERENCES IN ONE DIMENSION DUE: WEDS FEB 21ST 2018 1. Theory Beam bendng s a classcal engneerng analyss. The tradtonal soluton technque makes smplfyng assumptons such as a constant
More informationSENSITIVITY APPROACH TO OPTIMAL CONTROL FOR AFFINE NONLINEAR DISCRETE-TIME SYSTEMS
448 Asan Journal o Control Vol. 7 No. 4 pp. 448-454 December 25 Bre Paper SENSIIVIY APPROACH O OPIMAL CONROL FOR AFFINE NONLINEAR DISCREE-IME SYSEMS Gong-You ang Nan Xe Peng Lu ABSRAC hs paper deals wth
More informationwhere the sums are over the partcle labels. In general H = p2 2m + V s(r ) V j = V nt (jr, r j j) (5) where V s s the sngle-partcle potental and V nt
Physcs 543 Quantum Mechancs II Fall 998 Hartree-Fock and the Self-consstent Feld Varatonal Methods In the dscusson of statonary perturbaton theory, I mentoned brey the dea of varatonal approxmaton schemes.
More information= = = (a) Use the MATLAB command rref to solve the system. (b) Let A be the coefficient matrix and B be the right-hand side of the system.
Chapter Matlab Exercses Chapter Matlab Exercses. Consder the lnear system of Example n Secton.. x x x y z y y z (a) Use the MATLAB command rref to solve the system. (b) Let A be the coeffcent matrx and
More informationOpen Systems: Chemical Potential and Partial Molar Quantities Chemical Potential
Open Systems: Chemcal Potental and Partal Molar Quanttes Chemcal Potental For closed systems, we have derved the followng relatonshps: du = TdS pdv dh = TdS + Vdp da = SdT pdv dg = VdP SdT For open systems,
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 informationSUMMARY OF STOICHIOMETRIC RELATIONS AND MEASURE OF REACTIONS' PROGRESS AND COMPOSITION FOR MULTIPLE REACTIONS
UMMAY OF TOICHIOMETIC ELATION AND MEAUE OF EACTION' POGE AND COMPOITION FO MULTIPLE EACTION UPDATED 0/4/03 - AW APPENDIX A. In case of multple reactons t s mportant to fnd the number of ndependent reactons.
More informationTranslational Equations of Motion for A Body Translational equations of motion (centroidal) for a body are m r = f.
Lesson 12: Equatons o Moton Newton s Laws Frst Law: A artcle remans at rest or contnues to move n a straght lne wth constant seed there s no orce actng on t Second Law: The acceleraton o a artcle s roortonal
More informationErrors for Linear Systems
Errors for Lnear Systems When we solve a lnear system Ax b we often do not know A and b exactly, but have only approxmatons  and ˆb avalable. Then the best thng we can do s to solve ˆx ˆb exactly whch
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 informationDepartment of Electrical and Computer Engineering FEEDBACK AMPLIFIERS
Department o Electrcal and Computer Engneerng UNIT I EII FEEDBCK MPLIFIES porton the output sgnal s ed back to the nput o the ampler s called Feedback mpler. Feedback Concept: block dagram o an ampler
More informationCHAPTER 7 CONSTRAINED OPTIMIZATION 2: SQP AND GRG
Chapter 7: Constraned Optmzaton CHAPER 7 CONSRAINED OPIMIZAION : SQP AND GRG Introducton In the prevous chapter we eamned the necessary and suffcent condtons for a constraned optmum. We dd not, however,
More informationA Simple Research of Divisor Graphs
The 29th Workshop on Combnatoral Mathematcs and Computaton Theory A Smple Research o Dvsor Graphs Yu-png Tsao General Educaton Center Chna Unversty o Technology Tape Tawan yp-tsao@cuteedutw Tape Tawan
More informationFormal solvers of the RT equation
Formal solvers of the RT equaton Formal RT solvers Runge- Kutta (reference solver) Pskunov N.: 979, Master Thess Long characterstcs (Feautrer scheme) Cannon C.J.: 970, ApJ 6, 55 Short characterstcs (Hermtan
More information6.3.4 Modified Euler s method of integration
6.3.4 Modfed Euler s method of ntegraton Before dscussng the applcaton of Euler s method for solvng the swng equatons, let us frst revew the basc Euler s method of numercal ntegraton. Let the general from
More informationGoal Programming Approach to Solve Multi- Objective Intuitionistic Fuzzy Non- Linear Programming Models
Internatonal Journal o Mathematcs rends and echnoloy IJM Volume Number 7 - January 8 Goal Prorammn Approach to Solve Mult- Objectve Intutonstc Fuzzy Non- Lnear Prorammn Models S.Rukman #, R.Sopha Porchelv
More informationMoments of Inertia. and reminds us of the analogous equation for linear momentum p= mv, which is of the form. The kinetic energy of the body is.
Moments of Inerta Suppose a body s movng on a crcular path wth constant speed Let s consder two quanttes: the body s angular momentum L about the center of the crcle, and ts knetc energy T How are these
More informationHidden Markov Models & The Multivariate Gaussian (10/26/04)
CS281A/Stat241A: Statstcal Learnng Theory Hdden Markov Models & The Multvarate Gaussan (10/26/04) Lecturer: Mchael I. Jordan Scrbes: Jonathan W. Hu 1 Hdden Markov Models As a bref revew, hdden Markov models
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 informationNON-CENTRAL 7-POINT FORMULA IN THE METHOD OF LINES FOR PARABOLIC AND BURGERS' EQUATIONS
IJRRAS 8 (3 September 011 www.arpapress.com/volumes/vol8issue3/ijrras_8_3_08.pdf NON-CENTRAL 7-POINT FORMULA IN THE METHOD OF LINES FOR PARABOLIC AND BURGERS' EQUATIONS H.O. Bakodah Dept. of Mathematc
More informationPHYS 1443 Section 004 Lecture #12 Thursday, Oct. 2, 2014
PHYS 1443 Secton 004 Lecture #1 Thursday, Oct., 014 Work-Knetc Energy Theorem Work under rcton Potental Energy and the Conservatve Force Gravtatonal Potental Energy Elastc Potental Energy Conservaton o
More informationLinear Momentum. Equation 1
Lnear Momentum OBJECTIVE Obsere collsons between two carts, testng or the conseraton o momentum. Measure energy changes durng derent types o collsons. Classy collsons as elastc, nelastc, or completely
More informationA new Approach for Solving Linear Ordinary Differential Equations
, ISSN 974-57X (Onlne), ISSN 974-5718 (Prnt), Vol. ; Issue No. 1; Year 14, Copyrght 13-14 by CESER PUBLICATIONS A new Approach for Solvng Lnear Ordnary Dfferental Equatons Fawz Abdelwahd Department of
More informationLecture 16 Statistical Analysis in Biomaterials Research (Part II)
3.051J/0.340J 1 Lecture 16 Statstcal Analyss n Bomaterals Research (Part II) C. F Dstrbuton Allows comparson of varablty of behavor between populatons usng test of hypothess: σ x = σ x amed for Brtsh statstcan
More informationTransfer Functions. Convenient representation of a linear, dynamic model. A transfer function (TF) relates one input and one output: ( ) system
Transfer Functons Convenent representaton of a lnear, dynamc model. A transfer functon (TF) relates one nput and one output: x t X s y t system Y s The followng termnology s used: x y nput output forcng
More informationPropagation of error for multivariable function
Propagaton o error or multvarable uncton ow consder a multvarable uncton (u, v, w, ). I measurements o u, v, w,. All have uncertant u, v, w,., how wll ths aect the uncertant o the uncton? L tet) o (Equaton
More informationLesson 16: Basic Control Modes
0/8/05 Lesson 6: Basc Control Modes ET 438a Automatc Control Systems Technology lesson6et438a.tx Learnng Objectves Ater ths resentaton you wll be able to: Descrbe the common control modes used n analog
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 information