{параметри ланцюга}
E0=10.0
R1=1.0
R2=1.0
L=1.0
C=1.0
k=1.5
{початкові умови}
tнач=0.0
Il=0.0
Uc=0.0
t=tнач
{масштабні коефіцієнти}
mi=0.2
mu=0.2
Собственние_Значенія (J, 300, 1E - 15) -> (Re, Im)
{параметри обчислювального експерименту}
h=0.1
tкон=10
{установка початкової точки графіків Il (t) і Uc (t)}
Ізменіть_Абсціссу (t) -> (Точка)
Ізменіть_Ордінату (Il * mi) -> (Точка)
Добавіть_Елемент (Точка) -> (Il_t)
Ізменіть_Ордінату (Uc * mu) -> (Точка)
Добавіть_Елемент (Точка) -> (Uc_t)
Оновити (графіка)
t=tнач
{розрахунковий цикл інтегрування}
поки t < tкон
початок
{розрахунок методом Рунге-Кутта 4-го порядку}
Iln=Il Ucn=Uc
{1-й подшаг}
tt=t=h * dIldt du1=h * dUcdt
{2-й подшаг}=t + h / 2.=Iln + di1 / 2. Uc=Ucn + du1 / 2.=H * dIldt du2=h * dUcdt
{3-й подшаг}=t + h / 2.=Iln + di2 / 2. Uc=Ucn + du2 / 2.=H * dIldt du3=h * dUcdt
{4-й подшаг}=t + h=Iln + di3 Uc=Ucn + du3=h * dIldt du4=h * dUcdt=Iln + (di1 + 2. * di2 + 2. * di3 + di4) / 6.0=Ucn + (du1 + 2. * du2 + 2. * du3 + du4) / 6.0=t + h
Ізменіть_Абсціссу (t) -> (Точка)
Ізменіть_Ордінату (Il * mi) -> (Точка)
Добавіть_Елемент (Точка) -> (Il_t)
Ізменіть_Ордінату (Uc * mu) -> (Точка)
Добавіть_Елемент (Точка) -> (Uc_t)
кінець
Чісло_Елементов (Il_t) -> (N)
цикл i=1 до n
початок
Ордината (Il_t [i]) -> (Y)
Ордината (Uc_t [i]) -> (X)
Добавіть_точку (y, x) -> (IlUc)
кінець
Оновити (ГРАФІКА)
КІНЕЦЬ
Пірімери побудови фазових портретів для нелінійних систем
Дана система диференціальних уравнеій:
Модель для системи Тополог:
ІМ'Я «Рішення системи диференціальних рівнянь»
ДАНІ: вещ: вещs: вещs: вещ_мін: вещ=- 1.0_макс: вещ=3.5_мін: вещ=- 1.0_макс: вещ=1.5_1: вещ_2: вещ: вещ
Шаг_Сеткі: вещ=0.2
{параметри обчислювального експерименту}: вещ: вещнач: вещкон: вещ: вещ
{робочі переменнние методу інтегрування} n: вещ X2n: вещ: вещ di2: вещ di3: вещ di4: вещ: вещ du2: вещ du3: вещ du4: вещ: функція (x1, x2)=(3 - x1 -5 * x2) * x1: функція (x1, x2)=(2 - 4 * x1 - 2 * x2) * x2
{вектор похідних}: вектор=(dX1, dX2) |
{матриця Якобі}: якобиан (V)
...
