我想确定音乐文件中的频率。它将从音乐文件中检索一大块数据并打印其中的频率。然后它将选择另一个数据块。有没有可能使用python实现它?我是这个领域的新人。如果有人能帮我做到这一点,我将不胜感激他。我的目标设备是树莓派。在import pyaudio
import wave
import numpy as np
chunk = 2048
wf = wave.open('/home/pi/music.wav', 'rb')
swidth = wf.getsampwidth()
RATE = wf.getframerate()
window = np.blackman(chunk)
p = pyaudio.PyAudio()
stream = p.open(format =
p.get_format_from_width(wf.getsampwidth()),
channels = wf.getnchannels(),
rate = RATE,
output = True)
data = wf.readframes(chunk)
while len(data) != '':
stream.write(data)
indata = np.array(wave.struct.unpack("%dh"%(len(data)/swidth),\
data))*window
fftdata=abs(np.fft.rfft(indata))**2
fftData=abs(np.fft.rfft(indata))**2
# find the maximum
which = fftData[1:].argmax() + 1
# use quadratic interpolation around the max
if which != len(fftData)-1:
y0,y1,y2 = np.log(fftData[which-1:which+2:])
x1 = (y2 - y0) * .5 / (2 * y1 - y2 - y0)
# find the frequency and output it
thefreq = (which+x1)*RATE/chunk
print "The freq is %f Hz." % (thefreq)
else:
thefreq = which*RATE/chunk
print "The freq is %f Hz." % (thefreq)
# read some more data
data = wf.readframes(chunk)
if data:
stream.write(data)
stream.close()
p.terminate()
我上面的程序对单频采样效果很好。但对于一首歌,它在indata=np.数组(wave.struct.解包(“%dh”%(len(data)/swidth),data))*窗口。error msg is操作数无法与形状一起广播。如何解决这个问题??在