[exestrain]

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

Код:

1 2 3 4 5 6 7 8 9 10 11 12 13

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]

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

Код:

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96

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]

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