MFCC功率谱图的绘制

Question

有没有人能帮我把mfcc特征画成功率谱图！我只能绘制mfcc系数，它不是在时域中表示的。我想在时域中绘制mfcc特征图。我还尝试了应用DCT将其转换为时域！

from python_speech_features import mfcc
#from python_speech_features import logfbank `enter code here`
import scipy.io.wavfile as wav
import pandas as pd
import matplotlib.pyplot as plt
from python_speech_features import logfbank
import scipy.fftpack

(rate,sig) = wav.read("Voice0003.wav")
mfcc_feat = mfcc(sig,rate)


(rate2,sig2) = wav.read("Voice0004.wav")
mfcc_feat2 = mfcc(sig2,rate2)
yf = scipy.fftpack.fft(sig,rate)
#fbank_feat = logfbank(yf)
#dct=scipy.fftpack.dct(fbank_feat      tried converting into Time Domain            didnt help


print(len(mfcc_feat))
print(len(mfcc_feat2))

pd.DataFrame(mfcc_feat2).T.plot()
plt.show()
pd.DataFrame(mfcc_feat).T.plot()
plt.show()

Answer 1

Pxx，freqs，bin，im = plt.specgram(signal，NFFT=NFFT，Fs=2，noverlap=100，cmap=None) #signal是.wav文件的一维数组