- •Нелінійні моделі та аналіз складних систем Навчальний посібник
- •Глава 1. Модель, моделювання, прогнозування і управління 20
- •Глава 2. Елементи теорії систем і системний аналіз. Синергетичний підхід 66
- •Глава 3. Рівняння і аналіз складних систем 115
- •Глава 4. Елементи теорії стійкості систем 167
- •Глава 5. Моделювання і аналіз соціально-економічних систем 195
- •Глава 6. Синергетичне моделювання і управління складними системами 234
- •Глава 7. Елементи теорії хаосу і хаотичної динаміки. Фрактали 260
- •Глава 1. Модель, моделювання, прогнозування і управління 20
- •Глава 1. Модель, моделювання, прогнозування і управління 20
- •Глава 2. Елементи теорії систем і системний аналіз. Синергетичний підхід 66
- •Глава 2. Елементи теорії систем і системний аналіз. Синергетичний підхід 66
- •Глава 3. Рівняння і аналіз складних систем 115
- •Глава 3. Рівняння і аналіз складних систем 115
- •Глава 4. Елементи теорії стійкості систем 167
- •Глава 4. Елементи теорії стійкості систем 167
- •Глава 5. Моделювання і аналіз соціально-економічних систем 195
- •Глава 5. Моделювання і аналіз соціально-економічних систем 195
- •Глава 6. Синергетичне моделювання і управління складними системами 234
- •Глава 6. Синергетичне моделювання і управління складними системами 234
- •Глава 7. Елементи теорії хаосу і хаотичної динаміки. Фрактали 260
- •Глава 7. Елементи теорії хаосу і хаотичної динаміки. Фрактали 260
- •Передмова
- •Розділ і моделювання та системний аналіз динамічних процесів
- •Глава 1. Модель, моделювання, прогнозування і управління
- •1.1. Поняття моделей і моделювання
- •1.2. Класифікація засобів моделювання
- •1.3. Поняття економічної системи і принципи її моделювання
- •1.4. Етапи економіко-математичного моделювання
- •1.5. Основні принципи опису виробничо-технологічного рівня економічних систем
- •1.6. Загальний вид математичних моделей та основні напрямки їхнього аналізу
- •1.7. Класифікація економіко-математичних моделей
- •1.8. Моделювання еколого-економічного управління виробничою системою в умовах нестабільності
- •1.9. Деякі принципи моделювання складних систем
- •1.10. Новий підхід до прогнозування поведінки складних систем і катастрофічних процесів (русла і джокери)
- •1.11. Моделювання та управління ризиком
- •Питання для самоперевірки
- •Глава 2. Елементи теорії систем і системний аналіз. Синергетичний підхід
- •2.1. Історія розвитку теорії систем і системного аналізу
- •2.3. Наукові напрямки системного аналізу
- •2.4. Система, її структура і функціонування
- •Типи системних задач залежно від ситуації
- •2.5. Синергетичний підхід в аналізі складних систем
- •Розділ іі математичні основи нелінійної динаміки і аналізу складних систем
- •Глава 3. Рівняння і аналіз складних систем
- •3.1. Моделі і аналіз нелінійних динамічних систем
- •3.1.1. Системи звичайних диференційних рівнянь
- •3.1.1.1. Основні визначення і теореми
- •3.1.1.2. Особливі точки та їхнє інваріантне різноманіття
- •Стійкість визначається нелінійними членами
- •3.1.1.3. Періодичні та неперіодичні розв’язки, граничні цикли та інваріантні тори
- •3.1.2. Аналіз нелінійної системи з дискретним часом
- •3.1.3. Використання теорії різницевих рівнянь для моделювання процесу мобілізації
- •3.1.4. Основи теорії диференційних рівнянь
- •3.2. Нелінійні моделі складних систем з хаотичною динамікою (стислий огляд)
- •Глава 4. Елементи теорії стійкості систем
- •4.1. Аналіз нелінійних економічних систем, що розвиваються
- •4.1.1. Основні дослідження стійкості нелінійних динамічних систем
- •4.1.2. Якісний аналіз економічної системи, що знаходиться під впливом новітніх інформаційних технологій (ніт)
- •4.2. Елементи теорії структурної динаміки
- •4.2.1. Основи теорії катастроф
- •Збурення Збурення Збурення Збурення Збурення Збурення Збурення Збурення збурення збурення
- •4.2.2. Біфуркації на фазовій площині
- •4.2.3. Дисипативні структури і. Пригожина
- •Розділ ііі моделі та аналіз складних процесів і систем
- •Глава 5. Моделювання і аналіз соціально-економічних систем
- •5.1. Деякі базові математичні моделі та їхнє застосування в дослідженні соціально-економічних процесів
- •5.1.1. Модель Мальтуса
- •5.1.2. Логістична модель
- •5.1.3. Експоненційна модель з відловом
- •5.1.4. Логістична модель з відловом
- •5.1.5. М’яка логістична модель з відловом
- •5.1.6. Модель Лотки–Вольтерра
- •5.2. Приклад моделювання і аналізу соціально-економічних процесів
- •5.2.1. Стійкість ринкових механізмів
- •5.2.2. Народження хаосу
- •5.3. Елементи теорії м’якого моделювання
- •5.3.1. Модель війни або битви
- •5.3.2. Оптимізація як шлях до катастрофи
- •5.3.3. Жорсткі моделі як шлях до помилкових передбачень
- •5.3.4. Небезпека багатоступеневого управління
- •5.3.5. Математичні моделі «перебудови» в срср
- •5.3.6. Математика і математична освіта у сучасному світі
- •5.4. Моделі нелінійної економічної динаміки з урахуванням швидкості встановлення змінних
- •5.4.1. Окремі моделі нелінійної економічної динаміки
- •5.4.2. Узагальнена модель динаміки економіки
- •Глава 6. Синергетичне моделювання і управління складними системами
- •6.1. М’яке нелінійне управління: синергетичний підхід в управлінні
- •6.2. Глобальне моделювання і аналіз світової динаміки
- •6.2.1. Модель світової динаміки Форрестера
- •6.2.1. Глобальна модель динаміки Месаровіча–Пестеля (м–п-модель)
- •6.2.2. Феноменологічна макромодель світової динаміки і стійкого розвитку
- •Глава 7. Елементи теорії хаосу і хаотичної динаміки. Фрактали
- •7.1. Теорія динамічного хаосу та її застосування
- •7.1.1. Динамічний хаос
- •7.1.2. Економічний хаос у детермінованих системах
- •7.1.2.1. Хаос у детермінованих системах
- •7.1.2.2. Економічний хаос у дискретній системі
- •7.1.2.3. Аперіодичне оптимальне економічне зростання
- •7.1.2.4. Динаміка міст – система Лоренца (приклад застосування)
- •7.1.2.5. Хаос у моделі міжнародної економіки
- •7.1.2.6. Хаос і економічне прогнозування
- •7.1.2.7. Деякі критерії класифікації атракторів
- •7.1.3. Дивні атрактори
- •7.1.4. Динамічний хаос і обмеження області прогнозу
4.1.2. Якісний аналіз економічної системи, що знаходиться під впливом новітніх інформаційних технологій (ніт)
Для прикладу візьмемо двовимірну макроекономічну структуру, яка добре моделює взаємодію двох суміжних галузей завдяки НІТ (газова промисловість, чорна металургія і т. д.).
Як відомо, рівняння, що описують таку макроекономічну систему, будуть мати такий вигляд:
(2.66)
де
Будемо
вважати, що коефіцієнти
мало змінюються в часі, і їх можна
прийняти на інтервалі спостереження
постійними. Тоді можна виділити праві
частини:
(2.67)
Знайдемо
стан рівноваги. Прирівнюючи ліві частини
рівняння (2.67) до нуля
одержимо
Очевидно, загальна кількість станів рівноваги буде дорівнювати комбінації отриманих станів, тобто їхня загальна кількість дорівнювати 4:
1)
2)
3)
4)
Обчислимо, у загальному виді, матрицю лінійного наближення:
Для складання дисперсійного (характеристичного) рівняння маємо:
звідки
дисперсійне (характеристичне) рівняння
отримує такий вигляд:
(2.68)
де
1.
Розглянемо перший з чотирьох станів
рівноваги
У
цьому випадку маємо
(2.69)
Оскільки
то ми отримаємо корені дійсні, одного
знака, і отже, стан
є нестійким вузлом. Фізично це означає,
що система не може знаходитися при
нульовому випуску, і за найменших змін
параметрів обов’язково спостерігається
інтенсивне зростання випуску як у
першій, так і в другій макросистемах.
2. Розглянемо
другий стан рівноваги:
Для цього випадку маємо
тому дисперсійне рівняння набуде
такого вигляду:
(2.70)
Тут
маємо критичну точку
в якій корені стають нульовими. При
система буде стійкою, стан рівноваги є
стійким вузлом. При
стан рівноваги є нестійким вузлом.
Очевидно, критична точка
– це стан нейтральної стійкості.
3. Тепер розглянемо третій стан рівноваги. У цьому випадку
Дисперсійне рівняння набуде такого вигляду:
(2.71)
(2.72)
Тут
маємо критичну точку
в якій корені стають нульовими і
система в ній нейтрально стійка. При
система буде стійкою, стан рівноваги –
стійким вузлом. При
стан рівноваги є нестійким вузлом.
4. Розглянемо четвертий стан рівноваги:
У цьому випадку
Дисперсійне рівняння буде виглядати так:
(2.73)
загальний розв’язок:
(2.74)
Якщо підкореневий вираз від’ємний, то стан рівноваги буде стійким вузлом, при зміні знака підкореневого виразу стан рівноваги перетвориться в стійкий фокус. У критичній точці, коли
система стрімко переходить зі стійкого фокуса в стійкий вузол.
4.2. Елементи теорії структурної динаміки
4.2.1. Основи теорії катастроф
На початку 70-х рр. став популярним термін «катастрофа», що позначає стрибкоподібні зміни, які виникають при плавних змінах значень параметрів. У популярних виданнях теорія катастроф рекламувалась як переворот у математиці, порівняно з винаходом диференційного й інтегрального обчислення. За останні двадцять років з’явилися сотні наукових публікацій з природознавства і техніки, в яких теорія катастроф успішно застосовувалась. Оприлюднено також праці про використання моделей теорії катастроф у економіці, психології, лінгвістиці, соціології.
Після періоду ейфорії, викликаною широкою саморекламою, з’явилися більш тверезі оцінки застосування теорії катастроф, більш того, з’ясувалося, що багато серйозних результатів було отримано до «проголошення» теорії катастроф.
Один з провідних радянських математиків, В. І. Арнольд, відзначає, що обґрунтованість використання теорії катастроф істотно залежить від обґрунтованості початкових передумов: «Наприклад, у теорії ударів пружних конструкцій і в теорії перекидання кораблів прогнози теорії цілком підтверджуються експериментом. З іншого боку, в біології, психології та соціальних науках (скажімо, в додатках до теорії поведінки біржових гравців або до вивчення нервових хвороб) як початкові передумови, так і висновки мають швидше евристичне значення».
Розглянемо основні положення теорії катастроф на прикладі катастрофи «зборка», якій відповідає диференційне рівняння
(2.76)
При
варіюванні значень параметрів
і
поведінка системи (кількість
стаціонарних точок, їхнє розташування)
буде також змінюватися. Для вивчення
якісного характеру цих змін розглянемо
потенційну функцію
(2.77)
Помітимо,
що
На рис. 2.39 наведено двовимірні
графіки, які характеризують поведінку
функції
На
рис. 2.39(а) зображена крива, що називається
біфуркаційною:
Ця
крива розділяє площину (
)
на дві частини. Всередині кривої
функція
має два мінімуми (рис. 2.39(b)). За межами
цієї кривої функція
має тільки один мінімум (рис. 2.39(с)). Як
відомо, екстремальні значення функції
можна визначити, прирівнявши до
нуля першу похідну:
(2.78)
але
це ж рівняння визначає стаціонарні
точки диференційного рівняння
(2.76). Отже, крива
розділяє площину (
)
на дві області, у першій з яких
є одна особлива точка, а у другій – дві.
Рис. 2.39. Графіки потенційної функції
Якщо
в тривимірному просторі на вертикальній
осі відкласти положення особливих
точок рівняння (2.76) х
(теорія
катастроф вивчає питання структурної
стійкості спеціального класу диференційних
рівнянь, права частина яких може бути
представлена у вигляді градієнтної
системи, тобто як рух у полі потенційних
сил з потенціалом
:
),
а на двох інших осях – значення
параметрів
і
,
то одержимо поверхню катастроф
(рис. 2.40). Проекція на площину параметрів
(
)
точок поверхні, в яких є вертикальна
дотична, є біфуркаційною кривою.
Припустимо,
що безперервній зміні значень параметрів
і
на рис. 2.40 відповідає рух по кривій
,
тоді в точці
відбувається
катастрофа – система стрибком переходить
з одного стану рівноваги в інший (
).
Відзначимо, що кожному значенню параметрів і , що знаходяться всередині біфуркаційної кривої, відповідають два різні стани системи (бімодальність).
На
поверхні катастроф можна спостерігати
явище гістерезису, коли різка зміна в
поведінці системи істотно залежить від
передісторії процесу. Наприклад, при
зміні стану системи вздовж кривої
,
як було наведено вище, відбувається
стрибок з верхнього аркуша на нижній,
з точки
в точку
.
Але при русі вздовж кривої
стрибок з нижнього аркуша на верхній
відбудеться не в точці
,
а в точці
.
Параметр одержав назву нормального фактора, параметр – розщеплюючого фактора.
Монографія Постона і Стюарта була присвячена використанню теорії катастроф у проблемах соціологічного моделювання, зокрема в ній аналізуються дослідження порушень режиму у в’язниці Гартрі протягом 1972 р. Застосовуючи факторний аналіз, автори дослідження розділили фактори, що впливають на хаос, на дві загалом незалежні групи: напруженість (почуття розчарування і безвиході, тяжкий стан) та роз’єднаність (взаємне відчуження, відсутність спілкування, розподіл на два табори).
тихо
тихо
Рис. 2.40. Катастрофа «зборка» |
Рис. 2.41. Модель збурень у в’язниці |
збурення
збурення
збурення
збурення
збурення
збурення
збурення
збурення
збурення
збурення
збурення
збурення
збурення
збурення
збурення
збурення
збурення
збурення
збурення
збурення
збурення
збурення
тихо
тихо
тихо
тихо
тихо
тихо
тихо
тихо
тихо
А
– теорія
Аналіз
показав, що зі зростанням напруженості
підвищується імовірність збурень,
а збільшення роз’єднаності пов’язане
з характером з
збурення
збурення
збурення
збурення
збурення
Динаміка системи може бути проаналізована за допомогою катастрофи типу «зборки», з віссю напруженості, ортогональної до вістря, і віссю роз’єднаності, що йде вздовж вістря (рис. 2.41).
З рис. 2.41 видно, що при низьких значеннях роз’єднаності система прагне до стійкого стану помірного хвилювання, а при високому рівні роз’єднаності – змінює свій стан стрибком з нижнього аркуша на верхній і назад. Система робить коливання всередині біфуркаційної множини катастрофи «зборка», що підтверджується реальними даними.
В. І. Арнольд вважає, що застосування теорії катастроф за умов, коли не тільки не відомі досліджувані функції, але й саме їхнє існування є проблематичним, має спекулятивний характер. Яскравий приклад такого твердження належить англійському математикові К. Зіману, який запропонував модель творчої особистості.
Вважатимемо, що творчу особистість (вченого) можна описати трьома параметрами: «техніка», «захопленість», «досягнення». Припустимо, що катастрофа «зборки», зображена на рис. 2.42, задовільно описує досягнення вчених, умовно розділених на дві групи: «генії» та «маніяки».
Дослідимо,
як залежать досягнення вченого від його
техніки і захопленості. Якщо
захопленість (
)
мала, то досягнення (
)
монотонно зростають разом зі
зростанням техніки (
).
Біфуркації відбуваються при досить
великій захопленості. Наприклад, якщо
техніка і захопленість змінюються
вздовж кривої 1 на рис. 2.42, тоді в точці
досягнення зростають стрибком. При
цьому ми переходимо в область високих
досягнень (генії).
Однак
якщо зростання захопленості не
супроводжується ростом техніки, то рух
відбувається по кривій 3 і досягнення
стрибком потрапляють у точку
.
Ми опиняємося в області малих досягнень
(«маніяки»). Зауважимо, що «маніяк»
і «геній» у цій моделі можуть мати рівні
техніку і захопленість, але при цьому
їхні досягнення різко відрізняються –
внаслідок того, що стрибки зі стану
«геній» у стан «маніяк» відбуваються
на різних лініях.
Моделі теорії катастроф є математичною інтерпретацією діалектики Гегеля – законів переходу кількості в якість, єдності та боротьби протилежностей.
