What is the general algorithm for computing the null space of a 3x3 matrix?

Click For Summary
SUMMARY

The general algorithm for computing the null space of a 3x3 matrix involves augmenting the matrix with a zero vector, performing Gauss-Jordan elimination, and then extracting the null space from the resulting matrix. A specific example illustrates that the method requires adjustment in step three to accurately derive the null space. Additionally, a closed-form solution for computing the eigenvector corresponding to an eigenvalue of a 3x3 matrix is presented, although it is noted that this method does not apply to non-symmetric matrices. The discussion emphasizes the need for a generalized approach to effectively compute the null space across various matrix types.

PREREQUISITES
  • Understanding of Gauss-Jordan elimination
  • Familiarity with null space concepts in linear algebra
  • Knowledge of eigenvalues and eigenvectors
  • Basic proficiency in matrix operations
NEXT STEPS
  • Research advanced techniques for computing the null space of non-symmetric matrices
  • Explore the Strassen algorithm for matrix multiplication
  • Learn about Singular Value Decomposition (SVD) and its applications
  • Study the implications of matrix exponentiation on algorithm efficiency
USEFUL FOR

Students and professionals in mathematics, computer science, and engineering who are focused on linear algebra, matrix computations, and algorithm optimization.

junglebeast
Messages
514
Reaction score
2
There are 2 issues I want to talk about in this post.
(1) General algorithm for gauss-jordan elimination computation of null space
(2) Closed form solution to 3x3 null space

Following the example here,

https://en.wikipedia.org/wiki/Kernel_(linear_algebra)
I thought a general algorithm to compute the null space would be to

1) augment with 0 vector on the right
2) compute gauss-jordan elimination
3) take 2nd to last column, and fill in extra elements with 1's to get the null space

This works in the example provided there. However, on the next example, step 3 needs to be changed...

1, 0, 1
2, 1, 3
1, 1, 2

which has a null space of 1, 1, -1

using gauss-jordan elimination, the closest I can get is

1, 0.5, 1.5, 0
0, 1, 1, 0
0, 0, 0, 0

x1 = -0.5 x2 - 1.5 x3
x2 = -x3
[1, 1, -1]

This gives me the right null space, but step #3 of my above method clearly wasn't right..how can I generalize step 3 into a straight-forward algorithm?

Now onto my second issue.

I found a method to compute the eigenvector corresponding to an eigenvalue of a 3x3 matrix closed form. It is simply:

a1*a5 - a2*(a4-e)
a1*a2 - a5*(a0-e)
(a0-e)*(a4-e) - a1*a1

Removing the 'e's, this is essentially a short cut to get the null space. However it doesn't seem to work for non-symmetric matrices. I feel like there should be a similar method that works for non-symmetric 3x3's...which could be used to avoid the SVD method in this case
 
Last edited by a moderator:
Physics news on Phys.org
Algorithms to solve linear equation systems are well known and are part of basic computer science courses. The fact that there might be examples with faster algorithms is well known, too, and we cannot reason upon examples. This doesn't allow a generalization which would be needed to talk about an algorithm.

If you are interested in the subject, then you might want to read about the improvements on the matrix exponent: https://en.wikipedia.org/wiki/Strassen_algorithm#Asymptotic_complexity
 

Similar threads

  • · Replies 1 ·
Replies
1
Views
3K
  • · Replies 4 ·
Replies
4
Views
16K
  • · Replies 5 ·
Replies
5
Views
3K
  • · Replies 2 ·
Replies
2
Views
2K
  • · Replies 1 ·
Replies
1
Views
4K
  • · Replies 5 ·
Replies
5
Views
2K
  • · Replies 10 ·
Replies
10
Views
2K
  • · Replies 1 ·
Replies
1
Views
2K
  • · Replies 5 ·
Replies
5
Views
12K
  • · Replies 2 ·
Replies
2
Views
2K