C#移植MATLAB代码时double类型数值小数点后结果不准确求助

该问题由二进制浮点数的固有精度缺陷导致：

• 你用到的0.025、0.075这类十进制小数，无法使用基于二进制存储的double类型完全精确表示，每次执行+= incxy累加操作时，微小的舍入误差会持续累积，多次累加后误差就会放大到可观测的程度。
• MATLAB默认对浮点数输出做了截断处理，你看到的0.025这类值只是展示层面做了简化，内部实际也存在微小误差；而C#默认会输出浮点数的全部有效位，因此误差会直接暴露。

可根据业务场景选择以下方案：

• 方案1：改用整数做步长计数，单次计算浮点值（最稳妥，无累加误差）
  放弃累加浮点值的写法，改用整数记录循环次数，每个坐标值直接通过初始值和步长倍数单次计算，误差仅产生一次不会累积，示例代码：

  // 计算x方向总步数，用四舍五入避免边界值误差
  int xTotalSteps = (int)Math.Round((x2 - x1) / incxy);
  for (int xStep = 0; xStep <= xTotalSteps; xStep++)
  {
      double xgriglia = x1 + incxy / 2 + xStep * incxy;
      // 同理处理y方向
      int yTotalSteps = (int)Math.Round((y2 - y1) / incxy);
      for (int yStep = 0; yStep <= yTotalSteps; yStep++)
      {
          double ygriglia = y1 + incxy / 2 + yStep * incxy;
          contagriglia++;
          griglia0[contagriglia - 1] = xgriglia;
          griglia1[contagriglia - 1] = ygriglia;
      }
  }
  

• 方案2：按业务需求保留固定小数位
  如果你的场景仅需要精确到3位小数，可在赋值时对结果做截断：

  griglia0[contagriglia - 1] = Math.Round(xgriglia, 3);
  griglia1[contagriglia - 1] = Math.Round(ygriglia, 3);
  

• 方案3：改用decimal类型计算和存储
  如果对精度要求极高且性能压力不大，可将所有相关变量替换为decimal类型，decimal为十进制浮点数，可精确表示0.025这类10的负幂次相关的小数，不会出现二进制浮点的精度问题。

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