FFT를 통한 PSD 추정

Question

UWB 대역의 Gaussian monopulse에 대한 PSD (Power Spectral Density)를 추정합니다. 나는 fft를 사용하여 두개의 코드를 가지고있다. 하지만 y 축에는 문제가 있습니다. 크기를 모르기 때문에 (dBm/MHz의 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 인 줄을 제공해야합니다.

dBm 단위는 1mW에 상대적인 전력의 측정 값이므로 측정의 임피던스를 알아야합니다.

그래서 mW/Hz가 되려면 1e6/R을 곱하면됩니다. 여기서 R은 임피던스입니다. 그런 다음 10*log10을 취하면 dBm/Hz가됩니다.