模拟脉冲响应“声量子”算法.pdfVIP

模拟脉冲响应的“声量子”算法 马剑，李振格 天津大学建筑系，天津（300072 ） E-mail：mjtju@ 摘 要：本文通过实现模拟脉冲响应的“声量子”算法，解决了对房间各个位置的瞬时声压脉 冲响应的计算。它是在几何声学的基础上，利用波动声学的基本原理，根据声线追踪法的基 本原则，利用声波叠加原理计算该位置的声压，来模拟房间的瞬时声压脉冲响应，实现了对 任意复杂体形厅堂的脉冲响应模拟，并得到初步验证。 关键词：“声量子”追踪法，脉冲响应，厅堂音质 1. 概述 通常以几何声学为基础的声场模拟算法，仅仅适用于房间几何尺寸大于声波波长且房间 壁面粗糙度小于声波波长的情况，在模拟过程当中几何声波的干涉和衍射现象均被忽略，声 场的描述参量减少为能量、传播时间和传播方向，因此不能模拟计算房间各个位置的瞬时声 压脉冲响应，只能模拟某一位置的有效声压脉冲响应[1] 。鉴于此因，我们在声线追踪法的基 础之上，探索了一种模拟脉冲响应的新方法。 设想声源向四周无规则发射声粒子，其声压为： P P cos(ωt −kr) i A 其中，Pi——某时刻的瞬时声压； PA——声波振幅； ω——振动圆频率 ω=2πf ; k—— 圆波数k=ω/c ; r——传播距离。 根据声线追踪法的基本原理，计算到达房间指定位置的声粒子数，利用声波叠加原理计 算该位置的声压，来模拟房间的瞬时声压脉冲响应。此种模拟脉冲响应的方法不仅局限于高 频、大尺度房间，还可以扩充到低频、小尺度房间。此方法模拟的声压脉冲响应更加逼近现 实声场特性。 2. “声量子”追踪法的基本原理 [2] 这是在几何声学的基础上，利用波动声学的基本原理 ，根据声线追踪法的基本原则， 将声线看成具有波动性能的声粒子在封闭空间内的传播路径，计算到达指定位置的声粒子的 能量、相位和传播距离，再利用声波叠加原理计算该位置的声压，来模拟房间的瞬时声压脉 冲响应，实现了更加逼真的声场模拟。其具体特点如下所述： 1. 声源向四周无规则发射声粒子束，它们在与壁面碰撞以前都以平面声波的传播特性 沿直线向前运动。 2. 每个声粒子开始时都携带相同的能量，其大小取决于声源辐射的总能量和声粒子 数。 3. 声粒子遇到壁面作镜面反射后（其出射方向由反射定律决定），将沿新的方向继续 向前传播，遵循反射定律即入射角与反射角相等并且反射声波与入射声波在碰撞平面的 同侧。 -1- 4. 每次碰撞后能量减少为原来的（1-a）倍（a 为所碰壁面的吸声系数）。 5. 每次碰撞后相位角由它传播的距离及圆频率确定（φ= ωt + kr ）。 6. 事先设定的阈值，当其能量低于此值时，计算机就停止对该声粒子的跟踪，继而跟 踪下一声粒子。如此重复进行，直至所有声粒子均被跟踪完毕。 7. 在追踪过程当中，声粒子穿过预先设定的计算面积或计算体积时，计算机计 算并 记录下它到达的时间、能量、已传播的距离和新的传播方向。 8. 确定脉冲响应的样本宽度，计算在此时间间隔内的声粒子数，根据声波的叠加原理， 就可以得到指定位置的瞬时声压的时间分布图。 3. 用此算法求声场的脉冲响应 3.1 平面声波的叠加原理 平面声波的叠加原理是指声波在性状不甚规则的大房间内传播时，由于声波在壁面上无 数次反射的结果，对于房间内的任意一个位置，在某时刻的声压有可能是从各个方向传来的 反射声波的叠加，而且它们的相位都是随时间无规律变化的。为简单起见，先以两列具有相 同频率、但都有无规律变化位相的声波为例，讨论其合成声场的性质。 设某瞬时到达声场中某位置的两列声波分别为：