CN103916810A - 一种时域声能量对比控制方法和*** - Google Patents
一种时域声能量对比控制方法和*** Download PDFInfo
- Publication number
- CN103916810A CN103916810A CN201310001166.0A CN201310001166A CN103916810A CN 103916810 A CN103916810 A CN 103916810A CN 201310001166 A CN201310001166 A CN 201310001166A CN 103916810 A CN103916810 A CN 103916810A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- centerdot
- msup
- math
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 71
- 230000004044 response Effects 0.000 claims abstract description 100
- 230000014509 gene expression Effects 0.000 claims abstract description 67
- 238000005070 sampling Methods 0.000 claims description 10
- 238000012546 transfer Methods 0.000 claims description 9
- 239000011159 matrix material Substances 0.000 claims description 6
- 238000005457 optimization Methods 0.000 claims description 6
- 238000004088 simulation Methods 0.000 claims description 6
- 238000012360 testing method Methods 0.000 claims description 6
- 238000004364 calculation method Methods 0.000 claims description 3
- 238000013461 design Methods 0.000 description 5
- 238000010586 diagram Methods 0.000 description 4
- 230000000694 effects Effects 0.000 description 4
- 230000000052 comparative effect Effects 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 230000009467 reduction Effects 0.000 description 2
- 238000005316 response function Methods 0.000 description 2
- 238000001228 spectrum Methods 0.000 description 2
- 235000018185 Betula X alpestris Nutrition 0.000 description 1
- 235000018212 Betula X uliginosa Nutrition 0.000 description 1
- 238000003491 array Methods 0.000 description 1
- 230000015556 catabolic process Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000006731 degradation reaction Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 239000003292 glue Substances 0.000 description 1
- 230000010355 oscillation Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
Landscapes
- Circuit For Audible Band Transducer (AREA)
Abstract
本发明涉及一种时域声能量对比控制方法,包括:布放扬声器阵列,设定明区和暗区的控制点;其中,所述的明区是指设定为具有高声能量的区域,所述的暗区是指设定为具有低声能量的区域;选定参考频率,列出明区中的频率响应一致性约束表达式;分别列出明区和暗区的平均声能量表达式;给合频率响应一致性约束表达式和平均声能量表达式,按照频率响应一致性约束的时域声能量对比控制准则,计算出各个通道的时域冲激响应滤波器信号。
Description
技术领域
本发明涉及声学领域,特别涉及一种时域声能量对比控制方法和***。
背景技术
声能量对比控制方法被广泛应用于电视机、电脑和手机等个人音频***领域。声能量对比控制方法能够在不打扰其它用户(暗区)的情况下将大部分声能量聚集在用户周围(明区),从而为用户打造个性化聆听空间,提高生活品质。
现有技术中的声能量对比控制方法主要分为两大类:频域设计和时域设计。无论是频域声能量对比控制方法还是时域声能量对比控制方法都有各自的缺陷。
(1)、现有技术中的频域声能量对比控制方法通常只计算一系列离散控制频率点上的最优系数,然后通过直接对频域系数进行逆傅里叶变换得到时域冲激响应滤波器信号,因此时域冲激响应滤波器信号的因果性无法得到保证。这使得现有技术中的频域声能量对比控制方法在非选定频率点上的对比聚焦效果下降,并且在不同的频率点上有不同的响应从而无法在明区得到较好的音质。
(2)、现有技术中的时域声能量对比控制方法直接在时域中设计,避免了时域冲激响应滤波器信号的非因果性问题,因此可以部分解决频域声能量对比控制方法在非控制频率点对比聚焦效果下降的问题。但是现有技术中的时域声能量对比控制方法仍然没有考虑频率响应一致性的问题,通常会导致在明区的音质下降。
现有技术中的声能量对比控制方法所存在的这些问题造成了在个人音频***中并不适合处理宽带信号输入,降低了用户区域的音质,在暗区中泄露更多的声能量,无法为用户营造个性化的私人聆听空间。针对现有声能量对比控制方法中存在的音质和宽带对比聚焦能力下降的问题,需要寻找更为简便有效的方法来改善其音质和宽带对比聚焦能力。
发明内容
本发明的目的在于克服现有技术中的声能量对比控制方法所存在的音质和宽带对比聚焦能力下降的问题,从而提供一种能够改善其音质和宽带对比聚焦能力的时域声能量对比控制方法。
为了实现上述目的,本发明提供了一种时域声能量对比控制方法,包括:
步骤1)、布放扬声器阵列,设定明区和暗区的控制点;其中,所述的明区是指设定为具有高声能量的区域,所述的暗区是指设定为具有低声能量的区域;
步骤2)、选定参考频率,列出明区中的频率响应一致性约束表达式;
步骤3)、分别列出明区和暗区的平均声能量表达式;
步骤4)、给合步骤2)的频率响应一致性约束表达式和步骤3)中的平均声能量表达式,按照频率响应一致性约束的时域声能量对比控制准则,计算出各个通道的时域冲激响应滤波器信号。
上述技术方案中,在所述的步骤1)中,所布放的扬声器阵列为线性阵列或圆形阵列或随机阵列。
上述技术方案中,在所述的步骤1)中,所述明区或暗区的区域形状是方形或圆形或线型。
上述技术方案中,所述步骤2)包括:
步骤2-1)、获得扬声器第l通道到明区的第k个控制点在频率f处的传递函数gBlk(f),其中,1≤l≤L,1≤k≤K,K和L分别为明区的控制点个数和扬声器通道个数;
步骤2-2)、根据步骤2-1)所得到的传递函数gBlk(f),列出扬声器阵列到明区第k个控制点的频率响应pBk(f),其表达式如下:
其中Ts是采样时间间隔,也是采样频率fs的倒数,wl(m)是对应第l通道扬声器的时域冲激响应滤波器的待求的第m个系数,M是时域冲激响应滤波器的长度;
步骤2-3)、将步骤2-2)所得到的频率响应pBk(f)写成向量形式:
pBk(f)=wTsBk(f)
其中w为待求的ML×1向量,定义为:
w=[w1(0),…,w1(M-1),…,wL(0),…,wL(M-1)]T
sBk(f)为ML×1向量,定义为:
步骤2-4)、选定参考频率fr,根据之前得到的频率响应,定义明区的频率响应一致性约束RV,该频率响应一致性约束的表达式为:
其中是取元素的实数部分,Ω是所有约束频率点的集合 是取不超过元素的最大整数,fc是集合Ω中的一个约束频率点,BΩ是约束频率点的个数;Q的表达式为:
上述技术方案中,在所述的步骤2-1)中,采用诸如B&K PULSE的音频测试仪器或者通过建模仿真获得所述扬声器第l通道到明区的第k个控制点在频率f处的传递函数gBlk(f)。
上述技术方案中,所述步骤3)包括:
步骤3-1)、将扬声器第l通道到明区的第k个控制点之间的冲激响应建模成一个长度为I的FIR滤波器,并采用诸如B&K PULSE的音频测试仪器或者通过建模仿真获得其系数,定义该FIR滤波器第i个系数为hBlk(i);
步骤3-2)、假设输入给扬声器***的信号为x(n),则在明区第k个控制点的输出信号为
其中,rBk是ML×1向量,可定义如下:
rBk(n)=[rBlk(n),...rBlk(n-M+1),…,rBLk(n),…,rBLk(n-M+1)]T
步骤3-3)、假设输入信号x(n)为稳态信号,则在步骤3-2)中输出信号yBk(n)表达式的基础上,定义明区的平均声能量为:
其中E{·}代表取元素的期望值,RB(n)表达式如下:
步骤3-4)、将扬声器第l通道到暗区的第kD个控制点之间的响应建模成一个长度为I的FIR滤波器,则其第i个系数定义为假设输入给扬声器***的信号为x(n),则在暗区第kD个控制点的输出信号为
其中,1≤kD≤KD,KD为暗区控制点个数,是ML×1向量:
定义暗区的平均声能量为:
其中
上述技术方案中,所述的步骤4)包括:
步骤4-1)、结合步骤2)给出的频率响应一致性约束的表达式和步骤3)给出平均声能量表达式,按照频率响应一致性约束的时域声能量对比控制准则,列出其优化问题为:
步骤4-2)、对步骤4-1)所得到的优化问题求出其解:
其中Pmax{·}是求解矩阵对应最大特征值的单位特征向量,U为单位矩阵,δ为鲁棒性参数,β为权重参数;参数δ和β都取正数;
步骤4-3)、将步骤4-2)所得到的向量w按每隔M个元素进行切割,得到各个通道的时域冲激响应滤波器信号。
本发明还提供了一种时域声能量对比控制***,包括:
扬声器阵列布放模块,用于布放扬声器阵列,并设定明区和暗区的控制点;其中,所述的明区是指设定为具有高声能量的区域,所述的暗区是指设定为具有低声能量的区域;
频率响应一致性约束表达式获取模块,用于选定参考频率,列出明区中的频率响应一致性约束表达式;
平均声能量表达式获取模块,用于分别列出明区和暗区的平均声能量表达式;
时域冲激响应滤波器信号计算模块,用于给合所述频率响应一致性约束表达式获取模块所得到的频率响应一致性约束表达式和平均声能量表达式获取模块所得到的平均声能量表达式,按照频率响应一致性约束的时域声能量对比控制准则,计算出各个通道的时域冲激响应滤波器信号。
本发明的优点在于:
1、本发明直接在时域中设计,避免了常规频域声能量对比控制设计方法中由逆傅里叶变换得到的时域冲激响应滤波器信号的非因果性。
2、相比于常规声能量对比控制设计方法,在相同配置情况下,本发明的宽带信号对比聚焦能力会远大于常规声能量对比控制方法的宽带信号比聚焦能力。
3、本发明加入频率响应一致性约束,使得明区的频率响应平坦度加强,从而改善了明区的音质。
4、本发明的声能量对比控制方法能够广泛应用于电视机,电脑和手机等个人音频领域,在宽带信号输入时,能够提高明区的音质,并取得较好的聚焦能力,从而营造出较好的私人聆听空间。
附图说明
图1是本发明的基于频率响应一致性约束的时域声能量对比控制方法的流程图;
图2是在一个实施例中,线性扬声器阵列的布放示意图;
图3(a)为扬声器第三通道到明区中心控制点的频率响应的示意图;
图3(b)为扬声器第三通道到明区中心控制点的冲激响应的示意图;
图4(a)为本发明方法与现有频域方法在宽带信号上的对比聚焦性能的比较示意图;
图4(b)为本发明方法与现有时域方法在宽带信号上的对比聚焦性能的比较示意图;
图5为本方法和现有技术中的声能量对比控制方法在明区中心控制点处的频率响应的示意图。
具体实施方式
现结合附图对本发明作进一步的描述。
本发明的基本思想是直接在时域中进行设计,并利用频率一致性约束,改善明区控制点的频率响应平坦度,从而使声能量对比控制方法在明区能够改善音质,同时提高了宽带信号的对比聚焦能力。基于上述思想设计的本发明的方法消除了现有技术中的声能量控制方法存在的音质下降和宽带信号聚焦能力的下降问题。
参考图1,本发明的基于频率响应一致性约束的时域声能量对比控制方法包括以下步骤:
步骤1)、布放扬声器阵列,设定明区和暗区的控制点;其中,所述的明区是指设定为具有高声能量的区域,所述的暗区是指设定为具有低声能量的区域;
步骤2)、选定参考频率,列出明区中的频率响应一致性约束表达式;
步骤3)、分别列出明区和暗区的平均声能量表达式;
步骤4)、给合步骤2)的频率响应一致性约束表达式和步骤3)中的平均声能量表达式,按照频率响应一致性约束的时域声能量对比控制准则,计算出各个通道的时域冲激响应滤波器信号。
下面对本发明方法中的各个步骤做进一步的描述。
在所述的步骤1)中,所布放的扬声器阵列为线性阵列或圆形阵列,也可以为随机阵列。所述明区或暗区的区域形状可以是方形或圆形,也可以是线型。
所述的步骤2)具体包括以下步骤:
步骤2-1)、采用B&K PULSE等音频测试仪器或者通过建模仿真获得扬声器第l(1≤l≤L)通道到明区的第k(1≤k≤K)个控制点在频率f处的传递函数gBlk(f),其中,K和L分别为明区的控制点个数和扬声器通道个数;
步骤2-2)、根据步骤2-1)所得到的传递函数gBlk(f),列出扬声器阵列到明区第k个控制点的频率响应pBk(f),其表达式如下:
其中Ts是采样时间间隔,也是采样频率fs的倒数,wl(m)是对应第l通道扬声器的时域冲激响应滤波器的待求的第m个系数,M是时域冲激响应滤波器的长度;
步骤2-3)、将步骤2-2)所得到的频率响应pBk(f)写成向量形式:
pBk(f)=wTsBk(f)
其中w为待求的ML×1向量,定义为:
w=[w1(0),…,w1(M-1),…,wL(0),…,wL(M-1)]T
sBk(f)为ML×1向量,定义为:
步骤2-4)、选定参考频率fr,根据之前得到的频率响应,定义明区的频率响应一致性约束RV,该频率响应一致性约束的表达式为:
其中是取元素的实数部分,Ω是所有约束频率点的集合 是取不超过元素的最大整数,fc是集合Ω中的一个约束频率点,BΩ是约束频率点的个数。Q的表达式为:
所述的步骤3)具体包括以下步骤:
步骤3-1)、将扬声器第l通道到明区的第k个控制点之间的冲激响应建模成一个长度为I的FIR滤波器,并采用B&K PULSE等音频测试仪器或者通过建模仿真获得其系数,定义该FIR滤波器第i个系数为hBlk(i)。
步骤3-2)、假设输入给扬声器***的信号为x(n),则在明区第k个控制点的输出信号为
其中,rBk是ML×1向量,可定义如下:
rBk(n)=[rBlk(n),…rBlk(n-M+1),…,rBLk(n),…,rBLk(n-M+1)]T
步骤3-3)、假设输入信号x(n)为稳态信号,则在步骤3-2)中输出信号yBk(n)表达式的基础上,明区的平均声能量定义为:
其中E{·}代表取元素的期望值,RB(n)表达式如下:
步骤3-4)、与明区类似,将扬声器第l通道到暗区的第kD个控制点之间的响应同样建模成一个长度为I的FIR滤波器,则其第i个系数定义为假设输入给扬声器***的信号为x(n),则在暗区第kD(1≤kD≤KD,KD为暗区控制点个数)个控制点的输出信号为
其中,是ML×1向量,与rBk定义类似:
同理,暗区的平均声能量可以定义为
RD(n)的定义与RB(n)类似。
所述的步骤4)具体包括以下步骤:
步骤4-1)、结合步骤2)给出的频率响应一致性约束的表达式和步骤3)给出平均声能量表达式,按照频率响应一致性约束的时域声能量对比控制准则,列出其优化问题为:
步骤4-2)、对步骤4-1)所得到的优化问题求出其解:
其中Pmax{·}是求解矩阵对应最大特征值的单位特征向量,U为单位矩阵,δ为鲁棒性参数,β为权重参数;参数δ和β都取正数。
步骤4-3)、将步骤4-2)所得到的向量w按每隔M个元素进行切割,得到各个通道的时域冲激响应滤波器信号。
为了更好的理解本发明,下面结合附图和具体实施例对本发明的方法作进一步详细描述。
在一个实施例中,如图2所示,在中国科学院声学研究所全消声室内,放置一线性扬声器阵列,明区和暗区处在扬声器阵列中垂线上左右各45度方向,距离扬声器阵列间距离都为1m,并与扬声器阵列处于同一水平面上;其中扬声器阵列由8个动圈式惠威扬声器组成,间距为12cm,扬声器单元最大功率为15W,额定阻抗8Ω;明区和暗区的控制点都由5元传声器阵列定义,间隔为8cm,传声器为丹麦B&K公司的4189传声器;实验选用的硬件装置为功放、PULSE测量仪和PC机。
本实施例的具体实施过程包括以下几步:
(1)***采样率fs设为8kHz,时域冲激响应滤波器信号的长度M设为100。首先测量扬声器阵列到明区控制点和暗区控制点的传递函数,这可以通过PULSE仪器直接测量扬声器阵列到各个控制点的频响应函数得到,频率间隔为10Hz。扬声器阵列到明区和暗区控制点的冲激响应可以直接通过对频响应函数做逆傅里叶变换得到,其长度为800。图3(a)为扬声器第三通道到明区中心控制点的频率响应的示意图,图3(b)为扬声器第三通道到明区中心控制点的冲激响应的示意图。从图中可以看出,在考虑的频带频响应比较平坦。这里***采样率fs设为8kHz,只考虑人声频段,因此***里接了一个固定的低通滤波器模块,截止频率设在了3.5kHz,这也保证了直接从频响应获得冲激响应的因果性。
(2)选定参考频率fr为1kHz,约束频率点fc为[80,80×2,…80×49]Hz,列出频率响应一致性的约束表达式。
(3)假定扬声器阵列***的输入信号x(n)为白噪声,即满足:
基于上述白噪声假设,列出在明区和暗区的声能量表达式。
(4)按照频率响应一致性约束的时域声能量对比控制,计算出权向量w,其中δ取0.05,β取0.05。
(5)将向量w按每隔M个元素进行切割,得到各个通道的时域冲激响应滤波器信号。
图4给出了本发明方法在宽带信号上的对比聚焦性能,并与现有技术中的方法做了比较(为了看的更清楚,图中分开进行比较,图4(a)为本发明方法与现有频域方法的比较,图4(b)为本发明方法与现有时域方法的比较)。从图中可以看出,现有技术中的频域声能量对比控制方法的宽带聚焦效果最差,在某些频率点对比聚焦性能急剧下降,只有有限的控制点取得较好的效果,而现有技术中的时域声能量对比控制方法可以部分解决对比聚焦性能下降的问题,经过对比,可以看出,本发明的方法具有最好的聚焦效果。这表明相比于现有技术中的声能量对比控制方法,本方法更适合处理宽带信号输入。
图5分别给出了本方法和现有技术中的声能量对比控制方法在明区中心控制点处的频率响应。从图中可以看出,现有技术中的频域声能量控制方法振荡幅度大,而现有技术中的时域声能量控制方法则最差,整个频谱趋向于线谱,而采用本发明方法的频率响应最平坦,因此能极大改善音质。
改善音质可以通过PESQ(Perceptual evaluation of speech quality)分数来评估,PESQ的分值为0到4.5分之间,0代表最差,4.5代表最好,处理后的语音与原语音只有纯延迟。这里采用的两个语音信号采样率为8kHz,分别命名为e1和e2,e1和e2的内容分别是“the birch canoe slid on the smooth planks”和“glue the sheet to thedark blue background”,由男声朗读。表1给出了明区中心控制点处产生的语音音质PESQ分数,从表中可以看出,现有技术中的声能量对比控制方法对语音音质有损害,而相比于现有技术中的声能量对比控制方法,本发明的方法极大改善了语音音质,几乎可以将语音音质恢复成原始语音的音质。
表1明区中心控制点处产生的语音音质评价,采用PESQ分数
本发明方法 | 现有技术方法(时域) | 现有技术方法(频域) | |
e1 | 4.270 | 2.907 | 3.040 |
e2 | 4.305 | 2.824 | 3.030 |
本实施例中虽然限定了采样频率为8kHz,并选定明区和暗区为线型区域,但这仅仅是对本发明所提供方法的一个举例说明,并不限定本发明所提供方法仅适用人说话声频率范围内,或者明区,暗区只能选择线型。事实上,本发明所提供方法能够拓展到整个可听声频率段的宽带信号并实现声能量对比聚焦。
本发明还提供了一种与方法相对应的***,包括:
扬声器阵列布放模块,用于布放扬声器阵列,并设定明区和暗区的控制点;其中,所述的明区是指设定为具有高声能量的区域,所述的暗区是指设定为具有低声能量的区域;
频率响应一致性约束表达式获取模块,用于选定参考频率,列出明区中的频率响应一致性约束表达式;
平均声能量表达式获取模块,用于分别列出明区和暗区的平均声能量表达式;
时域冲激响应滤波器信号计算模块,用于给合所述频率响应一致性约束表达式获取模块所得到的频率响应一致性约束表达式和平均声能量表达式获取模块所得到的平均声能量表达式,按照频率响应一致性约束的时域声能量对比控制准则,计算出各个通道的时域冲激响应滤波器信号。
最后所应说明的是,以上实施例仅用以说明本发明的技术方案而非限制。尽管参照实施例对本发明进行了详细说明,本领域的普通技术人员应当理解,对本发明的技术方案进行修改或者等同替换,都不脱离本发明技术方案的精神和范围,其均应涵盖在本发明的权利要求范围当中。
Claims (8)
1.一种时域声能量对比控制方法,包括:
步骤1)、布放扬声器阵列,设定明区和暗区的控制点;其中,所述的明区是指设定为具有高声能量的区域,所述的暗区是指设定为具有低声能量的区域;
步骤2)、选定参考频率,列出明区中的频率响应一致性约束表达式;
步骤3)、分别列出明区和暗区的平均声能量表达式;
步骤4)、给合步骤2)的频率响应一致性约束表达式和步骤3)中的平均声能量表达式,按照频率响应一致性约束的时域声能量对比控制准则,计算出各个通道的时域冲激响应滤波器信号。
2.根据权利要求1所述的时域声能量对比控制方法,其特征在于,在所述的步骤1)中,所布放的扬声器阵列为线性阵列或圆形阵列或随机阵列。
3.根据权利要求1所述的时域声能量对比控制方法,其特征在于,在所述的步骤1)中,所述明区或暗区的区域形状是方形或圆形或线型。
4.根据权利要求1所述的时域声能量对比控制方法,其特征在于,所述步骤2)包括:
步骤2-1)、获得扬声器第l通道到明区的第k个控制点在频率f处的传递函数gBlk(f),其中,1≤l≤L,1≤k≤K,K和L分别为明区的控制点个数和扬声器通道个数;
步骤2-2)、根据步骤2-1)所得到的传递函数gBlk(f),列出扬声器阵列到明区第k个控制点的频率响应pBk(f),其表达式如下:
其中Ts是采样时间间隔,也是采样频率fs的倒数,wl(m)是对应第l通道扬声器的时域冲激响应滤波器的待求的第m个系数,M是时域冲激响应滤波器的长度;
步骤2-3)、将步骤2-2)所得到的频率响应pBk(f)写成向量形式:
pBk(f)=wTsBk(f)
其中w为待求的ML×1向量,定义为:
w=[w1(0),…,w1(M-1),…,wL(0),…,wL(M-1)]T
sBk(f)为ML×1向量,定义为:
步骤2-4)、选定参考频率fr,根据之前得到的频率响应,定义明区的频率响应一致性约束RV,该频率响应一致性约束的表达式为:
其中是取元素的实数部分,Ω是所有约束频率点的集合fc是集合Ω中的一个约束频率点,是取不超过元素的最大整数,BΩ是约束频率点的个数;Q的表达式为:
5.根据权利要求4所述的时域声能量对比控制方法,其特征在于,在所述的步骤2-1)中,采用诸如B&K PULSE的音频测试仪器或者通过建模仿真获得所述扬声器第l通道到明区的第k个控制点在频率f处的传递函数gBlk(f)。
6.根据权利要求4所述的时域声能量对比控制方法,其特征在于,所述步骤3)包括:
步骤3-1)、将扬声器第l通道到明区的第k个控制点之间的冲激响应建模成一个长度为I的FIR滤波器,并采用诸如B&K PULSE的音频测试仪器或者通过建模仿真获得其系数,定义该FIR滤波器第i个系数为hBlk(i);
步骤3-2)、假设输入给扬声器***的信号为x(n),则在明区第k个控制点的输出信号为
其中,rBk是ML×1向量,可定义如下:
rBk(n)=[rBlk(n),…rBlk(n-M+1),…,rBLk(n),…,rBLk(n-M+1)]T
步骤3-3)、假设输入信号x(n)为稳态信号,则在步骤3-2)中输出信号yBk(n)表达式的基础上,定义明区的平均声能量为:
其中E{·}代表取元素的期望值,RB(n)表达式如下:
步骤3-4)、将扬声器第l通道到暗区的第kD个控制点之间的响应建模成一个长度为I的FIR滤波器,则其第i个系数定义为假设输入给扬声器***的信号为x(n),则在暗区第kD个控制点的输出信号为
其中,1≤kD≤KD,KD为暗区控制点个数,是ML×1向量:
定义暗区的平均声能量为:
其中
7.根据权利要求6所述的时域声能量对比控制方法,其特征在于,所述的步骤4)包括:
步骤4-1)、结合步骤2)给出的频率响应一致性约束的表达式和步骤3)给出平均声能量表达式,按照频率响应一致性约束的时域声能量对比控制准则,列出其优化问题为:
步骤4-2)、对步骤4-1)所得到的优化问题求出其解:
其中Pmax{·}是求解矩阵对应最大特征值的单位特征向量,U为单位矩阵,δ为鲁棒性参数,β为权重参数;参数δ和β都取正数;
步骤4-3)、将步骤4-2)所得到的向量w按每隔M个元素进行切割,得到各个通道的时域冲激响应滤波器信号。
8.一种时域声能量对比控制***,其特征在于,包括:
扬声器阵列布放模块,用于布放扬声器阵列,并设定明区和暗区的控制点;其中,所述的明区是指设定为具有高声能量的区域,所述的暗区是指设定为具有低声能量的区域;
频率响应一致性约束表达式获取模块,用于选定参考频率,列出明区中的频率响应一致性约束表达式;
平均声能量表达式获取模块,用于分别列出明区和暗区的平均声能量表达式;
时域冲激响应滤波器信号计算模块,用于给合所述频率响应一致性约束表达式获取模块所得到的频率响应一致性约束表达式和平均声能量表达式获取模块所得到的平均声能量表达式,按照频率响应一致性约束的时域声能量对比控制准则,计算出各个通道的时域冲激响应滤波器信号。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310001166.0A CN103916810B (zh) | 2013-01-05 | 2013-01-05 | 一种时域声能量对比控制方法和*** |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310001166.0A CN103916810B (zh) | 2013-01-05 | 2013-01-05 | 一种时域声能量对比控制方法和*** |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103916810A true CN103916810A (zh) | 2014-07-09 |
CN103916810B CN103916810B (zh) | 2016-03-02 |
Family
ID=51042128
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310001166.0A Active CN103916810B (zh) | 2013-01-05 | 2013-01-05 | 一种时域声能量对比控制方法和*** |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103916810B (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104469595A (zh) * | 2014-10-30 | 2015-03-25 | 苏州上声电子有限公司 | 一种基于误差模型的多区域声重放方法和装置 |
CN108476373A (zh) * | 2016-01-27 | 2018-08-31 | 华为技术有限公司 | 一种处理声场数据的方法和装置 |
CN109275084A (zh) * | 2018-09-12 | 2019-01-25 | 北京小米智能科技有限公司 | 麦克风阵列的测试方法、装置、***、设备和存储介质 |
CN110050471A (zh) * | 2016-12-07 | 2019-07-23 | 迪拉克研究公司 | 相对于亮地带和暗地带优化的音频预补偿滤波器 |
CN116684784A (zh) * | 2023-06-29 | 2023-09-01 | 中国科学院声学研究所 | 一种基于参量阵扬声器阵列的声重放方法及*** |
CN117119092A (zh) * | 2023-02-22 | 2023-11-24 | 荣耀终端有限公司 | 一种音频处理方法及电子设备 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2005004533A1 (en) * | 2003-07-01 | 2005-01-13 | University Of Southampton | Sound reproduction systems for use by adjacent users |
CN101426171A (zh) * | 2007-10-31 | 2009-05-06 | 株式会社东芝 | 声场控制方法以及*** |
CN102611962A (zh) * | 2012-02-09 | 2012-07-25 | 苏州上声电子有限公司 | 基于无线通信网络的扬声器阵列空域响应均衡与声场控制方法和装置 |
-
2013
- 2013-01-05 CN CN201310001166.0A patent/CN103916810B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2005004533A1 (en) * | 2003-07-01 | 2005-01-13 | University Of Southampton | Sound reproduction systems for use by adjacent users |
CN101426171A (zh) * | 2007-10-31 | 2009-05-06 | 株式会社东芝 | 声场控制方法以及*** |
CN102611962A (zh) * | 2012-02-09 | 2012-07-25 | 苏州上声电子有限公司 | 基于无线通信网络的扬声器阵列空域响应均衡与声场控制方法和装置 |
Non-Patent Citations (2)
Title |
---|
JOUNG-WOO ETC.: "Generation of an acoustically bright zone with an illuminated region using multiple sources", 《2002 ACOUSTICAL SOCIETY OF AMERICA》 * |
YONG ZHAO ETC.: "Robust Wideband Beamforming With Frequency Response Variation Constraint Subject to Arbitrary Norm-Bounded Error", 《IEEE TRANSACTIONS ON ANTENNAS AND PROPAGATION》 * |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104469595A (zh) * | 2014-10-30 | 2015-03-25 | 苏州上声电子有限公司 | 一种基于误差模型的多区域声重放方法和装置 |
WO2016065719A1 (zh) * | 2014-10-30 | 2016-05-06 | 苏州上声电子有限公司 | 一种基于误差模型的多区域声重放方法和装置 |
US9992604B2 (en) | 2014-10-30 | 2018-06-05 | Suzhou Sonavox Electronics Co., Ltd. | Error model-based multi-zone sound reproduction method and device |
CN108476373B (zh) * | 2016-01-27 | 2020-11-17 | 华为技术有限公司 | 一种处理声场数据的方法和装置 |
CN108476373A (zh) * | 2016-01-27 | 2018-08-31 | 华为技术有限公司 | 一种处理声场数据的方法和装置 |
CN110050471A (zh) * | 2016-12-07 | 2019-07-23 | 迪拉克研究公司 | 相对于亮地带和暗地带优化的音频预补偿滤波器 |
CN110050471B (zh) * | 2016-12-07 | 2022-01-21 | 迪拉克研究公司 | 相对于亮地带和暗地带优化的音频预补偿滤波器 |
US11246000B2 (en) | 2016-12-07 | 2022-02-08 | Dirac Research Ab | Audio precompensation filter optimized with respect to bright and dark zones |
CN109275084A (zh) * | 2018-09-12 | 2019-01-25 | 北京小米智能科技有限公司 | 麦克风阵列的测试方法、装置、***、设备和存储介质 |
CN117119092A (zh) * | 2023-02-22 | 2023-11-24 | 荣耀终端有限公司 | 一种音频处理方法及电子设备 |
CN117119092B (zh) * | 2023-02-22 | 2024-06-07 | 荣耀终端有限公司 | 一种音频处理方法及电子设备 |
CN116684784A (zh) * | 2023-06-29 | 2023-09-01 | 中国科学院声学研究所 | 一种基于参量阵扬声器阵列的声重放方法及*** |
CN116684784B (zh) * | 2023-06-29 | 2024-03-12 | 中国科学院声学研究所 | 一种基于参量阵扬声器阵列的声重放方法及*** |
Also Published As
Publication number | Publication date |
---|---|
CN103916810B (zh) | 2016-03-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Mertins et al. | Room impulse response shortening/reshaping with infinity-and $ p $-norm optimization | |
CN103916810B (zh) | 一种时域声能量对比控制方法和*** | |
CN101622669B (zh) | 用于信号分离的***、方法及设备 | |
CN102947685B (zh) | 用于减少环境噪声对收听者的影响的方法和装置 | |
CN101883304B (zh) | 用于声复制的补偿*** | |
CN103856866B (zh) | 低噪微分麦克风阵列 | |
CN106535076B (zh) | 一种立体声音响***的空间校准方法及其移动终端设备 | |
CN110557710B (zh) | 具有语音控制的低复杂度多声道智能扩音器 | |
CN103634726A (zh) | 一种扬声器自动均衡方法 | |
CN104170408A (zh) | 应用组合的或混合的声场控制策略的方法 | |
WO2016065719A1 (zh) | 一种基于误差模型的多区域声重放方法和装置 | |
US9264812B2 (en) | Apparatus and method for localizing a sound image, and a non-transitory computer readable medium | |
Møller et al. | A moving horizon framework for sound zones | |
CN106162482A (zh) | 一种用于抑制啸叫的扬声器阵列扩声***及方法 | |
CN103916730B (zh) | 一种能够改善音质的声场聚焦方法及*** | |
CN103916733B (zh) | 基于均方误差最小准则的声能量对比控制方法及*** | |
Clifford et al. | Proximity effect detection for directional microphones | |
CN113763984B (zh) | 一种用于分布式多说话人的参数化噪声消除*** | |
JP6567216B2 (ja) | 信号処理装置 | |
Chang et al. | Distributed parameterized topology-independent noise reduction in acoustic sensor networks | |
Betlehem et al. | Temperature robust active-compensated sound field reproduction using impulse response shaping | |
Tuna et al. | Data-driven local average room transfer function estimation for multi-point equalization | |
CN105792075B (zh) | 一种串声消除滤波器的生成方法及三维声音重放方法 | |
JP5948259B2 (ja) | 収音装置及び再生装置 | |
CN114915875B (zh) | 一种可调波束形成方法、电子设备及存储介质 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |