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格式做处理，不过，只要稍微加一更改，就可以处理别的了。具体改代码的事情，还是谁用谁做吧。做个程序员，总要付出一些代价的。

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