2018版matlab心理声学参数计算

在MATLAB中，可以使用PsySound3工具箱来计算心理声学参数。其中，常用的心理声学参数包括：

1. 声压级（Sound Pressure Level，SPL）
2. 频率加权声压级（A-weighted Sound Pressure Level，SPL(A)）
3. 声音的清晰度（Speech Clarity Index，SCI）
4. 声音的可懂度（Speech Intelligibility Index，SII）
5. 声音的响度（Loudness）
6. 声音的音调（Pitch）

以下是一个计算SPL和SPL(A)的示例代码：

% 读取音频文件
[y,Fs] = audioread('test.wav');

% 计算SPL
P_ref = 20e-6; % 参考声压级
SPL = 20*log10(rms(y)/P_ref);

% 计算SPL(A)
freq = (0:2047)/2048*Fs/2; % 频率向量
A_weight = db2mag(-39.4 + 12.5*log10(freq.^2./(freq.^2+12194^2)) + ...
    19.9*log10(freq.^2./(freq.^2+20.6^2))); % A权重
SPL_A = 10*log10(sum((y.*A_weight).^2)/sum(A_weight.^2)/P_ref^2);

其中，`test.wav`是要计算的音频文件，`P_ref`是参考声压级，`rms(y)`是音频信号的RMS值，`db2mag`是将分贝转换为线性幅度，`sum`是对向量求和。