<?xml version="1.0" encoding="UTF-8"?>
<rss version="2.0"
	xmlns:content="http://purl.org/rss/1.0/modules/content/"
	xmlns:wfw="http://wellformedweb.org/CommentAPI/"
	xmlns:dc="http://purl.org/dc/elements/1.1/"
	xmlns:atom="http://www.w3.org/2005/Atom"
	xmlns:sy="http://purl.org/rss/1.0/modules/syndication/"
	xmlns:slash="http://purl.org/rss/1.0/modules/slash/"
	>

<channel>
	<title>dspLog &#187; Diversity</title>
	<atom:link href="http://www.dsplog.com/category/diversity/feed/" rel="self" type="application/rss+xml" />
	<link>http://www.dsplog.com</link>
	<description>Signal Processing for Communication</description>
	<lastBuildDate>Tue, 20 Jul 2010 16:36:56 +0000</lastBuildDate>
	<generator>http://wordpress.org/?v=2.9.2</generator>
	<language>en</language>
	<sy:updatePeriod>hourly</sy:updatePeriod>
	<sy:updateFrequency>1</sy:updateFrequency>
			<item>
		<title>Transmit beamforming</title>
		<link>http://www.dsplog.com/2009/04/13/transmit-beamforming/</link>
		<comments>http://www.dsplog.com/2009/04/13/transmit-beamforming/#comments</comments>
		<pubDate>Mon, 13 Apr 2009 15:53:26 +0000</pubDate>
		<dc:creator>Krishna Sankar</dc:creator>
				<category><![CDATA[Diversity]]></category>
		<category><![CDATA[Transmitter]]></category>
		<category><![CDATA[beamforming]]></category>
		<category><![CDATA[Rayleigh]]></category>

		<guid isPermaLink="false">http://www.dsplog.com/?p=535</guid>
		<description><![CDATA[In this post lets discuss a closed-loop transmit diversity scheme, where the transmitter has the knowledge of the channel. As there is a feedback path required from the receiver, to communicate the channel seen by the receiver to the transmitter, the scheme is called closed-loop transmit diversity scheme. Recall that the transmit diversity using Space [...]


Related posts:<ol><li><a href='http://www.dsplog.com/2008/09/19/equal-gain-combining/' rel='bookmark' title='Permanent Link: Equal Gain Combining (EGC)'>Equal Gain Combining (EGC)</a></li><li><a href='http://www.dsplog.com/2009/03/15/alamouti-stbc-2-receive-antenna/' rel='bookmark' title='Permanent Link: Alamouti STBC with 2 receive antenna'>Alamouti STBC with 2 receive antenna</a></li><li><a href='http://www.dsplog.com/2008/09/28/maximal-ratio-combining/' rel='bookmark' title='Permanent Link: Maximal Ratio Combining (MRC)'>Maximal Ratio Combining (MRC)</a></li></ol>]]></description>
			<content:encoded><![CDATA[<p></p><p>In this post lets discuss a <strong>closed-loop transmit diversity scheme</strong>, where the transmitter has the knowledge of the channel. As there is a feedback path required from the receiver, to communicate the channel seen by the receiver to the transmitter, the scheme is called closed-loop transmit diversity scheme. Recall that the <a href="http://www.dsplog.com/2008/10/16/alamouti-stbc/">transmit diversity using Space Time Coding (Alamouti STBC)</a> does not require the knowledge of the channel.  In this post, we will restrict our discussion to a 2 transmit, 1 receive case. We will assume that the channel is a flat fading <a title="Discussion on Rayleigh multipath channel" href="http://www.dsplog.com/2008/07/14/rayleigh-multipath-channel/" target="_self">Rayleigh multipath channel</a> and the modulation is BPSK.</p>
<p><span id="more-535"></span></p>
<h2>Channel Model</h2>
<p>1. We have 1 receive antennas and two transmit antenna.</p>
<p>2. The channel is flat fading &#8211; 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 <a title="Digital Communications: Fundamentals and Applications (2nd Edition) - Bernard Sklar (Author)" href="http://www.amazon.com/gp/redirect.html?ie=UTF8&amp;location=http%3A%2F%2Fwww.amazon.com%2FDigital-Communications-Fundamentals-Applications-Technologies%2Fdp%2F0130847887&amp;tag=dl04-20&amp;linkCode=ur2&amp;camp=1789&amp;creative=9325">[DIGITAL COMMUNICATIONS: SKLAR]</a><img style="border: medium none  ! important; margin: 0px ! important;" src="http://www.assoc-amazon.com/e/ir?t=dl04-20&amp;l=ur2&amp;o=1" border="0" alt="" width="1" height="1" /></p>
<p>3. The channel experienced by each receive antenna is randomly varying in time. For the <img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?i^{th}" border="0" alt="" align="absmiddle" /> receive antenna, each transmitted symbol gets multiplied by a randomly varying complex number <img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?h_i" border="0" alt="" align="absmiddle" />. As the channel under consideration is a Rayleigh channel, the real and imaginary parts of <img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?h_i" border="0" alt="" align="absmiddle" /> are Gaussian distributed having mean <img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?\mu_{h_i}=0" border="0" alt="" align="absmiddle" /> and variance <img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?\sigma^2_{h_i}=\frac{1}{2}" border="0" alt="" align="absmiddle" />.</p>
<p>4. The channel experience by each transmit antenna to receive antenna is independent from the channel experienced by other transmit antennas.</p>
<p>5. On each receive antenna, the noise<img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?n" border="0" alt="" align="absmiddle" /> has the Gaussian probability density function with</p>
<p><img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?p%28n%29%20=%20%5Cfrac%7B1%7D%7B%5Csqrt%7B2%5Cpi%5Csigma%5E2%7D%7De%5E%7B%5Cfrac%7B-%28n-%5Cmu%29%5E2%7D%7B2%5Csigma%5E2%7D" border="0" alt="" align="absmiddle" /> with <img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?%5Cmu=0" border="0" alt="" align="absmiddle" /> and <img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?%5Csigma%5E2%20=%20%5Cfrac%7BN_0%7D%7B2%7D" border="0" alt="" align="absmiddle" />.</p>
<p>6. At each transmit antenna, the channel <img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?h_i" border="0" alt="" align="absmiddle" /> is known.</p>
<p><img class="alignnone" title="2 Trasmit 1 Receive beam steering" src="http://www.dsplog.com/db-install/wp-content/uploads/2009/04/2tx_1rx_beamforming.png" alt="" /></p>
<p><strong>Figure: 2 transmit 1 receive beam steering</strong></p>
<h2>Transmit Beamforming</h2>
<p>On the receive antenna, the received signal is,</p>
<p><img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?y = \begin{eqnarray}[h_1 &amp; h_2]\end{eqnarray} \left[ \begin{eqnarray}x \\ \\ x \end{eqnarray}\right]+n = \underbrace{(h_1+h_2)}x + n" border="0" alt="" align="absmiddle" /> where,</p>
<p><img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?y" border="0" alt="" align="absmiddle" /> is the received symbol,<br />
<img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?h_i" border="0" alt="" align="absmiddle" /> is the channel on the <img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?i^{th}" border="0" alt="" align="absmiddle" /> transmit antenna,<br />
<img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?x" border="0" alt="" align="absmiddle" /> is the transmitted symbol and<br />
<img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?n" border="0" alt="" align="absmiddle" /> is the noise on the receive antenna.</p>
<p>When transmit beamforming is applied, we <strong>multiply the symbol from each transmit antenna with a complex number corresponding to the inverse of the phase of the channel</strong> so as to ensure that the signals add constructively at the receiver. In this scenario, the received signal is,</p>
<p><img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?y = \begin{eqnarray}[h_1 &amp; h_2]\end{eqnarray} \left[ \begin{eqnarray}e^{-j\theta_1} \\ \\ e^{-j\theta_2} \end{eqnarray}\right]x+n" border="0" alt="" align="absmiddle" />,</p>
<p>where,</p>
<p><img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?h_1 = |h_1|e^{j\theta_1}" border="0" alt="" align="absmiddle" /> and</p>
<p><img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?h_2 = |h_2|e^{j\theta_2}" border="0" alt="" align="absmiddle" />.</p>
<p>In this case, the signal at the receiver is,</p>
<p><img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?y = \underbrace{\left(|h_1| + |h_2|\right)}x + n" border="0" alt="" align="absmiddle" />.</p>
<p>For equalization, we need to divide the received symbol <img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?y" border="0" alt="" align="absmiddle" />with the new effective channel, i.e,</p>
<p><img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?\hat{y} = \frac{y}{\left(|h_1| + |h_2|\right)}=x + \frac{n}{\left(|h_1| + |h_2|\right)}" border="0" alt="" align="absmiddle" />.</p>
<h2>BER Simulation Model</h2>
<p>The Matlab/Octave script performs the following</p>
<p>(a) Generate random binary sequence of +1&#8217;s and -1&#8217;s.</p>
<p>(b) Multiply the symbols with the beam steering matrics &#8211; corresponding to the phase of the channel</p>
<p>(c) Perform equalization at the receiver</p>
<p>(d) Perform hard decision decoding and count the bit errors</p>
<p>(e) Repeat for multiple values of <img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?\frac{E_b}{N_0}" border="0" alt="" align="absmiddle" /> and plot the simulation and theoretical results.</p>
<p>Click here to download <a href="http://www.dsplog.com/db-install/wp-content/uploads/2009/04/script_ber_2transmit_1receive_beamforming_bpsk_rayleigh_channel.m">Matlab/Octave script for simulting BER for BPSK in flat fading Rayleigh channel with and without beamforming</a></p>
<p><a href="http://www.dsplog.com/db-install/wp-content/uploads/2009/04/script_ber_2transmit_1receive_beamforming_bpsk_rayleigh_channel.m"></a><img class="alignnone" title="BER plot for 2 transmit 1 receive beamforming for BPSK in Rayleigh channel" src="http://www.dsplog.com/db-install/wp-content/uploads/2009/04/ber_plot_2transmit_1receive_transmit_beamforming_with_bpsk_in_rayleigh_channel.png" alt="" width="448" height="336" /></p>
<p><strong>Figure: BER plot for 2 transmit 1 receive beamforming for BPSK in Rayleigh channel</strong></p>
<h2>Observations</h2>
<p>1. Sending the same information on multiple transmit antenna does not provide diversity gain. Intuituvely, this is due to the fact that the effective channel <img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?h_1 + h_2" border="0" alt="" align="absmiddle" /> in a 2 transmit antenna case is again a Rayleigh channel; hence the BER performance is identical to 1 transmit 1 receive Rayleigh channel case.</p>
<p>2. If the transmit symbols are multiplied by a complex phase to ensure that the phases align at the receiver, there is diversity gain. However, the BER performance seems to be slighly poorer than the 1 transmit 2 receive MRC case. I guess, the noise is scaled by <img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?|h_1| + |h_2|" border="0" alt="" align="absmiddle" /> in the case of transmit beamforming, whereas the noise scaling is different in the case of <a title="maximal ratio combining on dsplog.com" href="http://www.dsplog.com/2008/09/28/maximal-ratio-combining/">Maximal Ratio Combining</a>. I need to study bit more for a precise answer.</p>
<h2>Reference</h2>
<p><a title="Digital Communication: Third Edition, by John R. Barry, Edward A. Lee, David G. Messerschmitt" href="http://www.amazon.com/gp/redirect.html?ie=UTF8&amp;location=http%3A%2F%2Fwww.amazon.com%2FDigital-Communication-John-R-Barry%2Fdp%2F0792375483&amp;tag=dl04-20&amp;linkCode=ur2&amp;camp=1789&amp;creative=9325">[DIG-COMM-BARRY-LEE-MESSERSCHMITT] </a><a title="Digital Communication: Third Edition, by John R. Barry, Edward A. Lee, David G. Messerschmitt" href="http://www.amazon.com/gp/redirect.html?ie=UTF8&amp;location=http%3A%2F%2Fwww.amazon.com%2FDigital-Communication-John-R-Barry%2Fdp%2F0792375483&amp;tag=dl04-20&amp;linkCode=ur2&amp;camp=1789&amp;creative=9325">Digital Communication: Third Edition, by John R. Barry, Edward A. Lee, David G. Messerschmitt</a><img style="border-style: none ! important; margin: 0px" src="http://www.assoc-amazon.com/e/ir?t=dl04-20&amp;l=ur2&amp;o=1" border="0" alt="" width="1" height="1" /></p>
<p><a title="Note on Receive diversity by Prof. RaviRaj Adve" href="http://www.comm.utoronto.ca/~rsadve/Notes/DiversityReceive.pdf"></a><a title="IEEExplore link for the Alamouti STBC paper" href="http://ieeexplore.ieee.org/iel4/49/15739/00730453.pdf" target="_self"> </a></p>


<p>Related posts:<ol><li><a href='http://www.dsplog.com/2008/09/19/equal-gain-combining/' rel='bookmark' title='Permanent Link: Equal Gain Combining (EGC)'>Equal Gain Combining (EGC)</a></li><li><a href='http://www.dsplog.com/2009/03/15/alamouti-stbc-2-receive-antenna/' rel='bookmark' title='Permanent Link: Alamouti STBC with 2 receive antenna'>Alamouti STBC with 2 receive antenna</a></li><li><a href='http://www.dsplog.com/2008/09/28/maximal-ratio-combining/' rel='bookmark' title='Permanent Link: Maximal Ratio Combining (MRC)'>Maximal Ratio Combining (MRC)</a></li></ol></p>]]></content:encoded>
			<wfw:commentRss>http://www.dsplog.com/2009/04/13/transmit-beamforming/feed/</wfw:commentRss>
		<slash:comments>51</slash:comments>
		</item>
		<item>
		<title>Alamouti STBC with 2 receive antenna</title>
		<link>http://www.dsplog.com/2009/03/15/alamouti-stbc-2-receive-antenna/</link>
		<comments>http://www.dsplog.com/2009/03/15/alamouti-stbc-2-receive-antenna/#comments</comments>
		<pubDate>Sun, 15 Mar 2009 16:01:18 +0000</pubDate>
		<dc:creator>Krishna Sankar</dc:creator>
				<category><![CDATA[Diversity]]></category>
		<category><![CDATA[Alamouti]]></category>
		<category><![CDATA[STBC]]></category>

		<guid isPermaLink="false">http://www.dsplog.com/?p=501</guid>
		<description><![CDATA[In the past, we had discussed two transmit, one receive antenna Alamouti Space Time Block Coding (STBC) scheme. In this post, lets us discuss the impact of having two  antennas at the receiver. For the discussion, we will assume that the channel is a flat fading Rayleigh multipath channel and the modulation is BPSK.
Alamouti [...]


Related posts:<ol><li><a href='http://www.dsplog.com/2008/10/16/alamouti-stbc/' rel='bookmark' title='Permanent Link: Alamouti STBC'>Alamouti STBC</a></li><li><a href='http://www.dsplog.com/2008/10/24/mimo-zero-forcing/' rel='bookmark' title='Permanent Link: MIMO with Zero Forcing equalizer'>MIMO with Zero Forcing equalizer</a></li><li><a href='http://www.dsplog.com/2009/04/13/transmit-beamforming/' rel='bookmark' title='Permanent Link: Transmit beamforming'>Transmit beamforming</a></li></ol>]]></description>
			<content:encoded><![CDATA[<p></p><p>In the past, we had discussed <a title="Alamouti STBC in dsplog.com" href="http://www.dsplog.com/2008/10/16/alamouti-stbc/">two transmit, one receive antenna Alamouti Space Time Block Coding (STBC)</a> scheme. In this post, lets us discuss the impact of having two  antennas at the receiver. For the discussion, we will assume that the channel is a flat fading <a title="Discussion on Rayleigh multipath channel" href="http://www.dsplog.com/2008/07/14/rayleigh-multipath-channel/" target="_self">Rayleigh multipath channel</a> and the modulation is BPSK.</p>
<h2><span id="more-501"></span>Alamouti STBC with two receive antenna</h2>
<p>The principle of space time block coding with 2 transmit antenna and one receive antenna is explained in the post on <a title="2 transmit, 1 receive Alamouti STBC on dsplog.com" href="http://www.dsplog.com/2008/10/16/alamouti-stbc/">Alamouti STBC</a>. With two receive antenna&#8217;s the system can be modeled as shown in the figure below.</p>
<p><img class="alignnone" title="2 Transmit 2 Receive Alamouti STBC" src="http://www.dsplog.com/db-install/wp-content/uploads/2009/03/2transmit_2receive_alamouti_stbc.png" alt="" width="378" height="430" /></p>
<p><strong>Figure: 2 Transmit 2 Receive Alamouti STBC</strong></p>
<p>For discussion on the channel and noise model, please refer to the post on <a title="Alamouti STBC in dsplog.com" href="http://www.dsplog.com/2008/10/16/alamouti-stbc/">two transmit, one receive antenna Alamouti Space Time Block Coding (STBC)</a> scheme.</p>
<p>The received signal in the first time slot is,</p>
<p><img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?\left[\begin{eqnarray}y_1^1 \\ y_2^1\end{eqnarray}\right] = {\left[\begin{array}{rr}\ h_{11}\   h_{12} \\h_{21}\ h_{22}\end{enarray}\right]}\left[\begin{eqnarray}x_1 \\ x_2 \end{eqnarray}\right]+\left[\begin{array}n_1^1\\n_2^1 \end{eqnarray}\right]" border="0" alt="" align="absmiddle" />.</p>
<p>Assuming that the channel remains constant for the second time slot, the received signal is in the second time slot is,</p>
<p><img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?\left[\begin{eqnarray}y_1^2 \\ y_2^2\end{eqnarray}\right] = {\left[\begin{array}{rr}\ h_{11}\   h_{12} \\h_{21}\ h_{22}\end{enarray}\right]}\left[\begin{eqnarray}-x_2^* \\ x_1^* \end{eqnarray}\right]+\left[\begin{array}n_1^2\\n_2^2 \end{eqnarray}\right]" border="0" alt="" align="absmiddle" /></p>
<p>where</p>
<p><img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?\left[\begin{eqnarray}y_1^1 \\ y_2^1\end{eqnarray}\right] " border="0" alt="" align="absmiddle" /> are the received information at time slot 1 on receive antenna 1, 2 respectively,</p>
<p><img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?\left[\begin{eqnarray}y_1^2 \\ y_2^2\end{eqnarray}\right] " border="0" alt="" align="absmiddle" /> are the received information at time slot 2 on receive antenna 1, 2 respectively,</p>
<p><img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?h_{ij}" border="0" alt="" align="absmiddle" /> is the channel from <img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?i^{th}" border="0" alt="" align="absmiddle" /> receive antenna to <img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?j^{th}" border="0" alt="" align="absmiddle" /> transmit antenna,</p>
<p><img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?x_1" border="0" alt="" align="absmiddle" />, <img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?x_2" border="0" alt="" align="absmiddle" />are the  transmitted symbols,</p>
<p><img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?\left[\begin{eqnarray}n_1^1 \\ n_2^1\end{eqnarray}\right] " border="0" alt="" align="absmiddle" /> are the noise at time slot 1 on receive antenna 1, 2 respectively and</p>
<p><img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?\left[\begin{eqnarray}n_1^2 \\ n_2^2\end{eqnarray}\right] " border="0" alt="" align="absmiddle" /> are the noise at time slot 2 on receive antenna 1, 2 respectively.</p>
<p>Combining the equations at time slot 1 and 2,</p>
<p><img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?\left[\begin{array}{ll}y_1^1 \\ y_2^1 \\ y_1^2^* \\y_2^2^*\end{array}\right] = {\left[\begin{array}{rr}h_{11} &amp; h_{12}\\h_{21} &amp; h_{22}\\h_{12}^* &amp; -h_{11}^*\\h_{22}^* &amp; -h_{21}^* \end{enarray}\right]}\left[\begin{eqnarray}x_1 \\ x_2 \end{eqnarray}\right]+\left[\begin{array}n_1^1\\n_2^1 \\ n_1^2^*\\n_2^2^*\end{array}\right]" border="0" alt="" align="absmiddle" />.</p>
<p>Let us define <img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?\mathbf{H}= {\left[\begin{array}{rr}h_{11} &amp; h_{12}\\h_{21} &amp; h_{22}\\h_{12}^* &amp; -h_{11}^*\\h_{22}^* &amp; -h_{21}^* \end{enarray}\right]}" border="0" alt="" align="absmiddle" />.</p>
<p>To solve for <img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?\left[\begin{eqnarray}x_1 \\ x_2 \end{eqnarray}\right]" border="0" alt="" align="absmiddle" />, we know that we need to find the inverse of   <img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?\mathbf{H}" border="0" alt="" align="absmiddle" />.</p>
<p>We know, for a general m x n matrix, the <a title="Pseduo Inverse on Planet Math" href="http://planetmath.org/encyclopedia/Pseudoinverse.html">pseudo inverse</a> is defined as,</p>
<p><img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?\mathbf{H^+}=(H^HH)^{-1}H^H" border="0" alt="" align="absmiddle" />.</p>
<p>The term,</p>
<p><img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?(H^HH) =\left[\begin{array}{cc}|h_{11}|^2+|h_{21}|^2+ |h_{12}|^2 + |h_{22}|^2 &amp; 0 \\ 0 &amp; |h_{11}|^2+|h_{21}|^2+ |h_{12}|^2 + |h_{22}|^2\end{array}\right]" border="0" alt="" align="absmiddle" /></p>
<p>Since this is a diagonal matrix, the inverse is just the inverse of the diagonal elements, i.e</p>
<p><img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?(H^HH)^{-1} =\left[\begin{array}{cc}\frac{1}{|h_{11}|^2+|h_{21}|^2+ |h_{12}|^2 + |h_{22}|^2} &amp; 0 \\ 0 &amp; \frac{1}{|h_{11}|^2+|h_{21}|^2+ |h_{12}|^2 + |h_{22}|^2}\end{array}\right]" border="0" alt="" align="absmiddle" /></p>
<p>The estimate of the transmitted symbol is,</p>
<p><img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?\hat{\left[\begin{eqnarray}x_1 \\ x_2^*\end{eqnarray}\right]} = (H^HH)^{-1}H^H\left[\begin{array}{ll}y_1^1 \\ y_2^1 \\ y_1^2^* \\y_2^2^*\end{array}\right]" border="0" alt="" align="absmiddle" />.</p>
<h2>Simulation Model</h2>
<p>The Matlab/Octave script performs the following</p>
<p>(a) Generate random binary sequence of +1&#8217;s and -1&#8217;s.</p>
<p>(b) Group them into pair of two symbols</p>
<p>(c) Code it per the Alamouti Space Time code, multiply the symbols with the channel and then add white Gaussian noise.</p>
<p>(d) Equalize the received symbols</p>
<p>(e) Perform hard decision decoding and count the bit errors</p>
<p>(f) Repeat for multiple values of <img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?\frac{E_b}{N_0}" border="0" alt="" align="absmiddle" /> and plot the simulation and theoretical results.</p>
<p>Click here to download <a href="http://www.dsplog.com/db-install/wp-content/uploads/2009/03/script_ber_2x2_alamouti_stbc_code_bpsk_rayleigh_channel.m">Matlab/Octave script for computing BER for 2 transmit 2 receive Alamouti STBC for BPSK modulation in flat fading Rayleigh channel</a></p>
<p><img class="alignnone" title="BER plot for 2 transmit 2 receive alamouti STBC " src="http://www.dsplog.com/db-install/wp-content/uploads/2009/03/ber_plot_2transmit_2receive_alamouti_stbc_with_bpsk_in_rayleigh_channel.png" alt="" width="448" height="336" /></p>
<p><strong>Figure: BER plot for 2 transmi 2 receive Alamouti STBC</strong></p>
<h2>Observations</h2>
<p>1. Can observe that the BER performance is much better than 1 transmit 2 receive MRC case. This is because the effective channel concatenating the information from 2 receive antennas over two symbols results in a diversity order of 4.</p>
<p>2. In general, with <img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?m" border="0" alt="" align="absmiddle" /> receive antennas, the diversity order for 2 transmit antenna Alamouti STBC is  <img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?2m" border="0" alt="" align="absmiddle" />.</p>
<p>3. As with the case of <a title="alamouti STBC in dsplog.com" href="http://www.dsplog.com/2008/10/16/alamouti-stbc/">2 transmit, 1 receive Alamouti STBC</a>, the fact that <img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?(H^HH) " border="0" alt="" align="absmiddle" /> is a diagonal matrix ensured that there is no cross talk between <img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?x_1" border="0" alt="" align="absmiddle" />, <img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?x_2" border="0" alt="" align="absmiddle" /> after the equalizer and the noise term is still white.</p>
<h2>Reference</h2>
<p><a title="Digital Communication: Third Edition, by John R. Barry, Edward A. Lee, David G. Messerschmitt" href="http://www.amazon.com/gp/redirect.html?ie=UTF8&amp;location=http%3A%2F%2Fwww.amazon.com%2FDigital-Communication-John-R-Barry%2Fdp%2F0792375483&amp;tag=dl04-20&amp;linkCode=ur2&amp;camp=1789&amp;creative=9325">[DIG-COMM-BARRY-LEE-MESSERSCHMITT] </a><a title="Digital Communication: Third Edition, by John R. Barry, Edward A. Lee, David G. Messerschmitt" href="http://www.amazon.com/gp/redirect.html?ie=UTF8&amp;location=http%3A%2F%2Fwww.amazon.com%2FDigital-Communication-John-R-Barry%2Fdp%2F0792375483&amp;tag=dl04-20&amp;linkCode=ur2&amp;camp=1789&amp;creative=9325">Digital Communication: Third Edition, by John R. Barry, Edward A. Lee, David G. Messerschmitt</a><img style="border-style: none ! important; margin: 0px" src="http://www.assoc-amazon.com/e/ir?t=dl04-20&amp;l=ur2&amp;o=1" border="0" alt="" width="1" height="1" /></p>
<p><a title="Note on Receive diversity by Prof. RaviRaj Adve" href="http://www.comm.utoronto.ca/~rsadve/Notes/DiversityReceive.pdf">A</a><a title="IEEExplore link for the Alamouti STBC paper" href="http://ieeexplore.ieee.org/iel4/49/15739/00730453.pdf" target="_self"> Simple Transmit Diversity Technique for Wireless Communication </a>Siavash M Alamouti, IEEE Journal on selected areas in Communication, Vol 16, No, 8, October 1998</p>


<p>Related posts:<ol><li><a href='http://www.dsplog.com/2008/10/16/alamouti-stbc/' rel='bookmark' title='Permanent Link: Alamouti STBC'>Alamouti STBC</a></li><li><a href='http://www.dsplog.com/2008/10/24/mimo-zero-forcing/' rel='bookmark' title='Permanent Link: MIMO with Zero Forcing equalizer'>MIMO with Zero Forcing equalizer</a></li><li><a href='http://www.dsplog.com/2009/04/13/transmit-beamforming/' rel='bookmark' title='Permanent Link: Transmit beamforming'>Transmit beamforming</a></li></ol></p>]]></content:encoded>
			<wfw:commentRss>http://www.dsplog.com/2009/03/15/alamouti-stbc-2-receive-antenna/feed/</wfw:commentRss>
		<slash:comments>79</slash:comments>
		</item>
		<item>
		<title>Alamouti STBC</title>
		<link>http://www.dsplog.com/2008/10/16/alamouti-stbc/</link>
		<comments>http://www.dsplog.com/2008/10/16/alamouti-stbc/#comments</comments>
		<pubDate>Thu, 16 Oct 2008 01:16:32 +0000</pubDate>
		<dc:creator>Krishna Sankar</dc:creator>
				<category><![CDATA[Channel]]></category>
		<category><![CDATA[Diversity]]></category>
		<category><![CDATA[Error Rate]]></category>
		<category><![CDATA[Alamouti]]></category>
		<category><![CDATA[BPSK]]></category>
		<category><![CDATA[Rayleigh]]></category>
		<category><![CDATA[STBC]]></category>

		<guid isPermaLink="false">http://www.dsplog.com/?p=247</guid>
		<description><![CDATA[In the recent past, we have discussed three receive diversity schemes &#8211; Selection combining, Equal Gain Combining and Maximal Ratio Combining. All the three approaches used the antenna array at the receiver to improve the demodulation performance, albeit with different levels of complexity. Time to move on to a transmit diversity scheme where the information [...]


Related posts:<ol><li><a href='http://www.dsplog.com/2009/03/15/alamouti-stbc-2-receive-antenna/' rel='bookmark' title='Permanent Link: Alamouti STBC with 2 receive antenna'>Alamouti STBC with 2 receive antenna</a></li><li><a href='http://www.dsplog.com/2008/10/24/mimo-zero-forcing/' rel='bookmark' title='Permanent Link: MIMO with Zero Forcing equalizer'>MIMO with Zero Forcing equalizer</a></li><li><a href='http://www.dsplog.com/2009/04/13/transmit-beamforming/' rel='bookmark' title='Permanent Link: Transmit beamforming'>Transmit beamforming</a></li></ol>]]></description>
			<content:encoded><![CDATA[<p></p><p>In the recent past, we have discussed three receive diversity schemes &#8211; <a title="post on selection diversity in dsplog.com" href="http://www.dsplog.com/2008/09/06/receiver-diversity-selection-diversity/">Selection combining</a>, <a title="post on EGC in dsplog.com" href="http://www.dsplog.com/2008/09/19/equal-gain-combining/">Equal Gain Combining</a> and <a title="Post on MRC in dsplog.com" href="http://www.dsplog.com/2008/09/28/maximal-ratio-combining/">Maximal Ratio Combining</a>. All the three approaches used the antenna array at the receiver to improve the demodulation performance, albeit with different levels of complexity. Time to move on to a <strong>transmit diversity</strong> scheme where the information is spread across multiple antennas at the transmitter. In this post, lets discuss a popular transmit diversity scheme called <strong>Alamouti Space Time Block Coding</strong> <strong>(STBC)</strong>. For the discussion, we will assume that the channel is a flat fading <a title="Discussion on Rayleigh multipath channel" href="http://www.dsplog.com/2008/07/14/rayleigh-multipath-channel/" target="_self">Rayleigh multipath channel</a> and the modulation is BPSK.</p>
<p><span id="more-247"></span></p>
<h2>Alamouti STBC</h2>
<p>A simple Space Time Code, suggested by Mr. Siavash M Alamouti in his landmark October 1998 paper &#8211; A<a title="IEEExplore link for the Alamouti STBC paper" href="http://ieeexplore.ieee.org/iel4/49/15739/00730453.pdf" target="_self"> Simple Transmit Diversity Technique for Wireless Communication</a>, offers a simple method for achieving spatial diversity with two transmit antennas. The scheme is as follows:</p>
<p>1. Consider that we have a transmission sequence, for example <img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?\{x_1, x_2, x_3, \ldots, x_n \}" border="0" alt="" align="absmiddle" /></p>
<p>2. In normal transmission, we will be sending <img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?x_1" border="0" alt="" align="absmiddle" />in the first time slot, <img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?x_2" border="0" alt="" align="absmiddle" />in the second time slot, <img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?x_3" border="0" alt="" align="absmiddle" /> and so on.</p>
<p>3. However, Alamouti suggested that we group the symbols into groups of two. In the first time slot, send <img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?x_1" border="0" alt="" align="absmiddle" />and <img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?x_2" border="0" alt="" align="absmiddle" />from the first and second antenna. In second time slot send <img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?-x_2^*" border="0" alt="" align="absmiddle" /> and <img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?x_1^*" border="0" alt="" align="absmiddle" />from the first and second antenna. In the third time slot send <img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?x_3" border="0" alt="" align="absmiddle" /> and <img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?x_4" border="0" alt="" align="absmiddle" />from the first and second antenna.In fourth time slot, send <img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?-x_4^*" border="0" alt="" align="absmiddle" /> and <img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?x_3^*" border="0" alt="" align="absmiddle" />from the first and second antenna and so on.</p>
<p>4. Notice that though we are grouping two symbols, we still need two time slots to send two symbols. Hence, there is no change in the data rate.</p>
<p>5. This forms the simple explanation of the transmission scheme with Alamouti Space Time Block coding.</p>
<p><img class="alignnone size-full wp-image-249" title="2Transmit, 1Receive Alamouti STBC coding" src="http://www.dsplog.com/db-install/wp-content/uploads/2008/10/2tx_1rx_stbc_coding.png" alt="" width="400" height="472" /></p>
<p><strong>Figure: 2-Transmit, 1-Receive Alamouti STBC coding</strong></p>
<h2>Other Assumptions</h2>
<p>1. The channel is flat fading &#8211; 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 <a title="Digital Communications: Fundamentals and Applications (2nd Edition) - Bernard Sklar (Author)" href="http://www.amazon.com/gp/redirect.html?ie=UTF8&amp;location=http%3A%2F%2Fwww.amazon.com%2FDigital-Communications-Fundamentals-Applications-Technologies%2Fdp%2F0130847887&amp;tag=dl04-20&amp;linkCode=ur2&amp;camp=1789&amp;creative=9325">[DIGITAL COMMUNICATIONS: SKLAR]</a><img style="border: medium none  ! important; margin: 0px ! important;" src="http://www.assoc-amazon.com/e/ir?t=dl04-20&amp;l=ur2&amp;o=1" border="0" alt="" width="1" height="1" /></p>
<p>2. The channel experience by each transmit antenna is independent from the channel experienced by other transmit antennas.</p>
<p>3. For the <img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?i^{th}" border="0" alt="" align="absmiddle" /> transmit antenna, each transmitted symbol gets multiplied by a randomly varying complex number <img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?h_i" border="0" alt="" align="absmiddle" />. As the channel under consideration is a Rayleigh channel, the real and imaginary parts of <img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?h_i" border="0" alt="" align="absmiddle" /> are Gaussian distributed having mean <img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?\mu_{h_i}=0" border="0" alt="" align="absmiddle" /> and variance <img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?\sigma^2_{h_i}=\frac{1}{2}" border="0" alt="" align="absmiddle" />.</p>
<p>4. The channel experienced between each transmit to the receive antenna is randomly varying in time. However, the channel is <strong>assumed to remain constant over two time slots</strong>.</p>
<p>5. On the receive antenna, the noise<img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?n" border="0" alt="" align="absmiddle" /> has the Gaussian probability density function with</p>
<p><img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?p%28n%29%20=%20%5Cfrac%7B1%7D%7B%5Csqrt%7B2%5Cpi%5Csigma%5E2%7D%7De%5E%7B%5Cfrac%7B-%28n-%5Cmu%29%5E2%7D%7B2%5Csigma%5E2%7D" border="0" alt="" align="absmiddle" /> with <img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?%5Cmu=0" border="0" alt="" align="absmiddle" /> and <img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?%5Csigma%5E2%20=%20%5Cfrac%7BN_0%7D%7B2%7D" border="0" alt="" align="absmiddle" />.</p>
<p>7. The channel <img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?h_i" border="0" alt="" align="absmiddle" /> is known at the receiver.</p>
<h2>Receiver with Alamouti STBC</h2>
<p>In the first time slot, the received signal is,</p>
<p><img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?y_1 =h_1x_1 + h_2x_2 + n_1 = [h_1\  h_2] \left[\begin{eqnarray}x_1 \\ x_2 \end{eqnarray}\right]+n_1" border="0" alt="" align="absmiddle" />.</p>
<p>In the second time slot, the received signal is,</p>
<p><img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?y_2 =-h_1x_2^* + h_2x_1^* + n_2 = [h_1\  h_2] \left[\begin{eqnarray}-x_2^* \\ x_1^*\end{eqnarray}\right]+n_2" border="0" alt="" align="absmiddle" />.</p>
<p>where</p>
<p><img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?y_1" border="0" alt="" align="absmiddle" />, <img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?y_2" border="0" alt="" align="absmiddle" /> is the received symbol on the first and second time slot respectively,<br />
<img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?h_1" border="0" alt="" align="absmiddle" /> is the channel from <img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?1^{st}" border="0" alt="" align="absmiddle" /> transmit antenna to receive antenna,<br />
<img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?h_2" border="0" alt="" align="absmiddle" /> is the channel from <img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?2^{nd}" border="0" alt="" align="absmiddle" /> transmit antenna to receive antenna,<br />
<img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?x_1" border="0" alt="" align="absmiddle" />, <img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?x_2" border="0" alt="" align="absmiddle" />are the  transmitted symbols and<br />
<img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?n_1,\ n_2" border="0" alt="" align="absmiddle" /> is the noise on <img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?1^{st}, 2^{nd}" border="0" alt="" align="absmiddle" /> time slots.</p>
<p>Since the two noise terms are independent and identically distributed,</p>
<p><img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?E\left{\left[\begin{eqnarray}n_1\\n_2^*\end{eqnarray}\right]\left[\begin{eqnarray}n_1^*\ n_2\end{eqnarray}\right]\right}=\left[\begin{eqnarray}|n_1|^2\ \ \ \  0 \\0\ \ \ \ |n_2|^2\end{eqnarray}\right]" border="0" alt="" align="absmiddle" />.</p>
<p>For convenience, the above equation can be represented in matrix notation as follows:</p>
<p><img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?\left[\begin{eqnarray}y_1 \\ y_2^*\end{eqnarray}\right] = \underbrace{\left[\begin{eqnarray}\ h_1\ \ \  h_2 \\ h_2^*\ -h_1^*\end{enarray}\right]}\left[\begin{eqnarray}x_1 \\ x_2 \end{eqnarray}\right]+\left[\begin{eqnarray}n_1\\n_2^* \end{eqnarray}\right]" border="0" alt="" align="absmiddle" />.</p>
<p>Let us define <img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?\mathbf{H}= \left[\begin{eqnarray}\ h_1\ \ \  h_2 \\ h_2^*\ -h_1^*\end{enarray}\right]" border="0" alt="" align="absmiddle" />. To solve for <img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?\left[\begin{eqnarray}x_1 \\ x_2 \end{eqnarray}\right]" border="0" alt="" align="absmiddle" />, we know that we need to find the inverse of   <img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?\mathbf{H}" border="0" alt="" align="absmiddle" />.</p>
<p>We know, for a general m x n matrix, the <a title="Pseduo Inverse on Planet Math" href="http://planetmath.org/encyclopedia/Pseudoinverse.html">pseudo inverse</a> is defined as,</p>
<p><img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?\mathbf{H^+}=(H^HH)^{-1}H^H" border="0" alt="" align="absmiddle" />.</p>
<p>The term,</p>
<p><img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?(H^HH) = \left[\begin{eqnarray}\ h_1^*\ \ \  h_2 \\ h_2^*\ -h_1\end{eqnarray}\right]\left[\begin{eqnarray}\ h_1\ \ \  h_2 \\ h_2^*\ -h_1^*\end{eqnarray}\right]=\left[\begin{eqnarray}|h_1|^2+|h_2|^2\ \ \ \ \ \ \ \ \ \ \ \ \ \    0 \\ 0\ \ \ \ \ \ \ \ \ \ \ \ \|h_1|^2+|h_2|^2\end{eqnarray}\right]" border="0" alt="" align="absmiddle" />. Since this is a digonal matrix, the inverse is just the inverse of the diagonal elements, i.e</p>
<p><img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?(H^HH)^{-1} = \left[\begin{eqnarray}\frac{1}{|h_1|^2+|h_2|^2}\ \ \ \ \ \ \ \ \ \ \ \ \ \    0 \\ 0\ \ \ \ \ \ \ \ \ \ \ \ \frac{1}{\|h_1|^2+|h_2|^2}\end{eqnarray}\right]" border="0" alt="" align="absmiddle" />.</p>
<p>The estimate of the transmitted symbol is,</p>
<p><img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?\begin{eqnarray}\hat{\left[\begin{eqnarray}x_1 \\ x_2\end{eqnarray}\right]} &amp; = &amp;(H^HH)^{-1}H^H\left[\begin{eqnarray}y_1 \\ y_2^* \end{eqnarray}\right]\\&lt;br /&gt; &amp; = &amp; (H^HH)^{-1}H^H\left(H\left[\begin{eqnarray}x_1 \\ x_2 \end{eqnarray}\right]+\left[\begin{eqnarray}n_1\\n_2^* \end{eqnarray}\right]\right)&lt;br /&gt; \\&amp;=&amp;\left[\begin{eqnarray}x_1 \\ x_2 \end{eqnarray}\right] + (H^HH)^{-1}H^H\left[\begin{eqnarray}n_1\\n_2^* \end{eqnarray}\right]\\\end{eqnarray}" border="0" alt="" align="absmiddle" />.</p>
<p>If you compare the above equation with the estimated symbol following equalization in <a title="post on Maximal Ratio Combining " href="http://www.dsplog.com/2008/09/28/maximal-ratio-combining/#MRC">Maximal Ratio Combining</a>, you can see that the equations are identical.</p>
<h2>BER with Almouti STBC</h2>
<p>Since the estimate of the transmitted symbol with the Alamouti STBC scheme is identical to that obtained from MRC, the BER with above described Alamouti scheme should be same as that for MRC. However, there is a small catch.</p>
<p>With Alamouti STBC, we are transmitting from two antenna&#8217;s. Hence the total transmit power in the Alamouti scheme is twice that of that used in MRC. To make the comparison fair, we need to make the total trannsmit power from two antennas in STBC case to be equal to that of power transmitted from a single antenna in the MRC case. With this scaling, we can see that <strong>BER performance of 2Tx, 1Rx Alamouti STBC case has a roughly 3dB poorer performance that 1Tx, 2Rx MRC case</strong>.</p>
<p>From the post on <a title="post on MRC in dsplog.com" href="http://www.dsplog.com/2008/09/28/maximal-ratio-combining">Maximal Ratio Combining</a>, the bit error rate for BPSK modulation in Rayleigh channel with 1 transmit, 2 receive case is,</p>
<p><img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?P_{e,MRC} =p_{MRC}^2\left[1+2(1-p_{MRC})\right]" border="0" alt="" align="absmiddle" />, where</p>
<p><img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?p_{MRC}=\frac{1}{2}-\frac{1}{2}\left(1+\frac{1}{E_b/N_0}\right)^{-1/2}" border="0" alt="" align="absmiddle" />.</p>
<p>With <strong>Alamouti 2 transmit antenna, 1 receive antenna STBC</strong> case,</p>
<p><img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?p_{STBC}=\frac{1}{2}-\frac{1}{2}\left(1+\frac{2}{E_b/N_0}\right)^{-1/2}" border="0" alt="" align="absmiddle" /> and <strong>Bit Error Rate is</strong></p>
<p><img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?P_{e,STBC} =p_{STBC}^2\left[1+2(1-p_{STBC})\right]" border="0" alt="" align="absmiddle" />.</p>
<h2>Key points</h2>
<p>The fact that <img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?(H^HH) " border="0" alt="" align="absmiddle" /> is a diagonal matrix ensured the following:<br />
1. There is no cross talk between <img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?x_1" border="0" alt="" align="absmiddle" />, <img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?x_2" border="0" alt="" align="absmiddle" /> after the equalizer.</p>
<p>2. The noise term is still white.</p>
<p><img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?E\left{H^H\left[\begin{eqnarray}n_1\\n_2^*\end{eqnarray}\right]\left[\begin{eqnarray}n_1^*\ n_2\end{eqnarray}\right]H\right}=H^H\left[\begin{eqnarray}|n_1|^2\ \ \ \  0 \\0\ \ \ \ |n_2|^2\end{eqnarray}\right]H=\left[\begin{eqnarray}|n_1|^2\ \ \ \  0 \\0\ \ \ \ |n_2|^2\end{eqnarray}\right]\left[\begin{eqnarray}{|h_1|^2+|h_2|^2}\ \ \ \ \ \ \ \ \ \ \ \ \ \    0 \\ 0\ \ \ \ \ \ \ \ \ \ \ \ {\|h_1|^2+|h_2|^2}\end{eqnarray}\right]" border="0" alt="" align="absmiddle" />.</p>
<h2>Simulation Model</h2>
<p>The Matlab/Octave script performs the following</p>
<p>(a) Generate random binary sequence of +1&#8217;s and -1&#8217;s.</p>
<p>(b) Group them into pair of two symbols</p>
<p>(c) Code it per the Alamouti Space Time code, multiply the symbols with the channel and then add white Gaussian noise.</p>
<p>(d) Equalize the received symbols</p>
<p>(e) Perform hard decision decoding and count the bit errors</p>
<p>(f) Repeat for multiple values of <img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?\frac{E_b}{N_0}" border="0" alt="" align="absmiddle" /> and plot the simulation and theoretical results.</p>
<p>Click here to download <a href="http://www.dsplog.com/db-install/wp-content/uploads/2008/10/script_ber_alamouti_stbc_code_bpsk_rayleigh_channel.m">Matlab/Octave script for simulating BER for 2 transmit, 1 receive Alamouti STBC coding for BPSK modulation in Rayleigh fading channel</a></p>
<p><img class="alignnone size-full wp-image-251" title="BER plot for 2 transmit, 1 receive Alamouti STBC coding for BPSK modulation in Rayleigh fading channel" src="http://www.dsplog.com/db-install/wp-content/uploads/2008/10/ber_plot_alamouti_stbc_2tx_1rx_bpsk_rayleigh_channel.png" alt="" width="448" height="336" /></p>
<p><strong>Figure: BER plot for BPSK in Rayleigh channel with 2 Transmit and 1 Receive Alamouti STBC<br />
</strong></p>
<h2><strong>Observations</strong></h2>
<p>Compared to the BER plot for nTx=1, nRx=2 Maximal ratio combining, we can see the <strong>Alamouti Space Time Block Coding</strong> has around <strong>3dB poorer</strong> performance.</p>
<h2>Reference</h2>
<p><a title="Digital Communication: Third Edition, by John R. Barry, Edward A. Lee, David G. Messerschmitt" href="http://www.amazon.com/gp/redirect.html?ie=UTF8&amp;location=http%3A%2F%2Fwww.amazon.com%2FDigital-Communication-John-R-Barry%2Fdp%2F0792375483&amp;tag=dl04-20&amp;linkCode=ur2&amp;camp=1789&amp;creative=9325">[DIG-COMM-BARRY-LEE-MESSERSCHMITT] </a><a title="Digital Communication: Third Edition, by John R. Barry, Edward A. Lee, David G. Messerschmitt" href="http://www.amazon.com/gp/redirect.html?ie=UTF8&amp;location=http%3A%2F%2Fwww.amazon.com%2FDigital-Communication-John-R-Barry%2Fdp%2F0792375483&amp;tag=dl04-20&amp;linkCode=ur2&amp;camp=1789&amp;creative=9325">Digital Communication: Third Edition, by John R. Barry, Edward A. Lee, David G. Messerschmitt</a><img style="border-style: none ! important; margin: 0px" src="http://www.assoc-amazon.com/e/ir?t=dl04-20&amp;l=ur2&amp;o=1" border="0" alt="" width="1" height="1" /></p>
<p><a title="Note on Receive diversity by Prof. RaviRaj Adve" href="http://www.comm.utoronto.ca/~rsadve/Notes/DiversityReceive.pdf">A</a><a title="IEEExplore link for the Alamouti STBC paper" href="http://ieeexplore.ieee.org/iel4/49/15739/00730453.pdf" target="_self"> Simple Transmit Diversity Technique for Wireless Communication </a>Siavash M Alamouti, IEEE Journal on selected areas in Communication, Vol 16, No, 8, October 1998</p>


<p>Related posts:<ol><li><a href='http://www.dsplog.com/2009/03/15/alamouti-stbc-2-receive-antenna/' rel='bookmark' title='Permanent Link: Alamouti STBC with 2 receive antenna'>Alamouti STBC with 2 receive antenna</a></li><li><a href='http://www.dsplog.com/2008/10/24/mimo-zero-forcing/' rel='bookmark' title='Permanent Link: MIMO with Zero Forcing equalizer'>MIMO with Zero Forcing equalizer</a></li><li><a href='http://www.dsplog.com/2009/04/13/transmit-beamforming/' rel='bookmark' title='Permanent Link: Transmit beamforming'>Transmit beamforming</a></li></ol></p>]]></content:encoded>
			<wfw:commentRss>http://www.dsplog.com/2008/10/16/alamouti-stbc/feed/</wfw:commentRss>
		<slash:comments>164</slash:comments>
		</item>
		<item>
		<title>Maximal Ratio Combining (MRC)</title>
		<link>http://www.dsplog.com/2008/09/28/maximal-ratio-combining/</link>
		<comments>http://www.dsplog.com/2008/09/28/maximal-ratio-combining/#comments</comments>
		<pubDate>Sun, 28 Sep 2008 06:11:26 +0000</pubDate>
		<dc:creator>Krishna Sankar</dc:creator>
				<category><![CDATA[Diversity]]></category>
		<category><![CDATA[Error Rate]]></category>
		<category><![CDATA[BPSK]]></category>
		<category><![CDATA[PSK]]></category>
		<category><![CDATA[Rayleigh]]></category>

		<guid isPermaLink="false">http://www.dsplog.com/?p=241</guid>
		<description><![CDATA[This is the third post in the series discussing receiver diversity in a wireless link. Receiver diversity is a form of space diversity, where there are multiple antennas at the receiver. The presence of receiver diversity poses an interesting problem &#8211; how do we use &#8216;effectively&#8216; the information from all the antennas to demodulate the [...]


Related posts:<ol><li><a href='http://www.dsplog.com/2008/09/19/equal-gain-combining/' rel='bookmark' title='Permanent Link: Equal Gain Combining (EGC)'>Equal Gain Combining (EGC)</a></li><li><a href='http://www.dsplog.com/2008/09/06/receiver-diversity-selection-diversity/' rel='bookmark' title='Permanent Link: Selection Diversity'>Selection Diversity</a></li><li><a href='http://www.dsplog.com/2008/08/19/receive-diversity-in-awgn/' rel='bookmark' title='Permanent Link: Receive diversity in AWGN'>Receive diversity in AWGN</a></li></ol>]]></description>
			<content:encoded><![CDATA[<p></p><p>This is the third post in the series discussing <strong>receiver diversity</strong> in a wireless link. <strong>Receiver diversity</strong> is a form of space diversity, where there are multiple antennas at the receiver. The presence of receiver diversity poses an interesting problem &#8211; how do we use &#8216;<em>effectively</em>&#8216; the information from all the antennas to demodulate the data.  In the previous posts, we discussed <strong><a title="Discussion on Selection Diversity" href="http://www.dsplog.com/2008/09/06/receiver-diversity-selection-diversity/">selection diversity</a></strong> and <a title="Equal Gain Combining in dspLog.com" href="http://www.dsplog.com/2008/09/19/equal-gain-combining/"><strong>equal gain combining</strong> (EGC)</a>.</p>
<p>In this post, we will discuss<strong> Maximal Ratio Combining (MRC)</strong>. For the discussion, we will assume that the channel is a flat fading <a title="Discussion on Rayleigh multipath channel" href="http://www.dsplog.com/2008/07/14/rayleigh-multipath-channel/" target="_self">Rayleigh multipath channel</a> and the modulation is BPSK.<br />
<span id="more-241"></span></p>
<h2>Background</h2>
<p>We use the same constraints as defined in the <a title="Discussion on Selection Diversity" href="http://www.dsplog.com/2008/09/06/receiver-diversity-selection-diversity/">Selection Diversity</a> and  <a title="Equal Gain Combining in dspLog.com" href="http://www.dsplog.com/2008/09/19/equal-gain-combining/">Equal Gain Combining (EGC)</a> post. Let me repeat the same.</p>
<p>1. We have N receive antennas and one transmit antenna.</p>
<p>2. The channel is flat fading &#8211; 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 <a title="Digital Communications: Fundamentals and Applications (2nd Edition) - Bernard Sklar (Author)" href="http://www.amazon.com/gp/redirect.html?ie=UTF8&amp;location=http%3A%2F%2Fwww.amazon.com%2FDigital-Communications-Fundamentals-Applications-Technologies%2Fdp%2F0130847887&amp;tag=dl04-20&amp;linkCode=ur2&amp;camp=1789&amp;creative=9325">[DIGITAL COMMUNICATIONS: SKLAR]</a><img style="border: medium none  ! important; margin: 0px ! important;" src="http://www.assoc-amazon.com/e/ir?t=dl04-20&amp;l=ur2&amp;o=1" border="0" alt="" width="1" height="1" /></p>
<p>3. The channel experienced by each receive antenna is randomly varying in time. For the <img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?i^{th}" border="0" alt="" align="absmiddle" /> receive antenna, each transmitted symbol gets multiplied by a randomly varying complex number <img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?h_i" border="0" alt="" align="absmiddle" />. As the channel under consideration is a Rayleigh channel, the real and imaginary parts of <img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?h_i" border="0" alt="" align="absmiddle" /> are Gaussian distributed having mean <img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?\mu_{h_i}=0" border="0" alt="" align="absmiddle" /> and variance <img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?\sigma^2_{h_i}=\frac{1}{2}" border="0" alt="" align="absmiddle" />.</p>
<p>4. The channel experience by each receive antenna is independent from the channel experienced by other receive antennas.</p>
<p>5. On each receive antenna, the noise<img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?n" border="0" alt="" align="absmiddle" /> has the Gaussian probability density function with</p>
<p><img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?p%28n%29%20=%20%5Cfrac%7B1%7D%7B%5Csqrt%7B2%5Cpi%5Csigma%5E2%7D%7De%5E%7B%5Cfrac%7B-%28n-%5Cmu%29%5E2%7D%7B2%5Csigma%5E2%7D" border="0" alt="" align="absmiddle" /> with <img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?%5Cmu=0" border="0" alt="" align="absmiddle" /> and <img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?%5Csigma%5E2%20=%20%5Cfrac%7BN_0%7D%7B2%7D" border="0" alt="" align="absmiddle" />.</p>
<p>The noise on each receive antenna is independent from the noise on the other receive antennas.</p>
<p>6. At each receive antenna, the channel <img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?h_i" border="0" alt="" align="absmiddle" /> is known at the receiver.</p>
<p>7. In the presence of channel <img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?h_i" border="0" alt="" align="absmiddle" />, the instantaneous bit energy to noise ratio at <img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?i^{th}" border="0" alt="" align="absmiddle" /> receive antenna is <img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?\frac{|h_i|^2E_b}{N_0}" border="0" alt="" align="absmiddle" />. For notational convenience, let us define,</p>
<p><img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?\gamma_i = \frac{|h_i|^2E_b}{N_0}" border="0" alt="" align="absmiddle" />.<br />
<a name="MRC"></a></p>
<h2>Maximal Ratio Combining (MRC)</h2>
<p>On the <img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?i^{th}" border="0" alt="" align="absmiddle" /> receive antenna, the received signal is,</p>
<p><img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?y_i = h_ix+n_i" border="0" alt="" align="absmiddle" /> where</p>
<p><img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?y_i" border="0" alt="" align="absmiddle" /> is the received symbol on the <img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?i^{th}" border="0" alt="" align="absmiddle" /> receive antenna,<br />
<img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?h_i" border="0" alt="" align="absmiddle" /> is the channel on the <img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?i^{th}" border="0" alt="" align="absmiddle" /> receive antenna,<br />
<img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?x" border="0" alt="" align="absmiddle" /> is the transmitted symbol and<br />
<img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?n_i" border="0" alt="" align="absmiddle" /> is the noise on <img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?i^{th}" border="0" alt="" align="absmiddle" /> receive antenna.</p>
<p>Expressing it in matrix form, the received signal is,</p>
<p><img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?\mathbf{y} = \mathbf{h}x+\mathbf{n}" border="0" alt="" align="absmiddle" />, where</p>
<p><img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?\mathbf{y} = [y_1 y_2 \ldots y_N]^T" border="0" alt="" align="absmiddle" /> is the received symbol from all the receive antenna</p>
<p><img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?\mathbf{h} = [h_1 h_2 \ldots h_N]^T" border="0" alt="" align="absmiddle" /> is the channel on all the receive antenna</p>
<p><img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?x" border="0" alt="" align="absmiddle" /> is the transmitted symbol and</p>
<p><img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?\mathbf{n} = [n_1 n_2 \ldots n_N]^T" border="0" alt="" align="absmiddle" /> is the noise on all the receive antenna.</p>
<p>The equalized symbol is,</p>
<p><img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?\hat{\mathbf{x}} = \frac{\mathbf{h}^H\mathbf{y}}{\mathbf{h}^H\mathbf{h}} = \frac{\mathbf{h}^H\mathbf{h}x}{\mathbf{h}^H\mathbf{h}}+\frac{\mathbf{\mathbf{h}^Hn}}{\mathbf{h}^H\mathbf{h}} = x+\frac{\mathbf{\mathbf{h}^Hn}}{\mathbf{h}^H\mathbf{h}}" border="0" alt="" align="absmiddle" />.</p>
<p>It is intuitive to note that the term,</p>
<p><img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?\mathbf{h}^H\mathbf{h} = \sum_{i=1}^{N}|h_i|^2" border="0" alt="" align="absmiddle" /> i.e sum of the channel powers across all the receive antennas.</p>
<p>Note: The equations in the post refers the note on <a title="Note on Receive diversity by Prof. RaviRaj Adve" href="http://www.comm.utoronto.ca/~rsadve/Notes/DiversityReceive.pdf">Receive diversity</a> by Prof. RaviRaj Adve.</p>
<h2>Effective Eb/No with Maximal Ratio Combining (MRC)</h2>
<p>Earlier, we noted that in the presence of channel <img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?h_i" border="0" alt="" align="absmiddle" />, the instantaneous bit energy to noise ratio at <img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?i^{th}" border="0" alt="" align="absmiddle" /> receive antenna is</p>
<p><img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?\gamma_i = \frac{|h_i|^2E_b}{N_0}" border="0" alt="" align="absmiddle" />.</p>
<p>Given that we are equalizing the channel with <img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?\mathbf{h}^H" border="0" alt="" align="absmiddle" />, with the <img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?N" border="0" alt="" align="absmiddle" /> receive antenna case, the effective bit energy to noise ratio is,</p>
<p><img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?\begin{eqnarray}\gamma &amp;= &amp;\sum_{i=1}^{N}\frac{|h_i|^2E_b}{N_0}\\ &amp; &amp; \\&amp; = &amp; N\gamma_i\end{eqnarray}" border="0" alt="" align="absmiddle" />.</p>
<p><strong>E</strong><strong>ffective bit energy to noise ratio in a N receive antenna case is N times the bit energy to noise ratio for single antenna case</strong>. Recall, this gain is same as the improvement which we got in <a title="Receive Diversity in AWGN" href="http://www.dsplog.com/2008/08/19/receive-diversity-in-awgn/">Receive diversity for AWGN</a> case <img src='http://www.dsplog.com/db-install/wp-includes/images/smilies/icon_smile.gif' alt=':)' class='wp-smiley' /> </p>
<p>Click here to download <a href="http://www.dsplog.com/db-install/wp-content/uploads/2008/09/script_maximal_ratio_combining_effective_snr.m">Matlab/Octave script for plotting effective SNR with Maximal Ratio Combining in Rayleigh channel</a></p>
<p><img class="alignnone size-full wp-image-243" title="Plot of Effective SNR with Maximal Ratio Combining in Rayleigh fading channel" src="http://www.dsplog.com/db-install/wp-content/uploads/2008/09/plot_effective_snr_maximal_ratio_combining_bpsk_rayleigh_channel.png" alt="" width="448" height="336" /></p>
<p><strong>Figure: Effective SNR with Maximal Ratio Combining in Rayleigh fading channel</strong></p>
<h2>Error rate with Maximal Ratio Combining (MRC)</h2>
<p>From the discussion on <a title="Chi-square random variable with two degrees of freedom" href="http://www.dsplog.com/2008/07/28/chi-square-random-variable/#twodegree">chi-square random variable</a>, we know that, if <img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?h_i" border="0" alt="" align="absmiddle" /> is a Rayleigh distributed random variable, then <img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?h_i^2" border="0" alt="" align="absmiddle" /> is a chi-squared random variable with two degrees of freedom. The pdf of <img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?\gamma_i" border="0" alt="" align="absmiddle" /> is<br />
<img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?p(\gamma_i)=\frac{1}{(E_b/N_0)}e^{\frac{-\gamma_i}{(E_b/N_0)}" border="0" alt="" align="absmiddle" />.</p>
<p>Since the effective bit energy to noise ratio <img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?\gamma" border="0" alt="" align="absmiddle" /> is the sum of <img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?N" border="0" alt="" align="absmiddle" />such random variables, the pdf of <img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?\gamma" border="0" alt="" align="absmiddle" /> is a <a title="Chi Square Random variable with m degrees of freedom" href="http://www.dsplog.com/2008/07/28/chi-square-random-variable/#mdegree">chi-square random variable</a> with <img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?2N" border="0" alt="" align="absmiddle" />degrees of freedom. The pdf of <img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?\gamma" border="0" alt="" align="absmiddle" />is,<br />
<img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?p(\gamma) = \frac{1}{(N-1)!(E_b/N_0)^N}\ \gamma^{N-1}e^{\frac{-\gamma}{(E_b/N_0)}},\ \ \ \  \gamma \ge 0" border="0" alt="" align="absmiddle" />.</p>
<p>If you recall, in the post on <a title="BER for BPSK in AWGN" href="http://www.dsplog.com/2007/08/05/bit-error-probability-for-bpsk-modulation/">BER computation in AWGN</a>, with bit energy to noise ratio of <img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?\frac{E_b}{N_0}" border="0" alt="" align="absmiddle" />, the bit error rate for BPSK in AWGN is derived as</p>
<p><img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?P_b=%5Cfrac%7B1%7D%7B2%7Derfc%5Cleft%28%7B%5Csqrt%7B%5Cfrac%7BE_b%7D%7BN_0%7D%7D%7D%5Cright%29" border="0" alt="" align="absmiddle" />.</p>
<p>Given that the effective bit energy to noise ratio with <strong>maximal ratio combining</strong> is<img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?\gamma" border="0" alt="" align="absmiddle" />, the total bit error rate is the integral of the conditional BER integrated over all possible values of <img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?\gamma" border="0" alt="" align="absmiddle" />.</p>
<p><img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?\begin{eqnarray}P_{e}&amp;=&amp;\int_0^{\infty}\frac{1}{2}erfc\left(\sqrt{\gamma}\right)p(\gamma)d\gamma\\&lt;br /&gt; &amp;=&amp;\int_0^{\infty}\frac{1}{2}erfc\left(\sqrt{\gamma}\right) \frac{1}{(N-1)!(E_b/N_0)^N}\ \gamma^{N-1}e^{\frac{-\gamma}{(E_b/N_0)}}d\gamma\end{eqnarray}" border="0" alt="" align="absmiddle" />.</p>
<p>This equation reduces to</p>
<p><img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?\begin{eqnarray}P_{e}&amp;=p^N\sum_{k=0}^{N-1}\left(\begin{eqnarray}N-1+k\\k\end{eqnarray}\right)\left(1-p\right)^k\end{eqnarray}" border="0" alt="" align="absmiddle" />, where</p>
<p><img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?p=\frac{1}{2}-\frac{1}{2}\left(1+\frac{1}{E_b/N_0}\right)^{-1/2}" border="0" alt="" align="absmiddle" />.</p>
<p>Refer Equation 11.12 and Equation 11.13 in Section 11.3.1 Performance with Maximal Ratio Combining in <a title="Digital Communication: Third Edition, by John R. Barry, Edward A. Lee, David G. Messerschmitt" href="http://www.amazon.com/gp/redirect.html?ie=UTF8&amp;location=http%3A%2F%2Fwww.amazon.com%2FDigital-Communication-John-R-Barry%2Fdp%2F0792375483&amp;tag=dl04-20&amp;linkCode=ur2&amp;camp=1789&amp;creative=9325">[DIG-COMM-BARRY-LEE-MESSERSCHMITT].</a> Again, I do not know the proof <img src='http://www.dsplog.com/db-install/wp-includes/images/smilies/icon_sad.gif' alt=':(' class='wp-smiley' /> </p>
<h2>BER Simulation Model</h2>
<p>The Matlab/Octave script performs the following</p>
<p>(a) Generate random binary sequence of +1&#8217;s and -1&#8217;s.</p>
<p>(b) Multiply the symbols with the channel and then add white Gaussian noise.</p>
<p>(c) Chose that receive path, equalize the received symbols per maximal ratio combining</p>
<p>(d) Perform hard decision decoding and count the bit errors</p>
<p>(e) Repeat for multiple values of <img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?\frac{E_b}{N_0}" border="0" alt="" align="absmiddle" /> and plot the simulation and theoretical results.</p>
<p>Click here to download <a title="Matlab/Octave script for simulating BER for BPSK in Rayleigh channel with Maximal Ratio Combining" href="http://www.dsplog.com/db-install/wp-content/uploads/2008/09/script_ber_bpsk_rayleigh_channel_maximal_ratio_combining.m">Matlab/Octave script for simulating BER for BPSK in Rayleigh channel with Maximal Ratio Combining</a></p>
<p><img class="alignnone size-full wp-image-245" title="BER plot for BPSK in Rayleigh channel with Maximal Ratio Combining" src="http://www.dsplog.com/db-install/wp-content/uploads/2008/09/ber_plot_maximal_ratio_combining_bpsk_rayleigh_channel.png" alt="" width="448" height="336" /></p>
<p><strong>Figure: BER plot for BPSK in Rayleigh channel with Maximal Ratio Combining</strong></p>
<p><strong><br />
</strong></p>
<h2>Reference</h2>
<p><a title="Digital Communication: Third Edition, by John R. Barry, Edward A. Lee, David G. Messerschmitt" href="http://www.amazon.com/gp/redirect.html?ie=UTF8&amp;location=http%3A%2F%2Fwww.amazon.com%2FDigital-Communication-John-R-Barry%2Fdp%2F0792375483&amp;tag=dl04-20&amp;linkCode=ur2&amp;camp=1789&amp;creative=9325">[DIG-COMM-BARRY-LEE-MESSERSCHMITT] </a><a title="Digital Communication: Third Edition, by John R. Barry, Edward A. Lee, David G. Messerschmitt" href="http://www.amazon.com/gp/redirect.html?ie=UTF8&amp;location=http%3A%2F%2Fwww.amazon.com%2FDigital-Communication-John-R-Barry%2Fdp%2F0792375483&amp;tag=dl04-20&amp;linkCode=ur2&amp;camp=1789&amp;creative=9325">Digital Communication: Third Edition, by John R. Barry, Edward A. Lee, David G. Messerschmitt</a><img style="border-style: none ! important; margin: 0px" src="http://www.assoc-amazon.com/e/ir?t=dl04-20&amp;l=ur2&amp;o=1" border="0" alt="" width="1" height="1" /></p>
<p><a title="Note on Receive diversity by Prof. RaviRaj Adve" href="http://www.comm.utoronto.ca/~rsadve/Notes/DiversityReceive.pdf">Receive diversity &#8211; Notes by Prof. Raviraj Adve</a></p>


<p>Related posts:<ol><li><a href='http://www.dsplog.com/2008/09/19/equal-gain-combining/' rel='bookmark' title='Permanent Link: Equal Gain Combining (EGC)'>Equal Gain Combining (EGC)</a></li><li><a href='http://www.dsplog.com/2008/09/06/receiver-diversity-selection-diversity/' rel='bookmark' title='Permanent Link: Selection Diversity'>Selection Diversity</a></li><li><a href='http://www.dsplog.com/2008/08/19/receive-diversity-in-awgn/' rel='bookmark' title='Permanent Link: Receive diversity in AWGN'>Receive diversity in AWGN</a></li></ol></p>]]></content:encoded>
			<wfw:commentRss>http://www.dsplog.com/2008/09/28/maximal-ratio-combining/feed/</wfw:commentRss>
		<slash:comments>90</slash:comments>
		</item>
		<item>
		<title>Equal Gain Combining (EGC)</title>
		<link>http://www.dsplog.com/2008/09/19/equal-gain-combining/</link>
		<comments>http://www.dsplog.com/2008/09/19/equal-gain-combining/#comments</comments>
		<pubDate>Fri, 19 Sep 2008 17:53:26 +0000</pubDate>
		<dc:creator>Krishna Sankar</dc:creator>
				<category><![CDATA[Diversity]]></category>
		<category><![CDATA[Error Rate]]></category>
		<category><![CDATA[BPSK]]></category>
		<category><![CDATA[PSK]]></category>
		<category><![CDATA[Rayleigh]]></category>

		<guid isPermaLink="false">http://www.dsplog.com/?p=236</guid>
		<description><![CDATA[This is the second post in the series discussing receiver diversity in a wireless link. Receiver diversity is a form of space diversity, where there are multiple antennas at the receiver. The presence of receiver diversity poses an interesting problem &#8211; how do we use &#8216;effectively&#8216; the information from all the antennas to demodulate the [...]


Related posts:<ol><li><a href='http://www.dsplog.com/2008/09/28/maximal-ratio-combining/' rel='bookmark' title='Permanent Link: Maximal Ratio Combining (MRC)'>Maximal Ratio Combining (MRC)</a></li><li><a href='http://www.dsplog.com/2008/09/06/receiver-diversity-selection-diversity/' rel='bookmark' title='Permanent Link: Selection Diversity'>Selection Diversity</a></li><li><a href='http://www.dsplog.com/2008/08/19/receive-diversity-in-awgn/' rel='bookmark' title='Permanent Link: Receive diversity in AWGN'>Receive diversity in AWGN</a></li></ol>]]></description>
			<content:encoded><![CDATA[<p></p><p>This is the second post in the series discussing <strong>receiver diversity</strong> in a wireless link. <strong>Receiver diversity</strong> is a form of space diversity, where there are multiple antennas at the receiver. The presence of receiver diversity poses an interesting problem &#8211; how do we use &#8216;<em>effectively</em>&#8216; the information from all the antennas to demodulate the data.  In the previous post, we discussed <strong><a title="Discussion on Selection Diversity" href="http://www.dsplog.com/2008/09/06/receiver-diversity-selection-diversity/">selection diversity</a></strong>. In this post, we will discuss <strong>equal gain combining</strong> (EGC). For the discussion, we will assume that the channel is a flat fading <a title="Discussion on Rayleigh multipath channel" href="http://www.dsplog.com/2008/07/14/rayleigh-multipath-channel/" target="_self">Rayleigh multipath channel</a> and the modulation is BPSK.</p>
<p><span id="more-236"></span></p>
<h2>Background</h2>
<p>We use the same constraints as defined in the <a title="Discussion on Selection Diversity" href="http://www.dsplog.com/2008/09/06/receiver-diversity-selection-diversity/">Selection Diversity</a> post. Let me repeat the same.</p>
<p>1. We have N receive antennas and one transmit antenna.</p>
<p>2. The channel is flat fading &#8211; 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 <a title="Digital Communications: Fundamentals and Applications (2nd Edition) - Bernard Sklar (Author)" href="http://www.amazon.com/gp/redirect.html?ie=UTF8&amp;location=http%3A%2F%2Fwww.amazon.com%2FDigital-Communications-Fundamentals-Applications-Technologies%2Fdp%2F0130847887&amp;tag=dl04-20&amp;linkCode=ur2&amp;camp=1789&amp;creative=9325">[DIGITAL COMMUNICATIONS: SKLAR]</a><img style="border: medium none  ! important; margin: 0px ! important;" src="http://www.assoc-amazon.com/e/ir?t=dl04-20&amp;l=ur2&amp;o=1" border="0" alt="" width="1" height="1" /></p>
<p>3. The channel experienced by each receive antenna is randomly varying in time. For the <img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?i^{th}" border="0" alt="" align="absmiddle" /> receive antenna, each transmitted symbol gets multiplied by a randomly varying complex number <img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?h_i" border="0" alt="" align="absmiddle" />. As the channel under consideration is a Rayleigh channel, the real and imaginary parts of <img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?h_i" border="0" alt="" align="absmiddle" /> are Gaussian distributed having mean <img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?\mu_{h_i}=0" border="0" alt="" align="absmiddle" /> and variance <img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?\sigma^2_{h_i}=\frac{1}{2}" border="0" alt="" align="absmiddle" />.</p>
<p>4. The channel experience by each receive antenna is independent from the channel experienced by other receive antennas.</p>
<p>5. On each receive antenna, the noise<img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?n" border="0" alt="" align="absmiddle" /> has the Gaussian probability density function with</p>
<p><img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?p%28n%29%20=%20%5Cfrac%7B1%7D%7B%5Csqrt%7B2%5Cpi%5Csigma%5E2%7D%7De%5E%7B%5Cfrac%7B-%28n-%5Cmu%29%5E2%7D%7B2%5Csigma%5E2%7D" border="0" alt="" align="absmiddle" /> with <img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?%5Cmu=0" border="0" alt="" align="absmiddle" /> and <img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?%5Csigma%5E2%20=%20%5Cfrac%7BN_0%7D%7B2%7D" border="0" alt="" align="absmiddle" />.</p>
<p>The noise on each receive antenna is independent from the noise on the other receive antennas.</p>
<p>6. At each receive antenna, the channel <img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?h_i" border="0" alt="" align="absmiddle" /> is known at the receiver.</p>
<p>7. In the presence of channel <img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?h_i" border="0" alt="" align="absmiddle" />, the instantaneous bit energy to noise ratio at <img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?i^{th}" border="0" alt="" align="absmiddle" /> receive antenna is <img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?\frac{|h_i|^2E_b}{N_0}" border="0" alt="" align="absmiddle" />. For notational convenience, let us define,</p>
<p><img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?\gamma_i = \frac{|h_i|^2E_b}{N_0}" border="0" alt="" align="absmiddle" />.</p>
<p>From the discussion on <a title="Chi-square random variable with two degrees of freedom" href="http://www.dsplog.com/2008/07/28/chi-square-random-variable/#twodegree">chi-square random variable</a>, we know that, if <img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?h_i" border="0" alt="" align="absmiddle" /> is a Rayleigh distributed random variable, then <img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?h_i^2" border="0" alt="" align="absmiddle" /> is a chi-squared random variable with two degrees of freedom. The pdf of <img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?\gamma_i" border="0" alt="" align="absmiddle" /> is<br />
<img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?p(\gamma_i)=\frac{1}{(E_b/N_0)}e^{\frac{-\gamma_i}{(E_b/N_0)}" border="0" alt="" align="absmiddle" />.</p>
<h2>Equal Gain Combining (EGC)</h2>
<p>On the <img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?i^{th}" border="0" alt="" align="absmiddle" /> receive antenna, equalization is performed at the receiver by dividing the received symbol <img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?y_i" border="0" alt="" align="absmiddle" /> by the apriori known phase of <img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?h_i" border="0" alt="" align="absmiddle" />. The channel <img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?h_i" border="0" alt="" align="absmiddle" /> is  represented in polar form as <img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?|h_i|e^{j\theta_i}" border="0" alt="" align="absmiddle" />. The decoded symbol is the sum of the phase compensated channel from all the receive antennas.</p>
<p><img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?\begin{eqnarray}\hat{y} &amp; = &amp;\sum_i\frac{y_i}{e^{j\theta_i}} \\ &amp; = &amp;\sum_i\frac{|h_i|e^{j\theta_i}x+n_i}{e^{j\theta_i}} \\ &amp; = &amp;\sum_i |h_i|x+\tilde{n_i}\end{eqnarray}" border="0" alt="" align="absmiddle" /></p>
<p>where<br />
<img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?\tilde{n_i} = \frac{n_i}{e^{j\theta_i}}" border="0" alt="" align="absmiddle" />is the additive noise scaled by the phase of the channel coefficient.</p>
<p>For demodulation, we use the classical definition i.e.</p>
<p><img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?\hat y &gt;0 \rightarrow 1" border="0" alt="" align="absmiddle" /> and <img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?\hat y \le 0 \rightarrow 0" border="0" alt="" align="absmiddle" /></p>
<p><strong>Note:</strong></p>
<p>For PSK modulation schemes, the equalization by the phase of the channel coefficients suffice. However, for QAM case, we need to compensate for the amplitude also when equalizing. <strong>We are NOT discussing QAM case in this post</strong>.</p>
<h2>Effective Eb/N0 with Equal Gain Combining</h2>
<p>The equations listed below obtained from the article <a title="Note on Receive diversity by Prof. RaviRaj Adve" href="http://www.comm.utoronto.ca/~rsadve/Notes/DiversityReceive.pdf">Receive diversity &#8211; Notes by Prof. Raviraj Adve</a>.</p>
<p>In the presence of channel <img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?h_i" border="0" alt="" align="absmiddle" />, the instantaneous bit energy to noise ratio at <img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?i^{th}" border="0" alt="" align="absmiddle" /> receive antenna is <img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?\frac{|h_i|^2E_b}{N_0}" border="0" alt="" align="absmiddle" />. For notational convenience, let us define,</p>
<p><img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?\gamma_i = \frac{|h_i|^2E_b}{N_0}" border="0" alt="" align="absmiddle" />.</p>
<p>The effective Eb/N0 with equal gain combining is the channel power accumulated over all receive chains, i.e.</p>
<p><img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?\begin{eqnarray}E(\gamma_i) &amp; = &amp;\frac{E_b}{N_0}\frac{1}{N}\left(\sum_{i=1}^{N}|h_i|^2\right)\\&amp;=&amp;\frac{E_b}{N_0}\frac{1}{N}\left(\sum_{i=1}^{N}\sum_{k=1}^{N}|h_i||h_k| \right)\\&amp;=&amp;\frac{E_b}{N_0}\frac{1}{N}\left(\sum_{i=1}^{N}|h_i|^2 + \sum_{i=1}^{N}\sum_{k=1,k\ne n}^{N}|h_i||h_k| \right)\end{eqnarray}" border="0" alt="" align="absmiddle" />.</p>
<p>The first term is chi-square random variable with <img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?2N" border="0" alt="" align="absmiddle" /> degrees of freedom having mean value of <img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?2N\sigma^2_{h_i} " border="0" alt="" align="absmiddle" />. Hence the first term reduces to,</p>
<p><img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?\begin{eqnarray}\sum_{i=1}^{N}|h_i|^2=N\end{eqnarray}" border="0" alt="" align="absmiddle" />.</p>
<p>The second term is a product of two Rayleigh random variables. The mean of Rayleigh random variable with variance <img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?\sigma^2_{h_i}" border="0" alt="" align="absmiddle" /> is <img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?\sigma_{h_i}\sqrt{\frac{\pi}{2}}" border="0" alt="" align="absmiddle" />. Hence the second term is,</p>
<p><img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?\begin{eqnarray}\sum_{i=1}^{N}\sum_{k=1,k\ne n}^{N}|h_i||h_k|)=N\sqrt{\frac{\pi}{4}}(N-1)\sqrt{\frac{\pi}{4}}=N(N-1)\frac{\pi}{4}\end{eqnarray}" border="0" alt="" align="absmiddle" />.</p>
<p>Simplifying, <strong>the effective Eb/N0 with equal gain combining is</strong>,</p>
<p><img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?\begin{eqnarray}E(\gamma_i) &amp; = &amp;\frac{E_b}{N_0}\frac{1}{N}\left[N+ N(N-1)\frac{\pi}{4} \right]\\&amp; = &amp; \frac{E_b}{N_0}\left[1+(N-1)\frac{\pi}{4}\right]\end{eqnarray}" border="0" alt="" align="absmiddle" />.</p>
<p>Click here to download <a href="http://www.dsplog.com/db-install/wp-content/uploads/2008/09/script_equal_gain_combining_effective_snr.m">Matlab/Ocatve script for computing effective SNR with equal gain combining</a></p>
<p><img class="alignnone size-full wp-image-240" title="Gain in Eb/N0 with equal gain combining" src="http://www.dsplog.com/db-install/wp-content/uploads/2008/09/snr_gain_equal_gain_combining.png" alt="" width="494" height="404" /></p>
<p><strong>Figure : Gain in Eb/N0 with equal gain combining</strong></p>
<h2>Bit Error Rate with Equal Gain Combining</h2>
<p>The IEEE paper <a title="BER with EGC in Rayleigh channel" href="http://ieeexplore.ieee.org/iel1/26/12187/00558680.pdf?tp=&amp;arnumber=558680&amp;isnumber=12187">[ZHANG97] </a>discuss the BER computation with Equal Gain Combining. Since the proof is tedious (and I did not understand) I am just noting the final results.</p>
<p>With two receive antennas, the BER with equal gain combining is,</p>
<p><img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?P_e = \frac{1}{2}\left[1-\frac{\sqrt{E_b/N_0(E_b/N_0+2)}}{E_b/N_0+1}\right]" border="0" alt="" align="absmiddle" />.</p>
<h2>Matlab model for simulating BER with Equal Gain Combining</h2>
<p>The Matlab/Octave script performs the following</p>
<p>(a) Generate random binary sequence of +1&#8217;s and -1&#8217;s.</p>
<p>(b) Multiply the symbols with the channel and then add white Gaussian noise.</p>
<p>(c) At the receiver, for each receive path, equalize by compensating with the known channel phase</p>
<p>(d) Accumulate the equalized symbols from all the receive paths</p>
<p>(d) Perform hard decision decoding and count the bit errors</p>
<p>(e) Repeat for multiple values of <img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?\frac{E_b}{N_0}" border="0" alt="" align="absmiddle" /> and plot the simulation and theoretical results.</p>
<p>Click here to download <a title="Matlab/Octave code for simulating BER for BPSK in Rayleigh fading channel with Equal Gain Combining" href="http://www.dsplog.com/db-install/wp-content/uploads/2008/09/script_ber_bpsk_rayleigh_channel_equal_gain_combining.m">Matlab/Octave code for simulating BER for BPSK in Rayleigh fading channel with Equal Gain Combining</a></p>
<p><img class="alignnone size-full wp-image-237" title="BER plot Receive Diversity with Equal Gain Combining" src="http://www.dsplog.com/db-install/wp-content/uploads/2008/09/ber_plot_receive_diversity_equal_gain_combining.png" alt="" width="500" height="397" /></p>
<p><strong>Figure: BER plot Receive Diversity with Equal Gain Combining</strong></p>
<p>Can observe that the simulation results are in good agreement with the theoretical results. <img src='http://www.dsplog.com/db-install/wp-includes/images/smilies/icon_smile.gif' alt=':)' class='wp-smiley' /> </p>
<h2>Reference</h2>
<p><a title="BER with EGC in Rayleigh channel" href="http://ieeexplore.ieee.org/iel1/26/12187/00558680.pdf?tp=&amp;arnumber=558680&amp;isnumber=12187">[ZHANG97] Probability of error for equal-gain combiners over Rayleigh channels: some closed-form solutions</a> Zhang, Q.T. Communications, IEEE Transactions on Volume 45,  Issue 3,  Date: Mar 1997,  Pages: 270 &#8211; 273</p>
<p><a title="Note on Receive diversity by Prof. RaviRaj Adve" href="http://www.comm.utoronto.ca/~rsadve/Notes/DiversityReceive.pdf">Receive diversity &#8211; Notes by Prof. Raviraj Adve</a></p>


<p>Related posts:<ol><li><a href='http://www.dsplog.com/2008/09/28/maximal-ratio-combining/' rel='bookmark' title='Permanent Link: Maximal Ratio Combining (MRC)'>Maximal Ratio Combining (MRC)</a></li><li><a href='http://www.dsplog.com/2008/09/06/receiver-diversity-selection-diversity/' rel='bookmark' title='Permanent Link: Selection Diversity'>Selection Diversity</a></li><li><a href='http://www.dsplog.com/2008/08/19/receive-diversity-in-awgn/' rel='bookmark' title='Permanent Link: Receive diversity in AWGN'>Receive diversity in AWGN</a></li></ol></p>]]></content:encoded>
			<wfw:commentRss>http://www.dsplog.com/2008/09/19/equal-gain-combining/feed/</wfw:commentRss>
		<slash:comments>32</slash:comments>
		</item>
		<item>
		<title>Selection Diversity</title>
		<link>http://www.dsplog.com/2008/09/06/receiver-diversity-selection-diversity/</link>
		<comments>http://www.dsplog.com/2008/09/06/receiver-diversity-selection-diversity/#comments</comments>
		<pubDate>Sat, 06 Sep 2008 09:44:11 +0000</pubDate>
		<dc:creator>Krishna Sankar</dc:creator>
				<category><![CDATA[Diversity]]></category>
		<category><![CDATA[Error Rate]]></category>
		<category><![CDATA[BPSK]]></category>
		<category><![CDATA[PSK]]></category>
		<category><![CDATA[Rayleigh]]></category>

		<guid isPermaLink="false">http://www.dsplog.com/?p=228</guid>
		<description><![CDATA[This is the first post in the series discussing receiver diversity in a wireless link. Receiver diversity is a form of space diversity, where there are multiple antennas at the receiver. The presence of receiver diversity poses an interesting problem &#8211; how do we use &#8216;effectively&#8216; the information from all the antennas to demodulate the [...]


Related posts:<ol><li><a href='http://www.dsplog.com/2008/09/28/maximal-ratio-combining/' rel='bookmark' title='Permanent Link: Maximal Ratio Combining (MRC)'>Maximal Ratio Combining (MRC)</a></li><li><a href='http://www.dsplog.com/2008/08/19/receive-diversity-in-awgn/' rel='bookmark' title='Permanent Link: Receive diversity in AWGN'>Receive diversity in AWGN</a></li><li><a href='http://www.dsplog.com/2008/09/19/equal-gain-combining/' rel='bookmark' title='Permanent Link: Equal Gain Combining (EGC)'>Equal Gain Combining (EGC)</a></li></ol>]]></description>
			<content:encoded><![CDATA[<p></p><p>This is the first post in the series discussing <strong>receiver diversity</strong> in a wireless link. <strong>Receiver diversity</strong> is a form of space diversity, where there are multiple antennas at the receiver. The presence of receiver diversity poses an interesting problem &#8211; how do we use &#8216;<em>effectively</em>&#8216; the information from all the antennas to demodulate the data. There are multiple ways to approach the problem. The three typical approaches to be discussed are &#8211; <strong>selection diversity</strong>, <strong>equal gain combining</strong> and <strong>maximal ratio combining</strong>. In this post we will discuss <strong>selection diversity</strong>. For the discussion, we will assume that the channel is a flat fading <a title="Discussion on Rayleigh multipath channel" href="http://www.dsplog.com/2008/07/14/rayleigh-multipath-channel/" target="_self">Rayleigh multipath channel</a> and the modulation is BPSK.</p>
<p><span id="more-228"></span></p>
<h2>Background</h2>
<p>1. We have N receive antennas and one transmit antenna.</p>
<p>2. The channel is flat fading &#8211; 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 <a title="Digital Communications: Fundamentals and Applications (2nd Edition) - Bernard Sklar (Author)" href="http://www.amazon.com/gp/redirect.html?ie=UTF8&amp;location=http%3A%2F%2Fwww.amazon.com%2FDigital-Communications-Fundamentals-Applications-Technologies%2Fdp%2F0130847887&amp;tag=dl04-20&amp;linkCode=ur2&amp;camp=1789&amp;creative=9325">[DIGITAL COMMUNICATIONS: SKLAR]</a><img style="border: medium none  ! important; margin: 0px ! important;" src="http://www.assoc-amazon.com/e/ir?t=dl04-20&amp;l=ur2&amp;o=1" border="0" alt="" width="1" height="1" /></p>
<p>3. The channel experienced by each receive antenna is randomly varying in time. For the <img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?i^{th}" border="0" alt="" align="absmiddle" /> receive antenna, each transmitted symbol gets multiplied by a randomly varying complex number <img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?h_i" border="0" alt="" align="absmiddle" />. As the channel under consideration is a Rayleigh channel, the real and imaginary parts of <img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?h_i" border="0" alt="" align="absmiddle" /> are Gaussian distributed having mean 0 and variance 1/2.</p>
<p>4. The channel experience by each receive antenna is independent from the channel experienced by other receive antennas.</p>
<p>5. On each receive antenna, the noise<img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?n" border="0" alt="" align="absmiddle" /> has the Gaussian probability density function with</p>
<p><img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?p%28n%29%20=%20%5Cfrac%7B1%7D%7B%5Csqrt%7B2%5Cpi%5Csigma%5E2%7D%7De%5E%7B%5Cfrac%7B-%28n-%5Cmu%29%5E2%7D%7B2%5Csigma%5E2%7D" border="0" alt="" align="absmiddle" /> with <img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?%5Cmu=0" border="0" alt="" align="absmiddle" /> and <img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?%5Csigma%5E2%20=%20%5Cfrac%7BN_0%7D%7B2%7D" border="0" alt="" align="absmiddle" />.</p>
<p>The noise on each receive antenna is independent from the noise on the other receive antennas.</p>
<p>6. At each receive antenna, the channel <img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?h_i" border="0" alt="" align="absmiddle" /> is known at the receiver. For example, on the <img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?i^{th}" border="0" alt="" align="absmiddle" /> receive antenna, equalization is performed at the receiver by dividing the received symbol <img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?y_i" border="0" alt="" align="absmiddle" /> by the apriori known <img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?h_i" border="0" alt="" align="absmiddle" /> i.e.</p>
<p><img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?\hat{y} = \frac{y_i}{h_i} = \frac{h_ix+n_i}{h_i} = x+\tilde{n_i}" border="0" alt="" align="absmiddle" /></p>
<p>where<br />
<img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?\tilde{n_i} = \frac{n_i}{h_i}" border="0" alt="" align="absmiddle" />is the additive noise scaled by the channel coefficient.</p>
<p>7. In the presence of channel <img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?h_i" border="0" alt="" align="absmiddle" />, the instantaneous bit energy to noise ratio at <img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?i^{th}" border="0" alt="" align="absmiddle" /> receive antenna is <img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?\frac{|h_i|^2E_b}{N_0}" border="0" alt="" align="absmiddle" />. For notational convenience, let us define,</p>
<p><img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?\gamma_i = \frac{|h_i|^2E_b}{N_0}" border="0" alt="" align="absmiddle" />.</p>
<p>From the discussion on <a title="Chi-square random variable with two degrees of freedom" href="http://www.dsplog.com/2008/07/28/chi-square-random-variable/#twodegree">chi-square random variable</a>, we know that, if <img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?h_i" border="0" alt="" align="absmiddle" /> is a Rayleigh distributed random variable, then <img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?h_i^2" border="0" alt="" align="absmiddle" /> is a chi-squared random variable with two degrees of freedom. The pdf of <img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?\gamma_i" border="0" alt="" align="absmiddle" /> is<br />
<img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?p(\gamma_i)=\frac{1}{(E_b/N_0)}e^{\frac{-\gamma_i}{(E_b/N_0)}" border="0" alt="" align="absmiddle" />.</p>
<h2>What is selection diversity?</h2>
<p>Consider a scenario where we have a single antenna for transmission and multiple antennas at the receiver (as shown in the figure below).</p>
<p><img class="alignnone size-full wp-image-229" title="Receive Diversity" src="http://www.dsplog.com/db-install/wp-content/uploads/2008/09/receive_diversity.png" alt="" width="415" height="355" /></p>
<p><strong>Figure: Receive diversity in a wireless link</strong></p>
<p>At the receiver we have now N copies of the same transmitted symbol. Which then poses the problem &#8211; how to effectively combine them to reliably recover the data.</p>
<p><strong>Selection diversity</strong> approach is one way out &#8211; With <strong>selection diversity</strong>, the receiver selects the antenna with the highest received signal power and ignore observations from the other antennas. The chosen receive antenna is one which gives <img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?\max_i (\gamma_i)" border="0" alt="" align="absmiddle" />.</p>
<h2>Outage probability in Selection Diversity</h2>
<p>The equations in the post refers the note on <a title="Note on Receive diversity by Prof. RaviRaj Adve" href="http://www.comm.utoronto.ca/~rsadve/Notes/DiversityReceive.pdf">Receive diversity</a> by Prof. RaviRaj Adve.</p>
<p>To analyze the bit error rate, let us first find the outage probability on the <img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?i^{th}" border="0" alt="" align="absmiddle" /> receive antenna. Outage probability is the probability that the bit energy to noise ratio falls below a threshold. The probability of outage on <img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?i^{th}" border="0" alt="" align="absmiddle" /> receive antenna is,</p>
<p><img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?\begin{eqnarray}P_{out,\gamma_i}&amp;=&amp;P[\gamma_i &lt; \gamma_s]\\&amp; = &amp;\int_0^{\gamma_s} \frac{1}{(E_b/N_0)}e^{\frac{-\gamma_i}{(E_b/N_0)}}d\gamma_i  \\&amp;=&amp; 1-e^{-\frac{\gamma_s}{(E_b/N_0)}}\end{eqnarray}" border="0" alt="" align="absmiddle" />.</p>
<p><img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?\gamma_s" border="0" alt="" align="absmiddle" /> is the defined threshold for bit energy to noise ratio.</p>
<p>In N reveive antenna case, the probability that all bit energy to noise ratio on all the receive antenna are below the threshold <img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?\gamma_s" border="0" alt="" align="absmiddle" />is,</p>
<p><img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?P_{out}=P[\gamma_1, \gamma_2,\cdots,\gamma_N &lt; \gamma_s]" border="0" alt="" align="absmiddle" />,</p>
<p>where</p>
<p><img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?\gamma_1, \gamma_2,\cdots,\gamma_N" border="0" alt="" align="absmiddle" /> are the bit energy to noise ratio on the 1st, 2nd and so on till the Nth receive antenna.</p>
<p>Since the channel on each antenna is assumed to independent, the joint probability is the product of individual probabilities.</p>
<p><img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?\begin{eqnarray}P_{out}&amp;=&amp;P[\gamma_1 &lt; \gamma_s]P[\gamma_2 &lt; \gamma_s]\cdots P[\gamma_N &lt; \gamma_s]\\&lt;br /&gt; &amp; = &amp;\prod_{i=1}^{N}P[\gamma_i &lt; \gamma_s]\\&lt;br /&gt; &amp; = &amp; \left[1-e^{-\frac{\gamma_s}{(E_b/N_0)}} \right]^N\end{eqnarray}" border="0" alt="" align="absmiddle" />.</p>
<p>Note that the equation above defines the probability that the effective bit energy to noise ratio with N receive antennas (lets call <img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?\gamma" border="0" alt="" align="absmiddle" />) is lower than the threshold <img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?\gamma_s" border="0" alt="" align="absmiddle" />.  This is infact the cumulative distribuition function (CDF) of <img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?\gamma" border="0" alt="" align="absmiddle" />. The probability density function (PDF) is then the deriviate of the CDF. <img src='http://www.dsplog.com/db-install/wp-includes/images/smilies/icon_smile.gif' alt=':)' class='wp-smiley' /> </p>
<p><img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?\begin{eqnarray}p(\gamma) &amp;=&amp; \frac{dP_{out}}{d\gamma}\\&amp;= &amp;\frac{N}{(E_b/N_0)}e^{-\frac{\gamma}{(E_b/N_0)}}\left[1-e^{-\frac{\gamma}{(E_b/N_0)}} \right]^{N-1}\\\end{eqnarray}" border="0" alt="" align="absmiddle" />.</p>
<p>Given that we know the PDF of <img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?\gamma" border="0" alt="" align="absmiddle" />, the average output bit energy to noise ratio is,</p>
<p><img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?\begin{eqnarray}E(\gamma) &amp; = &amp; \int_0^{\infty}\gamma p(\gamma)d\gamma \\&lt;br /&gt; &amp; = &amp; \int_0^{\infty}\gamma \frac{N}{(E_b/N_0)}e^{-\frac{\gamma}{(E_b/N_0)}}\left[1-e^{-\frac{\gamma}{(E_b/N_0)}} \right]^{N-1} \\&lt;br /&gt; &amp; = &amp; \frac{E_b}{N_0}\sum_{i=1}^N\frac{1}{i}\end{eqnarray}" border="0" alt="" align="absmiddle" />.</p>
<address>I do not know how to reduce the above integral to this simple sum. <img src='http://www.dsplog.com/db-install/wp-includes/images/smilies/icon_sad.gif' alt=':(' class='wp-smiley' /> </address>
<p>This means that,</p>
<p>- with two receive antennas the effective bit energy to noise ratio is 1.5 times <img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?\frac{E_b}{N_0}" border="0" alt="" align="absmiddle" />,</p>
<p>- with three receive antennas, the effective bit energy to noise ratio is 1.833 times <img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?\frac{E_b}{N_0}" border="0" alt="" align="absmiddle" />,</p>
<p>- with four receive antennas, the effective bit energy to noise ratio is 2 times <img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?\frac{E_b}{N_0}" border="0" alt="" align="absmiddle" /> and so on.</p>
<p>If you recall the results from the <a title="BER analysis in AWGN with receive diversity" href="http://www.dsplog.com/2008/08/19/receive-diversity-in-awgn/">AWGN with receive diversity case</a>,</p>
<blockquote><p>Effective bit energy to noise ratio in a N receive antenna case is N times the bit energy to noise ratio for single antenna case.</p></blockquote>
<p>With <strong>selection diversity</strong> we are seeing that the effective SNR improvement is <strong>not a linear improvement</strong> with increasing the number of receive antennas. The returns diminish.</p>
<p><img class="alignnone size-full wp-image-232" title="SNR gain with selection diversity" src="http://www.dsplog.com/db-install/wp-content/uploads/2008/09/effective_snr_selection_diversity.png" alt="" width="500" height="399" /></p>
<p><strong>Figure: SNR gain with selection diversity</strong></p>
<p>Click here to download<a href="http://www.dsplog.com/db-install/wp-content/uploads/2008/09/script_selection_diversity_effective_snr.m"><span style="text-decoration: none; color: #000000;"> </span></a><a title="Matlab/Octave script for computing the effective SNR in Ralyeigh channel with selection diversity" href="http://www.dsplog.com/db-install/wp-content/uploads/2008/09/script_selection_diversity_effective_snr.m">Matlab/Octave script for computing the effective SNR in Ralyeigh channel with selection diversity</a></p>
<h2>Bit Error probability with selection diversity</h2>
<p>If you recall, in the post on <a title="BER for BPSK in AWGN" href="http://www.dsplog.com/2007/08/05/bit-error-probability-for-bpsk-modulation/">BER computation in AWGN</a>, with bit energy to noise ratio of <img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?\frac{E_b}{N_0}" border="0" alt="" align="absmiddle" />, the bit error rate for BPSK in AWGN is derived as</p>
<p><img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?P_b=%5Cfrac%7B1%7D%7B2%7Derfc%5Cleft%28%7B%5Csqrt%7B%5Cfrac%7BE_b%7D%7BN_0%7D%7D%7D%5Cright%29" border="0" alt="" align="absmiddle" />.</p>
<p>Given that the effective bit energy to noise ratio with selection diversity is<img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?\gamma" border="0" alt="" align="absmiddle" />, the total bit error rate is the integral of the conditional BER integrated over all possible values of <img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?\gamma" border="0" alt="" align="absmiddle" />.</p>
<p><img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?\begin{eqnarray}P_{e}&amp;=&amp;\int_0^{\infty}\frac{1}{2}erfc\left(\sqrt{\gamma}\right)p(\gamma)d\gamma\\&lt;br /&gt; &amp;=&amp;\int_0^{\infty}\frac{1}{2}erfc\left(\sqrt{\gamma}\right)\frac{N}{(E_b/N_0)}e^{-\frac{\gamma}{(E_b/N_0)}}\left[1-e^{-\frac{\gamma}{(E_b/N_0)}} \right]^{N-1}d\gamma\end{eqnarray}" border="0" alt="" align="absmiddle" />.</p>
<p>This equation reduces to</p>
<p><img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?\begin{eqnarray}P_{e}&amp;=&amp;\frac{1}{2}\sum_{k=0}^{N}(-1)^k\left(\begin{eqnarray}N\\k\end{eqnarray}\right)\left(1+\frac{k}{(E_b/N_0)}\right)^{-1/2}\end{eqnarray}" border="0" alt="" align="absmiddle" />.</p>
<address>Refer Equation 11.24 in Section 11.3.2 Performance with Selection combining in <a title="Digital Communication: Third Edition, by John R. Barry, Edward A. Lee, David G. Messerschmitt" href="http://www.amazon.com/gp/redirect.html?ie=UTF8&amp;location=http%3A%2F%2Fwww.amazon.com%2FDigital-Communication-John-R-Barry%2Fdp%2F0792375483&amp;tag=dl04-20&amp;linkCode=ur2&amp;camp=1789&amp;creative=9325">[DIG-COMM-BARRY-LEE-MESSERSCHMITT].</a> Again, I do not know the proof <img src='http://www.dsplog.com/db-install/wp-includes/images/smilies/icon_sad.gif' alt=':(' class='wp-smiley' /> </address>
<p><strong><br />
</strong></p>
<h2>BER Simulation Model</h2>
<p>The Matlab/Octave script performs the following</p>
<p>(a) Generate random binary sequence of +1&#8217;s and -1&#8217;s.</p>
<p>(b) Multiply the symbols with the channel and then add white Gaussian noise.</p>
<p>(c) At the receiver, find the receive path with maximum power</p>
<p>(d) Chose that receive path, equalize (divide) the received symbols with the known channel</p>
<p>(d) Perform hard decision decoding and count the bit errors</p>
<p>(e) Repeat for multiple values of <img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?\frac{E_b}{N_0}" border="0" alt="" align="absmiddle" /> and plot the simulation and theoretical results.</p>
<p>Click here to download <a href="http://www.dsplog.com/db-install/wp-content/uploads/2008/09/script_ber_bpsk_rayleigh_channel_selection_diversity.m">Matlab/Octave script for simulating BER for BPSK in Rayleigh channel with selection diversity</a></p>
<p><img class="alignnone size-full wp-image-233" title="BER plot for BPSK in Rayleigh channel with Selection Diversity" src="http://www.dsplog.com/db-install/wp-content/uploads/2008/09/ber_plot_bpsk_rayleigh_channel_selection_diversity.png" alt="" width="500" height="386" /></p>
<p><strong>Figure: BER plot for BPSK in Rayleigh channel with Selection Diversity</strong></p>
<h2>Observations</h2>
<p>Around 16dB improvement at <img src="http://www.dsplog.com/cgi-bin/mimetex.cgi?10^{-4}" border="0" alt="" align="absmiddle" /> BER point by with two receive antenna selection diversity <img src='http://www.dsplog.com/db-install/wp-includes/images/smilies/icon_smile.gif' alt=':)' class='wp-smiley' /> </p>
<h2>References</h2>
<p><a title="Note on Receive diversity by Prof. RaviRaj Adve" href="http://www.comm.utoronto.ca/~rsadve/Notes/DiversityReceive.pdf">Receive diversity &#8211; Notes by Prof. Raviraj Adve</a></p>
<p><a title="Digital Communication: Third Edition, by John R. Barry, Edward A. Lee, David G. Messerschmitt" href="http://www.amazon.com/gp/redirect.html?ie=UTF8&amp;location=http%3A%2F%2Fwww.amazon.com%2FDigital-Communication-John-R-Barry%2Fdp%2F0792375483&amp;tag=dl04-20&amp;linkCode=ur2&amp;camp=1789&amp;creative=9325">[DIG-COMM-BARRY-LEE-MESSERSCHMITT] </a><a title="Digital Communication: Third Edition, by John R. Barry, Edward A. Lee, David G. Messerschmitt" href="http://www.amazon.com/gp/redirect.html?ie=UTF8&amp;location=http%3A%2F%2Fwww.amazon.com%2FDigital-Communication-John-R-Barry%2Fdp%2F0792375483&amp;tag=dl04-20&amp;linkCode=ur2&amp;camp=1789&amp;creative=9325">Digital Communication: Third Edition, by John R. Barry, Edward A. Lee, David G. Messerschmitt</a><img style="border-style: none ! important; margin: 0px" src="http://www.assoc-amazon.com/e/ir?t=dl04-20&amp;l=ur2&amp;o=1" border="0" alt="" width="1" height="1" /></p>
<p>Happy learning.</p>


<p>Related posts:<ol><li><a href='http://www.dsplog.com/2008/09/28/maximal-ratio-combining/' rel='bookmark' title='Permanent Link: Maximal Ratio Combining (MRC)'>Maximal Ratio Combining (MRC)</a></li><li><a href='http://www.dsplog.com/2008/08/19/receive-diversity-in-awgn/' rel='bookmark' title='Permanent Link: Receive diversity in AWGN'>Receive diversity in AWGN</a></li><li><a href='http://www.dsplog.com/2008/09/19/equal-gain-combining/' rel='bookmark' title='Permanent Link: Equal Gain Combining (EGC)'>Equal Gain Combining (EGC)</a></li></ol></p>]]></content:encoded>
			<wfw:commentRss>http://www.dsplog.com/2008/09/06/receiver-diversity-selection-diversity/feed/</wfw:commentRss>
		<slash:comments>38</slash:comments>
		</item>
	</channel>
</rss>
