{робочі масиви}: масив з вещ: масив з вещ
Рішення: масив з точка_вещ
{Опис змінних і об'єктів для графічного виведення}
Чісло_Решеній: цілий
Чісло_Собственних_Значеній: цілий
Текущее_Решеніе: цілий
Фазовий_Портрет: масив з ломаная_вещ
Текущая_Ломаная: ломаная_вещ
Устойчівие_Узли: масив з точка_вещ
Неустойчівие_Узли: масив з точка_вещ
Устойчівие_Фокуси: масив з точка_вещ
Неустойчівие_Фокуси: масив з точка_вещ
Центри: масив з точка_вещ
Сідла: масив з точка_вещ
МОДЕЛЬ
{початкові умови} поч=0.0
{параметри обчислювального експерименту}=0.05кон=2
Найті_Все_Решенія (V, 1E - 8, x1_мін, x1_макс, _мін, x2_макс, Шаг_Сеткі, Рішення)
Чісло_Елементов (Рішення) -> (Чісло_Решеній)
цикл Текущее_Решеніе=1 до Чісло_Решеній
початок
Абсциса (Рішення [Текущее_Решеніе]) -> (X1)
Ордината (Рішення [Текущее_Решеніе]) -> (X2)
Собственние_Значенія (J, 300, 1E - 15) -> (Re, Im)
Чісло_Елементов (Re) -> (Чісло_Собственних_Значеній)
якщо Чісло_Собственних_Значеній=2 то
початок
Re_1=Re [1] _2=Re [2]=Re_1 * Re_2
якщо Im [1]=0.0 то {речові СЗ}
якщо R < 0.0 то
Добавіть_Елемент (Рішення [Текущее_Решеніе]) ->
(Сідла)
інакше
якщо Re_1 < 0.0 то
Добавіть_Елемент (Рішення [Текущее_Решеніе]) ->
(Устойчівие_Узли)
інакше
Добавіть_Елемент (Рішення [Текущее_Решеніе]) ->
(Неустойчівие_Узли)
інакше {комплексні СЗ}
якщо Re_1 < 0.0 то
Добавіть_Елемент (Рішення [Текущее_Решеніе]) ->
(Устойчівие_Фокуси)
інакше
якщо Re_1=0.0 то
Добавіть_Елемент (Рішення [Текущее_Решеніе]) ->
(Центри)
інакше
Добавіть_Елемент (Рішення [Текущее_Решеніе]) ->
(Неустойчівие_Фокуси)
кінець
конецs=0
поки x1s < x1_макс
началоs=0
поки x2s < x2_макс
початок
гудок=x1s=x2s
Очістіть_Массів (Текущая_Ломаная)
Добавіть_Точку (x1, x2) -> (Текущая_Ломаная)=tнач
{розрахунковий цикл інтегрування}
поки t < tкон
початок
{розрахунок методом Рунге-Кутта 4-го порядку} n=x1 X2n=x2
{1-й подшаг}
tt=t=h * dX1 du1=h * dX2
{2-й подшаг}=t + h / 2.0=X1n + di1/2.0 x2=X2n + du1/2.0=h * dX1 du2=h * dX2
{3-й подшаг}=t + h / 2.0=X1n + di2/2.0 x2=X2n + du2/2.0=h * dX1 du3=h * dX2
{4-й подшаг}=t + h=X1n + di3 x2=X2n + du3=h * dX1 du4=h * dX2=X1n + (di1 + 2.0 * di2 + 2.0 * di3 + di4) / 6.0= X2n + (du1 + 2.0 * du2 + 2.0 * du3 + du4) / 6.0
t=t + h
Добавіть_Точку (x1, x2) -> (Текущая_Ломаная)
кінець
Добавіть_Елемент (Текущая_Ломаная) ->
(Фазовий_Портрет) s=x2s + Шаг_Сеткі
конецs=x1s + Шаг_Сеткі
кінець
Оновити (ГРАФІКА)
КІНЕЦЬ
Фазовий портрет:
Дана система:
Модель для системи Тополог взята з попереднього прикладу з параметрами: _мін: вещ=- 1.0_макс: вещ=5.5_мін: вещ=- 1.0_макс: вещ=2.5
: функція (x1, x2)=(4 - 2.5 * x2) * x1: функція (x1, x2)=(- 2 + x1) * x2
=0.01кон=2
Отримано фазовий портрет:
<...
