Aмплитудно-Частотная Характеристика, delphi 7

построить АЧХ характеристику сигнала, используя ряд фурье 
 
procedure TForm1.Button1Click(Sender: TObject);
const
M=15;
Var
a,c,b,f: array [1..M] of real;
t,n:integer;
w,dt: real;
begin
form1.Series1.Clear;
dt:=StrtoFloat(Edit1.text);
 
 
for n:=1 to M do
  begin
 a[n]:=0;
  b[n]:=0;
 
 
  for t:=1 to M do
 begin
 w:=2*pi/t+dt;
    f[t]:=sin((5/20)*t)+sin((2/20)*t);
  a[n]:=a[n]+(f[t]*cos(n*w*(t+dt)))*dt;
     b[n]:=b[n]+(f[t]*sin(n*w*(t+dt)))*dt;
 
 
 
 
   c[n]:=sqrt(a[n]*a[n]+b[n]*b[n]);
    form1.Series1.AddXY(w,c[n]);
end;
 end;
   end;

procedure TForm1.FormPaint(Sender: TObject);
var x0, y0: Integer;
    x,y: Real;
    m: Integer;
    A: Array of TPoint;
    k,i: Integer;
begin
  x0 := 0;
  y0 := ClientHeight div 2;
  x := 0;
  m := 40;
  y := 0;
  k := 0;
 
  while x <= 20 do
    begin
      y := Cos(x)+1;
      //Canvas.Pixels[Round(x0+x*m), Round(y0-y*m)] := clGreen;
      SetLength(A,k+1);
      A[k].X := Round(x0+x*m);
      A[k].Y := Round(y0-y*m);
      x := x+0.1;
      k := k+1;
    end;
 
  for i := 0 to Length(A)-1 do
    begin
      Canvas.MoveTo(A[i].X, y0);
      Canvas.LineTo(A[i].X, A[i].Y);
    end;
end;

 for n:=1 ti 20 do
 begin 
  for t:=1 to N do

unit Unit1;
 
interface
 
uses
 Windows, Messages, SysUtils, Variants, Classes, Graphics, Controls, Forms,
 Dialogs, TeEngine, Series, ExtCtrls, TeeProcs, Chart, StdCtrls;
 
type
 TSingleArray = array[0..10000] of extended; //Single
 TForm1 = class(TForm)
 Button1: TButton;
 Chart1: TChart;
 Chart2: TChart;
 Series3: TLineSeries;
 Series4: TLineSeries;
 Series1: TLineSeries;
 Series2: TLineSeries;
 Procedure DTF(var Array_Re,Array_Im:TSingleArray;size:LongInt);
 procedure Button1Click(Sender: TObject);
 private
 { Private declarations }
 public
 { Public declarations }
 end;
 
var
 Form1: TForm1;
 Array_Re,Array_Im:TSingleArray;
 implementation
 
{$R *.dfm}
 
procedure TForm1.Button1Click(Sender: TObject);
 var
 i,n: integer;
 
begin
 n:=512;
 Series1.Clear;
 Series2.Clear;
 Series3.Clear;
 Series4.Clear;
 for i := 0 to n - 1 do
 begin
 Array_Re[i]:=sin(i)+1;
 Array_Im[i]:=0;
 Series1.AddY(Array_Re[i]);
 Series2.AddY(Array_Im[i]);
 end;
 DTF(Array_Re,Array_Im,n);
 for i := 0 to n - 1 do
 begin
 Series3.AddY(Array_Re[i]);
 Series4.AddY(Array_Im[i]);
 end;
 end;
 
Procedure TForm1.DTF(var Array_Re,Array_Im:TSingleArray;size:LongInt);
 Var
 t,f : LongInt;
 k,norm : Single;
 Temp_Array_Re,Temp_Array_Im : Array Of Single;
 Begin
 k:=2*pi/size;
 norm:=1/sqrt(size);
 SetLength(Temp_Array_Re,size);
 SetLength(Temp_Array_Im,size);
 For t:=0 To size-1 Do
 begin
 Temp_Array_Re[t]:=Array_Re[t];
 Temp_Array_Im[t]:=Array_Im[t];
 end;
 For t:=0 To size-1 Do
 Begin
 Array_Re[t]:=0;
 Array_Im[t]:=0;
 For f:=0 To size-1 Do
 begin
 Array_Re[t]:=Array_Re[t]+Temp_Array_Re[f]*cos(k*f*t);
 Array_Im[t]:=Array_Im[t]+Temp_Array_Re[f]*sin(k*f*t);
 end;
 Array_Re[t]:=norm*Array_Re[t];
 Array_Im[t]:=norm*Array_Im[t];
 End;
 End;
 end.

procedure TForm1.Button1Click(Sender: TObject);
const
M=20;
 
Var
a,c,b,f: array [1..M] of real;
 t,n:integer;
  w,dt:real;
begin
   N:=200;
form1.Series1.Clear;
 dt:=StrtoFloat(Edit1.text);
  w:=2*pi/M;
 
 
 for  t:= 0 to N  do
 begin
  f[t]:=Sin(t)+1;
 end;
 
 
 
{ for n:=0 to m-1 do
 f[t]:=fy[t] }
 
 for n:=0 to M do
  begin
  a[n]:=0;
  b[n]:=0;
 for t:=0 to M do
   begin
 
  a[n]:=a[n]+(f[t]*cos(n*w*t))*dt;
  b[n]:=b[n]+(f[t]*sin(n*w*t))*dt;
  c[n]:=sqrt(a[n]*a[n]+b[n]*b[n]);
  form1.Series1.AddXY(w,c[n]);
  end;
   End;
 End;

Form1.Canvas.Pixels[w,c[n]] := clBlack;

procedure TForm1.Button1Click(Sender: TObject);
const
size=50;
 
var
 f,a,b,c: array [1..size] of real;
  w: real ;
   dt:real;
   i,n: integer;
    t,o : integer;
 
begin
 n:=200;
  Series1.Clear;
  dt:=StrtoFloat(Edit1.text);
  for i := 0 to n - 1 do
 begin
 f[i]:=Sin(1*i*pi/180)+1;
 end;
 
  w:=2*pi/size;
 
 
 For t:=0 To size-1 Do
  Begin
   a[t]:=0;
    b[t]:=0;
 For i:=0 To size-1 Do
   begin
    a[t]:=a[t]+(f[i]*cos(w*i*t))*dt;
     b[t]:=b[t]+(f[i]*sin(w*i*t))*dt;
      c[t]:=sqrt(a[t]*a[t]+b[t]*b[t]);
       Series1.AddY(w,c[t]);--- выдилина красным
 
   end;
  end;
end;
 end.

Series1.AddY(Round(w),Round(c[t]));

Series1.AddY(Round(w),Round(c[t]));

Series1.AddY

Series1.AddХY