转:使用Python对音频进行频谱分析
2017-12-01 08:34
141 查看
#转自 http://blog.sina.com.cn/s/blog_40793e970102w3m2.html[/code]import wave import pyaudio import numpy import pylab #打开WAV文档,文件路径根据需要做修改 wf = wave.open("test.wav", "rb") #创建PyAudio对象 p = pyaudio.PyAudio() stream = p.open(format=p.get_format_from_width(wf.getsampwidth()), channels=wf.getnchannels(), rate=wf.getframerate(), output=True) nframes = wf.getnframes() framerate = wf.getframerate() #读取完整的帧数据到str_data中,这是一个string类型的数据 str_data = wf.readframes(nframes) wf.close() #将波形数据转换为数组 # A new 1-D array initialized from raw binary or text data in a string. wave_data = numpy.fromstring(str_data, dtype=numpy.short) #将wave_data数组改为2列,行数自动匹配。在修改shape的属性时,需使得数组的总长度不变。 wave_data.shape = -1,2 #将数组转置 wave_data = wave_data.T #time 也是一个数组,与wave_data[0]或wave_data[1]配对形成系列点坐标 #time = numpy.arange(0,nframes)*(1.0/framerate) #绘制波形图 #pylab.plot(time, wave_data[0]) #pylab.subplot(212) #pylab.plot(time, wave_data[1], c="g") #pylab.xlabel("time (seconds)") #pylab.show() # # 采样点数,修改采样点数和起始位置进行不同位置和长度的音频波形分析 N=44100 start=0 #开始采样位置 df = framerate/(N-1) # 分辨率 freq = [df*n for n in range(0,N)] #N个元素 wave_data2=wave_data[0][start:start+N] c=numpy.fft.fft(wave_data2)*2/N #常规显示采样频率一半的频谱 d=int(len(c)/2) #仅显示频率在4000以下的频谱 while freq[d]>4000: d-=10 pylab.plot(freq[:d-1],abs(c[:d-1]),'r') pylab.show()
相关文章推荐
- 使用Python对音频进行频谱分析
- 使用JMF进行视频、音频的采集与实时传送
- 使用ropemacs对python代码进行补全
- 使用 Python 进行线程编程(转)
- 使用SWIG对Python进行扩展
- 使用 Python 进行线程编程
- 使用python+cron对php状态进行定时监控
- 使用 Google API 进行翻译的 Python 小工具
- 使用NetBeans IDE进行python编程: ImportError: No module named pygtk
- 使用python进行下载地址转换
- 使用JMF进行视频、音频的采集与实时传送
- 在WinXP平台使用Directshow进行音频录制
- 在WinCE平台使用Directshow进行音频录制
- 使用Python进行二进制文件读写
- 使用 Python 进行线程编程
- 使用SWIG和Python对C/C++进行单元测试(一)
- 使用python对txt格式的小说进行处理
- 由于Python本身自带的界面库功能并不强大,我们使用Python+wxPython作为界面开发平台,在公司新的项目中进行界面开发。开发过程使用Eclipse+PyDev 作为集成开发环境。产品发布时使用Py2exe进行打包。如今,在我的团队中,Pytho
- windows下配置eclipse + Python + PyDev + tornado,并使用pylint进行静态检测
- 使用SWIG和Python对C/C++进行单元测试(二)