ПРИМЕНЕНИЕ МЕТОДА ЦЕНТРА НЕОПРЕДЕЛЕННОСТИ ДЛЯ ОЦЕНИВАНИЯ ПАРАМЕТРОВ ЛИНЕЙНЫХ ФУНКЦИЙ

При решении прикладных задач необходимо по набору экспериментальных данных, содержащих систематические и случайные ошибки, выбрать такую математическую модель, которая наиболее адекватно описывает исследуемое явление или процесс. Если изучаемый объект описывается линейной по параметрам функцией вида
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
263 KB...
21 07 2026 19:32:22
Статья в формате PDF
252 KB...
20 07 2026 21:17:26
Статья в формате PDF
320 KB...
19 07 2026 20:40:43
В статье представлены актуальные данные о проблеме урогeнитaльного xлaмидиоза. Рассмотрены современные вопросы эпидемиологии, патогенеза и терапии инфекции.
...
18 07 2026 21:12:43
Статья в формате PDF
317 KB...
17 07 2026 6:40:25
Статья в формате PDF
104 KB...
15 07 2026 15:29:34
14 07 2026 1:50:57
Статья в формате PDF
241 KB...
13 07 2026 4:17:14
Статья в формате PDF
254 KB...
12 07 2026 6:42:29
Статья в формате PDF
115 KB...
11 07 2026 0:35:40
Статья в формате PDF
122 KB...
10 07 2026 18:48:30
Статья в формате PDF
121 KB...
09 07 2026 3:51:32
Статья в формате PDF
101 KB...
08 07 2026 23:17:42
Статья в формате PDF
163 KB...
07 07 2026 21:13:59
Статья в формате PDF
134 KB...
06 07 2026 15:39:10
Статья в формате PDF
101 KB...
04 07 2026 11:34:43
Статья в формате PDF 120 KB...
03 07 2026 20:24:45
Статья в формате PDF
122 KB...
02 07 2026 23:24:33
Статья в формате PDF
120 KB...
01 07 2026 9:10:57
Статья в формате PDF
218 KB...
30 06 2026 19:49:18
В работе рассмотрены термодинамические аспекты люминесцентного газового анализа. Молекулы красителя, адсорбированные на поверхности пористого вещества или внедренные в полимерную пленку, рассматриваются как система невзаимодействующих частиц, погруженная в термостат. Для относительной интенсивности флюоресценции молекул красителя получена связь с основной термодинамической хаpaктеристикой термостата – энергией Гиббса. Определены термодинамические ограничения точности газового анализа. Показано, что оптимальной основой для люминесцентного анализатора является полимерная пленка с наименьшим значением поверхностного натяжения.
...
29 06 2026 14:50:35
Статья в формате PDF
253 KB...
27 06 2026 10:31:45
Статья в формате PDF
103 KB...
25 06 2026 12:31:46
Статья в формате PDF
112 KB...
24 06 2026 19:38:44
Статья в формате PDF
108 KB...
23 06 2026 23:12:29
Статья в формате PDF
109 KB...
22 06 2026 18:36:39
Обсуждены методика и некоторые результаты моделирования вероятных конфигураций межфазных границ на поверхности композиционных материалов, полученные методом итерации прямоугольных генераторов на определенных сетках Кеплера-Шубникова.
...
21 06 2026 20:22:33
Статья в формате PDF
111 KB...
20 06 2026 0:40:18
Статья в формате PDF
110 KB...
19 06 2026 10:30:14
Статья в формате PDF
224 KB...
18 06 2026 20:33:40
Статья в формате PDF
107 KB...
17 06 2026 9:52:24
Статья в формате PDF
162 KB...
16 06 2026 0:19:28
Статья в формате PDF
111 KB...
14 06 2026 15:48:22
Статья в формате PDF
118 KB...
13 06 2026 19:36:40
Статья в формате PDF
119 KB...
12 06 2026 20:28:10
Еще:
Поддержать себя -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 ::