System.Windows.Media.Matrix双精度疑问：行列式非零却无法求逆？

这个问题其实是WPF框架在数值稳定性上的设计选择，不是double精度范围的锅，听我给你拆解：

核心原因：数值稳定性优先于严格数学定义 #

数学上只要行列式不为零，矩阵就可逆，但在计算机浮点数运算里，极小的行列式会导致逆矩阵的元素变得极大，带来严重的数值误差——甚至可能超出double的有效精度范围，得到完全不可靠的结果。

WPF的Matrix类为了避免这种情况，并没有严格判断行列式是否为零，而是用了一个阈值来判定：当行列式的绝对值小于这个阈值时，就认为矩阵“不可逆”，拒绝执行逆变换。

你的两个案例对比 #

• 案例1：行列式是1e-14，这个值大于WPF内部使用的阈值（具体是1e-10，可以看DoubleUtil.IsZero的实现），所以HasInverse返回true，Invert()可以正常执行。
• 案例2：行列式是~1e-15，这个值小于1e-10的阈值，所以HasInverse返回false，调用Invert()就会抛出异常。

为什么double范围足够还是有这个限制？ #

double的范围确实能容纳很大的数，但问题不在范围，而在精度。比如案例2里，逆矩阵的元素会是原矩阵元素除以行列式：1e-7 / 1e-15 = 1e8，看起来没问题，但如果行列式更小（比如1e-30），逆矩阵元素会到1e30，接近double的上限，而且浮点数的相对误差会被放大，最终的变换结果完全不可信。WPF作为UI框架，需要保证变换的稳定性，所以用阈值过滤了这种边缘情况。

如果你一定要处理这种矩阵 #

如果你的场景能接受数值误差，可以手动计算逆矩阵，绕过WPF的阈值判断：

double det = mat.Determinant;
// 这里可以用更小的阈值，比如double.Epsilon，根据你的需求调整
if (Math.Abs(det) > double.Epsilon)
{
    var invMat = new Matrix(
        mat.M22 / det, -mat.M12 / det,
        -mat.M21 / det, mat.M11 / det,
        (mat.M21 * mat.OffsetY - mat.M22 * mat.OffsetX) / det,
        (mat.M12 * mat.OffsetX - mat.M11 * mat.OffsetY) / det
    );
    // 使用invMat做后续操作
}

但要注意，这种方式得到的逆矩阵可能存在较大误差，一定要根据你的场景评估风险。

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