Signaling method in an OFDM multiple access system

Abstract

A method for reducing the peak-to-average ratio in an OFDM communication signal is provided. The method includes defining a constellation having a plurality of symbols, defining a symbol duration for the OFDM communication signal, and defining a plurality of time instants in the symbol duration. A plurality of tones are allocated to a particular communication device, and a discrete signal is constructed in the time domain by mapping symbols from the constellation to the time instants. A continuous signal is generated by applying an interpolation function to the discrete signal such that the continuous signal only includes sinusoids having frequencies which are equal to the allocated tones.

Claims

1. A method of wireless communication, comprising: mapping received data symbols to a discrete signal of complex data symbols; receiving the discrete signal of complex data symbols; calculating, based on the received discrete signal of complex data symbols, frequency responses corresponding to tones through an M-point DFT; outputting a frequency response vector of length M; expanding the frequency response vector of length M to a new vector of length N by inserting zero-value symbols corresponding to tones other than tones allocated to a transmitter; performing an N-point inverse discrete Fourier transform operation on the new vector of length N after zero insertion to obtain a vector of digital signal samples; and pre-pending a cyclic prefix to the vector of digital signal samples.

2. The method of claim 1, further comprising generating a digital-to-analog converted OFDM-based signal using the vector of digital signal samples with the pre-pended cyclic prefix for transmission.

3. The method of claim 1, wherein the tones are allocated tones available for transmission.

4. The method of claim 1, wherein the allocated tones are contiguous, and wherein the new vector of length N is larger than the frequency response vector of length M, and wherein the allocated tones are equally spaced.

5. The method of claim 1, wherein the new vector of length N includes entries corresponding to the tones other than the tones allocated to the transmitter.

6. The method of claim 1, wherein the discrete signal of complex data symbols includes odd numbered symbols and even numbered symbols, the mapping further comprising phase rotating each even numbered symbol by π/4.

7. An apparatus for wireless communication, comprising: a symbol-to-symbol mapping circuit for mapping received data symbols to a discrete signal of complex data symbols; a digital signal processor (DSP), the DSP configured for: receiving the discrete signal of complex data symbols; calculating, based on the received discrete signal of complex data symbols, frequency responses corresponding to tones through an M-point DFT; outputting a frequency response vector of length M; expanding the frequency response vector of length M to a new vector of length N by inserting zero-value symbols corresponding to tones other than tones allocated to a transmitter; performing an N-point inverse discrete Fourier transform operation on the new vector of length N after zero insertion to obtain a vector of digital signal samples; and pre-pending a cyclic prefix to the vector of digital signal samples.

8. The apparatus of claim 7, further comprising a digital-to-analog converter for generating an OFDM-based signal using the vector of digital signal samples with the pre-pended cyclic prefix for transmission.

9. The apparatus of claim 7, wherein the tones are allocated tones available for transmission.

10. The apparatus of claim 7, wherein the allocated tones are contiguous, and wherein the new vector of length N is larger than the frequency response vector of length M, and wherein the allocated tones are equally spaced.

11. The apparatus of claim 7, wherein the new vector of length N includes entries corresponding to the tones other than the tones allocated to the transmitter.

12. The apparatus of claim 7, wherein the discrete signal of complex data symbols includes odd numbered symbols and even numbered symbols, and wherein symbol-to-symbol mapping circuit is further adapted for phase rotating each even numbered symbol by π/4.

13. An apparatus for wireless communication, in a system comprising: a processor; memory in electronic communication with the processor; and instructions stored in the memory and operable, when executed by the processor, to cause the apparatus to: map received data symbols to a discrete signal of complex data symbols; receive the discrete signal of complex data symbols; calculate, based on the received discrete signal of complex data symbols, frequency responses corresponding to tones through an M-point DFT; output a frequency response vector of length M; expand the frequency response vector to a new vector of length N by inserting zero-value symbols corresponding to tones other than tones allocated to a transmitter; perform an N-point inverse discrete Fourier transform operation on the new vector of length M after zero insertion to obtain the vector of digital samples; and pre-pend a cyclic prefix to the vector of the digital signal samples.

14. The apparatus of claim 13, wherein the instructions are further executable by the processor to generate an OFDM-based signal using the vector of digital signal samples with the pre-pended cyclic prefix for transmission.

15. The apparatus of claim 13, wherein the tones are allocated tones available for transmission.

16. The apparatus of claim 13, wherein the allocated tones are contiguous, and wherein the new vector of length N is larger than the frequency response vector of length M, and wherein the allocated tones are equally spaced.

17. The apparatus of claim 13, wherein the new vector of length N includes entries corresponding to the tones other than the tones allocated to the transmitter.

18. The apparatus of claim 13, wherein the discrete signal of complex data symbols includes odd numbered symbols and even numbered symbols, and wherein the instructions are further executable by the processor to adapt for phase rotating each even numbered symbol by π/4.

Description

DESCRIPTION OF DRAWINGS

(1) FIG. 1 is a block diagram of an OFDM system.

(2) FIG. 2A is a block diagram of an interpolation system used by the OFDM system of FIG. 1.

(3) FIG. 2B is a block diagram of another interpolation system used by the OFDM system of FIG. 1.

(4) FIG. 3A is a graph showing symbols mapped to prescribed time instants in the time domain according to the OFDM technique implemented by the system of FIG. 1.

(5) FIG. 3B is a graph showing the frequency domain response of the graph of FIG. 3B.

(6) FIG. 4A shows an implementation technique for producing a digital signal sample vector using time domain symbol mapping in the case where the allocated tones are contiguous.

(7) FIG. 4B is a block diagram showing a communication system for producing a digital signal sample vector in the case where the allocated frequency tones are contiguous.

(8) FIG. 4C is a graph showing the mapping of the symbols to the prescribed time instants, the expansion of the mapped symbols, and the use of a sinc function to interpolate the expanded symbols.

(9) FIG. 4D is a graph showing the large peak-to-average ratio of the resulting digital signal sample vector when the symbols are mapped in the frequency domain in the prior OFDM systems.

(10) FIG. 4E is a graph showing the reduced peak-to-average ratio of the resulting digital signal sample vector when the symbols are mapped in the time domain using the technique of FIGS. 4A-4C.

(11) FIG. 5A shows another implementation technique for producing the digital signal sample vector using time domain symbol mapping in the case where the allocated tones are equally spaced in frequency.

(12) FIG. 5B is a block diagram showing a communication system for producing a digital signal sample vector in the case where the allocated frequency tones are equally spaced.

(13) FIG. 5C is a graph showing the mapping of the symbols to the prescribed time instants, the expansion of the mapped symbols, and the use of a sinc function to interpolate the symbols.

(14) FIG. 5D is a graph showing the reduced peak-to-average ratio of the resulting digital signal sample vector when the symbols are mapped in the time domain using the technique of FIGS. 5A-5C.

(15) FIG. 6 is a graph showing π/4 symbol rotation.

(16) FIG. 7 shows the use of a cyclic shift of the real and imaginary signal components.

(17) FIG. 8A is a graph showing application of a windowing function in the frequency domain to further reduce the peak-to-average ratio.

(18) FIG. 8B is a block diagram showing a technique using more tones than the number of symbols to be transmitted for producing a digital signal sample vector.

(19) FIG. 8C is a graph showing the use of an interpolation function corresponding to the window function of FIG. 8B to the symbols mapped to the prescribed time instants.

(20) FIG. 8D is a graph showing the reduced peak-to-average ratio of the resulting digital signal sample vector when the symbols are mapped in the time domain using the technique of FIGS. 8A-8C.

(21) Like reference symbols in the various drawings indicate like elements.

DETAILED DESCRIPTION

(22) Referring to FIG. 1, an orthogonal frequency division multiplexing (OFDM) communication system 10 is shown. OFDM communication system 10 receives a first constellation of symbols {B.sub.i} 12 and provides the symbols to a symbol-to-symbol mapping circuit 14, that produces a second constellation of complex symbols {C.sub.i} 16. The complex symbols 16 represent data or a stream of data to be transmitted by the OFDM communication system, and may be chosen from a variety of symbol constellations including, but not limited to phase shift keying (PSK) and quadrature amplitude modulation (QAM) symbol constellations. The symbol-to-symbol mapping performed by the mapping circuit 14 is an optional step performed by the OFDM communication system 10.

(23) Next, a time instant mapping circuit 18 maps each complex symbol 16 to a prescribed time instant within a given OFDM symbol duration. The mapping operation is performed in the time domain such that the mapping circuit 18 generates a discrete signal of mapped symbols within the time domain symbol duration. The output of the mapping circuit 18 is provided to an interpolation circuit 20, that produces a series of digital signal samples {S.sub.i} 22. The digital signal samples 22 are formed by sampling a continuous signal, which is constructed by applying one or more predetermined continuous interpolation functions to the mapped complex symbols 19. Alternatively, the digital signal samples 22 are formed by directly applying one or more predetermined discrete interpolation functions to the mapped complex symbols 19. When using the technique of applying discrete interpolation functions, no intermediate continuous signal is generated and the step of sampling the continuous signal is not necessary. The operation of the interpolation circuit 20 is described in greater detail below. A cyclic prefix circuit 24 receives the series of digital signal samples 22 from the interpolation circuit 20 and prepends a cyclic prefix to the digital signal samples 22. The cyclic prefix circuit 24 operates to copy and prepend the last portion of the digital signal sample vector S 22 to the beginning of the OFDM symbol duration. The resulting digital signal samples 22 with the prepended cyclic prefix are converted to an analog signal by a digital to analog converter 28. The resulting analog signal is further processed by a pulse shaping filter 30, the output of which is modulated to a carrier frequency, and amplified by a power amplifier unit 32 for transmission through an antenna 34.

(24) In one implementation of the OFDM communication system 10, the symbol-to-symbol mapping circuit 14, the time instant mapping circuit 18, the interpolation circuit 20, and the cyclic prefix circuit 24 are implemented in a digital signal processor (DSP) 26, and may include a combination of hardware modules and/or software modules. These circuits 14, 18, 20, and 24 can also be implemented as separate discrete circuits within the OFDM communication system 10.

(25) The details of the interpolation circuit 20 are shown in FIG. 2A. The interpolation circuit 20 includes an interpolation function module 21 that applies one or more continuous interpolation functions to the discrete signal of mapped symbols 19 to generate a continuous signal in which signal variation between adjacent symbols is minimized. Thus, the continuous signal has a low peak-to-average ratio. The interpolation functions may be precomputed and stored in an interpolation function memory 23 connected to the interpolation function module 21. A frequency tone and time instant allocation circuit 27 is connected to the interpolation function memory 23 and defines an allocated tone set selected from frequency tones distributed over a predetermined bandwidth associated with the OFDM communication system 10. The allocated tone set is then provided to the interpolation function memory 23. The frequency tone and time instant allocation circuit 27 also defines the prescribed time instants distributed over the time domain symbol duration, which can also be stored in the interpolation function memory 23 for use by the interpolation function module 21 as well as other modules within the DSP 26. The interpolation circuit 20 also includes a sampling circuit 25 for receiving and sampling the continuous signal at discrete time instants distributed over the time domain symbol duration to generate the vector of digital signal samples 22. Alternatively, in FIG. 2B the interpolation function module 21 applies one or more discrete interpolation functions to the discrete signal of mapped symbols 19 to directly generate the digital signal sample vector 22, in which case the sampling circuit 25 (of FIG. 2A) is not needed. Through applying the discrete interpolation functions, the interpolation function module 21 effectively combines the processing steps of applying the continuous interpolation functions and sampling the intermediate continuous signal.

(26) FIG. 3A graphically depicts the signal processing steps performed by the various circuits of the DSP 26. More specifically, FIG. 3A shows the construction of the signal to be transmitted in a given OFDM time domain symbol duration 40. The time domain symbol duration 40 is a time interval from 0 to T. For purposes of the following description, the OFDM symbol duration T does not include the cyclic prefix. The signal to be transmitted in the symbol duration 40 is represented by complex symbols C.sub.1, C.sub.2, C.sub.3, . . . , C.sub.M 16 that are mapped to the prescribed time instants, where M denotes the number of symbols to be transmitted in the symbol duration 40.

(27) In one implementation, the OFDM communication system 10 is a multiple access communication system where the entire bandwidth available to all transmitters within the system is divided into F orthogonal frequency tones, f.sub.1, f.sub.2, . . . , f.sub.F. In the given symbol duration 40, a particular transmitter operating within a multiple access communication system is allocated M frequency tones f.sub.i(1), f.sub.i(2), . . . , f.sub.i(M), which is a subset of f.sub.1, f.sub.2, . . . , f.sub.F, (the total number of frequency tones) in order to transmit the signal. As part of this implementation, the number of tones allocated to a particular transmitter is equal to the number of symbols to be transmitted by that transmitter. Later in FIG. 8A, the number of allocated tones can be greater than the number of symbols to be transmitted. The remaining frequency tones can be used by other transmitters within the communication system. This technique allows OFDM communication system 10 to operate as a multiple access communication system.

(28) The complex data symbols C.sub.1, C.sub.2, C.sub.3, . . . , C.sub.M 16 are first mapped to t.sub.1, t.sub.2, t.sub.3, . . . , t.sub.M, respectively, where t.sub.1, t.sub.2, t.sub.3, . . . , t.sub.M are M prescribed time instants within the time domain symbol duration 40. The mapping operation generates a discrete signal of mapped symbols. It should be noted that the number of prescribed time instants is equal to the number of symbols M to be transmitted. As described above, the symbol mapping occurs in the time domain. Continuous interpolation functions 42 are then applied to the discrete signal of mapped symbols 16 to generate a continuous function CF(t) for tin the time interval from 0 to T.

(29) The interpolation functions 42 are constructed such that the values of the continuous function CF(t) at time instants t.sub.1, t.sub.2, t.sub.3, . . . , t.sub.M are respectively equal to C.sub.1, C.sub.2, C.sub.3, . . . , C.sub.M and the frequency response of the continuous function CF(t) contains only sinusoids at the allocated tones. Therefore, CF(t) is constructed as

(30) $\mathrm{CF}\left(t\right)=\underset{k=1}{\overset{M}{.Math.}}{A}_k{e}^{J2\pi f_{i\left(k\right)}t}$
where J=√{square root over (−1)} and coefficients A.sub.k are given by

(31) $\left[\begin{array}{c}{A}_1\\ \mathrm{.Math.}\\ A_M\end{array}\right]={\left[\begin{array}{ccc}{e}^{J2\pi f_{i\left(1\right)}{t}_1}& \mathrm{.Math.}& e^{J2\pi f_{i\left(M\right)}{t}_1}\\ \mathrm{.Math.}& & \mathrm{.Math.}\\ e^{J2\pi f_{i\left(1\right)}{t}_M}& \mathrm{.Math.}& e^{J2\pi f_{i\left(M\right)}{t}_M}\end{array}\right]}^{-1}\left[\begin{array}{c}{C}_1\\ \mathrm{.Math.}\\ C_M\end{array}\right]$
Thus, each coefficient A.sub.k is generated by multiplying a matrix of predetermined sinusoids with the single column of data symbols C.sub.1, C.sub.2, C.sub.3, . . . , C.sub.M 16.

(32) FIG. 3B shows the frequency response of the continuous function CF(t).

(33) More specifically, FIG. 3B shows that the frequency response of the continuous function is non-zero only at the allocated frequency tones f.sub.i(1), f.sub.i(2), . . . , f.sub.i(M), and is zero at all other frequency tones.

(34) The output of the DSP 26 is a vector of digital signal samples S 22, which are the samples of the continuous function CF(t) at discrete time instants 0, T/N, 2T/N, . . . , T(N−1)/N, that is, S.sub.1=CF(0), S.sub.2=CF(T/N), S.sub.3=CF(2T/N), . . . , S.sub.N=CF(T(N−1)/N), where N is the number of discrete time instants in the vector of digital signal samples 22. In a general form, t.sub.1, . . . , t.sub.M may not necessarily be equal to any of the time instants 0, T/N, 2T/N . . . , T(N−1)/N. Therefore, while the digital signal samples S 22 may occur at the time instants t.sub.1, . . . , t.sub.M, the OFDM communication system 10 does not require that the time instants 0, T/N, 2T/N . . . , T(N−1)/N be equal to t.sub.1, . . . , t.sub.M.

(35) In another implementation of OFDM communication system 10, the digital signal samples S 22 may be generated by the DSP 26 by directly multiplying a matrix of precomputed sinusoidal waveforms Z, operating as discrete interpolation functions, with the discrete signal of mapped symbols C in order to satisfy the transformation function S=ZC according to the following:

(36) $S=\left[\begin{array}{c}{S}_1\\ \mathrm{.Math.}\\ S_N\end{array}\right]=\left[\begin{array}{ccc}{e}^{J2\pi f_{i\left(1\right)}0}& \mathrm{.Math.}& e^{J2\pi f_{i\left(M\right)}0}\\ \mathrm{.Math.}& & \mathrm{.Math.}\\ e^{J2\pi f_{i\left(1\right)}T\frac{N-1}{N}}& \mathrm{.Math.}& e^{J2\pi f_{i\left(M\right)}T\frac{N-1}{N}}\end{array}\right]\left[\begin{array}{c}{A}_1\\ \mathrm{.Math.}\\ A_M\end{array}\right]={\left[\begin{array}{ccc}{e}^{J2\pi f_{i\left(1\right)}0}& \mathrm{.Math.}& e^{J2\pi f_{i\left(M\right)}0}\\ \mathrm{.Math.}& & \mathrm{.Math.}\\ e^{J2\pi f_{i\left(1\right)}T\frac{N-1}{N}}& \mathrm{.Math.}& e^{J2\pi f_{i\left(M\right)}T\frac{N-1}{N}}\end{array}\right]\left[\begin{array}{ccc}{e}^{J2\pi f_{i\left(1\right)}{t}_1}& \mathrm{.Math.}& e^{J2\pi f_{i\left(M\right)}{t}_1}\\ \mathrm{.Math.}& & \mathrm{.Math.}\\ e^{J2\pi f_{i\left(1\right)}{t}_M}& \mathrm{.Math.}& e^{J2\pi f_{i\left(M\right)}{t}_M}\end{array}\right]}^{-1}\hspace{1em}\left[\begin{array}{c}{C}_1\\ \mathrm{.Math.}\\ C_M\end{array}\right]=\mathrm{ZC}$
where C represents the symbol vector, and the matrix Z represents the product of the two matrices in the second line of the above equation. Each column (i) of matrix Z represents the interpolation function 42 of a corresponding symbol C.sub.i to generate the digital signal samples S 22. As such, the matrix Z can be pre-computed and stored in the interpolation function memory 23 of the interpolation circuit 20 (FIG. 2B). The interpolation circuit 20 then applies the discrete interpolation functions 42 defined by the matrix Z to the discrete signal of mapped complex symbols C 16 in order to satisfy the criteria of S=ZC and to generate the vector of digital signal samples 22.

(37) The purpose of constructing the signal in the time domain is to directly map the symbols 16, which have a desirable low peak-to-average ratio property, to the prescribed time instants within the symbol duration 40. Appropriate interpolation functions 42 are selected to obtain the continuous function CF(t) and the digital signal samples 22 such that the desirable low peak-to-average ratio property of the symbols 16 is substantially preserved for the continuous function and for the digital signal samples 22. The peak-to-average ratio property of the resulting (interpolated) continuous function CF(t) and the digital signal samples 22 is dependent upon the interpolation functions 42, the choice of allocated frequency tones f.sub.i(1), f.sub.i(2), . . . , f.sub.i(M) from the set of tones, and the prescribed time instants t.sub.1, . . . , t.sub.M.

(38) Referring to FIG. 4A, one implementation of the OFDM communication system 10 allocates tones f.sub.i(1), f.sub.i(2), . . . , f.sub.i(M) to the transmitter associated with the communication system that are a subset of contiguous tones in the tone set f.sub.1, f.sub.2, . . . , f.sub.F. Therefore, f.sub.i(k)=f.sub.0+(k−1)Δ, for k=1, . . . , M, where M is the number of symbols. If the OFDM communication system 10 is a multiple access system, each transmitter associated with the communication system is allocated a non-overlapping subset of frequency tones. For purposes of description, let f.sub.0=0. The construction for the other cases where f.sub.0≠0 can be similarly obtained.

(39) Complex symbols C.sub.1, . . . , C.sub.M 16 are mapped in the time domain to the following time instants t.sub.k=(k−1)T/M, for k=1, . . . , M. As part of this implementation, the prescribed time instants t.sub.1, . . . , t.sub.M are equally-spaced time instants uniformly distributed over the entire OFDM symbol duration 40 as shown in the first time domain graph of FIG. 4A. Given the choice of the allocated frequency tones and prescribed time instants, the matrix Z, which is used to generate the digital signal samples S as discussed in FIGS. 3A-3B, can be simplified to

(40) $Z=\frac{1}{M}\left[\begin{array}{ccc}{e}^{J2\pi f_{i\left(1\right)}0}& \mathrm{.Math.}& e^{J2\pi f_{i\left(M\right)}0}\\ \mathrm{.Math.}& & \mathrm{.Math.}\\ e^{J2\pi f_{i\left(1\right)}T\frac{N-1}{N}}& \mathrm{.Math.}& e^{J2\pi f_{i\left(M\right)}T\frac{N-1}{N}}\end{array}\right]\left[\begin{array}{ccc}{e}^{-J2\pi f_{i\left(1\right)}{t}_1}& \mathrm{.Math.}& e^{-J2\pi f_{i\left(1\right)}{t}_M}\\ \mathrm{.Math.}& & \mathrm{.Math.}\\ e^{-J2\pi f_{i\left(M\right)}{t}_1}& \mathrm{.Math.}& e^{-J2\pi f_{i\left(M\right)}{t}_M}\end{array}\right].$

(41) The second time domain graph of FIG. 4A shows the resulting digital signal sample vector S 22 after the interpolation circuit 20 applies the interpolation functions 42 defined by the matrix Z to the complex symbols 16 according to the expression S=ZC. As part of this implementation, the sampling module 25 is not generally used as the digital signal sample vector S 22 is directly generated from the discrete signal of mapped symbols using the transformation function S=ZC.

(42) Turning to FIG. 4B, a digital processing system 50 provides another technique for obtaining the vector of digital signal samples S. A DFT circuit 52 receives a discrete signal of complex data symbols C.sub.i, and calculates the frequency responses A.sub.1, . . . , A.sub.M, at tones f.sub.i(1), f.sub.i(2), . . . , f.sub.i(M), through an M-point discrete Fourier transform (DFT). The vector [A.sub.1, . . . , A.sub.M] 54 output by the DFT circuit 52 is then expanded to a new vector of length N (the total number of time instants in the discrete signal vector S) by zero insertion at block 56. More specifically, this process involves putting the k.sup.th symbol A.sub.k to the i(k).sup.th element of the new vector, for k=1, . . . , M, where f.sub.i(k) is the k.sup.th tone allocated to the transmitter, and inserting zeros in all the remaining elements. Finally, an IDFT circuit 58 performs an N-point inverse discrete Fourier transform on the resulting vector (after zero insertion) to obtain the digital signal sample vector S. The collective procedure of DFT, zero insertion and IDFT is one way of implementing the discrete interpolation functions.

(43) Turning to FIG. 4C, another technique for obtaining the digital signal samples S is shown. For simplicity of description, it is assumed that the allocated contiguous tones f.sub.i(1), f.sub.i(2), . . . , f.sub.i(M) are centered at frequency 0. The construction for the other cases where the allocated tones are not centered at frequency 0 can be similarly obtained. As with FIG. 4A, the prescribed time instants t.sub.1, . . . , t.sub.M are equally-spaced time instants uniformly distributed over the entire OFDM symbol duration 40.

(44) The complex symbols C.sub.1, . . . , C.sub.M are first mapped in the time domain to time instants t.sub.1, . . . t.sub.M respectively. Next, the mapped symbols C.sub.1, . . . , C.sub.M are leftwards and rightwards shifted and replicated to an expanded set of prescribed time instants, which is a superset of t.sub.1, . . . t.sub.M and consists of an infinite number of equally-spaced time instants covering the time interval from −∞ to +∞. This technique creates an infinite series of mapped symbols C. The continuous function CF(t) is then constructed by interpolating the infinite series of mapped symbols using a sinc interpolation function 60. Mathematically, the above steps construct the continuous function CF(t) as

(45) $\mathrm{CF}\left(t\right)=\underset{i=1}{\overset{M}{.Math.}}\left\{C_i\underset{k=-\infty }{\overset{\infty }{.Math.}}\sin c\left(t-t_i-\mathrm{kT},\frac{T}{M}\right)\right\}.$
where sinc(a,b)=sin(πa/b)/(πa/b). The sine interpolation function 60 can also be precomputed and stored in the interpolation function memory 23. As discussed in FIG. 3A the digital signal samples S 22 are the samples of the continuous function CF(t) at time instants 0, T/N, . . . , T(N−1)/N. In FIGS. 4A-4C, if N is a multiple of M, then S.sub.1+(k−1)N/M=C.sub.k, for k=1, . . . , M. It should be noted that the continuous function CF(t) only applies to the symbol duration 40 from 0 to T. The use of time interval from −∞ to +∞ is solely for the purpose of mathematically constructing CF(t). The discrete interpolation functions, which combine the continuous interpolation functions and the sampling function, can be derived easily from the above description.

(46) For comparison purposes, FIG. 4D illustrates the resulting peak-to-average ratio for a digital signal sample vector S 62 and its associated transmitted OFDM signal 64 produced by symbols 16 where the signal is constructed in the frequency domain. As described above, this known technique of mapping the symbols 16 in the frequency domain produces a large signal variation in the transmitted OFDM signal 64 and results in a large peak-to-average ratio.

(47) FIG. 4E illustrates the resulting small signal variation and low peak-to-average ratio of the digital signal sample vector S 66 associated with the transmitted OFDM signal 68. As will be appreciated by comparing FIGS. 4D and 4E, mapping the constellation of complex symbols 16 in the time domain produces an OFDM signal 68 having a significantly reduced peak-to-average ratio.

(48) FIG. 5A shows a second implementation of the OFDM communication system 10, and serves to further generalize the system shown in FIGS. 4A-4C. As part of OFDM system 10, tones, f.sub.i(1), f.sub.i(2), . . . , f.sub.i(M), allocated to the transmitter associated with the communication system, are a subset of equally-spaced tones in the tone set f.sub.1, f.sub.2, . . . , f.sub.F. Therefore, f.sub.i(k)=f.sub.0+(k−1)LΔ, for k=1, . . . , M, and L is a positive integer number representing the spacing between two adjacent allocated frequency tones. When L=1, this implementation is equivalent to the implementation technique described in FIGS. 4A-4C. For the sake of description, let f.sub.0=0. The construction for the other cases where f.sub.0≠0 can be similarly obtained.

(49) In this case where the allocated tones are equally-spaced tones, the constructed continuous function CF(t) is identical in each of the L time intervals, [0,T/L), [T/L,2T/L), . . . , and [(L−1)T/L, T/L). As part of this technique, symbols C.sub.1, . . . , C.sub.M 16 are mapped to the following time instants t.sub.k=(k−1)T/M/L, for k=1, . . . , M. In this implementation, the prescribed time instants t.sub.1, . . . , t.sub.M are equally-spaced time instants uniformly distributed over a fraction (1/L) of the symbol duration 70. As a comparison, in the case of allocated contiguous tones (FIG. 4A), the prescribed time instants are equally-spaced and distributed over the entire symbol duration, as discussed with respect to FIG. 4A.

(50) The procedure for obtaining the digital signal samples S 22 described in FIG. 4A can also be applied with respect to FIG. 5A. More specifically, the digital signal sample vector S is the product of matrix Z (defining the discrete interpolation functions) and the symbol vector C. Given the choice of the allocated frequency tones and prescribed time instants, the matrix Z, which is used to generate the digital signal samples 22 from the discrete signal of mapped symbols, can be simplified to the same formula as in FIG. 4A with the only change in the definition off f.sub.i(1), f.sub.i(2), . . . , f.sub.i(M) and t.sub.1, . . . , t.sub.M.

(51) In FIG. 5B, the procedure of obtaining the digital signal sample vector S 22 described in FIG. 4B can also be applied to the case of allocated frequency tones that are equally spaced tones. More specifically, a digital processing system 100 provides another technique for obtaining the vector of digital signal samples S. A DFT circuit 102 receives a discrete signal of complex data symbols C.sub.i and calculates the frequency responses A.sub.1, . . . , A.sub.M, at tones f.sub.i(1), f.sub.i(2), . . . , f.sub.i(M), through an M-point discrete Fourier transform (DFT). The vector [A.sub.1, . . . A.sub.M] 104 output by the DFT circuit 102 is then expanded to a new vector of length N (the total number of time instants in the digital signal sample vector S) by zero insertion at block 106. More specifically, this process involves putting the k.sup.th symbol A.sub.k to the i(k)th element of the new vector, for k=1, . . . , M, where f.sub.i(k) is the k.sup.th tone allocated to the transmitter, and inserting zeros in all the remaining elements. Finally, an IDFT circuit 108 performs an N-point inverse discrete Fourier transform on the resulting vector (after zero insertion) to obtain the time domain digital signal sample vector S. The collective procedure of DFT, zero insertion and IDFT is one way of implementing the discrete interpolation functions.

(52) FIG. 5C is the counterpart of FIG. 4C, where symbols C.sub.1, . . . , C.sub.M are first mapped to t.sub.1, . . . , t.sub.M respectively over a fraction (1/L) of the symbol duration 70. The symbol mapping is also performed in the time domain. Next the mapped symbols C.sub.1, . . . , C.sub.M are leftwards and rightwards shifted and replicated to an expanded set of prescribed time instants from −∞ to +∞ which creates an infinite series of symbols. The continuous function CF(t) is then constructed by interpolating the infinite series of mapped symbols with a sinc interpolation function 72. Thus, the continuous function CF(t) includes the digital signal samples mapped to the prescribed time instants as well as digital sample points between the prescribed time instants. Mathematically, the above steps construct the continuous function as

(53) $\mathrm{CF}\left(t\right)=\underset{i=1}{\overset{M}{.Math.}}\left\{C_i\underset{k=-\infty }{\overset{\infty }{.Math.}}\sin c\left(t-t_i-\mathrm{kT}\frac{1}{L},\frac{T}{M}\frac{1}{L}\right)\right\}.$
With continued reference to FIG. 5C, each sinc interpolation function 72 is narrower and therefore decays faster than the sinc interpolation function 60 shown in FIG. 4C. The sinc interpolation function 72 can also be precomputed and stored in the interpolation function memory 23 for use by the interpolation function module 21. The digital sample vector S 22 can be obtained in the same technique shown in FIG. 4C. In FIGS. 5A and 5C, if N is a multiple of ML, then S.sub.1+(k−1)N/M/L+(j−1)N/L=C.sub.k, for k=1, . . . , M, and j=1, . . . , L. The discrete interpolation functions, which combine the continuous interpolation functions and the sampling function, can be derived easily from the above description.

(54) FIG. 5D illustrates the resulting small signal variation and low peak-to-average ratio of the digital signal sample vector S 74 associated with the transmitted OFDM signal 76. As will be appreciated by comparing FIGS. 4D and 5D, mapping the constellation of complex symbols 16 in the time domain produces an OFDM signal 76 having a significantly lower peak-to-average ratio.

(55) Referring now to FIG. 6, a π/4 symbol rotation technique is used to further reduce the peak-to-average ratio of the transmitted OFDM signal. At an OFDM symbol duration, if symbols B.sub.1, . . . , B.sub.M of the constellation are to be transmitted, symbols B.sub.1, . . . , B.sub.M are mapped to another block of complex symbols C.sub.1, . . . , C.sub.M, where each odd number symbol remains unchanged and each even number symbol is phase rotated by π/4. For example, if symbols B.sub.1, . . . , B.sub.M belong to a QPSK constellation {0, π/2, π, π3/2}, the odd number symbols C.sub.k still belong to the same QPSK constellation, while after being phase rotated the even number symbols C.sub.k belong to another QPSK constellation {π/4, π3/4, π5/4, π7/4}. Symbols C.sub.1, . . . , C.sub.M are then used to construct the digital signal samples 22 in the time domain as described above with respect to FIGS. 3A-5C.

(56) With reference to FIG. 7, another technique for reducing the peak-to-average ratio is shown, which introduces a cyclic offset of the real and imaginary signal components. This technique involves a first step of offsetting the imaginary components of the digital signal samples S 22, which have been generated using the technique of FIGS. 3A-5C, by an integer number of samples. If necessary, the technique then involves a second step of adjusting the timing by a fraction of a sample period between the real and the imaginary signal components in the transmit path.

(57) At an OFDM symbol duration, if the digital signal samples S.sub.1, S.sub.2, . . . , S.sub.N have been obtained using the method as described in FIGS. 3A-5C, the digital signal sample vector S is then mapped to another vector S′ as follows. The real component of digital signal sample S′.sub.k is equal to that of digital signal sample S.sub.k. The imaginary component of digital signal sample S′.sub.k is equal to that of digital signal sample S.sub.j where index j=(k+d−1)mod N+1, for k=1, . . . , N, with mod representing a module operation. The parameter d is an integer representing the cyclic offset, in terms of number of samples, between the real and imaginary components.

(58) In one implementation, the value of d is determined by

(59) $\frac{N}{2\mathrm{LM}},$
where L is discussed in FIG. 5A. In one aspect of this technique, d is chosen to be close to

(60) $\frac{N}{2\mathrm{LM}}.$
For example, d can be the integer closest to

(61) $\frac{N}{2\mathrm{LM}},$
the largest integer not greater than

(62) 0$\frac{N}{2\mathrm{LM}},$
or the smallest integer not smaller than

(63) $\frac{N}{2\mathrm{LM}}.$
In one example, d is chosen to be the largest integer not greater than

(64) $\frac{N}{2\mathrm{LM}}.$
This example can be easily extended for other choices of d.

(65) The digital signal sample vector S′ is then passed to the cyclic prefix prepender circuit 24, as shown in FIG. 1. Therefore, the operation of half symbol cyclic shifting is carried out before the operation of prepending the cyclic prefix, such as that performed by the cyclic prefix circuit 24 of FIG. 1.

(66) Not specifically shown in FIG. 7, when or after the sample vector S′ and the cyclic prefix are outputted to the digital to analog converter 28, the imaginary components are further delayed by an amount of

(67) $\left(\frac{N}{2\mathrm{LM}}-d\right)\frac{T}{N},$
which is a fraction of a sample period T/N.

(68) As a variation of the technique shown in FIG. 7 (not specifically shown), another technique for achieving a similar result can be used to eliminate the second step of adjusting timing by a fraction of a sample period between the real and the imaginary signal components in the transmit path. As part of this technique, the real and the imaginary components of the desired digital signal samples S 22 are generated separately as described by the following.

(69) A first series of digital signal samples 22 are generated using the technique of FIGS. 3A-5C. The real components of the desired digital signal samples 22 are equal to those of the first series of samples. A second series of digital signal samples 22 are generated using the technique of FIGS. 3A-5C except for the following changes. The imaginary components of the desired digital signal samples are equal to those of the second series of samples. In the general method described in FIGS. 3, 4A, and 5A, the matrix

(70) $\hspace{1em}\left[\begin{array}{ccc}{e}^{J2\pi f_{i\left(1\right)}0}& \mathrm{.Math.}& e^{J2\pi f_{i\left(M\right)}0}\\ \mathrm{.Math.}& & \mathrm{.Math.}\\ e^{J2\pi f_{i\left(1\right)}T\frac{N-1}{N}}& \mathrm{.Math.}& e^{J2\pi f_{i\left(M\right)}T\frac{N-1}{N}}\end{array}\right]$
is changed to

(71) $\hspace{1em}\left[\begin{array}{ccc}{e}^{J2\pi f_{i\left(1\right)}\left(0-\frac{T}{2\mathrm{LM}}\right)}& \mathrm{.Math.}& e^{J2\pi f_{i\left(M\right)}\left(0-\frac{T}{2\mathrm{LM}}\right)}\\ \mathrm{.Math.}& & \mathrm{.Math.}\\ e^{J2\pi f_{i\left(1\right)}\left(T\frac{N-1}{N}-\frac{T}{2\mathrm{LM}}\right)}& \mathrm{.Math.}& e^{J2\pi f_{i\left(M\right)}\left(T\frac{N-1}{N}-\frac{T}{2\mathrm{LM}}\right)}\end{array}\right].$

(72) In the block diagram method described with regard to FIG. 4B, an additional operation is required after zero insertion (block 56) and before N-point IDFT (block 58), where each element k in the expanded length N vector is phase rotated by

(73) $e^{-J2\pi f_k\frac{T}{2\mathrm{LM}}}.$

(74) Referring to FIGS. 8A-8D, another technique for further reducing the peak-to-average ratio is implemented by allocating more frequency tones than the number of complex symbols to be transmitted in a symbol duration 40. In FIGS. 3-7, the number of tones allocated to the transmitter associated with the communication system is equal to the number of symbols to be transmitted in a given OFDM symbol duration. Compared with the other techniques described with respect to the previous figures, the technique of FIGS. 8A-8D requires additional overhead of bandwidth to transmit the same number of complex symbols.

(75) For example, if the communication system 10 is allocated M+M.sub.ex contiguous frequency tones, f.sub.i(1), f.sub.i(2), . . . , f.sub.i(M+Mex), and M symbols C.sub.1, . . . , C.sub.M of the constellation are to be transmitted at an OFDM symbol duration, from the comparison of FIGS. 4A and 5A, the case of allocated contiguous tones can be easily extended to the case of allocated equally-spaced tones. As part of this implementation of the OFDM communication system 10, M.sub.ex is a positive number representing the number of excess tones to be used and is assumed to be an even number. Therefore, the allocated tone

(76) $f_{i\left(k\right)}=f_0+\left(k-\frac{M_{\mathrm{ex}}}{2}-1\right)\Delta ,$
for k=1, . . . M+M.sub.ex. For purposes of description, let f.sub.0=0. The construction for the other cases where f.sub.0≠0 can be similarly obtained.

(77) As with the technique described with respect to FIG. 4A, the prescribed time instants are t.sub.k=(k−1)T/M, for k=1, . . . , M, that is, the prescribed time instants t.sub.1, . . . , t.sub.M are equally-spaced time instants in the symbol duration 40.

(78) As part of this technique shown in FIG. 8A, P(f) is a smooth windowing function 90 in the frequency domain, which is non-zero only over interval [f.sub.i(1), f.sub.i(M+Mex)]. In addition, P(f) 90 also satisfies the Nyquist zero intersymbol interference criterion, i.e.,

(79) $\underset{k=-\infty }{\overset{\infty }{.Math.}}P\left(f-\mathrm{kM}\Delta \right)=1$
for any frequency f, where □ is the spacing between adjacent tones.

(80) FIG. 8B shows the block diagram of the technique. As described above, a symbol-to-symbol mapping is optionally performed to generate a discrete signal of mapped complex symbols C.sub.1, . . . , C.sub.M, 16. The frequency responses A.sub.1, . . . , A.sub.M 84 are calculated through an M-point discrete Fourier transform (DFT) of the complex symbols 16 at block 82. At block 86, vector [A.sub.1, . . . , A.sub.M] 84 is cyclically expanded to a new vector A′ of length N and windowed with a windowing function 90 as follows:
A′.sub.kA.sub.g(k)*P((k−1)□+f.sub.1)
where index g(k)=mod(k−i(1)−M.sub.ex/2,M)+1, for k=1, . . . , N.

(81) At block 88, the digital signal sample vector S is obtained by taking an N-point inverse discrete Fourier transform (IDFT) of the new vector A′. Finally, the cyclic prefix is added by cyclic prefix circuit 24 as described above with regard to FIG. 1.

(82) To provide additional insight to the above signal construction technique, assume that the allocated tones f.sub.i(1), f.sub.i(2), . . . , f.sub.i(M+Mex) are centered at frequency 0. In FIG. 8C (as with FIG. 4C), symbols C.sub.1, . . . , C.sub.M are first mapped to equally-spaced time instants in the symbol duration 40, and are then leftwards and rightwards shifted and replicated from −∞ to +∞. What is different from FIG. 4C is that a different interpolation function 92, which is determined by the windowing function 90, is used to generate the continuous function,
CF(t)=Σ.sub.i=1.sup.MC.sub.iΣ.sub.k=−∞.sup.∞p(t−t.sub.i−kT)
where p(t) 92 is the time domain response of P(f) 90. As with FIG. 4C, the digital signal samples are obtained by letting t=0, . . . , T(N−1)/N.

(83) In one exemplary aspect of this technique, if a raised cosine windowing function is used, i.e.,

(84) $P\left(f\right)=\{\begin{array}{cc}\frac{T}{M}& if.Math.f.Math.<\left(1-\beta \right)\frac{M}{2T}\\ \frac{T}{2M}\left\{1+\cos \left[\frac{\pi T}{\beta M}\left(.Math.f.Math.-\frac{\left(1-\beta \right)M}{2T}\right)\right]\right\}& \mathrm{if}\left(1-\beta \right)\frac{M}{2T}\le .Math.f.Math.\le \left(1+\beta \right)\frac{M}{2T}\\ 0& \mathrm{if}.Math.f.Math.>\left(1+\beta \right)\frac{M}{2T}\end{array}$
where β=(M.sub.ex+2)/M represents the percentage of excess tone overhead, then, the interpolation function p(t) 92 is given by

(85) 0$p\left(t\right)=\frac{\sin \left(\pi \mathrm{tM}/T\right)}{\pi \mathrm{tM}/T}\frac{\cos \left(\pi \beta \mathrm{tM}/T\right)}{1-4{\beta }^2{t}^2{M}^2/T^2}.$
As β increases, the interpolation function p(t) 92 decays faster, thereby reducing the probability of having large peak at samples between t.sub.i.

(86) FIG. 8D shows the resulting small signal variation and low peak-to-average ratio of the digital signal sample vector S 94 associated with the transmitted OFDM signal 96. As will be appreciated, mapping the constellation symbols 16 in the time domain produces an OFDM signal 96 having a significantly lower peak-to-average signal ratio.

(87) A number of embodiments of the invention have been described. Nevertheless, it will be understood that various modifications may be made without departing from the spirit and scope of the invention. Accordingly, other embodiments are within the scope of the following claims.