% plot_DSF - plot Doppler Spread function for isotropic scatterer

%% parameters
v  = 10;              % velocity m/s
fc = 1e9;             % carrier frequency
c  = 299792458;       % speed of light in m/s

fd = v/c*fc;          % doppler

%% Bessel function
ff = (-0.99:0.01:0.99)*fd;
Psi = 1./sqrt(1-(ff/fd).^2);

%% plot
plot(ff, Psi)
xlabel( 'Doppler Frequency (Hz)')
ylabel( '\Psi_d(f_d)')
grid
axis([-40 40 0 7])

%% Store name of script for automated figure generation
ScriptName = mfilename;

