3. Моделирующая программа.
Моделирующая программа составлена на языке программирования Turbo Pascal 7.0 [2]. Текст программы приведен в приложении 1.
3.1. Последовательность вычислений.
3.1.1. Ввод исходных данных для расчета:
Сист – истинное значения измеряемой константы;
Emax – максимальная величина ошибки (по модулю);
N – количество замеров;
Окончание ввода при нажатии клавиши < Enter >.
3.1.2. Получение массива Х зашумленных замеров по формуле ;
Xi = Сист + 2* Еmax*(Random (1)-0.5)
3.1.3. Расчет оценки Сср замеров константы по формуле (2) и среднеквадратической ошибки Ссрerr по формуле (3).
3.1.4. Расчет границ FH (i), FВ (i) – интервала неопределенности Hi для минимаксного метода по формуле (4);
3.1.5. Сортировка двух массивов FH (i) и FВ (i) для определения максимального FH (i) и минимального FВ (i) (сортировка методом «пузырька»), которым присваиваются значения IH и IB.
3.1.6. Проверка на равенство IH и IB, если они равны, то оценке Cminmax присваивается значение Cminmax = IH = IB.
3.1.7. Определение вероятности попадания оценок Сср и Cminmax в заданный диапазон +-5% от Сист, т.е. определение интервалов наложения Нср и Hminmax и отношений (Нср/ Сср) и (Нminmax/ Сminmax) которые и будут являться вероятностью.
3.1.11. Вывод всех рассчитанных значений в файл result.txt;
3.1.12. Выход из программы.
3.2. Блок-схема программы
Начало
0
1
1
0
1
0
0
1
4. Результаты моделирования.
Контрольный расчет произведен для следующих исходных данных:
Сист =50
Emax =15
N = 20
I |
Emax |
Сист |
Xi |
Fн(I) |
Fв(I) |
Cср |
Сerr |
1 |
15 |
50 |
63,83 |
48,83 |
78,83 |
49,04 |
1,88 |
2 |
15 |
50 |
56,85 |
41,85 |
71,85 |
|
|
3 |
15 |
50 |
42,08 |
27,08 |
57,08 |
|
|
4 |
15 |
50 |
47,20 |
32,20 |
62,20 |
|
|
5 |
15 |
50 |
37,67 |
22,67 |
52,67 |
|
|
6 |
15 |
50 |
56,12 |
41,12 |
71,12 |
|
|
7 |
15 |
50 |
47,67 |
32,67 |
62,67 |
|
|
8 |
15 |
50 |
54,68 |
39,68 |
69,68 |
|
|
9 |
15 |
50 |
52,50 |
37,50 |
67,50 |
|
|
10 |
15 |
50 |
44,36 |
29,36 |
59,36 |
|
|
11 |
15 |
50 |
55,60 |
40,60 |
70,60 |
|
|
12 |
15 |
50 |
38,05 |
23,05 |
53,05 |
|
|
13 |
15 |
50 |
51,84 |
36,84 |
66,84 |
|
|
14 |
15 |
50 |
55,24 |
40,24 |
70,24 |
|
|
15 |
15 |
50 |
63,02 |
48,02 |
78,02 |
|
|
16 |
15 |
50 |
36,84 |
21,84 |
51,84 |
|
|
17 |
15 |
50 |
47,47 |
32,47 |
62,47 |
|
|
18 |
15 |
50 |
39,05 |
24,05 |
54,05 |
|
|
19 |
15 |
50 |
51,02 |
36,02 |
66,02 |
|
|
20 |
15 |
50 |
39,70 |
24,70 |
54,70 |
|
|
Iн |
Iв |
Сminmax |
CerrH |
CerrB |
+5%Сист |
-5%Сист |
Pcp |
Pminmax |
48,83 |
51,84 |
50.34 |
47,16 |
50,92 |
47,5 |
52,5 |
91,2% |
100% |
Приложение 1