Lecture 3. In this lecture, we will discuss algorithms for solving systems of linear equations.

Similar documents
How do we solve these things, especially when they get complicated? How do we know when a system has a solution, and when is it unique?

Matrix Algebra. Matrix Addition, Scalar Multiplication and Transposition. Linear Algebra I 24

How do we solve these things, especially when they get complicated? How do we know when a system has a solution, and when is it unique?

Bases for Vector Spaces

dx dt dy = G(t, x, y), dt where the functions are defined on I Ω, and are locally Lipschitz w.r.t. variable (x, y) Ω.

Quadratic Forms. Quadratic Forms

Lecture Solution of a System of Linear Equation

Things to Memorize: A Partial List. January 27, 2017

Lecture 2e Orthogonal Complement (pages )

set is not closed under matrix [ multiplication, ] and does not form a group.

Lecture 08: Feb. 08, 2019

Torsion in Groups of Integral Triangles

September 13 Homework Solutions

2. VECTORS AND MATRICES IN 3 DIMENSIONS

Infinite Geometric Series

Section 4: Integration ECO4112F 2011

Chapter 3 MATRIX. In this chapter: 3.1 MATRIX NOTATION AND TERMINOLOGY

Coalgebra, Lecture 15: Equations for Deterministic Automata

Farey Fractions. Rickard Fernström. U.U.D.M. Project Report 2017:24. Department of Mathematics Uppsala University

Chapter 0. What is the Lebesgue integral about?

Lecture 3: Equivalence Relations

Chapter 14. Matrix Representations of Linear Transformations

Improper Integrals. The First Fundamental Theorem of Calculus, as we ve discussed in class, goes as follows:

Intermediate Math Circles Wednesday, November 14, 2018 Finite Automata II. Nickolas Rollick a b b. a b 4

The Regulated and Riemann Integrals

Determinants Chapter 3

2.4 Linear Inequalities and Interval Notation

Linear Systems with Constant Coefficients

Bridging the gap: GCSE AS Level

Parse trees, ambiguity, and Chomsky normal form

378 Relations Solutions for Chapter 16. Section 16.1 Exercises. 3. Let A = {0,1,2,3,4,5}. Write out the relation R that expresses on A.

The Algebra (al-jabr) of Matrices

p-adic Egyptian Fractions

Here we study square linear systems and properties of their coefficient matrices as they relate to the solution set of the linear system.

SCHOOL OF ENGINEERING & BUILT ENVIRONMENT. Mathematics

308K. 1 Section 3.2. Zelaya Eufemia. 1. Example 1: Multiplication of Matrices: X Y Z R S R S X Y Z. By associativity we have to choices:

SCHOOL OF ENGINEERING & BUILT ENVIRONMENT

20 MATHEMATICS POLYNOMIALS

Rudimentary Matrix Algebra

Discrete Mathematics and Probability Theory Spring 2013 Anant Sahai Lecture 17

Math 4310 Solutions to homework 1 Due 9/1/16

MTH 505: Number Theory Spring 2017

CS 301. Lecture 04 Regular Expressions. Stephen Checkoway. January 29, 2018

5. (±±) Λ = fw j w is string of even lengthg [ 00 = f11,00g 7. (11 [ 00)± Λ = fw j w egins with either 11 or 00g 8. (0 [ ffl)1 Λ = 01 Λ [ 1 Λ 9.

Linear Algebra 1A - solutions of ex.4

1B40 Practical Skills

HW3, Math 307. CSUF. Spring 2007.

Lecture 7 notes Nodal Analysis

I1 = I2 I1 = I2 + I3 I1 + I2 = I3 + I4 I 3

Review of Gaussian Quadrature method

Linear Inequalities. Work Sheet 1

TOPIC: LINEAR ALGEBRA MATRICES

AUTOMATA AND LANGUAGES. Definition 1.5: Finite Automaton

1 Nondeterministic Finite Automata

Unit #9 : Definite Integral Properties; Fundamental Theorem of Calculus

Bernoulli Numbers Jeff Morton

fractions Let s Learn to

DEFINITION The inner product of two functions f 1 and f 2 on an interval [a, b] is the number. ( f 1, f 2 ) b DEFINITION 11.1.

CS 373, Spring Solutions to Mock midterm 1 (Based on first midterm in CS 273, Fall 2008.)

ARITHMETIC OPERATIONS. The real numbers have the following properties: a b c ab ac

Geometric Sequences. Geometric Sequence a sequence whose consecutive terms have a common ratio.

Section 6.1 INTRO to LAPLACE TRANSFORMS

Chapters Five Notes SN AA U1C5

Chapter 6 Techniques of Integration

Matrices. Elementary Matrix Theory. Definition of a Matrix. Matrix Elements:

M344 - ADVANCED ENGINEERING MATHEMATICS

Lecture 3. Introduction digital logic. Notes. Notes. Notes. Representations. February Bern University of Applied Sciences.

MAA 4212 Improper Integrals

Lecture 09: Myhill-Nerode Theorem

1 Online Learning and Regret Minimization

ECON 331 Lecture Notes: Ch 4 and Ch 5

The Wave Equation I. MA 436 Kurt Bryan

MATH 573 FINAL EXAM. May 30, 2007

13: Diffusion in 2 Energy Groups

and that at t = 0 the object is at position 5. Find the position of the object at t = 2.

If we have a function f(x) which is well-defined for some a x b, its integral over those two values is defined as

The Dirichlet Problem in a Two Dimensional Rectangle. Section 13.5

CMPSCI 250: Introduction to Computation. Lecture #31: What DFA s Can and Can t Do David Mix Barrington 9 April 2014

Convert the NFA into DFA

Matrices and Determinants

7.2 The Definite Integral

8. Complex Numbers. We can combine the real numbers with this new imaginary number to form the complex numbers.

Lecture Note 9: Orthogonal Reduction

221B Lecture Notes WKB Method

Operations with Matrices

1 PYTHAGORAS THEOREM 1. Given a right angled triangle, the square of the hypotenuse is equal to the sum of the squares of the other two sides.

State space systems analysis (continued) Stability. A. Definitions A system is said to be Asymptotically Stable (AS) when it satisfies

Homework Solution - Set 5 Due: Friday 10/03/08

Lecture 6. Notes. Notes. Notes. Representations Z A B and A B R. BTE Electronics Fundamentals August Bern University of Applied Sciences

Theoretical foundations of Gaussian quadrature

CM10196 Topic 4: Functions and Relations

QUADRATIC EQUATIONS OBJECTIVE PROBLEMS

Designing Information Devices and Systems I Spring 2018 Homework 7

Week 10: Line Integrals

Lecture 2: January 27

UNIT 5 QUADRATIC FUNCTIONS Lesson 3: Creating Quadratic Equations in Two or More Variables Instruction

Review of Calculus, cont d

CS 311 Homework 3 due 16:30, Thursday, 14 th October 2010

Matrix Solution to Linear Equations and Markov Chains

5.7 Improper Integrals

Transcription:

Lecture 3 3 Solving liner equtions In this lecture we will discuss lgorithms for solving systems of liner equtions Multiplictive identity Let us restrict ourselves to considering squre mtrices since one cn form the product of two squre mtrices in either order nd e le to compre the results In this sitution there is specil mtrix 1 0 0 0 I n := ; 0 0 1 we will sometimes lso write dig(1 1 for the resulting digonl mtrix Computing directly using the definition of mtrix multipliction one just checks tht for ny n n squre mtrix A we hve AI n = I n A = A Definition 31 An n n mtrix B is multiplictive identity (for rel squre mtrices if for ny squre mtrix A the identity AB = BA = A is stisfied Question 32 Are there other multiplictive identity mtrices? Multiplictive inverses All of wht we sid ove ws well nd good ut very forml Mtrices were invented with the gol of solving systems of liner equtions so let s get ck to tht The simplest liner eqution is of the form x = How do we solve this eqution? Assuming 0 we just divide y nd get x = If we hd system of liner equtions of the form Ax = B where A is n m n-mtrix x is n n 1 column vector nd B is n m 1-column vector then it might e nice to find mtrix tht plys the role of 1 Assume m = n so oth sides of the eqution cn e esily compred We would like to find mtrix C such tht CA = I n If we could do tht then we d hve C(Ax = CB 1

2 3 Solving liner equtions nd thus x = CB Writing our equtions like this ws just n ritrry choice: we could just s well hve written our eqution s the product of row vector with mtrix (on the right In this reversed order we would wnt to multiply mtrices on the right We cn now turn this drem into property Definition 33 An multiplictive inverse for n n n-mtrix A is mtrix B such tht AB = BA = I n If A hs n inverse we sy it is invertile We write A 1 for multiplictive inverse of squre mtrix A when it exists However so fr this is just drem We don t know if given mtrix A ctully hs multiplictive inverse In fct mny mtrices do not (consider the mtrix 0 n n Even if we know strctly tht n inverse exists we do know how to compute the inverse! Nevertheless given mtrix A we cn recognize the inverse if it is presented to us y explicit mtrix multipliction Question 34 If CA = I n wht cn we sy out AC? Question 35 If A is invertile cn it hve more thn one inverse? Question 36 Under wht conditions re ech of the following mtrices invertile nd wht re their inverses when they exist? ( 0 0 ( 0 1 0 1 1 0 In clss we worked this out in detil since it will e of some importnce lter Row opertions Given system of liner equtions sy x 1 + x 2 = 21 x 1 + 22 x 2 = 2 we would like to if possile solve the system for x 1 nd x 2 The generl procedure is s follows: first we solve the first eqution for one vrile in terms of the others To e s systemtic s possile: if is not zero we cn divide the first eqution y without chnging its solutions In other words we cn write x 1 + x 2 = nd our old system of equtions ecomes the new system of equtions: x 1 + x 2 = 21 x 1 + 22 x 2 = 2 Let us now keep trck of this mnipultion in terms of coefficient mtrices

3 3 Solving liner equtions We chnge the originl coefficient mtrix ( 11 21 The solution vector chnges similrly: ( 1 2 21 ( 1 2 There is third description of this trnsition: it corresponds to multiplying the originl coefficient mtrix nd the solution mtrix on the left y the digonl mtrix 1 Let us lso oserve tht since 0 this mtrix is invertile nd this encodes the fct tht we could return to the originl system of equtions from our new one y multiplying y the inverse Next we cn solve the new first eqution for x 1 in terms of x 2 ie we cn write x 1 = x 2 Hving performed this opertion we cn sustitute the given expression for x 1 into the other expression In this cse: 21 ( x 2 + 22 x 2 = 2 After some reorgnizing we cn write this s: ( 22 21 x 2 = 2 21 In other words we hve otined the system of equtions x 1 + x 2 = ( 22 21 x 2 = 2 21 Agin let us keep trck of the chnge in coefficient mtrices The new coefficient mtrix for the originl eqution nd the new eqution involving only x 2 is encoded in trnsition of the form 21 0 22 21 while the trnsition for the solution vector now tkes the form: ( 1 ( 2 2 1 21

4 3 Solving liner equtions String t this we hve dded multiple ( 21 to the second row This trnsition lso cn e gin encoded in terms of mtrix products: oserve tht ( ( 1 = 21 1 21 0 22 21 nd similrly ( 0 1 ( = 21 1 2 2 1 21 ( 1 0 To sy it differently multipliction y the mtrix hs the the following effect 21 1 on equtions: we dd multiple of the first eqution to the second (in this cse the multiple is 21 Agin the invertiility of this mtrix encodes the fct tht the new system of equtions is completely equivlent to the originl system of equtions At the end of this step we lmost hve the vlue of x 2 Indeed we hve the eqution ( 22 21 x 2 = 2 21 If the expression ( 22 21 is non-zero we cn divide y it nd get x 2 = 2 21 ( 22 21 As sidenote let us oserve tht so fr we hve ssumed 0 nd ( 22 21 0 Since the product of two non-zero numers is gin non-zero we lso see tht ( 22 21 0 As in the first step this cn e encoded in the following trnsition of the coefficient mtrix: 0 22 21 nd the corresponding trnsition for the solution vector: ( 2 1 21 2 21 ( 22 21 While this looks cumersome gin this trnsition cn e effected y multiplying y digonl mtrix of the form dig(1 ( 22 21 1 on the left Under the hypothesis tht ( 22 21 0 this mtrix is lso invertile To finish we hve to solve for x 1 Since now we know the vlue of x 2 we could plug this ck into the eqution for x 1 In other words we tke the eqution x 1 = ( 2 21 ( 22 21

5 3 Solving liner equtions The coefficient mtrix hs undergone the trnsition: ( 1 0 while the solution vector hs undergone the trnsition 11 2 21 ( 22 21 ( 2 1 21 ( 22 21 2 21 ( 22 21 While this looks very complicted it is esier to how to effect this chnge: we multiply oth the coefficient nd solution vectors on the left y the mtrix ( Once more this mtrix is invertile Along the wy we mde vrious ssumptions out certin things eing non-zero Under these hypotheses we cn keep trck of ll the mtrix multiplictions we hve performed on coefficient mtrices: ( ( ( ( 1 0 1 ( 0 ( 22 21 1 1 = 21 1 21 ( 1 0 Grouping the four mtrices on the left together oserve tht this shows tht our originl coefficient mtrix ws invertile under the two hypotheses: 0 nd 22 21 0 Ws the condition 0 relly necessry? Suppose = 0 In tht cse we cn simply chnge the order of the equtions which of course does not lter the solutions! After tht we get coefficient mtrix tht hs 0 in the lower-left hnd coefficient nd this looks similr to the plce we were prtwy through our nlysis Agin this trnsition in coefficient mtrices cn e effected y left multipliction y mtrix nmely left multipliction y the mtrix ( 1 0 hs the effect of swpping the two rows The fct tht this mtrix is its own inverse corresponds to the fct tht if we swp two rows nd then swp them ck we hve done nothing to the originl system of equtions Next time we will finish the nlysis in this cse