TI-Nspire CX CAS Lua二次方程程序在电脑与计算器上结果不一致问题排查

我之前也碰到过类似的浮点精度坑！你遇到的这个问题其实是浮点数抵消误差加上计算器和电脑的浮点运算环境差异导致的，咱们一步步拆解：

为什么会出现差异？ #

电脑和TI-Nspire计算器的处理器架构、Lua解释器的浮点实现细节不一样，当计算两个非常接近的数值做加减（比如这里的-s + sqrt(s²-4ft)，如果s是正数且判别式很小，sqrt(s²-4ft)会非常接近s，两者相加后其实是一个极小的数），就会触发灾难性抵消——浮点数的有效位数不够，导致结果变成类似-4.44089...e-16的舍入误差，这个值本质上是0的近似，但计算器的显示逻辑把它当成了非零的极小值。

而电脑的浮点运算单元精度更高，或者刚好在计算时的舍入方向不同，所以没有显示出这个误差。

解决办法 #

有两个核心思路可以解决这个问题：

1. 使用避免抵消的二次方程公式变形 #

利用韦达定理（二次方程两根乘积为c/a），先计算绝对值较大的那个根，再用乘积关系算出另一个根，彻底避免相近数相减的场景：

• 如果b是正数，先算(-b - sqrt(b²-4ac))/(2a)（这个根的绝对值更大），再用c/(a*root)得到另一个根；
• 如果b是负数，先算(-b + sqrt(b²-4ac))/(2a)，再用同样的韦达定理推导式得到另一个根。

2. 对极小值做近似处理 #

判断计算结果的绝对值是否小于某个极小阈值（比如1e-12），如果是就直接返回0，避免显示科学计数法的极小值。

改进后的代码示例 #

-- 改进后的二次方程求解函数，规避浮点抵消误差
function quadraticRoots(a, b, c)
    a = tonumber(a)
    b = tonumber(b)
    c = tonumber(c)
    local discriminant = b^2 - 4*a*c
    local sqrtDisc = math.sqrt(discriminant)
    
    local root1, root2
    -- 根据b的符号选择先计算的根，避免抵消
    if b >= 0 then
        root1 = (-b - sqrtDisc) / (2*a)
        -- 用韦达定理变形计算第二个根，避免重复计算可能带来的误差
        root2 = (2*c) / (-b - sqrtDisc)
    else
        root1 = (-b + sqrtDisc) / (2*a)
        root2 = (2*c) / (-b + sqrtDisc)
    end
    
    -- 把极小值近似为0，优化显示效果
    local epsilon = 1e-12
    if math.abs(root1) < epsilon then root1 = 0 end
    if math.abs(root2) < epsilon then root2 = 0 end
    
    return root1, root2
end

function on.paint(gc)
    local r1, r2 = quadraticRoots(a, b, c)
    formula:setExpression("0s: "..r1)
    formulaB:setExpression(r2)
end

效果验证 #

用这个改进后的函数，不管是在电脑还是TI-Nspire计算器上，都会得到更准确的结果，不会再出现那个恼人的极小科学计数法数值了。

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