如何在单个图形上绘制声波、mfcc 和 mel 频谱图？

Question

我有一个包含多个 wav 文件的文件夹（目前说 4 个）。我想连续绘制 wav、它的 mfcc 和 mel 频谱图，所以最后是一个有 12 个图的图（每行有三个图，因此有四行）。我无法绘制图表，只能提取特征。请有人帮忙解决这个for循环。我的意思是如何使用 subplot 命令以及如何将每个图形存储在循环中。

问候

path=glob.glob('the path having four wav files/*.wav')

for p in path:
    y, sr = librosa.load(p, sr=16000)
    mfcc=librosa.feature.mfcc(y)
    S = librosa.feature.melspectrogram(y, sr)
    fig, ax = plt.subplot(4,3,.....)
    librosa.display.waveplot(y, sr=sr)
    librosa.display.specshow(librosa.power_to_db(S, ref=np.max))
    librosa.display.specshow(mfcc, x_axis="time",y_axis="mel")
    plt.show()

Answer 1

import matplotlib.pyplot as plt
Arrays=np.random.rand(6,10,10)
N_figures=Arrays.shape[0]
for Id in range(1,N_figures+1):
    Array=Arrays[Id-1]
    plt.subplot(N_figures/2,N_figures/3,Id)
    plt.imshow(Array)

您只需要更改 imshow 所需的绘图。因为我没有频谱图文件，所以我使用了随机创建的 NumPy 数组。您的实现不起作用，因为fig, ax = plt.subplot(4,3,.....) 不适合使用。 Fig 和 axs 在函数 plt.subplots 而不是 plt.subplot 中返回。您也可以使用该方法，但在 for 循环之外，然后访问 axs 元素。

Answer 2

最终代码为：

import matplotlib.pyplot as plt
import librosa
import librosa.display
import glob

path=glob.glob('E:/Python_On_All_Dataset/Four emotion_for plot/*.wav') 

fig, ax =  plt.subplots(nrows=4, ncols=3, sharex=True)

    
for i in range(4) :
   
    y, sr = librosa.load(path[i], sr=16000)
    librosa.display.waveplot(y, sr, ax=ax[i, 0])  # put wave in row i, column 0
    plt.axis('off')
    
    
   
    mfcc=librosa.feature.mfcc(y) 
    librosa.display.specshow(mfcc, x_axis='time', ax=ax[i, 1]) # mfcc in row i, column 1
   

    S = librosa.feature.melspectrogram(y, sr)
    librosa.display.specshow(librosa.power_to_db(S), x_axis='time', y_axis='log', ax=ax[i, 2])  # spectrogram in row i, column 2
   

尝试放置这个轴（在每次绘图后关闭，但不知何故它不起作用）