Фильтрация цифрового сигнала

 procedure TForm1.Button1Click(Sender: TObject);
  var s_tek, j: integer;
  begin
  i:=0;
  Series1.Clear;
  with OpenDialog1 do
  if Execute then
    begin
      assignfile(t, FileName);
      Reset(t);
      while not EOF(t) do
             begin
             Readln(t, s_tek);
             try
             Series1.AddXY(i/200,s_tek,'',clRed);
             except end;
             i:=i+1;
             end;
      closefile(t);
      end;
      n:=i-1;
      SetLength(a,i);
      for j:=0 to n do
      begin
        a[j]:=Series1.YValues[j] // собственно, вот этот массив который надо отфильтровать
      end;
      end;

procedure TForm1.Button3Click(Sender: TObject); // фильтрация - адаптивный фильтр
var i:integer;
S,S1,t,out,w1,w2:array of extended;
mu: double;
begin
setlength(S, Length(a));
setlength(t, Length(a));
setlength(Out, Length(a));
setlength(S1, Length(a));
setlength(w1, Length(a));
setlength(w2, Length(a));
Series2.Clear;
    mu:=1e-11;    // коэффициент определяющий глубину адаптации
    w1[0]:=0.5;   // изначальное значение весового коэффициента 1
    w2[0]:=0.5;   // изначальное значение весового коэффициента 2
    Out[0]:=a[0];
    for i:=0 to length(a)-1 do
    begin
    S[i]:=1500*sin(2*3.14*50*i);     // эталонный сигнал
    S1[i]:=1500*cos(2*3.14*50*i);   // эталонный сигнал сдвинутый на 90
    w1[i]:=w1[i]+2*mu*Out[i]*S[i]; // весовой коэффициент 1 по текущему отсчету
    w2[i]:=w2[i]+2*mu*Out[i]*S1[i]; // весовой коэффициент 2 по текущему отсчету
    Out[i]:=a[i]-(w1[i]*S[i]+w2[i]*S1[i]); // выходной сигнал
    w1[i+1]:=w1[i]+2*mu*Out[i]*S[i];    // перерасчет весового коэффициента 1
    w2[i+1]:=w2[i]+2*mu*Out[i]*S1[i];   // перерасчет весового коэффициента 2
    Series2.AddXY(i/200,Out[i],'',clRed);
    end;