|
Коэффициенты Ьнс вычисляются |
по |
формулам: |
|
|
|
||||||||||||||||
|
|
|
|
|
т |
|
|
|
|
|
|
гп |
|
|
|
|
т |
|
|
|
||
|
|
|
>п = |
Е |
|
|
|
|
|
|
1—1 |
|
|
* |
» |
= |
Е “ в! |
|
<4- |8> |
|||
|
|
|
|
|
;=1 |
т |
|
|
|
|
|
|
|
|
/=> |
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
т |
|
|
|
|
|||
|
|
|
&12 5= ^21 = |
5 ] |
а /1 а / 2» |
^13 = |
^31 = |
^ 1 |
а / 1а / 3» |
|
|
|||||||||||
|
|
|
|
|
|
/=1 |
|
|
|
|
|
т |
|
|
|
/=1 |
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
^23 = |
|
^32 = |
2 |
Л/2а/3* |
|
|
|
|
(4*19) |
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
/=1 |
|
|
|
|
|
|
|
|
|
|
Свободные |
плены |
|
с* |
|
вычисляются |
|
из |
соотношений: |
|
|
|||||||||||
|
|
|
|
т |
|
|
|
|
|
т |
|
|
|
|
|
т |
|
|
|
|||
|
|
|
С1 = |
2 |
аПу П |
|
с2 = |
|
Ц |
а/2у ; |
с3 = 5] а/зУг |
(4*20) |
||||||||||
|
|
|
|
/=1 |
|
|
|
|
|
у=1 |
|
|
|
|
1=1 |
|
|
|
||||
ч Ч„ |
|
|
|
|
|
% |
|
|
V |
/ |
|
|
|
|
|
|
|
|
|
|
||
|
г . ' |
|
|
|
Оц |
|
V |
Р а |
|
|
|
|
|
|
|
|
|
|
||||
|
ак |
\ X/ |
|
|
|
|
|
\ |
|
|
> |
Ч/ |
'“ О пО ыО зз'О '.РиОз^Оцйг&г'О пО гхО ц'О иОиРи'ОиОдОл |
|||||||||
> |
|
т% \й п ~ а21®!г |
У |
Он |
/ |
аи |
р а |
|||||||||||||||
|
ч |
|
|
V |
\ |
V |
|
|
|
|
|
|
|
|
|
|||||||
|
Оц |
4га |
|
|
|
/р31 /Р *2 |
/Р *3\ |
Ч/ |
*7» |
|
|
|
|
|
|
|
|
|||||
|
|
|
ч+ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Рис. 24. К вопросу вычисления |
определителей |
второго и |
третьего |
|||||||||||||||||||
порядков |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
Решение системы |
нормальных |
уравнении. |
При |
числе |
неизвест- |
||||||||||||||||
'ных п < 4 удобно их значения находить, решая систему |
нормаль |
|||||||||||||||||||||
ных уравнений с помощью определителей. |
Вычисление определите |
|||||||||||||||||||||
лей 2-го и 3-го порядка показано на рис. 24. |
|
|
2. |
|
||||||||||||||||||
|
Рассмотрим |
|
вначале |
решение |
системы |
(4.13) для п = |
|
|||||||||||||||
|
•Составляем и вычисляем главный определитель |
данной системы |
||||||||||||||||||||
уравнений: |
|
^ |
|
|
|
|
|
^11 |
^12 |
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
(4.21) |
||||
|
|
|
|
|
|
|
|
|
|
|
|
&21 |
Ь22 |
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
Составляем и вычисляем частные определители |
и ^ 2, |
заме |
|||||||||||||||||||
нив в системе (4.13) коэффициенты |
Ьы при соответствующих |
неиз |
||||||||||||||||||||
вестных |
на свободные |
|
члены |
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
0 2 = |
|
С1 |
|
Ь12 |
|
Д2 = |
|
^11 |
С1 |
|
|
|
(4.22) |
||||
|
|
|
|
|
|
^2 |
|
^22 |
|
|
Ь21 |
С2 |
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
Вычисляем |
наиболее |
вероятные |
значения |
неизвестных: |
|
||||||||||||||||
|
|
|
|
|
хх = |
Вг/0; |
х2 = |
^ 2/0* |
|
|
|
|
|
(4.23) |
||||||||
п = |
Аналогично |
решается |
|
и |
система |
нормальных |
уравнений |
для |
||||||||||||||
|
3. Вычисляем главный |
определитель |
системы: |
|
|
|
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
Ьп |
Ьп |
Ьп |
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
о |
|
= |
&21 |
^22 |
^23 |
• |
|
|
|
|
(4.24) |
||||
|
|
|
|
|
|
|
|
|
|
|
У31 |
У32 |
&33 |
|
|
|
|
|
|
|||
Вычисляем частные определители этой системы:
^1 ь» С2 А22
сз ^32
^13
^23 *ососо
|
А ц |
*1 *13 |
; |
Д 2 = ^21 |
А23 |
|
^31 С* А3З |
|
Ьп |
ь12 С1 |
|
Ь21 А,2 с2 |
(4.25> |
|
|
А32 % |
|
Наиболее вероятные значения неизвестных Г
х2 = |
х$ == 0^/5 . |
(4.26) |
Среднеквадратичное отклонение результатов совокупных (сов местных) ^измерений. Подставив вычисленные наиболее вероятные значения XI неизвестных в условные уравнения (4.04), можно найти
V^^ затем получить V* и сумму квадратов остаточных погрешностей
т
Среднеквадратичное отклонение результатов совокупных (со
вместных) измерений
Е
/= 1
5 - (4.27)
т — п Б *
Здесь /л — число условных уравнений; п — число неизвестных; АН1— адъюнкты элементов главной диагонали определителя О
(при А« / ) , полученные вычеркиванием строки А и столбца и, соот ветствующих данному элементу ЬН(> и последующим умножением на
Для |
п = |
2 |
адъюнкты |
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
А ц .— |
А22; |
Л22 — |
АА1. |
|
|
|
(4.28) |
||||
При числе |
неизвестных |
п == 3 |
адъюнкты |
Ли, |
Л 22 |
И Л33 |
вы- |
||||||||
числяются следующим |
образом: |
|
|
|
|
|
|
|
|
|
|||||
а |
_ |
А22 Ао3 |
в |
- |
_1Ап |
А131 |
|
л |
1Ах! А121. |
(4.29) |
|||||
11 ” “ /> |
Ъ |
1 |
Л22 — |
|
к |
и |
|
» |
л 33 =* |
||||||
|
|
"32 "33 |
|
|
|
1"31 |
"33 1 |
|
|
1А21 |
А22| |
|
|||
Доверительные границы случайной составляющей погрешности результатов совокупных (совместных) измерений. Задавшись дове рительной вероятностью у, из прил. 7 находим соответствующее значение коэффициента доверия В этом случае число степеней
свободы
к = т — п. |
(4.30) |
Тогда можно найти доверительные границы в^ случайной составляющей погрешности результата совокупных (совместных) измерений:
10. ПРОГРАММЫ ДЛЯ ОБРАБОТКИ СОВОКУПНЫХ И СОВМЕСТНЫХ ИЗМЕРЕНИЙ С ПОМОЩЬЮ микро ЭВМ «ЭЛЕКТРОНИКА БЗ-34»
Программа № 18 (обработка совокупных или совместных изме рений для любого т при п = 2). Программа обеспечивает вычисле ние наиболее вероятных значений неизвестных х? и доверительных границ гс случайных погрешностей этих результатов; число неиз вестных п = 2, количество условных уравнений т — любое. Не которые исходные данные вводятся в #регистры памяти до начала вычислений; значения параметров условных уравнений (а^, а-2
и //у) вводятся в машину при выполнении вычислений дважды* Программа № 18 приведена в табл. 57.
Инструкция для работы с |
программой: |
||||
1. |
Ввести и проверить программу. |
||||
2. |
Ввести |
исходные |
данные |
и |
очистить некоторые регистры |
памяти: |
т = |
ПО; 0 = |
П1 = П2 = |
ПЗ = П4 = П5 = П6; т — |
|
-л =■ ПД.
3. Нажав С/П, запустить вычисления по программе, на инди каторе / = 1. Ввести в регистры прмяти параметры первого услов ного уравнения
#и П7 |
П8 у± П9. |
|
|
Нажав С/П, получить на индикаторе / = 2. |
Ввести параметры |
||
второго условного уравнения |
|
|
|
#21 П7 #22 П8 У2-П9. |
|
|
|
Снова С/П и т. д. |
|
|
|
4. После ввода параметров последнего, т-го |
условного |
урав |
|
нения нажатием на С/П запустить вычисления, |
в результате |
чего |
|
на индикаторе высвечивается^значение хг. Нажатием на клавишу
ХУ вызвать |
|
на |
индикатор х2.. |
|
|
|
|
|
|
|
|||
5. |
Найти |
в |
прил. |
7 для к = т — п при выбранном у значение |
|||||||||
коэффициента |
доверия |
^ и ввести |
его в память, |
/ = |
П2. |
|
|||||||
Установить число циклов и очистить повторно один из регист |
|||||||||||||
ров памяти,— т = |
ПО; 0 = П6. |
|
|
|
|
|
должен |
||||||
6. |
Снова С/П (перед этим на индикаторе обязательно |
||||||||||||
быть нуль!), на индикаторе высвечивается |
/ = |
1. |
Повторно |
вве |
|||||||||
сти параметры |
условных уравнений,— |
|
|
|
|
|
|
||||||
|
|
|
|
|
Оц П7 012 П8 |
У1 П9 С/П |
|
|
|
|
|
||
И т. д. |
После |
ввода |
параметров последнего условного |
уравнения |
|||||||||
7. |
|||||||||||||
и С/П — на |
индикаторе появляется |
значение |
Нажатием |
на |
XV |
||||||||
вызвать на |
индикатор |
значение е2. |
|
|
|
|
памяти: |
||||||
8. |
Все искомые величины записаны также и в регистрах |
||||||||||||
|
|
ИПА = |
хг\ |
ИПВ = х2; |
ИПС = |
ех; |
ИПД = |
е2. |
|
|
|||
9. |
Для |
обработки |
новой совокупности |
данных |
повторить |
опе |
|||||||
рации по пп. 2—8. |
|
регистров |
адресуемой |
памяти |
показано |
||||||||
10. |
Использование |
||||||||||||
в табл. |
58. |
|
|
|
|
|
|
|
|
|
|
|
|
Адрес
|
|
|
о |
|
|
|
О |
|
|
Коман |
|
Инди |
о |
К оман- |
|
Инди |
и |
|
Инди |
Код |
ЕС |
Код |
а ■ Коман |
К од |
|||||
да |
кация |
Да |
кация |
да |
кация |
||||
|
|
< |
|
|
|
00 |
КИП6 |
Гб |
33 |
ИП2 |
62 |
|
01 |
ИП6 |
66 |
34 |
Рх2 |
22 |
|
02 |
с / п |
50 |
/ |
|
11 |
|
|
|
|
35 |
— |
||
03 |
X |
12 |
36 |
ПС |
4С |
|
-04 |
ИП5 |
65 |
37 |
ИП4 |
64 |
|
|
|
|
38 |
ИПЗ |
-63 |
|
05 |
+ |
10 |
39 |
X |
12 |
|
06 |
П5 |
45 |
40 |
ИП5 |
65 |
|
07 |
XV |
14 |
||||
08 |
ИП9 |
69 |
41 |
ИП2 |
62 |
|
09 |
X |
12 |
42 |
X |
12 |
|
|
|
|
43 |
— |
11 |
|
10 |
ИП4 |
64 |
44 |
ИП1 |
61 |
|
11 |
+ |
10 |
45 |
ИП5 |
65 |
|
12 |
П4 |
44 |
||||
13 |
ИП8 |
68 |
46 |
X |
12 |
|
14 |
Р х 2 |
22- |
47 |
ИП2 |
62 |
|
|
|
|
48 |
ИП4 |
64 |
|
15 |
ИПЗ |
63 |
49 |
X |
12 |
|
16 |
+ |
10 |
50 |
|
11 |
|
17 |
п з |
43 |
____ |
|||
18 |
ИП8 |
68 |
51 |
ипс |
6С |
|
19 |
ИП7 |
67 |
52 |
: |
13 |
|
|
|
|
53 |
пв |
41 |
|
20 |
X |
12 |
54 |
XV |
14 |
|
21 |
ИП2 |
62 |
55 |
и пс |
|
|
22 |
+ |
10 |
6С |
|||
23 |
П2 |
42 |
56 |
; |
13 |
|
24 |
ИП7 |
67 |
57 |
ПА |
4- |
|
25 |
Рха |
22 |
58 |
С/П |
50 |
|
|
|
|
||||
26 |
ИП1 |
61 |
59 |
П4 |
44 |
|
27 |
+ |
10 |
||||
|
|
|
||||
28 |
П1 |
41 |
60 |
П5 |
45 |
|
29 |
РЬО |
5Г |
||||
61 |
КИП6 |
Гб |
||||
30 |
00 |
00 |
62 |
ИП6 |
66 |
|
63 |
С/П |
50 |
||||
31 |
ИПЗ |
63 |
||||
32 |
X |
12 |
64 |
ИП7 |
67 |
|
|
|
|
*1
/
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
ИПА |
6- |
X |
12 |
ИП8 |
68 |
ИПВ |
6Ь |
X |
12 |
-4~ |
10 |
ИП9 |
69 |
— |
11 |
Рха |
22 |
ИП4 |
64 |
+ |
10 |
П4 |
44 |
РЬО |
5Г |
61 |
61 |
ИПД |
6Г |
: |
13 |
ИПС |
6С |
• |
13 |
П5 |
45 |
ИП1 |
61 |
X |
12 |
Р * |
21 |
ИП2 |
62 |
X12
пд 4Г
ИП565
ИПЗ 63
X12
Р/ 21
ИП2 62
X |
12 |
|
ПС |
4С |
г 1 |
С/П |
50 |
Пример 23 (вычисления по программе М 18). |
При исследова |
|||||||||
нии зависимости сопротивления медной |
пройолоки |
от |
температуры |
|||||||
получено ‘семь |
пар |
значений |
/? = Р (/), |
внесенных в табл. 59. |
||||||
Полагая, |
что |
зависимость между |
/? |
и 1 описывается |
формулой |
|||||
(1.06), |
найти значения /?0 и а |
в этой формуле, |
а также |
оценить точ |
||||||
ность |
полученных |
результатов. |
получаем |
выражение |
п/-1л:1 + |
|||||
Преобразовав, |
формулу |
(1.06), |
||||||||
•+- а/2*г = У{, где |
= 1; хг ^ Л„; «/2 = |
*’> *г = |
а#о'> У/ |
|
|
|||||
0 т
1 0
20
30
40
50
6 0
Вводим
Организация циклов |
7 |
ап |
|
|
|
|
1 |
|
|
|
|
|
|
|
Е «},; Ьп |
|
8 |
°/2 |
|
Е ац а& ь п> |
9 |
У/ |
|
|
|
|
|||
^ Я/2> ^22 |
|
А |
|
Ч |
Е апУп |
Е °? |
В |
|
Ч |
Е а12У,’ ч , -{~ - п) ь |
С |
|
Ч |
|
Косвенная |
индикация |
д |
т— п |
бо |
исходные данные,—
т — 7 = ПО; т — п — 7 — 2 = 5 = ПД.
Очищаем регистры памяти, — 0 = П1 = .. = П6. Запускаем вычисления (С/П). Вводим данные для / = 1.
1 П7 10,0 П8 55,28 П9 С/П. |
|
Таблица 59 |
||||||||
Время |
вычисления — около |
10 с. |
|
1, *С |
|
|||||
Вводим |
данные |
для |
/ = |
2: |
|
/ |
Я. Ом |
|||
1 П7 15,1 П8 |
56,42 П9 С/П и т. д. |
|
|
|
||||||
После |
ввода |
параметров |
|
последнего |
1 |
10,0 |
55,28 |
|||
условного уравнения время вычисления око |
2 |
15,1 |
56,42 |
|||||||
ло 18 с; на индикаторе |
х\ — 53,014 577. |
3 |
20,0 |
57,4Ь |
||||||
Находим также Тг — 0,225 373 49. |
|
|
|
|||||||
Снова |
вводим |
некоторые |
исходные |
4 |
25,0 |
58,65 |
||||
данные: |
|
|
|
|
|
|
||||
т = |
7 = |
ПО; |
(у = |
2,571 |
= |
5 |
30,0 |
59,80 |
||
6 |
35,1 |
60,92 |
||||||||
= |
П2 (для т — п = |
5); |
7 |
40,1 |
62,05 |
|||||
Очищаем регистр 6: |
0 = |
П6.- |
||||||||
|
|
|
||||||||
Нажимаем С/П — на индикаторе/ = 1; вводим повторно параметры условных урав нений
1П7 10,0 П8 55,28 П9 С/П
ит. д. Время вычисления — 7 с.