The present invention relates broadly to a method of digitally filtering an audio signal. The invention relates particularly although not exclusively to digitally filtering an audio signal in audio equalisation (EQ). The invention extends to other digital filtering including filtering images and other signals including signals associated with digital communications and processing.
In digital recording and playback an analog signal representative of audio is converted into a digital signal which lends itself to manipulation and storage. The conversion is performed in an analog to digital converter (ADC). The stored digital signal can be converted back to an analog signal in a digital to analog converter (DAC). The analog signal is played back using conventional audio equipment such as amplifiers and speakers. The digital signal can be manipulated prior to the DAC to improve its quality before playback. This manipulation includes audio EQ where selected parts of the frequency spectrum of the audio are filtered to, for example, compensate for irregularities in the frequency response. The digital signal may also be filtered to resolve problems from its conversion into a digital signal or back to an analog signal.
According to a first aspect of the present invention there is provided a method of digitally filtering an audio signal, said method comprising the steps of:
Preferably the step of increasing the sample rate of the audio filter and/or the other audio filter includes:
Preferably the nominated neighbouring audio signals are shifted in the time domain substantially midway between the neighbouring sample point and the intermediate sample point.
Preferably the step of increasing the sample rate of the audio filter and/or the other audio filter includes:
Preferably the shifted hypothetical audio signal is expanded in the time domain by a factor of substantially two (2).
Preferably the step of combining the audio filters is performed at an adjusted sampling rate wherein the other audio filter includes one or more intervening sample points between adjacent of its neighbouring sample points. More preferably the adjusted sampling rate for applying the audio filter to the other audio filter is inversely proportional to the number of intervening sample points relative to the number of neighbouring sample points for the other filter.
Preferably the step of increasing the sample rate of the audio filter and/or the other audio filter includes:
Alternatively the step of increasing the sample rate of the audio filter and/or the other audio filter includes:
Preferably weighting is applied across a predetermined number of said neighbouring sample points.
Preferably the composite audio filter includes a bank of filters. More preferably the bank of filters together define a frequency bandwidth generally representative of the audio signal to be filtered.
Preferably the composite audio filter is a lowpass or bandpass filter which approaches the Nyquist frequency.
Preferably the audio filter and/or the other audio filter are represented by an impulse response produced by an impulse fed to said filter More preferably the method also comprises the step of applying an averaging curve to the impulse response. Still more preferably the averaging curve is adjusted to a width proportional to the frequency of the impulse response to which it is applied.
Preferably the impulse response is in the time domain represented by a sinc function. Alternatively the impulse response is in the time domain represented by a sine function of absolute values.
According to a second aspect of the invention there is provided a computer or device-readable medium including instructions for digitally filtering an audio signal using a plurality of audio filters each of a predetermined sample rate, said instructions when executed by a processor cause said processor to:
According to a third aspect of the invention there is provided a system for digitally filtering an audio signal, the system comprising:
According to a fourth aspect of the invention there is provided a method of digitally filtering a signal, said method comprising the steps of:
Preferably the signal is an electronic signal derived from displacement of a transducer or measurement device.
According to a fifth aspect of the invention there is provided a method of digitally filtering an image, said method comprising the steps of:
Preferably the image may include a matrix of pixels to which the composite image filter is applied.
According to a sixth aspect of the invention there is provided a method of digitally filtering an audio signal, said method comprising the step of:
In order to achieve a better understanding of the nature of the present invention an embodiment of a method of digitally filtering an audio signal will now be described, by way of example only, with reference to the accompanying drawings in which:
The present invention in a preferred embodiment is directed to a method of digitally filtering an audio signal at a predetermined sample rate by applying a composite audio filter derived at an increased sample rate. The composite audio filter is obtained by combining one audio filter with another audio filter at the increased sample rate. In this embodiment the sample rates of the audio filters may be increased from their predetermined to the increased sample rate by various techniques which involve weighting intermediate sample points.
It will be understood that the various embodiments of the present disclosure can be applied:
Some embodiments of the present disclosure may be embodied in computer program code or software. The digital filter of the digital signal processor 14 is represented by a particular frequency response. The particular frequency response is generally dependent on the impulse response of the filter which is characterised by the software or techniques of the various embodiment of this invention. The present embodiment may cover the basic types of frequency response by which digital filters are classified including lowpass, highpass, bandpass and bandreject or notch filters. The digital filters are broadly categorised as Finite Impulse Response (FIR) or Infinite Impulse Response (IIR) filters.
In order to understand this embodiment of audio filtering with a virtual sample rate increase the composite audio filter is for simplicity derived from two (2) audio filters although it will be appreciated that any number of filters may be used. The composite audio filter generally includes a bank of the filters. The bank of filters together define a frequency bandwidth representative of the audio signal or spectrum to be filtered. In this embodiment an impulse response is produced by an impulse fed to the respective filter. The impulse response for each of the filters may be represented by a sinc function according to the equation:
where lpf is the corner frequency for the lowpass filter, x is the time variable on the x-axis, and e−(qx)2 represents an averaging curve with q representing the aspect ratio of the averaging curve. It is to be understood that the sinc function is the sum of cos components.
In this embodiment the filters are combined by convolution to obtain one of the composite audio filters. This convolution of impulse responses a and b is represented by an array of samples which can also be mathematically defined by the equation:
where N is the number of samples for each of impulse responses a and b, and k is from 0 to N−1 for each of the samples for impulse response b. The array of samples thus includes 2N−1 rows and columns. The sum of the sample values for each row of the array represents the composite audio filter. It is also possible that the composite audio filter is represented mathematically by integrating the impulse responses across an infinite number of samples.
The composite audio filter is in this example a lowpass filter which approaches the Nyquist frequency. The Nyquist frequencies and above are substantially removed in performing the sample rate increase on the various impulse responses. The composite filter or other composite filters may also function as bandpass or bandreject filter depending on the application.
The composite audio filter is “constructed” with the benefit of increased accuracy at the increased sample rate. As shown in
In some embodiments, the sample rate increase on each of the audio filters may be performed by the following two techniques involving:
In weighting values of the impulse response using the shifted neighbouring audio signals, neighbouring impulse responses are nominated from either side of the intermediate sample point to be determined. Each of the nominated neighbouring samples is then shifted in the time domain substantially midway between the neighbouring sample point and the intermediate sample point. In this example the relevant weighting is calculated by summing values which each of the shifted neighbouring impulse responses contribute at the relevant intermediate sample point. This technique is schematically illustrated in
In using this weighting technique, combining of the audio filters is performed at the adjusted sampling rate so that neighbouring sample points for the audio filter align or correspond with at least each of the intervening sample points of the other audio filter to which it is applied. This involves shifting the audio filter at the adjusted sampling rate relative to the other audio filter. For example, if the other audio filter includes intervening sample points located substantially midway between adjacent of its neighbouring sample points, the adjusted sampling rate for applying the filters to one another is substantially half the predetermined sample rate.
The sampling rate is adjusted in this embodiment by convolving every other impulse response. This means the uppermost impulse response of
For a predetermined sample rate of 44.1 kHz the adjusted sampling rate in this example is 22.05 kHz. If the other audio filter includes nine (9) intervening sample points between adjacent of its neighbouring sample points the adjusted sampling rate will be one tenth of the predetermined sample rate. This equates to an adjusted sampling rate of 4.41 kHz for a predetermined sample rate of 44.1 kHz. It is understood that adjusting the sampling rate “corrects” for shifting of the nominated neighbouring sample points in calculating weightings for each of the intermediate sample points. The shift in the nominated neighbouring samples in the time domain is generally proportional to the adjustment in the sampling rate in convolving the audio filters. Thus, a shift in the nominated neighbouring samples midway between neighbouring sample point and the intermediate sample point means an adjustment in the sampling rate by a factor of one-half.
In weighting values of the impulse response using the expanded hypothetical impulse response, the relevant impulse response is effectively replicated as a hypothetical impulse response with its time domain shifted to align with the intermediate sample point to be determined. In some embodiments, the hypothetical and shifted impulse response is expanded in its time domain by factor of substantially 2. In this example the relevant weighting is calculated by summing values for the expanded impulse response at the neighbouring sample points. This technique is schematically illustrated in
In some embodiments, the sample rate increase on each of the audio filters may alternatively be performed by using the following two (2) techniques involving i) a hypothetical audio signal, and/or ii) neighbouring audio signals.
In weighting values of the impulse response using the hypothetical audio signal, the relevant impulse response is effectively replicated with its time domain shifted to align with the intermediate sample point to be determined. The weighting is calculated by summing values for the hypothetical audio signal at the neighbouring sample points and the weighting is a factor inversely proportional to the sum of these values. The relevant weighting or factor is applied to the impulse response of the filter at respective of the intermediate sample points. This technique is schematically illustrated in
where n is the sample number, q represents the aspect ratio of the averaging curve, and lpf is the corner frequency for the lowpass filter.
In weighting values of the impulse response using the neighbouring audio signals, neighbouring impulse responses are nominated from either side of the intermediate sample point to be determined. In this example the relevant weighting is calculated by summing values which each of the nominated neighbouring impulse responses contribute at the relevant intermediate sample point. This technique is schematically illustrated in
In another embodiment the averaging curve applied to the impulse response may be adjusted to a width proportional to the frequency of the impulse response to which it is applied.
In another aspect of the invention an audio filter is provided at an increased sample rate and applied to an audio signal at its predetermined sample rate. The sample rate increase on the audio filter is provided using any one of the weighting techniques described wherein the intermediate sample points are weighted according to the influence of the neighbouring sample points.
Now that several embodiments of the present disclosure have been described it will be apparent to those skilled in the art that the method of digitally filtering an audio signal has at least the following advantages over the prior art:
Those skilled in the art will appreciate that the invention described herein is susceptible to variations and modifications other than those specifically described. For example, the impulse response may be of practically any waveform. If represented by a mathematical equation, the impulse response is not limited to a sinc function but includes other waveforms such as:
The processing of audio signals need not be limited to acoustics but extends to other sound applications including ultrasound and sonar. The invention also extends beyond audio signals to other signals including signals derived from a physical displacement such as that obtained from measurement devices, for example a strain gauge or other transducer which generally converts displacement into an electronic signal. The invention also covers digital filtering of signals associated with digital communications.
The invention in another embodiment is applied to imaging and image filters where, for example, the matrix of pixels in an image is filtered with a composite image filter. In some embodiments, the composite image filter is obtained by combining two (2) or more image filters at an increased sample rate. In increasing the sample rate to include intermediate sample points, these intermediate points are weighted depending on the influence of neighbouring sample points.
All such variations and modifications are to be considered within the scope of the present invention the nature of which is to be determined from the foregoing description.
This application claims priority from U.S. patent application No. 61/805,432 filed on 26 Mar. 2013, the contents of which are to be taken as incorporated herein by this reference. This application is related to and if required claims priority from US patent application Nos. 61/805,406, 61/805,466, 61/805,469, 61/805,449 and 61/805,463 all filed on 26 Mar. 2013, the contents of which are to be taken as incorporated herein by these references. This application is also related to and if required claims priority from U.S. patent application No. 61/819,630 filed on 5 May 2013 and U.S. patent application No. 61/903,225 filed on 12 Nov. 2013, the contents of which are to be taken as incorporated herein by these references.
Filing Document | Filing Date | Country | Kind |
---|---|---|---|
PCT/AU2014/000319 | 3/26/2014 | WO | 00 |
Number | Date | Country | |
---|---|---|---|
61903225 | Nov 2013 | US | |
61819630 | May 2013 | US | |
61805406 | Mar 2013 | US | |
61805466 | Mar 2013 | US | |
61805469 | Mar 2013 | US | |
61805449 | Mar 2013 | US | |
61805463 | Mar 2013 | US | |
61805432 | Mar 2013 | US |