白噪声的功率谱密度估计

《白噪声的功率谱密度估计》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、randn('state',0)NFFT=1024;%NFFT为取样点数Fs=10000;%Fs为取样频率t=(0:NFFT-1)/Fs;%时间y=randn(NFFT,1);%产生高斯白噪声，2*pi为其功率。figure(1);subplot(3,1,1);plot(t,y);gridon;title('白噪声波形');%计算白噪声的自相关系数[cory,lags]=xcorr(y,200,'unbiased');subplot(3,1,2);plot(lags,cory);%自相关函数(无偏差的)，cory为要求的自相关函数，lag为自相关函

2、数的长度。title('白噪声相关数');gridon;f=fft(cory);%对自相关系数进行傅里叶变换k=abs(f);fl=(0:length(k)-1)*Fs/length(k);%fl为fc的长度。subplot(3,1,3)plot(fl,k);gridon;title('白噪声功率谱谱');%自相关函数的傅里叶变换：即功率谱密度。其中。k是cory付里叶变换后的幅值。