%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);