±³¼ö´Ô²²¼ ÁֽŠ±¸¹®À¸·Î ¸ÅÆ®·¦À» »ç¿ëÇغ¸¾Ò´Âµ¥
°á°ú°¡ ÀÌ·¸°Ô ³ª¿É´Ï´Ù. Àú±â 4¹ø° ±×·¡ÇÁ°¡ ±³¼ö´Ô²²¼ Á÷Á¢ º¸¿©ÁÖ½Ç ¶§¶û °ªÀÌ ´Ù¸£°Ô ³ª¿À´Âµ¥
Ȥ½Ã Input signal¿¡ ³ëÀÌÁî°¡ ¾ø¾î¼ CrossCorrelation ±×·¡ÇÁ¿¡µµ ³ëÀÌÁî°¡ ¾øÀÌ °á°ú°¡ ³ª¿À´Â °ÇÁö ¾Æ´Ï¸é ±¸¹®ÀÌ À߸øµÈ°ÇÁö
Áú¹® µå¸³´Ï´Ù..
Ã¥¿¡´Â
Àú·± ³ëÀÌÁî°¡ ÀÖ½À´Ï´Ù.
%
corrX.m CrossCorrelation (Pig4.9)
clear all;
L = 200; % delay
M = 10; % pulse width
N = 1000; % # of data
n = 0:1:N-1; % time index
f = [ones(1, M), zeros(1, N-M)] * 2;
fL = [zeros(1, L), f(1:N-L)];
w = randn(1, length(fL));
g = fL + w;
Rfg = xcorr(f, g);
subplot (2, 2, 1)
plot (n, f);
title('Input Signal: f(n)' );
axis ([0, N, -0.2, 2.2]);
ylabel('f(n)')
subplot(2, 2, 2)
plot(n, w);
title( 'Noise Signal: w(n)');
axis ( [0, N, -4, 4]);
ylabel('w(n)');
subplot (2, 2, 3)
plot (n, g);
title( 'Function: g(n) = f(n-100)+w(n)');
axis ([0, N, -4, 4]);
ylabel('g(n)');
subplot (2, 2, 4)
plot(n, Rfg(N:2*N-1));
title( 'CrossCorrelation of f(n)&g(n)' );
axis ([0, N, -20, 40]);
xlabel( 'Lag Index' ) ;
ylabel ('Rfg(1)');