Azzera filtri
Azzera filtri

How do I return something only if all steps of an conditional statement is true?

2 visualizzazioni (ultimi 30 giorni)
I am attempting to generate a random 4x4 matrix and if all the real parts of the eigenvalues are negative, then compute the inverse matrix of the original random matrix.
m=4;
G=randn(m);
H=eig(G);
J=real(H);
[rowsJ,colsJ]=size(J);
for r=1:rowsJ
if J(r,1)<1
K=G^-1
else
printf("Unable to Compute")
end
end
  4 Commenti

Accedi per commentare.

Risposta accettata

Walter Roberson
Walter Roberson il 12 Set 2019
if all(J(:)<1)
K=inv(G);
else
K = [];
printf("Unable to Compute")
end
However... you should rarely be computing the inverse of a matrix. In most cases you should be using other approaches than inv(), such as using the \ operator.
  4 Commenti
Stephen23
Stephen23 il 13 Set 2019
Modificato: Stephen23 il 13 Set 2019
"but why should I not be computing the inverse of a matrix?"
inv is slower and less accurate than the recommended methods of solving systems of linear equations.
The inv documentation explains: "It is seldom necessary to form the explicit inverse of a matrix. A frequent misuse of inv arises when solving the system of linear equations Ax = b. One way to solve the equation is with x = inv(A)*b. A better way, from the standpoint of both execution time and numerical accuracy, is to use the matrix backslash operator x = A\b."

Accedi per commentare.

Più risposte (0)

Categorie

Scopri di più su Linear Algebra in Help Center e File Exchange

Community Treasure Hunt

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

Start Hunting!

Translated by