フィルターのクリア

How to tell if a random 3x3 Matrix is invertible

15 ビュー (過去 30 日間)
Kyle Donk
Kyle Donk 2020 年 1 月 10 日
編集済み: Meg Noah 2020 年 1 月 11 日
I have to produce a random 3x3 matrix A that is invertible and display it. I have a couple questions:
  1. How do I know when a matrix is invertible? I used the command (inv) on the random 3x3 matrix that I had created and I got a 3x3 matrix with different numbers. Does this mean that the matrix is invertible?
  2. I also got a hint with the question: Use a while-loop until you get one with non-zero determinant. I am confused by this because I used the determinants command (det) on my 3x3 matrix and got a nonzero determinant. I feel like I might be missing something here.
  2 件のコメント
Walter Roberson
Walter Roberson 2020 年 1 月 10 日
rank()
Most random matrices with floating point entries are invertible. Not all, but most. So unless you are using integer random values, do not be surprised if the first one you generate works.

サインインしてコメントする。

採用された回答

Meg Noah
Meg Noah 2020 年 1 月 10 日
編集済み: Meg Noah 2020 年 1 月 11 日
Back to your question, I have to produce a random 3x3 matrix A that is invertible and display it. One way could be to start with a matrix that you know will have a determinant of zero and then add random noise to each element. It worked for me to generate random matrices that are invertable.
for MC = 1:10000
% first create a matrix that you know has a low rcond value:
A = double(uint32(1000.*rand(3,1)).*uint32(1000.*rand(1,3)));
% then add noise
C = A + 100.0*rand(3,3);
if (rcond(C)<1e-20)
disp('algorithm fails');
C
inv(C)
end
end
There were objections to this suggestion about checking the determinant value. I had said: If the determinant of a square matrix is 0, it can't be inverted. I'd suggestion to test with - using your tolerance on the last argument. See comments below.
  7 件のコメント
Steven Lord
Steven Lord 2020 年 1 月 11 日
Better is simply to not invert a matrix. If you're trying to invert the matrix to solve a system of equations, use the backslash operator (\).
Meg Noah
Meg Noah 2020 年 1 月 11 日
編集済み: Meg Noah 2020 年 1 月 11 日
But that wasn't the question. He has a task to produce a matrix that can be inverted.

サインインしてコメントする。

その他の回答 (0 件)

カテゴリ

Help Center および File ExchangeSwitches and Breakers についてさらに検索

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by