CN101902680B - 基于有源控制的虚拟耦合空间的设计方法 - Google Patents
基于有源控制的虚拟耦合空间的设计方法 Download PDFInfo
- Publication number
- CN101902680B CN101902680B CN2010102168282A CN201010216828A CN101902680B CN 101902680 B CN101902680 B CN 101902680B CN 2010102168282 A CN2010102168282 A CN 2010102168282A CN 201010216828 A CN201010216828 A CN 201010216828A CN 101902680 B CN101902680 B CN 101902680B
- Authority
- CN
- China
- Prior art keywords
- sound
- space
- msub
- coupling
- sound pressure
- 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
- 238000000034 method Methods 0.000 title claims abstract description 16
- 230000008878 coupling Effects 0.000 claims description 52
- 238000010168 coupling process Methods 0.000 claims description 47
- 238000005859 coupling reaction Methods 0.000 claims description 47
- 238000011217 control strategy Methods 0.000 claims description 3
- 238000004613 tight binding model Methods 0.000 claims description 2
- 238000005259 measurement Methods 0.000 claims 2
- 230000000694 effects Effects 0.000 abstract description 4
- 239000011159 matrix material Substances 0.000 description 6
- 230000014509 gene expression Effects 0.000 description 3
- 238000010521 absorption reaction Methods 0.000 description 2
- 238000013459 approach Methods 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 230000008595 infiltration Effects 0.000 description 1
- 238000001764 infiltration Methods 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 238000005192 partition Methods 0.000 description 1
- 230000000737 periodic effect Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 239000002699 waste material Substances 0.000 description 1
Images
Landscapes
- Building Environments (AREA)
- Soundproofing, Sound Blocking, And Sound Damping (AREA)
Abstract
本发明公开了一种利用有源控制实现虚拟耦合空间的设计方法。由该种方法设计的虚拟耦合空间***,可以在一个普通的矩形空间中,实现某种目标耦合空间中的声场分布。该方法能用较小的成本实现耦合空间声场,并且有较好的效果。
Description
一.技术领域
本发明涉及一种基于有源控制的虚拟耦合空间的设计方法。
二.背景技术
可变室内声学***自问世已经有超过50年的历史。比较有代表性的有受援共振(AR)***(Parkin P H and Morgan K.″Assisted resonance″in the royal festival hall,london.J.Sound Vib.,1965;2(1):74-85.),混响室立体混响***(Harold B-M and Vincent M.Sound control techniquesfor the legitimate theatre and opera.Journal of the Audio Engineering Society,1961;9(3):184-186.),多通道混响(MCR)***(Franssen N V.Sur i′amplification des champs acoustiques.Acustica,1968;20:315-323.),声学控制***(ACS),可变室内声学***(VRAS;Poletti,M.,andEllison,S.Control of early and late energy in rooms with the variable room acoustics system.J.Acoust.Soc.Am.114,(2003).2341.),扩展多通道(EMCR)***。而迄今为止国内外尚未查到相关技术的发明专利。
以上这些***都仅能够改变空间内的混响时间,并不能准确地模拟出一个复杂结构的声学空间中所特有的声场,比如耦合空间。由于耦合空间能提供一种兼顾清晰度和混响感的声场,所以现代的厅堂设计中越来越多地引入了耦合空间的元素。但是引入的耦合空间利用率很低甚至完全不可利用,带来了较严重的空间资源浪费问题。
本发明提出的可变室内声学***,能够比较准确地在一个普通的空间中近似实现出耦合空间中所具有的声场分布,既能够在设计时节省空间,又能够获得较好品质的声场。
三.发明内容
本发明的目的是提供一种可以在普通的矩形空间中,产生接近于耦合空间中的声场的***。该方法先通过测量或者求解目标耦合空间中的声场来获得有源控制的目标量pc,下标c意为耦合(coupled)。然后利用有源控制策略,在普通的矩形空间中,布放控制声源和误差传声器,通过优化控制声源的强度,使得矩形空间中原声源、本发明所引入的控制声源产生的总声压逼近于原声源在目标耦合空间中产生的声压pc。
本发明是通过以下技术方案来实现的:
a)构建目标耦合空间的模型,其中包含声源的子空间的结构与控制的对象空间相同,利用模态展开的方法联合两个子空间在耦合开口上声压相等的条件,利用变分计算法,求解获得主子空间中的声压表达(即pc),
其中为耦合开口上的质点速度分布,通过改变其中的Gb(与图1中空间Vb的结构直接相关)可达到改变耦合子空间结构的目的,以产生不同耦合声场的效果;ρ是空气密度,c是空气中的声波速度,k为波数,Sω(r)是分布式声源的强度,Sa和Sb分别为两个子空间中的内表面积。为含源子空间(图1中的Va空间)未与无源子空间(图1中的Vb空间)耦合时其中的声压。
Λn=LxLyLz (4)
其中,
nx,ny,nz为偶数时,公式(8)中Lx,Ly和Lz的表达中分别取加号;nx,ny,nz为奇数时,分别取负号。βx0和βxl为x方向两平行平面(分别对应x=0平面和x=l平面)上的比声导纳,βy0和βyl为y方向两平行平面上的比声导纳,βz0和βzl为z方向两平行平面上的比声导纳。
b)在一个与目标耦合空间中的主子空间相同结构的普通矩形空间中,布放控制声源和误差传声器。M×N个控制声源以M行N列的形式布置于一个与目标耦合空间的耦合开口平行且大小相同,距离λ/4的平面内;M×N个传声器以M行N列的形式布置于一个与目标耦合空间的耦合开口平行且大小相同,距离λ/8的平面内,其中的λ为所控制的声压在最高频率的波长。
控制的目标为该总声压逼近目标耦合空间中相同位置的声压pc,代价函数为,
其中, c=(Zpqp-pc)*(Zpqp-pc);pp为初级声源产生声压的矩阵表达,ps为控制声源产生声压的矩阵表达,pc为当把初级声源至于虚拟耦合空间中时产生声压的矩阵表达。qp为初级声源的强度矩阵,qs为控制声源的强度矩阵。Zp为初级声源至误差传声器的传递函数,Zs为控制声源到误差传声器的传递函数,它们有如下表达,
qs=-b/a (12)
本发明的特点是:通过理论和数值求解目标耦合空间声场的声压,或者通过在现有耦合空间中直接测量的方法获得需要布放误差传声器处的声压,再利用有源控制的方法,能够较准确地在普通空间中实现耦合空间中的声场分布。本发明相比于以前的混响控制***,目标更明确,实现准确,且可通过改变控制参数容易地实现不同结构目标耦合空间的声场。
四.附图说明
图1目标耦合空间的示意图,Va为有源子空间,Vb为无源子空间,S0为耦合开口
图2普通矩形空间中有源控制***的布放图,该空间与图1中的Va尺度相同,取高度z=1.5m的一个xy剖面,re为误差传声器布放面,rs为控制声源布放面,C1为y=2,z=1.5的观察线,C2为x=4,z=1.5的观察线。
图3沿着如图2中的观察线C1方向的控制效果预测图,横轴为x∈[0,5]m,纵轴为C1线上的相对声压级(SPL)差。
五.具体实施方式
下面结合实施例参照附图对本发明进行详细说明:
本发明实现的目标是在普通矩形空间产生为如图1所示耦合空间中Va中的声场。
首先,按照公式(1-8)求得目标声压pc。本例中,Va为主空间,声源位于其中,其3维尺寸为5×4×3m3,其内壁比声导纳设定为βa=0.07663+0.06425i,换算成吸声系数约为0.4。Vb为次空间,其内无声源,一般用于调节Va中的声学特性,比如声场分布,混响时间长短等。Vb的3维尺寸为7×4×3m3,其内壁比声导纳设定为βb=0.0322-0.0383i,换算成吸声系数约为0.2。两个子空间中间的分隔上,在正中央开有面积为1×1m2的耦合开口S0,初级声源位置固定在(1.0,1.0,1.0)m处。
求得目标耦合空间中的声压分布后,取若干点布置误差传声器,本例中选择如图2中re平面上的9点,re平面与目标耦合空间中耦合开口大小相同,平行,且距离λ/8。控制声源分布于rs平面,距离耦合开口λ/4。其中的λ为所控制的声压在最高频率的波长。图2显示为高度z=1.5m的xy平面,S0为虚拟目标耦合空间中耦合开口所在的位置,C1为沿x方向穿越空间正中的观察线。利用公式(12),求得控制声源的强度矩阵。
按照如上条件,数值仿真获得的控制效果预测图如图3所示。曲线为250Hz频率下,声压级沿着图2中C1曲线的分布图,横轴为沿着x方向的C1线,纵轴为标准声压级差,方块线为普通矩形空间中C1线上的相对声压分布,三角线为目标耦合空间中C1线上的相对声压分布,空心圆线为在普通矩形空间中使用有源控制后的相对声压分布。可以看到耦合空间的存在较大地抑制了空间中某些模态(或者叫振型)过于强烈而导致的声压分布周期性起伏,对均匀化声压空间分布有帮助。故在普通空间中引入有源控制***,构造虚拟耦合空间,可以有效地控制矩形空间中的声场,逼近存在耦合空间时的声压分布。
Claims (2)
1.一种虚拟耦合空间***的设计方法,其特征是通过在普通的矩形空间中引入有源控制***来实现一种新的可变室内声学***,使普通的矩形空间中产生类似于耦合空间中所具有的声场,并且可通过改变控制参数模拟不同结构的耦合空间,包括以下步骤:
(1)使用矩形空间中的格林函数和两空间的耦合条件,求解获得目标耦合空间中的声场的声压,或者通过测量获得一个已有耦合空间中的声场的声压,并以此声压作为有源控制的目标;
(2)使用闭空间中的有源控制策略,布置控制声源和误差传声器,优化控制声源强度,具体的的有源控制策略如下:
1)M×N个控制声源以M行N列的形式布置于一个与目标耦合空间的耦合开口平行且大小相同的平面内,两者相距λ/4,λ为所要控制的声压的最高频率对应的声波波长;
2)M×N个传声器以M行N列的形式布置于一个与目标耦合空间的耦合开口平行且大小相同的平面内,两者相距λ/8;
3)以原矩形空间声场的声压和控制声源产生的总声压与如权利要求2所述的pc之差的平方值作为代价函数对次级源强度进行优化,即令J=|pp+ps-pc|2,其中pp为原矩形空间中的原有声源产生的声压,ps为该矩形空间中控制声源产生的声压。
2.如权利要求1所述的虚拟耦合空间***的设计方法,其特征在于,确定***控制所要达到的目标,即目标耦合空间中的声场的声压,可通过理论模型预测获得,也可通过在一个实际的耦合空间中测量获得,设为pc,下标c意为耦合(coupled)。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2010102168282A CN101902680B (zh) | 2010-07-05 | 2010-07-05 | 基于有源控制的虚拟耦合空间的设计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2010102168282A CN101902680B (zh) | 2010-07-05 | 2010-07-05 | 基于有源控制的虚拟耦合空间的设计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101902680A CN101902680A (zh) | 2010-12-01 |
CN101902680B true CN101902680B (zh) | 2011-11-16 |
Family
ID=43227808
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2010102168282A Expired - Fee Related CN101902680B (zh) | 2010-07-05 | 2010-07-05 | 基于有源控制的虚拟耦合空间的设计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101902680B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104103267A (zh) * | 2014-06-19 | 2014-10-15 | 广东电网公司汕头供电局 | 用于变压器低频噪声的虚拟声屏障 |
CN104778321B (zh) * | 2015-04-14 | 2017-11-17 | 哈尔滨工程大学 | 一种复杂形状声学空间的分析预报方法 |
CN110111765B (zh) * | 2019-05-21 | 2022-06-14 | 东南大学 | 一种一维声场条件下的反射声有源控制方法 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2299735B1 (en) * | 2000-07-19 | 2014-04-23 | Koninklijke Philips N.V. | Multi-channel stereo-converter for deriving a stereo surround and/or audio center signal |
JP4658968B2 (ja) * | 2004-01-19 | 2011-03-23 | コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ | 広い領域にわたってステレオサウンド感を提供する点サウンド生成手段及び空間サウンド生成手段を有する装置 |
DE102005016021A1 (de) * | 2005-04-07 | 2006-10-12 | Airbus Deutschland Gmbh | Aktives Gegenschallsystem mit spezieller Anordnung der Sekundäraktuatoren zur Reduzierung eines Schalldurchgangs an einer offenen Grenzfläche zweier Volumina, aktive Gegenschallsystem-Anordnung, Verfahren zur aktiven Schallwellenreduzierung und Verwendung eines aktiven Gegenschallsystems zur aktiven Schallwellenreduzierung für ein zumindest teilweise geöffnetes Kontrollvolumen |
-
2010
- 2010-07-05 CN CN2010102168282A patent/CN101902680B/zh not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
CN101902680A (zh) | 2010-12-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Yu et al. | On the sound insulation of acoustic metasurface using a sub-structuring approach | |
CN105427853A (zh) | 宽频带微穿孔板吸声体及其性能预测方法、结构设计方法 | |
CN101902680B (zh) | 基于有源控制的虚拟耦合空间的设计方法 | |
CN107944108A (zh) | 基于统计能量分析的船舶舱室噪声预报方法 | |
CN106096184A (zh) | 一种基于声强法的变压器多点声源模型的噪声计算方法及*** | |
CN109307711B (zh) | 基于混响室及消声室的列车贯通道隔声性能测试方法 | |
CN110543669A (zh) | 一种声学超材料板的隔声仿真计算方法 | |
CN109635486B (zh) | 一种高分辨率三维成像声纳换能器阵列稀疏优化方法 | |
GB2473968A (en) | Loudspeaker array configuration with 3D plot/display | |
Greschner et al. | Measurements and wall modeled LES simulation of trailing edge noise caused by a turbulent boundary layer | |
Zhang et al. | Secondary source and error sensing strategies for the active control of sound transmission through a small opening | |
CN109543227B (zh) | 大型反射面天线的风致指向误差实时测量方法、雷达天线 | |
CN106153176A (zh) | 一种用于对单相油浸式变压器噪声进行预估的***及方法 | |
CN104008287B (zh) | 基于pso‑mvdr的声场重构与鬼影抑制方法 | |
CN111638272B (zh) | 一种在役水工混凝土结构内损监测预警***及其方法 | |
Chen et al. | Study on the Optimization of the Distribution of Absorbing Material on a Noise Barrier | |
Zou et al. | Applying an active noise barrier on a 110 KV power transformer in Hunan | |
Schoenwald | Flanking sound transmission through lightweight framed double leaf walls: Prediction using statistical energy analysis | |
CN102589676B (zh) | 一种基于室内空间剖分的声线追踪方法 | |
Zhang et al. | Numerical investigation on the transmission loss of skin panels based on the intelligent pso-cga algorithm | |
CN104778321B (zh) | 一种复杂形状声学空间的分析预报方法 | |
CN103742892A (zh) | 用于余热锅炉的共振式***的设计方法、结构及其应用 | |
SheikhMozafari | Enhancing Sound Absorption in Micro-Perforated Panel and Porous Material Composite in Low Frequencies: A Numerical Study Using FEM | |
Oyelade et al. | Sound transmission through triple plates separated by air cavities in the low-frequency range | |
CN110866333A (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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20111116 Termination date: 20150705 |
|
EXPY | Termination of patent right or utility model |