Дослідження методів інтерполяції
Курсовой проект - Компьютеры, программирование
Другие курсовые по предмету Компьютеры, программирование
-------------------------------------------------------------------------}
procedure zastavka;
begin
clrscr; GoToXY (14,2);
Writeln ( (C) CopyRight by Alexander Kyrsta, BHTY, 1AM-03);
GoToXY (29,5); Writeln (PROGRAM INTERPOL, 2006);
Writeln; GoToXY (32,13); Writeln (VUBERITb METOD); Writeln;
GoToXY (24,15); Writeln (1-METOD ROZPODILENUX RIZNUCb. );
GoToXY (24,17); Writeln ( 2-METOD LAGRANJA. );
GotoXY (24,25); Writeln ( Press Esc+Enter for exit );
end;
{-------------------------------------------------------------------}
{------------------OSNOVNA PROGRAMA------------------------------------}
begin
zastavka;
repeat
ch: =readkey;
case ch of
{-------------METOD ROZPODILENUX RIZNUCb----------------------------}
#49: begin
{ROZPODILENI RIZNUCI 1 PORIADKY}
for i: =0 to 4 do begin
a1 [i]: = ( (u [i+1] -u [i]) / (t [i+1] -t [i]));
{writeln (a1 [i]: 3: 6); } end;
{ROZPODILENI RIZNUCI 2 PORIADKY}
for i: =0 to 3 do begin
a2 [i]: = ( (a1 [i+1] -a1 [i]) / (t [i+2] -t [i]));
{writeln (a2 [i]: 3: 7); } end;
{ROZPODILENI RIZNUCI 3 PORIADKY}
for i: =0 to 2 do begin
a3 [i]: = ( (a2 [i+1] -a2 [i]) / (t [i+3] -t [i]));
{writeln (a3 [i]: 3: 16); } end;
{ROZPODILENI RIZNUCI 4 PORIADKY}
for i: =0 to 1 do begin
a4 [i]: = ( (a3 [i+1] -a3 [i]) / (t [i+4] -t [i]));
{writeln (a4 [i]: 3: 16); } end;
{ROZPODILENI RIZNUCI 5 PORIADKY}
for i: =0 to 0 do begin
a5 [i]: = ( (a4 [i+1] -a4 [i]) / (t [i+5] -t [i]));
{writeln (a5 [i]: 3: 16); } end;
clrscr; tabluca;
P1: =u [0] +a1 [0] * (t1-t [0]) +a2 [0] * (t1-t [0]) * (t1-t [1]) +a3 [0] * (t1-t [0]) * (t1-t [1]) * (t1-t [2]);
P2: =a4 [0] * (t1-t [0]) * (t1-t [1]) * (t1-t [2]) * (t1-t [3]) +a5 [0] * (t1-t [0]) * (t1-t [1]) * (t1-t [2]) * (t1-t [3]) * (t1-t [4]);
PP: =P1+P2;
Z1: =u [0] +a1 [0] * (t2-t [0]) +a2 [0] * (t2-t [0]) * (t2-t [1]) +a3 [0] * (t2-t [0]) * (t2-t [1]) * (t2-t [2]);
Z2: =a4 [0] * (t2-t [0]) * (t2-t [1]) * (t2-t [2]) * (t2-t [3]) +a5 [0] * (t2-t [0]) * (t2-t [1]) * (t2-t [2]) * (t2-t [3]) * (t2-t [4]);
ZZ: =Z1+Z2;
writeln; GoToXY (25,17);
writeln (Pru T=75C U=,PP: 2: 6,mV);
Writeln; GoToXY (25,18);
writeln (Pru T=93C U=,ZZ: 2: 6,mV);
GoToXY (30,22);
writeln (PRESS ENTER);
readln;
graphik (x,x1,dx,y,y1,dy0,dy1,dy2,dy3,dy4,dy5,dy6,dy7,driver,mode);
end;
{-------METOD LAGRANJA-------------------------------------------------------}
#50: begin
clrscr; tabluca;
e: =t1;
LL0: =u [0] * (e-40) * (e-60) * (e-80) * (e-100) * (e-120) / ( (20-40) * (20-60) * (20-80) * (20-100) * (20-120));
LL1: =u [1] * (e-20) * (e-60) * (e-80) * (e-100) * (e-120) / ( (40-20) * (40-60) * (40-80) * (40-100) * (40-120));
LL2: =u [2] * (e-20) * (e-40) * (e-80) * (e-100) * (e-120) / ( (60-20) * (60-40) * (60-80) * (60-100) * (60-120));
LL3: =u [3] * (e-20) * (e-40) * (e-60) * (e-100) * (e-120) / ( (80-20) * (80-40) * (80-60) * (80-100) * (80-120));
LL4: =u [4] * (e-20) * (e-40) * (e-60) * (e-80) * (e-120) / ( (100-20) * (100-40) * (100-60) * (100-80) * (100-120));
LL5: =u [5] * (e-20) * (e-40) * (e-60) * (e-80) * (e-100) / ( (120-20) * (120-40) * (120-60) * (120-80) * (120-100));
PP: =LL0+LL1+LL2+LL3+LL4+LL5;
writeln; GoToXY (25,17);
writeln (Pru T=75C U=,PP: 3: 6,mV);
e: =t2;
LL0: =u [0] * (e-40) * (e-60) * (e-80) * (e-100) * (e-120) / ( (20-40) * (20-60) * (20-80) * (20-100) * (20-120));
LL1: =u [1] * (e-20) * (e-60) * (e-80) * (e-100) * (e-120) / ( (40-20) * (40-60) * (40-80) * (40-100) * (40-120));
LL2: =u [2] * (e-20) * (e-40) * (e-80) * (e-100) * (e-120) / ( (60-20) * (60-40) * (60-80) * (60-100) * (60-120));
LL3: =u [3] * (e-20) * (e-40) * (e-60) * (e-100) * (e-120) / ( (80-20) * (80-40) * (80-60) * (80-100) * (80-120));
LL4: =u [4] * (e-20) * (e-40) * (e-60) * (e-80) * (e-120) / ( (100-20) * (100-40) * (100-60) * (100-80) * (100-120));
LL5: =u [5] * (e-20) * (e-40) * (e-60) * (e-80) * (e-100) / ( (120-20) * (120-40) * (120-60) * (120-80) * (120-100));
ZZ: =LL0+LL1+LL2+LL3+LL4+LL5;
writeln; GoToXY (25, 19);
writeln (Pru T=93C U=,ZZ: 3: 6,mV);
GoToXY (30,22);
writeln (PRESS ENTER);
readln;
graphik (x,x1,dx,y,y1,dy0,dy1,dy2,dy3,dy4,dy5,dy6,dy7,driver,mode);
end;
end;
zastavka;
until ch=#27;
readln
end.
Додаток Г. Результати роботи програми