|
|
Регистрация | << Правила форума >> | FAQ | Пользователи | Календарь | Поиск | Сообщения за сегодня | Все разделы прочитаны |
|
Опции темы | Поиск в этой теме | Опции просмотра |
|
#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)*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; АЧХ - a[n,b] n-частота c<1 (0,75...0,95) Заранее спасибо. Если поможете, очень меня выручите |
#2
|
|||
|
|||
Вот, что я сам смог написать:
Код:
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. |
#3
|
|||
|
|||
Вообще-то cos(2*pi*64/128) = cos(pi) = -1
sin(2*pi*64/128) = cos(pi) = 0 |
#4
|
|||
|
|||
Это считается максимальное значение функции, по которому она нормируется
|
#5
|
|||
|
|||
Ну, да я уже исправил ошибки. спасибо
|