Vectors and Matrices

A vector is

[Graphics:../Images/mma_primer_gr_19.gif]
[Graphics:../Images/mma_primer_gr_20.gif]

Another is

[Graphics:../Images/mma_primer_gr_21.gif]
[Graphics:../Images/mma_primer_gr_22.gif]

Take the "dot" product of these.

[Graphics:../Images/mma_primer_gr_23.gif]
[Graphics:../Images/mma_primer_gr_24.gif]

This looks good for the scalar multiplication of vectors.  We can also do:

[Graphics:../Images/mma_primer_gr_25.gif]
[Graphics:../Images/mma_primer_gr_26.gif]

One kind of vector product is:

[Graphics:../Images/mma_primer_gr_27.gif]
[Graphics:../Images/mma_primer_gr_28.gif]

The normal cross product can also be obtained from the determinant of a matrix as you will see below.

You can also do

[Graphics:../Images/mma_primer_gr_29.gif]
[Graphics:../Images/mma_primer_gr_30.gif]

But this multiplication does not arise in any physical problems that we will be considering this semester.  It could be useful for numerical solution techniques.    

We would like to make a matrix that Mathematica can understand.  Use an opening {{ then close each row with } and then open the next with { .......  after the last row you obviously need a closing }.  The one that is defined here has the name "a"

[Graphics:../Images/mma_primer_gr_31.gif]
[Graphics:../Images/mma_primer_gr_32.gif]
[Graphics:../Images/mma_primer_gr_33.gif]
[Graphics:../Images/mma_primer_gr_34.gif]

The matrix a is a 3 by 3 square matrix.

A matrix could correspond to many different physical entities.  For example a 3 by 3 could be the elements in a second order tensor which corresponds to the stresses in a fluid at a point and some time.  Of course, if this were the case, the matrix would be symmetric for all simple fluids. (Yes, this course will be complicated at times!)

How can we get a symmetric matrix from a and what is "symmetric".

First define transpose, which is a switching of a[i,j] --> a[j,i]

[Graphics:../Images/mma_primer_gr_35.gif]
[Graphics:../Images/mma_primer_gr_36.gif]

You see that we left the diagonal elements unchanged and "reflected" the off diagonal elements across the diagonal.  For example elements a[i,j] , where i is row and j is column were transformed to a[ j,i]

[Graphics:../Images/mma_primer_gr_37.gif]
[Graphics:../Images/mma_primer_gr_38.gif]
[Graphics:../Images/mma_primer_gr_39.gif]
[Graphics:../Images/mma_primer_gr_40.gif]

Now we can easily get the symmetric part of a (such that a[i,j]=a[j,i] and (of course) its antisymmetric piece.

[Graphics:../Images/mma_primer_gr_41.gif]
[Graphics:../Images/mma_primer_gr_42.gif]
[Graphics:../Images/mma_primer_gr_43.gif]
[Graphics:../Images/mma_primer_gr_44.gif]

We see that the antisymmetric part does not have any diagonal elements because a[i,i] cannot be equal to -a[i,i].  Thus our definition of an anti symmetric matrix is that a[i,j]=-a[j,i].  

Symmetric and antisymmetric matrixes are useful, for example, in defining the basic elements of fluid motion in a study of "kinematics".  It is convenient that

[Graphics:../Images/mma_primer_gr_45.gif]
[Graphics:../Images/mma_primer_gr_46.gif]

thus we have recovered a again.

We need to clear the value of a because we use it again below.

[Graphics:../Images/mma_primer_gr_47.gif]
[Graphics:../Images/mma_primer_gr_48.gif]
[Graphics:../Images/mma_primer_gr_49.gif]

The determinant is a standard operation.  You recall

[Graphics:../Images/mma_primer_gr_50.gif]
[Graphics:../Images/mma_primer_gr_51.gif]

As mentioned above, the cross product can be defined with the determinant

[Graphics:../Images/mma_primer_gr_52.gif]
[Graphics:../Images/mma_primer_gr_53.gif]
[Graphics:../Images/mma_primer_gr_54.gif]
[Graphics:../Images/mma_primer_gr_55.gif]

Where we had to work a little too hard with the {ex,ey,ez} being a unit vector and the Coefficient being used to recover the components of each direction.  

[Graphics:../Images/mma_primer_gr_56.gif]
[Graphics:../Images/mma_primer_gr_57.gif]

From these operations we see that addition and subtraction of matrices is done by simple addition and subtraction of corresponding elements.


Converted by Mathematica      August 20, 2000