R语言中提取.wav文件dB声压级数值的相关函数有哪些？

核心原理 #

WAV文件本身存储的是音频采样的相对振幅值，无法直接输出dB SPL（声压级），需要先计算指定窗口内的振幅均方根（RMS），再基于参考值转换为dB单位的结果，你已安装的seewave包可直接完成该流程，无需额外安装工具。

具体实现代码 #

# 加载依赖包（seewave依赖tuneR，已安装可直接加载）
library(seewave)
library(tuneR)

# 1. 读取本地WAV文件
wav_obj <- readWave("你的音频文件路径.wav")
# 提取采样率
samp_rate <- wav_obj@samp.rate

# 2. 配置分析参数
window_ms <- 10 # 分析窗口长度，单位毫秒，可按需调整
overlap_pct <- 50 # 窗口重叠比例，单位%，可按需调整
window_length <- round(window_ms / 1000 * samp_rate)

# 3. 计算各窗口的RMS振幅
rms_series <- rms(
  wave = wav_obj,
  f = samp_rate,
  wl = window_length,
  ovlp = overlap_pct
)

# 4. 转换为dB值
# 仅需相对音量对比时ref设为1即可；需要真实物理SPL时ref替换为校准得到的20μPa对应振幅值
dB_series <- dB(rms_series, ref = 1)

# 5. 生成对应时间轴，合并为结果表
time_series <- seq(
  from = window_length/2 / samp_rate,
  to = length(wav_obj@left) / samp_rate,
  by = (window_length * (100 - overlap_pct)/100) / samp_rate
)
result_df <- data.frame(
  time_s = time_series,
  sound_pressure_level_dB = dB_series
)

输出的result_df即为每个时间点对应的dB值结果。

注意事项 #

• 若需真实物理意义的声压级，需提前用声级校准器录制已知dB值的校准音，计算对应RMS振幅后代入dB()函数的ref参数，否则结果为相对dB值，仅可用于同一段音频内的音量大小对比。
• 窗口长度越小时间精度越高，但dB值波动越大，可根据你的分析需求调整参数：需要毫秒级时间精度可设1-10ms窗口，需要平滑的整体音量变化趋势可设50-100ms窗口。
• 双声道音频默认处理左声道，如需处理右声道可将wav_obj@left替换为wav_obj@right，也可先对双声道振幅做平均后再计算。

内容的提问来源于stack exchange，提问作者stressedstudent