The disclosure relates to a system and method (generally referred to as a “system”) for adaptive beamforming.
Far field sound capturing allows, for example, recording sounds from a desired sound source positioned at a greater distance (e.g., several meters) from a far field microphone. The greater the distance between sound source and the far field microphone, however, the lower the desired sound to noise ratio will be. The term “noise” in the instant case includes sound that carries no information, ideas or emotions, e.g., no speech or music. Noise is commonly undesired and may be also referred to as interfering noise. When speech or music is introduced into a noise-filled environment such as a home or office interior, the noise present in the interior can have an undesired interfering effect on a desired speech communication or music presentation. Noise reduction is commonly the attenuation of undesired signals but may also include the amplification of desired signals. Desired signals may be speech or music signals, whereas undesired signals can be any sounds in the environment which interfere with the desired signals. There have been three main approaches used in connection with noise reduction: Directional beamforming, spectral subtraction, and pitch-based speech enhancement. Systems designed to receive spatially propagating signals often encounter the presence of interference signals. If the desired signal and interferers occupy the same temporal frequency band, then temporal filtering cannot be used to separate the desired signal from interference. There is a desire to improve noise reduction.
An adaptive beamforming system is configured to process at least two input signals to provide an output signal, wherein a first input signal of the at least two input signals includes a desired signal as a main component and a second input signal of the at least two input signals includes an undesired signal as main component. The system includes an error extraction block configured to adaptively process the second input signal and at least one of the first input signal and the output signal to provide an estimated undesired signal representative of an estimate of undesired signal components included in the first input signal. The system further includes a subtractor configured to take the difference between the estimated undesired signal and the first input signal to provide the output signal.
An adaptive beamforming method is configured to process at least two input signals to provide an output signal, wherein a first input signal of the at least two input signals includes a desired signal as a main component and a second input signal of the at least two input signals includes an undesired signal as main component. The method includes adaptive error processing the second input signal and at least one of the first input signal and the output signal to provide an estimated undesired signal representative of an estimate of undesired signal components included in the first input signal. The method further includes taking the difference between the estimated undesired signal and the first input signal to provide the output signal.
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. In the Figures, like referenced numerals designate corresponding parts throughout the different views.
The Figures describe concepts in the context of one or more structural components. The various components shown in the figures can be implemented in any manner including, for example, software or firmware program code executed on appropriate hardware, software and any combination thereof. In some examples, the various components may reflect the use of corresponding components in an actual implementation. Certain components may be broken apart into plural sub-components and certain components can be implemented in an order that differs from that which is illustrated herein, including a parallel manner.
It has been found that desired signals and interfering signals often originate from different spatial locations. Therefore, beamforming techniques may be used to improve signal-to-noise ratios in audio applications. Common beamforming techniques include delay-sum techniques, adaptive finite impulse response (FIR) filtering techniques using algorithms such as the Griffiths-Jim algorithm, and techniques based on models of human binaural hearing.
Beamformers can be classified as either data independent or statistically optimum, depending on how the weights are chosen. The weights in a data independent beamformer do not depend on the array data and are chosen to present a specified response for all signal/interference scenarios. Statistically optimum beamformers select the weights to optimize the beamformer response based on statistics of the data. The data statistics are often unknown and may change with time so adaptive algorithms may be used to obtain weights that converge to the statistically optimum solution. Due to computational considerations, partially adaptive beamformers with arrays composed of large numbers of sensors may be used. Many different approaches have been proposed for implementing optimum beamformers. For example, the statistically optimum beamformer places nulls in the directions of interfering sources in an attempt to maximize the signal to noise ratio at the beamformer output.
In many applications the desired signal may be of unknown strength and may always be present. In such applications, a correct estimation of signal and noise covariance matrices in a maximum signal-to-noise ratio (SNR) is not possible. Lack of knowledge about the desired signal may prevent utilization of a reference signal approach. These limitations can be overcome through the application of linear constraints to the weight vector. Use of linear constraints is an approach that allows for extensive control over the adapted response of the beamformer. A universal linear constraint design approach does not exist and in many applications a combination of different types of constraint techniques may be effective. However, attempting to find either a single best way or a combination of different ways to design the linear constraint limits the use of techniques that rely on linear constraint design for beamforming applications.
Generalized sidelobe canceller (GSC) technology presents an alternative to a linear constraint design technique for beamforming applications. Essentially, GSC is a mechanism for changing a constrained minimization problem into unconstrained form. GSC leaves the desired signals from a certain direction undistorted, while, at the same time, undesired signals radiating from other directions are suppressed. However, GSC uses a two path structure; an upper path to realize a fix beamformer pointing to the direction of the desired signal, and a lower part that adaptively generates an ideally pure noise estimate, which can then be subtracted from the output signal of the fix beamformer, thus increasing its SNR by suppressing noise.
The estimation of the noise, is usually realized in a two-step approach. The first stage of the lower path is configured to remove remaining signal parts of the desired signal from the input signals of this stage. The second stage of the lower path further comprises M (M is an integer) adaptive interference cancellers (AIC) in order to generate a single-channel, estimated noise signal, which can possibly then be subtracted from the output signal of the upper path, which is a time delayed output signal of the fix beamformer. Thus, the noise contained in the output signal of the time delayed fix beamformer can be suppressed, leading to a better SNR, as the desired signal component ideally would not be affected by this processing. This holds true if and only if all desired signal parts within the noise estimation could successfully be united, which is rarely the case in practice, and which thus represents one of the major drawbacks related to current adaptive beamforming algorithms.
Acoustic echo cancellation may be attained, e.g., by subtracting an estimated echo signal from the total sound signal. To provide an estimate of the actual echo signal, algorithms may operate in the time domain that may employ adaptive digital filters processing time-discrete signals. Such adaptive digital filters operate in such a way that the network parameters defining the transmission characteristics of the filter are optimized with reference to a preset quality function. Such a quality function may be implemented, for example, by minimizing the average square errors of the output signal of the adaptive network with reference to a reference signal.
Referring to
As can be seen from
A FIR filter 202 that has a transfer function ĥ(n) of order L−1, wherein L is a length of the FIR filter, may be used to model the echo path transfer function vector h(n). The transfer function vector ĥ(n)=h(n) can be described as
[h(0,n), . . . h(L−1,n),]T
The desired microphone signal d(n) at block 203 for the adaptive filter is given as
d(n)=sT(n)h(n)+v(n),
wherein s(n)=[s(n)s(n−1) . . . s(n−L+1)]T is a real-valued vector containing L (L is an integer) most recent time samples of the input signal, s(n), and v(n), i.e., the near-end signal.
Using the previous notations, the feedback/echo error signal x(n) is given as
x(n)=d(n)−sT(n−1)ĥ(n)=sT(n)[h(n)−ĥ(n)]+v(n),
wherein the vector ĥ(n) contains adaptive filter coefficients at time n. Vector ĥ(n) is estimated using e.g. a Least Mean Square (LMS) algorithm or any state-of the art recursive algorithm. The LMS update process using a step size of μ(n) of the LMS-type algorithm can be expressed as
ĥ(n)=ĥ(n−1)+μ(n)s(n)x(n).
Referring again to
wherein M is again the number of microphones. The FS beamformer block 104 may include a summer 301 which receives signals xi(n) from AEC block 103 via filter paths 302 having transfer functions w(n) and a coefficient element 303 which divides the output signal of summer 302 by M to generate an output signal of the FS beamformer block 104. Further output signals of FS beamformer block 104 are derived accordingly with different transfer functions w(n). A DS beamformer can be implemented by choosing wi(n) according to wi(n)=e−j2πf·τ
In
Referring to
The smoothing block 402 outputs a smoothed signal that may still contain a considerable level of noise and thus, may cause noticeable sharp discontinuities as described above. The level of voice signals typically differs distinctly from the variation of the level of background noise, particularly due to the fact that the dynamic range of the level changes of voice signals is significantly greater and the level changes occur in much shorter intervals than the level changes of background noise. Linear smoothing in a noise estimation block 403 following the smoothing block 402 would therefore smear out the sharp variation in the desired music or voice signal as well as filter out the noise. Such smearing of a music or voice signal is not desired in many applications, therefore, alternatively or additionally, non-linear smoothing (not shown) in noise estimation block 403 may be applied to the smoothed signal to overcome the artifacts mentioned above. The noise estimation block 403 outputs signals representative of estimates of the noise contained in the corresponding output signals of smoothing block 402.
In an SNR calculation block 404, based on the smoothed signals from smoothing block 402 and the estimated background noise signals from noise estimation block 403, variations of the respective SNR are evaluated and output as, e.g., a distribution of each of SNR values SNR1(n) . . . SNRB(n) over time n for a corresponding beam 1 . . . B. The evaluation and, thus, the distribution may extend over a certain time period, i.e., a frame or time window, such as from n to n−100 or the like. By evaluating variations in the SNR, a noise source can be differentiated from a desired speech or music signal. For example, low SNR values may represent a variety of noise sources such as an air-conditioner, fan, an open window, or an electrical device such as a computer etc. The SNR values SNR1(n) . . . SNRB(n) may be determined in a time domain (as shown in
In a subsequent comparator block 405, the SNR values SNR1(n) . . . SNRB(n) output by block 404 are compared to pre-determined thresholds. The thresholds may be fix or controllable, frequency dependent or independent etc. If the current SNR value exceeds a corresponding pre-determined threshold, a flag indicating a desired signal such as speech will be set to a first logic value, e.g., ‘1’. Alternatively, if the current SNR value is less than a corresponding pre-determined threshold, a flag indicating an undesired signal such as noise from an air-conditioner, fan, an open window, or an electrical device such as a computer will be set to a second logic value, e.g., ‘0’. In the present example, all thresholds have the same value SNRTH.
SNR values SNR1(n) . . . SNRB(n) from block 404 and the flags from comparator block 405 are passed on to a controller block 406 via paths #1 . . . # B. The controller block 406 counts the number of logic values “1”, i.e. when SNR signals exceed the threshold SNRTH, over a certain time period (frame, window) which may be adjustable by a parameter “TimeFrame” for all B beams. Thereby, a histogram is created for those SNR values per beam over the certain time period defined by the tunable parameter “TimeFrame”. For each frame, the sum or count over all logic values “1” for each of the B beams, which represent distinct (horizontal) directions of the sound, is determined. The direction of the B beams that has the highest (maximum or “Max”) count of logic values “1” in this frame will be considered as a positive beam, pointing to the desired signal source, e.g., a speaker. The corresponding indices identifying a time variant direction of the positive beam (its steering vector or look direction) to a desired source will be output, referred to as signal S(n). The signal S(n) is fed into a fading block 407 which prevents acoustical artifacts, such as clicks that occur when switching between beams, by fading (soft switching) between beams.
The concept described above (counting “0” values) may have some drawbacks when applied to the determination of indices of a negative beam, which ideally should point in the direction of a potential noise source, in the event of low SNR values of all B beam signals. In such a case, each of the B beam signals that have the lowest SNR may also be indicated with a logic value “1” (by inverting the logic value “0” used so far), and a histogram, defined by the same, adjustable “TimeFrame” parameter, may again be compiled over time. This means that the controller block 406 counts the number of logic values “1”, i.e. the number of indices (look directions) that have, over a certain time period (frame, window), the lowest SNR values which may be adjustable by a parameter “TimeFrame” for all B beams. Then again, the maximum of this histogram would also indicate the indices (steering vector, look direction) of the negative beam that point to a noise source.
If no distinct solution can be found or none of the B beams meet the given requirements, then the previous indices (look direction) for the positive and/or negative beam are applied. Similarly to signal S(n), the indices of the negative beam are also supplied to the fading block 407 as time varying signal I(n). Again, fading block 407 ensures that the final output signal of the identified negative beam does not suffer from acoustical artifacts, similarly to the positive beam signal. For example, a histogram of the maximum and minimum values may be compiled for a pre-determined time duration. The minimum and maximum values in a histogram are representative of at least two different output signals: at least one signal is directed towards a desired source denoted by signal S(n) and at least one signal is directed towards an interference source denoted by signal I(n).
If indices of low and high SNR values in controller block 406 change over time, a fading process is initiated in a fading block 407 that allows a smooth transition from one to the other output signal, without generating acoustic artifacts. Referring again to
The outputs of BS block 105 provide a signal with a high SNR, i.e., a positive beam signal b(n), which is used as a reference signal by the adaptive blocking filter (ABF) block 106 and one with a low SNR, i.e., a negative beam bn(n), used as a desired noise signal (further input signal) by the ABF block 106. The ABF filter block 106 may use filters controlled according to an LMS algorithm to adaptively subtract the signal of interest, represented by the reference signal b(n) (corresponding to the positive beam) from the desired noise signal bn(n) (corresponding to the negative beam) and provides an error signal e(n). The error signal e(n) obtained from the ABF block 106, which ideally is a pure noise signal also included in the positive beam signal b(n), is passed to the adaptive interference canceller (AIC) block 107 which adaptively removes the signal components (noise) contained in the positive beam signal b(n) that are correlated to the error signals, i.e., to the positive beam signal b(n).
AIC block 107 computes an interference signal employing adaptive filtering. The output of this adaptive filter is subtracted (e.g., by way of a subtractor 109) from the optionally delayed (e.g., by way of delay 108) reference signal b(n) to eliminate the remaining interference and noise components in the reference signal b(n). Finally, an adaptive post filter 110 may be connected to an output of subtractor 109 for the reduction of statistical noise components (not having a distinct autocorrelation). As in the block 106, the filter coefficients in the AIC block 107 may be updated using the adaptive LMS algorithm. The norm of the filter coefficients in the ABF block 106 and/or AIC block 107 may be constrained to prevent them from growing excessively large.
The adaptive filter controller 501 and the controllable filter 502 constitute an exemplary adaptive blocking filter which may employ, as in the present example, a constraint C so that a magnitude frequency response, e.g., the magnitude of the transfer function a(n), of the adaptive blocking filter does not exceed a predetermined value at any frequency and at any time. For example, the constraint C may limit, in the time domain, filter coefficients of the adaptive blocking filter to a predetermined value such as +/−1. Alternatively, a constraint may be implemented in the frequency domain according to
[U,V]=find(abs(A)>Max ALim),
A(U)=A(U)/abs(A(U)),
Max ALim=1,
wherein U is a parameter representing frequencies at which the constraint is to be applied, V is a parameter representing amplitudes at the frequencies U, A is the transfer function of the adaptive blocking filter in the frequency domain, and MaxALim is the predetermined maximum value of the transfer function.
Referring to
Frequency range and resolution on a frequency axis of a spectrum graph depend on a sampling rate fSAMPLE and the size of the data record (the number NRECORD of acquisition points). The number of frequency points or lines or bands in the (power) spectrum is NRECORD/2, where NRECORD is the number of signal points captured in the time domain. The first frequency line in the power spectrum represents a frequency of zero (DC). The last frequency line can be found at fSAMPLE/2. Frequency lines are spaced at even intervals of fSAMPLE/NRECORD. They are commonly referred to as frequency bins or FFT bins.
In the example shown in
The (optional) spectral smoothing of each of the beam signals B1(ω) . . . BB(ω) results in a number of frequency bins that may deteriorate the total spectral resolution. In order to reduce this effect, each spectrally smoothed signal, i.e., each output signal of the spectral smoothing block 602, may further be smoothed by way of a temporal smoothing block 603 for every bin of each of the beam signals B1(ω) . . . BB(ω). The temporal smoothing block 603 is employed to smear out impulsive distortions such as speech in the spectrally smoothed signals, and may combine frequency bins over time for each output signal of the spectral smoothing block 602 to reduce the temporal dynamics in these signals.
The temporal smoothing block 603 outputs for each of the output signals of spectrally smoothing block 602 one or more additionally temporally smoothed signals, (herein referred to as smoothed signals) that may still contain impulsive distortions as well as background noise. A noise estimation block 604, which is coupled to outputs of the temporal smoothing block 603, is employed to smear out remaining impulsive distortions and to estimate the current background noise in each output signal of the temporal smoothing block 603 (smoothed signal) derived from each output signal of the spectral smoothing block 602 (spectrally smoothed signal). In order to reduce or avoid smearing of a desired signal such as music or voice signals, non-linear smoothing (not shown) may be employed in noise estimation block 604.
Based on the smoothed signals from temporal smoothing block 603 and the estimated quasi stationary background noise signal from noise estimation block 604, variations in SNR can be calculated (e.g., as a distribution of SNR values over frequency). In an SNR calculation block 605, for each output signal of noise estimation block 604 in connection with the corresponding output signal from temporal smoothing block 603 an SNR value is computed. Due to variations in the SNR, a noise source can be differentiated from a desired speech or music signal. For example, a low SNR value may represent a variety of noise sources such as an air-conditioner, fan, an open window, or an electrical device such as a computer etc. The SNR may be evaluated in the time domain or in the frequency domain or in the sub-band domain.
In a comparator block 606, each spectral line of each SNR value SNR1(n) . . . SNRB(n) provided by SNR calculation block 605 is compared to a mutual (or individual) pre-determined threshold SNRTH. If the respective current SNR value SNR1(n) . . . SNRB(n) is greater than the pre-determined threshold SNRTH, a flag SnrFlag, which may indicate a desired speech signal, will be set to a first logic value of, e.g., ‘1’. If the respective current SNR value SNR1(n) . . . SNRB(n) is less than the pre-determined (alternatively, controllable and/or spectrally dependent) threshold SNRTH, the flag SnrFlag will be set per spectral line to a second logic value of, e.g., ‘0’ to indicate an undesired signal such as noise from an air-conditioner, fan, an open window, or an electrical device such as a computer.
The flags SnrFlag from block 606 are passed on to a min max block 607. The min max block 607 determines, for each of beam signals B1(ω) . . . BB(ω), from all flags SnrFlag related to the respective one of beach signals B1(ω) . . . BB(ω) and the corresponding SNR values the corresponding minimum SNR values Min1(k) . . . MinB(k) and maximum SNR values Max1(k) . . . MaxB(k), in which k represents a discrete time. The min max block 607 further counts beam indices 1 . . . B corresponding to the minimum SNR values Min1(k) . . . MinB(k) and beam indices 1 . . . B corresponding to the maximum SNR values Max1(k) . . . MaxB(k) and, if counting is not possible, uses a previous distribution of beam indices. For this, a histogram of the maximum and minimum SNR values may be compiled for a pre-determined time period. The minimum and maximum SNR values in a histogram are representative of at least two different output signals. At least one signal is directed towards a desired source denoted by a time variant steering vector S(n) and at least one signal is directed towards an interference source denoted by a time variant steering vector I(n).
A fader block 608 receives the steering vectors S(n) and I(n) and the signals b1(n) . . . bn(n) and provides the reference signal, i.e., the positive beam signal b(n), and the desired signal, i.e., the negative beam signal bn(n). If the indices of low and high SNR values in the min max block 607 change over time, a fading process is initiated in a fader block 608 that allows a smooth transition between its output signal(s), without generating acoustic artifacts. The output signals of the min max block 607, the reference signal b(n) and the desired signal bn(n), are the positive beam signal b(n) and the negative beam signal bn(n) over discrete time n.
Referring back to
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.
The embodiments of the present disclosure generally provide for a plurality of circuits, electrical devices, and/or at least one controller. All references to the circuits, the at least one controller, and other electrical devices and the functionality provided by each, are not intended to be limited to encompassing only what is illustrated and described herein. While particular labels may be assigned to the various circuit(s), controller(s) and other electrical devices disclosed, such labels are not intended to limit the scope of operation for the various circuit(s), controller(s) and other electrical devices. Such circuit(s), controller(s) and other electrical devices may be combined with each other and/or separated in any manner based on the particular type of electrical implementation that is desired.
It is recognized that any computer, processor and controller as disclosed herein may include any number of microprocessors, integrated circuits, memory devices (e.g., FLASH, random access memory (RAM), read only memory (ROM), electrically programmable read only memory (EPROM), electrically erasable programmable read only memory (EEPROM), or other suitable variants thereof) and software which co-act with one another to perform operation(s) disclosed herein. In addition, any controller as disclosed utilizes any one or more microprocessors to execute a computer-program that is embodied in a non-transitory computer readable medium that is programmed to perform any number of the functions as disclosed. Further, any controller as provided herein includes a housing and the various number of microprocessors, integrated circuits, and memory devices (e.g., FLASH, random access memory (RAM), read only memory (ROM), electrically programmable read only memory (EPROM), electrically erasable programmable read only memory (EEPROM)) positioned within the housing. The computer(s), processor(s) and controller(s) as disclosed also include hardware based inputs and outputs for receiving and transmitting data, respectively from and to other hardware based devices as discussed herein.
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.
Number | Date | Country | Kind |
---|---|---|---|
16199773 | Nov 2016 | EP | regional |
Filing Document | Filing Date | Country | Kind |
---|---|---|---|
PCT/EP2017/079576 | 11/17/2017 | WO | 00 |
Publishing Document | Publishing Date | Country | Kind |
---|---|---|---|
WO2018/091648 | 5/24/2018 | WO | A |
Number | Name | Date | Kind |
---|---|---|---|
20050141731 | Hamalainen | Jun 2005 | A1 |
20110274291 | Tashev | Nov 2011 | A1 |
20120027218 | Every | Feb 2012 | A1 |
20120224715 | Kikkeri | Sep 2012 | A1 |
20120239385 | Hersbach et al. | Sep 2012 | A1 |
20130287069 | Su et al. | Oct 2013 | A1 |
20140037100 | Giesbrecht | Feb 2014 | A1 |
20150317999 | Gao | Nov 2015 | A1 |
20170076720 | Gopalan et al. | Mar 2017 | A1 |
20180033447 | Ramprashad et al. | Feb 2018 | A1 |
20190104371 | Ballande et al. | Apr 2019 | A1 |
Number | Date | Country |
---|---|---|
2014154325 | Oct 2014 | WO |
Entry |
---|
Sissi Xiaoxiao Wu et al: “Rank-two transmit beamformed Alamouti space-time coding for physical-layer multicasting”, 2012 IEEE International Conference on Acoustics, Speech and Signal Processing (ICASSP 2012) : Kyoto, Japan, Mar. 25-30, 2012 ; [Proceedings), IEEE, Piscataway, NJ, Mar. 25, 2012 (Mar. 25, 2012), pp. 2793-2796, XP032227734, DOI: 10.1109/ICASSP.2012.6288497 ISBN: 978-1-4673-0045-2 the whole document. |
Abdel-Sainad A et al: Robust transmit eigen-beamforming based on imperfect channel state information, Smart Antennas, 2004. ITG Workshop on Munich, Germany Mar. 18-19, 2004, Piscataway, NJ, USA,IEEE Mar. 18, 2004 (Mar. 18, 2004), pp. 200-207, XP010780100, 001: 10. 1109/WSA.2004. 1407669 ISBN: 978-1-7803-8327-2 p. 201-p. 206. |
Number | Date | Country | |
---|---|---|---|
20190320261 A1 | Oct 2019 | US |