1 Row elimination and row-equivalent matrices.

Transcription

1 Geometria Lingotto LeLing: Echelon form and solutions C ontents: Row elimination and row-equivalent matrices Echelon form and solved systems Gauß-Jordan elimination and echelon form The general solution The theorem of Steinitz R ecommended exercises: GeoLing 1, 1 Row elimination and row-equivalent matrices Two matrices A, B are said row-equivalent when either is obtained from the other by successive EROs Theorem 11 Two homogeneous linear systems are equivalent if and only if (iff) the associated matrices are row-equivalent Therefore, to decide whether two systems are equivalent we should check if their matrices are equivalent under row operations For this we need the so-called echelon form 11 Echelon form and solved systems The basic idea to keep in mind is that a matrix is in echelon form if it represents a system which is essentially already solved In other words we needn t proceed further because the solution is self-evident Definition 1 A matrix A is in echelon form (or is an echelon matrix) when: (i) The zero rows (if present) are grouped together below the nonzero ones (ii) In each nonzero row the first element different from zero is 1, which is also the only nonzero element in its own column (ie above and below the 1 only 0s appear) The number 1 and the column containing it are called special Ingegneria dell Autoveicolo, LeLing 1 Geometry

2 1 Gauß-Jordan elimination and echelon form Geometria Lingotto (iii) The special element of a row lies at the left of any special element below it The matrices ( On the other hand The matrix Also the matrix ), , is not an echelon matrix for (i) fails Here is an important example of an echelon matrix: Example 13 The matrix is not echelon because (ii) does not hold 0 0 are echelon is not echelon, in fact (iii) is not satisfied is echelon 1 Gauß-Jordan elimination and echelon form The Gauß-Jordan elimination method is an algorithm that takes an arbitrary matrix A and uses EROs to produce and echelon matrix E equivalent to A In order to remember the procedure easily, we divide it into two steps, the Gauß step and the Jordan step: Gauß step: it consists essentially in obtaining zeroes below the elements equal to 1 Starting from the first column of A we identify (moving from left to right) the first nonzero column 1 We find in the first nonzero column the first (from the top) nonzero 1 Usually it is the first column to be already nonzero, ie containing a nonzero element a i 1 0; if not, the unknown x 1 is not explicitly involved in the system to be solved Ingegneria dell Autoveicolo, LeLing Geometry

3 1 Gauß-Jordan elimination and echelon form Geometria Lingotto element a i0 j 0 0 Then we use 0 0 a i0 j 0 R i0 a i0 j 0 = a i0 j 0 a i0 j 0 a i0 j 0 From the 1 on row i 0 we obtain 0s below it using the ERO R i + cr i0 This gives a matrix a i0 j 0 a i0 j 0 a i0 j 0 Now we swapp the i 0 row with the first row by using the ERO R i0 1 to get: a i0 j 0 a i0 j 0 a i0 j b b 1 Repeating the procedure from the start on the (sub-)matrix: Ingegneria dell Autoveicolo, LeLing 3 Geometry

4 1 Gauß-Jordan elimination and echelon form Geometria Lingotto b b 1 will lead eventually to a matrix where: (i) The zero rows (if present) are all written below the nonzero ones If not, Gauß s method could be used from the position of the last column (ii) In each nonzero row the first element different from zero is 1; below it only 0s are present (iii) The first 1 of a row lies at the left of the 1 s below it This ends the Gauß step and is where Jordan begins: Jordan step: This works like a backwards gear: starting from the 1 in the last nonzero row we obtain zeroes above it, this time, using the ERO R i + cr j : = The last column is now done with, and the 1 has now become the only nonzero element in the whole column, precisely as required by condition (iii) in the definition of echelon matrix Once the engine is started, we repeat over and over until we get to the 1 of the fist column, the one in position i 0, j 0 This concludes the Jordan step and yields an echelon matrix Ingegneria dell Autoveicolo, LeLing 4 Geometry

5 13 The general solution Geometria Lingotto 13 The general solution When the matrix is echelon the system is practically already solved In fact, the unknowns corresponding to the special columns (those containing just a 1) are variables that depend upon the other unknonws, in turn called independent variables, or parameters Hence x 1, x,, x n can be divided into two groups, one of which is written in terms of the other For example: Example 14 Let S be Its matrix reads: S = x 1 6x + 4 = 0 x + x 3 = 0 3x 1 + 7x = 0 3x 1 + x + x 3 = 0 The Gauß-Jordan method produces the echelon matrix: x 1 + 7x 4 = 0 so the system S is equivalent to: x + 3 x 4 = 0 x 3 + 3x 4 = 0 It is then clear that x 1, x, x 3 depend on The solutions of S read: Note all solutions are multiples of the single column , Example 15 Let S be the system: Ingegneria dell Autoveicolo, LeLing 5 Geometry

6 Geometria Lingotto S = { x 1 + 4x = 0 x 1 + x x 3 + = 0 ( ) with matrix: 1 ( Gauß-Jordan yields the echelon matrix: 5 { x1 + 4x so that S is equivalent to: = 0 x 3 + 5x 4 = 0 Now x 1, x 3 are functions of x, The system s solutions are : 4x + x 5 The solutions are linear combinations of two columns 4x + 4 x 5 x = x The theorem of Steinitz ), Recall that nontrivial solution of a system means nonzero Gauß-Jordan follows this result of Steinitz From the elimination of Theorem 1 A homogeneous system of linear equations with more unknowns than equations admits always a non-trivial solution Ingegneria dell Autoveicolo, LeLing 6 Geometry