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

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