БИОТЕХНИЧЕСКИЙ ЗАКОН, ЭВРОРИТМ И АЛГОРИТМ ПОИСКА ПАРАМЕТРОВ > Полезные советы
Тысяча полезных мелочей    

БИОТЕХНИЧЕСКИЙ ЗАКОН, ЭВРОРИТМ И АЛГОРИТМ ПОИСКА ПАРАМЕТРОВ

БИОТЕХНИЧЕСКИЙ ЗАКОН, ЭВРОРИТМ И АЛГОРИТМ ПОИСКА ПАРАМЕТРОВ

Мазуркин П.М. Статья в формате PDF 182 KB

В основе вычислений параметров у исходной конструкции идентифицируемой модели [1] лежат методы случайного поиска. По мнению проф. Л.А. Растригина [2, с.405] методы случайного поиска являются прямым развитием метода «проб и ошибок», когда решение принимается случайно. При удаче приближения выходных результатов математической модели к результатам объекта исследования решение оставляется, а в ином случае - отвергается. Случайность поиска становится источником возможностей.

Наше уточнение состоит лишь в том, что сама случайность является проявлением неизвестной биотехнической закономерности. Поэтому поиск при случайных вариациях значений и состава параметров модели в программной среде типа CurveExpert-1.3 внутри ПЭВМ ведется закономерно итерациями (сеансами) от простого к сложному в конструкции модели. Постепенно, в ходе структурно-параметрической идентификации и эвристического осмысления получаемых результатов, происходит наращивание модели. Она структурируется по типу - от простого фрагмента биотехнического закона к сложной аддитивной структуре волновых составляющих с переменными амплитудами и частотами колебаний у готовой статистической модели.

Пусть функционирование объекта описывается математической моделью вида

,  (1)

где y - показатель (критерий) функционирования по выходному результату математической модели, x - множество объясняющих факторов (переменных), A - параметры модели, требующие определения.

Для модели (1) пусть известна исходная информация   и  по регистрации данных функционирования объекта исследования. Тогда задача идентификации сводится к задаче нелинейного программирования. Критерием качества (параметрической верификации) принимается сумма квадратов остатков, то есть

,  (2)

где K - целевая функция минимизации остатков, n - общее число наблюдений.

По модели (2) параметры А исходной модели принимаются за независимые переменные. При этом имеем множества

, (2а)

. (2б)

При этом справедливо соотношение, что ,(3) где  - остатки между выходными результатами объекта исследования и его модели.

Сущность поиска значений множества  сводится к тому, чтобы свести остатки   к минимально возможным значениям пори постепенном наращивании фрагмента биотехнического закона до полной конструкции. Затем происходит переход к следующей аддитивной составляющей.

Адекватность модели оценивается по остаткам, которые принимаются как единицы статистической совокупности. Средняя арифметическая остатков будет равна

.                                                                                    (4)

При  получим улучшение сходимости модели к фактическим данным.

Процедурная модель параметрической идентификации в специализированной программной среде РЕК (до 1600 точек, 25 факторов и 40 параметрах модели) приведена на рис. 1.

 

Рис. 1. Процедурная модель параметрической идентификации математической модели случайным поиском значений параметров

Из схемы видно, что исходная модель предполагается заданной, то есть до параметрической идентификации проведены этапы концептуальной и структурной идентификации моделей в виде эвроритмов. Процедурой 17 завершается создание готовой модели, которая в дальнейшем направляется для исследования влияния объясняющих переменных x  на показатели y . После такого исследования проводится обратная идентификация объекта исследования по модели, то есть выводы и рекомендации, полученные после моделирования и исследования готовой модели, внедряются в производство и учебный процесс с целью улучшения процесса функционирования объекта исследования. Так происходит итерация исследований. Процедуры 3, 8 и 12 выполняются исследователем на ЭВМ. Для этого была разработана программа ПЭК(РЕК), полностью отражающая содержание процедуры 3. Все остальные процедуры выполняются вручную, то есть в неавтоматизированном режиме. Причем процедуры 5, 6, 7, 9, 12, 13, 16 требуют применения эвристик. Например, процеДypa 13 предполагает два случая: а) изъятие малозначимых факторов и соответственно их массивов данных; б) трaнcформация (чаще всего в сторону упрощения) частных функций, что не приводит к изменению массива исходных данных. В первом случае далее выполняется процеДypa 15, а во втором - деятельность исследователя направляется к процедуре 7.

Сходимость модели по процедуре 7 оценивается по снижению статистических показателей остатков, сравненных относительно выходных результатов объекта по данным блока 4.

Сходимость считается удовлетворительной, если происходит «рысканье» статистического показателя точности относительно некоторого минимально достигнутого значения. Если такой тремор наблюдается, то исследователь переходит к выполнению процедуры 8, иначе к процедуре 12, в котором корректируются параметры гипер-параллелепипеда выбора стартовых значений .

Уравнения для расчета статистических показателей остатков приведены далее. При высокой адекватности возможно исключение тех наблюдений, которые отклоняются от границ доверительного интервала наблюдений. Причем в процедуре 16 выполняется содержательный анализ причин «аномальных» отклонений. Если дополнительные факторы выявить невозможно, то исследователь переходит к процедуре 11.

От процедуры 16, если найдено объяснение отклоняющимся точкам, возможно два случая:

а) отбросить резко отклоняющиеся точки и перейти к 15;

б) найденное содержательное объяснение выразить в виде дополнительного фактора с его частной функцией и затем перейти к процедуре 13 с целью усложнения структуры математической модели.

Из схемы на рис. 1 видно, что чем лучше подготовлена конструкция исходной модели, то тем быстрее происходит процесс перехода к готовой модели. Чем больше число факторов х, то тем труднее мысленно охватывать конструкцию модели, поэтому здесь необходимо выполнять принципы поискового иерархического проектирования конструкции (в данной книге не рассматривается), например, по методике [3].

Алгоритм поиска параметров модели. Инструмент параметрической идентификации находится в процедуре 3 по схеме на рис. 1, реализованный в виде программы ПЭК. В процеДypaх 8 и 11 нами использовались стандартные известные программы.

Пpaктически сеанс автоматического поиска содержит 5-15 циклов, а в каждом цикле принимается 300-2000 шагов поиска (в зависимости от числа факторов, числа наблюдений, а также заданного лимита машиносчетного времени). Программа ПЭК используется в итеративном режиме расчетов на ПЭВМ, что облегчает использование вычислительной машины и позволяет, в перерывах между сеансами, без спешки обдумывать эвристические приемы в соответствующих процеДypaх по схеме на рис. 1. В диалоговом режиме пpaктически на ПЭВМ возможна идентификация только относительно небольших моделей, сущность которых и выходной результат (готовая модель) уже априори предугаданы. Диалоговый режим эффективен для идентификации функционирования различных объектов. В математической среде «Эврика» для ПЭВМ возможна идентификация регрессионных моделей с числом наблюдений до 24 и числом факторов до 4 (версия 1.0 Eureka).

После завершения каждого цикла на печать или на дисплей выводятся значения: . Здесь: - улучшенные значения параметров модели, приведшие к приближению  ;  - число заданных шагов поиска; - число удачных шагов поиска; - сумма квадратов отклонений (остатков); - среднеквадратичное отклонение; - ошибка среднеквадратичного отклонения; - показатель изменчивости остатков относительно выходных результатов объекта исследования; - показатель точности (риска) сходимости модели относительно выходных результатов объекта исследования (100 -   -доверительная вероятность математической модели); - относительное отклонение (относительная погрешность) выходных результатов объекта   и модели y .

От цикла к циклу происходит изменение   при  . Первые (стартовые) значения выбираются случайно в некотором гипер-параллелепипеде.

,   (5)

где a,b - границы интервалов, задаваемые ориентировочно, N - общее число параметров модели. В программе РЕК для ПЭВМ принимались просто некоторые значения  , которые исследователю кажутся правдоподобными.


Внутри цикла поиск продолжается при снятых ограничениях (5). В последующих циклах в целях расширения стартовой области поиска принимается

 и  . (6)

Например, для автоматического переключения конфигурации гипер-параллелепипеда в стартах последующих циклов (после первого) поиска принимались

 и .  (6а)

После коррекции гипер-параллелепипеда новые значения   снова случайно принимаются по условию (5). С учетом итераций по формулам (5) и (6) от цикла к циклу меняется область задания старта исходных значений параметров модели. Поэтому происходит псевдоглобальный поиск значений параметров математической модели (1).

Для решения задачи (2) совместно с В.Г. Грудачевым был предложен алгоритм случайного поиска в «несвязных» подпрострaнcтвах, имеющий линейную тактику поиска. Случайный шаг выполняется после предыдущего неудачного, а в случае удачи шаги повторяются в одном и том же направлении.

Отличие этого алгоритма от известных алгоритмов случайного поиска [4, с.236-239] заключается в способе выбора приращений. Случайно выбирается целое число S, распределенное на отрезке [1,N]  и определяющее количество одновременно изменяемых параметров модели. Номера которых определяются также случайно. Предложенный алгоритм случайного поиска можно записать в виде следующих рекуррентных выражений:


,   (7)

 при  ,   (8)

где h - число последовательно неудачных шагов поиска.

Синтез приращений  происходит следующим образом:

                                (9)

где: , ,  - векторы приращений соответственно на t-1-м, t -м, t+1 -м шаге поиска; K(A t-1), K(A t) - значения критериев идентификации после осуществления t-1  и  t  шагов поиска; c - вектор максимального рабочего шага поиска, c={c j} ;

 

-  (10)

вектор случайных чисел, равномерно распределенных на отрезке [-1;-u;u;1] ; u -коэффициент, учитывающий минимально допустимое значение длины рабочего шага, в нашем случае принималось u =0.001; p - случайный номер параметра математической модели, упорядоченный в интервале [1;N]  соотношением, при s=1,N ,

.    (11)

Первые (стартовые) значения параметров математической модели задаются или случайным образом выбираются в заданной стартовой области поиска (5) с помощью соотношения

, (12)

 

где V j - вектор случайных чисел в интервале [u,1] .

Массив исходных данных (процеДypa 2 по схеме на рис. 2) включает: массив значений объясняющих факторов   и выходных результатов  объекта исследования, массивы границ стартовой области поиска a j и b j , массив вектора максимального рабочего шага поиска c j , а также значения числа наблюдений n , границ структуры математической модели N  и m , границ сеанса поиска l  (число циклов) и цикла h 0  (число принятых шагов поиска).

О второго цикла поиска значения вектора максимального рабочего шага вычисляются по формуле

.  (13)


для которого мы принимали α  = 0,1.

Анализ сходимости и адекватности модели. Случайность признана жизненным атрибутом любого явления [5, с.12]. Регрессия - это односторонняя стохастическая зависимость, выраженная с помощью функции [5, с.15]. В регрессионном анализе исследуется форма связи. А в корреляционном анализе - сила стохастической связи [5, с.18]. Если структура регрессионной модели подобрана верно, то необходимость в корреляционном анализе в методике МЭРА отпадает.

Сила связей между факторами определяется на основе исследования готовой модели.

Сходимость - это свойство приближаемости теоретической линии регрессии к выходным результатам объекта исследования. Критерием для оценки сходимости является время поиска параметров модели, или общее число шагов поиска по всем сеансам. Поиск прекращается после того, как сходимость модели по значениям её параметров прекращается, что можно оценить уменьшением остатков. Если остатки после нескольких циклов поиска начинают колeбaться около одного и того же значения, то поиск прекращается. Далее по конечным остаткам оценивается адекватность модели процессу функционирования объекта исследования.

Такой подход эффективен для решения многих экономических задач. Как отмечается в работе [5, с.154], множественная нелинейная регрессия лучше отражает многообразие связей в экономике. Применение ПЭВМ с высоким быстродействием и большой памятью снимает все вычислительные проблемы, которые раньше были препятствием для создания нелинейных многофакторных моделей.

При применении метода наименьших квадратов (МНК) к оценке адекватности готовой модели к функционированию объекта исследования (а промежуточные значения показателей МНК - к оценке сходимости) будем считать, что стратегия, выбирающая наиболее эффективную проверку, всегда оптимальна. Такая эффективная проверка существует в общей теории статистики [68], в соответствии с которой ниже приведены расчетные показатели анализа адекватности выходных результатов готовой математической модели и объекта исследования.

Из сходных предпосылок получим формулу для расчета остатков

,(14)

где  - выходной результат объекта (фактические значения показателя), y i - выходной результат математической модели (теоретические значения показателя), i - номер наблюдения, n - общее число наблюдений.

Причем теоретическое значение показателя оценивается по формуле (1). Все значения остатков   образуют статистическую выборку для следующей составляющей искомой модели при условии, что данная составляющая доведена по сложности конструкции до полной структуры по биотехническому закону как по амплитуде, так и по половине периода колебательного возмущения.

Статья опубликована при поддержке гранта 3.2.3/4603 МОН РФ


СПИСОК ЛИТЕРАТУРЫ:

  1. Мазуркин, П.М. Статистическое моделирование. Эвристико-математический подход / П.М. Мазуркин. - Научное издание. - Йошкар-Ола: МарГТУ, 2001. - 100 с.
  2. Растригин, Л.А. Системы экстремального управления / Л.А. Растригин. - М.: Наука. 1974. - 632 с.
  3. Капустян, В.М. Конструктору о конструировании атомной техники / В.М. Капустян, Ю.А. Махотенко. - М.: Атомиздат, 1981. - 190 с.
  4. Автоматизация поискового конструирования (искусственный интеллект в машинном проектировании) / Под ред. А.И. Половинкина. - М.: Радио и связь, 1981. - 344 с.
  5. Фёрстер, Э. Методы корреляционного и регрессионного анализа: Руководство для экономистов / Э. Фёрстер, Б. Рёнц. - М.: Экономика и статистика, 1983. - 302 с.


Продажа товаров в кредит

Продажа товаров в кредит Статья в формате PDF 113 KB...

21 04 2024 21:17:23

ГИГИЕНИЧЕСКАЯ ОЦЕНКА ОПАСНОСТИ ЗАГРЯЗНЕНИЯ ПОЧВЫ СВИНЦОМ

ГИГИЕНИЧЕСКАЯ ОЦЕНКА ОПАСНОСТИ ЗАГРЯЗНЕНИЯ ПОЧВЫ СВИНЦОМ Повышение уровня свинца в крови у детей дошкольного возраста на 1 мкг/дл ведет к снижению интеллектуального развития ребенка на 1/4–1/2 балла, причем негативные последствия обнаруживаются и через 10 лет после воздействия свинца в раннем возрасте. Целью данного исследования было дать гигиеническую оценку загрязнения почвы в качестве депонирующей системы свинцом в городе Шымкент. Для исследования почвы на содержание свинца был произведен забор проб согласно ГОСТу 17.4.02-84. Определение свинца проводили на атомно-абсорбционном спектрометре МГА-915 с электротермической атомизацией. В результате исследования установлено превышение содержания свинца в почве по отношению к ПДК во всех отобранных пробах. Причем, по мере удаления от завода концентрация свинца в почве хотя и уменьшалось, но превышало ПДК в 3–4 раза. При исследовании овощей произрастающих на загрязненной территории, максимальное содержание свинца установлено в картофеле (в среднем 3 ПДК в 1 зоне). Таким образом, полученные результаты показали, что наибольшее загрязнение наблюдается на расстоянии 500–1000 м от завода. Вместе с тем обнаружено загрязнение почвы по всей изучаемой территории, где складывается нeблагоприятная санитарная ситуация по свинцу. ...

19 04 2024 3:15:21

ИЗУЧЕНИЕ ВЕКТОРОВ В ШКОЛЕ. ВОПРОСЫ И ОТВЕТЫ

ИЗУЧЕНИЕ ВЕКТОРОВ В ШКОЛЕ. ВОПРОСЫ И ОТВЕТЫ Статья в формате PDF 250 KB...

17 04 2024 18:28:27

ПРОГНОЗИРОВАНИЕ ЭКОНОМИЧЕСКОЙ ДЕЯТЕЛЬНОСТЬЮ ПРЕДПРИЯТИЯ И УПРАВЛЕНИЕ ЕГО РАЗВИТИЕМ

ПРОГНОЗИРОВАНИЕ ЭКОНОМИЧЕСКОЙ ДЕЯТЕЛЬНОСТЬЮ ПРЕДПРИЯТИЯ И УПРАВЛЕНИЕ ЕГО РАЗВИТИЕМ На основе системного анализа функционирования экономической деятельности промышленного предприятия введена его теоретическая кривая прогнозирования бизнеса и разработан алгоритм выхода на данную кривую в процессе стратегического управления развитием предприятия. ...

08 04 2024 21:23:38

ИЗВЛЕЧЕНИЕ ФЛАВОНОИДОВ ИЗ ПИЖМЫ ОБЫКНОВЕННОЙ

ИЗВЛЕЧЕНИЕ ФЛАВОНОИДОВ ИЗ ПИЖМЫ ОБЫКНОВЕННОЙ Статья в формате PDF 193 KB...

02 04 2024 2:31:53

ОЦЕНКА ВЛИЯНИЯ НЕБЛАГОПРИЯТНОЙ ЭКОЛОГО-ГИГИЕНИЧЕСКОЙ ОБСТАНОВКИ ЮЖНО-КАЗАХСТАНСКОЙ ОБЛАСТИ НА ЗАБОЛЕВАЕМОСТЬ МИГРАНТОВ-РЕПАТРИАНТОВ

ОЦЕНКА ВЛИЯНИЯ НЕБЛАГОПРИЯТНОЙ ЭКОЛОГО-ГИГИЕНИЧЕСКОЙ ОБСТАНОВКИ ЮЖНО-КАЗАХСТАНСКОЙ ОБЛАСТИ НА ЗАБОЛЕВАЕМОСТЬ МИГРАНТОВ-РЕПАТРИАНТОВ На здоровье населения особое влияние оказывают экологические, гигиенические, социально-медицинские причины. В работе была реализована специально созданные социологические карты. Результаты социологического исследования показали, что к причинам, сильно влияющим на здоровье мигрантов-репатриантов относятся экологически нeблагоприятные условия окружающей среды. Заболеваемость мигрантов-репатриантов, проживающих в высокой степени опасности экологически наблагоприятных районах достигает от 2227,9 до 3010,9 ‰. Этот показатель указывает на значительное повышение показателей мигрантов, проживающих районах, где экологическая обстановка средняя, низкая и неопасная .Между загрязнением атмосферного воздуха и почвы и патологиями иммунной системы, минерализацией воды и заболеваниями мочепoлoвoй системы, загрязнением атмосферного воздуха и патологиями дыхательных путей есть прямая и в высокой степени связь. ...

01 04 2024 18:14:19

АНОМАЛИИ ПЛАНЕТ СОЛНЕЧНОЙ СИСТЕМЫ

АНОМАЛИИ ПЛАНЕТ СОЛНЕЧНОЙ СИСТЕМЫ С целью проверки космологических и геологических теорий всё больший интерес вызывают измерения аномалий: увеличение радиусов орбит планет, увеличение радиусов планет, замедление вращения планет. Технические возможности таких измерений имеются. Эмпирическая Теория Вселенной позволяет легко вычислять указанные аномалии. В статье показан метод расчёта аномалий и некоторые результаты для планет Солнечной системы. Сравнение расчёта с уже имеющимися измерениями (удаление Луны от Земли, удаление Земли от Солнца, замедление вращения Земли) показывает хорошее согласие расчёта и измерения. ...

30 03 2024 17:40:56

СОСТОЯНИЕ ОКРУЖАЮЩЕЙ СРЕДЫ В АРКТИКЕ

СОСТОЯНИЕ ОКРУЖАЮЩЕЙ СРЕДЫ В АРКТИКЕ В Арктике масштабы деградации окружающей среды приобретают опасные тенденции, нарушение хрупкой арктической природы может иметь необратимый хаpaктер. Анализ данных официальных источников показал, что к территориям «риска» по загрязнению питьевой воды относятся Ямало-Ненецкий автономный округ и Республика Саха. Высокий уровень загрязнения атмосферного воздуха зарегистрирован в Красноярском крае, а самые высокие показатели загрязнения почвы показаны в Мурманской области. ...

29 03 2024 13:33:45

ХЕМОКИНЕТИКОТРОПНЫЕ СВОЙСТВА АТФ

ХЕМОКИНЕТИКОТРОПНЫЕ СВОЙСТВА АТФ Статья в формате PDF 88 KB...

22 03 2024 12:31:10

МЕДИЦИНСКАЯ ГЕОГРАФИЯ БИОГЕЛЬМИНТОЗОВОПИСТОРХОЗА

МЕДИЦИНСКАЯ ГЕОГРАФИЯ БИОГЕЛЬМИНТОЗОВОПИСТОРХОЗА Статья в формате PDF 227 KB...

21 03 2024 12:46:29

ЛИТВИНА ЛИДИЯ АЛЕКСЕЕВНА

ЛИТВИНА ЛИДИЯ АЛЕКСЕЕВНА Статья в формате PDF 283 KB...

18 03 2024 2:51:40

ЭЛЕКТРОМАГНИТНЫЙ ФАКТОР ПЛАНЕТЫ

ЭЛЕКТРОМАГНИТНЫЙ ФАКТОР ПЛАНЕТЫ Статья в формате PDF 190 KB...

16 03 2024 10:56: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 ::