Материал: 4572

Внимание! Если размещение файла нарушает Ваши авторские права, то обязательно сообщите нам

 

 

Операции

 

Наименование кла-

 

 

 

 

 

 

 

виш

 

 

 

 

1

 

 

2

 

 

 

 

 

 

 

 

 

 

1.

Включить микрокалькулятор

 

 

 

 

ВКЛ

 

 

 

 

 

 

 

 

 

2.

Перейти в режим программирования

 

 

F ПРГ

 

 

 

 

 

 

 

 

3.

Ввести программу

 

 

 

(см. табл. 4.2)

 

 

 

 

 

 

 

 

4.

Очистить программный счетчик

 

 

 

F ABT В/О

 

 

 

 

 

 

 

 

 

 

5.

Ввести значения: y1 в регистр Р7

 

 

 

y1

Х-П

7

 

 

y2 в регистр Р8

 

 

 

y2

Х-П

8

 

 

y3 в регистр Р9

 

 

 

y3

Х-П

9

6.

Ввести количество циклов, равное - 2 в регистр памяти ПО

n- 2

Х-П

0

 

 

 

 

 

 

 

 

 

7.

Определить численное значение

 

 

 

 

С/П

 

8.

Определить численное значение У2

 

 

С/П

 

9.

Ввести yi

 

 

 

 

 

 

 

 

 

 

 

 

 

 

10. Определить значение

 

 

 

 

С/П

 

11. Ввести следующее значение уi

 

 

 

 

 

 

12. Определить значение

 

 

 

 

С/П

 

13. Провести вычисления остальных значений yi включая

(см. пункты

 

 

 

11 и 12)

 

 

 

 

 

 

 

14. Определить значение

 

 

 

 

С/П

 

 

 

 

 

 

 

 

Таблица 4.4

Номера

Экспериментальные значе-

 

Выровненные значе-

Время выполнения расчетов по

 

узлов

ния

 

ния

циклам

 

 

 

 

1

y0=0,9

 

=0,97

 

4 с.

 

 

 

2

у1=2,12

 

=1,98

 

3 с.

 

 

 

3

у2=2,92

 

=3,0633333

 

7 с.

 

 

 

4

y3 =4,15

 

=4,0783333

 

4 с.

 

 

n= 4

 

 

 

 

 

 

 

Контрольные вопросы

1.Дайте определение промаха измеряемой величины и поясните порядок ее выбраковки.

2.Какими способами осуществляется сглаживание экспериментальных данных. 3.Порядок сглаживания экспериментальных данных по трем ординатам на программируемом микрокалькуляторе.

5 ОПРЕДЕЛЕНИЕ КОЭФФИЦИЕНТОВ УРАВНЕНИЯ РЕГРЕССИИ ПЕРВОГО ПОРЯДКА

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

5.1 Применение метода наименьших квадратов для составления уравнения регрессии

В результате выполнения экспериментальных исследований интересующая нас функция обычно оказывается заданной в табличном виде, когда значениям аргумента (варьируемым факторам)x1; х2…хn соответствуют конкретные значения функции отклика y1 ;y2 ...уn . Для получения аналитического выражения, являющегося функцией отклика, проводят регрессионный анализ, позволяющий получить уравнение регрессии в виде приближенной эмпирической формулы, с помощью которой можно с заданной точностью оценить степень и характер влияния каждого из факторов на выходную величину.

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

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

у = f (х).

Далее находят основные параметры аналитической функции, которая с минимальной погрешностью приближалась бы к исходной функции

у = f (x). Таковому условию в наилучшей степени отвечает требование о том, чтобы сумма квадратов отклонений ∑∆ откликов (y1;y2 ;...уn) от линий функции у = f (x) была бы минимальной. Указанный метод и получил наименование метода наименьших квадратов. С его помощью могут быть определены коэффициенты уравнений регрессии линейного, параболического, степенного, гиперболического, экспоненциального и т.п. видов.

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

На первых этапах исследований чаще всего обращаются к регрессионным моде-

лям первого порядка, линейным уравнениям вида

 

y = b0+b1x,

(5.1)

где b0 и b1 - коэффициенты линейного приближения, определяемые по формулам:

)

(5.2)

 

(5.3)

Применение уравнений 5.1-5.3 становится наиболее целесообразным, когда и ходе эксперимента наблюдается при возрастании факторов близкое к пропорциональному возрастание или убывание откликов.

Ниже приводится программа НИП - 6 (табл.5.2) и инструкции по ее пользованию(табл. 5.3),предназначенная для определения коэффициентов линейного уравнения регрессии. В разделе 6 приведена программа для определения коэффициентов уравнения регрессии параболического вида. Для отыскания другого типа можно воспользоваться программами, приводимыми в специальной литературе (4).

5.2 Задание. Вычислить коэффициенты и найти линейное уравнение регрессии

Порядок выполнения 1.Получить у преподавателя таблицу парных значений изучаемых эксперимен-

тальных величин хi и уi.

Примечание: При выполнении задания с элементами УИРС или при выполнении научных исследований студенты используют значения хi и уi полученные в результате проведенных опытов.

2.Составить блок-схему алгоритма программы вычислений.

3.Проанализировать таблицу 5.1 распределения регистров памяти предстоящих расчетов по программе.

4.Проанализировать программу (табл.5.2) и инструкцию (5.3) по работе с программой.

5.Ввести программу в микрокалькулятор и, пользуясь инструкцией (табл.5.2), поочередно ввести в него:

-количество парных данных n;

-значения хi; -значения уi.

Примечание. Программа составлена с использованием диалоговой формы введения данных. Чтобы не допускать ошибок, каждый раз перед вводом в программу значений хi или yi в индикаторе высвечивает их порядковый номер (1, 2,3...n).

6.После завершения счета по программе в соответствующих регистрах микрокалькулятора отыскать значения b0а) и bi(Pb) и занести их в соответствующую строку табл.5.1.

7.Составить уравнение регрессии

Пример 8. Найти коэффициенты b0 и b1| и составить уравнение регрессии для экспериментальных данных, деформации (мм) древесных образцов в зависимости от давления прижима (мПа). Результаты эксперимента (1) сведены в таблицу 5.4.

Таблица 5.1

Распределение регистровой памяти микрокалькулятора

Наименование информации

Номера адресных регистров

 

ресных регист

 

 

 

ров

 

 

 

 

 

 

 

 

информации

РО

Р1

 

Р2

РЗ

Р4

Р5

Р6

Р7

Р8

Ра

Рb

Исходные данные

 

 

 

 

 

 

 

 

 

 

 

 

Результаты вычислений

000001

n

 

xi

yi

∑xi

Xi2

∑yi

∑xiyi

(∑xi)2

b1

b2

Команда вызова

 

ХП1

 

 

 

 

 

 

 

 

ХПа

ХПв

Таблица 5.2 НИП - 6. Программа вычисления коэффициентов b0 и b1, методом наименьших

квадратов

Адрес

Команда

Код

Комментарий

Адрес

Команда

Код

Комментарий

00

ХП1

41

n→ в Р1

32

+

10

∑ хi уi → в Р7

01

ХПО

40

т→ в Р0

33

ХП7

47

 

 

02

0

00

 

34

Flo

5F

Повторение цикла

03

ХП4

44

 

35

7

07

для ввода очеред-

04

ХП5

45

Обнуление

36

ПХ4

64

ных хi

и yi

05

ХП6

46

Р4, Р5, Р6, Р7

37

Fx2

22

 

 

06

ХР7

47

 

38

ПХ1

61

 

 

07

ПХ1

61

 

39

ПХ5

65

 

 

08

ПХО

60

 

40

x

12

 

 

09

ПХ0

11

Счет в режиме

41

-

11

( )2→n

→в

10

1

01

диалога

42

ХП8

48

Р8

 

11

+

10

 

43

ПХ4

64

 

 

12

С/П

50

Ввод хi → в Р2

44

ПХ6

66

 

 

13

ХП2

42

 

45

X

12

 

 

14

С/П

50

Ввод yi → в Р2

46

ПХ1

61

 

 

15

ХПЗ

43

Расчеты по ф. 6.2 и

47

ПХ7

67

 

 

16

ПХ6

66

6.3

48

x

12

 

 

17

+

10

 

49

-

11

 

 

18

ХП6

46

∑уi→ вР6

50

ПХ8

68

 

 

19

ПХ2

62

 

51

÷

13

 

 

20

ПХ4

64

 

52

Хпа

4-

b1 →в Ра

21

+

10

 

53

ПХ6

66

 

 

22

ХП4

44

∑ хi → в Р4

54

Пха

6-

 

 

23

ПХ2

62

 

55

ПХ4

64

 

 

24

Fx2

22

 

56

x

12

 

 

25

ПХ5

65

 

57

-

11

 

 

26

+

10

 

58

ПХ1

61

Индикация bо

27

ХП5

45

( )2→ в Р5

59

÷

13

 

 

28

ПХ2

62

 

60

ХПb

4L

bo→ в РЬ

29

ПХ3

63

 

61

С/П

50

Конец счета

30

x

12

 

 

 

 

 

 

31

ПХ7

67

 

 

 

 

 

 

После обработки данных по программе НИЛ - 6 (табл. 5.2) получены следующие результаты:

b0 =0,01439186, b1=0,12256756

Округлив значения b0 и b1, составим уравнение регрессии y = -0,0144 + 0,123х

Таблица 5.3

Инструкция работы с программой НИЛ – 6

 

 

Операции

 

Наименование клавиш

 

 

1.

Включите микрокалькулятор

 

ВКЛ

 

 

2. Перейдите в режим «Пpограммирование»

 

F ПРГ

 

 

3.

Введите программу НИП - 5

 

(Табл. 5.2)

 

 

4.

Очистите программный счетчик

 

F АВТ В/О

 

 

5.

Наберите значение n

 

 

 

 

 

6.

Введите п в программу

 

 

С/П

 

 

7.

Наберите значение x1

 

 

 

 

 

8.

Введите х1 в программу

 

С/П

 

 

9.

Наберите значение y1

 

 

 

 

 

10. Введите y1 в программу

 

С/П

 

 

11. Повторите вычисления ввода последовательно:

 

 

 

 

 

- х2 3 ...xi

 

 

 

 

 

 

-y2;y3…yi

(см. п. п. 7-10)

 

 

 

 

12. По окончании счета - в индикаторе высвечивает

 

 

 

 

значение коэффициента b0.

 

 

 

 

Это же значение b0 → в РЬ, значение коэффициента

 

 

 

 

b1→ в Ра

 

 

ПХb

 

 

 

 

 

 

ПХа

 

 

 

 

Таблица 5.4

 

 

 

 

 

Результаты испытаний образцов древесины(4)

 

 

 

 

 

 

 

 

 

 

Номер опыта

Р (МПа)

 

мм

 

 

 

1

1,5

 

0,16

 

 

 

2

1,0

 

0,11

 

 

 

3

2,0

 

0,215

 

 

 

4

2,5

 

0,33

 

 

 

5

3,5

 

0,4

 

 

 

n=5

-

 

-

 

Контрольные вопросы

1.В чем заключается сущность метода наименьших квадратов?

2.Дайте анализ регрессионной модели первого порядка.

3.По каким формулам рассчитываются коэффициенты линейного уравнения регрессии?

4.Порядок расчета коэффициентов уравнения регрессии с использованием программируемого калькулятора.

6 СОСТАВЛЕНИЕ УРАВНЕНИЯ РЕГРЕССИИ ВТОРОГО ПОРЯДКА