Міністерство освіти і науки України Національний університет „Львівська політехніка” Кафедра ЕОМ Звіт з лабораторної роботи №4 на тему: „ Швидкі алгоритми обчислень дискретних тригонометричних перетворень.” [ім’я автора приховано] Прийняв: Кузьо М. М. Львів 2007


Результат:
uses Windows, Messages, SysUtils, Variants, Classes, Graphics, Controls, Forms, Dialogs, StdCtrls, Buttons, jpeg, ExtCtrls; type TForm1 = class(TForm) BitBtn1: TBitBtn; ListBox1: TListBox; ListBox2: TListBox; BitBtn2: TBitBtn; Image1: TImage; procedure BitBtn1Click(Sender: TObject); procedure BitBtn2Click(Sender: TObject); procedure FormCreate(Sender: TObject); private { Private declarations } public { Public declarations } end; var Form1: TForm1; y:real; implementation {$R *.dfm} procedure TForm1.BitBtn1Click(Sender: TObject); const pi = 3.1415926535897932384626433832795; N = 512; var x: array [0..300] of real; i,n1,k: integer; tt: string; begin //for i:=0 to 128 do begin k := 0; n1 := 0; x[0] := 1/(N/2) * (sin(2*pi*k*n1/N)); tt := FloatToStr(x[0]); ListBox1.Items.Add(' | При x(0) | ' + tt); k := 1; n1 := 1; x[1] := 1/(N/2)*( x[0] + (sin(2*pi*k*n1/N))); tt := FloatToStr(x[1]); ListBox1.Items.Add(' | При x(1) | ' + tt); k := 2; n1 := 2; x[2] := 1/(N/2)*( x[0] + x[1] + (sin(2*pi*k*n1/N))); tt := FloatToStr(x[2]); ListBox1.Items.Add(' | При x(2) | ' + tt); k := 3; n1 := 3; x[3] := 1/(N/2)*( x[0] + x[1] + x[2] + (sin(2*pi*k*n1/N))); tt := FloatToStr(x[3]); ListBox1.Items.Add(' | При x(3) | ' + tt); k := 4; n1 := 4; x[4] := 1/(N/2)*( x[0] + x[1] + x[2] + x[3] + (sin(2*pi*k*n1/N))); tt := FloatToStr(x[4]); ListBox1.Items.Add(' | При x(4) | ' + tt); k := 5; n1 := 5; x[5] := 1/(N/2)*( x[0] + x[1] + x[2] + x[3] + x[4] + (sin(2*pi*k*n1/N))); tt := FloatToStr(x[5]); ListBox1.Items.Add(' | При x(5) | ' + tt); k := 6; n1 := 6; x[6] := 1/(N/2)*( x[0] + x[1] + x[2] + x[3] + x[4] + x[5] + (sin(2*pi*k*n1/N))); tt := FloatToStr(x[6]); ListBox1.Items.Add(' | При x(6) | ' + tt); k := 7; n1 := 7; x[7] := 1/(N/2)*( x[0] + x[1] + x[2] + x[3] + x[4] + x[5] + x[6] + (sin(2*pi*k*n1/N))); tt := FloatToStr(x[7]); ListBox1.Items.Add(' | При x(7) | ' + tt); k := 8; n1 := 8; x[8] := 1/(N/2)*(x[0]+x[1]+x[2]+x[3]+x[4]+x[5]+x[6]+x[7]+(sin(2*pi*k*n1/N))); tt := FloatToStr(x[8]); ListBox1.Items.Add(' | При x(8) | ' + tt); k := 9; n1 := 9; x[9] := 1/(N/2)*(x[0]+x[1]+x[2]+x[3]+x[4]+x[5]+x[6]+x[7]+x[8]+(sin(2*pi*k*n1/N))); tt := FloatToStr(x[9]); ListBox1.Items.Add(' | При x(9) | ' + tt); k := 10; n1 := 10; x[10] := 1/(N/2)*(x[0]+x[1]+x[2]+x[3]+x[4]+x[5]+x[6]+x[7]+x[8]+x[9]+(sin(2*pi*k*n1/N))); tt := FloatToStr(x[10]); ListBox1.Items.Add(' | При x(10) | ' + tt); k := 11; n1 := 11; x[11] := 1/(N/2)*(x[0]+x[1]+x[2]+x[3]+x[4]+x[5]+x[6]+x[7]+x[8]+x[9]+x[10]+(sin(2*pi*k*n1/N))); tt := FloatToStr(x[11]); ListBox1.Items.Add(' | При x(11) | ' + tt); k := 12; n1 := 12; x[12] := 1/(N/2)*(x[0]+x[1]+x[2]+x[3]+x[4]+x[5]+x[6]+x[7]+x[8]+x[9]+x[10]+x[11]+(sin(2*pi*k*n1/N))); tt := FloatToStr(x[12]); ListBox1.Items.Add(' | При x(12) | ' + tt); end; procedure TForm1.BitBtn2Click(Sender: TObject); const pi = 3.1415926535897932384626433832795; N = 512; var x: array [0..300] of real; i,n1,k: integer; tt: string; begin BitBtn2.Enabled := False; for i:=0 to 127 do begin x[i] := 1/(N/2) * (y + cos(2*pi*i*i/N)); y := x[i]; tt := FloatToStr(y); ListBox2.Items.Add(' | При ('+ IntToStr(i) +') | ' + tt); end; end; procedure TForm1.FormCreate(Sender: TObject); begin y:=0; end; end. Висновок: на даній лабораторній роботі я дослідив швидкі алгоритми дискретних тригонометричних перетворень та порівняв їх з алгоритмами безпосереднього обчислення тригонометричних перетворень за допомогою даної програми.