Оценка PSD с помощью FFT

Question

Я оцениваю спектральную плотность мощности (PSD) для гауссова монопульса в диапазоне UWB. У меня есть два кода, используя fft. Но есть проблема с у-осями, так как я не знаю его размерности (мне нужен PSD в дБм/МГц). А также это должно быть ошибкой в ​​первом коде, поскольку оно показывает только одно значение в осях y.

Code1 

fs=1e11; 

g=0.1e-9; 

tmax=1e-9; 

fftl=2048; 


t=(-tmax:1./fs:tmax)'; 

s=t./(g.^3.*sqrt(2.*pi)).*exp(-t.^2./(2.*g.^2))./2.5e19; 

figure(1) 

plot(t,s); 

xlabel('Time, s'); 

ylabel('Amplitude, V'); 

ffts=abs(fft(s,fftl)); 

ffts=2.*ffts./fftl; 

fftp=abs(ffts.*conj(ffts))./2; 

fftps=(fftp-30)./1e-6; 

f=0:fs./fftl:fs./2-fs./fftl; 

figure(2) 

plot(f,fftps(1:length(f))),grid; 

xlabel('Frequency, Hz'); 

Code2 

fs=1e11; 

g=0.1e-9; 

tmax=1e-9; 

t=(-tmax:1./fs:tmax)'; 

s=t./(g.^3.*sqrt(2.*pi)).*exp(-t.^2./(2.*g.^2))./2.5e19; 

figure(1) 

plot(t,s); 

xlabel('Time, s'); 

ylabel('Amplitude, V'); 

S=fft(s,8192); 

f=fs.*(0:4095)./8192; 

Pss=S.*conj(S)./8192; 

figure(2) 

plot(f,Pss(1:4096)); 

Благодарим за помощь!

Answer 1

Ваш второй участок должен предоставить вам участок с блоками V^2/Hz.

дБм единиц измерения мощности относительно 1 мВт, поэтому вам необходимо знать сопротивление вашего измерения.

Итак, чтобы перейти на mW/Hz, вы хотели бы умножить на 1e6/R, где R - ваш импеданс. Затем возьмите 10*log10 результата, и у вас есть дБм/Гц.