E S C U E L A P O L I T É C N I C A N A C I O N A L

FACULTAD DE INGENIERÍA ELÉCTRICA

PROGRAMA DIGITAL PARA ANÁLISIS DE CIRCUITOS ELÉCTRICOS

EN EL DOMINIO DE LA FRECUENCIA

TESIS DE G R A D O

QUITO, MARZO DE 1983

\H DEL ROSARIO RODRÍGUEZ SALAZAR

CERTIFICO QUE EL PRESENTE TRABAJO
DE TESIS HA SIDO. REALIZADO EN SU
TOTALIDAD POR LA SEÑORITA ELIZABETH
RODRÍGUEZ S.

Ing. Efraín Del Pino

DIRECTOR DE TESIS
DEDICATORIA

A mis padres

y
A Delita
 Í N D I C E

INTRODUCCIÓN

CAPITULO PRIMERO

1.1. A n á l i s i s de sistemas y circuitos eléctricos ' 1

1.2 Lineal i dad. 3
1.3 Análisi-s sinusoidal de estado permanente 4
1.4 Desarrollo sistemático de ecuaciones eléctricas
para el análisis de circuitos 6
1.4.1 Matriz de incidencia 16
1.4.2 Análisis de producto [Link] 20
1.5 Ejemplos de aplicación
1.5.1 Circuito'DC 22
1.5.2 Circuito AC 26
CAPITULO SEGUNDO
DESCRIPCIÓN DE ALGORITMOS EMPLEADOS Y SOLUCIÓN DE LAS
ECUACIONES QUE DESCRIBEN EL CIRCUITO
2.1 Introducción • 30
2.2 Estructuración de la matriz incidencia (A) y ma
triz de admitancia ramal (Yb)
2.2.1 Matriz incidencia (Estructuración) 32
2.2.2 Estructuración de matriz de admitancia ramal Yb 39
2.2.3 Estructuración de Yb cuando existen inductancias
mutuas 42
.2.3 Producto ([Link]) AT - Yn 45
2.4- Solución numérica de ecuaciones lineales simultá
neas _ 46
2.5 Inversión de matrices
2.5.1 Método de Gauss . . 54
2.5.2 Método de Shipley 57

CAPITULO TERCERO
3.1 Estructura general de los programas 62
3.2 Programa maestro CIRCUITO 64
3.3 Programa PR/1NGRESO 67
3.4 Programas: PR/FREQ, PR/OPTIOC 70
3.5 Programa PR/CALCULO 75
3.6 Programas: PR/YE1, PR/YE2 78
3.7 Análisis para varias frecuencias 82
CAPITULO CUARTO
1, Ejemplo de aplicación
2. Conclusiones y recomendaciones

ANEXO 1
Listado de programas

ANEXO 2
MANUAL DE USO DEL PROGRAMA CIRCUITO
Descripción del circuito
Algunas consideraciones sobre la descripción del
circuito
Modo de empleo, ejemplo
 I N T R O D U C C I Ó N

El presente trabajo-constituye un conjunto de programas para el

a n á l i s i s de circuitos eléctricos l i n e a l e s usando un computador. Estos
programas han sido desarrollados e implementados en el mini - computador
TEXTRONIX 4051 de la Facultad de Ingeniería Eléctrica y formarán parte
de la biblioteca de programas que se dispone en el. Departamento de Con-
trol .

El espíritu o m o t i v a c i ó n para el desarrollo de este tema es la

creciente necesidad de disponer de un programa digital de uso sencillo,
dinámico y de disponibilidad inmediata, lo que lo'convierte en un instrj¿
mentó útil y puesto al día para f a c i l i t a r el diseño y a n á l i s i s de circui_
tos eléctricos l i n e a l e s en el dominio de la frecuencia para excitaciones
sinusoidales en estado estable.

Como datos se tienen los valores de los elementos que constitu-

yen el circuito y la topología del mismo. Los elementos del circuito
son elementos pasivos (resistencias, i n d u c t a n c i a s , inductancias mutuas,
capacidades, c o n d u c t a n c i a s ) , fuentes i n d e p e n d i e n t e s de voltajes y/o co-
rrientes y fuentes controladas que son particularmente útiles en la es -
tructuración de modelos de dispositivos electrónicos. A base de estos
datos el programa procede a estructurar las ecuaciones que describen to-
das las relaciones entré voltajes de nodo y corrientes de rama..

El método de análisis que se emplea se basa en un desarrollo si_s_

temático de ecuaciones eléctricas fasoriales (se obtienen de un análisis
nodal). Estas ecuaciones implican operaciones y resoluciones matricia -
les, apropiadas para implementarlas directamente en un computador digi-
tal. Sin embargo este tema tiene una orientación escencialmente prácti-
ca por lo que-se ha puesto especial atención en la estructuración de es-
tas ecuaciones con el fin de economizar localizaciones de memoria, para
emplear de mejor manera el área disponible en el procesamiento y análi-
sis de circuitos -largos o complicados.

Este programa puede realizar el análisis para corriente continua

(da) y/o alterna (una o varias frecuencias). Como soluciones se obtie -
nen todos los voltajes de la red,, pudiéndose luego calcular las corrien-
tes de rama. Para un rango de frecuencias puede obtenerse especifícame^
te el voltaje de un nodo (Vni), una corriente de rama (Ibi), la diferen-
cia entre dos nodos cualesquiera (Vni - Vnj) o relaciones (Vni/Vnj),
(Ibi/Ibj), (Vni/Ibj), (Ibi/Vnj).

Los resultados se los presenta en forma numérica y/o gráfica del

módulo o la fase en función de la frecuencia.

El tema del trabajo se desarrolla de la siguiente forma:

En el Primer Capítulo se expone a manera de introducción ideas

sobre el análisis de sistemas y circuitos eléctricos, conceptos básicos
sobre la linealidad y el análisis sinusoidal de estado permanente, para
continuar con la exposición sistemática de las ecuaciones que describen
la red eléctrica.

En el Segundo Capítulo contiene la descripción de los algorit -

mos empleados en la elaboración y solución de las ecuaciones que descri_
ben el circuito.

En el Tercer Capítulo se explica la estructura general de la bj_

blioteca de programas y su descripción de una manera clara y sencilla..

En el Cuarto Capitulo como aplicación se desarrollan [Link].s -.

de utilidad práctica. Se obtienen conclusiones deducidas de la [Link].e
experimental y se hacen los comentarios pertinentes.

Se incluyen en un anexo el listado de los programas y un manual

de uso .al que se lo presenta de una manera consisa y lógica para fac.ilj_
tar la consulta rápida y permitir una utilización eficaz,

Para concluir deseo expresar mi especial agradecimiento al

Ing. Efraín del Pino por su acertada dirección, sus múltiples y útiles
sugerencias han permitido la culminación de este trabajo y a la Escuela
Politécnica Nacional.
 C A P I T U L O P R I M E R O

1.1 Análisis de sistemas y circuitos eléctricos

1.2 Linealidad

1.3 Análisis sinusoidal de estado permanente

1.4 ' Desarrollo sistemático de ecuaciones eléctricas para el análisis

de circuitos

1.4.1 Matriz de incidencia

1.4.2 Análisis del producto A,Ib

1.5 Ejemplos de aplicación

1.5.1 Circuito DC

1.5.2 -Circuito AC
 C A P I T U L O P R I M E R O

l'.l ANÁLISIS DE SISTEMAS Y CIRCUITOS ELÉCTRICOS

• Un -Sistema en general de ingeniería, es un conglomerado de comp_o_

nentes físicos que ejecutan en forma conjunta una función útil. AsT los

sistemas electrónicos por ejemplo se pueden definir como combinaciones -

de circuitos que manejan señales. Las señales pueden ser voltajes o co-

rrientes, pero si existen transductores de energía las señales pueden

ser otras variables físicas.

La distinción entre circuitos y sistemas depende más del punto

de vista del sistema que- de complejidad. El punto de vista del sistema

enfatiza en las funciones operacionales o relaciones entre la entrada y

salida. El punto de vista del circuito trata de la determinación de co-

rrientes y voltajes para ramales o mallas de un circuito.

Un conocimiento_ de las posibilidades y limitaciones de los cir-

cuitos s e - h a c e necesario para el diseño práctico de sistemas, por lo ta£

to la piedra angular del diseño es el análisis sistemático de estos cir-

cuitos.
 - 2-

El análisis y diseño de sistemas eléctricos, electrónicos se fa-

cilita por el uso de modelos de circuito aproximados del sltema físico,
ya que el objetivo consiste en substituir los dispositivos o sistemas fí_
slcos con modelos de circuito apropiados y luego completar el análisis
del sistema de acuerdo con el modelo general.

El modelo particular que se utiliza para un dispositivo depende-

c- 5
rá de la presiden requerida por lo tanto este modelo no es único» exis-
ten muchos modelos posibles y antes de iniciar el análisis es necesario
hacer una elección de entre aquellos que se disponga.

Asi por ejemplo se ilustra para un transistor bipolar tres mode-

los que representan la complejidad cada vez mayor de acuerdo con el nüme_
ro de elementos que se utilizan, aparecen en (b), (c), y (d) de la
fig. 1.1.
 - 3 -

— cb
3-1

o-
(b)

0 II nc

un E ^ Q *- L C b ^ Ag
V
^
m 6e
v«

CcJ
11
Ce.

o II11 1

-C 0 <S| -
>,„,„, =
O
Cd)
1-1-

1.2 LINEALIDAD

Un sistema se define como lineal en términos -de su excitación y

respuesta. En el caso de las redes eléctricas en general se puede defi-
nir una corriente o un voltaje de respuesta Y y una corriente o un
taje de excitación x . En tal caso, si Y = f (x) la red es lineal
si :
 a Y - f (a x )

es decir Y es proporcional a x. (Propiedad de homogeneidad)

Además una condición necesaria para determinar una red lineal es

que la respuesta producida por la aplicación simultánea de varias funci_o_
nes excitadoras distintas, es la suma de las respuestas individuales.
Por lo tanto, para redes lineales se puede calcular la respuesta a diver
sas entradas, tratando una entrada a la vez, y añadiendo o sumando los
resultados.

Muchos circuitos eléctricos prácticos contienen elementos elec -

irónicos no lineales como los transistores pero para ciertos análisis se
puede considerar que son lineales para una gama restringida de operación
y proceder a su análisis según su modelo circuital aproximado.

1.3 ANÁLISIS SINUSOIDAL DE ESTADO PERMANENTE

El análisis de circuitos en régimen permanente sinusoidal tiene

una gran importancia, no solo porque el senoide es la forma de onda que
se usa comunmente en la transmisión de información y energía. Sino que
hay otras señales periódicas que se pueden representar por una suma de
senoides y el análisis de circuitos con fuentes periódicas no sinusoida-
les pueden realizarse convenientemente empleando la superposición.

La función sinusoidal es la única función periódica real que

 _ 5 -

tiene una integral y una derivada con la misma onda sinusoidal. En cuaj_
quier circuito lineal ello significa que la respuesta, en régimen perma-
nente, a una fuente sinusoidal será también sinusoidal prescindiendo del
tipo o número de elementos de circuitos que se conecten.

Mediante el empleo de la notación fasorial en la representación

de las tensiones' e intensidades de corriente se facilita el estudio en
régimen permanente senoidal. Sobre todo puede reducirse el cálculo a la
solución de un conjunto de'ecuaciones algebraicas con coeficientes de va_
lores complejos.

La relación de estado estacionario entre los fasores I y V de

elementos simples son:

1.) Resistencia VKn = R I

Vn
impedancia Ze = — =R
I

2) Inductancía V, = j w L I

impedancia Ze = — =jwL
I

admitancia - Ve = — -
w L
 - 6 -

3) Capacidad Ve = —¿
w C

impedancia - i
Zc = Ve = —**—
I w C

admitancia .Ye = j w C

El método de análisis que se utilizará es un médodo nodal con

orientación .al uso del computador digital.

Es interesante observar que los métodos de análisis nodal se pre

fieren en forma casi universal debido a su generalidad y a la facilidad
con que se programa.

1.4 DESARROLLO SISTEMÁTICO DE ECUACIONES ELÉCTRICAS PARA EL ANÁLISIS

DE CIRCUITOS

Para el análisis consideramos a la red constituida por m • núme-

ro de ramas y n números de nodos independientes.

Se [Link] un NODO como el punto en el cual dos o más elementos

tienen conexión común y una RAMA como el camino que conecta un nodo con
otro nodo cualquiera.

En el análisis nodal se hace uso de los voltajes de nodo como va_

Hables auxiliares. Llamamos voltaje de nodo, los voltajes de los nodos
de red referido a un nodo denominado de referencia, que generalmente se
 - 7-

lo considera a cero voltios. De una manera general la rama se considera^

rá constituida por:

1) Un elemento pasivo que puede ser una (R, G, L, C) de admitan-

cia Ye.
2) Fuentes independientes de voltaje y de corriente, de valor

VQ, 13-
3) Fuentes dependientes o controladas de voltaje y corriente de
valor vd, i¿,

Estos elementos constitutivos de una rama generalizada se mués -

tran en la figura 1.2 con los sentidos y polaridades que se emplearán
para el análisis,- y obtención de las ecuaciones eléctricas.

ve
ib -O
Nodo A Nodo B

-vb-

FIG. 1.2.
 Aplicando la primera ley de Kirchhoff a uno de los nodos A o B
se tiene:

ib + ig = ie + id - (1.1)

donde ib = corriente de rama

ie = corriente que circula por el elemento
ig = fuente independiente
id'= fuente controlada

El voltaje de rama ub será la suma algebraica de:

ub = ve + vd - vg (1-2)

ve = diferencia de potencial de elemento

vg = fuente independiente
vd = fuente controlada
y. - ie = Ye ve (1.3)

Los problemas en análisis de circuitos lineales (descriptible en

términos de ecuaciones algebraicas lineales simultáneas) pueden organizajr
se en forma compacta con la ayuda del algebra de matrices. La notación
de matrices no elimina el trabajo computacional requerido para la solu -
ción numérica pero permite la formulación del problema en un estilo uni-
ficado y ofrece un programa muy sistemá[Link] las soluciones aritméti_
cas.

Asi para una red ramal arbitraria con m ramas se define:

 - 9-

vector corriente de rama 16

Ib = [ ibia ib2, . . . ibm ]

—> ib-c corriente de la rama -I ¿ = 1, 2, 3 ...m

simultáneamente se definen los vectores:

Ig = [ igi> 192* . . . igm ]

=> ig¿ corriente de fuente independiente en la-rama

le = [ iei, ie2, - - . iem ]

—> ie¿ corriente de elemento en la rama -I

vector de voltaje de rama Vb

Vb = [ ubi, vbz, . . . u b m ] vb-¿ = voltaje de rama 4,

Ji - 1, 2, 3 . . . m

de la misma manera Vg, Vd, Ve vectores de voltajes de fuentes indepen -

dientes, de fuentes dependientes y de elementos.

Adicionalmente le = Ye Ve
donde Ye es una matriz simétrica de orden m

Ye = Ze"1
donde la diagonal principal esta constituida por las admitancias propias
de la rama yi¿ =' rama -i/ es resistiva
 - 10 -

= j w C Á, rama ¿/ es capacitiva

rama ¿/ es inductiva
j w U/

La matriz Ye tiene nulos todos los elementos excepto los de

la diagonal principal cuando no existe acoplamientos mutuos.

Al existir inductancias mutuas

significa que existe acoplamiento inductivo entre las ramas -i y

FUENTES CONTROLADAS

Una fuente controlada es una fuente cuyo valor no es independiej]_

te del resto del circuito, sino que es proporcional a la corriente o al
voltaje de otras ramas cualesquiera.

Existe cuatro tipos de fuentes controladas.

1. F C C C Fuente de Corriente controlado por Corriente

2. F V C C Fuente de Voltaje controlado por Corriente
3. F C C V Fuente de Corriente controlado por Voltaje
4. F V C V Fuente de Voltaje controlado por Voltaje
K - rama controlada
 - 11 -

j - rama controlante

1. F C C C ife= pfe.y iy

3fey fac-
tor de ampl if.
de corriente.

F V C C

vfe = rfe./ iy
rkj'-

3. F C CV

ik = gkj

= concJuctancia

4. F V C V

+* == factor
vj
de amplifica-
ción de volta-
je.
FIG. 1.3
 - 12 -

Igualmente se puede definir para una red ramal de m ramas las

matrices de tamaño ( m x m )

B = [ efe/ ]

R = [ rfe/ ]

G = [ gfe/ ]

M - [

Por consiguiente las corrientes controladas

Idi = [ Bfe/ ] le = B le
Id2 - [ gfe/ ] Ve - G Ve
Id. = B le + G Ve pero Ve = Ye -"-le

Id = B le + G Ye~1Ie
Id = ( B + G Ye"1 ) le
Id = F le con: F - B + G Ye a (1.4)

Los voltajes controlados:

Vd! = [ r/fe ] le - R le
Vd 2 - [ y/fe ] Ve = M Ve

Vd = R le + II Ve pero le - Ye Ve

Vd = R Ye Ve + M Ve = ( R Ye + M ) Ve

Vd = D Ve con D - R Ye + M (1,5)

De la ecuación (1.2) se tiene que:

 - 13 -

en términos matriciales tendremos:

Vb + Vg - Ve + Vd pero Vd - D Ve ecuación (1.5)

Ve + D Ve - Vb + Vg
( U + D ) Ve = Vb + Vg
Ve - ( U + D )-1 (Vb + Vg) (1.6)

De la ecuación (1.1) en términos matriciales tenemos que:

Ib + .Ig =- le + Id pero Id = F le
Ib + Ig = le + F le
Ib + Ig = ( U + F ) le *
e
pero Ib' ~ Ye Ve . reemplazando la ecuación (1.6)
se tiene:
le - Ye ( U + D )~1 ( Vb + Vg)
por lo tanto:
Ib = ( U + F ) Ye (U + D)'1 ( Vb + Vg ) - Ig
donde Yb = ( U + F ) Ye ( .U + D J'1 matriz admitancia de rama (1.7)
Ib - Yb (Vb + Vg) - Ig (1.8)

INDUCTANCIAS MUTUAS

Un componente adicional que debe ser considerado para el anal i -

sis, es la inductancia mutua. En general puede existir cualquier número
de ramas acopladas entre sí, donde la relación volt-amperimétrica para -
otros elementos (R, 6, L, C) incluidos las inductancias mutuas queda de-
 - 14 -

finida por:

Ve = Ze le

Ze es la matriz de impedancias de elementos, tiene m filas y m colum_

ñas ( m es el número de ramas del circuito), por lo tanto:

i = Rx si la rama -i* contiene una resis_

de valor fU.
.¿ = l/(j w¿) Si la rama -i/ contiene un capacj_
tor de valor C¿
I-u. = jw LL Si la rama AJ contiene una in -
ductancia de valor L¿
I ~ jw Mcfe- para toda inductancia mutua que
acopla la rama -I/ con la fe.
Los e'lementos de la matriz Ye pueden obtenerse por la ecuación
•
de Ze

Ye - Ze'1

ACOPLAMIENTO INDUCTIVO ENTRE DOS RAMAS

o-

roma p 'rama S
)Lp Ls
\

Sean las ramas P y S, ramas con acoplamiento inductivo

 - 15 -

Las relaciones básicas de voltaje y corriente para astas dos i_n

4,

ductancias acopladas por Mps = Msp = M son las siguientes: .

Ep = jw Lp Ip + jw Mps Is . (1.9)

Es = jw Ls Is + jw Msp Ip (1.10)

Por lo tanto las corrientes en cada rama serán:

Tn - -EP Lp Ls Es M Ls (1.11)
jw Lp (Lp Ls ~,M2) jw Ls (Lp Ls - M2)'

Is = Es LP Ls EP M LP (1.12)
jw Ls (Lp Ls.- M2) jw Lp (Lp-Ls - M2)

La interrelación entre las inductancias propias Lp, Ls y la i_n

ductancia mutua, se expresa asf:

M
= K (1.13)
/Lp Ls

donde K es denominado, coeficiente de acoplamiento entre las bobinas.

Físicamente hablando, K puede acercarse a la unidad, pero nunca sobre -
pasarla, porque K es sencillamente la fracción del flujo de una bobina
que se acopla a la segunda.

Todas las .bobinas reales, tienen un K < 1 y pueden acercarse a

valores 0,915 a 0,998 (mediante^un diseño apropiado). No se considera
 - 16 -

en este análisis al "transformador ideal", comúnmente utilizado en ana

Tisis de circuitos teóricos. .

1.4.1 Matriz de Incidencia A

a) Relaciones entre voltajes de ramas y voltajes de nodos

En una red de . m ramas y n nodos independientes definimos:

Vb ~ [ubi, ub 2) - - . vbm] vector de voltajes de las ramas

Vn = [uní, unn] vector de voltajes nodal

El voltaje de cada rama puede escribirse como una diferencia

de potencial entre dos nodos así:

n - número de nodos = 3
m = número de ramas = 5

1.4.

ubi
 - 17 -

De una manera general para la rama ft

rama

vbk

Sean v¿, vy. los voltajes en los nodos L,- / respectivame_n_

te.

El voltaje de la rama ft es:

notando al nodo ¿ como nodo "desde"

/ corno nodo "hasta"

.-Considerando todos los voltajes de rama del ejemplo y expresa^

dos en forma matricial se tiene:
 - 18 -

"Vbx. 1 -1 0 "Vn,"

Vb 2 o - i - i Víl2

Vb 3 0 0 - 1 Vn 3
rr X
Vb 4 0 1 0

Vb s - 1 0 0
— ~_ i _ _

5X1 5X3 3X1

Generalizando para una red de m número de ramas y n número

de nodos será:

Rama
"vb, .1 bu f r k z - b j Á by .. b m "vm"
Vb 2 2 bzi b 2 2 . b 2 i b 2 j .. b 2n Vn 2

Vb 3 3 bai b32 - b3i b 3 j . . b 3tt Vn^

• . . , . . . .
Vnj
• — - . . -. x

Vbfc ft bfei bk2 . bfe^ bfe/ . . bfen

Vbm m bmi bm¿ -. b-m.£ bm/ bmn Vn n

m x 1 m x n íl X 1

Donde la matriz [bfe¿] es de orden m x n y para l a rama fe

(o f i l a fr) se tiene que:
 - 19 -

bk¿ - I si ¿ es un nodo "desde"

bkj = -1 si / es un nodo "hasta"

Se define la transpuesta de esta matriz como matriz incidencia

A = ]T

Vb = AT. Vn (1-14)

Si A = [ b/¿¿ ] entonces la matriz incidencia A es de orden

n x m y para el ejemplo dado Fig. 1.4.

2 3 4 5
1
I 1 0 0 0 -1

A -2 -1 1 0 1 0

3 0 -1 -1 0 0

A - [ . L ~ 1,2, ... n
k = 1,2, ... m

y
1 Si la rama ft sale "desde" nodo x-
-1 Si la rama ^ entra "hasta" nodo -i
O Si l a rama I? no encide al nodo -I
 - 20 -

1.4.2 Análisis del producto A,Ib

Consideremos el siguiente ejemplo:

2

b5

b3
b1 b7

2 3 4 5 6 7

- 1 0 0 1 0 0

0 0 1 - 1 1 0 Ib 2
A Ib-
1 1 - 1 0 - 1 1 Ib 3 .

Ibi*

Ib 5

Ib 6

Ib7

n x m m x 1

Tenemos el producto de la Sera, fila de matriz incidencia por Ib

Fila 3 x Ib - [ O 1 1 -1 O -1 1 ] X Ibi
Ib 2
Ib3
Ib*
Ib s
Ib6
Ib7
 - 21 -

= Ib2 + Ib3 - Ibtj - Ib6 + Ib7

= (Ib2 + Ib3 + Ib7) - (Ib«> - Ib6)

Este producto representa:

Suma de corríe_n_ Suma de corrientes

tes de rama que - de rama que entran
sale de nodo 3 a nodo 3

De acuerdo con las leyes de Kirchoff la expresión anterior debe

ser cero entonces se tiene que para todos los nodos

[Link] = O • (1.15)

Reemplazando 1.8 en 1.15 tenemos:

A [ Yb (Vb + Vg) - Ig ] = O

A Yb Vb + A Yb Vg - A Ig

como Vb = A Vit

entonces

A Yb AT \¡n = A. (Ig.- [Link])

donde se define Yn = A Yb A matriz de admitancia nodal

Is = A (Ig - Yb Vg) matriz de corrientes equi
val entes.
 - 22 -

Yn Vn - Is (1-16)

Por consiguiente al determinar los voltajes de nodos de la red

Vn (resolviendo este sistema de ecuaciones simultaneas) también
se podra determinar los voltajes y corrientes de rama Vb, Ib res_
pectivamente.

1.5 EJEMPLOS DE APLICACIÓN

1.5.1 Circuito DC

[Link].

zov 20 SL

Para iniciar el análisis, el circuito debe estar completamente,

definido, por lo tanto a más de los datos del valor de los elementos que
lo componen se debe:

1. Seleccionar el nodo de referencia

2. A cada nodo del circuito se le asigna un número (1, 2, ...n )
n & de nodos
3. Se enumera las ramas (1, 2, . . . m)
4. Se asigna los sentidos de las corrientes, o de las ramas
 - 23 -

6 ramas,m = 6

3 nodos, n - 3

independ.

-(o)
Nodo referencia.

Matn'z de incidencia

I 2 3 4 5 6
1 0 -1 0 1 1 0~

A = 2 0 0 -1 -1 0 1
3 -1 1 1 0 0 0

Matriz de admitancia de elementos Matriz de Volt, inde-

pendiente.

1/5 0 0 0 0 0 20
0 1/10 0 0 0 0 O
0 0 1/6 0 0 0 O
Ye = Vg
0 0 0 1/2 0 0 O
0 0 0 0 1/20 0 O
0 0 0 0 0 1/4 O

C á l c u l o de matriz admitancia de rama

Yb - ( U + F ) Ye ( U + D )
 -24 -

pero F = O

D = O

entonces Yb - Ye

T
Cálculo de Yn - A Yb A

O -1 O 1 .1 O 1/5 0 0 0 0 0
O 0 - 1 - 1 0 1 0 1/10 0 0 0 0
A Yb =
- 1 1 1 0 0 0 0 0 1/6 0 0 0
0 0 0 1/2 0 0
0 0 0 0 1/20 0
ó 0 0 '0 0 1/4

• O -1/10 O 1/2 1/20' O

A Yb - O O -1/6 -1/2 O 1/4

-1/5 1/10 1/6 O ' O O

0 -1/10 0 1/2 1/20 0 0 0 -1

T
Yn = (A Yb) A 0 . 0 -1/6 -1/2 0 1/4 -1 0 1
-1/5 1/10 1/6 0 0 0 0 -1 1
1 -1 0
1 0 0
'o 1 0
_
 - 25 -

13 -1/2 t -l/10

20
Yn = -1/2 11/12 -1/6
-1/10 -1/6 7/15

Cálculg : Is = A (Ig - Yb Vg), pero ya que Ig = O

entonces:
Is - - A Yb Vg

0 -1/10 0 1/2. 1/20 0 20

Is- - 0 0 -1/6 -1/2 0 1/4 O

-1/5 1/10 1/6 0 0 0 O

O
O
O

0 0

Is- - 0 = 0

_4 4

Por lo tanto se tiene que:

 - 26 -

YÍI Mu = Is

13 1 -1/10 Vru 0
20 2

__! 11/12 -1/6 Vn 2 = 0

2

_JL -1/6 7/15 Vn 3 +4

10

Vni = 5,948

Vn 2 = 5,377

Vji 3 - 11,776

1.5.2 C i r c u i t o AC

bn
ti
1'
1f
T\lv.
•i-' 1 r
rb, b*

- 1 1 0 0 • o 10

A = 0 - 1 1 1 0 Ig = 0

0 0 0 - 1 1 0

0
 - 27 -

Yb = Ye -

0.5

AxYb -

- j 0.5 1/3 0

0 - 1 / 3

(A x Yb) AT = O O o -i 0 ó

0.5 1/3 O I -1 0

/I 1
0 1 0
0 -1/3 wj
4, 0 1 -1
0 0 1
 - 28 -

(A Yb)AT =

4 -1/3

I-s - A (Ig - Yb/Vg)

Is = A Ig = -1 1 0 0 0 10 -10
-1 1 1 0 0
0 o
0
0 0 0 -1 1 o
0
_o _

Para w = 2

1 - j J O Vil! -10
j 0,833-j -0,333 Vílo = 0
O -0,333 0,333 + 0.5J Vita 0

donde

- 6,06 - l,24j
4,83 - 2,71j

Vn 1,53 + 0,123
 - 29 -

Como se puede observar en estos ejemplos, particularmente en los

casos de corriente alterna (AC), los cálculos se vuelven engorrosos y
más aún si se debe hallar o resolver el circuito para varias frecuencias.
Por consiguiente la imperiosa necesidad de utilizar un computador digi -
tal, más aun cuando se tiene que analizar circuitos complejos que son -
los que tienen utilidad práctica.
 - 30 -

C A P I T U L O S E G U N D O

2.1 INTRODUCCIÓN

Como se expreso anteriormente, en el análisis nodal se hace uso

de los voltajes de nodos como variables auxiliares, para luego determinar

los voltajes y corrientes de rama.

Según las ecuaciones desarrolladas en el capítulo anterior la se_

cuencia del cálculo es la siguiente.

(1) YÍI - A Yb AT

(2) Is = A (Ig - Yb Vg)

(3) [Y)t I Is] - Vit

(4) Yb = AT Vn

(5) Ib = Yb [Vb + Vg] - Ig

La primera fase (1), el cálculo de Yn (Matriz de admitancia no 1 -

dal), se la puede dividir en tres subfases; la formación de las matrices
de incidencia y de-admitancia ramal Yb, el producto de estas A . Y b , y
 C A P I T U L O S E G U N D O

DESCRIPCIÓN DE ALGORITMOS EMPLEADOS Y SOLUCIÓN DE LAS ECUACIONES

QUE. DESCRIBEN EL CIRCUITO

2.1 Introducción
2.2 Estructuración de la matriz incidencia (A) y matriz de admitan-
cia ramal (Yb)
2.2.1 Matriz incidencia (Estructuración)
2.2.2 Estructuración de matriz de admitancia rama Yb
2.2.3 .Estructuración de Yb cuando existen inductancias mutuas
2.3 Producto ([Link]) AT - Yn
2.4 Solución numérica de ecuaciones lineales simultáneas
2.5 Inversión de matrices
2.5.1 Método de Gauss
2.5.2 Método de Shipley
 - 31 -

finalmente la posmultiplicación,por la matriz de incidencia transpuesta

(AT).

La segunda fase (2), consiste en la obtención de la matriz de co_

rrientes equivalentes Is. Se lo hace igualmente en tres pasos, primero,
'el producto [Link], seguido de la substracción de este producto de Ig y
la premultiplicación por A. De este modo se obtienen la matriz de coe-
ficientes (Y¡i) y el vector de términos independientes (Is) de la ecuación
matricial Yn Vn = Is que representa un sistema de n ecuaciones line^
les simultáneas de donde se obtienen los voltajes de nodos de la red.
Obtenidos estos voltajes r la primera parte del análisis esta completa.
Cuando se requieren las corrientes de rama previamente se debe obtener
los voltajes de rama.

Esta secuencia es indistintamente utilizada tanto para análisis

en corriente continua (d.c) como en corriente alterna (a.e.), tomando en
cuenta que para el análisis en (a.c.) las ecuaciones formales del circuí^
to serán:

(A Yb(w) (AT))Vn - A [Ig(w) - Yb(w) Vg(w)]

A continuación se explica en detalle "los algoritmos empleados en

la formación de las matrices de incidencia, de admitancia ramal, los prp_
ductos [Link], (A Yb).A , la reducción y solución del' sistema de ecuacio-
nes para obtener los voltajes de nodos, inversión de matrices, la métod_o_
logia utilizada en el almacenamiento de resultados cuando se analiza el
circuito para un rango de frecuencias.
 - 32 -

Los algoritmos que se han empleado, están en función de la me-

jor utilización de las localidades de memoria disponibles.

2.2. ESTRUCTURACIÓN DE LA MATRIZ DE INCIDENCIA (A) Y MATRIZ DE

ADMITANCIA RAMAL Yb

2.2.1 Matriz Incidencia

La matriz de incidencia en sí, es una matriz de interconexión,

que define como se conectan las ramas a los nodos independientes. Tie-
ne tantas filas como nodos independientes y tantas columnas como número
de ramas tenga la red por lo tanto su dimensión es (n x m).

Para la formación de la matriz incidencia se debe conocer la

descripción de cada rama (clase de elementos que la conforman), se debe
indicar también, para cada rama el nodo "desde" y el nodo "hasta".

El procedimiento usual para la formación de la matriz inciden -

cia es el siguiente:

!)• Se debe inicializar la matriz de incidencia con ceros.

2) Analizar uno por uno los diferentes elementos si es R, L, C,

G encontrar el número de rama en que se localiza, el núme
ro del nodo "desde" y el número de nodo "hasta".

3) En la columna indicada por el número de rama y en la fila

indicada por el número del nodo "desde" poner un 1 igno -
 - 33 -

rando si el nodo "desde" es cero.

4) En la columna indicada por el número de rama y en la fila

indicada por el número del nodo."hasta" poner -1 , ignorar
si el nodo "hasta" es cero.

Por lo tanto, en todo caso los elementos de la matriz inciden -.

cia A serán O, 1, -1. Caso contrario hay error en datos.

En el minicomputador TEXTRONIX 4051 todas las cantidades numérj_

cas son tratadas como números reales y un número es representado en ocho
bytes, de ocho bits cada uno; en cambio cada carácter literal se repre
senta en un byte. Lo que significa que una red que tenga veinte nodos
y treinta ramas, la-matriz incidencia tendrá seiscientos elementos y
ocupará cuatro mil ochocientos bytes. No así, si estos elementos son
tratados como caracteres alfa o literales ya que utilizan únicamente
seiscientos bytes esto implica un ahorro de memoria y más aún teniéndo-
se en cuenta que éstos, son un conjunto de 1, O, o' -1.

Por lo tanto la matriz de incidencia se formará en una variable

literal A$ y tiene n x tn caracteres, y los

1 2 K rnxm

FIG. 2.1.

caracteres a almacenarse serán "O", "1" ó "2" que corresponden a

 - 34 -

-1, O, 1, respectivamente. Como se observa la variable literal A$ t1_e

ne tantos caracteres como elementos tiene la matriz incidencia A, por
consiguiente a un elemento aij de A le corresponderá un carácter en
la posición K de A$. (Fig. 2.1), por ejemplo un circuito con tres np_
dos independientes y cuatro ramas.

n =3
m = 4
A (n x m) - A (3 x 4)

0-11 0-12 0-13 «-1 f

A = 2 0-21 0-22 0-23 «-2 f

0-31 «-32 «-33 0-3 if

0-11 0-12 0-13 O-l^ «-21 0-22 0-23 0-24 0-31 0-32 0-33 0-34
A$ =

K ->• 1 2 3 4 5 6 7 S 9 10 11 12

Posición en A$
 - 35 -

Al elemento a.23 le corresponde un carácter que se encuentra en

la séptima posición de la variable A$. Por lo tanto generalizando la
posición K se determinará.

K = (i- 1 )m +j

En la matriz A( 3 x 4 ) si: i = 3, j = 2, m == 4 entonces

K = 10 y corresponderá al elemento a32-

El procedimiento para la estructuración de la matriz de inciden-

cia, en la variable literal A$ es el siguiente:

1.- Se inicializa la variable literal .con caracteres "1".

2.- Según la clase de elementos, sea R, 0 , 1 , 6 : se colocaría

un carácter "2" en la posición K determinada por:

/
K = (nodo "desde"'- lUnúmeros de ramas \\o de la
\de la red J \rama

se colocará un carácter "O" (cero) en la posición K .

K = (nodo "hasta" - 1)/número de ramas\ número de la

,de la red / V rama

Esto significa que la matriz incidencia se ha transformado en

una sucesión de caracteres literales "1", "O", "2" y cualquier opera -
ción que se realice con A, como los productos [Link], (A Yb).A se
 - 36 -

simplifican a una simple decisión' según sea el carácter y a la adi-

ción del valor correspondiente.

El procedimiento para la formación de A se resume en el dia-

grama de flujo de la Fig. 2.3.

EJEMPLO DE APLICACIÓN

Dado el circuito Fig. 2.2 estructurar la matriz de incidencia

en caracteres literales según la metodología antes expuesta.

. b2 b3
-t>—AAAAA-

b5 b6
V

FIG. 2.2.
 - 37 -

Clase de Nodo Nodo

N^ elemento # Rama "desde" u.,d ii

1 R 1 0 1

2 R 2 1 2
3 C 5 2 0
4 L 3 2 3

5 C ' 4 3 0

6 Vg 1 - -

A$ tendrá 15 elementos inicial izados con "I 1

A$ = 111111111111-111

1. Se pondrá "O" en K (1-1) 5 + 1 = 1

2. Se pondrá "2" en K (1-1) 5 + 2 - 2
"O" en K (2-1) 5 + 2 = 7
3. Se pondrá "2" en K (2-1) 5.+ 5 =10
4. Se pondrá "2" en K (2-1) 5 + 3 - 8
"O" en K (3-1) 5 + 3 =13
5. Se pondrá "2" en K (3-1) 5 + 4 =14

A$ .= 0 2 1 1 1 1 0 2 1 2 1 1 0 2 1

Con fines de comprobación se formaría la matriz i n c i d e n c i a ,

 - 38 -

leer Ji = número de nodos

m = número de ramas
AS (n x w)

I n j c i a l U a r AS
con caracteres "1"

I = 1, TO MI

Leer la clase de elemento

Leer
SI Húmero de rama
nodo "desde"
nodo "hasta"

f)0

SI

NO

poner un "z 1 . en A5 en la posición

K = (n."desde"-l)*n +• $ de rama

poner un "0" en la posición

K
K = (n."hasta"~l)* n + í de rama

FIG. 2.3.
Diagrama de f l u j o de la estructuración de la matriz de i n c i d e n c i a
 - 39 -

normalmente-

±3 :

A = -1 1 O O O

0 - 1 1 0 1
0 0 - 1 1 0

Para cada elemento de A: -1 se reemplaza "O"

Ó se reemplaza "1"
1 se reemplaza "2"

A$ = 0 2 1 1 1 1 0 2 1 2 1 1 0 2 1 .

2.2.2 Estructuración de matriz de admitancia ramal Yb

Previamente a la obtención de la matriz de admitancia Yb, se de-

be obtener la matriz de admitancia de elementos Ye donde

Ye = [ Y¿j ] ([Link])
YA] - O

Y-i-i- - Ye - - admitancia de elemento en la ra_

ma A,

Generalizando:

Ye.. = + -j w CU i- - ([Link])
-ÍX.
w
 Nodo A
o

Rama i
r
Nodo &

FIG. 2.4.

Significa que elementos en paralelo se los puede considerar inte_

grantes de una misma rama que va, del nodo A, al nodo B.

Ye(U

donde se definió

F - B + G Ye

D = R Ye + M

La existencia de fuentes dependientes de corriente sean controla

das por voltaje o corriente determinan la presencia de F, igualmente -
las fuentes de voltaje controladas por corriente o voltaje definen la pr
sencia de D. El caso más general será cuando existan los cuatro tipos
de fuentes dependientes, por lo tanto

Yb - [ U + B + G Ye x] Ye [ U + R Ye + M ]
 •- 41 -'

jSi existen una clase, dos o tres los cálculos a realizarse disminu-
yen. Se considera para el análisis 21* = 16 posibilidades de combina -
clones de fuentes, que pueden existir en el circuito, determinando para
cada caso la ecuación de Yb a ser estructurada.

Estas combinaciones se indican en la Fig. 2.5 en donde, 1a

existencia o no, del tipo de fuente se representa por los números lógi -
eos 160 respectivamente.

Resumiendo, se ha definido un código binario que indica la pre-

sencia de las fuentes dependientes y el tipo de las mismas, éste sé alnia_
cena en un vector al que se lo transforma a su equivalente y constituye
un indicador para tomar la decisión de la clase de ecuación de Yb a ser
estructurada. (Fig. 2.6}.

Por ejemplo si:

CÓDIGO BINARIO D0= O O O O

INDICADOR D!= O

el indicador señala la ausencia de fuentes controladas. Por lo tanto:

Yb - Ye
 - 42 -

INDICADOR
CÓDIGO BINARIO ECUACIÓN Yb A ESTRUCTURARSE
DECIMAL

0 0 0 0 0 Ye

1 0 0 0 Í (U + B) Ye

2 0 0 1 0 Ye (U + R Ye)" 1

3 0 0 1 1 (U + B) Ye (U + R Y e ) " 1

4 0 1 0 0 Ye + G

5 0 1 0 1 Ye + G + B Ye

: 6
0 1 1 0 (Ye + G) (U + R Ye)" 1

7 0 1 1 1 (Ye + G + B Ye) (U + R Ye)' 1

8 1 0 0 0 Ye (U + M)" 1

9 1 0 0 0 (U + B) Ye (U + M)' 1

10 1 0 0 0 Ye (U + R Ye + M)'1

11 1 .0 1 1 (U + B) Ye (U + R Ye + M)" 1

12 1 1 0 0 (Ye + G) (U + M)" 1

13 1 1 0 1 (Ye + G + B Ye) (U + M)"1

14 1 1 1 0 (Ye + G) (U + R Ye + M)' 1

15 1 1 1 1 (Te + G + B Ye) (U + R Ye + M)

FIG. 2.5.

2.2.3 Estructuración de Yb cuando existen inductancias mutuas

Para el análisis de varias ramas con acoplamiento inductivo en-

 Inicial ización de DO,
variable que almacena
código binario

Formación del código binario

según existan fuentes controladas
de Volt, y/o corr.

Transformación del código binario a

un indicador decimal
. DI

Formación de matriz
Admitancia de"elementos Ye

Estructuración de Yb -
Matriz de Admitancia Ramal
según indicador decimal

F!G Z.6

Áe. Flujo 6e>naral Je \a "Estructuración

¿la la. M a t r í * ' d e AJvnil3ncia Ram¿») Yb

 44 -

tre sí, se debe obtener la matriz de impedanciá Ze. Por consiguiente 3

se podría considerar [Link] rama a elementos pasivos conectados en se

ríe, donde la impedancia equivalente de la rama será la suma de las imp_e_
dancias de los elementos. Esto puede dificultar o provocar ambigüedades
en la utilización del programa, pues se definió que una rama puede for-
marse de varios elementos pasivos en paralelo. Para mantener consisten-
cia con esta definición, se analizó el acoplamiento inductivo entre dos.
ramas y sus admitancias fueron obtenidas en las ecuaciones ([Link])y(1.12)
Las admitancias de las ramas p y s se encuentran en la diagonal princi-
pal en Ybpp, Yb44 estos elementos de la matriz Yb, al igual que Ybp4
y Yb¿p deben ser modificados según (1.11) y (1.12)

Formar Yb
considerando
Ye diagonal

Modificación de los elementos

Ybpp, Ybp¿, Y8¿¿, Yb¿p
según Ees. (1.11) y (1.12)

Matriz Yb
se ha estructurado

T(6 2.3-

Díaqrama de Flujo para haflcfv Vb cuando

exi'J
xie Mutuas
 - 45 -

2.3 PRODUCTO (A. Yb) A T = Yp.

El producto A Yb se define:
m
A . Y b = Y = [Yjtj] = £1 CUfe^bfej (2.3.1)
k - í

¿ = 1, 2, 3 . . . n

j = 1, 2, 3 . . . m

Y es una matriz de (n x m). Sin embargo A se encuentra como l_j_

teral en A$ con caracteres "O", "1", "2" que corresponden -1, O, +1.
Por lo tanto si :

fe = "O" i m p l i c a , que (cUfe) i/bfej = -

i = "1" implica que '(CUfe) ^bkj = O

QMi = "2" i m p l i c a que (cafe ) ^jbfej - Vbfej

La posición de CLtfe. en A$ está dada por;

K l = (¿ - 1)* m - fe (2.3.2)

Se define A como la matriz transpuesta de A

AT = [ ¿J ¿j] ^. - 1, . . . m
= 1 . - - n
 - 46 -

y el producto:
n
Y AT = Yn - [ 5n¿j] - £ ^Ufe ó! fej (2.3,3)
= 7

¿ = 1, . . . n
j = 1, . . . n

la posición de O-fej está dada por

-K 1 = (¿ - 1} . m + fe . (2.3.4)

La prernul tiplicación de Yb por A y la posmultiplicación de

Y por A se reduce a la ejecución de simples adiciones.

2.4 SOLUCIÓN NUMÉRICA DE ECUACIONES LINEALES SIMULTANEAS

Utilizando el análisis de nodo se han estructurado n ecuacio -

nes lineales simultáneas, de cuya solución se determinará el vector, de
voltaje de nodos (Vn).

Consideremos el problema de resolver conjuntos de ecuaciones si-

multáneas lineales expresadas en las formas equivalentes:

i + 0.12 X 2 + . . . + O-m X n =

£L 21 Xi + 0-22 X 2 + . . .

Q-Yi-2 y.2 + . . . +&nnXn = bn

 - 47 -

o bien»

Xi bi

X2 b2

(2-4.1)

Xn bn

o en su forma de representación compacta

Ax = b (2.4.2)

donde A = matriz real de orden n x n

x = vector incógnita de longitud n

b = vector independiente de longitud n

Existirá solución única para el sistema, si el d e t l A ) ? 5 O, signj_

fica que el rango de la matriz A es igual al rango de la matriz ampliada
( AÍ b).

Si b = O el sistema es homogéneo y la solución única es la tri-

vial * = 0. Analíticamente cuando b * O la solución está dada por:

x - A x b (2.4.3)
 - 48 -

El problema real consiste en el cálculo de la solución puesto

que, aun cuando u es comparativamente pequeña, y si se resuelve el s1s_
tema por algún método conocido (p.e. la regla de Cramer) y cada determi_
nante se expande por separado se necesitarán mucho más, que n 1 muí ti -
plicaciones y divisiones, razón, que su aplicación es numéricamente ine-
ficaz. A continuación se examina un método sencillo y, con todo, de
gran utilidad para encontrar x a partir de A y B.

PROCEDIMIENTO DE ELIMINACIÓN DE GAUSS

Esta eliminación es simplemente una versión sistemática de la

técnica algebraica de "eliminar o deshacerse de las incógnitas", y trans_
forma al sistema en otro equivalente a través de varias operaciones; Las
operaciones posibles, denominadas operaciones elementales, que se pueden
realizar con las filas son:

1. Sumar una fila a otra.

2. Multiplicar una fila por un escalar
3. Sumar a una fila el producto de otra fila por un escalar.

También, se puede intercambiar las filas sin alterar la solución

del sistema.

La resolución de las YL ecuaciones lineales simultáneas se pue-

de resumir en dos etapas;
 - 49 -

1. Reducción de la matriz de coeficientes a una matriz triangular

superior.

= o

2. Evaluación de las incógnitas

A continuación se describirá el procedimiento para la resolución

de las n ecuaciones lineales simultáneas.

1. El proceso se inicia con ¿ = 1

.2. Se debe buscar al elemento pivote, elemento de mayor valor

absoluto de. la columna ¿ - ésima, seVá buscado desde •

<X¿¿ hasta o~ói

3. Si el elemento pivote se encuentra en una posición diferente

a a¿¿, colocarlo en esta, realizando intercambio de filas

4. Introducir ceros en la columna ¿ a partir de d¿¿, éste se

exceptúa. Se consigue esto obteniéndose para cada fila K
( & * ¿ ) , el factor multiplicador F definido como:

(2.4.4)
¿ -í- elemento pivote

Este factor se m u l t i p l i c a por la f i l a Jiy desde la columna L

hasta n, para luego sumar a cada elemento de la f i l a fc-
 • - 50 -

5. Se incrementa el valor de t y se repite el proceso hasta

obtener la matriz triangular superior esto es L ~ n - 1.

6. Terminada la reducción de Gaus-s se procede a la evaluación re

gresiva-de las incógnitas es decir se obtiene primero Xn, lue_
go Xn - 1 . . . hasta Xi.

El diagrama de flujo de la eliminación de Gauss en la (Fig. 2.8)

Ilustraremos esto con un ejemplo:

- 2x2 - x 3 - 2 (2.4.5)
+ —X2 - 2x3 = O (2.4.6)
- 2
- 2x2 + Ilx3 = -2 (2.4.7)

4 -2 -2 Xi 2

-2 +17/2 -2 X2 = 0

-1 • -2 11 X3 -2

El elemento pivote de la primera columna es 4 y está en la posi-

ción (1,1) no necesitamos intercambio de filas. Añadiremos múltiplos de
la Ec. (2.4.5) a las ecuaciones (2.4.6) y (2.4.7) para eliminar Xi ' de
estas ecuaciones. Los múltiplos o factores apropiados son anotados con-
venientemente en la columna de la izquierda.'. Haciendo esto se tiene:
 - 51 -

4 -2 -1 "x7 2

0 — ~- X2 1
2 2 '

-1 0 ~- — X3
2 4 2

El pivote de la segunda columna, buscado a partir de la posición

(2,2) hacia abajo es 15 .

- 2x2 - x 3 = 2 (2.4.5)

_15_ (2.4.8)
2

-3
(2.4.9)
2

.Para eliminar X2 .de la ecuación (2.4.9). se debe añadir el pro -

ducto de _!_ por la ecuación (2.4.8) a la ecuación (2.4.9) y se tiene:
3

- 2x2 - x 3 = 2 (2.4.5)

15
~ (2.4.10)
2

119
X3 = - (2.4.11)
12
 - 52 -

"X!~
4 -2 -1 2

X2 1
2 2

uo un 119
-7
X3
12 6

Se ha obtenido la matriz triangular superior. La resolución aho_

ra es directa y se obtiene de la Ecuación (2.4,11) x 3 = -2
j , substitu -
o
yendo ésta en la Ecuación (2.4.10) se tiene xz = or y X], de la E.'cua.-
ción (2.4.5) es — .
85

44

85 -10

El escoger como pivote al elemento de mayor valor absoluto ase-

gura una mayor precisión en la solución.

La eliminación gaussiana puede ser llevado a cabo con pivota -

nriento completo, ésto es, en cada paso usamos como pivote al elemento
de magnitud más grande de la matriz en cuestión.

La práctica demuestra que el utilizar el pivotaje simple o el

elemento de mayor valor absoluto en la columna correspondiente de la ma_
triz reducida con sistemas de orden alto o pequeño se obtienen resulta -
dos satisfactorios.
 - 53 -

í = L Til u - 1

Búsqueda del elemento pivote

ABS (pivote)\I

Si es necesario intej^
cambio de filas.

/• 1, tt

Obtención del
Factor F Ec. (2.4.1)

Transfonnación de los elementos

de A y del vector b de la fe-Ssina
fila

Evaluación Regresiva
de incógnitas

FIG. 2.8,

Diagrama de f l u j o de evaluación de Ecuaciones Lineales Simultá-

neas (Reducción Gauss J o r d á n ) .
2.5 EVALUACIÓN DE LA INVERSA DE UNA MATRIZ

2.5.1 Método de Gauss

La eliminación gaussiana con pivoteamiento parcial puede ser

adaptada al problema de evaluar la inversa de una matriz dada. Para
ilustrar la idea tornemos A una matriz de 3 x 3, digamos:

A - 0-23

«-31 0-33

Y la matriz inversa de A es la matriz X donde:

Xl3

X21 X22 X23

X32

entonces

A .X - (2.5.1)
 - 55 -

I matriz identidad

O-l 1 O-i 2 O,i 3 Xn Xl2 Xl 3 1 0 0

O-21 O.22 O.23 X2X X22 X23 — 0 1 0 (2.5.2)
0.31 0.32 0.33 X 31 X 32 X33 0 0 1

Significa que se tiene tres sistemas de ecuaciones lineales si

multáneas con tres incógnitas.

ler. Sistema: O-l 1 O-i 2 O_i 3 Xn 1

0-21 0-22 0-23 X21 = 0 (2.5.3)

0-31 0.32 .0.33 X3l .0

2do. Sistema O-l l O-! 2 O_! 3 A]_ 2 0

0-21 O.22 0,23 A2 2 - 1 (2.5.4)

•0-31 ít.32 0,33 ^32 0

3er. Sistema O-l l 0-12 0-13 Xl3 0

0-21 O.22 O-23 ^23 = 0 (2.5.5)
O-3 l 0.3 2 0.3 3 X33 1
 - 56 -

Los tres sistemas tienen igual la matriz de coeficientes (matriz

A) por lo tanto se pueden resolver al mismo tiempo usando el proceso de
eliminación gaussiiana, con la consideración que los términos independie_n_
tes constituyen ahora una matriz de tercer orden.

«-11 «-12 «-1 3 1 0 0 1 0 0 Xn Xi2 X13

«-21 «-22 0-2 3 0 4 0 -J-

0 1 0 X21 X22 X2 3

«-31 «-32 «-3 3 0 0 1 0 0 1 XB i Xa 2 X33

- —
A I
FIG. 2.9. FIG. 2.10.

Utilizando el esquema inicial de la matriz ampliada Fig. 2.9.

se obtiene la matriz de incógnitas Fig. 2.10. donde:

A. X - I

A"1 A.X = A"1

X '= A"1 (2.5.6)

Se utiliza el mismo esquema para la inversión de una matriz de

orden( n x n ).

Este método tiene la particularidad que para obtener la solución

se ha utilizado como punto de partida dos matrices, la matriz a invertir
y la matriz identidad, en cuya localización terminado el proceso se ob-
tendrá la matriz inversa.
 - 57 -

2.5.2 Método de Shipley

Este método tiene la ventaja que para obtener la matriz inversa

A 1, únicamente requiere la matriz a invertirse; pues A'1 se obtiene en
la misma l o c a l i z a c i ó n de la matriz dada.

El método a seguirse es el siguiente:

1. Iniciamos con K = I
I
2. Reemplazamos al elemento u^ftft por

3. Multiplicamos en la ¿ - ésima f i l a ( ¿ * h) al elemento

por - G-ftfr, este nuevo CLcft será el factor para

transformar los restantes elementos de Va <L - ésima f i l a .

4. La transformación de estos elementos de la ¿ - ésima f i l a

- excepto CL¿/f se obtiene sumando a e'sta el producto del
factor por los elementos de la h- ésima r matemáticamente:

5. .Se. repite el procedimiento para las filas restantes por lo

tanto se debe incrementar en 1 el valor de ¿ y repetir
desde el numeral 3 hasta que ¿ ~ n

6. Terminado estas reducciones se multiplica los elementos

a fio por d#í? " j # ft
 - 58 -

7. Se Incrementa en 1 el valor de K y se repite desde 2. F1nalj_

za cuando K = n.

Obtensión del
elemento [Link]

Transformación
de los elementos
de la ¿-ésima fila

Transformación de
los elementos de
k-ésima fila excepto Qfefe.

FIG. 2.11.

Diagrama de flujo de inversión de matrices (Método de Shipley)

 - 59 -

Ejemplo Numérico (Método de Shipley)

Sea A- 2 3

4 5

Hallar A 1 (matriz inversa)

1. Iniciamos con J? =
1
1 3
1 2
3.11 2
4 5

2. Transformación de los elementos de la ¿ - ésima fila

Ji = 2 Ji * K
1

. «-21 = -("14 - - 2 —
2.
-2 5

éste será el factor para obtener

<x 22 = a a i . 0.12 + 0-22 = (-2) (3) + 5 = -

-2 -1
 - 60 -

3: Modificamos los elementos de la fr - ésima f i l a excepto

«-12 - «-12 • 3 —

Se repite de i g u a l forma para fr = 2

«-22 =
«-22

—> 0.12 - - = - (-1) í -\ - factor

= «-12 «-21 + «-11 - —

-5 _3
2 2

-2 -1
 - 61 -

= (-2) (-1) = 2

-A ti
2 2
= A"
i2 -1

Como se observa a través de operaciones elementales y ocupando

únicamente las localidades que ocupa una matriz r se ha obtenido, la matriz
inversa.
 C A P I T U L O T E R C E R O

3.1 Estructura general de los diagramas

3.2 Programa maestro CIRCUITO

3.3 Programa PR/INGRESO

3.4 Programas: PR/FREQ, PR/OPTIOC

3.5 Programa PR/CALCULO

3.6 Programas: PR/YE1, PR/YE2

3.7 Análisis para varias frecuencias

 - 62 -

C A P I T U L O T E R C E R O

3.1 ESTRUCTURA GENERAL DE LOS PROGRAMAS

El programa digital en sí, consta de un conjunto de programas

los que se ejecutan individualmente, están regidos por un programa maes-

tro, el que conserva en memoria únicamente al programa que interesa ser

procesado en un determinado momento, ésto permite una mejor utilización

del área de memoria para el análisis de circuitos eléctricos más "largos

y complicados.

En el minicomputador TEXTRONIX 4051 se dispone de 30K bytes de

memoria y el lenguage de programación que emplea es el BASIC, con la pa_r

ticularidad que las variables tanto numéricas como literales no tienen

independencia entre un programa y otro. Esto obliga a mantener consiste_n_

cia en la utilización de los nombres de las variables durante el desarro_
lio global del programa.

En la figura 3.1. se presenta la estructura general de los pro-

gramas, los que han sido organizados y almacenados en una biblioteca, y

se ejecutarán, según requiera la lógica del proceso.

Programa nuestro: CIRCUITO

en
co

FIG. 3.1.
Estructura General de los Programas.
 - 64 -

A continuación se describe la función, características, y los

nombres de las variables que se utilizan en cada programa con su respec-
tivo contenido.

3.2 PROGRAMA MAESTRO: CIRCUITO

El programa maestro CIRCUITO conduce automáticamente el funci_o

namiento de los demás programas de la biblioteca (Fig. 3.1).

Este programa se limita a;

1. Verificar si el programa seleccionado .automáticamente, ,o re -

querido por el usuario, se encuentra en la memoria.

2. Si está lo ejecuta.

3. Si no está, limpia la memoria, carga al programa deseado y lo

ejecuta,

Al iniciar el programa, se presentan dos opciones de trabajo,

analizar un nuevo circuito o el proceder a grafizar un problema que ha
sido analizado con anterioridad. Esto tiene como finalidad facilitar al
usuario la ejecución de circuitos grandes, que requieran realizarse en -
una o más sesiones.

El funcionamiento del programa maestro se resume en el diagrama

de flujo de la Fig. 3.2 . Las variables utilizadas y su contenido son:
 - 65 -

-NOMBRE , CONTIENE

D9 Número de la unidad de discos utili-

zados.
P9 Dirección de la unidad de impresión
01 Número del programa que se desea
ejecutar.
02 Número del programa que está en memo-
ria.
03 Indicador controla la secuencia de
ejecución en el programa.
PR/INGRESO

06 , Bandera en el programa PR/INGRESO

H7 Bandera que controla la impresión de
resultados parciales, en pantalla
con el siguiente significado:
H7 = O no hay impresión
K7 = 1 si hay impresión.

El programa maestro permanecerá en memoria en el transcurso de

toda la ejecución, de modo que se podrá cargar otro programa,según se re_
quiera o según el usuario lo desee, en este caso se mostrará en pantalla
el índice y posibilidades de redefinir el circuito y/o cambiar las condj_
ciones iniciales antes definidas.
 - 66 -

Mostrar en la pantalla:
1. índice
2. Analizar un nuevo problema
6. Grafizar

ÍNDICE
En pantalla
1. índice
2. Ingreso/Corrección de datos
3. Ingresar
4. Operaciones de solución
6. Grafizar

Esta el progra_
ma deseado en menoría?

Limpiar la nercaria y

cargar el programa deseado

i
Ejecutar el programa

FIG. 3.2.
Diagrama de flujo del funcionamiento del programa maestro
circuito.
 - 67 -

3.3 PROGRAMA: PR/INGRESO

PR/INGRESO es el programa inicial del proceso, permite introdu-

cir datos de un nuevo circuito, o 1eer información de un problema antes
almacenado.

Este programa admite las siguientes opciones:

1. Ingreso y almacenamiento de datos en archivo de disco.

2. Lectura de archivo de datos (DISCO) e impresión en pantalla.
3. Impresión en papel.
4. Corrección / Actualización de datos.
5. Fin de ingreso / verificación / corrección.

Como se ve las opciones 3, 4, 5 no podrán ejecutarse sin antes

haber realizado las opciones 1 o 2.

Con la primera opción se realiza el ingreso inicial de los datos

del circuito, éste debe estar completamente definido. Se requerirá del
usuario una identificación del problema con un máximo de 72 caracteres,
el número de las ramas existentes, número de nodos independientes, para
cada rama indicar el número del nodo "desde" y el número del nodo "hasta1
la clase y valor de los componentes. A cada componente se asignará un
número según el código indicado en la siguiente tabla (FI6. 3.3 ) y
 . . .C:O.M-P O.N.E.N-T.E . . . .
NUMERO
SÍMBOLO. . '. . . N . O M-B R.E

1 . . . . V g . : . Fuente. Indep. de Voltaje

2 . ' ' ' .'la' '

3 " ' . R" " " Resistencia

4 ' * " . . ' 11 . ' . ' " ' Iridüctáfióiá
5 ' ; ' : c .' ." . ; Capacidad . .
.6. .. . . . ' . G . : . . Conductancia
7 . . . : ; M . . . - Inductancia. Mutua
' 8 . F . C C ' C . Fuente de Corr. controlada por Corr.

9 T V C C • Fuente de Volt, controlada por Corr.

. 10 -. .F.C C.V Fuente de Corr. controlada por Volt.
.11. . . . F . V . - C . V . Fuente de Volt, controlada por Volt.

12 . F.I.N A L.I.Z A C I 0 N DE I N G R E S O

FIG. 3.3.
Código utilizado para la identificación de la clase de. elemento

de acuerdo a éste, se preguntará al usuario, el número de la rama en que

está situado el elemento, el número de las ramas acopladas inductivamen-
te, o las ramas controlada y controlante si se trata de fuentes de voltea
je y/o corriente dependientes.

Concluido el ingreso de la información referente a un componente,

la que ha sido condensada en memoria en un vector, se almacena en un ar-
 - 69 -

chivo del disco, para poder reutilizar estas localidades de memoria en

el ingreso de un nuevo componente. Al finalizar la entrada de toda la i_n_
formación del circuito, esta ha sido almacenada en el disco y estará dis-
ponible, cada vez que se requiera en el proceso de cálculo.

El almacenamiento de la información ha sido organizada de manera

que si se necesitan los datos de un componente permanecerán únicamente
e'stos en memoria, durante el tiempo necesario.

Las variables utilizadas son:

VARIABLE CONTIENE

M Numero de ramas del circuito

N , Número de nodos independientes
NI ..,. Numero de componentes
V Vector de dimensión 7 3 contiene toda
la información referente a un compo -
nente.
P9 ., Dirección de l a u n i d a d de impresión
D9 Unidad de disco u t i l i z a d a
 - 70 -

VARIABLE CONTENIDO

03 Controla la lógica de ejecución

04 Bandera en el ingreso inicial de da-
tos
06 Bandera que permite corrección de un
componente o el aumento de varios.
U5 Bandera indica si el archivo que se
lee es el deseado
D$ Contiene el nombre del archivo en el
que se almacenan los datos.
Fl, F2 'Contienen la frecuencia o las fre-
cuencias .inferior y superior para
un rango.
Si Fl '= O y F2 = O significa que
el análisis se realiza en corriente
continua.

En el diagrama de flujo de la FIG. 3.4. se explica el funciona-

miento de este programa.

3.4 . PROGRAMAS: PR/FREQ, PR/OPTIOC

Estos son programas adicionales en el ingreso de la información

PR/FREQ permite la elección de realizar el análisis en corrien-

te continua, o corriente alterna. Si se opta por la alternativa de co-
 íf<

Escoger opción:
1. ¡nnreso y almacenamiento do datos en 3. Impresión en papel
archivo de disco.
4. Corrección / Actualización de datos
2. Lectura de archivo de datos e impre- $> Fln dg ingreso/i;eri flcactón/Corrección
sión en pantal la.

2 3 ¡i .)
_XZI~ *
I j
Ingresar el nombre Ingresar el nombre impresión en
papel Transferir el control
del archivo del archivo r>^<(í)>jü-> <J
i
l^n al programa maestro
1 ^ 1 1 HO
Ingresar: Debe leer o ingre ¡no para que cargue al
- Identificación del Lectura de archivo sor datos (T) pragraiM Pfi/FREQ
problema 2Í de datos

Corrección o
actualización

Ingresar datos del

componente en V ( 7 )

SI I Grabar NI
{# componentes)

Ó
FIG. 3.4.

Diagrama de flujo del programa PR/ INGRESO

 - 72 -

rriente alterna, se puede realizar el análisis a una frecuencia o para

un rango de frecuencias.. Para un intervalo se debe indicar las frecue_n_
cias límites: inferior, superior y el número de valores de frecuencia.
Si la relación entre las frecuencias límites es superior a 10 automáti-
camente el intervalo de incremento es logarítmico, si es inferior a este
factor, se puede escoger entre éste, o el incremento lineal.

Cuando se realiza el análisis en (D.C.) o (A.C.) a una sola fre-

cuencia las soluciones posibles a obtenerse son: los voltajes de nodos
o, los voltajes de nodos y las corrientes de rama; por el contrario si
el análisis se realiza para un rango de frecuencias las posibles soluci_o_
nes a obtenerse son;

1. Algún voltaje de nodo Vn¿

2. Cualquier diferencia de voltajes de nodo
Vn/ - Vnj
-t-

3. Cualquier corriente de rama I tu..

4. Relación VM¿ / Vnj
5. Relación Vwx / Ibj
6. Relación Ib¿ / Ibj
7. Relación Ib¿ / Vnj

quedando estas alternativas a elección del usuario. Esto se realiza en

el programa PR/OPTIOC.

A estos programas se los ha dividido en dos,diferentes,para faci_

litar al usuario reanaliza'r un mismo circuito a 'varias frecuencias
 - 73 -

y/o obtener como solución diferentes funciones de salida según se desee.

Este funcionamiento se esquematiza en los diagramas de las figu-

ras 3.5 y 3.6. Las variables empleadas en estos programas son:

si
Leer las frecuencias
FI y n

Debe leer o ingresar datos frecuencias

tienen valor /-I

Transfiere control a
programa Maestro

Transferir el control a programa Transferir control a programa

Maestro para cargar Maestro para cargar

PR/CALCULO PR/OPTIOC

FIG. 3.5.'

Diagrama, de flujo de PR/FREQ.

 Debe leer Análisis se 1i
o ingresar 'mita a [Link]-
datos ,corr. cíe ramas

Vaya a prog.
maestro

Escoger opción de
cálculo

Transferir el con-
trol al prog. maestro
para cargar PR/CALCU
LO.

FIG. 3.6.

Diagrama de flujo de PR/OPTIOC.

 - 75 -

VARIABLE CONTIENE
HS Bandera: H8 = 1 sólo voltajes de
nodo H8 - 2 voltaje de nodo y co_
rrientes de rama.
N2 Número de valores de Frecuencia.
Z Z = 1 frecuencia lineal
Z = 2 frecuencia logarítmica
D Incremento o factor de frecuencia.
O Clase de solución deseada para un
rango de frecuencias.
15 Almacena el valor ¿ de la opción.
J5 Almacena el valor j de la opción.
Fl, F2 Contiene información sobre la clase
de análisis . D.C. o A.C.

En PR/FREQ se decide cargar PR/OPTIOC si el análisis es para

un rango de frecuencias y si analiza para una frecuencia se transfiere
el control al programa maestro para que cargue PR/CALCULO.

3.5 PROGRAMA PR/CALCULO

En este programa se inicia el proceso de cálculo, se estructura

la matriz de incidencia, la que es almacenada en el disco en un archivo
auxiliar de trabajo (aleatorio). Determina luego la clase de fuentes
controladas existentes, formando el código binario y su equivalente de-
cimal .
 - 76 -

La matriz de admitancia- de elementos Ye es una matriz diagonal

de orden ( m x m ), cuando no existen inductancias mutuas. Si el circuj_
to es puramente resistivo se trata de una matriz diagonal real, pudiéndo_
se tratar a esta matriz como un vector de dimensión m que contendrá -
los elementos diagonales de Ye, Por el contrario, si el circuito contie_
ne elementos 'L,C , la matriz de admitancia de elementos es una matriz
diagonal compleja, que igualmente puede almacenarse en un vector comple-
jo.

Estos dos casos pueden ser tratados paralelamente, para el prime_

ro se desarrollará un programa que operará con números reales, en cambio
en el segundo caso la matriz Ye no puede ser almacenada en un vector com_
piejo, ya que el sistema de computación utilizado no hay posiblidad de
definir variables o cantidades complejas, de manera que se puede conside_
rar Ye constituida por el conjunto de dos vectores reales, que contie-
ne el uno la parte real y el otro la parte imaginaria, tomando en cuenta
que se debe evaluarlos considerando la aritmética compleja.

La presencia de inductancias mutuas será tratado incluyéndolo en

el segundo caso, considerando las variaciones explicadas en el Cap. I y
Cap. II.

En consecuencia en este, se determina el nombre del programa .a

ser cargado y estos son PR/YE1 o PR/YE2
 - 77 -

c io

1. Estructuración de la matriz de inductancia

2. Almacenamiento de la matriz Incidencia

3. Determinación de la existencia y clase de

fuentes controladas.

El circuito es s o t o
resistivo

Transferir el control al

programa maestro para que

cargue al programa •

PR/VE1

Banderas F y FO = O

Transferir el control al
programa maestro para que

cargue al programa:

PÍÍ/YE2
Banderas F » 1

FIG. 3.7.
Diagrama de flujo del programa PR/CALCULO
 - 78 -

VARIABLES CONTIENE

A$ Los (n x m) elementos de la matriz de

incidencia
DO . Código binario
(es un vector de dimensión (4))
DI Indicador decimal de la existencia y
clase de las fuentes controladas
FO
F Indicadores para seleccionar los pro-
gramas PR/YE-lo PR/YE2
F3 Valor de frecuencia al cual se reali-
za el cálculo.
Fl < F3 < F2

3,6 PROGRAMA' PR/YE1 y PR/YE2

En estos programas se estructuran las ecuaciones de nodos descri-

tas en el Capitulo I, Puede hab-erse desarrollado un programa integro, -
que ejecute todas las fases indicadas en la sección 2.1 del Capítulo II,
es decir la estructura de Yb, luego AYb y asi sucesivamente, pero gui_a_
dos por eV criterio de la optimización de la memoria se ha fragmentado en
programas cortos que realizan una función específica, así PR/YE1 se en-
carga de la estructuración de la matriz de admitancia de ramas Yb, y caj^
ga automáticamente al programa PR/YE1INC1D -, en éste se obtiene la matriz
de admitancia de nodos Yn. Este programa carga a su vez a PR/[Link]
donde se obtiene Is y los 'voltajes de nodos Vn, utilizando eliminación
P rcg rair-a Programa
Programa
PR/VE1 PR/YE11NCIO
PR/YE1SOLU

Estructura Yb (reales) Estructura Yti (reales)

Estructura Is^calcula
VH(Véales)

ls es complejo
Programa Programa
PR/ALKS PR/ESCRIBE .
Calcula Ib s1
se. desea

Prograna Programa
Programa
PR/YE2 PR/YE3IKCID
PR/YE3SOLU

Estructura Vb{ccí[Link]) Estructura Yb(complejos)

Estructur I s j
Vü (complejos)

FI6. 3.8.
Esquema de funcionamiento de los Programas PR/YE1 y PR/YE2 a una frecuencia
 - 80 -

de Gauss. Igual ocurre con PR/YE2 con las variantes, que este, tiene
que trabajar con números complejos, carga a los programas PR/YE31NC1D ,
PR/YE3SOLU, que desempeñan funciones similares que en el caso anterior.
Se debe anotar que al tener un circuito puramente resistivo, que conten-
ga fuentes con desfasaje, Is será un vector complejo, y se hará la trans_
ferencia al programa PR/YE350LU. Fig. (3.8).

Los capacitores e inductores son elementos que para el análisis

en (D.C.) pueden ser considerados como resistencias de 10 inega-ohmlos
y 0.10 ohmios respectivamente. En este programa se ha optado por tomar •
una w = 27rí suficientemente pequeña w = 10 10.

Los nombres de las variables utilizados en estos programas man -

tienen una relación correspondiente, como se ve en el cuadro de la
Fig. 3.9.

Debe anotarse que cada vez que se ejecuta un programa se tiene

en memoria sólo los datos que son necesarios para su procesamiento, alma_
cenando y leyéndolos de un archivo común-de trabajo.
 V A R I A B L E S rARürTFPT^TTPñ DF C O N T E N I D O
PR/YE1 PR/YE2
LA VARIABLE D E L A S V A R I A B L E S
R e a l I m a g.

El El E2 Vectores dim (M) elementos diagonales Ye

Gl Gl G2 Matrices dim (M,M) matriz de admitancia ramal ( Y b )

Yl * Yl * y2 ** Matrices dim (M,M) producto [Link]

y * YO * Matrices dim (N,N) matriz de admitancia nodal Yit
y ** I
CO
Cl + Cl ++ C2 ++ Vectores dim (M) términos independientes Is.

terminado evaluación de Gauss contienen volta

jes de nodos

E + E ++ EO ++
Vectores dim (M) contiene fuentes de Volt./corr- ind.
C C CO

C3 C3 C4 Vectores dim (M) contiene vector de corrientes de rama Ib

* PR/YE1 INCID ** PR/YE3INCID

+ PR/YE1 SOLU ++ PR/YE3 SOLU
PR/ALMS
FIG. 3.9
 - 82 -

3.7 ANÁLISIS PARA VARIAS FRECUENCIAS

Es Interesante observar que los caminos que en el principio se -

bifurcaron siguiendo caminos diferentes de ejecución, dentro de las dos j
trayectorias principales (PR/YE1, PR/YE2), convergen en el programa
PR/ALMS (Fig. 3.8). .Este programa toma varias decisiones; calcula las
corrientes de rama [Link] requerido y carga al programa PR/ESCRIBE si el
análisis es para una sola frecuencia.

Cuando el análisis es para un rango o intervalo de frecuencias ,•

el almacenar para cada frecuencia todos los voltajes de nodos y corrien-
tes de rama es prácticamente imposible, por lo que se almacena únicamen-
te la función de salida deseada y la frecuencia a "13 que se realiza el •
cálculo, posteriormente se debe ejecutar todo el procedimiento anterior
para la siguiente frecuencia. Se observa que existen vectores y matri -
ees que no requieren nuevamente ser estructurados, por ejemplo, en el
vector de admitancia de elementos la parte real es independiente de la
frecuencia y sólo varia la parte imaginaria, ésta es la razón por la que
se organizó programas adicionales que. conducen a un mejoramiento del tieni
po de ejecución.

Este funcionamiento se esquematiza en el diagrama de flujo de la

Fig. 3.10.
 - 83 -
I N I C I O

In1c1al1zac1<5n:
1. De la frecuencia
F3 - Fl
2. Del contador de I de fre
cuencias
II - 1

C á l c u l o de Vb en
la frecuencia F3.

Cálculo de "in

Cargar PR/YE3SOLU

Cálculo Is

Evaluación Vn
Gauss Jordán

Cálculo del valor real/imag.

de la opción de solución
y almacenamiento

Obtención de la
frecuencia de operación
F3

FIG. 3.10.

Diagrama de flujo del Procesamiento para varias frecuencias

 - 84 -

Terminado el cálculo para varias frecuencias el paso siguiente es

la grafización de los resultados. El programa de GRÁFICOS comienza grafi_
zando los puntos obtenidos en el análisis; si el gráfico y las necesida -
des lo requieren se puede realizar una interpolación con un número de pun_
tos mayor, al número de puntos evaluados. Esto constituye una gran venta,
ja, ya que en circuitos grandes es posible realizar el análisis para un
número pequeño, razonable, de frecuencias y a través de una interpolación
obtener la respuesta de frecuencia para dicha función, ahorrando asi,
tiempo de ejecución, en el análisis del circuito.

La técnica de interpolación utilizada es la seccional cúbica y

consiste en encontrar un polinomio de tercer grado para cada intervalo de
dos puntos, con continuidad de la primera y segunda derivada en los pun -
tos de quiebre.
 - 85 -

Escoger clasa de nrJflco

1. Parte ííeal
2. Parte Imaginarla
3. Magnitud
4. Magnitud db
5. Magnitud db con referen-
cia
6. Ángulo (Grados)

Ingresar parámetros
para el gráfico
Escalas
Dispositivo de grafiza-
ción. etc.

Presentar gráfico de
los puntos calculados

SI

Real Jzar
interpolación
HO

1
Gráfico de
puntos con
interpolación

Gráfico de ejes
y/o reticulado

Impresión de
leyendas

F I H

FI6. 3.11.

Diagrama de flujo de la grafización

 C A P I T U L O C U A R T O

1. Ejemplo de aplicación

2. Conclusiones y recomendaciones
 - 86

C A P I T U L O C U A R T O

EJEMPLO DE APLICACIÓN

Como aplicación de la utilización se presenta el segundo circuito

equalizador de audio.

41

ImV ¿4100?
500K

6BK

50-w.F

FIG. 4.1.
 - 87 -

La función de Interés es obtener la respuesta de frecuencia entre

los 10 Hz y los 50 KHz, de la ganancia de voltaje en "decibeles. Esto es
la relación de V2 Para Ingresar los datos es necesario tener pre -
VI dB
viamente definido el circuito equivalente.

Vi

-t- bl
< 47 K 100 V

-ION 500K
-1000

A b6

V3

FIG. 4.2.

A continuación se presentan:

1. Listado de los datos del circuito TABLA 4,1.

2. Gráfico de la respuesta de frecuencia de V2 , Fig. 4 . 3 .

VI dB
3. Igualmente ésta se puede referir a un valor en-una determinada
frecuencia (f - 1000 Hz) Fig. 4.4.
4. Se puede también obtener los voltajes de nodos y corrientes de
rama. En este ejemplo se ha realizado para la frecuencia de
ESCUELA POLITÉCNICA NACIONAL
FACULTAD DE INGENIERÍA ELÉCTRICA
o
ANÁLISIS HE CIRCUITOS ELÉCTRICOS o
o

N
NOMBRE DEL PROBLEMA; CIRCUITO ECUALIZADOR (EJEMPLO HE LIBRO DE TEXTO)

ARCHIVO! ND3 í-t-

NO* DE RAMAS: 9 o.
o
NO, DE NODOS INBEPENB, ; 5
N O » DE ELEMENTOS? 11
CD

NO, CLASE RAMA DESDE HASTA VALOR FACTOR FASE

OD
~."í ro CO
1 FUENTE VOLT» 1 0 0 1*000 0,000
2 RESISTENCIA 1 0 1 47,000 ó
3 RESISTENCIA 9 1 7 10,000 6
4 RESISTENCIA 3 3 4 500,000 0
5 CAPACIDAD 4 4 0 50,000 ~A
6 RESISTENCIA 5 0 o
A-, 100,000 0
7 FVCV 0 o
A" 5 100000,000 0
'-j ~!
& RESISTENCIA ó 500,000 3
O
9 CAPACIDAD 7 5 4700,000 -12 .
10 RESISTENCIA Q 5 3 08,000 3
11 CAPACIDAD 9 o ¿ 1000,000 -12

T A B L f l ¿1.1
NOMBRE DEL PROBLEMA: CIRCUITO ECUALIZADGR (EJEMPLO DE LIBRO BE TEXTO)
VOLTAJES HE NODO Cf = 1000 Hz>

NODO MAGNITUD/VOLTIOS ÁNGULO/GRADOS

0*00100000 0*000
O»109¿S901 -40,732
0,00100034 -0,041
0,00000037 ~89,Ó7Ó
0*09347942 -14,450

co
CORRIENTES DE RAMA LO

RAMA MAGNITUD/AMPERIO; ÁNGULO/GRADOS

1 0 ,00000000 139, 337

15
0 ,00000000 O w/
139, -r77
T
0 .00000200 0,324
4 0 ,00000200 0,324
5 0,00000200 0, 324
ó 0 * 00000022 -41 , 075
7 0 ,00000143 -14,002
S 0 ,00000143 -14. 002
9 0,00000008 48, 925

TABLfl
 MAGNITUD EN dB DE
rvj no co w 4^ oí en -j
en G> en es tn co en C£> en en en

n
H
o
c:
H
-H
o
Zl

es
r~
H

m O
o XJ
c: ͙1
m n
'z. c_
o rn
H
ITJ
jr"
GD
:o
ni
H

m
H
rn
x
H
O
<S> ••_-'

^«

O)

G)
CS)
 ¿.V
CIRCUITO ECUALIZADOR CEJEMPLO DE LIBRO DE TEXTO)
r\
— -— -^
X,

1S X,
* \
G>
G> V
18
H \
m
-o c5
O X
w - •^ --
0
/N ^\
xy
>
X -5°
r\0 \
^< \ 18 180
>-te
Ld
Q
cn-15
"u
'

S -20
Q
D

g-25
C£>
<
SI
-38
c
FRECUENCIA / Hz
FIG. 4.4.
 CONCLUSIONES Y RECOMENDACIONES

Los resultados que pueden obtenerse del programa de la presente

tesis se aplican fundamentalmente al diseño y análisis de circuitos -
eléctricos, permitiendo determinar su comportamiento (respuesta de fre_
cuencia), conocer impedancias de entrada, de salida, alguna admitancia
o impedancia de transferencia etc., sin necesidad de implementar físi-
camente el circuito.

La efectividad y precisión del programa se ha verificado con la

utilización de .ejemplos simples, los que fueron .analizados previamente
por medio de procedimientos convencionales. Al menos se ha procesado
un ejemplo por cada una de las posibilidades de bifurcación que se tie_
ne, y los resultados obtenidos son totalmente satisfactorios y confia-
bles.

El desarrollo del trabajo se lo hizo tomando en consideración

las limitaciones del sistema de computación y optimizando lo disporii -
ble. El señalar un criterio específico en cuanto a la capacidad del
.programa para procesar circuitos extremadamente grandes, de un número
elevado de ramas y nodos es relativo, igual sucede con el tiempo de eje-
cución, ya que son dieciseis posibilidades de cálculo por cada bifurca -
ción principal.

Se sugiere que para el análisis de circuitos que no puedan ser

procesados en .el .computador dado que su magnitud está limitada por la ca_
gacidad_de memoria del equipo de computación,, analizarlos por partes o a
su vez se -recomienda para trabajos posteriores desarrollar técnicas 'de
-programación con acceso al disco para el almacenamiento matricial. "Esto
s&r-ía a ;eqsta de la [Link] un mayor tiempo de ejecución.

_Es de esperars.e que en el futuro puedan desarrollarse-e implemen_

•tarse programas de análtsis de circuitos no lineales, que constituirán al
igual que .éste, un instrumento didáctico a disposición del maestro y el
estudiante de la Facultad -de Ingeniería £léctric'a.

En el anexo 2 se indica de una manera explícita, en base a un

[Link], la manera de utilizar el programa.
 A N E X-0 1

LISTADO DE PROGRAMAS
 ANEXO 1 CIRCUITO PAG* 1
o !0*= " *
3 Ü9=~l
4 GO TO 100
8 01 = 1
9 03=0
10 06=0
11 GO TO 1000
12 01=2
13 GO TO 1000
14 GO TO 920
16
17 GO TO 1000
24
2 5 GO TO 1000
40 01=6
41 GO TO 1000
100 REM CIRCUITO
110 REM
111 REM ANÁLISIS DE CIRCUITOS ELÉCTRICOS LINEALES EN EL DOMINIO
112 REM DE LA FRECUENCIArPARA EXITACIOHES SINUSOIDALES EN ESTADO
113 REM ESTABLE
114 REM
115 REM ESCUELA POLITÉCNICA NACIONAL
120 REM TESIS DE GRADO DE ELIZABETH RODRÍGUEZ SALAZAR
121 REM
130 REM 5 DE MARZO DE 1982
140 REM
.150 PRINT LANALIS1S DE CIRCUITOS ELÉCTRICOS
151 IF D9< -1 THEN ISO
154 PRINT TECLA 1 - ÍNDICE DE PROGRAMAS"
156 PRINT TECLA ANÁLISIS DE UN PROBLEMA NUEVO"
158 PRINT TECLA ó GRÁFICOS DE UN PROBLEMA YA CALCULADO
160. G'O TO 10
180 PRINT TECLA ÍNDICE DE PROGRAMAS"
190 PRINT TECLA INGRESO/VERIFICACIÓN/CORRECCIÓN Y/O IHPRES
191 PRINT ION, DE DATOS DEL CIRCUITO
200 PRINT TECLA INGRESO O CAMBIO DE RECUENCIA"
201 PRINT TECLA OPCIONES DE SOLUCIÓN
203 PRINT TECLA GRÁFICOS"
PRINT ESCOJA TECLAS"
DELETE
X$(140)

INEA 800 , ' SELECCIÓN DE NOMBRE DE PROGRAMA

THEN 1030
-1 THEN 850
H7-1 IMPRESIÓN DE RESULTADOS INTERMEDIOS
H7=0 NO HAY IMPRESIÓN

IUNIDAD DONDE ESTA EL DISCO!

D9
IF NOT<D9=0 OR D9=l) THEN 820
CALL "UNIT'í09
DATA "INGRESO" FREQ" f "OPTIOC CALCULO YE! YE2" > "ALMS'y
S¿>2 ItATA " G R Á F I C O S
870 K'iZBTORE 860
nC'> t::'OR J=l TU 01
$90 fri^AD X$
' -

•••'^r/:t I'CL ETC 1001^20000

ÍV'/ÁÍ; 4 P p E ;\n XsilOOO
;j. O •"•••"' tt£M #&'* L ÍNEA 1.000 'SELECCIÓN AUTOMÁTICA DE PROGRAMAS
;!_ i) ". .\O
".O.?0 (30 TÜ 300
:!. 030 REM
': :"!4:^ ::LND
 ANEXO 1 SPR/INGRESO PAG, 3

REM ### LINEA 1000 SPR/INGRESO (TECLA 2)

G2~l
GO TO 800
REM
Y$=°INGRESO/ VERIFICACIÓN Y CORRECCIÓN DE DATOS"
GLOSE
PAGE
PRINT USING 'SXjFA'íY*
PRINT 'JJJJ 1 ~ INGRESO Y ALMACENAMIENTO DE DATOS EN•í
PRINT * ARCHIVO DE DISCO1
PRINT B J 2 - LECTURA DE ARCHIVO DE DATOS (DISCO) E IM 1 ?
PRINT "PRESIÓN EN PANTALLA'
IF 03=0 THEN 1160
PRINT "J 3 - IMPRESIÓN EN PAPEL"
PRINT 'J 4 - CORRECCIÓN/ACTUALIZACIÓN DE DATOS'
PRINT "J 5 - FIN DE INGRESO/VERIFICACION/CORRECCION•
PRINT "JJJJ ESCOGER OPCIÓN 'í
INPUT K
GO TO K OF 1200?3030?3970?4100?5420
GO TO 1160
PRINT "LINGRESO DE DATOS"
PRINT 'JNOMBRE DEL ARCHIVO PARA ALMACENAR LOS DATOS CO=FIN>:
INPUT C*
IF C$="0" THEN 1000

CALL 'FILE"-D9 ? D$ ? X$
Xp x$="* THEN 1320
PRINT VJYA EXISTE EL ARCHIVO a ?C$?' (BIBLIOTECA DAT)'
PRINT "JDESEA DESTRUIR SU CONTENIDO? (SI O NO)i G'í
INPUT X*
IF NQTÍX$="S" OR X*='SI") THEN 1200 '
KILL D$
PRINT "JJIDENTIFICACION DEL PROBLEMA"
INPUT Z*
PRINT -J NO * DE RAMAS: "?
INPUT M
PRINT "JNO* DE NODOS INDEPENDIENTES I Ȓ
INPUT N
REM INICIALIZACION DE FRECUENCIA

*M>?72+5
I A

'VECTOR DE ALMACENAMIENTO TEMPORAL DE DATOS DEL

COMPONENTE
CLASE DE! COMPONENTE
N U M f D E RAMA EN LA QUE ES EL ELEMENTO
[Link] NODO DESDE O RAMAS CONTROLANTE , o ACOPLADA1
N U M > I J E NODO HASTA O RAMAS CONTROLADA o ACOPLADA 2
VALOR DEL COMPONENTE
 ANEXO 1 @F'R/INGRESO PAG

1540 REM y (6)

15 5 O REM y í 7) -
1560 DELETE y
:! 570 DIM y < 7 >
:!. 580 NI-O
:!.590 NI:~Nl-fl
.1.600 y-O
:!.óio ü^--0
1620 RCM DESDC AQUÍ INGRESO DE NI-AVO COMPONENTE C6 NÚMEROS)
i!. 6 3 O P RIM T n LIN O R E S C) í'1 E. D A T Ü S 0 E L O S C O H p O N E N T E S ir EL. CIR C UIT O J "
•i. £ 4 O I: A T A a F ü E N T E IN D E F' > D ii U O L.. T A J E "
;i. 650 D A T A "FUENTE INDEP* DE C O R R I E N T E "
16 6 O D A "í" A " P E SIS T E N CIA u
•1.670 IWA " I N D U C T A N C I A "
::.6f;?o DATA "CAPACIDAD"? -CONDUCTANCIA" ? U INDUCTAWCIA MUTUA'
:!. 6 9 O D A 7' A fl F U til i'-! T E D E C U R R * C O N T R O L.. A 0 A ' P O R C O R R * 9
1700 DATA "FUENTE HE yQLT* CONTROLADA POR CQRR.*
;;.710 DATA "FUENTE DE CQRR > CONTROLADA POR VOLT."
Ü.720 DATA "FUENTE 0E '-/'OLT» CONTROLADA POR VOLT,
::. 730 DATA "FIN DE INGRESO DE DATOS"
17-O REST Ü RE 164 O
1750 FOR J-l TQ 12
1760 READ X$
1 77 O P RIN T U SIN O " 5 0 u " n " F A n í J ? X$
1700 NEXT j
17«o DATA "PICO" f "NANO 3 , "MICRO" r 'MILI' ? " " j "KILO" í "MEG'A11 ? "GIBA"
1800 DATA "TERA U
1810 HOME
1820 PRINT a J *
1 B3 O R E S T O R E 17 9 O
1840 FOR J=-12 Tu 12 STEP 3
1B50 READ X$
Í.GÓO PRINT USINO a 55X M 5A ? 2X ? 4D " í X £ r J
1870 NEXT j
1880 PRINT V J J °
1890 PRINT US1NG VSX"'COMPONENTE NO * u * 3 D ? 3 X n " C L A S E ; a n ? S n : N l
1900 INPU"í' y (1 )
1910 PEM SE CONTROLARA ÜUE SE INGRESE DATOS SEGÚN LAS ESPECIFICACIONES
19 2 O IF V (1) < -• 6 T H E N 2 J. O O
1930 GO TQ y(l)~6 ÜF 2510 - 2630.2630?2630?2630»2750
1940 PRINT USING a 5 0 X " n K E R R O R í CLASE MAX - 1 2 G • D a í
1950 QO TO 1390
1960 PACE
1970 PRINT USING " 7X " a NO * DE. [Link] DESDE C MAX H " 3D ü n ) \B ? S ' : N
19SO INPUT yí3)
; 1.990 04-1
2000 IF V(3)=>0 AND y<3X=N THEN 2040
2010 PRINT USING "50X""KERRORí NODOS MAX - G ' " F H " í N
2020 GOSUB 2860
2030 GO TO 1970
2040 PRTNT USING " 7 X n n N O * DE NODO HASTA ( M A X u n 3 D " B ) í U D :'S B :N
2050 INPUT y(4)
2 O 6 O IF y C 4 ) ~ > O A N D y C 4 ) < -- N T H E N 21S O
2070 PRINT USING U 50X""KERROR; NODOS MAX - G'"FU'íN
 ANEXO 1 8PR/INGRESQ PAG, 3

PRINT USING "7X NO( ti £ RAMA í MAX " " 3D "" ) í ' ' , 3 n J M

I !-' y ( 2 5 :.:- O A N D y < 2 ) < = M T I--I E N 2 1 7 O

PRINT USING "50XH "KERRORÍ R A M A MAX = G^FD"
GOSUB 2860
60 TO 2100
i! i"" y<'15>2 THEN 1.970
D A T A " V O L T I O S " • " A K P E R I O S " -u O H M I O S * -:l H E N RÍOS
DATA a P ARABIOS" ? "MOHOS" :• "HEHR10S " . '
BATA "OHMIOS " ? "MOHOS
':;: E S T O R E 2180
FOR ]>1. TO V(l )
REAB X*
N£XT I
f- R t M T U S I N G " 2 8 X " " V A L O R í * H ? S " :
INPUT V(5)
p'rÜÍNT IJSING "26X" " PREF:'UOí " " S " í
ÍNPUT 2$
22^0
2"íOO rp 2$=" u THEN 2400
fc'ESTQRE 1790

I F Y $ -' Z * "í H E N 238 O

Ir J<12 THEN 2330
GO TO 2270

USING "50X B H K u "FA" ÍX*

1 ) >2 THEN 244O
INGRESE LA FASE <£N GRAÜOS)! " í

> 1 T !••! E N 2 4 7 O

GO TO
WRITE
FOR J =l TO 70
NEXT J
QO TO 1590
2 S 3.0 PRINT USING a 5 X n s N Q , R A M A 1 ACOPLADA? ( MAX H H 3DU " )? ? su; H
2'";^'O INPUT y (3)
2530 IF y( »0 AND VÍ3X-M THEN 2570
2540 F1 R 1 N T U S I N G B 5 O X " " K E R R t) R í R A H A S M A X « G a a F B n t M
2550 GüSUB 2860
2560 GO TO 2510
2570 PRINT USING " 5 X n n N O , R A M A 2 ACOPLADA < MAX " " 3D H " ) ; ;3at M
INPUT
IF y<
PRINT
GOSUB
2620 GÜ TO 2570
2 a- 3 O P R1N T U S1N G • 2 X " a N O * D E R A H A C O N T R O L A N T E ( M A X " " 3 D ' D ) í " " > S " ; ¡
:^40 INPUT V < 3 )
.'¿ASO IF V(3»0 AND V(3X=M THEN 2690
2660 PRINT USING u 50X" "KERR'OR í RAMAS MAX = G ' - F B ' í M
2¿>70 GÜSUB 2860
2680 GO TO 2630
vA9O p R IN T U 8'!.' N G " 3 X " ' N O * D E P. A M A C O N T R O L A D A ( M A X " " 31'' " ) í " " « S ' i M
2700 INPUT V <4>

2720 PRINT UB1NG a50X" °KERROR; R A M A S MAX = G a " F D a S M

2730 GOSUB 2S60
;27¿o GO TO 269O
27ÍÍO IF Q6O1 THEN 2SOO
2760 F-'Fíi'NT "J INGRESE CORRECTAMENTE LA CLASE DE ELEMENTO »
2770 PRINT "J PAPA CONTINUAR PRESIONE ÍRETURN)'?
2780 INPUT ;¿$
2790 GO TO 1000
2800 Nl^Nl-l
2810 URITE *l?3;M/NfNl
:¿Q20 IF 03 = 1 THEN 4620
2S30 03=1
2040 Y$='DATOS FUERON INGRESADOS0
28ÜO GO TO 1060
2860 IF 04=0 THEN 2900
2670 l' = N
2B80 IÍ-" NODOS-"
2890 GO TQ 2920
2900 I-"H
2910 I^^-RAMAS"
2920 PRINT USING 2930íl$íl
2930 IMAGE /"NO, M A X * DE u FA' ES"'FD"r ES CORRECTO? (SI O N O > I G
2940 INPUT XT>
2950 IF X$="S" OR X$™ 9 SI" THEN 3020
.2960 PRINT USING "12X' B NO* CORRECTO DE B . « F A " a ? D Í ? S w : i $
2970 INPUT I
29BO SO TO 04 OF 3010
2990 M-l
3000 GO TO 3020
3010 N=I
3020 RETURN
3030 REM*«* LECTURA HE ARCHIVO E IMPRESIÓN EN PANTALLA £X<£
3040 PRINT B L LECTURA E IMPRESIÓN EN PANTALLA "
3050 PPXNT "4 ARCHIVO CON DATOS DEL PROBLEMA DESEADO (0=FIN)t "•
3060 INPUT Cí>
JO/O IF C$"S0" THEN 1000
3080 GLOSE
3 O 9 O D $ ~ " 'í1 D A T / ' S C $
3100 DELETE V
3110 DIM VÍ7>
VL20 CALL u FILE " r D9 .• D* r X*
5Í30 IF X*='' THEN 3230
3140 P9=s32
3150 80SUB 3280
3160 00 Tu U5 OF 3170?3050
'317'O Q3--1
31^0. rii=" DATOS LEÍDOS DE ARCHIVO; " SC4
'.'•íl0'? '30 TO tJ6 OF 3200*3220
;>20C PtflNT "J PARA CONTINUAR APLASTE (RETURN) "
3210 INPUT X'f>
3220 GO TO 1060
3230 PSINT "J ARCHIVO? "?C$ÍB NO EXISTE*
3240 -'RINT ' JJ DESEA LEER OTRO ARCHIVO ? * í '
3250 l'WP'UT Xí>
32.::-> ir NOT(x-t>~= H s n OR x*="si") THEN 1040
3270' GO TO 3050
'.12x0 "¡rri*** SUBRUTINA DE LECTURA E IMPRESIÓN **#

33-0 *EM U5Ȓ ARCH, ES CORRECTO.- U5^2 CAMBIAR DE ARCHIVO

33!. O U ¿--I
332'"1 PRINT (2P9? USINO 3330;
3330 TMAGE P"ESCUELA POLITÉCNICA NACIONAL"
3340 --RINT £P9í 'FACULTAD DE INGENIERÍA ELÉCTRICA*
33 '•; O P RIN T @ P 9 í " J A N A L131S I' E CIR C UIT O S E L E C T RIC O S '
3 3 •.:':- O P RIN T @ P 9 í U SIN S 3 3 7 O í
3370 IHAGE 72C"-")
3 3 S O C A1- L,. " FIL E u ? D 9 .- D $ < X $ •
3 3 9 O IF PO 3(X $ P " O OP E N'?13 O)=O T HEN 3410
3400 OPEN Díy 1 r ' F " # X $
3410 READ *l?lí2$
3420 PRINT @P9í USING 3430tZÍ
:í430 II^AíSE /"NOMBRE DEL PROBLEMA; fl ? FA
V5 -¿ A o x í=S E G < Ü $ ? ¿ ? L E N ( D * ) - 5)
3450 PRINT SP9Í USINO 3400VX$
3 4 i- O IM A G E /15 X ' A R C HIV O í ° 3 5 T F A
3470 READ # 1 ? 3 Í M : N ? N I -
3 -A £ O ^ r-;: IN T © P' 9 i U SIN G 3 4 9 O í M
3490 IHAGE lOX'NO- DE RAMASí u 35T3D
3500 PRINT EP9Í US1NG 3510iM
3 51 O IM AGE 10X"NO, HE N O D OS IUDEPE N D,í"35 T3D
3 5 2 O !::' RIN T í¡? P 9 ? U SIN G 35 3 O ? N1
3530 1MA6E lOX'NO» DE ELEMENTOS I"35T3D
35^0 F'RINT f?P9í USINO 3370;
3550 1F P9O32 THEN 3070
3500 PRINT USING 3570í
3570 1MAC3E //5X"ES EL PROBLEMA DESEADO? (SI O NO)? G " * S
3580 INPUT X*
3 5 í:> O 1F N C? T í X $ ~n S " O R X $ -a SI " ) T H E N 3 9 5 O
3000 IF N1O10 THEN 3070
3610 PRINT USING 3620;
3620 IHAGE /5X"DESEA UER (PANTALLA) L03 VALORES DE LOS COMPONENTES? 'S
3630 P-RINT "(SI O NO) D ?
3640 INPUT X$
3650 T.F NOTCX*="S n . OR Xí=" SI" ) THEN 3930
3660 -PAGE
3 6 7 O L'i A T A " F U E N TE 'v' O L T * " t ' F U E N T E C O R R , u
3 6 8 O D A T A ' R E SIS T E N CIA a .- * IN D U C T A N CIA "
3 6 9 O D A T A " C A F1 A CID A D " ? ' C O N D U C T A N C JA'
 ANEXO 1 ePR/INGRESO PAG, 8
3 7 OO í A T A ' IN D U C * M U T U A " > " B E T A ( F C C C ) '
3710 DATA °FVCC 'r'FCCV
3720 DATA Q FVCV "
3730 PRINT GP9Í USING 3740?
'3740 IMAGE" NO, CLASE"¿X a RAMA'IX a DESDE"IX1 HASTA"4X« VALOR",s
3"oO s-'RINT @P9í USING 3760 í
3;>>C IMAGE IX" FACTOR'3X" FASE'
37 >'0 K'RINT £ P 9 í USING 3370 í
3/80 FOR £1-5 TO Nl-l-4
37^0 -:EAD -ti - II I V
'><;;•'O RESTÓLE 3670
38:0 FOR J«l TO i; < j. >
382:./ P EA D X $
3830 NEXT J
3^40 IF Vfi)=l OR VCl}-2 THEN 3880
3850 PRINT H?P'9í USVNG 3860 111-4 * XT> ? V ( 2) •• V < 3 ) t V < 4 ) r V C 5 > ? V < 6 )
386~' IMAGEN 3D2XFA?2X3i:!2X?2C3X3D2X) j7D»3D?5D
3870 30 TU 3900
3Í-Í80 P'RINT @ P 9 Í USING 3890 ÍI1-4 ? X$ ? V (.2 ) ? V (3 ) ? V C 4 ) • V < 5 ) .- V ( 6 ) • V (7 )
3890 IMAGE 3D2XFA s- 2X3D2X r 2 í 3X3D2X ) i 7D » 3D * SD , 9D » 3D
3t;;'00 [Link] II
3910 REM
3910 GQ TO 3960

3940 60 TO 3960
3950 U5-2
3960 RE TUPÍN
3 9 7 O ft £ M # * Sí IM P R E 310 N E N F' A P E L * £ #
3930 IF 03^1 THEN 4030
3990 PRINT ' DEBE INGRESAR DATOS O LEER ARCHIVO '?
4000 PRINT ' (RETURN PARA CÜNTINUAR)G"
4010 INPUT X$
4020 00 TO 1060
4030 P9-51
4040 PRINT 3 ALISTE EL IMPRESOR (RETURN PARA CONTINUARA"?
^050 INPUT X$
4060 GOSUB 3280
4070 Y$=s" DATOS DE ARCHIVO* * SC$
4080 Y$=Y$£ a IMPRESOS EN PAPEL"
4090 GO TO 1060
4100 REh***COn:RECCION DE HATOS***
4110 IF 03-1 THEN 4130
4L20 GO TU 3990
4130 REM
4140 Y*="CORRECCION -DE DATOS DE ARCHIVO: 'SC$
4ISO PAGE
¿ .1.6 O f F: IN T U SIN G " 5 X ? F A " í Y $
4170 PRINT "JJJJ i - CORRECCIÓN DE IDENTIFICACIÓN DEL PROBLEMA 0
4180 PRINT 'J 2 - CORRECCIÓN DE NO* RAMAS y NODOS; AÑADIR " ?
41.9 O P RIN T " C O M F' O N E N T E S °
4 2 O O F' RIN T " J 3 — °
.4210 PRINT "J 4 - CORRECIÓN DE DATOS DE LOS COMPONENTES'
4220 PRINT "J 5 - FIN DE CORRECCIÓN"
4230 PRINT "JJJJ ESCOJA OPCIÓN 'í
 ANEXO 1 gPR/INGRESQ P A G . 9

4240 INPUT Z
4230 GO TO 2 OF 42SO ? 4370 • 5410 ? 4650 r 5370 '
42-50 00 Tu 4230
¿2 'O tfEMttüí* MODIFICACIÓN DE IDENTIFICACIÓN ***
42 ?0 PRINT " J NOMBRE A C T U A L ES *?Z<&
42 >0 F'RINT " J J INGRESE NOMBRE C O R R E C T O ; • ;
4 3 "i O 1 N F' U T Z '£
4310 IF Z1K>" B

4'5¿-: r$"' IDENTIFICACIÓN FUE MODIFICADA ÍAROI-Uí

4350 y1í-"T''$£ " ) a
4360 GO Tu) 41.50
4370 fl'Eh*** MODIFICACIÓN DE M ? N? NI***
/i 3 80 O ó™ O
4 3 9 O P R I N T " J N O t RA H A B ~ a ? M ? ' N O * NODO S I N D » = ' ? N ? " N O * C O M P - «¡NI
4400 F'RINT " J DESEA CORREGIR NO* DE RAMAS (SI O NQ) Q ?
4410 INPUT Xí
4 4 2 O I F N 0 T < X $ = H S " O R X $ -:J S I * ) T H E N 4490
4 4 3 O P R I N T fl J I N G R E S E N O * DE R A MAS ' ?
4440 INPUT X
4430 IF X=>M THEN 44SO
44.60 PRINT aJ NO PUEDE DISMINUIR EL MUMERQ DE LAS RAMAS '
4470 GO TQ 4140
448O N~X
4 4 9 O P i:;: I N T a J n E SEA C O R R E G I R N O * DE- N O D O S (SI O N O ) " ?
4500 INPUT XS
4510 IF NOTÍX^'S" OR X^'SI") THEN 4580
4520 PRINT "J INGRESE N O » DE NODOS "•
4 5 30 I N P U T X
4540 IF X~>N THEN 4570
4550 PRINT -VJ NO PUEDE DISMINUIR EL NUMERO DE NQÜQSGG"
4560 GO TQ 4140 -
4570 N~X
4580 PRINT "J DESEA AÑADIR COMPONENTES (SI O NO) •?
4590 INPUT X$
4 6 O O I F N 0 T ( X * ~- " S a O R X * = " S I " > T H E N 2 310
4610 GO TO 1590 '
4620 Y*=" DATOS DE H ? N - N 1 FUERON CORREGIDOS < A R C H » Í " " SC$
4630 Y$~Y*SP ) "
4640 GO TO 4150
4650 REM*^:¿í MODIFICACIÓN DE DATOS 'DE COMPONENTES***
4660 '"'AGE
4670 PRINT " NO» DEL ELEMENTO QUE DESEA CORREGIR (MAX " Í N 1 ? " 0=FIN) "
4680 INPUT X
4690 IF X<=0 THEN 5340
4700 IF X>N1 THEN 4670
4710 F'RINT
4720 PRINT IJSING 4730:
4730 I MAOE " NO t 1 > CLASE u 6X H 2 ) RAMA " IX B 3 ) DESDE " IX " 4 ) HASTA ' 4X " 5 ) VALORfl ?
4740 PRINT USING 4750:
4750 1MAGE IX ' 6 ) FACTOR • 3X" 7 ) FASE '
4 7 6 O F' R I N T U S I N G 3 3 7 O t
4770 I'ELETE U
 A N E X O 1 G ? P R / I N G R E S O P A G * 10
v

47SO DIN y (7)

47^0 READ #l.-X4-4:y
*i«'.>0RESTORE 3670
•••M^í.O
FOR J--1 TO y < i )
46?O READ Xt
¿UíK1 Ni::.XT J
*u>'0 PRXNT USING 3890 \? X$? V (2) • y < 3 > .- V (4) ? V C 5 ) • V ( 6 ) ? U < 7 )
*- 'ú:: 'í O
F' F! IN T U S 3! N G 3 3 7 O *
48oO PRINT "J NCU DE LA COLUMNA A CORREGIR (MAX 7 0=FIN) 'í
•¿;:-/'0
1NPUT 1
-«íiííSO
¿F I-• O THEN 5340
QR I>7 THEN 4860
THEN 4940
'¿910 IF y<1)=1 OR y < 1 )=2 THEN 4940
4 « 2 O P RIN T " J J S O L O SE F' U E D E C O R R E G í R F A S E D E F U E N T E S IN D E P E N DIE N T E S
4^30 GO TQ 4860
4 9 4 O P F: IN T " J E L V A L O R A L M A C E' N A D O E S í " ? M C I )
49SO IF 1-1 THEN 4980
4 9 ó O F' RIN T " IN G R E S E V A L O R C O R R E C T O í " ?
4970 INPUT Vi!)
^980 GO TQ I OF 5000 ? 5060 ? 51.50 ? 5150 ? 5320? 5270 - 5320
4990 GO TO 4860
5000 Al-NI
5010 Nl-X
5020 06=1
5030 PRI "J INGRESE 'DATOS PARA ESTE COMPONENTE (PARA . CONTINUAR RETURN)"
5040 INPUT Z$
5050 GO TO 1600
5000 IF y( 1X7 THEN 5120
5070 IF y<I)~0 THEN 5110
5 O 8 O F' RIN T H J E L y A L O P . IN G R E S A D O DEBE SER < C E R O ) . - ?
5090 PRINT ' ES UN COMPONENTE DE ACOPLAMIENTO '
5100 G O• TO 4 9 6 O
5110 GO TO 4660
5120 IF yci)=:>l AND y<IK=H THEN 5320
5130 PRINT USING D 50XH n KERROR RAMAS HAX = G""F;D"íM
5140 GÜ TQ 4960
5.1 SO I F' y < 1 ) > 2 T H E N 5200
5160 PRINT vjj ERROR, pftRA FUENTES IND t NO SE NECESITA ACOPLAR RAMASG"
5170 PRINT "J PARA CONTINUAR PRESIONE CRETURN)"?
5180 INPUT Z$
5190 GO TO 4650
5200 • IF y<lX7 THEN 5240
5210 IF yCDOO AND y (I ><=M THEN 5320
5220 PRINT USING •" 50Xn ' KERROR R A M A S MAX = G"FD'íM
5230 GO TO 4960
5240 IF ya)~>0 AND y (IX-N THEN 5320
5250 PRINT USING a50X " KERROR: NODOS MftX = G""FD'iN
5260 GO TO 4960
5270 A~y(I)/3
5 2 S O Q=y CI) ••-1N T ( A ) * 3
5290 IF Q=0 THEN 5320
5300 PRINT "JJ ERROR G'
5310 GO TO 4960
 ANEXO 1 (2PR/INGRESG PAG» 11
5320 WRITE *1
5330 GO TO 4800
5340 Y*=' VALORES CORREGIDOS (ARCHIVOí "£C$
5350 Y$-Y$S')'
5300 GO TO 4150
5370 REM
5330 Y$=" FIN DE CORRECCIÓN (ARCHIVO: 'SC$
5390 Y$=Y$&')'
5400 GO TO 1060
5410 GO TO 4140
5420 REM CARGAR PROGRAMA INGRESO DE FRECUENCIA
5430 01=2
5440'GO TO SOO
 ANEXO 1 (2FR/FREQ PAG. 12

REM ********** LINEA 1000 @PR/FREQ *********

02 = 2
GO TU 800
REM INGRESO DE FRECUENCIA /CORRECCIÓN O CAMBIO HE FRECUENCIA
IF 03- 1 THEN 1090
PRINT a L DEBE PRIMERO INGRESAR O LEER DATOS DEL CIRCUITO"
PRINT "J PARA CONTINUAR APLASTE (RETURN) '?
INPUT X*
GO TO 1000
REM
READ # l í 4 í F l ? F 2
IF F1O-1 THEN 1440
REM INGRESO DE FRECUENCIAS

F2=0
PRINT "LJJJJ DESEA ANALIZAR EL PROBLEMA EN ¡ '
PRINT *JJ i,- CORRIENTE CONTINUA CEU C + > B
PRINT 'J 2*-- CORRIENTE ALTERNA C A . C O "
PRINT "JJ ESCOJA OPCIÓN " í
INPUT K
O O T O K O F 1 4 2 Of 1 2 2 O
GO TO 1190
PRINT °L ANÁLISIS EN Í A » C * > PARA í "
PRINT "JJ :L*- UNA FRECUENCIA'
PRINT B J 2»-- RANGO DE FRECUENCIAS"
PRINT "JJ ESCOJA OPCIÓN; "?
INPUT K
GO TO K OF 1390 y 1290
GO TO 1250
PAGE
PRINT "J INGRESE FRECUENCIA INFERIOR .EN < Hz * > " í
INPUT Fl
IF F1OO THEN 1300
PRINT "j EN EL ANÁLISIS PARA UN RANGO HE FRECUENCIAS- NO SE CON n
PRINT 'SIDERA'f'J EL ANÁLISIS A CORRIENTE CONTINUA**
GO TO 1300
PRINT "J INGRESE FRECUENCIA SUPERIOR EN (Hz*) 'í
INPUT F2
GO TO 1420
PRINT '!_ INGRESE LA FRECUENCIA EN <HZ ( ) '?
INPUT Fl
F2 = F1
REM
WRITE *1,4ÍF1 :F2
REM AERIFICACIÓN / CAMBIO O CORRECCIÓN DE FRECUENCIA
IF F1OF2 THEN 1660
IF Fl=0 THEN 17 SO
PRINT "LJJ EL ANÁLISIS ES PARA í ."?Flf" Hs*"
PRINT "J ES CORRECTO (SI O NO) "?
INPUT X$
IF X*='SI" OR X$-'S° THEN 1790
•X*=' 1 CAMBIO O CORRECCIÓN DE FRECUENCIA"
Y$="ANALISIS PARA UN RANGO DE FRECUENCIAS'
PRINT "jj DESEA; *
 ANEXO 1 GíPR/FREQ P A G . 1

1540 PRINT 'JJ 1 — ANÁLISIS EN C P . C . ) "

PRINT 'J 2 — ?X«
PRINT 'J 3 —
:!.':< 7 O PRINT "JJ ESCOJA OPCIÓN:
INPUT K
r^o GO TO K OF 1610; 1640,- 1650
.1.6^0 GO TU 1570
1.610 r"l=0
;L620 F2--0
GO TO 1420
:i.640 GO TO 1390
j. 6 ?; o GO TO 12*? O
1660 PRjNT "U E.-.L ANÁLISIS ES PARA í
1670 P R 3 N T u J Fí;' E C i J £ N CIA INFERID R J Hz,
:í 680 PRjNT 3 J FRECUENCIA SUPERIOR! Hs t
PRINT 'J ES CORRECTO (SI O NO)
1700 INPUT X$
.!. 7 10 IF Xfc-'S]C " OR X$~ D S' THEN 1790
Xía"ANALISIS PARA UNA FRECUENCIA 1
:¡.7^0 YT^'CAMBIAR O CORREGIR FL RANGO DE FRECUENCIAS'
1 740 GO TO 1530
1750 PRINT U L EL ANÁLISIS ES EN (D,C)a
1760 PRINT "J ES CORRECTO <SI O NO) B?
1770 INPUT X$
1780 IF wQT(Xí*="SI* OR Xí-'S") THEN 1220
17--? O REM
:[Link] rp F1-F2 THEN 2050
1810 REM INGRESO DEL NUMERO DE FRECUENCIAS Y TIPO DE INCREMENTO
1820 REM LINEAL O LOGARÍTMICO
1830 PRINT J L EL ANÁLISIS SE HARÁ PARA UN RANGO COMPRENDIDO ENTRE
18^0 PRINT "JJ FRECUENCIA INFERIOR '-Fl.*' Hs."
1850 PRINT 'J FRECUENCIA SUPERIOR 'jFS.'' Hs * "
PRINT "J INGRESE NUMERO DE FRECUENCIAS "?
1870 INPUT N2
1880
IF F2 10 THEN 1960
.1900 PRINT DESEA :
:L9io PRINT JJ l t - INCREMENTO LINEAL"
1920 PRINT i 2»- INCREMENTO LOGARÍTMICO'
1930 PRINT °J ESCOJA OPCIONGGG 'í
1940 INPUT 2
1950 IF NQTCZ- OR z=2> THEN 1930
1960 R E A D * 11 2 [Link]'JS
1970 WRITE 4-1? O r 1 5 ? J 5 ? Z
1980 GO TO Z OF 1990.2010
.1.990 D-(F2-Fi.V(N2-l¿
2000 GÜ TO 2020
2010 D=CF2/F1)"Í1/ÍN2-1))
2020 REM CARGA AL PROGRAMA fl£ OPCIONES DE SOLUCIÓN
2030 01 = 3
2040 GO TO SOO
2050 PRINT 'J 1 - CALCULO DE VOLTAJES DE NODO"
2060 PRINT "J 2 ~ CALCULO HE VOLTAJES HE NODO Y CORRIENTES DE R A M A
2070 PRINT B J CLASE DESEADA! G " í
 ANEXO 1 ffPR/FRKQ PAO* 14

2080 INPUT H8
2090 IF H8=l ÜR H8=2 THEN 2110
2100 GO TU 2070
21:!.O REM CARGA AL PROGAMA DE CALCULO
21.10 01-4
.?J.:--0 GO TO 800
 ANEXO 1 2PR/ÜPTIOC PAG* 13
1000 REM ********* LINEA 1000 (3PR/OPTIÜC *********
1010 02=3
1020 GO TO 800
1030 REM OPCIONES DE SOLUCIÓN
1040 IF 03-1 THEN .1090
1050 PRINT "L DEBE PRIMERO INGRESAR O LEER DATOS DEL CIRCUITO"
1060 PRINT "J PARA CONTINUAR APLASTE (RETURN) '?
1070 INPIJT' X$
:!.OBO GO TO 100
1090 1T F1O-F2 THEN 1160
1100 P'RINT "LJ PARA ANÁLISIS DEL PROBLEMA EN UNA SOLA FRECUENCIA"
1110 PRINT "J EL CALCULO SE LIMITA A i 1 — LOS VOLTAJES DE NODOS"
PRINT "J 2 — VOLTAJES DE NODOS Y CORRÍ *
"ENTES DE RAMAS"
"JJ PARA CUALQUIER OTRO CAMBIO PRESIONE TECLA #1 "

"L PARA UN RANGO DE FRECUENCIAS LAS SOLUCIONES POSIBLES SONÍ"

'JJ 1 — ALGÚN VOLTAJE DE NODO"
2 — CUALQUIER DIFERENCIA DE VOLTAJES DE NODO'
1210 PñINT B J ([Link] i - VOLT*NüDQ J)"
1220 PRINT "J 3 — CUALQUIER CORRIENTE DE RAMA'
1.230 PRINT "J 4 —- RELACIÓN VOLT, NODO i / VOLT, NODO J'
1240 PRINT "J 5 — a" VOLT+NQDO i / C O R R * R A M A J "
1250 PRINT "J ó — °" [Link] i / CORR * RAMA J"
1260 PRINT 9 J 7 — C O R R * R A M A i / VOLT* NODO J"
1270 PRINT VJ • OPCIÓN DESEADA: G " > '
1280 1NPUT O
1290 GO. TO O QF 1310 ? 1400 1 1350 • 1400? 1400 - 1400 ? 1400
1300 GO TO 1270
1310 PRINT " INGRESE EL NUMERO DEL NODO DESEADO (MAX - ? N ? m > : '?
132.0 INPUT 15
1330 IF !5=>1 AND I3<=N THEN 1670
1340 GO TO 1310
1350 PRINT " INGRESE EL NUMERO DE LA CORRIENTE DE RAMA DESEADA (MAX "?
1360 PRINT Mí " ) í ' 5
1370 INPUT 15
1380 IF I3=>1 AND I5<=M THEN 1670
1390 GO TO 1350
1400 PRINT " INGRESE LOS NÚMEROS DE LOS NODOS Y/O CORRIENTES DESEADOS?."
1410 GO TO O OF 1310? 1420.- 1330? 1440? 1490? 1550 • 1620
1420 PRINT a J INGRESE EL VALOR DE i (MAX " ? N í " > í 'í
1430 INPUT 15
1440 IF I5<1 OR I5>N THEN 1420
1450 PRINT "J INGRESE EL VALOR DE J (MAX 'ÍN?'); H ?
1460 INPUT J5
1470 IF J5<1 OR J5>N THEN 1450
1480 GO TO 1670
1490 PRINT "J INGRESE EL VALOR DE i (MAX " ? N ? ' > I ' í
1500 IF .I5<1 OR I5>N THEN 1490
1510 PRINT "J INGRESE EL VALOR DE J (MAX " ? M í " ) t D í
1520 INPUT J5
1530 IF J5<1 OR J5>M THEN 1510
 ANEXO 1 SPR/OPTIOG PAG, 16
1540 GO TO 1070
1550 PRINT 'J INGRESE EL VALOR DE i (MAX 'JM?'):
15AO PRINT M? " ) '•?
1.570 IF I5<1 OR I5>M THEN 1550
1580 PRINT 'J INGRESE EL UALOR DE J (MAX 'ÍMÍ'K
1590 INPUT J5
1600 IF J5<1 OR J5>M THEN 1580
1610 GO Tu 1670
1020 PRINT "J INGRESE EL VALOR HE i (MAX " ? M í ' > í
1630 IF I5<1 OR I5>M THEN 1620
1640 P'RINT "J INGRESE EL 'UALOR DE J (MAX ' ? N ? " > ;
1650 INF'UT J5
1660 IF J5<1 OR J5>N THEN 1640
16^0 REH
1680 REAÜ *l?2;Z:ZíZ:2
1690 WRITE *l ? 2ÍOrI5/J5?Z
1700 01-4
1710 GQ TO 800
 ANEXO 1 Í2PR/CALCULO PAG, 17
REM ### LINEA 1000 GPR/CALCULO (TECLA 5)
O 2 «4
G O TO SO O
REM **# CALCULO DE LA MATRIZ INCIDENCIA ***
PAGE
CALI. "FILE"? D?.-D£íX$
IF POS (x^^o UPEN- ? 130 >~o THEN ioeo
UPEN D$?l? "F 1 ?X$
READ * l ? 3 Í M j N y N l
DELETE A$?U
DIH A$Omi) jVÍ7)
X$=" "
FOR J=l TO M
X$ = Xí&S'lu
NEXT J
A$^ H "
FÜR 1=1 TO N
A$-A$SX^
NEXT I
TO NI +4
READ *lfi;t;
IF V (1X3 Oft y < l » ó THEN 12¿0
IF V(3)<=0 THEN 1240
A$ = REP ( ' 2 ' * < V < 3 > - 1 ) *:M-K' C 2 )? 1 )
IF V<4)<=0 THEH 1200
A$ = REP í D O B ? < M ( 4 ) - 1 ) #M + V ( 2 ) - 1 )
NEXT I
REM### ALMACENAMIENTO DE MATRIZ INCIDENCIA ***
REM*** PARA M(MAXIMO) DE 30 RAMAS Y N í M A X I M O ) DE 30 NODOS
Oí= "@DAT/ AUXILIAR '
KILL 0$
CRÉATE 0$?9.-H*M*10
OREN Q$?2j. a F " ? X*
FOR j=i TO 9
UJRITE *2fj:0
NEXT J
[Link] *2? 1ÍA*
REM *** ANÁLISIS DE CASO DE CALCULO
DELETE DO ? A$
DIM D0(4>
DO -O
FOR 1=5 TO NI +4
READ *l?i;V
GO TQ y(l)-7 OF 1450íl470rl490íl510
GO TO 1520
nO(4)=l
GO TO 1520
nO(3)=l
GO TO 1520
DOÍ2)-1
GO TO 1520

NEXT I
D 1 = DO < 1 ) *8 + DO C 2 (4
 ANEXO 1 e P R / C A L C U L O P A G * 18

540 BELETE DO
L 550
!. 560 F=0
570 POR 1 = 5 TO NI "1-4
i 580 READ tfl.-IIV
.590 IF V C 1 )*4 OR V < D = 5 THEN 1630
i 600 1F M (1 )<>7 THEN 1640
A 0 FO---1
1620 GQ ro 1640
.630 F=l
-640 NEXT I
A 1 50
o» I:F Fo*.0 AND F=0 THEN 1630
.660 01 =ó
!670 00 TO 1700
i 680 OI- 5
1 690 GO T0 1700
.700 REM
.710 F3"-''::'l
,7;ÍO DEL.E Tu A$
730 GO T0 SOO
tD
<L

o
><
lü
o
OJ

CD
<r
Cu

LU
G)

O
X
UJ

IK Iff
tH
10
o*. ».

^r Kt
II vH
rH O **
V—i O O N
O O-
CN i-t tn
TH *. 11
CD -O vH o
••_•• O CO M O CO
OJ «í- ~¡ oj <r
•••s U íil OJ m oj
T—1 R^
---- rH •M •- rH
O <!
"3 •<! LU »• O Q
•s»* -w- ». O
^> O O rH X- O fx <-
1-1 LL 2> ^-s
Lü C O O- w ^>
rx o «
O O ^ O rH '"-. C ~J O OJ
*t Lü •^
<-s !— OJ rH M i—i "• " oo 01 '• U
CO TH M *, O rH
~J 2: CU & *• '• C> Lü
CN #• ^^- O CN Ci
*. Lü rH *-- M O !> 11
"i O rH 1—1 -^~*' ÍG <í O
'rH O -J Í3l —' W
~f ••-' TH CN C'J C) OJ Ú1 CN OJ LU'
-0 b_ *~, Oí O
OJ <3- Lü > CD o-í LV: + -^r !*H *•-'
££ U. O rH -í- Lü > ~J LL LE:
~f- OJ -CC -r X 11 -H
.*"s rH j_ ^ ^N Lu O ü£ CD r-f X II Lu O L¿
21 í- •-. 1C n O Lü CO 1C
^_
su V
-U. X- J --. s: M U C¿ x;
O
^ Z> *"•> >^- Oi Ü
hJ -—• > O M ^1 tH LÜ
o O rH f rH CD tH •"• O
C- X 0 LU XI O ** TH *- O rH tO rH CD _Jí ^ o *•- co ro
o O O M [_ H CD 1 ¡ C3 c, _ ^: H
i— ¡— ••_' ÍX }- LTJ !— rH >. !— n A. => í~ Lü I ! CD ZD •x. i- fH A. -^
i fEj -rH -rH CD —i *>
iH L¿ N í-9 ü"; CD 31 V '- U vH TH CD O CCi 11
^J TH ,—1 .-—• !-! !—! n 3Ü t~i
t]"j rH ^-v í-s TH --•• V-1 V-t H ^J 3: ü"j
rH f] -r-í rH !! OJ Oí OJ I-H *-'
'X ÜJ — ii ft rH -~v 1-1 11 ü ~l ~J M LÜ Lü ^< II
Ü ~ M !! 11 ~ -J rH LÜ Ü X LÜ TH 11 * -rH <- 1-H ü Í-H 5—1
M •v^ *í" H -J =, í— O O h- I— _J I—1
H- CQ TH O Lü i— CD hH •^ ^ H r. O O i- U í- ffi
rH í-H M "J J
V- 1- h- 2: Lü í-t 1> •"-' h- rH h- H Lü H h~ 2: Lü
h- ». ¡_ Í— Lü ID C? 1— Lü O W > Z> f- i-i
>< •w* X üi ^-" M k-I li 1 x: o L¿ <E > x c¿ ££ •- X X _l w Q 211 ,J líH 0 L¿ X
X C¿ tV t-l X x _i en E: _j :s: 11 L¿ 'X
Lü O O üi 21 Lü Lü M ¡| 0 Lü Lu --- Lü 0 O rH LÜ LU LU LJ Q ü¿ 2: Lü Lü M 11 0
o •"• Lü 0 0 -^ LULU Lü O Lu O Lü Lü VH rH O LÜ LU rH LÜ O rH
uj ÍT) Lu L¿ M PJ 2: LL Lu CD 21 21 W! CD CD ü_ U.1 Cu' (-.! V-1 r> iL.
í- 1 21 Lu Lu =i 2! 21 W CD M CD ££ í-i <-) CO b_ íj£ rH CD 2: Lu CD 2: CD CD ü_ Lü L¿ Q
O O O O O O O O O O O O O c- o 0 O0 o ooooo oooo o 0
oo 0 C' O O O o ooo o opo O OO o oo oo O rH OJ EO -í LO ^O N 03 CN O rH c>! ro
Oí Í"'J •e if; ••o ¡x o?. G-- O TH O-l í'"J "3" i* J *>"-> ¡x 00 O- O t-l M '<"1 ^ LD s3 ¡X 03 0-- O rH O-í M M- LO -í! f^ OJ CN Ci O O O O O O rH rH
•j-i -.j-f -n -Ti iX ¡X FX ÍX ÍX ÍX ÍX EX ¡x rx o. =>
03 CO CO O:' CO CO CO 03 CTj CO CN CN CN CN 0*- 0'- O1- 0^ o oo
_i _i _:r -J T-l
 ANEXO 1 ePR/YEi PAG, 21
FOR J=5 TO Nl + 4
READ #1?J:V
IF V U ) < > 1 1 THEN 2200

2770
GÜ TO 01-7 OF 2260 r 2350
GO TO Dl-11 OF 2350*2350
RR1NT "GGGERROR Dl= H ?D1
EHD
DELETE Gl
IHM G l í M - M )
CU -O
2290 FOR 1=1 TO H
2300 FOR J=l TO M
2310 G 1 < I ? J ) =E 1 < I ) * U í 1 ? J )
2320 NEXT J
2330 NEXT I
2340 GO TU 2940
2350 DELETE G2
2300 DIM G2(H?M)
2370 G2-0
2330 FOR 1=1 TO M
2390 FOR J=l TU M
2400 FOR K=l TO M
2410 G2 ( I y J ) =G2 ( I f J ) -i-G 1 ( I r K ) *U C K r J )
2420 NEXT K
2430 NEXT J
2440 NEXT I
2450 G1=G2
2400 DELETE G2?U
2470 GO TO 2940
2480 REH Dl=10 .
2490 HELETE U
2500 DIM UCH.-M)
2510 U-0
2520 FOR 1=1 TO M
2530 U<IíI)=l
2540 NEXT I
2550 DELETE R
2500 DIH R(M/M)
2570 R=0
2580 FÜR ]>5 TO N1Í4
2590 READ * l j I t V
2000 IF V(l)<>9 THEN 2030

Có)
 ANEXO 1 ePR/YEl PAG, 2
2700 NEXT I
2710 GÜSUB 2770
2720 GO TO ül-9 OF 226092350
2730 GO TO Bl-13 OF 231 0 50
2740 PRINT "GGGERROR fl:
2750 ENB
2700 REM
2770 REM INVERSIÓN
2780 FOR K=l TO M
2800 FOR 1 = 1 TO h"
2810 IF K-I THEN 2S70
2830 FOR J=l TO M
2840 IF J •K THEN 2800
2850 U(I? K rJ
2800 NEXT
2870 NEXT I
2880 FOR J=l TO H
2890 IF J =K THEN 2910
K)
2910 NEXT J
2920 NEXT K
2930 RETURN
2940 X$""§PR/YI
2950 01=28
2960 GO TO 920
 ANEXO 1 EPR/YEIINCIÜ PAG, 2
1000 REM ************ LINEA 1000 ©PR/YEIINCID
1010 02-2S
1020 GO TO SOO
1030 REM CALCULO DEL A*YB ( M INCIDENCIA #YB CALCULADA)
1040 REM IMPRESIÓN Ye> Yb
1041 IF H7=0 THEN 1060
1050 PRINT @F'9: " JVECTOR DE ADMITANCIA
A DE ELEMENTOS' * El
1060 IF Dl=0 THEN 10SO
.1001 IF H7=0 THEN 1080
1070 PRINT (3P9Í ' JMATRIZ. YBa Gl
1080 DELETE A $ ? Y l j Y
1090 DIM A*<M*N) ? Y 1 C N ? M )
1100' Yl~ 0
1120 READ *2?líft$
1150 IF D1OO THEN 1270
11 00 FOR 1=1 TQ N
1170 FOR J=l TO M
USO X* = SEG ( A$ ? C 1-1) *M-f J ? 1 >
1190 T = yAL (X$)-f 1
1200 GO TÓ T OF 121051240^1 O
1210 YKI? J)--El(J)
1220 GO TQ 1240
1230 YKIf J)-EKJ)
1240 NEXT J
1250 NEXT I
1200 GO TO 1410
1270 REM PRODUCTO (A*YB)=Yl
1280 FOR 1=1 TO N
1290 FOR J=l TQ M
1300 Y 1 C I ? J)=0
1310 FOR K=l TO M
1320 X$ = SEG< A$? C 1-1 >*M-fK? 1 )
1330 T"'v'AL CX*)-f 1
1340 GO TO T OF 1350? 1380?!-%
1350 Y 1 C I ? J)=YKI? J)-G1(K? J)
1360 GO TO 1380
1370 Yl ( I .- J ) =Y1 ( I ? J ) iGl (
)K * J
1380 NEXT K
1390 NEXT J
NEXT I
REM
REM ALMACENAMIENTO Yb
IF D1OO THEN 1490
j ¡ T**I ~r ~t* F— .s. *"^ í~^ *• r** ^
W K o. 1 ti í x; * V f 11
GO TO 1501
WRITE #2?9ÍG1
DELETE G1?E1
DIM Y C N í N )
Y=0'
REM Y1*A' DONDE A TRANSPUESTA DE A
FOR I=i TO N
FOR J=l TO N
YCI?J)=0
FOR K=l TO M
 ANEXO 1 0PR/YE.1.INC1D PAG* 24
1560 X*=SEG(A*> (J-l)#M-fK, 1)
1570 T=VAL(X*)-H
:i.S80 GO TO T OF 1590, 1620 / 1610
1590 Y (I- J)=Y<I> J)~Y1(I?K)
1000 GO TO 1620
1610 Y(IrJ)=Y<IjJ)ÍY1CI?K)
1620 NEXT K
:LÓ30 NEXT J
:l.640 NEXT 1
:¡.Ó50 IF H7-0 TREN 1080
:tó5i PRINT @P9:" JMATRIZ YI=A*EICGD ' T Y I
1060 PRINT @P9t " JMATRX7. DE ADMITANCIA DE NODOS A*Y1*A' =Y' , Y
1670 PRINT (?P9Í" JMATRIZ DE INCIDENCIA ALMACENADA EN DISCO " ? A $
1680 WRITE *2í4íY
1730 DELETE Y
.1740 GO TO 1790
1790 X*="t?F1R/YElSOLU"
1800 01=30
1810 GO TO 920
 ANEXO 1 OPR/YE1SOLU
1000 REM ************ LINEA 1000 SPR/YEISOLU
1010 02 = 30 *************
1.020 GO TU 800
L030 REM CALCULO DE TÉRMINOS INDEPENDIENTES
1040 IF F1OF2 THEN 1070
1050 IF Fl=0 THEN 1130
.000 GO TO 1480
.070 PRINT 'L PARA UN CIRCUITO SOLO RESISTIVO EL CALCULO SE LIMITA
loso PRINT 'U A UNA SOLA FRECUENCIA "
1090 GO TO 1480
1.100 REM ****** ANÁLISIS PARA CD,C) ******
!. 110 REM CALCULO DE OOG1*E
1120
1.130 DELETE E,C?.C1
1.140 DIM E(M) Í C(M) ? C1(N)
.150 E=0
U 60
L170
.180 TO Nl+4
1190 V
1200 THEN 1240
L210 THEN 1250
220
1230
.240 E C V Í 2 )
L250 NEXT I
.260 FOR 1=1 TO N
.270 FOR J=l TO M
L2SO Cl i I ) =C1 ( I ) -Yl C I? J ) *E C J >
90 NEXT J
L300 NEXT I
.310 IF H7=0 THEN 1340
.320 PRINT OP9: "JVECTOR DE FUENTES DE VOLTAJES INDEPENDIENTES
1330 REM NECESITAMOS PARA ALMACENAMIENTO
1340 FOR 1=1 TO N
350 FOR J=l TO M
300 XÍ» = SEG í A* ? C 1-1 ) *M+ J ? 1 )
370
.380 00 TO T OF 1390-1420?1410
>90 CICI)=C1ÍI)-CCJ)
L400 GO TO 1420
.4.1.0 C1<I)=CICI)+C<J)
.420 NEXT J
.430 NEXT I
.440 IF H7=0 THEN 1460
.450 PRINT £P9; a JMATRIZ DE TÉRMINOS INDEPENDIENTES (Is) CÍ
i 460 DELETE C,Y1
470 GO TO 1960
.480 REM ****** ANÁLISIS PARA (A.O ******
490 REM E CONTIENE PARTE REAL' DE Va
00 REM EO CONTIENE PARTE IMAGINARIA DE V<á
. 510 REM C REAL DE la
520 REM co IMAGINARIA DE
1 ÍO DELETE C ? COyE?EOíC:UC2
 ANEXO 1 0PR/YE1SOLU PAG. 26
1540 E C M ) r E O C M ) .CCM) ,CO(M) >C1(N) r C 2 < N >
1550
.1560
1570
.1580
1590 Cl = 0
;L600
1 ó 10 REM *** FORMACIÓN DE la y Va (COMPLEJOS) ***
1020 SET DEGREES
1030 FOR 1=5 TO Nl+4
1040 READ #l?i:y
1650 IF V U > ~ 1 THEN 1700
IF V(l)<>2 THEN 1720
16 ó O'
1670 C (V(2))=yC 5)*COS(V C 7))*10~V(6)
:LÓSO CO(V C 2))=V(5)#SIN(V < 7))*10"y(ó)
1690 GO TO 1720
1700 E < y ( 2) )< 5 > *-cos í u < 7 > > <ó>
:L 710 EOCVC2)
1720 NEXT I
1730 FOR 1=1 TO N
1740 FOR J-l TO M
:i.750 Cl ( I > -Cl (I > --Y1 < IJ) #E ( J )
1760 C2 < I ) =C2 < I) - Yl (JI) UCEO ( J )
.1770 NEXT' j
1780 NEXT I
1790 IF H7-0 THEN 1820
1800 PRINT "J Vsí" í E í E O ? "IG
.1810 REM £#«## FORMACIÓN DEL VECTOR Is
1820 FOR 1=1 TO N
1330 FOR J=l TO M
1840 X* = SEG í A$ , < I--1 ) *H + J > 1
1850
1860 GO TO T OF 1870rl920f 1900
1870 CKI)=C1(I)-C(J)
1880 C2<I)"C2(I)~CO(J)
1890 GO TO 1920 .
1900 C1(I)"C1CI)+C(J)
1910 C2(I)==C2CI)+COCJ)
1920 NEXT J
1930 NEXT I
1940 IF H7-0 THEN 1970
1950 PRINT @P9Í' TÉRMINOS INDEPENDIENTES SON ' C1-C2
1960 REM LECTURA DE MATRICES DE COEFICIENTES
1970 DELETE Y
1980 HIM YCN.'N)
1990 REAÜ *2f4íY
2000 IF F1OF2 THEN 2020
2010 IF Fl = 0 THEN. 2100
2020 REM PARA ANÁLISIS C A > C ) NECESITAMOS GAUSS GORDAN COMPLEJA
2030 DELETE YO
2040 DELETE YO
2050 DIM YO(NfN)
2060 YO-0
2070
 ANEXO 1 (?Pr</YElSQLU PAG, 27
20SO 01=34
2090 GO TO 920
2100 REM SOLUCIÓN DEL SITEMA DE ECUACIONES (ELIMINACIÓN DE GAUSS)
2110 REM Y#Vn=Is
2120 IF N=l THEN 2580
2130 REM CON PILOTAJE SIMPLE
2140 DELETE R
2150 R==0
2100 POR 1=1 TO N-l
2170 R=I
21SO T=ABS<Y(I?I))
2190 POR J = I4-1 TO N
2200 IF T=>ABS(Y<JjI» THEN 2230
2210 T=ABS(Y<J?I))
2220 R=J
2230 NEXT J
2240 IF TO1,OE-12 THEN 2550
2250 IF R=I THEN 2360
2200 REM INTERCAMBIAR FILAS R y I
2270 POR K=I TO N
*"* o o f~) ™r — v í T k** ^t

2300 Y(R?K)~T
2310 NEXT K
2320 T=C1(I)
2330 CKI)=CKR)
2340 C1(R)= T
2350 REM CEROS BAJO PIVOTE
2300 POR K=I+1 TO N

23SO POR J=I+1 TO N

2390 Y í K í J ) =Y ( K, J ) +T)KY CI, J )
2400 NEXT J
2410 C1(K)=C1<K>+T#C1<I)
2420 NEXT K
2430 NEXT I
2440 REM EVALUACIÓN DE INCÓGNITAS Vn
2450 IF ABS(YÍN.'N) )< = !-> OE-12 THEN 2550
2460 Cl<N)=C1CN)/Y<N?N)
2470 POR I=N-1 TO 1 STEP -1
24SO S-0
2490 POR K=I-K TO N
2500 S=S+Y(IíK)*Cl(K)
2510 NEXT K . .

2530 NEXT I
2540 GO TO 2590
2550 PRINT "LMATRIZ HE ADMITANCIAS DE NODOS DEL CIRCUITO ES SINGULAR*
2560 PRINT "JPROBLEMA CON DATOS MAL ESTRUCTURADOS, REVISARGG'
2570 END
2580 CKl)=Cl(l)/Y(ljl)
2590 IF H7=0 THEN 2030
2600 PRINT g'P9: ' JVECTQR COLUMNA DE NODOS (Vn) ",C1
2610 PRINT 'JJ PAPA CONTINUAR APLASTE (RETURN) " í
 ANEXO i EF'R/YEISOLU PAG. 28

2020 INPUT X$
2030 X*="@F'
2640 01=42
2650 GO TO 920
 ANEXO 1 0PR/YE2 PAG, 29
1000 REH ************ LINEA 1000 PR/YE2 ****
1010 Q2--6
GO TO SOO
REM OCUPA 095? b DE 30530 b
DELETE E1 ? E2
DIM E1ÍH) 7E2(M>
E1 = 0
E2-0
11 = 1
IF F1OF2 THEN 1110
IF Fl=0 THEM 1130
W=2*PI*F3
Gü TU 1140
W=1*OE-10
REM
FOR 1=5 TO NI +4
REAÜ #l*i:V
GO TO y < l > - 2 OF 1190.1210,123071250
GO TO 1200
U E K V C 2 ) >=E1(V<2) )+l/ÍV<5)*10~V(ó>)
GO TO 1260
E2<V(2) > = E 2 < V < 2 ) ) ~l/< U*( V (5) *10"Y (ó) ) )
GO TO 1200
E2CVÍ2) )=E2(V<2) >+W* C '-,' (5) *10"V (6) )
GO TO 1260
Eli V C 2) )=E1 W2»+V(5>*iO~V<6)
NEXT I
IF H7=0 THEN 1290
PRINT eP9t'JVECTORES DE AÜMITANCIA DE ELEMENTOS" , El , E2
REH ALMACENAMIENTO DE Ye (ADMITANCIA DE ELEMENTOS)
WRITE *2 ? ÓÍE1
WR1TE *2?7ÍE2
GO TO DI OF 1370? 1030? 1370? 1670? 1670
GO TO Dl~5 OF 1670rló70íS050y 1370? 1630
GO TO Dl-10 ÜF 1370flo70ilo70ile.70i.1070
GQ TO 2070
REM Dl=l FCCC <B?
DELETE U
DIM U ( M f M )

FOR M

H£XT I
FOR J=5 JO Nl+4
READ * l T J t y
IF 'v'(l)<>8 THEN 1470
U < V < 4 ) ? V < 3 » = U < V < 4 > ? M < 3 » + V < 5 ) *10"V C ó )
NEXT J
DELETE G1 ? G2
DIM GKM.'M) ? G2(MíM)
Gl = 0
62 = 0
FOR 1=1 TO M
FOR J=l TO M
 ANEXO 1 GPR/YE2 PAG*' 30

540 G l C I r J)=U(I* J)*EKJ)

550 G2CI*J)=UCIjJ)*E2<J>
SóO NEXT J
570 NEXT I
80 DELETE U
1590 GO TO DI OF 26707l61Q.'1630
L600 GO TO Dl-8 OF 2050,1610?1630
.010 PRINT "GGG ERROR Dl= 'JD1
620 END
.630 REM CARGAR PROGRAMA @PR/YE2CONT
.640
-¿50 01-45
660 GO TO 920
.670 REM EXISTE UNA FUENTE FUCC (0) Dl=
_680 DELETE GlíG2
.690 DIM G 1 C M . M ) ? G 2 ( H f M )
.700 Gl-O
L710 G2~0
-720 FOR 1=5 TO Nl+4
730 READ * 1 ? I Í V
_740 IF UCDO10 THEN 1760
,750 GKVC4) 7^(3) )=VC5>#10~V(Ó>
.760
1770 FOR 1=1 TO M
r^H f ~r T \• r~ -! / T \r / T T \ \ \' i -4* i j ( J * i
.730 *— X \. J J i ^ l ' M L _ J _ \ _ ¡ , . ' T L ? - í - k , J - 7 -L J
,790 G2<I jI)=E2CI)'
LSOO NEXT I
810 GO TO EI 1-3 OF 2670fl850j 1630*1850
820 GO TO Dl-11 OF 2050?1850?1030*1850
.330 PRINT 'GGGERRORD1== 'ÍD1
.840 END
.850 REM CALCULO G14-B*YE IH-5,7.
_860 DELETE B
.870 DIM B C M . M )
LSSO
890 FOR 1-5 TO Nl+4
.900 READ * l í I t V
.910 IF y<l)O8 THEN 1930
.920
.930 NEXT I
.940 FOR 1=1 TO M
950 FOR J=l TO M
.960 GlCI,J)=G1(I *J)+B C 1 7 J)*£1(J)
.970 G2 < I * J ) =G2 < 1 9 J ) 4-B ( I, J ) :ííE2 < J )
980 NEXT J
.990 NEXT I
'000 DELETE B
J010 00 TO 01-4 OF 2670,203071630
020 GO TO Dl-12 OF 2050*2030*1030
2030 PRINT "GGG ERROR Dl= "?01
1040 END
2050 REM DI-=8
ÍOóO DELETE U
DIM U ( M * M >
 ANEXO i (2PR/YE2 P A G , 31

2080 u=o
2090 FOR 1=1 TO M
2100 U < I? I) - 1
2110 NEXT I
2120 FOR J=5 TO Nl+4
2130 REAU *i?j:v
2140 IF U(1)<>11 THEH 2160
2150 U ( V ( 4 ) - V C 3.) ) =U ( M ( 4 ) ? V C 3 ) ) + V
2160 NEXT J
2 170 REM SIMPLE INVERSIÓN
2 ISO FOR K=l TO M
2190 U(K?K)=i/U(K7K>
2200 FOR 1=1 TO h
2210 IF K=I THEN 2270
'JIIQ U(IíK)=-U(I;K)*U(KíK)
2230 FOR J=l TO M
2240 IF J = K THEN 2200
2250 U ( I ? J ) =U ( I r J ) + U ( I ? K ) #U ( K ? J )
2260 NEXT J
2270 NEXT I
2280 FQR J=l TO M
2290 IF J = K THEH 2310
2300 LKK? J)=U(K-J)*UCK-K)
2310 NEXT J
2320 NEXT K
2330 REM ALMACENAMIENTO DE U
2340 WRITE #2? Sí U
2350 GO TO IH-7 QF 2390? 2510
2360 GO TO Dl-11 OF 2510^2510
2370 PRIHT 'GGGERROR 01= 'íüi
2380 END
2390 DELETE G1 ? G2
2400 HIM G l ( M í M ) í G 2 C M 7 M )
2410 Gl™0
2420 G2 = 0
2430 FOR 1=1 TO M
2440 FOR J=l TO M
2-450 SKI'? J)=E1 (l)^UCI.. J)
2400 G2 < 1 7 J ) = E2 (I ) £U C I } J )
2470 NEXT J
2480 NEXT I
2490 HELETE U
2500 GO TO 2670
2310 HELETE G3jG4
2520 DELETE El ? Eli
2530 DIM G 3 < M ? M ) . G 4 C M ? M )
2540 G3=0
2550 G4 = 0
2560 FOR 1=1 TO M
2570 FOR J=l TO M
2580 FOR K=l TO M
2590 G3 ( I ! J ) =G3 í I : J ) -i-Gl C 1 7 K* ) $U
J )( K
2600 G4 ( 11 J ) =04 C 1 7 J ) +G2 C 1 7 K ) #UJ() K
2010 NEXT K
 ANEXO .1 8PR/YE2 P A G - 32

2620 NEXT J
2630 NEXT I
2640 (31 = 63
2050 G2 = G4
'2660 [Link] G3rG4 ? U
2070 REM
2080 X$="Í?PR/YE3INCID"
2690 01=32
2700 GO TO 920
 ANEXO 1 (3PR/YE2CQNT P A G » 33

1000 REM ************ LINEA 1000 0PR/YE2CONT

1010 02=45
1020 GO TO 800
1030 REM
1040 DELETE R r U l 7 U 2
1050 DIM R C M > M ) ? U K M í M ) jU2(M?M)
'1060 R=0
1070 Ul = 0
1080 U 2=0
1090 FOR I 7= 1 TO M
:L100 Ul < I i I )= 1
1110 NEXT J
1120 GO TO Dl-9 OF 1200*1200
1130 GO TO D 1 - 1 3OF 1200,1200
1140 FOR 1 = 5 TO NI + 4
1 150 REAÜ *1 1 1 : v
1100 IF y< 1 )<>9 THEN 1180
1170 R ( V ()4y U ( 3 )=V ) < 5)#10"'v!(6)
USO NEXT I
1190 GO TO 1 2SO
1200 FOR I = 5 TO NI + 4
:L210 READ dt 1 ? i í y
1220 IF V < 1 )<>9 THEN 1250
1230 R C V C 4 ),V ( 3 )=V<5)*10"y<6)
)
1240 GO TO 1270
1250 IF V < 1 )Olí THEN 1270
:L2¿0 U 1 C V (4)j V < 3 ))=UKV(4) ? V (3) )+V<5)*10"V(ó-)
1270 NEXT I
IL2SO FOR I = 1 TO M
1290 FOR J =z 1 TO M
1300 U l C I í J )= U 1 C I j J)+R-<Ií J)*E1CJ)
1 3 1 0U2 < 1 1 J )= R<I .• J)*E2<J)
1320 NEXT J
1330 NEXT I
1340 REM ALMACENAMIENTO DE Ul
1350 WRITE *2?8Í Ul
:1360 DELETE R
1370 REM S UB , INVERSIÓN DE MATRIZ COMPLEJA — HET» DE SHIF'LEY C S H I P L E Y C )
1330 GOSUB 1 750
1390 IF DI z= o OR rn=10 THEN 1410
1-400 GO TO 1 530
1410 DELETE G l í G 7'
1 420 DIM G 1 (M ? H ) j 62 ( M f M)
1430 61=0
1440 62=0
1450 FOR I =3 1 TO M
1460 FOR J=l TO M
1470 GKIí J )= E1( I) *U1 < Ir J ) --E2 < I ) *U2 ( 1 7 J )
1480 G2 C I j J )= E1( I)*U2ÍI? J)+E2(I)*U1(I7 J)
1490 NEXT J
1500 NEXT I
1510 REM DEL Ul? U2;E1 ? E2
1520 GO TO ii 470
1530 REM DEL El r E 2
 ANEXO 1 (2PR/YE2CONT P A G , 34

1540 DELETE G3><G4

1550 D I M G 3 ( M , M ) , G 4 Í M ? M )
1560
1570
1580 REM ALMACENAMIENTO DE Gl
1590 WRITE #2.'9ÍG1
1600 FOR 1=1 TO M
1610 FOR J=l TO M
1620 FOR K=l TO M
1630 G3(I?J)=G3(I J)+GKI?K>*UlCKí ) - G 2 C I f K ) * U 2 C K ? J)
1640 G4(l"y J)=G4<I J)+Gl(I?K)*U2(Ky
1650 NEXT K
1660 NEXT J
1670 NEXT I
1680
1690 02=04
1700 DELETE G3íG4?UÍ?U2
1710 REM
1720 X*=="@Pft/YE3lNCID"
1730 01 = 32
1740 GO TO 920
1750 REM SUB+ INVERSIÓN DE MATRIZ COMPLEJA — MET > DE SHIPLEY (SHIPLEYC)
1760 REM
1770 REM M — ORDEN DE LA MATRIZ
1780 REM Ul — MATRIZ DE LAS PARTES REALES
1790 REM U2 — MATRIZ DE LAS PARTES .IMAG,
1800 REM INVERSA QUEDA EN LAS MISMAS LOCALIZACIONES Ul -f JU2
1810 REM
1820 FOR Q0=l TO M
1830
1840 U i(QO-QO > =U1í QO,QO)*G
1850 U2(QO ? QO)=-U2< QO fQO)#Q
1860 FOR Qi=:i TO M
18 7 O IF -01=50 THEN 1970
Q=~U1í QO y QO > #U1<Ql ? QO)+U2C QOfQO)#U2(Ql,QO)
U2(Q1j QO)=-U1<GO? Q O > #U2(Q11QO)-U2(QO,QO)*U1(Q1fQO)
Ul(QlyQO)~Q
FOR Q2=l TO M
IF (32 = 00 THEN 1960
Q™U1(QlYQ2)+U1(Ql,QO)*UIC QO? Q2 >-U2 <Q1-QO)#U2 C QO tQ2)

UKQlíQ2)=Gl
NEXT Q2
NEXT Ql
FOR Q2-1' TO M
IF Q2=QO THEN 2030
Q = U1 ( QOf QO ) XíUl C QO t Q2 ) -U2 ( QO t QO ) *U2 ( QO ? Q2 )
U2 < QO , Q2 ) =U 1 < QO f QO ) *U2 C QO 7 Q2 > -f U2 ( GO j QO ) *UQ2
1 () QO
UlCQOyQ2)=Q
NEXT G2
NEXT QO
RETURN
 ANEXO 1 (2PR/YE3INCIB PAG, 3!

1000 REM ************* LINEA 1000 0PR/YE3INCIB

1010 02-32
IL020 GO TO 300
.1030 REM CORRECION DE MATRIZ DE ADMITANCIA R A M A L SI EXISTE I ÑU* MUTUAS
1040 REM Y ALMACENAMIENTO
1050 REM PRINT £P9t'MATRIZ. YB SIN M O D I F I C A R a r ' G l ' , G l • 'G2', G2
1060 IF FOO1 THEN 1680
1070 IF D1OO THEN 1450
.1080 REM CORRECION DE YB CUANDO ES MATRIZ DIAGONAL COMPLEJA
1090 DELETE G1 ? G2
1100 DIM G K M / M ) r G 2 ( M » M )
1110 Gl = 0
1120
1130 FOR 1=1 TO M
1140 GKIí I)=E1(I)
1150 G2<IjI)=E2(I)
1160 NEXT I
1170 DELETE ElíE2
USO GOSUB 1200
1190 GO TO 1290
1200 DELETE L
1210 DIM L C M )
1220
1230 FOR 1-5 TO Nl+4
1240 READ *l?i:V
1250 IF U(l)<>4 THEN 1270
1260 LCVÍ2))=VC5)*10~Y<Ó>
1270 NEXT I
1280 RETURN
1290 REM CORRECCIÓN
1300 FDR 1=5 TO Nl+4
13 í O READ #lfi;V
1320 IF y<l)<>7 THEN 1430
1330 T2-G2CUÍ3) ? VC3))
1340
1350 T5=LCV(4>>*L<VC3)>-M<5)*V<5)*10~<2*V(ó))
1360 IF ABS(T5)>1(OE-16 THEN 1380
1370 GO TO 1560
1380 G2 ( V < 3 ) .- U < 4 ) ) — C T4*V C 5 ) *10"V ( 6 ) *L ( V ( 4 ) ) ) /T5
1390 G2 C V C 4),V(3))=-(T2*V <5)*10"MC 6 > *L(V < 3)»/T5
1400
1410 G2(y(3)? VC3))=T2*T5
1420 G2( y (4) y V C4) )=T4*"T5
1430 NEXT I
1440 GQ TO 1690
1450 REM CORRECCIÓN DE YB CUANDO ES UNA MATRIZ COMPLEJA
1460 GOSUB 1200
1470 FOR 1=5 TO Nl+4
1480 READ *l?i:V
1490 IF V C D O 7 THEN 1670
1500
1510
1520
1530 T4-G2CVC4)
 ANEXO 1 8PR/YE3INCin PAG» 36

1540 T5=L<V(4) )#L(V(3) )-V(5)£'v'(5) *1Q"<2#V(6) >

1550 IF ABS(T5»1 .OE-ló THEN 1580
1560 PRINT 'J EXISTE ERROR EN LOS DATOS DE INDUCTACIAS (REVISE)
1570 END
:L580 Gia'(3) : V C 4 > )=-( V(ó))*L(V(4))/T5
1590 V ( 6 ) # L < V ( 4 ) ) )/T5
100 O G1(V(4);U<3))=-( V(Ó)#L(V(3)))/T5
1610
1620 T5=L(V(3)>#LCV<4))/T5
i~M f \ / ~Z \ I f "T \T -\
1030 b 1 ', V \ ) .' V (. o ) ) — I 1 # I O
1040
1650 G K V ( 4 ) 7 V ( 4 ) )=T3*T5
1600 G2(V(4)?V(4))=T4*T5
1670 NEXT I
1ÓSO REM
1690 REM CALCULO DEL A*YB
.1.700 DELETE A $ ? Y l i Y 2
:i7io DIM A*(M*N) -YKN'M) /Y2(N ? M)
1720 Yl=0
.1730 Y2=0
1740 READ *2?ltA*
1750 IF F0=l THEN 1910
1760 IF D1OO THEN 1910
1770 FOR 1=1 TO N
1780 FOR J=l TO M
1790
1SOO
1810 GO TO T OF 1820 187071850
1S20 Y1U j J)=~E1(J)
1830 Y2(I? J)=-E2(J)
1840 GO TO 1870
1850 Y l í l r J)=E1(J)
1300 Y2(I? J)=E2(J)
1870 NEXT J
1880 NEXT I
1390 0ELETE E1.E2
.1900 GO TO 2130
1910 REM PRODUCTO A*G2=Y2
1920 DELETE E1 ? E2
1930 FOR 1=1 TO N
1940 FOR J=l TO M
1950 FOR K=l TO H
1960 X*=SEG(A$y (1-1 )
1970
1980 GO TU T OF 19907204072020
1990 Yl í I ? J ) = Y 1 < I f J ) ~ G 1 ( K , J )
2000 Y2CI?J)=Y2(I>J)-G2(K»J)
2010 GO TO 2040
2020 YKIr J)«YKIj J)ÍG1(K/ J)
2030 Y2 < I j J)=Y2(I *J)+G2< K 7 J)
2040 NEXT K
2050 NEXT J
2060 NEXT 1
2070 REM ALMACENAMIENTO Yb
 ANEXO 1 (?PR/YKolNCID PAG. 37

2080 WRITE -42>2;G1

2090 WRITE *2r3;G2
2100 IF H7=0 THEN 2120
2110 PRINT @P9Í"JMATRIZ YB " r G i ? G 2
2120 DELETE G1*G2
2130 REM
2140 REM Y1#A' DONDE A' = TRANSPUESTA DE A
2150 DELETE YíYO
2100 tUM Y C N r N ) í Y O ( N ; N )
2170 Y=0
2180 Y0=0
2190 FOR 1=1 TO N
2200 FOR J=l TO N
2210 FOR K=l TO M
2220 X*=SEG(A$><J-1>#M+K?1)
2230 T=VAL(X*)+1
2240 GO TO T OF 2250,2300j2280
2250 Y<IíJ)=YCI?J)-YlCIíK)
2200 YO(I ? J)=YO<I7J)-Y2 C I > K )
2270 GO TO 2300
2280 Y(l7 J)=YU7 Jm'KI.'K)
2290 YO(I,J)-YO(I> J)+Y2(IjK)
2300 NEXT K
2310 NEXT J
2320 NEXT I
2330 IF H7=0 THEN 2360
2340 PRINT @P9Í"JMATRIZ Y1=A*E1ÍG1) i ? Y l ? Y 2 .
2350 PRINT (3P9Í " JMATRIZ DE INCIDENCIA ALMACENADA EN DISCO
2360 WRITE *2?4iY
2370 WRITE *2-5íYO
2330 IF H7=0 THEN 2400
2390 PRINT @P9:>JMATRIZ DE ADMITANCIA DE NODOS A*Y1*A / -Y"
2400 DELETE Y 7 YO
2410 X$= i ePR/YE3SOLU 1
2420 01=34
2430 GO TO -920
 ANEXO 1 SPR/YE3SGLU PAG* 33

REM ********** LINEA 1000 (3PR/YE3SOLU

02=34
GO TO 800
REM ********* CALCULO DE TÉRMINOS INDEPENDIENTES *********
REM ********* Y VOLTAJES DE NODOS *********
IF F0=0 AND F O THEN 2260
IF F1OO THEN 1080
IF F1=F2 THEN 1420
REM ****** ANÁLISIS PARA Í A * C )
REH E CONTIENE PARTE REAL DE Va
REM EO CONTIENE PARTE .IMAGINARIA DE Va
REM C REAL DE Isí
REM CO IMAGINARIA DE I sí
CÉLETE C ? CO,E ? EO ? C1 ? C2
DIM E(M) ? EO(M) ? C<M)jCOCM)

EO-0

CO-0
REM COMPROBACIÓN DE EXISTENCIA DE FUENTES
IF F0=l THEN 1310
IF D1OO THEN 1310
GOSUB 1910
GOSUB 1780
REH PRODUCTO l£*Vá <E1+E2J>*<£+EOJ)
FOR 1=1 TO M
C CI)=C<I)-El(I)*E <I)+E2 <I)&EO <I)
CO<I>=CGíI>~E1<I>#EO(I>-EíI>#E2<I>
NEXT I
BELETE E1?E2
GO TO 2030
GOSUB 1970
GOSUB 1780
REM PRODUCTO DE Yb*Vá (G1+G2J)*(E+EOJ)
FOR 1=1 TO M
FOR J=l TO M
C<I)~-=Cf IJ-G1CI J ) *E ( J ) + G2 ( If J ) *EO < J )
I 3 J)*£0(J')-62CI? J)*E(
NEXT J
NEXT I
BELETE G1 ? G2
GO TO 030
REM ****** ANÁLISIS CD.C)
REM Vá=E
BELETE E » C
DIM E(H)?C O(M) CN)

C=0
co=o
REM COMPROBACIÓN DEL PROGRAMA
IF F0=l THEN 1590
IF B1OO THEN 1590
GOSUB 1910
GOSUB 1680
 ANEXO i 6PR/YE3SOLU PAG, 39

1540 FOR 1=1 TO M

1550 C(I)=C<I)-E1CI)*E(I>
1500 CO(I)=CO<I)-E2<I)*E<I)
1570 HEXT I
.1580 GO TO 2030
1590 GOSUB 1970
1400 GOSUB 1680
1610 FOR 1 a TO M
1020 FOR J TO M
1630 )
1640 CO CI) =1O ( I ) -G2 C I > J ) #E C J >
1650 NEXT J
1660'NEXT I
1670 GO TO 2030
1080 REM *** SUB* FORMACIÓN DE la y Ms (REALES)
1690 FOR 1=5 TO Nl + 4
1700 READ #l ? i:y
1710 IF y<l)=l THEN 1750
1720 IF y(l)<>2 THEN 1760
1730 CCVÍ2))=V<5
1740 GO TO 1760
1750 E< y<2))=y(5
1760 NEXT I
1770 RETURN
1780 REM *** SUB FORMACIÓN DE (COMPLEJOS)
1790 SET DEGREES
1800 FOR J-5 TO Nl+4
1810 READ *l*i:y
1820 IF y<i)=l THEN 1870
1830 IF U(l)<>2 THEN 1890
1840 C(V<2))=y<5)*COS<VC7» (6 )

GO TO 1960
E < V < 2 >)=y < 5 > #CQS C y < 7 > > *1Q^y < ó >
E0<y<2) }~'v'í5}#SIN(y (7) )*10"y(6)
NEXT I
RETURN
REM *** SUB* LECTURA DE E1.E2 *^'*'
DELETE E1?E2 '
DIM E1(M)?E2<M>
READ *2.-6íEl
READ *2?7:E2
RETURN
REM *** SUB* LECTURA DE G1.G2 ***
DELETE G1 ? G2
DIM G 1 ( M ? M ) ? G 2 < M y M >
READ *2?2ÍG1
REAU *2?3;G2
RETURN
REM ****** FORMACIÓN DEL UECTOR Is ******
DIM CKN) íC2CN)

C2~0
FOR 1=1 TO N
 ANEXO 1 @PR/YE3SOLU PAG. 40
2080 FOR J=l TO M
.2090 X$ = SEG(A$> (I-im-í-fJ.l)
2100 T=VAL<X$)+1
2110 GO TO T OF 2120.2170.2150
2120 C1(I)=C1CI)-C(J)
2130 C2(I)=C2(I)-CO(J)
2140 GO TO 2170
2150 C1(J)=CKI)+C<J)
2100 C2(I)=C2(I)-fCO(J>
2170 NEXT J
2180 [Link] I
2190 IF H7=Q TREN 2220
2200 F'RINT @P9: " TÉRMINOS INDEPENDIENTES SON "-C1.C2
2210 REM LECTURA DE MATRICES DE COEFICIENTES
2220 DELETE Y - Y O
2230 DIM Y ( N r N ) í Y O < N r N )
2240 READ *2.-4íY
2250 READ *2r5ÍYO
2260 REM ****** SOLUCIÓN DEL SISTEMA DE ECUACIONES ******
2270 REM
2280 REM N — NUMERO DE ECUACIONES
2290 REM [Link] — MATRICES DE COEFICIENTES < N , N > REAL. IMAG
2300 REM ClfC2 — VECTORES DE TÉRMINOS INDEP, <N) — SOLUCIONES RE?IM
2310 REM
2320 Q2=1*OE-12
2330 FOR Q6-1 TO N-l
2340 REM **# BUSCAR RIMÓTE
2350 G5=SQR(Y<GórQÓ)*Y<GÓ,QÓ>4-YQ<Gó r GÓ)*YO<QóyQ6))
2360 Q4 = Q6
2370 FOR Q3=G6+1 TO N
23SO Q1=SQR< Y(Q3 ? Q6)#Y(Q3,Q6)+YO(03•Qó)*YO(Q3 , Q6 ) )
2390 IF Q5=>Q1 THEN 2420
2400 Q5=Q1
2410 Q4=Q3
2420 NEXT Q3
2430 IF GJ5=>Q2 THEN 2400
2440 PRINT 'MATRIZ DE COEFICIENTES ES SINGULARGGGG'
2450 END
2460 IF Q4-QÓ THEN 2630
2470 REM *** INTERCAMBIAR FILAS 04 Y .06
2480 FOR Q3=Q6 TO N . -
2490 Q5=Y(Q4.'Q3)
2500 Y(Q4jQ3)=Y(Qój03)
2510 Y<QÓjQ3?=G5
2520 Q5=YO(Q4.Q3)'
2530 YO <Q4,Q3)=YOC Q6 ?Q3)
2540 YOCQ6.Q3)-Q5
2550 NEXT Q3
2560 Q5-CKQ4)
2570 C1CQ4)~C1CQ6)
25SO C1CÜ6)=Q5
2590 Q5=C2(Q4)
2600 C2(Q4)=C2<Qó)
2610 C2(QÓ)=Q5
 ANEXO 1 (2PR/YE3SQLU P A G » 41

2620 REM #*# CEROS BAJO PIVOTE

2630 FOR 05=06+1 TO N
2640
G3 = Ql!fí<YCQ5íQÓ )#Y ( 06 < Q6 ) + YO ( 05 , 06 ) #YO < 06 , 06 > )
2660 Gl=Ql*(YOCQ5íG 6) #Y ( 06 , 06 ) -Y ( Q5 > Gó ) ^YO ( 06 * Q6 ) )
2670 FOR G4-QÓ+1 TO N
2680 QO=Y< 05,04) +03£Y(Q6rQ4)~Ql#YO(QÓ 7 Q4)
2690 YO(G5íG4)=YO(0 5 íG4)+G3!*:YOCQó>Q4)+Ql*Y(Q6iQ4)
2700 Y<Q5íG4)=QO
2710 NEXT G4
2720 00-C1 (05 )+G3*C1 (06)- Ql^C2(G6)
2730 C2Í05)"~* l^ 2(05 ) + Q 3*C2C Q6)+Q1*'C1(Q6)
2740 Cl (05)~Q 0
2730 NEXT 05
2760 NEXT 0 6
2770 IF SOR (Y ( N f N)# Y <N ? N)+ YOÍN ? N )#YO < N ? N ) )<02 THEN 2440
2780 REM *££ CALCULO DE INCÓGNITAS
2790 00=1 /( Y( N f N )*Y < N? N)+YO(N • N ) Y0(
* N* N) )
2800 01=00* (CKN)#Y (N r N)+C2CN)*YO (N.- N) )
2810 C 2 < N > = 00 *(C2 (N )*Y (Nr N)-Cl(N)*YO (HrN) )
2820 C1(N)= 01
2830 FOR GÓ=N -1 TO 1 STEP _ -i
2840 Q3 = 0
2850 Gl = 0
2860 FOR Q5= Q6 + Í TO N
2870 03-03+ Y (06 y 0 5))K Cl (05 )-YO (06 y05)*C2CQ5)
2880 Q1-Q1+ Y( 06-0 5)*C2 (05 )+YO (Q6j 05)*C1CQ5)
2890 NEXT 0 c- J
2900 07~1 /( Y( 06 y 06)*Y( Q6? 06) +YO ( ó? 0 G6)#YQ(Qóf Gó> )
2910 QO=Q7* < C CKQ ó)-03 )*Y (06r 06) +(C2 <Q6)-Ql)*YOCQ6jQÓ) )
2920 C2(Q6)~Q 7#< ( C2 (06 )-Ql)#Y (Q6? 06)•~(C.l(Q6)-Q3)*YO(Q6rQÓ) )
2930 CKQ6)= QO
2940 NEXT 06
2950 IF H7 = 0 THEN o 990
2960 F'RIN.T OP 9t ' JSOLUCIONVEC TOR COMPLEJO DE NODOS ( Vn ) 'jC
2970 PRINT 1 J J PARA CONTINUAR APLASTE (RETURN) ' ?
29SO INPUT X$
2990 X$='@PR/ALMS *
3000 01-42
3010 GO TO 920
 ANEXO 1 gPR/ALMS PAG, 42

1000 REM ****** LINEA 1000 2PR/ALMS T^ 'T- 'T- ** -

1010 02 = 42
¡020 GO TO 800
1.030 REM CALCULO DE CORRIENTES DE R A M A /ALMACENAMIENTO
1040
LOSO IF F1OF2 THEN 1690
.060 GO TO H8 OF 1630j1070
.070 IF FO-0 AND F=0 THEN 1090
LOBO GO TO 1520
L090 IF F1OO THEN 153O
U O O DELETE MI
.110 DIM VKÍ-S)
L120
L130 FOR 1=1 TO M
L140 FOR J=l TO N
LISO
.160
.170 GO TO T OF 1180-1210yl200
1180 V K D - V K D - C I C J )
.190 GO TO 1210
200 yiCI)=yi(I)-fCl< J)
210 NEXT J
220' NEX'T I
.230 Vl=yi+E
.240 IF D1OO THEN 1300
.250 DELETE El
.200 DIM El(h)
.270 READ *2í9ÍEl
.280 FOR 1=1 TO h
.290 yici)=Eia)*yi(i)
300 NEXT I
310 DELETE C3,E1
1320' DIM C3(M)
L330 C3=V1
L340 DELETE VI
L330 GO TO 1450
.300 DELETE G1 ? C3
>70 DIM G K M f M ) ? C3CM)
.380 READ *2?9íGl
L390 C3=0
.400 FOR 1=1 TO M
.410 FOR J=l TO M
.420 C3(I)=C3 íI)+G1<IfJ)*y1(J)
L430 NEXT J
',40 NEXT I
1450 REM
460 FOR 1=5 TO Nl+4
.470 READ *1,IíV
.480 IF U C D O 2 THEN 1500
490 C3 < y < 2 ) ) =C3 < V < 2 ) ) - V < S ) *10'
ISOO NEXT I
LSI O GO TO 1610
L520 IF Fl=0 THEN 1540
530
 ANEXO 1 (?PR/ALMS P A G * 43

1540 GOSUB 2320

1550 REM LLAMADA A SUBRUTINA CLCULO BE Ib PARA PR/YE2 PR/YE3
1500 IF H7^0 THEN 1630
.1570 IF NOTCF^O AND F^0> THEN 1620
1580 IF F1OF2 THEN 1620
1590 IF F1OO THEN 1620
:!.¿oo PRINT @P9: 'CORRIENTES DE RAMAS"-C3
1010 GO TO 1630
1620 PRINT @P9Í"CORRIENTES DE RAMAS' ? C3>C4
1030 REM X*** FIN DE CALCULO PARA UNA SOLA FREC*
:LÓ40 REM HS=1 SOLO Vn . H8—2 Vn f Ib
1650 REM CARGAR PROGRAMA DE IMPRESIÓN
loó O X$="í?pR/ESCRIBE*
1670 01-43
1680 GO TO 920
1690 REM CALCULO DE Ib PARA VARIAS FRECUENCIAS
170 O IF F0=0 AND F=0 THEN 1720
1710 GO TO 1750
1720 PRINT flL PROBLEMA CUYA ADMITANCIA DE ELEMENTOS ES SOLO REAL
1730 PRINT •J NO REQUIERE CÁLCULOS PARA VARIAS FRECUENCIAS '
1740 GO TO 1060
1750 REM
1700 READ
1770 GO TO O OF 1780* 1810.*1380í1340í1880rl880*1880
1780 A8=C1(I5>
1790 A9=C2CI5)
1300 GO TO 2070
1810 A8=C1<I5)-CKJ
1820 A9-C2(Í5)-C2<J
1830 [Link] 2070
1840 J5)+C2<- »
1850 AS=(C1(I5)*C1( ;;> + C2(Ii •)*C2<j:
1800 )#T
1870 GO TO 2070
1880 REM CALCULO DE LA CORRIENTE DE RAMA
1890 GOSUB 2320
1900 IF OO3 THEN 1940
1910 Mí? C3<I5>
1920 A9 C4CI5)
1930 GO TO 2070
1940 GO TO ü-4 OF 19SO?1960,2040
1950 REM
1960 T = l/ < C3 ( J5) #C3 (J5 ) +C4 ( J5 ) *C4 < J5) )
1970 IF QO5 THEN 2010
1980 A8=T*(ClC15)*C3C J5>+C2(
1990 A9=T*(C2(15 >*C3<J5)-C1C
2000 GO TO 2070
2010 A ci = T '¿i f n ~\i S ) ^í C 3 ( ü 5 ) -}" C 4 C15 "* & C 4 (-J
2020 A9==T*<C4(I5)*C3(J5)-C3(I5)#C4<J
203Q GO TO 2070
2040 T=1/<C1(J5)*C1CJ5)+C2CJ5)*C2<J5))
2050 A8=T*(C3<I5)*Ci(J5)+C4CI5)*C2<J5))
2060 A9==T* C C4 C15 ) #C1 ( J5 ) -C3 (15 ) *C2 C J5) >
2070 U$=D$&'G'
 ANEXO 1 E'PR/ALMS PAG, 44
2080 IF I1O1 THEN- 2180
2090 CALL 'FILE 1 ,D9 Í U$ ? X$
2100 IF XÍO1 " THEN 2130
2110 CRÉATE U*?2007 O
2120 OPEN U $ P 8 7 ' U B íX$
2130 IF P Q S < X * > " 0 OPEN'?130)=0 THEN 2150
2140 OPEN U$;Sj'U'rX$
2150 CALL -REWIND 1 t 8
.2160 WRITE *8:N27F3íA8*A9
2170 GQ TO 2200
2180 REM
2190 WRITE +8íF3 ? A8fA9
2200
2210 IF I1>N2 THEN 2270
2^20 IF ZO2 THEN 2240
2230 F3=F3#D
2240 X$='£PR/YE2VAF2'
2250 CU = 36
2200 GO TO 920
2270 PRINT "L VALORES DESEADOS CALCULADOS Y ALMACENADOS
2230 CLOSE 8
2290 03-1
2300 PRINT 'JJ PARA GRÁFICOS TECLA 6G'
2310 END
2320 REM SUBRUTINA DE CALCULO DE Ib (GPR/YE2? 0PR/YE3)
,¿. w*1~> O DELETE V l í V 2
2340 DIM V1CM)rV2(M>
2350 y 1=0
2300
2370 FOR 1=1 TO M ' -
2380 FOR J-l TO N
2390
2400 T=UALCX$)+1
2410 GO TO T OF 2420; 2470r 2450
2420
2430
2440 GO TO 2470
2450
2460
2470 NEXT J
2480 NEXT I
2490
2500 IF L8=0 THEN 2520
2510
2520 IF FOO1 AND Dl THEN 2770
2530 DELETE Glf(?2;C3
2540 DIM G 1 ( M ? M ) 7G2C ) íC3(M) ?C4(M)
2550
2500 C4-0
2570 IF F0=0 AND F=0 THEN 2680
2580 READ #2;2:G1
2590 READ *2j3ÍG2
2600 FOR 1=1 TO M
2610 FOR J=l TO M
 A N E X O 1 GJPR/ALMS P A G , 43

2620 C3 ( I ) ~C3 ( I ) -i-Gl < I ? J ) *V1 ( J )J >-(32 # l - >(2( I

J)?
2630 C4 < I ) =C4 < I ) +G1 < I y J > #V2 C J ) +G2J ( I
>*V1U) >
2640 NEXT J
2650 NEXT I
2660 DELETE G l > G 2 ? V l 7 V 2
2670 GO TO 2930
2680 READ *2;9:G1
2690 FOR 1=1 TO M
2700 FOR J=l TO M
2710
2720
2730 NEXT J
2740 NEXT 1
2750 UT"I ti
r** i r" "t* i"* r*^ H
L lü f ti G? 1 y iV* -i1 ? Tví /**!j¿
L

2760 GO TO 2930
2770 DELETE C37C4ÍE1..E2
2780 IF FOO THEN 2860
2790 DIM E 1 ( H ) ? C3 ( M ) y C4 ( M )
2800 READ #2í9¡£l
2810 FOR 1=1 TO M
2820
2830
2840 NEXT I
2850 GQ TO 2930
2860 DIM El CM)? E2CM) ? C3 < H ) >C4(M)
2870 READ *2í¿ÍEl
2880 READ *2?7tE2
2890 FOR 1=1 TO M
2900 C3 C I ) «El < I ) #V1 C I ) ~E2 ( I ) *V2 C I )
2910 C4 í I ) ~E1 ( I ) JKV2 ( I ) +E2 < I ) #V1 C I )
2920 NEXT I
2930 REM
2940 IF L8«0 THEN 3030
2950 SET DEGREES
2960 FOR 1=5 TO Nl+4
2970 READ # 1 ? I Í M
2980 IF y<l)<>2 THEN 3010
2990 C3CVC2) )=C3 C V C 2 ) )-y(5)#CQS<•y ( 7 )
3000 C4 CU( 2 ) )~C4í V C 2) )-V (5)*SINÍ y < 7 )
3010 NEXT I
3020 GO TO 3090
3030 FOR 1=5 TO NH-4
"3 0 4 0 READ #l?i;V
3050 IF V(l)<>2 THEN 3080
3060 w PX
vi ( \U \«(¿ 1 a'(2> )-V(S)^10"tvJ < 6 )
J "1 W)O=:P"X
3070 REM C4-C4
3080 NEXT I
3090 RETL1RM
 ANEXO 1 @PR/YE2VAF2 PAG, 46
1000 REM ************ LINEA 1000 @PR/YE2VAF2 *************
1010 02-30
1020 GO TO SOO
1030 REM
1040 DJ.M E 2 C M )
1050
1060 IF ZOl THEN 1080
.1070 F3=F1+<I1-1)*C
:l.080 W = 2*F'I*F3
Ü.090 FOR 1 = 5 TO Nl + 4
;L100 READ * l í I t V
1110 GO TO V C D - 3 OF 1130-11 0
1120 GO TO 11 00
1130 £ 2 < Y ( 2 > ) = E 2 < V ( 2 > > - ! / < U * V < 5 ) * 1 0 ~ V í ó; )
1140 GO TO 11 00
1150 E2(U<2> )-E2(V<2»+U*Y(5)#10~VCó)
11 ¿O NEXT I
1170 WRITE # 2 r 6 : E 2
1180 IF H7 = 0 THEN 1200
1190 PRINT @P9: "UECTOR DE ADMITANCIA IMAGINARIO 1 ? E2
1200 GO TO DI OF 1250 , 1490 , 1250? 1690 , 1690
1210 GO TO Ül-5 OF 1690, 1690 • 1990 7 1250, 1490
1220 GO TO Dl-10 OF 1250,1690,1090,1690*1090
1230 GO TO 2050
1240 REM Dl=l FCCC CB)
1250 CÉLETE U
1260 DIM U ( M , M >
1270 U"0
;L280 FOR 1 = 1 TO M
1290 IKI, I)
1300 NEXT 1
1310 FOR j= TO Nl+4
1320 READ #l j:V
1330 IF V(l) >8 THEN 1350
1340 U (U (4) , <3) }-UCU(4)? (5 ) íó)
1350 NEXT J
1300 DELETE G2
1370 DIM G2(M,M>
1380 G2-0
1390 FOR .1 = 1 TO M
1400 FOR J=l TO M
1410 G2<I, J)=U(I, )*E2(J)
1420 NEXT J
1430 NEXT I
;l.44Q' CÉLETE U
1450 GO TO DI OF 205071470^1490
1460 GO TO Dl 8 OF 1990,1470,1490
1470 PRINT 'GGG ERROR Dl= 'ÍD1
1480 END
1490 CÉLETE U1?U2,R
1500 DIM U 2 C M / M ) r R ( M r M ) í U K M í M
1510 U2=0
1520 R=0
1530 FOR 1=5 TO Nl+4
 ANEXO 1 (2PR/YE2VAF2 P A G » 47
1540 READ *•1; ir y
1550 IF V i l )<>9 THEN 1570
1560 R ( U C 4r )V (3 ))= l v l C5) #10~U(6) '
1570 NEXT I
1580 FOR 1 =1 TO M
1590 FOR J= 1 TO M
1600 U2 C I • J > = R < i f J)#E2(J)
1610 NEXT J
1620 NEXT I
1630 DELETE R
1640 REM LECTURA Ul
1650 , READ *2í 81 u 1
1660 GOSUB 2710
1670 IF Dl = 2 ÜR 01=10 THEN 2340
1680 GO TO 2480
1690 REM EXISTE UNA FUENTE FUCO (G) 01=4, , 6 ^ 7
1700 DELETE G2
1710 DIM G2 ( M ? M)
1720 G2 = 0
1730 FOR 1 =1 TO M
1740 G2CI, I)=E2 < I)
1750 NEXT I
1760 GO TO ni -3 OF 2650,1800,149071300
1770 GO TO 01--1 i OF 1990? 1800r 1490; 1800
1780 PRINT •GGGERRORD1= " ÍD1
1790 END
1800 REM CALCULO Gl+BKYE ni=5?7f 13í
1810 DELETE B
1820 DIM B( MÍ M)
1830 B=0
1840 FOR 1 =5 TO Nl+4
1850 READ * i.-i: y
ISóO IF U C 1 )<>8 THEN 1830
1870 B C U ( 4 r)V < 3 > >=y<5):*:iO~V<6)
1880 NEXT I
1890 FOR 1 =1 TO M
1900 FOR J = 1 TO M
1910 G2ÍI- J )=G2 (Ir J ) + B C I - J)#E2<J>
1920 NEXT J
1930 NEXT I
1940 DELETE B
1950 GO TO 01 -4 OF 2650-1970?1490
1960 GO TO IU-1 2 OF 1 990 7 1 970 y 1490
1970 PRINT •GGG ERROR Dl= ' í DI
1980 END
1990 REM ni — 8
2000 DELETE U
2010 DIM U C M i M )
2020 GO TO ni -7 OF 2060r2170
2030 GO TO Dl-1 1 OF 217072170
2040 PRINT •GGGERRQR Dl= ' J D 1
2050 END
2060 DELETE 02.. u
2070 DIM G2 C M r M )j U (M i M)
 ANEXO 1 8PR/YE2VAF2 PAG* 48

2080 G2=0
2090 GOSUB 2310
2100 FOR 1 =1 TO H
2110 FOR J =1 TO M
2120 G2 < I * J )=:£2(I)*U(Ií J)
2130 NEXT J
2140 NEXT I
2150 DELETE U
2100 GO TO 2650
2170 DELETE G'4? U
2180 D I M G 4(M .'M )? U C M ? M )
2190 G4~0
2200 GOSUB 23 10
2210 FOR 1 =1 TO H
2220 FOR j =1 TO M
2 230 FOR K=l TO M
2240 G4 <1' 7 J )= G4 (I ? J ) +G2 ( I y K ) #U ( K .- J )
2250 NEXT K
22éO NEXT J
2270 NEXT I
22SO G2=G4
2290 DELETE G 4 * U
2300 GO TO 26 50
2310 REM LEC TURA U
2320 READ * o }8Í U
2330 RETURN
2340 DELETE Gl- G2, El
2350 niM Gl < H - M) .- G2(M? H) ? El (M)
2300 Gl = 0
2370 G2~0
2380 READ *2: 61 El
2390 READ *2? 71 E2
2400 FOR 1 =1 TO M
2410 FOR j =1 TO M
2420 Gl ( I í J) =El í I ) #01 < I ? ü ) -E2 C I ) #U2 C 1 1 J ).
2430 G2CI, J )z^ El < I ) #U2 < I í J ) + E2 ( I ) #U1 < I ? J )
2440 NEXT J
2450 NEXT I
2460 DELETE U 1 1 U 2r El jE2
2470 GO TO 2080
2480 DELETE G•7 G 4 .' Gl ? El f E2
2490 DIH G3 (M ?M) f G4(MrM> jGKMjM)
2500 G3-0
2510 G4~0
2520 READ *2í9: Gl
2530 FOR j =1 TO M
2540 FOR J = 1 TO H
2550 FOR K= 1 TO M
2500 G3 C I , J)= G3 (1 7 J ) +G1 ( I ? K ) #U1 ( K , J ) -G2 C I > K ) #U2 ( K r J )
2570 G4 ( I - J)a:G4 (1 7 J ) +G1 C I r K ) *U2 ( K j J ) 4-G2 ( I , K ) #U1 C K , J )
25SO NEXT K
2590 NEXT J
2000 NEXT I
2010 G1-G3
 ANEXO 1 0PR/YE2VAF2 PAG, A9
2020
2630 DELETE G3 ? G4rUl,U2
2040 GO TO 2080
2050 X*='@PR/YE3JNCini'
2660 01 = 33
2670 GO TG 920
2080
.2090 01=32
2700 GO TO 920
2710 REM *M SUB, GJORUC
2720 REM
2730 REM N — NUMERO DE ECUACIONES
2740 REM YíYO — MATRICES DE COEFICIENTES (N>N) REAL.* IMAG
2750 REM ClíC2 — VECTORES DE TÉRMINOS INDEP. (N) --- SOLUCIONES
2700 REM
2770 Q2=ltOE-12
2730 FOR Qó = l TO N-l
2790 REM ### BUSCAR PIVOTE
2300 G5=SGR(Y<QóyGÓ>#YCG6,GÓ>4-YO<GÓjQÓ)*YO(G6 ? QÓ> )
2810 G4-QÓ
2320 FOR Q3 = GÓ-hl TO N

2340 1F G5=>G1 THEN 2870

2850 G5=Q1
2800
2370 NEXT Q3
2880 IF G5=>G2 THEN 2910
2890 PRINT 'MATRIZ DE COEFICIENTES ES SINGULARGGGG'
2900 END
2910 IF. Q4=QÓ THEN 3080
2920 REM *** INTERCAMBIAR FILAS Q4 Y Qó
2930 FOR Q3=Qó TO N
2940
2950 YCG4íQ3)=Y(QófQ3)
2900
2970 Q5=YO(Q4jQ3)
2930 YO(Q4jQ3)-YO(GójQ3)
2990 YO(GóíQ3)=Q5
3000 NEXT Q3
3010
3020 CKQ4)=Cl(Qó)
3030
3040
3050
3000
3070 REM W& CEROS B A J O P I V O T E
3080 FOR Q5-Q6+1 TO N
3090 Ql=-l / ( Y < G6 . 0.6 ) )KY C Qó r Qó ) +YO ( Q6 - Qó ) ^:YO ( Qé • Qó ) )
3100 {33=01* C Y < Q5 f Qó ) ücY ( Qó r Qó ) +YO ( Q5 , Qó ) * YO C Qó , Qó ) )
3110 Gl-OU*: C YO C Q5 7 Qó ) #Y < Qó j Qó ) -Y C Q5 ? Qó ) #YO < Ció T Qó ) )
3120 FOR Q4=QÓ+1 TO N
3130 QO=YCQ5íQ4)+Q3*Y(QójQ4)-Ql*YO<QÓjQ4)
3140 YOÍQ5?Q4)=YO<Q5íQ4).+Q3*YO(QÓ/Q4)+Ql*Y(QófQ4)
3150
 ANEXO 1 Í?PR/YE2UAF'2 PAG* 30

3160 NEXT Q4
3170 QO=C1(Q5
3180 C2 < 05 > =02(05)+03*02 < 06 )+01*CKG6)
3190 C1 ( G 5 ) ~ Q O
3200 NEXT 05
3210 NEXT Qó
3220 IF S0R(Y(NrN)*Y(N?N)+YO(NíN)*YO<NjN))<Q2 THEN 2890
3230 REM ##* CALCULO DE INCÓGNITAS
3240 00=1/(Y(NfN)#Y(NfN)+YO(N-N)*YO(N,N))

3260 C2 C N)=00* í C2 < N)#YC N,N)-C1C N)#YO(N,N))

3270 C1(N)=01
3280 POR 0Ó=N-1 TO 1 STEP -1
3290 03=0
3300 Ql=0
3310 FOR 05=00+1 TO N

3330 01 = 01+Y(Gó yQ5)*C2(05)+YO< Q6 ? 05)#01< 05 >

3340 NEXT Q5
3350 07=1/C Y(Qó,00)*Y(QótQó)+YOC Gó,00)#YO<Qó? Qó) )

3370 C2CQÓ)=Q7#f(C2<Gó)-01)*YíGó•Qó)~(Cl<Qó)-03)*YO(06 / Qó)

3380 CKQ6)=QO
3390 NEXT Qó
3400 RETURN
 ANEXO 1 @PR/YE3JNCID1 PAG, 51
1000 REM ********: LINEA 1000 0PR/YE3INCID1
1010 02-33
'.020 GO TO 800
L030 REM CORRECION DE MATRIZ DE ADMITANCIA RAMAL SI EXISTE IND, MUTUAS
1040 REM Y ALMACENAMIENTO
1050 REM PRINT @P9:' MATRIZ YB SIN MODIFICAR ' , ' Gl V- Gl, ' G2 V- G2
L060 IF FOO1 THEN 1430
.070 IF D1OO THEN 1260
LOSO REM CORRECION DE YB CUANDO ES MATRIZ DIAGONAL COMPLEJA
1090 DELETE G2
L100 IHM G2(M,M)
L110 G2=0
U 20 FOR 1=1 TO M
L130 G2<IrI)=E2(I)
U 40 NEXT I
1.150 DELETE E2
L160 GO TO 1270
U 70 DELETE L
U 80 DIM L(M)
..190
1.200 FOR 1=5 TO NI+4
1210 READ *líi:V
L220 IF UCDO4 THEN 1240
1230
.240 NEXT I
1.250 RETURN
L260 REM CORRECCIÓN DE YB CUANDO ES UNA MATRIZ COMPLEJA
1270 GOSUB 1170
..280 FOR 1=5 TO Nl+4
L290 READ *l.i;y
.300 IF U(!)<>/ THEN 1420
L310
..320
L330 T5=L(V C 4)? #L(V(3))- V(5)*V< 5)*10"í2*V(ó > >
1340 IF ABSCT5)>1 ( OE~16 THEN 1370
1350 PRINT 'J EXISTE ERROR EN LOS DATOS DE INDUCTANCIAS (REVISE)
1300 END
570 G2 í V(3)fV(4))=-<T4*V(5)*10~V(ó)*L C V(4)))/T5
1380
1390 T5~L < V(3 » *L(V(4 »/T5
1400 G2(V<3)yVC3))=T2*T5
L 410 G2(V(4)rV(4))=T4*T5
NEXT I
1430 REM
1.440 REM CALCULO DEL A*YB
1450 DELETE A$ ? Y2
i. 460 DIM A $ ( M * N > y Y 2 ( N T M )
L470
1480 READ #2.'i:A$
.490 IF FO-1 THEN 1630
'.500 IF D1OO THEN 1630
;1510 FOR 1=1 TO N
1520 FOR J=l TO M
1530
 ANEXO 1 <?PR/YE3INCini PAG* 52
1540
1550 GO TO T OF 15607Í590f 1580
1560
1570 GO TO 1590
1580 Y2CI, J)=E2(J>
1590 NEXT J
1600 NEXT I
1610 DELETE E2
1620 GO TO 1820
1630 REM PRODUCTO A*G1=Y1 A*G2=Y2
1640 DELETE E2
1650 FOR 1=1 TO N
1660' FOR J=l TO M
1670 FOR K=l TO M
1680
1690
1700 GO TO T OF 1710. 1740r 1730
1710 Y2<I? J)=Y2< I 7 J)-G2<K f J)
1720 GO TO 1740
1730 Y2(I? J)=Y2í I, J)+G2<Kj J)
1740 NEXT K
1750 NEXT J
1760 NEXT I
1770 REM ALMACENAMIENTO Yb
1780 WRITE #2>3:G2
1790 IF H?=0 THEN 1810
1800 PRINT £P9: ' JMATRIZ YB PARTE IMAGINARIA ' , G2
1810 DELETE G2
1820 REM
1830 REM Yi*A' DONDE A'= TRANSPUESTA DE A
1S40 DELETE YO
1850 DIM Y O C N í N )
1860 Y0=0
1S70 FOR 1=1 TO N
1880 FOR J=l TO N
1890 FOR K=l TO M
1900 X$=SEG< A$Í < J-DfcM+Ky 1)
1910 T-UAL (X$) +1
1920 GO TO T OF 1930? 1960? 1950
1930 YOÍIj J)=YO< Ií J)-Y2<IjK)
1940 GO TO 1960
1950 Y O < I ? J)=YOC I/ J)+Y2CIíK)
1960 NEXT K
1970 NEXT J
1930 NEXT I
1990 IF H7=0 THEN 2020
2000 PRINT @P9:'JMATRIZ Y2=A*E2<G2) ' Y2
2010 PRINT @P9í'JMATRIZ DE INCIDENCIA ALMACENADA EN DISCO
2020 REM
2030 URITE *2í5íYO
2040 IF H7=0 THEN 2060
2050 PRINT (3P9: ' JMATRIZ DE ADMITANCIA DE NODOS PARTE IMAGINAR!A',YO
2060 DELETE YO
2070 Xí>-D(3PR/YE3SOLU'
 OEÓ 01 03 06O£
fr£ = TO OSOS

T OX3NM
 ANEXO 1 0PR/GRAFICOS PAG»

1000 REM #*# LINEA 1000 0PR/GRAFICOS (TECLA ó)

1010 02 = 9
:L020 GO TO 800
loso REM #*# GRÁFICOS DE RESULTADOS EN FUNCIÓN DE FRECUENCIA ##*
1040 DELETE XOíYO
1050 X0=l,792
ILOÓO Y0=2*816
1070 PRINT -L GRÁFICOS EN FUNCIÓN DE LA FRECUENCIA'
10SO IF D*="* THEN 1130
1090
1100 PRINT 'JARCHIVO ESI ' í X í í 1 (SI O N0)í G'í
1110 INPUT X*
.1.120 IF X*='S' OR X$="SI' THEN 1230
1130 PRINT "JJ NOMBRE DEL ARCHIVO DE DATOS: '?
1140 INPUT B$
1150
11 00 CALL 'FILE' jO,Ii*jX*
1170 IF X$<>"" THEN 1230
1180 PRINT 'J NO EXISTE ARCHIVO 1
1190 PRINT 'J DESEA EJECUTAR OTRO ARCHIVO? (SI O NO) *í
1200 INPUT X$
1210 IF XÍ>=-S' OR X*='SI' THEN 1130
1220 GQ TO 100
.1230 OPEN n $ í l í ' R w ? X*
1240 READ *líliZ*
1250 READ *lf2íOfI5íJ5?Z
.1200 READ *lí4:Fl^F2
1270 GLOSE 1
1280 IF F2>F1 THEN 1340
1290 PRINT -JERROR; NO SE PUEDEN HACER GRÁFICOS SOLO .PARA UN
1300 PRINT 'VALOR DE FRECUENCIAG'
:L310 PRINT 'JJRETURN PARA CONTINUARG"
1320 INPUT XS
1330 GO TO 100
1340
1350 OPEN G $ í 1 r 'R"rX*
1360 READ *i:N2
1370
1380 KS=1/LGT(F3)
1390
1400 REM POSICIÓN PARA FREC, LOG * J = K8 # LGTCF) -f K9
1410
1420 IF NOT(0=3 OR 0=0 OR 0=7) THEN 1440
1430

1450 IF NOT<0=5 OR 0=6) THEN 1470

1460
1470
1480
1490
1500 IF 0=1 OR 0=3 THEN 1600
1510 IF 0<>2 THEN 1540
1520
1530 GO TO 1550
 ANEXO 1 0PR/GRAFICGS PAG, 55

1540 C$~C$& 1 / '

1550 C$~C$£ G
1560 Xí-STR (J5)
1570 X$=REP < 1 * ? 1 ? 1 )
1580 C*-C$S X $
1590 C$s=C$& • ) "
1600 REM c$ = NOMBRE DEL GRÁFICO
1610 PRINT USING 1620:¿$
1620 IMAGE /4X'PR OBLEMA: •EA
1630 PRINT " J 1 - GR iFICO DE PARTE REAL 1
1640 PRINT • 2 - GRA 'ICO DE PARTE IMAGINARIA'
1 ASO PRINT • 3 - GRA "ICO DE MAGNITUD'
1660' PRINT • 4 - GRA "ICO DE VALOR DE MAGNITUD EN DB* '
1670 PRINT • 5 - GRA 'ICO DE MAGNITUD EN BB * r CON RESPECTO
1680 PRINT • ALGÚN V A L O R 1
1690 PRINT • 6 - GRA. ICO DE ÁNGULO DE FASE 1
1700 PRINT J 9
CLASE DESEADA DE GRÁFICO DE "?C$Í G'
1710 INPUT Kl
1720 IF Kl< 1 OR Kl >6 THE. ! 1700
1730 DATA " PARTE REAL 1 ,a PARTE I M A G I N A R I A 1 ? MAGNITUD'
1740 DATA 'MAGNITUD EN. d •MAGNITUD EN dB -'ÁNGULO DE FASE1
1750 RESTORE 1730
1760 FOR J =1 TO Kl
1770 READ X$
1780 NEXT J

IF K1O5 THEN 1900

PRINT 'J FRECUENCIA DE REFERENCIA PARA Q DB,
PRINT "('ÍF1.; 1 - " Í F 2 Í 1 Hz>: 'í
1840 INPUT K3
350 IF K3<F1 OR K3>F2 THEN 1820
1860 C$=C$&< (O dB = '
L870 X$=STR(K3)
LSSO
1.890 C$=C*S' Hz)"
1900 PAGE
L910 GO TO 2090
L920 REM #*# SUB, DE WINDOW Y DE VIEWPORT
1.930 WINDOW W l f W 2 í U 3 í U 4
1.940 IF P8"32 THEN 1970
L950 VIEWPORT 127145-10.96
1.960 GQ TO 1980
1970 VIEWPORT 12?1257lOy9ó
1980 RETURN
..990 REM #*# SUB* HE DISPOSITIVO DE GRAFIZACION
2000 PRINT 'JDESEA GRÁFICO EN PAPEL? (SI O NO): 'í
203.0 INPUT X$
2020 IF X$='S" OR X*="SI' THEN 2050
20 o O P8-32
2040 GO TO 2080
2050
2060 PRINT ALISTE GRAFIZADOR (RETURN PARA CONTINUAR)G Í
2070 INPUT
 ANEXO 1 RPR/GRAFICGS PAG» 5¿
2080 RETURN
2090 REM
2100 N-N2-1
2110 REM N = NUMERO DE INTERVALOS
2120 DELETE X l ? Y l j X 9 í Y 9
2130 DIM XKN2) fYlíN2>
2140 SET DEGREES
2150 FOR J=l TO N2
2100 READ # 1 ; X 1 ( J ) > X 9 í Y 9
2170 GO TO Z OF 2190,2130
2180 X1(J)=J
2190 GO TO Kl OF 2200y2220y2240i2240í2240?22óO
2200 YKJ)=X9
2210 GO TO 2310
2220 Yl<J)=Y9
2230 GO TO 2310
2240 Y1CJ)=SQR(X9#X9+Y9*Y9>
2250 GO TO 2310
2200 Y1(J)=0
2270 IF X9>0 AND Y9=0 THEN 2310
22SO Y1(J)=180
2290 IF X9<0 AND Y9=0 THEN 2310
2300 Y1< J)=SGN(Y9)*ACS(X9/SQR C X9*X9 + Y9*Y9))
2310 NEXT J
2320 GO TO Kl OF 2420j2420?2420.'2400j233072420
2330 GOSUB 4960
2340 GO TO Z OF 2350?2370
2350 X=K3
2360 GO TO ;o
2370
2380
2390 Y1=Y1/P
2400 Y1=LGT(Y1)
2410 Y1=20*Y1
2420 REM ^^* XI, Yl DATOS PARA GRÁFICO (SIN INTERPOLACIÓN)
2430 REM ##& CALCULO DE MÍNIMO Y MÁXIMO HE ORDENADAS
2440 REM MI = MÍNIMO
2450 REM M2 = MÁXIMO
2460 Ml^[Link]+100
2470 M2=-M1
2480 FOR J=l TO N2
2490 MI-MI MIN Y1CJ)
2500 M2=M2 MAX Y 1 C J )
2510 NEXT J
2520 PRINT LGRAFICO DE
2530 PRINT JPARAMETROS PARA EL GRAFICOJ"
2540 PRINT F m i n == •; F i
2550 PRINT FfTlBX = ' ? F2
2560 PRINT USING 30X"' KKFmin PARA GRÁFICO
2570 INPUT F4
2580 PRINT USING :ox*mF(ri3x PARA GRÁFICO
2590 INPUT F5
2600 PRINT
2<610 PRINT Y ni i n M1
 ANEXO 1 0PR/GRAFICOS PAG. 57
2620 PRINT ' YITIBX - ' ÍM2
2030 PRINT USING -30X" •30X''KKYmin
"KKYiriin PARA GRÁFICO: a " , S " :
2040 INPUT U3
2050 PRINT USING "30X a "Ymsx PARA GRÁFICO: '
2600 INPUT W4
2661 PRINT 'J 1 - EJES CON MARCAS'
2662 PRINT ' 2 - RETICULADO"
2663 PRINT ' 3 - SOLO CURVA (SIN EJES NI RETICULADO)'
2664 PRINT 'J CLASE DESEADA: •;
2665 INPUT W9
2 6 6 6IF NOKW9-1 OR W9=-2 OR THEN 2664
2667 IF W9=3 THEN 2720
2670 GO TO Z OF 2680.'2700
2680 PRINT USING VÓX 1 "ESPACIO ENTRE MARCAS PARA EJE HORIZONTAL:
2090 INPUT W7
2700 PRINT USING VÓX'"ESPACIO ENTRE MARCAS PARA EJE VERTICAL:
2710 INPUT W5
2720 PAGE
2730 RE'M XXX SIGUE GRÁFICO DE PUNTOS SIN INTERPOLACIÓN
2740 PS-32
2750 GO TO Z OF 2760*2800
2760 REM #$# FREC, LINEAL
.2770 W1 = F4
2780 W2-F5
2790 GO TO 2S30
2800 REM ### FREC, LOG
2310 W1 = KS'*LGT(F4)+K9
2820 W2=K8*LGT(F5)+K9
2830 GOSUB 1920
28 4 O MOVE gpeiXICl),Y1(1)
2850 DRAW {?F'8:XlrYl
2S60 IF P3O32 THEN 2920
£870 HOME
2880 GOSUB 1990
2S9Ü IF P8=32 THEN 2920
2900 PAGE
2910 GO TO 2750
.2920 PRINT 'DESEA MAS PUNTOS (SI O NO): 'í
O O "i A INPUT X*
¿- / •uv
2940 IF NOTCX*='S' OR X*="SI a > THEN 3310
2950 REM *** SIGUE INTERPOLACIÓN Y LUEGO GRÁFICO
2960 PRINT 'JNO+ DE PUNTOS: '?
2970 INPUT N9
2980 PAGE
2990 P8-32
3000 DELETE X2-Y2
3010 DIM X2CN9),Y2CN9)
3020 GOSUB 4960
3030 GO TO Z OF 3040j3070
3040 U1 = F1
3050
3060 GO TO 3090
3070
3080 U2-N2
 ANEXO 1 SPR/GRAFICGS PAG» 58

3090 U3=(U2-U1>/CN9-1)
3100 1= 1
3110 X2<1>=X1(1)
3120 X2(N9)=X1(N+1)
3130 FOR J=l TO N9
3140 IF J=l OR J-N9 THEN 3160
3150 X2(J)=U1+<J-1)#U3
3160 IF X1CIX = X2(J> AND X2 ( J ) <=X1 C 1 + 1 ) THEN3190
3170 1=1 + 1
3180 GO TO 3160
3190 n=X2CJ)-XlCI)
3200 Y2C J)=C( 1 y I )+D*(CC2í I )+fl#<C(3, 1 >+B#C(4 ,
3210 NEXT J
3220 GOSUB 1920
3230 hOVE @P8ÍX2<1) ,Y2(1)
3240 IlRAW @P8:X2,Y2
3250 IF P8O32 THEN 3310
3260 HOME
3270 GOSUB 1990
3280 IF P8=32 THEN 3310
3290 PAGE
3300 GO TO 3220
3310 GO TO W9 OF 3340-3340*4690
3340 GO TO Z OF 3790-3350
3350 REM #*### EJE HORIZONTAL CLOG)
3360 E4=LGT(F4+l*OE-ó)
3370 IF ABSCE4-INT(E4»<1>OE-12 THEN 3390
3380 E4=INT(E4)
3390 E5=LGT(F5)
3400 IF ABS(E5-INT(E5»<1.0E~12 THEN 3420
3410 E5=INTCE5>+1
3420 E6=E5-E4
3430 GOSUB 1920
3440 HOME @P8ÍW1,W3
3450 FOR J-E4 TO E5~l
3460 J4~10"J
3470 J2=2
3480 IF JOE4 THEN 3500
3490 J2=l
3500 FOR J1=J2 TO 10
3510 X=J1*J4
3520 IF X<F4 OR X>F5 THEN 3760
3530 GOSUB 1920
3540 Q1=K8#LGT(X)+K9
3550 GO TO W9* OF 3560.3680
3560 HRAW @P8:Q1,W3
3570 SCALE 1,1
3580 RURAW (?P8:Oj 1
3590 RMOVE @PS:O,-I
3600 IF NOTCX=F4 OR X=F5 OR Jl=10) THEN 3760
3610 X$"STR(X)
3620 X*=REP( ' ' > 1,1)
3630 RMOVE @P8:-LEN<X*)/2#XOí-l*l*YO
3640 PRINT SPSÍX*?
 ANEXO 1 GPR/GRAFICQS PAG, 59
3650 'GOSUB 1920
3660 MOVE GP8tQljW3
3670 GO TO 3760
3680 MOVE @ P 8 ¡ C U , W 4
3690 DRAIO G P 8 : Q l j U 3
3700 SCALE l.'l
3710 IF NOT(X=F4 OR X=F5 OR Jl=iO) THEN 3760
3720 X*=STR(X)
3730 X*=REP< " tlr 1)
3740 RMOVE GP8;-LEN(X*)/2*XOf-l»l*YO
3750 PRINT GF'8;X$Í
3760 NEXT Jl
3770 NEXT J
37SO GO TO 4110
3790 REM ***** EJE HORIZONTAL (LINEAL)
3800 GOSUB 1920
3810 IF Wl=>0 THEN 3840
3820 IF U2OO THEN 3840
3S30 GO TO 3890
3840 MOVE G P 8 í W l í W 3
3850 FOR J=W1 TO W2 STEP W7
3860 GOSUB 3930
3870 NEXT J
3880 GO TO 4100
3890 MOVE GP8tOfW3
3900 FOR J=0 TO Wl STEP -W7
3910 GOSUB 3980
3920 NEXT J
3930 MOVE GP8ÍOjU3
3940 FOR J=0 TO U2 STEP W7
3950 GOSUB 3980
3960 NEXT J
3970 GO TO 4100
3980 REM ##* SUB. GRÁFICO DE MARCAS EN EJE HORIZv
3990 DRAW GP8;JíU3
4000 SCALE l í l
4010 RDRAW GPSIOíl
4020
4030 X*=STR(J)
4040 X*~REP('"flrl)
4050 RMOVE GP8í-LEN(X$)/2*XOj-l*l*YO
4060 PRINT @PS¡X*Í
4070 GOSUB 1920
4080 MOVE GPStJiWS
4090 RETURN
4100 REM
4110 REM ***** EJE VERTICAL (LINEAL)
4120 W6=0
4130 GOSUB 1920
4140 IF W3=>0 THEN 4170
4150 IF W4<=0 THEN 4170
4160 GO TO 4220
4170 MOVE @P8:WlíW3
4180 FOR J=W3 TO W4 STEP W5
 ANEXO 1 SPR/GRAFICOS PAG. 60
4190 GOSUB W9 OF 4370^4510
4200 NEXT J
4210 GO TO 4660
4220 MQVE @P8:UlfO
4230 J = W5
4240 J=J-W5
4250 IF J=>W3 THEN 4270
4260
4270 GOSUB U9 OF 4370 .* 4510
4280 IF J>W3 THEN 4240
4290 MOVE @P8:U1,Q
4300 J--W5
4310 J-J+W5
4320 IF J<=W4 THEN 4340
4330 J-W4
4340 GQSUB. W9 OF 4370-4510
4350 IF J<W4 THEN 4310
4360 GO TO 4660
4370 REM ### SUB, GRÁFICO MARCAS EJE VERTICAL
4380 DRAW e P Q í W l y J
4390 SCALE 1..1
4400 RDRAW @P8ílfO
4410 RMÜVE @P8:~1-
4420
4430 Xí>=REPC ' ' -1 ? 1)
4440 IF-L£N<X$)=Ó THEN 4470
4450 X$= a 'SX$
4460 GG TO 4440
4470 PRINT GPS: USING 'FA'tXS
4480 GOSUB 1920
4490 MOVE G P S r U l i J
4500 RETURN
451.0 REM W* SUB* LINEAS HORIZONTALES
4520
4530 DRAW ©PSÍ
4540 SCALE 1 - 1
4550 RMOVE GPS
4560 X*=STR<J)
4570
4 58 O WÓ=WÓ MAX LEN(X*>
4590 IF LEN(X$)=Ó THEN 4620

4010 GO TG 4590
4620 PRINT @P8: USING 'FA^íX*
4630 GOSUB 1920
4640
4650 RETURN
4660 REM
4670 HOME
4680 PRINT ' JJJJ'
4690 IF P8=32 THEN 4910
4700 PRINT 'DESEA LEYENDAS? (SI O NO): '?
4710 1NPUT X^
4720 IF NOT<X*= 1 'S' OR X$='SI") THEN 4910
 ANEXO 1 epR/GRAFICOS PAG. 61
4730 X$='FRECUENCIA / Hz1
4740 GOSUB 1920
4750 MOVE @P8;(W1+W2)/2?W3
4760 SCALE 1,1
4770 RMOVE @P8í-LENC X*)/2#XO,~2 * 2#YO
4780 PRINT @P8:X*
4790 PRINT gP8,25:90
4300 GOSUB 1920
4S10 MOVE @P8:U1,(W3+W4)/2
4S20 SCALE l í l
4830 RMOVE @P8:-l*5-Wó*XOj-LEN<C*>/2#XO
4840 PRINT @F8ÍC*
4350 PRINT @P8,25:0
4SéO GOSUB 1920
4S70 MOVE @P8:(W1+U2)/2,W4
4880 SCALE 1,1
4890 R M Ü V E @P8:-LEN(Z$)/2#XO»O *2*YO
4900 PRINT @P8:Z*
4910 CLOSE
4920 IF P8O32 THEN 4950
4930 PRINT -RETURN PARA CONTINUARG'
4940 INPLIT X$
4950 GO TO 100
4960 REM ### SUB DE CALCULO DE POLINOMIOS DE INTERPOLACIÓN
4970 DELETE C
4980 1UM C<4,Nil)
4990 C-0
5000 REM ALMACENAR Y l C I ) EN C C 1 r I )
5010 FOR 1=1 TO N+l
5020
5030 NEXT
5040 REM SE INICIALIZA DERIVADA EN PUNTO INICIAL Y EN PTO* FINAL
5050 REM 1=1 APROXIMACIÓN LINEAL PARA DERIVABAS
5000 REM 1-2 APROXIMACIÓN CUADRÁTICA PARA DERIVADAS
5070 I ~~ *1-O
M

5080 GO TO I OF 5090, 120

5090 Y1<2)-Y1C1) )/<XK2)-Xl(D)
5100 = CYKN-M)-Yl(N) ) / CXI ( N-f 1) -XI (N) )
5110 GO TO 5180
5120 B=CY1C2)-Y1C1) )/(Xl(2)-Xl(l»
5130 D0=( <Y1(3)-Y1C2) )/(Xl(3)~Xl(2) ) ~D) / ( XI (3) - 1) )
XI (
5140 C < 2 í 1>=B+BO#(X1(1)-X1(2))
5150 D=<Y1(N)-YKN-1) )/<Xl(N)"Xl(N~D)
5160 D0=í (Y1(N+1)-Y1CN> )/(XKN-M)~Xl(N) ) -D )/CXl N-f 1 )- N-l
5170 c ( 2 .- N+i ) =n+rio# < 2*xi c N-M > »xi < N- D -xi ( N»
5180 REM CALCULAR DERIVADAS Y ALMACENAR EN C(2,I>
5190 DELETE D,DO
5200 DIM BíN+1) fDO(N+l)
5210 REM D AUXILIAR PARA ALMACENAR DELTA (Xi)
5220 REM DO AUXILIAR
5230 HO<1)=1
5240 B(l)=0
5250 FOR J=2 TO N-M
5260 ri(J)=Xl(J)-Xl<J-l)
 ANEXO 1 0PR/GRAFICOS PAG, 02
5270 n O < J ) = < C < l , J)~C(1 7 J-1> )/D(J)
5280 NEXT J
5290 FOR J=2 TO N
5300 C < 2 r J ) =3* ( D C J ) #DO ( J + i ) -f-D < J + l ) #DO < J ) )
5310 DO(J)=2*(D(J)-KKJ + 1) )
5320 NEXT J
5330 REM #*# SIGUE REDUCCIÓN DE MATRIZ DE COEFICIENTES
5340 FOR J=2 TO N
5350 T=-D<J+1)/DOCJ-1)
5300 BOCJ)=DO<J>+T*D<J-1)
5370 C<2, J)=C(2f J)+T*C<2. J-l)
5330 NEXT J
5390 REM XXX SIGUE EVALUACIÓN REGRESIVA DE INCÓGNITAS
5400 FOR J=N TO 2 STEP -1
5410 C(2? J)=(C(2rJ)-D(J)#C(2r J+l) )/ÜO<J) •
5420 NEXT J
5430 ÜELETE D ? DO
5440 REM ### CALCULO DE COEFICIENTES C(3?I)r CÍ4,I)
5450 FOR 1=1 TO N
5400 D = X1ÍI + 1>-XKI)
5470 D1 = (C<1 .I + D-Cíl.*!) )/D
5480 D3=CC27l)+C(2jI+l>-2*Dl
5490 C(3?I)=(Dl-CC2rI)-D3)/D
5500 C(4íI)=D3/D/D
5510 NEXT I
5520 REM *** F'OL* DE INTERPOLACIÓN LISTOS (COEF EN C<4,N-hl>
5530 RETURN
5540 REM *#* SUB DE EVALUACIÓN PARA X ALEATORIO ENTRE XI (1) Y X1CN2)
5550 10=1
5500 JO=N
5570 IF JO-IOM THEN 5040
5580 1=10-1
5590 1= 1 + 1
5600 IF X=>X1CI) AND X <=X1 CI+1) THEN
5610 IF KJO THEN 5590
5620 PRINT 'X ESTA FUERA DE RANGOG'
5630 STOP
5640 KO=INT< (JO+IO)/2)
5650 IF X<XKKO) THEN 5680
5660 IO = KO
5670 GO TO 5570
5680 JO = KO
5690 GO TO 5570
5700 D=X-X1(I>
5710 P = C(l7l)+imC(2<I >+D# (C(3f I)+D*
720 RETURN
 ANEXO 1 @F'R/ESCRIBE PAG» 63
1000 REM ** 0PR/ESCRIBE
1010 02=43
1020 GO TO 800
1030 REM IMPRESIÓN
IM DE VOLTAJES DE NODO Y CORRIENTES DE RAMA
1040 Hé~0
1050 PRINT •LIMPRESIQN DE RESULTADOS 1
1000 P9=32
1070 PRINT "JDESEA LA IMPRESIÓN EN PAPEL? (SI O NO)I G"?
1080 INPUT X$
1090 IF NOT CX$='S" OR X$=°SI"> THEN 1130
1100 P9=51
ano PRINT "JALISTE EL IMPRESOR' (RETURN)G'
1120 INPUT :X*
1130 IF P9 = 32 AND Hó=l THEN 1870
1140 IF F=0 AND F0=0 AND Fl = 0 THEN 1570
1150 REM **:* IMPRESIÓN VALORES COMPLEJOS
1160 PRINT @P9Í USING ll/OtZ*
1170 IMAGE .P/'NOMBRE DEL PROBLEMA: 'FA
1180 PRINT i@P9: m JVOLTAJES DE NODO Cf = Flí
1190 PRINT (?P9: USING 1200;
1200 IMAGE 72('-'>
1210 PRINT '(?P9-; USING 1220 í
1220 IMAGE 10X'NODO"1QXMMAGNITUD/VOLTIOS' 8X'ÁNGULO/GRADOS
1230 PRINT i@F'9Í USING 1200;
1240 PRINT i(?P9i
1250 SET DEi :ES
1200 FOR J= 1 TO N
1270 J2=SQR
1280 IF C2(.J)<>0 THEN 1330
1290 J3=0
1300 IF C1C J»0 THEN 1340
1310 J3=180
1320 GO TO 1340
1330 J3=SGN (C2CJ))*ACS(C1CJ)/J2)
1340 PRINT @P9Í USING 1350íJ ? J2fJ3
1350 IMAGE 15D*8D?15D*3D
1300 NEXT J
1370 IF HS = 1 THEN 1830
1380 PRINT @P9:"JJCORRIENTES DE R A M A '
1390 PRINT SP9: USING 1400:
1400 IMAGE 72<l-')
1410 PRINT @P9 USING 1420t
1420 IMAGE 10X RAMA"lOX'MAGNITUD/AMPERIOS'8X 1 ÁNGULO/GRADOS'
1430 PRINT (?P9 USING 14001
.1.440 PRINT :@P9Í
1450 FOR J= 1 TO M
1460 J2-SQR (C3(J)"2+C4(J)^2)
1470 IF C4(.J)<>0 THEN 1520
1480 J3~0
1490 IF C3< J»0 THEN 1530
1500 J3=180
1510 GO TO 1530
1520 J3 = SGN C C4 C J ) ) *ACS CCCJ ) / J2 )
1530 • PRINT 1@P9t USING 15 13
 ANEXO 1 GPR/ESCRIBE PAG, 04

1540 IMAGE 13D,15D,8D ,1511,30

1550 NEXT J
1560 GO TO 1830
1570 REM #*# IMPRESIÓN SOLO UALORES REALES (CIRCUITO REISITIUO A D.C.)
1580 PRINT @P9: USING 159012$
1590 IMAGE P/'NOMBRE DEL PROBLEMA: 'FA
1600 PRINT @P9:'JVOLTAJES DE NODO (f = 'íFlf' Hz)1
1610 PRINT [?P9t USING 1620Í
1620 IMAGE 72C-1 )
1.630 PRINT GP9Í USING 1640t
1640 IMAGE 10X" NODO'15X'MAGNITUD/VOLTIOS'
1650 PRINT @P9: USING 1620Í
1060 PRINT (?P9Í
1670 FOR j=i TO N
1630 PRINT @P9: USING 1690:J?C1(J)
1690 IMAGE [Link]
1700 NEXT J
1710 IF H8=l THEN 1830
1720 PRINT 6P9Í'JJCORRIENTES DE RAMA'
1730 PRINT (?P9: USING 1740:
1740 IMAGE 72('-')
1750 PRINT @F9: USING 1760:
:L760 IMAGE lOX'RAMA'lSX'MAGNITUD/AMPERIOS'
1770 PRINT (3P9: USING 1740t
.1.780 PRINT @P9Í
1790 FOR J=l TO M
1800 PRINT 0P9; USING 1810:JrC3<J)
1S10 IMAGE 13Dí20D*8n
1820 NEXT J
1830 REM
1340 IF P9O32 THEN 1870
1850
1860 GO TO 1070
END
 A N E X O 2

MANUAL DE USO DEL PROGRAMA : CIRCUITO

Descripción del circuito

Algunas consideraciones sobre la descripción del circuito

Modo de empleo, ejemplo

 A N E X O 2

DESCRIPCIÓN DEL CIRCUITO

Conocido el circuito y el valor de los componentes se debe:

1- Enumerar arbitrariamente los nudos del circuito empezando por

O (cero), éste está reservado para el nado de-.r-eferencia'
(^1, U 2 . . . U R ) .

2- Se enumera las ramas a partir del numero 1 en adelante (b l s

b2i ...bm).

3- Se asignan los sentidos para las corrientes en las ramas, és-

to determina cual es el 'nodo "desde" y el nodo "hasta".

4- Para cuando existen fuentes dependientes debe conocerse el nu_

mero de la rama "controlante" y "controlada", si son inducta_n_
cías mutuas los números de las ramas acopladas.

Los sentidos de referencia son importantes para la determinación

correcta del signo de lasfuentes independientes y/o dependientes de vol-
taje y corriente. (Ver Fig. 1.2 , Cap. I)
'92

ALGUNAS CONSIDERACIONES SOBRE LA DESCRIPCIÓN DEL CIRCUITO

El programa es capaz de tratar solo circuitos "completos", por

ejemplo circuitos en los que la diferencia de potencial -entre dos
nodos arbitrarlos está bien definida. En otras palabras no debe
existir el caso de que el potencial en una parte del circuito
"Flote" en relación al resto. Así el programa no es capaz de de-
sarrollar cálculos correctos en el siguiente circuito.

La razón es que los nudos 4,5,6 flotan con relación al resto del
circuito. Se soluciona cortoclrcultando los nudos O y 6.
2. No se pueden especificar fuentes de voltaje o corriente ideales
como constitutivos únicos de una rama, debe especificarse con su
resistencia y/o conductancia interna.

Fuentes ideales Fuentes

En un caso dado se puede poner Rg - 10 6Q o bien Gs = 1

3. Varios elementos pasivos en paralelo pueden ser considerados con-

formando una rama.
-2,

b,
(D 3 ramas
2 nudos ind
 Cuando se trata de inductancias mutuas y la auto-inductancia en
paralelo con otro elemento, deben considerarse como ramas di fe -
rentes la causa, es que, sólo la corriente que circula por la
inductancia es lo que provoca el flujo magnético que acopla a la
segunda inductancia

5 ramas
2 nudos

MODO DE USO DEL PROGRAMA

. Antes de explicar el modo de utilización del programa, es necesji_

rio que el usuario defina su circuito y organice sus datos para un ingre_
so ágil. A continuación se presenta un ejemplo.
Se tiene el circuito:

. Jí.

El circuito con indicación de ramas y nodos es:

b,

"3

bs
Se debe disponer de un modelo o circuito equivalente del transi_
tor, según los requerimientos de preci sion.

Enumerando las ramas y señalando su sentido en estos,se tiene:

b
vs

La organización de los datos cuando se trata de circuitos gran-

des es importante. Esto evita cometer errores u omisiones, que
si bien pueden ser corregidos, no dejan de ser incómodos.

La información a ingresarse en el computador para el ejemplo

dado es:
 F A C T O R FASE(GRADOS)
CLASE R A M A # RAMA CONTROLAN RAMA CONTROLADA V A L O R
0 PREFIJO
i TE 0 NODO "DES 0 NODO "HASTA"
DE.
Fuente Volt. 1 1 MILI 0
Resistencia 1 0 1 1 MICRO -
Capacidad 2 1 2 47.00 MICRO -
desistencia 3 2 0 27,00 MILI -
Desistencia 4 3 4 56.00 - -
Desistencia 5 4 0 940.00 - -
Capacidad 5 ' 4 0 47.00 MICRO -
Capacidad 6 6 0 47.00 MICRO •
-
Desistencia ' 7 0' ' 7 2.20 KILO -
Desistencia 8 0 6 56.00 KILO —
Capacidad 9 8 540.00 MICRO —
[_—— —1—
Resistencia 10 ' 8' 0 10.00 KILO -
Desistencia •11 6' 2 51.00 KILO —
Desistencia 12 2 3 2.00 KILO —
Beta (FCCC) - 12 13' 100.00 -
—
Conductancia 13 5. 3 0.25 MICRO -
Resistencia 14 6 5 2.525 KILO -
_
Beta (FCCC) 14 15 100.00 - ~~\
Conductancia 15 " T ' 5 0.25 MICRO
Con los datos debidamente preparados

Prenda el computador
Colocar el disco en una de las unidades
I n i c i a l i z a r el. reloj interno del computador mediante la i n s t r u c
ción:
CALL"SETTIM" 5 ."DD-MMM-AA / K H H : MM:SS"

luego presionar la tecla RETURN

DD : día
MMM: mes
AA : año

.tf : espacio en blanco

HH : hora
MM : minutos
SS .: segundos

Montar el disco en el sistema mediante la i n s t r u c c i ó n :

CALL"MOUNT",#ÍX$

# numero de la u n i d a d en la que se ha colocado el disco.

Luego p r e s i o n a r la tecla RETURN
I n d i c a r al sistema la u n i d a d donde está el disco
UNIT # (Esto es i n d i s p e n s a b l e si se coloca en la u n i d a d 1).

Cargar a la memoria del computador el programa maestro

C I R C U I T O , mediante l a i n s t r u c c i ó n :
OLD"CIRCUITQ"
Presionar RETURN
 - Comenzar la ejecución del programa con la Instrucción RUN y
luego presionar RETURN.
- Continuar adelante de una manera interactiva, siguiendo las i_n_
dicaciones que da el computador por medio de la pantalla.

En caso de emergencia se puede interrumpir la ejecución aplasta_n_

do dos veces la tecla BREAK, pero ésto debe ser evitado.
ESCUELA POLITÉCNICA NACIONAL
FACULTAD DE INGENIERÍA ELÉCTRICA C
v*
ANÁLISIS'DE CIRCUITOS ELÉCTRICOS • £
LJ—

NOMBRE DEL PROBLEMA! CASCUDO n

Q-
ARCH-IVOÍ CASC <~t-
NO, BE RAMAS: 15 o»
NO, DE NODOS INDEPENDA S ^
NO, DE ELEMENTOS; 19 S
NO, CLASE RAMA DESDE HASTA MALOR FACTOR FASE ?
—r

1 FUENTE VOLT. 1 0 0 1*000 -3 0,000. %•

?. RESISTENCIA 1 O 1 i , 000 ~6 J,
3 CAPACIDAD 2 1 2 47 * 000 -ó •
4 RESISTENCIA 3 2 0 27,000 ¿
5 RESISTENCIA 4 ' 3 4 56,000 O
ó RESISTENCIA . 5 4 O 940,000 O
7 CAPACIDAD r 5 4 O 47*000 -ó
8 CAPACIDAD ó ó O 47,000 -ó
~ ? ' RESISTENCIA 7 0 7 2,200 3
10 RESISTENCIA S O • 6 56,000 3
1 1 CAPACIDAD 9 7 8 540,000 -6
12 RESISTENCIA 10 8 O 10,000 3
13 RESISTENCIA 11 . ó 2 51,000 3
14 RESISTENCIA 12 2 3 2,000 3
15 BETA (FCCC) O 12 13 100,000 O
16 CONDUCTANCIA 13 5 3 0,250 -6
•17 RESISTENCIA 14 ó 5 2,525 3
18 BETA (FCCC) O 14 15 100,000 O
19 CONDUCTANCIA 15 7 5 • 0,250 -ó
NOMBRE BEL PROBLEMA i CASCOD.O

VOLTAJES DE NODO <f ~ 1000 Hz)

NODO MAGNlTUn/UOLTl'GS ÁNGULO/GRADOS

1 0» 00100000 0 f 000 eí'

o 0,00099998 0*036 O-
-V—
_, o
0, 00073940 ~0 * 8 6 6 Crt
4 0* 00004462 -87,200
5 -177,472 o
0,00032630 O"
í> 0, 00000037 93,058 íí
3
7 0. 02329292 -177,407
o 0-
0,02329292 -177*405 0

ORRIENTES HE RAMA

RAMA MAGN1 TUD/AMPERIGS ÁNGULO/GRADOS

1 0,00000019 1*817
? 0 * OÜOOOOÍ9 1,817
3 0* 00000004 0*036
4 0, 00001318 2.594
5 0» 00001318- 2«594
o/ 0. OOOOOOIl -176*942
7 0, 00001059 2*593
S 0, 00000000 -80*942
9 0,00000233 -177,405
10 0}00000233 -177.405
11 0» 00000002 -179.> 985
12 Q, 00000013 2*594
13 0. 00001305 2*594
14 0, 00000013 2*594
15 ' 0» 00001292 2,594
U

I I B L I 0 6 R A F I A

1. STANDHAMMER 0., "Circuit Analysis by Digital Computer", Prentice

Hall, 1975.

2. CALAHAN D., Macnee A . B . , Me Mahon L., "Análisis moderno de circuitos

Interamericana, 1977.

3. CONTÉ S . D . y de BOOR C . , "Análisis numérico", Me Graw Hill Book C o . 3

1972.

4. "Introduction to BASIC programming", Tektronix.

5. "File Mananger 4907, Operator's Manual", Tektronix.