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格式做處理，不過，只要稍微加一更改，就可以處理別的了。具體改代碼的事情，還是誰用誰做吧。做個程序員，總要付出一些代價的。

每日一言：君子役物，小人役於物。——《勸學》