The field of application of the invention is the transmission of a digital signal, the invention relates especially to the block coding of a signal for transmission through more than one antenna.
As is well known, the transmission path used for the transmission of signals in telecommunications connections causes interference to telecommunications. This occurs regardless of the physical form of the transmission path, whether it is a radio link, optical fibre or copper cable. Especially in radio communications, there are often situations, in which the quality of the transmission path varies from one connection to another and even during a connection. Fading of the radio path is one typical phenomenon that causes changes in the transmission channel. Other simultaneous connections can also cause interference that may change with time and location. In a typical radio communications environment, signals between a transmitter and receiver propagate along several routes. This multipath propagation occurs mainly due to the fact that the signal reflects from the surrounding surfaces. The signals that have propagated along different routes arrive at the receiver at different times due to their different propagation delays. Different methods have been developed to compensate for the fading caused by multipath propagation.
One of the most effective ways to compensate for the fading on the radio path is controlling the transmission power of the transmitter. If the properties of the radio path are known, the power of the transmitter can be controlled in such a manner that the effect of the fading can be eliminated. In practice, such a solution is, however, very difficult to implement, because firstly, the transmitter should know the quality of the channel, and transmitting this information to the transmitter in real time is difficult. Secondly, transmission power restrictions set on transmitters and the dynamic range of transmitters set their own limitations. In addition, power control in itself may lead to ineffective transmission by increasing power to a high level in fading dips. A second solution to the problem is using diversity in the transmitter. Time diversity employs interleaving and coding to provide temporal diversity to the transmitted signal. However, this has the drawback that it causes delays in the transmission especially when the channel is a slow-fading one. In frequency diversity, the signal is transmitted on several frequencies at the same time. This is, however, an ineffective method when the coherence bandwidth of the channel is wide.
In transmission antenna diversity, the same signal or different parts of the same signal are transmitted to the receiver by using two or more antennas. The signal components that multipath-propagated through different channels will then probably not suffer from interference from a simultaneous fading. Publication WO 99/14871 discloses a diversity method especially suited for two antennas, in which the symbols to be transmitted that are made up of bits are coded in blocks of given lengths and in which each block is coded to contain a given number of channel symbols according to formula (1).
In the formula, the horizontal lines of the matrix show transmission time instants so that the upper horizontal line shows information transmitted at time instant t and the lower horizontal line shows information transmitted at time instant t+T, wherein T refers to a symbol cycle. The vertical lines of the matrix in turn show antennas so that the first vertical line shows information transmitted through antenna 1 and the second vertical line shows information transmitted through antenna 2. However, a block code of complex modulation as shown in formula (1) only exists for a maximum of two antennas. A maximum data transmission rate, i.e. coding rate, transmitted to more than two antennas is calculated according to formula (2), where N is the number of transmission antennas and the square brackets indicate the smallest integer that is greater than or equal to the expression in brackets. It should be noted that herein the coding rate refers to a symbol-coding rate, i.e. the number of symbols transmitted during a symbol cycle T. For a third and fourth antenna, the maximum coding rate for an orthogonal code is ¾.
Publication Tirkkonen, Boariu, Hottinen, IEEE 6th Symposium on Spread-Spectrum Tech. & Appli., NJIT, New Jersey, USA, September 2000, discloses some solutions of a full coding rate 1. The publication describes a coding method, in which a non-orthogonal block code according to formula (3) having four antennas is formed by utilizing an orthogonal Alamout matrix according to formula (1).
The matrix is of form ABBA according to presentation method (4), in which the A matrix follows the Alamout matrix for symbols z1 and z2, whereas B follows the Alamout matrix for symbols z3 and z4.
Known non-orthogonal block codes do, however, contain a significant drawback. Said block codes do not have full diversity, i.e. the number of independently decoded channels is smaller than the number of antennas, whereby transmission capacity provided by the antennas is lost. The diversity of the block code is the smallest number of nonzero eigenvalues.
Hce=DceHDce, (5)
where Dce is defined by formula (6).
Dce=Cc−Ce (6)
In the above, Cc is a transmitted code matrix and Ce a detected defective code matrix. In the case of a channel symbol pair, in which the same error Δ occurs with symbols z1 and z3 and no error occurs with symbols z2 and z4, a matrix according to formula (7) is obtained as the difference matrix Dce.
The matrix according to formula (7) is singular, i.e. the matrix does not have an inverse matrix. The matrix has only two nonzero eigenvalues, 2Δ and 2Δ*. Thus, the diversity degree of the ABBA block code according to formula (3) is only 2. Low diversity begins to show in the coding performance as a decreasing bit error rate when the bit energy to interference level ratio exceeds 5 dB.
It is thus an object of the invention to provide a method and an apparatus implementing the method in such a manner that coding achieves full diversity when two or more transmission antennas are used. This is achieved by a method for transmitting digital symbols, in which method at least two complex symbols belonging to the same transmission block are read in the transmitter. In the method, a non-orthogonal block code comprising channel symbols is formed from the symbols belonging to the transmission block by performing at least one of the following operations to each symbol: repetition, conjugation, multiplying by a weighting coefficient, summing two or more symbols, and multiplying at least one symbol belonging to the transmission block by a first nonzero weighting coefficient and at least one symbol belonging to the transmission block by a second nonzero weighting coefficient, the division ratio of the first and second weighting coefficients differing from values ±1, ±j and ±1/sqrt(2), and transmitting the formed channel symbols via two or more transmission antenna routes.
The invention also relates to an arrangement for transmitting digital symbols, comprising a coder for coding digital symbols into channel symbols, and one or more transmission antennas connected to the coder for transmitting the channel symbols. The coder is arranged to read at least two complex symbols belonging to the same transmission block and, by means of the non-orthogonal block code defining the coding, to form channel symbols of the symbols belonging to the transmission block by performing at least one of the following operations to each symbol: repetition, conjugation, multiplying by a weighting coefficient, summing two or more symbols, the coder being arranged to multiply at least one symbol belonging to the transmission block by a first nonzero weighting coefficient and at least one symbol belonging to the transmission block by a second nonzero weighting coefficient, the division ratio of the first and second weighting coefficients differing from values ±1, ±j and ±1/sqrt(2), and the arrangement being arranged to transmit the formed channel symbols via two or more transmission antenna routes formed by means of said one or more antennas.
Preferred embodiments of the invention are disclosed in the dependent claims.
The invention thus relates to a method and apparatus for performing a space-time block coding in a radio transmitter. A certain number of symbols are read during block coding and coded into channel symbols for transmission at different time instants and through at least two transmission antennas. The number of reception antennas is not essential for the invention. The inventive solution is based on having at least two symbols in the transmission block that are multiplied by nonzero weighting coefficients in such a manner that the division ratio of the coefficients differs from values ±1 and ±j.
The code matrix of the invention fulfils the requirement of maximum symbol-wise diversity. This means that if all other symbols except one were removed from the code matrix, the code matrix would be unitary with respect to the remaining symbol. The above-mentioned condition applies to all code matrix symbols. Thus, equation (8) applies to the symbol-wise code matrix Ck of code matrix C. The symbol-wise code matrix Ck refers herein to a code matrix, in which all other symbols except zk are replaced by zeros.
CkHCk≈|zk|2IN (8)
In equation (8), CH is a Hermitian conjugate of matrix C, i.e. a transpose of a complex conjugate, |z| is the absolute value of z, and IN is an N-dimensional identity matrix. If M=N, the simplest way of detecting symbol-wise unitarity is to form a symbol-wise matrix Ck, in which case exactly one symbol on each horizontal and vertical line is zk or its complex conjugate while the other elements are zero.
A general code matrix that fulfils the full symbol-wise diversity can be formed as follows. An integer nk between 1≦nk≦N, where N is the number of transmission antennas, is selected for each symbol zk. An N×N matrix is formed, in which nk zk's and N−nk z
Ck=UkγkVk (9)
The general code matrix C that fulfils the full symbol-wise diversity is of the form shown in formula (10)
where Uk is a unitary M×M matrix and Vk is a unitary N×N matrix. In the inventive solution, matrices Uk, Vk have at least one element not belonging to the value range {0,±1,±j}.
The code matrix of the invention should preferably be such that the non-orthogonality matrix formed by means of the Hermitian square of the code matrix is traceless. Otherwise, bits coded into some of the symbols would be non-homogeneously coded—for instance an error in the detection might be more probable if bit 0 is transmitted instead of bit 1. The code matrix is also preferably required to have symbol-level homogeneity, i.e. that all symbols have an equal position with respect to each other in the code matrix.
In a preferred embodiment of the invention, all symbols of the block code are formed using the same modulation constellation, but for at least one symbol, the modulation constellation is phase-rotated in relation to the other symbols in the block. For QPSK symbols, for instance, the phase rotation in codes according to equation (3) is preferably between 25 and 65 degrees. In a second preferred embodiment of the invention, two or more symbols are summed when forming the channel symbol of a block-coding matrix. In one preferred embodiment, one symbol is multiplied by a weighting coefficient and then summed with another symbol from one or more blocks to form the channel symbol.
Preferably, at least two symbols belonging to the transmission block are coded orthogonal with each other, even though as a whole, the block code of the invention is non-orthogonal. For instance, in the case of a four-symbol block, one preferred embodiment is a non-orthogonal 2+2 matrix, in which two symbols are coded orthogonal with each other and the two remaining symbols are correspondingly coded orthogonal with each other. Another preferred embodiment is a 3+1 block code, in which three symbols are coded orthogonal with each other. The invention is not limited to whether the block matrix is power-balanced or not.
The invention provides significant advantages. The coding rate of the code matrix can be made higher than the coding rate allowed by orthogonality, and by weighting one or more symbols by a complex weighting coefficient, the invention enables the presentation of a code matrix that achieves full diversity. The solution of the invention provides significant advantages in the level of the bit error rate (BER) in comparison with the prior art solution, especially at high signal-to-noise ratios (SNR).
The invention will now be described in more detail by means of preferred embodiments and with reference to the attached drawings, in which
The invention can be used in radio systems, in which it is possible to transmit at least a part of a signal by using two or more transmission antennas or two or more lobes provided by two or more transmission antennas. The transmission channel can be established using a time-division, frequency-division, or code-division symbol multiplexing or multiple access method. Systems using combinations of different multiple access methods are also systems of the invention. The examples describe the use of the invention in a UMTS (Universal Mobile Telephony System) system employing a wide-band code-division multiple-access method implemented by a direct sequence technique, without limiting the invention to it, however.
The structure of a mobile system is described by way of example with reference to
The description in
The base station 204 further has a control unit 210 that controls the operation of the transceivers 208 and the multiplexer 212. The multiplexer 212 places the traffic and control channels used by several transceivers 208 on one transmission connection 214. The transmission connection 214 forms an interface Iub.
The transceivers 208 of the base station 204 are connected to an antenna unit 218 implementing a bi-directional radio connection 216 to the user equipment 202. The structure of frames transmitted in the bi-directional radio connection 216 is defined specifically for each system and is called an air interface Uu.
The radio network controller 206 comprises a group-switching field 220 and a control unit 222. The group-switching field 220 is used to switch speech and data and to connect signaling circuits. A radio network subsystem 224 formed by the base station 204 and the radio network controller 206 also comprises a transcoder 226. The transcoder 226 is usually located as close to a mobile switching center 228, as possible, because speech can then be transmitted in cellular radio network format between the transcoder 226 and the radio network controller 206, saving transmission capacity. The transcoder 226 transforms the different digital speech coding formats used between a public telephone network and a radio telephone network to be compatible with each other, for instance from a fixed network format to a cellular radio network format and vice versa. The control unit 222 takes care of call control, mobility management, collection of statistics and signaling.
The coder 308 is connected to the transmission antennas 314 to 318 through radio frequency parts 312. In the case of
Next, the invention will be described by means of a non-orthogonal block code designed for N antennas and transmitting K symbols in M symbol cycles, where K is more than M times the maximum orthogonal coding rate allowed by equation (2). The code matrix is required to have full symbol-wise diversity that is described by equations (8) and (9). This means that the code provides a full N-fold diversity protection against bit errors in one symbol, which can be considered a basic requirement for the code providing at least an approximately full diversity protection.
The coding and modulation methods can be taken as selecting points in a point space. Depending on the used method and channels, there always is a metric defining the geometry of the alphabet. Some points of the alphabet are closer to each other than others. The performance depends greatly on the distances between the closest neighbouring points. Optimally, the points have equal positions, so the distance to the closest neighbours should not be dependent on the selected points. This is called the homogeneity principle. From the homogeneity principle follows that traceless non-orthogonality can be required of a block code. The Hermitian square of the block code fulfilling the one-symbol diversity requirement described above can be shown as equation (11).
It is now required of the preferred code matrix that the non-orthogonality matrix X shown in equation (12) is traceless (13).
TrX=0. (13)
The above can be argued as follows. The form of matrix X shows that it is a Hermitian matrix, which means that its diagonal elements are real. Because X is a bilinear combination of matrices Ck that are linear with respect to symbols zk, the real diagonal values of X must be real linear combinations of the following: Re[z
The performance of block codes is greatly dependent on the properties of the squared difference matrix (5) of code words. The trace of matrix Hce is the Euclidian distance between the symbol constellations Ce and Cc. In a linearly coded code, the difference matrices (6) of code words use the equations (11, 12) filled by the code matrices, where symbols zk are replaced by symbol differences
Δk=zk(c)−zk(e). (14)
Here, zk(c) refers to the transmitted symbols and zk(e) to symbols possibly received with errors. This way, the difference matrix of the code words complies with form (14), in which case the squared matrix Hce is of form (15).
The trace of matrix Hce has two parts,
from the first term and Tr(X(Δ
Let us now concentrate on describing preferred embodiments when N=4, M=4 and K=4, i.e. the coding rate is 1. In the preferred embodiment in question, symbol-level homogeneity is required, i.e. that all symbols have an equal position with respect to each other in the code matrix. Further, in the preferred embodiment, the code matrix also has a maximum residual orthogonality. In a four-antenna 2+2 embodiment according to these requirements, two symbols z1 and z2 are coded orthogonally with each other in a pair, as are the remaining two symbols z3 and z4. QPSK, for instance, is used as the modulation method. In the embodiments of the invention, the matrices U and V are selected in such a manner that the code has full diversity, contrary to the prior-art 2+2 code that is according to equation (3).
In this preferred embodiment, the code matrices according to equation (10) are according to equation (16)
where A and B are orthogonal 2*2 matrices in the form of the Alamout matrix. 0 is a 2*2 zero matrix and U and V are unitary 4×4 matrices in such a manner, however, that at least one matrix element in U or V is unequal to 0, ±1 and ±j. Let us now concentrate on the exemplary case, in which matrix V is a 4×4 singular matrix.
In the preferred embodiment, traceless non-orthogonality can be achieved for instance by setting
where W is a unitary 2×2 matrix whose determinant is 1 and q and p are real numbers that fulfil the condition
q2+p2=1. (18)
The codes having different values of q are fully equivalent to each other and the codes have the same eigenvalues of the squared difference matrix of the code words. To optimise performance, it is enough to examine one example of the above mentioned coefficients by setting q=0, whereby the code becomes
From the form of matrix (19), it is possible to form both a power-balanced and a non-power-balanced block code. The power-balanced code means herein that each antenna transmits at a constant power at all time instants, if constant-power modulation, such as QPSK, is used, whereas in non-power-balanced code, the transmission power is not constant at all time instants. In a power-balanced block code, W is expected to be a diagonal matrix, for instance
where φ is selected so as to optimise performance. In the above, W is unitary and its determinant is 1. The code matrix is then in the form of equation (3), except that symbols z3 and z4 are multiplied by phase coefficient λjφ. Each φ≠0 makes the code a non-singular full diversity code of the invention. A second way of interpreting these inventive solutions is to think that in a code according to equation (3), the modulation points of symbols z1 and z2 are selected from {±1,±j}, whereas symbols z3 and z4 are selected from the phase rotated constellation λjφ{1, j,−1,−j}, where each φ≠0 makes the code a non-singular full diversity code of the invention. In a partly non-power-balanced block code, matrix W can be selected to be a general unitary 2×2 matrix whose determinant is 1. These are of form
where α and β belong to complex numbers and the sum of the squares of their absolute values is 1. By combining the W described above with the Alamout block B described in formula (19)
the code matrix can be presented by means of pseudosymbols
{tilde over (z)}3=ejφ(αz
{tilde over (z)}4=ejφ(αz
The code matrix can then be presented as
In the code matrices of the invention, either α or β or both are unequal to 0, ±1 and ±j. In non-power-balanced embodiments, both α and β are nonzero. The non-orthogonality matrix is now
where I2 is a 2×2 singular matrix, and it is clearly traceless.
At a high signal-to-noise ratio, the eigenvalues of the squares Hce of the code matrix difference matrices according to equation (5) determine explicitly the performance of the code. The most important property of the difference matrices Hce is their rotatability, i.e. that all eigenvalues are nonzero. The code then obtains full diversity. This criterion is what is known as the rank criterion of space-time codes, namely that all code matrix pairs Cc and Ce have a difference matrix square Hce of maximum rank. The preferred embodiments of the present invention, for instance those according to equation (16), where U and V are according to the invention, fulfil the rank criterion; no square Hce of a code word difference matrix is singular, and the codes provide full diversity.
When the ranks of the matrices Hce are maximized, the distribution of their eigenvalues can next be maximized so as to provide as good a performance as possible. At high signal-to-noise ratios, this is achieved by maximizing the determinants of the matrices Hce. According to the determinant criterion of the space-time codes, the smallest possible determinant of Hce should be maximized, when taking all code matrix pairs Cc and Ce into consideration. This is what is known as the MAX MIN DET criterion. In the preferred embodiments of the invention, the angle φ in equation (20), or in the non-power-balanced embodiments, the complex numbers α and β in equation (21) can now be selected to be in accordance with the MAX MIN DET criterion.
If a lower signal-to-noise ratio is used, it is possible to minimize the bit error rate, for instance, with code matrices of a given form instead of using the determinant criterion.
The eigenvalues of the square Hce of a four-antenna block code difference matrix and thus also the determinant and bit error rate of Hce can be calculated exactly. The eigenvalues can be calculated from the traces of the first, second, third and fourth power of Hce. The invariants t1 to t4 thus obtained for the code matrices according to the preferred embodiments of equation (16) are presented in equations (27). The thus formed eigenvalues of matrix Hce according to equation (5) are presented in equations (28) and the determinant is presented in equation (29).
t
t
t
In the above, Δi=zic−zie are symbol differences. In a non-power-balanced case according to equation (22), Δ3 and Δ4 should be interpreted as differences of pseudo-symbols {tilde over (z)}3 and {tilde over (z)}4. Equation (29) shows that the determinants of symbol difference pairs Δ1, Δ3 and Δ2, Δ4 can be maximized separately. If Δ2=Δ4=0, it results in
det[H]=|Δ
According to the rank criterion, the modulation constellations are selected such that the determinant is never zero. The determinant is zero only if Δ1=Δ3.
Thus, it is clear that the modulation constellations for symbols z1 and z3 should be selected in such a manner that they do not overlap. This can be done for instance by selecting the constellations with different phase factors φ, as done in the preferred embodiments shown in equations (19, 20). In the case of QPSK modulation presented by way of example, the results of different optimisations can be seen in
In the phase factor φ shown in equation (31), the minimum distance between points in a squared symbol difference alphabet for symbols z1 and z3 is at most as in the presentation 410 of
A second way of ensuring that the determinant (29) is never lost is to select a non-power-balanced embodiment (19, 21). If both α and β are nonzero, possible errors occurring in symbols z1 and the pseudo-symbols {tilde over (z)}3 of equation (23) can never be the same. The complex variables α and β in W can be presented by means of three Euler angles.
α=ei(φ+Φ)cos θ
β=ei(φ−Φ)sin θ. (32)
The union bound limit of the bit error rate is minimized by selecting
φ=π/8
Φ=π/8
θ=π/5. (33)
There are also some other angles that provide the same performance, but listing them herein is not essential.
It is clear that the method presented herein is not limited to the QPSK modulation method. For instance the MAX MIN DET calculated for the 16QAM modulation method is achieved by phase factor φ≈0.172π. A second advantageous phase rotation for the 16QAM modulation method is about φ≈π/4≈45°, which is again technically easier to implement than the optimal φ≈0.172π.
The above describes the phase rotation of one symbol in the case of a 2+2 ABBA code. Equation (34) shows a non-orthogonal 3+1-form block code, in which symbols z1, z2 and z3 are coded orthogonally with respect to each other, but symbol z4 is not coded orthogonally with respect to any of the other symbols.
The code in formula (34) is thus not homogeneous at symbol level. In code matrix (34), the modulation points of symbols z1, z2 and z3 preferably belong to the value range {+/−1, +/−j}. The optimal modulation points of symbol z4 depend on the desired SNR, when the bit error rate is minimized. It has been noted that by minimizing the union bound limit of the bit error rate, the optimal phase shift for symbol z4 is +/−29 degrees at an average signal-to-noise ratio of 10 dB. In practice, the phase shift of the alphabet of symbol z4 can be selected freely between 25 and 65 degrees without a significant impact on the performance.
A non-orthogonal 2+2 code optimized on the basis of the results obtained from simulations is always better with respect to bit error rate minimization than the known ABBA code shown in formula (3).
Let us now examine the receiver shown in
From the combiner, the signal is taken to a detector that detects symbols using detection methods. It is, for instance, possible to calculate the Euclidian distance of combined symbol estimates from possible symbol states or define the a posteriori probabilities of received symbols or bits. In the latter case, information on the channel needs to be obtained from the channel estimator. From the detector, the signal is taken to a channel decoder and on to other parts of the receiver. Symbol or bit interleaving and deinterleaving performed in both the transmitter and receiver and a possible channel coding are missing from the above description. These can be performed by known methods, if necessary. The above description is only one example of a possible receiver. The calculation and use of channel estimates, for instance, can be implemented in many ways, as is apparent to a person skilled in the art.
The facilities executing the steps of the invention in the coder and other parts of equipment belonging to the transmitter and arrangement are preferably implemented by program by means of a processor and suitable software in both the transmitting and receiving end. The facilities can also be implemented by means of separate components or circuits, for instance.
Even though the invention has been explained in the above with reference to examples in accordance with the accompanying drawings, it is apparent that the invention is not restricted to them but can be modified in many ways within the scope of the inventive idea disclosed in the attached claims.
Number | Date | Country | Kind |
---|---|---|---|
20002845 | Dec 2000 | FI | national |
Filing Document | Filing Date | Country | Kind | 371c Date |
---|---|---|---|---|
PCT/FI01/01133 | 12/19/2001 | WO | 00 | 11/13/2003 |
Publishing Document | Publishing Date | Country | Kind |
---|---|---|---|
WO02/058311 | 7/25/2002 | WO | A |
Number | Name | Date | Kind |
---|---|---|---|
5848103 | Weerackody | Dec 1998 | A |
6088408 | Calderbank et al. | Jul 2000 | A |
6097771 | Foschini | Aug 2000 | A |
6178196 | Naguib et al. | Jan 2001 | B1 |
6307851 | Jung et al. | Oct 2001 | B1 |
6317411 | Whinnett et al. | Nov 2001 | B1 |
6317466 | Foschini et al. | Nov 2001 | B1 |
6359874 | Dent | Mar 2002 | B1 |
6445730 | Greszczuk et al. | Sep 2002 | B1 |
6631168 | Izumi | Oct 2003 | B2 |
6741658 | Ionescu | May 2004 | B1 |
6760388 | Ketchum et al. | Jul 2004 | B2 |
6865237 | Boariu et al. | Mar 2005 | B1 |
6922447 | Ionescu | Jul 2005 | B1 |
7006848 | Ling et al. | Feb 2006 | B2 |
7061854 | Tarokh et al. | Jun 2006 | B2 |
20010040928 | Sakoda | Nov 2001 | A1 |
20030147343 | Onggosanusi et al. | Aug 2003 | A1 |
20040071240 | Betts | Apr 2004 | A1 |
20040243904 | Kim et al. | Dec 2004 | A1 |
20050190853 | Tirkkonen | Sep 2005 | A1 |
20050201481 | Calderbank et al. | Sep 2005 | A1 |
Number | Date | Country |
---|---|---|
0 905 920 | Mar 1999 | EP |
2237706 | May 1991 | GB |
WO 9741670 | Nov 1997 | WO |
WO 9914871 | Mar 1999 | WO |
WO 9923766 | May 1999 | WO |
WO 0011806 | Mar 2000 | WO |
WO 0018056 | Mar 2000 | WO |
WO 0049780 | Aug 2000 | WO |
WO 0051265 | Aug 2000 | WO |
WO 0105060 | Jan 2001 | WO |
WO 0119013 | Mar 2001 | WO |
WO 0163826 | Jun 2001 | WO |
WO 0154305 | Jul 2001 | WO |
WO 0156218 | Aug 2001 | WO |
WO 0169814 | Sep 2001 | WO |
WO 0176094 | Oct 2001 | WO |
WO 0243313 | May 2002 | WO |
Number | Date | Country | |
---|---|---|---|
20040081254 A1 | Apr 2004 | US |