МИНИСТЕРСТВО ОБРАЗОВАНИЯ РЕСПУБЛИКИ БЕЛАРУСЬ
Белорусский национальный технический университет
Кафедра «Электрические станции»
Группа 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 |
Блок –схема: |
|||||
|
|||||
Текст программы:
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
График
Уважаемый посетитель!
Чтобы распечатать файл, скачайте его (в формате Word).
Ссылка на скачивание - внизу страницы.