[exestrain]

Пожалуйста, помогите с программой для построения графиков характеристик адаптивного фильтра.
Сам я далек от програмирования.
Нужно, чтобы программа строила графики АЧХ и ФЧХ фильтра.
Формулы для ФЧХ и АЧХ я для удобства разбил:
Код:

Код:

am:=(SQRT(SQR(1-COS(2*(Pi)*64/128))+SQR(SIN(2*(PI)*64/128))))/
(SQRT(SQR(1-c*COS(2*(PI)*64/128))+SQR(c*SIN(2*(PI)*64/128))));
q1:=1-COS(2*(Pi)*n/128);
q2:=1-c*COs(2*(Pi)*n/128);
s1:=SiN(2*(PI)*n/128);
s2:=s1*c;
a1 :=SQRT(SQR(q1)+SQR(s1));
a2:=SQRT(SQR(q2)+SQR(s2));
f1:=ArcTan(s1/q1);
f2:=ArcTan(s2/q2);

a[n,b]:=(a1/a2)/am;
f[n,b]:=f1-f2;

ФЧХ - это f[n,b]
АЧХ - a[n,b]
n-частота
c<1 (0,75...0,95)
Заранее спасибо. Если поможете, очень меня выручите

[exestrain]

Вот, что я сам смог написать:

Код:

var
 Form1: TForm1;
    d,am,n,c,q1 ,q2,s1 ,s2,a1 ,a2,f1 ,f2,z:real;
  canvas: TCanvas;

implementation

{$R *.DFM}

Function f(x:real):real;
begin
c:=0;
for d:=1 to 10 do
begin
c:=c+0.1;
am:=(SQRT(SQR(1-COS(2*(Pi)*64/128))+SQR(SIN(2*(PI)*64/128))))/
(SQRT(SQR(1-c*COS(2*(PI)*64/128))+SQR(c*SIN(2*(PI)*64/128))));
q1:=1-COS(2*(Pi)*x/128);
q2:=1-c*COS(2*(Pi)*x/128);
s1:=SiN(2*(PI)*x/128);
s2:=s1*c;
if n<>64 then
a1 :=SQRT(SQR(q1)+SQR(s1));
a2:=SQRT(SQR(q2)+SQR(s2));
f:=(a1/a2)/am;
end;
end;

procedure GrOfFunc;
var
    n1,n2:real;
 y1,y2:real;
 n:real;
 y:real;
 dn:real;
 l,b:integer;
 w,h:integer;
 mx,my:real;
 x0,y0:integer;

begin

 l:=10;
 b:=Form1.ClientHeight-20;
 h:=Form1.ClientHeight-40;
 w:=Form1.Width-40;

 n1:=0;
 n2:=25;
 dn:=0.01;


 y1:=f(n1);
 y2:=f(n1);
 n:=n1;
 repeat
  y := f(n);
  if y < y1 then y1:=y;
  if y > y2 then y2:=y;
  n:=n+dn;
 until (n>=n2);
 my:=h/abs(y2-y1);
 mx:=w/abs(n2-n1);
 x0:=l;
 y0:=b-Abs(Round(y1*my));

 with form1.Canvas do
 begin

 MoveTo(l,b);LineTo(l,b-h);
 MoveTo(x0,y0);LineTo(x0+w,y0);
 TextOut(l+5,b-h,FloatToStrF(y2,ffGeneral,6,3));
 TextOut(l+5,b,FloatToStrF(y1,ffGeneral,6,3));

 n:=n1;
 repeat
  y:=f(n);
  Pixels[x0+Round(n*mx),y0-Round(y*my)]:=clRed;
  n:=n+dn;
 until (n>=n2);
end;
end;

procedure TForm1.FormPaint(Sender: TObject);
begin
GrOfFunc;
end;


procedure TForm1.FormResize(Sender: TObject);
begin
form1.Canvas.FillRect(Rect(0,0,ClientWidth,ClientHeight));
GrOfFunc;
end;

end.

[Kapitoshka438]

Вообще-то cos(2*pi*64/128) = cos(pi) = -1
sin(2*pi*64/128) = cos(pi) = 0

[exestrain]

Это считается максимальное значение функции, по которому она нормируется

[exestrain]

Ну, да я уже исправил ошибки. спасибо