EP1502475A1 - Procede et systeme de representation d un champ acoustique - Google Patents

Procede et systeme de representation d un champ acoustique

Info

Publication number
EP1502475A1
EP1502475A1 EP03749929A EP03749929A EP1502475A1 EP 1502475 A1 EP1502475 A1 EP 1502475A1 EP 03749929 A EP03749929 A EP 03749929A EP 03749929 A EP03749929 A EP 03749929A EP 1502475 A1 EP1502475 A1 EP 1502475A1
Authority
EP
European Patent Office
Prior art keywords
representative
parameters
acquisition means
sensors
determining
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.)
Granted
Application number
EP03749929A
Other languages
German (de)
English (en)
Other versions
EP1502475B8 (fr
EP1502475B1 (fr
Inventor
Rémy BRUNO
Arnaud Laborie
Sébastian Montoya
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.)
Trinnov Audio
Original Assignee
Individual
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 Individual filed Critical Individual
Publication of EP1502475A1 publication Critical patent/EP1502475A1/fr
Application granted granted Critical
Publication of EP1502475B1 publication Critical patent/EP1502475B1/fr
Publication of EP1502475B8 publication Critical patent/EP1502475B8/fr
Anticipated expiration legal-status Critical
Expired - Lifetime legal-status Critical Current

Links

Classifications

    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04R—LOUDSPEAKERS, MICROPHONES, GRAMOPHONE PICK-UPS OR LIKE ACOUSTIC ELECTROMECHANICAL TRANSDUCERS; DEAF-AID SETS; PUBLIC ADDRESS SYSTEMS
    • H04R3/00—Circuits for transducers, loudspeakers or microphones
    • H04R3/005—Circuits for transducers, loudspeakers or microphones for combining the signals of two or more microphones
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04R—LOUDSPEAKERS, MICROPHONES, GRAMOPHONE PICK-UPS OR LIKE ACOUSTIC ELECTROMECHANICAL TRANSDUCERS; DEAF-AID SETS; PUBLIC ADDRESS SYSTEMS
    • H04R3/00—Circuits for transducers, loudspeakers or microphones
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04R—LOUDSPEAKERS, MICROPHONES, GRAMOPHONE PICK-UPS OR LIKE ACOUSTIC ELECTROMECHANICAL TRANSDUCERS; DEAF-AID SETS; PUBLIC ADDRESS SYSTEMS
    • H04R3/00—Circuits for transducers, loudspeakers or microphones
    • H04R3/04—Circuits for transducers, loudspeakers or microphones for correcting frequency response

Definitions

  • the present invention relates to a method and a device for representing an acoustic field from signals delivered by acquisition means.
  • the methods and systems for acquiring and representing existing sound environments use models based on physically impracticable acquisition means, in particular as regards the electro-acoustic and / or structural characteristics of these acquisition means.
  • the acquisition means are, for example, made up of a set of elementary measurement elements or sensors arranged in specific places in space and having intrinsic electro-acoustic acquisition characteristics.
  • the systems grouped under the term "Ambisonic" only consider the directions of origin of the sounds relative to the center of the acquisition means formed by a plurality of elementary sensors, which leads to assimilate the acquisition means to a point microphone.
  • these systems represent the sound environment by modeling virtual sources whose angular distribution around the center theoretically makes it possible to obtain such a sound environment.
  • the object of the invention is to solve this problem by providing a method and a device delivering a representation of the acoustic field substantially independent of the characteristics of the acquisition means.
  • the subject of the present invention is a method of representing an acoustic field comprising a step of acquiring measurement signals delivered by acquisition means formed by one or more elementary sensors exposed to said acoustic field, characterized in that 'it comprises :
  • Said structural characteristics include at least position characteristics of said elementary sensors with respect to a predetermined reference point of said acquisition means;
  • - Said encoding filters are also representative of electro-acoustic characteristics of said acquisition means; - Said electro-acoustic characteristics include at least characteristics linked to the electro-acoustic capacities for intrinsic acquisition of said elementary sensors;
  • said step of determining the encoding filters comprises: a sub-step of determining a sampling matrix representative of the acquisition capacities of said acquisition means; a substep for determining an intercorrelation matrix representative of the resemblance between said measurement signals delivered by the elementary sensors forming said acquisition means; and a sub-step of determining an encoding matrix from said sampling matrix, from said intercorrelation matrix, and from a parameter representative of a desired compromise between the fidelity of representation of the acoustic field. and minimizing the background noise induced by the acquisition means, which matrix is representative of said encoding filters;
  • said step of determining the sampling matrix is carried out, for each of said elementary sensors forming said acquisition means, from:
  • step of determining the sampling matrix is further carried out on the basis of at least one of the parameters from:
  • parameters representative of the orientations of all or part of the sensors namely of their direction of maximum sensitivity;
  • said calibration step comprises, for at least one of said elementary sensors forming said acquisition means:
  • said calibration step further comprises:
  • a substep for transmitting a specific acoustic field to said at least one sensor said acquisition substep corresponding to the acquisition of the signals delivered by this sensor when exposed to said specific acoustic field;
  • said calibration step comprises a substep for receiving a finite number of signals representative of the electro-acoustic and structural characteristics of said sensors forming said acquisition means, which signals are directly used during said substep for determining the electro-acoustic and / or structural characteristics of said acquisition means;
  • the invention also relates to a computer program comprising program code instructions for executing the steps of the method, as described above when said program is executed on a computer.
  • the invention also relates to a mobile support of the type comprising at least one processing processor and a non-volatile memory element, characterized in that said memory comprises a program comprising code instructions for the execution of the steps of the method as described previously when said processor executes said program.
  • the invention also relates to a device for representing an acoustic field connectable to acquisition means formed by one or more elementary sensors delivering measurement signals when they are exposed to said acoustic field, characterized in that it comprises a module for processing the measurement signals by applying encoding filters representative of at least structural characteristics of said acquisition means to these measurement signals to deliver a signal which comprises a finite number of coefficients representative over time and in the three dimensions of the space of said acoustic field, said coefficients making it possible to obtain a representation of said acoustic field substantially independent of the characteristics of said acquisition means,
  • - Said encoding filters are also representative of electro-acoustic characteristics of said acquisition means
  • said means for determining encoding filters receive as input at least one of the following parameters: - parameters representative of the positions relative to the center of said means for acquiring all or part of the sensors;
  • said means for determining the parameters received by said means for determining the encoding filters comprising at least one of the following elements:
  • - Fig.1 is a representation of a spherical coordinate system
  • - Fig.2 is a representative diagram of the acquisition means used
  • - Fig.3 is a general flowchart of the method of the invention: - Fig.4 is a flow diagram of the detail of an embodiment of the calibration step of the method of the invention;
  • - Fig.5 is a flow diagram of the detail of an embodiment of the step of determining the encoding filters of the method of the invention
  • - Fig.6 is a detail diagram of an embodiment of the step of applying the encoding filters
  • - Fig.7 is a block diagram of a device suitable for implementing the method of the invention.
  • Figure 1 there is shown a conventional spherical coordinate system, so as to specify the coordinate system to which reference is made in the text.
  • This coordinate system is an orthonormal coordinate system, of O origin and comprising three axes (QX), (OY) and (02).
  • a position denoted x is described by means of its spherical coordinates (r, ⁇ ,), where r denotes the distance from the origin O, ⁇ orientation in the vertical plane and orientation in the horizontal plane.
  • an acoustic field is known if we define at any point at each instant t the acoustic pressure noted p (r, ⁇ , ⁇ , t), whose Fourier transform is noted P (r, ⁇ , ⁇ J) where / designates the frequency.
  • the method of the invention is based on the use of spatiotemporal functions making it possible to describe any sound field in time and in the three dimensions of space.
  • these functions are so-called spherical Fourier-Bessel functions of the first kind, hereinafter called Fourier-Bessel functions.
  • the Fourier-Bessel functions correspond to the solutions of the wave equation and constitute a base which generates all the acoustic fields produced by sources located outside this zone. Any three-dimensional acoustic field can therefore be expressed by a linear combination of the Fourier-Bessel functions, according to the expression of the inverse Fourier-Bessel transform which is expressed:
  • J v (x) is the Bessel function of pre- first species of order v
  • yf ( ⁇ , ⁇ ) is the real spherical harmonic of order / and of term m, with m going from - / to /, defined by:
  • the Fourier-Bessel coefficients are also expressed in the time domain by the coefficients p !> M (t) corresponding to the inverse temporal Fourier transform of the coefficients P ⁇ , m f).
  • the acoustic field is broken down on the basis of functions, where each of the functions is expressed by an optionally infinite linear combination of Fourier-Bessel functions.
  • each elementary sensor is arranged at specific points in space around a predetermined point 4 designated as the center of the acquisition means 1.
  • the position of each elementary sensor can be expressed in space in a spherical coordinate system such as that described with reference to FIG. 1, centered on the center 4 of the acquisition means 1.
  • each sensor 2 n of the acquisition means 1 When exposed to an acoustic field P, each sensor 2 n of the acquisition means 1 delivers a measurement signal c sauwhich corresponds to the measurement made by this sensor in the acoustic field P.
  • the acquisition means 1 deliver a plurality of signals Ci to c / v which are the signals for measuring the acoustic field P by the acquisition means 1.
  • These measurement signals ci to cw delivered by the acquisition means 1 are therefore directly linked to the acquisition capacities of the elementary sensors
  • FIG 3 there is shown a general flowchart of the method of the invention.
  • the method begins with a step 10 for entering parameters and a step 20 for calibrating the acquisition means, which make it possible to define a set of parameters representative of the structural and / or electro-acoustic characteristics of the acquisition means 1.
  • Certain parameters and in particular parameters representative of electro-acoustic characteristics are dependent on the frequency.
  • the input step 10 and the calibration step 20, which is described in more detail with reference to FIG. 4, can be carried out simultaneously or in any order.
  • the method of the invention may include only the input step 10.
  • the steps 10 of input and 20 of calibration make it possible to determine, for one or more sensors, all or part of the following parameters:
  • H n (f) representative of the frequency response of the corresponding sensor 2 n , for each frequency /; to the sensitivity of the sensor 2 n in the direction aJJ);
  • - parameters cf n (f) representative of the power spectral density of the background noise of the sensor 2 ⁇ ;
  • ⁇ (f) specifying a compromise between the fidelity of representation of the acoustic field P and the minimization of the background noise provided by the sensors 2 ⁇ to n and which can take all the values between 0 and 1:
  • all or part of the parameters described is considered to be independent of the frequency.
  • the parameters ⁇ (f), L (f) and ⁇ (k, m k ) ⁇ (f) are representative of the optimization strategies making it possible to control the extraction of spatiotemporal information from the acoustic field P from the measurement signals.
  • ci at C / v and are entered during the entry step 10.
  • the other parameters can be entered during the input step 10 or determined during the calibration step 20.
  • the method of the invention is carried out only with the parameters ⁇ (f), L (f) and the set of parameters x precedeor the set of parameters B n j > m (f) or a combination of parameters x propeland Bn j), so as to have at least one parameter per elementary sensor 2 n .
  • the method comprises a step 30 for determining encoding filters representative of the at least structural and advantageously electro-acoustic characteristics of the acquisition means 1.
  • This step 30, described in more detail with reference to FIG. 5, makes it possible to take into account all the parameters determined during the input steps 10 and / or calibration.
  • These encoding filters are therefore representative at least of the position characteristics of the elementary sensors 2 n relative to the reference point 4 of the acquisition means 1.
  • these filters are also representative of other structural characteristics of the acquisition means 1, such as the orientation of the elementary sensors 2 ⁇ to 2N OR their mutual influences, as well as their electro-acoustic acquisition capacities and in particular their noise of background, their directivity diagram, their frequency response, ...
  • the encoding filters obtained at the end of step 30 can be memorized, so that steps 10, 20 and 30 are only repeated in the event of modification of the acquisition means 1 or of the optimization strategies. .
  • These encoding filters are applied during a step 40 of processing the signals ci to c / v coming from the elementary sensors 2 ⁇ to 2 / v.
  • This processing corresponds to filtering of the signals and to combinations of the filtered signals.
  • this step 40 of processing the measurement signals by applying the encoding filters to these signals a finite number of coefficients representative in time and in the three dimensions of the space of the acoustic field P is issued. These coefficients are so-called Fourier-Bessel coefficients, denoted P ⁇ , m (f) and correspond to a representation of the acoustic field P substantially independent of the characteristics of the acquisition means 1.
  • the calibration step 20 makes it possible to directly determine the coefficients B ni f representative of the acquisition capacities of the acquisition means 1.
  • This step 20 begins with a sub-step 22 of emission of a specific acoustic field towards the acquisition means 1 and with a sub-step 24 of acquisition of measurement signals by the acquisition means 1 exposed to the sound field emitted.
  • step 20 of calibration is implemented using means of generating an acoustic field which comprise only a fixed speaker, assumed to be punctual and with a flat frequency response, the speaker and the acquisition means 1 being placed in an anechoic environment.
  • the loudspeaker emits the same sound field and the acquisition means 1 are placed in the same position but they are oriented in different and known directions.
  • the loudspeaker is in a different position (] ⁇ ⁇ , ⁇ q p , ⁇ q p ) for each field q generated.
  • the acquisition means 1 are thus exposed to an acoustic field q whose Fourier-Bessel coefficients P ⁇ , m , q (f), in the reference frame of the acquisition means 1, are known up to a given order, noted E 3 .
  • the measurement signals delivered following the acquisition sub-step 24 are a finite number of coefficients represented sensitive to the acoustic field q generated, as well as the acquisition capacities of the acquisition means 1.
  • the method includes a modeling sub-step 26 making it possible to determine a representation of the Q acoustic fields emitted during sub-step 22.
  • a modeling matrix P representative of all of the Q known fields to which the acquisition means 1 are successively exposed is determined.
  • This matrix P is a matrix of size (E 3 + l) 2 over Q consisting of the elements P ⁇ , m , q (f), the indices (l, m) designating the line l 2 + l + m and the index q designating column q.
  • the matrix P therefore has the following form:
  • the sound field produced by the loudspeaker is modeled by spherical radiation, thus, in the reference of the acquisition means 1, the coefficients P ⁇ , m , q (f) of each sound field q thus generated are known, thanks to the relation:
  • the coefficients obtained during sub-step 26 are then used during a sub-step 28 in order to determine parameters representative of the structural and / or acoustic characteristics of the acquisition means 1.
  • this sub-step 28 also uses the modeling matrix R determined during sub-step 26.
  • This sub-step 28 begins with the determination of a matrix C representative of all the signals c n ⁇ (t) collected at the output of the N sensors in response to the Q known fields.
  • This matrix is a matrix of size N over Q consisting of the elements C n , q (f), the index n designating the row n and the index q designating the column q.
  • the elements C are deduced from the signals c voyage > q (t) by Fourier transform.
  • the matrix C therefore has the following form:
  • the matrix C is representative of the acquisition capacities of the acquisition means 1 and of the Q acoustic fields emitted.
  • the coefficients B provokej tJtt (f) are determined from the matrices C and P using conventional methods of generalized matrix inversion, applied to the relation which binds C to P.
  • the coefficients R,., /, m (/) are placed in a matrix B determined by the following relation:
  • the matrix B is a matrix of size N on (E 3 + l) 2 consisting of the coefficients B nm if, the index n designating the row n and the indices (l, m) designating the column l 2 + l + m.
  • the matrix B therefore has the following form:
  • the substeps 26 and 28 of the calibration step 20 can be carried out in different ways, depending on the parameters to be determined.
  • the sub-steps 26 and 28 use the propagation times of the waves emitted by the loudspeakers to reach the sensors 2zen.
  • the position of each sensor 2 n is determined using at least three propagation time measurements according to triangulation methods.
  • sub-steps 26 and 28 make it possible to determine, from the signals c n> q (t), the impulse responses of each sensor 2 n when the loudspeaker emits a given pulse.
  • the calibration step 20 allows the determination of the electro-acoustic characteristics of the sensors. It then begins with the determination of the directivity diagram of each sensor 2 n for each frequency / considered, for example, by determining the frequency response of each sensor 2 n for several directions. Then, all or part of the following parameters are determined:
  • This parameter d n (f) can be determined using the usual methods of parameter estimation, for example by applying a method to the least dres squares providing the value of d n (f) which minimizes the error between the real directivity diagram and the modeled directivity diagram.
  • the calibration step 20 also makes it possible to determine the parameter ⁇ A n (J) corresponding to the power spectral density of the background noise of the sensors.
  • the signal delivered by the sensor 2 n is collected in the absence of an acoustic field.
  • the parameter ⁇ A n (f) is determined by means of power spectral density estimation methods, for example the so-called periodogram method.
  • all or part of the sub-steps 22 to 28 is repeated, for example to allow the determination of several types of parameters, certain sub-steps possibly being common to the determination of different types of parameters.
  • the calibration step 20 can also be carried out by means other than those described, such as direct measurements, for example using optical means for measuring the position of each elementary sensor 2 ⁇ relative to the center 4 of the means of acquisition 1.
  • the calibration step 20 can implement a simulation, for example using a computer, of signals representative of the acquisition capacities of the elementary sensors 2 n . It therefore appears that this calibration step 20 makes it possible to determine all or part of the parameters representative of the structural and / or electro-acoustic characteristics of the acquisition means 1, which are used during step 30 of determining the encoding filters. .
  • FIG. 5 a flow diagram of an embodiment of the step 30 for determining the encoding filters is shown.
  • Step 30 includes a sub-step 32 of determining a matrix B representative of the acquisition capacities of the acquisition means 1 or sampling matrix.
  • the matrix B is determined from the parameters ⁇ nj H Mone(f), d n (f), cc ll (f) and B n f) and is a matrix of size N over
  • Certain elements of the matrix R can be directly determined during steps 10 or 20.
  • the matrix B is then supplemented with elements determined from a modeling of the sensors.
  • each sensor n is modeled by a point sensor placed at the position ⁇ n] having a directivity composed of a combination of omnidirectional and bidirectional diagrams of proportion d n (f), oriented in the direction aAf) and having a frequency response H Tar(j).
  • Step 30 then comprises a sub-step 34 for determining an intercorrelation matrix A representative of the resemblance between the signals ci to cw delivered by the sensors 2 ⁇ to 2A / due to the fact that these sensors 2 ⁇ to 2 / v perform measurements on the same sound field P.
  • the matrix A is determined from the sampling matrix B.
  • A is a matrix of size N over N obtained by means of the relation:
  • the matrix A is determined more precisely by using a matrix B completed to an order L 2 according to the method of the previous step.
  • the sub-step 34 of determining the intercorrelation matrix A can be considered as an intermediate calculation step and can therefore be integrated into another sub -step of step 30.
  • Step 30 then comprises a sub-step 36 for determining an encoding matrix E (J) representative of the encoding filters for a given frequency.
  • the matrix E (j) is determined from the matrices A and B and the parameters L (f), ⁇ (f), ⁇ (h, nik) ⁇ (f) e ⁇ admir 2 (f).
  • the matrix E (f) is a matrix of size (L (f) +1) 2 over N consisting of elements E ⁇ > m , n (f), the indices (l, m) designating the line l 2 + l + m and the index n designating the column n.
  • the matrix E (f) therefore has the following form:
  • the matrix E (f) is determined line by line. For each operating frequency / each line E /,nd, of index (l, m) of the matrix E (f) takes the following form:
  • Substeps 32, 34 and 36 for determining the matrices A, B and E (f) are repeated for each frequency / of operation.
  • the parameters are independent of the frequency and the substeps 32, 34 and 36 are carried out only once.
  • Sub-step 36 then directly allows the determination of a matrix E independent of the frequency.
  • parameters FD representative of the encoding filters are determined from the matrix E (f).
  • Each element E ⁇ ⁇ m , n if) of the matrix E (f) represents the frequency response of an encoding filter.
  • Each encoding filter can be described by the FD parameters in different forms.
  • the parameters FD representative of the filters E / , m , " (/) are: - frequency responses, the parameters FD are then directly the E ⁇ , m> f) calculated for certain frequencies /;
  • step 30 of determining the encoding filters delivers parameters FD describing encoding filters representative of the at least structural and / or electro-acoustic capacities of the acquisition means 1.
  • these filters are representative of the following characteristics:
  • step 40 the coefficients p i ⁇ m ( ⁇ representative of the acoustic field P are deduced from the signals Ci to CN coming from the elementary sensors 2 ⁇ to 2A /, by the application of the encoding filters for response in French accordingly E ,, m, .C0 as follows.:
  • coefficients p hm are a finite number of coefficients representative in time and in the three dimensions of the space of the acoustic field and constitute a faithful representation of this acoustic field.
  • other filtering by E7, w , “(/) can be carried out according to different filtering methods, such as for example:
  • the filtering is carried out by means of filtering methods in the frequency domain such as, for example, convolution techniques by blocks;
  • the invention makes it possible to faithfully represent an acoustic field by a substantially independent representation of the characteristics of the acquisition means in the form of Fourier-Bessel coefficients.
  • the method of the invention can be implemented using only the knowledge of the parameters x announcedrepresentative of the position of the sensors 2 n relative to the center 4 of the acquisition means 1 and of the parameters ⁇ and L relating to the optimization strategy. Furthermore, in this simplified embodiment, it is considered that the parameters are independent of the frequency.
  • the matrices A and B are calculated simultaneously or sequentially in any order during sub-steps 32 and 34.
  • the elements B n ⁇ m (f) of the matrix B are then organized as follows: 22
  • the matrix A is obtained from the matrix B by means of the relation:
  • the elements A nn (j) of the matrix A are determined with better precision by the relation:
  • E 2 is the order in which the determination of the matrix ⁇ is carried out and is an integer greater than L.
  • the encoding matrix E representative of the encoding filters is determined from the matrices A and B and from the parameter ⁇ according to the expression:
  • Sub-steps 32, 34 and 36 for determining the matrices A and B then E are repeated for all of the frequencies / of operation.
  • Each element E., m , .. (H corresponds to an encoding filter which integrates the spatial distribution of the sensors 2 n as well as the optimization strategy.
  • phase 40 the signals to CM from sensors 2- ⁇ to 2N are filtered using the encoding filters described by the parameters FD.
  • Each coefficient p l> m (t) delivered is deduced from the signals c- to c ⁇ by applying the filters as follows:
  • the coefficients p lm ( ⁇ are determined by means of filtering methods in the frequency domain, such as for example block convolution techniques.
  • the representation of the sound field therefore takes into account the position of the sensors and the optimization parameters chosen and constitutes a faithful estimate of the sound field.
  • FIG 7 there is shown a block diagram of a device suitable for implementing the method of the invention.
  • a device 50 for representing the acoustic field P is connected to the acquisition means 1 as described with reference to FIG. 2.
  • the device 50 or encoding device is also connected at input to means 60 for determining parameters representative of the structural and / or electro-acoustic characteristics of the acquisition means 1.
  • means 60 include in particular means 62 for entering parameters and calibration means 64 which are suitable for implementing steps 10 and 20 respectively of the method of the invention as described above.
  • the encoding device 50 receives, means 60 for determining the parameters, a plurality of parameters representative of the characteristics of the acquisition means 1 distributed between a signal CL for defining the structural characteristics and a signal CP for configuring the structural characteristics and / or electro-acoustic.
  • the device also receives parameters relating to the representation strategies in a signal OS for optimizing the representation.
  • ⁇ (j) specifying the compromise between the fidelity of representation of the acoustic field and the minimization of the background noise provided by the sensors
  • - a parameter L (f) specifying the order in which the representation is carried out
  • this device 50 includes means 51 for shaping the input signals adapted to deliver from the signals Ci to CN, a signal SI put into corresponding shape.
  • the means 51 comprise analog-to-digital converters, amplifiers or even filtering systems.
  • the device 50 further comprises means 52 for determining the encoding filters which comprise a module 55 for calculating the sampling matrix B, a module 56 for calculating the intercorrelation matrix A, which are both connected to a module 57 for calculating the encoding matrix E (j).
  • This encoding matrix E (f) is used by a module 58 for determining encoding filters which delivers a signal SFD which contains the parameters FD representative of the encoding filters.
  • This SFD signal is used by a processing module 59 which applies the encoding filters to the signal SI in order to deliver a SIFB signal which includes the Fourier-Bessel coefficients representative of the acoustic field P.
  • the device 50 includes a non-volatile memory in which the parameters which constitute the SFD signal which have been determined beforehand are stored.
  • the acquisition means 1 are tested and calibrated by their manufacturer in order to directly supply a memory comprising all the parameters of the SFD signal which it is necessary to integrate into an encoding device in order to carry out the acquisition. of the acoustic field P and to deliver a faithful representation of the latter.
  • this memory only includes the matrices B and possibly A and the device 50 includes means for entering the parameters constituting the optimization signal OS in order to implement the determination of the matrix E (f) d encoding and determination of the FD parameters representative of the encoding filters.
  • the different modules described can be envisaged as required.

Description

Procédé et système de représentation d'un champ acoustique.
La présente invention concerne un procédé et un dispositif de représentation d'un champ acoustique à partir de signaux délivrés par des moyens d'acquisition.
Les procédés et systèmes d'acquisition et de représentation d'environnements sonores existants utilisent des modélisations basées sur des moyens d'acquisition physiquement irréalisables, notamment en ce qui concerne les caractéristiques électro-acoustiques et/ou structurelles de ces moyens d'acquisition.
Les moyens d'acquisition sont, par exemple, constitués d'un ensemble d'éléments de mesure ou capteurs élémentaires disposés en des endroits spécifiques de l'espace et présentant des caractéristiques électro-acoustiques d'acquisition intrinsèques.
Les systèmes existants sont limités par les caractéristiques structurelles des moyens d'acquisition, telles que la disposition physique des capteurs élé- mentaires ainsi que par leurs caractéristiques électro-acoustiques, et délivrent des représentations dégradées de l'environnement sonore à acquérir.
Par exemple, les systèmes regroupés sous le terme « Ambisonic » ne considèrent que les directions de provenance des sons par rapport au centre des moyens d'acquisition formés d'une pluralité de capteurs élémentaires, ce qui conduit à assimiler les moyens d'acquisition à un microphone ponctuel.
Cependant, l'impossibilité de positionner l'ensemble des capteurs élémentaires en un même point limite les performances de ces systèmes.
En outre ces systèmes représentent l'environnement sonore par une modélisation de sources virtuelles dont la distribution angulaire autour du centre permet théoriquement l'obtention d'un tel environnement sonore.
Cependant, l'absence de disponibilité de capteurs élémentaires de caractéristiques de directivité élevées restreint ces systèmes à un niveau de précision de représentation couramment appelé ordre un sur une base mathématique dite base des harmoniques sphériques. Dans d'autres systèmes, tels que celui mettant en œuvre le procédé et le dispositif d'acquisition décrits dans la demande de brevet WO-01 -58209, l'acquisition est basée sur la mesure, dans un plan, d'informations représentatives de l'environnement sonore à acquérir. Cependant, ces systèmes utilisent des modélisations basées sur des capteurs élémentaires parfaits disposés nécessairement sur un cercle et aboutissent à une amplification importante des bruits de fond des capteurs.
Ces systèmes nécessitent donc des capteurs dont le bruit de fond in- trinsèque est extrêmement faible, et sont donc irréalisables en pratique.
En outre, dans ces systèmes, l'environnement sonore est décrit uniquement par une modélisation bidimensionnelle, ce qui correspond à une approximation importante et réductrice des caractéristiques sonores réelles.
Il apparaît donc que les représentations faites des environnement so- nores par les systèmes existants sont incomplètes et dégradées et qu'aucun système ne permet d'en obtenir une représentation fidèle.
Le but de l'invention est de résoudre ce problème en fournissant un procédé et un dispositif délivrant une représentation du champ acoustique sensiblement indépendante des caractéristiques des moyens d'acquisition. La présente invention a pour objet un procédé de représentation d'un champ acoustique comprenant une étape d'acquisition de signaux de mesure délivrés par des moyens d'acquisition formés d'un ou de plusieurs capteurs élémentaires exposés audit champ acoustique, caractérisé en ce qu'il comporte :
- une étape de détermination de filtres d'encodage représentatifs de caractéristiques au moins structurelles desdits moyens d'acquisition ; et
- une étape de traitement desdits signaux de mesure par l'application desdits filtres d'encodage à ces signaux pour déterminer un nombre fini de coefficients représentatifs dans le temps et dans les trois dimensions de l'espace dudit champ acoustique, lesdits coefficients permettant d'obtenir une représentation dudit champ acoustique sensiblement indépendante des caractéristiques desdits moyens d'acquisition.
Suivant d'autres caractéristiques :
- lesdites caractéristiques structurelles comportent au moins des caractéristiques de position desdits capteurs élémentaires par rapport à un point de référence prédéterminé desdits moyens d'acquisition ;
- lesdits filtres d'encodage sont en outre représentatifs de caractéristiques électro-acoustiques desdits moyens d'acquisition ; - lesdites caractéristiques électro-acoustiques comportent au moins des caractéristiques liées aux capacités électro-acoustiques d'acquisition intrinsèques desdits capteurs élémentaires ;
- lesdits coefficients permettant d'obtenir une représentation du champ acoustique sont des coefficients dits de Fourier-Bessel et/ou des combinaisons linéaires de coefficients de Fourier-Bessel ;
- ladite étape de détermination des filtres d'encodage comprend : une sous-étape de détermination d'une matrice d'échantillonnage représentative des capacités d'acquisition desdits moyens d'acquisition ; une sous-étape de détermination d'une matrice d'intercorrélation représentative de la ressemblance entre lesdits signaux de mesure délivrés par les capteurs élémentaires formant lesdits moyens d'acquisition ; et - une sous-étape de détermination d'une matrice d'encodage à partir de ladite matrice d'échantillonnage, de ladite matrice d'intercorrélation, et d'un paramètre représentatif d'un compromis souhaité entre la fidélité de représentation du champ acoustique et la minimisation du bruit de fond induit par les moyens d'acquisition, laquelle matrice est représentative desdits filtres d'encodage ;
- lesdites sous-étapes de détermination des matrices sont réalisées pour un nombre fini de fréquences de fonctionnement ;
- ladite étape de détermination de la matrice d'échantillonnage est réalisée, pour chacun desdits capteurs élémentaires formant lesdits moyens d'acquisition, à partir :
- de paramètres représentatifs de la position dudit capteur par rapport au centre desdits moyens d'acquisition ; et/ou
- d'un nombre fini de coefficients représentatifs des capacités d'acquisition dudit capteur ; - ladite étape de détermination de la matrice d'échantillonnage est réalisée en outre à partir au moins d'un des paramètres parmi :
- des paramètres représentatifs des réponses en fréquence de tout ou partie des capteurs ; - des paramètres représentatifs des diagrammes de directivité de tout ou partie des capteurs ;
- des paramètres représentatifs des orientations de tout ou partie des capteurs, à savoir de leur direction de sensibilité maximale ; - des paramètres représentatifs des densités spectrales de puissance du bruit de fond de tout ou partie des capteurs ;
- d'un paramètre spécifiant l'ordre auquel est conduite la représentation ; et
- d'un paramètre représentatif d'une liste de coefficients dont on exige que la puissance soit égale à la puissance du coefficient correspondant dans le champ acoustique à représenter ;
- il comporte une étape de calibrage permettant de délivrer tout ou partie des paramètres utilisés dans ladite étape de détermination des filtres d'encodage ; - ladite étape de calibrage comporte, pour au moins l'un desdits capteurs élémentaires formant lesdits moyens d'acquisition :
- une sous-étape d'acquisition de signaux représentatifs des capacités d'acquisition dudit au moins un capteur ; et
- une sous-étape de détermination de paramètres représentatifs de caractéristiques électro-acoustiques et ou structurelles dudit au moins un capteur ;
- ladite étape de calibrage comporte en outre :
- une sous-étape d'émission d'un champ acoustique spécifique vers ledit au moins un capteur, ladite sous-étape d'acquisition correspondant à l'acquisition des signaux délivrés par ce capteur lorsque exposé audit champ acoustique spécifique ; et
- une sous-étape de modélisation dudit champ' acoustique spécifique en un nombre fini de coefficients afin de permettre la réalisation de ladite sous-étape de détermination de paramètres représentatifs de caractéristiques électro-acoustiques et ou structurelles du capteur ;
- ladite étape de calibrage comporte une sous-étape de réception d'un nombre fini de signaux représentatifs des caractéristiques électro-acoustiques et structurelles desdits capteurs formant lesdits moyens d'acquisition, lesquels signaux sont directement utilisés lors de ladite sous-étape de détermination des caractéristiques électro-acoustiques et/ou structurelles desdits moyens d'acquisition ; et
- il comporte une étape de saisie permettant de déterminer tout ou partie des paramètres utilisés lors de ladite étape de détermination des filtres d'encodage.
L'invention concerne également un programme d'ordinateur comprenant des instructions de code de programme pour l'exécution des étapes du procédé, tel que décrit précédemment lorsque ledit programme est exécuté sur un ordinateur. L'invention concerne également un support mobile du type comportant au moins un processeur de traitement et un élément de mémoire non volatile, caractérisé en ce que ladite mémoire comprend un programme comportant des instructions de code pour l'exécution des étapes du procédé tel que décrit précédemment lorsque ledit processeur exécute ledit programme. L'invention concerne également un dispositif de représentation d'un champ acoustique connectable à des moyens d'acquisition formés d'un ou plusieurs capteurs élémentaires délivrant des signaux de mesure lorsqu'ils sont exposés audit champ acoustique, caractérisé en ce qu'il comporte un module de traitement des signaux de mesure par l'application de filtres d'encodage repré- sentatifs de caractéristiques au moins structurelles desdits moyens d'acquisition à ces signaux de mesure pour délivrer un signal qui comporte un nombre fini de coefficients représentatifs dans le temps et dans les trois dimensions de l'espace dudit champ acoustique, lesdits coefficients permettant d'obtenir une représentation dudit champ acoustique sensiblement indépendante des caractéristiques desdits moyens d'acquisition,
Selon d'autres caractéristiques du dispositif de l'invention :
- lesdits filtres d'encodage sont en outre représentatifs de caractéristiques électro-acoustiques desdits moyens d'acquisition ;
- il comporte en outre des moyens de détermination desdits filtres d'encodage représentatifs de caractéristiques structurelles et/ou électroacoustiques desdits moyens d'acquisition ;
- lesdits moyens de détermination de filtres d'encodage reçoivent en entrée au moins l'un des paramètres parmi les paramètres suivants : - des paramètres représentatifs des positions par rapport au centre desdits moyens d'acquisition de tout ou partie des capteurs ;
- un nombre fini de coefficients représentatifs des capacités d'acquisition de tout ou partie des capteurs ; - des paramètres représentatifs des réponses en fréquence de tout ou partie des capteurs ;
- des paramètres représentatifs des diagrammes de directivité de tout ou partie des capteurs ;
- des paramètres représentatifs des orientations de tout ou partie des capteurs, à savoir de leur direction de sensibilité maximale ;
- des paramètres représentatifs des densités spectrales de puissance du bruit de fond de tout ou partie des capteurs ;
- un paramètre représentatif du compromis souhaité entre la fidélité de représentation du champ acoustique et la minimisation du bruit de fond induit par les moyens d'acquisition ;
- un paramètre spécifiant l'ordre auquel est conduit l'encodage ; et
- un paramètre représentatif d'une liste de coefficients dont on exige que la puissance soit égale à la puissance du coefficient correspondant dans le champ acoustique à représenter ;
- il est associé à des moyens de détermination de tout ou partie des paramètres reçus par lesdits moyens de détermination des filtres d'encodage, lesdits moyens comportant au moins l'un des éléments suivants :
- des moyens de saisie des paramètres ; et/ou - des moyens de calibrage ;
- il est associé à des moyens de mise en forme desdits signaux de mesure afin de délivrer un signal mis en forme correspondant.
L'invention sera mieux comprise à la lecture de la description qui va suivre, donnée uniquement à titre d'exemple et faite en se référant aux dessins annexés, sur lesquels :
- la Fig.1 est une représentation d'un repère sphérique ;
- la Fig.2 est un schéma représentatif de moyens d'acquisition utilisés ;
- la Fig.3 est un organigramme général du procédé de l'invention : - la Fig.4 est un organigramme du détail d'un mode de réalisation de l'étape de calibrage du procédé de l'invention ;
- la Fig.5 est un organigramme du détail d'un mode de réalisation de l'étape de détermination des filtres d'encodage du procédé de l'invention ; - la Fig.6 est un schéma du détail d'un mode de réalisation de l'étape d'application des filtres d'encodage ; et
- la Fig.7 est un schéma synoptique d'un dispositif adapté pour la mise en œuvre du procédé de l'invention.
Sur la figure 1 , on a représenté un repère sphérique classique, de manière à préciser le système de coordonnées auquel il est fait référence dans le texte.
Ce repère est un repère orthonormal, d'origine O et comportant trois axes (QX), (OY) et (02).
Dans ce repère, une position notée x est décrite au moyen de ses coordonnées sphériques (r,θ, ), où r désigne la distance par rapport à l'origine O, θ l'orientation dans le plan vertical et l'orientation dans le plan horizontal.
Dans un tel repère, un champ acoustique est connu si l'on définit en tout point à chaque instant t la pression acoustique notée p(r, θ, φ,t), dont la transformée de Fourier est notée P(r,θ,φJ) où /désigne la fréquence. Le procédé de l'invention se base sur l'utilisation de fonctions spatiotemporelles permettant de décrire un champ acoustique quelconque dans le temps et dans les trois dimensions de l'espace.
Dans les modes de réalisation décrits, ces fonctions sont des fonctions dites de Fourier-Bessel sphériques de première espèce, appelées par la suite fonctions de Fourier-Bessel.
Dans une zone vide de sources et vide d'obstacles, les fonctions de Fourier-Bessel correspondent aux solutions de l'équation des ondes et constituent une base qui engendre tous les champs acoustiques produits par des sources situées à l'extérieur de cette zone. Tout champ acoustique tridimensionnel peut donc s'exprimer par une combinaison linéaire des fonctions de Fourier-Bessel, selon l'expression de la transformée de Fourier-Bessel inverse qui s'exprime :
Dans cette équation, les termes Pi f) sont définis comme les coefficients de Fourier-Bessel du champ p(r,θ,φ,t), k=~^- , c est la célérité du son dans
l'air (340 ms*1), j (kr) est la fonction de Bessel sphérique de première espèce
d'ordre / définie par où Jv(x) est la fonction de Bessel de pre- mière espèce d'ordre v, et yf(θ,φ) est l'harmonique sphérique réelle d'ordre / et de terme m, avec m allant de -/ à /, définie par :
^(^^)=R/ |m|(cost?)trgm(^) avec :
Dans cette équation, les Pj m(x) sont les fonctions de Legendre associées définies par :
avec Pι(x) les polynômes de Legendre, définis par :
Les coefficients de Fourier-Bessel s'expriment aussi dans le domaine temporel par les coefficients p!>m(t) correspondant à la transformée de Fourier temporelle inverse des coefficients Pι,mf).
Dans d'autres modes de réalisation, le champ acoustique est décomposé sur une base de fonctions, où chacune des fonctions s'exprime par une combinaison linéaire éventuellement infinie de fonctions de Fourier-Bessel.
Sur la figure 2, on a représenté schématiquement des moyens d'acquisition 1 formés de N capteurs élémentaires 2-ι à 2N.
Ces capteurs élémentaires sont disposés en des points spécifiques de l'espace autour d'un point 4 prédéterminé désigné comme le centre des moyens d'acquisition 1. Ainsi, la position de chaque capteur élémentaire peut s'exprimer dans l'espace dans un repère sphérique tel que celui décrit en référence à la figure 1 , centré sur le centre 4 des moyens d'acquisition 1.
Lorsqu'il est exposé à un champ acoustique P, chaque capteur 2n des moyens d'acquisition 1 délivre un signal de mesure c„qui correspond à la mesure faite par ce capteur dans le champ acoustique P.
Ainsi, les moyens d'acquisition 1 délivrent une pluralité de signaux Ci à c/v qui sont les signaux de mesure du champ acoustique P par les moyens d'acquisition 1. Ces signaux de mesure c-i à cw délivrés par les moyens d'acquisition 1 sont donc directement liés aux capacités d'acquisition des capteurs élémentaires
Sur la figure 3, on a représenté un organigramme général du procédé de l'invention. Le procédé débute par une étape 10 de saisie de paramètres et une étape 20 de calibrage des moyens d'acquisition, qui permettent de définir un ensemble de paramètres représentatifs des caractéristiques structurelles et/ou électro-acoustiques des moyens d'acquisition 1.
Certains paramètres et notamment des paramètres représentatifs de caractéristiques électro-acoustiques sont dépendants de la fréquence.
L'étape 10 de saisie et l'étape 20 de calibrage, laquelle est décrite plus en détail en référence à la figure 4, peuvent être réalisées simultanément ou dans un ordre quelconque.
De même, le procédé de l'invention peut ne comporter que l'étape 10 de saisie.
Les étapes 10 de saisie et 20 de calibrage permettent de déterminer, pour un ou plusieurs capteurs, tout ou partie des paramètres suivants :
- des paramètres x„ représentatifs de la position du capteur 2n par rapport au centre 4 des moyens d'acquisition 1 , qui s'écrivent en coordonnées sphériques (rn, θ„, φn) ;
- des paramètres d„(f) représentatifs du diagramme de directivité du capteur 2„ qui peut prendre toutes les valeurs entre 0 et 1 et permet de décrire la directivité du capteur 2n par une combinaison de diagrammes omnidirectionnels et de diagrammes bidirectionnels : si d„(f) ≈ 0, le capteur est omnidirectionnel si dn(f) = V_, le capteur est cardioïde si dn( ) = 1, le capteur est bidirectionnel ;
- des paramètres a ) représentatifs de l'orientation du capteur 2n c'est-à-dire de sa direction de sensibilité maximale qui est donnée par le couple d'angles (θfï,φn)(f) ;
- des paramètres Hn(f) représentatifs de la réponse en fréquence du capteur 2n correspondant, pour chaque fréquence/; à la sensibilité du capteur 2n dans la direction aJJ) ; - des paramètres cfn(f) représentatifs de la densité spectrale de puissance du bruit de fond du capteur 2π ;
- des paramètres Bn>ι>m(f) représentatifs des capacités d'acquisition du capteur 2n, c'est-à-dire de la façon dont le capteur 2n prélève des informations sur le champ acoustique P. Ainsi chaque 5„,/jW( est représentatif des capacités d'acquisition d'un capteur et notamment de sa position dans l'espace et l'ensemble des Rn,.,m(/) est représentatif de l'échantillonnage du champ acoustique P réalisé par les moyens d'acquisition 1 ;
- un paramètre μ(f) spécifiant un compromis entre la fidélité de représentation du champ acoustique P et la minimisation du bruit de fond apporté par les capteurs 2ι à n et pouvant prendre toutes les valeurs entre 0 et 1 :
- si μ(f) ≈ 0, le bruit de fond est minimal ;
- si μ(J) = 1, la qualité spatiale est maximale ;
- un paramètre L(f) spécifiant l'ordre auquel est conduite la représentation ; et - un paramètre {(k,mk)}(f) représentatif d'une liste de coefficients dont on exige que la puissance soit égale à la puissance du coefficient correspondant dans le champ acoustique à représenter.
Dans des modes de réalisation simplifiés, tout ou partie des paramètres décrits est considéré comme indépendant de la fréquence. Les paramètres μ(f), L(f) et {(k,mk)}(f) sont représentatifs des stratégies d'optimisation permettant de maîtriser l'extraction d'informations spatiotemporelles du champ acoustique P à partir des signaux de mesure c-i à C/v et sont saisis lors de l'étape 10 de saisie. Les autres paramètres peuvent être saisis lors de l'étape 10 de saisie ou déterminés lors de l'étape 20 de calibrage. Dans des modes de réalisation simplifiés, le procédé de l'invention est réalisé uniquement avec les paramètres μ(f), L(f) et l'ensemble des paramètres x„ ou l'ensemble des paramètres Bnj>m(f) ou une combinaison de paramètres x„ et Bn j), de manière à disposer d'au moins un paramètre par capteur élémen- taire 2n.
Bien entendu, tout ou partie des paramètres utilisés peut être délivré par des mémoires ou des dispositifs dédiés, ces techniques étant assimilables à l'étape 10 de saisie directe par un opérateur telle que décrite.
A l'issue des étapes 10 de saisie et/ou 20 de calibrage, le procédé comporte une étape 30 de détermination de filtres d'encodage représentatifs des caractéristiques au moins structurelles et avantageusement électro-acoustiques des moyens d'acquisition 1.
Cette étape 30, décrite plus en détail en référence à la figure 5, permet de prendre en compte tous les paramètres déterminés lors des étapes 10 de sai- sie et/ou 20 de calibrage.
Ces filtres d'encodage sont donc représentatifs au moins des caractéristiques de position des capteurs élémentaires 2n par rapport au point de référence 4 des moyens d'acquisition 1.
Avantageusement, ces filtres sont également représentatifs d'autres caractéristiques structurelles des moyens d'acquisition 1, telles que l'orientation des capteurs élémentaires 2ι à 2N OU leurs influences mutuelles, ainsi que leurs capacités électro-acoustiques d'acquisition et notamment leur bruit de fond, leur diagramme de directivité, leur réponse en fréquence, ...
Les filtres d'encodage obtenus à l'issue de l'étape 30 peuvent être mémorisés, de sorte que les étapes 10, 20 et 30 ne sont répétées qu'en cas de modification des moyens d'acquisition 1 ou des stratégies d'optimisation.
Ces filtres d'encodage sont appliqués lors d'une étape 40 de traitement des signaux c-i à c/v issus des capteurs élémentaires 2ι à 2/v.
Ce traitement correspond à un filtrage des signaux et à des combinai- sons des signaux filtrés.
A l'issue de cette étape 40 de traitement des signaux de mesure par l'application des filtres d'encodage à ces signaux, un nombre fini de coefficients représentatifs dans le temps et dans les trois dimensions de l'espace du champ acoustique P est délivré. Ces coefficients sont des coefficients dits de Fourier-Bessel, notés Pι,m(f) et correspondent à une représentation du champ acoustique P sensiblement indépendante des caractéristiques des moyens d'acquisition 1.
Il apparaît donc que grâce au procédé de l'invention, on obtient une représentation fidèle du champ acoustique dont on transcrit des caractéristiques temporelles et spatiales quels que soient les moyens d'acquisition utilisés.
Sur la figure 4, on a représenté un organigramme d'un mode de réalisation de l'étape 20 de calibrage.
Dans ce mode de réalisation, l'étape 20 de calibrage permet de dé- terminer directement les coefficients Bni f représentatifs des capacités d'acquisition des moyens d'acquisition 1.
Cette étape 20 commence par une sous-étape 22 d'émission d'un champ acoustique spécifique vers les moyens d'acquisition 1 et par une sous- étape 24 d'acquisition de signaux de mesure par les moyens d'acquisition 1 ex- posés au champ acoustique émis.
Ces sous-étapes 22 et 24 sont répétées pour une pluralité Q de champs spécifiques différents et requièrent des moyens de génération de champs acoustiques spécifiques et des moyens de déplacement et/ou de rotation des moyens d'acquisition 1. Par exemple, l'étape 20 de calibrage est mise en œuvre à l'aide de moyens de génération d'un champ acoustique qui ne comportent qu'un haut- parleur fixe, supposé ponctuel et de réponse en fréquence plate, le haut-parleur et les moyens d'acquisition 1 étant placés dans un environnement anéchoïque.
A chaque sous-étape 22 de génération, le haut-parleur émet le même champ acoustique et les moyens d'acquisition 1 sont placés à la même position mais ils sont orientés selon des directions différentes et connues.
Bien entendu, il est également possible de déplacer le haut-parleur.
Ainsi, dans le repère des moyens d'acquisition 1 , le haut-parleur est dans une position (]{ψ ,θqp ,φqp) différente pour chaque champ q généré. Les moyens d'acquisition 1 sont ainsi exposés à un champ acoustique q dont les coefficients de Fourier-Bessel Pι,m,q{f), dans le repère des moyens d'acquisition 1 , sont connus jusqu'à un ordre donné, noté E3.
Dans le mode de réalisation décrit, les signaux de mesure délivrés à la suite de la sous-étape 24 d'acquisition sont un nombre fini de coefficients repré- sentatifs du champ acoustique q généré, ainsi que des capacités d'acquisition des moyens d'acquisition 1.
Les paramètres E3 et Q sont choisis pour respecter la condition : Q≥ iLAlf
Avantageusement, par la suite, le procédé comporte une sous-étape 26 de modélisation permettant de déterminer une représentation des Q champs acoustiques émis lors de la sous-étape 22.
Ainsi, au cours de la sous-étape 26, une matrice P de modélisation représentative de l'ensemble des Q champs connus auxquels sont successivement exposés les moyens d'acquisition 1 est déterminée. Cette matrice P est une matrice de taille (E3+l)2 sur Q constituée des éléments Pι,m,q(f), les indices (l,m) désignant la ligne l2+l+m et l'indice q désignant la colonne q. La matrice P a donc la forme suivante :
I3,i3,β (f)
Dans le mode de réalisation décrit, le champ acoustique produit par le haut-parleur est modélisé par un rayonnement sphérique, ainsi, dans le repère des moyens d'acquisition 1 , les coefficients Pι,m,q(f) de chaque champ acoustique q ainsi généré sont connus, grâce à la relation :
-B,« ,,(/)= r. !Ψ ~ξ>(rJφ,f)yr(θïp,Φ$p)
« «^ u t ko . 1 ( +*)! f j2 πrq hp f avec ξ,(r9*,f) =∑ „> 2A**! (i/--*k)! I— ~
Les coefficients obtenus lors de la sous-étape 26 sont ensuite utilisés lors d'une sous-étape 28 afin de déterminer des paramètres représentatifs de caractéristiques structurelles et/ou acoustiques des moyens d'acquisition 1.
Dans le mode de réalisation décrit, cette sous-étape 28 utilise également la matrice R de modélisation déterminée lors de la sous-étape 26. Cette sous-étape 28 débute par la détermination d'une matrice C représentative de l'ensemble des signaux cnΛ(t) recueillis en sortie des N capteurs en réponse aux Q champs connus. Cette matrice est une matrice de taille N sur Q constituée des éléments Cn,q(f), l'indice n désignant la ligne n et l'indice q désignant la colonne q. Les éléments C„ιgf) se déduisent des signaux c„>q(t) par transformée de Fourier. La matrice C a donc la forme suivante :
' Cl (f) Cltl(f)-Cv,Q(fi C2,}(/)C_,2if)-C2>g(f)
CNλ(f)C (f -C f)_
La matrice C est représentative des capacités d'acquisition des moyens d'acquisition 1 et des Q champs acoustiques émis.
Dans le mode de réalisation décrit, au cours de la sous-étape 28, les coefficients B„jtJtt(f) sont déterminés à partir des matrices C et P en utilisant des méthodes classiques d'inversion matricielle généralisée, appliquées à la relation qui lie C à P. Par exemple, les coefficients R,.,/,m(/) sont placés dans une matrice B déterminée par la relation suivante :
B = cpτ(ppτ l
B est une matrice de taille N sur (E3+l)2 constituée des coefficients Bn mif , l'indice n désignant la ligne n et les indices (l,m) désignant la colonne l2+l+m. La matrice B a donc la forme suivante :
BnMt r-Wm NM • - N^ if)- • - ^if)- • -BN^L if}_
Ces sous-étapes 26 et 28 sont réalisées pour chaque fréquence de fonctionnement et les coefficients ainsi déterminés constituent directement les paramètres représentatifs des capacités d'acquisition des moyens d'acquisition 1.
Les sous-étapes 26 et 28 de l'étape 20 de calibrage peuvent être réalisées de différentes manières, en fonction des paramètres devant être déterminés.
Par exemple, dans le cas où l'étape 20 de calibrage permet de déterminer la position χ„ de chaque capteur 2n, les sous-étapes 26 et 28 exploitent les temps de propagation des ondes émises par les haut-parleurs pour atteindre les capteurs 2„. La position de chaque capteur 2n est déterminée à l'aide d'au moins trois mesures de temps de propagation selon des méthodes de triangulation.
Dans un autre cas, les sous-étape 26 et 28 permettent de déterminer, à partir des signaux cn>q(t), les réponses impulsionnelles de chaque capteur 2n lorsque le haut-parleur émet une impulsion donnée.
Par exemple, on utilise dans ce cas les techniques usuelles de détermination de réponses impulsionnelles, telles que MLS (en anglais : Maximum Lenght Séquence). Avantageusement, l'étape 20 de calibrage permet la détermination de caractéristiques électro-acoustiques des capteurs. Elle débute alors par la détermination du diagramme de directivité de chaque capteur 2n pour chaque fréquence / considérée, par exemple, en déterminant la réponse en fréquence de chaque capteur 2n pour plusieurs directions. Dans un second temps, tout ou partie des paramètres suivants est déterminé :
- des paramètres an(f) représentatifs de l'orientation de chaque capteur 2n, c'est-à-dire de sa direction de sensibilité maximale donnée par les angles (θi,φ )(j) pour lesquels le diagramme de directivité admet un maximum à la fré- quence courante/;
- des paramètres H„(f représentatifs de la réponse en fréquence de chaque capteur 2n dans la direction de sensibilité maximale qui correspond donc à la valeur du diagramme de directivité pour la direction ; et
- des paramètres d„ f) représentatifs du diagramme de directivité de chaque capteur qui permet de décrire la directivité de chaque capteur par un modèle constitué d'une combinaison de diagrammes omnidirectionneis et de diagrammes bidirectionnels orientés selon la direction „( ), à l'aide du modèle de directivité suivant :
1 - dn( ) + dn(J) cos( α..( ).(6> 0 ) où θn(f).(θ, ) désigne le produit scalaire entre les directions ocnif) et
(θ, ).
Ce paramètre dn(f) peut être déterminé à l'aide des méthodes usuelles d'estimation de paramètres, par exemple en appliquant une méthode aux moin- dres carrés fournissant la valeur de dn(f) qui minimise l'erreur entre le diagramme de directivité réel et le diagramme de directivité modélisé.
Avantageusement, l'étape 20 de calibrage permet également de déterminer le paramètre <An(J) correspondant à la densité spectrale de puissance du bruit de fond des capteurs. Ainsi, au cours de cette étape 20, on recueille le signal délivré par le capteur 2n en l'absence de champ acoustique. Le paramètre <An(f) est déterminé au moyen de méthodes d'estimation de densité spectrale de puissance, par exemple la méthode dite du périodogramme.
En fonction des modes de réalisation, tout ou partie des sous-étapes 22 à 28 est répété, par exemple pour permettre la détermination de plusieurs types de paramètres, certaines sous-étapes pouvant être communes à la détermination de différents types de paramètres.
L'étape 20 de calibrage peut également être réalisée par d'autres moyens que ceux décrits tels que des mesures directes, par exemple à l'aide de moyens de mesure optiques de la position de chaque capteur élémentaire 2π par rapport au centre 4 des moyens d'acquisition 1.
En outre, l'étape 20 de calibrage peut mettre en œuvre une simulation, par exemple à l'aide d'un ordinateur, de signaux représentatifs des capacités d'acquisition des capteurs 2n élémentaires. II apparaît donc que cette étape 20 de calibrage permet de déterminer tout ou partie des paramètres représentatifs des caractéristiques structurelles et/ou électro-acoustiques des moyens d'acquisition 1 , qui sont utilisés lors de l'étape 30 de détermination des filtres d'encodage.
Sur la figure 5, on a représenté un organigramme d'un mode de réali- sation de l'étape 30 de détermination des filtres d'encodage.
L'étape 30 comporte une sous-étape 32 de détermination d'une matrice B représentative des capacités d'acquisition des moyens d'acquisition 1 ou matrice d'échantillonnage.
Dans le mode de réalisation décrit, la matrice B est déterminée à partir des paramètres χn j H„(f), dn(f), ccll(f) et Bn f) et est une matrice de taille N sur
(L(f)+lf constituée d'éléments Bnj,m(f), l'indice n désignant la ligne n et les indices (l,m) désignant la colonne l2+l+m. La matrice B a donc la forme suivante : .0.0O7 ,ι,-ι( ) .1.0O7 BlΛΛ(f)- • • B L,_L(f)- • • BhLM • • BlJJf) ' £2,o.o(/) B2Λ l(f) B2Λi0(f) B2Λ (f)- • • Bu -M- • • BlrL/f)> • • B2>LiL(f)
BNAoif)BNΛ>- NΛfii )BN ιif)- • -BN^L(f)- • .5Λii0( > • -5W(/).
Certains éléments de la matrice R peuvent être directement déterminés lors des étapes 10 ou 20. La matrice B est ensuite complétée avec des éléments déterminés à partir d'une modélisation des capteurs.
Dans ce mode de réalisation, chaque capteur n est modélisé par un capteur ponctuel placé à la position χn ] présentant une directivité composée d'une combinaison de diagrammes omnidirectionnels et bidirectionnels de proportion dn(f), orienté dans la direction aAf) et possédant une réponse en fréquence H„(j).
Les éléments Bnχm(j) complémentaires sont alors déterminés selon la relation
f. (kt )yr(θn,φ,)ur - yf" θ ,φ, ou :
fiΦA - l Ji- "») - (*+!) JiAkr»)
21+1
et où ur ≈sinθn sinθ,f(f)cos(φn -φ,?(f)) + ∞sθr, cosθ(ï(f) uθ=cosθn siιι<9.?( )cos( φn -φ>?(f)) ~ sinθn cosθfï(f)
Dans le cas où les capteurs sont orientés radialement, la relation admet une expression plus simple :
lωfl≈ H'M Jl ( -) ( W-C ) f -n ) - j dnif) ^frθ^+l)^*»^ L'étape 30 comporte ensuite une sous-étape 34 de détermination d'une matrice A d'intercorrélation représentative de la ressemblance entre les signaux ci à cw délivrés par les capteurs 2ι à 2A/ due au fait que ces capteurs 2\ à 2/v réalisent des mesures sur un même champ acoustique P. La matrice A est déterminée à partir de la matrice B d'échantillonnage. A est une matrice de taille N sur N obtenue au moyen de la relation :
A ≈B BT
Avantageusement, la matrice A est déterminée plus précisément en utilisant une matrice B complétée jusqu'à un ordre L2 selon la méthode de l'étape précédente.
La matrice A pouvant s'exprimer en fonction uniquement de la matrice B, la sous-étape 34 de détermination de la matrice A d'intercorrélation peut être considérée comme une étape intermédiaire de calcul et peut à ce titre, être intégrée dans une autre sous-étape de l'étape 30.
L'étape 30 comprend ensuite une sous-étape 36 de détermination d'une matrice E(J) d'encodage représentative des filtres d'encodage pour une fréquence donnée. La matrice E(j) est déterminée à partir des matrices A et B et des paramètres L(f), μ(f), {(h,nik)}(f) e σ„2(f). La matrice E(f) est une matrice de taille (L(f)+1)2 sur N constituée d'éléments Eι>m,n(f), les indices (l,m) désignant la ligne l2+l+m et l'indice n désignant la colonne n. La matrice E(f) a donc la forme suivante :
La matrice E(f) est déterminée ligne par ligne. Pour chaque fréquence de fonctionnement / chaque ligne E/,„, d'indice (l,m) de la matrice E(f) prend la forme suivante :
[E m (f)E m,2(f) EhmιN(f)]
Les éléments E/,m,„(/) de la ligne Eι>m s'obtiennent par les expressions suivantes : - si (l,m) appartient à la liste {(k,nik)}(f) alors :
El>m=μ(j)Blm((μ(j)-λ)A+(l-μ(fj)ΣN)-1 où λ vérifie la relation :
( )2^ (^ω-^+(i-^ )^)~^((^ -^ +(i-^ω)^w)_IR =ι et où la valeur de λ est déterminée au moyen de méthodes analytiques ou numériques de recherche de racines d'équations, en utilisant éventuellement des méthodes de diagonalisation de matrices ; et
- si (l,m) n'appartient pas à la liste {(k,mù}if)> alors :
El>n ~μ(f)Blm(μ( )A+(l-μ(J))ΣNγ Dans ces expressions, R/jW est la colonne (l,m) de la matrice B et ΣN est une matrice diagonale de taille N sur N représentative du bruit de fond des capteurs où l'élément n de la diagonale est σ„z(f).
Les sous-étapes 32, 34 et 36 de détermination des matrices A, B et E(f) sont répétées pour chaque fréquence/de fonctionnement. Bien entendu, dans des modes de réalisation simplifiés, les paramètres sont indépendants de la fréquence et les sous-étapes 32, 34 et 36 sont effectuées une seule fois. La sous-étape 36 permet alors directement la détermination d'une matrice E indépendante de la fréquence.
Lors d'une sous-étape 38 suivante, des paramètres FD représentatifs des filtres d'encodage sont déterminés à partir de la matrice E(f). Chaque élément Eιιm,nif) de la matrice E(f) représente la réponse en fréquence d'un filtre d'encodage. Chaque filtre d'encodage peut être décrit par les paramètres FD sous différentes formes.
Par exemple, les paramètres FD représentatifs des filtres E/,m,«(/) sont : - des réponses en fréquence, les paramètres FD sont alors directement les Eι,m> f) calculés pour certaines fréquences/;
- des réponses impulsionnelles finies e/,m,„(t) calculées par transformée de Fourier inverse de E/,m,„( ), chaque réponse impulsionnelle eι,m,„(t) est échantillonnée puis tronquée à une longueur propre à chaque réponse ; et - des coefficients de filtres récursifs à réponses impulsionnelles infinies calculées à partir des E!ιm,n(f) avec des méthodes d'adaptation. Ainsi, l'étape 30 de détermination des filtres d'encodage délivre des paramètres FD décrivant des filtres d'encodage représentatifs des capacités au moins structurelles et/ou électro-acoustiques des moyens d'acquisition 1.
Notamment, ces filtres sont représentatifs des caractéristiques suivan- tes :
- position des capteurs 2ι à 2N ,'
- caractéristiques électro-acoustiques intrinsèques des capteurs 2ι à 2N, notamment densité spectrale de puissance du bruit de fond et capacités d'acquisition du champ acoustique ; et - stratégies d'optimisation, notamment le compromis entre la fidélité spatiale d'acquisition du champ acoustique et la minimisation du bruit de fond apporté par les capteurs.
Sur la figure 6, on a représenté le détail d'un mode de réalisation de l'étape 40 de traitement des signaux de mesure délivrés par les moyens d'acquisition 1 par l'application des filtres d'encodage à ces signaux et par sommation des signaux filtrés.
Lors de l'étape 40, les coefficients piιm(ή représentatifs du champ acoustique P sont déduits à partir des signaux Ci à CN issus des capteurs élémentaires 2ι à 2A/, par l'application des filtres d'encodage de réponse en fré- quence E,,m,..C0 de la façon suivante :
où Pm(f) est la transformée de Fourier de β m(t) et Cn(J) est la transformée de Fourier de cn(t).
Dans l'exemple, on a décrit le cas d'un filtrage par réponse impulsion- nelle finie. Ce filtrage requiert la détermination initialement d'un paramètre Tnχm, correspondant au nombre d'échantillons propre à chaque réponse enχ1n(t), ce qui conduit à l'expression de convolution suivante :
Ces coefficients phm sont un nombre fini de coefficients représentatifs dans le temps et dans les trois dimensions de l'espace du champ acoustique et constituent une représentation fidèle de ce champ acoustique. Selon la nature des paramètres FD, d'autres filtrages par E7,w,«(/) peuvent être réalisés selon des méthodes de filtrage différentes, telles que par exemple :
- si les paramètres FD fournissent directement les réponses en fré- quence Eι!in>n(f), le filtrage est effectué au moyen de méthodes de filtrage dans le domaine fréquentiel comme par exemple, des techniques de convolution par blocs ;
- si les paramètres FD fournissent la réponse impulsionnelle finie e/,m,«( » le filtrage est effectué dans le domaine temporel par convolution ; et - si les paramètres FD fournissent les coefficients d'un filtre récursif à réponse impulsionnelle infinie, le filtrage est effectué dans le domaine temporel au moyen de la relation de récurrence.
Il apparaît donc que l'invention permet de représenter fidèlement un champ acoustique par une représentation sensiblement indépendante des carac- téristiques des moyens d'acquisition sous la forme de coefficients de Fourier- Bessel.
Par ailleurs, ainsi que cela a été dit précédemment, le procédé de l'invention peut être mis en œuvre dans des modes de réalisation simplifiés.
Par exemple, si tous les capteurs 2ι à 2 / sont sensiblement omnidirec- tionnels et sensiblement identiques en sensibilité et en niveau de bruit de fond, le procédé de l'invention peut être mis en œuvre à l'aide uniquement de la connaissance des paramètres x„ représentatifs de la position des capteurs 2n par rapport au centre 4 des moyens d'acquisition 1 et des paramètres μ et L relatifs à la stratégie d'optimisation. En outre, dans ce mode de réalisation simplifié, on considère que les paramètres sont indépendants de la fréquence.
Ainsi, à l'aide de ces paramètres; les matrices A et B sont calculées simultanément ou séquentiellement dans un ordre quelconque lors des sous- étapes 32 et 34. Les éléments BnÂm(f) de la matrice B sont alors organisés de la manière suivante : 22
avec :
De même, les éléments n π f) de la matrice -t sont alors organisés de la manière suivante :
_
Dans ce mode de réalisation, la matrice A est obtenue à partir de la matrice B au moyen de la relation :
A =B BT
Avantageusement, les éléments An n (j) de la matrice A sont déterminés avec une meilleure précision par la relation :
- φ où E2 est l'ordre auquel est conduite la détermination de la matrice^ et est un entier supérieur à L. Plus E2 sera choisi grand, plus le calcul des A n ,nτ (f) sera précis mais long.
Lors de la sous-étape 36, la matrice d'encodage E représentative des filtres d'encodage est déterminée à partir des matrices A et B et du paramètre μ selon l'expression :
E = μBτ( μA + (l-μ)IN
Les éléments E/,?„,„(/) de la matrice E sont organisés de la manière suivante :
Les sous-étapes 32, 34 et 36 de détermination des matrices A et B puis E sont répétées pour l'ensemble des fréquences/de fonctionnement.
Chaque élément E.,m,..(H correspond à un filtre d'encodage qui intègre la répartition spatiale des capteurs 2n ainsi que la stratégie d'optimisation.
Lors de la phase 40, les signaux à CM issus des capteurs 2-ι à 2N sont filtrés à l'aide des filtres d'encodage décrits par les paramètres FD. Chaque coefficient pl>m(t) délivré est déduit à partir des signaux c- à c^ par l'application des filtres de la façon suivante :
où Phm(f) est la transformée de Fourier de plm(ή et Cn(f) est la transformée de Fourier de c„(t).
Dans ce mode de réalisation, les coefficients pl m(ή sont déterminés au moyen de méthodes de filtrage dans le domaine fréquentiel, comme par exemple des techniques de convolution par bloc.
La représentation du champ acoustique prend donc en compte la position des capteurs et les paramètres d'optimisation choisis et constitue une estimation fidèle du champ acoustique.
Sur la figure 7, on a représenté un schéma synoptique d'un dispositif adapté pour la mise en œuvre du procédé de l'invention.
Sur cette figure, un dispositif 50 de représentation du champ acoustique P est relié aux moyens d'acquisition 1 tels que décrits en référence à la figure 2.
Le dispositif 50 ou dispositif d'encodage est également relié en entrée à des moyens 60 de détermination des paramètres représentatifs des caractéristiques structurelles et/ou électro-acoustiques des moyens d'acquisition 1. Ces moyens 60 comportent notamment des moyens 62 de saisie de paramètres et des moyens 64 de calibrage qui sont adaptés pour mettre en œuvre respectivement les étapes 10 et 20 du procédé de l'invention tel que décrit précédemment. Le dispositif d'encodage 50 reçoit, des moyens 60 de détermination des paramètres, une pluralité de paramètres représentatifs des caractéristiques des moyens d'acquisition 1 répartis entre un signal CL de définition des caractéristiques structurelles et un signal CP de paramétrage des caractéristiques structurelles et/ou électro-acoustiques. Le dispositif reçoit également des paramètres relatifs aux stratégies de représentation dans un signal OS d'optimisation de la représentation.
Dans ces signaux, les paramètres se répartissent de la manière suivante :
- dans le signal CL de définition : - des paramètres x„ représentatifs de la position du capteur 2n ;
- dans le signal CP de paramétrage :
- des paramètres Hn(J) représentatifs de la réponse en fréquence du capteur 2n ;
- des paramètres dn(f) représentatifs du diagramme de directivité du capteur 2n ;
- des paramètres an(f) représentatifs de l'orientation du capteur 2n ;
- des paramètres <An(j) représentatifs de la densité spectrale de puissance du bruit de fond du capteur 2n ; et
- des paramètres R„,/,,„C représentatifs des capacités d'acquisition du capteur 2n ; et
- dans le signal OS d'optimisation :
- un paramètre μ(j) spécifiant le compromis entre la fidélité de représentation du champ acoustique et la minimisation du bruit de fond apporté par les capteurs ; - un paramètre L(f) spécifiant l'ordre auquel est conduite la représentation ; et
- un paramètre {(lk,mk)}(f) représentatif de la liste des coefficients dont on exige que la puissance soit égale à la puissance du coefficient correspondant dans le champ acoustique à représenter P. Avantageusement, ce dispositif 50 comporte des moyens 51 de mise en forme des signaux d'entrée adaptés pour délivrer à partir des signaux Ci à CN, un signal SI mis en forme correspondant.
Par exemple, les moyens 51 comportent des convertisseurs analogi- que-numérique, des amplificateurs ou encore des systèmes de filtrage.
Le dispositif 50 comporte en outre des moyens 52 de détermination des filtres d'encodage qui comportent un module 55 de calcul de la matrice B d'échantillonnage, un module 56 de calcul de la matrice A d'intercorrélation, lesquels sont tous deux reliés à un module 57 de calcul de la matrice E(j) d'encodage.
Cette matrice E(f) d'encodage est utilisée par un module 58 de détermination de filtres d'encodage qui délivre un signal SFD qui contient les paramètres FD représentatifs des filtres d'encodage.
Ce signal SFD est utilisé par un module 59 de traitement qui applique les filtres d'encodage au signal SI afin de délivrer un signal SIFB qui comporte les coefficients de Fourier-Bessel représentatifs du champ acoustique P.
Eventuellement, le dispositif 50 comporte une mémoire non volatile dans laquelle sont mémorisés les paramètres qui constituent le signal SFD qui ont été déterminés préalablement. Par exemple, les moyens d'acquisition 1 sont testés et calibrés par leur constructeur afin de fournir directement une mémoire comportant l'ensemble des paramètres du signal SFD qu'il convient d'intégrer dans un dispositif d'encodage afin de réaliser l'acquisition du champ acoustique P et de délivrer une représentation fidèle de ce dernier. De même, en variante, cette mémoire ne comporte que les matrices B et éventuellement A et le dispositif 50 comporte des moyens de saisie des paramètres constituant le signal OS d'optimisation afin de mettre en œuvre la détermination de la matrice E(f) d'encodage et la détermination des paramètres FD représentatifs des filtres d'encodage. Bien entendu, d'autres répartitions entre les différents modules décrits peuvent être envisagées en fonction des besoins.

Claims

REVENDICATIONS
1. Procédé de représentation d'un champ acoustique comprenant une étape d'acquisition de signaux de mesure (cn) délivrés par des moyens d'acquisition (1) formés d'un ou de plusieurs capteurs (2n) élémentaires exposés audit champ acoustique (P), caractérisé en ce qu'il comporte :
- une étape (30) de détermination de filtres d'encodage représentatifs de caractéristiques au moins structurelles desdits moyens d'acquisition (1) ; et
- une étape (40) de traitement desdits signaux de mesure (cn) par l'application desdits filtres d'encodage à ces signaux (cn) pour déterminer un nombre fini de coefficients représentatifs dans le temps et dans les trois dimensions de l'espace dudit champ acoustique (P), lesdits coefficients permettant d'obtenir une représentation dudit champ acoustique (P) sensiblement indépendante des caractéristiques desdits moyens d'acquisition (1).
2. Procédé selon la revendication 1, caractérisé en ce que lesdites ca- racteristiques structurelles comportent au moins des caractéristiques de position desdits capteurs élémentaires (2n) par rapport à un point de référence (4) prédéterminé desdits moyens d'acquisition (1).
3. Procédé selon l'une quelconque des revendications 1 ou 2, caractérisé en ce que lesdits filtres d'encodage sont en outre représentatifs de caracté- ristiques électro-acoustiques desdits moyens d'acquisition (1).
4. Procédé selon la revendication 3, caractérisé en ce que lesdites caractéristiques électro-acoustiques comportent au moins des caractéristiques liées aux capacités électro-acoustiques d'acquisition intrinsèques desdits capteurs élémentaires (2n).
5. Procédé selon l'une quelconque des revendications 1 à 4, caractérisé en ce que lesdits coefficients permettant d'obtenir une représentation du champ acoustique (P) sont des coefficients dits de Fourier-Bessel et/ou des combinaisons linéaires de coefficients de Fourier-Bessel.
6. Procédé selon l'une quelconque des revendications 1 à 5, caractéri- se en ce que ladite étape (30) de détermination des filtres d'encodage comprend :
- une sous-étape (32) de détermination d'une matrice (B) d'échantillonnage représentative des capacités d'acquisition desdits moyens d'acquisition (1) ; 03/096742
27
- une sous-étape (34) de détermination d'une matrice (A) d'intercorrélation représentative de la ressemblance entre lesdits signaux de mesure (cn) délivrés par les capteurs (2n) élémentaires formant lesdits moyens d'acquisition (1) ; et - une sous-étape (36) de détermination d'une matrice (E(f) ; E) d'encodage à partir de ladite matrice (B) d'échantillonnage, de ladite matrice (A) d'intercorrélation, et d'un paramètre (μ(f)) représentatif d'un compromis souhaité entre la fidélité de représentation du champ acoustique et la minimisation du bruit de fond induit par les moyens d'acquisition (1), laquelle matrice est représentative desdits filtres d'encodage.
7. Procédé selon la revendication 6, caractérisé en ce que lesdites sous-étapes de détermination des matrices sont réalisées pour un nombre fini de fréquences de fonctionnement.
8. Procédé selon l'une quelconque des revendications 6 ou 7, caracté- risé en ce que ladite étape (32) de détermination de la matrice (B) d'échantillonnage est réalisée, pour chacun desdits capteurs élémentaires (2n) formant lesdits moyens d'acquisition (1 ), à partir :
- de paramètres (χn ) représentatifs de la position dudit capteur (2n) par rapport au centre (4) desdits moyens d'acquisition (1) ; et/ou - d'un nombre fini de coefficients (Bn.ι,m(O) représentatifs des capacités d'acquisition dudit capteur (2n).
9. Procédé selon la revendication 8, caractérisé en ce que ladite étape de détermination de la matrice (B) d'échantillonnage est réalisée en outre à partir au moins d'un des paramètres parmi : - des paramètres (Hn(f)) représentatifs des réponses en fréquence de tout ou partie des capteurs (2n) ;
- des paramètres (dn(f)) représentatifs des diagrammes de directivité de tout ou partie des capteurs (2n) ;
- des paramètres (αn(f)) représentatifs des orientations de tout ou par- tie des capteurs (2n), à savoir de leur direction de sensibilité maximale ;
- des paramètres (σ2 n(f)) représentatifs des densités spectrales de puissance du bruit de fond de tout ou partie des capteurs (2n) ;
- d'un paramètre (L(f)) spécifiant l'ordre auquel est conduite la représentation ; et - d'un paramètre ({(lk,mk)}(f)) représentatif d'une liste de coefficients dont on exige que la puissance soit égale à la puissance du coefficient correspondant dans le champ acoustique à représenter (P).
10. Procédé selon l'une quelconque des revendications 1 à 9, caracté- risé en ce qu'il comporte une étape (20) de calibrage permettant de délivrer tout ou partie des paramètres utilisés dans ladite étape (30) de détermination des filtres d'encodage.
11. Procédé selon la revendication 10, caractérisé en ce que ladite étape (20) de calibrage comporte, pour au moins l'un desdits capteurs élémentai- res (2n) formant lesdits moyens d'acquisition (1 ) :
- une sous-étape (24) d'acquisition de signaux représentatifs des capacités d'acquisition dudit au moins un capteur (2n); et
- une sous-étape (28) de détermination de paramètres représentatifs de caractéristiques électro-acoustiques et ou structurelles dudit au moins un cap- teur (2n).
12. Procédé selon la revendication 11 , caractérisé en ce que ladite étape (20) de calibrage comporte en outre :
- une sous-étape (22) d'émission d'un champ acoustique spécifique vers ledit au moins un capteur (2π), ladite sous-étape (24) d'acquisition corres- pondant à l'acquisition des signaux délivrés par ce capteur (2n) lorsque exposé audit champ acoustique spécifique ; et
- une sous-étape (26) de modélisation dudit champ acoustique spécifique en un nombre fini de coefficients afin de permettre la réalisation de ladite sous-étape (28) de détermination de paramètres représentatifs de caractéristi- ques électro-acoustiques et/ou structurelles du capteur (2n).
13. Procédé selon l'une quelconque des revendications 10 à 12, caractérisé en ce que ladite étape (20) de calibrage comporte une sous-étape de réception d'un nombre fini de signaux représentatifs des caractéristiques électroacoustiques et structurelles desdits capteurs (2n) formant lesdits moyens d'acquisition (1), lesquels signaux sont directement utilisés lors de ladite sous- étape de détermination des caractéristiques électro-acoustiques et/ou structurelles desdits moyens d'acquisition (1).
14. Procédé selon l'une quelconque des revendications 1 à 13, caractérisé en ce qu'il comporte une étape (10) de saisie permettant de déterminer tout ou partie des paramètres utilisés lors de ladite étape (30) de détermination des filtres d'encodage.
15. Programme d'ordinateur comprenant des instructions de code de programme pour l'exécution des étapes du procédé selon l'une quelconque des revendications 1 à 14, lorsque ledit programme est exécuté sur un ordinateur.
16. Support mobile du type comportant au moins un processeur de traitement et un élément de mémoire non volatile, caractérisé en ce que ladite mémoire comprend un programme comportant des instructions de code pour l'exécution des étapes du procédé selon l'une quelconque des revendications 1 à 14 lorsque ledit processeur exécute ledit programme.
17. Dispositif de représentation d'un champ acoustique connectable à des moyens d'acquisition (1) formés d'un ou plusieurs capteurs (2n) élémentaires délivrant des signaux de mesure (cn) lorsqu'ils sont exposés audit champ acoustique (P), caractérisé en ce qu'il comporte un module (59) de traitement des si- gnaux de mesure (cn) par l'application de filtres d'encodage représentatifs de caractéristiques au moins structurelles desdits moyens d'acquisition (1) à ces signaux de mesure (cn) pour délivrer un signal (SIFB) qui comporte un nombre fini de coefficients représentatifs dans le temps et dans les trois dimensions de l'espace dudit champ acoustique (P), lesdits coefficients permettant d'obtenir une représentation dudit champ acoustique (P) sensiblement indépendante des caractéristiques desdits moyens d'acquisition (1 ).
18. Dispositif selon la revendication 17, caractérisé en ce que lesdits filtres d'encodage sont en outre représentatifs de caractéristiques électroacoustiques desdits moyens d'acquisition (1).
19. Dispositif selon l'une quelconque des revendications 17 ou 18, caractérisé en ce qu'il comporte en outre des moyens (52) de détermination desdits filtres d'encodage représentatifs de caractéristiques structurelles et/ou électroacoustiques desdits moyens d'acquisition (1).
20. Dispositif selon la revendication 19, caractérisé en ce que lesdits moyens (52) de détermination de filtres d'encodage reçoivent en entrée au moins l'un des paramètres parmi les paramètres suivants :
- des paramètres (x„ ) représentatifs des positions par rapport au centre desdits moyens d'acquisition (1) de tout ou partie des capteurs (2π) ; - un nombre fini de coefficients (Br,.ι,m(f)) représentatifs des capacités d'acquisition de tout ou partie des capteurs (2n) ;
- des paramètres (Hn(f)) représentatifs des réponses en fréquence de tout ou partie des capteurs (2n); - des paramètres (dn(f)) représentatifs des diagrammes de directivité de tout ou partie des capteurs (2n) ;
- des paramètres (αn(f)) représentatifs des orientations de tout ou partie des capteurs (2n), à savoir de leur direction de sensibilité maximaie ;
- des paramètres (σ2 n(f)) représentatifs des densités spectrales de puissance du bruit de fond de tout ou partie des capteurs (2n);
- un paramètre μ(f) représentatif du compromis souhaité entre la fidélité de représentation du champ acoustique et la minimisation du bruit de fond induit par les moyens d'acquisition (1 )
- un paramètre (L(f)) spécifiant l'ordre auquel est conduit l'encodage ; et
- un paramètre ({(lk,ι ik)}(f)) représentatif d'une liste de coefficients dont on exige que la puissance soit égale à la puissance du coefficient correspondant dans le champ acoustique à représenter (P).
21. Dispositif selon la revendication 20, caractérisé en ce qu'il est as- socié à des moyens (60) de détermination de tout ou partie des paramètres reçus par lesdits moyens (52) de détermination des filtres d'encodage, lesdits moyens (60) comportant au moins l'un des éléments suivants :
- des moyens (62) de saisie des paramètres ; et/ou
- des moyens (64) de calibrage.
22. Dispositif selon l'une quelconque des revendications 17 à 21 , caractérisé en ce qu'il est associé à des moyens (51) de mise en forme desdits signaux de mesure (ci à CN) afin de délivrer un signal (SI) mis en forme correspondant.
EP03749929A 2002-05-07 2003-05-06 Procede et systeme de representation d un champ acoustique Expired - Lifetime EP1502475B8 (fr)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
FR0205741 2002-05-07
FR0205741A FR2839565B1 (fr) 2002-05-07 2002-05-07 Procede et systeme de representation d'un champ acoustique
PCT/FR2003/001410 WO2003096742A1 (fr) 2002-05-07 2003-05-06 Procede et systeme de representation d'un champ acoustique

Publications (3)

Publication Number Publication Date
EP1502475A1 true EP1502475A1 (fr) 2005-02-02
EP1502475B1 EP1502475B1 (fr) 2005-07-27
EP1502475B8 EP1502475B8 (fr) 2005-09-28

Family

ID=29286370

Family Applications (1)

Application Number Title Priority Date Filing Date
EP03749929A Expired - Lifetime EP1502475B8 (fr) 2002-05-07 2003-05-06 Procede et systeme de representation d un champ acoustique

Country Status (12)

Country Link
US (1) US7212956B2 (fr)
EP (1) EP1502475B8 (fr)
JP (1) JP4293986B2 (fr)
KR (1) KR100972419B1 (fr)
CN (1) CN1659926B (fr)
AT (1) ATE300852T1 (fr)
AU (1) AU2003255562B2 (fr)
CA (1) CA2484588C (fr)
DE (1) DE60301146T2 (fr)
DK (1) DK1502475T3 (fr)
FR (1) FR2839565B1 (fr)
WO (1) WO2003096742A1 (fr)

Families Citing this family (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8073157B2 (en) * 2003-08-27 2011-12-06 Sony Computer Entertainment Inc. Methods and apparatus for targeted sound detection and characterization
US7809145B2 (en) * 2006-05-04 2010-10-05 Sony Computer Entertainment Inc. Ultra small microphone array
US7783061B2 (en) * 2003-08-27 2010-08-24 Sony Computer Entertainment Inc. Methods and apparatus for the targeted sound detection
US8947347B2 (en) * 2003-08-27 2015-02-03 Sony Computer Entertainment Inc. Controlling actions in a video game unit
US9174119B2 (en) 2002-07-27 2015-11-03 Sony Computer Entertainement America, LLC Controller for providing inputs to control execution of a program when inputs are combined
US8233642B2 (en) 2003-08-27 2012-07-31 Sony Computer Entertainment Inc. Methods and apparatuses for capturing an audio signal based on a location of the signal
US7803050B2 (en) 2002-07-27 2010-09-28 Sony Computer Entertainment Inc. Tracking device with sound emitter for use in obtaining information for controlling game program execution
US8139793B2 (en) * 2003-08-27 2012-03-20 Sony Computer Entertainment Inc. Methods and apparatus for capturing audio signals based on a visual image
US8160269B2 (en) 2003-08-27 2012-04-17 Sony Computer Entertainment Inc. Methods and apparatuses for adjusting a listening area for capturing sounds
US20070223732A1 (en) * 2003-08-27 2007-09-27 Mao Xiao D Methods and apparatuses for adjusting a visual image based on an audio signal
GB0523946D0 (en) * 2005-11-24 2006-01-04 King S College London Audio signal processing method and system
US20140167972A1 (en) * 2012-12-13 2014-06-19 General Electric Company Acoustically-responsive optical data acquisition system for sensor data
EP2765791A1 (fr) 2013-02-08 2014-08-13 Thomson Licensing Procédé et appareil pour déterminer des directions de sources sonores non corrélées dans une représentation d'ambiophonie d'ordre supérieur d'un champ sonore
CN104935913B (zh) * 2014-03-21 2018-12-04 杜比实验室特许公司 处理多个装置采集的音频或视频信号
CN105898668A (zh) * 2016-03-18 2016-08-24 南京青衿信息科技有限公司 一种声场空间的坐标定义方法

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4060850A (en) * 1977-04-25 1977-11-29 The United States Of America As Represented By The Secretary Of The Navy Beam former using bessel sequences
JPH0728470B2 (ja) * 1989-02-03 1995-03-29 松下電器産業株式会社 アレイマイクロホン
US5216640A (en) * 1992-09-28 1993-06-01 The United States Of America As Represented By The Secretary Of The Navy Inverse beamforming sonar system and method
US6201872B1 (en) * 1995-03-12 2001-03-13 Hersh Acoustical Engineering, Inc. Active control source cancellation and active control Helmholtz resonator absorption of axial fan rotor-stator interaction noise
US6216540B1 (en) * 1995-06-06 2001-04-17 Robert S. Nelson High resolution device and method for imaging concealed objects within an obscuring medium
US7348181B2 (en) * 1997-10-06 2008-03-25 Trustees Of Tufts College Self-encoding sensor with microspheres
JP3584800B2 (ja) * 1999-08-17 2004-11-04 ヤマハ株式会社 音場再現方法およびその装置

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
See references of WO03096742A1 *

Also Published As

Publication number Publication date
DE60301146D1 (de) 2005-09-01
US7212956B2 (en) 2007-05-01
EP1502475B8 (fr) 2005-09-28
CA2484588C (fr) 2013-03-12
US20050177606A1 (en) 2005-08-11
ATE300852T1 (de) 2005-08-15
JP4293986B2 (ja) 2009-07-08
FR2839565A1 (fr) 2003-11-14
CA2484588A1 (fr) 2003-11-20
JP2005531016A (ja) 2005-10-13
CN1659926A (zh) 2005-08-24
KR100972419B1 (ko) 2010-07-27
AU2003255562B2 (en) 2009-04-23
KR20050010784A (ko) 2005-01-28
WO2003096742A1 (fr) 2003-11-20
FR2839565B1 (fr) 2004-11-19
EP1502475B1 (fr) 2005-07-27
DE60301146T2 (de) 2006-06-01
CN1659926B (zh) 2010-05-12
DK1502475T3 (da) 2005-11-28
AU2003255562A1 (en) 2003-11-11

Similar Documents

Publication Publication Date Title
EP1546916B1 (fr) Procede et systeme de traitement d&#39;une representation d&#39;un champ acoustique
EP1479266B1 (fr) Procede et dispositif de pilotage d&#39;un ensemble de restitution d&#39;un champ acoustique
EP1586220B1 (fr) Procede et dispositif de pilotage d&#39;un ensemble de restitution a partir d&#39;un signal multicanal
EP1502475A1 (fr) Procede et systeme de representation d un champ acoustique
EP3348079B1 (fr) Procédé et système d&#39;élaboration d&#39;une fonction de transfert relative à la tête adaptée à un individu
EP1563485B1 (fr) Procede de traitement de donnees sonores et dispositif d&#39;acquisition sonore mettant en oeuvre ce procede
EP1836876B1 (fr) Procédé et dispositif d&#39;individualisation de hrtfs par modélisation
EP2258119B1 (fr) Procede et dispositif pour la determination de fonctions de transfert de type hrtf
EP2898707A1 (fr) Calibration optimisee d&#39;un systeme de restitution sonore multi haut-parleurs
EP1946612A1 (fr) Individualisation de hrtfs utilisant une modelisation par elements finis couplee a un modele correctif
EP1652406B1 (fr) Systeme et procede de determination d&#39;une representation d&#39;un champ acoustique
EP3895446B1 (fr) Procede d&#39;interpolation d&#39;un champ sonore, produit programme d&#39;ordinateur et dispositif correspondants.
WO2018050292A1 (fr) Dispositif et procede de captation et traitement d&#39;un champ acoustique tridimensionnel
EP3384688B1 (fr) Décompositions successives de filtres audio
EP3934282A1 (fr) Procédé de conversion d&#39;un premier ensemble de signaux représentatifs d&#39;un champ sonore en un second ensemble de signaux et dispositif électronique associé
WO2022184800A1 (fr) Procédé de détermination d&#39;une direction de propagation d&#39;une source sonore par création de signaux sinusoïdaux à partir des signaux sonores reçus par des microphones
FR3072820A1 (fr) Procede de modelisation et de simulation en temps reel de sources acoustiques dans un environnement par decomposition en sources elementaires
FR3107982A1 (fr) Procédé et système de détermination de filtres d’égalisation sonore d’un système audio
Otshudi et al. Synthesis of principal techniques of measuring impulse response in room acoustics
WO2014102199A1 (fr) Dispositif et procede d&#39;interpolation spatiale de sons

Legal Events

Date Code Title Description
PUAI Public reference made under article 153(3) epc to a published international application that has entered the european phase

Free format text: ORIGINAL CODE: 0009012

GRAP Despatch of communication of intention to grant a patent

Free format text: ORIGINAL CODE: EPIDOSNIGR1

17P Request for examination filed

Effective date: 20041101

AK Designated contracting states

Kind code of ref document: A1

Designated state(s): AT BE BG CH CY CZ DE DK EE ES FI FR GB GR HU IE IT LI LU MC NL PT RO SE SI SK TR

GRAS Grant fee paid

Free format text: ORIGINAL CODE: EPIDOSNIGR3

GRAA (expected) grant

Free format text: ORIGINAL CODE: 0009210

AK Designated contracting states

Kind code of ref document: B1

Designated state(s): AT BE BG CH CY CZ DE DK EE ES FI FR GB GR HU IE IT LI LU MC NL PT RO SE SI SK TR

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: IT

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT;WARNING: LAPSES OF ITALIAN PATENTS WITH EFFECTIVE DATE BEFORE 2007 MAY HAVE OCCURRED AT ANY TIME BEFORE 2007. THE CORRECT EFFECTIVE DATE MAY BE DIFFERENT FROM THE ONE RECORDED.

Effective date: 20050727

Ref country code: EE

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20050727

Ref country code: SK

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20050727

Ref country code: FI

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20050727

Ref country code: IE

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20050727

Ref country code: SI

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20050727

Ref country code: RO

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20050727

Ref country code: TR

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20050727

Ref country code: AT

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20050727

Ref country code: NL

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20050727

Ref country code: CZ

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20050727

REG Reference to a national code

Ref country code: GB

Ref legal event code: FG4D

Free format text: NOT ENGLISH

REG Reference to a national code

Ref country code: CH

Ref legal event code: EP

REG Reference to a national code

Ref country code: IE

Ref legal event code: FG4D

Free format text: LANGUAGE OF EP DOCUMENT: FRENCH

REF Corresponds to:

Ref document number: 60301146

Country of ref document: DE

Date of ref document: 20050901

Kind code of ref document: P

RAP2 Party data changed (patent owner data changed or rights of a patent transferred)

Owner name: TRINNOV AUDIO

RIN2 Information on inventor provided after grant (corrected)

Inventor name: MONTOYA, SEBASTIAN

Inventor name: BRUNO, REMY

Inventor name: LABORIE, ARNAUD

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: SE

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20051027

Ref country code: BG

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20051027

Ref country code: GR

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20051027

NLT2 Nl: modifications (of names), taken from the european patent patent bulletin

Owner name: TRINNOV AUDIO

Effective date: 20050914

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: ES

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20051107

REG Reference to a national code

Ref country code: DK

Ref legal event code: T3

GBT Gb: translation of ep patent filed (gb section 77(6)(a)/1977)

Effective date: 20051104

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: PT

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20051227

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: HU

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20060128

NLV1 Nl: lapsed or annulled due to failure to fulfill the requirements of art. 29p and 29m of the patents act
REG Reference to a national code

Ref country code: IE

Ref legal event code: FD4D

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: MC

Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES

Effective date: 20060531

PLBE No opposition filed within time limit

Free format text: ORIGINAL CODE: 0009261

STAA Information on the status of an ep patent application or granted ep patent

Free format text: STATUS: NO OPPOSITION FILED WITHIN TIME LIMIT

26N No opposition filed

Effective date: 20060428

PGFP Annual fee paid to national office [announced via postgrant information from national office to epo]

Ref country code: BE

Payment date: 20060712

Year of fee payment: 4

BERE Be: lapsed

Owner name: MONTOYA, SEBASTIAN

Effective date: 20070531

Owner name: LABORIE, ARNAUD

Effective date: 20070531

Owner name: BRUNO, REMY

Effective date: 20070531

REG Reference to a national code

Ref country code: CH

Ref legal event code: PL

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: LI

Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES

Effective date: 20070531

Ref country code: CH

Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES

Effective date: 20070531

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: BE

Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES

Effective date: 20070531

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: LU

Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES

Effective date: 20060506

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: CY

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20050727

REG Reference to a national code

Ref country code: FR

Ref legal event code: PLFP

Year of fee payment: 14

REG Reference to a national code

Ref country code: FR

Ref legal event code: PLFP

Year of fee payment: 15

REG Reference to a national code

Ref country code: FR

Ref legal event code: PLFP

Year of fee payment: 16

PGFP Annual fee paid to national office [announced via postgrant information from national office to epo]

Ref country code: GB

Payment date: 20220519

Year of fee payment: 20

Ref country code: FR

Payment date: 20220411

Year of fee payment: 20

Ref country code: DK

Payment date: 20220426

Year of fee payment: 20

Ref country code: DE

Payment date: 20220511

Year of fee payment: 20

REG Reference to a national code

Ref country code: DE

Ref legal event code: R071

Ref document number: 60301146

Country of ref document: DE

REG Reference to a national code

Ref country code: DK

Ref legal event code: EUP

Expiry date: 20230506

REG Reference to a national code

Ref country code: GB

Ref legal event code: PE20

Expiry date: 20230505

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: GB

Free format text: LAPSE BECAUSE OF EXPIRATION OF PROTECTION

Effective date: 20230505