autocorrelation.m

AR模型的源程序。

% File: autocorrelation.m
% ------------------------------------
% This function calculate the autocorrelation of x(n)

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