频域声压级 matlab

频域声压级是指在一定频率范围内的声压级。在Matlab中，可以通过以下步骤计算频域声压级：

1. 读取音频文件并进行FFT变换，得到频域数据。

[y, Fs] = audioread('audio.wav'); % 读取音频文件
N = length(y); % 信号长度
Y = fft(y); % 进行FFT变换
f = (0:N-1)*(Fs/N); % 计算频率

2. 计算声压级。

p = abs(Y).^2/N/Fs; % 计算功率谱密度
p(2:end-1) = 2*p(2:end-1); % 修正功率谱密度
SPL = 10*log10(sum(p)/2e-5^2); % 计算声压级

3. 对于多个文件的情况，可以使用循环读取并计算每个文件的声压级。

folder = 'path/to/folder'; % 文件夹路径
files = dir(fullfile(folder, '*.wav')); % 获取所有wav文件
for i = 1:length(files)
    filename = fullfile(folder, files(i).name); % 获取文件名
    [y, Fs] = audioread(filename); % 读取音频文件
    N = length(y); % 信号长度
    Y = fft(y); % 进行FFT变换
    f = (0:N-1)*(Fs/N); % 计算频率
    p = abs(Y).^2/N/Fs; % 计算功率谱密度
    p(2:end-1) = 2*p(2:end-1); % 修正功率谱密度
    SPL = 10*log10(sum(p)/2e-5^2); % 计算声压级
    disp(['文件', num2str(i), '的声压级为', num2str(SPL), 'dB']); % 显示结果
end