Додаток Б. Лістинг програми
Program interpol;
Uses Crt,Graph;
Const t1=75; t2=93;
T: array [0. .5] of integer = (20,40,60,80,100,120);
U: array [0. .5] of real = (-0.67,-0.25,-0.17,0.61,1.06,1.52);
Var
A1: array [0. .5] of real;
A2: array [0. .5] of real;
A3: array [0. .5] of real;
A4: array [0. .5] of real;
A5: array [0. .5] of real;
driver,mode, i,j: integer;
x,x1,dx,y,y1,dy0,dy1,dy2,dy3,dy4,dy5,dy6,dy7: integer;
e,PP,P1,P2,Z1,Z2,ZZ: real; Ch: char;
LL0,LL1,LL2,LL3,LL4,LL5: real;
{----------------------------------------------------------------------------------------}
{TABLUCIA}
procedure tabluca;
begin
TextBackGround (white);
ClrScr;
TextColor (red);
Writeln ('');
GoToXY (20,10);
Writeln ('TARUROVO4NA TABLUCIA DLIA TERMOPARU');
GoToXY (15,12);
Writeln ('------------------------------------------------');
GoToXY (15,13);
Writeln ('| T,C | 20 | 40 | 60 | 80 | 100 | 120 |');
GoToXY (15,14);
GoToXY (15,15);
Writeln ('|U,mV | - 0.67| - 0.25| - 0.17| 0.61 | 1.06 | 1.52|');
GoToXY (15,16);
end;
{-------------------------Graphik-----------------------------------------------------}
procedure graphik (driver,mode,x,x1,dx,y,y1,dy0,dy1,dy2,dy3,dy4,dy5,dy6,dy7: integer);
driver: =detect;
initgraph (driver,mode,'C: \progra~1\TP\bgi');
setviewport (20, 20,getmaxx-59,getmaxy-59,true);
clearviewport;
setbkcolor (white);
x: =20; y: =20;
dy0: = (getmaxy div 2 - (trunc (80*u [0])));
dy1: = (getmaxy div 2 - (trunc (80*u [1])));
dy2: = (getmaxy div 2 - (trunc (80*u [2])));
dy3: = (getmaxy div 2 - (trunc (80*PP)));
dy4: = (getmaxy div 2 - (trunc (80*u [3])));
dy5: = (getmaxy div 2 - (trunc (80*ZZ)));
dy6: = (getmaxy div 2 - (trunc (80*u [4])));
dy7: = (getmaxy div 2 - (trunc (80*u [5])));
dx: =t [1] -t [0] ;
setcolor (red);
line (4*dx,dy0,6*dx,dy1);
line (6*dx,dy1,8*dx,dy2);
line (8*dx,dy2,6*dx+t1,dy3);
line (6*dx+t1,dy3,10*dx,dy4);
line (10*dx,dy4,7*dx+t2,dy5);
line (7*dx+t2,dy5,12*dx,dy6);
line (12*dx,dy6,14*dx,dy7);
repeat
setcolor (green);
line (20,y,getmaxx,y);
line (x, 20,x,getmaxy);
x: =x+20; y: =y+20;
y1: =getmaxy div 2;
setcolor (blue);
outtextxy (2*20,y1,'0 20 40 60 80 100 120');
outtextxy (20,y1-60*2,' 1.5');
outtextxy (20,y1-40*2,' 1.0');
outtextxy (20,y1-20*2,' 0.5');
outtextxy (20,y1+60*2,'-1.5');
outtextxy (20,y1+40*2,'-1.0');
outtextxy (20,y1+20*2,'-0.5');
until keypressed;
if readkey=#0 then driver: =ord (readkey);
closegraph
{---------------------ZASTAVKA--------------------------------------------------}
{---------------------------------------------------------------------------------------}
procedure zastavka;
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 ');
{-------------------------------------------------------------------}
{------------------OSNOVNA PROGRAMA------------------------------------}
zastavka;
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);
{-------METOD LAGRANJA-------------------------------------------------------}
#50: begin
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 ('Pru T=75C U=',PP: 3: 6,'mV');
e: =t2;
ZZ: =LL0+LL1+LL2+LL3+LL4+LL5;
writeln; GoToXY (25, 19);
writeln ('Pru T=93C U=',ZZ: 3: 6,'mV');
until ch=#27;
readln
end.
Додаток Г. Результати роботи програми
Страницы: 1, 2
