FFT- Algorithm
X=fft(eirp,N); % conversion from MATLAB spectrum to Fourier spectrum X=X/N
% DOUBLE-SIDED ESD of the waveform E = fftshift(abs(X).^2/(df^2));
% SINGLE-SIDED ESD of the waveform Ess = 2.*E((N/2+1):N);
% PSD of the combination in dBm/MHz PSD = 10 * log10 ((1/Ts) * Ess/377)+ 90; %30 for dBm, 60 for /MHz
frequency=linspace(0,df*size(PSD,2),size(PSD,2)); figure(); plot(frequency,PSD); axis([0 10e9 -120 0]); grid on; title('PSD of the Time Domain Pulse UWB Pulse'); xlabel('Frequency [Hz]'); ylabel('PSD [dBm/MHz]');