rcond

Arguments

X

matrice carrée réelle ou complexe

r

nombre réel positif

Description

rcond(X) est une estimation de l'inverse du conditionnement de X pour la norme L_1.

Si X est bien conditionnée, rcond(X) est proche 1. Sinon, rcond(X) est proche de 0.

L'estimation de l'inverse du conditionnement L-1 avec rcond est beaucoup (plus de 200x) plus rapide que cond(X). En contrepartie, rcond peut être moins fiable.

Algorithme : la norme L-1 de X est calculée avec Lapack/DLANGE. La décomposition LU de X est effectuée avec Lapack/DGETRF. Finalement, le conditionnement est estimé avec Lapack/DGECON.

rcond([]) vaut %inf.

Exemples

A = diag([1:10]);
rcond(A)
A(1,1) = 0.000001;
rcond(A)

Bancs-tests comparés cond() / rcond() :

A = ones(1000, 1000);
timer(); cond(A);    timer()
timer(); 1/rcond(A); timer()

Voir aussi

• svd — décomposition en valeurs singulières
• cond — conditionnement
• inv — inverse d'une matrice

Historique