[那年下雪]

为什么我通过相关和fft变换来仿真高斯白噪声功率谱不是一条近似直线，相关函数仿真是冲击响应，且冲击值为1，但功率谱不是波动在1附近，
请指教 万分感激
clear all;
fs=1000;
nfft=1024;
idx=0:round(nfft/2-1);
k=idx*fs/nfft;
t=0:1/fs:1;
x1=randn(1,length(t));
[rx1 lag1]=xcorr(x1,'unbiased');%求自相关函数
fx1=fft(rx1,nfft);%求功率谱密度;
我还利用多次循环也不可以？谢谢了