CN106446532B - 一种虚拟源强的配置方法 - Google Patents
一种虚拟源强的配置方法 Download PDFInfo
- Publication number
- CN106446532B CN106446532B CN201610813656.4A CN201610813656A CN106446532B CN 106446532 B CN106446532 B CN 106446532B CN 201610813656 A CN201610813656 A CN 201610813656A CN 106446532 B CN106446532 B CN 106446532B
- Authority
- CN
- China
- Prior art keywords
- sound source
- source
- point sound
- barycenter
- point
- 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.)
- Expired - Fee Related
Links
Classifications
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16Z—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS, NOT OTHERWISE PROVIDED FOR
- G16Z99/00—Subject matter not provided for in other main groups of this subclass
Landscapes
- Stereophonic System (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
本发明旨在提供一种虚拟源强的配置方法:A、在目标区域建立坐标系,建立目标区域的声压图像,将声压图像中对应的灰度值作为对应点声源的声压值;B、将各点声源的声压值作为权值通过重心法求出目标区域声场系统的质心球虚拟源强的球心坐标;C、将目标区域划分为多个单位体积,计算出各个单位体积的声源聚集度,根据声源聚集度进行归类分析,根据分类得到多个非质心虚拟源强;D、结合各点声源的声压值作为权值通过重心法求出各个非质心球虚拟源强的球心坐标;E、自定义质心球虚拟源强与非质心球虚拟源强的半径,完成虚拟源强的配置。该配置方法克服现有技术声场重建精度不高的缺陷,具有精度高、更适于结构形状振动体适用性更强的特点。
Description
技术领域
本发明涉及声场重建领域,具体涉及一种虚拟源强的配置方法。
背景技术
在声场重建领域,波叠加法是常用的手段。波叠加法是由Koopmann等最先提出的,其基本思想是:物体辐射的声场可表示为,置于该辐射体内部一系列简单源(单极子、偶极子等)产生的声场线性叠加,来拟合实际声场。现有文献《基于多球域波叠加法的Patch近场声全息》中提出了基于多球域波叠加法的PNAH,该方法采用多球形虚拟源强配置的灵活性,改进了共形问题影响重建精度的困难,提高了重建精度和计算效率
根据所采用等效源的不同,虚拟源强源配置方法可以分为以下两类:
(1)采用简单源(单极子或偶极子等)作为等效源的NAH技术。采用此种方法时,等效源分布在实际声源表面背离分析域一定距离的虚源面上,且等效源的数量等于表面节点数时数据仿真结果最佳。
(2)采用球面波作为等效源的NAH技术。通过一系列不同阶次的球面波源加权组合来近似实际声源所辐射的外部声场。
现有声源分析中,采用的声源分布并不是随机的,是按照一定的规律进行分布的;然而对于任何结构形状振动体,由于声源产生是随机分布的,那么就需要根据不同的随机分布特性,进行实际的虚拟声源的布置。
发明内容
本发明旨在提供一种虚拟源强的配置方法,该配置方法克服现有技术声场重建精度不高的缺陷,具有精度高、更适于结构形状振动体适用性更强的特点。
本发明的技术方案如下:一种虚拟源强的配置方法,包括以下步骤:
A、在需要进行虚拟源强配置的目标区域建立坐标系,基于坐标系获得各点声源的坐标,建立目标区域的声压图像,将声压图像中对应的灰度值作为对应点声源的声压值;所述的声压图像通过结合目标区域的声场情况建模模拟测量得到;
B、将各点声源的声压值作为权值通过重心法求出目标区域声场系统的重心坐标,作为质心球虚拟源强的球心坐标;
C、将目标区域划分为多个单位体积,计算出各个单位体积的声源聚集度,根据声源聚集度进行归类分析,将相邻区域聚集度比较高的进行归类,每一类点声源作为一个非质心虚拟源强;将区域聚集度不高的点声源进行分散处理,每一个点声源作为一个非质心虚拟源强;
D、将各个非质心虚拟源强中各点声源的声压值作为权值通过重心法求出各个非质心虚拟源强的重心坐标,作为非质心球虚拟源强的球心坐标;
E、自定义质心球虚拟源强与非质心球虚拟源强的半径,完成虚拟源强的配置。
优选地,所述的目标区域为待重建声场对应的不规则振动体区域。
优选地,所述的步骤C具体为:将以各点声源坐标为中心,以预设的搜索半径为半径的球形区域作为单位体积,计算各点声源对应的单位体积的声源聚集度;将最大声源聚集度对应的单位体积内的点声源划分为一类;对剩余各点声源对应的单位体积的声源聚集度计算,将最大声源聚集度对应的单位体积内的点声源划分为新的一类;循环进行上述对剩余各点声源的声源聚集度的计算及分类,直至全部点声源都完成分类。
优选地,步骤B具体如下:
利用式(1)求出整个系统的重心作为质心球虚拟源强球心坐标;
其中(xi,yi,zi)为各点声源的坐标,f(xi,yi,zi)为点声源(xi,yi,zi)的声压值,v为点声源的数量。
优选地,步骤C具体如下:
C1、预设搜索半径为d,将以各点声源坐标为中心,以d为半径的球形区域作为单位体积,利用式(2)计算各点声源对应的单位体积的声源聚集度;
式中,F为单位体积的声源聚集度,a为单位体积内的点声源数量,A为声场系统内总的点声源数量,b为单位体积内点声源的声压值之和,B为声场系统内点声源的声压值之和;
C2、将最大声源聚集度对应的单位体积内的点声源划分为一类;对剩余各点声源对应的单位体积的声源聚集度计算,将最大声源聚集度对应的单位体积内的点声源划分为新的一类;重复进行上述对剩余的声源聚集度的计算及分类,直至全部点声源都完成分类。
优选地,步骤D具体如下:
将属于同一类别点声源构成一个子系统,将各个子系统中各点声源的声压值作为权值通过重心法求出各个子系统的重心坐标,作为非质心球虚拟源强的球心坐标;
公式如下:
式中为各个子系统的重心坐标,j为子系统的数量,1≤j<v;w为对应子系统中点声源的数量。
本发明采用多球形虚源的简单源进行研究,有效解决了球面波有限阶数截断的问题,克服了采用简单源配置中与声源表面的共形问题;同时,引入了声源空间聚集度的聚类分析法,有限地针对声源分布的随机性来配置虚拟源强,更适用于任何结构形状振动体应用,为后续声场重建的高精确度打下基础。
附图说明
图1为本发明提供的虚拟源强的配置方法的流程图
图2为实施例1的点声源及坐标建立图
图3为实施例1对应声场重建生成的L曲线
图4为实施例1对应声场重建结果与误差图。
具体实施方式
下面结合附图和实施例具体说明本发明。
实施例1
如图1所示,本实施例提供的一种虚拟源强的配置方法,包括以下步骤:
A、在需要进行虚拟源强配置的目标区域建立坐标系,基于坐标系获得各点声源的坐标,建立目标区域的声压图像,将声压图像中对应的灰度值作为对应点声源的声压值;
B、将各点声源的声压值作为权值通过重心法求出目标区域声场系统的重心坐标,作为质心球虚拟源强的球心坐标;
具体如下:
利用式(1)求出整个系统的重心作为质心球虚拟源强球心坐标;
其中(xi,yi,zi)为各点声源的坐标,f(xi,yi,zi)为点声源(xi,yi,zi)的声压值,v为点声源的数量;
C、预设搜索半径,将以各点声源坐标为中心,以搜索半径为半径的球形区域作为单位体积,计算各点声源对应的单位体积的声源聚集度;将最大声源聚集度对应的单位体积内的点声源划分为一类;对剩余各点声源对应的单位体积的声源聚集度计算,将最大声源聚集度对应的单位体积内的点声源划分为新的一类;循环进行上述对剩余各点声源的声源聚集度的计算及分类,直至全部点声源都完成分类;
具体如下:
C1、预设搜索半径为d,将以各点声源坐标为中心,以d为半径的球形区域作为单位体积,利用式(2)计算各点声源对应的单位体积的声源聚集度;
式中,F为单位体积的声源聚集度,a为单位体积内的点声源数量,A为声场系统内总的点声源数量,b为单位体积内点声源的声压值之和,B为声场系统内点声源的声压值之和;
C2、将最大声源聚集度对应的单位体积内的点声源划分为一类;对剩余各点声源对应的单位体积的声源聚集度计算,将最大声源聚集度对应的单位体积内的点声源划分为新的一类;重复进行上述对剩余的声源聚集度的计算及分类,直至全部点声源都完成分类;
D、将各个非质心虚拟源强中各点声源的声压值作为权值通过重心法求出各个非质心虚拟源强的重心坐标,作为非质心球虚拟源强的球心坐标;
步骤D具体如下:
将属于同一类别点声源构成一个子系统,将各个子系统中各点声源的声压值作为权值通过重心法求出各个子系统的重心坐标,作为非质心球虚拟源强的球心坐标;
公式如下:
式中为各个子系统的重心坐标,j为子系统的数量,1≤j<v;w为对应子系统中点声源的数量;
E、预设质心球虚拟源强与非质心球虚拟源强的半径,完成虚拟源强的配置。
本实施例预设目标区域的不规则振动体位于测量面的一侧,其中S为25个点声源组成的任意形状的不规则振动体,分布在0.4mх0.8mх0.8m的长方体内,以测量面的中心为坐标原点,测量面指向点声源S方向为x轴正方向,如图2所示;最初的测量距离d=0.8m,测量面采用的是8×8的声阵列布置方式;质心球虚拟源强的半径预设为0.025m,非质心球虚拟源强的半径预设为0.05m;
通过步骤C将25个点声源分为6类,其中类5与类6均只包含一个点声源,求出各个非质心球虚拟源强的球心坐标后,完成虚拟源强的配置;声场重建的后续步骤参照发明专利:声场重建方法(专利号:201310733451.1)中的相应步骤进行,声场重建生成的L曲线图如图3所示,重建声场的幅值及相位误差分析结果如图4所示,其中L曲线在获得最大曲率点时的曲线拟合程度高,这样更易于选择正则化参数,也证明了本实施例虚源源强配置方法的有效性。
Claims (3)
1.一种虚拟源强的配置方法,其特征在于包括以下步骤:
A、在需要进行虚拟源强配置的目标区域建立坐标系,基于坐标系获得各点声源的坐标,建立目标区域的声压图像,将声压图像中对应的灰度值作为对应点声源的声压值;
B、将各点声源的声压值作为权值通过重心法求出目标区域声场系统的重心坐标,作为质心球虚拟源强的球心坐标;
步骤B具体如下:
利用式(1)求出整个系统的重心作为质心球虚拟源强球心坐标;
其中(xi,yi,zi)为各点声源的坐标,f(xi,yi,zi)为点声源(xi,yi,zi)的声压值,v为点声源的数量;
C、将目标区域划分为多个单位体积,计算出各个单位体积的声源聚集度,根据声源聚集度进行归类分析,将相邻区域聚集度比较高的进行归类,每一类点声源作为一个非质心虚拟源强;将区域聚集度不高的点声源进行分散处理,每一个点声源作为一个非质心虚拟源强;
所述的步骤C具体为:将以各点声源坐标为中心,以预设的搜索半径为半径的球形区域作为单位体积,计算各点声源对应的单位体积的声源聚集度;将最大声源聚集度对应的单位体积内的点声源划分为一类;对剩余各点声源对应的单位体积的声源聚集度计算,将最大声源聚集度对应的单位体积内的点声源划分为新的一类;循环进行上述对剩余各点声源的声源聚集度的计算及分类,直至全部点声源都完成分类;
D、将各个非质心虚拟源强中各点声源的声压值作为权值通过重心法求出各个非质心虚拟源强的重心坐标,作为非质心球虚拟源强的球心坐标;
步骤D具体如下:
将属于同一类别点声源构成一个子系统,将各个子系统中各点声源的声压值作为权值通过重心法求出各个子系统的重心坐标,作为非质心球虚拟源强的球心坐标;
公式如下:
式中为各个子系统的重心坐标,j为子系统的数量,1≤j<v;w为对应子系统中点声源的数量;
E、自定义质心球虚拟源强与非质心球虚拟源强的半径,完成虚拟源强的配置。
2.如权利要求1所述的虚拟源强的配置方法,其特征在于:所述的目标区域为待重建声场对应的不规则振动体区域。
3.如权利要求1所述的虚拟源强的配置方法,其特征在于:
步骤C具体如下:
C1、预设搜索半径为d,将以各点声源坐标为中心,以d为半径的球形区域作为单位体积,利用式(3)计算各点声源对应的单位体积的声源聚集度;
式中,F为单位体积的声源聚集度,a为单位体积内的点声源数量,A为声场系统内总的点声源数量,b为单位体积内点声源的声压值之和,B为声场系统内点声源的声压值之和;
C2、将最大声源聚集度对应的单位体积内的点声源划分为一类;对剩余各点声源对应的单位体积的声源聚集度计算,将最大声源聚集度对应的单位体积内的点声源划分为新的一类;重复进行上述对剩余的声源聚集度的计算及分类,直至全部点声源都完成分类。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610813656.4A CN106446532B (zh) | 2016-09-09 | 2016-09-09 | 一种虚拟源强的配置方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610813656.4A CN106446532B (zh) | 2016-09-09 | 2016-09-09 | 一种虚拟源强的配置方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106446532A CN106446532A (zh) | 2017-02-22 |
CN106446532B true CN106446532B (zh) | 2018-09-21 |
Family
ID=58167567
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610813656.4A Expired - Fee Related CN106446532B (zh) | 2016-09-09 | 2016-09-09 | 一种虚拟源强的配置方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106446532B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10609485B2 (en) | 2017-09-29 | 2020-03-31 | Apple Inc. | System and method for performing panning for an arbitrary loudspeaker setup |
CN112216113B (zh) * | 2020-10-16 | 2021-07-27 | 温州职业技术学院 | 基于密度峰值的交通流量数据聚类方法 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP3986785B2 (ja) * | 2001-09-20 | 2007-10-03 | 日本放送協会 | 音源分離収音マイクロホン装置および方法 |
EP2051540A1 (en) * | 2007-10-19 | 2009-04-22 | Weistech Technology Co., Ltd. | Three-dimension array structure of surround-sound speaker |
EP2297556A1 (en) * | 2008-07-08 | 2011-03-23 | Brüel & Kjaer Sound & Vibration Measurement A/S | Method for reconstructing an acoustic field |
JP2011527422A (ja) * | 2008-07-08 | 2011-10-27 | ブリュエル アンド ケアー サウンド アンド ヴァイブレーション メジャーメント エー/エス | 音場の再構成 |
CN102901559A (zh) * | 2012-09-27 | 2013-01-30 | 哈尔滨工程大学 | 一种基于单面测量和局部声全息法的分离声场方法 |
CN103167373A (zh) * | 2011-12-09 | 2013-06-19 | 现代自动车株式会社 | 定位声源的方法和系统 |
CN103712684A (zh) * | 2013-12-25 | 2014-04-09 | 广西科技大学 | 声场重建方法 |
-
2016
- 2016-09-09 CN CN201610813656.4A patent/CN106446532B/zh not_active Expired - Fee Related
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP3986785B2 (ja) * | 2001-09-20 | 2007-10-03 | 日本放送協会 | 音源分離収音マイクロホン装置および方法 |
EP2051540A1 (en) * | 2007-10-19 | 2009-04-22 | Weistech Technology Co., Ltd. | Three-dimension array structure of surround-sound speaker |
EP2297556A1 (en) * | 2008-07-08 | 2011-03-23 | Brüel & Kjaer Sound & Vibration Measurement A/S | Method for reconstructing an acoustic field |
JP2011527422A (ja) * | 2008-07-08 | 2011-10-27 | ブリュエル アンド ケアー サウンド アンド ヴァイブレーション メジャーメント エー/エス | 音場の再構成 |
JP2011527421A (ja) * | 2008-07-08 | 2011-10-27 | ブリュエル アンド ケアー サウンド アンド ヴァイブレーション メジャーメント エー/エス | 音響場を再構成するための方法 |
CN103167373A (zh) * | 2011-12-09 | 2013-06-19 | 现代自动车株式会社 | 定位声源的方法和系统 |
CN102901559A (zh) * | 2012-09-27 | 2013-01-30 | 哈尔滨工程大学 | 一种基于单面测量和局部声全息法的分离声场方法 |
CN103712684A (zh) * | 2013-12-25 | 2014-04-09 | 广西科技大学 | 声场重建方法 |
Non-Patent Citations (2)
Title |
---|
基于多球域波叠加法的 Patch 近场声全息;何 伟等;《广 西 工 学 院 学 报》;20130331;全文 * |
基于多球域波叠加法的 PNAH 的虚源配置研究;王玉江等;《沈 阳 理 工 大 学 学 报》;20141231;全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN106446532A (zh) | 2017-02-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105069423B (zh) | 一种人体姿态检测方法及装置 | |
CN104063903B (zh) | 三维实体模型的四面体网格生成方法和装置 | |
Zhou et al. | Random generation of natural sand assembly using micro x-ray tomography and spherical harmonics | |
CN106940439B (zh) | 基于无线声传感器网络的k均值聚类加权声源定位方法 | |
CN106777822B (zh) | 二维高含量颗粒增强复合材料三相细观模型及建立方法 | |
CN103425815B (zh) | 用于高聚合物材料的模拟方法 | |
CN106446532B (zh) | 一种虚拟源强的配置方法 | |
CN104596636B (zh) | 声场分离方法 | |
CN103488806A (zh) | 用于模拟填料在高聚物材料中的分散性的方法 | |
CN104066178A (zh) | 一种基于人工神经元网络的室内无线定位指纹生成方法 | |
Barecasco et al. | Simple free-surface detection in two and three-dimensional SPH solver | |
CN103438985A (zh) | 一种用于声场合成的声场信息采集方法 | |
CN105067641A (zh) | 基于模板粒子群寻优的复杂体异物的微波检测系统 | |
Decker et al. | Geometric accuracy prediction for additive manufacturing through machine learning of triangular mesh data | |
CN104992035B (zh) | 一种太赫兹频段表面粗糙目标电磁散射的快速计算方法 | |
US10319132B2 (en) | Method and system for representing objects with velocity-dependent particles | |
US11500197B2 (en) | Efficient optical system design and components | |
CN110390124A (zh) | 确定物体形状的方法和设备、制造计算装置的方法 | |
CN111080790B (zh) | 一种基于凸包算法的凹体生成方法和装置 | |
CN105184033A (zh) | 基于阶数步进金属目标的宽频带电磁特性快速预估方法 | |
CN110763328B (zh) | 一种半空间声场重建方法和装置 | |
Menzel | Evolvable free-form deformation control volumes for evolutionary design optimization | |
Liu et al. | A novel stochastic generation method of shale blocks in S-RM considering the particle size effect of morphological features | |
CN116843862B (zh) | 一种三维薄壁模型网格表面纹理合成方法 | |
CN104777329A (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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
TR01 | Transfer of patent right |
Effective date of registration: 20200518 Address after: Room 501-1, Jingxuan building, Liuzhou Vocational and technical college, No. 30, Shewan Road, Liuzhou City, 545000 Guangxi Zhuang Autonomous Region Patentee after: Liuzhou Daoyuan Technology Co.,Ltd. Address before: 545006 No. 268 East Ring Road, the Guangxi Zhuang Autonomous Region, Liuzhou Patentee before: GUANGXI UNIVERSITY OF SCIENCE AND TECHNOLOGY |
|
TR01 | Transfer of patent right | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20180921 |
|
CF01 | Termination of patent right due to non-payment of annual fee |