ПРИМЕНЕНИЕ МЕТОДА ЦЕНТРА НЕОПРЕДЕЛЕННОСТИ ДЛЯ ОЦЕНИВАНИЯ ПАРАМЕТРОВ ЛИНЕЙНЫХ ФУНКЦИЙ
При решении прикладных задач необходимо по набору экспериментальных данных, содержащих систематические и случайные ошибки, выбрать такую математическую модель, которая наиболее адекватно описывает исследуемое явление или процесс. Если изучаемый объект описывается линейной по параметрам функцией вида
y = F(x, a), (1)
то необходимо определить точечные и интервальные оценки параметров, а этой функции. Пусть ошибка измерения в каждом опыте составляет εi , то тогда истинные значения параметров функции (1) а* должны удовлетворять системе неравенств:
(2)
В работе [1, 2, 3] показано, что область Ω возможных значений параметров а при интервальной ошибке измерения представляет собой выпуклый многогранник.
В качестве точечной ошибки параметров можно использовать либо геометрический центр множества неопределенности, либо решение экстремальной задачи вида
(3)
При пpaктическом использовании метода центра неопределенности [МЦН] представляют интерес аппроксимации множества Ω простыми геометрическими фигурами сверху и снизу. В качестве таких фигур использовали прямоугольник и эллипс. В этом случае точечной оценкой параметров может служить геометрический центр выбранных фигур. Необходимо подчеркнуть, что разработанные алгоритмы обладают свойством рекуррентности, что позволяет по мере поступления информации об изучаемом объекте последовательно уточнять точечные оценки и уменьшать интервальные оценки.
В данной работе приведены результаты расчетов начальной концентрации исходного вещества и констант скорости реакции обратимой
g реакции первого порядка превращения γ-оксимасляной кислоты в соответствующей ланктон и воду [4]. В таблице 1 приведены исходные кинетические параметры этой реакции.
Таблица 1
i |
1 |
2 |
3 |
4 |
5 |
6 |
7 |
8 |
9 |
А(i) моль/л |
0.1582 |
0.1450 |
0.1324 |
0.1213 |
0.1115 |
0.1012 |
0.0923 |
0.0788 |
0.0668 |
ti, мин |
21 |
36 |
50 |
65 |
80 |
100 |
120 |
160 |
220 |
При проведении вычислительных экспериментов оценки погрешности измерений текущих концентраций вещества и времени варьировали в широких пределах. В таблице 2 приведены результаты расчетов параметров a, b линейной зависимости y = a + bx, где b = -(k1 + k2). Значение начальной концентрации определили по формуле где = 0,0495 моль/л - точное значение концентрации γ- оксимасляной кислоты. Расчеты проводили по методикам прямоугольника в МЦН, рекуррентного прямоугольника в МЦН (для двух случаев: при точном измерении входной переменной ti и неточном измерении входной переменной ti), взвешенного эллипса неопределенности в МЦН. Данные методики описаны в [5-7].
Таблица 2
|
A0 |
Размеры области неопределенности |
Метод |
|
2.016872 |
0.009414 |
0.1827896 |
1,34 • 10-7 |
Прямоугольник в МЦН |
2.038648 |
0.009217 |
0.182355 |
1,42 10-7 |
Рекуррентный прямоугольник в МЦН |
2.01807 |
0.00940 |
0.18241 |
1 • 10-7 |
Взвешенный эллипс в МЦН |
2.027299 |
0.009306 |
0.018119 |
4,17 • 10-7 |
МНК |
Непосредственно из экспериментов получены значения параметров: b = 0,0094, (-а) = 2,01807, которому соответствует начальная концентрация γ -оксимасляной кислоты А0 = 0,18241 моль/л. Таким образом
оценки параметров, полученные по предложенным методикам пpaктически совпадают с экспериментальным значениям. Все расчеты выполнены по программе «Аппроксимация экспериментальных данных».
СПИСОК ЛИТЕРАТУРЫ
- Вощинин А.П. Решение задач оптимизации по интервальным оценкам критерия // Зав. лаб. - 1987. - Т.53, № 7. С. 45-48.
- Вощинин А.П., Дывак Н.П. Планирование оптимального насыщенного эксперимента в задачах анализа интервальных данных // Завод. Лаб. - 1991. - Т.57, № 7. - С. 56-59.
- Вощинин А.П. Метод анализа данных с интервальными ошибками в задачах проверки гипотез и оценивания параметров неявных линейно-параметризованных функций // Завод. лаб.-2000.- Т.66, №3. - С. 51-65.
- Эмануэль Н.М., Кнорре Д.Г. Курс химической кинетики. М.: Высш. шк., 1984. - 463 с.
- Белов В.М., Гончаров С.А., Пролубников В.И., Унгер Ф.Г., Лукьянцева М.В. Алгоритмы прямоугольника в методе центра неопределенности для оценивания параметров линейных функций. -Томск, 2001. - 36 с. /Препринт ТНЦ СО РАН/.
- Белов В.М., Евстигнеев В.В., Суханов В.А., Лагуткина Е.В. Интервальная кинетика химических реакций. Методика взвешенного эллипса неопределенности в кинетике обратимых реакций первого порядка // Изв. вузов. Химия и хим. технология. - 1999. - Т.42, № 1.- С. 41 - 45.
- Белов В.М., Гончаров С.А., Гончарова Н.Л. Рекуррентный алгоритм оценивания параметров линейной двухпараметрической функции // Междунар. конференции «Физико-химические процессы в неорганических материалах», Ч.2. - Кемерово: Кузбасвузиздат, 2001. - С. 134-135.
Статья в формате PDF 126 KB...
28 03 2024 7:30:55
Исследованы особенности взаимосвязи увеличения массы и продольных размеров тела и его частей у здоровых новорожденных (535), детей с задержкой внутриутробного развития (938) и 221 детей с отставанием в росте одной из нижних конечностей. Использованы методы антропометрии детей, ультразвуковое исследование плода и оценка психического развития дошкольников. В этих группах обследуемых выявлены хаpaктерные различия динамики продольных размеров тела, отнесенных к его массе. Показано, что церебральный тип конституции появляется у детей при отсутствии диспластической задержки роста тела. Для успешного психического развития дошкольников благоприятны не максимальные, а средние размеры тела детей. ...
26 03 2024 16:55:23
Статья в формате PDF 122 KB...
24 03 2024 2:25:44
Статья в формате PDF 106 KB...
23 03 2024 10:19:24
Статья в формате PDF 120 KB...
21 03 2024 6:41:28
20 03 2024 0:35:37
Статья в формате PDF 253 KB...
19 03 2024 18:59:44
Статья в формате PDF 105 KB...
18 03 2024 5:22:57
Статья в формате PDF 101 KB...
17 03 2024 3:33:22
Статья в формате PDF 789 KB...
16 03 2024 1:57:26
Статья в формате PDF 124 KB...
15 03 2024 23:26:43
Статья в формате PDF 118 KB...
14 03 2024 1:59:55
13 03 2024 2:28:31
Статья в формате PDF 334 KB...
12 03 2024 17:27:33
Статья в формате PDF 230 KB...
11 03 2024 1:53:52
Статья в формате PDF 129 KB...
10 03 2024 22:33:17
09 03 2024 22:12:21
Статья в формате PDF 112 KB...
08 03 2024 6:31:34
Статья в формате PDF 120 KB...
07 03 2024 14:40:36
Статья в формате PDF 93 KB...
06 03 2024 10:14:44
Статья в формате PDF 146 KB...
05 03 2024 11:59:36
Построена октетная электродинамика. Обсуждена возможность объединения механики и электродинамики. Выявлена дальнодействующая структуризация октетного прострaнcтва. Исследуются свойства интервала. ...
04 03 2024 1:54:42
Статья в формате PDF 99 KB...
03 03 2024 22:28:40
Статья в формате PDF 312 KB...
02 03 2024 11:22:52
Статья в формате PDF 359 KB...
01 03 2024 10:56:13
Статья в формате PDF 120 KB...
28 02 2024 6:43:52
27 02 2024 10:12:12
26 02 2024 11:10:48
Статья в формате PDF 283 KB...
25 02 2024 17:23:43
Статья в формате PDF 112 KB...
23 02 2024 13:20:59
Статья в формате PDF 251 KB...
22 02 2024 2:44:55
Статья в формате PDF 112 KB...
20 02 2024 22:58:31
Статья в формате PDF 294 KB...
19 02 2024 21:57:37
Статья в формате PDF 346 KB...
18 02 2024 15:49:26
Еще:
Поддержать себя -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 ::