DEMO #8 - GAUSSIAN ELIMINATION USING MATHEMATICA. 1. Matrices in Mathematica

Size: px
Start display at page:

Download "DEMO #8 - GAUSSIAN ELIMINATION USING MATHEMATICA. 1. Matrices in Mathematica"

Transcription

1 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 Gauss-Jordan elmnaton on a matrx Bran Knowlton ChE 50B, Chemcal Engneerng Department, UC Santa Barbara created: 1/7/98 revsed: 1/7/98 revsed b G. Fredrckson for Eng 5A: 10/27/98 1. Matrces n Mathematca In Mathematca, matrces are called Lsts. A Lst s ust a collecton of numbers, separated b commas, and surrounded b curl braces { }. For example, the followng defnes a 1x matrx. u = 8 1, 2, < 81, 2, < You can also use varables n matrces. v = 8x,, < 8x,, < The elements of a Lst don't have to be numbers or varables; the can also be other lsts. For example, a 2x5 matrx can be thought as a 2x1 Lst, where both elements are 1x5 Lsts. Therefore, to specf a 2x5 matrx n Mathematca, specf a two element Lst, where each element s a fve element Lst. Ths tpe of Lst s also referred to as a nested Lst, snce the fve element Lsts resde are 'nested' n another Lst. A = 8 8a11, a12, a1, a14, b1<, 8a21, a22, a2, a24, b2< < 88a11, a12, a1, a14, b1<, 8a21, a22, a2, a24, b2<< The matrx A mght be an augmented matrx derved from a par of lnear equatons n 4 unknowns, a 11 x 1 + a 12 x 2 + a 1 x + a 14 x 4 = b 1 a 21 x 1 + a 22 x 2 + a 2 x + a 24 x 4 = b 2

2 demo8.nb 2 where a are coeffcents multplng the unknowns x. The matrces u and v are called one-dmensonal matrces, where A s a two-dmensonal matrx. You can also specf three or hgher dmensonal matrces n Mathematca b further nestng of Lsts, and although we wll not use these n Eng 5A, the are useful mathematcal structures n varous branches of phscs and engneerng. 2. Eas-To-Read Matrx Output The format that Mathematca uses for matrx output s a lttle dffcult to read. Mathematca provdes a functon, MatrxForm, that produces more readable output. MatrxForm@AD a11 a12 a1 a14 b1 J a21 a22 a2 a24 b2 N. Applng Lnear Algebra to a Real Problem In Fgure 1 s a crcut dagram, representng a collecton of tpcal household applances, connected (perhaps unsafel) to a sngle wall outlet. The power that each applance uses s lsted n the fgure, along wth the voltage provded b the outlet. The problem s to predct the amount of current that s requred to run all the applances at once. Ths nformaton can be mportant, snce f too much current s drawn from the outlet, a fuse wll fal and ou wll lose to power to all applances. Fgure 1: Crcut Dagram of Common Household Applances In order to fnd the current passng through the outlet, we need to fnd the current runnng through each applance and the voltage at varous places n the crcut. The frst step n analng ths crcut s to fnd the resstance of each applance. Usng the nformaton suppled b the manufacturer, the resstance for each applance can be determned. Fgure 2 shows the resstance of each applance, and also the remanng unknown currents and voltages.

3 demo8.nb Fgure 2: Applance Resstances and Labelng Notaton of Unknowns To fnd the unknowns, we use two phscal laws to wrte down lnear equatons descrbng the behavor of ths crcut. Frst, at each pont where the crcut branches or ons (ponts A or C), we emplo a "balancng equaton" of currents. Because current s ust the movement of electrons, the total current enterng an gven pont must equal the total current flowng out. Usng the notaton n Fgure 2, we can wrte the followng two (redundant) equatons. I 1 = I 2 + I I 2 + I = I 1 Equatons can be wrtten usng Ohm's Law, whch descrbes how the electrc potental drops when current flows through a resstor. Ohm's Law sas that the voltage drop across a resstor s proportonal to the current flowng through the resstor, where the constant of proportonalt s the resstance. As an equaton, Ohm's Law can be wrtten as, V = IR Applng ths law to each of the applances n Fgure 2, we can wrte fve more equatons, V 1 = 12 I 1 V 1 - V 2 = 9 I V 2 - V = 60 I V 1 - V = (2/) I 2 V - 0 = (2/) I 1 In total, we have seven equatons n sx unknowns. Wrtten n a more standard form, the equatons look lke the followng. I 1 - I 2 - I = 0 - I 1 + I 2 + I = 0 12 I 1 + V 1 = I + V 1 - V 2 = 0-60 I + V 2 - V = 0 - (2/) I 2 + V 1 - V = 0

4 demo8.nb 4 - (2/) I 1 + V = 0 NOTE: We have chosen here to nclude the redundant equaton. In man applcatons, the fact that two equatons are redundant wll not be so obvous. Also, ou wll see that ncludng a redundant equaton wll not affect the soluton to the problem (see below). 4. Usng Mathematca To Perform Gauss-Jordan Elmnaton Usng Mathematca, we can wrte the equatons nto an augmented matrx: A = 8 8 1, 1, 1, 0, 0, 0, 0<, 8 1, 1, 1, 0, 0, 0, 0<, 812, 0, 0, 1, 0, 0, 120<, 8 0, 0, 9, 1, 1, 0, 0<, 8 0, 0, 60, 0, 1, 1, 0<, 8 0, 2, 0, 1, 0, 1, 0<, 8 2, 0, 0, 0, 0, 1, 0< <; MatrxForm@AD k { To perform Gauss-Jordan Elmnaton on ths augmented matrx, we can call the alread prepared Mathematca functon RowReduce, whch performs Gauss-Jordan elmnaton, producng a matrx n reduced row-echelon form. B = RowReduce@AD; MatrxForm@BD k {

5 demo8.nb 5 N@ MatrxForm@BDD k The fnal answer s: I 1 =.6; I 2 =.51; I = 0.102; V 1 = ; V 2 = ; V = 8.748; The total amperage requred to run all the applances at once s I 1 or.6 Amperes. { NOTE: The matrx B, the reduced matrx, has a row of eros as the last row. Ths row s the result of the redundant equaton that we ncluded n our sstem of equatons. An equaton that s not ndependant of the others wll show up as a row of eros, and can be neglected. Appendx 1: An Even Better Wa of Formattng Matrces The use of the MatrxForm functon s a rather cumbersome wa to format matrx output, snce ou must call the functon ever tme ou output a matrx. However, there s a wa ou can tell Mathematca to alwas prnt Lsts n matrx form. The followng command, when evaluated, wll cause Mathematca, b default, to use the MatrxForm functon for output of Lsts. $Post := If@MatrxQ@#1D, MatrxForm@#1D, #1D& Tpcall, ths command would come at the begnnng of the notebook. In order to restore the orgnal formattng, smpl clear the varable $Post. See the followng example for a demonstraton.

6 demo8.nb 6 Clear@$PostD A $Post := If@MatrxQ@#1D, MatrxForm@#1D, #1D& A 981, 1, 1, 0, 0, 0, 0<, 8 1, 1, 1, 0, 0, 0, 0<, 812, 0, 0, 1, 0, 0, 120<, 80, 0, 9, 1, 1, 0, 0<, 80, 0, 60, 0, 1, 1, 0<, 90, 2, 0, 1, 0, 1, 0=, 9 2, 0, 0, 0, 0, 1, 0== k {

MAE140 - Linear Circuits - Fall 13 Midterm, October 31

MAE140 - Linear Circuits - Fall 13 Midterm, October 31 Instructons ME140 - Lnear Crcuts - Fall 13 Mdterm, October 31 () Ths exam s open book. You may use whatever wrtten materals you choose, ncludng your class notes and textbook. You may use a hand calculator

More information

MAE140 - Linear Circuits - Winter 16 Midterm, February 5

MAE140 - Linear Circuits - Winter 16 Midterm, February 5 Instructons ME140 - Lnear Crcuts - Wnter 16 Mdterm, February 5 () Ths exam s open book. You may use whatever wrtten materals you choose, ncludng your class notes and textbook. You may use a hand calculator

More information

The Fundamental Theorem of Algebra. Objective To use the Fundamental Theorem of Algebra to solve polynomial equations with complex solutions

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

Introduction to Vapor/Liquid Equilibrium, part 2. Raoult s Law:

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

Advanced Circuits Topics - Part 1 by Dr. Colton (Fall 2017)

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

MAE140 - Linear Circuits - Fall 10 Midterm, October 28

MAE140 - Linear Circuits - Fall 10 Midterm, October 28 M140 - Lnear rcuts - Fall 10 Mdterm, October 28 nstructons () Ths exam s open book. You may use whatever wrtten materals you choose, ncludng your class notes and textbook. You may use a hand calculator

More information

n α j x j = 0 j=1 has a nontrivial solution. Here A is the n k matrix whose jth column is the vector for all t j=0

n α j x j = 0 j=1 has a nontrivial solution. Here A is the n k matrix whose jth column is the vector for all t j=0 MODULE 2 Topcs: Lnear ndependence, bass and dmenson We have seen that f n a set of vectors one vector s a lnear combnaton of the remanng vectors n the set then the span of the set s unchanged f that vector

More information

Physics 4B. A positive value is obtained, so the current is counterclockwise around the circuit.

Physics 4B. A positive value is obtained, so the current is counterclockwise around the circuit. Physcs 4B Solutons to Chapter 7 HW Chapter 7: Questons:, 8, 0 Problems:,,, 45, 48,,, 7, 9 Queston 7- (a) no (b) yes (c) all te Queston 7-8 0 μc Queston 7-0, c;, a;, d; 4, b Problem 7- (a) Let be the current

More information

Lecture 10 Support Vector Machines II

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

Formulation of Circuit Equations

Formulation of Circuit Equations ECE 570 Sesson 2 IC 752E Computer Aded Engneerng for Integrated Crcuts Formulaton of Crcut Equatons Bascs of crcut modelng 1. Notaton 2. Crcut elements 3. Krchoff laws 4. ableau formulaton 5. Modfed nodal

More information

Key component in Operational Amplifiers

Key component in Operational Amplifiers Key component n Operatonal Amplfers Objectve of Lecture Descrbe how dependent voltage and current sources functon. Chapter.6 Electrcal Engneerng: Prncples and Applcatons Chapter.6 Fundamentals of Electrc

More information

Selected Student Solutions for Chapter 2

Selected Student Solutions for Chapter 2 /3/003 Assessment Prolems Selected Student Solutons for Chapter. Frst note that we know the current through all elements n the crcut except the 6 kw resstor (the current n the three elements to the left

More information

1 Matrix representations of canonical matrices

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

Transfer Functions. Convenient representation of a linear, dynamic model. A transfer function (TF) relates one input and one output: ( ) system

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

KIRCHHOFF CURRENT LAW

KIRCHHOFF CURRENT LAW KRCHHOFF CURRENT LAW ONE OF THE FUNDAMENTAL CONSERATON PRNCPLES N ELECTRCAL ENGNEERNG CHARGE CANNOT BE CREATED NOR DESTROYED NODES, BRANCHES, LOOPS A NODE CONNECTS SEERAL COMPONENTS. BUT T DOES NOT HOLD

More information

CinChE Problem-Solving Strategy Chapter 4 Development of a Mathematical Model. formulation. procedure

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

Electrical Circuits 2.1 INTRODUCTION CHAPTER

Electrical Circuits 2.1 INTRODUCTION CHAPTER CHAPTE Electrcal Crcuts. INTODUCTION In ths chapter, we brefly revew the three types of basc passve electrcal elements: resstor, nductor and capactor. esstance Elements: Ohm s Law: The voltage drop across

More information

Section 3.6 Complex Zeros

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

C/CS/Phy191 Problem Set 3 Solutions Out: Oct 1, 2008., where ( 00. ), so the overall state of the system is ) ( ( ( ( 00 ± 11 ), Φ ± = 1

C/CS/Phy191 Problem Set 3 Solutions Out: Oct 1, 2008., where ( 00. ), so the overall state of the system is ) ( ( ( ( 00 ± 11 ), Φ ± = 1 C/CS/Phy9 Problem Set 3 Solutons Out: Oct, 8 Suppose you have two qubts n some arbtrary entangled state ψ You apply the teleportaton protocol to each of the qubts separately What s the resultng state obtaned

More information

Complex Numbers, Signals, and Circuits

Complex Numbers, Signals, and Circuits Complex Numbers, Sgnals, and Crcuts 3 August, 009 Complex Numbers: a Revew Suppose we have a complex number z = x jy. To convert to polar form, we need to know the magntude of z and the phase of z. z =

More information

Polynomial Regression Models

Polynomial 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

Inner Product. Euclidean Space. Orthonormal Basis. Orthogonal

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

SUMMARY OF STOICHIOMETRIC RELATIONS AND MEASURE OF REACTIONS' PROGRESS AND COMPOSITION FOR MULTIPLE REACTIONS

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

Structure and Drive Paul A. Jensen Copyright July 20, 2003

Structure and Drive Paul A. Jensen Copyright July 20, 2003 Structure and Drve Paul A. Jensen Copyrght July 20, 2003 A system s made up of several operatons wth flow passng between them. The structure of the system descrbes the flow paths from nputs to outputs.

More information

8.4 COMPLEX VECTOR SPACES AND INNER PRODUCTS

8.4 COMPLEX VECTOR SPACES AND INNER PRODUCTS SECTION 8.4 COMPLEX VECTOR SPACES AND INNER PRODUCTS 493 8.4 COMPLEX VECTOR SPACES AND INNER PRODUCTS All the vector spaces you have studed thus far n the text are real vector spaces because the scalars

More information

Resource Allocation and Decision Analysis (ECON 8010) Spring 2014 Foundations of Regression Analysis

Resource Allocation and Decision Analysis (ECON 8010) Spring 2014 Foundations of Regression Analysis Resource Allocaton and Decson Analss (ECON 800) Sprng 04 Foundatons of Regresson Analss Readng: Regresson Analss (ECON 800 Coursepak, Page 3) Defntons and Concepts: Regresson Analss statstcal technques

More information

Pre-Calculus Summer Assignment

Pre-Calculus Summer Assignment Pre-Calculus Summer Assgnment Dear Future Pre-Calculus Student, Congratulatons on our successful completon of Algebra! Below ou wll fnd the summer assgnment questons. It s assumed that these concepts,

More information

Robert Eisberg Second edition CH 09 Multielectron atoms ground states and x-ray excitations

Robert Eisberg Second edition CH 09 Multielectron atoms ground states and x-ray excitations Quantum Physcs 量 理 Robert Esberg Second edton CH 09 Multelectron atoms ground states and x-ray exctatons 9-01 By gong through the procedure ndcated n the text, develop the tme-ndependent Schroednger equaton

More information

Formulas for the Determinant

Formulas for the Determinant page 224 224 CHAPTER 3 Determnants e t te t e 2t 38 A = e t 2te t e 2t e t te t 2e 2t 39 If 123 A = 345, 456 compute the matrx product A adj(a) What can you conclude about det(a)? For Problems 40 43, use

More information

Salmon: Lectures on partial differential equations. Consider the general linear, second-order PDE in the form. ,x 2

Salmon: Lectures on partial differential equations. Consider the general linear, second-order PDE in the form. ,x 2 Salmon: Lectures on partal dfferental equatons 5. Classfcaton of second-order equatons There are general methods for classfyng hgher-order partal dfferental equatons. One s very general (applyng even to

More information

2.3 Nilpotent endomorphisms

2.3 Nilpotent endomorphisms s a block dagonal matrx, wth A Mat dm U (C) In fact, we can assume that B = B 1 B k, wth B an ordered bass of U, and that A = [f U ] B, where f U : U U s the restrcton of f to U 40 23 Nlpotent endomorphsms

More information

EE 2006 Electric Circuit Analysis Fall September 04, 2014 Lecture 02

EE 2006 Electric Circuit Analysis Fall September 04, 2014 Lecture 02 EE 2006 Electrc Crcut Analyss Fall 2014 September 04, 2014 Lecture 02 1 For Your Informaton Course Webpage http://www.d.umn.edu/~jngba/electrc_crcut_analyss_(ee_2006).html You can fnd on the webpage: Lecture:

More information

The Exact Formulation of the Inverse of the Tridiagonal Matrix for Solving the 1D Poisson Equation with the Finite Difference Method

The Exact Formulation of the Inverse of the Tridiagonal Matrix for Solving the 1D Poisson Equation with the Finite Difference Method Journal of Electromagnetc Analyss and Applcatons, 04, 6, 0-08 Publshed Onlne September 04 n ScRes. http://www.scrp.org/journal/jemaa http://dx.do.org/0.46/jemaa.04.6000 The Exact Formulaton of the Inverse

More information

Econ107 Applied Econometrics Topic 3: Classical Model (Studenmund, Chapter 4)

Econ107 Applied Econometrics Topic 3: Classical Model (Studenmund, Chapter 4) I. Classcal Assumptons Econ7 Appled Econometrcs Topc 3: Classcal Model (Studenmund, Chapter 4) We have defned OLS and studed some algebrac propertes of OLS. In ths topc we wll study statstcal propertes

More information

Physics 5153 Classical Mechanics. D Alembert s Principle and The Lagrangian-1

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

MAE140 - Linear Circuits - Winter 16 Final, March 16, 2016

MAE140 - Linear Circuits - Winter 16 Final, March 16, 2016 ME140 - Lnear rcuts - Wnter 16 Fnal, March 16, 2016 Instructons () The exam s open book. You may use your class notes and textbook. You may use a hand calculator wth no communcaton capabltes. () You have

More information

ELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM

ELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM ELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM An elastc wave s a deformaton of the body that travels throughout the body n all drectons. We can examne the deformaton over a perod of tme by fxng our look

More information

G = G 1 + G 2 + G 3 G 2 +G 3 G1 G2 G3. Network (a) Network (b) Network (c) Network (d)

G = G 1 + G 2 + G 3 G 2 +G 3 G1 G2 G3. Network (a) Network (b) Network (c) Network (d) Massachusetts Insttute of Technology Department of Electrcal Engneerng and Computer Scence 6.002 í Electronc Crcuts Homework 2 Soluton Handout F98023 Exercse 21: Determne the conductance of each network

More information

PHYSICS - CLUTCH CH 28: INDUCTION AND INDUCTANCE.

PHYSICS - CLUTCH CH 28: INDUCTION AND INDUCTANCE. !! www.clutchprep.com CONCEPT: ELECTROMAGNETIC INDUCTION A col of wre wth a VOLTAGE across each end wll have a current n t - Wre doesn t HAVE to have voltage source, voltage can be INDUCED V Common ways

More information

Chapter 12. Ordinary Differential Equation Boundary Value (BV) Problems

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

Kinematics in 2-Dimensions. Projectile Motion

Kinematics in 2-Dimensions. Projectile Motion Knematcs n -Dmensons Projectle Moton A medeval trebuchet b Kolderer, c1507 http://members.net.net.au/~rmne/ht/ht0.html#5 Readng Assgnment: Chapter 4, Sectons -6 Introducton: In medeval das, people had

More information

Fundamental loop-current method using virtual voltage sources technique for special cases

Fundamental loop-current method using virtual voltage sources technique for special cases Fundamental loop-current method usng vrtual voltage sources technque for specal cases George E. Chatzaraks, 1 Marna D. Tortorel 1 and Anastasos D. Tzolas 1 Electrcal and Electroncs Engneerng Departments,

More information

LECTURE 9 CANONICAL CORRELATION ANALYSIS

LECTURE 9 CANONICAL CORRELATION ANALYSIS LECURE 9 CANONICAL CORRELAION ANALYSIS Introducton he concept of canoncal correlaton arses when we want to quantfy the assocatons between two sets of varables. For example, suppose that the frst set of

More information

Chapter Newton s Method

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

Linear Feature Engineering 11

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

Chapter 6 Electrical Systems and Electromechanical Systems

Chapter 6 Electrical Systems and Electromechanical Systems ME 43 Systems Dynamcs & Control Chapter 6: Electrcal Systems and Electromechancal Systems Chapter 6 Electrcal Systems and Electromechancal Systems 6. INTODUCTION A. Bazoune The majorty of engneerng systems

More information

Finite Differences, Interpolation, and Numerical Differentiation

Finite Differences, Interpolation, and Numerical Differentiation 4 Fnte Dfferences, Interpolaton, and Numercal Dfferentaton 4. INTRODUCTION Lnear nterpolaton s dscussed n the precedng chapter as a method for fndng a partcular root of a polynomal, or, transcendental

More information

Sections begin this week. Cancelled Sections: Th Labs begin this week. Attend your only second lab slot this week.

Sections begin this week. Cancelled Sections: Th Labs begin this week. Attend your only second lab slot this week. Announcements Sectons begn ths week Cancelled Sectons: Th 122. Labs begn ths week. Attend your only second lab slot ths week. Cancelled labs: ThF 25. Please check your Lab secton. Homework #1 onlne Due

More information

Important Instructions to the Examiners:

Important Instructions to the Examiners: Summer 0 Examnaton Subject & Code: asc Maths (70) Model Answer Page No: / Important Instructons to the Examners: ) The Answers should be examned by key words and not as word-to-word as gven n the model

More information

U.C. Berkeley CS294: Beyond Worst-Case Analysis Luca Trevisan September 5, 2017

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

Analytical Chemistry Calibration Curve Handout

Analytical Chemistry Calibration Curve Handout I. Quck-and Drty Excel Tutoral Analytcal Chemstry Calbraton Curve Handout For those of you wth lttle experence wth Excel, I ve provded some key technques that should help you use the program both for problem

More information

Lectures - Week 4 Matrix norms, Conditioning, Vector Spaces, Linear Independence, Spanning sets and Basis, Null space and Range of a Matrix

Lectures - Week 4 Matrix norms, Conditioning, Vector Spaces, Linear Independence, Spanning sets and Basis, Null space and Range of a Matrix Lectures - Week 4 Matrx norms, Condtonng, Vector Spaces, Lnear Independence, Spannng sets and Bass, Null space and Range of a Matrx Matrx Norms Now we turn to assocatng a number to each matrx. We could

More information

Chapter 11: Simple Linear Regression and Correlation

Chapter 11: Simple Linear Regression and Correlation Chapter 11: Smple Lnear Regresson and Correlaton 11-1 Emprcal Models 11-2 Smple Lnear Regresson 11-3 Propertes of the Least Squares Estmators 11-4 Hypothess Test n Smple Lnear Regresson 11-4.1 Use of t-tests

More information

1 Derivation of Rate Equations from Single-Cell Conductance (Hodgkin-Huxley-like) Equations

1 Derivation of Rate Equations from Single-Cell Conductance (Hodgkin-Huxley-like) Equations Physcs 171/271 -Davd Klenfeld - Fall 2005 (revsed Wnter 2011) 1 Dervaton of Rate Equatons from Sngle-Cell Conductance (Hodgkn-Huxley-lke) Equatons We consder a network of many neurons, each of whch obeys

More information

Lecture 3. Ax x i a i. i i

Lecture 3. Ax x i a i. i i 18.409 The Behavor of Algorthms n Practce 2/14/2 Lecturer: Dan Spelman Lecture 3 Scrbe: Arvnd Sankar 1 Largest sngular value In order to bound the condton number, we need an upper bound on the largest

More information

The optimal delay of the second test is therefore approximately 210 hours earlier than =2.

The optimal delay of the second test is therefore approximately 210 hours earlier than =2. THE IEC 61508 FORMULAS 223 The optmal delay of the second test s therefore approxmately 210 hours earler than =2. 8.4 The IEC 61508 Formulas IEC 61508-6 provdes approxmaton formulas for the PF for smple

More information

Foundations of Arithmetic

Foundations of Arithmetic Foundatons of Arthmetc Notaton We shall denote the sum and product of numbers n the usual notaton as a 2 + a 2 + a 3 + + a = a, a 1 a 2 a 3 a = a The notaton a b means a dvdes b,.e. ac = b where c s an

More information

β0 + β1xi. You are interested in estimating the unknown parameters β

β0 + β1xi. You are interested in estimating the unknown parameters β Revsed: v3 Ordnar Least Squares (OLS): Smple Lnear Regresson (SLR) Analtcs The SLR Setup Sample Statstcs Ordnar Least Squares (OLS): FOCs and SOCs Back to OLS and Sample Statstcs Predctons (and Resduals)

More information

Homework Notes Week 7

Homework Notes Week 7 Homework Notes Week 7 Math 4 Sprng 4 #4 (a Complete the proof n example 5 that s an nner product (the Frobenus nner product on M n n (F In the example propertes (a and (d have already been verfed so we

More information

e i is a random error

e i is a random error Chapter - The Smple Lnear Regresson Model The lnear regresson equaton s: where + β + β e for,..., and are observable varables e s a random error How can an estmaton rule be constructed for the unknown

More information

Poisson brackets and canonical transformations

Poisson brackets and canonical transformations rof O B Wrght Mechancs Notes osson brackets and canoncal transformatons osson Brackets Consder an arbtrary functon f f ( qp t) df f f f q p q p t But q p p where ( qp ) pq q df f f f p q q p t In order

More information

Linear Approximation with Regularization and Moving Least Squares

Linear Approximation with Regularization and Moving Least Squares Lnear Approxmaton wth Regularzaton and Movng Least Squares Igor Grešovn May 007 Revson 4.6 (Revson : March 004). 5 4 3 0.5 3 3.5 4 Contents: Lnear Fttng...4. Weghted Least Squares n Functon Approxmaton...

More information

Numerical Solution of Ordinary Differential Equations

Numerical Solution of Ordinary Differential Equations Numercal Methods (CENG 00) CHAPTER-VI Numercal Soluton of Ordnar Dfferental Equatons 6 Introducton Dfferental equatons are equatons composed of an unknown functon and ts dervatves The followng are examples

More information

EE 2006 Electric Circuit Analysis Spring January 23, 2015 Lecture 02

EE 2006 Electric Circuit Analysis Spring January 23, 2015 Lecture 02 EE 2006 Electrc Crcut Analyss Sprng 2015 January 23, 2015 Lecture 02 1 Lab 1 Dgtal Multmeter Lab nstructons Aalable onlne Prnt out and read before Lab MWAH 391, 4:00 7:00 pm, next Monday or Wednesday (January

More information

NUMERICAL DIFFERENTIATION

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

p 1 c 2 + p 2 c 2 + p 3 c p m c 2

p 1 c 2 + p 2 c 2 + p 3 c p m c 2 Where to put a faclty? Gven locatons p 1,..., p m n R n of m houses, want to choose a locaton c n R n for the fre staton. Want c to be as close as possble to all the house. We know how to measure dstance

More information

x yi In chapter 14, we want to perform inference (i.e. calculate confidence intervals and perform tests of significance) in this setting.

x yi In chapter 14, we want to perform inference (i.e. calculate confidence intervals and perform tests of significance) in this setting. The Practce of Statstcs, nd ed. Chapter 14 Inference for Regresson Introducton In chapter 3 we used a least-squares regresson lne (LSRL) to represent a lnear relatonshp etween two quanttatve explanator

More information

APPENDIX 2 FITTING A STRAIGHT LINE TO OBSERVATIONS

APPENDIX 2 FITTING A STRAIGHT LINE TO OBSERVATIONS Unversty of Oulu Student Laboratory n Physcs Laboratory Exercses n Physcs 1 1 APPEDIX FITTIG A STRAIGHT LIE TO OBSERVATIOS In the physcal measurements we often make a seres of measurements of the dependent

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

Chapter 13: Multiple Regression

Chapter 13: Multiple Regression Chapter 13: Multple Regresson 13.1 Developng the multple-regresson Model The general model can be descrbed as: It smplfes for two ndependent varables: The sample ft parameter b 0, b 1, and b are used to

More information

Indeterminate pin-jointed frames (trusses)

Indeterminate pin-jointed frames (trusses) Indetermnate pn-jonted frames (trusses) Calculaton of member forces usng force method I. Statcal determnacy. The degree of freedom of any truss can be derved as: w= k d a =, where k s the number of all

More information

UNIVERSITY OF UTAH ELECTRICAL & COMPUTER ENGINEERING DEPARTMENT. 10k. 3mH. 10k. Only one current in the branch:

UNIVERSITY OF UTAH ELECTRICAL & COMPUTER ENGINEERING DEPARTMENT. 10k. 3mH. 10k. Only one current in the branch: UNIERSITY OF UTH ELECTRICL & COMPUTER ENGINEERING DEPRTMENT ECE 70 HOMEWORK #6 Soluton Summer 009. fter beng closed a long tme, the swtch opens at t = 0. Fnd (t) for t > 0. t = 0 0kΩ 0kΩ 3mH Step : (Redraw

More information

PHYSICS - CLUTCH 1E CH 28: INDUCTION AND INDUCTANCE.

PHYSICS - CLUTCH 1E CH 28: INDUCTION AND INDUCTANCE. !! www.clutchprep.com CONCEPT: ELECTROMAGNETIC INDUCTION A col of wre wth a VOLTAGE across each end wll have a current n t - Wre doesn t HAVE to have voltage source, voltage can be INDUCED V Common ways

More information

Physics 1202: Lecture 11 Today s Agenda

Physics 1202: Lecture 11 Today s Agenda Physcs 122: Lecture 11 Today s Agenda Announcements: Team problems start ths Thursday Team 1: Hend Ouda, Mke Glnsk, Stephane Auger Team 2: Analese Bruder, Krsten Dean, Alson Smth Offce hours: Monday 2:3-3:3

More information

FE REVIEW OPERATIONAL AMPLIFIERS (OP-AMPS)( ) 8/25/2010

FE REVIEW OPERATIONAL AMPLIFIERS (OP-AMPS)( ) 8/25/2010 FE REVEW OPERATONAL AMPLFERS (OP-AMPS)( ) 1 The Op-amp 2 An op-amp has two nputs and one output. Note the op-amp below. The termnal labeled l wth the (-) sgn s the nvertng nput and the nput labeled wth

More information

Physics 114 Exam 2 Fall 2014 Solutions. Name:

Physics 114 Exam 2 Fall 2014 Solutions. Name: Physcs 114 Exam Fall 014 Name: For gradng purposes (do not wrte here): Queston 1. 1... 3. 3. Problem Answer each of the followng questons. Ponts for each queston are ndcated n red. Unless otherwse ndcated,

More information

Rockefeller College University at Albany

Rockefeller College University at Albany Rockefeller College Unverst at Alban PAD 705 Handout: Maxmum Lkelhood Estmaton Orgnal b Davd A. Wse John F. Kenned School of Government, Harvard Unverst Modfcatons b R. Karl Rethemeer Up to ths pont n

More information

SIMPLE LINEAR REGRESSION

SIMPLE LINEAR REGRESSION Smple Lnear Regresson and Correlaton Introducton Prevousl, our attenton has been focused on one varable whch we desgnated b x. Frequentl, t s desrable to learn somethng about the relatonshp between two

More information

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

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

More information

EEE 241: Linear Systems

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

Supplemental document

Supplemental document Electronc Supplementary Materal (ESI) for Physcal Chemstry Chemcal Physcs. Ths journal s the Owner Socetes 01 Supplemental document Behnam Nkoobakht School of Chemstry, The Unversty of Sydney, Sydney,

More information

DO NOT DO HOMEWORK UNTIL IT IS ASSIGNED. THE ASSIGNMENTS MAY CHANGE UNTIL ANNOUNCED.

DO NOT DO HOMEWORK UNTIL IT IS ASSIGNED. THE ASSIGNMENTS MAY CHANGE UNTIL ANNOUNCED. EE 539 Homeworks Sprng 08 Updated: Tuesday, Aprl 7, 08 DO NOT DO HOMEWORK UNTIL IT IS ASSIGNED. THE ASSIGNMENTS MAY CHANGE UNTIL ANNOUNCED. For full credt, show all work. Some problems requre hand calculatons.

More information

Expected Value and Variance

Expected Value and Variance MATH 38 Expected Value and Varance Dr. Neal, WKU We now shall dscuss how to fnd the average and standard devaton of a random varable X. Expected Value Defnton. The expected value (or average value, or

More information

DC Circuits. Crossing the emf in this direction +ΔV

DC Circuits. Crossing the emf in this direction +ΔV DC Crcuts Delverng a steady flow of electrc charge to a crcut requres an emf devce such as a battery, solar cell or electrc generator for example. mf stands for electromotve force, but an emf devce transforms

More information

e - c o m p a n i o n

e - c o m p a n i o n OPERATIONS RESEARCH http://dxdoorg/0287/opre007ec e - c o m p a n o n ONLY AVAILABLE IN ELECTRONIC FORM 202 INFORMS Electronc Companon Generalzed Quantty Competton for Multple Products and Loss of Effcency

More information

CHEMICAL REACTIONS AND DIFFUSION

CHEMICAL REACTIONS AND DIFFUSION CHEMICAL REACTIONS AND DIFFUSION A.K.A. NETWORK THERMODYNAMICS BACKGROUND Classcal thermodynamcs descrbes equlbrum states. Non-equlbrum thermodynamcs descrbes steady states. Network thermodynamcs descrbes

More information

An introduction to least-squares fitting

An introduction to least-squares fitting atonale An ntroducton to least-squares fng.p. Pla hp://www.sppla.co.uk Ths paper provdes a mnmall mathematcal ntroducton to least-squares fng, ntended to e of some modest value to engneerng students needng

More information

Nodal analysis of finite square resistive grids and the teaching effectiveness of students projects

Nodal analysis of finite square resistive grids and the teaching effectiveness of students projects 2 nd World Conference on Technology and Engneerng Educaton 2 WIETE Lublana Slovena 5-8 September 2 Nodal analyss of fnte square resstve grds and the teachng effectveness of students proects P. Zegarmstrz

More information

6.01: Introduction to EECS I Lecture 7 March 15, 2011

6.01: Introduction to EECS I Lecture 7 March 15, 2011 6.0: Introducton to EECS I Lecture 7 March 5, 20 6.0: Introducton to EECS I Crcuts The Crcut Abstracton Crcuts represent systems as connectons of elements through whch currents (through arables) flow and

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.

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

Module 2. Random Processes. Version 2 ECE IIT, Kharagpur

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

Negative Binomial Regression

Negative Binomial Regression STATGRAPHICS Rev. 9/16/2013 Negatve Bnomal Regresson Summary... 1 Data Input... 3 Statstcal Model... 3 Analyss Summary... 4 Analyss Optons... 7 Plot of Ftted Model... 8 Observed Versus Predcted... 10 Predctons...

More information

Prof. Dr. I. Nasser Phys 630, T Aug-15 One_dimensional_Ising_Model

Prof. Dr. I. Nasser Phys 630, T Aug-15 One_dimensional_Ising_Model EXACT OE-DIMESIOAL ISIG MODEL The one-dmensonal Isng model conssts of a chan of spns, each spn nteractng only wth ts two nearest neghbors. The smple Isng problem n one dmenson can be solved drectly n several

More information

% & 5.3 PRACTICAL APPLICATIONS. Given system, (49) , determine the Boolean Function, , in such a way that we always have expression: " Y1 = Y2

% & 5.3 PRACTICAL APPLICATIONS. Given system, (49) , determine the Boolean Function, , in such a way that we always have expression:  Y1 = Y2 5.3 PRACTICAL APPLICATIONS st EXAMPLE: Gven system, (49) & K K Y XvX 3 ( 2 & X ), determne the Boolean Functon, Y2 X2 & X 3 v X " X3 (X2,X)", n such a way that we always have expresson: " Y Y2 " (50).

More information

College of Computer & Information Science Fall 2009 Northeastern University 20 October 2009

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

9 Derivation of Rate Equations from Single-Cell Conductance (Hodgkin-Huxley-like) Equations

9 Derivation of Rate Equations from Single-Cell Conductance (Hodgkin-Huxley-like) Equations Physcs 171/271 - Chapter 9R -Davd Klenfeld - Fall 2005 9 Dervaton of Rate Equatons from Sngle-Cell Conductance (Hodgkn-Huxley-lke) Equatons We consder a network of many neurons, each of whch obeys a set

More information

INDUCTANCE. RC Cicuits vs LR Circuits

INDUCTANCE. RC Cicuits vs LR Circuits INDUTANE R cuts vs LR rcuts R rcut hargng (battery s connected): (1/ )q + (R)dq/ dt LR rcut = (R) + (L)d/ dt q = e -t/ R ) = / R(1 - e -(R/ L)t ) q ncreases from 0 to = dq/ dt decreases from / R to 0 Dschargng

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

Introduction to circuit analysis. Classification of Materials

Introduction to circuit analysis. Classification of Materials Introducton to crcut analyss OUTLINE Electrcal quanttes Charge Current Voltage Power The deal basc crcut element Sgn conventons Current versus voltage (I-V) graph Readng: 1.2, 1.3,1.6 Lecture 2, Slde 1

More information

CHAPTER 4. Vector Spaces

CHAPTER 4. Vector Spaces man 2007/2/16 page 234 CHAPTER 4 Vector Spaces To crtcze mathematcs for ts abstracton s to mss the pont entrel. Abstracton s what makes mathematcs work. Ian Stewart The man am of ths tet s to stud lnear

More information