2017-05-05から1日間の記事一覧

高速フーリエ変換

# coding:utf-8 import numpy as np from pylab import * import pdb f = np.matrix([0,1,1,2,3,2,3,2]).transpose() plot(f) show() w8 = e**(pi/4*1j) #1の原始8乗根 w4 = e**(pi/2*1j) #1の原始4乗根 w2 = e**(pi*1j) #1の原始2乗根 a0 = f[0] a1 = f[1] …

離散フーリエ変換

離散フーリエ変換を1の原始N乗根による表現で計算する1の原始N乗根の定義:N=8の場合 import numpy as np from pylab import * import pdb f = np.matrix([0,1,1,2,3,2,3,2]).transpose() plot(f) show() w8 = e**(pi/4*1j) A = 1/8.0 * np.matrix([ [1,1,1,…