tmp602.m

来自「信号分析与处理的一些源代码」· M 代码 · 共 43 行

%MATLAB PROGRAM 6-2
clf

load dat1;
N=4000;
n=0:N-1;
Fs=1000;
t=n/Fs;
Lag=200;

%signal x
%========
x=sin(2*pi*10*t)+0.6*randn(1,length(t));
[c,lags]=xcorr(x,Lag,'unbiased');

subplot(221)
plot(t,x)
xlabel('t');ylabel('x(t)');
title('Original signal x');
grid

subplot(222)
plot(lags/Fs,c)
xlabel('t');ylabel('Rx(t)');
title('Autocorrelation ');
grid

%signal x1
%==========
x1=randn(1,length(x));
[c,lags]=xcorr(x1,Lag,'unbiased');

subplot(223)
plot(t,x1)
xlabel('t');ylabel('x(t)');
title('Original signal x1');
grid

subplot(224)
plot(lags/Fs,c)
xlabel('t');ylabel('Rx1(t)');
title('Autocorrelation ');
grid