一种基于最大似然的语音信号混响时间(reverberation time)估计方法的MATLAB实现 文中所用音频文件 博客链接： https://blog.csdn.net/wlwdecs_dn/article/details/126328067

3．混响时间的优化 我们以中频为例对混响时间进行优化，方法是将部分吊顶材料由Plaster-Joists Suspended(五合板)改为Suspended Absorber(水泥膨胀珍珠岩板)。这两种材料不仅在做法上有所不同，它们的声学性质也有着很大的差异。 我们可以通过材质的吸声系数曲线对它们进行比较，在主工具条中点击“Material Properties”按钮启动 材质管理器对话框， 然后进入到 “Acoustic Data(声 学数据)”选项卡，这 里显示了材质在倍频 程中心频率上的吸声 系数及其曲线。通过 对比，我们可以发现 它们的差别主要集中 在中高频范围，五合 板对于中高频的的吸 声作用较弱，大部分 中高频声都被反射回 室内空间，而水泥膨 胀珍珠岩板的中高频 吸声作用相对于前者 来说要强一些。 五合板吸声系数曲线

该算法允许在 0.2 秒到 1.2 秒的范围内估计混响时间（RT 或 T60），并假设源和接收器不在临界距离内。 此功能不执行去噪，必须提前完成。 使用的算法在Heinrich W. Löllmann、Emre Yilmaz、Marco Jeub 和 Peter Vary： “一种改进的盲混响时间估计算法” 声学回声和噪声控制国际研讨会（IWAENC）， 以色列特拉维夫，2010 年 8 月。 （可参考www.ind.rwth-aachen.de/~bib/loellmann10a ） 笔记： 通过秒快速跟踪变化的 RT 的方法此版本中未实现直方图以进一步减少算法的复杂性。该程序的参数设置与用于模拟示例的参数设置不同参考文件。

用于计算室内混响时间，文件整合了计算理论依据，不同吸声结构的吸声系数表等内容，可方便设置吸声结构，调整混响时间。

rt_script.m 是主程序。 它生成一个文本文件和一个 pdf 文件来记录混响时间估计。 可以使用两种测量方法：1) 扬声器开扬声器关闭法，和 2) 气球爆破法。 自述文件中给出了测量两种方法的基本程序。 已发现扬声器开扬声器关闭方法比气球爆破方法准确得多。 气球爆破方法在 1000 Hz 以下通常会有 50% 以上的误差。 reverb_time.m 计算 1/3 倍频程时间记录的混响时间。 使用来自 makelNHANESNoisesm_ed.m 的随机时间测试信号的时间记录（也在 Matlab Central File Exchange 上）是使用扬声器开-扬声器关方法测量混响时间的理想选择。 使用气球爆破法，可以多次处理同一个文件，粗略估计每个 1/3 倍频程的混响时间。

适用于水下声学实验中，对于一个采集信号的音频信号，进行分析混响时间。