Engineering Analysis/Linear Transformations

From Wikibooks, open books for an open world
< Engineering Analysis
Jump to: navigation, search

Linear Transformations[edit]

A linear transformation is a matrix M that operates on a vector in space V, and results in a vector in a different space W. We can define a transformation as such:

T:V\to W

In the above equation, we say that V is the domain space of the transformation, and W is the range space of the transformation. Also, we can use a "function notation" for the transformation, and write it as:

M(x) = Mx = y

Where x is a vector in V, and y is a vector in W. To be a linear transformation, the principle of superposition must hold for the transformation:

M(av_1 + bv_2) = aM(v_1) + bM(v_2)

Where a and b are arbitrary scalars.

Null Space[edit]

The Nullspace of an equation is the set of all vectors x for which the following relationship holds:

Mx = 0

Where M is a linear transformation matrix. Depending on the size and rank of M, there may be zero or more vectors in the nullspace. Here are a few rules to remember:

  1. If the matrix M is invertable, then there is no nullspace.
  2. The number of vectors in the nullspace (N) is the difference between the rank(R) of the matrix and the number of columns(C) of the matrix:
N = R - C

If the matrix is in row-eschelon form, the number of vectors in the nullspace is given by the number of rows without a leading 1 on the diagonal. For every column where there is not a leading one on the diagonal, the nullspace vectors can be obtained by placing a negative one in the leading position for that column vector.

We denote the nullspace of a matrix A as:

\mathcal{N}\{A\}

Linear Equations[edit]

If we have a set of linear equations in terms of variables x, scalar coefficients a, and a scalar result b, we can write the system in matrix notation as such:

Ax = b

Where x is a m × 1 vector, b is an n × 1 vector, and A is an n × m matrix. Therefore, this is a system of n equations with m unknown variables. There are 3 possibilities:

  1. If Rank(A) is not equal to Rank([A b]), there is no solution
  2. If Rank(A) = Rank([A b]) = n, there is exactly one solution
  3. If Rank(A) = Rank([A b]) < n, there are infinitely many solutions.

Complete Solution[edit]

The complete solution of a linear equation is given by the sum of the homogeneous solution, and the particular solution. The homogeneous solution is the nullspace of the transformation, and the particular solution is the values for x that satisfy the equation:

A(x) = b
A(x_h + x_p) = b

Where

x_h is the homogeneous solution, and is the nullspace of A that satisfies the equation A(x_h) = 0
x_p is the particular solution that satisfies the equation A(x_p) = b

Minimum Norm Solution[edit]

If Rank(A) = Rank([A b]) < n, then there are infinitely many solutions to the linear equation. In this situation, the solution called the minimum norm solution must be found. This solution represents the "best" solution to the problem. To find the minimum norm solution, we must minimize the norm of x subject to the constraint of:

Ax - b = 0

There are a number of methods to minimize a value according to a given constraint, and we will talk about them later.

Least-Squares Curve Fit[edit]

If Rank(A) doesnt equal Rank([A b]), then the linear equation has no solution. However, we can find the solution which is the closest. This "best fit" solution is known as the Least-Squares curve fit.

We define an error quantity E, such that:

E = Ax - b \ne 0

Our job then is to find the minimum value for the norm of E:

\|E\|^2 = \|Ax - b\|^2 = <Ax -b, Ax-b>

We do this by differentiating with respect to x, and setting the result to zero:

\frac{\partial \|E\|^2}{\partial x} = 2A'(Ax - b) = 0

Solving, we get our result:

x = (A^TA)^{-1}A^Tb