matlab计算混响时间t60代码

计算混响时间T60的Matlab代码如下：

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

% 将音频文件转换为单声道信号
y = mean(y, 2);

% 计算信号的短时能量
win_len = round(fs * 0.02); % 窗口长度为20ms
hop_len = round(fs * 0.01); % 帧移为10ms
E = st_energy(y, win_len, hop_len);

% 计算信号的自相关函数
r = xcorr(E);

% 找到自相关函数中第一个峰值的位置及其值
[~, locs] = findpeaks(r);
p = locs(1);

% 计算混响时间T60
T60 = -p / log(0.001);

其中，`st_energy`函数用于计算信号的短时能量，`xcorr`函数用于计算信号的自相关函数，`findpeaks`函数用于找到自相关函数中的峰值位置及其值。需要注意的是，该代码仅适用于单声道信号。