注意:此功能现在可从IOSR MATLAB Toolbox作为IOSR.Coustics.Irstats获得。
------------------------
计算RT,DRR,CTE和EDT用于脉冲响应文件
RT = IR_STATS(文件名)返回混响时间(至-60 dB)
使用基于ISO 3382-1:2009的方法。该功能使用反向
累积梯形整合以估计衰减曲线,以及一个
线性最小二乘拟合估计0 dB和-60 dB之间的斜率。
估计在八度乐队中采取,整体数字是一个
500 Hz和1 kHz带的平均值。
文件名应该是音频文件的完整路径或名称
Matlab搜索路径上的音频文件。该文件可以是任何格式
万博1manbetx由Audioread功能支持,并有任何数量的渠道;
估计(和图标)将为每个频道返回。
该函数返回1xn向量的RTS,其中n是数量
音频文件中的通道。
该功能确定直接声音作为平方的峰值
脉冲响应。
[RT,DRR] = IR_STAT(文件名)返回直响混合率
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(文件名)返回早期索引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(文件名)返回Octave频段中心
频率CFS用于计算RT。
[RT,DRR,CTE,CFS,EDT] = IR_STATS(文件名)返回早期衰减
时间EDT,大小与RT相同。衰减曲线的斜率
从0到10 dB之间的拟合确定。衰减时间是
根据60 dB衰减所需的时间从斜率计算。
... = ir_stats(...,'参数',值)允许众多
要指定的参数。这些参数是:
'图':{false} |真的
控制是否绘制衰减曲线。具体而言,图形
被脉冲响应,衰减曲线和线性绘制
最小二乘适用于音频的每个八度频段和通道
文件。如果指定了EDT输出,则EDT拟合也将是
绘制。
'te':{0.05} |标量子
指定早期限制(以秒为单位)。
'spec':{'均值'} |'满的'
确定RT和EDT输出的性质。规格='卑鄙'
(默认)报告的RT和EDT是500 Hz的平均值
和1 kHz的乐队。使用spec ='full',函数返回
RT和EDT按返回的每个八度频带计算
CFS;RT和EDT具有大小[M n],其中m =长度(CFS)。
'y_fit':{[0 60]} |双元素矢量
指定衰减曲线应该的分贝范围
评估。例如,'y_fit'可以是[-5 -25]或[-5 -35]
对应于RT20和RT30。
'校正':{0.0025} |标量子
指定上面给出的校正参数C(以秒为单位)
用于DRR和CTE计算。最多10毫秒的值
建议在文献中。
八度带滤波器根据ANSI S1.1-1986和IEC计算 作者要感谢Feifei Xiong的意见 参考 [1] Zahorik,P.,2002:'直接反向的能量比 另见Audioread,Octdsgn。
标准。请注意,OCTDSGN功能推荐中心频率
FS / 200
校正参数。
敏感度'美国声学学会,
112,2110-2117。
Christopher Hummersone(2021)。脉冲响应声学信息计算器(https://github.com/iosr-surrey/matlabtoolbox),github。检索到。