Chapter System of Equations

Similar documents
Statistics for Financial Engineering Session 1: Linear Algebra Review March 18 th, 2006

Lesson 4 Linear Algebra

Handout #2. Introduction to Matrix: Matrix operations & Geometric meaning

Vectors. Vectors in Plane ( 2

Autar Kaw Benjamin Rigsby. Transforming Numerical Methods Education for STEM Undergraduates

Section 7.3, Systems of Linear Algebraic Equations; Linear Independence, Eigenvalues, Eigenvectors (the variable vector of the system) and

The total number of permutations of S is n!. We denote the set of all permutations of S by

Chapter 2 Infinite Series Page 1 of 9

Lecture 2: Matrix Algebra

DETERMINANT. = 0. The expression a 1. is called a determinant of the second order, and is denoted by : y + c 1

PROGRESSIONS AND SERIES

Elementary Linear Algebra

lecture 16: Introduction to Least Squares Approximation

For all Engineering Entrance Examinations held across India. Mathematics

Unit 1. Extending the Number System. 2 Jordan School District

Taylor Polynomials. The Tangent Line. (a, f (a)) and has the same slope as the curve y = f (x) at that point. It is the best

Notes 17 Sturm-Liouville Theory

1. (25 points) Use the limit definition of the definite integral and the sum formulas to compute. [1 x + x2

M3P14 EXAMPLE SHEET 1 SOLUTIONS

INFINITE SERIES. ,... having infinite number of terms is called infinite sequence and its indicated sum, i.e., a 1

Matrix Algebra Notes

Section 6.3: Geometric Sequences

( a n ) converges or diverges.

MTH 146 Class 16 Notes

Surds, Indices, and Logarithms Radical

MATRIX ALGEBRA, Systems Linear Equations

SOLUTION OF SYSTEM OF LINEAR EQUATIONS. Lecture 4: (a) Jacobi's method. method (general). (b) Gauss Seidel method.

Fast Fourier Transform 1) Legendre s Interpolation 2) Vandermonde Matrix 3) Roots of Unity 4) Polynomial Evaluation

Sequence and Series of Functions

ALGEBRA. Set of Equations. have no solution 1 b1. Dependent system has infinitely many solutions

Chapter 7 Infinite Series

=> PARALLEL INTERCONNECTION. Basic Properties LTI Systems. The Commutative Property. Convolution. The Commutative Property. The Distributive Property

0 otherwise. sin( nx)sin( kx) 0 otherwise. cos( nx) sin( kx) dx 0 for all integers n, k.

Schrödinger Equation Via Laplace-Beltrami Operator

ECE 102 Engineering Computation

Basic Maths. Fiorella Sgallari University of Bologna, Italy Faculty of Engineering Department of Mathematics - CIRAM

A GENERAL METHOD FOR SOLVING ORDINARY DIFFERENTIAL EQUATIONS: THE FROBENIUS (OR SERIES) METHOD

[ 20 ] 1. Inequality exists only between two real numbers (not complex numbers). 2. If a be any real number then one and only one of there hold.

In an algebraic expression of the form (1), like terms are terms with the same power of the variables (in this case

Review of Sections

1.3 Continuous Functions and Riemann Sums

Closed Newton-Cotes Integration

SUTCLIFFE S NOTES: CALCULUS 2 SWOKOWSKI S CHAPTER 11

General properties of definite integrals

Eigenfunction Expansion. For a given function on the internal a x b the eigenfunction expansion of f(x):

MA123, Chapter 9: Computing some integrals (pp )

Math 3B Midterm Review

Week 13 Notes: 1) Riemann Sum. Aim: Compute Area Under a Graph. Suppose we want to find out the area of a graph, like the one on the right:

INTEGRATION TECHNIQUES (TRIG, LOG, EXP FUNCTIONS)

CHAPTER 2d. MATRICES

Inverse Matrix. A meaning that matrix B is an inverse of matrix A.

M.A. (ECONOMICS) PART-I PAPER - III BASIC QUANTITATIVE METHODS

10.5 Power Series. In this section, we are going to start talking about power series. A power series is a series of the form

2.1.1 Definition The Z-transform of a sequence x [n] is simply defined as (2.1) X re x k re x k r

Chapter Vectors

Chapter Real Numbers

GRAPHING LINEAR EQUATIONS. Linear Equations. x l ( 3,1 ) _x-axis. Origin ( 0, 0 ) Slope = change in y change in x. Equation for l 1.

Approximations of Definite Integrals

Addendum. Addendum. Vector Review. Department of Computer Science and Engineering 1-1

Dynamics of Structures

( ) k ( ) 1 T n 1 x = xk. Geometric series obtained directly from the definition. = 1 1 x. See also Scalars 9.1 ADV-1: lim n.

RADICALS. Upon completion, you should be able to. define the principal root of numbers. simplify radicals

Review of the Riemann Integral

is an ordered list of numbers. Each number in a sequence is a term of a sequence. n-1 term

Numerical Methods (CENG 2002) CHAPTER -III LINEAR ALGEBRAIC EQUATIONS. In this chapter, we will deal with the case of determining the values of x 1

BC Calculus Review Sheet

Discrete Mathematics I Tutorial 12

FREE Download Study Package from website: &

Fig. 1. I a. V ag I c. I n. V cg. Z n Z Y. I b. V bg

Inner Product Spaces (Chapter 5)

RULES FOR MANIPULATING SURDS b. This is the addition law of surds with the same radicals. (ii)

Student Success Center Elementary Algebra Study Guide for the ACCUPLACER (CPT)

Remarks: (a) The Dirac delta is the function zero on the domain R {0}.

Interpolation. 1. What is interpolation?

BRILLIANT PUBLIC SCHOOL, SITAMARHI (Affiliated up to +2 level to C.B.S.E., New Delhi)

Definite Integral. The Left and Right Sums

SM2H. Unit 2 Polynomials, Exponents, Radicals & Complex Numbers Notes. 3.1 Number Theory

A GENERALIZATION OF GAUSS THEOREM ON QUADRATIC FORMS

Numbers (Part I) -- Solutions

EVALUATING DEFINITE INTEGRALS

is continuous at x 2 and g(x) 2. Oil spilled from a ruptured tanker spreads in a circle whose area increases at a

Chapter Real Numbers

We will begin by supplying the proof to (a).

Infinite Series Sequences: terms nth term Listing Terms of a Sequence 2 n recursively defined n+1 Pattern Recognition for Sequences Ex:

Exponents and Radical

Crushed Notes on MATH132: Calculus

Numerical Solution of Fuzzy Fredholm Integral Equations of the Second Kind using Bernstein Polynomials

Prior distributions. July 29, 2002

Physics of Semiconductor Devices Vol.10

Solving Systems of Equations

 n. A Very Interesting Example + + = d. + x3. + 5x4. math 131 power series, part ii 7. One of the first power series we examined was. 2!

(3) If you replace row i of A by its sum with a multiple of another row, then the determinant is unchanged! Expand across the i th row:

Zeros of Polynomials

Linear regression. Daniel Hsu (COMS 4771) (y i x T i β)2 2πσ. 2 2σ 2. 1 n. (x T i β y i ) 2. 1 ˆβ arg min. β R n d

Matrix Algebra 2.2 THE INVERSE OF A MATRIX Pearson Education, Inc.

NICK DUFRESNE. 1 1 p(x). To determine some formulas for the generating function of the Schröder numbers, r(x) = a(x) =

POWER SERIES R. E. SHOWALTER

Frequency-domain Characteristics of Discrete-time LTI Systems

Similar idea to multiplication in N, C. Divide and conquer approach provides unexpected improvements. Naïve matrix multiplication

Theorem: Let A n n. In this case that A does reduce to I, we search for A 1 as the solution matrix X to the matrix equation A X = I i.e.

Transcription:

hpter 4.5 System of Equtios After redig th chpter, you should be ble to:. setup simulteous lier equtios i mtrix form d vice-vers,. uderstd the cocept of the iverse of mtrix, 3. kow the differece betwee costet d icostet system of lier equtios, d 4. ler tht system of lier equtios c hve uique solutio, o solutio or ifiite solutios. Mtrix lgebr used for solvig systems of equtios. you illustrte th cocept? Mtrix lgebr used to solve system of simulteous lier equtios. I fct, for my mthemticl procedures such s the solutio to set of olier equtios, iterpoltio, itegrtio, d differetil equtios, the solutios reduce to set of simulteous lier equtios. Let us illustrte with exmple for iterpoltio. Exmple The upwrd velocity of rocket give t three differet times o the followig tble. Tble 5.. Velocity vs. time dt for rocket Time, t Velocity, v (s) (m/s) 5 6.8 8 77. 79. The velocity dt pproximted by polyomil s v t t bt c, 5 t. Set up the equtios i mtrix form to fid the coefficiets, b, c of the velocity profile. Solutio The polyomil goig through three dt poits, v, t, v, d t v tble 5.. t 5, v t, v 8 6.8 77. t where from 3, 3 4.5.

4.5. hpter 4.5 t 3, v 3 79. Requirig tht v t t bt c psses through the three dt poits gives vt v t bt c vt v t bt c vt 3 v3 t3 bt3 c Substitutig the dt, v, t, v, d t, v 5 b 5 c 6. 8 8 b 8 c 77. b c 79. t gives 3 3 or 5 5b c 6.8 64 8b c 77. 44 b c 79. Th set of equtios c be rewritte i the mtrix form s 5 5b c 6.8 64 8b c 77. 44 b c 79. The bove equtio c be writte s lier combitio s follows 5 5 6.8 64 b 8 c 77. 44 79. d further usig mtrix multiplictio gives 5 5 6.8 64 8 b 77. 44 c 79. The bove illustrtio of why mtrix lgebr eeded. The complete solutio to the set of equtios give lter i th chpter. A geerl set of m lier equtios d ukows, x x x c x x x c. m x m x... m x cm c be rewritte i the mtrix form s

System of Equtios 4.5.3.. x c.. x c m m.. m x cm Deotig the mtrices by A, X, d, the system of equtio AX, where A clled the coefficiet mtrix, clled the right hd side vector d X clled the solutio vector. A X systems of equtios re writte i the ugmeted form. Tht Sometimes... c... c A m m... m c A system of equtios c be costet or icostet. Wht does tht me? A system of equtios AX costet if there solutio, d it icostet if there o solutio. However, costet system of equtios does ot me uique solutio, tht, costet system of equtios my hve uique solutio or ifiite solutios (Figure ). [A][X]= [B] ostet System Icostet System Uique Solutio Ifiite Solutios Figure 5.. ostet d icostet system of equtios flow chrt. Exmple Give exmples of costet d icostet system of equtios. Solutio ) The system of equtios

4.5.4 hpter 4.5 4x 6 3y 4 costet system of equtios s it hs uique solutio, tht, x. y b) The system of equtios 4x 6 y 3 lso costet system of equtios but it hs ifiite solutios s give s follows. Expdig the bove set of equtios, x 4y 6 x y 3 you c see tht they re the sme equtio. Hece, y combitio of x, y tht stfies x 4y 6 solutio. For exmple x, y, solutio. Other solutios iclude x, y (.5,.5), x, y (,.5), d so o. c) The system of equtios 4x 6 y 4 icostet s o solutio exts. How c oe dtiguh betwee costet d icostet system of equtios? A system of equtios X ugmeted mtrix A A system of equtios X the ugmeted mtrix A. A costet if the rk of A equl to the rk of the A icostet if the rk of A less th the rk of But, wht do you me by rk of mtrix? The rk of mtrix defied s the order of the lrgest squre submtrix whose determit ot zero. Exmple 3 Wht the rk of 3 A 5? 3 Solutio The lrgest squre submtrix possible of order 3 d tht [A] itself. Sice det( A ) 3, the rk of [ A ] 3.

System of Equtios 4.5.5 Exmple 4 Wht the rk of 3 A 5? 5 7 Solutio The lrgest squre submtrix of [A] of order 3 d tht [A] itself. Sice det( A ), the rk of [A] less th 3. The ext lrgest squre submtrix would be mtrix. Oe of the squre submtrices of [A] 3 B d det( B ). Hece the rk of [A]. There o eed to look t other submtrices to estblh tht the rk of [A]. Exmple 5 How do I ow use the cocept of rk to fid if 5 5 x 6.8 64 8 x 77. 44 x 79. 3 costet or icostet system of equtios? Solutio The coefficiet mtrix 5 5 A 64 8 44 d the right hd side vector 6.8 77. 79. The ugmeted mtrix 5 5 6.8 B 64 8 77. 44 79. Sice there re o squre submtrices of order 4 s [B] 34 mtrix, the rk of [B] t most 3. So let us look t the squre submtrices of [B] of order 3; if y of these squre submtrices hve determit ot equl to zero, the the rk 3. For exmple, submtrix of the ugmeted mtrix [B]

4.5.6 hpter 4.5 5 5 [D] 64 8 44 hs det( D ) 84. Hece the rk of the ugmeted mtrix [B] 3. Sice [ A] [ D], the rk of [A] 3. Sice the rk of the ugmeted mtrix [B] equls the rk of the coefficiet mtrix [A], the system of equtios costet. Exmple 6 Use the cocept of rk of mtrix to fid if 5 5 x 6.8 64 8 x 77. 89 3 x 84. 3 costet or icostet? Solutio The coefficiet mtrix give by 5 5 A 64 8 89 3 d the right hd side 6.8 77. 84. The ugmeted mtrix 5 5 :6.8 B 64 8 :77. 89 3 : 84. Sice there re o squre submtrices of order 4 s [B] 43 mtrix, the rk of the ugmeted [B] t most 3. So let us look t squre submtrices of the ugmeted mtrix [B] of order 3 d see if y of these hve determits ot equl to zero. For exmple, squre submtrix of the ugmeted mtrix [B] 5 5 D 64 8 89 3 hs det( D ). Th mes, we eed to explore other squre submtrices of order 3 of the ugmeted mtrix [B] d fid their determits. Tht,

System of Equtios 4.5.7 5 E 8 3 det( E ) 5 F 64 89 det( F ) 6.8 77. 84. 5 6.8 8 77. 3 84. 5 6.8 G 64 77. 89 84. det( G ) All the squre submtrices of order 33 of the ugmeted mtrix [B] hve zero determit. So the rk of the ugmeted mtrix [B] less th 3. Is the rk of ugmeted mtrix [B] equl to?. Oe of the submtrices of the ugmeted mtrix [B] 5 5 H 64 8 d det( H ) So the rk of the ugmeted mtrix [B]. Now we eed to fid the rk of the coefficiet mtrix [B]. 5 5 A 64 8 89 3 d det( A ) So the rk of the coefficiet mtrix [A] less th 3. A squre submtrix of the coefficiet mtrix [A] 5 J 8 det( J ) 3 So the rk of the coefficiet mtrix [A]. Hece, rk of the coefficiet mtrix [A] equls the rk of the ugmeted mtrix [B]. So the system of equtios [ A][ X ] [ ] costet.

4.5.8 hpter 4.5 Exmple 7 Use the cocept of rk to fid if 5 5 x 6.8 64 8 x 77. 89 3 x3 8. costet or icostet. Solutio The ugmeted mtrix 5 5 :6.8 B 64 8 :77. 89 3 : 8. Sice there re o squre submtrices of order 44 s the ugmeted mtrix [B] 43 mtrix, the rk of the ugmeted mtrix [B] t most 3. So let us look t squre submtrices of the ugmeted mtrix (B) of order 3 d see if y of the 33 submtrices hve determit ot equl to zero. For exmple, squre submtrix of order 33 of [B] 5 5 D 64 8 89 3 det(d) = So it mes, we eed to explore other squre submtrices of the ugmeted mtrix [B] 5 6.8 E 8 77. 3 8. det( E.. So the rk of the ugmeted mtrix [B] 3. The rk of the coefficiet mtrix [A] from the previous exmple. Sice the rk of the coefficiet mtrix [A] less th the rk of the ugmeted mtrix [B], the system of equtios icostet. Hece, o solutio exts for [ A][ X ] [ ]. If solutio exts, how do we kow whether it uique? I system of equtios [ A][ X ] [ ] tht costet, the rk of the coefficiet mtrix [A] the sme s the ugmeted mtrix [ A ]. If i dditio, the rk of the coefficiet mtrix [A] sme s the umber of ukows, the the solutio uique; if the rk of the coefficiet mtrix [A] less th the umber of ukows, the ifiite solutios ext.

System of Equtios 4.5.9 [A] [X] = [B] ostet System if rk (A) = rk (A.B) Icostet System if rk (A) < rk (A.B) Uique solutio if rk (A) = umber of ukows Ifiite solutios if rk (A) < umber of ukows Figure 5.. Flow chrt of coditios for costet d icostet system of equtios. Exmple 8 We foud tht the followig system of equtios 5 5 x 6.8 64 8 x 77. 44 x 79. 3 costet system of equtios. Does the system of equtios hve uique solutio or does it hve ifiite solutios? Solutio The coefficiet mtrix 5 5 A 64 8 44 d the right hd side 6.8 77. 79. While fidig out whether the bove equtios were costet i erlier exmple, we foud tht the rk of the coefficiet mtrix (A) equls rk of ugmeted mtrix A equls 3. The solutio uique s the umber of ukows = 3 = rk of (A). Exmple 9 We foud tht the followig system of equtios 5 5 x 6.8 64 8 x 77. 89 3 x 84. 3 costet system of equtios. Is the solutio uique or does it hve ifiite solutios.

4.5. hpter 4.5 Solutio While fidig out whether the bove equtios were costet, we foud tht the rk of the coefficiet mtrix [A] equls the rk of ugmeted mtrix A equls Sice the rk of [ A ] < umber of ukows = 3, ifiite solutios ext. If we hve more equtios th ukows i [A] [X] = [], does it me the system icostet? No, it depeds o the rk of the ugmeted mtrix A d the rk of [A]. ) For exmple 5 5 6.8 x 64 8 77. x 44 79. x 3 89 3 84. costet, sice rk of ugmeted mtrix = 3 rk of coefficiet mtrix = 3 Now sice the rk of (A) = 3 = umber of ukows, the solutio ot oly costet but lso uique. b) For exmple 5 5 6.8 x 64 8 77. x 44 79. x 3 89 3 8. icostet, sice rk of ugmeted mtrix = 4 rk of coefficiet mtrix = 3 c) For exmple 5 5 6.8 x 64 8 77. x 5 3.6 x 3 89 3 8. costet, sice rk of ugmeted mtrix = rk of coefficiet mtrix = But sice the rk of [A] = < the umber of ukows = 3, ifiite solutios ext. ostet systems of equtios c oly hve uique solutio or ifiite solutios. system of equtios hve more th oe but ot ifiite umber of solutios? No, you c oly hve either uique solutio or ifiite solutios. Let us suppose [ A][ X ] [ ] hs two solutios [Y ] d [Z] so tht

System of Equtios 4.5. [ A][ Y ] [ ] [ A][ Z] [ ] If r costt, the from the two equtios ray r raz r Addig the bove two equtios gives ray raz r r Ar Y rz Hece ry rz solutio to AX Sice r y sclr, there re ifiite solutios for [ A][ X ] [ ] of the form r Y r Z you divide two mtrices? If [ A][ B] [ ] defied, it might seem ituitive tht B [ A], but mtrix divio ot defied like tht. However iverse of mtrix c be defied for certi types of squre mtrices. The iverse of squre mtrix [A], if extig, deoted by [ A ] such tht [ A][ A] [ I] [ A] [ A] Where [I] the idetity mtrix. I other words, let [A] be squre mtrix. If [B] other squre mtrix of the sme size such tht [ B][ A] [ I], the [B] the iverse of [A]. [A] the clled to be ivertible or osigulr. If [ A ] does ot ext, [A] clled oivertible or sigulr. If [A] d [B] re two mtrices such tht [ B][ A] [ I], the these sttemets re lso true [B] the iverse of [A] [A] the iverse of [B] [A] d [B] re both ivertible [A] [B]=[I]. [A] d [B] re both osigulr ll colums of [A] d [B]re lierly idepedet ll rows of [A] d [B] re lierly idepedet. Exmple Determie if 3 [B] 5 the iverse of 3

4.5. hpter 4.5 3 [A] 5 Solutio 3 3 3 [ B][ A] 5 3 5 3 [I] Sice [ B][ A] [ I], [B] the iverse of [A] d [A] the iverse of [B]. But, we c lso show tht 3 3 [ A][ B] 5 3 5 3 [I] to show tht [A] the iverse of [B]. I use the cocept of the iverse of mtrix to fid the solutio of set of equtios [A] [X] = []? Yes, if the umber of equtios the sme s the umber of ukows, the coefficiet mtrix [A] squre mtrix. Give [ A][ X ] [ ] The, if [ A ] exts, multiplyig both sides by [ A ]. [ A] [ A][ X ] [ A] [ ] [ I][ X ] [ A] [ ] [ X ] [ A] [ ] Th implies tht if we re ble to fid [ A ], the solutio vector of [ A][ X ] [ ] simply multiplictio of [ A ] d the right hd side vector, []. How do I fid the iverse of mtrix? If [A] mtrix, the [ iverse of mtrix [ A][ A] [ I] Deotig A ] mtrix d ccordig to the defiitio of

System of Equtios 4.5.3 A ] [ ] [ A ] [I Usig the defiitio of mtrix multiplictio, the first colum of the ] [ A mtrix c the be foud by solvig Similrly, oe c fid the other colums of the ] [ A mtrix by chgig the right hd side ccordigly. Exmple The upwrd velocity of the rocket give by Tble 5.. Velocity vs time dt for rocket Time, t (s) Velocity, v (m/s) 5 6.8 8 77. 79. I erlier exmple, we wted to pproximte the velocity profile by 5, t c bt t t v We foud tht the coefficiets c b d,, i t v re give by

4.5.4 hpter 4.5 79. 77. 6.8 c b 44 8 64 5 5 First, fid the iverse of 44 8 64 5 5 A d the use the defiitio of iverse to fid the coefficiets.,d, c b Solutio If 33 3 3 3 3 A the iverse of ] [A, the 44 8 64 5 5 33 3 3 3 3 gives three sets of equtios 44 8 64 5 5 3 44 8 64 5 5 3 44 8 64 5 5 33 3 3 Solvig the bove three sets of equtios seprtely gives 3 = 4.57.954.476 3 = 5..47.8333

System of Equtios 4.5.5 Hece Now where [ A ] 3 3 33.357.4643.49.476.954 4.57 AX X b c 6.8 77. 79. Usig the defiitio of A, A AX A X A Hece So.476.954 4.57.95 b 9.69 c.86.8333.47 5..8333.47 5..357.4643.49.357 6.8.4643 77..49 79. v t.95t 9.69t.86, 5 t Is there other wy to fid the iverse of mtrix? For fidig the iverse of smll mtrices, the iverse of ivertible mtrix c be foud by A dja deta where

4.5.6 hpter 4.5 where dj A ij re the cofctors of ij T. The mtrix itself clled the mtrix of cofctors from [A]. ofctors re defied i hpter 4. Exmple Fid the iverse of 5 5 A 64 8 44 Solutio From Exmple 4.6 i hpter 4, we foud deta 84 Next we eed to fid the djoit of [A]. The cofctors of A re foud s follows. The mior of etry 5 5 M 64 8 44 8 4 The cofctors of etry M M 4 The mior of etry 5 5 M 64 8 44

System of Equtios 4.5.7 64 44 8 The cofctor of etry M M (8) 8 Similrly 3 384 7 9 3 4 3 3 3 39 33 Hece the mtrix of cofctors of [A] 4 8 384 7 9 4 3 39 T The djoit of mtrix [A] [], dja T 4 7 3 8 9 39 384 4 Hece A dja det A 4 7 3 8 9 39 84 384 4.476.8333.357.954.47.4643 4.57 5..49

4.5.8 hpter 4.5 If the iverse of squre mtrix [A] exts, it uique? Yes, the iverse of squre mtrix uique, if it exts. The proof s follows. Assume tht the iverse of [A] [B] d if th iverse ot uique, the let other iverse of [A] ext clled []. If [B] the iverse of [A], the [ B][ A] [ I] Multiply both sides by [], [ B][ A][ ] [ I][ ] [ B][ A][ ] [ ] Sice [] iverse of [A], [ A][ ] [ I] Multiply both sides by [B], [ B][ I] [ ] [ B] [ ] Th shows tht [B] d [] re the sme. So the iverse of [A] uique. Key Terms: ostet system Icostet system Ifiite solutios Uique solutio Rk Iverse