The disclosure relates to a system and method (generally referred to as a “system”) for controlling noise, for example, road noise.
Sound is a pressure wave which consists of alternating periods of compression and expansion. For noise-cancellation a sound wave is emitted with the same amplitude but with inverted phase (also known as antiphase) to the original sound. The waves combine to form a new wave, in a process called interference, and effectively cancel each other out—an effect which is called destructive interference. Modern active noise control (ANC) is commonly achieved through the use of analog and/or digital signal processing. Adaptive algorithms can be designed to analyze the waveform of the background aural or non-aural noise, and, based on the specific algorithm, can generate a signal that will either phase shift or invert the polarity of the original signal. This inverted signal (antiphase signal) is then amplified and a transducer creates a sound wave directly proportional to the amplitude of the original waveform, creating destructive interference. This effectively reduces the loudness of the perceivable noise.
A noise-cancellation transducer may be co-located with the sound source to be attenuated. In this case it should have the same audio power level as the source of the unwanted sound. Alternatively, the transducer emitting the cancellation signal may be located at the location where sound attenuation is wanted (e.g. a user's ear). This requires a much lower power level for cancellation but is effective only for a single user. Noise cancellation at other locations is more difficult as the three-dimensional wave fronts of the unwanted sound and the cancellation signal could match and create alternating zones of constructive and destructive interference, reducing noise in some spots while increasing noise in others. In small enclosed spaces (e.g. the passenger compartment of a vehicle) global noise reduction can be achieved via multiple speakers and error microphones, and through measurement of the modal responses of the enclosure.
Land based vehicles, when driven upon roads and other surfaces, generate low frequency noise known as road noise. As the wheels are driven over the road surface, the road noise is at least in part structure borne, i.e., it is transmitted through vehicle components such as tires, wheels, hubs, chassis components, suspension components such as suspension control arms or wishbones, dampers, anti-roll or sway bars and the vehicle body, and can be heard in the vehicle cabin. In order to reduce the vibrations in the vehicle components and hence road noise experienced by cabin occupants, ANC systems of the kind described above may be employed.
A widely used adaption algorithm with ANC systems is the Normalized Filtered X Least Mean Square (NFX-LMS) algorithm, which is used because of its known advantage of speedy convergence and therefore quick adaption to new boundary conditions. To achieve additional speed in the convergence, the goal of the algorithm may be defined so as to increase its step-size to the biggest values possible, thereby running the risk of creating an instable system. Selecting a static step-size will always be a trade-off between speed and stability. As a consequence there is a demand for new techniques allowing accelerated normalized convergence without compromising on stability. It is desirable to achieve a fast but robust ANC system e.g. for Road Noise Cancellation (RNC), without compromising performance and without taking additional risks involving instability.
An example active noise control filter arrangement with an adaptive filter structure includes a controllable filter module configured to process, according to a controllable K×M filter matrix with K≥1 and M≥1, K input signals to provide M output signals, the K×M filter matrix having variable filter coefficients and being controlled by updating the filter coefficients. The filter arrangement further includes a filter control module configured to update the filter coefficients dependent on the K input signals and L≥1 error signals, the L error signals being representative of at least one performance criterion of the filter module. The filter arrangement further includes an update control module configured to apply a leakage functionality and a convergence functionality to the updated filter coefficients. At least one of the following applies: The leakage functionality is controlled by at least one of a flush functionality, freeze functionality, spatial freeze functionality and leakage threshold, and the convergence functionality is controlled by at least one of freeze functionality and spatial freeze functionality.
An example active noise control filtering method using an adaptive filter structure includes processing, according to a controllable K×M filter matrix with K≥1 and M≥1, K input signals to provide M output signals, the K×M filter matrix having variable filter coefficients and being controlled by updating the filter coefficients. The method further includes updating the filter coefficients dependent on the K input signals and L≥1 error signals, the L error signals being representative of at least one performance criterion of the filter module. The method further includes applying a leakage functionality and a convergence functionality to the updated filter coefficients. At least one of the following applies: The leakage functionality is controlled by at least one of a flush functionality, freeze functionality, spatial freeze functionality and leakage threshold, and the convergence functionality is controlled by at least one of freeze functionality and spatial freeze functionality.
Other systems, methods, features and advantages will be, or will become, apparent to one with skill in the art upon examination of the following detailed description and appended figures. It is intended that all such additional systems, methods, features and advantages be included within this description, be within the scope of the invention, and be protected by the following claims.
The system may be better understood with reference to the following drawings and description. The components in the figures are not necessarily to scale, emphasis instead being placed upon illustrating the principles of the invention. Moreover, in the figures, like referenced numerals designate corresponding parts throughout the different views.
Referring to
The reference signals 110 are representative of the disturbing sound waves. The reference signals 111 in the frequency domain are (optionally) filtered with a filter module 112 with transfer functions that model the secondary path transfer functions to provide filtered reference signals 113 (also referred to as filtered input signals) in the frequency domain. The filtered reference signals 113 in the frequency domain and the error (microphone) signals 108 in the frequency domain (the error signals represent performance criterions of the system, e.g., the cancellation performance) are supplied to a control module 114 which generates control signals 115 in the frequency domain. The control signals 115 in the frequency domain are transformed by an inverse fast Fourier transform (IFFT) module 119 into control signals 116 in the time domain which are used to update a controllable filter module 117 (also referred to as w-filter) connected upstream of the loudspeaker 101 to supply loudspeaker signals 118 (also referred to as a filter output signals) thereto and supplied with the reference signals 110. The controllable filter module 117 provides, for example, a controllable w-filter matrix (with controllable w-filter transfer functions). Although no distinction is made in
The exemplary ANC multichannel system shown in
The primary and secondary paths may have a spectral behavior that changes over time. For example, the secondary paths may be modified whenever something is impacting or changing the acoustic chamber geometry. Thus the primary and secondary paths can also be described as Pk,l[n] with k=1 . . . K and l=1 . . . L, which is a (K×L) matrix of transfer functions representing the time dependent primary paths in the frequency domain, and Sm,l[k, n] with m=1 . . . M and l=1 . . . L, which is a (M×L) matrix of transfer functions representing the time dependent secondary paths in the frequency domain. The measured secondary paths are only “snapshots” of a given set-up so that they are treated as estimations representing a significant contribution to the adaptation process. The contribution to the adaptation process can be described by the “Summed-Cross-Spectrum”. The “Summed-Cross-Spectrum” SCSk,m[k,n] for each m and k combination may be as follows:
Taking this into account, the w-filter matrix update (coefficients wk,m[n], and updated coefficients wk,m[n+1]) can be described as below:
in which μglobal is the global defined adaptation step size (μ) and Ck,m[k,n] with k=1 . . . K and m=1 . . . M is a (K×M) matrix of time dependent convergence values (also referred to as Convergencek,m[k,n]) in the frequency domain.
The update is performed, in this example, according to a Filtered X Least Mean Square (FX-LMS) algorithm, in which X represents an input signal (e.g., one of the reference signals 111) filter update routine. However, any other appropriate algorithm may be used as well. The stability of the FX-LMS algorithm is highly dependent on the secondary path estimation accuracy and level of disturbance within the reference signals. The baseline (or background) may additionally include reference signal normalization, e.g., by way of an NFX-LMS algorithm. One normalization option is:
So the w-Filter matrix update can be rewritten as:
w
k,m[n+1]=wk,m[n]−IFFT{{NCk,m[k]}}
in which NCk,m[k,n] with k=1 . . . K and m=1 . . . M is a (K×M) matrix of normalized and time dependent convergence values in the frequency domain.
The system described below will not distinguish between different Normalized Filtered X Least Mean Square (NFX-LMS) variants. It is further assumed that the previously proposed normalization is used. The normalization applies a reciprocal, frequency dependent scaling to the summed cross spectrum by the energy of the reference signal. Hence the convergence step size automatically adjusts to the reference signal's spectral energy, leading to an adaptation rate which will be as fast as possible, independent from the spectral energy content of the reference signals.
Reference signal normalization does, by no means, eliminate the need of introducing a reference signal threshold definition to control the update process, known as Modified Filtered X Least Mean Square (MFX-LMS) algorithm. Nevertheless, the introduced system can further be enhanced by including such an algorithm. Although the normalization already improves ANC systems, additional techniques may be applied to further enhance stability and/or performance.
For example, the baseline assumes predefined, frequency dependent step size (μ) values which may be defined as:
Here the w-filter update process can be rewritten as:
in which μk,m[k] with k=1 . . . K and m=1 . . . M is a (K×
It is noted that independent of the applied convergence method further such methods can be applied without any restrictions. Therefore, the w-filter update process can be rewritten as:
w
k,m[n+1]=wk,m[n]−FFT{∇Jk,m[k,n]}
This implies that each convergence method can be substituted by another method without affecting the proposed improvements.
As can be seen, in this example the step-sizes are shaped over all frequency bins for each w-filter matrix index ‘m’ and I′, which represent one step size tuning set. Additionally, the baseline assumes a leakage factor that is already introduced within the w-filter update process along with the above described normalized convergence step-size, as shown in
In the flow chart shown in
The leakage value (in the following also referred to as Lk,m[k]) can be regarded as the w-filter's “oblivion” factor, with which the currently adapted w-filter coefficient values will be “forgotten”, i.e. slowly driven to zero. The value may be tunable over frequency for each individual w-filter matrix element. If the leakage shall be used as an individual multiplication factor, the w-filter update may be performed in the frequency domain in order to avoid an otherwise required, complicated convolution.
Thus the w-filter matrix update can be described as follows:
in which
W
k,m[k,n]=FFT{wk,m[n]}
W
k,m[k,n+1]=FFT{wk,m[n+1]}
L
k,m[k]=Leakagek,m[k]]
However, by definition, introduction of a leakage factor reduces the system performance because leakage and the update term act against each other. Therefore, in the following, leakage is only used as an instrument for protection against instability due to changes in the secondary paths. Furthermore, basic control features which provide control over the w-filter update via leakage and the update term are introduced. The basic control features allow for enhancing the flush mechanism, freeze mechanism, spatial freeze mechanism, and leakage threshold. Those basic control features further stabilize the system without requiring additional memory and central processing module (CPU) capacity. Introduction of basic control features within the w-filter update process in the frequency domain may be performed, for example, with basic logic modules that control the update process as shown in
Referring to
in which
FLk,m[k]=FlushLeakagek,m[k].
It is assumed that the update term contribution is weak compared to the leakage factor weighted by the flush effect and therefore the w-filter coefficients start to fade out.
Referring to
The freeze control module 401 is implemented to immediately freeze the current adaption process by bypassing the leakage factor (205) and to zero the matrix of update terms 206 in the frequency domain:
Referring to
The update process may be disabled by the freeze mechanism. The spatial freeze module 501 may toggle a spatial freeze flag and change the adaption process as follows:
in which SFLk,m[k] with k=1 . . . K and m=1 . . . M is a matrix of spatial freeze leakage values in the frequency domain, and SF∇Jk,m[k, n] with k=1 . . . K and m=1 . . . M is a matrix of spatial freeze update terms that are time dependent in the frequency domain.
Protection is achieved by the spatial freeze module 501 as it temporary limits the bandwidth of both the leakage and the update term in the frequency domain. Once the spatial freeze applies, only the upper frequency bins of the update term and the leakage are frozen, while the lower frequency bins stay as tuned:
in which SFBin is the spatial freeze limit/boundary bin. This method need not be limited to a sharp transmission between non spatial frozen and spatial frozen values, also variations of smooth transmissions techniques may be applied.
Vehicle information such as vehicle chassis information 302 and/or reference signal evaluation information such as reference signal evaluation information 402 is used to provide feedback to the flush control module 301, freeze control module 401 and/or spatial freeze control module 501. The vehicle information and/or reference signal evaluation information may execute common debounce algorithms, e.g., including hysteresis techniques, in order to avoid unwanted on/off feedback behavior to consecutive modules.
The flush control module 301 provides a flush detection that may be triggered, for example, by the reference signal and/or a vehicle information in case an already adapted w-filter has an invalid w-filter matrix and may cause hearable artifacts, because the primary path is expected to be permanently changing or one or more system components (e.g. sensors or loudspeakers) are detected as permanently being offline. Here the regular adaption process of applying the update term 206 and leakage factor 205 is insufficient or slow. Therefore, in order to ensure a safe re-adaption of the w-filter within a given new situation to an optimal w-filter setup, the w-filters become partly or completely flushed within a defined fading time.
The flush mechanism may be suitable in special scenarios in which a permanent significant and rapid change of the road noise and/or primary path is expected such as, for example, when using retractable tire studs, changing tires (summer to winter and vice versa), modifying suspension or acoustically relevant chassis components, applying dynamic driving modes as (e.g., sport and comfort mode), and in off-road suspension stiffness setups, and car-lift setups.
Also, if one or more peripheral sub-systems permanently fail, the remaining system may continue successfully with normal operation after a complete flush and re-adaptation. The vehicle on-board or on-system diagnostic may detect such permanent failures. According to a decision matrix it may be evaluated whether an operation on the remaining system can successfully continue. The term successfully is understood herein to mean that a sufficient attenuation is expected based on real measurements or simulations of such scenarios. For example, sub-systems such as error microphones, accelerometers and loudspeakers may fail.
Freeze trigger evaluation may be used to trigger the freeze module in order to preventthe already adapted system from becoming instable and/or losing performance as the w-filter coefficients could adapt to an un-desired target during ramp-up. It is assumed that the freeze control module 401 will be active only temporally, for example, in case of non-road related disturbances, high reference signal impacts, and/or low reference signal levels.
Regarding the non-road related disturbances, the impact of wind noise, for example, increases with increasing vehicle speed and at a certain level the wind noise drowns out the internal cabin noise. In such a scenario, further w-filter adaptation may be disabled by defining a maximum vehicle speed threshold to trigger the freeze mechanism. Non-road related disturbances may include at least one of wind noise, fan noise (e.g. air conditioning or other compressor modules using ventilators), audio signals from infotainment and/or entertainment systems, passenger speech and other vehicle interior disturbances.
Regarding the high reference signal impacts, adequate evaluation of the reference signal (e.g., reference signal and/or a vehicle information evaluation 303) may detect roads with too many excessively high impacts. In order to protect the adaption process to a high number of such unusual broadband impacts and an absence of stationary ones, the freeze mechanism may be triggered. Vehicle off-road information may also be used to enhance the detection process.
Regarding the low reference signal levels, another suitable scenario for freezing the adaptation entails defining a lower threshold limit for the reference signal level, so that the freeze control module 401 is triggered if the reference signal level is below a minimum value. For example, one of the two ways described below may be advantageous over simply detecting an excessively low reference signal level. One is to permanently evaluate the reference signal and to trigger the freeze control module 401 once the signal is below a certain threshold level. The other is to define a vehicle speed range, e.g. 0-15 [km/h] in which the reference signal level is known to be below a certain threshold level.
In order to evaluate the spatial freeze trigger, the spatial freeze control module 501 is employed which improves the robustness and stability of the system, e.g., in situations in which the secondary path is expected to change such as when a door or window, or the roof, sunroof or trunk is opened or closed, seats are modified, shifted or folded, and sunblinds are used. As some changes may not lead to a complete invalidation of the secondary paths and, respectively, of the estimations, in such cases the adaption process may partly continue with restrictions. The lower spectral components of an estimated secondary path may be still valid and may be used by the adaptation process. Here the spatial freeze bandwidth limit may be individually set to the last valid secondary path spectral component.
Referring to
In leakage threshold module 601, a threshold may be defined for enabling leakage so that the w-filters could first deploy to a certain level at the beginning of an adaption or in case they have been flushed. The leakage threshold module 601 distinguishes between already adapted systems and systems in the ramp-up phase of the adaptation. It is assumed that during ramp-up, the leakage factors should be less pronounced compared to the leakage applied once the system is fully deployed:
in which NBins is the number of frequency bins, RLk,m[k] with k=1 . . . K and m=1 . . . M is a matrix of ramping leakage values, and LTHk,m with k=1 . . . K and m=1 . . . M is a matrix of leakage threshold values in the frequency domain, in which:
RL
k,m[k]=RampingLeakagek,m[k] and
1.0≥RLk,m[k]≥Lk,m[k].
Leakage freeze may be applied once the ramping leakage values equal one, which may be a valid setup for fast adaptation. For example, the ramping leakage values needs to be greater than the tuned leakage values to allow an accelerated deployment of the w-filter coefficients. Instead of a single threshold value, several threshold values (LTHi,k,m) may be used to gradually change the applied leakage value, but the used leakage values may always comply with the following inequality, in which NThreshold is the number of threshold boundaries, LTHi,k,m with i=1 . . . NThreshold, k=1 . . . K and m=1 . . . M is a matrix of leakage threshold values in the frequency domain, and RLi,k,m[k] with i=1 . . . NThreshold, k=1 . . . K and m=1 . . . M is a matrix of leakage values, in which:
RL
i,k,m[k]=RampingLeakagei,k,m[k] and
1.0≥RL1,k,m[k]≥RL2,k,m[k]≥. . . ≥RLn,k,m[k].
so that
Referring to
The flush functionality may detect the validity of the updated filter coefficients and set to a given value or ramp down the updated filter coefficients within a defined time period if the updated filter coefficients are detected to be invalid. The freeze functionality may withhold the updated filter coefficients so that the updating of the filter coefficients is disabled. The spatial freeze functionality may lower spectral parts of the filter coefficients with either a hard spectral limit or a smooth spectral transition. The leakage threshold may detect whether the active noise control filter is in an adapting or re-adapting state (e.g., after a flush process) or adapted state and adjusts the leakage functionality dependent on the detected state. The filter control module and the update control module may be operated in the frequency domain, wherein, in the frequency domain, the leakage functionality may be applied to the updated filter coefficients by multiplying a leakage factor with the updated filter coefficients and the convergence functionality may be applied to the updated filter coefficients by subtracting a convergence value from the updated filter coefficients. The at least one of flush functionality, freeze functionality and spatial freeze functionality may be controlled dependent on at least one of ambient information or the input signal. Ambient information may be, for example, information provided by a vehicle on its conditions and ambient conditions in case the method is applied in a road noise control system, an engine order control system, or any other noise control system in the vehicle.
Referring to
The spatial freeze control module 501 receives, for example, vehicle/chassis information 302. The spatial freeze control module 501 provides the spatial freeze request signal to the detection module 502. If no spatial freeze request is detected by the detection module 502, the leakage factor 205 is kept unchanged (Lk,m[k] or RLk,m[k]) and the update term 206 is set to ∇Jk,m[k,n]. If, however, a spatial freeze request is detected by the detection module 502, the leakage factor 205 is set to matrix SFLk,m[k], k=1 . . . K and m=1 . . . M, which represents the spatial freeze leakage values in the frequency domain, and the update term 206 is set to matrix SF∇Jk,m[k,n], k=1 . . . K and m=1 . . . M, which represents a frequency dependent spatial freeze update term in the frequency domain.
The freeze control module 401 receives, for example, vehicle/chassis information 302, reference signal evaluation information 402 and error signal evaluation information 403. The freeze control module 401 provides the freeze request signal to detection module 404. If no freeze request is detected by the detection module 404, the leakage factor 205 and the update term 206 are kept unchanged. If, however, a freeze request is detected by the detection module 404, the leakage factor 205 is set to 1 and the update term 206 is set to 0. The update term 206 is used for the subtraction 207.
The flush control module 301 receives, for example, vehicle/chassis information 302 and diagnostic information 303. The flush control module 301 provides a flush request signal to detection module 304. If no flush request is detected by the detection module 304, the leakage factor 205 is kept unchanged. If, however, a flush request is detected by the detection module 304, the leakage factor 205 is multiplied (e.g., by multiplier 306) with a flush leakage matrix 305, i.e., FLk,m[k], and the product of the two is used in the multiplication operation 204.
As can be seen, the leakage factor 205 and the update terms 206 shown in
The description of embodiments has been presented for purposes of illustration and description. Suitable modifications and variations to the embodiments may be performed in light of the above description or may be acquired from practicing the methods. For example, unless otherwise noted, one or more of the described methods may be performed by a suitable device and/or combination of devices. The described methods and associated actions may also be performed in various orders in addition to the order described in this application, in parallel, and/or simultaneously. The described systems are exemplary in nature, and may include additional elements and/or omit elements.
As used in this application, an element or step recited in the singular and proceeded with the word “a” or “an” should be understood as not excluding plural of said elements or steps, unless such exclusion is stated. Furthermore, references to “one embodiment” or “one example” of the present disclosure are not intended to be interpreted as excluding the existence of additional embodiments that also incorporate the recited features. The terms “first,” “second,” and “third,” etc. are used merely as labels, and are not intended to impose numerical requirements or a particular positional order on their objects.
While various embodiments of the invention have been described, it will be apparent to those of ordinary skilled in the art that many more embodiments and implementations are possible within the scope of the invention. In particular, the skilled person will recognize the interchangeability of various features from different embodiments. Although these techniques and systems have been disclosed in the context of certain embodiments and examples, it will be understood that these techniques and systems may be extended beyond the specifically disclosed embodiments to other embodiments and/or uses and obvious modifications thereof.
Filing Document | Filing Date | Country | Kind |
---|---|---|---|
PCT/IB2016/056305 | 10/20/2016 | WO | 00 |