Как добраться из системы уравнений MNA в пространство состояний [A, B, C, D] mimo в MATLAB? - PullRequest
Новая
       58

Как добраться из системы уравнений MNA в пространство состояний [A, B, C, D] mimo в MATLAB?

0 голосов
/ 24 января 2019

Я хочу программно сгенерировать систему пространства состояний [A B C D] из сгенерированной системы уравнений MNA (модифицированного узлового анализа). Я не знаю, какие функции MATLAB вызывать.

Я использую скрипт SCAM link от Эрика Чивера, чтобы сгенерировать систему уравнений.

Я создал файл example3.cir link , представляющий собой RLC-схему Spice с источником входного сигнала Vs.

Я также пробовал netlist2linss ссылка github, но он не находится в домене Лапласа и не работает правильно. 

example3.cir:
R1 2 1 R
C1 3 0 C
L1 2 3 L
Vs 1 0 V

>> fname = 'example3.cir'; scam; system_of_eqns = [Z == A*X]
Solved variables:
v_1
v_2
v_3
I_Vs

system_of_eqns =
                      0 == I_Vs + v_1/R1 - v_2/R1
 0 == v_2*(1/R1 + 1/(L1*s)) - v_1/R1 - v_3/(L1*s)
          0 == v_3*(C1*s + 1/(L1*s)) - v_2/(L1*s)
                                        Vs == v_1

>> [A_,b_] = equationsToMatrix(system_of_eqns, X)

A_ =
[ -1/R1,              1/R1,                 0, -1]
[  1/R1, - 1/R1 - 1/(L1*s),          1/(L1*s),  0]
[     0,          1/(L1*s), - C1*s - 1/(L1*s),  0]
[    -1,                 0,                 0,  0]

b_ =
0
0
0
-Vs

ожидается (напряжение на конденсаторе в качестве выхода): 

A =
[     0,  1/C1 ]
[ -1/L1, -R/L1 ]

B =
[ 0, 1/L1 ]

C =
[ 1, 0 ]

D = 
[ 0 ]

states = [sVC1 sIL1]

input =  [ Vs ]

1 Ответ

0 голосов
/ 26 июня 2019

Я нашел решение в python под названием lcapy: https://github.com/mph-/lcapy

...