unit Unit1;
interface
uses
Windows, Messages, SysUtils, Variants, Classes, Graphics, Controls, Forms,
Dialogs, StdCtrls, Grids;
type
TForm1 = class(TForm)
StringGrid1: TStringGrid;
Edit1: TEdit;
Edit2: TEdit;
Button1: TButton;
Button2: TButton;
Label2: TLabel;
lbl1: TLabel;
Edit3: TEdit;
procedure Button2Click(Sender: TObject);
procedure Button1Click(Sender: TObject);
private
{ Private declarations }
public
{ Public declarations }
end;
var
Form1: TForm1;
a:array [0..21,0..21] of real;
implementation
{$R *.dfm}
procedure sort();
var
i:integer;
buf:real;
flag:boolean;
begin
flag:=true;
while flag do
begin
flag:=false;
for i:=0 to 20 do
if a[i,0]>a[i+1,0] then
begin
buf:=a[i+1,0];
a[i+1,0]:=a[i,0];
a[i,0]:=buf;
buf:=a[i+1,1];
a[i+1,1]:=a[i,1];
a[i,1]:=buf;
flag:=true;
end;
end;
end;
function interpol(p1,p2:integer; x0:real):real;
var
Y1,Y2,X1,X2,y0,k,c:real;
begin
X1:=a[p1,0];
Y1:=a[p1,1];
X2:=a[p2,0];
Y2:=a[p2,1];
k:=(Y2-Y1)/(X2-X1);
c:= Y1-k*X1;
y0:=k*x0+c;
interpol:=(y0);
end;
function interpolkv(p1,p2,p3:integer; x0:real):real;
var
Y1,Y2,Y3,X1,X2,X3,yk,d,h,v:real;
begin
X1:=a[p1,0];
Y1:=a[p1,1];
X2:=a[p2,0];
Y2:=a[p2,1];
X3:=a[p3,0];
Y3:=a[p3,1];
d:=((Y3-Y1)*(X2-X1)-(Y2-Y1)*(X3-X1))/(Sqr(X3)-Sqr(X1))*(X2-X1)-(Sqr(X2)-Sqr(X1))*(X3-X1);
h:=(Y2-Y1-d*(Sqr(X2)-Sqr(X1)))/(X2-X1);
v:=Y1-(d*Sqr(X1)+h*X1);
yk:=d*Sqr(x0)+h*x0+v;
interpolkv:=(yk);
end;
Function find(x0:Real):integer;
var
i:integer;
begin
i:=1;
while a[i,0]>x0 do
i:=i+1;
find:=(i);
end;
procedure TForm1.Button2Click(Sender: TObject);
var
i:integer;
begin
{Randomize;}
StringGrid1.Cells[0,0]:='Npk';
StringGrid1.Cells[0,1]:='X';
StringGrid1.Cells[0,2]:='y';
StringGrid1.ColCount:=23;
For i:=1 to StringGrid1.ColCount do
begin
StringGrid1.Cells[i,1]:=FloatToStr(5-(i-1)*0.5);
StringGrid1.Cells[i,0]:=IntToStr(i);
StringGrid1.Cells[i,2]:=FloatToStr(Random(10));
end;
begin
Edit1.Text:='2,5';
Edit2.Text:='';
Edit3.Text:='';
end;
end;
procedure TForm1.Button1Click(Sender: TObject);
var
i,j,b:integer;
x0,y0,yk:real;
begin
for i:=0 to 20 do
begin
for j:=0 to 1 do
a[i,j]:=StrToFloat(StringGrid1.Cells[i+1,j+1]);
end;
sort;
for i:=0 to 20 do
begin
for j:=0 to 1 do
StringGrid1.Cells[i+1,j+1]:=FloatToStr(a[i,j]);
end;
x0:=StrToFloat(Edit1.Text);
b:=find(x0);
y0:=interpol(b-1,b,x0);
Edit2.Text:=FloatToStr(y0);
yk:=interpolkv(b-1,b,x0);
Edit3.Text:=FloatToStr(yk);
end;
end.