How to check if a matrix is symmetric in MATLAB

A matrix is a two-dimensional array of numbers arranged in rows and columns. In MATLAB, matrices are fundamental data structures used for various numerical computations.

Let's see an example of how to create a matrix:

A = [1, 2, 3; 4, 5, 6; 7, 8, 9]; % Defines a 3x3 matrix

Symmetric matrix

A symmetric matrix is a square matrix that is equal to its transpose. In other words, for a symmetric matrix $A$ , $A = A'$ .

Let's see an example of how to create a symmetric matrix:

B = [1, 2, 3; 2, 4, 5; 3, 5, 6]; % Defines a 3x3 symmetric matrix

In the illustration below notice that in the symmetric matrix, the values are symmetric across the main diagonal, whereas in the nonsymmetric matrix, this symmetry property is not satisfied.

Symmetric matrix vs. nonsymmetric

Transpose

Transpose of a matrix is obtained by interchanging its rows and columns. In MATLAB, the transpose of a matrix $A$ is obtained using the $'$ operator or the transpose function.

Let's see an example of taking the transpose of a matrix:

A_transpose = A'; % Using the ' operator
A_transpose_alt = transpose(A); % Using the transpose function

Symmetry check in MATLAB

If we need to check whether a matrix is symmetric, we can compare it to its transpose. We can check it in MATLAB using the following ways:

Using `issymmetric()`

The issymmetric() function in MATLAB checks whether a given matrix is symmetric. It returns 1 if the matrix is symmetric, meaning it is equal to its own transpose, and 0 otherwise.

Syntax

The basic syntax of the function is as follows:

issymmetric(A)

Here, A is the matrix for which we want to check if it is symmetric or not.

Using `if-else` statement

The if-else statement works with the isequal() function comparing the matrix with its transpose.
If a matrix, $A$ , is symmetric, then $A$ equals its transpose, i.e., $A = A'$ .

Example

Let’s implement these ways in the code widget below:

A = [1 0; 0 1];
B = issymmetric(A); % Check if the matrix is symmetric
if B
    disp('A is symmetric.');
else
    disp('A is not symmetric.');
end
C = [1 2 3; 4 5 6];
% Check if the matrix is symmetric
if isequal(C, C')
    disp('C is symmetric.');
else
    disp('C is not symmetric.');
end