Автор: Пользователь скрыл имя, 27 Ноября 2011 в 22:20, курсовая работа
Во второй половине 20 столетия мир вступил в эпоху вычислительной математики. На сцену вышла новая наука – кибернетика, которая сразу обрела бытовую практику и в ее силу уверовали практически все.
Джей Форрестер, один из крупнейших специалистов в области теории управления, был профессором в Школе управления Альфреда П. Слоуна в Массачусетсом технологическом институте (МТИ). С 1939 г. до конца Второй Мировой войны он занимался разработкой сервомеханизмов в МТИ, а позднее — цифровых ЭВМ.
ВВЕДЕНИЕ………………………………………………………………………….……………….2
1 СТРУКТУРА МИРОВОЙ СИСТЕМЫ………………………………………………………...…4
1.1 ПОЛНАЯ СХЕМА МИРОВОЙ СИСТЕМЫ…………………………………………………...3
2 МИРОВАЯ МОДЕЛЬ: СТРУКТУРА И ПРЕДПОЛОЖЕНИЯ………………………………....6
2.1 ФОРМАЛИЗАЦИЯ МОДЕЛИ ДЖ.ФОРРЕСТЕРА ……………………………………….…..7
2.2 СТОХАСТИЧЕСКАЯ МОДЕЛЬ………...……………………………………………………..11
3 ПОСТАНОВКА ЗАДАЧИ ОПТИМАЛЬНОГО УПРАВЛЕНИЯ………………………….…...12
3.1 ЧИСЛЕННОЕ ИССЛЕДОВАНИЕ МОДЕЛИ МИРОВОЙ ДИНАМИКИ ФОРРЕСТЕРА………………………………………………………………………………………13
ВЫВОДЫ……………………………………………………………………………………………27
СПИСОК ИСПОЛЬЗУЕМОЙ ЛИТЕРАТУРЫ…………………………………………………...29
[4]
В качестве цели управления можно выбрать стабилизацию на каком-либо заданном уровне одной или нескольких переменных (переменных состояния, переменных-темпов или переменных-параметров). В этом случае можно выбрать квадратичный критерий управления. Синтез оптимального управления нелинейной системой, тем более при неаналитическом задании нелинейностей, в общем случае является крайне сложным делом даже при таком простом критерии оптимальности, как квадратичный. Однако для линеаризованной подходящим образом системы в отсутствие ограничений на управление найти структуру оптимального управления не составляет труда, поскольку решение задачи оптимального управления линейной системой.
3.1 .ЧИСЛЕННОЕ ИССЛЕДОВАНИЕ
МОДЕЛИ МИРОВОЙ ДИНАМИКИ ФОРРЕСТЕРА
Прежде всего, проверим, выполняется ли для модели Форрестера предположение
[ 3]
Где матрица A(k, x(k)) в любой момент времени до и после текущего момента.
вычислим матрицы A(k, x(k)) на интервале времени с t(1) = 1900 г. по t(N) = 2100 г. с шагом h = 1/ 4 года. на неуправляемой стохастической траектории. При этом случайные возмущения оставим только в темпах рождаемости–смертности. Матрицы A(k, x(k)) найдем численным дифференцированием. Вычислим эмпирические средние A и эмпирические дисперсии S2 всех элементов этой матрицы на выбранной траектории:
Для некоторой конкретной реализации траектории неуправляемого движения стохастической динамической системы, описываемой моделью Форрестера, получаем результат, представленный в табл. 1 и 2
Таблица 1
Ai, j | j = 1 | 2 | 3 | 4 | 5 |
i = 1 | 9,473
E–001 |
4,005
E–003 |
5,412
E–002 |
3,071
E–006 |
7,787
E+007 |
2 | 1,250
E–002 |
4,527
E–001 |
2,375
E–001 |
0 | 0 |
3 | 1,250
E–003 |
0 | 9,988
E–001 |
2,848
E–006 |
–2,540
E+007 |
4 | –4,261
E–015 |
0 | –1,112
E–001 |
9,999
E–001 |
5,644
E+008 |
5 | 4,405
E–013 |
7,783
E–015 |
–2,766
E–013 |
1,536
E–016 |
9,805
E–001 |
Таблица 2
Si, j | j = 1 | 2 | 3 | 4 | 5 |
i = 1 | 2,8555
E–002 |
1,2276
E–003 |
9,9051
E–003 |
2,7107
E–006 |
1,9208
E+007 |
2 | 1,0290
E–012 |
5,3336
E–002 |
2,9151
E–012 |
0 | 0 |
3 | 4,2597
E–014 |
0 | 1,6008
E–003 |
2,4816
E–006 |
1,0204
E+007 |
4 | 5,7536
E–013 |
0 | 3,5573
E–002 |
5,5147
E–005 |
2,2675
E+008 |
5 | 8,7312
E–014 |
3,1843
E–016 |
1,8775
E–013 |
1,1860
E–016 |
3,2152
E–004 |
Эти же данные
приведены в логарифмическом
масштабе на рис.1 для абсолютного
значения матрицы A , вытянутой по столбцам,
т.е. представленной в виде массива-столбца,
составленного путем последовательного
присоединения к первому столбцу остальных
ее столбцов. Для наглядности соседние
по номерам значения массива соединены
линиями. с матричными
индексами (3,2), (4,2), (2,4), (2,5), имеющие
соответственно номера 8, 9, 17, 22, равны 0, т.е.
в логарифмическом масштабе равны −∞
, поэтому на рис. 1 они не представлены
(пропущены). Эти пропущенные значения
надо понимать просто как равные нулю.
На этом же рисунке приведены матрицы
A m S в том же представлении. Как видно из
сравнения табл. 1 и 2 и представленных
на рис. 1 массивов, элементы матрицы
относительно незначительно флуктуируют
на траектории движения (эволюции) системы.
Это дает основание считать систему слабо
нелинейной и предположение справедливым.[3]
Изменчивость
матрицы А на траектории
Рисунок 1- Изменчивость матрицы А на траектории
Рассмотрим
в качестве примера реализацию и
детерминированной и
x1(1)=1, 65E+9 человек. Результаты
численного моделирования на двухсотлетнем
интервале времени (с 1900 по 2100 г.) неуправляемой
и управляемой систем мировой динамики
Форрестера в отсутствие и при наличии
стохастических возмущений приведены
на рис. 2. естественно, в силу ограниченности
места мы приводим ход только некоторых
переменных модели.
Рисунок2-Субоптимальное
ограниченное управление
На рис.
2 приведен ход субоптимальных ограниченных
управлений с обратной связью для
детерминированной и
Рисунок 3- Динамика численности населения
На приведенных
здесь же кривых роста численности
населения в неуправляемой
Рисунок 4- Относительная
величина фондов в сельском хозяйстве
На рис. 5 и 6 приведены темпы рождаемости и смертности населения. Видно, что в управляемой системе вместе со стабилизацией численности населения стабилизируются
и темпы рождаемости и смертности, даже при наличии весьма значительных собственных (естественных) флуктуаций этих темпов. В неуправляемой системе темпы рождаемости и
смертности растут
с ростом численности населения.
[5]
Рисунок 5- Темп
рождаемости
Рисунок 6-Темп смертности
Обратим внимание, что колебания темпов рождаемости и смертности в стохастической системе существенно различаются, хотя темп прироста населения равняется просто разности темпов рождаемости и смертности. Здесь проявляются достаточно сложные влияния на темп смертности материального уровня жизни, уровня питания, уровня загрязнения, плотности населения, учитываемые моделью Форрестера. На рис.7 приведен ход безразмерного относительного уровня питания населения. Этот показатель определяет количество пищи на душу населения в единицах среднего мирового уровня питания на душу населения в каком-нибудь выбранном году (в данном случае в неуправляемой модели это соответствует уровню питания населения в 1940 г.). Видно, что в начальном периоде управляемого процесса благодаря интенсивным капиталовложениям в сельское хозяйство уровень питания населения быстро растет, но с ростом населения уровень питания неизбежно начинает падать, причем падение его в управляемой системе происходит быстрее, чем в неуправляемой (требование стабилизации уровня населения). По достижении некоторого минимального для системы значения уровень питания стабилизируется. [5]
Рисунок 7- Относительный
уровень питания
Очевидно, перераспределение капиталовложений в сельское хозяйство должно привести к уменьшению доли фондов в другие сферы человеческой деятельности, в том числе в материальную сферу. Поэтому материальный уровень жизни населения должен снизиться, тем более что капиталовложения в материальную сферу резко снижаются в начальный период интервала управления, затрудняя поддержание и развитие материальной сферы. На рис.8 приведены кривые хода материального уровня жизни при наличии и отсутствии управления. Материальный уровень жизни в модели Форрестера есть безразмерная величина, которая описывает степень изменения эффективности относительной величины фондов на душу населения в сравнении с некоторым стандартом – фиксированным значением.