EP1733381B1 - System for actively reducing sound - Google Patents

System for actively reducing sound Download PDF

Info

Publication number
EP1733381B1
EP1733381B1 EP05729401A EP05729401A EP1733381B1 EP 1733381 B1 EP1733381 B1 EP 1733381B1 EP 05729401 A EP05729401 A EP 05729401A EP 05729401 A EP05729401 A EP 05729401A EP 1733381 B1 EP1733381 B1 EP 1733381B1
Authority
EP
European Patent Office
Prior art keywords
array
loudspeakers
microphone
microphones
error signal
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.)
Not-in-force
Application number
EP05729401A
Other languages
German (de)
French (fr)
Other versions
EP1733381A1 (en
Inventor
Arthur Perry Berkhoff
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.)
Nederlandse Organisatie voor Toegepast Natuurwetenschappelijk Onderzoek TNO
Original Assignee
Nederlandse Organisatie voor Toegepast Natuurwetenschappelijk Onderzoek TNO
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 Nederlandse Organisatie voor Toegepast Natuurwetenschappelijk Onderzoek TNO filed Critical Nederlandse Organisatie voor Toegepast Natuurwetenschappelijk Onderzoek TNO
Priority to EP05729401A priority Critical patent/EP1733381B1/en
Publication of EP1733381A1 publication Critical patent/EP1733381A1/en
Application granted granted Critical
Publication of EP1733381B1 publication Critical patent/EP1733381B1/en
Not-in-force legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G10MUSICAL INSTRUMENTS; ACOUSTICS
    • G10KSOUND-PRODUCING DEVICES; METHODS OR DEVICES FOR PROTECTING AGAINST, OR FOR DAMPING, NOISE OR OTHER ACOUSTIC WAVES IN GENERAL; ACOUSTICS NOT OTHERWISE PROVIDED FOR
    • G10K11/00Methods or devices for transmitting, conducting or directing sound in general; Methods or devices for protecting against, or for damping, noise or other acoustic waves in general
    • G10K11/16Methods or devices for protecting against, or for damping, noise or other acoustic waves in general
    • G10K11/175Methods or devices for protecting against, or for damping, noise or other acoustic waves in general using interference effects; Masking sound
    • G10K11/178Methods or devices for protecting against, or for damping, noise or other acoustic waves in general using interference effects; Masking sound by electro-acoustically regenerating the original acoustic waves in anti-phase
    • G10K11/1785Methods, e.g. algorithms; Devices
    • G10K11/17857Geometric disposition, e.g. placement of microphones
    • GPHYSICS
    • G10MUSICAL INSTRUMENTS; ACOUSTICS
    • G10KSOUND-PRODUCING DEVICES; METHODS OR DEVICES FOR PROTECTING AGAINST, OR FOR DAMPING, NOISE OR OTHER ACOUSTIC WAVES IN GENERAL; ACOUSTICS NOT OTHERWISE PROVIDED FOR
    • G10K11/00Methods or devices for transmitting, conducting or directing sound in general; Methods or devices for protecting against, or for damping, noise or other acoustic waves in general
    • G10K11/16Methods or devices for protecting against, or for damping, noise or other acoustic waves in general
    • G10K11/175Methods or devices for protecting against, or for damping, noise or other acoustic waves in general using interference effects; Masking sound
    • G10K11/178Methods or devices for protecting against, or for damping, noise or other acoustic waves in general using interference effects; Masking sound by electro-acoustically regenerating the original acoustic waves in anti-phase
    • G10K11/1781Methods or devices for protecting against, or for damping, noise or other acoustic waves in general using interference effects; Masking sound by electro-acoustically regenerating the original acoustic waves in anti-phase characterised by the analysis of input or output signals, e.g. frequency range, modes, transfer functions
    • G10K11/17821Methods or devices for protecting against, or for damping, noise or other acoustic waves in general using interference effects; Masking sound by electro-acoustically regenerating the original acoustic waves in anti-phase characterised by the analysis of input or output signals, e.g. frequency range, modes, transfer functions characterised by the analysis of the input signals only
    • G10K11/17825Error signals
    • GPHYSICS
    • G10MUSICAL INSTRUMENTS; ACOUSTICS
    • G10KSOUND-PRODUCING DEVICES; METHODS OR DEVICES FOR PROTECTING AGAINST, OR FOR DAMPING, NOISE OR OTHER ACOUSTIC WAVES IN GENERAL; ACOUSTICS NOT OTHERWISE PROVIDED FOR
    • G10K11/00Methods or devices for transmitting, conducting or directing sound in general; Methods or devices for protecting against, or for damping, noise or other acoustic waves in general
    • G10K11/16Methods or devices for protecting against, or for damping, noise or other acoustic waves in general
    • G10K11/175Methods or devices for protecting against, or for damping, noise or other acoustic waves in general using interference effects; Masking sound
    • G10K11/178Methods or devices for protecting against, or for damping, noise or other acoustic waves in general using interference effects; Masking sound by electro-acoustically regenerating the original acoustic waves in anti-phase
    • G10K11/1787General system configurations
    • G10K11/17879General system configurations using both a reference signal and an error signal
    • G10K11/17881General system configurations using both a reference signal and an error signal the reference signal being an acoustic signal, e.g. recorded with a microphone

Definitions

  • the invention relates to a system for actively reducing sound from a primary noise source, such as traffic noise, comprising: a loudspeaker connector for connecting to at least one loudspeaker for generating anti-sound for reducing said noisy sound; a microphone connector for connecting to at least a first microphone placed adjacent to said loudspeaker; a control unit coupled to said first microphone connector, for providing an error signal, based on the output of said first microphone; and a control unit for outputting a signal to said loudspeaker connector, for controlling said loudspeaker based on said error signal of said control unit.
  • a primary noise source such as traffic noise
  • the invention has as one of its goals to provide a system that at the one time offers a compact and robust way of implementing a noise barrier but where the hereabove described disadvantages are mitigated.
  • D1 is indeed similar to the subject matter of the invention. However, D1 assumes that the primary source near the real microphone has identical characteristics as the source near the virtual microphone.
  • a simulated error signal of a virtual microphone is a result of a combination of a first transfer function expressing the far field character of the primary source near said first microphone and a second transfer function expressing the near field character of the secondary source near said first microphone.
  • the noise suppression is much more successful than conventional active noise barriers and offers reductions even up to 10 dB or more compared to a sound reduction of a conventional system with adjacently placed microphones.
  • the primary and secondary sources may have a very different nature they may lead to rather different requirements for a transfer function.
  • One big difference, for example, is that the near field sensors are located rather far away from the primary sources while the near field sensors are, by definition, in the near field of the secondary sources.
  • the loudspeaker may be arranged to produce anti-sound in a direction away from said primary noise source. Further, said microphone is preferably placed adjacent in front of said loudspeaker.
  • a second microphone may be placed between said loudspeaker and said primary source.
  • Said second microphone may be arranged to provide a second error signal that is used as a feed forward error signal in order to achieve an expected sound pressure level of the primary source near said loudspeaker.
  • Said at least one loudspeaker and said at least first microphone may be part of an array of loudspeakers and an array of microphones respectively, wherein loudspeakers and said microphones are placed at a distance less than 5 times an interspacing between two adjacent loudspeakers.
  • Such an array may serve as an active noise barrier on the side of a road etc.
  • said at least one loudspeaker and said at least one microphone are part of an array of loudspeakers and an array of microphones respectively, wherein said loudspeakers and said microphones are placed relative to each other in a range between 10% and 100% of an interspacing between two adjacent loudspeakers. It has been found that a position of said microphone at such a general "close" distance of the loudspeaker, is still able to derive a sufficient error signal in order to be able to calculate a far field sound pressure.
  • said microphone is formed integral with said loudspeaker in a panel to be placed on the side of a road.
  • the invention is also related to a panel for a noise screen to be placed on the side of the road, for actively reducing noise from a primary source, comprising: a loudspeaker to be directed away from the primary source a microphone attached to said panel and placed on a distance away from said panel; and a system according to any of the preceding aspects.
  • said panel may further comprise a second microphone placed opposite to said first microphone, viewed in a direction away from said panel. Furthermore, said loudspeaker is placed on top of said panel. In this way the influence of wind is more or less the same for the primary sources and the secondary sources. This can lead to a more robust system with respect to fluctuations of the wind.
  • the invention relates to providing a a system according to any of the preceding aspects, further comprising an input for inputting data representing environmental conditions such as wind speed and/or temperature, said processor arranged to provide a transfer function that is dynamically adapted to said environmental conditions.
  • the invention is related to a method as described in claim 12 and 13.
  • FIG. 1 A proposed configuration of the system according to the invention can be found in Figure 1 for an active noise barrier 6, for instance, for blocking primary traffic noise of passing vehicles 1.
  • reference sensors 2 microphones
  • near secondary sources 3 ladspeakers
  • nearfield error sensors 4 microphones
  • farfield error sensors 5 can be used-
  • One section (n) of the active noise barrier 6 is considered.
  • the sensors and sources of adjacent sections should be taken into account.
  • virtual error signals are derived from the nearfield error sensors. These virtual error signals should represent the degrees of freedom in the farfield of the angular sector of interest.
  • the primary disturbance d reference signals x
  • actuator control signals u nearfield error signals e _ y
  • Lower case variables denote vectorial quantities that are a function of time
  • upper case variables denote matrice s that operate on the time-dependent vectors.
  • the assumption is that G ux is compensated in the controller, leading to a control structure based on Internal Model Control. Furthermore, aspects related to decentralization are not taken into account yet but will be discussed in the remainder.
  • Figure 2 shows a system using a single transfer function H yz . Such systems are used for actively controlling noise from panels. However, the performance for the present application is insufficient, as can be demonstrated in simulations.
  • Figure 3 shows a system where the basic strategy is to minimize the nearfield pressure, using a feedforward controller, as in Figure 2 .
  • Figure 4 shows a system in which the nearfield pressure is directly fed back to the controller, where the basic strategy is also to minimize the nearfield pressure. The latter configuration has been tested by Japanese researchers. In the present application, reduction of the nearfield pressure is very different from minimizing the farfield pressure.
  • the simulations described below will show that the configuration of Figure 5 performs substantially better than the configurations of Figure 2 - Figure 4 .
  • FIG. 5 in a first approach for the reduction of farfield sound pressure the sound pressure is reduced near the noise barrier.
  • This is the approach that is used in, for example, a Japanese active noise barrier as illustrated in Figure 4 .
  • a disadvantage of this system is that a reduction of the nearfield pressure not necessarily leads to a reduction of the farfield pressure. This holds for both feedback configurations and for feedforward configurations.
  • an improved technique is based on the estimation of the farfield error signals e_z (which are assumed to be unavailable when the system is operational) from the nearfield error signals e_y.
  • a transfer function H yz is designed which makes a real-time estimate of e_z from measured data e_y.
  • H yz the transfer function
  • the nearfield sensors are located rather far away from the primary sources while the nearfield sensors are, by definition, in the nearfield of the secondary sources. So, although a single transfer function H yz may have acceptable performance in some systems, such as active panels, the performance for active noise barriers can be improved further.
  • Stability robustness and performance robustness are governed by the sensitivity of the control system to changes in the secondary path.
  • the propagation distances involved in the secondary path are considerably reduced if nearfield sensors are used instead of the farfield sensors. Therefore, it is expected that the reduced change of the propagation path will have a positive effect on the stability robustness of the system.
  • the primary sources are effectively diffraction sources on top of a noise barrier and if the secondary sources are close to these diffraction sources then we obtain a system in which a large part of the possibly varying transmission path is common to both the primary signals and the secondary signals. This results in an expected improved performance robustness with respect to changes in wind, temperature, etc.
  • an initial calibration phase is assumed in which sources are placed at random positions in the region where the primary noise sources are to be expected.
  • the estimate of the primary signal on one of the farfield sensors is shown in Figure 8 .
  • the estimate of one of the secondary source transfer functions to the farfield as estimated from the secondary transfer functions to the nearfield sensors is shown in Figure 9 .
  • calibration sources have been used, which are indicated as green asterixes in Figure 12 .
  • Three independent primary sources are used, also shown in Figure 12 .
  • the identification procedure is based on the solution of the equations resulting from an assumed multi-input, multi-output Finite Impulse Response model leading to a block-Toeplitz structure.
  • This value of ⁇ was obtained by reducing the value of ⁇ until the prediction error for primary signal validation data did not reduce anymore.
  • the latter prediction error was obtained with validation data from an independent set of primary sources at random positions in the same region as for the solution of the system of equations.
  • the system can be adapted to conditions such as wind speed and/or temperature.
  • the transfer function is dynamically adapted to said environmental conditions.
  • the objective is to make modifications to the J ⁇ I dimensional matrix transfer functions H yz d t and H yz s t using measurements of the wind speed vector U such that these functions more accurately predict farfield sound pressure from measured near field sound pressure.
  • a single effective wind speed vector is assumed for the whole region of interest.
  • This interpolations functions have the effect of modifying the sound of speed due to wind speed.
  • the transformation can be derived from measurements at various wind speed magnitudes and directions but also from a physical model.
  • an I ⁇ I dimensional static transfer function matrix V is defined. This matrix represents a dependance of windspeed and direction.
  • the example of Figure 10 is an active noise barrier in which the error signal equals the pressure as measured on microphones near the secondary sources.
  • this distance i.e. the difference between the z- coordinate of the 5 secondary sources and z-coordinate of the 5 error sensors is 0.5 m.
  • 5 reference sensors are used, which are positioned, as seen from the secondary sources, 1 m towards the primary sources.
  • the primary noise sources are three independent broadband noise sources which are positioned at a distance of 4.5 m from the secondary sources.
  • Figure 10 shows the configuration and the resulting sound pressure in the x,z-plane. It can be seen that sound pressure reductions are mainly obtained near the error sensors. The reduction at the error sensors is 24.2 dB but the reduction in the farfield is much less, being 1.9 dB.
  • Figure 12 shows the results in case the error signal equals virtual farfield signals, using the same microphone positions as in Figure 10 .
  • the virtual sensor signals are obtained by processing the nearfield error signals e _ y with a fixed operator which is determined in a calibration phase.
  • the sources as used in this calibration phase are positioned in the region where the primary noise sources are to be expected (see Figure 12 ). It can be seen that the resulting farfield sound pressure reductions are slightly less than in Figure 11 , being 9.0 dB. However, the farfield sound pressure reductions are considerably higher than obtained with the nearfield error sensors as used in Figure 10 .
  • the reduction of the error signals is approximately equal in case of minimizing the virtual error signals, being 11.6 dB, and in case of minimizing the true farfield error signals (11.3 dB).
  • the numerical results for the three cases are collected in Table 2. [Table 2] Error signal error signal av. reduction [dB] farfield pressure av. reduction [dB] Figure 10 24.2 1.9 Figure 11 11.3 11.3
  • the influence of a simultaneous change of the positions of the primary sources and a changing primary spectrum was investigated using a modified primary spectrum that was obtained by applying a fourth-order bandpass filter with Butterworth characteristic with lower and upper cutoff frequencies of 0.15 f s and 0.35 f s , respectively, to the primary source signals.
  • a slight improvement was found of the performance on the farfield evaluation microphones from 10.6 dB to 11.4 dB.
  • the performance for a moving source was investigated.
  • the radiation characteristics of the moving sources were simulated by taking into account the doppler shift as well as the changing radiation characteristics depending on the direction of the movement. Also, simulations were performed in order to study the performance variations for a change in the height of the primary sources.

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Acoustics & Sound (AREA)
  • Multimedia (AREA)
  • Soundproofing, Sound Blocking, And Sound Damping (AREA)

Description

  • The invention relates to a system for actively reducing sound from a primary noise source, such as traffic noise, comprising: a loudspeaker connector for connecting to at least one loudspeaker for generating anti-sound for reducing said noisy sound;
    a microphone connector for connecting to at least a first microphone placed adjacent to said loudspeaker;
    a control unit coupled to said first microphone connector, for providing an error signal, based on the output of said first microphone; and
    a control unit for outputting a signal to said loudspeaker connector, for controlling said loudspeaker based on said error signal of said control unit.
  • Recently, such systems have been considered as active noise barriers for reduction of traffic noise. The goal of the active noise barriers is to reduce the noise on positions distant from the noise barrier, such as on the facades of houses. Microphones at these positions can be used in order to create error signals to be minimized by the controller. A problem with this is the vulnerability of connections and microphones and a therefore expected short economic life. In another configuration, microphones were placed directly near the active noise barrier. This leads to obvious advantages due to the absence of sensors at inconvenient positions. This can lead to a possibly more robust system with respect to fluctuations of the wind. However, it has been found that the minimization of sound near the loudspeakers is not always appropriate for reducing the sound pressure at more distant locations.
  • The invention has as one of its goals to provide a system that at the one time offers a compact and robust way of implementing a noise barrier but where the hereabove described disadvantages are mitigated.
  • This is achieved by providing a system according to the features of claim 1.
  • In Garcia-Bonito et al: "Generation of zones of quiet using a virtual microphone arrangement", Journal of the Acoustical Society of America, American Institute of Physics. New York, US, vol. 101, no. 6, 1 June 1997, an arrangement is described of a local active noise control system for use near the head of a travelling person. The description mentions that a secondary source "projects" zones of quit, using a so called virtual microphone arrangement. In this respect D1 is indeed similar to the subject matter of the invention. However, D1 assumes that the primary source near the real microphone has identical characteristics as the source near the virtual microphone.
  • In particular, according to the invention, a simulated error signal of a virtual microphone is a result of a combination of a first transfer function expressing the far field character of the primary source near said first microphone and a second transfer function expressing the near field character of the secondary source near said first microphone.
  • Conveniently, it has been found that by deriving a simulated error signal of a virtual microphone placed in the far-field of the loudspeaker, the noise suppression is much more successful than conventional active noise barriers and offers reductions even up to 10 dB or more compared to a sound reduction of a conventional system with adjacently placed microphones. Since the primary and secondary sources may have a very different nature they may lead to rather different requirements for a transfer function. One big difference, for example, is that the near field sensors are located rather far away from the primary sources while the near field sensors are, by definition, in the near field of the secondary sources.
  • The loudspeaker may be arranged to produce anti-sound in a direction away from said primary noise source. Further, said microphone is preferably placed adjacent in front of said loudspeaker.
  • In a still further preferred embodiment, a second microphone may be placed between said loudspeaker and said primary source. Said second microphone may be arranged to provide a second error signal that is used as a feed forward error signal in order to achieve an expected sound pressure level of the primary source near said loudspeaker.
  • Said at least one loudspeaker and said at least first microphone may be part of an array of loudspeakers and an array of microphones respectively, wherein loudspeakers and said microphones are placed at a distance less than 5 times an interspacing between two adjacent loudspeakers. Such an array may serve as an active noise barrier on the side of a road etc. Still more preferably, said at least one loudspeaker and said at least one microphone are part of an array of loudspeakers and an array of microphones respectively, wherein said loudspeakers and said microphones are placed relative to each other in a range between 10% and 100% of an interspacing between two adjacent loudspeakers. It has been found that a position of said microphone at such a general "close" distance of the loudspeaker, is still able to derive a sufficient error signal in order to be able to calculate a far field sound pressure.
  • As a practical embodiment, preferably, said microphone is formed integral with said loudspeaker in a panel to be placed on the side of a road. To this end, the invention is also related to a panel for a noise screen to be placed on the side of the road, for actively reducing noise from a primary source, comprising: a loudspeaker to be directed away from the primary source
    a microphone attached to said panel and placed on a distance away from said panel; and a system according to any of the preceding aspects.
  • Furthermore, said panel may further comprise a second microphone placed opposite to said first microphone, viewed in a direction away from said panel. Furthermore, said loudspeaker is placed on top of said panel. In this way the influence of wind is more or less the same for the primary sources and the secondary sources. This can lead to a more robust system with respect to fluctuations of the wind.
  • In another aspect, the invention relates to providing a a system according to any of the preceding aspects, further comprising an input for inputting data representing environmental conditions such as wind speed and/or temperature, said processor arranged to provide a transfer function that is dynamically adapted to said environmental conditions.
  • In another aspect, the invention is related to a method as described in claim 12 and 13.
  • The invention will be further clarified with reference to the figures. In the figures:
    • Figure 1 shows a basic configuration according to an embodiment of the invention;
    • Figure 2, 3 and 4 show prior art system diagrams wherein noise control is actively used;
    • Figure 5 is a schematic system diagram of the system according to the invention;
    • Figure 6 shows a sectional diagram of Figure 5 explaining the invention in more detail;
    • Figure 7 shows a configuration of an active noise control system for reduction of noise at the position of the farfield sensors by using information derived from nearfield error sensors showing the region for identification related to the primary noise sources;
    • Figure 8 shows a comparison of a primary farfield signal with an estimate obtained from the nearfield sensor signals;
    • Figure 9 shows a comparison of a secondary transfer function to the farfield with an estimate obtained from the nearfield sensor signals;
    • Figure 10 shows a reduction of the broadband sound pressure for three independent broadband primary noise sources at z= -4 m; the error signal equals the nearfield pressure;
    • Figure 11 shows a reduction as in; the error signal equals a measured farfield pressure;
    • Figure 12 shows a reduction as in Figure 10; the error signal equals a virtual farfield pressure.
  • In the description, like or corresponding elements in will be referenced with the same reference numerals.
  • A proposed configuration of the system according to the invention can be found in Figure 1 for an active noise barrier 6, for instance, for blocking primary traffic noise of passing vehicles 1. In this figure reference sensors 2 (microphones) are present near secondary sources 3 (loudspeakers), and nearfield error sensors 4 (microphones). In an initial calibration phase also farfield error sensors 5 (microphones) can be used- One section (n) of the active noise barrier 6 is considered. In order to assess the stability and performance of the active noise barrier 6, including the interactions between different subsystems(n-1, n+1), also the sensors and sources of adjacent sections should be taken into account.
  • One of the aspects of the system is that virtual error signals are derived from the nearfield error sensors. These virtual error signals should represent the degrees of freedom in the farfield of the angular sector of interest. In the configuration, we distinguish the following components: the primary disturbance d, reference signals x, actuator control signals u, nearfield error signals e_y, and farfield error signals e_z such that x ϵ y ϵ z = G dx G ux G dy G uy G dz G uz d u
    Figure imgb0001
  • The action of the control operator W is described by u = Wx
    Figure imgb0002
  • Lower case variables denote vectorial quantities that are a function of time, whereas upper case variables denote matrice s that operate on the time-dependent vectors. The assumption is that G ux is compensated in the controller, leading to a control structure based on Internal Model Control. Furthermore, aspects related to decentralization are not taken into account yet but will be discussed in the remainder.
  • For purposes of comparison, some existing systems will be described with reference to the Figure 2-Figure 4. Figure 2 shows a system using a single transfer function Hyz. Such systems are used for actively controlling noise from panels. However, the performance for the present application is insufficient, as can be demonstrated in simulations. Figure 3 shows a system where the basic strategy is to minimize the nearfield pressure, using a feedforward controller, as in Figure 2. Figure 4 shows a system in which the nearfield pressure is directly fed back to the controller, where the basic strategy is also to minimize the nearfield pressure. The latter configuration has been tested by Japanese researchers. In the present application, reduction of the nearfield pressure is very different from minimizing the farfield pressure. The simulations described below will show that the configuration of Figure 5 performs substantially better than the configurations of Figure 2- Figure 4.
  • Turning to Figure 5, in a first approach for the reduction of farfield sound pressure the sound pressure is reduced near the noise barrier. This is the approach that is used in, for example, a Japanese active noise barrier as illustrated in Figure 4. A disadvantage of this system is that a reduction of the nearfield pressure not necessarily leads to a reduction of the farfield pressure. This holds for both feedback configurations and for feedforward configurations. According to the invention, an improved technique is based on the estimation of the farfield error signals e_z (which are assumed to be unavailable when the system is operational) from the nearfield error signals e_y.
  • In a first approach a transfer function Hyz is designed which makes a real-time estimate of e_z from measured data e_y. However, it was found that the maximum performance with such a transfer function is limited by the fact that we are dealing with two sets of sources, the primary sources and the secondary sources. These two sets of sources may have a very different nature and may lead to rather different requirements for the transfer function Hyz. One big difference, for example, is that the nearfield sensors are located rather far away from the primary sources while the nearfield sensors are, by definition, in the nearfield of the secondary sources. So, although a single transfer function Hyz may have acceptable performance in some systems, such as active panels, the performance for active noise barriers can be improved further.
  • In Figure 5 and Figure 6 is illustrated, that, for further improving the noise control in the far field, a new strategy was devised in which the primary sources and the secondary sources were treated separately. This estimation proceeds as follows. The actual contributions on e_y and e_z consist of a part due to the primary signal and a part due to the secondary signal. As the transfer functions from the primary source and secondary source to both e_y and e_z are quite different, these contributions are taken into account separately. Therefore, we define Gdy as the transfer function between primary source and nearfield error sensor, Gdz as the transfer function between primary source and farfield error sensor, Guy as the transfer function between secondary source and nearfield error sensor, and Guz as the transfer function between secondary source and farfield error sensor. The farfield signal e_z can be expressed as e_z = d_z + s_z
    Figure imgb0003

    in which the primary signal d_z is given by d_z = G dz d
    Figure imgb0004

    and the secondary signal s_z is given by s_z = G uz u
    Figure imgb0005
  • The estimate d^z of the signal d_z is derived from an estimate d^y of the primary nearfield error signals as follows d z = H yz d d y
    Figure imgb0006

    where Hyz {d} is obtained from system identification. The estimate d^y is obtained by subtracting the secondary signal from the measured signal e_y: d y = e_y - s y
    Figure imgb0007

    where s^y is an estimate of the secondary signal on the nearfield error sensors, which is given by s y = G uy u
    Figure imgb0008
  • An estimate of the contribution of the secondary sources on the farfield sensors is given by s z = H yz s s y
    Figure imgb0009

    where Hyz {s} is also obtained from system identification. Finally, the estimate e^z of the farfield sensor can be obtained from the summation of the estimates of the primary sources and secondary sources: e z = d_z + s_z
    Figure imgb0010
  • The estimate e^z can be written as e z = H yz d d y + H yz s s y = H yz d e_y + H yz s - H yz d G uy u
    Figure imgb0011
  • A block diagram of the above procedure is given in Figure 6. The nearfield error signal is given by e_y = d_y + s_y
    Figure imgb0012
  • Therefore, using d_y = G dy d
    Figure imgb0013

    and s_y = G uy u
    Figure imgb0014

    the estimate e^z can be written as e z = H yz d G dy d + G uy u + H yz s - H yz d G uy u
    Figure imgb0015
  • In case of perfect modeling we have Guy = G^uy. Then the above expression reduces to e z = H yz d G dy d + H yz s G uy u
    Figure imgb0016
  • If the secondary sources are close to the diffraction sources and if they radiate in a similar way then we have Hyz{d} = Hyz{s} = Hyz, leading to e z = H yz G dy d + G uy u = H yz e_y
    Figure imgb0017
  • Stability robustness and performance robustness are governed by the sensitivity of the control system to changes in the secondary path. The propagation distances involved in the secondary path are considerably reduced if nearfield sensors are used instead of the farfield sensors. Therefore, it is expected that the reduced change of the propagation path will have a positive effect on the stability robustness of the system. Furthermore, if the primary sources are effectively diffraction sources on top of a noise barrier and if the secondary sources are close to these diffraction sources then we obtain a system in which a large part of the possibly varying transmission path is common to both the primary signals and the secondary signals. This results in an expected improved performance robustness with respect to changes in wind, temperature, etc.
  • Eq. (17) suggests that minimizing e_y would be sufficient to reduce e^z. This only holds if the primary sources exactly coincide with the secondary sources. In practical situations it is very difficult to exactly realize the condition Hyz{d} = Hyz{s} = Hyz. Differences are caused by, for example, the distributed nature of the primary (diffraction) sources while the secondary sources are concentrated sources. These different behaviors are especially seen on the nearfield sensors. Also the exact position of the primary sources may be difficult to determine. Therefore, for a practical implementation a more general approach is required, such as based on the scheme of Figure 6 where the estimation of the primary path and the estimation of the secondary path are treated separately. For the estimation of the primary transfer functions we define a region 7 with identification sources 8 as shown in Figure 7.
  • Furthermore, only those (and precisely those) nearfield pressure distributions should be controlled that contribute to the farfield, as governed by the transfer functions Hyz. If e_y is minimized directly then usually only sound pressure reductions in the nearfield are obtained. An additional reduction of the order of Hyz can lead to systems that are potentially more robust, have less spillover and require less actuator effort.
  • In order to be able to compute the virtual error signals an initial calibration phase is assumed in which sources are placed at random positions in the region where the primary noise sources are to be expected. The estimate of the primary signal on one of the farfield sensors is shown in Figure 8. The estimate of one of the secondary source transfer functions to the farfield as estimated from the secondary transfer functions to the nearfield sensors is shown in Figure 9. In an example 50 calibration sources have been used, which are indicated as green asterixes in Figure 12. Three independent primary sources are used, also shown in Figure 12. The identification procedure is based on the solution of the equations resulting from an assumed multi-input, multi-output Finite Impulse Response model leading to a block-Toeplitz structure. In order to improve the stability of the solution, the mean-square value of the coefficients of the FIR-model are weighted with a normalized coefficient weighting factor β = 10-2 was compared to the mean-square prediction error of the farfield pressure. This value of β was obtained by reducing the value of β until the prediction error for primary signal validation data did not reduce anymore. The latter prediction error was obtained with validation data from an independent set of primary sources at random positions in the same region as for the solution of the system of equations.
  • Although the transfer functions Hyz{d} and Hyz{s} relate the signals from identical nearfield sensors and identical farfield sensors, they should be obtained from separate calibration procedures, as mentioned previously. The necessity of this is demonstrated in Table 1, which gives the estimation results on farfield sensor for different combinations of calibration sources and the two transfer functions Hyz{d} and Hyz{s}. The configuration is identical to that in Fig. 12. It can be seen that the estimator Hyz{d} only gives accurate results if the calibration sources are placed in the same region as where the primary sources are to be expected. The same holds for Hyz{s}, which is only accurate if the calibration sources are taken to be the secondary sources. [Table 1]
    calibration sources Hyz{d} estimation error[dB] Hyz{s} estimation error[dB]
    sources in primary source region -14.6 10.5
    secondary sources -2.6 -36.3
  • In one embodiment, the system can be adapted to conditions such as wind speed and/or temperature. To this end the transfer function is dynamically adapted to said environmental conditions. The objective is to make modifications to the J × I dimensional matrix transfer functions H yz d t
    Figure imgb0018
    and H yz s t
    Figure imgb0019
    using measurements of the wind speed vector U such that these functions more accurately predict farfield sound pressure from measured near field sound pressure. A single effective wind speed vector is assumed for the whole region of interest. The wind speed vector U is transformed to a matrix of temporal interpolation functions τ ij (t) having the dimension time, which are a function of the argument time t, which also has the dimension time; the index i =1,..., I denotes a near field sensor and the index j =1,...,J denotes a farfield sensor. This interpolations functions have the effect of modifying the sound of speed due to wind speed.
  • The transformation can be derived from measurements at various wind speed magnitudes and directions but also from a physical model. In addition, an I × I dimensional static transfer function matrix V is defined. This matrix represents a dependance of windspeed and direction. The modified transfer function matrix yz d
    Figure imgb0020
    (t) for the primary field is obtained by the operation yz d t = H y i z j d τ ij t V ,
    Figure imgb0021

    where the square brackets are used to denote that the j -th row and i-th column of the resulting J × I dimensional matrix are the transfer functions H y i z j d τ ij t .
    Figure imgb0022
    . Similarly, the modified transfer function matrix yz s
    Figure imgb0023
    for the secondary field is obtained by the operation yz s t = H y i z j s τ ij t V .
    Figure imgb0024
  • The example of Figure 10 is an active noise barrier in which the error signal equals the pressure as measured on microphones near the secondary sources. In the present case this distance, i.e. the difference between the z-coordinate of the 5 secondary sources and z-coordinate of the 5 error sensors is 0.5 m. In addition, 5 reference sensors are used, which are positioned, as seen from the secondary sources, 1 m towards the primary sources. The primary noise sources are three independent broadband noise sources which are positioned at a distance of 4.5 m from the secondary sources. The position of the passive noise barrier could be at z=0 but for reasons of simplicity there is no such noise barrier. The sound pressure in the farfield is evaluated with 5 microphones at z = 20 m. The distance between these microphones is 1.5 m. The sampling frequency is f_s = 1 kHz, the number of controller coefficients equals 128, the distance between the secondary sources is 0.3 m, as are the distances between the reference sensors and the error sensors. All transducers and sources are at a height y = 2 m. Figure 10 shows the configuration and the resulting sound pressure in the x,z-plane. It can be seen that sound pressure reductions are mainly obtained near the error sensors. The reduction at the error sensors is 24.2 dB but the reduction in the farfield is much less, being 1.9 dB.
  • If the error sensors are moved to the farfield, viz. to the positions of the evaluation microphones, then the reductions in the farfield become much higher. The average reduction of the error signal and consequently the reduction at the evaluation microphones now becomes 11.3 dB. The results are shown in Figure 11.
  • Figure 12 shows the results in case the error signal equals virtual farfield signals, using the same microphone positions as in Figure 10. The virtual sensor signals are obtained by processing the nearfield error signals e_y with a fixed operator which is determined in a calibration phase. The sources as used in this calibration phase are positioned in the region where the primary noise sources are to be expected (see Figure 12). It can be seen that the resulting farfield sound pressure reductions are slightly less than in Figure 11, being 9.0 dB. However, the farfield sound pressure reductions are considerably higher than obtained with the nearfield error sensors as used in Figure 10. The reduction of the error signals is approximately equal in case of minimizing the virtual error signals, being 11.6 dB, and in case of minimizing the true farfield error signals (11.3 dB). The numerical results for the three cases are collected in Table 2. [Table 2]
    Error signal error signal av. reduction [dB] farfield pressure av. reduction [dB]
    Figure 10 24.2 1.9
    Figure 11 11.3 11.3
    Figure 12 11.6 9.0
  • Robustness of the controllers was evaluated for the following cases: for a different order of the controller, for a new set of primary sources in the same region as for the nominal case, for a change of the primary spectrum as compared to the nominal case, and for a moving source. It was found that a reduction of the number of controller coefficients leads to less reduction for the nominal case but also to less performance degradation for changing primary source distributions. Similar trends were found for the case of 32 controller coefficients. Hence, as compared to the controllers with 128 coefficients, the performance of the systems with 32 and 64 coefficients is less dependent on the actual primary source distribution and is therefore more stable and predictable. The influence of a simultaneous change of the positions of the primary sources and a changing primary spectrum was investigated using a modified primary spectrum that was obtained by applying a fourth-order bandpass filter with Butterworth characteristic with lower and upper cutoff frequencies of 0.15 fs and 0.35 fs, respectively, to the primary source signals. With the modified primary source positions and spectrum, a slight improvement was found of the performance on the farfield evaluation microphones from 10.6 dB to 11.4 dB. Furthermore, the performance for a moving source was investigated. The radiation characteristics of the moving sources were simulated by taking into account the doppler shift as well as the changing radiation characteristics depending on the direction of the movement. Also, simulations were performed in order to study the performance variations for a change in the height of the primary sources.

Claims (13)

  1. Active noise barrier system for attenuating sound from a primary noise source, such as traffic noise comprising:
    - an array of loudspeakers to be directs away from the primary source;
    - an array of microphones as a secondary source placed adjacent to said array of loudspeakers;
    - a loudspeaker connector for connecting to the array of loudspeakers for generating anti-sound for reducing said noisy sound;
    - a microphone connector for connecting to the array of microphones placed adjacent to said loudspeaker;
    - a control unit coupled to said microphone connector, for providing an error signal, based on the output of said microphone; and
    - a control unit for outputting a signal to said loudspeaker connector, for controlling said array of loudspeakers based on said error signal of said control unit;
    wherein said control unit is provided with a processor adapted to provide said error signal as a simulated error signal of a virtual microphone placed in the far-field of the array of loudspeakers,
    characterized in that said simulated error signal of said virtual microphone is a result of a combination of a first transfer function expressing the far field character of the primary source near said first microhone array and a second transfer function expressing the near field character of the secondary source near said array of microphones, both transfer functions being obtained from separate calibration procedures.
  2. A system according to claim 1, wherein said array of microphones is placed adjacent in front of said array of loudspeakers.
  3. A system according to any of the preceding claims, wherein a second array of microphones is placed between said array of loudspeakers and said primary source.
  4. A system according to claim 3, wherein said second array of microphones is arranged to provide a reference signal that is used as a feed forward signal in order to achieve an expected sound pressure level of the primary source near first array of microphones.
  5. A system according to claim 1, wherein said array of loudspeakers and said array of microphones are placed at a distance less than 5 times an interspacing between two adjacent loudspeakers.
  6. A system according to claim 1, wherein said array of loudspeakers and said array of microphones are placed relative to each other in a range between 10% and 100% of an interspacing between two adjacent loudspeakers.
  7. A system according to claim 1, wherein said array microphone is formed integral with said array of loudspeakers in a panel to be placed on the side of a road.
  8. A system according to any of the preceding claims, further comprising an input for inputting data representing environmental conditions such as wind speed and/or temperature, said processor arranged to provide a transfer function that is dynamically adapted to said environmental conditions.
  9. A system according to claim 1 further comprising a second array of microphones placed opposite to said first array of microphones, viewed in a direction away from said noise barrier.
  10. A system according to claim 1 wherein said loudspeaker is placed on top of said barrier.
  11. A method of calibrating a system according to any of the claims 1-10 comprising:
    - placing said system near the side of the road
    by placing a far field microphone at a distance away from said panel and
    said road
    placing an array of loudspeakers for generating anti-sound for reducing said traffic noise and
    placing an array of microphones adjacent said array of loudspeakers;
    - measuring an error signal derived from a first microphone placed adjacent to said array of loudspeakers;
    - measuring a reference signal derived from the far field microphone;
    - calibrating said system by a comparison between said error signal and said reference signal.
  12. A method for actively reducing sound from a primary noise source, such as traffic noise, comprising:
    - deriving an error signal, based on the output of a first microphone placed adjacent to an array of loudspeakers as a secondary source; and
    - outputting a driving signal to said array of loudspeakers based on said error signal for generating anti-sound for reducing said noisy sound;
    - wherein said error signal is a simulated error signal of a virtual microphone placed in the far-field of the array of loudspeakers and wherein said simulated error signal of said virtual microphone is obtained by combining a first transfer function expressing the far field character of the primary source near said first microphone array and a second transfer function expressing the near field character of the secondary source near said first microphone, both transfer functions being obtained from separate calibration procedures.
  13. A method according to claim 12, further comprising receiving data representing environmental conditions such as wind speed and/or temperature, and providing a transfer function that is dynamically adapted to said environmental conditions.
EP05729401A 2004-03-31 2005-03-31 System for actively reducing sound Not-in-force EP1733381B1 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
EP05729401A EP1733381B1 (en) 2004-03-31 2005-03-31 System for actively reducing sound

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
EP04076031A EP1583075A1 (en) 2004-03-31 2004-03-31 System for actively reducing sound
PCT/NL2005/000234 WO2005096269A1 (en) 2004-03-31 2005-03-31 System for actively reducing sound
EP05729401A EP1733381B1 (en) 2004-03-31 2005-03-31 System for actively reducing sound

Publications (2)

Publication Number Publication Date
EP1733381A1 EP1733381A1 (en) 2006-12-20
EP1733381B1 true EP1733381B1 (en) 2012-12-19

Family

ID=34878274

Family Applications (2)

Application Number Title Priority Date Filing Date
EP04076031A Withdrawn EP1583075A1 (en) 2004-03-31 2004-03-31 System for actively reducing sound
EP05729401A Not-in-force EP1733381B1 (en) 2004-03-31 2005-03-31 System for actively reducing sound

Family Applications Before (1)

Application Number Title Priority Date Filing Date
EP04076031A Withdrawn EP1583075A1 (en) 2004-03-31 2004-03-31 System for actively reducing sound

Country Status (2)

Country Link
EP (2) EP1583075A1 (en)
WO (1) WO2005096269A1 (en)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9704509B2 (en) * 2015-07-29 2017-07-11 Harman International Industries, Inc. Active noise cancellation apparatus and method for improving voice recognition performance
EP3185240A1 (en) * 2015-12-22 2017-06-28 Helmut-Schmidt-Universität System and method for actively reducing noise passing through an opening in a sound barrier

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5315661A (en) * 1992-08-12 1994-05-24 Noise Cancellation Technologies, Inc. Active high transmission loss panel
US5519637A (en) * 1993-08-20 1996-05-21 Mcdonnell Douglas Corporation Wavenumber-adaptive control of sound radiation from structures using a `virtual` microphone array method
US5715320A (en) * 1995-08-21 1998-02-03 Digisonix, Inc. Active adaptive selective control system
US5917919A (en) * 1995-12-04 1999-06-29 Rosenthal; Felix Method and apparatus for multi-channel active control of noise or vibration or of multi-channel separation of a signal from a noisy environment
GB9927131D0 (en) * 1999-11-16 2000-01-12 Royal College Of Art Apparatus for acoustically improving an environment and related method

Also Published As

Publication number Publication date
EP1583075A1 (en) 2005-10-05
WO2005096269A1 (en) 2005-10-13
EP1733381A1 (en) 2006-12-20

Similar Documents

Publication Publication Date Title
CN105405438A (en) Adaptive noise control system with improved robustness
EP1856948B1 (en) Position-independent microphone system
Haneda et al. Common acoustical pole and zero modeling of room transfer functions
US7778425B2 (en) Method for generating noise references for generalized sidelobe canceling
CN104158990A (en) Method for processing an audio signal and audio receiving circuit
JP4376902B2 (en) Voice input system
US20170251301A1 (en) Selective audio source enhancement
Nakatani et al. Blind dereverberation of single channel speech signal based on harmonic structure
KR20190126069A (en) Signal processing apparatus and method, and program
DE102005037034B4 (en) Method and system for controlling the energy density using a two-dimensional energy density sensor
Shi et al. Feedforward multichannel virtual-sensing active control of noise through an aperture: Analysis on causality and sensor-actuator constraints
Saito et al. Influence of modeling error on noise reduction performance of active noise control systems using filtered-x LMS algorithm
CN104781880A (en) Apparatus and method for providing informed multichannel speech presence probability estimation
Chen et al. An active noise barrier with unidirectional secondary sources
US11276387B2 (en) Proximity compensation for remote microphone ANC algorithm
CN111213359A9 (en) Echo canceller and method for echo canceller
JP4735319B2 (en) Active vibration noise control device
EP1733381B1 (en) System for actively reducing sound
Elliott et al. Superposition of the uncertainties in acoustic responses and the robust design of active control systems
US6959092B1 (en) Noise reduction panel arrangement and method of calibrating such a panel arrangement
Lombard et al. Multidimensional localization of multiple sound sources using blind adaptive MIMO system identification
Jung et al. The effect of remote microphone technique and head-tracking on local active sound control
De Diego et al. Performance evaluation of multichannel adaptive algorithms for local active noise control
Comminiello et al. A novel affine projection algorithm for superdirective microphone array beamforming
CN115604627A (en) Audio signal processing method and device, electronic equipment and readable storage medium

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

17P Request for examination filed

Effective date: 20061013

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 IS IT LI LT LU MC NL PL PT RO SE SI SK TR

DAX Request for extension of the european patent (deleted)
17Q First examination report despatched

Effective date: 20081216

GRAP Despatch of communication of intention to grant a patent

Free format text: ORIGINAL CODE: EPIDOSNIGR1

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 IS IT LI LT LU MC NL PL PT RO SE SI SK TR

REG Reference to a national code

Ref country code: GB

Ref legal event code: FG4D

REG Reference to a national code

Ref country code: CH

Ref legal event code: EP

REG Reference to a national code

Ref country code: AT

Ref legal event code: REF

Ref document number: 589743

Country of ref document: AT

Kind code of ref document: T

Effective date: 20130115

REG Reference to a national code

Ref country code: DE

Ref legal event code: R096

Ref document number: 602005037486

Country of ref document: DE

Effective date: 20130221

REG Reference to a national code

Ref country code: CH

Ref legal event code: NV

Representative=s name: PATENTANWAELTE SCHAAD, BALASS, MENZL AND PARTN, CH

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

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: 20121219

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: 20121219

Ref country code: LT

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: 20121219

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: 20130330

REG Reference to a national code

Ref country code: NL

Ref legal event code: T3

REG Reference to a national code

Ref country code: AT

Ref legal event code: MK05

Ref document number: 589743

Country of ref document: AT

Kind code of ref document: T

Effective date: 20121219

REG Reference to a national code

Ref country code: LT

Ref legal event code: MG4D

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

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: 20121219

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: 20130320

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

Ref country code: IS

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: 20130419

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: 20121219

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: 20121219

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: 20121219

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: 20121219

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: 20130319

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: 20130419

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: 20121219

Ref country code: PL

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: 20121219

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

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: 20130331

Ref country code: DK

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: 20121219

26N No opposition filed

Effective date: 20130920

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: 20121219

REG Reference to a national code

Ref country code: IE

Ref legal event code: MM4A

REG Reference to a national code

Ref country code: DE

Ref legal event code: R097

Ref document number: 602005037486

Country of ref document: DE

Effective date: 20130920

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

Ref country code: IE

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

Effective date: 20130331

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

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: 20121219

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; INVALID AB INITIO

Effective date: 20050331

Ref country code: LU

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

Effective date: 20130331

REG Reference to a national code

Ref country code: FR

Ref legal event code: PLFP

Year of fee payment: 12

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

Ref country code: NL

Payment date: 20160321

Year of fee payment: 12

Ref country code: CH

Payment date: 20160323

Year of fee payment: 12

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

Ref country code: GB

Payment date: 20160321

Year of fee payment: 12

Ref country code: FR

Payment date: 20160321

Year of fee payment: 12

Ref country code: BE

Payment date: 20160321

Year of fee payment: 12

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

Ref country code: DE

Payment date: 20160330

Year of fee payment: 12

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

Ref country code: IT

Payment date: 20160324

Year of fee payment: 12

REG Reference to a national code

Ref country code: DE

Ref legal event code: R119

Ref document number: 602005037486

Country of ref document: DE

REG Reference to a national code

Ref country code: CH

Ref legal event code: PL

REG Reference to a national code

Ref country code: NL

Ref legal event code: MM

Effective date: 20170401

GBPC Gb: european patent ceased through non-payment of renewal fee

Effective date: 20170331

REG Reference to a national code

Ref country code: FR

Ref legal event code: ST

Effective date: 20171130

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

Ref country code: DE

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

Effective date: 20171003

Ref country code: FR

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

Effective date: 20170331

Ref country code: NL

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

Effective date: 20170401

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 NON-PAYMENT OF DUE FEES

Effective date: 20170331

Ref country code: GB

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

Effective date: 20170331

Ref country code: LI

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

Effective date: 20170331

Ref country code: CH

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

Effective date: 20170331

REG Reference to a national code

Ref country code: BE

Ref legal event code: MM

Effective date: 20170331

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: 20170331