prog21.m

Cross-correlation of several signals

clear all;
clg;
fp=2000;
N=1000;
f1=100;
f2=101;
f3=340;
t=linspace(0,(N-1)/fp,N);
x1=sin(2*pi*t*f1);
x2=sin(2*pi*t*f2);
x3=sin(2*pi*t*f3);
x4=rand(1,N)-0.5;
x5=randn(1,N);
x6=randn(1,N);
x7=square(2*pi*t*f1);
x8=sawtooth(2*pi*t*f1);
l=700;
tc=linspace(-l/fp,l/fp,2*l+1);
%dwa sinusy o r罂nych cz阺totliwo渃iach
%im wi阫sza r罂nica cz阺totliwo渃i tym mniejsza warto滄 max korelacji
c1=xcorr(x1,x2,l,'coeff');
subplot(4,2,1);
plot(tc,c1);
title('korelacja wzajemna sin 100Hz i 101Hz')
xlabel('czas[s]')
ylabel('wartosc')
c2=xcorr(x1,x3,l,'coeff');
subplot(4,2,2);
plot(tc,c2);
title('korelacja wzajemna sin 100Hz i 340Hz')
xlabel('czas[s]')
ylabel('wartosc')
%szum r體nomierny z gausowskim
c3=xcorr(x4,x5,l,'coeff');
subplot(4,2,3);
plot(tc,c3);
title('korelacja wzajemna szumu rownomiernego z gausowskim')
xlabel('czas[s]')
ylabel('wartosc')
%dwie realizacje szumu gausowskiego
c4=xcorr(x5,x6,l,'coeff');
subplot(4,2,4);
plot(tc,c4);
title('korelacja wzajemna 2 realizacji szumu gausowskiego')
xlabel('czas[s]')
ylabel('wartosc')
%sinus i szum gausowski
c5=xcorr(x1,x5,l,'coeff');
subplot(4,2,5);
plot(tc,c5);
title('korelacja wzajemna sinusa 100Hz i szumu gausowskiego')
xlabel('czas[s]')
ylabel('wartosc')
%sinus i szum r體nomierny
c6=xcorr(x1,x4,l,'coeff');
subplot(4,2,6);
plot(tc,c6);
title('korelacja wzajemna sinusa 100Hz i szumu rownomiernego')
xlabel('czas[s]')
ylabel('wartosc')
%sinus z prostok箃em
c7=xcorr(x1,x7,l,'coeff');
subplot(4,2,7);
plot(tc,c7);
title('korelacja wzajemna sin 100Hz z prostokatem 100Hz')
xlabel('czas[s]')
ylabel('wartosc')
%sinus z tr骿k箃em o tej samej cz阺totliow渃i
c8=xcorr(x1,x8,l,'coeff');
subplot(4,2,8);
plot(tc,c8);
title('korelacja wzajemna sin 100Hz z trojkatem 100Hz')
xlabel('czas[s]')
ylabel('wartosc')