%Samp9_10
Fs=10000;   %信號(hào)采樣頻率
randn('state',0);   %設(shè)置隨機(jī)狀態(tài)初始值
x = randn(16384,1);    %產(chǎn)生第一個(gè)白噪聲信號(hào)
y= randn(16384,1);   %產(chǎn)生第二個(gè)白噪聲信號(hào)
z=2*sin(2*pi*1000*[1:16384]'/Fs)+randn(16384,1);%產(chǎn)生含周期成分的噪聲信號(hào)
[Pxy,f]=csd(x,y,1024,Fs,triang(500),0);%白噪聲信號(hào)互相關(guān)譜,無重疊三角窗
[Pxz,f]=csd(x,z,1024,Fs,triang(500),0);%噪聲與周期信號(hào)互相關(guān)譜無重疊三角窗
subplot(2,1,1),plot(f,10*log10(Pxy)),grid on; %繪制x,y信號(hào)互相關(guān)譜
ylabel('互相關(guān)譜/dB'),title('噪聲信號(hào)的互相關(guān)譜')
subplot(2,1,2),plot(f,10*log10(Pxz)) ,grid on; %繪制x,z信號(hào)的互相關(guān)譜
ylabel('互相關(guān)譜/dB'),xlabel('頻率/Hz')
title('帶噪聲周期信號(hào)的互相關(guān)譜')