Hayth example 13.

2
0 0 0 0 0 0 𝑒1̇ 𝐺1 −𝐺1 0 1 0 0 𝑒1 0
0 0 0 0 0 0 𝑒2̇ −𝐺1 𝐺1 0 0 1 0 2 𝑒 0
0 0 0 0 0 0 𝑒3̇ 0 0 𝐺2 0 0 1 𝑒3 0
+ +
0 0 0 0 0 0 𝑖𝑠̇ 1 0 0 0 0 0 𝑖𝑠 −𝑉𝑠
0 0 0 0 −𝐿1 −𝑀12 𝑖1̇ 0 1 0 0 0 0 𝑖1 0
[0 0 0 0 −𝑀21 −𝐿2 ] [ 𝑖2̇ ] [ 0 0 1 0 0 0] [ 𝑖2 ] [ 0 ]
Discretizando C con Euler inverso:
𝐿1 𝑑𝐼𝐿1 𝑀12 𝑑𝐼𝐿2
𝜐𝐿1 = +
𝑑𝑡 𝑑𝑡
𝐿1 𝑛+1 𝑀
𝜐𝐿𝑛+1 = ( 𝑖1 − 𝑖1𝑛 ) + 12 ( 𝑖2𝑛+1 − 𝑖2𝑛 )
1 ℎ ℎ
Y por lo tanto para L2:
𝐿2 𝑑𝐼𝐿2 𝑀21 𝑑𝐼𝐿1
𝜐𝐿2 = +
𝑑𝑡 𝑑𝑡
𝐿2 𝑛+1 𝑀
𝜐𝐿𝑛+1 = ( 𝑖2 − 𝑖2𝑛 ) + 21 (𝑖1𝑛+1 − 𝑖1𝑛 )
2 ℎ ℎ
Los valores por obtener en la ecuación son los de subíndice “n+1” y los valores con subíndice “n”
son conocidos, que van al vector f, similar a la estampa para una fuente de voltaje independiente.
La nueva matriz queda de la forma:
𝐺1 −𝐺1 0 1 0 0 𝑒1 𝑛+1 0 𝑛
−𝐺1 𝐺1 0 0 1 0 𝑒2 0
0 0 𝐺2 0 0 1 𝑒3 0
1 0 0 0 0 0 𝑖𝑠 + −𝑉 𝑠
=0
0 1 0 0 −𝐿𝑒𝑞 −𝑀𝑒𝑞12 𝑖1 𝑉𝑒𝑞1
1
[ 0 0 1 0 −𝑀𝑒𝑞 [
−𝐿𝑒𝑞 ] 2 𝑖 ] 𝑉
[ 𝑒𝑞 2 ]
21 2

Donde:
𝐿1
𝐿𝑒𝑞 =
1 ℎ
𝐿2
𝐿𝑒𝑞 2 =
ℎ
𝑀12 𝑀21
𝑀𝑒𝑞 = 𝑀𝑒𝑞 21 = =
12 ℎ ℎ
𝐿1 𝑛 𝑀12 𝑛
𝑉𝑒𝑞1 = 𝑖 + 𝑖
ℎ 1 ℎ 2
𝐿2 𝑛 𝑀21 𝑛
𝑉𝑒𝑞2 = 𝑖 + 𝑖
ℎ 2 ℎ 1