samp9_10.m

MATLAB7.x数字信号处理 光盘内容

%Samp9_10
Fs=10000;   %信号采样频率
randn('state',0);   %设置随机状态初始值
x = randn(16384,1);    %产生第一个白噪声信号
y= randn(16384,1);   %产生第二个白噪声信号
z=2*sin(2*pi*1000*[1:16384]'/Fs)+randn(16384,1);%产生含周期成分的噪声信号
[Pxy,f]=csd(x,y,1024,Fs,triang(500),0);%白噪声信号互相关谱,无重叠三角窗
[Pxz,f]=csd(x,z,1024,Fs,triang(500),0);%噪声与周期信号互相关谱无重叠三角窗
subplot(2,1,1),plot(f,10*log10(Pxy)),grid on; %绘制x,y信号互相关谱
ylabel('互相关谱/dB'),title('噪声信号的互相关谱')
subplot(2,1,2),plot(f,10*log10(Pxz)) ,grid on; %绘制x,z信号的互相关谱
ylabel('互相关谱/dB'),xlabel('频率/Hz')
title('带噪声周期信号的互相关谱')