Кратковременное преобразование Фурье в python

Я хочу получить частоту с максимальной мощностью в каждый момент в файле wav. Итак, я написал STFT на Python, используя fft из scipy. Я использовал оконную функцию kaiser от scipy. Все выглядит отлично, но мои результаты выглядят странно. У него очень маленькие числа, а некоторые очень высокие.

вот результат для одного файла wav: http://pastebin.com/5Ryd2uXj и вот код на Python:

import scipy, pylab
import wave
import struct
import sys

def stft(data, cp, do, hop):
    dos = int(do*cp)
    w = scipy.kaiser(dos,12) //12 is very high for kaiser window
    temp=[]
    wyn=[]
    for i in range(0, len(data)-dos, hop):
        temp=scipy.fft(w*data[i:i+dos])
        max=-1
        for j in range(0, len(temp),1):
            licz=temp[j].real**2+temp[j].imag**2
            if( licz>max ):
                max = licz
                maxj = j
        wyn.append(maxj)
    #wyn = scipy.array([scipy.fft(w*data[i:i+dos])
        #for i in range(0, len(data)-dos, 1)])
    return wyn

file = wave.open(sys.argv[1])
bity = file.readframes(file.getnframes())
data=struct.unpack('{n}h'.format(n=file.getnframes()), bity)
file.close()

cp=44100 #sampling frequency
do=0.05 #window size
hop = 5

wyn=stft(data,cp,do,hop)
print len(wyn)
for i in range(0, len(wyn), 1):
    print wyn[i]

5
задан user1226419 22 February 2012 в 17:10
поделиться