n = 1024;
t = 1:n;
noisey1 = rand(1,n);
noisey = noisey1 - mean(noisey1);
data = 0.2*sin(0.1*t);
total = noisey + data;
data2 = 0.2*sin(0.1*t + 10);
total2 = noisey + data2;
M = moviein(50);
for j = 1:50;
 dt = 1*(j-1);
subplot(2,1,1);
plot(t,total,'o', dt + t(1:n-dt),total2(1:n-dt),'x');
subplot(2,1,2);
crosscor = sum(total(1:n-dt).*total2(dt:n-1))/sum(total(1:n).*total2(1:n));
plot(dt,crosscor,'x');
axis([0 200 -1 1]);
  M(:,j) = getframe;
end;
movie(M);