Моделирование систем |
||
4.1. Основные понятия моделирования
Модель - объект или описание объекта, системы для замещения (при определенных условиях предложениях, гипотезах) одной системы (т.е. оригинала) другой системы для изучения оригинала или воспроизведения его каких - либо свойств. Модель - результат отображения одной структуры на другую. Отображая физическую систему (объект) на математическую систему (например, математический аппарат уравнений) получим физико - математическую модель системы или математическую модель физической системы. В частности, физиологическая система - система кровообращения человека, подчиняется некоторым законам термодинамики и описав эту систему на физическом (термодинамическом) языке получим физическую, термодинамическую модель физиологической системы. Если записать эти законы на математическом языке, например, выписать соответствующие термодинамические уравнения, то получим математическую модель системы кровообращения. Эту модель можно назвать физиолого - физико - математической моделью или физико - математической моделью. Модели, если отвлечься от областей, сфер их применения, бывают трех типов: познавательные, прагматические и инструментальные. Познавательная модель - форма организации и представления знаний, средство соединение новых и старых знаний. Познавательная модель, как правило, подгоняется под реальность и является теоретической моделью. Прагматическая модель - средство организации практических действий, рабочего представления целей системы для ее управления. Реальность в них подгоняется под некоторую прагматическую модель. Это, как правило, прикладные модели. Инструментальная модель - является средством построения, исследования и/или использования прагматических и/или познавательных моделей. Познавательные отражают существующие, а прагматические - хоть и не существующие, но желаемые и, возможно, исполнимые отношения и связи. По уровню, "глубине" моделирования модели бывают эмпирические - на основе эмпирических фактов, зависимостей, теоретические - на основе математических описаний и смешанные, полуэмпирические - использующие эмпирические зависимости и математические описания.
Математическая модель М описывающая ситему
S (x Основные требования к модели: наглядность построения; обозримость основных его свойств и отношений; доступность ее для исследования или воспроизведения; простота исследования, воспроизведения; сохранение информации, содержавшиеся в оригинале (с точностью рассматриваемых при построении модели гипотез) и получение новой информации. Проблема моделирования состоит из трех задач:
Модель М называется статической, если среди xi нет временного параметра t. Статическая модель в каждый момент времени дает лишь "фотографию" сиcтемы, ее срез. Модель - динамическая, если среди xi есть временной параметр, т.е. она отображает систему (процессы в системе) во времени. Модель - дискретная, если она описывает поведение системы только в дискретные моменты времени. Модель - непрерывная, если она описывает поведение системы для всех моментов времени из некоторого промежутка времени. Модель - имитационная, если она предназначена для испытания или изучения, проигрывания возможных путей развития и поведения объекта путем варьирования некоторых или всех параметров xi модели М. Модель - детерминированная, если каждому входному набору параметров соответствует вполне определенный и однозначно определяемый набор выходных параметров; в противном случае - модель недетерминированная, стохастическая (вероятностная). Можно говорить о различных режимах использования моделей - об имитационном режиме, о стохастическом режиме и т. д.
Модель включает в себя: объект О, субъект (не обязательный) А, задачу
Z, ресурсы B, среду моделирования С: М= Свойства любой модели таковы:
Жизненный цикл моделируемой системы:
Основными операциями используемыми над моделями являются:
Модели и моделирование применяются по следующим основным и важным направлениям.
В базовой четверке информатики: "модель - алгоритм - компьютер - технология" при компьютерном моделировании главную роль играют уже алгоритм (программа), компьютер и технология (точнее, инструментальные системы для компьютера, компьютерные технологии). Например, при имитационном моделировании (при отсутствии строгого и формально записанного алгоритма) главную роль играют технология и средства моделирования; аналогично и в когнитивной графике. Основные функции компьютера при моделировании систем:
Компьютерное моделирование - основа представления знаний в ЭВМ (построения различных баз знаний). Компьютерное моделирование для рождения новой информации использует любую информацию, которую можно актуализировать с помощью ЭВМ. Разновидностью компьютерного моделирования является вычислительный эксперимент. Компьютерное моделирование, вычислительный эксперимент становится новым инструментом, методом научного познания, новой технологией также из-за возрастающей необходимости перехода от исследования линейных математических моделей систем . Компьютерное моделирование, от постановки задачи - до получения результатов, проходит следующие этапы.
4.2. Модели и моделирование систем Математическое и компьютерное моделирование. Данный период характерен необходимостью моделирования различных социально- экономических процессов и систем и принятия решений на основе результатов моделирования, ибо такие системы и процессы достаточно сложны, многогранны, динамичны, подвержены случайным воздействиям. Достаточно интенсивно моделируются сейчас такие социально-экономические процессы, как демографические (например, эволюция и цикличность), социальные (например, поведение социальных групп, социальных последствий тех или иных решений), экономические (например, рыночные отношения, налоговые сборы, риски), гуманитарные (например, воздействие на человека информационного потока) и др. Приведём пример математического моделирования некоторой системы (полный жизненный цикл моделирования).Решение поставленной задачи разобьем на этапы, в соответствии с этапами жизненного цикла моделирования, объединяя для удобства некоторые этапы для удобства и краткости изложения. Этап 1. Содержательная постановка задачи Современное производство характеризуется тем, что некоторая часть производимой продукции (в стоимостном выражении) возвращается в виде инвестиций (т.е. части конечной продукции, используемой для создания основных фондов производства) в производство. При этом время возврата, ввода в оборот новых фондов может быть различной для различного рода производства. Необходимо промоделировать данную эту ситуацию и выявить динамику изменения величины основных фондов производства (капитала). Сложность и многообразие, слабая структурированность и плохая формализуемость основных экономических механизмов, определяющих работу предприятий не позволяют преобразовать процедуры принятия решений в экономической системе в полностью эффективные математические модели и алгоритмы прогнозирования. Поэтому часто эффективно использование простых, но гибких и надёжных процедур принятия решения. Рассмотрим одну такую простую модель. Эта модель будет полезна для прогноза событий и связанных с ними социально-экономических процессов. Этап 2. Формулировка гипотез, построение, исследование модели Структура производства и сбыта часто зависит от изменений в окружающей среде (социально-экономических условий). Динамика изменения величины капитала определяется, в основном, в нашей модели, простыми процессами производства и описывается так называемыми обобщенными коэффициентами амортизации (расхода фондов) и потока инвестиций (часть конечного продукта, используемого в единицу времени для создания основных фондов). Эти коэффициенты - относительные величины (за единицу времени). Необходимо разработать и исследовать модель динамики основных фондов. Считаем при этом допустимость определённых гипотез, определяющих рассматриваемую систему производства. Пусть x(t) - величина основных фондов (капитала) в момент времени t, где 0 £ t £ N. Через промежуток времени Dt она будет равна x(t+Dt). Абсолютный прирост равен Относительный прирост будет равен Примем следующие гипотезы:
Считая Dt ® 0, а также учитывая определение производной, получим из предыдущего соотношения следующее математическое выражение закона изменения величины капитала - математическую модель (уравнение) динамики капитала (такие уравнения называются дифференциальными): где х(0) - начальное значение капитала в момент времени t=0. Эта величина х0 везде в дальнейшем будет считаться заданной. Эта простейшая модель динамики величины капитала. Эта простейшая модель не отражает того факта, что социально- экономические ресурсы производства таковы, что между выделением инвестиций и их введением и использованием в выпуске новой продукции проходит некоторое время - лаг. Учитывая это можно записать модель (1) в виде: Данной непрерывной, дифференциальной, динамической модели можно поставить в соответствие простую дискретную модель: где n - предельное значение момента времени при моделировании. Эта дискретная модель получается из непрерывной при Dt=1, а также заменой производной x'(t) на относительное приращение Dt (замена, как это следует из определения производной, справедлива при малых Dt). Этап 3. Построение алгоритма и программы моделирования Рассмотрим для простоты режим моделирования когда m, c, y - известны и постоянны, а также рассмотрим наиболее простой алгоритм моделирования в укрупнённых шагах.
Приведём программу на Паскале для имитационного моделирования (компактность программы «принесена в жертву» структурированности).
PROGRAM MODFOND; {Исходные данные находятся в файле in.dat текущего каталога} {Результаты записываются в файл out.dat текущего каталога} Uses Crt, Graph, Textwin; Type Vector = Array[0..2000] of Real; Mas = Array[0..2000] of LongInt; Var Time,Lag,t,dv,mv,i,yi,p : Integer; tmax,tmin : LongInt; a,b,m,X0,maxx,minx,aa,bb,cc,sx,tk : Real; x : Vector; ax,ay : Mas; ch : Char; f1,f2 : Text; {----------------------------------------------------------------------------} Procedure InputKeyboard; { Ввод с клавиатуры } Begin OpenWindow(10,5,70,20,' Ввод данных ',14,4); ClrScr; WriteLn; WriteLn('Введите время Т прогнозирования системы:'); Repeat Writeln('Для удобства построения графика введите Т не меньше 2'); Write('Т ='); ReadLn(Time); until Time>=2; WriteLn('Введите лаг:'); Repeat Write('Лаг должен быть строго меньше Т - '); ReadLn(Lag); until Lag Этап 4. Проведение вычислительных экспериментов Эксперимент 1. Поток инвестиций постоянный и в каждый момент времени равен 111. В начальный момент капитал - 1000 руб. Коэффициент амортизации - 0.0025. Построить модель динамики (посуточно) и найти величину основных фондов через 50 суток, если лаг равен 10 суток. Эксперимент 2. Основные фонды в момент времени t=0 была равны 50000. Через какое время общая их сумма превысит 1200000 руб., если поток инвестиций постоянный и равен 200, а m=0.02, T=5? Этап 5. Модификация (развитие) модели
Для модели динамики фондов с переменным законом потока инвестиций:
Компьютерное моделирование и вычислительные эксперименты Рассмотрим следующую модель течения эпидемии. Пусть существует группа из n контактирующих индивидуумов, в которой в момент времени t имеется x восприимчивых индивидуумов, y - источников инфекции и z - удаленных (изолированных или выздоровевших и ставших невосприимчивыми к инфекции) индивидуумов. Таким образом, имеем: x+y+z=n. Определим исходные параметры системы: b — частота контактов между членами группы; g — частота случаев удаления; m — скорость пополнения восприимчивых индивидуумов извне (она же и есть скорость гибели индивидуумов, удаленных из популяции); y* — критическое значение, при котором начинается эпидемия. После какой — либо одной вспышки эпидемии, в результате которой плотность восприимчивых индивидуумов упадет ниже критического значения y*, наступает период относительного затишья, длящийся до тех пор, пока снова не будет достигнуто критическое значение y* и не возникнет новая вспышка. За время Dt группа восприимчивых индивидуумов, с одной стороны, уменьшается на b xy Dt за счет заражения части из них, а с другой — увеличивается на n Dt. При t=0 заданы x(0)=x0, y(0)=y0, z(0)=z0 - число восприимчивых, инфекционных и удаленных индивидуумов соответственно. Рассматривается промежуток времени t (0 £ t £ T). При малых Dt=1 (например, минута) получим соотношения (дискретную модель): Приведём результаты проведённых вычислительных экспериментов.
Первые два случая демонстрируют периодичность эпидемии. В последнем примере имеет место случай стационарного состояния, т.е. после некоторого момента времени t количества индивидуумов в группах x, y, z остаются постоянными. Программу компьютерного моделирования и вычислительные эксперименты реализовали студентки математического факультета КБГУ Буздова Асият и Мильман Евгения. Компьютерное моделирование. Рассмотрим проблему расчета влажности почвы с учетом накапливаемой биомассы и прогнозирования урожайности сельхозкультур по заданной (экологически обоснованной) влагообеспеченности корнеобитаемого слоя почвы. Разработать соответствующую компьютерную моделирующую среду, которая позволяет решать задачи прогноза влажности корнеобитаемого слоя почвы и урожайности (биомассы) сельхозкультур на заданный момент времени с развитыми интерфейсными средствами, рассчитанными на неподготовленного пользователя-агронома, эколога и др. Опишем одну такую программную среду, реализованную реально в среде Delphi 2.0 Windows 95 автором и студентами КБГУ Кирьязевой С.К. и Кирьязевым Д.А. для расчета влажности почвы и определения урожайности сельхозкультур. Общение с пользователем осуществляется с помощью диалогового окна “Расчет влажности и урожайности”, содержащего 5 страниц: “Эксперимент”, “С/х культура”, “Регион”, “Рабочая” и “Результат”. Страница “С/х культура”- для ввода входной информации по культуре. Страница “Регион” - для ввода информации по региону эксперимента. Страница “Эксперимент” выглядит следующим образом. Данная страница предназначена для ввода названия эксперимента, выбора названия культуры из списка имеющихся (которые были введены на странице “С/х культура”), региона проведения эксперимента из списка регионов, введенных на странице “Регион”, ввода даты посева культуры и даты снятия урожая, для ввода данных по величинам, зависящим от вида культуры, типа почв (по фазам вегетации). Предусмотрена возможность выбора: вычислять ли значения коэффициентов e и l или же рассматривать их как постоянные величины? После заполнения страницы “Эксперимент”, можно произвести расчет влажности почвы и прогноз урожайности культуры. Для этого достаточно лишь нажать кнопку “Произвести расчет” на странице “Эксперимент”. После этого автоматически раскрывается страница “Результат” с таблицей рассчитанных величин и выводится график зависимости влажности почвы от времени (1 - синим цветом), накопления биомассы растений от времени при вычисленной влажности (2 - зеленым цветом) и оптимального развития растений - по экспериментальным данным за прошлый год (3 - красным цветом). Страница “Рабочая” - для визуального анализа расчётных величин. Были проведены вычислительные эксперименты для двух сельхозкультур - кукурузы и ярового ячменя с использованием общедоступных данных (это также можно отнести к достоинствам системы). Данные по температуре воздуха, величине осадков, уровню грунтовых вод и относительной влажности воздуха представлены с интервалом в 10-15 суток за весь период вегетационного цикла растения. Программа отображает результаты расчета в таблице и на графике. График оптимального развития рассматриваемой культуры имеет “ступенчатый” характер ввиду того, что экспериментально полученные значения за прошлый год вводятся по фазам вегетации, а для межфазных периодов программно рассчитываются по соответствующим математическим моделям. Результаты экспериментов приведены ниже. Эксперимент 1
С/х культура: Кукуруза "Луч-300"; Fmax = 20 Дж/(м2сут.); s = 0,6; m = 10-8; а = 0,8; Время проведения посева с 01.04. по 20.04. Тип почвы: Черноземные почвы.
где Р - величина осадков (мм); Н - уровень грунтовых вод (м3/га);
С/х культура: Кукуруза "Луч-300". Тип почвы: Черноземные почвы. Дата посева: 02.04.97. Дата снятия: 10.07.97. e = 0,0370; l = 0,0002. Результаты расчетов - в виде графиков, таблица расчётов не выводится. Экспертная система. Спроектируем одну гипотетическую базу знаний и экспертную систему. Структура базы знаний и экспертной системы изображена ниже. Экспертная система может рассуждать (имитировать рассуждения) и настраиваться на предметную область. Рассматривается база знаний и экспертная система “Социально- экономико - экологическая система”, которая построена (Казиевым В. М. и студентом КБГУ Тебуевым М. Д.) с использованием аппарата нечётких множеств и логики и позволяет оценивать (качественно) социально - экономико - экологическое состояние некоторой среды по задаваемым пользователям (экспертом) количественным оценкам тех или иных параметров среды (выбираемых из базы знаний системы). Для каждого входного фактора в диалоговом режиме задаются относительные (от 0 до 1) оценки влияния этого фактора (вес фактора). После анализа этих данных (этой экологической обстановки) система принимает на основе базы знаний решение о состоянии социально-экологической среды, используя удельную количественную оценку (от 0 до 1) и десятибалльную качественную систему оценок. Функции и работу системы характеризует нижеприведённый сценарий и протокол диалога с этой системой.
Экспертная система
(23.02.1998 - Понедельник, 11: 23: 37) Входные данные:
Принятие решения о социально-экономико-экологической обстановке:
Пакеты прикладных программ (ППП). ППП - комплекс программ, имеющих следующие особенности (отличающие его от "большого программного комплекса"):
Пример интегрированного ППП - простой и универсальный пакет статистического анализа данных SPSS. Интерфейс пользователя с SPSS для Windows реализуется с помощью простых меню и диалоговых окон т.е. как и предыдущая разработка, SPSS свободна от использования специально изучаемого командного языка пакета. Имеется редактор Data Editor для визуального контроля вводимых данных, функционально аналогичный, например, Excel. По столбцам отображаются варьируемые переменные, а по строкам - наборы их вариации, причем с каждой из переменных можно ознакомиться, вызвав её имя. Ввод данных - аналогичен вводу данных, например, в Excel. В диалоговых окнах можно определять и сложных выражений арифметического или логического типа, используемых далее в расчётах. Опишем расчёты с использованием ППП по анализу систем и динамики задолженности, например, множественный регрессионный анализ; результаты одного такого анализа с помощью авторского ППП приведены ниже. Смысл переменных: х(1) - коэффициент абсолютной ликвидности; х(2) - коэффициент текущей ликвидности; х(3) - дебиторская задолженность; х(4) - кредиторская задолженность; х(5) - превышение кредиторской задолженности над дебиторской; х(6) - коэффициент финансовой зависимости; х(7) - коэффициент соотношения привлечённых и собственных средств; х(8) - кредиторская задолженность перед бюджетом; х(9) - кредиторская задолженность по социальному страхованию и внебюджетным платежам; х(10) - коэффициент собираемости налоговых платежей; х(11) - коэффициент собираемости налога на добавленную стоимость. Были проведены различные вычислительные эксперименты, например, если: х(1) - коэффициент абсолютной ликвидности, х(2) - текущей ликвидности, х(6) - финансовой зависимости, х(7) - коэффициент привлечения собственных средств, а y = х(10) - коэффициент собираемости налогов, то находится зависимость коэффициента собираемости налогов от коэффициентов абсолютной ликвидности, текущей ликвидности, финансовой зависимости и привлечения собственных средств, т.е. зависимость вида Типы экспериментов определяются экономическими соображениями, например, с целью выявления факторов, наиболее влияющих на собираемость налогов (на коэффициент собираемости налогов). БАКСАНСКИЙ РАЙОН КБР
(205 предприятия(й)) Регрессионная модель вида Y = a(0) + a(1)*x(1) + ... + a(7)*x(7)
Таблица 1. Таблица коэффициентов модели
Примечание: "+" - коэффициент значим, "-" - коэффициент не значим.
Коэффициент множественной корреляции значим и равен: 0.98. Таблица 2.
Таблица корреляции y и x(i), i=1,2,...,7
Таблица 3.
Таблица адекватности модели
Рис. 27. Фрагменты результатов вычислительных экспериментов. В результате анализа проведенных экспериментов по каждому району, городу найдена регрессионная зависимость с очень высокой степенью адекватности; коэффициент множественной корреляции равен 0.99 - 1.0, а относительная погрешность в среднем порядка 5 - 8 процентов (для таких зависимостей такая погрешность считается очень низкой). При этом:
Вопросы для самоконтроля
| ||
Copyright © 2001, В.М. Казиев Веб-мастер: Артур Балкаров Дизайн: Феликс Джегутанов Разрешается использование в некоммерческом образовании. Учреждениям коммерческого использования - разрешается только с письменного разрешения автора. |