Time domain aliasing cancellation apparatus and signal processing method thereof

• US 5,857,000 A
• Filed: 12/06/1996
• Issued: 01/05/1999
• Est. Priority Date: 09/07/1996
• Status: Expired due to Term

• Current Assignee: National Cheng KUNG University (Government of The Republic of China)
• Original Assignee: National Science Council
• Inventors: Jar-Ferr, Yang, Din-Yuen, Chan
• Primary Examiner(s): Chin, Stephen
• Assistant Examiner(s): GLUCK, JEFFREY WILLIAM
• Application Number:
  US08/759,672
• Time in Patent Office: 760 Days
• Field of Search: 375/240, 370/477, 704/203, 704/204, 704/500, 364/725.03, 364/400.01, 348/395, 348/403, 382/250
• US Class Current: 375/240
• CPC Class Codes: G10L 19/0212 using orthogonal transformation
• CPC Class Codes: G10L 19/022 Blocking, i.e. grouping of samples in time; Choice of analysis windows; Overlap factoring
• CPC Class Codes: G10L 21/0364 for improving intelligibility

Abstract

A time domain aliasing cancellation (TDAC) apparatus and its signal processing method to be used with the AC-3 high-fidelity audio signal compression system of the MPEG-2 international video standard. This invention proposes two preferred embodiments to realize the compression encoding and decoding processes of the TDAC apparatus. The first preferred embodiment employs a data reordering technique to change the TDAC encoding to a discrete cosine transform (DCT), and furthermore, it changes the TDAC decoding to a inverse discrete cosine transform (IDCT). This implementation has the least computational complexity. The second preferred embodiment utilizes data reordering to change the TDAC encoding and decoding into a type IV discrete cosine transformation, and then converts the DCT transformation into a 2nd order infinite impulse filter. The multiplication coefficients in this filter can be fixed to improve the precision and also to reduce the amount of computations. This implementation of the TDAC apparatus has the simplest hardware structure. Both preferred embodiments are suitable for implementation using VLSI technology.

Claims

1. A time domain aliasing cancellation apparatus comprising an encoding device and a decoding device, said encoding device performing time-domain aliasing cancellation encoding for transforming an m-th input signal frame x_m (n) in the time domain into an m-th signal frame X_m (k) in the frequency domain, said decoding device performing time-domain aliasing cancellation decoding for transforming said input frequency-domain signal frame X_m (k) back to a time-domain signal frame x'"'"'_m (n), said time-domain signal frame x_m (n) and x'"'"'_m (n) and frequency-domain signal frame X_m (k) having N terms wherein N is a positive integer number and n, k, and m are integers, said encoding device of the time-domain aliasing cancellation apparatus comprising:

a modified analysis window function unit for multiplying individual terms of the input time-domain signal frame x_m (n) with a modified analysis window function w_E (n), thereby generating a first time sequence s(n) expressed as;

s(n)=x_m (n)w_E (n), wherein ##EQU56## wherein h(n) is the original analysis window function of a coding system;

an encoding unit, rearranging said first time sequence s(n) for generating a second time sequence y(n) with length N, the first N/4 terms of said second time sequence y(n) being composed of the last N/4 terms of said first time sequence s(n), and the last 3N/4 terms of y(n) being composed of the first 3N/4 terms of s(n);

a subtraction unit for subtracting from the N/2 terms of said second time sequence y(n) the terms in the second half of said time sequence in reversed order and generating a third time sequence u(n) expressed as;

u(n)=y(n)-y(N-1-n);

a discrete cosine transform unit, performing a discrete cosine transformation on said third time sequence u(n) for generating a first frequency sequence U(k), wherein k is an integer, and the transformation equation being expressed as ##EQU57## a frequency sequence adder, utilizing said first frequency sequence U(k) for generating a second frequency sequence Y(k) with length N/2 expressed as Y(k)=U(k+1)+U(k); and

an output unit, using said second frequency sequence Y(k) for generating an output encoded frequency sequence X_m (k) with length N, the first N/2 terms of X_m (k) being Y(k) multiplied by a phase factor (-1)_m k and expressed as X_m (k)=(-1)^mk Y(k), and the last N/2 terms of X_m (k) being Y(k) multiplied by a phase factor (-1)^mk+1 in reversed order and expressed as X_m (k)=(-1)^mk+1 Y(N-k-1),

and said decoding device of the time-domain aliasing cancellation apparatus comprising;

an input sign adjustment unit, adding a phase factor of (-1)^mk to the frequency sequence X_m (k) and generating a third frequency sequence Y(k) expressed as Y(k)=(-1)^mk X_m (k);

an adder, utilizing said third frequency sequence Y(k) for generating a fourth frequency sequence Z(k) with length N/2 wherein Z(k)=2Y(k-1)+2Y(k) when k is between 1 to N/2-1 and Z(k)=2Y(0) when k equals zero;

an inverse discrete cosine transformation unit, performing an inverse discrete cosine transformation on said fourth frequency sequence Z(k) for generating a fourth time sequence z(n) with length N/2, the transformation equation being ##EQU58## decoding unit, rearranging said fourth time sequence z(n) for generating a fifth time sequence q_m (n), wherein the first N/4 terms of q_m (n) are composed of the second half of z(n), the second N/4 terms of q_m (n) are composed of the second half of z(n) in reversed order, the third N/4 terms of q_m (n) are the first half of z(n) in reversed order, and the last N/4 terms of q_m (n) are the first half of the fourth time sequence z(n); and

a modified synthesis window function unit for multiplying said fifth time sequence q_m (n) and the previous input time sequence q_m-1 (n) by a modified synthesis window function w_D (n) for generating an output time sequence x'"'"'_m (n) expressed as ##EQU59## with ##EQU60## wherein f(n) is an original synthesis window function.

2. The time domain aliasing cancellation apparatus according to claim 1, wherein said original analysis window function h(n) is equal to said original synthesis window function f(n).

3. A time domain aliasing cancellation method comprising an encoding method and a decoding method, said encoding method performing time domain aliasing cancellation encoding for transforming an m-th signal frame x_m (n) in the time domain into an m-th signal frame X_m (k) in the frequency domain, said decoding method performing time domain aliasing cancellation decoding for transforming said input frequency-domain signal frame X_m (k) back to a time-domain signal frame x'"'"'_m (n), said time-domain signal frames x_m (n) and x'"'"'_m (n) and said frequency-domain signal frame X_m (k) having N terms wherein N is a positive integer number, and n, k, and m are integers, said encoding method of the time-domain aliasing cancellation method comprising the steps of

multiplying individual terms of said input time-domain signal frame x_m (n) by a modified analysis window function w_E (n) for generating a first time sequence s(n) with length N expressed as:

s(n)=x_m (n)w_E (n) with ##EQU61## wherein h(n) is the original analysis window function;

rearranging said first time sequence s(n) for generating a second time sequence y(n) with length N, the first N/4 terms of said second time sequence y(n) being composed of the last N/4 terms of said first time sequence s(n), and the last 3N/4 terms of y(n) being composed of the first 3N/4 terms of s(n);

subtracting from the first N/2 terms of said second time sequence y(n) the terms in the second half of said time sequence in reversed order for generating a third time sequence u(n) expressed as;

u(n)=y(n)-y(N-1-n);

performing a discrete cosine transformation on said third time sequence u(n) for generating a first frequency domain signal U(k) expressed as;

##EQU62## adding said first frequency sequence U(k) with neighboring terms thereof for generating a second frequency sequence Y(k) with length N/2 expressed as;

Y(k)=U(k+1)+U(k); and

rearranging said second frequency sequence Y(k) for generating said encoded output frequency sequence X_m (k) with length N, the first N/2 terms of X_m (k) being Y(k) multiplied by a phase factor (-1)^mk expressed as;

X_m (k)=(-1)^mk Y(k), the last N/2 terms of X_m (k) being Y(k) in reversed order and multiplied by a phase factor (-1)^mk+1 expressed as;

X_m (k)=-(1)^mk+1 Y(N-k-1),

and said decoding method comprising the steps of

adding a phase factor of (-1)^mk to said frequency sequence X_m (k) for generating a third frequency sequence Yr(k) expressed as Yr(k)=(-1)^mk X_m (k);

utilizing said third frequency sequence Yr(k) for generating a fourth frequency sequence Z(k) with length N/2, wherein Z(k)=2Yr(k-1)+2Yr(k) when k is between 1 to N/2-1, and Z(k)=2Yr(0) when k is equal to zero;

performing an inverse discrete cosine transform on said fourth frequency sequence Z(k) for generating a fourth time sequence z(n) expressed as ##EQU63## rearranging said fourth time sequence z(n) for generating a fifth time sequence q_m (n) with length N, wherein the first 1/4 of the terms of q_m (n) being composed of the second half of z(n), the second N/4 terms of q_m (n) being composed of the second half of z(n) in reversed order, the third N/4 terms of q_m (n) being the first half of z(n) in reversed order, and the last N/4 terms of q_m (n) being the first half of the fourth time sequence z(n); and

multiplying said fifth time sequence q_m (n) and the previous frame input time sequence q_m-1 (n) by a modified synthesis window function W_D (n) for generating the output time sequence x (n) expressed as ##EQU64## wherein f(n) being the original analysis window function of a coding system.

4. A time domain aliasing cancellation apparatus comprising an encoding device and a decoding device, said encoding device performing time domain aliasing cancellation encoding for transforming an m-th signal frame x_m (n) in the time domain into an m-th signal frame X_m (k) in the frequency domain, said decoding device performing time domain aliasing cancellation decoding for transforming said input frequency-domain signal frame X_m (k) back to a time-domain signal frame x (n), said time-domain signal frames x_m (n) and x'"'"'_m (n) and said frequency-domain signal frame X_m (k) having N terms wherein N is a positive integer number and n, k, and m are integers, said encoding device of the time-domain aliasing cancellation apparatus comprising:

a modified analysis window function unit for multiplying term by term the input time-domain signal frame x_m (n) with a modified analysis window function w_E (n) for generating a first time sequence s(n) expressed as s(n)=x_m (n)w_E (N-1-n) with ##EQU65## wherein h(n) is the original analysis window function of a coding system 0<

=J<

=N/2-1;

an encoding unit, rearranging said first time sequence s(n) for generating a second time sequence y(n) with length N, the first N/4 terms of said second time sequence y(n) being composed of the last N/4 terms of said first time sequence s(n), and the last 3N/4 terms of y(n) being composed of the first 3N/4 terms of s(n);

a subtraction unit, subtracting from the N/2 terms of said second time sequence y(n) the terms in the second half of said time sequence in reversed order for generating a third time sequence v(n);

a first buffer register with N/2 random access memory registers for storing said third time sequence signal v(n);

a first address unit, using a first parameter as the address to select said third time sequence v(n) from said first buffer, rearranging said third time sequence v(n) for generating a fourth time sequence v'"'"'(n);

a first sign adjustment unit, using a second parameter for adjusting the sign of each term of said fourth time sequence v'"'"'(n);

a first digital filter for transforming said sign-adjusted fourth time sequence v'"'"'(n) into a first encoded frequency sequence Y(k), said first encoded frequency sequence Y(k) being the discrete cosine transformation of said fourth time sequence v'"'"'(n); and

an output unit, using said first frequency sequence Y(k) for generating said output encoded frequency X_m (k), wherein the first N/2 terms of X_m (k) are Y(k) multiplied by a phase factor (-1)^mk expressed as X_m (k)=(-1)^mk Y(k), and the last N/2 terms of X_m (k) being Y(k) multiplied by a phase factor (-1)^mk+1 expressed as X_m (k)=(-1)^mk+1 Y(k),

said decoding apparatus comprising

an input sign adjusting unit, using said input frequency sequence X_m (k) with length N by a phase factor (-1)^mk, and shifting the result to the left by 1 bit to perform a multiplication by 2 for generating a second frequency sequence 2Y(k) expressed as 2Y(k)=2(-1)^mk X_m (k);

a second buffer with N/2 random access memory registers for storing said second frequency sequence 2Y(k);

a second address unit, using said first parameter as the address to select each term of said second frequency sequence 2Y(k) from said second buffer, rearranging for generating a third frequency sequence Y'"'"'(k);

a second sign adjustment unit, using said second parameter for correcting the sign of each term of said third frequency sequence Y'"'"'(k);

a second digital filter, transforming said sign-adjusted third frequency sequence Y'"'"'(k) into a fifth time sequence y(n), said fifth time sequence y(n) being the inverse discrete cosine transformation of said third frequency sequence Y'"'"'(k);

a reordering encoder, rearranging said fifth time sequence y(n) for generating a sixth time sequence q_m (n) with length N, the first 3N/4 terms of said sixth time sequence q_m (n) being composed of the last 3N/4 terms of said fifth time sequence y(n), and the last N/4 terms of said sixth time sequence q_m (n) being composed of the negative of the first N/4 terms of said fifth time sequence y(n); and

a modified synthesis window function unit for multiplying said fifth time sequence q_m (n) and the previous input time sequence q_m-1 (n) by a modified synthesis window function w_D (n) for generating the output time sequence x'"'"'_m (n) expressed as ##EQU66## wherein ##EQU67## and f(n) is an original synthesis window function.

5. The time domain aliasing cancellation apparatus according to claim 4, wherein each of said first and second digital filters receives an input signal for generating an output signal, and each digital filter includes:

a first adder, adding said received input signal and a fifth internal signal for generating a first internal signal;

a first delay unit, receiving said first internal signal and inserting a time delay for generating a second internal signal;

a second delay unit, receiving said second internal signal and inserting a time delay for generating a third internal signal;

a constant coefficient multiplier, receiving said second internal signal and multiply with a constant coefficient for generating a fourth internal signal;

a second adder, adding said third internal signal and said fourth internal signal together for generating said fifth internal signal; and

a third adder, adding said first internal signal and said second internal signal for generating said output signal.

6. A time domain aliasing cancellation signal processing method comprising an encoding method and a decoding method, said encoding method performing time domain aliasing cancellation encoding for transforming an m-th sequence x_m (n) in the time domain into an m-th sequence X_m (k) in the frequency domain, said decoding method performing time domain aliasing cancellation decoding for transforming said input frequency-domain sequence X_m (k) back to an output sequence x'"'"'_m (n), said time-domain sequences x_m (n) and x_m '"'"'(n) and said frequency-domain sequence X_m (k) having N terms wherein N is a positive integer number, and n, k, and m are integers, said encoding method of the time domain aliasing cancellation method comprising the steps of:

multiplying a modified analysis window function w_E (n) term by term with said input time domain signal frame x_m (n) for generating a first time sequence s(n) expressed as s(n)=x_m (n)w_E (N-1-n) with ##EQU68## wherein h(n) being an original analysis window function, 0<

=J<

=N/2-1;

rearranging said first time sequence s(n) for generating a second time sequence y(n) with length N, the first N/4 terms of said second time sequence y(n) being composed of the last N/4 terms of said first time sequence s(n), and the last 3N/4 terms of y(n) being composed of the first 3N/4 terms of s(n);

subtracting from the first N/2 terms of said second time sequence y(n) the terms in the second half of said time sequence in reversed order for generating a third time sequence v(n);

storing said third time sequence signal v(n) in a first buffer with N/2 registers;

using a first parameter as the address to select said third time sequence v(n) from said first buffer, rearranging said third time sequence v(n) for generating a fourth time sequence v'"'"'(n);

using a second parameter for adjusting the sign of each term of said fourth time sequence v'"'"'(n);

filtering and transforming said sign-adjusted fourth time sequence v'"'"'(n) into a first encoded frequency sequence Y(k), said first encoded frequency sequence Y(k) being the discrete cosine transformation of said fourth time sequence v'"'"'(n); and

reordering said first frequency sequence Y(k) for generating said output encoded frequency X_m (k), wherein the first N/2 terms of X_m (k) being Y(k) multiplied by a phase factor (-1)^mk expressed as

space="preserve" listing-type="equation">X.sub.m (k)=(-1).sup.mk Y(k)

wherein the last N/2 terms of X_m (k) being Y(k) multiplied by a phase factor (-1)^mk+1 expressed as

space="preserve" listing-type="equation">X.sub.m (k)=(-1).sup.mk+1 Y(k),

said decoding procedure comprising the steps of

multiplying said input frequency sequence X_m (k) with length N by a phase factor (-1)^mk, and shifting the result to the left by 1 bit to perform a multiplication by 2 for generating a second frequency sequence 2Y(k) expressed as 2Y(k)=2(-1)^mk X_m (k);

storing said 2nd frequency sequence 2Y(k) in a second buffer with N/2 RAM registers;

using said first parameter as the address to select each term of said second frequency sequence 2Y(k) from said second buffer registers and rearrange order thereof for generating a third frequency sequence Y'"'"'(k);

using said second parameter for correcting the sign of each term of said third frequency sequence Y'"'"'(k);

filtering and transforming said sign-adjusted third frequency sequence Y'"'"'(k) into a fifth time sequence y(n), said fifth time sequence y(n) being the inverse discrete cosine transformation of said third frequency sequence Y'"'"'(k);

rearranging said fifth time sequence y(n) for generating a sixth time sequence q_m (n) with length N, the first 3N/4 terms of said sixth time sequence q_m (n) being composed of the last 3N/4 terms of said fifth time sequence y(n), and the last N/4 terms of said sixth time sequence q_m (n) being composed of the negative of the first N/4 terms of said 5th time sequence y(n);

multiplying said fifth time sequence q_m (n) and the previous input time sequence q_m-1 (n) by a modified synthesis window function W_D (n) for generating the output time sequence x_m '"'"'(n) expressed as ##EQU69## wherein ##EQU70## and f(n) is an original synthesis window function.

7. The signal processing method according to claim 6, wherein each of the filtering steps further comprises the steps of:

receiving said input signal, adding it to a fifth internal signal for generating a first internal signal;

receiving said first internal signal and inserting a time delay for generating a second internal signal;

receiving said second internal signal and inserting a time delay for generating a third internal signal;

receiving said second internal signal and multiplying with a constant coefficient for generating a fourth internal signal;

receiving said third and fourth internal signals, add these two signals together for generating said fifth internal signal; and

receiving said first and second internal signal, add these two signals together for generating said output signal.

8. The signal processing method according to claim 6, further comprising a procedure for generating said first parameter, comprising the steps of:

letting n denote said first parameter, wherein n satisfies the condition expressed as

space="preserve" listing-type="equation">(2k+1)n|.sub.mod N =(2J+1) n+J-k |.sub.mod N

wherein J is the selected index corresponding to said constant multiplicative coefficient (2K+1) associated with said first parameter;

storing the value of the left-hand-side and the right-hand-side of the above equation in a left accumulator and right a accumulator respectively;

fixing the value of n in said right accumulator while maintaining the value of said right accumulator positive;

incrementing n from zero, until the value of the content of said right accumulator equals the value of the content of said left accumulator;

if said intermediate solution n being less than or equal to N/2-1, then said first parameter is equal to n; and

if said intermediate solution n is greater than N/2-1, then said first parameter is set to N-1-n.

9. The signal processing method according to claim 8, further comprising a procedure for generating said second parameter, comprising the steps of:

using a left counter and a right counter to register the carry bits of said left accumulator and said right accumulator, respectively; and

if said intermediate solution n is less than or equal to N/2-1, then said second parameter is the result of the XOR of the contents of said left counter and said right counter; and

if said intermediate solution n is greater than N/2-1, then said second parameter is the result of the XOR of the contents of said right counter and said left counter.