声压频域matlab实现三分之一倍频程

以下是在Matlab中实现声压频域三分之一倍频程的方法：

% 假设已经获取了声压信号的频域数据，存储在变量P中
% 假设采样率为Fs，频率范围为f

% 计算中心频率
f_c = [1, 1.25, 1.6, 2, 2.5, 3.15, 4, 5, 6.3, 8, 10]; % 中心频率数组
f_c_third = f_c(1:3:end); % 三分之一倍频程的中心频率数组

% 计算三分之一倍频程的上下限频率
f_l_third = f_c_third ./ sqrt(2); % 下限频率
f_u_third = f_c_third .* sqrt(2); % 上限频率

% 找到频率范围内所有符合条件的频率点
idx_third = find(f >= f_l_third(1) & f <= f_u_third(end));

% 计算三分之一倍频程的平均幅值
P_third = zeros(size(f_c_third));
for i = 1:length(f_c_third)
    idx = find(f >= f_l_third(i) & f <= f_u_third(i));
    P_third(i) = mean(P(idx));
end

% 绘制三分之一倍频程谱
semilogx(f_c_third, 20*log10(P_third), 'b.-');
xlabel('Frequency (Hz)');
ylabel('Sound Pressure Level (dB)');
title('1/3 Octave Band Spectrum');