Аппроксимация экспериментальных данных функциональными зависимостями

Страницы работы

Содержание работы

МИНИСТЕРСТВО ОБРАЗОВАНИЯ РЕСПУБЛИКИ БЕЛАРУСЬ

Белорусский  национальный технический университет

Кафедра «Электрические станции»

Группа 106129

Лабораторная работа №22

АППРОКСИМАЦИЯ ЭКСПЕРИМЕНТАЛЬНЫХ ДАННЫХ ФУНКЦИОНАЛЬНЫМИ ЗАВИСИМОСТЯМИ

Исполнитель:  Чижевский А.М.

Руководитель: Климкович П.И.

Минск, 2010

Цель работы: Получение практических навыков составления программ, реализующих численные методы аппроксимации экспериментальных данных функциональными зависимостями. Сравнительный анализ рассмотренных методов.

Задание:

По заданным экспериментальным данным предыдущей лабораторной работы определить вид нелинейной аппроксимирующей зависимости, дающее наилучшее приближение к данному набору точек. В соответствии с выбранной зависимостью построить схему алгоритма и составить программу. Выполнить программу на ЭВМ. Нанести заданные точки на координатную плоскость. Построить в тех же координатах рассчитанную зависимость.   

№

варианта

Значение

Аргумента

Значение

Функции

Пределизменения

аргумента

Шаг

интегрирования

28

0,88

0,98

1,22

1,63

2,22

2,68

2,93

3,02

0,76

1,649

2,41

2,986

3,135

3,06

2,654

2,52

[0,7;3,1]

0,1

Таблица идентификаторов:

N

Xi

A

YI

XN

XK

X

B

C

D

B1

B0

N

X(I)

A

Y(I)

XN

XK

DX

B

C

D

B1

B0

Блок –схема:

началоA=,B=,C=,D=ввод  NВвод x(i),y(i),xkxn,xk,dxDo 2 i=1,Ny1=B1=,B0=Вывод B1,B0конецIf 100< N<1If x(i)=0A=,B=,C=,D=Do  x1=xk,xn,dxВывод x,y

Текст программы:                                                                                                   

D=0

DO 2 I=1,N

IF(X(I).EQ.0) GO TO 2

A=A+alog(X(I))

B=B+alog(Y(I))

C=C+alog(X(I))**2

D=D+alog(Y(I))*alog(X(I))

2   CONTINUE

B1=(A*B-N*D)/(A*A-N*C)

B0=exp(1/N*(B-B1*A))

WRITE (6,*)'B1=',B1

WRITE (6,*)'B0=',B0

DO 1 X1=XN,XK,DX

Y1=B0*X1**B1

WRITE(6,*)'PRI X=',X1,'Y=',Y1

1   CONTINUE

STOP

END

Результаты:


B1=  0.70969260

 B0=   1.0000000


 PRI X=  0.69999999     Y=  0.77636725

 PRI X=  0.80000001     Y=  0.85353965

 PRI X=  0.90000004     Y=  0.92795360

 PRI X=   1.0000000     Y=   1.0000000

 PRI X=   1.1000000     Y=   1.0699811

 PRI X=   1.2000000     Y=   1.1381365

 PRI X=   1.3000001     Y=   1.2046608

 PRI X=   1.4000001     Y=   1.2697142

 PRI X=   1.5000001     Y=   1.3334315

 PRI X=   1.6000001     Y=   1.3959262

 PRI X=   1.7000002     Y=   1.4572965

 PRI X=   1.8000002     Y=   1.5176270

 PRI X=   1.9000002     Y=   1.5769917

 PRI X=   2.0000002     Y=   1.6354557

 PRI X=   2.1000001     Y=   1.6930768

 PRI X=   2.2000000     Y=   1.7499065

 PRI X=   2.3000000     Y=   1.8059909

 PRI X=   2.3999999     Y=   1.8613716

 PRI X=   2.4999998     Y=   1.9160862

 PRI X=   2.5999997     Y=   1.9701689

 PRI X=   2.6999996     Y=   2.0236511

 PRI X=   2.7999995     Y=   2.0765610

 PRI X=   2.8999994     Y=   2.1289251

 PRI X=   2.9999993     Y=   2.1807673


График

Похожие материалы

Информация о работе

Предмет:
Информатика
Тип:
Отчеты по лабораторным работам
Размер файла:
34 Kb
Скачали:
0