Estimation PSD via FFT

Question

J'estime la densité spectrale de puissance (PSD) pour la monopulse gaussienne dans la bande UWB. J'ai deux codes, en utilisant fft. Mais il y a un problème avec les axes y, puisque je ne connais pas sa dimension (j'ai besoin de PSD en dBm/MHz). Et il devrait également être une erreur dans le premier code, car il montre une seule valeur sur les axes 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)); 

Merci beaucoup pour toute aide!

Answer 1

Votre deuxième parcelle devrait vous donner un terrain avec des unités V^2/Hz. Les unités dBm mesurent la puissance par rapport à 1 mW. Vous devez donc connaître l'impédance de votre mesure. Donc pour arriver à mW/Hz, vous voudriez multiplier par 1e6/R, où R est votre impédance. Ensuite, prenez 10*log10 du résultat et vous avez dBm/Hz.