matlab麦克风计算混响时间t60代码

以下是MATLAB计算混响时间（T60）的代码示例：

% 输入音频文件路径
audio_file = 'path/to/audio.wav';

% 读取音频文件和采样率
[x, Fs] = audioread(audio_file);

% 预设参数
N = 4096; % FFT窗口大小
n1 = 2000; % 分析起始点
n2 = n1 + N - 1; % 分析终止点

% 计算短时能量
E = sum(x(n1:n2).^2);

% 计算初始瞬时声压级
Lp0 = 20*log10(max(abs(x(n1:n2)))/0.00002);

% 重复计算短时能量和瞬时声压级，直到瞬时声压级下降至初始值的-60dB
Lp = Lp0;
while Lp-Lp0 > -60
    n1 = n1+N/2;
    n2 = n1+N-1;
    E = [E sum(x(n1:n2).^2)];
    Lp = [Lp 20*log10(max(abs(x(n1:n2)))/0.00002)];
end

% 计算信号能量下降到原始能量的0.001倍所需的时间（即混响时间T60）
E_ratio = E/E(1);
T60 = find(E_ratio <= 0.001, 1, 'first')/Fs;

% 打印T60值
fprintf('T60 = %.2f s\n', T60);

这个代码根据输入的音频文件计算混响时间（T60），并输出T60值。需要注意的是，这个代码只适用于单声道音频文件，如果需要处理多声道音频文件，需要对代码进行适当修改。