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

В основе вычислений параметров у исходной конструкции идентифицируемой модели [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 МОН РФ
СПИСОК ЛИТЕРАТУРЫ:
- Мазуркин, П.М. Статистическое моделирование. Эвристико-математический подход / П.М. Мазуркин. - Научное издание. - Йошкар-Ола: МарГТУ, 2001. - 100 с.
- Растригин, Л.А. Системы экстремального управления / Л.А. Растригин. - М.: Наука. 1974. - 632 с.
- Капустян, В.М. Конструктору о конструировании атомной техники / В.М. Капустян, Ю.А. Махотенко. - М.: Атомиздат, 1981. - 190 с.
- Автоматизация поискового конструирования (искусственный интеллект в машинном проектировании) / Под ред. А.И. Половинкина. - М.: Радио и связь, 1981. - 344 с.
- Фёрстер, Э. Методы корреляционного и регрессионного анализа: Руководство для экономистов / Э. Фёрстер, Б. Рёнц. - М.: Экономика и статистика, 1983. - 302 с.
Статья в формате PDF
124 KB...
02 05 2026 8:23:55
Статья в формате PDF
300 KB...
30 04 2026 13:14:46
Статья в формате PDF
122 KB...
28 04 2026 16:39:41
Выбрать оптимальный метод введения больных в период реабилитации после черепно-мозговой травмы. Материалы и методы: За 2011 год в Новокуйбышевской центральной городской больницы пролечено 960 пострадавших с черепно-мозговой травмой, из них 780 пострадавших с сотрясением головного мозга. Все пациенты с сотрясением головного мозга, первых семь дней находились на стационарном лечении в условиях травматологического отделения. Под наблюдением врачей нейрохирурга, травматолога, невролога и окулиста, проводилась дегидратационная и симптоматическая терапия. После первой недели стационарного лечения данных пациентов разделили на три равных группы по 260 человек и в дальнейшем их вели по- разному. Результаты: Удовлетворительные результаты лечения получены в первой группе у 252 пациентов (97%), у второй группы 243 пациентов(93%), а в третьей 156 пациентов (60%). Один день дневного стационара в травматологическом отделение в Новокуйбышевской центральной городской больницы НЦГБ стоит 360 рублей, а один день дневного стационара, стоит 190 рублей. Таким образом стоимость лечения пациентов первой группы = (7 + 7)·360 = 5040 рублей, стоимость лечения пациентов второй группы = 7·360 + 7·190 = 2520 + 1330 = 3850 рублей, стоимость лечения пациентов третьей группы = 7·360 = 2520 рублей. Из данных расчетов видно, что пациенты третьей группы, требует меньше расходов, но к сожалению, у них намного хуже результаты лечения. Результаты лечения пациентов первой и второй группы пpaктически одинаковы, а стоимость пациентов второй группы намного меньше.
...
27 04 2026 11:50:36
26 04 2026 20:57:51
Статья в формате PDF
119 KB...
25 04 2026 3:11:33
Статья в формате PDF
204 KB...
23 04 2026 15:42:31
Статья в формате PDF
154 KB...
22 04 2026 14:25:17
Статья в формате PDF
154 KB...
21 04 2026 1:56:24
Статья в формате PDF
129 KB...
20 04 2026 19:48:55
Статья в формате PDF
148 KB...
19 04 2026 3:36:36
Статья в формате PDF
111 KB...
18 04 2026 5:50:57
Статья в формате PDF
137 KB...
17 04 2026 14:52:36
Статья в формате PDF
198 KB...
15 04 2026 3:34:43
Статья в формате PDF
262 KB...
14 04 2026 17:26:53
Статья в формате PDF
124 KB...
13 04 2026 20:11:39
Статья в формате PDF
114 KB...
12 04 2026 15:25:22
Статья в формате PDF
124 KB...
11 04 2026 18:39:51
Статья в формате PDF
115 KB...
10 04 2026 19:10:51
Статья в формате PDF
108 KB...
09 04 2026 2:43:48
Статья в формате PDF
119 KB...
08 04 2026 19:41:32
Изучены онтогенез и возрастная структура ценопопуляций многолетних травянистых поликарпических видов, относящихся к различным типам экобиоморф: стержнекорневых – дягиль лекарственный (Angenica archangelica L.) и цикорий обыкновенный (Cichorium intybus L.) и длиннокорневищных – вязель разноцветный (Coronilla varia L.).В онтогенезе выбранных видов выделены следующие 4 периода и 9 возрастных состояний: 1). период первичного покоя (покоящиеся семена); 2). виргинильный период (проростки, ювенильное, имматурное, виргинильное); 3). генеративный (молодое, средневозрастное, старое генеративное); 4). сенильный (сенильное).
Изучение возрастной структуры ценопопуляций данных видов было проведено в сравнительно-георафическом аспекте с учетом приуроченности к определенным типам растительных сообществ. Установлено наличие полночлeнных возрастных спектров, представленых прегенеративными, генеративными и сенильными растениями с преобладанием молодых вегетирующих особей. Преобладающим типом самоподдержания дягиля и цикория является семенное, а вязеля – вегетативное размножение.
Отмечено, что возрастные спектры ценопопуляций выбранных видов имеют адаптивный хаpaктер, заметно меняются в зависимости от условий внешней среды и антропогенного воздействия и отражают флуктуационный хаpaктер динамических процессов в фитоценозах.
...
06 04 2026 22:57:25
Статья в формате PDF
110 KB...
05 04 2026 0:36:37
Статья в формате PDF
129 KB...
04 04 2026 5:36:29
Статья в формате PDF
111 KB...
03 04 2026 17:16:18
Статья в формате PDF
127 KB...
02 04 2026 15:30:11
Статья в формате PDF
111 KB...
01 04 2026 22:22:26
Статья в формате PDF
276 KB...
30 03 2026 23:13:32
Статья в формате PDF
241 KB...
29 03 2026 15:19:14
Статья в формате PDF
266 KB...
28 03 2026 7:44:12
Статья в формате PDF
141 KB...
27 03 2026 19:28:22
Статья в формате PDF
552 KB...
26 03 2026 22:15:48
В работе дана экологическая оценка возможных последствий на окружающую среду, жизнь и здоровье населения. Показано, что при решении выбора источника энергии необходимо учитывать не только экономические, но и экологические последствия возможного влияния объектов энергетики при строительстве и эксплуатации. Комбинированное производство энергии двух видов на мини–ТЭЦ способствуют гораздо более экологичному использованию топлива по сравнению с раздельной выработкой электроэнергии и тепловой энергии на котельных установках, но и повышению чистоты воздушного бассейна, улучшению общего экологического состояния окружающей среды. Интенсивное шумовое воздействие на организм человека нeблагоприятно влияет на протекание нервных процессов, способствует развитию утомления, изменениям в сердечно-сосудистой системе и появлению шумовой патологии, среди многообразных проявлений которой ведущим клиническим признаком является медленно прогрессирующее снижение слуха. В работе рассмотрено воздействие мини-ТЭЦ с дизельными и газопоршнеыми двигателями мощностью 1000 кВт на окружающую среду. Установлено что, шум, создаваемый электростанцией, состоящей их 4 газопоршневых двигателей мощностью 1000 кВт, будет ниже допустимого для территории, непосредственно прилегающей к жилым домам. Поэтому специальных мероприятий по снижению шума не требуется.
...
25 03 2026 18: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 ::