发明名称 一种声场定量重现控制系统的扬声器回放阵列控制信号的处理方法
摘要 本发明涉及一种声场定量重现的控制系统及其方法,属于声场重现及控制技术领域;该系统包括:扬声器回放阵列、采样阵列、多通道数模转换器、存储有信号处理程序的计算机和多通道数据采集卡;该方法包括:对目标声场进行采样:对采样信号进行频谱分析:通过傅里叶信号处理方法对目标声压信号进行频谱分析,得到该目标声压信号的频谱结构,提取频谱结构中每个频率成分的幅值和相位,将该目标声压信号转换到频域量:测量扬声器回放阵列对采样阵列的电声传递函数,基于传递函数识别方法求解扬声器回放阵列控制信号,将该控制信号输入回放扬声器回放阵列执行回放;该方法使得回放声场与目标声场误差较小。
申请公布号 CN102857852B 申请公布日期 2014.10.22
申请号 CN201210337089.1 申请日期 2012.09.12
申请人 清华大学 发明人 郑四发;李申广;彭博;连小珉
分类号 H04R29/00(2006.01)I 主分类号 H04R29/00(2006.01)I
代理机构 北京清亦华知识产权代理事务所(普通合伙) 11201 代理人 廖元秋
主权项 一种声场定量重现控制系统的扬声器回放阵列控制信号的处理方法,该系统包括:扬声器回放阵列、采样阵列、多通道数模转换器、存储有信号处理程序的计算机和多通道数据采集卡;其中,采样阵列由多个传声器组成,位于扬声器回放阵列的中心区域;每个传声器为一个采样测点,每个传声器的输出端均与多通道数据采集卡的输入端相连,多通道数据采集卡的输出端与计算机的输入端相连,计算机的输出端与多通道数模转换器相连,多通道数模转换器的输出端同时与扬声器回放阵列中的各扬声器的输入端和多通道数据采集卡的另一个输入端相连;其特征在于,该方法包括以下步骤:1)对目标声场进行采样:设在限定区域内布置M个采样测点,分别用x<sub>1</sub>,x<sub>2</sub>,x<sub>3</sub>,…,x<sub>M</sub>表示其空间位置,由此可获得目标时段内M个采样测点的目标声压信号[P<sub>d</sub>(t)]<sub>M×1</sub>,表示为:[P<sub>d</sub>(t)]<sub>M×1</sub>=[P<sub>d</sub>(x<sub>1</sub>;t),P<sub>d</sub>(x<sub>2</sub>;t),...,P<sub>d</sub>(x<sub>M</sub>;t)]<sup>T</sup>   (1)式(1)中t为时间变量,[·]<sup>T</sup>表示矩阵的转置;2)对采样信号进行频谱分析:通过傅里叶信号处理方法对目标声压信号进行频谱分析,得到该目标声压信号的频谱结构,提取频谱结构中每个频率成分的幅值和相位,将该目标声压信号[P<sub>d</sub>(t)]<sub>M×1</sub>转换到频域量[P<sub>d</sub>(f)]<sub>M×1</sub>:[P<sub>d</sub>(f)]<sub>M×1</sub>=[P<sub>d</sub>(x<sub>1</sub>;f),P<sub>d</sub>(x<sub>2</sub>;f),...,P<sub>d</sub>(x<sub>M</sub>;f)]<sup>T</sup>   (2)式(2)中f为频率变量;3)测量扬声器回放阵列对采样阵列的电声传递函数:设扬声器回放阵列由L个扬声器组成,首先给1号扬声器一个已知信号激励,测量扬声器回放阵列11中的1号扬声器对采样阵列12中的第1个采样测点的电声传递函数,然后同步测量L通道数模转换器的输出端对应通道的电信号e(t)以及同时刻采用采样阵列的采样测点1处的声压信号p(t),通过所述电信号e(t)和声压信号p(t)两个信号计算出1号扬声器对采样测点1的电声传递函数,依此类推,分别测量扬声器回放阵列上的L个扬声器对采样阵列上M个采样测点的电声传递函数,从而得到传递函数矩阵[H(x;f)]<sub>M×L</sub>:<img file="FDA0000526298580000021.GIF" wi="1541" he="328" />采用分频段激励和测量,激励电信号e(t)中每个频率成分幅值均为A<sub>0</sub>,如式(4):<img file="FDA0000526298580000022.GIF" wi="1246" he="139" />其中f<sub>k</sub>为频率变量、<img file="FDA0000526298580000023.GIF" wi="58" he="59" />为相位变量;对激励电压信号e(t)进行傅里叶变换时进行整周期截取,计算传递函数H<sub>lm</sub>(x;f);依此,对每个扬声器分别进行传递函数的测量,可得到传递函数矩阵[H(x;f)]<sub>M×L</sub>;4)基于传递函数识别方法求解扬声器回放阵列控制信号,设扬声器回放阵列输入控制信号中f频率分量由<img file="FDA0000526298580000024.GIF" wi="184" he="88" />表示,则声场控制系统模型如式(7)所示:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mrow><mo>[</mo><mi>H</mi><mrow><mo>(</mo><mi>x</mi><mo>;</mo><mi>f</mi><mo>)</mo></mrow><mo>]</mo></mrow><mrow><mi>M</mi><mo>&times;</mo><mi>L</mi></mrow></msub><mo>&CenterDot;</mo><msub><mrow><mo>[</mo><mover><mi>e</mi><mo>^</mo></mover><mrow><mo>(</mo><mi>f</mi><mo>)</mo></mrow><mo>]</mo></mrow><mrow><mi>L</mi><mo>&times;</mo><mn>1</mn></mrow></msub><mo>=</mo><msub><mrow><mo>[</mo><msub><mi>P</mi><mi>d</mi></msub><mrow><mo>(</mo><mi>f</mi><mo>)</mo></mrow><mo>]</mo></mrow><mrow><mi>M</mi><mo>&times;</mo><mn>1</mn></mrow></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>7</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000526298580000025.GIF" wi="1511" he="101" /></maths>使用最小误差平方法求解<img file="FDA0000526298580000026.GIF" wi="184" he="90" />使得<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mo>|</mo><mo>|</mo><msub><mrow><mo>[</mo><mi>H</mi><mrow><mo>(</mo><mi>x</mi><mo>;</mo><mi>f</mi><mo>)</mo></mrow><mo>]</mo></mrow><mrow><mi>M</mi><mo>&times;</mo><mi>L</mi></mrow></msub><mo>&CenterDot;</mo><msub><mrow><mo>[</mo><mover><mi>e</mi><mo>^</mo></mover><mrow><mo>(</mo><mi>f</mi><mo>)</mo></mrow><mo>]</mo></mrow><mrow><mi>L</mi><mo>&times;</mo><mn>1</mn></mrow></msub><mo>-</mo><msub><mrow><mo>[</mo><msub><mi>P</mi><mi>d</mi></msub><mrow><mo>(</mo><mi>f</mi><mo>)</mo></mrow><mo>]</mo></mrow><mrow><mi>M</mi><mo>&times;</mo><mn>1</mn></mrow></msub><mo>|</mo><mo>|</mo></mrow>]]></math><img file="FDA0000526298580000027.GIF" wi="910" he="110" /></maths>最小,得到正则化最小误差平方的扬声器回放阵列控制信号f频率分量如式(8)所示:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><mrow><mo>[</mo><mover><mi>e</mi><mo>^</mo></mover><mrow><mo>(</mo><mi>f</mi><mo>)</mo></mrow><mo>]</mo></mrow><mrow><mi>L</mi><mo>&times;</mo><mn>1</mn></mrow></msub><msup><mrow><mo>[</mo><msup><msub><mrow><mo>[</mo><mi>H</mi><mrow><mo>(</mo><mi>x</mi><mo>;</mo><mi>f</mi><mo>)</mo></mrow><mo>]</mo></mrow><mrow><mi>M</mi><mo>&times;</mo><mi>L</mi></mrow></msub><mi>T</mi></msup><msub><mrow><mo>[</mo><mi>H</mi><mrow><mo>(</mo><mi>x</mi><mo>;</mo><mi>f</mi><mo>)</mo></mrow><mo>]</mo></mrow><mrow><mi>M</mi><mo>&times;</mo><mi>L</mi></mrow></msub><mo>+</mo><mi>&lambda;I</mi><mo>]</mo></mrow><mrow><mo>-</mo><mn>1</mn></mrow></msup><msup><msub><mrow><mo>[</mo><mi>H</mi><mrow><mo>(</mo><mi>x</mi><mo>;</mo><mi>f</mi><mo>)</mo></mrow><mo>]</mo></mrow><mrow><mi>M</mi><mo>&times;</mo><mi>L</mi></mrow></msub><mi>T</mi></msup><msub><mrow><msub><mrow><mo>[</mo><mi>P</mi></mrow><mi>d</mi></msub><mrow><mo>(</mo><mi>f</mi><mo>)</mo></mrow><mo>]</mo></mrow><mrow><mi>M</mi><mo>&times;</mo><mn>1</mn></mrow></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>8</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000526298580000028.GIF" wi="1715" he="121" /></maths>式(8)中λ为正则化参数,I为单位矩阵,[·]<sup>‑1</sup>表示矩阵的逆;将各个频率分量的扬声器回放阵列控制信号<img file="FDA0000526298580000029.GIF" wi="188" he="85" />叠加即可得到最后总的扬声器回放阵列控制信号<img file="FDA00005262985800000210.GIF" wi="129" he="85" />5)将步骤4)中求解得到的控制信号<img file="FDA00005262985800000211.GIF" wi="92" he="86" />输入回放扬声器回放阵列执行回放,即实现对目标声场的重现。
地址 100084 北京市海淀区清华园1号