jacobi’s formula: a -1 = adj (a) where adj (a) i,j = (-1) i+j det (a[ j,i ])
DESCRIPTION
Jacobi’s formula: A -1 = Adj (A) where Adj (A) i,j = (-1) i+j det (A[ j,i ]) and A[ j,i ] is A with row j and column i deleted. Use this to find the inverse of :. Computing the Adjoint matrix:. Assignment #3 is available. Due : Monday Nov. 5, beginning of class. Initial dictionary: - PowerPoint PPT PresentationTRANSCRIPT
![Page 1: Jacobi’s formula: A -1 = Adj (A) where Adj (A) i,j = (-1) i+j det (A[ j,i ])](https://reader035.vdocument.in/reader035/viewer/2022062521/5681666d550346895dda094b/html5/thumbnails/1.jpg)
1
Jacobi’s formula:A-1 = Adj(A) where Adj(A)i,j = (-1)i+j det(A[j,i]) and A[j,i] is A with row j and column i deleted.Use this to find the inverse of:
![Page 2: Jacobi’s formula: A -1 = Adj (A) where Adj (A) i,j = (-1) i+j det (A[ j,i ])](https://reader035.vdocument.in/reader035/viewer/2022062521/5681666d550346895dda094b/html5/thumbnails/2.jpg)
2
Computing the Adjoint matrix:
![Page 3: Jacobi’s formula: A -1 = Adj (A) where Adj (A) i,j = (-1) i+j det (A[ j,i ])](https://reader035.vdocument.in/reader035/viewer/2022062521/5681666d550346895dda094b/html5/thumbnails/3.jpg)
3
Assignment #3 is available.Due: Monday Nov. 5, beginning of class
![Page 4: Jacobi’s formula: A -1 = Adj (A) where Adj (A) i,j = (-1) i+j det (A[ j,i ])](https://reader035.vdocument.in/reader035/viewer/2022062521/5681666d550346895dda094b/html5/thumbnails/4.jpg)
Initial dictionary:X4 =5 -1 X1 -1 X2 +1 X3X5 =3 -1 X1 +0 X2 -2 X3X6 =4 +0 X1 -1 X2 -6 X3------------------------z = 0 +1 X1 +3 X2 +6 X3
Final dictionary:X2 =4 -6 X3 +0 X4 -1 X6 X1 =1 +7 X3 -1 X4 +1 X6 X5 =2 -9 X3 +1 X4 -1 X6 -----------------------z =13 -5 X3 -1 X4 -2 X6
What are B and B-1 for the final dictionary?
4
![Page 5: Jacobi’s formula: A -1 = Adj (A) where Adj (A) i,j = (-1) i+j det (A[ j,i ])](https://reader035.vdocument.in/reader035/viewer/2022062521/5681666d550346895dda094b/html5/thumbnails/5.jpg)
5
The initial matrix A: X1 X2 X3 X4 X5 X6 [ 1 1 -1 1 0 0 ][ 1 0 2 0 1 0 ][ 0 1 6 0 0 1 ]
The final matrix: X1 X2 X3 X4 X5 X6[ 0 1 6 0 0 1 ] [ 1 0 -7 1 0 -1 ] [ 0 0 9 -1 1 1 ]
![Page 6: Jacobi’s formula: A -1 = Adj (A) where Adj (A) i,j = (-1) i+j det (A[ j,i ])](https://reader035.vdocument.in/reader035/viewer/2022062521/5681666d550346895dda094b/html5/thumbnails/6.jpg)
6
The initial matrix A: X1 X2 X3 X4 X5 X6 [ 1 1 -1 1 0 0 ][ 1 0 2 0 1 0 ][ 0 1 6 0 0 1 ]The final matrix: X1 X2 X3 X4 X5 X6[ 0 1 6 0 0 1 ] [ 1 0 -7 1 0 -1 ] [ 0 0 9 -1 1 1 ] B = X2 X1 X5[ 1 1 0 ][ 0 1 1 ][ 1 0 0 ]
B-1 = X4 X5 X6 [ 0 0 1 ][ 1 0 -1 ][ -1 1 1 ]
![Page 7: Jacobi’s formula: A -1 = Adj (A) where Adj (A) i,j = (-1) i+j det (A[ j,i ])](https://reader035.vdocument.in/reader035/viewer/2022062521/5681666d550346895dda094b/html5/thumbnails/7.jpg)
7
Initial dictionary:X4 =5 -1 X1 -1 X2 +1 X3X5 =3 -1 X1 +0 X2 -2 X3X6 =4 +0 X1 -1 X2 -6 X3------------------------z = 0 +1 X1 +3 X2 +6 X3
What of this information do we really need to know to decide where to pivot?
![Page 8: Jacobi’s formula: A -1 = Adj (A) where Adj (A) i,j = (-1) i+j det (A[ j,i ])](https://reader035.vdocument.in/reader035/viewer/2022062521/5681666d550346895dda094b/html5/thumbnails/8.jpg)
8
Initial dictionary:X4 =5 -1 X1 -1 X2 +1 X3X5 =3 -1 X1 +0 X2 -2 X3X6 =4 +0 X1 -1 X2 -6 X3------------------------z = 0 +1 X1 +3 X2 +6 X3
1. Look at z row to choose the pivot column.
With smallest subscript rule, we would choose X1 to enter.
![Page 9: Jacobi’s formula: A -1 = Adj (A) where Adj (A) i,j = (-1) i+j det (A[ j,i ])](https://reader035.vdocument.in/reader035/viewer/2022062521/5681666d550346895dda094b/html5/thumbnails/9.jpg)
9
Initial dictionary:X4 =5 -1 X1 -1 X2 +1 X3X5 =3 -1 X1 +0 X2 -2 X3X6 =4 +0 X1 -1 X2 -6 X3------------------------z = 0 +1 X1 +3 X2 +6 X3
2. Look at entering column and b column to find tightest constraint (pivot row).
![Page 10: Jacobi’s formula: A -1 = Adj (A) where Adj (A) i,j = (-1) i+j det (A[ j,i ])](https://reader035.vdocument.in/reader035/viewer/2022062521/5681666d550346895dda094b/html5/thumbnails/10.jpg)
10
After 1 pivot:
X4 = 2- 1 X2 + 3 X3 + 1 X5 X1 = 3+ 0 X2 - 2 X3 - 1 X5 X6 = 4- 1 X2 - 6 X3 + 0 X5 --------------------------------z = 3+ 3 X2 + 4 X3 - 1 X5
Idea behind the revised Simplex method:
just keep track of what we need to know instead of the entire tableau.
![Page 11: Jacobi’s formula: A -1 = Adj (A) where Adj (A) i,j = (-1) i+j det (A[ j,i ])](https://reader035.vdocument.in/reader035/viewer/2022062521/5681666d550346895dda094b/html5/thumbnails/11.jpg)
11
Note: A standard convention in matrix algebra is that a vector x is by default a column vector and xT is a row vector. The text does not follow this standard convention, but I am using it. Hence, the notation here may differ a bit from the text. x= xT= [x1 x2 x3 ... xn]
![Page 12: Jacobi’s formula: A -1 = Adj (A) where Adj (A) i,j = (-1) i+j det (A[ j,i ])](https://reader035.vdocument.in/reader035/viewer/2022062521/5681666d550346895dda094b/html5/thumbnails/12.jpg)
12
In the Simplex Method, what information is needed to find the next basic feasible solution?
1. Coefficients of non-basic variables in the z row to determine the pivot column.
2. The pivot column and the current solution to determine the pivot row. The new solution is then determined from the pivot column and the old solution.
![Page 13: Jacobi’s formula: A -1 = Adj (A) where Adj (A) i,j = (-1) i+j det (A[ j,i ])](https://reader035.vdocument.in/reader035/viewer/2022062521/5681666d550346895dda094b/html5/thumbnails/13.jpg)
HBT = indices of basic columns.
HNT = indices of non-basic columns.
Split the matrix A into two parts:AB = columns corresponding to HB
T.
AN = columns corresponding to HNT.
Create cB, cN, and xB, xN the same way.
What are these for this initial dictionary?X4 =5 -1 X1 -1 X2 +1 X3X5 =3 -1 X1 +0 X2 -2 X3X6 =4 +0 X1 -1 X2 -6 X3------------------------z = 0 +1 X1 +3 X2 +6 X3
13
![Page 14: Jacobi’s formula: A -1 = Adj (A) where Adj (A) i,j = (-1) i+j det (A[ j,i ])](https://reader035.vdocument.in/reader035/viewer/2022062521/5681666d550346895dda094b/html5/thumbnails/14.jpg)
HBT = indices of basic columns
HNT = indices of non-basic columns
Split the matrix A into two parts:AB = columns corresponding to HB
T
AN = columns corresponding to HNT
Create cB, cN, and xB, xN the same way. What are these for this finaldictionary?
X2 =4 -6 X3 +0 X4 -1 X6 X1 =1 +7 X3 -1 X4 +1 X6 X5 =2 -9 X3 +1 X4 -1 X6 -----------------------z =13 -5 X3 -1 X4 -2 X6
14
![Page 15: Jacobi’s formula: A -1 = Adj (A) where Adj (A) i,j = (-1) i+j det (A[ j,i ])](https://reader035.vdocument.in/reader035/viewer/2022062521/5681666d550346895dda094b/html5/thumbnails/15.jpg)
15
Another problem:X4 = 5 - 2 X1 - 3 X2 - 1 X3 X5 = 11 - 4 X1 - 1 X2 - 2 X3 X6 = 8 - 3 X1 - 4 X2 - 2 X3 ------------------------------z = 0 + 5 X1 + 4 X2 + 3 X3
X1 enters. X4 leaves.
X1 = 2.5 - 1.5 X2 -0.5 X3 - 0.5 X4 X5 = 1.0 + 5.0 X +0.0 X3 + 2.0 X4 X6 = 0.5 + 0.5 X2 -0.5 X3 + 1.5 X4 --------------------------------------z = 12.5 - 3.5 X2+ 0.5 X3 - 2.5 X4
![Page 16: Jacobi’s formula: A -1 = Adj (A) where Adj (A) i,j = (-1) i+j det (A[ j,i ])](https://reader035.vdocument.in/reader035/viewer/2022062521/5681666d550346895dda094b/html5/thumbnails/16.jpg)
16
For this problem we have: Maximize cT xcT= [5 4 3 0 0 0] x=[x1 x2 x3 x4 x5 x6]T
subject to A x = bA = [ 2 3 1 1 0 0 ] [ 5] [ 4 1 2 0 1 0 ] b= [11] [ 3 4 2 0 0 1 ] [ 8]Initially, the basis corresponds to the slack variables x4, x5, x6: HB
T= [ 4 5 6] // subscripts for basisThe non-basic variables are x1, x2, and x3: HN
T= [1 2 3]// subscripts of non-basic variables.
![Page 17: Jacobi’s formula: A -1 = Adj (A) where Adj (A) i,j = (-1) i+j det (A[ j,i ])](https://reader035.vdocument.in/reader035/viewer/2022062521/5681666d550346895dda094b/html5/thumbnails/17.jpg)
17
Split the matrix A into two parts: AB = columns corresponding to HB
T
AN = columns corresponding to HNT
HBT =
[4 5 6]AB =[1 0 0][0 1 0][0 0 1]Create cB, cN, and xB, xN the same way: cB = [0 0 0]T
cN= [5 4 3]T
xB= [x4 x5 x6]T
xN= [x1 x2 x3]T
HNT =
[1 2 3]AN =[2 3 1][4 1 2][3 4 2]
![Page 18: Jacobi’s formula: A -1 = Adj (A) where Adj (A) i,j = (-1) i+j det (A[ j,i ])](https://reader035.vdocument.in/reader035/viewer/2022062521/5681666d550346895dda094b/html5/thumbnails/18.jpg)
18
The system of equations can then be expressed as:
(a) A x = b A⟹ B xB + AN xN = b
(b) z= cBT xB + cN
T xN
From (a), xB = AB-1 b - AB
-1 AN xN
= AB-1 [b - AN xN]
Plugging this into the formula for z:
z= cBT AB
-1 [ b - AN xN ] + cNT xN
So the dictionary with basis AB reads:
xB = AB -1 * [ b - AN * xN ] -------------------------------------------z= cBT * AB -1 * [ b - AN * xN ] + cNT * xN
![Page 19: Jacobi’s formula: A -1 = Adj (A) where Adj (A) i,j = (-1) i+j det (A[ j,i ])](https://reader035.vdocument.in/reader035/viewer/2022062521/5681666d550346895dda094b/html5/thumbnails/19.jpg)
19
From the previous slide:
xB = AB-1 [b - AN xN]
z= cBT AB
-1 [ b - AN xN ] + cNT xN
So the dictionary with basis AB reads:
xB = AB-1 [ b – AN xN ]
-----------------------------z = cB
T AB-1 [ b - AN xN ] + cN
T xN
![Page 20: Jacobi’s formula: A -1 = Adj (A) where Adj (A) i,j = (-1) i+j det (A[ j,i ])](https://reader035.vdocument.in/reader035/viewer/2022062521/5681666d550346895dda094b/html5/thumbnails/20.jpg)
20
The dictionary with basis AB reads: xB = AB
-1 [ b – AN xN ]-----------------------------z = cB
T AB-1 [ b - AN xN ] + cN
T xN
xB= =
x4 = 5 - 2 x1 - 3 x2 - 1 x3 x5 = 11 - 4 x1 - 1 x2 - 2 x3 x6 = 8 - 3 x1 - 4 x2 - 2 x3 ------------------------------z = 0 + 5 x1 + 4 x2 + 3 x3
![Page 21: Jacobi’s formula: A -1 = Adj (A) where Adj (A) i,j = (-1) i+j det (A[ j,i ])](https://reader035.vdocument.in/reader035/viewer/2022062521/5681666d550346895dda094b/html5/thumbnails/21.jpg)
21
z = cBT AB
-1 [ b - AN xN ] + cNT xN
z=
(0 0 0) +
(5 4 3 )
x4 = 5 - 2 x1 - 3 x2 - 1 x3 x5 = 11 - 4 x1 - 1 x2 - 2 x3 x6 = 8 - 3 x1 - 4 x2 - 2 x3 ------------------------------z = 0 + 5 x1 + 4 x2 + 3 x3
![Page 22: Jacobi’s formula: A -1 = Adj (A) where Adj (A) i,j = (-1) i+j det (A[ j,i ])](https://reader035.vdocument.in/reader035/viewer/2022062521/5681666d550346895dda094b/html5/thumbnails/22.jpg)
22
Step 1: Find coefficients of non-basic variables in the z row. z = cB
T AB-1 [ b - AN xN ] + cN
T xN
Let yT= cBT * AB
-1
Find yT by solving yT AB = cBT
or equivalently ABT y = cB.
Set z= yT b + [cNT - yT AN ] xN
![Page 23: Jacobi’s formula: A -1 = Adj (A) where Adj (A) i,j = (-1) i+j det (A[ j,i ])](https://reader035.vdocument.in/reader035/viewer/2022062521/5681666d550346895dda094b/html5/thumbnails/23.jpg)
23
Find yT by solving ABT y = cB.
Set z= yT b + [cNT - yT AN ] xN
=
z = [0 0 0]+ z = 0 + 5 x1 + 4 x2 + 3 x3
The yT b gives the constant term in the z row and the rest gives the terms corresponding to the non-basic variables.
![Page 24: Jacobi’s formula: A -1 = Adj (A) where Adj (A) i,j = (-1) i+j det (A[ j,i ])](https://reader035.vdocument.in/reader035/viewer/2022062521/5681666d550346895dda094b/html5/thumbnails/24.jpg)
Step 2: Determine the leaving variable. Solve for entering column d in current dictionary: d= AB
-1 a where a is the entering column taken from the initial problem.
Or equivalently, solve for d: AB d = a
If tightest constraint corresponds to
xleaving= v – t * xentering
then the new value of xentering will bes= v/t.
24
![Page 25: Jacobi’s formula: A -1 = Adj (A) where Adj (A) i,j = (-1) i+j det (A[ j,i ])](https://reader035.vdocument.in/reader035/viewer/2022062521/5681666d550346895dda094b/html5/thumbnails/25.jpg)
25
Solving AB d = a: Recall: z = 0 + 5 x1 + 4 x2 + 3 x3
Choose x1 to enter. Which variable leaves?
Solving AB d= a:
=
![Page 26: Jacobi’s formula: A -1 = Adj (A) where Adj (A) i,j = (-1) i+j det (A[ j,i ])](https://reader035.vdocument.in/reader035/viewer/2022062521/5681666d550346895dda094b/html5/thumbnails/26.jpg)
26
Which equation imposes the tightest constraint? 5 - 2x1 ≥ 0 ⟹ x1 ≤ 5/2 (*)11 - 4x1 ≥ 0 ⟹ x1 ≤ 11/4 8 - 3x1 ≥ 0 ⟹ x1 ≤ 8/3 The first basis variable x4 leaves because the first equation is the tightest. The value of the entering variable is 5/2 because this is the tightest constraint. Plug in the value of the entering variable: x4= 5 - 2(5/2) = 0 x1= 5/2x5= 11 - 4(5/2) = 1x6= 8 - 3(5/2) = 1/2
![Page 27: Jacobi’s formula: A -1 = Adj (A) where Adj (A) i,j = (-1) i+j det (A[ j,i ])](https://reader035.vdocument.in/reader035/viewer/2022062521/5681666d550346895dda094b/html5/thumbnails/27.jpg)
27
Updating all the variables:HB
T =[1 5 6]AB =[2 0 0][4 1 0][3 0 1] cB = [5 0 0]T
cN= [0 4 3]T
xB= [x1 x5 x6]T
xN= [x4 x2 x3]T
Current solution:[5/2 1 1/2]
HNT =
[4 2 3]AN =[1 3 1][0 1 2][0 4 2]
Z is:(previous value of Z) + (coeff. of entering var. in Z row) * (new value of entering variable)= 0 + 5*(5/2) = 25/2
![Page 28: Jacobi’s formula: A -1 = Adj (A) where Adj (A) i,j = (-1) i+j det (A[ j,i ])](https://reader035.vdocument.in/reader035/viewer/2022062521/5681666d550346895dda094b/html5/thumbnails/28.jpg)
28
The Revised Simplex Method uses more work to determine the
• z row coefficients,
• the column that should enter,
• and the pivot row number (exiting variable)
but then it takes less work to pivot.
![Page 29: Jacobi’s formula: A -1 = Adj (A) where Adj (A) i,j = (-1) i+j det (A[ j,i ])](https://reader035.vdocument.in/reader035/viewer/2022062521/5681666d550346895dda094b/html5/thumbnails/29.jpg)
29
Summary of the steps of the Revised Simplex Method: Maintain for each step: 1. HB- the current basis header. 2. HN- the header for the non-basic variables (optional). 3. the current solution. 4. AB- the columns from A which correspond to the basis in the same order as in HB. [Actually, most programs maintain some factorization of this matrix].
![Page 30: Jacobi’s formula: A -1 = Adj (A) where Adj (A) i,j = (-1) i+j det (A[ j,i ])](https://reader035.vdocument.in/reader035/viewer/2022062521/5681666d550346895dda094b/html5/thumbnails/30.jpg)
30
5. AN- the columns from A which correspond to the basis in the same order as in HN. [Actually, you can get these columns from A when you need them so you don't really have to store this.] 6. cB- the costs of the basic variables in the same order as HB. 7. cN- the costs of the non-basic variables in the same order as HN. 8. z- the current value of the objective function.
![Page 31: Jacobi’s formula: A -1 = Adj (A) where Adj (A) i,j = (-1) i+j det (A[ j,i ])](https://reader035.vdocument.in/reader035/viewer/2022062521/5681666d550346895dda094b/html5/thumbnails/31.jpg)
The Revised Simplex Algorithm Step 1: Determine pivot column.
Solve ABT y = cB for y.
compute [cNT - yT AN] * xN
to get coefficients of non-basic variables.
Look for a positive coefficient, say r corresponding to non-basic xj.
31
![Page 32: Jacobi’s formula: A -1 = Adj (A) where Adj (A) i,j = (-1) i+j det (A[ j,i ])](https://reader035.vdocument.in/reader035/viewer/2022062521/5681666d550346895dda094b/html5/thumbnails/32.jpg)
Step 2: Determine the leaving variable. Solve for entering column d in current dictionary: d= AB
-1 a where a is the entering column taken from the initial problem.
Or equivalently, solve for d: AB d = a
If tightest constraint corresponds to
xleaving= v – t * xentering
then the new value of xentering will be s= v/t.32
![Page 33: Jacobi’s formula: A -1 = Adj (A) where Adj (A) i,j = (-1) i+j det (A[ j,i ])](https://reader035.vdocument.in/reader035/viewer/2022062521/5681666d550346895dda094b/html5/thumbnails/33.jpg)
Step 3: Update variables (xj enters, xk leaves). Update basic variables headers HB by replacing k with j. Update HN by replacing j with k. Set xj = s in the new solution. Plug this value for xj into the other equations to update the values of the other basic variables. The leaving variable will be 0.
Recall r = the coefficient of xj in the z row:Set z = z + r s.
Update AB, AN, cB, cN, xB, xN to match basis headers.
33