% plot_PSDRect - plot PSD with rectangular pulse

%% frequency axis: (normalized to 1/T)
ffT      = -10:0.01:10;

%% PSD
Ps = sinc( ffT ).^2;

%% plot
plot( ffT, 10*log10(Ps) )
grid
xlabel('Normalized Frequency (fT)')
ylabel('PSD (dB)')
axis( [-10 10 -70 0])

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