使用Maxima的mat_norm计算矩阵2范数报错，求问题原因及解决方法

矩阵存在的问题 #

你定义的矩阵是数值奇异矩阵，计算它的行列式：

determinant(A);

结果约为-1e-8，在浮点数精度下可视为0，说明矩阵的秩为1，属于奇异矩阵。Maxima内置的mat_norm函数计算2范数时，内部依赖特征值求解逻辑，而奇异矩阵的数值特征值计算会因精度问题触发报错，导致无法完成计算。

Maxima中计算矩阵2范数的推荐方法 #

方法1：利用奇异值分解（最可靠） #

矩阵的2范数等于其最大奇异值，直接使用singular_values函数获取：

A : matrix([0.780, 0.563],[0.913, 0.659]);
sv_list : singular_values(A);
max(sv_list);

运行结果即为矩阵的2范数。

方法2：手动按定义计算 #

根据2范数的定义，它是矩阵转置与自身乘积（A^T·A）的最大特征值的平方根，步骤如下：

A : matrix([0.780, 0.563],[0.913, 0.659]);
A_trans = transpose(A);
A_T_A = A_trans . A;
eigen_vals = eigenvalues(A_T_A);
sqrt(max(abs(eigen_vals)));

加入abs是为了处理浮点数精度导致的微小负特征值，避免开平方错误。

方法3：尝试加载linalg包使用norm函数 #

部分Maxima版本中，加载linalg扩展包后可直接调用norm函数：

load(linalg);
A : matrix([0.780, 0.563],[0.913, 0.659]);
norm(A, 2);

若此方法仍报错，优先选择前两种方案。

内容的提问来源于stack exchange，提问作者user32330751