CN101902680B - Design method of virtual coupled space based on active control - Google Patents

Design method of virtual coupled space based on active control Download PDF

Info

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
acoustic pressure
control
coupled
sound
sound source
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
Application number
CN2010102168282A
Other languages
Chinese (zh)
Other versions
CN101902680A (en
Inventor
浦宏杰
邱小军
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Nanjing University
Original Assignee
Nanjing University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Nanjing University filed Critical Nanjing University
Priority to CN2010102168282A priority Critical patent/CN101902680B/en
Publication of CN101902680A publication Critical patent/CN101902680A/en
Application granted granted Critical
Publication of CN101902680B publication Critical patent/CN101902680B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Soundproofing, Sound Blocking, And Sound Damping (AREA)
  • Building Environments (AREA)

Abstract

The invention discloses a design method for realizing a virtual coupled space based on active control. The virtual coupled space system designed by using the method can, in a common rectangular space, realize sound field distribution in a certain target coupled space. The method can realize sound field in the coupled space at small cost, and achieves good effect.

Description

Method for designing based on the virtual coupled space of active control
One. technical field
The present invention relates to a kind of method for designing of the virtual coupled space based on active control.
Two. background technology
Sound system has had the history above 50 years in the variable cells from coming out.More representational resonance (AR) system of receiving aid (Parkin P H and Morgan K. " Assisted resonance " in the royal festival hall, london.J.Sound Vib., 1965; 2 (1): 74-85.), reverberation room ambiophonic system (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.), multichannel reverberation (MCR) system (Franssen N V.Sur i ' amplification des champs acoustiques.Acustica, 1968; 20:315-323.), sound system (VRAS in the acoustics control system (ACS), variable cells; Poletti, M., andEllison, S.Control of early and late energy in rooms with the variable room acoustics system.J.Acoust.Soc.Am.114, (2003) expansion multichannel (EMCR) system .2341.).And the domestic and international up to now patent of invention of not finding correlation technique as yet.
More than these systems all only can change reverberation time in the space, can not simulate peculiar sound field in the acoustic space of a labyrinth exactly, such as coupled room.Because coupled room can provide a kind of sound field of taking into account definition and reverberation sense, so introduced the element of coupled room in the design of the modern hall more and more.But the coupled room utilance of introducing is very low even can not utilize fully, has brought more serious space resources waste problem.
Sound system in the variable cells that the present invention proposes, can be more exactly in a common space the approximate sound field that realizes out in the coupled room being had distribute, can either can obtain the sound field of better quality again at design time save space.
Three. summary of the invention
The purpose of this invention is to provide a kind of can be in common coffin, produce the system that approaches the sound field in the coupled room.This method obtains the aim parameter p of active control by measuring or find the solution sound field in the target coupled room earlier c, subscript c means coupling (coupled).Utilize active control strategy then, in common coffin, lay control sound source and error microphone, by the intensity of optimal control sound source, the total acoustic pressure that makes acoustic source in the coffin, the control sound source that the present invention introduced produce approaches the acoustic pressure p that acoustic source produces in the target coupled room c
The present invention is achieved through the following technical solutions:
A) model of establishing target coupled room, the structure of subspace that wherein comprises sound source is identical with the object space of control, the method of utilizing mode to launch is united the two sub spaces condition that acoustic pressure equates on coupling opening, utilizes the Variational Calculation method, and finding the solution the acoustic pressure expression that obtains in the principal subspace (is p c),
p c ( t ) = - jkρc ∫ ∫ ∫ S ω / G a ( r | r 0 ) S ω ( r 0 ) dr 0 + jkρc ∫ ∫ S 0 G a ( r | r 0 ) u 0 ( r 0 ) dr 0 - - - ( 1 )
Wherein
Figure BSA00000167955400012
For the particle velocity on the coupling opening distributes, by changing G wherein b(with space V among Fig. 1 bStructure directly related) can reach the purpose that changes coupling subspace structure, to produce the effect of different coupling sound fields; ρ is an atmospheric density, and c is airborne acoustic velocity, and k is a wave number, S ω(r) be the intensity of distributed sound source, S aAnd S bBe respectively the internal surface area in two sub spaces.
Figure BSA00000167955400021
For containing subspace, the source (V among Fig. 1 aThe space) not with the passive subspace (V among Fig. 1 bThe space) when coupling acoustic pressure wherein.
G aAnd G bBe respectively the Green's function in two sub spaces, have expression formula respectively
Figure BSA00000167955400022
Figure BSA00000167955400023
Eigenfunction φ n, eigenvalue η nAverage Λ with eigenfunction nBe respectively,
Figure BSA00000167955400024
× cos ( q zn z πz l z + j β z 0 kl z πq zn z )
η n 2 = ( πq xn x l x ) 2 + ( πq yn y l y ) 2 + ( πq zn z l z ) 2 - - - ( 3 )
Λ n=L xL yL z (4)
Wherein,
Figure BSA00000167955400028
Figure BSA00000167955400029
Figure BSA000001679554000211
Figure BSA000001679554000212
L x = 1 2 ± 1 2 πq xn x sin ( πq xn x ) , L y = 1 2 ± 1 2 πq yn y sin ( πq yn y ) , L z = 1 2 ± 1 2 π q zn z sin ( πq zn z ) - - - ( 8 )
n x, n y, n zDuring for even number, L in the formula (8) x, L yAnd L zExpression in get plus sige respectively; n x, n y, n zDuring for odd number, get negative sign respectively.β X0And β XlBe the ratio acoustic admittance on x direction two parallel planes (corresponding x=0 plane and x=l plane respectively), β Y0And β YlBe the ratio acoustic admittance on y direction two parallel planes, β Z0And β ZlBe the ratio acoustic admittance on z direction two parallel planes.
B) one with the target coupled room in the ordinary rectangular space of principal subspace same structure in, lay control sound source and error microphone.M * N control sound source is parallel and big or small identical in the coupling opening with the target coupled room with the arranged in form of the capable N row of M, in the plane of distance lambda/4; M * N microphone is parallel and big or small identical in the coupling opening with the target coupled room with the arranged in form of the capable N row of M, and in the plane of distance lambda/8, λ wherein is the acoustic pressure the controlled wavelength at highest frequency.
The target of control is approached the acoustic pressure p of same position in the target coupled room for this total acoustic pressure c, cost function is,
J = | p p + p s - p c | 2
= q s * aq s + q s * b + q s b * + c - - - ( 9 )
Wherein,
Figure BSA00000167955400034
C=(Z pq p-p c) *(Z pq p-p c); p pBe the expression matrix of elementary sound source generation acoustic pressure, p sBe the expression matrix of control sound source generation acoustic pressure, p cFor when elementary sound source as for virtual coupled space in the time produce the expression matrix of acoustic pressure.q pBe the intensity matrix of elementary sound source, q sIntensity matrix for the control sound source.Z pBe the transfer function of elementary sound source to error microphone, Z sBe the transfer function of control sound source to error microphone, they have following expression,
Figure BSA00000167955400035
Figure BSA00000167955400036
Wherein
Figure BSA00000167955400037
Through optimizing, obtaining the optimal control strength of sound source be
q s=-b/a (12)
Characteristics of the present invention are: by the acoustic pressure of theoretical and numerical solution target coupled room sound field, perhaps obtain to lay the acoustic pressure at error microphone place by the method for in existing coupled room, directly measuring, utilize the method for active control again, can realize in general spatial more exactly that the sound field in the coupled room distributes.The present invention than before the reverberation control system, target is clearer and more definite, realizes accurately, and can be by changing the sound field that Control Parameter easily realize different structure target coupled room.
Four. description of drawings
The schematic diagram of Fig. 1 target coupled room, V aBe active subspace, V bBe passive subspace, S 0Be coupling opening
Active guidance system lays figure, the V among this space and Fig. 1 in Fig. 2 ordinary rectangular space aYardstick is identical, gets the xy section of height z=1.5m, r eFor error microphone lays face, r sFor the control sound source lays face, C1 is y=2, and the observation line of z=1.5, C2 are x=4, the observation line of z=1.5.
Fig. 3 is along the control effect prediction figure as the observation line C1 direction among Fig. 2, and transverse axis is x ∈ [0,5] m, and the longitudinal axis is that the relative sound pressure level (SPL) on the C1 line is poor.
Five. embodiment
The present invention is described in detail with reference to accompanying drawing below in conjunction with embodiment:
The target that the present invention realizes is to be produced as shown in Figure 1 V in the coupled room in the ordinary rectangular space aIn sound field.
At first, try to achieve target acoustic pressure p according to formula (1-8) cIn this example, V aBe the principal space, sound source is positioned at wherein, and its 3 dimension is of a size of 5 * 4 * 3m 3, its inwall is set at β than acoustic admittance a=0.07663+0.06425i is converted into acoustic absorptivity and is about 0.4.V bBe inferior space, no sound source in it generally is used to regulate V aIn acoustic characteristic, distribute reverberation time length etc. such as sound field.V b3 dimensions be of a size of 7 * 4 * 3m 3, its inwall is set at β than acoustic admittance b=0.0322-0.0383i is converted into acoustic absorptivity and is about 0.2.In the separation in the middle of two sub spaces, having area in centre is 1 * 1m 2Coupling opening S 0, elementary sound source position is fixed on (1.0,1.0,1.0) m place.
After trying to achieve the acoustic pressure distribution in the target coupled room, get some somes placement error microphones, select as r among Fig. 2 in this example eOn the plane 9 point, r eThe plane is identical, parallel with coupling opening size in the target coupled room, and distance lambda/8.Guide sound source distribution is in r sThe plane is apart from coupling opening λ/4.λ wherein is the acoustic pressure the controlled wavelength at highest frequency.Fig. 2 is shown as the xy plane of height z=1.5m, S 0Be the position at coupling opening place in the virtual target coupled room, C1 is for passing through the observation line of center, space along the x direction.Utilize formula (12), try to achieve the intensity matrix of control sound source.
According to as above condition, the control effect prediction figure that numerical simulation obtains as shown in Figure 3.Curve is under the 250Hz frequency, the distribution map of sound pressure level C1 curve in Fig. 2, transverse axis is the C1 line along the x direction, the longitudinal axis is that standard sound pressure is differential, the square line is that the relative acoustic pressure on the C1 line distributes in the ordinary rectangular space, triangle line is that the relative acoustic pressure on the C1 line distributes in the target coupled room, and the open circles line is that the relative acoustic pressure after the active control of use distributes in the ordinary rectangular space.The existence that can see coupled room has suppressed the acoustic pressure distribution periodic undulations that some mode (perhaps being the vibration shape) causes too strongly in the space significantly, and is helpful to homogenizing acoustic pressure spatial distribution.So introduce active guidance system in general spatial, the constructing virtual coupled room can be controlled the sound field in the coffin effectively, and the acoustic pressure of approaching when having coupled room distributes.

Claims (2)

1. the method for designing of a virtual coupled space system, it is characterized in that realizing sound system in a kind of new variable cells by in common coffin, introducing active guidance system, make in the common coffin and to produce the sound field that is similar in the coupled room to be had, and can may further comprise the steps by changing the coupled room of Control Parameter simulation different structure:
(1) Green's function in the use coffin and the coupling condition in two spaces are found the solution the acoustic pressure that obtains the sound field in the target coupled room, perhaps obtain the acoustic pressure of the sound field in the existing coupled room by measurement, and with the target of this acoustic pressure as active control;
(2) use the active control strategy close in the space, arrange control sound source and error microphone, the optimal control strength of sound source, concrete active control strategy is as follows:
1) in coupling opening with the target coupled room in the parallel and big or small identical plane, both are at a distance of λ/4 with the arranged in form of the capable N of M row for M * N control sound source, and λ is the wave length of sound of the highest frequency correspondence of the acoustic pressure that will control;
2) in coupling opening with the target coupled room in the parallel and big or small identical plane, both are λ/8 apart with the arranged in form of the capable N row of M for M * N microphone;
3) with the acoustic pressure of former coffin sound field and the total acoustic pressure and the p as claimed in claim 2 of control sound source generation cThe square value of difference as cost function secondary source intensity is optimized, even J=|p p+ p s-p c| 2, p wherein pBe the acoustic pressure that the original sound source in the former coffin produces, p sAcoustic pressure for control sound source generation in this coffin.
2. the method for designing of virtual coupled space system as claimed in claim 1, it is characterized in that, determine the target that system's control will reach, it is the acoustic pressure of the sound field in the target coupled room, can obtain by the theoretical model prediction, also can obtain, be made as p by in the coupled room of a reality, measuring c, subscript c means coupling (coupled).
CN2010102168282A 2010-07-05 2010-07-05 Design method of virtual coupled space based on active control Expired - Fee Related CN101902680B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN2010102168282A CN101902680B (en) 2010-07-05 2010-07-05 Design method of virtual coupled space based on active control

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN2010102168282A CN101902680B (en) 2010-07-05 2010-07-05 Design method of virtual coupled space based on active control

Publications (2)

Publication Number Publication Date
CN101902680A CN101902680A (en) 2010-12-01
CN101902680B true CN101902680B (en) 2011-11-16

Family

ID=43227808

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2010102168282A Expired - Fee Related CN101902680B (en) 2010-07-05 2010-07-05 Design method of virtual coupled space based on active control

Country Status (1)

Country Link
CN (1) CN101902680B (en)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104103267A (en) * 2014-06-19 2014-10-15 广东电网公司汕头供电局 Virtual sound barrier applied to low frequency noise of transformer
CN104778321B (en) * 2015-04-14 2017-11-17 哈尔滨工程大学 A kind of analysis and prediction method of complicated shape acoustic space
CN110111765B (en) * 2019-05-21 2022-06-14 东南大学 Reflected sound active control method under one-dimensional sound field condition

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
ES2461167T3 (en) * 2000-07-19 2014-05-19 Koninklijke Philips N.V. Multi-channel stereo converter to derive a stereo surround signal and / or audio center
CN1910958B (en) * 2004-01-19 2011-03-30 皇家飞利浦电子股份有限公司 Device having a point and a spatial sound generating-means for providing stereo sound sensation over a large area
DE102005016021A1 (en) * 2005-04-07 2006-10-12 Airbus Deutschland Gmbh Active counter sound system has secondary actuators arranged at edge of open depressing surface to be penetrable by primary sound wave

Also Published As

Publication number Publication date
CN101902680A (en) 2010-12-01

Similar Documents

Publication Publication Date Title
CN104464739A (en) Audio signal processing method and device and difference beam forming method and device
CN101902680B (en) Design method of virtual coupled space based on active control
CN110543669A (en) Sound insulation simulation calculation method of acoustic metamaterial plate
CN203136135U (en) High-pitch horn implanted with phase plug for sound amplification
Abhayapala et al. Spherical harmonic analysis of wavefields using multiple circular sensor arrays
CN103616071B (en) Patch near field acoustic holography-sound quality objective parameter distributed in three dimensions method for visualizing
EP3847719A1 (en) Open cavity system for directed amplification of acoustic signals
CN104733850A (en) Holographic modulation based artificial tensor impedance surfaced antenna and implementation method thereof
CN104323794A (en) Split-type rank addressing method for three-dimensional ultrasonic imaging
CN109100685A (en) A kind of passive acoustic direction blending algorithm of two-sided quaternary cross battle array
Gao et al. Broadband ventilated sound insulation in a highly sparse acoustic meta-insulator array
CN110095755B (en) Sound source positioning method
CN104023304A (en) Method for simplifying five-loudspeaker system to four-loudspeaker system
CN101718869A (en) Planar acoustic array with wideband and super-wide covering property
CN106454648B (en) A kind of acoustic waveguide
CN102984642A (en) Three-dimensional translation method for five loudspeakers
CN105098370B (en) The method that low-frequency sound wave auto-collimation transmits extremely is realized using sound metamaterial
CN110736959B (en) Planar co-prime array design method based on sum-difference cooperative array construction
WO2022104529A1 (en) Generating method and generating system for acoustic tweezers
Xing et al. A fusion algorithm of passive sound source localization based on the two-plane four-element cross array
CN112887874A (en) Strong sound system based on full-band sound wave control
CN1110701C (en) Method for processing sonar wave beams and its data processing system
Reeder et al. Short-range acoustic propagation variability on a shelf area with strong nonlinear internal waves
CN110334480A (en) Curve design method is extended for reducing the double offset antenna minor face of noise temperature
CN104363548A (en) Piezoelectric network method for improving sound insulating performance of plate and shell structure

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