C言語の計算結果が疑問です

Question

Cで作成されたソースで計算結果がよくわかりません

バイナリファイルを読み込んでR、c、d、E、Dを求めeとfを計算してます

printf("Rは %f\n",R);
printf("Eは %f\n",E);
printf("Dは %f\n",D);
printf("cは %d\n",c);
printf("dは %d\n",d);
e = (R + c * E) / D;
f = (R + d * E) / D;
printf("eは %f\n",e);
printf("fは %f\n",f);

出力結果
Rは 839.896912
Eは 0.250000
Dは 1.000000
cは 2480
dは 2474
eは 1459.896973
fは 1458.396973

電卓でもEXCELでも
eは1459.896912
fは1458.396912となりなぜ上記の計算結果になるのでしょうか？

Comments

[takabosoft]

各変数の「型」が判るようにソースを修正してください。

Answer 1

単精度浮動小数点数（float型）の有効数値は十進で約7桁です。つまり、8桁目以降はゴミです。
倍精度浮動小数点数（double型）にすれば十進で約15桁の有効数値を持ちます。

Comments

[marymills]

詳しくありがとうございます。桁数の違いで計算結果が異なると勉強になりました

Answer 2

C言語では、２進数により小数点以下を表現しているため、ニンゲンが計算する１０進数の計算結果とはずれが発生します
「浮動小数点 誤差」でぐぐってみればいいかと。

Answer 3

e と f が float で定義していませんかね。

gcc で試したらこうなりました。

  double R = 839.896912;
  double E = 0.250000;
  double D = 1.000000;
  int c = 2480;
  int d = 2474;

  printf("R is %f\n", R);
  printf("E is %f\n", E);
  printf("D is %f\n", D);
  printf("c is %d\n", c);
  printf("d is %d\n", d);

  double e = (R + c * E) / D;
  double f = (R + d * E) / D;
  printf("e is %f\n", e);
  printf("f is %f\n", f);

  float e2 = (R + c * E) / D;
  float f2 = (R + d * E) / D;
  printf("e2 is %f\n", e2);
  printf("f2 is %f\n", f2);

R is 839.896912
E is 0.250000
D is 1.000000
c is 2480
d is 2474
e is 1459.896912
f is 1458.396912
e2 is 1459.896973
f2 is 1459.396973

Comments

[marymills]

float とdoubleで計算結果が違うのですね。わかりました。
ありがとうございました。

Answer 4

どこかの変数がfloatになっている、ということはありませんか？（doubleとは精度が違うので、結果も異なってきます）