通過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。