Maximum Ratio Combining (MRC) architecture simulation - GaussianWaves (2024)

In the previous post on Single Input Multiple Output (SIMO) models for receive diversity, various receiver diversity techniques were outlined. One of them is maximum ratio combining, the focus of the topic here.

Channel model

Assuming flat slow fading channel, the received signal model is given by

where, is the channel impulse response, is the received signal, is the transmitted signal and is the additive Gaussian white noise.

Assuming small scale Rayleigh fading, the channel impulse response is modeled as complex Gaussian random variable with zero mean and variance

Therefore, the instantaneous channel power is exponentially distributed

In the context of AWGN channel, the signal-to-noise ratio (SNR) for a given channel condition, is a constant. But in the case of fading channels, the signal-to-noise ratio is no longer a constant as the signal is fluctuating when passed through a fading channel. Therefore, for fading channel, the SNR has a random variable component built into it. Hence, we just don’t call it SNR, instead it is called instantaneous SNR which depends on the current conditions of the channel (or equivalently, the value of the random variable at that instant). Since the SNR is a random variable, we can also talk about its expected (average) value, which is called average SNR. Denoting the average SNR as and for convenience, let’s assume that the average power of the channel is unity, i.e,

The instantaneous SNR is given by

Therefore, like the channel impulse response, the instantaneous SNR is also exponentially distributed

Maximum Ratio Combining (MRC)

The selection combining technique is the simplest technique, where in, the received signal from the antenna that experiences the highest SNR (i.e, the strongest signal from N received signals) is chosen for processing at the receiver. Therefore this technique throws away of observations. Whereas, in maximum ratio combining (MRC) all observations are used.

MRC works on the signal in spatial domain and is very similar to what a matched filter in frequency domain does to the incoming signal. MRC maximizes the inner product of the weights and the signal vector.

The maximum ratio combining technique, uses all the received signal elements (Figure 1), it weighs them and combines the weighted signals so that the output SNR is maximized. Requiring the knowledge of the individual channels , the weights are chosen as

With the weights set as , the output of the MRC combiner is

Therefore, the output SNR after MRC processing is

MRC processing results in the weighted average of the received signals and hence the overall output SNR is equal to the sum of the SNRs of all individual receive signals, which yields the maximum possible diversity gain of . This is the maximum achievable SNR for all possible receive diversity schemes (selection combining, equal gain combining, etc..,).

Generally, two figures of merits are used to gauge the performance of the diversity schemes – outage probability and error rate performance for PSK modulation.

Outage probability

As we know, fading channels are characterized by deep fades, i.e, the period when the signal level falls below a certain threshold or certain noise level. During such fades, the user experiences signal outage. We would like to compute the probability, in certain fading channel, that a user will experience signal outage. This is called outage probability. Outage probability can be easily computed if we know the probability distribution characteristics of the fading.

The outage probability with which the instantaneous output SNR of MRC falls below a given SNR target is

For high average SNR conditions , the outage probability can be approximated as

Python code

import numpy as npimport matplotlib.pyplot as pltfrom scipy.special import factorialgamma_ratio_dB = np.arange(start=-10,stop=40,step=2)Ns = [1,2,3,4,10,20] #number of received signal pathsgamma_ratio = 10**(gamma_ratio_dB/10) #Average SNR/SNR threshold in dBfig, ax = plt.subplots(1, 1)for N in Ns: n = np.arange(start=0,stop=N,step=1) P_outage = 1 - np.exp(-1/gamma_ratio)*np.sum(((1/gamma_ratio)**n[:,None])/factorial(n[:,None]),axis=0) ax.semilogy(gamma_ratio_dB,P_outage,label='N='+str(N))ax.legend()ax.set_xlim(-10,40);ax.set_ylim(0.0001,1.1)ax.set_title('MRC outage probability (Rayleigh fading channel)')ax.set_xlabel(r'$10log_{10}\left(\Gamma/\gamma_t\right)$')ax.set_ylabel('Outage probability');fig.show()

Figure 2, plots the outage probability against (the ratio of average SNR and the SNR threshold) for different values – the number of received signals received over an Rayleigh flat fading channel. For example, the outage probability dramatically improves when going from branch to branches. At outage probability of 0.01% (projected y-value in the graph at ), there is an approximate reduction in the required SNR.

Error rate performance

In the case of receive diversity schemes with antennas, the received signal vector is given by

Considering the QPSK modulated symbols that are transmitted (denoted as ), the maximum likelihood detection criterion for detecting the transmitted symbols by the equalizer block at the receiver is given by,

The solution to this problem can be obtained using the least squares method (refer equation 8.17 given in chapter 8 of this book)

The solution can be re-written as

As an example, the symbol error rate performance of a QPSK modulated transmission over a Rayleigh flat fading SIMO channel, for a range of values for the number of receive antennas () is simulated here. Maximum ratio combining is used in the receiver.

The code utilizes the modem class discussed in the book here. The modem class incorporates modulation and demodulation techniques for PSK,PAM,QAM and FSK modulation schemes. It uses the object oriented programming method for implementing the various modems.

The addition of Gaussian white noise needs to be multidimensional. The method discussed in this article is extended here for computing and adding the required amount of noise across the branches of signals.

"""Eb/N0 Vs SER for PSK over Rayleigh flat fading with MRC@author: Mathuranathan ViswanathanCreated on Jan 16, 2020"""import numpy as np # for numerical computingimport matplotlib.pyplot as plt # for plotting functions#from matplotlib import cm # colormap for color palettefrom numpy.random import standard_normalfrom DigiCommPy.modem import PSKModemfrom DigiCommPy.channels import awgn#---------Input Fields------------------------nSym = 10**6 # Number of symbols to transmitEbN0dBs = np.arange(start=-20,stop = 36, step = 2) # Eb/N0 range in dB for simulationN = [1,2,4,8,10] # [1,2,3,4,10,20] #number of diversity branchesM = 4 #M-ary PSK modulationk=np.log2(M)EsN0dBs = 10*np.log10(k)+EbN0dBs # EsN0dB calculationfig, ax = plt.subplots(nrows=1,ncols = 1) #To plot figurefor nRx in N: #simulate for each # of received branchs #Random input symbols to modulator inputSyms = np.random.randint(low=0, high = M, size=nSym) modem = PSKModem(M) s = modem.modulate(inputSyms) #modulated PSK symbols #nRx signal branches s_diversity = np.kron(np.ones((nRx,1)),s); ser_sim = np.zeros(len(EbN0dBs)) # simulated symbol error rates for i,EsN0dB in enumerate(EsN0dBs): #Rayleigh flat fading channel as channel matrix h = np.sqrt(1/2)*(standard_normal((nRx,nSym))+1j*standard_normal((nRx,nSym))) signal = h*s_diversity #effect of channel on the modulated signal #Computing the signal power and adding noise gamma = 10**(EsN0dB/10) #converting EsN0dB to linear scale P = np.sum(np.abs(signal)**2,axis=1)/nSym #calculate power in each branch of signal N0 = P/gamma #required noise spectral density for each branch #Scale each row of noise with the calculated noise spectral density noise = (standard_normal(signal.shape)+1j*standard_normal(signal.shape))*np.sqrt(N0/2)[:,None] r = signal+noise #received signal branches #Receiver processing equalized = np.sum(r*np.conj(h),axis=0) #equalized signal detectedSyms = modem.demodulate(equalized) #demodulation decisions ser_sim[i] = np.sum(detectedSyms != inputSyms)/nSym #ax.grid(True,which='both'); ax.semilogy(EbN0dBs,ser_sim,label='N='+str(nRx))#plot simulated error ratesax.set_xlim(-20,35);ax.set_ylim(0.0001,1.1);ax.grid(True,which='both');ax.set_xlabel('Eb/N0(dB)');ax.set_ylabel('Symbol Error Rate($P_s$)')ax.set_title('SER performance for QPSK over Rayleigh fading channel with MRC')ax.legend();fig.show()

Rate this article: (12 votes, average: 3.92 out of 5)

References

[1] Andrea Goldsmith, “Wireless Communications”, ISBN: 978-0521837163, Cambridge University Press, 1 edition.↗

[2] Barry-Lee-Messerschmitt, “Digital Communication”, ISBN: 978-0792375487 , Springer, 3rd edition, September 30, 2003.↗

Articles in this series

Books by the author