MAT 232 Chapters 3 and 4

1613 days ago by kcrisman

Chapter 3 - Determinants

Wouldn't it be nice to have a function that told us something about the matrix in one number?

@interact() def _(A=matrix(RDF,[[1,0],[0,1]])): def maketriangle(M): vertex1 = M*vector((0,0)) vertex2 = M*vector((1,3)) vertex3 = M*vector((-1,2)) edges = line([vertex1,vertex2,vertex3,vertex1])+point(vertex1,size=40,color='black')+point(vertex2,size=40,color='red')+point(vertex3,size=40,color='green')+point((0,0),size=0) return edges html('triangle, transformed by $%s$'%latex(A)) G = maketriangle(A) + maketriangle(identity_matrix(2)) G.show(aspect_ratio=1,figsize=[3,3]) html('Notice the new triangle is $%s$ times as big'%latex(det(A))) 
       

Click to the left again to hide and once more to show the dynamic interactive window

We'll return to this.  The function we need is called the determinant.

How might one calculate the determinant?  With Sage, all you need is one thing.

M = matrix([[1,2],[3,4]]) 
       

Just use "det".

det(M) 
       
-2
-2

(Alternately, you can use the dot notation.)

M.determinant() 
       
-2
-2

Use square matrices!

det(matrix([[1,2]])) 
       
Traceback (click to the left of this block for traceback)
...
ValueError: self must be a square matrix
Traceback (most recent call last):
  File "<stdin>", line 1, in <module>
  File "_sage_input_9.py", line 10, in <module>
    exec compile(u'open("___code___.py","w").write("# -*- coding: utf-8 -*-\\n" + _support_.preparse_worksheet_cell(base64.b64decode("ZGV0KG1hdHJpeChbWzEsMl1dKSk="),globals())+"\\n"); execfile(os.path.abspath("___code___.py"))
  File "", line 1, in <module>
    
  File "/tmp/tmp4XvcVm/___code___.py", line 3, in <module>
    exec compile(u'det(matrix([[_sage_const_1 ,_sage_const_2 ]]))
  File "", line 1, in <module>
    
  File "/usr/local/sage-5.6-linux-64bit-ubuntu_8.04.4_lts-x86_64-Linux/local/lib/python2.7/site-packages/sage/misc/functional.py", line 337, in det
    return x.det()
  File "matrix2.pyx", line 1213, in sage.matrix.matrix2.Matrix.det (sage/matrix/matrix2.c:9041)
  File "matrix_integer_dense.pyx", line 3289, in sage.matrix.matrix_integer_dense.Matrix_integer_dense.determinant (sage/matrix/matrix_integer_dense.c:25942)
ValueError: self must be a square matrix
det(matrix([[-1,2,3,0],[3,4,3,0],[5,4,6,6],[4,2,4,3]])) 
       
114
114

We can visualize that $\det(AB)=\det(A)\det(B)$.

@interact() def _(A=matrix(RDF,[[1,0],[0,1]]),B=matrix(RDF,[[1,0],[0,1]])): def maketriangle(M): vertex1 = M*vector((0,0)) vertex2 = M*vector((1,3)) vertex3 = M*vector((-1,2)) edges = line([vertex1,vertex2,vertex3,vertex1])+point(vertex1,size=40,color='black')+point(vertex2,size=40,color='red')+point(vertex3,size=40,color='green')+point((0,0),size=0) return edges html('triangle, transformed by $%s$'%latex(A)) html('and then transformed further by $%s$'%latex(B)) G = maketriangle(A) H = maketriangle(identity_matrix(2)) I = maketriangle(B*A) show(H+G+I,aspect_ratio=1,figsize=[3,3]) html('Notice the second triangle is $%s$ times as big'%latex(det(A))) html('And the third triangle is another $%s$ times as big'%latex(det(B))) html('For a total of $%s\cdot %s=%s$ times as big'%(latex(det(A)),latex(det(B)),latex(det(B*A)))) 
       

Click to the left again to hide and once more to show the dynamic interactive window

How long does it take to compute a determinant?

@interact def _(n=5): if n>100: html("Don't get too big and hog Sage!") else: M = random_matrix(ZZ, n) show(M) time a = det(M) html("determinant is $%s$"%a) 
       

Click to the left again to hide and once more to show the dynamic interactive window

Chapter 4 - Vector Spaces

We didn't use Sage this chapter.  But we could have!

M = matrix([[1,2,3],[1,2,3],[1,2,3]]); show(M) 
       
M.kernel() 
       
Free module of degree 3 and rank 2 over Integer Ring
Echelon basis matrix:
[ 1  0 -1]
[ 0  1 -1]
Free module of degree 3 and rank 2 over Integer Ring
Echelon basis matrix:
[ 1  0 -1]
[ 0  1 -1]
M.column_space() 
       
Free module of degree 3 and rank 1 over Integer Ring
Echelon basis matrix:
[1 1 1]
Free module of degree 3 and rank 1 over Integer Ring
Echelon basis matrix:
[1 1 1]
M.row_space() 
       
Free module of degree 3 and rank 1 over Integer Ring
Echelon basis matrix:
[1 2 3]
Free module of degree 3 and rank 1 over Integer Ring
Echelon basis matrix:
[1 2 3]
M.rank() 
       
1
1
dimension(M.kernel()) + M.rank() 
       
3
3