Replace pwrsd.m

pwrsd.m

@@ -17,7 +17,7 @@ dt = 1/freq;
 dF = freq/nt;
 f = -freq/2:dF:freq/2-dF;
 
-p = abs(fftshift(fft(s)))./nt;
+p = abs(fftshift(fft(s))).*(dt/nt);
 
 if ifplot
     figure; set(gcf,'Color','w');