PRG-FmModulation

FM変調をプログラムで

http://ism1000ch.hatenablog.com/entry/2014/06/17/205452

class Vibrato(ef.Effecter):
   '''
   member
       depth : 変調具合[frame]
       freq  : 変調周波数[hz]
       rate  : サンプリングレート[hz]
       n     : 現在フレーム[frame]
   '''

   def __init__(self,depth=2, freq=2, rate=44100):
       self.depth = int(rate * depth / 1000) # input:[ms]
       self.freq  = freq
       self.rate  = rate 
       self.n     = 0

   def process(self,data):
       # 時間軸をゆがめる
       frames = self.calc_frames(np.arange(data.size))
      
       # 対応するシグナルを線形補完する
       data = self.calc_signal(frames,data)
       return data

   # 時間軸をゆがめる. N(n)の計算
   # input  : 時間軸[ 0,   1,   2, 3,..]
   # output : 時間軸[ 0, 0.5, 2.5, 3,..]
   def calc_frames(self,frames):
       vfunc = np.vectorize(self.calc_frame)
       return vfunc(frames)
       #return frames

   def calc_frame(self,n):
       # n = n ~ n + 2*depth
       n = n + self.depth * (1 + np.sin(self.n * (2 * np.pi * self.freq /　self.rate)))
       self.n += 1
       return n

   # N(n)を与え，線形補完してy(n)を計算
   def calc_signal(self,frames,data):
       # frames: N(n)のリスト
       # framesのlimit番目以降に対し処理を行う
       limit = self.depth * 2
       calc_data = [self.calc_interp(frame,data) for frame in frames[limit:]]
       calc_data = np.hstack([data[:limit],calc_data])
       return np.array(calc_data)

   def calc_interp(self,frame,data):
       x = int(np.floor(frame))
       d = np.interp(frame,[x,x+1],data[x:x+2]) #indexオーバーする可能性アリ
       return d

Tips

http://aidiary.hatenablog.com/entry/20110514/1305377659