ПРИМЕНЕНИЕ МЕТОДА ЦЕНТРА НЕОПРЕДЕЛЕННОСТИ ДЛЯ ОЦЕНИВАНИЯ ПАРАМЕТРОВ ЛИНЕЙНЫХ ФУНКЦИЙ
При решении прикладных задач необходимо по набору экспериментальных данных, содержащих систематические и случайные ошибки, выбрать такую математическую модель, которая наиболее адекватно описывает исследуемое явление или процесс. Если изучаемый объект описывается линейной по параметрам функцией вида
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
314 KB...
20 09 2026 4:56:44
Статья в формате PDF 110 KB...
18 09 2026 8:54:24
Статья в формате PDF
134 KB...
17 09 2026 7:50:49
Статья в формате PDF
361 KB...
16 09 2026 5:52:16
Статья в формате PDF
102 KB...
15 09 2026 1:11:26
Статья в формате PDF
131 KB...
14 09 2026 5:31:54
Статья в формате PDF
140 KB...
13 09 2026 23:40:47
Статья в формате PDF
102 KB...
12 09 2026 10:20:36
Статья в формате PDF
116 KB...
10 09 2026 8:10:27
Статья в формате PDF
131 KB...
09 09 2026 16:50:56
Статья в формате PDF
208 KB...
06 09 2026 0:23:54
Гравитационные силы обусловлены тем, что в материальные тела поступает энергия из космического прострaнcтва, которая создает давление и увеличивает массу тел. Гипотеза находит подтверждение в виде космологического красного смещения. Возникновение инерционных сил (вопреки теории относительности А. Эйнштейна) наступает вследствие взаимодействия элементарных частиц с эфиром. Проанализирована структура электрона, и на ее основе проведена оценка скорости гравитационных волн, которая оказалась равной 4.7∙108 м/с.
...
05 09 2026 7:23:57
Статья в формате PDF
92 KB...
04 09 2026 15:40:40
Статья в формате PDF
108 KB...
03 09 2026 17:29:54
Статья в формате PDF
125 KB...
02 09 2026 17:20:20
Статья в формате PDF
110 KB...
31 08 2026 20:59:33
Статья в формате PDF
92 KB...
30 08 2026 4:28:39
Статья в формате PDF
268 KB...
29 08 2026 16:56:36
Статья в формате PDF
125 KB...
28 08 2026 0:22:58
Статья в формате PDF
217 KB...
27 08 2026 12:56:48
Статья в формате PDF
250 KB...
26 08 2026 4:16:26
Статья в формате PDF
261 KB...
25 08 2026 20:34:15
Статья в формате PDF
248 KB...
24 08 2026 9:37:48
Статья в формате PDF
127 KB...
23 08 2026 20:40:44
Статья в формате PDF
236 KB...
20 08 2026 15:40:33
Статья в формате PDF
117 KB...
19 08 2026 8:33:20
Статья в формате PDF
254 KB...
18 08 2026 23:53:13
На основе анализа s-d обменного взаимодействия в структурах типа NiAs с частично вакантными катионными позициями, моделировались различного рода зависимости результирующей намагниченности от температуры нестехиометрических ферримагнетиков. На основе исследований пирротина методами ЯГР и РФА доказано, что двухподрешеточный ферримагнетик, содержащий в структуре катионные вакансии, должен рассматриваться, при определенном типе распределения вакансий, как ферримагнетик с четырьмя магнитными подрешетками. В данном случае, дополнительные магнитные подрешетки можно рассматривать как подрешетки, индуцированные хаpaктером распределения катионных вакансий в структуре. Квантово-механические расчеты в рамках модели молекулярного поля температурных изменений намагниченности отдельно для каждой из подрешеток, а также анализ результирующей термокривой намагниченности, объясняют ряд экспериментально полученных кривых зависимости намагниченности от температуры нестехиометрического пирротина с различной плотностью вакансий в структуре.
...
17 08 2026 19:55:51
Статья в формате PDF
127 KB...
16 08 2026 19:59:29
Статья в формате PDF
286 KB...
15 08 2026 0:16:13
Статья в формате PDF
256 KB...
14 08 2026 6:18:18
Статья в формате PDF
123 KB...
13 08 2026 10:49:42
Еще:
Поддержать себя -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 ::
