How to calculate the rank of a matrix using Python?

by | Oct 3, 2023 | Featured, Linear Algebra

What is the rank of a matrix?

The rank of a matrix A is the maximal number of linearly independent columns of A. This is, in turn, equal to the maximal number of linearly independent rows of a matrix. For example, let’s say A is a matrix that has 3 rows and 3 columns. And the third column of the matrix is a linear combination of the first two rows. So, the maximal number of independent columns of the matrix is 2. So, the rank of the matrix will be 2.

For example, let’s consider the following matrix:

B=\begin{bmatrix}  1 & 0 & 1 \\  2 & 3 & -1 \\  3 & 3 & 0  \end{bmatrix}

Now, if we subtract the second column from the first column of the matrix, we will get the third column. So, the third column is a linear combination of the first two columns of the matrix. Hence, the rank of the matrix is 2.

Now, a matrix A is said to have full rank if the rank of the matrix is the maximum for the given dimension of the matrix. And that will be the minimum of the number of rows and columns of the matrix. For example, if the number of rows is more than the number of columns, then the full rank will be equal to the number of columns of the matrix.

If a matrix does not have full rank, then the matrix is said to be a rank-deficient matrix. And the rank deficiency of a matrix is equal to the difference between the lesser of the number of rows or columns of the matrix and the rank of the matrix.

How to calculate the rank of a matrix using Python?

We can use the following Python code to calculate the rank of a matrix.

import numpy
from numpy.linalg import matrix_rank

A = numpy.random.randint(low=1, high=10, size=(3,3))
print("Matrix A: \n", A)
r = matrix_rank(A)
print("Rank of the matrix A: ", r)

Here, we are first creating a matrix A that contains random integers. We are using the numpy.random.randint() function for that purpose. The generated random integers will be within the range [low, high).

Now, we are using the matrix_rank() function from the numpy.linalg module to calculate the rank of the matrix using NumPy.

The output of the mentioned program will be:

Matrix A: 
 [[8 8 4]
 [2 7 9]
 [4 5 5]]
Rank of the matrix A:  3
Facebooktwitterredditpinterestlinkedinmail

Calculate the pseudoinverse of a matrix using Python

What is the pseudoinverse of a matrix? We know that if A is a square matrix with full rank, then A-1 is said to be the inverse of A if the following condition holds: $latex AA^{-1}=A^{-1}A=I $ The pseudoinverse or the Moore-Penrose inverse of a matrix is a...

Cholesky decomposition using Python

What is Cholesky decomposition? A square matrix A is said to have Cholesky decomposition if it can be written as a product of a lower triangular matrix and its conjugate transpose. $latex A=LL^{*} $ If all the entries of A are real numbers, then the conjugate...

Tensor Hadamard Product using Python

In one of our previous articles, we already discussed what the Hadamard product in linear algebra is. We discussed that if A and B are two matrices of size mxn, then the Hadamard product of A and B is another mxn matrix C such that: $latex H_{i,j}=A_{i,j} \times...

Perform tensor addition and subtraction using Python

We can use numpy nd-array to create a tensor in Python. We can use the following Python code to perform tensor addition and subtraction. import numpy A = numpy.random.randint(low=1, high=10, size=(3, 3, 3)) B = numpy.random.randint(low=1, high=10, size=(3, 3, 3)) C =...

How to create a tensor using Python?

What is a tensor? A tensor is a generalization of vectors and matrices. It is easily understood as a multidimensional array. For example, in machine learning, we can organize data in an m-way array and refer it as a data tensor. Data related to images, sounds, movies,...

How to combine NumPy arrays using horizontal stack?

We can use the hstack() function from the numpy module to combine two or more NumPy arrays horizontally. For example, we can use the following Python code to combine three NumPy arrays horizontally. import numpy A = numpy.array([[1, 2, 3], [4, 5, 6], [7, 8, 9]]) B =...

How to combine NumPy arrays using vertical stack?

Let’s say we have two or more NumPy arrays. We can combine these NumPy arrays vertically using the vstack() function from the numpy module. For example, we can use the following Python code to combine three NumPy arrays vertically. import numpy A = numpy.array([[1, 2,...

Singular Value Decomposition (SVD) using Python

What is Singular Value Decomposition (SVD)? Let A be an mxn rectangular matrix. Using Singular Value Decomposition (SVD), we can decompose the matrix A in the following way: $latex A_{m \times n}=U_{m \times m}S_{m \times n}V_{n \times n}^T $ Here, U is an mxm matrix....

Eigen decomposition of a square matrix using Python

Let A be a square matrix. Let’s say A has k eigenvalues λ1, λ2, ... λk. And the corresponding eigenvectors are X1, X2, ... Xk. $latex X_1=\begin{bmatrix} x_{11} \\ x_{21} \\ x_{31} \\ ... \\ x_{k1} \end{bmatrix} \\ X_2=\begin{bmatrix} x_{12} \\ x_{22} \\ x_{32} \\ ......

How to calculate eigenvalues and eigenvectors using Python?

In our previous article, we discussed what eigen values and eigenvectors of a square matrix are and how we can calculate the eigenvalues and eigenvectors of a square matrix mathematically. We discussed that if A is a square matrix, then $latex (A- \lambda I) \vec{u}=0...

Amrita Mitra

Author

Ms. Amrita Mitra is an author, who has authored the books “Cryptography And Public Key Infrastructure“, “Web Application Vulnerabilities And Prevention“, “A Guide To Cyber Security” and “Phishing: Detection, Analysis And Prevention“. She is also the founder of Asigosec Technologies, the company that owns The Security Buddy.

0 Comments

Submit a Comment

This site uses Akismet to reduce spam. Learn how your comment data is processed.

Not a premium member yet?

Please follow the link below to buy The Security Buddy Premium Membership.

Featured Posts

Translate »