이 번역 페이지는 최신 내용을 담고 있지 않습니다. 최신 내용을 영문으로 보려면 여기를 클릭하십시오.

xcorr

상호상관

구문

r = xcorr(x,y)

r = xcorr(x)

r = xcorr(___,maxlag)

r = xcorr(___,scaleopt)

[r,lags] = xcorr(___)

설명

r = xcorr(x,y)는 2개의 이산시간 시퀀스의 상호상관을 반환합니다. 상호상관은 벡터 x와 벡터 y의 이동된(지연된) 복사본 사이의 유사성을 그 지연 시간의 함수로 측정합니다. x와 y의 길이가 다르면 이 함수는 더 짧은 벡터의 끝부분에 0을 추가하여 더 긴 벡터와 동일한 길이를 가지도록 합니다.

r = xcorr(x)는 x의 자기상관 시퀀스를 반환합니다. x가 행렬이면 r도 행렬이 되며, 반환된 행렬의 열에는 x의 열의 모든 조합에 대한 자기상관 및 상호상관 시퀀스가 포함됩니다.

r = xcorr(___,maxlag)는 위에 열거된 구문에 대한 지연 범위를 -maxlag에서 maxlag까지로 제한합니다.

r = xcorr(___,scaleopt)는 상호상관 또는 자기상관에 대한 정규화 옵션도 지정합니다. 'none'(디폴트 값)이 아닌 다른 옵션을 사용하려면 x와 y의 길이가 같아야 합니다.

[r,lags] = xcorr(___)은 상관이 계산되는 지연값도 반환합니다.

예제

두 벡터의 상호상관

벡터 x와, x를 오른쪽으로 5개 요소만큼 이동한 벡터 y를 만듭니다. x와 y의 상호상관 추정값을 계산하고 플로팅합니다. 가장 큰 스파이크는 x와 y의 요소가 정확히 일치하는 지연값 (-5)에서 발생합니다.

n = 0:15;
x = 0.84.^n;
y = circshift(x,5);
[c,lags] = xcorr(x,y);
stem(lags,c)

벡터의 자기상관

벡터 x의 자기상관 추정값을 계산하고 플로팅합니다. 가장 큰 스파이크는 x가 자신과 정확하게 일치하는 지연 0에서 발생합니다.

n = 0:15;
x = 0.84.^n;
[c,lags] = xcorr(x);
stem(lags,c)

정규화된 상호상관

단위 피크를 갖는 벡터 x와 벡터 y의 정규화된 상호상관을 계산하고 플로팅한 다음 최대 지연값 10을 지정합니다.

n = 0:15;
x = 0.84.^n;
y = circshift(x,5);
[c,lags] = xcorr(x,y,10,'normalized');
stem(lags,c)

입력 인수

x — 입력 배열

벡터 | 행렬 | 다차원 배열

입력 배열로, 벡터, 행렬, 다차원 배열 중 하나로 지정됩니다. x가 다차원 배열인 경우 xcorr은 모든 차원에서 열 방향 연산을 수행하고 각 자기상관과 상호상관을 행렬의 열로 반환합니다.

데이터형: single | double
복소수 지원 여부: 예

y — 입력 배열

벡터

입력 배열로, 벡터로 지정됩니다.

데이터형: single | double
복소수 지원 여부: 예

maxlag — 최대 지연값

정수 스칼라

최대 지연값으로, 정수 스칼라로 지정됩니다. maxlag를 지정할 경우 반환되는 상호상관 시퀀스의 범위는 -maxlag에서 maxlag까지입니다. maxlag를 지정하지 않으면 지연값 범위가 2N – 1이 되며, 여기서 N은 x와 y 중 더 긴 쪽의 길이입니다.

데이터형: single | double

scaleopt — 정규화 옵션

'none' (디폴트 값) | 'biased' | 'unbiased' | 'normalized' | 'coeff'

정규화 옵션으로, 다음 중 하나로 지정됩니다.

• 'none' — 스케일링되지 않은 원시 상호상관입니다. 'none'은 x와 y의 길이가 다른 경우 유일하게 유효한 옵션입니다.

• 'biased' — 상호상관의 편향 추정값입니다.

  $${\widehat{R}}_{xy,\text{biased}}(m)=\frac{1}{N}{\widehat{R}}_{xy}(m).$$

• 'unbiased' — 상호상관의 무편향 추정값입니다.

  $${\widehat{R}}_{xy,\text{unbiased}}(m)=\frac{1}{N-\left|m\right|}{\widehat{R}}_{xy}(m).$$

• 'normalized' 또는 'coeff' — 지연 0에서 자기상관이 1이 되도록 시퀀스를 정규화합니다.

  $${\widehat{R}}_{xy,\text{coeff}}(m)=\frac{1}{\sqrt{{\widehat{R}}_{xx}(0){\widehat{R}}_{yy}(0)}}{\widehat{R}}_{xy}(m).$$

출력 인수

r — 상호상관 또는 자기상관

벡터 | 행렬

상호상관 또는 자기상관으로, 벡터 또는 행렬로 반환됩니다.

x가 M × N 행렬인 경우 xcorr(x)는 x의 열의 자기상관과 상호상관을 포함하는 (2M – 1) × N² 행렬을 반환합니다. maxlag를 지정할 경우 r의 크기는 (2 × maxlag + 1) × N²입니다.

예를 들어, S에 3개의 열 $$\text{S}=\left(\begin{array}{ccc}{x}_1& x_2& x_3\end{array}\right)$$이 있다면 R = xcorr(S)의 결과는 다음과 같이 구성됩니다.

$$\text{R}=\left(\begin{array}{lllllllll}{R}_{x_1{x}_1}\hfill & R_{x_1{x}_2}\hfill & R_{x_1{x}_3}\hfill & R_{x_2{x}_1}\hfill & R_{x_2{x}_2}\hfill & R_{x_2{x}_3}\hfill & R_{x_3{x}_1}\hfill & R_{x_3{x}_2}\hfill & R_{x_3{x}_3}\hfill \end{array}\right).$$

lags — 지연값 인덱스

벡터

지연값 인덱스로, 벡터로 반환됩니다.

세부 정보

상호상관과 자기상관

xcorr의 결과는 두 랜덤 시퀀스 간의 상호상관 추정값으로 해석되거나 두 결정적 신호 간의 결정론적 상관으로 해석될 수 있습니다.

결합 정상 확률 과정 x_n과 y_n의 실제 상호상관 시퀀스는 다음과 같이 지정됩니다.

$$R_{xy}(m)=E\left\{x_{n+m}{y}_n^{*}\right\}=E\left\{x_n{y}_{n-m}^{*}\right\},$$

여기서 −∞ < n < ∞입니다. 별표는 켤레 복소수를 나타내고, E는 기대값 연산자입니다. xcorr은 실제로 무한 길이 확률 과정의 한 실현에서 유한 부분만 사용 가능하기 때문에 시퀀스만 추정할 수 있습니다.

기본적으로, xcorr은 정규화 없이 원시 상관을 계산합니다.

$${\widehat{R}}_{xy}(m)=\{\begin{array}{ll}{\displaystyle \sum _{n=0}^{N-m-1}{x}_{n+m}{y}_n^{\ast },}\hfill & m\ge 0,\hfill \\ {\widehat{R}}_{yx}^{*}(-m),\hfill & m<0.\hfill \end{array}$$

출력 벡터 c는 다음으로 지정되는 요소를 가집니다.

$$\text{c}(\text{m})={\widehat{R}}_{xy}(m-N),m=1,2,\dots ,2N-1.$$

일반적으로 상관 함수가 정확한 추정값을 생성하려면 정규화가 필요합니다. 입력 인수 scaleopt를 사용하여 상관의 정규화를 제어할 수 있습니다.

참고 문헌

[1] Buck, John R., Michael M. Daniel, and Andrew C. Singer. Computer Explorations in Signals and Systems Using MATLAB^®. 2nd Edition. Upper Saddle River, NJ: Prentice Hall, 2002.

[2] Stoica, Petre, and Randolph Moses. Spectral Analysis of Signals. Upper Saddle River, NJ: Prentice Hall, 2005.

확장 기능

tall형 배열

메모리에 담을 수 없을 정도로 많은 행을 가진 배열을 계산할 수 있습니다.

xcorr 함수는 tall형 배열을 지원하지만 다음과 같은 사용법 관련 참고 및 제한 사항이 있습니다.

• x 입력값은 tall형 열 벡터여야 합니다.

• y 입력값은 tall형이 아닌 벡터여야 합니다.

• 구문 xcorr(x)는 지원되지 않습니다.

• 옵션 scaleopt는 지원되지 않습니다.

• maxlag를 지정할 경우 maxlag <= max(numel(x),numel(y))-1을 충족해야 합니다.

• lags 출력값은 tall형 열 벡터로 반환됩니다.

자세한 내용은 tall형 배열 항목을 참조하십시오.

C/C++ 코드 생성

MATLAB® Coder™를 사용하여 C 코드나 C++ 코드를 생성할 수 있습니다.

사용법 관련 참고 및 제한 사항:

• size(x)에서 선행하는 1의 값(차원 길이 값이 1이 아닌 첫 번째 값 앞에 오는 차원 길이 값 1)은 모든 입력값 x에 대해 일정해야 합니다. x가 가변 크기이고 행 벡터인 경우 1×:여야 합니다. 런타임에 size(x,1) = 1인 경우에는 :×:이 될 수 없습니다.

스레드 기반 환경

MATLAB®의 backgroundPool을 사용해 백그라운드에서 코드를 실행하거나 Parallel Computing Toolbox™의 ThreadPool을 사용해 코드 실행 속도를 높일 수 있습니다.

이 함수는 스레드 기반 환경을 완전히 지원합니다. 자세한 내용은 스레드 기반 환경에서 MATLAB 함수 실행하기 항목을 참조하십시오.

GPU 배열

Parallel Computing Toolbox™를 사용해 GPU(그래픽스 처리 장치)에서 실행하여 코드 실행 속도를 높일 수 있습니다.

xcorr 함수는 GPU 배열을 완전히 지원합니다. GPU에서 이 함수를 실행하려면 입력 데이터를 gpuArray (Parallel Computing Toolbox)로 지정하십시오. 자세한 내용은 GPU에서 MATLAB 함수 실행하기 (Parallel Computing Toolbox) 항목을 참조하십시오.

예를 들어, 단일 x에서 gpuArray 객체를 만들고 정규화된 자기상관을 계산할 수 있습니다.

t = 0:0.001:10-0.001;
x = cos(2*pi*10*t) + randn(size(t));  
X = gpuArray(x);                       
[r,lags] = xcorr(X,200,'normalized');     
r = gather(r);

버전 내역

R2006a 이전에 개발됨

R2019a: Signal Processing Toolbox에서 MATLAB으로 이동됨

이전에는 xcorr 함수를 실행하려면 Signal Processing Toolbox™가 필요했습니다.

참고 항목

conv | corrcoef | cov | xcov