% plot_TCF - plot time cohrence 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
tt = 0:1/(20*fd):10/fd;
Psi = besselj(0, 2*pi*fd*tt);

%% plot
plot(tt/1e-3, Psi)
xlabel( 'Time \Deltat (ms)')
ylabel( '\Psi_H(\Deltat)')
grid

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

