crosscorrelation.m

AR模型的源程序。

% File: crosscorrelation.m
% ------------------------
% This function calculate the crosscorrelation of x(n) and s(n)

function[crosscorrelation] = crosscorrelation(m, signalvector_x, signalvector_s);
% L: Sample length
% m: Timeshift interval
% signalvector_x: Input sequence
% signalvector_s: Signal sequence
L = length(signalvector_x);
crosscorrelation = 0;
for i = 1: L - abs(m)
    crosscorrelation = signalvector_x(i) * signalvector_s(i + abs(m)) + crosscorrelation;
end
crosscorrelation = crosscorrelation / (L - abs(m));