![Page 1: ME451 Kinematics and Dynamics of Machine Systems Review of Matrix Algebra – 2.2 September 13, 2011 Dan Negrut University of Wisconsin-Madison © Dan Negrut,](https://reader036.vdocument.in/reader036/viewer/2022062300/56649dc55503460f94ab8b98/html5/thumbnails/1.jpg)
ME451 Kinematics and Dynamics
of Machine Systems
Review of Matrix Algebra – 2.2 September 13, 2011
Dan NegrutUniversity of Wisconsin-Madison
© Dan Negrut, 2011ME451, UW-Madison
![Page 2: ME451 Kinematics and Dynamics of Machine Systems Review of Matrix Algebra – 2.2 September 13, 2011 Dan Negrut University of Wisconsin-Madison © Dan Negrut,](https://reader036.vdocument.in/reader036/viewer/2022062300/56649dc55503460f94ab8b98/html5/thumbnails/2.jpg)
Before we get started…
Due next week: Problems: 2.2.5, 2.2.8. 2.2.10 out of Haug’s book (http://sbel.wisc.edu/Courses/ME451/2010/bookHaugPointers.htm)
Due on Th: In class, if pen and paper version is submitted At 23:59 PM if electronic form submitted
Moving to full electronic submission starting after September See Forum posting for some ideas on how to go about typing equations in your document in Windows
Last time: Covered Geometric Vectors & operations with them Justified the need for Reference Frames (using a vector basis) Introduced algebraic representation of a vector & related operations Rotation Matrix (for switching from one RF to another RF)
Today: Dealing with the kinematics of a body: rotation + translation Quick review of matrix/vector algebra Discuss concept of “generalized coordinates”
2
![Page 3: ME451 Kinematics and Dynamics of Machine Systems Review of Matrix Algebra – 2.2 September 13, 2011 Dan Negrut University of Wisconsin-Madison © Dan Negrut,](https://reader036.vdocument.in/reader036/viewer/2022062300/56649dc55503460f94ab8b98/html5/thumbnails/3.jpg)
The Rotation Matrix A
Very important observation ! the matrix A is orthonormal:
3
![Page 4: ME451 Kinematics and Dynamics of Machine Systems Review of Matrix Algebra – 2.2 September 13, 2011 Dan Negrut University of Wisconsin-Madison © Dan Negrut,](https://reader036.vdocument.in/reader036/viewer/2022062300/56649dc55503460f94ab8b98/html5/thumbnails/4.jpg)
Important Relation
Expressing a given vector in one reference frame (local) in a different reference frame (global)
4Also called a change of base.
![Page 5: ME451 Kinematics and Dynamics of Machine Systems Review of Matrix Algebra – 2.2 September 13, 2011 Dan Negrut University of Wisconsin-Madison © Dan Negrut,](https://reader036.vdocument.in/reader036/viewer/2022062300/56649dc55503460f94ab8b98/html5/thumbnails/5.jpg)
Example 1[Deals with the rotation of a body wrt a Global Reference Frame (GRF)]
5
Express the geometric vector in the local reference frame O’x’y’.
Express the same geometric vector in the global reference frame OXY
Do the same for the geometric vector
y’ x’ θ
O’
E
B
L
X
Y
O
![Page 6: ME451 Kinematics and Dynamics of Machine Systems Review of Matrix Algebra – 2.2 September 13, 2011 Dan Negrut University of Wisconsin-Madison © Dan Negrut,](https://reader036.vdocument.in/reader036/viewer/2022062300/56649dc55503460f94ab8b98/html5/thumbnails/6.jpg)
6
Express the geometric vector
in the local reference frame O’x’y’. Express the same geometric
vector in the global reference frame OXY
Do the same for the geometric vector
L
y’
x’
θ
P
G
Y
XO
O’
Example 2[Deals with the rotation of a body wrt a Global Reference Frame (GRF)]
![Page 7: ME451 Kinematics and Dynamics of Machine Systems Review of Matrix Algebra – 2.2 September 13, 2011 Dan Negrut University of Wisconsin-Madison © Dan Negrut,](https://reader036.vdocument.in/reader036/viewer/2022062300/56649dc55503460f94ab8b98/html5/thumbnails/7.jpg)
The Kinematics of a Rigid Body:
Handling both Translation + Rotation
What we just discussed dealt with the case when you are interested in finding the representing the location of a point P when you change the reference frame, but yet the new and old reference frames share the same origin
7
What if they don’t share the same origin (see picture at right)? How would you represent the position of the point P in this new reference frame?
![Page 8: ME451 Kinematics and Dynamics of Machine Systems Review of Matrix Algebra – 2.2 September 13, 2011 Dan Negrut University of Wisconsin-Madison © Dan Negrut,](https://reader036.vdocument.in/reader036/viewer/2022062300/56649dc55503460f94ab8b98/html5/thumbnails/8.jpg)
More on Body Kinematics
A lot of ME451 is based on the ability to look at the location of one point P in two different reference frames: a local reference frame (LRF) and a global reference frame (GRF)
Local reference frame is typically fixed (rigidly attached) to a body that is moving in space
Global reference frame is the “world” reference frame: it’s not moving, and serve as the universal reference frame
² In the LRF, the position of point P is described by s0P (sometimes the notationused is ¹sP )
² In theGRF, the position of point P is described by rP (see next slide)
8
![Page 9: ME451 Kinematics and Dynamics of Machine Systems Review of Matrix Algebra – 2.2 September 13, 2011 Dan Negrut University of Wisconsin-Madison © Dan Negrut,](https://reader036.vdocument.in/reader036/viewer/2022062300/56649dc55503460f94ab8b98/html5/thumbnails/9.jpg)
ME451 Important Slide
O’
X
Y
O
'Psx’
y’
frP
r
P
9
![Page 10: ME451 Kinematics and Dynamics of Machine Systems Review of Matrix Algebra – 2.2 September 13, 2011 Dan Negrut University of Wisconsin-Madison © Dan Negrut,](https://reader036.vdocument.in/reader036/viewer/2022062300/56649dc55503460f94ab8b98/html5/thumbnails/10.jpg)
Example
The location of point O’ in the OXY global RF is [x,y]T. The orientation of the bar is described by the angle 1. Find the location of C and D expressed in the global reference frame as functions of x, y, and 1.
1
D
y1′
x1′φ1
2
2Y
X
C
O
O′
10
![Page 11: ME451 Kinematics and Dynamics of Machine Systems Review of Matrix Algebra – 2.2 September 13, 2011 Dan Negrut University of Wisconsin-Madison © Dan Negrut,](https://reader036.vdocument.in/reader036/viewer/2022062300/56649dc55503460f94ab8b98/html5/thumbnails/11.jpg)
11
Matrix Review
![Page 12: ME451 Kinematics and Dynamics of Machine Systems Review of Matrix Algebra – 2.2 September 13, 2011 Dan Negrut University of Wisconsin-Madison © Dan Negrut,](https://reader036.vdocument.in/reader036/viewer/2022062300/56649dc55503460f94ab8b98/html5/thumbnails/12.jpg)
Recall Notation Conventions
A bold upper case letter denotes matrices Example: A, B, etc.
A bold lower case letter denotes a vector Example: v, s, etc.
A letter in italics format denotes a scalar quantity Example:
12
1,a b
![Page 13: ME451 Kinematics and Dynamics of Machine Systems Review of Matrix Algebra – 2.2 September 13, 2011 Dan Negrut University of Wisconsin-Madison © Dan Negrut,](https://reader036.vdocument.in/reader036/viewer/2022062300/56649dc55503460f94ab8b98/html5/thumbnails/13.jpg)
Matrix Review
What is a matrix? A tableau of numbers ordered by row and column.
Matrix addition:
Addition is commutative
13
11 12 1
21 22 2 22
1
1
1 2
T
T
T
n
nn
m m m mn
a a a
a a a
a a a
é ùé ùê úê úê úê úê úê ú é ù= = = ê úê
¼
¼¼
¼ ¼ ¼ ¼
¼
ú ê úë û ê úê úê úê úê úê úë û ë û
A a a a
a
a
a
L
+ = +A B B A
[ ], 1 , 1
[ ], 1 , 1
[ ],
ij
ij
ij ij ij ij
a i m j n
b i m j n
c c a b
A
B
C A B
![Page 14: ME451 Kinematics and Dynamics of Machine Systems Review of Matrix Algebra – 2.2 September 13, 2011 Dan Negrut University of Wisconsin-Madison © Dan Negrut,](https://reader036.vdocument.in/reader036/viewer/2022062300/56649dc55503460f94ab8b98/html5/thumbnails/14.jpg)
Matrix Multiplication
Recall dimension constraints on matrices so that they can be multiplied: # of columns of first matrix is equal to # of rows of second matrix
Matrix multiplication is not commutative
Distributivity of matrix multiplication with respect to matrix addition:
14
*
*
*
1
[ ],
[ ],
[ ]· ,
m nij
n pij
m pij
n
ij ik kjk
a
c
d
d a c=
= Î
= Î
= = Î
= å
A A
C C
D AC D
¡
¡
¡
![Page 15: ME451 Kinematics and Dynamics of Machine Systems Review of Matrix Algebra – 2.2 September 13, 2011 Dan Negrut University of Wisconsin-Madison © Dan Negrut,](https://reader036.vdocument.in/reader036/viewer/2022062300/56649dc55503460f94ab8b98/html5/thumbnails/15.jpg)
Matrix-Vector Multiplication
A column-wise perspective on matrix-vector multiplication (part of your HW)
Example:
A row-wise perspective on matrix-vector multiplication:
15
11 12 1 1 1
21 22 2 2 21 2
1
1 2
n
nn
n i ii
m m mn n n
a a a v v
a a a v vv
a a a v v=
é ùé ù é ùê úê ú ê úê úê ú ê úê úê ú ê úé ù= = =ê úê ú ê úê úë ûê úê ú ê úê úê ú ê úê úê ú ê úë ûë û ë û
¼
¼¼
¼ ¼ ¼ ¼
¼
åAv a a a aL L
1 1
2 2
T T
T T
T Tm m
é ù é ùê ú ê úê ú ê úê ú ê ú= =ê ú ê úê ú ê úê ú ê úê ú ê úë û ë û
v
A vv
v
va a
a a
a a
L L
1 4 2 0 1 1 4 2 0 7
2 3 1 1 2 2 3 1 1 8(1) (2) ( 1) (1)
1 0 1 1 1 1 0 1 1 3
0 1 1 2 1 0 1 1
· · · ·
2 1
Av
![Page 16: ME451 Kinematics and Dynamics of Machine Systems Review of Matrix Algebra – 2.2 September 13, 2011 Dan Negrut University of Wisconsin-Madison © Dan Negrut,](https://reader036.vdocument.in/reader036/viewer/2022062300/56649dc55503460f94ab8b98/html5/thumbnails/16.jpg)
Orthogonal & Orthonormal Matrices
Definition (Q, orthogonal matrix): a square matrix Q is orthogonal if the product QTQ is a diagonal matrix
Matrix Q is called orthonormal if it’s orthogonal and also QTQ=In
Note that in many cases people fail to make a distinction between an orthogonal and orthonormal matrix. We’ll observe this distinction
Note that if Q is an orthonormal matrix, then Q-1=QT
Example, orthonormal matrix:
16
![Page 17: ME451 Kinematics and Dynamics of Machine Systems Review of Matrix Algebra – 2.2 September 13, 2011 Dan Negrut University of Wisconsin-Madison © Dan Negrut,](https://reader036.vdocument.in/reader036/viewer/2022062300/56649dc55503460f94ab8b98/html5/thumbnails/17.jpg)
Exercise
Prove that the orientation A matrix is orthonormal
A =
2
4cosÁ ¡ sinÁ
sinÁ cosÁ
3
5
17
![Page 18: ME451 Kinematics and Dynamics of Machine Systems Review of Matrix Algebra – 2.2 September 13, 2011 Dan Negrut University of Wisconsin-Madison © Dan Negrut,](https://reader036.vdocument.in/reader036/viewer/2022062300/56649dc55503460f94ab8b98/html5/thumbnails/18.jpg)
Remark:
On the Columns of an Orthonormal Matrix
Assume Q is an orthonormal matrix
18
In other words, the columns of an orthonormal matrix have unit norm and are mutually perpendicular to each other
![Page 19: ME451 Kinematics and Dynamics of Machine Systems Review of Matrix Algebra – 2.2 September 13, 2011 Dan Negrut University of Wisconsin-Madison © Dan Negrut,](https://reader036.vdocument.in/reader036/viewer/2022062300/56649dc55503460f94ab8b98/html5/thumbnails/19.jpg)
Matrix Review [Cntd.]
Scaling of a matrix by a real number: scale each entry of the matrix
Example:
Transpose of a matrix A dimension m£n: a matrix B=AT of dimension n£m whose (i,j) entry is the (j,i) entry of original matrix A
19
· ·[ ] [ · ]ij ija a A
1 4 2 0 1.5 6 3 0
2 3 1 1 3 4.5 1.5 1.5(1.5)
1 0 1 1 1.5 0 1.5 1.5
0 1 1 2 0 1.
·
5 1.5 3
1 4 2 0 1 2 1 0
2 3 1 1 4 3 0 1
1 0 1 1 2 1 1 1
0 1 1 2 0 1 1 2
T
![Page 20: ME451 Kinematics and Dynamics of Machine Systems Review of Matrix Algebra – 2.2 September 13, 2011 Dan Negrut University of Wisconsin-Madison © Dan Negrut,](https://reader036.vdocument.in/reader036/viewer/2022062300/56649dc55503460f94ab8b98/html5/thumbnails/20.jpg)
Linear Independence of Vectors
Definition: linear independence of a set of m vectors, v1,…, vm :
The vectors are linearly independent if the following condition holds
If a set of vectors are not linearly independent, they are called dependent Example:
Note that v1-2v2-v3=0 20
1 1 1.... 0m m n m v v 0
1 2 3
1 0 1
1 1 3
2 4 6
v v v
v1; ::::;vm 2 Rn
![Page 21: ME451 Kinematics and Dynamics of Machine Systems Review of Matrix Algebra – 2.2 September 13, 2011 Dan Negrut University of Wisconsin-Madison © Dan Negrut,](https://reader036.vdocument.in/reader036/viewer/2022062300/56649dc55503460f94ab8b98/html5/thumbnails/21.jpg)
Matrix Rank
Row rank of a matrix Largest number of rows of the matrix that are linearly independent A matrix is said to have full row rank if the rank of the matrix is equal to
the number of rows of that matrix
Column rank of a matrix Largest number of columns of the matrix that are linearly independent
NOTE: for each matrix, the row rank and column rank are the same This number is simply called the rank of the matrix It follows that
21
![Page 22: ME451 Kinematics and Dynamics of Machine Systems Review of Matrix Algebra – 2.2 September 13, 2011 Dan Negrut University of Wisconsin-Madison © Dan Negrut,](https://reader036.vdocument.in/reader036/viewer/2022062300/56649dc55503460f94ab8b98/html5/thumbnails/22.jpg)
Matrix Rank, Example
What is the row rank of the matrix J?
What is the rank of J?
22
2 1 1 0
4 2 2 1
0 4 0 1
J
![Page 23: ME451 Kinematics and Dynamics of Machine Systems Review of Matrix Algebra – 2.2 September 13, 2011 Dan Negrut University of Wisconsin-Madison © Dan Negrut,](https://reader036.vdocument.in/reader036/viewer/2022062300/56649dc55503460f94ab8b98/html5/thumbnails/23.jpg)
Matrix Review [Cntd.]
Symmetric matrix: a square matrix A for which A=AT
Skew-symmetric matrix: a square matrix B for which B=-BT
Examples:
Singular matrix: square matrix whose determinant is zero
Inverse of a square matrix A: a matrix of the same dimension, called A-1, that satisfies the following:
23
2 1 1 0 1 2
1 0 3 1 0 4
1 3 4 2 4 0
A B
![Page 24: ME451 Kinematics and Dynamics of Machine Systems Review of Matrix Algebra – 2.2 September 13, 2011 Dan Negrut University of Wisconsin-Madison © Dan Negrut,](https://reader036.vdocument.in/reader036/viewer/2022062300/56649dc55503460f94ab8b98/html5/thumbnails/24.jpg)
Singular vs. Nonsingular Matrices
Let A be a square matrix of dimension n. The following are equivalent:
24
![Page 25: ME451 Kinematics and Dynamics of Machine Systems Review of Matrix Algebra – 2.2 September 13, 2011 Dan Negrut University of Wisconsin-Madison © Dan Negrut,](https://reader036.vdocument.in/reader036/viewer/2022062300/56649dc55503460f94ab8b98/html5/thumbnails/25.jpg)
Other Useful Formulas[Pretty straightforward to prove true]
If A and B are invertible, their product is invertible and
Also,
For any two matrices A and B that can be multiplied
For any three matrices A, B, and C that can be multiplied
25
![Page 26: ME451 Kinematics and Dynamics of Machine Systems Review of Matrix Algebra – 2.2 September 13, 2011 Dan Negrut University of Wisconsin-Madison © Dan Negrut,](https://reader036.vdocument.in/reader036/viewer/2022062300/56649dc55503460f94ab8b98/html5/thumbnails/26.jpg)
26
Absolute (Cartesian) Generalized Coordinates
vs.
Relative Generalized Coordinates