3DICA 1: Vector and Matrix Algebra

eZine lover (@eZine)

Published in 

3DICA

 · 1 year ago

1.1 Vectors

1.1.1 General

Vectors are lines which have a constant direction and length (or magnitude), and can be moved arbitrarily in a given coordinate system. A vector is drawn as in the following picture (a short perpendicular line (not necessary) at the starting point, and an arrowhead at the ending point):

Vectors are usually marked the same way as the coordinates of a point: a=(x_1, y_1, z_1) , where x_1 is the x-axis vector coefficient, y_1 and z_1 respectively. This vector \vec a can be drawn so that we draw a point into the coordinates (x_1,y_1,z_1) , and a line from the origin to that point. The x-, y-, and z-axis unit vectors (vectors that have a length of one) have their own special symbols: i, j, and k, respectively. So the x-axis unit vector is i=(1,0,0) , etc.

Vectors are marked as alphabetical symbols with a line over them, or the symbol being written with a bold font. I'm naturally using here the latter way.

1.1.2 Vector equation

Now we'll calculate the equation of a vector \overrightarrow {AB} . If the point A is (Ax, Ay, Az), and the point B (Bx, By, Bz), where x, y, and z are subscripts, the equation of \overrightarrow {AB} (the vector which starts in the point A and ends in the point B that is) is

Note! The starting point is always subtracted from the ending point!

1.1.3 Vector length

The length of the vector AB is marked as the vector's absolute value and calculated by adding the squares of its terms together and taking a square root of the result:

Calculating a vector's length is needed for example with phong shading and vertex normal calculation.

Problem: Let a = (1,1,1). Derive the unit vector which is parallel to a.

Solution: |a| = \sqrt {\left( 1^2 + 1^2 + 1^2 \right)} = \sqrt {\left(3\right)} , the unit vector which is parallel to a:

Vector addition is performed by adding the terms together. Subtraction respectively, but a negative vector points to the opposite direction than a positive vector.

Solution: 2a = (2*2,2*(-6)) = (4,-12)

2a-b = (4-1,-12-3) = (3,-15)

then

One vector can be multiplied by another by two ways, and the first one is called the dot product (or scalar product if you wish). A dot product is read like "a dot b" and it is calculated by multiplying the length of a by the length of b and by the cosine of the angle between the two vectors. Another way to calculate it is to multiply the terms of the two vectors together and to add the products together. The geometric meaning of dot product is the length of a projected to b.

Note! The angle between the two vectors means always the smaller one between them when the two vectors start from the same point.

a \cdot b = |a| |b| \cos(a,b) = A_x \cdot B_x + A_y \cdot B_y + A_z \cdot B_z
Cos(a,b) means the angle alpha.

Note! Vectors are perpendicular when the angle between them is 90 degrees (cos90 = 0 -> a dot b = 0)!

Dot product is needed for example in phong illuminating and many other light manipulation operations.

Problem: Calculate the angle between the vectors a+b and a-b, where a = (3,1) and b = (1,-7).

Solution:

1.1.6 Vector projection

We are trying to calculate the projection of a vector on another vector, or to calculate vector a's component which is parallel to a vector b.

Note! All vectors should be ensured to be unit vectors for the calculations to give the correct answer!

To clarify, all the vectors in the picture start from the same point, and the projection vector of a on b is the short vector which is parallel to b. All we need to do is to calculate the dot product of a and b and to multiply b by it:

a's component parallel to b = (a dot b)*b.

The result of a dot product is a constant, so (a dot b) just scales b which at least I find rational :) Vector projection is needed at least when converting a vector to a matrix.

1.1.7 Cross product

In the 3-space, we define the cross product (or vector product) axb ("a cross b") of two vectors a and b as a vector which is perpendicular to each of the vectors a and b and whose length |a x b| tells the area of the parallelogram defined by the two vectors:

|a \times b| = |a| \, |b| \sin(a , b)
(Note: this is the equation for the area of a parallelogram, and is used also in geometry)

Note! a and b need not to be perpendicular to each other even though you might think so by looking at the picture X)

Some other things to remember about cross product:

Note! Every second sign is a +, every second a -!

REMEMBER this! (the signs go like +-+-+...)

There's no restrictions in the size of a determinant, and all of them can be solved respectively. For example, from a four-row determinant we generate four three-row determinants etc.

Cross product is needed for example when calculating normals and plane equations.

Problem: Find the area of the parallelogram determined by the vectors (1,-1,2) and (2,3,-1).

Solution: (We use the hard way to demonstrate determinants)

A scalar triple product is marked as (axb) dot c, and its absolute value tells the volume of the parallelogramic object determined by the three vectors (cf. parallelogram in cross product).

Note! If this volume is zero, the three vectors are of course on the same plane!

Problem: Find a constant a so that the vectors (1,-1,2), (1,3,-1), and (a,-4,1) are on the same plane. Solution: On the same plane -> scalar triple product = 0. We first calculate the cross product of arbitrary two of the vectors and then the dot product of this vector and the third vector. This ought to be zero, -5a - 15 = 0 <=> -5a = 15 <=> a = -3.

1.2 Matrices

1.2.1 General

(Algebra ja Geometria, Simo K. Kivela, Otatieto 1989; translated into English)

A specially ordered table of real or complex numbers is caller a matrix:

The numbers a(y,x) are matrix elements. They form p rows and n columns. The type of the matrix is p x n; if needed, it can also be written together with the matrix symbol: A(pxn). The first index j of a matrix element a(j,k) shows the row on which the element is located; the latter index k shows respectively the column (the indices go thus like in C's tables: compared to the representation of a point (x,y), in matrices and C language tables the indices are (y,x)).

We use capital latin letters as names or symbols for matrices.

Matrices of the type 1 x n are called horizontal- or row vectors. Exceptionally, they are marked with small latin letters:

A matrix whose all elements are zero, is called a zero matrix. The type of a zero matrix is not restricted. The symbol is O (or o if the matrix is a vector).

A matrix is called orthogonal if its row or column vectors (that is, the vectors formed by the rows or columns of the matrix) are perpendicular to each other and they are unit vectors.

A matrix B defined like this

1.2.2 Matrix operations

1.2.2.1 Addition

Matrix addition is defined so that every element b(j,k) of the second matrix is added to every element a(j,k) of the first matrix:

(A scalar = a constant) This goes like addition, but instead of the +'s we put *'s, and instead of b(j,k) we put the scalar with which we want to multiply the matrix. Multiplying by a scalar scales a matrix.

1.2.2.3 Multiplication

Here's the algorithm:

Note! Only certain types of matrices can be multiplied by each other. That is, if the first matrix is of the size p x q, the second one should be of the size q x r, where p, q, and r are not restricted. When we multiply a p x q sized matrix by a q x r sized matrix, we get a p x r sized matrix as a result, as can easily be seen from the definition of matrix multiplication.

Matrix multiplication is needed in object transformations.

1.2.2.3.1 Multiplying by a vector

The case is easiest to understand by making the desired vector a horizontal vector:

1.2.2.4 Transposition

Definition: The transpose of a matrix A(pxn) is

NOTE! The transpose of an orthogonal matrix is also its inverse! (this piece of information is used in optimizing 3D routines -- see chapter 2.7.)