Main Content

Esta página aún no se ha traducido para esta versión. Puede ver la versión más reciente de esta página en inglés.

Correlación cruzada con entrada multicanal

Generar tres secuencias exponenciales de 11 muestras dadas por

<math display="block">
<mrow>
<mn>0</mn>
<mo>.</mo>
<msup>
<mrow>
<mn>4</mn>
</mrow>
<mrow>
<mi>n</mi>
</mrow>
</msup>
</mrow>
</math>
,
<math display="block">
<mrow>
<mn>0</mn>
<mo>.</mo>
<msup>
<mrow>
<mn>7</mn>
</mrow>
<mrow>
<mi>n</mi>
</mrow>
</msup>
</mrow>
</math>
Y
<math display="block">
<mrow>
<mn>0</mn>
<mo>.</mo>
<mn>9</mn>
<mn>9</mn>
<msup>
<mrow>
<mn>9</mn>
</mrow>
<mrow>
<mi>n</mi>
</mrow>
</msup>
</mrow>
</math>
Con
<math display="block">
<mrow>
<mi>n</mi>
<mo></mo>
<mn>0</mn>
</mrow>
</math>
. Se utiliza para trazar las secuencias una al lado de la otra.stem3

N = 11; n = (0:N-1)';  a = 0.4; b = 0.7; c = 0.999;  xabc = [a.^n b.^n c.^n];  stem3(n,1:3,xabc','filled') ax = gca; ax.YTick = 1:3; view(37.5,30)

Calcular las autocorrelaciones y correlaciones cruzadas mutuas de las secuencias. Salida de los retrasos para que no tenga que realizar un seguimiento de ellos. Normalizar el resultado para que las autocorrelaciones tengan el valor unitario en cero retardo.

[cr,lgs] = xcorr(xabc,'coeff');  for row = 1:3     for col = 1:3         nm = 3*(row-1)+col;         subplot(3,3,nm)         stem(lgs,cr(:,nm),'.')         title(sprintf('c_{%d%d}',row,col))         ylim([0 1])     end end

Restringir el cálculo a los retrasos entre

<math display="block">
<mrow>
<mo>-</mo>
<mn>5</mn>
</mrow>
</math>
Y
<math display="block">
<mrow>
<mn>5</mn>
</mrow>
</math>
.

[cr,lgs] = xcorr(xabc,5,'coeff');  for row = 1:3     for col = 1:3         nm = 3*(row-1)+col;         subplot(3,3,nm)         stem(lgs,cr(:,nm),'.')         title(sprintf('c_{%d%d}',row,col))         ylim([0 1])     end end

Calcular estimaciones imparciales de las autocorrelaciones y correlaciones cruzadas mutuas. De forma predeterminada, los retrasos se ejecutan entre

<math display="block">
<mrow>
<mo>-</mo>
<mo stretchy="false">(</mo>
<mi>N</mi>
<mo>-</mo>
<mn>1</mn>
<mo stretchy="false">)</mo>
</mrow>
</math>
Y
<math display="block">
<mrow>
<mi>N</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</math>
.

cu = xcorr(xabc,'unbiased');  for row = 1:3     for col = 1:3         nm = 3*(row-1)+col;         subplot(3,3,nm)         stem(-(N-1):(N-1),cu(:,nm),'.')         title(sprintf('c_{%d%d}',row,col))     end end

Consulte también

Funciones