我正在尝试绘制Amplitude (dBFS) vs. Time (s)音频图（.wav）文件使用matplotlib。我设法用以下代码做到了这一点：
def convert_to_decibel(sample):
ref = 32768 # Using a signed 16-bit PCM format wav file. So, 2^16 is the max. value.
if sample!=0:
return 20 * np.log10(abs(sample) / ref)
else:
return 20 * np.log10(0.000001)
from scipy.io.wavfile import read as readWav
from scipy.fftpack import fft
import matplotlib.pyplot as gplot1
import matplotlib.pyplot as gplot2
import numpy as np
import struct
import gc
wavfile1 = '/home/user01/audio/speech.wav'
wavsamplerate1, wavdata1 = readWav(wavfile1)
wavdlen1 = wavdata1.size
wavdtype1 = wavdata1.dtype
gplot1.rcParams['figure.figsize'] = [15, 5]
pltaxis1 = gplot1.gca()
gplot1.axhline(y=0, c="black")
gplot1.xticks(np.arange(0, 10, 0.5))
gplot1.yticks(np.arange(-200, 200, 5))
gplot1.grid(linestyle = '--')
wavdata3 = np.array([convert_to_decibel(i) for i in wavdata1], dtype=np.int16)
yvals3 = wavdata3
t3 = wavdata3.size / wavsamplerate1
xvals3 = np.linspace(0, t3, wavdata3.size)
pltaxis1.set_xlim([0, t3 + 2])
pltaxis1.set_title('Amplitude (dBFS) vs Time(s)')
pltaxis1.plot(xvals3, yvals3, '-')
给出以下输出：
我还绘制了Power Spectral Density (PSD, in dBm)使用下面的代码：
from scipy.signal import welch as psd # Computes PSD using Welch's method.
fpsd, wPSD = psd(wavdata1, wavsamplerate1, nperseg=1024)
gplot2.rcParams['figure.figsize'] = [15, 5]
pltpsdm = gplot2.gca()
gplot2.axhline(y=0, c="black")
pltpsdm.plot(fpsd, 20*np.log10(wPSD))
gplot2.xticks(np.arange(0, 4000, 400))
gplot2.yticks(np.arange(-150, 160, 10))
pltpsdm.set_xlim([0, 4000])
pltpsdm.set_ylim([-150, 150])
gplot2.grid(linestyle = '--')
输出如下：
上面的第二个输出使用韦尔奇方法绘制了更美观的输出。 dBFS 图虽然信息丰富，但在我看来并不是很美观。这是因为：
另外，有没有办法可以绘制我的dBFS输出为峰峰式的情节就像在我的PSD (dBm)情节而不是密集茎图?
将会很有帮助，并且会感谢这里专家的任何指示、答案或建议，因为我只是一个初学者matplotlib和情节python一般来说。
TLNR
pyplot.你应该构建数据，计算每帧的 RMS 或峰值，然后将该值转换为 dbFS而不是将这种变换应用于每个样本点。
当我们谈论振幅时，我们谈论的是周期信号。当我们从声音文件中读取一系列数据时，我们读取了一系列样本点信号的（可以是或不是周期性的）。每个采样点的值代表在特定时间采样的电压值或声压值。
We assume在很短的时间间隔内，例如 10 毫秒，信号是静止的。每个这样的区间称为frame.
通常会对每一帧应用一些特定的函数，以减少该帧边缘的突变，这些函数称为窗函数。如果您对每个帧不执行任何操作，则会向它们添加矩形窗口。
举个例子：当你的声音采样频率是44100Hz时，在10ms长的帧中，有44100*0.01=441样本点。这就是nperseg争论意味着在你的psd功能，但与 dbFS 无关。
有了上面的知识，现在我们可以讨论振幅了。
有两种方法可以获取每帧的幅度值：
之后，峰值或RMS值可以转换为dbFS值。
让我们开始编码：
import numpy as np
import matplotlib.pyplot as plt
from scipy.io import wavfile
# Determine full scall(maximum possible amplitude) by bit depth
bit_depth = 16
full_scale = 2 ** bit_depth
# dbFS function
to_dbFS = lambda x: 20 * np.log10(x / full_scale)
# Read in the wave file
fname = "01.wav"
fs,data = wavfile.read(fname)
# Determine frame length(number of sample points in a frame) and total frame numbers by window length(how long is a frame in seconds)
window_length = 0.01
signal_length = data.shape[0]
frame_length = int(window_length * fs)
nframes = signal_length // frame_length
# Get frames by broadcast. No overlaps are used.
idx = frame_length * np.arange(nframes)[:,None] + np.arange(frame_length)
frames = data[idx].astype("int64") # Convert to in 64 to avoid integer overflow
# Get RMS and peaks
rms = ((frames**2).sum(axis=1)/frame_length)**.5
peaks = np.abs(frames).max(axis=1)
# Convert them to dbfs
dbfs_rms = to_dbFS(rms)
dbfs_peak = to_dbFS(peaks)
# Let's start to plot
# Get time arrays of every sample point and ever frame
frame_time = np.arange(nframes) * window_length
data_time = np.linspace(0,signal_length/fs,signal_length)
# Plot
f,ax = plt.subplots()
ax.plot(data_time,data,color="k",alpha=.3)
# Plot the dbfs values on a twin x Axes since the y limits are not comparable between data values and dbfs
tax = ax.twinx()
tax.plot(frame_time,dbfs_rms,label="RMS")
tax.plot(frame_time,dbfs_peak,label="Peak")
tax.legend()
f.tight_layout()
# Save serval details
f.savefig("whole.png",dpi=300)
ax.set_xlim(1,2)
f.savefig("1-2sec.png",dpi=300)
ax.set_xlim(1.295,1.325)
f.savefig("1.2-1.3sec.png",dpi=300)
整个时间跨度看起来像（右轴单位为dbFS):
And the voiced part looks like:
您可以看到 dbFS 值变大，而元音起点处的振幅变大：