相同代码在MATLAB 2016与2021版本运行为何得到不同结果

你遇到的差异属于single精度下的正常浮点舍入误差，根源是2016和2021版MATLAB对混合精度运算的隐式转换逻辑、sum底层的累加实现（如并行累加顺序、SIMD指令优化策略）不同导致，两个结果的差值在single类型的精度允许范围内。可以通过以下方案避免这类差异：

• 显式统一运算精度，避免隐式混合精度转换。不要依赖MATLAB版本默认的混合精度转换规则，运算前手动将两个变量转为同一精度：
  统一为single运算：sum(single(x1) .* b)
  统一为double运算：sum(double(b) .* x1)
• 强制指定累加过程的精度，降低累加误差的版本差异。可以先将所有变量转为double完成累加计算后，再转换回single，全程累加过程精度一致，不受版本底层优化影响：
  示例代码：single(sum(double(b) .* x1))
• 关闭版本相关的浮点运算优化。如果不需要最高运算速度，可以关闭MATLAB的JIT加速和SIMD向量化优化，避免不同版本优化逻辑带来的运算顺序差异：
  执行命令：feature('Accelerate', 'off')
• 自定义固定顺序的累加逻辑。对可复现性要求极高的场景，可以放弃内置sum函数，手动按固定索引顺序循环累加，保证所有版本的累加顺序完全一致，浮点误差也会完全相同。

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