5-5.m

来自「经典数字信号处理滤波器的源代码 重点是利用巴特沃斯模拟滤波器转而设计其它数字滤」· M 代码 · 共 23 行

%例程5－5  计算随机信号的自相关函数
% e.g.5-5.m for example5-5;
% to test function of xcorr;
 
clear all;
N=10^3;
n=[0:N-1];
Fs=500;
t=n/Fs;
Mlag=100;
% Generate signal plus noise
x1=sin(2*pi*10*t)+0.7*randn(1,length(t));
figure(1);  plot(t,x1);
xlabel('n'); ylabel('x1(n)');
title('Original signal x1');
grid on;
% Computing the unbiased atuocorretaion of x1
[r1,lags]=xcorr(x1,Mlag,'unbiased');
figure(2); plot(lags/Fs,r1);
xlabel('n'); ylabel('Rx1(n)');
title('Autocorrelation of x1');
grid on;