Курс лекций по теоретическим основам электротехники. Исаев Ю.Н - 154 стр.

UptoLike

154
Определяем собственные числа матрицы состояния A => λ
λ eigenvals A():=
λ
1210.96
2454.41i
+
1210.96 2454.41i
244.75
=
Для проверки определяем корни характеристического уравнения через импеданс
схемы Z(p)
p
R2
1
C1 p
R2
1
C1 p
+
Lp+ R1+
R2
1
C2 p
R2
1
C2 p
+
+
solve p,
float 6,
1210.96 2454.41i
1210.96 2454.41i+
244.752
:=
Для проверки определяем принуждённые составляющие
i
Lпр
E
R
1
R
2
+ R
2
+
:=
u
C1
i
Lпр
R
1
R
2
+
(
)
:=
u
C2
i
Lпр
R
2
:=
u
C1
54.545=
i
Lпр
0.455=
u
C2
45.455=
A
1
B
54.545
45.455
0.455
=
Теперь обращаемся к одной из стандартных процедур решения
системы дифференциальных. уравнений
Dt x,() Ax
+:=
Расширенная матрица
τ
1
max Re λ
()()
⎯⎯⎯
:=
τ 0.0041=
T 6
τ
:=
y rkfixed
0
0
0
0, T, N, D,
:=
Метод Рунге Кутта
ty
0
:=
0 0.0041 0.0082 0.0123 0.0163 0.0204 0.0245
0
10
20
30
40
50
60
y
1
〈〉
t