为什么仿白噪声的功率谱密度得不到一条直线?
为什么仿白噪声的功率谱密度得不到一条平行于X轴的直线?我用了两种方法算, 一种是先求白噪声自相关, 再将它傅立叶变换;另一种是先傅立叶变换, 将它的模, 平方,得到总能量, 再除以长度N.
但都得不到一条直线.
下面是我用第一种方法算,写的代码,和仿出的图.
clear; clc;
fs = 1000; t = 0: 1/fs: (1-1/fs); noise = randn(1,1000);
maxlag = 1000; %最大延迟
= 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 编辑 ] 不用试了,结果就是这样,不会是理论的直线
回复 楼主 lihaoweiv 的帖子
个人认知是, 有限点randn是不可能如理论上的白噪声具体处理待讯号处理高手确定
[ 本帖最后由 ChaChing 于 2009-4-5 20:39 编辑 ] 那怎么样才能得到直线呢 mission impossible!?
[ 本帖最后由 ChaChing 于 2009-4-5 20:53 编辑 ] 所谓白噪声是统计意义上的,所以不会出现理论上的直线
以前我也和你一样,有同样的问题,问过老师,都是这样解释的,样本取多点
[ 本帖最后由 ChaChing 于 2009-4-11 21:57 编辑 ] 将采样频率fs取为1,就可以得到描述的白噪声 的功率谱图像 苏小贝 发表于 2015-3-24 19:50
将采样频率fs取为1,就可以得到描述的白噪声 的功率谱图像
我试了下,只取一个采样点,,采样频率取1,得到的功率谱就是平行于x轴的直线,不过这应该不是白噪声了吧,还是属于一个特例呢?
页:
[1]