为什么仿白噪声的功率谱密度得不到一条直线?

lihaoweiv

为什么仿白噪声的功率谱密度得不到一条平行于X轴的直线?
我用了两种方法算, 一种是先求白噪声自相关, 再将它傅立叶变换;  另一种是先傅立叶变换, 将它的模, 平方,  得到总能量, 再除以长度N.
但都得不到一条直线.
下面是我用第一种方法算,写的代码,和仿出的图.
clear; clc;
fs = 1000; t = 0: 1/fs: (1-1/fs); noise = randn(1,1000);
maxlag = 1000;                                               %最大延迟
[c, maxlags] = xcorr(noise,maxlag);                 %C为返回的自相关函数, maxlags为C的长度, 等于2*maxlag+1
n = length(c); py = abs(fft(c)); k = 0:floor(n/2-1);

figure(1)
subplot(2,2,1); plot(t,noise);  xlabel('t'); ylabel('x(t)'); title('white noise');
subplot(2,2,2); plot(maxlags/fs,c); xlabel('t'); ylabel('Rx(t)'); title('Autocorrelation of white noise');
subplot(2,2,3); plot(k, 10*log10(py(k+1))); xlabel('frequency');ylabel('power(dB)'); title('PSD of white noise');

[ 本帖最后由 ChaChing 于 2009-4-5 20:43 编辑 ]

mao

不用试了，结果就是这样，不会是理论的直线

ChaChing

个人认知是, 有限点randn是不可能如理论上的白噪声
具体处理待讯号处理高手确定

[ 本帖最后由 ChaChing 于 2009-4-5 20:39 编辑 ]

lihaoweiv

那怎么样才能得到直线呢

ChaChing

mission impossible!?

[ 本帖最后由 ChaChing 于 2009-4-5 20:53 编辑 ]

panxiefan

所谓白噪声是统计意义上的，所以不会出现理论上的直线
以前我也和你一样，有同样的问题，问过老师，都是这样解释的，样本取多点

[ 本帖最后由 ChaChing 于 2009-4-11 21:57 编辑 ]

苏小贝

将采样频率fs取为1，就可以得到描述的白噪声 的功率谱图像

苏小贝

苏小贝 发表于 2015-3-24 19:50
将采样频率fs取为1，就可以得到描述的白噪声 的功率谱图像

我试了下，只取一个采样点，，采样频率取1，得到的功率谱就是平行于x轴的直线，不过这应该不是白噪声了吧，还是属于一个特例呢？