|
|
Операции |
|
Наименование кла- |
||||
|
|
|
|
|
|
|
виш |
|
|
|
|
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 СОСТАВЛЕНИЕ УРАВНЕНИЯ РЕГРЕССИИ ВТОРОГО ПОРЯДКА