Спектральная плотность мощности гусинового белого шума

Question

Я пытаюсь генерировать белый шум Гуссина со стандартным отклонением 7% и частоту выборки 2000 Hz. Затем вычислите спектральную плотность мощности (PSD). Проблема в том, что когда я писал PSD, я знаю, что мне нужно получить плоский PSD. Однако, я не

Я использовал ниже код

clear; 
clc; 
x =7*randn(1,10000); 
N=10000; 
sample_frequency=2000; 
fax_bin=[0 :N/2-1]; 
fax_Hz=fax_bin*sample_frequency/N; 
FFT_A=fft(x); 
spectra=FFT_A.*conj(FFT_A); 
figure 
plot(fax_Hz,spectra(1,1:5000)); 

Answer 1

Это случайный процесс, так что вы никогда не получите это действительно плоский, он становится лучше с большим количеством образцов. Кроме того, построение графика на уровне дБ имеет больше смысла.

plot(fax_Hz,10*log10(spectra(1,1:N/2)));