I am not certain what the probme is with pwelch. It appears to give the correct result when I plot it.
Alternatively, try this:
s = readmatrix('Sumof3MIMOsDATASET1afterDACFs10G200K512eros.csv');
Fs=10e9;
L=numel(s);
t = linspace(0, L, L)/Fs;
Ts = 1/Fs;
Fn = Fs/2;
t = linspace(0, L, L)*Ts;
ms = mean(s);
FTs = fft(s)/L;
Fv = linspace(0, 1, fix(L/2)+1)*Fn;
Iv = 1:numel(Fv);
figure
plot(t, s)
grid
xlabel('Time (s)')
ylabel('Amplitude (Units)')
figure
semilogy(Fv, abs(FTs(Iv))*2)
grid
xlabel('Frequency (Hz)')
ylabel('Amplitude (Units)')
figure
plot(Fv, mag2db(abs(FTs(Iv))*2))
grid
xlabel('Frequency (Hz)')
ylabel('Power (Units^2)')
.
Best Answer