Ведаючи пачатковия Сћмови, разлічим сістему диференцияльних раСћнанняСћ (3). Лікави разлік сістеми (3) виканаем плиг дапамозе праграма разліку сістеми диферанцияльних раСћнанняСћ Метад Рунге-Кутта чацвертага парадку (праграма DIFRK4 у лабараторнай рабоце № 7). Неабходна зрабіць адаптапцию гета праграма для разліку сістеми (3). Змяненні праводзім у критим, што пішам нову падпраграму PRAV для разліку правих частак. p align="justify"> Файл уваходних дадзених для маделявання пераходнага режиму Сћ схемою пасли замикання ключа К травні виглядає:
1 0. 0.1 0.0001
2.7239 13.6196
Винікі разліку пераходнага режиму виводзяцца праграмай у виглядзе стаСћбцоСћ з дадзенимі Сћ файл з паширеннем REZ и GRF.
Графікі пераходнага режиму змешчани на малюнках.
Падпраграма правих частак для сістеми диферанцияльних раСћнанняСћ (3) у будз мець виглядає:
SUBROUTINE PRAV (X, Y, F)
REAL L1
DATA C2, C3, L1, R1, R2, R3, Em, fi/0.0001, 0.00002,0.2,25., 15., 25., 60.,
* 35./
DIMENSION Y (*), F (*)
! ЕЛЕКТРАРУХАЮЧАЯ Сіла У МОМАНТ x
e = Em * sin (314. * x + fi/180. * 3.14159)/sqrt (2.)
! ВИЛІЧЕННЕ КАЕФІЦИЕНТАРЋ, ШТО СТАЯЦЬ побачим Зх ПЕРАМЕННИМІ
quo11 = - (R1 * R2 + R2 * R3 + R1 * R3)/(L1 * (R2 + R3))
quo12 =-R2/(L1 * (R2 + R3))
quo13 =-R3/(L1 * (R2 + R3))
quo14 = e * R3/(L1 * (R2 + R3))
quo21 = R2/(C2 * (R2 + R3))
quo22 = -1./(C2 * (R2 + R3))
quo23 = 1./(C2 * (R2 + R3))
quo24 =-e/(C2 * (R2 + R3))
quo31 = R3/(C3 * (R2 + R3))
quo32 = 1./(C3 * (R2 + R3))
quo33 = -1./(C3 * (R2 + R3))
quo34 = e/(C3 * (R2 + R3))
! СІСТЕМА ДИФ.РАРЋНАННЯРЋ
F (1) = quo11 * Y (1) + quo12 * Y (2) + quo13 * Y (3) + quo14
F (2) = quo21 * Y (1) + quo22 * Y (2) + quo23 * Y (3) + quo24
F (3) = quo31 * Y (1) + quo32 * Y (2) + quo33 * Y (3) + quo34
RETURN
END
Ніжей приведзени файл з винікамі разліку сістеми диферинцияльних раСћнанняСћ (файл приведзени толькі часткова).
DIFRK4 - РАШЕННЕ СIСТЕМИ ДИФЕРЕНЦИЯЛЬНИХ РАУНАННЯУ
УВАХОДНИЯ Дано
1 0. 0.0299999993 9.99999975E-005
0. 2.72390008 13.6196003
РЕЗУЛЬТАТИ РАЗЛIКУ
xy (1) y (2) y (3) y (4) y (5) ...
