FREQUENCY-DOMAIN MIMO PROCESSING METHOD AND SYSTEM
FIELD OF THE INVENTION In general, the present invention relates to the field of multiple input- multiple output ("MIMO") communication systems. More specifically, the present invention relates to increasing the link capacity of broadband wireless single carrier systems operating in a high delay-spread channel.
BACKGROUND OF THE INVENTION
Both single data-stream array processing methods and multiple data stream methods as known in the art are utilized for purposes of increasing the link capacity of narrowband single carrier communication systems and orthogonal frequency division multiplexing ("OFDM") systems. Various methodologies require the transmit array to have knowledge of the channel response between each transmit antenna element and each receive antenna element, and are therefore referred to as "closed-loop" MIMO. Conversely, prior art transmission methods that do not require channel knowledge at the transmit array are referred to as "open-loop" MIMO. Theoretically, the benefits of having perfect channel knowledge versus not having any channel knowledge at the transmit array are obtained mainly at low signal-to-noise ratios and in correlated channels. The higher link capacities obtained when channel knowledge is present, however, translate into significant capacity gains on the system level. A Cyclic-Prefix single carrier ("CP-single carrier") system known in the art is characterized by the insertion of an OFDM-style cyclic prefix into a sequence of time-domain data symbols. It is well known that the cyclic prefix facilitates the application of frequency-domain equalization and receive array processing methodologies at a receiving device having one or more antennas. Many benefits of frequency-domain equalization and cyclic prefixes have been acquired. However, the idea of applying a MIMO methodology to cyclic
prefix-based single carrier systems is absent from the prior art. The present invention addresses this situation.
SUMMARY OF THE INVENTION One form of the invention is a method of operating a transmitting device having one or more transmission branches with each transmission branch receiving a time-domain signal. First, each time-domain signal is transformed into a frequency-domain signal. Second, each frequency-domain signal is weighted to form one or more weighted frequency-domain signals with each weighted frequency-domain signal being a function of each frequency-domain signal.
A second form of the invention is a method of operating a receiving device having one or more reception branches with each reception branch receiving a time-domain waveform and each time-domain waveform being a function of one or more time-domain signals. First, each time-domain waveform is transformed into a frequency-domain sample. Second, each frequency-domain sample is weighted to form one or more weighted frequency-domain samples with each weighted frequency-domain sample being a function of each frequency-domain sample. The foregoing and other features and advantages of the invention will become further apparent from the following detailed description of the presently preferred embodiment, read in conjunction with the accompanying drawings. The detailed description and drawings are merely illustrative of the invention rather than limiting, the scope of the invention being defined by the appended claims and equivalents thereof.
BRIEF DESCRIPTION OF THE DRAWINGS
FIG. 1 illustrates one embodiment of a timing diagram of a cyclic-prefix- based single-carrier system in accordance with the present invention;
FIG. 2 illustrates one embodiment of a closed-loop MIMO transmitting device for a cyclic-prefix single-carrier system in accordance with the present invention;
FIG. 3 illustrates one embodiment of a frequency-domain weighting module for the closed loop MIMO transmitting device of FIG. 1 ;
FIG. 4 illustrates one embodiment of a linear MIMO receiving device for a cyclic-prefix single-carrier system in accordance with the present invention; and
FIG. 5 illustrates one embodiment of a frequency-domain weighting module for the linear MIMO receiving device of FIG. 4.
DETAILED DESCRIPTION OF THE
PRESENTLY PREFERRED EMBODIMENTS
The present invention is a method for incorporating multiple-stream array processing techniques in single carrier systems comprising a transmitting device having Mr transmit antennas (e.g., a transmitting device 20 illustrated in FIG. 2 where Mτ= 3) and a receiving device having MR receive antennas (e.g., a receiving device 30 illustrated in FIG. 4 where MR = 3). Such methods can employ cyclic redundancy (i.e., cyclic prefixes, cyclic postfixes, or a mixture thereof).
FIG. 1 illustrates a time line wherein the time domain is divided up into a baud interval 11a, a baud interval 11b, and others baud intervals not shown. Each illustrated baud interval consists of a cyclic redundancy of Lcp symbols, such as a cyclic prefix 12a of Lc symbolsi and a cyclic prefix 12b of Lcp symbol. In other embodiments of baud intervals, other forms of cyclic redundancy can be used, such as, for example, cyclic postfix and a mixture of a cyclic prefix and a cyclic postfix. Each baud interval further consists of a data interval of N symbols (Λ/>1), such as a data interval 13a of N symbols-i and a data interval 13b of N symbol. The Lcp symbols of each cyclic prefix
are simply a repetition of the last Lc symbols of a corresponding data interval. Within a given baud interval, the goal of a transmitting device of the present invention (e.g., transmitting device 20) is to deliver Ns sequences of N symbols (e.g., Quadrature Amplitude Modulation ("QAM") symbols or Phase- Shift Keying ("PSK") symbols) within each baud interval. The mth time- domain symbol of they* stream is denoted Sj[m], fory'=1... Ns, and m=0...Λ/-1. The symbols of the Ns streams can belong to a traditional QAM or PSK constellation, time-domain chips in a CDMA system, or other weighted time- domain signals as known in the art. FIG. 2 illustrates transmitting device 20 operating in accordance with the present invention. Transmitting device 20 has three branches for processing a time-domain signal s^m], a time-domain signal S2[m], and a time-domain signal S3[m] (i.e., Mτ~ 3 and Λ/s=3). A following operational description on of transmitting device 20 involves baud interval 11a (FIG.1). As a result, in the following operational description of transmitting device 20, the dependencies of the mathematical quantities on a particular baud interval are omitted for clarity. Furthermore, from the following operational description of transmitting device 20, those having ordinary skill in the art will appreciate alternative embodiments of a transmitting device in accordance with the present invention for processing any number of streams (i.e., 1< v" s<°°) with any number of transmit antennas Mτ.
Transmitting device 20 includes a Time-Frequency Transformer ("TFT") module 21a, a TFT module 21b, and a TFT module 21c. TFT module 21a, TFT module 21 , and TFT module 21c transform time-domain signal stfm], time-domain signal s∑fm], and time-domain signal S
3[m], respectively, into frequency-domain signals of a continuous nature or a discrete nature. In one embodiment, TFT module 21a, TFT module 21 b, and TFT module 21c are Fast Fourier Transform ("FFT") based modules for implementing the following equation [1] in transforming time-domain signal sι[m], time-domain signal S
2[m], and time-domain signal s
3[m], respectively, into a frequency-domain signal Sι[k], a frequency-domain signal S
2[k], and a frequency-domain signal S
3/ΑJ, respectively:
[1]
where Sj[k] is the frequency-domain symbol on the /c"1 frequency bin (or subcarrier) for the h stream, where / =0...Λ/-1 , and N is the FFT size.
Once in the frequency-domain, the transmit array processing required by MIMO is carried out on each frequency bin by a frequency-domain weighting module 22. Specifically, frequency-domain weighting module 22 provides a weighted frequency-domain signal Xι[k], a weighted frequency- domain signal X∑[k], and a weighted frequency-domain signal X3[k], each being a function of frequency-domain signal Sι[k], frequency-domain signal S∑[k], and frequency-domain signal S3[k]. In one embodiment, as illustrated in FIG. 3, frequency-domain weighting module 22 receives frequency-domain signal Sι[k], frequency-domain signal S∑fk], and frequency-domain signal S3[k] on the Ac"1 subcarrier in an Λ/S χ1 vector according to the following equation [2]:
[2]
where =0.../V-1 , and where Λ/s=3.
Accordingly, TFT module 21a implements a serial-to-parallel conversion of the time-domain signal sι[m] (m=0...Λ ) to a frequency- domain signal Sι[0], a frequency-domain signal Sι[1], and a frequency- domain signal Sι[2]. TFT module 21b implements a serial-to-parallel conversion of the time-domain signal s2[m] (m=0...Λ/-1) to a frequency- domain signal S∑β], a frequency-domain signal S∑ft], and a frequency-
domain signal S2 2 . TFT module 21c implements a serial-to-parallel conversion of the time-domain signal S3[m] (/77=0.../V-1 ) to a frequency- domain signal S3[0], a frequency-domain signal S3[1], and a frequency- domain signal S3[2]. In this FFT embodiment, the FFT size N (which in most embodiments is equal to the number of symbols in a data interval 13a) is shown to be 3, although, as will be appreciated by those skilled in the art, N can take on other values such as 512 or 1024.
In response thereto, for each subcarrier k of the N subcarriers, frequency-domain weighting module 22 provides Mr weighted frequency- domain signals Xj[k], y'=1... Mr, which are collected into the /Wrχ1 vector X[k] in accordance with the following equation [3]:
[3]
As shown in FIG. 3, the number of antennas Mγ is shown to be 3, although it will be appreciated by those skilled in the art that Mr can take on other values that are, for example, not necessarily equal to the number of streams N
s. In most embodiments, Mγ is chosen to be greater than or equal to N
s. A Frequency-Time Transformer ("FTT") module 23a receives a weighted frequency-domain signal Xι[k] for k=0...Λ/-1 (i.e., Xι[0], X ] and Xι[2]). An FTT module 23b receives a weighted frequency-domain signal X
2[k] for / =0...Λ (i.e., X
2[0], X
2[1] and X
2[2]). An FTT module 23c receives a weighted frequency-domain signal X
3[k] for =0...Λ/-1 (i.e., 3/O/, X
3[1] and X
3[2J). In one embodiment, frequency-domain weighting module 22 includes a plurality of weight modules 22a-22i and a plurality of sum devices 22j-22r for providing the vector X[k] by a weighting of the frequency-domain signals Sι[0]-S
3[2] in accordance with the following equation [4]:
[4]
where the Mr χ1 vector X[/] is as in equation [3]; where the Λ/S χ1 vector S[/] is as in equation [2]; where the Mrχ1 transmit weight vector for the £ih stream is given by \e [k] in accordance with the following equation [5]:
[5]
and the Ns transmit vectors are arranged into a transmit matrix V[/] in accordance with the following equation [6]:
[6] where, as illustrated in FIG. 3, a plurality of signal paths SP1-SP27 are arranged between weight units 22a-22i and sum units 22j-22r in accordance with equations [4], [5], and [6].
FTT module 23a converts weighted frequency-domain signals Xι[0], Xι[1], and Xι[2J into a weighted time-domain signal xι[m], for m=0...Λ/-1. FTT module 23b converts weighted frequency-domain signals X∑fO], X2II], and X2/27 into weighted time-domain signal x∑ m], for m=0.../V-1. FTT module 23c converts weighted frequency-domain signals X3/0 * X3P] and X3[2] into a weighted time-domain signal X3[m], for m=0... V-1. In one embodiment, FTT module 23a, FTT module 23b, and FTT module 23c are inverse FFT based
modules whereby weighted time-domain signals Xι[m], x∑ϊm], and X3[m] are determined in accordance with the following equation [7]:
[7] At this point, the quantities vJt[m], m=0... Λ/-1 , are defined to be the time-domain equivalent representation of the frequency-domain transmit weights v [k], / =0.../V-1 ,y"=1...Mr, ^ =1...Λ/S> in accordance with the following equation [8]:
[8]
It can then be shown that the weighted time-domain signals xι[m], X2[m], and X3[m] of Equation [7] are given by the following equation [9]:
[9]
As shown in equation [9], the time-domain data is the convolution of the time-domain equivalent transmit weights and the Ns data streams.
A cyclic module 24a adds a cyclic redundancy in the form of a cyclic prefix to weighted time-domain signal xι[m] to form a weighted redundant time-domain signal x-|[i] . A cyclic module 24b adds a cyclic redundancy in the form of a cyclic prefix to weighted time-domain signal x [m] to form a weighted redundant time-domain signal X2[i] . A cyclic module 24c adds a cyclic redundancy in the form of a cyclic prefix to weighted time-domain signal
X3[m] to form a weighted redundant time-domain signal X3U] . In one embodiment, weighted redundant time-domain signals x<|[i] , X2_i] , 3D] are formed in accordance with the following equation [10]:
- fxj[' + N ~ LCp], O ≤ i ≤ LCp - 1
XJLIJ Xj[i - Lcp], Lcp < i ≤ N + Lcp - 1
[10] In an alternative embodiment, cyclic module 24a adds a cyclic redundancy in the form of a cyclic postfix to weighted time-domain signal Xι[m] to form weighted redundant time-domain signal x<|[i] . A cyclic module 24b adds a cyclic redundancy in the form of a cyclic postfix to weighted time- domain signal x∑fm] to form weighted redundant time-domain signal X2D] . A cyclic module 24c adds a cyclic redundancy in the form of a cyclic postfix to weighted time-domain signal X3[m] to form weighted redundant time-domain signal X3 [i] . In one embodiment, weighted redundant time-domain signals xι [i] , X2[i] , and X3[i] having cyclic postfixes are formed in accordance with the following equation [11]:
[11]
In another alternative embodiment, cyclic module 24a adds a cyclic redundancy in the form of a combined cyclic prefix and cyclic postfix to weighted time-domain signal xι[m] to form weighted redundant time-domain signal x-|[i] . A cyclic module 24b adds a cyclic redundancy in the form of a combined cyclic prefix and cyclic postfix to weighted time-domain signal X
2[m] to form weighted redundant time-domain signal X2[i] . A cyclic module 24c adds a cyclic redundancy in the form of a combined cyclic prefix and cyclic
postfix to weighted time-domain signal X
3[m] to form weighted redundant time-domain signal
.
A pulse shaping module 25a shapes weighted redundant time-domain signal x-][i] into a weighted time-domain waveform x-j(t) . A pulse shaping module 25b shapes weighted time-domain signal X2[i] into a weighted time- domain waveform X2(t). A pulse shaping module 25c shapes weighted time- domain signal X3 ] into a weighted time-domain waveform ><3(t) . In one embodiment, pulse shaping module 25a, pulse shaping module 25b, and pulse shaping module 25c provides weighted time-domain waveforms xι(t) , X2(t) . and X3(t), respectively, in accordance with the following equation [12]:
N+Lcp-1
*;( = ∑Zj[n]Psr(t - nTs) n=0
[12]
where psι{t) is the pulse shaping waveform (which typically has the square-root raised cosine spectrum) and Ts is the duration of a symbol.
Weighted time-domain waveforms x-|(t) , X2(t), and X3(t) propagate from antenna 26a, antenna 26b, and antenna 26c, respectively, to an antenna 31a, an antenna 31b, and an antenna 31c, respectively, of receiving device 30 through a multipath channel with an RF impulse response in accordance with the following equation [13]:
[13]
where h and τl are the complex gain for they* transmit antenna and the h receive antenna and time delay of the /*h path, respectively, and where δ{t) is a Dirac delta function. This channel impulse response is assumed to
be time-invariant during the transmitted block of N data symbols, although the extension of this equation to time-varying channels is straightforward.
FIG. 4 illustrates receiving device 30 operating in accordance with the present invention. Receiving device 30 has three branches (i.e., MR= 3) for processing three time-domain waveforms r (t) , r2( )and r3(f)which are a function of a propagation of the weighted time-domain waveforms x-ι(t) ,
X2(t), and X3(t), through the channel as would be appreciated by those having ordinary skill in the art. For purposes of describing receiving device 30, cyclic module 24a, cyclic module 24b, and cyclic module 24c as illustrated in FIG. 2 will be viewed as adding cyclic prefixes to weighted time-domain signals xι[m], x∑tm], and X3[m], respectively. Also, for clarity, the operation of receiving device 30 will be described with N symbols. However, other embodiments of a receiving device in accordance with the present invention can utilize any multiple of N (e.g., 2/V, 3/V, etc.). From such description, those having ordinary skill in the art will appreciate that alternative embodiments of a receiving device in accordance with the present invention can have any number of receive antennas MR.
In response to time-domain waveforms f (t) , r2(f)and r3{t), a matched filter 32a, a matched filter 32b, and a matched filter 32c (where in one embodiment, each has an impulse response pSr(-t), where psr{t) is the impulse response of each of the pulse shaping modules 25a, 25b, and 25c) provide a matched-filtered waveform rι(t), a matched-filtered waveform r2(t), and a matched-filtered waveform r3(t), respectively. In one embodiment, matched-filtered waveform n(t), matched-filtered waveform r2(t), and matched-filtered waveform 3(t) are determined in accordance with the following equation [14]:
Mγ Λ/+ -1
'7(0= ∑ ∑*j[n1 ij{t -nTs) + nj(t) 7=1 n=0
[14]
where hjj(t) is the overall impulse response between the weighted redundant time-domain signal Xj[ ] and the matched-filtered waveform r t) on the /* antenna as would be appreciated by those skilled in the art. In equation [14], n t) is the receiving device noise signal on the /* antenna branch, h^t) includes the transmit pulse shaping for transmit antenna y", the received matched filtering for receive antenna /', and the impulse response of the RF channel between they* transmit and the /* receive antenna is given by the following equation [15]:
-1 hjj{t) = ∑hjjWrc(t -τι)
1=0
[15]
where prc{t) is the convolution of psι t) with psr(-t), and prc(f) typically has the raised-cosine spectrum. After matched-filtering, a sampling module 33a, a sampling module
33b, and a sampling module 33c sample matched-filtered waveform n(t), matched-filtered waveform r (t), and matched-filtered waveform r3(t), respectively, at a symbol rate to form a received sequence rι[m], a received sequence r∑im], and a received sequence r3[m], respectively. In one embodiment, received sequence rι[m], received sequence r2[m], and received sequence ^Im] are formed in accordance with the following equation [16] (0<m<Λ/-1):
Mj N+Lcp- M-r / +L -l η[m] = η(mTs) = ∑ ∑*j[n ijj(mTs - nTs) = ∑ ∑xj[n ijj[m - n] y'=1 n=0 y'=1 π=0
[16]
where hij[m] is defined to be h mTs) and noise has been neglected.
At each receive branch, the first Lcp received samples for each baud interval are discarded, and the remaining N symbols of each block are transformed into the frequency-domain. Accordingly, a cyclic removal module 34a, a cyclic removal module 34b, and a cyclic removal module 34c provide a symbol-spaced sample yι[m], a symbol-spaced sample y2[m], and a symbol- spaced sample y^m], respectively, where /τ?=0...Λ/-1. In one embodiment, symbol-spaced sample yι[m], symbol-spaced sample y∑ m], and symbol- spaced sample y3[m] are determined in accordance with the following equation [17]:
MT Lcp-l VI M = ∑ ∑ hjj [n]xj [m - n] + n; [m] y'=1 n=0
[17]
where AT7=0.../V-1 , /7;[m] is the sampled noise waveform on the /* receive antenna, and for one embodiment, the following assumptions are made. First, the raised cosine waveform pκ{f) is assumed to be time limited for the interval -NpTs≤t≤NpTs, where Np is the number of significant "tails" of the raised cosine waveform. Second, in the RF channel response as defined in equation [13] above, the first multipath arrival time is NPTS, and the maximum multipath arrival time π is assumed to be no greater than {Lcp - 2/Vp) symbol times, where Lcp is the length of the prefix. These assumptions imply that h, ri = h^nTs) is nonzero only for 0</< cp-1. As will be appreciated by those skilled in the art, the extension of equation [17] to situations where these assumptions do not hold is straightforward.
In response to symbol-spaced sample yι[m], symbol-spaced sample y2[m], and symbol-spaced sample yβlm], respectively, where m=0...Λ/-1 , a TFT module 35a, a TFT module 35b, and a TFT module 35c provide a frequency-domain sample Yι[k], a frequency-domain sample Y∑fk], and a
frequency-domain sample Ys[k], respectively, for /c=0...Λ/-1. In one embodiment, TFT modules 35a-35c are FFT based modules whereby frequency-domain sample Yι[k], a frequency-domain sample Y∑[k], and a frequency-domain sample Y3[k] are determined in accordance with the following equation [18]:
[18]
where
[19]
[20]
Λ/-1 Ni[k] = ∑ nj[n]e-J2πnk l N n=0
[21]
Equation [18] can be concisely written as the following equation [22]:
[22]
where the MR*MT channel matrix on subcarrier k is in accordance with the following equations [23] and [24]:
[k] . H ] H2[ ] - HMτ[k]
[23] and
[24]
Using equation [4], equation [22] can be rewritten as the following equation [25]:
Y[/ ] = H[/ ]v[/ ]sM+N[/ ]
[25] For linear receive array processing, a frequency-domain weighting module 36 takes the vector of frequency-domain samples on each subcarrier (as represented in equation [22]) and applies a linear weighting matrix W[/] to form the vector of multiple outputs on each subcarrier in accordance with the following equation [26]:
Z[k] = WH[/ ]Y[/ ] = WH[/f l[ ]v[ ]s[ ]+ WH[/ ]N[/ ]
[26]
In one embodiment, as illustrated in FIG. 5, frequency-domain weighting module 36 includes a plurality of split modules 36a-36i and a plurality of weight and sum ("W&S") modules 36j-36r in communication via signal paths SP28-SP54.
A FTT module 37a converts weighted frequency-domain samples Zι[0], Zι[1], and Zι[2] into a weighted time-domain sample Zι[m], for A77=O...Λ/-1. A FTT module 37b converts weighted frequency-domain samples Z2/0 , Z2[1], and Z [2] into weighted time-domain sample z2[m], for m=0...N- 1. A FTT module 37c converts weighted frequency-domain samples Z3[0], Z3[1] and Z3[2] into a weighted time-domain sample Z3[m], for m=0...Λ/-1. In one embodiment, FTT module 37a, FTT module 37b, and FTT module 37c are inverse FFT based modules whereby weighted time-domain sample Zι[m], z∑ m], and z3[m] are determined in accordance with the following equation [27]:
[27]
A decoder 38a conventionally decodes weighted time-domain sample
Zι[m] to derive time-domain signal Sι[m] (FIG. 2). A decoder 38b conventionally decodes weighted time-domain sample z2[m] to derive time- domain signal s2[m] (FIG. 2). A decoder 38c conventionally decodes weighted time-domain sample Z3[m] to derive time-domain signal S3[m] (FIG. 2).
Each module of transmitting device 20 (FIG. 2) and each module of receiving device 30 (FIG. 4) as described may be implemented in hardware (analog or digital), software, or any combination of hardware and software. Those having ordinary skill in the art will appreciate a sequential operation of the modules of transmitting device 20 and the modules of receiving device 30 (e.g., in a software implementation) and a concurrent operation of the modules of transmitting device 20 and the modules of receiving device 30 (e.g., in a hardware implementation).
Numerous methods can be used for picking the transmit weight matrix implemented by frequency-domain weighting module 22 (FIGS. 2 and 3) and the receive weight matrix implemented by frequency-domain weighting
module 36 (FIGS. 4 and 5) for frequency flat channels. Such methods can be applied on a subcarrier in the frequency-domain of a cyclic-prefix-based single carrier system. For illustration purposes, several examples for picking the transmit and receive weights are given in the following discussion.
For an arbitrary transmit weight matrix V[/], the MMSE receive weights can be in accordance with the following equation [28]:
TJfc] = {lφ]v[ ]V* [ifcjH" [k]+ σ2ι}_1ιφ]v| ]
[28]
When channel knowledge is not available to transmitting device 20 (FIG. 2), the transmit weight matrix V[/c] can be chosen to be equal to the identity matrix. Another example of choosing the transmit weight matrix is as follows. First, take the singular value decomposition of H[ ] in accordance with the following equation [29]:
H[k]= U[k]∑[k]QH[k]
[29]
where ∑[k] is a diagonal matrix with the singular values of H[/ ] along the main diagonal. Setting W(/ ) = U[/]∑_1 and V[/] = Q[A], equation [26] becomes the following equation [30]:
Z[k] = S[k]+Σ~1[k]u[k)M[k] [30]
and the channel has effectively been "diagonalized" into parallel non- interfering spatial channels between the input to the transmit array and the input to the receive array. In this embodiment, the transmitter must first be provided channel information such as either H[/] or Q[ ]. In one embodiment,
the channel matrices H[/ ] are estimated at the receiving device and then transmitted back to the transmitting device.
The term "signal" was used in the operational description of the transmitting device 20 (FIG. 2) herein and term "sample" was used in the operational description of the transmitting device 20 (FIG. 4) herein for purposes of making a clear distinction between the operational description of transmitting device 20 and the operational description of receiving device 30. However, for purposes of claim construction, the terms "signal" and "sample" are interchangeable. The present invention may be embodied in other specific forms without departing from its spirit or essential characteristics. The described embodiments are to be considered in all respects only as illustrative and not restrictive. The scope of the invention is, therefore, indicated by the appended claims rather than by the foregoing description. All changes that come within the meaning and range of equivalency of the claims are to be embraced within their scope.