Sistemas Dinámicos
Sistemas Dinámicos
Sistemas Dinámicos
Circuito Chua.
1. Punto 1. Encontrar las soluciones de equilibrio, en las ecuaciones obtenidas del circuito Chua.
dI(3) −V2
=
dt L
dV(2) I(3) (V1 − V2 )
= +G
dt C2 C2
dV(1) (V2 − V1 ) f(VNR )
=G −
dt C1 C1
Escribiendolo de manera mas ordenada tenemos:
G 1
V˙1 = (V2 − V1 ) − f(V )
C1 C1 1
G 1
V˙2 = (V1 − V2 ) + I3
C2 C2
1
I˙3 = − V2
L
1
1 1 1 1 R5 Vsat
Ga = − − , Gb = − + ,B= .
R3 R6 R3 R4 R5 + R6
Con lo mencionado de nuestras ecuaciones presentaremos la forma normalizada o adimensional
de los voltajes V1 , V2 y la corriente I3 en el circuito de Chua, encontrando que:
V1 V2 R
x= ,y = , z = I, a = Ga R,
B B B
C2 R2 t
α= ,β = C2 , τ = , b = Gb R
C1 L RC2
Entonces, se tiene un sistema de tres ecuaciones diferenciales que en principio depende de
dos parámetros α y β. Pero, consideraremos la resistencia de la bobina como un valor muy
pequeño y despreciable, en lo cual no afectara nuestras ecuaciones.
dx
= α(y − x − f (x))
dt
,
dy
=x−y+z
dt
,
dz
= −βy
dt
La normalización de la función no lineal del diodo de Chua es f (x) = bx+0.5(a–b)(|x+1|–|x–1|).
Ahora consideraremos nuestro caso inicial en el instante donde no existe variación (0, 0, 0),
con lo mencionado tambien debemos establecer el comportamiento de nuestra bobina (induc-
tancia), condensador (capacitor) y resistencia, cuando no se presenta ningún tipo de carga al
circuito.
α(y − x − f (x)) = 0
,
x−y+z =0
,
−βy = 0
2
Ahora considerando nuestro parametros con el comportamiento del circuito llegamos a la con-
clusión de:
α = 1, no afectara el comportamiento del circuito. β = 1, no altera el circuito, f (x) en el t = 0
lo podemos tomar como 0 debido que no esta incidiendo en el sistema a causa de no haber
voltaje.
(y − x) = 0
,
x−y+z =0
,
−y = 0
3
como Ga y Gb son el valor de las pendientes probables de nuestro sistema, podemos considerar
los valores extremos sin tomarlos, encontrando dos puntos de equilibrio más en las cercanias
de (−1) y de (1).
como estos se aplican en f (x) tenemos que los puntos de equilibrio en nuestro sistema estan
en : (0, 0, 0), (1, 0, 0) y (−1, 0, 0)
Para aplicar la matriz jacobiana a nuestro sistema es necesario aclarar que esta es una matriz for-
mada por las derivadas parciales de primer orden de una función. Para nuestro caso como tenesmos
tres ecuaciones nos queda:
df df df
dx dy dz
dg dg dg
J =
dx dy dz
dm dm dm
dx dy dz
En la ecuación 1:
α(y − x − (bx + 0.5(a–b)(|x + 1|–|x–1|).))
Ẋ = α (y − x − bx − 0.5a |x + 1| + 0.5a |x − 1| + 0.5b |x + 1| − 0.5b |x − 1|)
Ẋ = α(−1 − b + −0.5a(x+1)+0.5b(x+1)
|1+x| + 0.5a(x−1)−0.5b(x−1)
|−1+x| )
Ẏ = α
Ż = 0
En la ecuación 2:
x−y+z
Ẋ = 1
Ẏ = −1
Ż = 1
En la ecuación 3:
−βy
Ẋ = 0
Ẏ = −β
Ż = 0
4
−0.5a(x+1)+0.5b(x+1) 0.5a(x−1)−0.5b(x−1)
α(−1 − b + |1+x| + |−1+x| ) α 0
1 −1 1
0 −β 0
Con lo anterior se remplaza cada punto de equilibrio en la matriz Jacobiana y se buscan los valores
y vectores propios
PUNTO (0,0,0)
−α(1 − 1.22) α 0
J(0,0,0) = 1 −1 1
0 −β 0
Valores Propios.
(A − λI)
α(0.22) − λ α 0
1 −1 − λ 1
0 −β −λ
Para simplificar nuestros valores, consisderaremos los valores de nuestro circuito, de acuerdo a
los más convenientes para que el sistema sea caotico en ello se toma C1 = 10nf, C2 = 100nf y
R = 2KΩ, L = 18mH, pérmitiendo determinar los valores de nuestros parametros α y beta.
por lo tanto los valores propios son (se omitiran los pasos del respectivo calculo):
5
Z(t) = (0.3013−0.1431i)k1 e3.086t +(−0.0394−0.1615i)k2 e(−0.943+3.867i)t +(−0.92)k3 e(−0.943+3.867i)t
PUNTO (1,0,0)
Para simplificar nuestros valores, consisderaremos los valores de nuestro circuito, de acuerdo a
los más convenientes para que el sistema sea caotico en ello se toma C1 = 10nf, C2 = 100nf y
R = 2KΩ, L = 18mH, pérmitiendo determinar los valores de nuestros parametros α y beta.
C2 R2 (2000)2
α= = 10, β = C2 = (100nf ) = 200/9
C1 L 0.018
Antes de evaluar nuestro sistema para x = 1 debemos considerar el limite del mismo para ver que
sucede en este punto, dado que nuestra f (x) presenta inconveniente.
0.5a(x − 1) − 0.5b(x − 1)
Limx→1
|x − 1|
. Al resolver el limite encontramos que el limite diverge, y se encuentra una cercania por derecha e
izquierda encontrando que:
0.5(a − b) o −0.5(a − b)
al resolverlo por izquierda encontramos el mismo sistema que en el punto (0, 0, 0)
−α(1 − 1.22) α 0
J(1,0,0) = 1 −1 1
0 −β 0
2, 2 10 0
J(1,0,0)
1
= −1 1
200
0 − 0
9
Pero al resolverlo por derecha encontramos que:
.
−α(1 − 0.728) α 0
J(1,0,0) = 1 −1 1
0 −β 0
Valores Propios.
(A − λI)
−α(0.272) − λ α 0
1 −1 − λ 1
0 −β −λ
Aplicamos el método de Montante o Regla de Sarrus o Regla de triángulo:
6
Reemplazazmos los valores de α y β y simplificamos:
−1
(225λ3 + 837λ2 + 3362λ + 13600)
225
por lo tanto los valores propios son, (se omitiran los pasos del respectivo calculo):
PUNTO (-1,0,0)
Como en los casos anteriores, para simplificar nuestros valores, consisderaremos los valores de
nuestro circuito, de acuerdo a los más convenientes:
α = 10, β = 200/9
Además debemos evaluar nuestro sistema para x = −1 considerando el limite del mismo para ver
que sucede en este punto, dado que nuestra f (x) presenta inconveniente.
−0.5a(x + 1) + 0.5b(x + 1)
Limx→−1
|x + 1|
. Al resolver el limite encontramos que el limite diverge, y se encuentra una cercania por derecha e
izquierda encontrando que:
0.5(a − b) o −0.5(a − b)
al resolverlo por izquierda y derecha encontramos el mismo sistema que en el punto (1, 0, 0), en lo
cual se determinaran los mismos valores y vectoers propios.
Con los valores y vectores propios encontrados, nos define un plano en el espacio fase, donde la parte
real la podriamos considerar como una trayectoria cercana a un punto de equilibrio y de acuerdo a
esto la trayectoria tiende a un espiral asintotico hacia el punto de equilibrio si consideramos el caso
< 0, a lo largo del vector correspondiente, luego por la trayectoria tendera a alejarse del punto de
equilibrio. En el caso que consideremos > 0 tendra una tendencia de alejarse en forma de espiral
asintotico al infinito y luego ingresa hacia el otro punto de equilibrio.
3. Punto 3 Consultar sobre punto fijo Hiperbólico y considerarlos en el modelo.
7
La hiperbolicidad representa un papel central ya que muestra el paradigma de los sistemas que se
consideran “caóticos” (sistemas de comportamiento impredecible) a pesar de lo cual se tiene una
descripción de su dinámica. Además presentan propiedades de estabilidad, lo que da a entender
que esta “caoticidad” no se destruye por perturbaciones ligeras que ocurren en el sistema.
Para poder considerar los puntos fijos hiperbolicos del sistema es necesario considerar si el espacio de
trabajo es un difeomorfismo, donde este se define como: un isomorfismo (tener la misma estructura)
en la categorı́a de las variedades diferenciables (es decir, un difeomorfismo es un homeomorfismo
(función biyectiva continua y cuya inversa es continua.) diferenciable con inversa diferenciable).
Como tal un difeomorfismo es una aplicación que posee aplicación inversa, por supuesto estas dos
aplicaciones son diferenciables.
dicho de otro modo un sistema plano ẋ = f (x) es hiperbólico si la matriz jacobiana en ese punto,
Df (ẋ) es hiperbólica (tiene valores propios todos con parte real no nula).
En puntos que no tengamos parte real los consideraremos no hiperbólicos y esto ocasionar la im-
posibilidad de conocer el comportamiento de las soluciones en vecindad del mismo a partir de la
linealización. El caso es contrario cuando se trata de puntos hiperbólicos y esto lo demuestra el
Teorema de Hartman-Grobman.
De acuerdo a lo mencionado anteriormente y a los resultados del Punto 2 podemos determinar que nuestro
sistema Chua tiene puntos hiperbolicos dados que todos sus valores propios son con parte real 6= 0
y analizando el comportamiento de ellos son espirales que se entrelazan en si.
Referencias
1. MATEMÁMTICAS MM (2010). DISCRIMINANTE O HESSIANO. Recuperado de:
https://blog.utp.edu.co/alejandropinto/files/2010/11/Discriminante-o-Hessiano.pdf
2. Flores, S. (2008). Modelo comportamental de una clase de circuitos caóticos mediante el uso de
SC-CNNs.
https://inaoe.repositorioinstitucional.mx/jspui/bitstream/1009/426/1/FloresGS.pdf
8
4. Nariño, J. (2002). Estudio de los efectos producidos al acoplar dos circuitos de chua y aplicaciones
del circuito chua a las telecomunicaciones. recuperado de:
https://repositorio.uniandes.edu.co/bitstream/handle/1992/15344/u234753.pdf?sequence=1
5. Zapata, A. Vidal, F. (). Diversas rutas alcaos en el circuito chua. Recuperado de:
file:///C:/Users/user/Downloads/Dialnet-DiversasRutasAlCaosEnElCircuitoDeChua-2986555.pdf
6. Jimenez, A. (1995). sobre los puntos fijos hiperbolicos de un difeomorfismo local del plano. Recu-
perado de:
https://dialnet.unirioja.es/descarga/articulo/7146005.pdf
7. (). Difeomorfismo.
https://es.wikipedia.org/wiki/Difeomorfismo: :text=En