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

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.*

{ 65 comments… read them below or add one }

Hey Krishna !

I am doing VBLAST detection… the program I am using is an OSIC detector function that implements VBLAST algorithm using MMSE (SINR based ordering), ZF (SNR based ordering) and column-norm based ordering. The function is as follows:

function [X_hat]=OSIC_detector(y,H,sigma2,NT,OSIC_type)

% Input

% y : Received signal, NRx1

% H : Channel matrix, NRxNT

% sigma2: Noise variance

% NT : Number of Tx antennas

% OSIC_type=1/2/3 for Post_detection_SINR/Column_max/Post_detection_SNR

% Output

% X_hat : Estimated signal, NTx1

if OSIC_type==1 % Post_detection_SINR

Order=[]; % Detection order

index_array=[1:NT]; % yet to be detected signal index

% V-BLAST

for stage=1:NT

Wmmse=inv(H’*H+sigma2*eye(NT+1-stage))*H’; % MMSE filter (Eq.(11.7))

WmmseH=Wmmse*H;

SINR=[];

for i=1:NT-(stage-1)

tmp= norm(WmmseH(i,[1:i-1 i+1:NT-(stage-1)]))^2+ sigma2*norm(Wmmse(i,:))^2;

SINR(i) = abs(Wmmse(i,i))^2/tmp; % SINR calculation Eq.(11.15)

end

[val_max,index_tmp] = max(SINR); % Ordering using SINR

Order = [Order index_array(index_tmp)];

index_array = index_array([1:index_tmp-1 index_tmp+1:end]);

x_tmp(stage) = Wmmse(index_tmp,:)*y; % MMSE filtering (Eq.(11.9))

X_hat(stage) = QPSK_slicer(x_tmp(stage)); % Slicing

y_tilde = y-H(:,index_tmp)*X_hat(stage); % Interference subtraction

H_tilde = H(:,[1:index_tmp-1 index_tmp+1:NT-(stage-1)]); % new H

H = H_tilde;

y = y_tilde;

end

X_hat(Order) = X_hat;

elseif OSIC_type==2 % Column_norm ordering detection

X_hat=zeros(NT,1);

G = inv(H); % Inverse of H

for i=1:NT % Column_norm calculation

norm_array(i) = norm(H(:,i));

end

[sorted_norm_array,Order_tmp] = sort(norm_array);

Order = wrev(Order_tmp);

% V-BLAST

for stage=1:NT

x_tmp = G(Order(stage),:)*y; % Tx signal estimation

X_hat(Order(stage)) = QPSK_slicer(x_tmp); % Slicing

y_tilde = y-H(:,Order(stage))*X_hat(Order(stage));

end

else % OSIC with Post_detection_SNR ordering

Order=[];

index_array=[1:NT]; % Set of indices of signals to be detected

% V-BLAST

for stage=1:NT

G = inv(H’*H)*H’;

norm_array=[];

for i=1:NT-(stage-1) % Detection ordering

norm_array(i) = norm(G(i,:));

end

[val_min,index_min] = min(norm_array); % Ordering in SNR

Order = [Order index_array(index_min)];

index_array = index_array([1:index_min-1 index_min+1:end]);

x_tmp(stage) = G(index_min,:)*y; % Tx signal estimation

X_hat(stage) = QPSK_slicer(x_tmp(stage)); % Slicing

y_tilde = y-H(:,index_min)*X_hat(stage); % Interference subtraction

H_tilde = H(:,[1:index_min-1 index_min+1:NT-(stage-1)]); % New H

H = H_tilde; y = y_tilde;

end

X_hat(Order) = X_hat;

end

Problem is its a function and I have to compare the BER of these three orderings using VBLAST… the code only gives the detected symbols X_hat..

I have used the program below, giving arguments of the function OSIC_detector(y,H,sigma2,NT,OSIC_type)

clc;clear all;close all;

M=2 ;

N=2 ;

EbN0 = 30;

OSIC_type=4;

NT=2;

total_bits = 1000;

j = 1;

H=(randn(N,M)+j*randn(N,M))/sqrt(2);

for EbN0=0:5:30

for i = 1:total_bits/2

N0 = 0.1.^(EbN0/10);

sigma2=(0.5/(10^(EbN0/10)));

x = randint(1,2,[0 3])’;

x_qam = qammod(x,4);

v=(randn(N,1)+j*randn(N,1))* sqrt(N0/2);

y=H*x+v;

xHat = OSIC_detector(y,H,sigma2,NT,OSIC_type);

xHat_det = qamdemod(xHat,4);

xbin = dec2bin(x,2);

xHat_det_bin = dec2bin(xHat_det,2);

z = bitxor(uint8(xbin),uint8(xHat_det_bin));

error_ite(i) = sum(sum(z));

end

snr(j) = EbN0;

bit_err(j) = sum(error_ite)./total_bits;

j = j+1;

end

semilogy(snr,bit_err);

The QPSK_slicer i am using is:

function [x_sliced] = QPSK_slicer(x,n)

sq05=1/sqrt(2); jsq05=j*sq05; n=length(x)

for i=1:length(x)

if imag(x(i))>real(x(i))

if imag(x(i))>-real(x(i)), x_sliced(i)=jsq05;

else x_sliced(i)=-sq05;

end

else

if imag(x(i))>-real(x(i)), x_sliced(i)=sq05;

else x_sliced(i)=-jsq05;

end

end

end

As such i haven’t found any problems in the logic of the above program but the results are extremely incorrect!… the BER I am getting is different every time!

most probably the error comes from the OSIC_function program. But i am not able to figure them out

PROBLEM: I am not getting any expected BER! results are too weird!

Please check them and let me know the possible errors in my program.

@Areeb: Sorry, I wont be able to help with the coding. Please ensure that you have zero error rate in the case where noise is zero.

Hi Krishna.

In articles about V-BLAST, it is always noted that “the number of receive antennas Nr must be at least as large as the number of transmit antennas Nt”.

Do you know that is why?

@Dar: Well, one can think of the V-BLAST story as trying to solve for N unknowns, and to solve for N unknowns we need atleast N equations. The number of receive antennas form the N equations part. Helps?

So you are saying, since we are implementing linear decoding like ZF and MMSE; we need N equations for N unknowns. So if we implement decoding using ML rather than ZF and MMSE; do we still have this constraint? What do you think?

@Dar: Yes, the constraint holds good.

So you are saying since we are using linear decoding algorithms we need N linearly independent equations for solving N unknowns. So if we rather use ML decoding algorithm does it mean we will not have such a constraint. What dou you think?

@Dar: I think the constraint holds good.

Sir

Please provide me the matlab code for Analysis and design of OFDM-IDMA systems

please help me it’s urgent

my mail id is phani08473@gmail.com

@phani: sorry, do not know about ofdm-idma

Hello sir

If it is possible, please send me the Matlab code of the above simulation.

Thank you very much

My email address: samin6182@gmail.com

@Samin: The individual posts referred in this article has the simulation code

Hello Sir,

can I use this ‘V-BLAST’ simulation code to model a 2×2 MIMO system using spatial multiplexing and investigate the effect of the signal processing at the transmitter on the SNR at the receiver?

Many Thanks

@Mary: Yes, please do so.

hi , exiuseme i want matlab code Six equalizers for V-BLAST

BER for BPSK modulation with 2×2 MIMO Rayleigh channel

pleaze send me way email

my email is ghahraman26@yahoo.com

@ghahraman: please refer to the individual posts in http://www.dsplog.com/tag/mimo

hi i want matlab code. pleaz send way e-mail for me. very tanks

my email is ghahraman26@yahoo.com

@ghahraman: emailed you the instructions

hello sir,

I am student, I am working on 2X2 MIMO . I need a source about spatial Multiplexinf and space time coding for QAM Modulation or instructions for programming from you.

Thank you

@Pattaraporn:

Hope the links

http://www.dsplog.com/tag/mimo/

http://www.dsplog.com/tag/alamouti/

helps.

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

@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.

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

@aydar:

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?

@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.

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?

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

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

thanks

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

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

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

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.

@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.

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.

@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.

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.

@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)

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

@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.

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

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

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.

hi good thank you blog

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

@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.

@benjamin u can find contact of krishna at

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

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

@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.

@krishna

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

@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.

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

@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.

Hi,

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

@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.

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

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

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

Thanks

@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.

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

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

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

@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.

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.

@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

{ 2 trackbacks }