{ **** UBPFD *********** by delphibase.endimus.com ****
>> Вычисление автокорреляционной функции
Зависимости: Math
Автор: lookin, <a href="mailto:lookin@mail.ru">lookin@mail.ru</a>, Екатеринбург
Copyright: lookin
Дата: 2 июня 2002 г.
***************************************************** }
//TDoubleArray = array of double
procedure AutoCorrelation(ValueArray: TDoubleArray; var ACArray: TDoubleArray;
 FromValue, ToValue: integer);
var
 i, j, N: integer;
 avr, dev, xxsum: double;
begin
 //ValueArray - массив типа double для которого вычисляется функция
 //FromValue - номер точки, начиная с которого выбираются элементы массива
 //ToValue - номер точки, на котором заканчивается выбор элементов массива
 //ACArray - массив возвращаемых значений автокорреляционной функции
 //для 5-и точек
 N := ToValue - FromValue;
 SetLength(ACArray, 5);
 if N < 5 then
 begin
  for i := 0 to 4 do
  ACArray[i] := 0;
  Exit;
 end
 else
 begin
  SetLength(rv, N);
  dev := 0;
  for i := 0 to N - 1 do
  rv[i] := ValueArray[i + FromValue];
  avr := Mean(rv);
  for i := 0 to N - 1 do
  dev := dev + Sqr(rv[i] - avr);
  dev := dev / N;
  for j := 0 to 4 do
  begin
  xxsum := 0;
  for i := 0 to (N - 1) - j do
  xxsum := xxsum + (rv[i] - avr) * (rv[i + j] - avr);
  ACArray[j] := xxsum / (dev * (N - j));
  end;
 end;
end;
//Пример использования:
var
 SourceArray, ACCoefs: TDoubleArray;
begin
 AutoCorrelation(SourceArray, ARCoefs, 0, Length(SourceArray) - 1);
 for i := 0 to Length(ACCoefs) - 1 do
  showmessage(FloatToStr(ACCoefs[i]));