Neverjetno Precision Problem je FFT Matlab

[qslazio]

preprost segment skript storiti FFT za ton plus dodatek belega šuma.

nbpts = 32768;
fs = 1e3 * nbpts;
nbprds = 7;
fi = fs / nbpts * nbprds;
amplituda = 1;
pnoise = 1e-6;

hrup = naključno ("norma", 0, sqrt (pnoise), 1, nbpts);
signal = amplitude * sin (2 * pi * fi / fs * [0: nbpts-1]) hrupa;

Sxx = abs (FFT (signal). / Nbpts). ^ 2;
SNN = abs (FFT (hrup). / Nbpts). ^ 2;
= sum(Sxx) - 0.5*amplitude^2

Pn_wi_signal
= vsota (Sxx) - 0,5 * amplitude ^ 2= sum(Snn)

Pn_wo_signal
= vsota (SNN)

Po vožnjo zgoraj scenarij, lahko vidimo, da Pn_wi_signal ni enako Pn_wo_signal.Pn_wo_signal je neposredna FFT, da hrup sama.
Pn_wi_signal (se spreminja veliko) je veliko večja kot Pn_wo_signal, ki ima prav.

Vsakdo ve, kaj je narobe z Pn_wi_signal?
Hvala veliko!