Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
Informatike.docx
Скачиваний:
1
Добавлен:
01.05.2025
Размер:
84.29 Кб
Скачать

34. Численные интегрирования.

35. Апроксимация данных.

Program lVl;

uses crt;

var a, b, c, d, a0, a1, a10, a11, p: real; i, n: integer;

x, y, x1, y1, f:array[1..12] of real;

begin

clrscr;

write('n='); readln (n);

write ('x(i) и y(i):');

for i:=1 to n do

begin

write ('x(',i:2,')='); readln (x[i]);

write('y(',i:2,')='); readln (y[i]);

x1[i]:=

y1[i]:=

end;

writeln ('числовые данные для аппроксимации:');

for i:=1 to n do begin

write('x(',i:2,')=',x[i]:2:4);

writeln ('y(',i:2,')=',y[i]:2:4);

end;

a:=0; b:=0; c:=0; d:=0;

for i:=1 to n do begin

a:=a+x1[i]; b:=b+y1[i]; c:= c+sqr(x1[i]); d:=d+x1[i]*y1[i];

end;

a10:= (b*c-a*d)/(n*c-sqr(a));

a11:=(n*d-a*b) /(n*c-sqr(a));

a0:=a10;

a1:=a11;

writeln ('параметры аппроксимирующей зависимости:');

writeln ('a0=', a0:2:4);

writeln ('a1=', a1:2:4);

writeln ('y[i]:');

for i:=1 to n do begin

f[i]:= 1/a0+a1*sqrt(x[i]));

write ('x(',i:2,')=', x[i]:2:4);

write ('f (',i:2,')=',f[i]:2:4);

end;

p:=0;

for i:=1 to n do p:=p+abs(y[i]-f[i])/y[i];

p:= p/n*100;

writeln ('Погрешность=',p:2:2,'%');

writeln(1st,’погрешность=’,p:2:2,’%’);

readln

Результат Y[i]:

x (1)= f (1)=

x (2)= f (2)=

x (3)= f (3)=

x (4)= f (4)=

x (5)= f (5)=

x (6)= f (6)=

x (7)= f (7)=

x (8)= f (8)=

x (9)= f (9)=

x (10)= f(10)=

x (11)= f (11)=

x (12)= f (12)=

Погрешность равна:

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]