注意:这个函数现在可以从IoSR的Matlab工具箱中获得:IoSR .acoustics. irstats。
-------------------------
计算脉冲响应文件的RT, DRR, Cte和EDT
返回混响时间(- 60db)
使用基于ISO 3382-1:2009的方法。函数使用reverse
累积梯形积分估计衰减曲线，以及
线性最小二乘拟合估计斜率在0 dB和-60 dB之间。
估计是在八度波段和总的数字是一个
500hz和1khz频带的平均值。
FILENAME应该是音频文件的完整路径或
音频文件的Matlab搜索路径。该文件可以是任何格式
万博1manbetx由AUDIOREAD功能支持，并且有任意数量的通道;
每个频道将返回估计值(和图)。
函数返回一个1xN的矢量RTs，其中N是
音频文件中的通道。
函数决定了直接声音的平方的峰值
脉冲响应。
[RT,DRR] = IR_STATS(FILENAME)返回直接混响比
DRR表示脉冲;DRR和rt大小相同，这是计算出来的
以下列方式:
DRR = 10 * log10 (X (T0-C: T0 + C) ^ 2 / X (T0 + C + 1:结束)^ 2)
其中X是脉冲的近似积分，T0是
C=2.5ms[1]。
[RT,DRR,CTE] = IR_STATS(FILENAME)返回早到晚的索引CTE
的冲动,CTE和rt的尺寸是一样的
以下方式:
CTE = 10 * log10 (X (T0-C: T0 + TE) ^ 2 / X (T0 + TE + 1:结束)^ 2)
其中TE为50毫秒。
[RT,DRR,CTE,CFS] = IR_STATS(FILENAME)返回倍频中心
频率CFS用于计算RT。
[RT,DRR,CTE,CFS,EDT] = IR_STATS(FILENAME)返回早期衰减
时间EDT，其大小与rt相同
是由0到-10 dB之间的配合确定的。衰变时间为
根据斜率计算60db衰减所需的时间。
．．．= IR_STATS(…，'PARAMETER'，VALUE)允许无数
要指定的参数。这些参数是:
'graph': {false} | true
控制是否绘制衰减曲线。具体来说,图
脉冲响应，衰减曲线和线性
最小二乘法适用于音频的每个八度频带和信道
文件。如果指定了EDT输出，则EDT匹配也将是
绘制。
'te':{0.05} |标量
指定早期时间限制(以秒为单位)。
'spec': {'mean'} | 'full'
确定RT和EDT输出的性质。规范=“的意思是”
(默认)报告的RT和EDT是500hz的平均值
和1千赫波段。对于spec='full'，函数返回
为返回的每个八度频带计算的RT和EDT
慢性疲劳综合症;RT和EDT有大小[M N]，其中M=长度(CFS)。
'y_fit':{[0 60]} |二元向量
指定衰减曲线应超过的分贝范围
被评估。例如，'y_fit'可能是[-5 -25]或[-5 -35]
分别对应RT20和RT30。
'correction':{0.0025} |标量
指定上面给出的校正参数C(秒)
用于DRR和CTE计算。值可达10毫秒
在文献中提到。
倍频带滤波器根据ANSI S1.1-1986和IEC计算
标准。注意，OCTDSGN函数推荐中心频率
取值范围为fs/200 < Fc < fs/5。
作者要感谢熊飞飞的意见
校正参数。
参考文献
Zahorik, P.， 2002:直接-混响能量比
《美国声学学会杂志》
112年,2110 - 2117。
请参见AUDIOREAD, OCTDSGN。
引用作为
克里斯托弗Hummersone(2021)。脉冲响应声学信息计算器GitHub (https://github.com/IoSR-Surrey/MatlabToolbox)。检索．
