Inversione matrici

Numero di condizionamento di una matrice

con numeric.js

function conditionNumberNumeric(A) {
    // numeric.svd restituisce { U, S, V }
    const svd = numeric.svd(A);
    const s = svd.S;
    
    const s_max = Math.max(...s);
    const s_min = Math.min(...s);
    
    if (s_min < 1e-15) {
        return {
            conditionNumber: Infinity,
            s_max: s_max,
            s_min: s_min,
            isSingular: true,
            message: 'Matrice singolare o mal condizionata'
        };
    }
    
    return {
        conditionNumber: s_max / s_min,
        s_max: s_max,
        s_min: s_min,
        isSingular: false,
        log10Condition: Math.log10(s_max / s_min)
    };
}
const A = [
            [4, 1, 2],
            [2, 5, 1],
            [1, 0, 6]
        ];
conditionNumberNumeric(A)