应用matlab语言实现房间脉冲响应合成

《应用matlab语言实现房间脉冲响应合成》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、!"#$%&’($’))*$’(声频工程应用!"#$"%语言实现房间·论文·脉冲响应合成彭健新&，吴硕贤’(&)华南理工大学应用物理系，广州*&+,-+’)华南理工大学建筑系，广州*&+,-+.【摘要】探讨了一种由倍频带能量脉冲响应合成房间脉冲响应的方法，给出了基于!"#/$"%语言的实现方案，并将采用该方法合成房间脉冲响应计算得到的一些声学参数与采用仿真软件和实测得到的结果进行对比。【关键词】能量脉冲响应；房间脉冲响应；可听化【!"#$%&’$】#01213045467893017:;:9<=442:2>?@71=17>4971A8191=<8:2>?@71=1

2、7>4971:94B3CD1/AC95:75:7B?7715:930:7>C>1=)#01>=4<=C2:7>=4D:515E:30!"#$"%)"3@C73F7421CB4?73:B>C=C2131=74A3C:915A830:7213045FC9530C3A821C7?=12193E:307:2?@C3:49C=1B42>C=15)【()*+,%-#】G91=<8:2>?@71=17>4971；H442:2>?@71=17>4971；"?=C@:;C3:49&&+引言#%($.R#%(&.&R$"R!#%(&.!(&S$".(&.%R&室内声场仿真是建筑声学领

3、域的重要研究课题式中%为对应中心频率为Q&)*P;T&,UP;（根据实际之一。目前许多室内声场仿真软件已商品化，如瑞典需要）的&+个倍频带。哥德堡大学的I"##、丹麦技术大学的JKGJL、德由于能量脉冲响应不含有相位信息，实验结果国"K"声学设计公司的G"MG和G"HM等N&O。室内表明在可听化滤波中其相位并不重要NQO，因而可以声场仿真的技术核心是应用各种计算方法获取房间用产生和房间能量脉冲响应长度相等、采样间隔相脉冲响应，并在此基础上实现聆听其仿真声场的效同的白噪声信号果（可听化）。为了获取房间声学参数(大部分为单耳’($.R’(&.&R$"R’(&.!(&S

4、$".(’.指标.，需将各倍频程的房间能量脉冲响应经过转换并将上述白噪声信号复制&+份，用各倍频带能量脉合成房间脉冲响应。一般做法是先将各倍频程能量冲响应的均方根分别调制对应频带的白噪声信号，脉冲响应转换至频域，用立方曲线插值的方法，求出即频域幅值传递函数，用P:@A1=3转换关系构造出该幅(%($.R"#%($.’($.(Q.值传递函数的最小相位函数，再将各频域复数函数然后，采用无限冲击脉冲响应(VVH.数字滤波器去掉叠加得到房间脉冲响应在频域的表示，最后经傅里各倍频带多余的频率成分。为了提高处理速度，在满叶反变换得到时域的房间脉冲响应N’ONQO，其实现过程足

5、采样定理的条件下，对各倍频带信号进行重采样。比较复杂。笔者介绍一种简单、快速的房间脉冲响应例如，当"R&27，对中心频率为Q&)*P;，,QP;，’*+合成方法，它利用!"#$"%语言中信号处理工具箱P;等倍频带进行)倍下采样，的一些功能函数来实现。最后比较了由合成的房间(%()+.+R+，&，⋯，,脉冲响应计算得到的一些声学参数与采用仿真软件(*%(+.R#(-.++!+，&，⋯，,和实测得到的结果。,为!W)的最大整数。对中心频率为*++P;T&,’房间脉冲响应合成方法UP;等倍频带信号进行$倍上采样设能量脉冲响应长度为!，采样时间间隔为"(%(-W..-R+

6、，.，’.，⋯，!.(*%(-.R#（*）(一般为*27或&27.，则能量脉冲响应#($.可表示+-!+，.，’.，⋯，!.为：同时，为了能将结果保存为一般的音频文件（/7R--)&X《电声技术》!""#年第$期总第#%&期声频工程!"#$%&’($’))*$’(!"#），需将各倍频带的脉冲响应进行!"倍上采样，%运用最小阶数&产生模拟低通滤波器原型。使采样频率为$$%&!"#。&运用固有频率把模拟低通滤波器原型转换成#%"($+!")$*,，!"，-!".⋯，&"!"模拟带通滤波器。#"’’($)*!(/),$!,，!"，-!".⋯，&"!"’运用冲击响应不变法

7、或双线性不变法把模拟&"为各倍频带房间脉冲响应的长度。最后，将得到滤波器转换成数字滤波器。的各倍频带房间脉冲响应相加得：采用637839语言实现房间脉冲响应合成的流&,程如图&所示。下面列出其中有关数字信号处理部#(’)’*$(’*#"($))"#"’’($)(0)分的一个具体的637839程序。"*&#(’)为合成的房间脉冲响应，式中(’*&+*1。最后，将房间脉冲响应归一化，保存为234文件。5基于637839语言的软件实现由于637839语言信号处理工具箱提供各种数字信号处理功能函数，可以很方便地直接调用，或者通过4:;;调用。下面对有关的主要函数作简单介绍

8、。(&)<