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');