发明名称 一种房间混响合成方法
摘要 本发明提供一种房间混响合成方法,通过利用房间混响镜像合成法合成早期反射,并以合成的早期反射为输入,在每个临界频带内利用参数方法计算出早期反射随时间衰减的曲线,并将每个频带内的曲线合并为一个所有频带范围内的曲面作为早期反射的EDR,利用该早期反射的EDR通过频域逼近方法求出逼近滤波器的系数作为FDN参数,通过反馈延时网络法对该FDN参数进行后期混响合成,进而合成房间混响,使本发明房间混响合成的方法具有了基于物理特性合成方法能够准确合成早期反射和基于感知合成方法的计算速度快的优点。
申请公布号 CN104703111B 申请公布日期 2016.09.28
申请号 CN201310661781.4 申请日期 2013.12.09
申请人 中国科学院声学研究所 发明人 李军锋;夏日升;颜永红
分类号 H04S3/00(2006.01)I;G10K15/12(2006.01)I 主分类号 H04S3/00(2006.01)I
代理机构 北京法思腾知识产权代理有限公司 11318 代理人 杨小蓉;杨青
主权项 一种房间混响合成方法,其特征在于,所述的方法包括:步骤1)设定房间声学环境的物理参数;步骤2)根据步骤1)中获得的物理参数计算空气对声音的吸收值;步骤3)利用步骤2)中获得的空气对声音的吸收值,通过下述计算公式求出声音在空气中传输时的声压衰减因子:p=exp(‑a)其中,p表示声压衰减因子,a表示空气对声音的吸收值;步骤4)根据房间混响镜像合成法计算每个声像到接收位置的距离,然后依据步骤3)获得的声压衰减因子求出空气对声音吸收后的声压,该声压作为合成的早期反射,该声压的计算公式表示为:p(d)=p<sub>0</sub>exp(‑ad)其中,p<sub>0</sub>为原声源在空气中经过传输距离d之后的声压;步骤5)以步骤4)中合成的早期反射为输入,在每个临界频带内利用参数方法计算出早期反射随时间衰减的衰减曲线,将每个频带内的衰减曲线合并为一个所有频带范围内的曲面作为早期反射的EDR,进而利用该早期反射的EDR通过频域逼近方法求出逼近滤波器的系数作为FDN参数;步骤6)利用反馈延时网络法对步骤5)中获得的FDN参数进行后期混响合成;步骤7)将测得的直达声、步骤4)中获得早期反射和步骤6)中获得的后期混响进行相加从而合成房间混响;所述步骤5)中的衰减曲线的计算公式表示为:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mover><mi>h</mi><mo>^</mo></mover><mi>p</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><mi>exp</mi><mrow><mo>(</mo><mo>-</mo><mi>a</mi><mi>d</mi><mo>)</mo></mrow></mrow><mrow><mn>8</mn><mi>d</mi><mi>r</mi></mrow></mfrac><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><mfrac><msub><mi>B</mi><mi>z</mi></msub><mrow><mi>log</mi><mrow><mo>(</mo><mfrac><msub><mi>B</mi><mi>y</mi></msub><msub><mi>B</mi><mi>x</mi></msub></mfrac><mo>)</mo></mrow></mrow></mfrac><mo>{</mo><mi>E</mi><mi>i</mi><mo>&lsqb;</mo><mi>log</mi><mrow><mo>(</mo><mfrac><msub><mi>B</mi><mi>z</mi></msub><msub><mi>B</mi><mi>x</mi></msub></mfrac><mo>)</mo></mrow><mo>&rsqb;</mo><mo>+</mo><mi>log</mi><mo>&lsqb;</mo><mi>log</mi><mrow><mo>(</mo><mfrac><msub><mi>B</mi><mi>z</mi></msub><msub><mi>B</mi><mi>x</mi></msub></mfrac><mo>)</mo></mrow><mo>&rsqb;</mo><mo>-</mo><mi>E</mi><mi>i</mi><mo>&lsqb;</mo><mi>log</mi><mrow><mo>(</mo><mfrac><msub><mi>B</mi><mi>z</mi></msub><msub><mi>B</mi><mi>y</mi></msub></mfrac><mo>)</mo></mrow><mo>&rsqb;</mo><mo>-</mo><mi>log</mi><mo>&lsqb;</mo><mi>log</mi><mrow><mo>(</mo><mfrac><msub><mi>B</mi><mi>z</mi></msub><msub><mi>B</mi><mi>y</mi></msub></mfrac><mo>)</mo></mrow><mo>&rsqb;</mo><mo>}</mo><mo>,</mo></mrow></mtd><mtd><mrow><mi>i</mi><mi>f</mi><mi> </mi><msub><mi>B</mi><mi>x</mi></msub><mo>&NotEqual;</mo><msub><mi>B</mi><mi>y</mi></msub><mo>&NotEqual;</mo><msub><mi>B</mi><mi>z</mi></msub></mrow></mtd></mtr><mtr><mtd><mrow><mfrac><msub><mi>B</mi><mi>z</mi></msub><mrow><mi>log</mi><mrow><mo>(</mo><mfrac><msub><mi>B</mi><mi>z</mi></msub><mi>B</mi></mfrac><mo>)</mo></mrow></mrow></mfrac><mo>{</mo><mi>E</mi><mi>i</mi><mo>&lsqb;</mo><mi>log</mi><mrow><mo>(</mo><mfrac><msub><mi>B</mi><mi>z</mi></msub><mi>B</mi></mfrac><mo>)</mo></mrow><mo>&rsqb;</mo><mo>+</mo><mi>log</mi><mo>&lsqb;</mo><mi>log</mi><mrow><mo>(</mo><mfrac><msub><mi>B</mi><mi>z</mi></msub><mi>B</mi></mfrac><mo>)</mo></mrow><mo>&rsqb;</mo><mo>+</mo><mi>r</mi><mo>}</mo><mo>,</mo></mrow></mtd><mtd><mrow><mi>i</mi><mi>f</mi><mi> </mi><msub><mi>B</mi><mi>z</mi></msub><mo>=</mo><msub><mi>B</mi><mi>y</mi></msub><mo>&NotEqual;</mo><msub><mi>B</mi><mi>x</mi></msub><mover><mo>=</mo><mi>&Delta;</mi></mover><mi>B</mi><mi> </mi><mi>o</mi><mi>r</mi><mi> </mi><msub><mi>B</mi><mi>z</mi></msub><mo>=</mo><msub><mi>B</mi><mi>x</mi></msub><mo>&NotEqual;</mo><msub><mi>B</mi><mi>y</mi></msub><mover><mo>=</mo><mi>&Delta;</mi></mover><mi>B</mi></mrow></mtd></mtr><mtr><mtd><mfrac><mrow><mi>B</mi><mo>-</mo><msub><mi>B</mi><mi>z</mi></msub></mrow><mrow><mi>log</mi><mrow><mo>(</mo><mfrac><mi>B</mi><msub><mi>B</mi><mi>z</mi></msub></mfrac><mo>)</mo></mrow></mrow></mfrac></mtd><mtd><mrow><mi>i</mi><mi>f</mi><mi> </mi><msub><mi>B</mi><mi>z</mi></msub><mo>&NotEqual;</mo><msub><mi>B</mi><mi>x</mi></msub><mo>=</mo><msub><mi>B</mi><mi>y</mi></msub><mover><mo>=</mo><mi>&Delta;</mi></mover><mi>B</mi></mrow></mtd></mtr><mtr><mtd><mrow><mi>B</mi><mo>,</mo></mrow></mtd><mtd><mrow><mi>i</mi><mi>f</mi><mi> </mi><msub><mi>B</mi><mi>x</mi></msub><mo>=</mo><msub><mi>B</mi><mi>y</mi></msub><mo>=</mo><msub><mi>B</mi><mi>z</mi></msub><mover><mo>=</mo><mi>&Delta;</mi></mover><mi>B</mi></mrow></mtd></mtr></mtable></mfenced></mrow>]]></math><img file="FDA0001017365260000011.GIF" wi="1710" he="599" /></maths>其中,r=0.57721,Ei表示一阶指数积分,d表示声音从声源到观测点的距离,t表示时间;且满足<img file="FDA0001017365260000021.GIF" wi="1195" he="110" />β为墙面的吸声系数,L<sub>x</sub>、L<sub>y</sub>和L<sub>z</sub>分别表示房间的每一边长度,B<sub>x</sub>、B<sub>y</sub>和B<sub>z</sub>表示不同墙面吸声系数的函数,<img file="FDA0001017365260000022.GIF" wi="35" he="47" />表示“定义为”。
地址 100190 北京市海淀区北四环西路21号