clear all
%*************Read data files***********
data =load('Z:\Myplaygroundd');
%*************Make it complex if necessary***********
data_complex = data(:,1) + i* data(:,2);
%*************PSD PWELCH FUNCTION***********
FFT_SIZE = 12000;
Fs_MHz = 12;
RBW = Fs_MHz * 1e6 ./ FFT_SIZE
[psd,freq] = pwelch(data_complex,[],[],FFT_SIZE,Fs_MHz,'centered','psd');
figure;
plot(freq,10*log10(psd))
xlabel('Frequency (MHz)')
ylabel('Magnitude (dB)')
grid on;
title('Double sided PSD');