US20040147229A1

US20040147229A1 - High frequency signal construction method and apparatus

Info

Publication number: US20040147229A1
Application number: US10/473,800
Authority: US
Inventors: David McGrath
Original assignee: Lake Technology Ltd
Current assignee: Dolby Laboratories Licensing Corp
Priority date: 2001-04-10
Filing date: 2002-04-10
Publication date: 2004-07-29
Also published as: WO2002084885A1; US7685218B2; AUPR433901A0; US20070098185A1

Abstract

A method of adding high frequency content to an input signal to form an augmented signal, the method comprising the steps of: (a) providing an initial signal having a first predetermined lower spectral range; (b) utilising the initial signal to form synthesised high frequency components of the initial signal; (c) filtering the initial signal with a low pass filter and filtering the synthesised high frequency components with a high pass filter (d) combining the filtered signals to form the augmented signal.

Description

BACKGROUND

1. Field of the Invention

The present invention relates to the synthesis of high frequency signals and, in particular, discloses a method and system for synthesising high frequency audio signals.

2. Background of the Invention

The digital recording of audio signals has become extremely popular. The most popular format for recording is the CD audio format which samples a signal at approximately 44.1 KHz. This is likely to produce a corresponding audio range of approximately 20 kHz which was thought to be adequate for reproducing the audio range that the human ear can detect. However, it is thought by some that the human ear is able to colour an audio signal through the utilization of portions of a signal above 20 kHz. Hence, recent standards have proposed either an 88.2 or a 96 kHz sampling rate. There is therefore the significant problem of how one takes, for example, a 44.1 kHz recorded signal and remasters the signal in say an 88.2 kHz format. One standard technique utilised is to use an interpolator that also uses some kind of linear filter to perform an anti alias filtering operation.

For the purposes of further discussion, the following terminology is defined:

The original signal is called the Original Audio Signal.

The original audio sample rate is called the Original Sample Rate.

The original audio signal is believed to be “accurate” up to a frequency known as the Original Frequency Range.

The Original Half Nyquist Frequency is defined as 0.5 times the Original Sample Rate.

The interpolated signal is called the Interpolated Audio Signal.

The new (higher) audio sample rate is called the Interpolated Sample Rate.

The Interpolated Half Nyquist Frequency is defined as 0.5 times the Interpolated Sample Rate.

The Oversampling Ratio is the Interpolated Sample Rate divided by the Original Sample Rate.

Typical values of the above defined quantities are:

For a CD player with 4× oversampling D/A converters

Original Sample Rate=44,100 Hz

Original Frequency Range=20,000 Hz

Original Half Nyquist Frequency=22,050 Hz

Interpolated Sample Rate=176, 400 Hz

Interpolated Half Nyquist Frequency=88, 200 Hz

In a system like this, the Original Audio Signal only contains reliable content up to 20 kHz, but it is assumed it may be desirable to synthesize new high frequency content up to say 88.2 kHz.

For a DVD player with 2× oversampling D/A converters:

Original Sample Rate=48,000 Hz ‘Original Frequency Range=20,000 Hz Original Half Nyquist Frequency=24,000 Hz

Interpolated Sample Rate=96,000 Hz Interpolated Half Nyquist Frequency=48,000 Hz

In a system like this, the Original Audio Signal only contains reliable content up to 20 kHz, but it may be desirable to synthesize new high frequency content up to 48 kHz. [ 0024] The standard prior art anti-aliasing approach to higher sampling rate extension operates on the principle that as no information about what audio content may have existed above the Original Half Nyquist Frequency is provided in the original audio material, it is necessary to ensure that an Interpolated Audio Signal has zero content in this upper frequency range.

The standard prior art method for producing an interpolated signal will now be described. Turning initially to FIG. 1, an original audio signal 1 is provided having samples e.g., 11,12. The samples are assumed to have been provided at a standard rate. The first step in forming the interpolated signal is to zero pad the audio signal as illustrated in FIG. 2. In zero padding, zero value signals e.g., 14,15 are added to the signal between samples. Next, as illustrated in FIG. 3, an interpolation process is provided where the signal e.g., 18 is formed from an interpolation of the two

signals

17,19. In the example provided, the interpolated sample rate is twice the original sample rate and hence the over sampling ratio is 2 with one zero sample inserted between each sample of the original audio signal. The zero-padding technique results in aliasing, meaning that the low frequency audio signal is duplicated in higher frequency bands. These higher frequency replicas (called aliases) are then filtered out (using a low-pass filter), to leave the Interpolated Audio Signal.

An example of aliasing is illustrated in FIG. 4 where an original audio signal having a

frequency spectrum

21 is zero padded resulting in the zero padded audio signal having a

frequency spectrum

23,24 with the lower frequency being replicated in high frequency bands. The interpolation process is equivalent to applying a low-pass filter 27 which results in the interpolated audio signal 29 which substantially reflects the original audio signal 21.

The arrangement of the prior art has a significant disadvantage in that none of the high frequency spectrum is utilized when a re-sampling occurs.

SUMMARY

It is an object of the present invention is to provide for alternative forms of high frequency signal extension of signals.

According to a first aspect of the invention there is provided method of adding high frequency content to an input signal to form an augmented signal, the method comprising the steps of:

(a) providing an initial signal having a first predetermined lower spectral range;

(b) utilising said initial signal to form synthesised high frequency components of said initial signal which extend beyond said lower spectral range;

(c) filtering said initial signal with a low pass filter and filtering said synthesised high frequency components with a high pass filter

(d) combining said high and low pass filtered signals to form said augmented signal.

Preferably, step (b) further comprises:

(i) for at least one portion of the input signal, determining the spectral content of said portion;

(ii) extrapolating a high frequency end portion of the spectral content to form said synthesised higher frequency components of said signal.

Conveniently, said portion is multiplied with a window function prior to determination of the spectral content and said synthesised higher frequency components are summed in an overlap-add fashion.

Advantageously, the method includes the step of dividing the input signal into a plurality of overlapping blocks, with each block being multiplied by a sliding window function to yield a series of windowed portions from which high frequency components are successively synthesised.

The window may be of a Gaussian or Hanning form.

The invention extends to a method of adding high frequency content to an input signal to form an augmented signal, the method comprising the steps of:

(b) utilising said initial signal to form synthesised high frequency components of said initial signal which extend beyond said lower spectral range, wherein step (b) further comprises:

(i) for at least one portion of the input signal, determining the spectral range of said portion;

(ii) extrapolating a high frequency end portion of the spectral range to form said synthesised higher frequency components of said signal.

Conveniently, said portion is multiplied with a window function prior to determination of the spectral content, and said synthesised high frequency components are summed in an overlap-add fashion.

Preferably, the method includes the steps of dividing the input signal into a plurality of overlapping blocks, with each block being multiplied by a sliding window function to yield a series of windowed portions from which high frequency components are successively synthesised.

Conveniently, at least some of the highest frequency components of said spectral content are discarded prior to the extrapolation of the remaining high frequency components.

Typically, the step of extrapolating said high frequency end portion comprises the steps of sampling the high frequency components, defining an extrapolation factor based on a geometric progression, and generating said geometric progression on the basis of the sampled high frequency components.

According to still further aspect of the invention there is provided apparatus for adding high frequency content to an input signal to form an augmented signal, the apparatus comprising:

(a) a synthesising processor for synthesising high frequency components from an initial signal having a first predetermined lower spectral range, said high frequency components extending beyond said lower spectral range;

(b) a low pass filter for filtering said initial signal;

(c) a high pass filter for filtering said synthesised high frequency components;

(d) a combiner for combining said high and low pass filtered signals to form said augmented signal.

Preferably, said synthesising processor comprises means for determining the spectral content of at least one portion of said input signal and means for extrapolating from a high frequency end portion of said spectral content to form said synthesised high frequency components of said signal.

The invention still further provides an apparatus for adding high frequency content to an input signal to form an augmented signal, the apparatus comprising:

(b) means for dividing the input signal into a plurality of overlapping portions;

(c) means for determining the spectral content of each of said overlapping portions; (d) means for extrapolating the high frequency end portion of the spectral content to form said synthesised high frequency components of said signal; (e) means for summing said synthesised high frequency components in a overlap-add fashion.

BRIEF DESCRIPTION OF THE DRAWINGS

Notwithstanding any other forms which may fall within the scope of the present invention, preferred forms of the invention will now be described, by way of example only, with reference to the accompanying drawings in which: [0060]
FIG. 1 illustrates a sampled original audio signal; [0061]
FIG. 2 illustrates a zero padded audio signal; [0062]
FIG. 3 illustrates an interpolated audio signal; [0063]
FIG. 4 illustrates the prior art process of forming an interpolated audio signal in the frequency domain; [0064]
FIG. 5 illustrates the basic process of forming an augmented audio signal in accordance with one embodiment; [0065]
FIG. 6 illustrates the frequency extension process of FIG. 5; [0066]
FIG. 7 illustrates a first embodiment of an apparatus for generating an augmented audio signal; and [0067]
FIG. 8 illustrates a second embodiment of an apparatus for generating an augmented audio signal.[0068]

DESCRIPTION OF PREFERRED AND OTHER EMBODIMENTS

In the preferred embodiment, there are provided various techniques for creating a reasonable estimate of the frequency response of the audio signal above the original frequency range. Further, the techniques are extended to include techniques for incorporating the extended frequency response signal into the interpolated audio signal whilst ensuring that the interpolated audio signal is an accurate match to the original signal in the more important lower frequency range. [0069]
All frequency extension techniques are, by definition, non-linear, because they cause the creation of new frequency content in the output signal that was not present in the input signal. Hence it is extremely difficult to ensure that a Frequency Extension Technique does not also introduce non-linear/distortion artifacts that are audible within the Original Frequency Range of the Original Audio Signal. Hence, the preferred embodiment proposes that the information from the Original Audio Signal, within the Original Frequency Range, should be preserved, by reinserting it into the Interpolated Audio Signal. [0070]
An example of this arrangement is shown at [0071] 30 in FIG. 5 where an original audio signal having a spectrum 21 is adapted utilizing a frequency extension technique 36 which is described in more detail below so as to provide for an extended audio signal having an extended frequency 32. The original signal is low-pass filtered using a low pass filter 37 and the extended audio signal is high-pass filtered using a high pass filter 38 before they are combined at 39 to produce the interpolated or augmented audio signal 34,35 which extends into the high frequency range. Hence, the Interpolated Audio Signal is composed of two signal components added together:
I. The low frequency part of the Original Audio Signal. [0072]
II. The high frequency part of the Extended Audio Signal. [0073]
In many cases, the Extended Audio Signal will be a very close approximation to the Interpolated Audio Signal, but the use of the low-[0074] pass 37 and high-pass 38 filters, and the summing element 39, ensure that any inaccuracies in the low frequency part of the Extended Audio Signal are removed, and replaced with the more accurate low-frequency components from the Original Audio Signal.
The use of the low-pass/high-pass technique, has the following benefits: [0075]
I. Low frequency information (from the Original Audio Signal) is preserved in an unaltered form; [0076]
II. High frequency information approximating the likely extension of the Original Audio Signal is added without affecting the low frequency information, but still using the Original Audio signal as a basis for the extension. [0077]
One method for [0078] high frequency extension 36 can operate by working on a sliding window on the Original Audio Signal, so that, with each iteration of the process, a windowed segment of the Original Audio Signal is analyzed, say in the Fourier domain, resulting in the method as is shown in FIG. 6, which shows a single iteration.
A segment or block [0079] 42.1 of the original audio signal 42 is multiplied with a window 43 (which can be of Gaussian form). The multiplied result of the two signals is then transferred into the Fourier domain using a Fast Fourier Transform (FFT) 44 or the like so as to produce a frequency response 41. The frequency response curve 41 will often include a ringing peak 46 that appears due to the anti-alias filter. However, the lower frequency points 47 and 48 adjacent the peak but still at the high frequency end of the response can be relied upon as truer indicators of the high frequency content of the original audio signal. Hence an extrapolation process 50 can be carried out so as to extend the representative high frequency components 47 and 48 of the audio signal. The components at points 47 and 48 are extrapolated, thereby yielding a reasonable estimate of the extended audio signal 49. The extended audio signal 49 then undergoes an inverse fast Fourier transform 60 before being multiplied by a Gaussian window 51 to yield a partially computed output audio signal 52 in the time domain.
Obviously, various other extrapolation techniques can be utilized. For example, in a 32 tap FFT filter, the FFT bins ranging from 13.5,15,16.5,18 to 19.5 kHz can be used. The 19.5 kHz bin may be adversely influenced by the [0080] peak 46. One form of extension can be made by extrapolating the difference between the samples 47 and 48 corresponding, say, to the 16.5 and 18 kHz bins, to higher frequencies, and by continuing them in a geometric series, as outlined in the mathematical summary below.
A high frequency audio signal augmentation system of the type schematically illustrated in FIG. 7 can be utilized. [0081]
Let x(k) be the original input signal and y(l), the Extrapolated Signal to be created. [0082]
The oversampling ratio can be S (typically, S=2 or S=4). This is implemented using an oversampler [0083] 54, into which the original audio signal is inputted.
Let the original FFT length be N; hence, the extended FFT length, N′=N.S [0084]
Let the two extrapolation FFT bins into which the representative high frequency samples are loaded (say 16.5 and 18 kHz) be defined as e[0085] ₁and e₂
Let the Overlap be L [0086]
Process the input signal blocks in overlapping blocks as follows, using the [0087] segmenting function 56. Each of the overlapping blocks are successively multiplied with a Gaussian window function 43 using multiplier 58. The forward shift in the window function essentially defines the blocks and their degree of overlap.
In iteration p: [0088]
Take the windowed input block of length N b[0089] _p(i)=x(L.p+I).w(i) for i=0 . . . N−1
Take the (real) FFT of this input block B[0090] _p=FFT[b_p] using FFT processor 44 (note, b_pis of length N, but B_pis of length N/2+1, because we are using the real FFT)
Define the Extrapolation factor f that tells us how the frequency response of B[0091] _pcan be extrapolated beyond the FFT bins e₁and e₂, by a geometric progression. The Extrapolation factor f is a complex number, constrained to lie on or within the unit circle:
f=B _p(e ₂)/B _p(e ₁) if |B _p(e ₂ |<|B _p(e ₁)|
f=0 if B _p(e ₁)=0
f=(B _p(e ₂)/B _p(e ₁))/|B _p(e ₂)/B _p(e ₁) otherwise
This is achieved using a [0092] frequency extrapolation processor 50.
Form the new, extended frequency response B′[0093] _p(which is of length N′/2+1) defined as:
B′ _p(i)=B _p(i) (0≦i≦e ₁)
B′ _p(i)=B _p(e ₁).f ^j−e1(e ₁ <i≦N′/2+1)
Transform this extended frequency response back to the time domain using [0094] Inverse FFT processor 60, creating a time-domain signal block of length N′:
b′_p=IFFT{B′_p}
This block of output is then summed into an output buffer after applying a suitable window w′ which can be in the form of the [0095] Gaussian window 51. The summer and buffer are shown at 62. Many other different window functions can be used with this method, with one desirable window including a Hanning window.
y _p(S.L.p+i)=y _p−1(S.L.p+i)+b′ _p(i).w′(i) for i=0 . . . N′−1
Following this summation operation, the first S.L samples are available to be output: [0096]
y(S.L.p+i)=yp(S.L.p+i) for i=0 . . . S.L−1
Part of the oversampled input audio signal is low pass filtered using the [0097] low pass filter 37 having a cut-off frequency of 19 kHz, towards the end of the audible frequency range. The summed extrapolated time domain samples are high pass filtered using the high pass filter 38, which has a cut-off frequency of 19 kHz, matching that the low pass filter 37 so as to prevent overlap of the low (audible) and high (ultrasonic) frequency parts of the signal, thereby to prevent lower frequency components of the ultrasonic signal interfering with the higher frequency components of the audible signal. These are then summed at summer 39 to yield an extrapolated or augmented output audio signal.
Referring now to FIG. 8, a further embodiment of a signal augmenting apparatus is shown which is specific to a doubling of the sampling rate. One part of the audio signal is processed through a sample rate converter or [0098] oversampler 65 which includes a 19 kHz low pass filter, with the low frequency output portion being fed to the summer 39. The other part of the input audio signal is fed to a segmenter 67 where it is broken into overlapping blocks. Each block is half as long as the equivalent blocks in FIG. 7 for the reason that the audio data has not been sample-rate converted or oversampled. Each block is in turn multiplied with the Gaussian window 43 and the result is converted into the frequency domain using a fast Fourier transform function 66 of length N, half that of the function 44 in FIG. 7. The frequency response is then processed using frequency extrapolator or interpolator 68, resulting in an extended audio signal of the type illustrated at 49 in FIG. 6. The augmented signal is inverse fast Fourier transformed back into the time domain using an inverse fast Fourier transform function 60. The resultant time domain signal has a length 2N. From then on, the processor is identical to that illustrated in FIG. 7.
It will be appreciated that in both FIGS. 7 and 8 each successive block defined by a shift in the window undergoes the process illustrated in FIG. 6 to yield a succession of partially computed output audio signals [0099] 52 which are then buffered and summed before being high pass filtered at 38. The main difference is that in FIG. 8 the frequency extrapolator 68 both doubles the length of the shorter frequency vector at the same time as extrapolating it.
In certain forms of the invention, the high and low pass filters may be eliminated. By way of example, with reference to FIG. 7 one embodiment may exclude [0100] low pass filter 37, high pass filter 38 and summer 39, with the extrapolation technique being sufficient to avoid corruption of the audible frequencies.
In a further possible embodiment, blocks [0101] 34, 37, 38 and 39 may be included, and the remaining blocks may be replaced by a crude interpolation function. In this embodiment, the high and low pass filters serve the primary function of preventing the high frequency portion of the signal from corrupting the low frequency portion.
The invention has numerous audio recordal and playback applications, including the following: [0102]
remastering of digital and analogue recordings having a relatively low sampled rate in the region of 44 kHz; [0103]
processing of recorded audio signals in CD, DVD and similar players having oversampling functions, where the audio content is coloured with but not altered by ultrasonic components; [0104]
particular application in audio playback devices such as CD and DVD players, and in any similar devices where oversampling is utilized. [0105]
It would be appreciated by a person skilled in the art that numerous variations and/or modifications may be made to the present invention as shown in the specific embodiment without departing from the spirit or scope of the invention as broadly described. The present embodiment is, therefore, to be considered in all respects to be illustrative and not restrictive. [0106]

Claims

1. A method of adding high frequency content to an input signal to form an augmented signal, the method comprising the steps of:

2. A method as claimed in claim 1 wherein said step (b) further comprises:

3. A method as claimed in claim 2 wherein said portion is multiplied with a window function prior to determination of the spectral content and said synthesised higher frequency components are summed in an overlap-add fashion.

4. A method as claimed in claim 3 which includes the step of dividing the input signal into a plurality of overlapping blocks, with each block being multiplied by a sliding window function to yield a series of windowed portions from which high frequency components are successively synthesised.

5. A method as claimed in claims 3 or 4 wherein said window is of a Gaussian or Hanning form.

6. A method of adding high frequency content to an input signal to form an augmented signal, the method comprising the steps of:

7. A method according to claim 6 wherein said portion is multiplied with a window function prior to determination of the spectral content, and said synthesised high frequency components are summed in an overlap-add fashion.

8. A method according to claim 7 which includes the steps of dividing the input signal into a plurality of overlapping blocks, with each block being multiplied by a sliding window function to yield a series of windowed portions from which high frequency components are successively synthesised.

9. A method as claimed in any one of claims 2 to 8 wherein at least some of the highest frequency components of said spectral content are discarded prior to the extrapolation of the remaining high frequency components.

10. A method as claimed in any one of claims 6 to 9 in which the step of extrapolating said high frequency end portion comprises the steps of sampling the high frequency components, defining an extrapolation factor based on a geometric progression, and generating said geometric progression on the basis of the sampled high frequency components.

11. A method as claimed in any previous claim wherein said signal is an audio signal, said lower spectral range corresponds to an audible component of said signal and said high frequency components correspond to ultrasonic components.

12. A method of adding high frequency content to an input signal substantially as hereinbefore described with reference to FIGS. 5 to 8 of the accompanying drawings.

13. Apparatus for adding high frequency content to an input signal to form an augmented signal, the apparatus comprising:

(b) a low pass filter for filtering said initial signal;

14. An apparatus according to claim 13 in which said synthesising processor comprises means for determining the spectral content of at least one portion of said input signal and means for extrapolating from a high frequency end portion of said spectral content to form said synthesised high frequency components of said signal.

15. An apparatus for adding high frequency content to an input signal to form an augmented signal, the apparatus comprising:

(c) means for determining the spectral content of each of said overlapping portions;

(d) means for extrapolating the high frequency end portion of the spectral content to form said synthesised high frequency components of said signal;

(e) means for summing said synthesised high frequency components in a overlap-add fashion.

16. Apparatus for adding high frequency content to an input signal substantially as herein before described with reference to FIGS. 5 to 8 of the accompanying drawings.