1 Star2 Stars3 Stars4 Stars5 Stars (5 votes, average: 3.8 out of 5)
Loading ... Loading ...
Print Print

Six equalizers for V-BLAST

by Krishna Sankar on April 21, 2009

In the past, we had discussed several posts on two transmit two receive MIMO communication, where the transmission was based on V-BLAST. The details about V-BLAST can be read from the landmark paper V-BLAST: An architeture for realizing very high data rates over the rich scattering wireless channel – P. W. Wolniansky, G. J. Foschini, G. D. Golden, R. A. Valenzuela. We will assume that the channel is a flat fading Rayleigh multipath channel and the modulation is BPSK.

Figure: 2 transmit 2 receive MIMO channel

V-BLAST transmission for 2×2 MIMO channel

In a 2×2 MIMO channel, probable usage of the available 2 transmit antennas can be as follows:

1. Consider that we have a transmission sequence, for example

2. In normal transmission, we will be sending in the first time slot, in the second time slot, and so on.

3. However, as we now have 2 transmit antennas, we may group the symbols into groups of two. In the first time slot, send and from the first and second antenna. In second time slot, send and from the first and second antenna, send and in the third time slot and so on.

4. Notice that as we are grouping two symbols and sending them in one time slot, we need only time slots to complete the transmission – data rate is doubled ! :)

5. This forms the simple explanation of a probable MIMO transmission scheme with 2 transmit antennas and 2 receive antennas.

Other Assumptions

1. The channel is flat fading – In simple terms, it means that the multipath channel has only one tap. So, the convolution operation reduces to a simple multiplication. For a more rigorous discussion on flat fading and frequency selective fading, may I urge you to review Chapter 15.3 Signal Time-Spreading from [DIGITAL COMMUNICATIONS: SKLAR]

2. The channel experience by each transmit antenna is independent from the channel experienced by other transmit antennas.

3. For the transmit antenna to receive antenna, each transmitted symbol gets multiplied by a randomly varying complex number . As the channel under consideration is a Rayleigh channel, the real and imaginary parts of are Gaussian distributed having mean and variance .

4. The channel experienced between each transmit to the receive antenna is independent and randomly varying in time.

5. On the receive antenna, the noise has the Gaussian probability density function with

with and .

7. The channel is known at the receiver.

System Model

The received signal on the first receive antenna is,

.

The received signal on the second receive antenna is,

.

where

, are the received symbol on the first and second antenna respectively,

is the channel from transmit antenna to receive antenna,

is the channel from transmit antenna to receive antenna,

is the channel from transmit antenna to receive antenna,

is the channel from transmit antenna to receive antenna,

, are the transmitted symbols and

is the noise on receive antennas.

We assume that the receiver knows , , and . The receiver also knows and . The unknown s are and . Two equations and two unknowns. Can we solve it? Answer is YES. :)

For convenience, the above equation can be represented in matrix notation as follows:

.

Equivalently,

Receiver structures

MIMO with Zero Forcing Equalization

The zero forcing approach tries to find a matrix which satisfies . The Zero Forcing (ZF) linear detector for meeting this constraint is given by,

.

MIMO with MMSE Equalization

The Minimum Mean Square Error (MMSE) approach tries to find a coefficient which minimizes the criterion,

.

Solving,

.

Zero Forcing Equalization with Successive Interference Cancellation

Using the Zero Forcing (ZF) equalization approach described above, the receiver can obtain an estimate of the two transmitted symbols , , i.e.

.

Take one of the estimated symbols (for example ) and subtract its effect from the received vector and , i.e.

.

Expressing in matrix notation,

,

The above equation is same as equation obtained for receive diversity case. Optimal way of combining the information from multiple copies of the received symbols in receive diversity case is to apply Maximal Ratio Combining (MRC).

The equalized symbol is,

.

This forms the simple explanation for Zero Forcing Equalizer with Successive Interference Cancellation (ZF-SIC) approach.

Zero Forcing Equalization with Optimally ordered Successive Interference Cancellation

In classical Successive Interference Cancellation, the receiver arbitrarily takes one of the estimated symbols, and subtract its effect from the received symbol and . However, we can have more intelligence in choosing whether we should subtract the effect of first or first. To make that decision, let us find out the transmit symbol (after multiplication with the channel) which came at higher power at the receiver. The received power at the both the antennas corresponding to the transmitted symbol is,

.

The received power at the both the antennas corresponding to the transmitted symbol is,

.

If then the receiver decides to remove the effect of from the received vector and and then re-estimate . Else if the receiver decides to subtract effect of from the received vector and , and then re-estimate .

MMSE equalization with optimaly ordered Successive Interference Cancellation

Using the Minimum Mean Square Error (MMSE) equalization, the receiver can obtain an estimate of the two transmitted symbols , , i.e.

.

If then the receiver decides to remove the effect of from the received vector and . Else if the receiver decides to subtract effect of from the received vector and , and then re-estimate .

Once the effect of either or is removed, the new channel becomes a one transmit antenna, 2 receive antenna case and the symbol on the other spatial dimension can be optimally equalized by Maximal Ratio Combining (MRC).

MIMO with ML equalization

The Maximum Likelihood receiver tries to find which minimizes,

Since the modulation is BPSK, the possible values of is +1 or -1 Similarly also take values +1 or -1. So, to find the Maximum Likelihood solution, we need to find the minimum from the all four combinations of and .

The estimate of the transmit symbol is chosen based on the minimum value from the above four values i.e

if the minimum is ,

if the minimum is ,

if the minimum is and

if the minimum is .

Simulation Results

The plot embedded below captures BER for 2 transmit 2 receive MIMO V-BLAST ttansmission/reception for BPSK modulation in a flat fading independent Rayleigh channel, for the different equalizer structures discussed above.

Figure: BER plot for 2 transmit 2 receive MIMO channel for BPSK modulation

Observations

1. The BER curve with ZF equalization for 2×2 MIMO channel is identical to BER plot for 1 transmit 1 receive system

2. Ordered variant of successive interference cancellation shows better performance than the simple successive interference cancellation

3. MMSE equalization with ordered successive interference cancellation provides performance which is slightly poorer than ML.

4. With ML equalization, we come close to the performance of 1 transmit 2 receive MRC case. We gain both throughput gain and diversity gain.

Reference

[DIG-COMM-BARRY-LEE-MESSERSCHMITT] Digital Communication: Third Edition, by John R. Barry, Edward A. Lee, David G. Messerschmitt

V-BLAST: An architeture for realizing very high data rates over the rich scattering wireless channel – P. W. Wolniansky, G. J. Foschini, G. D. Golden, R. A. Valenzuela.

Related posts

  1. MIMO with MMSE SIC and optimal ordering
  2. MIMO with ZF SIC and optimal ordering
  3. MIMO with Zero Forcing Successive Interference Cancellation equalizer
  4. MIMO with ML equalization
  5. MIMO with MMSE equalizer

D id you like this article? Make sure that you do not miss a new article by subscribing to RSS feed OR subscribing to e-mail newsletter. Note: Subscribing via e-mail entitles you to download the free e-Book on BER of BPSK/QPSK/16QAM/16PSK in AWGN.

{ 45 comments… read them below or add one }

1 D1EG0 April 21, 2009 at 9:38 pm

Can you put the matlab code of this graph for know how can I generate all the matrix.

Thanks

Reply

2 Krishna Pillai April 30, 2009 at 5:04 am

@D1EG0: Matlab code for the curve for each receiver structure is present in the individual articles (links are present in the post). I just made a top level simulation to generate all the curves in one go.

Reply

3 benjamin.frank May 6, 2009 at 8:05 am

after evaluating the bit error rate, can you consider the bit loading in ofdm scheme and in the receiver using V-blast or D-blast? thanks

Reply

4 Krishna Pillai May 12, 2009 at 5:10 am

@benjamin.frank: I will surely post on MIMO OFDM using V-BLAST. D-BLAST is bit more sophisticated and I need to stuby more.

Reply

5 saeed April 22, 2009 at 8:26 am

thank you so much for this file .and i hope to give me some code about quality of service and call admission control because i’am research in this topic.
thanks

Reply

6 Krishna Pillai April 30, 2009 at 5:09 am

@saeed: I have not tried modeling higher layer protocols like QoS, admission control etc. Good luck in your research.

Reply

7 ilhami gozcu April 23, 2009 at 8:03 pm

Hi,

Thanks for the useful post. Is it posible to add the Matlab code of the simulation?

Reply

8 Krishna Pillai April 30, 2009 at 5:10 am

@ilhami gozcu: Matlab code for the curve for each receiver structure is present in the individual articles (links are present in the post). I just made a top level simulation to generate all the curves in one go.

Reply

9 VC April 27, 2009 at 10:48 pm

Hi, can you please send me the source codes for the simulations. Thanks a lot.

Reply

10 Krishna Pillai April 30, 2009 at 5:30 am

@VC:Matlab code for the curve for each receiver structure is present in the individual articles (links are present in the post). I just made a top level simulation to generate all the curves in one go.

Reply

11 Alvina April 29, 2009 at 2:20 pm

@krishna
can you please give me the source code for the simulation. Thanks.

Reply

12 Krishna Pillai April 30, 2009 at 5:43 am

@Alvina: Matlab code for the curve for each receiver structure is present in the individual articles (links are present in the post). I just made a top level simulation to generate all the curves in one go.

Reply

13 benjamin.frank May 6, 2009 at 8:14 am

Krishna Pillai,
I am right now calculating the ber versus snr in qpsk-ofdm environement,but the curve is not what we expected, it is different from the theory,maybe the code has some problems, how can i contact you?my mail address:benjamin.frank@163.com

Reply

14 Krishna Pillai May 12, 2009 at 5:12 am

@benjamin.frank: Well, due to time constraints, I prefer not to debug the code. May I point you to a post on Symbol Error Rate with QPSK modulation in AWGN.
http://www.dsplog.com/2007/11/06/symbol-error-rate-for-4-qam/

Good luck.

Reply

15 Alvina May 7, 2009 at 4:55 pm

@benjamin u can find contact of krishna at

http://www.dsplog.com/contact-us/

Reply

16 Mark June 6, 2009 at 12:10 pm

can you please give me the source code for the simulation. Thanks.

Reply

17 Krishna Pillai June 7, 2009 at 2:37 pm

@Mark: The source code for each equalizer structure is provided as part of the individual articles. You may click the hyper link to see the articles in detail and obtain the source code.

Reply

18 sohbet July 2, 2009 at 6:11 pm

hi good thank you blog

Reply

19 Habib July 10, 2009 at 5:44 pm

Hi there,
I am M.Sc researcher in IIUM. I am working on 2X2 MIMO channel sounder. I need an urgent help for a simulink model of MIMO transmitter. can you help me, how to call M-file in simulink. I have created an M-function or code generator. I need to ad this as a source of transmitter in simulink.

Thank you

Reply

20 Krishna Pillai July 15, 2009 at 4:59 am

@Habib: Sorry, I have not worked with Simulink.

Reply

21 Ashok July 24, 2009 at 10:46 am

Hi Habib,
What kind of MIMO you need?
STBC MIMO demo is there in simulink if you have communication tool box. Check the WiMax Model.
For using matlab code use ‘Embedded MATLAB block’ in simulink.

Reply

22 raj September 19, 2009 at 11:54 am

hi sir,
I extremely thanks you to share this tremendous work in field of MIMO spatial multiplexing,
I am working in the same field , so i request you to share some details if you don’t mind
I have certain query about the V-BLAST programme,

1.traditionally BLAST scheme uses three different types of coding, V-BLAST, D-BLAST, H-BLAST, so this programme uses coding?

2. and one most important question is that What is the most recent possible research scope in this MIMO spatial multiplexing field?

Ones again thanks for showing interest in this field.

Reply

23 Krishna Sankar September 22, 2009 at 5:43 am

@raj: My replies:
1/ I used V-Blast where all transmit antennas transmit simultaneously.
2/ One topic which interests me is to simplify the equalizer structure and reach ML like performance (Sphere decoding etc)

Reply

24 raj September 23, 2009 at 11:07 am

Thank you for reply, but my question is that

In V-BLAST , is there any kind of coding (e.g. convolution) is used? b’cos in programs ,I am not found any kind of coding at transmitter, if I want to use coding at transmitter than what changes require at receiver? particularly in case of V-BLAST & D-BLAST.

In the programme of ZF-SIC with optimal ordering , you have assign channel with higher SNR to second stream (force-fully), but if I want to do SIC on the basis of original SNR values of symbol than would I get desired result?

thx

Reply

25 Krishna Sankar October 1, 2009 at 5:15 am

@raj: My replies
1/ No, I have not assumed presence of coding in transmitter. However, even if coding is present, it will just be that the output of the equalizer will go to the decoder and the decoder will generate the raw bits.
2/ I believe you meant optimal ordered successive interference cancellation. To do that select the stream having higher SNR. I have provided a post on Zero Forcing with optimal ordered successive interference cancellation.
http://www.dsplog.com/2008/11/29/mimo-zf-sic-optimal-order/

Hope this helps.

Reply

26 William October 6, 2009 at 12:29 pm

Hello Krishna,

My simulations show a little different results. In my results, the 2×2 BPSK ZF BER lacks the 1×1 BPSK BER by 3 dB. I believe this to be a sane result, because I have scaled the transmission power by 1/sqrt(2) to keep the total transmission power emitted by all antennas constant for a fare comparison between different schemes. Have you done this in your simulations? Similarly, my 2×2 ML BPSK transmission lacks the MRC by 5 dBs.

Reply

27 Krishna Sankar October 8, 2009 at 5:34 am

@William: If you look @ the simulation in post http://www.dsplog.com/2008/10/24/mimo-zero-forcing/, it can be seen that I have scaled the transmit voltage (not power) by 1/sqrt(2). Scaling the transmit voltage by 1/sqrt(2) ensures that power on each arm is 1/2 and totals to unity.

Further, from the Section 3.3 of book Fundamentals of Wireless Communication, David Tse, Pramod Viswanath, its mentioned that “BER for 2×2 MIMO channel in Rayleigh fading with Zero Forcing equalization is same as the BER derived for a 1×1 channel in Rayleigh fading.”

Hope this helps.

Reply

28 raj October 8, 2009 at 11:15 am

Thank you sir,
I have get the result for SIC using optimal ordering, but if I use V-BLAST with coding than the coding is perform before demultiplexing of data so data bits(symbols) get spreded spatially in channel and we have to decode them jointly at receiver which is not happened in H-BLAST (coding is done after demultiplexing) so for V-BLAST,

1. How can the bits are decoded at the receiver? Should it require some change in programme?

2. How many complexity operations are required in ML(maximum likelyhood), ZF, MMSE, ZF-SIC, MMSE-SIC, ZF-OSIC, MMSE-OSIC?

sir very very thanks for showing interest in my queries, bye.

Reply

29 Krishna Sankar October 12, 2009 at 5:20 am

@raj:
1/ In V-BLAST, there is interference between parallel transmissions. At the receiver, we try to null the interference and recover the required data symbol.
2/ I have not done precise calculations on the complexity difference between these algorithms.

Reply

30 Egerue Nnamdi October 8, 2009 at 2:58 pm

Hi all
Pls i am currently an MSC student in the UK and working on (OFDM) as my Dissertation. I am to investigate how to use Matlab to simulate the following channel impairment
1 Multipath Delay
2 Peak power Clipping
3 Channel Noise
4 Frequency selective fading
5 Intersymbol interference.
I am surpose to compare the performance of (OFDM) and BPSK, QPSK etc using the above channel impairments. currently the Matlab codes i have does not plot any thing. Pls i will appreciate if someone can help me out.

NNAMDI

Reply

31 Krishna Sankar October 12, 2009 at 5:27 am

@Egerue: Hope the posts under http://www.dsplog.com/tag/ofdm address some of your queries. Good luck.

Reply

32 raj October 26, 2009 at 4:06 pm

hi sir,
i want check the impact of MIMO (MMSE-OSIC) system on the higher-level protocol, wether it will be better idea for research?

And is there any other technique for MMSE detector, rathar than the mathamatical equation (H^H*H + I/snr)^-1*H^H), where can i find it?

Reply

33 Krishna Sankar November 7, 2009 at 7:10 am

@raj:
1/ What is higher level protocol?
2/ In general the equation for MMSE is as the one listed above. You may find the derivation in typical text books.

Reply

34 raj November 12, 2009 at 10:51 am

sir i want to check the of impact of MIMO (MMSE-OSIC) system on the transport control protocol (TCP) protocol.
and is there any different between MMSE detection and MMSE equalization, if yes than what is that?

Reply

35 Krishna Sankar November 13, 2009 at 5:39 am

@raj: MMSE-OSIC on TCP protocol, eh? How do you apply MMSE in the TCP framework. Can you please point me to some references.

Reply

36 raj November 18, 2009 at 12:52 pm

Sir i will send you some references regarding to that, please send me your e-mail address.

thanks

37 raj November 18, 2009 at 12:54 pm

Is there any different between MMSE detection and MMSE equalization, if yes than what is that?

38 aydar October 27, 2009 at 3:17 pm

wow, i didnt know that Reinaldo Valenzuela is one of the authors..

Reply

39 Krishna Sankar November 8, 2009 at 7:37 am

@aydar: :)

Reply

40 raj November 25, 2009 at 10:34 am

Sir i will send you some references regarding to that, please send me your e-mail address.

Is there any different between MMSE detection and MMSE equalization, if yes than what is that?

thanks

Reply

41 Krishna Sankar December 7, 2009 at 4:31 am

@raj: Well, am not sure. It depends on the problem which you are modeling. The MMSE algorithm can be used to ‘detect the arrival of the packet’ which is quite different from equalizing the packet which is corrupted by noise and interference.

Reply

42 Krishna Sankar December 6, 2009 at 3:51 pm

@raj: Thanks. You can find my email listed in http://www.dsplog.com/contact-us/

Reply

43 Krishna Sankar December 6, 2009 at 3:52 pm

@raj: Well, am not sure. It depends on the problem which you are modeling. The MMSE algorithm can be used to ‘detect the arrival of the packet’ which is quite different from equalizing the packet which is corrupted by noise and interference.

Reply

44 raj December 16, 2009 at 10:30 pm

hi sir,
thanks for reply,
but sir equalization and detection , both are same process and they are different from estimation, am i right?

Is there any difference between SVD (singular value decomposition), SQRD (sorted QR decomposition), and decorrelating process?, b’cos ultimately they all belongs to set diagonal components of any matrix.

What is the different between Sphere decoding (SD) and Lattice reduction?

i am curious to wait for your replay,
once again thanks for your great support.

Reply

45 Krishna Sankar December 23, 2009 at 5:27 am

@raj: my replies
a/ both equalization and detection can be different
b/ i think both svd and sqrd are different. have not studied much on this topic
c/ have not studied in detail sphere decoding and lattice reduction. hence unable to comment

Reply

Leave a Comment

{ 1 trackback }

Previous post:

Next post: