How to simulate two channel estimates within coherence time

Hello all,

So I am having trouble in simulating my research.

I need to estimate a wireless channel (rayleigh fading) at both the Tx and Rx within coherence time to obtain similar channnel estimates but not completely the same, rather, somewhat similar.

I tried code out there already, however, when I simulate it twice, it ends up producing uncorrelated channels.

I need a lot of help on this. Thank you!

Answer 1

[Muskan]

Hi Joshua,

I understand that you want to estimate a Rayleigh Fading at both the “Tx” and “Rx” within coherence time to obtain similar channel estimates.

As per my understanding of the question, you can try the following:

coherence_time = 1; % Coherence time in seconds 

sampling_rate = 1000; % Sampling rate in Hz 

doppler_frequency = 10; % Doppler frequency in Hz 

num_samples = coherence_time * sampling_rate; 

% Generate channel estimate at the Tx 

channel_tx = (randn(num_samples, 1) + 1i * randn(num_samples, 1)) / sqrt(2); 

% Apply Doppler shift to the channel at the Tx 

t = (0:num_samples-1) / sampling_rate; 

doppler_shift_tx = exp(1i * 2 * pi * doppler_frequency * t.'); 

channel_tx = channel_tx .* doppler_shift_tx; 

% Generate channel estimate at the Rx 

channel_rx = (randn(num_samples, 1) + 1i * randn(num_samples, 1)) / sqrt(2); 

% Apply Doppler shift to the channel at the Rx 

doppler_shift_rx = exp(1i * 2 * pi * doppler_frequency * t.'); 

channel_rx = channel_rx .* doppler_shift_rx; 

% Add small variations to the channel estimates 

variation = 0.1; % Adjust the variation level as desired 

channel_tx = channel_tx + variation * (randn(num_samples, 1) + 1i * randn(num_samples, 1)); 

channel_rx = channel_rx + variation * (randn(num_samples, 1) + 1i * randn(num_samples, 1)); 

% Plot the channel estimates 

time = (0:num_samples-1) / sampling_rate; 

figure; 

subplot(2, 1, 1); 

plot(time, abs(channel_tx)); 

xlabel('Time (s)'); 

ylabel('Magnitude'); 

title('Channel Estimate at Tx'); 

subplot(2, 1, 2); 

plot(time, abs(channel_rx)); 

xlabel('Time (s)'); 

ylabel('Magnitude'); 

title('Channel Estimate at Rx'); 

In this example, we generate channel estimates at both the “Tx” and “Rx” using complex Gaussian random processes. We then apply a Doppler shift to the channels to introduce a varying fading effect. Small variations are added to the channel estimates to make them similar but not identical. Finally, the channel estimates are plotted to visualize the fading characteristics.

You can adjust the coherence time, sampling rate, Doppler frequency, and variation level according to your requirements.

I hope the above information helps resolve your query.

Thanks