Matrix Inversion and Applications with NumPy

Welcome back! As we continue our journey through linear algebra with NumPy, we come to an essential concept: matrix inversion. Matrix inversion plays a crucial role in solving systems of linear equations and in various applications across science and engineering.

Before diving into code, it is important to remember what an invertible matrix is. An invertible matrix, also known as a non-singular matrix, is a square matrix that has an inverse. Only square matrices (matrices with the same number of rows and columns) that have non-zero determinants are invertible. This property is crucial for applications such as solving systems of linear equations, which we will explore soon.

Let's walk through the process of inverting a matrix using NumPy. We'll demonstrate this with a simple 2x2 matrix.

1. First, we define the matrix using np.array:

   Python

   import numpy as np
   
   # Defining an invertible matrix
   matrix = np.array([[4, 7], [2, 6]])

2. We now use NumPy's np.linalg.inv() function to calculate the inverse:

   Python

   # Matrix inversion
   inverse = np.linalg.inv(matrix)

3. To verify the inversion, we multiply the original matrix by its inverse. The result should be the identity matrix, showing that the inversion was successful:

   Python

   # Verifying the inversion (A * A_inv should yield the identity matrix)
   identity_check = np.matmul(matrix, inverse)

4. Finally, we can display the results:

   Python

   # Display results
   print("Matrix:\n", matrix)
   print("Inverse:\n", inverse)
   print("Identity Check (A * A_inv):\n", identity_check)
   
   # Output:
   # Matrix:
   #  [[4 7]
   #  [2 6]]
   # Inverse:
   #  [[ 0.6 -0.7]
   #  [-0.2  0.4]]
   # Identity Check (A * A_inv):
   #  [[1.00000000e+00 -1.11022302e-16]
   #  [1.11022302e-16 1.00000000e+00]]

Each step above confirms the validity of the inverse. Note that the identity matrix, while expected to be exact, may show very small numerical errors due to floating-point precision, as shown above.

Matrix inversion has significant practical applications. One of the classic uses is solving a system of linear equations. In this example, a, b, c, d, e, and f are known values, while x and y are the unknowns that we want to solve for. Consider the following system of equations:

⎡ a  b ⎤   ⎡ x ⎤   ⎡ e ⎤
⎢      ⎥ * ⎢   ⎥ = ⎢   ⎥
⎣ c  d ⎦   ⎣ y ⎦   ⎣ f ⎦

The solution [x, y] can be found using the inverse of the coefficient matrix:

# Example system of equations
coeff_matrix = np.array([[4, 7], [2, 6]])  # [a, b], [c, d] = [4, 7], [2, 6]
constants = np.array([1, 1])               # [e, f] = [1, 1]

# Solving for [x, y]
solution = np.linalg.inv(coeff_matrix).dot(constants)
print("Solution [x, y]:", solution)

# Output:
# Solution [x, y]: [-3.  2.]

Here, the matrix inversion helps us find the values of (x) and (y) that satisfy both equations.