El blog de Daniel Muñoz “No hay mayor riqueza que el conocimiento ni mayor pobreza que la ignorancia.” – Alí ibn Abi-Talib

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