Skip to main content

5G Channel Estimation using Orthogonal Matching Pursuit (OMP)



For millimeter wave massive MIMO communication in 5G, we observe that the number of available multipath that avails communication is much smaller than the maximum connections possible between the transmitter(TX) and receiver(RX). Only a few MPCs reach at receiver with good received signal strength. For example, the number of strong MPCs that reaches the receiver is L and there is N transmitter antenna on the transmitter side and N number of antennas on the receiver side.

So, from the channel matrix of the massive MIMO system, we can say the total number of available paths or connections between TX and RX is equal to, N X N or, N^(2)

Now, L << N^(2)

For simplicity, if the number of possible strong beams from the transmitter and receiver sides are NtBeams and NrBeams, then, L = NtBeams * NrBeams

If we look up the massive MIMO channel matrix, then, H=



Primarily, if the number of available MPCs to avail communication between TX and RX are much lesser than the maximum connections possible, N^(2), then the channel matrix will be sparse. Because MPCs other than stronger L number of paths are so weak that they cannot communicate between TX and RX.

Sparse channel matrix means where the elements are zeros. We have nothing to do with those elements at all. So, we need such an algorithm that can sense the zeros of the channel matrix and we can transmit our signal thru that selected first L number of stronger paths. The orthogonal matching Pursuit (OMP) algorithm helps us to do that. It can sense a sparse matrix. OMP is sometimes termed as the "Compressive Sensing Method".


Orthogonal Matching Pursuit Algorithm:

Let assume, NtBeams = NrBeams = 28

Total stronger connections or paths, L = 28 * 28 = 784 (it is observational)

Also assume, the number of antennas at both the transmitter and receiver side is equal to 32

So, the channel matrix contains 32 * 32 elements or 1024 elements

Now, the Equivalent sensing matrix, Q (say), will have an array dimension of 784 X 1024


Mathematically,

y = Qhb + n

Here, y = received vector

Q = Equivalent sensing matrix

hb = beamspace channel matrix where the number of elements is N X N = N^(2). It is a sparse matrix as there are only L-acceptable elements out of N^(2)

n = noise


Also assume, Matrix Q contains n number of columns, like that,

Q = [q1 q2 q3 ..... qn]

Now, we focus on finding out the maximum in the context of which column with a projectile with y generates maximum value, such that,

i(1) = argmax |(qjH / ||qj||)*y| where, j=1,2,...,n

Let's assume, qi1 column in Q matrix contributes maximum value. Then find out using the least square solution

min ||y - qi1*hb1|| (We try to make it zero)

Where, hb1 is any column of hb

Where, hb1= (qi1H * qi1)^(-1) * qi1H *y


Now, take a residue matrix as,

r1 = y - qi1*hb1

Now, we'll find the maximum correlation of residue matrix r1 by finding out which column of Q with a projectile with r1 generates maximum value, such that,

i(2) = argmax|(qjH / ||qj||)*r1| where, j=1,2,...,n

Let assume, qi2 column in Q generates the maximum value

Now form a matrix, Qn, like that,

Qn = [qi1 qi2]

Now find out hb2

hb2 = (QnH * Qn)^(-1) * QnH *y

Now update the residue matrix as,

r2 = y - Qn*hb2

If norm of ||r(n-1) - rn|| falls below a threshold value, then close the loop

In our case, if we consider r2-r1 falls below the threshold, then we close the loop.

If we find hb2=[3;2] then we do such operations in hb matrix, such that,

and if i(1) = 5 and i(2) = 2

Then hb will be [0;2;0;0;3;0; ... up to 28th row all are 0]

Here, we place the element value of hb2 in hb matrix in such row with row number matches with the value of i(1) and i(2)


Mathematical Example of Orthogonal Matching Pursuit (OMP):

Let's assume, for a MIMO communication system,

The size of the equivalent sensing matrix, Q is 4 X 6

And received signal matrix, y=




Now, y = Qhb

Or,

Where, hb =beamspace matrix =



Let assume, Q = [ q1 q2 q3 q4 q5 q6]

Here, q1 is the first column of Q, q2 is the second column of Q, and so on.


First Iteration of Orthogonal Matching Pursuit:

Now we find the maximum correlation of y by finding out which column in Q generates the maximum value with the projection of y,

Or,

QT*y =



Here, we can see the element in the 5th row is the maximum among all elements. So, we’ll select the 5th column of Q with which y has the maximum correlation value.

Now, hb1 = (q5'*q5)^(-1) * q5'*y

Where q5’ denotes the transpose of q5

Or,

hb1 = 4

Residue Matrix, r1 = y – q5* hb1

Or, r1 =


Here, we observe residual matrix r1 is not a zero matrix. So, we go for 2nd iteration.

Second Iteration of Orthogonal Matching Pursuit

Where we find the maximum correlation of r1 with respect to Q matrix.

Alternatively,

QT*r1=


Now, we see the element in the 2nd row of the above matrix generates the maximum value so r1 has a maximum correlation with the 2nd column of Q Matrix.

Now, we’ll form a new matrix, Qn = [q5 q2]

We find hb2 = (Qn'*Qn)^(-1)*Qn'*y;

Or, hb2 =


Now updated residue matrix, r2 = y – Qn* hb2

Or, r2=


Now we get the desired value in residue matrix r2 where all elements are zeros. So, beamspace matrix, hb will be


Here, we replace the elemental value of hb in that rows which are equal with the number of columns which generates maximum values with projection with y, then r1, and so on.

Now, from the mathematics we can say,

y = Q* hb

Or,



<< Back to Previous Page


Also read about

[1] 5G : Channel modelling for millimeter wave

[2] Time-delayed saleh valenzuala cluster model for UWB & mm-Wave

#beamforming

People are good at skipping over material they already know!

View Related Topics to







Admin & Author: Salim

s

  Website: www.salimwireless.com
  Interests: Signal Processing, Telecommunication, 5G Technology, Present & Future Wireless Technologies, Digital Signal Processing, Computer Networks, Millimeter Wave Band Channel, Web Development
  Seeking an opportunity in the Teaching or Electronics & Telecommunication domains.
  Possess M.Tech in Electronic Communication Systems.


Contact Us

Name

Email *

Message *

Popular Posts

BER vs SNR for M-ary QAM, M-ary PSK, QPSK, BPSK, ...

๐Ÿ“˜ Overview of BER and SNR ๐Ÿงฎ Online Simulator for BER calculation of m-ary QAM and m-ary PSK ๐Ÿงฎ MATLAB Code for BER calculation of M-ary QAM, M-ary PSK, QPSK, BPSK, ... ๐Ÿ“š Further Reading ๐Ÿ“‚ View Other Topics on M-ary QAM, M-ary PSK, QPSK ... ๐Ÿงฎ Online Simulator for Constellation Diagram of m-ary QAM ๐Ÿงฎ Online Simulator for Constellation Diagram of m-ary PSK ๐Ÿงฎ MATLAB Code for BER calculation of ASK, FSK, and PSK ๐Ÿงฎ MATLAB Code for BER calculation of Alamouti Scheme ๐Ÿงฎ Different approaches to calculate BER vs SNR What is Bit Error Rate (BER)? The abbreviation BER stands for bit error rate, which indicates how many corrupted bits are received (after the demodulation process) compared to the total number of bits sent in a communication process. It is defined as,  In mathematics, BER = (number of bits received in error / total number of transmitted bits)  On the other hand, SNR ...

RMS Delay Spread, Excess Delay Spread and Multi-path ...

๐Ÿ“˜ Overview of Delay Spread and Multi-path ๐Ÿงฎ Excess Delay spread ๐Ÿงฎ Power delay Profile ๐Ÿงฎ RMS Delay Spread ๐Ÿ“š Further Reading ๐Ÿ“‚ Other Topics on RMS Delay Spread, Excess Delay ... ๐Ÿงฎ Multipath Components or MPCs ๐Ÿงฎ Online Simulator for Calculating RMS Delay Spread ๐Ÿงฎ Why is there significant multipath in the case of very high frequencies? ๐Ÿงฎ Why RMS Delay Spread is essential for wireless communication? ๐Ÿงฎ Why the Power Delay Profile is essential? ๐Ÿงฎ MATLAB Codes for Calculating Different Types of delay Spreads Delay Spread, Excess Delay Spread, and Multipath (MPCs) The fundamental distinction between wireless and wired connections is that in wireless connections signal reaches at receiver thru multipath signal propagation rather than directed transmission like co-axial cable. Wireless Communication has no set communication path between the transmitter and the receiver. The line...

Constellation Diagrams of ASK, PSK, and FSK

๐Ÿ“˜ Overview of Energy per Bit (Eb / N0) ๐Ÿงฎ Online Simulator for constellation diagrams of ASK, FSK, and PSK ๐Ÿงฎ Theory behind Constellation Diagrams of ASK, FSK, and PSK ๐Ÿงฎ MATLAB Codes for Constellation Diagrams of ASK, FSK, and PSK ๐Ÿ“š Further Reading ๐Ÿ“‚ Other Topics on Constellation Diagrams of ASK, PSK, and FSK ... ๐Ÿงฎ Simulator for constellation diagrams of m-ary PSK ๐Ÿงฎ Simulator for constellation diagrams of m-ary QAM BASK (Binary ASK) Modulation: Transmits one of two signals: 0 or -√Eb, where Eb​ is the energy per bit. These signals represent binary 0 and 1.    BFSK (Binary FSK) Modulation: Transmits one of two signals: +√Eb​ ( On the y-axis, the phase shift of 90 degrees with respect to the x-axis, which is also termed phase offset ) or √Eb (on x-axis), where Eb​ is the energy per bit. These signals represent binary 0 and 1.  BPSK (Binary PSK) Modulation: Transmits one of two signals...

Comparisons among ASK, PSK, and FSK | And the definitions of each

๐Ÿ“˜ Comparisons among ASK, FSK, and PSK ๐Ÿงฎ Online Simulator for calculating Bandwidth of ASK, FSK, and PSK ๐Ÿงฎ MATLAB Code for BER vs. SNR Analysis of ASK, FSK, and PSK ๐Ÿ“š Further Reading ๐Ÿ“‚ View Other Topics on Comparisons among ASK, PSK, and FSK ... ๐Ÿงฎ Comparisons of Noise Sensitivity, Bandwidth, Complexity, etc. ๐Ÿงฎ MATLAB Code for Constellation Diagrams of ASK, FSK, and PSK ๐Ÿงฎ Online Simulator for ASK, FSK, and PSK Generation ๐Ÿงฎ Online Simulator for ASK, FSK, and PSK Constellation ๐Ÿงฎ Some Questions and Answers Modulation ASK, FSK & PSK Constellation MATLAB Simulink MATLAB Code Comparisons among ASK, PSK, and FSK    Comparisons among ASK, PSK, and FSK   Simulator for Calculating Bandwidth of ASK, FSK, and PSK The baud rate represents the number of symbols transmitted per second. Both baud rate and bit rate a...

MATLAB Code for Pulse Amplitude Modulation (PAM) and Demodulation

๐Ÿ“˜ Overview & Theory of Pulse Amplitude Moduation (PAM) ๐Ÿงฎ MATLAB Code for Pulse Amplitude Modulation and Demodulation of Analog Signal and Digital Signal ๐Ÿงฎ Simulation results for comparison of PAM, PWM, PPM, DM, and PCM ๐Ÿ“š Further Reading ๐Ÿ“‚ Other Topics on Pulse Amplitude Modulation ... ๐Ÿงฎ MATLAB Code for Pulse Amplitude Modulation and Demodulation of an Analog Signal (2) ๐Ÿงฎ MATLAB Code for Pulse Amplitude Modulation and Demodulation of Digital data ๐Ÿงฎ Other Pulse Modulation Techniques (e.g., PWM, PPM, DM, and PCM)   Pulse Amplitude Modulation (PAM) & Demodulation of an Analog Message Signal MATLAB Script clc; clear all; close all; fm= 10; % frequency of the message signal fc= 100; % frequency of the carrier signal fs=1000*fm; % (=100KHz) sampling frequency (where 1000 is the upsampling factor) t=0:1/fs:1; % sampling rate of (1/fs = 100 kHz) m=1*cos(2*pi*fm*t); % Message signal with per...

Difference between AWGN and Rayleigh Fading

๐Ÿ“˜ Introduction, AWGN, and Rayleigh Fading ๐Ÿงฎ Simulator for the effect of AWGN and Rayleigh Fading on a BPSK Signal ๐Ÿงฎ MATLAB Codes ๐Ÿ“š Further Reading Wireless Signal Processing Gaussian and Rayleigh Distribution Difference between AWGN and Rayleigh Fading 1. Introduction Rayleigh fading coefficients and AWGN, or additive white gaussian noise [↗] , are two distinct factors that affect a wireless communication channel. In mathematics, we can express it in that way.  Fig: Rayleigh Fading due to multi-paths Let's explore wireless communication under two common noise scenarios: AWGN (Additive White Gaussian Noise) and Rayleigh fading. y = h*x + n ... (i) Symbol '*' represents convolution. The transmitted signal  x  is multiplied by the channel coefficient or channel impulse response (h)  in the equation above, and the symbol  "n"  stands for the white Gaussian noise that is added to the si...

MATLAB Codes for Various types of beamforming | Beam Steering, Digital...

๐Ÿ“˜ How Beamforming Improves SNR ๐Ÿงฎ MATLAB Code ๐Ÿ“š Further Reading ๐Ÿ“‚ Other Topics on Beamforming in MATLAB ... MIMO / Massive MIMO Beamforming Techniques Beamforming Techniques MATLAB Codes for Beamforming... How Beamforming Improves SNR The mathematical [↗] and theoretical aspects of beamforming [↗] have already been covered. We'll talk about coding in MATLAB in this tutorial so that you may generate results for different beamforming approaches. Let's go right to the content of the article. In analog beamforming, certain codebooks are employed on the TX and RX sides to select the best beam pairs. Because of their beamforming gains, communication created through the strongest beams from both the TX and RX side enhances spectrum efficiency. Additionally, beamforming gain directly impacts SNR improvement. Wireless communication system capacity = bandwidth*log2(1+SNR)...

Mathematical Aspects of Beamforming in MIMO

๐Ÿ“˜ Overview ๐Ÿงฎ The Mathematics Behind Analog, Digital, and Hybrid Beamforming ๐Ÿ“š Further Reading ๐Ÿ“‚ Other Topics on Beamforming in MIMO ... MIMO / Massive MIMO Beamforming Techniques Beam steering, which permits strong directed beams towards the receiver to combat excessive pathloss, especially for higher frequency bands, immediately comes to mind when discussing mathematical aspects of Beamforming in MIMO antennas. On the other side, it also lessens signal interference and improves the effectiveness of spatial multiplexing in Massive MIMO communication. Let's go right to the mathematical parts of Beamforming, which will make it easier for you to code in Python and MATLAB. 1. Beam Steering (Analog Beamforming) In the first stage, the BS applies beam steering at the side of the mobile station (MS) while the MS enables omnidirectional transmission. In the following step, the MS uses beam steering while the BS is an...