Главная -> Книги

(0) (1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) (12) (13) (14) (15) (16) (17) (18) (19) (20) (21) (22) (23) (24) (25) (26) (27) (28) (29) (30) (31) (32) (33) (34) (35) (36) (37) (38) (39) (40) (41) (42) (43) (44) (45) (46) ( 47 ) (48) (49) (50) (51) (52) (53) (54) (55) (56) (57) (58) (59) (60) (61) (62) (63) (64) (65) (66) (67) (68) (69) (70) (71) (72) (73) (47)

(до, а,, аг.....а„ заносятся также в регистры РО, Р1, Р2.....Рп).

Для определения y(z) по формуле (5.74) для каждого г вводим п = РД, г = X - дсо = РХ и нажимаем клавиши БП 78 С/П.

ПА Сх ПД КИПД ПВ ИПД 1 + ПС ИПД

1 + ПД ИПВ КИПД ПВ XY - ИПС -г-

КПД ИПД ИПА - Jx=0 09 ИПС ПД ИПА -

Fx=0 03 ИПА 1 ис ПВ ИПВ 1 +

ПД кипв КИПД ИПС X - КПВ ипв 1 ч-

ПВ ИПА - Рлг=0 37 ИПС 1 - Рд:=0 35

ПД С/П пв КИПД ИПС н- КПД С/П ИПС ИПВ

X ПС ИПД I + ПД БП 63 t 1

f Сх F, X КИПД + ИПД 1 -

ПД Рх<0 83 F, С/П БП 78

Для указанного выше примера получим «о = 0,4860966, ki == е= 0,0035472167, «а «= -0,0003899, «з = 2,2583333-10-», а4 =. = -6•I0- Оо = 0,4860966, aj = 3,5472167.10- аг = -0,3899• 10-», as «= 0,022583333, 04 - -б-10-». Для х = 2,45 (г == 0,25) лолучаем г/(х)= 0,49285719,

5.14. Сглаживание эмпирических данных с ошибками

Нередко функция yi(xi) задается значениями yi прн равиоот-стоящих абсциссах xi. Значения yi могут содержать случайные отклонения. По ряду yt, используя усреднение на основе интерполяционных формул, можно уточнить значения ординат, т.е. провести сглаживание зависимости г/((х(). Сглаженные ординаты записываются в виде yi.

Программа 5.99. Линейное сглаживание по трем ординатам о помощью формул:

Уй = (5уо -1- 21 - У2)/6, г = О,

Ввод: уо = P7, г/1 = P8, г/а = P9 н далее (после нажатия клавиш В/О н С/П) Уу у, у в регистр X. Вывод р, у, у, --.,у в регистр X.

ИП7 5 X ИП8 2 X + ИП9 - 6 - С/П ИП7 ИПЗ ИП9 +4-3 -f- С/П Пб ИПЗ П7 ИПЭ П8 ИПб П9 БП 12 ИПЭ 5 X ИП8 2 X + ИП7 - 6 -f. С/П

После ввода г/о, г/i и y% нажимаем клавиши В/О и С/П и получаем Уй. Еще раз нажав клавишу С/П, получаем yi. Теперь вводим



уг И, нажав клавишу С/П, получаем уз а т. д. После ввода послед, него значения у и получения г/др , нажимаем клавиши БП 29 « С/П и получаем g.f. Результаты сглаживания для функции у т Юх при АГ = 0,1 -г- 1 даны в табл. 5.11.

Таблица 5.11

Сглаживание

Сглаживание

Готаое значение

по трем ординатам

по пяти ординатам

y=iOx

0,87

0,992

2.12

1,98

1.995

2,92

3,0633333

2,998

4,15

3.99

4,033

5,05

4.99S

£,9733333

6,044

6,92

7,056и663

7,024

8,15

8,04

8,004

9,05

8.957

9,323

9.91

Программа 5.100, Линейное сглаживание при пяти ординатах по формулам;

Уа ~ (Зг/о + 2у, + у2 yi)lo, i = О,

i/i = (4г/о + Зу1+2у2 + уз)/10, i = 1,

h = (у1--г + /-1 + «г + Ум + Ji+2)/5> 2<j < iV - 2,

5/., = [Ум-ъ + -N-2 + Ум-1 + 4у)/10. * -/V

h = (Зг/л/ + 2!/,v-i + Уы-2 - Уы-ъ)1 =

Ввод: = Р5, yi = Р6, уа == Р7, уа = Р8, у4 = Р9, у,, в регистр X (см. пояснения в контрольном примере).

N.

Проведем сглаживание значений у/, приведенных в табл. 5.!! Введя ya-i-yi и нажав клавиши В/О и трижды С/П (с перерывами на время счета), получим уа, Уи Уг- Далее вводим у, в получаем й



к т. д. Аналогично получаем у4, У, 4дг 2- После ввода у, нажимаем клавиши БП 56 С/П и получаем y, , затем еше рач нажимаем клавишу С/П и получаем у (числовые данные см. в табл. 6.11).

Для функций yi(Xi), заметно отличающихся от линейных, используется нелинейное сглаживание с помощью полиномов высоких степеней т. При m = 3 необходимое число ординат составляет 7.

Программа 5.101. Нелинейное сглаживание при семи орднната.ч по формуле

i = (Vi-3+ Vi-2+ "эг-! + "Ai/i + «вгг+. + чсг+а+ЗД+з)

где индексы у коэффициентов а соответствуют номерам регистров, нз которых берутся их численные значения. Ввод; г/о-г г/в = РО-Ь Р6

Таблица 5.12

Р8

Пуск программы

В/О С/П

В/О С/П

В/О С/П

В/О С/П

( • •

В/О С/П

В/О С/П

В/О С/П

Таблица 5.13

(Сглаженные)

у точные

0.70142857

0,69314717

1,08

1,0807142

1,0986122

1.39

1.385

1,3862944

1,64

1,6242857

1,6094379

>,7б

1,7995238

1,7917594

1,99

1,9509523

1,9459101

2,04

2,0733333

2,0794415

2,22

2,1914283

2,1972245

2,28

2,3028571

2,3025851

2,42

2,412619

2,3973952

И далее в соответствии с табл, 5.12 н указаниями контрольного при-мера.



(0) (1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) (12) (13) (14) (15) (16) (17) (18) (19) (20) (21) (22) (23) (24) (25) (26) (27) (28) (29) (30) (31) (32) (33) (34) (35) (36) (37) (38) (39) (40) (41) (42) (43) (44) (45) (46) ( 47 ) (48) (49) (50) (51) (52) (53) (54) (55) (56) (57) (58) (59) (60) (61) (62) (63) (64) (65) (66) (67) (68) (69) (70) (71) (72) (73)