Mejorando una matriz mal condicionada.
Tengo una matriz mal condicionada, cuyarcond()
es cercano a cero, y por lo tanto, la inversa de esa matriz no resulta ser correcta. He intentado usarpinv()
Pero eso no resuelve el problema. Así es como estoy tomando lo inverso:
X = (A)\(b);
Busqué una solución a este problema y encontréeste enlace (última solución) para mejorar la matriz. La solución allí sugiere usar esto:
A_new = A_old + c*eye(size(A_old));
Dóndec > 0
. Hasta ahora empleando esta técnica trabaja en hacer la matriz.A
Mejor acondicionado y la solución resultante se ve mejor. Sin embargo, investigué utilizando diferentes valores dec
y la solución resultante depende del valor elegido.c
.
Aparte de investigar manualmente el valor dec
, ¿hay una forma automática a través del cual puedo encontrar el valor dec
¿Para cuál obtengo la mejor solución?