samp9_15.m

数字信号处理的Matlab实现的随书光盘

%Samp9_15
clf;
%第一种情况：y与x相干
Fs=1000;    %采样频率
N=1024;Nfft=256;n=0:N-1;t=n/Fs; %数据长度和分段数据长度和时间序列
window=hanning(256);   %采用的窗口
noverlap=128;   %重叠数据个数
dflag='none';
randn('state',0);   %设置产生随机数的状态
xn=sin(2*pi*50*t)+randn(1,N);   %原始信号
h=fir1(30,0.2,boxcar(31));  %设计30阶截止频率为0.2的FIR滤波器
yn=filter(h,1,xn);   %运用上述的FIR滤波器进行滤波
subplot(2,1,1)
cohere(xn,yn,Nfft,Fs,window,noverlap,dflag);  %计算其相干函数
xlabel('频率/Hz'); ylabel('相干函数估计')
title('Yn与xn相干')
%第二种情况：y与x不相干
Fs=1000;   %采样点数
N=1024;Nfft=256;n=0:N-1;t=n/Fs; %数据长度和分段数据长度和时间序列
window=hanning(256);   %采用的窗口
noverlap=128;   %重叠数据长度
dflag='none';
randn('state',0);   %设置产生随机数的状态
xn=ones(1,N);         %第一个信号
yn=sin(2*pi*50*t);    %第二个信号
subplot(2,1,2)
cohere(xn,yn,Nfft,Fs,window,noverlap,dflag);  %两个信号的相干函数
title('Yn 与xn不相干');
xlabel('频率/Hz'); ylabel('相干函数估计')
legend('相干函数',4);