Método circuital o de las corrientes de mallas con matlab.
Comparto un script de matlab que hice para resolver circuitos en AC mediante el método de las corrientes de malla o método circuital (ver también método nodal).
Por ejemplo, para resolver el siguiente circuito:
Habría que ejecutar lo siguiente en matlab (muestro en negrita lo ingresado, el resto es generado por el script):
>> mallas
MÉTODO DE LAS MALLAS EN EL DOMINIO DE LAPLACE.
---------------------------------------------------------------------
¿Cuantas mallas tiene el circuito?: 3
¿Cuantos dígitos de precisión quiere ver en el resultado?: 3
---------------------------------------------------------------------
MATRIZ DE IMPEDANCIAS.
A continuación ingrese los elementos de la matriz de impedancias en el
dominio de Laplace utilizando "s" como variable generalizada de Laplace.
***Sepa que debe ingresar las coimpedancias con su SIGNO CORRESPONDIENTE**
Ingrese z[1,1]: 6
Ingrese z[1,2]: -2
Ingrese z[1,3]: -2
Ingrese z[2,2]: 6
Ingrese z[2,3]: -2
Ingrese z[3,3]: 6
Así quedó la matriz de impedancias:
[ 6, -2, -2]
[ -2, 6, -2]
[ -2, -2, 6]
¿Está bien?(S/N): s
---------------------------------------------------------------------
MATRIZ DE TENSIONES.
Ingrese V[1]: 10/s
Ingrese V[2]: 0
Ingrese V[3]: 0
Así quedó la matriz de Tensiones:
10/s
0
0
¿Está bien?(S/N): s
---------------------------------------------------------------------
DETERMINANTE PRINCIPAL
El determinante principal (Dp) es:
128
---------------------------------------------------------------------
PARA ENCONTRAR i1(t).
La matriz que origina el determinante sustituto 1 (Ds1) es:
[ 10/s, -2, -2]
[ 0, 6, -2]
[ 0, -2, 6]
El Ds1 es:
320
---
s
Ds1/Dp:
5
---
2 s
Expresado en fracciones parciales:
5
---
2 s
Antitransformando tenemos i1(t):
5/2
Expresado i1(t) con formato punto decimal:
2.5
---------------------------------------------------------------------
PARA ENCONTRAR i2(t).
La matriz que origina el determinante sustituto 2 (Ds2) es:
[ 6, 10/s, -2]
[ -2, 0, -2]
[ -2, 0, 6]
El Ds2 es:
160
---
s
Ds2/Dp:
5
---
4 s
Expresado en fracciones parciales:
5
---
4 s
Antitransformando tenemos i2(t):
5/4
Expresado i2(t) con formato punto decimal:
1.25
---------------------------------------------------------------------
PARA ENCONTRAR i3(t).
La matriz que origina el determinante sustituto 3 (Ds3) es:
[ 6, -2, 10/s]
[ -2, 6, 0]
[ -2, -2, 0]
El Ds3 es:
160
---
s
Ds3/Dp:
5
---
4 s
Expresado en fracciones parciales:
5
---
4 s
Antitransformando tenemos i3(t):
5/4
Expresado i3(t) con formato punto decimal:
1.25
