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.




Contact Us

Name

Email *

Message *

Popular Posts

MATLAB Code for MUSIC

  MATLAB Code clc; clear; close all ; %% Step 1: Define Parameters M = 8; % Number of array sensors d = 0.5; % Sensor spacing (lambda/2) K = 2; % Number of signals N = 200; % Number of snapshots theta = [-20 30]; % True signal angles (degrees) SNR = 10; % Signal-to-noise ratio (dB) fprintf( 'Step 1: Parameters Initialized\n' ); %% Step 2: Generate Signal Sources t = 1:N; s1 = exp(1j*2*pi*0.05*t); s2 = exp(1j*2*pi*0.1*t); S = [s1; s2]; figure; plot(real(S(1,:))) title( 'Signal 1 (Real Part)' ) xlabel( 'Samples' ) ylabel( 'Amplitude' ) figure; plot(real(S(2,:))) title( 'Signal 2 (Real Part)' ) xlabel( 'Samples' ) ylabel( 'Amplitude' ) fprintf( 'Step 2: Source Signals Generated\n' ); %% Step 3: Construct Steering Matrix A = zeros(M,K); for k = 1:K A(:,k) = exp(-1j*2*pi*d*(0:M-1)'*sin(theta(k)*pi/180)); end fprintf( 'Step 3: Steering Matr...

Direction of Arrival (DoA) Online Simulator (using MUSIC)

Interactive DOA Simulator X-axis XY angle (deg): 45 XZ angle (deg): 30 Noise: 0.05 Y-axis XY angle (deg): 60 YZ angle (deg): 45 Noise: 0.05 Z-axis XZ angle (deg): 60 YZ angle (deg): 30 Noise: 0.05 Estimated DOA (deg): 0 Simulation Workflow and Mathematical Background This simulator demonstrates Direction of Arrival (DOA) estimation using three-axis sensor signals (X, Y, Z), Maximal Ratio Combining (MRC) , and the MUSIC algorithm . It allows interactive control of signal angles and noise for teaching purposes. 1. Signal Generation A pure sinewave signal of frequency f is projected onto three axes using user-defined angles in different planes: X-axis: θ XY , θ XZ Y-axis: θ XY , θ YZ Z-axis: θ XZ , θ YZ Mathematically, for each time sample t : x(t) = s(t) * cos(θ_xy_x) * cos(θ_xz_x) + n_x(t) y(t) = s(t) * sin(θ_xy_y) * cos(θ_yz_y) + n_y(t) z(t) = s(t) * sin(θ_xz_z) * sin(θ_yz_z) + n_z(t) wh...

Theoretical BER vs SNR for BPSK

Theoretical Bit Error Rate (BER) vs Signal-to-Noise Ratio (SNR) for BPSK in AWGN Channel Let’s simplify the explanation for the theoretical Bit Error Rate (BER) versus Signal-to-Noise Ratio (SNR) for Binary Phase Shift Keying (BPSK) in an Additive White Gaussian Noise (AWGN) channel. Key Points Fig. 1: Constellation Diagrams of BASK, BFSK, and BPSK [↗] BPSK Modulation Transmits one of two signals: +√Eb or −√Eb , where Eb is the energy per bit. These signals represent binary 0 and 1 . AWGN Channel The channel adds Gaussian noise with zero mean and variance N₀/2 (where N₀ is the noise power spectral density). Receiver Decision The receiver decides if the received signal is closer to +√Eb (for bit 0) or −√Eb (for bit 1) . Bit Error Rat...

MATLAB code for BER vs SNR for M-QAM, M-PSK, QPSK, BPSK (with Simulation)

🧮 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; snr_db = -5:2:25; 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) ber_psk_results(i, :) = berawgn(snr_db, 'psk', psk_orders(i), 'nondiff'); end for i = 1:length(qam_orders) ber_qam_results(i, :) = berawgn(snr_db, 'qam', qam_orders(i)); end figure; semilogy(snr_db, ber_psk_results(1, :), 'o-', 'LineWidth', 1.5, 'DisplayName', 'BPSK'); hold on; for i = 2:length(psk_orders) semilogy(snr_db, ber_psk_results(i, :), 'o-', 'DisplayName', sprintf('%d-PSK', psk_or...

PSD Calculation with FFT: MATLAB Tutorial for Signal Analysis

  Implementation Steps 1. FFT Computes the Frequency Content of a Signal FFT converts a time-domain signal to the frequency domain. If: The signal is sampled at rate $f_s$ You compute an $N_{\text{FFT}}$-point FFT Then each FFT bin corresponds to a frequency resolution of: $$\Delta f = \frac{f_s}{N_{\text{FFT}}}$$ So the FFT gives you accurate frequency content, assuming the signal is stationary and adequately sampled (Nyquist criterion met).  2. Magnitude Squared Gives Power (Not Amplitude) $$P[k] = |X[k]|^2$$ This gives power at each frequency bin, not just amplitude. It represents how much energy is present at each frequency. It's a key step for PSD.  3. Normalization Makes the PSD Physically Meaningful The equation: $$\text{PSD}[k] = \frac{|X[k]|^2}{N_{\text{FFT}} \cdot f_s \cdot U}$$ is derived from first principles and ensures that the u...

Power Spectral Density Calculation Using FFT in MATLAB

📘 📘 Overview 🧮 🧮 Steps to calculate 💻 🧮 MATLAB Codes 📚 📚 Further Reading Power spectral density (PSD) tells us how the power of a signal is distributed across different frequency components, whereas Fourier Magnitude gives you the amplitude (or strength) of each frequency component in the signal. Steps to calculate the PSD of a signal Firstly, calculate the fast Fourier transform (FFT) of a signal. Then, calculate the Fourier magnitude (absolute value) of the signal. Square the Fourier magnitude to get the power spectrum. To calculate the Power Spectral Density (PSD), divide the squared magnitude by the product of the sampling frequency (fs) and the total number of samples (N). Formula: PSD = |FFT|^2 / (fs * N) Sampling frequency (fs): The rate at which the continuous-time signal is sampled (in Hz). ...

OFDM Symbols and Subcarriers Explained

This article explains how OFDM (Orthogonal Frequency Division Multiplexing) symbols and subcarriers work. It covers modulation, mapping symbols to subcarriers, subcarrier frequency spacing, IFFT synthesis, cyclic prefix, and transmission. Step 1: Modulation First, modulate the input bitstream. For example, with 16-QAM , each group of 4 bits maps to one QAM symbol. Suppose we generate a sequence of QAM symbols: s0, s1, s2, s3, s4, s5, …, s63 Step 2: Mapping Symbols to Subcarriers Assume N sub = 8 subcarriers. Each OFDM symbol in the frequency domain contains 8 QAM symbols (one per subcarrier): Mapping (example) OFDM symbol 1 → s0, s1, s2, s3, s4, s5, s6, s7 OFDM symbol 2 → s8, s9, s10, s11, s12, s13, s14, s15 … OFDM sym...

UGC NET Electronic Science Previous Year Question Papers with Solutions

Home / Engineering & Other Exams / UGC NET 2026 PYQ ⬇️ Download Papers and Solutions 📋 Exam Pattern 💡 Preparation Tips ❓ FAQs 📊 Exam Highlights: Electronic Science (88) Feature Details Junior Research Fellowship (JRF) ₹37,000 + HRA per month Eligibility M.Sc/M.Tech in Electronics (55%) Validity of Certificate JRF (3 Years) | Lectureship (Lifetime) 📥 Download UGC NET Electronics PDFs Complete collection of previous year question papers, answer keys and explanations for Subject Code 88. Start Downloading 📂 View All Question Papers June 2025 - Question Paper Download PDF June 2025 - Solved Paper + Explanation ...