ИССЛЕДОВАНИЕ ЗАВИСИМОСТИ КОЭФФИЦИЕНТОВ ЛИНЕЙНОЙ МОДЕЛИ ПАРНОЙ РЕГРЕССИИ ОТ ПАРАМЕТРОВ ФУНКЦИИ ВЕЙЕРШТРАССА-МАНДЕЛЬБРОТА ПРИ МОДЕЛИРОВАНИИ СРЕДНЕГО ЗНАЧЕНИЯ ФУНКЦИИ С ПОМОЩЬЮ ЕЕ ФРАКТАЛЬНОЙ РАЗМЕРНОСТИ
Введение
При моделировании микроускорений с помощью действительной части фpaктальной функции Вейерштрасса-Maндельброта (ФВМ) [1]:
(1)
важно грамотно провести отождествление параметров функции (1) и реальных условий процесса возникновения поля микроускорений внутри КА.
Рисунок 1. Динамика изменения среднего значения ФВМ
Этой теме, равно как и возможности использовании ФВМ в виде (1) при тождественно равной нулю случайной фазе, посвящен ряд работ [2-5], с которыми можно ознакомиться для понимания решаемой задачи. В работе [4] были получены корреляционные зависимости между средним значением ФВМ (1) и фpaктальной размерностью D (рис. 1).
Как видно из рис. 1, они пpaктически линейны. Этот факт доказывается с помощью коэффициента детерминации, который при оценке качества моделирования корреляционных зависимостей рис. 1 во всех случаях превышает 0,999 (рис. 2).
Как видно из рис. 1, они пpaктически линейны. Этот факт доказывается с помощью коэффициента детерминации, который при оценке качества моделирования корреляционных зависимостей рис. 1 во всех случаях превышает 0,999 (рис. 2).
Рисунок 2. Зависимость коэффициента детерминации от параметра b при моделировании корреляционных зависимостей рис.1 линейной моделью парной регрессии
Постановка задачи
Для построения функциональной зависимости между фpaктальной размерностью ФВМ D и средним значением ФВМ (1) с помощью линейной модели парной регрессии вида:
(2)
требуется исследовать влияние параметра b ФВМ на коэффициенты и правой части (2), исходя из полученных ранее корреляционных зависимостей рис. 1.
Основные результаты работы
Как видно из рис. 1, оба исследуемых коэффициента с ростом b изменяются. Для аппроксимации корреляционной зависимости и b (рис. 3)
Рисунок 3. Корреляционная зависимость коэффициента а1 от параметра b
Была сначала построена линейная модель:
(3)
и с помощью метода наименьших квадратов (МНК) оценены коэффициенты и , которые получились равными: ; . Таким образом, наилучшая с точки зрения МНК линейная зависимость коэффициента от b имеет вид:
(4)
Затем была предпринята попытка улучшить качество аппроксимации за счет учета квадратичного члeна, заменив зависимость (3) на квадратичную вида:
, (5)
для которой также с помощью МНК были оценены коэффициенты: ; ; . Таким образом, наилучшая с точки зрения МНК квадратичная зависимость коэффициента от b имеет вид:
(6)
И, наконец, был произведен учет кубического члeна с помощью зависимости:
, (7)
При этом коэффициенты, найденные с помощью МНК, оказались равными: ; ; ; .
Таким образом, наилучшая с точки зрения МНК кубическая зависимость коэффициента от b имеет вид:
(8)
Качество построенных зависимостей (4), (6) и (8) проверялось с помощью коэффициента детерминации (рис. 4) и критерия согласия -Пирсона (рис. 5).
Рисунок 4. Изменение коэффициента детерминации при усложнении формы модели зависимости коэффициента а1 от параметра b
Усложнение модели связано, прежде всего с высоким качеством аппроксимации корреляционных зависимостей рис. 1 моделью (2). Значения коэффициентов детерминации для линейной и квадратичной моделей составили 0,979 и 0,998 соответственно, что ниже значений коэффициентов детерминации, рассчитанных при моделировании (рис. 2). И только модель (8) позволяет с уверенностью сделать вывод о том, что потерь качества при замене коэффициента на правую часть (8) не будет, т.к. значение коэффициента детерминации составляет 0,9995.
Рисунок 5. Динамика изменения наблюдаемого значения критерия согласия при усложнении модели
Критерий согласия еще более наглядно показывает улучшение качества моделирования при усложнении формы модели. Критическое значение критерия для рассматриваемых двух степеней свободы равно 5,99147 (5 %-й уровень значимости). Наблюдаемое значение критерия для линейной модели (4) составляет 5,613, что очень близко к критическому, несмотря на высокое значение коэффициента детерминации. Для квадратичной модели (6) величина наблюдаемого значения критерия сокращается более, чем в десять раз: 0,505 и почти в четыре раза сокращается еще для кубической модели (8): 0,130.
Таким образом, проверка качества построенных моделей (4), (6) и (8) показала, что лишь последняя из них достаточно точно описывает динамику изменения коэффициента в модели (2).
Следующим этапом является исследование зависимости другого параметра модели (2) от b. При замене корреляционных зависимостей рис.1 на функциональные вида (2) была построена зависимость изменения при различных значениях b. Эта зависимость приведена на рис. 6.
Рисунок 6. Корреляционная зависимость коэффициента а0 от параметра b
Аналогично схеме построения зависимости от b, сначала была исследована линейная зависимость, подобная (3). С помощью МНК подобраны коэффициенты этой модели:
(9)
Затем были подобраны лучшие с точки зрения МНК коэффициенты квадратичной зависимости, подобной (5):
(10)
и кубической зависимости, подобной (7):
(11)
Качество построенных моделей (9), (10) и (11) оценивалось с помощью коэффициента детерминации (рис. 7) и критерия согласия -Пирсона (рис. 8).
Рисунок 7. Изменение коэффициента детерминации при усложнении формы модели зависимости коэффициента а0 от параметра b
Рисунок 8. Динамика изменения наблюдаемого значения критерия согласия при усложнении модели
Таким образом, в результате данной работы построена модель зависимости среднего значения ФВМ от параметров функции:
Причем, данная модель объясняет более 99,9% дисперсии корреляционных зависимостей рис.
СПИСОК ЛИТЕРАТУРЫ
- Седельников А.В., Бязина А.В., Антипов Н.Ю. Использование функции Вейерштрасса-Maндельброта для моделирования микроускорений на борту КА //Сборник научных трудов X Всероссийского научно-технического семинара по управлению движением и навигации ЛА. Самара. 2002. с. 124-128.
- Седельников А.В., Корунтяева С.С., Чернышева С.В. Анализ влияния параметров функции Вейерштрасса-Maндельброта на ее закон распределения //Современные наукоемкие технологии. - 2005 г. - № 9. - с. 43-46.
- Седельников А.В., Бязина А.В., Иванова С.А. Статистические исследования микроускорений при наличии слабого демпфирования колебаний упругих элементов КА //Научные чтения в Самарском филиале РАО. - Часть 1. Естествознание. - М.: Изд. УРАО. - 2003. - 137 - 158.
- Седельников А.В., Корунтяева С.С., Подлеснова Д.П. Исследование динамики изменения среднего значения фpaктальной функции Вейерштрасса-Maндельброта как случайной величины //Фундаментальные исследования. - № 4. - 2006. - с. 84-87.
- Седельников А.В., Корунтяева С.С., Чернышева С.В. Выявление коридора значений параметров фpaктальной функции Вейерштрасса-Maндельброта, при которых справедлив нормальный закон распределения функции //Современные наукоемкие технологии. - № 1. - 2006. - с. 85-87.
Статья в формате PDF 228 KB...
24 04 2024 16:40:56
Статья в формате PDF 107 KB...
22 04 2024 23:55:56
Статья в формате PDF 143 KB...
21 04 2024 21:43:57
Статья в формате PDF 268 KB...
20 04 2024 6:45:37
19 04 2024 3:56:10
Статья в формате PDF 576 KB...
18 04 2024 23:53:47
Статья в формате PDF 119 KB...
16 04 2024 5:29:11
Статья в формате PDF 104 KB...
15 04 2024 13:43:21
С помощью комплекса ядерно-физических методов, ЯМР-спектроскопии, выявлена неоднозначная степень насыщения связанной фазы воды молекулами воды и ряда химических элементов, где основу их специфической связи представляет многослойная поляризованная структура сыворотки крови и лимфы здоровых людей, пациентов с актуальными заболеваниями. Разработана иерархическая двухуровневая модель, согласно собственной концепции сопряженного действия и эффекта энергии, системного ЭМП, энергии биохимических цикловых процессов, объединенных потоком протонов, регулируемых буферной системой и гормонами стресса. ...
14 04 2024 23:42:44
Статья в формате PDF 111 KB...
13 04 2024 16:40:47
Статья в формате PDF 103 KB...
10 04 2024 9:11:56
Статья в формате PDF 106 KB...
09 04 2024 2:51:49
Статья в формате PDF 122 KB...
08 04 2024 12:47:33
Статья в формате PDF 119 KB...
07 04 2024 15:29:24
Рассмотрена экономико-математическая модель конкуренции двух фирм на однородном рынке сбыта. Приводится формулировка соответствующей задачи Коши для системы обыкновенных дифференциальных уравнений первого порядка, описывающей динамику развития системы, которая может быть легко обобщена на случай произвольного количества конкурирующих предприятий. Дана экономическая интерпретация полученных результатов. ...
06 04 2024 4:10:23
Статья в формате PDF 151 KB...
05 04 2024 0:47:46
Статья в формате PDF 257 KB...
04 04 2024 0:58:11
Статья в формате PDF 104 KB...
03 04 2024 20:22:55
Статья в формате PDF 121 KB...
02 04 2024 23:22:38
Статья в формате PDF 293 KB...
01 04 2024 22:15:26
Статья в формате PDF 294 KB...
30 03 2024 8:41:22
Статья в формате PDF 120 KB...
29 03 2024 13:39:38
Статья в формате PDF 280 KB...
28 03 2024 5:13:10
Статья в формате PDF 124 KB...
27 03 2024 1:59:55
Статья в формате PDF 100 KB...
25 03 2024 4:29:37
Статья в формате PDF 237 KB...
24 03 2024 8:42:54
Статья в формате PDF 111 KB...
23 03 2024 18:37:24
Статья в формате PDF 115 KB...
22 03 2024 7:21:17
Статья в формате PDF 131 KB...
20 03 2024 15:55:39
Статья в формате PDF 243 KB...
19 03 2024 0:44:57
Статья в формате PDF 732 KB...
18 03 2024 8:20:50
Статья в формате PDF 114 KB...
17 03 2024 3:51:48
Статья в формате PDF 120 KB...
16 03 2024 4:50:56
Еще:
Поддержать себя -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 ::