The present disclosure relates to sub-band processing and, in particular, to systems that reduce computational complexity and memory requirements.
Frequency domain based adaptive filtering is computationally intensive as it translates a time domain signal into multiple frequency components that are processed individually. Wider bandwidth networks provide higher throughput and better performance than narrowband networks but at the expense of an increase in processing loads and memory requirements—as the bandwidth increases, so too does the number of frequency components used to represent the signal. This increase in the number of components at higher sampling rates results in a proportional increase in processing loads and memory required.
To reduce processing load and memory requirements for ease of implementation, a frequency domain-based system compresses data through a lossy compression scheme and performs subsequent processing of the lossy signal before reconstructing back the data into the time domain with an inverse process to the one used to translate the time domain signal into multiple frequency components. It is desirable to perform the lossy compression in a way that ensures minimal perceptual distortion is introduced in the resulting reconstructed data.
Reference will now be made, by way of example, to the accompanying drawings which show example embodiments of the present application and in which:
Like reference numerals are used in the drawings to denote like elements and features.
In one aspect, the present disclosure describes a system comprising: a first logic stored in a computer-readable medium and executable by a processor that partitions and stores a frequency spectrum of bins of real and imaginary data into a smaller number of sub-bands; a second logic stored in the computer-readable medium and executable by the processor that executes a first lossy compression for a first set of the sub-bands, the first lossy compression compressing a designated magnitude of one bin in each sub-band of the first set that is representative of that sub-band and a designated phase of one bin in each sub-band of the first set that is representative of that sub-band, wherein the first set includes those sub-bands having indices that are greater than or equal to a first index; and a third logic stored in the computer-readable medium and executable by the processor that executes, subsequent to a frequency spectrum processing of the lossy compressed data rendered by the second logic, a second lossy compression for a second set of the sub-bands, the second lossy compression compressing a designated magnitude of one bin in each sub-band of the second set that is representative of that sub-band and a designated phase of one bin in each sub-band of the second set that is representative of that sub-band, wherein the second set includes those sub-bands having indices that are less than the first index and greater than or equal to a second index.
In another aspect, the present disclosure describes an echo cancellation system for processing a plurality of frames of data. The echo cancellation system includes: a processor that executes a computer readable medium comprising: computer program code that determines a first number of sub-bands for a module of the echo cancellation system to process for a first set of the plurality of frames and a second number of sub-bands for the module to process for a second set of the plurality of frames, such that a sub-band count alternates between the first number and the second number during processing of each pair of sequential frames.
Other example embodiments of the present disclosure will be apparent to those of ordinary skill in the art from a review of the following detailed descriptions in conjunction with the drawings.
A sub-band processing system processes data such that, after it is compressed and decompressed, it is restored to its original format. The sub-band processing system may compress video, audio, text, code, and/or numeric data such that little or no data is lost after a bin or file is decompressed. Some systems preserve the original data (or a representative data set) while compressing and decompressing operating data through a lossy compression. After further processing by an ancillary device or system, the sub-band processing system reconstructs and restores the data.
The sub-band processing system analysis may occur on frequency domain characteristics. To derive frequency domain properties, the signal may be broken into intervals though a multiplier function (retained in a local or a distributed computer readable medium) or multiplier device that multiplies the signal by a “window” function or a “frame” of fixed duration. To minimize spectral distortion, smooth window functions (such as Hann, Hamming, etc. retained in the local or the distributed computer readable medium) or a window filter may be used for the short-time spectral analysis. A time-to-frequency transform device, a Discrete Fourier Transform (DFT) device, or a Fast Fourier Transform (FFT) device may transform (or decompose) the short-time based signals into a complex spectrum. The spectrum may be separated into bins of magnitude and phase data or substantially equivalent complex (e.g. real and imaginary) data. A sub-band (or band) may be represented by a single bin of magnitude and phase spectra, or a collection of consecutive or successive bins represented by a common or single magnitude and phase spectra. Table 1 shows representative characteristics of an exemplary FFT device:
At a sample rate of about 8 kHz, an FFT device may transform the time domain signal into about 256 bins. Due to the complex symmetry, the FFT device may yield about 129 useful bins (e.g., 256/2+1). Each bin may represent a frequency resolution of about 31.25 Hz (e.g., 8 kHz/256). The frequency resolution of other sample rates (e.g., 16 kHz and 32 kHz) may be maintained by changing the FFT length. For example, at 16 kHz, the FFT length may be about double the FFT length of the 8 kHz sample rate. At 32 kHz, the FFT length may be about double the FFT length of the 16 kHz sample rate.
In some systems, the magnitude and phase spectra may be obtained from one or more signal processors that execute a Discrete Fourier Transform (DFT) stored in a local or a distributed memory. The output of the DFT may be represented by X(k):
for k=0 . . . N−1, where K is the frequency index for each bin, N is the time index for each sample, and N is the length of the DFT (or FFT).
The bins (R) of the FFT (or DFT) device may be partitioned into a fewer (or smaller R≥M) number of sub-bands (M). In some applications, the sub-band processing system may reduce M to a lowest possible integer that does not affect the performance or quality of a later process. In these applications, the system may generate a selected number of sub-bands that minimize perceptual error. The applications may exploit the sensitivity of the human auditory system or other systems that do not detect or process certain frequencies or are affected by certain signal distortions.
A lossy compression may compress the data such that some data is lost when the data is compressed into the sub-bands. Some sub-band processing systems compress q bins (q is an integer greater than 1) into individual sub-bands. Other systems apply a perceptual scale (through a processor or controller, for example) where the bins are grouped into sub-bands that match the frequency selectivity of the human auditory system such that the compression divides a variable sequence of uncompressed bins into a substantially equal sequence of compressed sub-bands. Perceptual distortions may be minimized by applying lower compression ratios at lower frequencies while applying higher compression ratios at higher frequencies. Table 2 describes an exemplary compression scheme in which each sub-band represents q bins for a 16 kHz sampling rate system:
3-4.6
The selected or designated magnitudes in each sub-band may be obtained by various schemes. When a maximum magnitude system is used, and a maximum magnitude is detected, the bin containing that magnitude is indexed, and stored in memory. Other systems may select the magnitude of the first bin within the sub-band that has the lowest frequency within that sub-band. Other systems may select a combination of both methods with lower compressed sub-bands making use of maximum magnitude as a selection criterion and higher compressed sub-bands making use of first bin within the sub-band.
While various lossy compression schemes may be used, the sub-band processing system may select or designate a representative phase for each sub-band. Some sub-band processing systems select the phase of a bin within the sub-band that has the lowest frequency within that sub-band. Other systems may select bins based on the index of the maximum magnitude found, and others may select some other qualitative measure.
Depending on the phase index selected, the maximum magnitude may be rotated or shifted (i.e. adjusted) to attain the selected or designated phase as the phase of index with maximum magnitude selected from its constituent bins, and the phase of the “preserved” or selected bin may be different.
In the sub-band processing system, the magnitude, |SBX(m)|, and phase, arg(SBX(m)), for each sub-band may be:
SBX(m)=|SBX(m)|,arg(SBX(m)) (2)
where |SBX(m)|=max(|X(jm)|,|X(jm+1)|, . . . ,|X(jm+Dm−1)|) (3)
arg(SBX(M))=arg(X(jm)) or arg(X(hm)) (4)
hm=arg max(|X(jm)|,|X(jm+1)|, . . . ,|X(jm+Dm−1)|) for m=0, . . . M−1 (5)
where
In some systems, common bins may be selected from the divided spectrum to attempt to preserve the phase of the sub-bands relative to each processed frame. In these systems jm and Dm may be constant (e.g. temporally invariant) while hm may change (e.g. time variant) from one aural or sound frame (or video, sound, text, code, and/or numeric data) to the next. Such systems may try to preserve the phase of the same bin within a sub-band on a frame-by-frame basis such as, for example, always the first bin of a sub-band in each frame or a common bin of a sub-band in each frame. Other systems may not try to preserve the phase from frame to frame, like when hm is the index of the selected phase that is preserved.
By maintaining magnitude and phase spectra through the adjusted sub-band spectrum, or “sparse spectrum”, the spectrum may be further processed in the frequency domain (or other domains). Adaptive filtering techniques or devices used by an acoustic echo canceller, noise cancellation, or a beam-former, for example, are sensitive to changes in phase and may need to process a consistent phase that does not change abruptly from frame to frame. In addition, an accurate approximation of the magnitude spectrum of the bins with the compressed sub-band representation is also critical. Abrupt phase changes may be identified as an impulse response that causes an acoustic echo canceller to diverge. When divergence occurs, a sub-optimal, reduced, or no echo cancellation may occur due to the mismatch between the filter coefficients and the echo path characteristics. When a divergence is declared, an adaptive filter may require time to achieve a convergence.
Making use of such systems may still result in the phase data not being consistently preserved from frame to frame for the compressed bins. In particular, any added perceptual distortion introduced by the compression may not be sufficiently minimized and, for example, the acoustic echo canceller may diverge. To ensure that subsequent processing (such as acoustic echo cancellation) is not impacted adversely by these compression techniques within a sub-band, some systems may have the magnitude, |SBX(m)|, and phase, arg(SBX(m)), compressed in two separate stages. In the first stage, only those sub-bands with indices m greater than or equal to Munc that are not susceptible to abrupt phase changes and inadequate magnitude representation in each frame are compressed. The sub-bands with indices less than Munc but greater than or equal to Munc2 are compressed afterward, following processing by a subsequent block such as the acoustic echo canceller. This can be done by selecting:
SBX(m)=|SBX(m)|,arg(SBX(m)) (6)
|SBX(m)|=|X(jm)| and (7)
arg(SBX(m))=arg(X(jm)) (8)
for m=Munc, . . . , M1−1, where
Munc2 has been previously obtained from Equations 2-5 above as the sub-band at which, if full compression had been performed, the bins in a sub-band are represented by a single bin and where Dm>1. In addition, though bands are not compressed between Munc2 and Munc in the first stage, their indices, as given by hm in Equation 5, are stored for subsequent use in the second stage of compression during the subsequent processing block.
This subsequent processing block requires as its input, processing bins that are not compressed below a threshold Runc to ensure adequate minimization of perceptual distortion being introduced by the compression process. After the subsequent processing step, such as acoustic echo cancellation for example, sub-bands less than Munc2 are not compressed and sub-bands greater than or equal to Munc2 are compressed in two stages, with bands greater than or equal to Munc being compressed in the initial stage and the remaining sub-bands being compressed in the subsequent processing step. Table 3 describes an exemplary two-stage compression scheme in which each sub-band represents q bins for a 16 kHz sampling rate system:
3-4.6
In Table 3 above, Munc2 is 96, Munc is 152, Runc is 251, M is 154, M1 is 253, and R is 257. The first stage compression will take R=257 bins and compress them down to M1=253 sub-bands, and then the subsequent processing module will take the M1=253 sub-bands and further compress them down to M=154 sub-bands. These compressed sub-bands can then be used in further subsequent processing after the acoustic echo cancellation stage, for example.
An acoustic echo cancellation processing system consists of various sub-systems that are processed for every input frequency band on a frame by frame basis. Some of the parameters computed in these sub-modules may vary significantly from frame to frame while others may vary seldomly. Some of the parameters may be sub-sampled, so that instead of having the parameters updated every frame, they may be updated every nth frame instead, slowing down the speed of adaptation, which in turn may negatively affect the performance of the echo canceller. Higher frequency sub-bands are much less negatively impacted by this sub-sampling and so some benefit could be obtained by sub-sampling unequally the different sub-bands in the system without any appreciable degradation in perceptual distortion in the reconstructed output speech.
In some systems, the number of sub-bands used in the various sub-modules within the acoustic echo cancellation module are alternated from frame to frame, with every odd processing frame making use of M1 number of sub-bands, and so some savings in complexity reduction could be obtained depending on how much smaller M1 is when compared to R. On the even processing frames, the number of bins used would be set to Munc2 which may be much smaller than M1. This may lead to processing load savings in the sub-module in question as the number of sub-bands would have dropped from M1 to Munc2. In an exemplary system, with a 16 kHz sampling rate as an example, the sub-band count could alternate between 253 and 96, thereby yielding similar savings in processing loads to one in which all compression was achieved in a single-stage. The impact of slowing down the adaptation of various sub-modules for the compressed sub-bands on the perceptual quality may be minimal. Some sub-modules may need to have their adaptation parameters changed to adapt faster for those sub-bands which were sub-sampled. This way, the net change in long-term adaptation is minimal across those sub-bands when compared to sub-bands that are not sub-sampled.
In some other systems, the compressed band index information hm that was computed by Equation 5 in the first stage could be used to decide which of the bins in a compressed sub-band to update in a frame. Only those bins represented by index hm are adapted. This implies that within each of the sub-modules, not all bins are processed, skipping over hm−1 bins in each sub-band. In addition, as not all bins are updated every frame, the non-updating bins may need to be adapted in some form to prevent the adaptation from resulting in poorer performance due to a mismatch in the parameter adaptation between the bins that are adapted and the other hm−1 bins that are not adapted in each sub-band. In some systems, this can be done by taking on the updated values of the non-adapted bins to be the same as the neighboring adapted bins.
In either of these systems, the output signal of the acoustic echo canceller module could be sub-sampled by taking from the M1 output samples, the sub-bands at only the indices specified by index hm from Equation 5, thereby achieving the second stage of compression.
This two-stage approach may provide equivalent memory savings and processing load reductions for any further processing downstream of the acoustic echo cancellation. In addition, equivalent processing load reductions may be achieved in the acoustic echo cancellation module in this exemplary scheme as well as memory savings within the acoustic echo cancellation module. This approach may also provide the added advantage of offering significant reduction in perceptual distortion in the reconstructed output signal.
When reconstructing the processed spectrum, the original spectral data (or a representative data set, or a data set of relative measures) is processed so that little or no data is lost when the decompression is complete. By processing the original spectral data (or the representative data or relative measure data set), the sub-band processing system may achieve a lossless or nearly lossless compression. Some systems may preserve almost the entire original spectrum to avoid generating perceivable artifacts when the spectrum is reconstructed.
An overlap-add synthesis may partially reconstruct the spectrum from the processed sparse spectrum. An overlap-add synthesis may avoid discontinuities in the reconstructed spectrum. For each sub-band, the system rotates the remaining bins that made up the sub-band by maintaining relative magnitudes and phases of the original spectrum (or representative data or relative measure data set). The magnitude and phase of the remaining reconstructed bins maintain the same relative magnitude and phase relationship with the restored peak magnitude bin, as the original spectral bins had with the original peak magnitude bin.
Because further processing (e.g., echo cancellation, noise reduction, beam former, signal attenuators, amplifiers, signal modifier, etc.) may alter the magnitude and phase of each sub-band, quantitatively each SBX(m) has been transformed into SBY(m). Equations 9-13 describe how the magnitude and phase for each sub-band may be expanded to its constituent bins. Equation (10) establishes that the magnitude of the restored selected bin is equal to the magnitude of the processed sub-band m. Equation (11) establishes that the phase of the restored selected is equal to the phase of the processed sub-band m after processing. Equations (12) and (13), respectively, establish how the remaining bins may be reconstructed.
for m=0, . . . , M−1 where
Once the complex spectrum is restored, a time domain signal may be generated by an Inverse Fourier Transform device (or function stored in a local or a distributed memory). If windows were used during system analysis, an overlap-add function may be used for synthesis.
Until the spectrum is restored, the original spectrum (or the representative data set) may be retained in a computer readable medium or memory so that the original relative magnitude and phase relationships may be maintained or restored in the decompressed spectrum. This retention potentially reduces audible artifacts that may be introduced by a compression scheme.
The system methods, and descriptions described may be programmed in one or more controllers, devices, processors (e.g., signal processors). The processors may comprise one or more central processing units that supervise the sequence of micro-operations that execute the instruction code and data corning from memory (e.g., computer readable medium) that generate, support, and/or complete an operation, compression, or signal modification. The dedicated applications may support and define the functions of the special purpose processor or general-purpose processor that is customized by instruction code (and in some applications may be resident to vehicles). In some systems, a front-end processor may perform the complementary tasks of gathering data for a processor or program to work with, and for making the data and results available to other processors, controllers, or devices.
The systems, methods, and descriptions may program one or more signal processors or may be encoded in a signal bearing storage medium, a computer-readable medium, or may comprise logic 402 stored in a memory that may be accessible through an interface and is executable by one or more processors 404 as shown in
The memory may retain an ordered listing of executable instructions in a processor, device, or controller accessible medium for implementing logical functions. A logical function may be implemented through digital circuitry, through source code, or through analog circuitry. The software may be embodied in any computer-readable medium or signal-bearing medium, for use by, or in connection with, an instruction executable system, apparatus, and device, resident to system that may maintain persistent or non-persistent connections. Such a system may include a computer system, a processor-based system, or another system that includes an input and output interface that may communicate with a publicly accessible or privately accessible distributed network through a wireless or tangible communication bus through a public and/or proprietary protocol.
A “computer-readable storage medium” “machine-readable medium,” “propagated-signal” medium, and/or “signal-bearing medium” may comprise a medium that stores, communicates, propagates, or transports software or data for use by or in connection with an instruction executable system, apparatus, or device. The machine-readable medium may selectively be, but not limited to, an electronic, magnetic, optical, electromagnetic, infrared, or semiconductor system, apparatus, device, or propagation medium. A non-exhaustive list of examples of a machine-readable medium would include: an electrical connection having one or more wires, a portable magnetic or optical disk, a volatile memory, such as a Random-Access Memory (RAM), a Read-Only Memory (ROM), an Erasable Programmable Read-Only Memory (EPROM or Flash memory), or an optical fiber. A machine-readable medium may also include a tangible medium, as the software may be electronically stored as an image or in another format (e.g., through an optical scan), then compiled, and/or interpreted or otherwise processed. The processed medium may then be stored in a computer and/or machine memory.
The various embodiments presented above are merely examples and are in no way meant to limit the scope of this application. Variations of the innovations described herein will be apparent to persons of ordinary skill in the art, such variations being within the intended scope of the present application. In particular, features from one or more of the above-described example embodiments may be selected to create alternative example embodiments including a sub-combination of features which may not be explicitly described above. In addition, features from one or more of the above-described example embodiments may be selected and combined to create alternative example embodiments including a combination of features which may not be explicitly described above. Features suitable for such combinations and sub-combinations would be readily apparent to persons skilled in the art upon review of the present application as a whole. The subject matter described herein and in the recited claims intends to cover and embrace all suitable changes in technology.
Number | Name | Date | Kind |
---|---|---|---|
5890125 | Davis | Mar 1999 | A |
5999899 | Robinson | Dec 1999 | A |
6501860 | Charrier | Dec 2002 | B1 |
7239253 | Denninghoff | Jul 2007 | B1 |
8457976 | Paranjpe | Jun 2013 | B2 |
8788277 | Vezyrtzis | Jul 2014 | B2 |
9076425 | Yoo | Jul 2015 | B2 |
9255318 | Ono | Feb 2016 | B2 |
10075802 | Kim | Sep 2018 | B1 |
20020006229 | Chao | Jan 2002 | A1 |
20050015249 | Mehrotra | Jan 2005 | A1 |
20050283370 | Mpr | Dec 2005 | A1 |
20060122825 | Oh | Jun 2006 | A1 |
20060140406 | Van DerVeen | Jun 2006 | A1 |
20060248135 | Cousineau | Nov 2006 | A1 |
20060271357 | Wang | Nov 2006 | A1 |
20090313009 | Kovesi | Dec 2009 | A1 |
20100198603 | Paranjpe | Aug 2010 | A1 |
20130257482 | Paranjpe | Oct 2013 | A1 |
20150046171 | Grancharov | Feb 2015 | A1 |
20170092282 | Choo | Mar 2017 | A1 |
20170142412 | Fuchs | May 2017 | A1 |
20180247662 | Fan | Aug 2018 | A1 |
20200013215 | Vosoughi | Jan 2020 | A1 |
20200294519 | Mariglio, III | Sep 2020 | A1 |
Number | Date | Country |
---|---|---|
2755205 | Jul 2014 | EP |
Entry |
---|
Extended European Search Report, EP Application No. 20177030.2 dated Nov. 12, 2020. |
Number | Date | Country | |
---|---|---|---|
20210012785 A1 | Jan 2021 | US |