fadingenv.m

OFDM multipath fading channel model jakes

% FadingEnv.m
% Draw fading envelope and pdf
function []=FadingEnv(r, name)
%******************************  variables *******************************
% r             Envelope of fading channel complex impulse response
% name          Title of the figure
%************************************************************************
% Fading Envelope: Time vs. Magnitude(dB)
fe = 20.*log10(r);
figure;clf;plot(fe);grid on
plot(fe);grid on
title(['Fading Envelope for ', name]);
xlabel('Time');
ylabel('Magnitude(dB)');

% Probability density function of fading envelope
step = 0:0.01:max(r);
r_pdf = histc(r, step);
[c,i] = max(r_pdf);
delta = step(i);
max_pdf = exp(-1/2)/delta;           % maximum pdf when x=delta
r_pdf = r_pdf./max(r_pdf).*max_pdf;  % normalization
figure;clf;plot(step,r_pdf);grid on
title (['pdf of envelope for ',name]);
xlabel('pdf');
ylabel('Magnitude');
%***************************** end of file *******************************