通过FFT进行PSD估计

Question

我估计了UWB频带中高斯单脉冲的功率谱密度（PSD）。我有两个代码，使用fft。但y轴有问题，因为我不知道它的尺寸（我需要PSD，单位为dBm/MHz）。而且它在第一个代码中应该是错误的，因为它在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。

dBm单位是相对于1 mW的功率测量值，因此您必须知道测量的阻抗。

所以要达到mW/Hz，你需要乘以1e6/R，其中R是你的阻抗。然后取10*log10的结果，你有dBm/Hz。