畳み込みリバーブ

convolution_reverb.py

import numpy as np
from scipy import signal
import matplotlib.pyplot as plt
import soundfile as sf

indata, infrq = sf.read("indata.wav")
impulse, impfrq = sf.read("impulse.wav")

outdata_0 = signal.convolve(indata[:, 0], impulse[:, 0], mode="full")
outdata_1 = signal.convolve(indata[:, 1], impulse[:, 1], mode="full")
outdata = np.array([outdata_0, outdata_1]).T

sf.write("outdata.wav", outdata, infrq)

fig = plt.figure(1)
plt.subplot(3, 1, 1)
plt.plot(indata[:, 0], "b-")

plt.subplot(3, 1, 2)
plt.plot(impulse[:, 0], "r-")

plt.subplot(3, 1, 3)
plt.plot(outdata[:, 0], "g-")

plt.savefig("grph.png")
plt.show()