python 实现wav的波形显示（时域和频域）

　　音频处理中，经常要看一下啊频域图是什么样子的，这里自己写了一个小程序，可以完美的同步显示时域和频域图，直接上代码：

 #wave data   -xlxw

 #import
 import wave as we
 import numpy as np
 import matplotlib.pyplot as plt
 import sys 

 def wavread(path):
     wavfile =  we.open(path,"rb")
     params = wavfile.getparams()
     framesra,frameswav= params[2],params[3]
     nchannels, sampwidth, framesra, frameswav = params[:4]
     print("nchannels:%d" % nchannels)
     print("sampwidth:%d" % sampwidth)
     datawav = wavfile.readframes(frameswav)
     wavfile.close()
     datause = np.fromstring(datawav,dtype = np.short)
     print(len(datause))
     if nchannels == 2:
         datause.shape = -1,2
     datause = datause.T
     time = np.arange(0, frameswav) * (1.0/framesra)
     return datause,time,nchannels

 def main():
     path = sys.argv[1]
     #path = input("The Path is:")
     print(path)
     wavdata,wavtime,nchannels = wavread(path)

     N=len(wavdata)
     framerate = 16000
     start=0
     df = 1
     freq = [df*n for n in range(0,len(wavdata))]
     print(len(wavdata))
     print(len(wavtime))

     c=np.fft.fft(wavdata)*nchannels
     d=int(len(c)/2)
     print(len(c))

     fig, ax = plt.subplots(2, 1)

     ax[0].plot(wavtime,wavdata,color = 'green')
     ax[0].set_xlabel('Time')
     ax[0].set_ylabel('Amplitude')

     ax[1].plot(freq,abs(c),color = 'red')
     ax[1].set_xlabel('Freq(HZ)')
     ax[1].set_ylabel('Y(freq)')

     plt.show()

　直接上图看结果：

　　这个只能对单声道１6k采样的wav格式做处理，不过，只要稍微加一更改，就可以处理别的了。具体改代码的事情，还是谁用谁做吧。做个程序员，总要付出一些代价的。

每日一言：君子役物，小人役于物。——《劝学》