pythonでバイノーラルビート

経緯

音声処理の練習にと思い挑戦して見ました。

wavファイルに書き込む方法がわからないという致命的な問題がありますが、なんとか再生まではできるようにはしました(汗)

ソースコード

#coding: utf-8
import struct, pyaudio as pa
from math import sin, pi

class StereoSound:

    def __init__(self, format=pa.paInt16, sampling=44100):
        self.form = format
        self.samp = int(sampling)

    def generateSinWave(self, left, right, volume=0.98, time=1.0):
        cast = {pa.paFloat32:(float, 1.0, 'f'),
                pa.paInt32:(int, 2147483647.0, 'i'),
                pa.paInt24:(int, 8388607.0, '3b'),
                pa.paInt16:(int, 32767.0, 'h'),
                pa.paInt8:(int, 127.0, 'b'),}
        length = int(time * self.samp)
        num = length << 1
        data = [0.0] * num
        left_theta = 2 * pi * left / self.samp
        right_theta = 2 * pi * right / self.samp
        for (n, l), r in zip(enumerate(range(0, num,2)), range(1, num,2)):
            data[l] = volume * sin(left_theta * n)
            data[r] = volume * sin(right_theta * n)
        data = [cast[self.form][0](x * cast[self.form][1]) for x in data]
        self.data = struct.pack(cast[self.form][2] * num, *data)

    def play(self, repeat=1):
        p = pa.PyAudio()
        stream = p.open(format=self.form,
                        channels=2,
                        rate=self.samp,
                        output=True)
        while repeat > 0:
            stream.write(self.data)
            repeat -= 1
        stream.close()
        p.terminate()

if __name__ == "__main__" :
    sound = StereoSound(pa.paFloat32, 192000)
    f1 = 440
    f2 = 450
    sound.generateSinWave(f1, f2, 0.05)
    sound.play(5)

github.com だいぶ見にくい気がしますが、どうやったらわかりやすくかけるんでしょうか？

課題

とりあえず、wavファイルとして保存できるようにしたいですね。waveモジュールの使い方がわからずに詰んでますが。

環境音とかを合成できるようにして、最終的にはheroku,djangoでウェブアプリケーションしたいところです。

どっちも触ったこと無いので、だいぶ先になりそうですがちょびちょび進めていくつもりです。

参考

aidiary.hatenablog.com この記事にだいぶ助けられました。これが原型になってます。