Linear equations and matrices

 We see that simultaneous linear equations can be concisely expressed by using a matrix. Furthermore, we can find a general way to solve such equations based on matrix algebra. Here, we mostly deal with 2×2 matrices, but the same principle applies to matrices of any size.



Definition (Kronecker's delta)

Kronecker's delta is the symbol δij such that

δij={1if i=j,0otherwise.

Definition (Identity matrix)

The n×n identity matrix is In=(δij).

The identity matrices have their diagonal elements that are all 1 and off-diagonal elements that are all 0.

Example

I1=(1),I2=(1001),I3=(100010001),

and so on. □

Theorem

Let X=(xij) be an n×m matrix. Then

InX=X=XIm.

Proof. Using the definition of δij, we have

(InX)ij=k=1nδikxkj=xij=(X)ij

and

(XIm)ij=k=1mxikδkj=xij=(X)ij.

Definition (Inverse matrix)

Let XMn. The matrix YMn is said to be an inverse matrix of X if XY=In. A matrix that has an inverse is said to be invertible.

It will be proved that an inverse matrix Y of X is unique. That is, if there were two inverse matrices of X, say Y and Z, then Y=Z. Based on this fact, we write X1 to denote the inverse matrix of X. Furthermore, it turns out that X1 is the unique matrix such that X1X=In.

Example. The inverse of (c)M1 is (1/c) if c0. (0) has no inverse.  □

Example. Consider M2. The ``zero'' matrix (0000) has no inverse. (Why?) □

Consider the 2×2 matrix

X=(abcd).

We want to find its inverse. Let Δ=adbc. If Δ=0, X is not invertible (i.e., X1 does not exist). If Δ0, then X is invertible and

X1=(d/Δb/Δc/Δa/Δ).

Exercise. Verify this claim by explicitly computing XX1. □

Definition (Determinant (2×2))

  The determinant of the 2×2 matrix X=(abcd), denoted as det(X) or |X|, is defined by

    det(X)=adbc.

We will study more general determinants in detail in later posts.

Consider the following simultaneous linear equations:

ax+by=p,cx+dy=q.

Using a matrix and vectors defined as

A=(abcd),x=(xy),p=(pq),

we can summarize the above equations as

Ax=p.

If |A|0, multiplying A1 on the both sides from the left, we have

A1(Ax)=A1p. But A1(Ax)=(A1A)x=I2x=x so that x=A1p,

and we have just solved the linear equations.

Example. Consider solving

2x+y=1,3xy=0.

Let

A=(2131), so det(A)=5 and A1=(1/51/53/52/5) so that (xy)=A1(10)=(1/51/53/52/5)(10)=(1535).
You should verify this is indeed the solution. □


Comments

Popular posts from this blog

Birth process

Branching processes: Mean and variance

Informal introduction to formal logic