python numpyにおけるarange関数のエラー

Question

print(np.arange(0.01,0.07,0.01))
print(np.arange(0.01,0.08,0.01))

を実行すると、
[ 0.01 0.02 0.03 0.04 0.05 0.06]
[ 0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.08]

のように出力されます。前者はいいとしても後者はなぜ0.08までの配列となっているのでしょうか。
[ 0.01 0.02 0.03 0.04 0.05 0.06 0.07] としたいので、解決方法があればご教授願いします。

Answer 1

コンピュータ上で浮動小数点数を扱う場合に出てくる一般的な動作仕様（問題）となります。
詳細はNumpyのarangeの動作はintとfloatで異なる？を参照ください。

今回のケースではnumpy.linspaceを用いることにより、たとえば２番目は
print(np.linspace(0.01,0.07,num = 7))または
print(np.linspace(0.01,0.08,num = 7, endpoint = False))
にて期待する結果を得られます。

Comments

[simasima]

ありがとうございます！arange のみではちょっと厳しそうですね。。。

[8524ba23]

他の回答にあるようにint型の数列をarangeで生成した後に浮動小数点数に変換してもよいです。
が、せっかく専用関数linspaceがあるので、そちらを使うほうが汎用/柔軟性は高いと思います。

Answer 2

結果だけであれば、

np.arange(0.01, 0.079, 0.01)

浮動小数の精度のせいでしょうか。

---

https://docs.scipy.org/doc/numpy-1.13.0/reference/generated/numpy.arange.html

When using a non-integer step, such as 0.1, the results will often not be consistent. It is better to use linspace for these cases.

Comments

[simasima]

ありがとうございます。できれば関数などを利用して一般的に解決をしたいのですが。。。

[mkgrei]

np.arange(1,8,1)/100.
でしょうか。

[mkgrei]

np.linspace(0.01, 0.07, 7)
でもよいみたいですね。

[mkgrei]

(0.08-0.01)/0.01が
7.00000000000000088817841970012523233890533447265625
になるせいですね。
これにceilがかかって8ステップ分計算されています。

from decimal import Decimal
Decimal(0.01)
などすると実際の計算精度が見られます。

[simasima]

勉強になります。ありがとうございます！

Answer 3

arangeのみで実装するならこんな形ではどうでしょうか？

import numpy as np
print(np.arange(1, 7) / 100)
print(np.arange(1, 8) / 100)

整数で処理すれば、arangeで問題ないことがわかっているので、それを100で割るだけのシンプル対応です。（mkgreiさんがコメントで書かれているものをさらに簡素化）

Comments

[simasima]

確かにそうですね！ありがとうございます！