Скачиваний:
10
Добавлен:
01.05.2014
Размер:
1.52 Кб
Скачать
uses sampler;

const
szUnit: String = 'MA_NAR.pas';
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,k:integer;
a,b:real;
begin
{creating new mass x and y}
for k:=1 to 80 do
begin
Sample ( szUnit,1 ) ;
randomize;
Sample ( szUnit,2 ) ;
writeln('Input Data: ');
Sample ( szUnit,3 ) ;
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;
Sample ( szUnit,4 ) ;
linfit1(x,y,y_calc,a,b,max);
Sample ( szUnit,5 ) ;
writeln('a: ',a:3:1,' b: ',b:3:1);
Sample ( szUnit,6 ) ;
for i:=1 to max do
writeln('y_calc: ',y_calc[i]:3:1);
Sample ( szUnit,7 ) ;
end;
end.
Соседние файлы в папке PAS