Скачиваний:
37
Добавлен:
01.05.2014
Размер:
1.29 Кб
Скачать
Program LineFit1;
const max = 20;
type ary = array[1..max] of real;

procedure linfit1(x,y: ary; var y_calc: ary; var a,b:real; n:integer);
{ fit a straight line (y_calc) through n sets of x and y pairs of points }
var i : integer;
sum_x,sum_y,sum_xy,sum_x2,
sum_y2,xi,yi,sxy,sxx,
syy: real;
begin { linfit }
sum_x:=0.0;
sum_y:=0.0;
sum_xy:=0.0;
sum_x2:=0.0;
sum_y2:=0.0;
for i:=1 to n do
begin
xi:=x[i];
yi:=y[i];
sum_x:=sum_x+xi;
sum_y:=sum_y+yi;
sum_xy:=sum_xy+xi*yi;
sum_x2:=sum_x2+xi*xi;
sum_y2:=sum_y2+yi*yi;
end;
sxx:=sum_x2-sum_x*sum_x/n;
sxy:=sum_xy-sum_x*sum_y/n;
syy:=sum_y2-sum_y*sum_y/n;
b:=sxy/sxx;
a:=((sum_x2*sum_y-sum_x*sum_xy)/n)/sxx;
for i:=1 to n do
y_calc[i]:=a+b*x[i]
end; { LINFIT }

var x,y,y_calc:ary;
i:integer;
a,b:real;
begin
{creating new mass x and y}
randomize;
writeln('Input Data: ');
for i:=1 to max do
begin
x[i] := i*10;
y[i] := i*10 + random(20) - 10;
writeln('x[',i:2,'] = ',x[i]:3:0,' y[',i:2,'] = ',y[i]:3:0);
end;
linfit1(x,y,y_calc,a,b,max);
writeln('a: ',a:3:1,' b: ',b:3:1);
for i:=1 to max do
writeln('y_calc: ',y_calc[i]:3:1);
end.
Соседние файлы в папке Лаб1