[发明专利]一种空间场主动降噪方法在审
申请号: | 201810725082.4 | 申请日: | 2018-07-04 |
公开(公告)号: | CN108877758A | 公开(公告)日: | 2018-11-23 |
发明(设计)人: | 徐银海;刘益帆 | 申请(专利权)人: | 北京安声科技有限公司 |
主分类号: | G10K11/16 | 分类号: | G10K11/16 |
代理公司: | 江苏致邦律师事务所 32230 | 代理人: | 闫东伟;毛禾枫 |
地址: | 100000 北京市朝阳区*** | 国省代码: | 北京;11 |
权利要求书: | 查看更多 | 说明书: | 查看更多 |
摘要: | |||
搜索关键词: | 降噪 噪声声场 声场模型 控制器参数 整体空间 基函数 三维空间 采集 连续函数 模型计算 频率声音 噪声水平 整体能量 主动降噪 空间场 构建 声场 噪声 抵消 | ||
1.一种空间场主动降噪方法,其特征在于,包括以下步骤:
S1、确定降噪区域,在需要降噪的空间内安装多个传入装置及输出装置;
S2、计算降噪系统基函数的系数,误差麦克风阵列采集空间的噪声信号,计算麦克风阵列采集到的对应监测点的统计声压数据阵列即噪声信号为e(x,f),x为在监测点该空间中的位置坐标,f为是监测点采集到信号的频率域统计声压值,e(x,f)用下列公式来表示:
式中φk(x,f)是根据降噪系统的使用环境选择的基函数,管道声场环境选用声音空腔模态,开放声场选用球面波,是一组固定的基函数,Ak(f)是基函数的系数,麦克风阵列采集到的噪声信号用来表示,的公式为:
即根据测得的噪声信号,通过上述公式推导出基函数系数的计算公式从而算出Ak(f),公式如下:
式中的H代表厄米特矩阵,-1代表逆矩阵,Ak(f)是通过上文的e(x,f)公式算出,因为基函数固定,通过测得的信号反向算出,计算出Ak(f)和B(f)后得到上述公式的傅里叶变换即
其中式中e(f)是噪声经过反向声波抵消后,误差麦克风采集到噪声信号,d(f)是噪声在误差麦克风处的信号,W(f)是ANCfilter的传递函数,v(f)为一个换算过程中的中间变量,F为傅里叶变换,是用时间的函数来表示该参数;
S3、启动降噪系统,对空间声场各点的能量进行计算,先对依据建模计算出的声场各点声音进行有针对性的主动降噪,此时信号如下:
上述e(f)是一组向量,而e(n)是实际执行降噪时第l处麦克风的值,e(f)包含e(n),式中l代表某个误差麦克风的编号,则el(n)表示在n时刻第l个误差麦克风处的信号,表示在n时刻第l个误差麦克风处的原始噪声源信号,加号后面为系统发出的反相声波信号,m为降噪扬声器的编号,其中gl,m,j中的g代表滤波器参数,j代表滤波器阶数,该系数通过前期测量测得,u代表驱动扬声器发出降噪声波的信号,um代表第m个扬声器的驱动信号,则该式表示的是,第l个误差麦克风处接收到的已经进行降噪后的噪声信号,本式代表一个麦克风附近的噪声信号经过公式计算出的声场环境能量,如果实际系统中有几个麦克风,需要将几个麦克风的信号计算出的声场环境能量进行相加并据此计算出声场环境的整体能量期望值,整体能量结合上述和的计算公式用以下公式表示:
式中的T指矩阵的转置,E代表的是随机变量的期望值;
S4、对步骤S3的J取最小值求得参数W,取得的系数W即为使该声场环境下J最小的系数,也就是使整体声场能量最低的控制器系数,得到W后即完成对控制器的设计,进而完成前期测量,正式开始降噪过程;
S5、将步骤S4求得的控制器参数W输入,主动降噪系统根据麦克风阵列采集的噪声信号计算出相应的反相声波,控制单元发出的反相声波信号结合扬声器到降噪区域的传递函数后进行降噪处理,麦克风阵列继续对经过降噪处理的噪声信号进行采集并由主动降噪系统继续进行降噪处理。
2.根据权利要求1所述的降噪方法,其特征在于,在步骤S1之后进行系统自检,降噪系统的自检模块接收并处理自身系统信息来判定出是否有传入装置和输出装置的增减,自检完成后降噪系统启动整体的计算过程,转入下一步骤开始计算最佳的控制器参数;系统自检模块内设有定时系统检测模块,检测到系统信息发生变化后会进行系统矫正,重新进行自检,如果有新增或新减少的传入装置或输出装置,系统会以新的装置数量为标准。
该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于北京安声科技有限公司,未经北京安声科技有限公司许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服】
本文链接:http://www.vipzhuanli.com/pat/books/201810725082.4/1.html,转载请声明来源钻瓜专利网。
- 上一篇:用于声学辐射控制的系统和方法
- 下一篇:一种用于家用电器的消音装置