%Samp9_2
clf; N=1000; Fs=500; %數據長度和采樣頻率
n=0:N-1;t=n/Fs;   %時間序列
Lag=100;   %延遲樣點數
randn('state',0);   %設置產生隨機數的初始狀態
x=sin(2*pi*10*t)+0.6*randn(1,length(t));   %原始信號
[c,lags]=xcorr(x,Lag,'unbiased');     %對原始信號進行無偏自相關估計
subplot(2,2,1),plot(t,x);           %繪原始信號x
xlabel('時間/s');ylabel('x(t)');title('帶噪聲周期信號');grid on;
subplot(2,2,2);plot(lags/Fs,c);     %繪x信號自相關,lags/Fs為時間序列
xlabel('時間/s');ylabel('Rx(t)');
title('帶噪聲周期信號的自相關');grid on;
%信號 x1
x1=randn(1,length(x));   %產生一與x長度一致的隨機信號
[c,lags]=xcorr(x1,Lag,'unbiased'); %求隨機信號x1的無偏自相關
subplot(2,2,3),plot(t,x1);  %繪制隨機信號x1 
xlabel('時間/s');ylabel('x1(t)');title('噪聲信號');
grid on;
subplot(2,2,4);plot(lags/Fs,c); %繪制隨機信號x1的無偏自相關
xlabel('時間/s');ylabel('Rx1(t)');
title('噪聲信號的自相關');
grid on