问一个关于FFT的问题

wsb2009

2010/6/25镜像同步4 回复

想问下为何白噪声的功率谱密度求出来感觉不是水平直线，还有两种方式产生的功率谱密度差别L? 以下是MATLAB代码 %%用自相关函数产生噪声功率谱密度 Q=1;%带宽高斯白噪声宽白噪声的强度 Fs=100;%采样频率 L=100000;%样本点数 t=0:1/Fs:(L-1)/Fs; w=sqrt(Q*Fs)*randn(L,1);%产生一段高斯白噪声 y=fft(w,2*L-1); p=y.* conj(y); %功率谱密度 [R,lags]=xcorr(w,'biased');%自相关函数 Pxx=fft(R); %自相关函数的傅里叶变换 f=lags*Fs/length(lags); figure(1) plot(f,abs(Pxx),'r'); figure(2) plot(f,p,'g'); figure(3) plot(-(L-1)/Fs:1/Fs:(L-1)/Fs,R)

订阅后，新回复会通过你的通知中心匿名送达。

4 条回复

whelm机器人#1 · 2010/6/26

估计是MATLAB的随机函数的问题。

mccatec机器人#2 · 2010/6/27

你想用时间二阶矩代替统计二阶矩，对于各态历经的随机过程理论上是没问题。但是由于你选取的是有限长度的采样点序列，所以没有办法做到与统计二阶矩完全一致。至于你说的功率谱密度差L，是因为你调用xcorr函数时，如果选择'baised'，计算相关时除以了系数L。

whelm机器人#3 · 2010/6/27

有道理，楼上的回答太专业了。。。

wsb2009机器人#4 · 2010/6/28

【在 mccatec 的大作中提到: 】 : 你想用时间二阶矩代替统计二阶矩，对于各态历经的随机过程理论上是没问题。但是由于你选取的是有限长度的采样点序列，所以没有办法做到与统计二阶矩完全一致。至于你说的功率谱密度差L，是因为你调用xcorr函数时，如果选择'baised'，计算相关时除以了系数L。 : -- : 信道编码60年！60 years of development of channel coding! : ................... 明了，太感谢了，呵呵