0% encontró este documento útil (0 votos)
50 vistas9 páginas

Sistemas Dinámicos

Descargar como pdf o txt
Descargar como pdf o txt
Descargar como pdf o txt
Está en la página 1/ 9

Assignment No 2 - SISTEMAS DINÁMICOS

Omar Granados (614181020)


Marzo 2021

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

f(V(1) ) = Gb V1 + 0.5(Ga + Gb )[|V1 + E| − [V1 − E]]


Buscamos el normalizar nuestras ecuaciones y para ello comenzaremos considerando la función
f(V(1) ) que tiene pendientes negativas en Ga , Gb y un punto de ruptura que depende del voltaje
de saturación del amplificador operacional B o número 2(detalles del fabricante). Estos elementos
están en función de los valores de componentes del circuito del modo siguiente:

(a) Diodo de Chua a (Circuito simulado)

(b) Diodo de Chua (Circuito integrado)


Figura a, b es el arreglo experimental para medir la curva caracteristica del diodo chua. Apli-
cando un Voltaje Vs de una función sinusoidal al circuito en serie compuesto por la resistencia.

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|).

BUSCANDO LOS PUNTOS DE EQUILIBRIO

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

(c) Circuito en estado de t = 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

Los puntos encontrados son −y = 0, x = −z y x = −y, aplicando transitividad concluimos


que x, y, z = 0, el punto de equilibrio en el instante 0 = 0, debido a que no tenemos influencia
de ningun elemento a parte de la resistencia del sistema, pero nuestro sistema partira desde el
origen enfocado al comportamiento que se desarrolle, por ello es necesario analizar los dos (2)
estados restantes de nuestro amplificador operacional (> 0, < 0), donde es dado por la función
f(x) (Los puntos de corte que se generan):

(d) Puntos de saturación

Gb x + (Gb − Ga )1 si (x < −1)


f (x) = Ga x si (−1 < x < 1)
Gb x + (Gb − Ga )1 si (x > 1)

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)

2. Punto 2 Calcular la matriz Jacobiana y evaluar los puntos de equilibrio

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

tomamos el jacobiano con las respectivas ecuaciones de nuestro sistema:


Ẋ = α(y − x − f (x)),
Ẏ = x − y + z,
Ż = −βy
NOTA. nuestra función f( x) tiene unos valores como son a, b estos se pueden considerar como:
a = –1, 22, b = –0, 728

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

Ahora escribimos nuestros resultados en la matriz jacobiana.

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.

Con esto procedemos a encontrar nuestros λ:


6 451 440 −1
−λ3 + λ2 − λ+ = (45λ3 − 54λ2 + 451λ − 2200)
5 45 9 45

por lo tanto los valores propios son (se omitiran los pasos del respectivo calculo):

λ1 = 3.08597391, λ2 = −0.94298695 + 3.8669192i, λ3 = −0.94298695 − 3.8669192i

Los vectores propios son:


V1 = (−0.84069685, 0.30130425 − 0.14306425i, 0.30130425 + 0.14306425i)
V2 = (−0.0744835, −0.03937774 + 0.16147683i, −0.03937774 − 0.16147683i)
V3 = (0.53635902, −0.92796713, −0.92796713)

con estos resultados


 planteamos la solucion
 general:  
−0.841 −0.0745
W (t) = k1 e3.086t  0.3013 + 0.1431i  + k2 e(−0.943+3.867i)t  −0.0394 + 0.1615i 
0.3013 − 0.1431i −0.0394 − 0.1615i
 
0.5364
+k3 e(−0.943+3.867i)t  −0.928 
−0.928

X(t) = (−0.841)k1 e3.086t + (−0.0745)k2 e(−0.943+3.867i)t + (0.5364)k3 e(−0.943+3.867i)t

Y (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

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:

34α + 125 2 91α − 125β 34αβ −1


−λ3 − λ + λ− = (125λ3 +(34α+125)λ2 −(91α−125β)λ+34αβ) = 0.
125 125 125 125

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):

λ1 = −3.88191717, λ2 = 0.08095859 + 3.94515107i, λ3 = 0.08095859 − 3.94515107i

Los vectores propios son:


V1 = (−0.82876512, 0.28108462 + 0.19101197i, 0.28108462 − 0.19101197i)
V2 = (0.09629564, 0.00337353 + 0.16439379i, 0.00337353 − 0.16439379i)
V3 = (0.55124906, −0.92599631, −0.92599631)

con estos resultados


 planteamos la solucion
 general:  
−0.829 0.096
W (t) = k1 e−3.88t  0.281 + 0.1910i  + k2 e(0.08096+3.945i)t  0.0034 + 0.1644i 
0.281 − 0.1910i
  0.0034 − 0.1644i
0.5512
+k3 e(0.08096−3.945i)t  −0.926 
−0.926

X(t) = (−0.829)k1 e−3.88t + (0.096)k2 e(0.08096+3.945i)t + (0.5512)k3 e(0.08096−3.945i)t

Y (t) = (0.281+0.1910i)k1 e−3.88t +(0.0034+0.1644i)k2 e(0.08096+3.945i)t +(−0.926)k3 e(0.08096−3.945i)t

Z(t) = (0.281−0.1910i)k1 e−3.88t +(0.0034−0.1644i)k2 e(0.08096+3.945i)t +(−0.926)k3 e(0.08096−3.945i)t

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.

Definición: Sea f : M →− M difeomorfismo y p un punto fijo de f . Decimos que p es hiperbólico si


Df p : T pM → − T pM es hiperbólico (no tiene valores propios de módulo uno). Un punto periódico
de perı́odo k se dice hiperbólico si es un punto fijo hiperbólico de f k

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.

TEOREMA (Hartman-Grobman): Sea f de clase C 1 y x− un punto crı́tico hiperbólico del sis-


tema plano ·x = f (x). Entonces hay una vecindad de x− en la cual ẋ = f (x) es topológicamente
equivalente a su linealización ẋ = Df (x− )x.
En otras palabras el teorema (Hartman-Grobman) garantiza que: para puntos hiperbólicos el com-
portamiento del sistema no lineal cerca del punto es “idéntico” al del sistema lineal correspondiente.

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

3. Sambarino, M. (2005). Tópicos de Sistemas Dinámicos. Recuperado de:


http://www.im.ufrj.br/ arbieto/ensino/2012/samba.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

También podría gustarte