Skip to main content

Adaptive Channel Estimation for Millimeter wave 5G


Adaptive Channel Estimation

7.1. Introduction

Channel estimation is to acquire site specific channel model of a particular environment. Millimetre wave cannot travel much distance in atmosphere due to its severe pathloss and high refraction properties. It also experiences high penetration loss. So, we need to equip MIMO system to avail beamforming. And it is also needed to deploy large antenna arrays in MIMO system to form narrower and stronger beam to fulfil the link margin gain at receiver side.

So, we need to train the beam to find the best beamforming vector at transmitter side and best combing vector pair at receiver side to maximize the signal to noise ratio. The beam training can be realized by analog beamformer but the limitation of analog beamforming is that only one data stream with same amplitude and different phases is possible. That cannot provide high data rates as required. That’s why we need to focus on MIMO hybrid architecture. In that architecture analog beam steering is realized by large antenna arrays and baseband precoding is realized by low dimensional digital precoder. It also helps us to maintain low overhead in network.

Especially, for vehicular communication channel varies fast, so, channel estimation plays a crucial role to find the best communicating beam from BS to find the best signal to noise ratio from beam steering.

7.2. System Model:

In MIMO channel estimation process, firstly we train the beam using beam steering method to pick the best beam from Tx and Rx side both. There is only one RF chain available in analog precoder with can enable single data stream between transmitter and receiver. So, we further go for hybrid architecture with equips lower dimensional baseband precoder which enables multiple simultaneous data stream and cancel interference between them as shown in figure 7.1.

Fig 7.1: Architecture of transceiver at Tx & Rx both which contains RF chains and baseband beamformers

On the other hand, hybrid precoder is used for low overhead in network because it uses relatively lesser RF chain as compare to digital optimal precoder without compromising much in performance.

Here I will discuss communication between a BS (base station) and vehicles for mm Wave communication channel. Both NBS and Nv have linear arrays (ULAs), with half-wavelength spaced antenna elements in a compact area uniformly. Let assume, fn is the beamforming vector at BS side. BS need to transmit a symbol, s, such that |s|2 = 1. On the other hand, vehicle assigns a unit-norm measurement vector wm. Now, the signal received at MS (user) side is written by

…………………………………………..(Appendix B)

Here in above equation BS and MS both enable beam steering. At a moment they find the best beam pairs to enable communication. Here ‘m’ indicates best beamforming vector from BS and ‘n’ from receiver side.

In the above figure 7.1, BS employs baseband precoder, FBB of size NRF X NS which is followed by RF precoder, FRF of size NBS X NRF. BS can communicate to MS thru NS data streams simultaneously and similarly it is applicable to MS side, where NS ≤ NRF ≤NBS and, NS ≤NRF≤NMS. Number of maximum possible simultaneous data streams between Tx and Rx is always less than or equal to number of RF chain.

If FT= FRF FBB is the combined then BS precoding matrix is of size NBS X NS, then transmitted signal is denoted by,

……………………………………………………………………….(7.1)

We know in case of mmWave channel estimation most of MPCs are so weak that makes mmWave channel matrix sparse. Only few MPCs are strong to avail communication between Tx and Rx. Let assume, the numbers of energy rich limited MPCs are L. Now, L MPCs contain AOA & AOD, and corresponding gain of each path. It is now essential to build the BS and MS training precoders and combiners in order to reduce the training overhead.

We’ll take advantage of the mmWave channel’s low scattering properties. Which indicates low rank property of the millimetre wave channel matrix or estimation problem as a sparse problem. It is shown here how adaptive compressed sensing (CS) works on certain objectives to find the training precoders and combiners for low overhead mmWave band communication. A novel hierarchical multi resolution codebook based on the hybrid analog/digital interface architecture is used to train beamforming vectors.

Fig 7.2: With beamforming vectors in each subset, architecture of a multi-resolution codebook with total beamforming vector N=8 and, total stage, K=2.

Fig 7.3: The beampatterns for the resulting beamforming vectors are shown in the first three codebook layers of an example hierarchical codebook.

In the above figure 7.2, let assume we want to cover the transmitted signal from BS in all directions in 0 to 360 degree azimuth angle range around the nearby area of BS. Then we divide the azimuth angle range of transmission in subparts, for example we divide it into two parts, and assign beamforming vector to each part, where one beamforming vector is supposed to cover 0 to 180 degree range and another for 180 to 360 degree range and it is only for the 1st stage. Then in second stage each azimuth angle range is further subdivided into different parts (shown in figure 7.4) and so on, only to find best beamformer and combiner pair to maximize the SNR.

Algorithm 2: Adaptive Estimation Algorithm for Single-Path mmWave Channels

Let assume for ‘P’ number of antenna elements in Tx and ‘Q’ number of antenna element in Rx, so there will be total P X Q possible paths between BS and MS. But due to high pathloss and poor scattering nature, only few stronger MPCs reach up to receiver. Here number of detectable strong paths or MPCs is ‘N’. The BS in the above algorithm uses the K training precoding vectors at the 1st level of the multi-resolution codebook, F, at the initial step. Similarly, the MS uses the measurement vectors of the first stage of W to combine the beamforming vectors from transmitter. Now at second stage each slot of previous stage is divided into K slots and we assign beamforming vector accordingly. At the final stage the total number of beamforming vector will be KS which is equal to total number of stronger MPCs or path, N. In adaptive estimation process we apply beam steering at different stages at find the best path with higher SNR (Alkhateeb, 2014).

7.3. Channel Estimation for Vehicular communication

The most promising technology for enabling ultra-fast and high-data-rate exchanges is mm wave communication. In case of vehicular communication BS can be placed in a particular place along road-side, or in simple word previously known location of road or infrastructure can improve initial access of communication. Millimeter wave connectivity has the advantage of supporting multi-Gbps throughput, but it suffers from substantial pathloss. On the other hand, as mm Wave wavelength lies between 1-10 millimetres so enable packing many antenna elements in a compact area.

Let assume, there is maximum L different beamforming and combining or measurement vectors is possible in different time slots at last stage (say, kth stage) of adaptive channel estimation process. So, there is total of L X L measurements among which we have to find which beamforming and combining or measurement vector pair has higher gain at receiver side. Now, arranging in matrix form, we can write

……………………………………………………………………..(7.2)

where W = [w1,…,wL], and F = [f1,…,fL]

N ~ Noise at different time slot with zero mean and variance σ2 , such that N ~ CN(0,σ2I)

Noise N are statistically independent for different time slots.

Fig 7.4: illustration of Adaptive channel estimation, where first stage contain ‘m’ number of beamforming vector at different time slot; then at stage 2, each previous time slot divided into 4 slots where each slot represents a angular range of a beamforming vector.

Here in figure, the coverage area in azimuth angular range is subdivided into four subranges and assigned a beam beamforming vector for each subrange or slot, similarly as mentioned for figure 7.3. & 7.4 above. In this adaptive channel estimation process required codebook for beamforming and measurement vector at Tx and Rx, respectively, are pre-computed if it is location aided. In this estimation process it is done in multiple stages. Our goal is to find best beamforming and measuring vector pair. At stage ‘k’ for received signal Y(k) is denoted as,

……………………………..(7.3)

At the beginning, at first stage we assume that suitable AOD and AOA for communication lies between 0 to 1800 range. And four beamforming vectors cover the whole angular range by applying corresponding beamwidth as shown in figure 11. At the next stage (stage 2) AOD and AOA range for a beam forming vector is further divided into four slots. In the above equations ^m(k-1) and ^n(k-1) subscripts denotes corresponding beamforming and combining vectors at stage k-1 and so forth. The Tx and Rx sound the channel in numerous AOD and AOA ranges corresponding to various stages in this operation, and the strongest SNR achieved is chosen. The previous ranges are divided into L sub-ranges in the next level. To achieve desired RSSI we apply the next stages in this channel estimation process. If 1st stage consists of L number of beamforming vectors then at stage k it will contain Lk independent slots. Here, it is now necessary that beamforming and combining or measurement vector will always choose from final stage k, it can be chosen from previous stages if we achieve stronger SNR at receiver side. Now, angle of arrival and departure (AOA & AOD) vector dependent gain can be denoted by,

……………………………………………………………….(7.10)

Let assume, at stage k, the range of the m-th weight vector is Rm(k) at the Tx or Rx, where the final stage k contains Lk disjoint sub-ranges which covers the AOA or AOD range from 00 to 1800 [0, π] . The mathematical expression for Rm(k) is given by,

……………………………….(7.4)

Fig 7.5: Here disc indicates priory known location information of the vehicle. Angle of departure (AOD) from BS and angle of arrival (AOA) at vehical side, which are expected to be priory known, represented as, uθ = [θdisk - arcsin( D/2ddisk), θdisk + arcsin( D/2ddisk)] and uφ = [Ï€ - θdisk - arcsin( D/2ddisk), Ï€ -θdisk +arcsin( D/2ddisk)] , respectively. Where, –Ï€/2 ≤ arcsin(x) ≤ Ï€/2 (Garcia, 2016).

By the way, channel estimation parameters totally depends on chosen number of total stages, how many beamforming vectors will be assigned at 1st stage, and codebooks at the Tx and Rx side (Garcia, 2016).

7.4. Mm wave channel estimation in the Presence of location Information:

It is assumed that the vehicle is located in the blue disk shown in figure 12 which has radius of D. The distance from BS to the centre of the disc is ddisk where angle of departure (AOD) from BS to the centre of disk is given by, θdisk. From the available location information, further we convert it to corresponding AOD/AOA by using basic trigonometric ratios. It has been said before that adaptive channel estimation method is multi-level process, so, which AOD/AOA range do not satisfy the link margin at receiver, is cancelled. Then we proceed to the next stage of channel estimation. Similarly, pathloss due to particular AOD/AOA can be calculated from the location information.

In the above figure 12, the priory known angle of departure (AOD) is indicated as, θ belongs to set u θ and angle of arrival φ belongs to set u φ ; such that sets u φ , u θ lies between [0, Ï€) . Suppose, the sequences of beamforming and measurement vectors are m(1),…,m(K) and n(1),…,n(K), respectively, where the indices indicates the stage. If P(m(K),n(K)) denotes the probability of beamforming and measuring vector pairs at stage k, then,

…………………………………………..(7.5)

All stages in adaptive channel estimation process is independent where individual stages happen at individual time intervals. We choose a beamforming and measuring vector pair which contributes a stronger RSSI at receiver at a particular stage. The indicator for mth beamforming vector nth measurement weight vectors, where each slot consist of L number of beamforming vector. The SNR value for a particular stage of adaptive channel estimation is defined as

……………………………..(7.6)

Where k= number of stage; α = channel gain; gv2 = Rx side antenna gain; gBS2 = Tx side antenna gain; P(m, n) is probability of mth beamformer vector and nth combiner vector for a particular stage.


People are good at skipping over material they already know!

View Related Topics to







Contact Us

Name

Email *

Message *

Popular Posts

Amplitude, Frequency, and Phase Modulation Techniques (AM, FM, and PM)

📘 Overview 🧮 Amplitude Modulation (AM) 🧮 Online Amplitude Modulation Simulator 🧮 MATLAB Code for AM 🧮 Q & A and Summary 📚 Further Reading Amplitude Modulation (AM): The carrier signal's amplitude varies linearly with the amplitude of the message signal. An AM wave may thus be described, in the most general form, as a function of time as follows .                       When performing amplitude modulation (AM) with a carrier frequency of 100 Hz and a message frequency of 10 Hz, the resulting peak frequencies are as follows: 90 Hz (100 - 10 Hz), 100 Hz, and 110 Hz (100 + 10 Hz). Figure: Frequency Spectrums of AM Signal (Lower Sideband, Carrier, and Upper Sideband) A low-frequency message signal is modulated with a high-frequency carrier wave using a local oscillator to make communication possible. DSB, SSB, and VSB are common amplitude modulation techniques. We find a lot of bandwi...

Analog vs Digital Modulation Techniques | Advantages of Digital ...

Modulation Techniques Analog vs Digital Modulation Techniques... In the previous article, we've talked about the need for modulation and we've also talked about analog & digital modulations briefly. In this article, we'll discuss the main difference between analog and digital modulation in the case of digital modulation it takes a digital signal for modulation whereas analog modulator takes an analog signal.  Advantages of Digital Modulation over Analog Modulation Digital Modulation Techniques are Bandwidth efficient Its have good resistance against noise It can easily multiple various types of audio, voice signal As it is good noise resistant so we can expect good signal strength So, it leads high signal-to-noise ratio (SNR) Alternatively, it provides a high data rate or throughput Digital Modulation Techniques have better swathing capability as compared to Analog Modulation Techniques  The digital system provides better security than the a...

Shannon Limit Explained: Negative SNR, Eb/No and Channel Capacity

Understanding Negative SNR and the Shannon Limit Understanding Negative SNR and the Shannon Limit An explanation of Signal-to-Noise Ratio (SNR), its behavior in decibels, and how Shannon's theorem defines the ultimate communication limit. Signal-to-Noise Ratio in Shannon’s Equation In Shannon's equation, the Signal-to-Noise Ratio (SNR) is defined as the signal power divided by the noise power: SNR = S / N Since both signal power and noise power are physical quantities, neither can be negative. Therefore, the SNR itself is always a positive number. However, engineers often express SNR in decibels: SNR(dB) When SNR = 1, the logarithmic value becomes: SNR(dB) = 0 When the noise power exceeds the signal power (SNR < 1), the decibel representation becomes negative. Behavior of Shannon's Capacity Equation Shannon’s channel capacity formula is: C = B log₂(1 + SNR) For SNR = 0: log₂(1 + SNR) = 0 When SNR becomes smaller (in...

Online Simulator for ASK, FSK, and PSK

Try our new Digital Signal Processing Simulator!   Start Simulator for binary ASK Modulation Message Bits (e.g. 1,0,1,0) Carrier Frequency (Hz) Sampling Frequency (Hz) Run Simulation Simulator for binary FSK Modulation Input Bits (e.g. 1,0,1,0) Freq for '1' (Hz) Freq for '0' (Hz) Sampling Rate (Hz) Visualize FSK Signal Simulator for BPSK Modulation ...

MATLAB code for BER vs SNR for M-QAM, M-PSK, QPSk, BPSK, ...(with Online Simulator)

🧮 MATLAB Code for BPSK, M-ary PSK, and M-ary QAM Together 🧮 MATLAB Code for M-ary QAM 🧮 MATLAB Code for M-ary PSK 📚 Further Reading MATLAB Script for BER vs. SNR for M-QAM, M-PSK, QPSK, BPSK % Written by Salim Wireless clc; clear; close all; num_symbols = 1e5; snr_db = -20:2:20; psk_orders = [2, 4, 8, 16, 32]; qam_orders = [4, 16, 64, 256]; ber_psk_results = zeros(length(psk_orders), length(snr_db)); ber_qam_results = zeros(length(qam_orders), length(snr_db)); for i = 1:length(psk_orders) psk_order = psk_orders(i); for j = 1:length(snr_db) data_symbols = randi([0, psk_order-1], 1, num_symbols); modulated_signal = pskmod(data_symbols, psk_order, pi/psk_order); received_signal = awgn(modulated_signal, snr_db(j), 'measured'); demodulated_symbols = pskdemod(received_signal, psk_order, pi/psk_order); ber_psk_results(i, j) = sum(data_symbols ~= demodulated_symbols) / num_symbols; end end for i...

BER vs SNR for M-ary QAM, M-ary PSK, QPSK, BPSK, ...(MATLAB Code + Simulator)

📘 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. BER = (number of bits received in error) / (total number of tran...

MATLAB Code for Pulse Width Modulation (PWM) and Demodulation

📘 Overview & Theory 🧮 MATLAB Code for Pulse Width Modulation and Demodulation 🧮 Generating a PWM Signal in detail 🧮 Other Pulse Modulation Techniques (e.g., PWM, PPM, DM, and PCM) 🧮 Simulation results for comparison of PAM, PWM, PPM, DM, and PCM 📚 Further Reading   MATLAB Code for Analog Pulse Width Modulation (PWM) clc; clear all; close all; fs=30; %frequency of the sawtooth signal fm=3; %frequency of the message signal sampling_frequency = 10e3; a=0.5; % amplitide t=0:(1/sampling_frequency):1; %sampling rate of 10kHz sawtooth=2*a.*sawtooth(2*pi*fs*t); %generating a sawtooth wave subplot(4,1,1); plot(t,sawtooth); % plotting the sawtooth wave title('Comparator Wave'); msg=a.*sin(2*pi*fm*t); %generating message wave subplot(4,1,2); plot(t,msg); %plotting the sine message wave title('Message Signal'); for i=1:length(sawtooth) if (msg(i)>=sawtooth(i)) pwm(i)=1; %is message signal amplitude at i th sample is greater than ...

Theoretical vs. simulated BER vs. SNR for ASK, FSK, and PSK (MATLAB Code + Simulator)

📘 Overview 🧮 Simulator for calculating BER 🧮 MATLAB Codes for calculating theoretical BER 🧮 MATLAB Codes for calculating simulated BER 📚 Further Reading BER vs. SNR denotes how many bits in error are received for a given signal-to-noise ratio, typically measured in dB. Common noise types in wireless systems: 1. Additive White Gaussian Noise (AWGN) 2. Rayleigh Fading AWGN adds random noise; Rayleigh fading attenuates the signal variably. A good SNR helps reduce these effects. Simulator for calculating BER vs SNR for binary ASK, FSK, and PSK Calculate BER for Binary ASK Modulation Enter SNR (dB): Calculate BER Calculate BER for Binary FSK Modulation Enter SNR (dB): Calculate BER Calculate BER for Binary PSK Modulation Enter SNR (dB): Calculate BER BER vs. SNR Curves MATLAB Code for Theoretical BER % The code is written by SalimWireless.Com clc; clear; close all; % SNR va...