Space-time coding/decoding method for frequency selective fading channel

<u>Field of the Technology</u>

0001The invention relates to mobile communication technology, more particularly to a method of space-time encoding and decoding for a frequency selective fading channel.

<u>Background of the Invention</u>

0002Space-time encoding is one of the important methods that can raise performance of a wireless communication system effectively. There are two branches of the space-time encoding, one is layered space-time encoding and the other is space-time encoding based on transmitting diversity. The space-time encoding based on transmitting diversity can be further divided into two kinds: space-time block encoding and space-time trellis encoding.

0003Along with the development of high-speed wireless communication technology, the signal transmission bandwidth and rate is increased continuously, which means that time for transmitting a data symbol is shorter and shorter. Therefore, time delay spread of a wireless transmission channel cannot be neglected, that is, the channel frequency selective fading is getting worse.

0004The prior space-time encoding and decoding schemes are designed based on a flat fading channel. When the schemes are used for a frequency selective fading channel, the performance becomes worse obviously, and it is too complex to implement when improving performance of the algorithm.

0005In a multi-user frequency selective fading channel, since Inter-Symbol Interference (ISI) and Multiple Access Interference (MAI) exist at the same time, implementation of space-time encoding and decoding has great difficulty. Some effective methods for counteracting ISI and MAI on a frequency selective channel, such as equalization, joint detection method etc., become very complex or even failure if space-time encoding is involved therein.

0006In the 3rd Generation Partnership Project - Time Division Duplex (3GPP-TDD) system, the prior space-time encoding and decoding scheme used for a multi-user frequency selective fading channel only performs space-time encoding for the basic common control channel, but in the 3rd Generation Partnership Project - Frequency Division Duplex (3GPP-FDD) system, space-time encoding is employed for most channels. The 3GPP-TDD space-time encoding scheme first takes a symbol as an encoding and decoding unit, later it takes a half of symbols of a data field as an encoding and decoding unit. Since these schemes have many changes that bring more complex compared with the schemes in which space-time encoding is not employed, the space-time transmission diversity (STTD) scheme is discarded in 3GPP-TDD v4.3 and space code transmission diversity (SCTD) scheme is applied. Although the reception processing procedure is simplified with the SCTD, more channel resource is occupied and more channels are employed; it is difficult to be applied to other kinds of channels, and also it is impossible to be applied on the multipath diversity situation.

0007There is a space-time encoding and decoding scheme that transforms high-speed data to multiple parallel low-speed data and transmits them on multiple channels, but really this is a space-time - orthogonal frequency division multiplexing (ST-OFDM) method that would thoroughly change the signal structure and system on the physical layer, so it is greatly limited in practical use.

0008In summary, the space-time encoding and decoding method is important to improve performance of a wireless communication system. Usually, a symbol or a string of symbols are used as a processing unit for encoding and decoding, but under the condition of a frequency selective fading channel, improvement of performance is limited by the calculation complexity.

0009<patcit id="pcit0001" dnum="EP1069707A"><text>EP-A-1 069 707</text></patcit> discloses a space time transmit diversity scheme using two antenna in which a data block of symbols (D22) is transmitted simultaneously with a complex conjugate of another block (D21*) With suitable receiver processing, the symbols can be recovered with only a small complexity increase.

0010"Computationally efficient multiuser detectors" discloses the use of iterative techniques to solve a symmetric system, which can reduce the computational complexity.

<u>Summary of the Invention</u>

0011An object of the invention is to design a space-time encoding and decoding method that is suitable to the frequency selective fading channel to obtain better performance with small amount of calculation complexity.

0012A space-time encoding and decoding method for a frequency selective fading channel according to the present invention comprises: <ul id="ul0001" list-style="none" compact="compact"><li>A. an encoder taking two independent data fields of a time slot in input data as a processing unit with space-time orthogonal encoding method, encoding them and generating two data vectors, thereby forming two diversity signals, and transmitting said two diversity signals simultaneously with each through one diversity antenna;</li><li>B. a terminal receiving said two diversity signals, and neglecting mutual interference between said two diversity signals caused by non-orthogonality; step B includes: <ul id="ul0002" list-style="none" compact="compact"><li>setting the upper right block and the lower left block of matrix <maths id="math0001"><math display="inline"><msup><mi mathvariant="normal">A</mi><mrow><mo mathvariant="normal">*</mo><mi mathvariant="normal">T</mi></mrow></msup><mo>⁢</mo><mi mathvariant="normal">A</mi><mo mathvariant="normal">=</mo><mfenced open="[" close="]"><mtable><mtr><mtd><msubsup><mi mathvariant="normal">A</mi><mn mathvariant="normal">1</mn><mrow><mo mathvariant="normal">*</mo><mi mathvariant="normal">T</mi></mrow></msubsup><mo>⁢</mo><msub><mi mathvariant="normal">A</mi><mn mathvariant="normal">1</mn></msub><mo mathvariant="normal">+</mo><msup><mfenced><msubsup><mi mathvariant="normal">A</mi><mn mathvariant="normal">2</mn><mrow><mo mathvariant="normal">*</mo><mi mathvariant="normal">T</mi></mrow></msubsup><mo>⁢</mo><msub><mi mathvariant="normal">A</mi><mn mathvariant="normal">2</mn></msub></mfenced><mo mathvariant="normal">*</mo></msup></mtd><mtd><msup><mtable><mtr><mtd><mfenced><msubsup><mi mathvariant="normal">A</mi><mn mathvariant="normal">1</mn><mrow><mo mathvariant="normal">*</mo><mi mathvariant="normal">T</mi></mrow></msubsup><mo>⁢</mo><msub><mi mathvariant="normal">A</mi><mn mathvariant="normal">2</mn></msub></mfenced></mtd></mtr></mtable><mi mathvariant="normal">T</mi></msup><mo mathvariant="normal">-</mo><msubsup><mi mathvariant="normal">A</mi><mn mathvariant="normal">1</mn><mrow><mo mathvariant="normal">*</mo><mi mathvariant="normal">T</mi></mrow></msubsup><mo>⁢</mo><msub><mi mathvariant="normal">A</mi><mn mathvariant="normal">2</mn></msub></mtd></mtr><mtr><mtd><msup><mfenced><msup><mfenced><msubsup><mi mathvariant="normal">A</mi><mn mathvariant="normal">1</mn><mrow><mo mathvariant="normal">*</mo><mi mathvariant="normal">T</mi></mrow></msubsup><mo>⁢</mo><msub><mi mathvariant="normal">A</mi><mn mathvariant="normal">2</mn></msub></mfenced><mi mathvariant="normal">T</mi></msup><mo mathvariant="normal">-</mo><msubsup><mi mathvariant="normal">A</mi><mn mathvariant="normal">1</mn><mrow><mo mathvariant="normal">*</mo><mi mathvariant="normal">T</mi></mrow></msubsup><mo>⁢</mo><msub><mi mathvariant="normal">A</mi><mn mathvariant="normal">2</mn></msub></mfenced><mrow><mo mathvariant="normal">*</mo><mi mathvariant="normal">T</mi></mrow></msup></mtd><mtd><msup><mtable><mtr><mtd><mfenced><msubsup><mi mathvariant="normal">A</mi><mn mathvariant="normal">1</mn><mrow><mo mathvariant="normal">*</mo><mi mathvariant="normal">T</mi></mrow></msubsup><mo>⁢</mo><msub><mi mathvariant="normal">A</mi><mn mathvariant="normal">1</mn></msub><mo mathvariant="normal">+</mo><msup><mfenced><msubsup><mi mathvariant="normal">A</mi><mn mathvariant="normal">2</mn><mrow><mo mathvariant="normal">*</mo><mi mathvariant="normal">T</mi></mrow></msubsup><mo>⁢</mo><msub><mi mathvariant="normal">A</mi><mn mathvariant="normal">2</mn></msub></mfenced><mo mathvariant="normal">*</mo></msup></mfenced></mtd></mtr></mtable><mo mathvariant="normal">*</mo></msup></mtd></mtr></mtable></mfenced></math><img file="EP1533928B1_D0001.tif" /></maths> to be null matrixes, and then calculating equation <b>d̂_t=(B)^-1A^-Tr</b> to obtain a simplified equation for joint detection; wherein <b>A₁</b> and <b>A₂</b> are system matrixes of signal transmission between first and second transmitting antennas and receiving antennas; <b>A</b> and <b>B</b> are matrixes; <b>d̂_t</b> is a value of continuous estimation of a receiving data field; <b>r</b> is a sample value of said receiving data field; T denotes a transpose operation; * denotes conjugate; wherein said matrix <b>B</b> is calculated by one of the following formulas: <maths id="math0002"><math display="block"><mi mathvariant="bold">B</mi><mo>=</mo><mrow><mo>{</mo><mtable><mtr><mtd><mi mathvariant="bold">I</mi></mtd><mtd><mi mathvariant="italic">MF</mi></mtd></mtr><mtr><mtd><msup><mi mathvariant="bold">A</mi><mrow><mo mathvariant="bold">*</mo><mi mathvariant="bold">T</mi></mrow></msup><mo>⁢</mo><mi mathvariant="bold">A</mi></mtd><mtd><mi mathvariant="italic">ZF</mi><mo mathvariant="italic">-</mo><mi mathvariant="italic">BLE</mi></mtd></mtr><mtr><mtd><msup><mi mathvariant="bold">A</mi><mrow><mo mathvariant="bold">*</mo><mi mathvariant="bold">T</mi></mrow></msup><mo>⁢</mo><mi mathvariant="bold">A</mi><mo mathvariant="bold">+</mo><msup><mi>σ</mi><mn mathvariant="bold">2</mn></msup><mo>⁢</mo><mi mathvariant="bold">I</mi></mtd><mtd><mi mathvariant="italic">MMSE</mi><mo mathvariant="italic">-</mo><mi mathvariant="italic">BLE</mi></mtd></mtr></mtable></mrow></math><img file="EP1533928B1_D0002.tif" /></maths></li><li>wherein MF represents a match filter scheme, ZF-BLE represents a zero-forcing block equalization scheme and MMSE - BLE represents a minimum mean-square-error block equalization scheme; wherein ² is noise power, and <b>I</b> is an identity matrix;</li></ul></li><li>C. said terminal performing joint detection only taking into account effect to said two diversity signals from multipath interference and multi-user interference, thereby obtaining a decoding result; said joint detection is calculated based on a simplified joint detection equation: <maths id="math0003"><math display="inline"><mrow><mo mathvariant="bold">{</mo><mtable columnalign="left"><mtr><mtd><mover><mi mathvariant="bold">d</mi><mo mathvariant="bold">^</mo></mover><mfenced><mn mathvariant="bold">1</mn></mfenced></mtd><mtd><mo mathvariant="bold">=</mo><msubsup><mi mathvariant="bold">B</mi><mi mathvariant="bold">S</mi><mrow><mo mathvariant="bold">-</mo><mn mathvariant="bold">1</mn></mrow></msubsup><mo>⁢</mo><mfenced><msubsup><mi mathvariant="bold">A</mi><mn mathvariant="bold">1</mn><mrow><mo mathvariant="bold">*</mo><mi mathvariant="bold">T</mi></mrow></msubsup><mo>⁢</mo><msub><mi mathvariant="bold">r</mi><mn mathvariant="bold">1</mn></msub><mo mathvariant="bold">+</mo><msup><mfenced><msubsup><mi mathvariant="bold">A</mi><mn mathvariant="bold">2</mn><mrow><mo mathvariant="bold">*</mo><mi mathvariant="bold">T</mi></mrow></msubsup><mo>⁢</mo><msub><mi mathvariant="bold">r</mi><mn mathvariant="bold">2</mn></msub></mfenced><mo mathvariant="bold">*</mo></msup></mfenced></mtd></mtr><mtr><mtd><mover><mi mathvariant="bold">d</mi><mo mathvariant="bold">^</mo></mover><mfenced><mn mathvariant="bold">2</mn></mfenced></mtd><mtd><mo mathvariant="bold">=</mo><msubsup><mi mathvariant="bold">B</mi><mi mathvariant="bold">S</mi><mrow><mo mathvariant="bold">-</mo><mn mathvariant="bold">1</mn></mrow></msubsup><mo>⁢</mo><mfenced><msubsup><mi mathvariant="bold">A</mi><mn mathvariant="bold">1</mn><mrow><mo mathvariant="bold">*</mo><mi mathvariant="bold">T</mi></mrow></msubsup><mo>⁢</mo><msub><mi mathvariant="bold">r</mi><mn mathvariant="bold">2</mn></msub><mo mathvariant="bold">-</mo><msup><mfenced><msubsup><mi mathvariant="bold">A</mi><mn mathvariant="bold">2</mn><mrow><mo mathvariant="bold">*</mo><mi mathvariant="bold">T</mi></mrow></msubsup><mo>⁢</mo><msub><mi mathvariant="bold">r</mi><mn mathvariant="bold">1</mn></msub></mfenced><mo mathvariant="bold">*</mo></msup></mfenced></mtd></mtr></mtable></mrow><mo mathvariant="bold">,</mo></math><img file="EP1533928B1_D0003.tif" /></maths>wherein <b>d̂(1)</b> and <b>d̂(2)</b> are values of continuous estimation of two receiving data fields, <b>B_s</b> is a matrix; r₁, and r₂ arc sample values of two receiving data fields; wherein said matrix <b>B_s</b> is calculated by one of the following formulas: <maths id="math0004"><math display="block"><msub><mi mathvariant="bold">B</mi><mi mathvariant="bold">S</mi></msub><mo>=</mo><mrow><mo>{</mo><mtable><mtr><mtd><mi mathvariant="bold">I</mi></mtd><mtd><mi mathvariant="italic">MF</mi></mtd></mtr><mtr><mtd><msubsup><mi mathvariant="bold">A</mi><mn mathvariant="bold">1</mn><mrow><mo mathvariant="bold">*</mo><mi mathvariant="bold">T</mi></mrow></msubsup><mo>⁢</mo><msub><mi mathvariant="bold">A</mi><mn mathvariant="bold">1</mn></msub><mo mathvariant="bold">+</mo><msup><mfenced><msubsup><mi mathvariant="bold">A</mi><mn mathvariant="bold">2</mn><mrow><mo mathvariant="bold">*</mo><mi mathvariant="bold">T</mi></mrow></msubsup><mo>⁢</mo><msub><mi mathvariant="bold">A</mi><mn mathvariant="bold">2</mn></msub></mfenced><mo mathvariant="bold">*</mo></msup></mtd><mtd><mi mathvariant="italic">ZF</mi><mo mathvariant="italic">-</mo><mi mathvariant="italic">BLE</mi></mtd></mtr><mtr><mtd><msubsup><mi mathvariant="bold">A</mi><mn mathvariant="bold">1</mn><mrow><mo mathvariant="bold">*</mo><mi mathvariant="bold">T</mi></mrow></msubsup><mo>⁢</mo><msub><mi mathvariant="bold">A</mi><mn mathvariant="bold">1</mn></msub><mo mathvariant="bold">+</mo><msup><mfenced><msubsup><mi mathvariant="bold">A</mi><mn mathvariant="bold">2</mn><mrow><mo mathvariant="bold">*</mo><mi mathvariant="bold">T</mi></mrow></msubsup><mo>⁢</mo><msub><mi mathvariant="bold">A</mi><mn mathvariant="bold">2</mn></msub></mfenced><mo mathvariant="bold">*</mo></msup><mo mathvariant="bold">+</mo><msup><mi>σ</mi><mn mathvariant="bold">2</mn></msup><mo>⁢</mo><mi mathvariant="bold">I</mi></mtd><mtd><mi mathvariant="italic">MMSE</mi><mo mathvariant="italic">-</mo><mi mathvariant="italic">BLE</mi></mtd></mtr></mtable></mrow></math><img file="EP1533928B1_D0004.tif" /></maths> wherein MF represents a match filter scheme, ZF-BLE represents a zero-forcing block equalization scheme and MMSE - BLE represents a minimum mean-square-error block equalization scheme; wherein σ² is noise power, and <b>I</b> is an identity matrix; and</li><li>D. performing interference counteraction based on result of joint detection to remove interference between two diversity signals, and then returning to step C to implement iteration for decoding processing; the step of performing interference couateraction based on result of joint detection comprises:. <ul id="ul0003" list-style="none" compact="compact"><li>D1. subtracting effect of a data field <b>d(1)</b> from received data signal based on the following formula, <maths id="math0005"><math display="block"><mrow><mo mathvariant="bold">{</mo><mtable columnalign="left"><mtr><mtd><msubsup><mi mathvariant="bold">r</mi><mn mathvariant="bold">1</mn><mi mathvariant="bold">ʹ</mi></msubsup></mtd><mtd><mo mathvariant="bold">=</mo><msub><mi mathvariant="bold">r</mi><mn mathvariant="bold">1</mn></msub><mo mathvariant="bold">-</mo><msub><mi mathvariant="bold">A</mi><mn mathvariant="bold">1</mn></msub><mo>⁢</mo><mover><mi mathvariant="bold">d</mi><mo>^</mo></mover><mfenced><mn mathvariant="bold">1</mn></mfenced></mtd></mtr><mtr><mtd><msubsup><mi mathvariant="bold">r</mi><mn mathvariant="bold">2</mn><mi mathvariant="bold">ʹ</mi></msubsup></mtd><mtd><mo mathvariant="bold">=</mo><msub><mi mathvariant="bold">r</mi><mn mathvariant="bold">2</mn></msub><mo mathvariant="bold">-</mo><msub><mi mathvariant="bold">A</mi><mn mathvariant="bold">2</mn></msub><mo>⁢</mo><msup><mover><mi mathvariant="bold">d</mi><mo>^</mo></mover><mo mathvariant="bold">*</mo></msup><mfenced><mn mathvariant="bold">1</mn></mfenced></mtd></mtr></mtable></mrow></math><img file="EP1533928B1_D0005.tif" /></maths> thereby obtaining <b>r'₁</b> and <b>r'₂</b> ; subtracting effect of another data field <b>d(2)</b> from received data signal based on the following formula: <maths id="math0006"><math display="block"><mrow><mo mathvariant="bold">{</mo><mtable columnalign="left"><mtr><mtd><msubsup><mi mathvariant="bold">r</mi><mn mathvariant="bold">1</mn><mi>ʺ</mi></msubsup></mtd><mtd><mo mathvariant="bold">=</mo><msub><mi mathvariant="bold">r</mi><mn mathvariant="bold">1</mn></msub><mo>+</mo><msub><mi mathvariant="bold">A</mi><mn mathvariant="bold">2</mn></msub><mo>⁢</mo><msup><mover><mi mathvariant="bold">d</mi><mo>^</mo></mover><mo>*</mo></msup><mfenced><mn mathvariant="bold">2</mn></mfenced></mtd></mtr><mtr><mtd><msubsup><mi mathvariant="bold">r</mi><mn mathvariant="bold">2</mn><mi>ʺ</mi></msubsup></mtd><mtd><mo mathvariant="bold">=</mo><msub><mi mathvariant="bold">r</mi><mn mathvariant="bold">2</mn></msub><mo mathvariant="bold">-</mo><msub><mi mathvariant="bold">A</mi><mn mathvariant="bold">1</mn></msub><mo>⁢</mo><mover><mi mathvariant="bold">d</mi><mo>^</mo></mover><mfenced><mn mathvariant="bold">2</mn></mfenced></mtd></mtr></mtable></mrow></math><img file="EP1533928B1_D0006.tif" /></maths> thereby obtaining <b>r"₁</b> and <b>r"₂</b> ; and</li><li>D2. substituting <b>r'₁</b> and <b>r'₂</b> for <b>r₁</b> and <b>r₂</b> in the second equation of said simplified joint detection formula used in step C, and substituting <b>r"₁</b> and <b>r"₂</b> for <b>r₁</b> and <b>r₂</b> in the first equation of said simplified joint detection formula used in step C, calculating said simplified joint detection formula, thereby obtaining iteration results of <b>d̂(1)</b> and <b>d̂(2).</b></li></ul></li></ul>

0013Two diversity signals in step A may be transmitted through two diversity beams of one smart antenna respectively and simultaneously.

0014The method may further comprise the step of predefining number of iteration times to determine execution times from step C to step D and from step D to step C again.

0015The method of the invention is a new space-time encoding and decoding method that is suitable to a frequency selective fading channel. The method takes an independent data field as a processing unit, and the decoding applies an iteration method based on joint detection and interference counteraction. When making joint detection, only effect to diversity signal from multipath interference and multi-user interference is taken into account, and then interference counteraction is implemented based on the joint detection result to remove interference between diversity signals.

0016With the space-time encoding and decoding scheme proposed in the invention, good performance can be reached with small amount of computation complexity in a frequency selective fading channel. The invention provides a simple and effective solution for space-time encoding and decoding on a frequency selective fading channel.

<u>Brief Description of the Drawings</u>

0017<ul id="ul0004" list-style="none"><li><figref idref="f0001">Figure 1</figref> is a schematic diagram illustrating a data structure for time-slot CDMA burst data.</li><li><figref idref="f0001">Figure 2</figref> shows the input and output of a space-time encoder when taking a data field as a processing unit.</li><li><figref idref="f0002">Figure 3</figref> is a flowchart illustrating the simplified decoding procedure of the invention.</li></ul>

<u>Detailed Description of the Invention</u>

0018The invention will be described in more detail hereinafter with reference to the accompanying drawings.

0019In a frequency selective fading channel, the ISI exists between neighbor symbols. Therefore, in a general space-time encoding method, there is not only interference between diversity signals, but also interference between neighbor encoding units and interference between neighbor encoding blocks. If the interference were neglected, the communication performance would be worse; and if the interference were considered, the calculation complexity would be increased greatly.

0020An independent data block is taken as the processing unit in the encoding and decoding method according to the present invention, so the interference between neighbor encoding units and neighbor encoding blocks does not exist, which simplifies the signal detection and decoding processing. This processing method is similar or compatible to the processing methods in which space-time encoding processing is not employed to a great extent. Many processing schemes for counteracting ISI and MAI can be used in the space-time encoding situation with small amount of change.

0021<figref idref="f0001">Figure 1</figref> shows data structure of a time-slot CDMA burst data in the 3GPP-TDD system. It can be seen from <figref idref="f0001">Fig.1</figref> that in parallel K burst data of a multi-user or multichannel, each time-slot includes two data fields (or referred as data blocks): DATA FIELD 1 and DATA FIELD 2, and there is a midamble between these two data fields. So these two data fields are independent and there is no mutual interference between them. The end of a time slot is a guard period (GP). Every data field has N symbols, and every symbol is consisted of Q chips.

0022In this invention, an independent data field is used as a processing unit for encoding, and at the receiving end the decoding also takes a data field as a processing unit. It is simple when this method is used in double diversity (theoretically it can be used in the multiple diversity, but the complexity is increased and therefore it is not practicable). Therefore, the encoding and decoding method according to the present invention can be employed with small amount of change for the signal structure and processing method. Furthermore, since the two data fields are independent, and there is no interference between them, the detection and decoding method can be simplified and improved, and better detection performance can be obtained with a small amount of calculation volume or a small quantity of iterative times.

0023<figref idref="f0001">Figure 2</figref> shows a space-time encoder 20 taking a data field as a processing unit. Input data of the encoder 20, i.e., information data vector d in a time slot, can be represented as <b>d=[d</b>^T(1) <b>d</b>^T(2)<b>]</b>^T<b>;</b> wherein <b>d</b>(1) and <b>d</b>(2) are two independent data fields, and in a 3GPP-TDD system they can be two data fields in one time slot; T denotes a transpose operation. The data which has been encoded by the encoder 20 is <b>[d</b>^T<b>(1)d</b>^T<b>(2)]</b>^T and [-<b>d</b>^*T(2)<b>d</b>^*T(1)]^T , wherein * represents conjugate. The two generated data vectors can be transmitted simultaneously by two different conventional diversity antennas: diversity antenna 1 and diversity antenna 2, or by two different diversity beams of one smart antenna: smart antenna diversity beam 1 and smart antenna diversity beam 2. In other words, the diversity antenna 1 (or diversity beam 1) transmits data vector [<b>d</b>^T(1)<b>d^T</b>(2)]^T and the diversity antenna 2 (or diversity beam2) transmits data vector [-<b>d</b>^*T(2)<b>d</b>^*T(1)]^T simultaneously, thereby the space-time encoding is realized. Specially, the input data [<b>d</b>^T(1)<b>d</b>^T(2)]^T is encoded by conventional space-time diversity orthogonal encoding method, and two output data vectors: [<b>d</b>^T(1)<b>d</b>^T(2)]^T and [-<b>d</b>*^T(2)<b>d</b>^*T(1)]^T are obtained and transmitted via two diversity antennas simultaneously and respectively.

0024In the following, the detection and decoding method for space-time encoding mode which takes a data field as a processing unit will be described, i.e. the detection and decoding scheme corresponding to the designed data field encoding scheme will be described. Here, it is called the original detection and decoding method.

0025Data corresponding to a data field can be expressed as the formula (1): <maths id="math0007" num="(1)"><math display="block"><mi mathvariant="bold">d</mi><mfenced><mi mathvariant="italic">i</mi></mfenced><mo mathvariant="italic">=</mo><msup><mfenced><msubsup><mi>d</mi><mn>1</mn><mfenced><mn>1</mn></mfenced></msubsup><mo>⋯</mo><msubsup><mi>d</mi><mn>1</mn><mfenced><mi>K</mi></mfenced></msubsup><msubsup><mi>d</mi><mn>2</mn><mfenced><mn>1</mn></mfenced></msubsup><mo>⋯</mo><msubsup><mi>d</mi><mn>2</mn><mfenced><mi>K</mi></mfenced></msubsup><mo>⋯</mo><msubsup><mi>d</mi><mi>N</mi><mfenced><mn>1</mn></mfenced></msubsup><mo>⋯</mo><msubsup><mi>d</mi><mi>N</mi><mfenced><mi>K</mi></mfenced></msubsup></mfenced><mi mathvariant="italic">T</mi></msup></math><img file="EP1533928B1_D0007.tif" /></maths> Wherein i is 1, 2; <i>K</i> is the number of terminal users working simultaneously; <i>N</i> is the number of symbols in a user data field.

0026Suppose <b>r</b>₁ and <b>r</b>₂ represent the sample values of two data fields received at the terminal, and there is no interference between <b>r</b>₁ and <b>r</b>₂ even on a frequency selective fading channel. The <b>r</b>₁ and <b>r</b>₂ can be expressed as the formula (2): <maths id="math0008" num="(2)"><math display="block"><mrow><mo mathvariant="bold">{</mo><mtable columnalign="left"><mtr><mtd><msub><mi mathvariant="bold">r</mi><mn mathvariant="bold">1</mn></msub></mtd><mtd><mo mathvariant="bold">=</mo><msub><mi mathvariant="bold">A</mi><mn mathvariant="bold">1</mn></msub><mo>⁢</mo><mi mathvariant="bold">d</mi><mfenced><mn mathvariant="bold">1</mn></mfenced><mo mathvariant="bold">-</mo><msub><mi mathvariant="bold">A</mi><mn mathvariant="bold">2</mn></msub><mo>⁢</mo><msup><mi mathvariant="bold">d</mi><mo mathvariant="bold">*</mo></msup><mfenced><mn mathvariant="bold">2</mn></mfenced><mo mathvariant="bold">+</mo><msub><mi mathvariant="bold">n</mi><mn mathvariant="bold">1</mn></msub></mtd></mtr><mtr><mtd><msub><mi mathvariant="bold">r</mi><mn mathvariant="bold">2</mn></msub></mtd><mtd><mo mathvariant="bold">=</mo><msub><mi mathvariant="bold">A</mi><mn mathvariant="bold">2</mn></msub><mo>⁢</mo><msup><mi mathvariant="bold">d</mi><mo>*</mo></msup><mfenced><mn mathvariant="bold">1</mn></mfenced><mo>+</mo><msub><mi mathvariant="bold">A</mi><mn mathvariant="bold">1</mn></msub><mo>⁢</mo><mi mathvariant="bold">d</mi><mfenced><mn mathvariant="bold">2</mn></mfenced><mo mathvariant="bold">+</mo><msub><mi mathvariant="bold">n</mi><mn mathvariant="bold">2</mn></msub></mtd></mtr></mtable></mrow></math><img file="EP1533928B1_D0008.tif" /></maths>

0027Wherein <b>n</b>_i is the noise vector of the ith data field; <b>A</b>_i is the system matrix of the signal transmission between the ith transmitting antenna and receiving antenna, and the system matrix is determined by the channel pulse response and the user transmission waveform. The formula (2) can be rewritten as formula (3): <maths id="math0009" num="(3)"><math display="block"><mi mathvariant="bold">r</mi><mo mathvariant="bold">=</mo><msub><mi mathvariant="bold">Ad</mi><mi mathvariant="bold">t</mi></msub><mo mathvariant="bold">+</mo><mi mathvariant="bold">n</mi></math><img file="EP1533928B1_D0009.tif" /></maths> Wherein <b>r</b>=[<b>r</b>₁^T,<b>r</b>₂^*T]^T,<b>d</b>_t =[<b>d</b>^T(1)<b>d</b>^*T(2)]^T,<b>n</b>=[<b>n</b>₁^T,<b>n</b>₂^*T]^T, and <maths id="math0010" num="(4)"><math display="block"><mi mathvariant="bold">A</mi><mo mathvariant="bold">=</mo><mfenced open="[" close="]"><mtable><mtr><mtd><msub><mi mathvariant="bold">A</mi><mn mathvariant="bold">1</mn></msub></mtd><mtd><mo mathvariant="bold">-</mo><msub><mi mathvariant="bold">A</mi><mn mathvariant="bold">2</mn></msub></mtd></mtr><mtr><mtd><msubsup><mi mathvariant="bold">A</mi><mn mathvariant="bold">2</mn><mo mathvariant="bold">*</mo></msubsup></mtd><mtd><msubsup><mi mathvariant="bold">A</mi><mn mathvariant="bold">1</mn><mo mathvariant="bold">*</mo></msubsup></mtd></mtr></mtable></mfenced><mspace width="1em" /><mn>.</mn></math><img file="EP1533928B1_D0010.tif" /></maths>

0028From formula (4), it can be obtained that: <maths id="math0011" num="(5)"><math display="block"><msup><mi mathvariant="bold">A</mi><mrow><mo mathvariant="bold">*</mo><mi mathvariant="bold">T</mi></mrow></msup><mo>⁢</mo><mi mathvariant="bold">A</mi><mo mathvariant="bold">=</mo><mfenced open="[" close="]"><mtable><mtr><mtd><msubsup><mi mathvariant="bold">A</mi><mn mathvariant="bold">1</mn><mrow><mo mathvariant="bold">*</mo><mi mathvariant="bold">T</mi></mrow></msubsup><mo>⁢</mo><msub><mi mathvariant="bold">A</mi><mn mathvariant="bold">1</mn></msub><mo mathvariant="bold">+</mo><msup><mfenced><msubsup><mi mathvariant="bold">A</mi><mn mathvariant="bold">2</mn><mrow><mo mathvariant="bold">*</mo><mi mathvariant="bold">T</mi></mrow></msubsup><mo>⁢</mo><msub><mi mathvariant="bold">A</mi><mn mathvariant="bold">2</mn></msub></mfenced><mo mathvariant="bold">*</mo></msup></mtd><mtd><msup><mfenced><msubsup><mi mathvariant="bold">A</mi><mn mathvariant="bold">1</mn><mrow><mo mathvariant="bold">*</mo><mi mathvariant="bold">T</mi></mrow></msubsup><mo>⁢</mo><msub><mi mathvariant="bold">A</mi><mn mathvariant="bold">2</mn></msub></mfenced><mi mathvariant="bold">T</mi></msup><mo mathvariant="bold">-</mo><msubsup><mi mathvariant="bold">A</mi><mn mathvariant="bold">1</mn><mrow><mo mathvariant="bold">*</mo><mi mathvariant="bold">T</mi></mrow></msubsup><mo>⁢</mo><msub><mi mathvariant="bold">A</mi><mn mathvariant="bold">2</mn></msub></mtd></mtr><mtr><mtd><msup><mfenced><mtable><mtr><mtd><msup><mfenced><msubsup><mi mathvariant="bold">A</mi><mn mathvariant="bold">1</mn><mrow><mo mathvariant="bold">*</mo><mi mathvariant="bold">T</mi></mrow></msubsup><mo>⁢</mo><msub><mi mathvariant="bold">A</mi><mn mathvariant="bold">2</mn></msub></mfenced><mi mathvariant="bold">T</mi></msup><mo mathvariant="bold">-</mo><msubsup><mi mathvariant="bold">A</mi><mn mathvariant="bold">1</mn><mrow><mo mathvariant="bold">*</mo><mi mathvariant="bold">T</mi></mrow></msubsup><mo>⁢</mo><msub><mi mathvariant="bold">A</mi><mn mathvariant="bold">2</mn></msub></mtd></mtr></mtable></mfenced><mrow><mo mathvariant="bold">*</mo><mi mathvariant="bold">T</mi></mrow></msup></mtd><mtd><msup><mfenced><msubsup><mi mathvariant="bold">A</mi><mn mathvariant="bold">1</mn><mrow><mo mathvariant="bold">*</mo><mi mathvariant="bold">T</mi></mrow></msubsup><mo>⁢</mo><msub><mi mathvariant="bold">A</mi><mn mathvariant="bold">1</mn></msub><mo mathvariant="bold">+</mo><msup><mfenced><msubsup><mi mathvariant="bold">A</mi><mn mathvariant="bold">2</mn><mrow><mo mathvariant="bold">*</mo><mi mathvariant="bold">T</mi></mrow></msubsup><mo>⁢</mo><msub><mi mathvariant="bold">A</mi><mn mathvariant="bold">2</mn></msub></mfenced><mo mathvariant="bold">*</mo></msup></mfenced><mo mathvariant="bold">*</mo></msup></mtd></mtr></mtable></mfenced></math><img file="EP1533928B1_D0011.tif" /></maths> The matrix <b>A</b> is not an orthogonal matrix under the condition of frequency selective fading channel, so the performance will become worse if a match filter is applied. In order to obtain better performance, the optimized linear joint detection scheme is provided. The continuously estimated value <b>d̂</b>_t of the receiving data <b>d̂</b>_t is: <maths id="math0012" num="(6)"><math display="block"><msub><mover><mi mathvariant="bold">d</mi><mo mathvariant="bold">^</mo></mover><mi mathvariant="bold">t</mi></msub><mo mathvariant="bold">=</mo><msup><mfenced><mi mathvariant="bold">B</mi></mfenced><mrow><mo mathvariant="bold">-</mo><mn mathvariant="bold">1</mn></mrow></msup><mo>⁢</mo><msup><mi mathvariant="bold">A</mi><mrow><mo mathvariant="bold">*</mo><mi mathvariant="bold">T</mi></mrow></msup><mo>⁢</mo><mi mathvariant="bold">r</mi></math><img file="EP1533928B1_D0012.tif" /></maths>

0029Wherein (<b>B</b>)^-1 has the function of interference suppression (implement inverse operation for B); <b>A</b>^<b>*</b>T<b>r</b> is the result of orthogonal match. The matrix <b>B</b> is shown in formula (7): <maths id="math0013" num="(7)"><math display="block"><mi mathvariant="bold">B</mi><mo>=</mo><mrow><mo>{</mo><mtable><mtr><mtd><mi mathvariant="bold">I</mi></mtd><mtd><mi mathvariant="italic">MF</mi></mtd></mtr><mtr><mtd><msup><mi mathvariant="bold">A</mi><mrow><mo mathvariant="bold">*</mo><mi mathvariant="bold">T</mi></mrow></msup><mo>⁢</mo><mi mathvariant="bold">A</mi></mtd><mtd><mi mathvariant="italic">ZF</mi><mo mathvariant="italic">-</mo><mi mathvariant="italic">BLE</mi></mtd></mtr><mtr><mtd><msup><mi mathvariant="bold">A</mi><mrow><mo mathvariant="bold">*</mo><mi mathvariant="bold">T</mi></mrow></msup><mo>⁢</mo><mi mathvariant="bold">A</mi><mo mathvariant="bold">+</mo><msup><mi>σ</mi><mn mathvariant="bold">2</mn></msup><mo>⁢</mo><mi mathvariant="bold">I</mi></mtd><mtd><mi mathvariant="italic">MMSE</mi><mo mathvariant="italic">-</mo><mi mathvariant="italic">BLE</mi></mtd></mtr></mtable></mrow></math><img file="EP1533928B1_D0013.tif" /></maths> wherein σ² is the noise power, and I is an identity matrix. The formula (7) shows three solutions among which MF represents matched filter scheme, ZF-BLE represents zero-forcing block equalization scheme and MMSE-BLE represents the minimum mean-square-error block equalization scheme.

0030The dimension of matrix <b>B</b> is 2NK × 2NK, so the calculation complexity of formula (6) is far greater than that under the situation of no space-time encoding. Therefore, the object of simplifying processing and reducing calculation volume proposed in the invention cannot be realized with the original joint detection algorithm.

0031Based on the above-mentioned design, the invention provides a simplified decoding procedure having a characteristic of a small amount of calculation volume.

0032The simplified algorithm is an iterative algorithm. It is divided into two steps. In the first step, only multi-path interference and multi-user interference to every diversity signal is taken into account, and the upper right block of the formula (5) and the lower left block of the formula (5) are set to null, then the equation (6) is calculated. In the second step, interference is counteracted with the result of the first step to remove interference between diversity signals. The procedure can be done by multiple iterations.

0033<figref idref="f0002">Figure 3</figref> is a flowchart of the simplified decoding procedure of the invention.

0034Step 31, the terminal receives data.

0035Step 32, mutual interference between diversity signals caused by their non-orthogonality is neglected, upper right block and lower left block of formula (5) are set to be null (they are null under the condition of a flat fading channel).

0036Step 33, joint detection is performed, wherein only multipath interference and multi-user interference to every diversity signal is taken account, in this case formula (6) is simplified to formula (8) shown in the following: <maths id="math0014" num="(8)"><math display="block"><mrow><mo mathvariant="bold">{</mo><mtable columnalign="left"><mtr><mtd><mover><mi mathvariant="bold">d</mi><mo mathvariant="bold">^</mo></mover><mfenced><mn mathvariant="bold">1</mn></mfenced></mtd><mtd><mo mathvariant="bold">=</mo><msubsup><mi mathvariant="bold">B</mi><mi mathvariant="bold">S</mi><mrow><mo mathvariant="bold">-</mo><mn mathvariant="bold">1</mn></mrow></msubsup><mo>⁢</mo><mfenced><msubsup><mi mathvariant="bold">A</mi><mn mathvariant="bold">1</mn><mrow><mo mathvariant="bold">*</mo><mi mathvariant="bold">T</mi></mrow></msubsup><mo>⁢</mo><msub><mi mathvariant="bold">r</mi><mn mathvariant="bold">1</mn></msub><mo mathvariant="bold">+</mo><msup><mfenced><msubsup><mi mathvariant="bold">A</mi><mn mathvariant="bold">2</mn><mrow><mo mathvariant="bold">*</mo><mi mathvariant="bold">T</mi></mrow></msubsup><mo>⁢</mo><msub><mi mathvariant="bold">r</mi><mn mathvariant="bold">2</mn></msub></mfenced><mo mathvariant="bold">*</mo></msup></mfenced></mtd></mtr><mtr><mtd><mover><mi mathvariant="bold">d</mi><mo mathvariant="bold">^</mo></mover><mfenced><mn mathvariant="bold">2</mn></mfenced></mtd><mtd><mo mathvariant="bold">=</mo><msubsup><mi mathvariant="bold">B</mi><mi mathvariant="bold">S</mi><mrow><mo mathvariant="bold">-</mo><mn mathvariant="bold">1</mn></mrow></msubsup><mo>⁢</mo><mfenced><msubsup><mi mathvariant="bold">A</mi><mn mathvariant="bold">1</mn><mrow><mo mathvariant="bold">*</mo><mi mathvariant="bold">T</mi></mrow></msubsup><mo>⁢</mo><msub><mi mathvariant="bold">r</mi><mn mathvariant="bold">2</mn></msub><mo mathvariant="bold">-</mo><msup><mfenced><msubsup><mi mathvariant="bold">A</mi><mn mathvariant="bold">2</mn><mrow><mo mathvariant="bold">*</mo><mi mathvariant="bold">T</mi></mrow></msubsup><mo>⁢</mo><msub><mi mathvariant="bold">r</mi><mn mathvariant="bold">1</mn></msub></mfenced><mo mathvariant="bold">*</mo></msup></mfenced></mtd></mtr></mtable></mrow></math><img file="EP1533928B1_D0014.tif" /></maths>

0037Formula (8) gives the result of simplified joint detection. The matrix <b>B_S</b> in formula (8) is given by formula (9): <maths id="math0015" num="(9)"><math display="block"><msub><mi mathvariant="bold">B</mi><mi mathvariant="bold">S</mi></msub><mo>=</mo><mrow><mo>{</mo><mtable><mtr><mtd><mi mathvariant="bold">I</mi></mtd><mtd><mi mathvariant="italic">MF</mi></mtd></mtr><mtr><mtd><msubsup><mi mathvariant="bold">A</mi><mn mathvariant="bold">1</mn><mrow><mo mathvariant="bold">*</mo><mi mathvariant="bold">T</mi></mrow></msubsup><mo>⁢</mo><msub><mi mathvariant="bold">A</mi><mn mathvariant="bold">1</mn></msub><mo mathvariant="bold">+</mo><msup><mfenced><msubsup><mi mathvariant="bold">A</mi><mn mathvariant="bold">2</mn><mrow><mo mathvariant="bold">*</mo><mi mathvariant="bold">T</mi></mrow></msubsup><mo>⁢</mo><msub><mi mathvariant="bold">A</mi><mn mathvariant="bold">2</mn></msub></mfenced><mo mathvariant="bold">*</mo></msup></mtd><mtd><mi mathvariant="italic">ZF</mi><mo mathvariant="italic">-</mo><mi mathvariant="italic">BLE</mi></mtd></mtr><mtr><mtd><msubsup><mi mathvariant="bold">A</mi><mn mathvariant="bold">1</mn><mrow><mo mathvariant="bold">*</mo><mi mathvariant="bold">T</mi></mrow></msubsup><mo>⁢</mo><msub><mi mathvariant="bold">A</mi><mn mathvariant="bold">1</mn></msub><mo mathvariant="bold">+</mo><msup><mfenced><msubsup><mi mathvariant="bold">A</mi><mn mathvariant="bold">2</mn><mrow><mo mathvariant="bold">*</mo><mi mathvariant="bold">T</mi></mrow></msubsup><mo>⁢</mo><msub><mi mathvariant="bold">A</mi><mn mathvariant="bold">2</mn></msub></mfenced><mo mathvariant="bold">*</mo></msup><mo mathvariant="bold">+</mo><msup><mi>σ</mi><mn mathvariant="bold">2</mn></msup><mo>⁢</mo><mi mathvariant="bold">I</mi></mtd><mtd><mi mathvariant="italic">MMSE</mi><mo mathvariant="italic">-</mo><mi mathvariant="italic">BLE</mi></mtd></mtr></mtable></mrow></math><img file="EP1533928B1_D0015.tif" /></maths>

0038Wherein σ² is the noise power, and <b>I</b> is the identity matrix. Calculating matrix B_s also has three ways: match filter MF, zero-forcing block equalization ZF-BLE and minimum mean-square-error block equalization MMSE - BLE.

0039The dimension of matrix <b>B_S</b> is NK × NK, so the calculation complexity of formula (8) is far less than that of formula (6) and is similar to that under the situation of no space-time encoding.

0040Step 34, whether the number of iteration times is equal to or greater than the predefined number <i>M</i> is judged. Usually <i>M</i> is set at 1 or 2.

0041Step 35, when the number of iteration times is equal to or greater than the predefined number <i>M</i>, the calculation result of formula (8) is outputted as the decoding result directly.

0042Step 36, when the number of iteration times is less than the predefined number <i>M</i>, the calculation result of formula (8) is used to implement interference counteraction so as to remove interference between diversity signals, that is, an interference counteraction method is used to counteract the remain interference between two diversity signals. This is done as follows: <ul id="ul0005" list-style="none" compact="compact"><li>Affect of <b>d(1)</b> is subtracted from the received data signal to obtain a 'clean' signal as shown in formula (10): <maths id="math0016" num="(10)"><math display="block"><mrow><mo mathvariant="bold">{</mo><mtable columnalign="left"><mtr><mtd><msubsup><mi mathvariant="bold">r</mi><mn mathvariant="bold">1</mn><mi mathvariant="bold">ʹ</mi></msubsup></mtd><mtd><mo mathvariant="bold">=</mo><msub><mi mathvariant="bold">r</mi><mn mathvariant="bold">1</mn></msub><mo mathvariant="bold">-</mo><msub><mi mathvariant="bold">A</mi><mn mathvariant="bold">1</mn></msub><mo>⁢</mo><mover><mi mathvariant="bold">d</mi><mo>^</mo></mover><mfenced><mn mathvariant="bold">1</mn></mfenced></mtd></mtr><mtr><mtd><msubsup><mi mathvariant="bold">r</mi><mn mathvariant="bold">2</mn><mi mathvariant="bold">ʹ</mi></msubsup></mtd><mtd><mo mathvariant="bold">=</mo><msub><mi mathvariant="bold">r</mi><mn mathvariant="bold">2</mn></msub><mo mathvariant="bold">-</mo><msub><mi mathvariant="bold">A</mi><mn mathvariant="bold">2</mn></msub><mo>⁢</mo><msup><mover><mi mathvariant="bold">d</mi><mo>^</mo></mover><mo mathvariant="bold">*</mo></msup><mfenced><mn mathvariant="bold">1</mn></mfenced></mtd></mtr></mtable></mrow></math><img file="EP1533928B1_D0016.tif" /></maths> and affect of <b>d(2)</b> is subtracted from the received signal to obtain another 'clean' signal as shown in formula (11): <maths id="math0017" num="(11)"><math display="block"><mrow><mo mathvariant="bold">{</mo><mtable columnalign="left"><mtr><mtd><msubsup><mi mathvariant="bold">r</mi><mn mathvariant="bold">1</mn><mi>ʺ</mi></msubsup></mtd><mtd><mo mathvariant="bold">=</mo><msub><mi mathvariant="bold">r</mi><mn mathvariant="bold">1</mn></msub><mo>+</mo><msub><mi mathvariant="bold">A</mi><mn mathvariant="bold">2</mn></msub><mo>⁢</mo><msup><mover><mi mathvariant="bold">d</mi><mo>^</mo></mover><mo>*</mo></msup><mfenced><mn mathvariant="bold">2</mn></mfenced></mtd></mtr><mtr><mtd><msubsup><mi mathvariant="bold">r</mi><mn mathvariant="bold">2</mn><mi>ʺ</mi></msubsup></mtd><mtd><mo mathvariant="bold">=</mo><msub><mi mathvariant="bold">r</mi><mn mathvariant="bold">2</mn></msub><mo mathvariant="bold">-</mo><msub><mi mathvariant="bold">A</mi><mn mathvariant="bold">1</mn></msub><mo>⁢</mo><mover><mi mathvariant="bold">d</mi><mo>^</mo></mover><mfenced><mn mathvariant="bold">2</mn></mfenced></mtd></mtr></mtable></mrow></math><img file="EP1533928B1_D0017.tif" /></maths></li></ul>

0043Iteration operation is implemented. In detail, <b>r</b>₁ and <b>r</b>₂ in the second formula of formula (8) are respectively substituted with the result of the formula (10) <b>r</b>₁' and <b>r</b>₂'; <b>r</b>₁ and <b>r</b>₂ in the first formula of formula (8) are respectively substituted with the result of the formula (11) <b>r</b>₁" and <b>r</b>₂"; and then formula (8) is recalculated to obtain a decoding result after one iteration operation.

0044A second iteration operation can be done by taking last result <b>d̂(1)</b> and <b>d̂(2)</b> of formula (8) in formula (10) and formula (11). In other words, the obtained results of <b>r</b>₁', <b>r</b>₂', <b>r</b>₁" and <b>r</b>₂" are substituted to formula (8) again, and then formula (8) is recalculated to obtain a decoding result after the second iteration operation.

0045The iteration procedure can be done for <i>M</i> times, and this is the simplified joint detection algorithm.

0046Experience shows that usually only several iteration times, such as one to two times, are necessary to reach performance of the original joint detection algorithm. Therefore, the simplified joint detection algorithm according to this invention can obtain better performance with small amount of calculation complexity.

0047The invention proposes a space-time encoding and decoding method that takes independent data field as a processing unit for encoding input, and this is different with the prior method that takes data symbols or a string of symbols as a processing unit. With the proposed method, a simplified decoding method is designed to obtain a better performance with small amount of calculation complexity.

0048The invention provides a simple and effective solution for space-time encoding and decoding on a frequency selective fading channel.