pythonの０での割り算

Question

発生している問題・エラーメッセージ

以下のようなエラーが出てしまっていて、ゼロで割ってしまっているとのことなのですがどこでゼロで割っているのかがわからなくて困っています。

Traceback (most recent call last):
  File "/Users/matsumura/PycharmProjects/wave/pinknoise_generater.py", line 29, in <module>
    a, b = synth_pinknoise(duration=10.0, samprate=48000)
  File "/Users/matsumura/PycharmProjects/wave/pinknoise_generater.py", line 17, in synth_pinknoise
    t = np.arange(start=0, stop=duration, step=(int(1/samprate)))
ZeroDivisionError: float division by zero

該当のソースコード

# coding: UTF-8
# 必要なライブラリ
import math
import numpy as np
import matplotlib
matplotlib.use('TkAgg')
import matplotlib.pyplot as plt

# ピンクノイズの関数定義
def synth_pinknoise(duration, samprate):
    """
    Generate pinknoise using Voss algorithm
    http://www.firstpr.com.au/dsp/pink-noise/
    """
    f_low = 10   # lowest frequency to keep pink (Hz)
    levels = math.ceil(np.log2(samprate/f_low))
    t = np.arange(start=0, stop=duration, step=1/samprate)
    out = np.zeros(len(t))
    x = np.random.normal(size=levels)
    for n in range(math.ceil(samprate * duration)):
        for m in range(levels):
            if n % 2**(m+1) == 0:
                x[m] = np.random.normal()
            out[n] = np.random.normal() + np.sum(x)
    return(out / max(abs(out)), t)


# 使い方（戻り値は、x：ピンクノイズ信号、t：各サンプルの時刻）
a, b = synth_pinknoise(duration=10.0, samprate=48000)
plt.plot(a,b)
plt.show()

Comments

[hayataka2049]

手元環境では再現しません。念の為、pythonおよびnumpyのバージョンをご提示ください。

[hayataka2049]

あ、と思ったけどエラーが出ているコードと貼り付けているコードは微妙に違いますね。エラーの記述に合わせたら再現できたので、回答します。

Answer 1

せっかくエラー・メッセージで行を指定してくれているんですからそれを見ようよ。

で、

t = np.arange(start=0, stop=duration, step=(int(1/samprate)))

このstepはいくつになるか考えてみよう

Answer 2

Traceback (most recent call last):
  File "/Users/matsumura/PycharmProjects/wave/pinknoise_generater.py", line 29, in <module>
    a, b = synth_pinknoise(duration=10.0, samprate=48000)
  File "/Users/matsumura/PycharmProjects/wave/pinknoise_generater.py", line 17, in synth_pinknoise
    t = np.arange(start=0, stop=duration, step=(int(1/samprate)))
ZeroDivisionError: float division by zero

step=(int(1/samprate))でstepに0を与えることになるので、np.arangeの内部でゼロ除算エラーになります。内部的にstepで割られるような計算があるのです（arangeの結果の出し方を想像するとすぐわかると思います）。

Comments

[sibazyun]

ご返事ありがとうございます。載せたコードとエラーが違い申し訳ございません、、、。凡ミスでした。
これはつまりstep=1/samprateとすればstepが小数点になり、0ではなくなるのでエラーではなくなるということでしょうか？

[hayataka2049]

そうですね。その場合はtは0から10.0の区間（厳密にはちょっと微妙ですが）を等分した480000要素の配列になることでしょう。

[sibazyun]

step=1/samprateとしたのですがやはり同じようなエラーが出てしまいました、、、。
```
Traceback (most recent call last):
File "/Users/matsumura/PycharmProjects/wave/pinknoise_generater.py", line 31, in <module>
a, b = synth_pinknoise(duration=10.0, samprate=48000)
File "/Users/matsumura/PycharmProjects/wave/pinknoise_generater.py", line 19, in synth_pinknoise
t = np.arange(start=0, stop=duration, step=1/samprate)
ZeroDivisionError: float division by zero
```

[sibazyun]

pythonのアップデートを行ったらうまく行きました！