Lévano Huamaccto If 2017
Lévano Huamaccto If 2017
Lévano Huamaccto If 2017
NATURALES Y MATEMATICA
SQ B 1 2013 f
{IE 6 030; 1 030)
:2 E 034»\§ )2
'l\> 030D:A\
55:5 (03
0 //$50
EM; 2018}
INFORME FINAL DEL PROYECTO DE INVESTIGACION
034CAOS
EN LA DINAMICA DE UN OSCILADOR NO LINEAL
(Periodo de Ejec}401ciénz
De 01 de enero del 2017 al 12 de Diciembre del 2017)
IJNDICE 1
1.1.|ndice de tablas 3
1.2.Indice de }401guras 4
1.3.1ndice de gra}401cas 5
lI.RESUMEN Y ABSTRACT
2.1. Resumen 6
III. INTRODUCCION
a) Seccién de Poincaré 15
b) Espectros de Potencia 15
c) Exponentes de Lyapunov 16
a) Caos 28
c)Linealidad y no linealidad 29
V. MATERIALES Y METODOS
5.1. Materiales 30
5.1.1. De ejecucién 30
5.1.2. De impresién 30
030n\
5.3. Las Técnicas de anélisis o Ios medios estadisticos 31
5.5. Metodologiaa 31
VI. RESULTADOS 34
VII. DISCUSION 41
VIII. REFERENCIALES 42
IX. APENDICES 43
X. ANEXOS 51
Matriz de Consistencia
2
TABLA DE CONTENIDO
3
TABLA DE FIGURAS
4 .
TABLA DE GRAFICAS
Gra}401ca
6.2 Espacio de fase sin Caos frecuencia 0.05 Pag.36
Gra}401ca
6.4 Espacio de fase con Caos frecuencia 0.1 Pag.38
Gra}401ca
6.5 Espacio de fase oon Caos frecuencia 0.3 Pag.39
6.6
Gra}401ca Espacio de fase oon Caos frecuencia 0.6 Pag.40
5
RESUMEN
fisico Oscilador no lineal amortiguado y forzado con una fuerza externa de tipo
' través de un programa desarrollado por Wolf et al. Se obtuvo valores del
trabajo de investigacion.
W
' 6
ABSTRACT
physical system with an external force of type sinusoidal through the indicator
oscillator. }401rst
determine the phase space where it is observed if it presents
Positive values of the Lyapunov exponents were obtained for disordered dynamics
-% 7
Ill. INTRODUCCION
Con el avance de la Ciencia se sabe que son mas sistemas no lineaies que los
\
031 Es importante determinar cuéndo un sistema presenta comportamiento de Caos
El problema identi}401cado
es la detenninacién de| estado de Coas en Ios
propuesto por el ruso matemético del mismo nombre que propone que el
a) Problema General
siguiente,
b) Problemas Especificos:
N 9
c) Objetivo General y Alcances de la investigacion
-Objetivo General
-Objetivos Especificos
dice que el Caos en ubicuo. Los trabajos de modelamiento que realizamos, |os
_ 030 que proponemos para estudiar otros fenémenos de la naturaleza debe ser
X 10
031 controlando Ios parémetros propios del sistema para saber |as condiciones
030 030 11
IV. MARCO TEORICO
' Zhang e Harne (2016) propone el estudio del Oscilador con efecto
un sistema. 030
12
i
Usualmente en }401sica
aprendimos que las leyes son descritas por ecuaciones
2013).
de la Mecanica por ejemplo las Leyes de Newton; por ello se dice también
Saletan E.,1998).
}401nales
del siglo XIX que eran contran'o al determinismo(Ra}401ada
A.,1990).
fuera ei caso de que las leyes naturales ya no tuvieran ningiin secreto para
. nosotros, pero no es asi en los resultados numéricos que nos describen que
& 13
solo .podn'amos conocer Ia situacién inicial aproximadamente. Si nos
E.,1998).
La teoria de| caos, es una teoria que trata sobre el fenomeno del
es por ello que en esta teoria se han desarrollado indicadores de| Caos para
\ 14
determinar si el sistema muestra el caos , entre los indicadores tenemos :
a) Seccién de Poincaré
2006).
FIGURA N° 4.1
3&1 ~
-}401?:§
mo» 2
53$
. _ _ _ $9; r
:<-on
Fuente: Ra}401ada.
Dinémica CIésica.1989
b) Espectros de Potencias
una transformada de Fourier continua, para ello se tiene que hacer variar
A.,1990)
FIGURA N°4.2
PERIODICO Y EN CAOS
030 ; I. q -_ .; ;:',-;2;": 030_ 031r::-_;'- .
o . .. - .- ,. 3 o _v ». '
Fuente: Ra}401ada,
Dinémica 030CIésica,
1989
c) Exponentes de Lyapunov
\ b 16
La evolucién temporal se rige por un espectro continuo de frecuencias
caos. Estudiemos el caso de sistemas con una sola variable que asumen dos
A.,1990)
los dos valores iniciales, Ia diferencia d,. entre los dos valores xn después de
aproximadamente.
«:1 = 034X0
+ s) ~ r(x.,) 2 5(: 024:)KO (2)
\ 17
Xn+1 = f(Xn)=-f(f(Xn-1)): f(f(f(Xn-2)))= = f(f(f(-»-(f(Xo)))))- (3)
es. 4
I d = f(x,1 + 2)
034 f(x,,) = f(x0 + 2)
024 f}401xo)
024 = sen" (4)
obtenemos I
In{[f n (X0 + 8) _
f S034
031 030_O
=
In[e 21}.] =
nix
-
(5)
Donde,
«u M, {T}{-- 024m(,,
}-..{ 024m
-nio_ 024dF1(X;1_1) dFl{Xn_7) dfutxo)
& 18
Reempiazando la ecuacién 7 en la ecuacién 6 para calcular el exponente
dzx dz 3
024,+B 024 024x+8x
=FCos(mt) (9)
dt dt
19
FiGURA N° 4.3
(3) (b)
6.0 6.0
030 = * if u 031
U
030(Swan
031 4035.0 u)6u.o
031 ""94.a 0.0 4,0
I x
FIGURA N 4.4
034 '
ESPACIO DE FASE
is .- J ' I xi}
4 _ I '7 O :.%
2004
20
La Dinémica esté compuesta por dos campos de estudio, Ia Dinémica Lineal
signi}401cativamente
el comportamiento dinémico del sistema convirtiendo los
sin rotar), suspendido entre dos muelles iguales, tal como se muestra en la
figura 5. Si ambos muelles estén sin de deformar(o sea, sin tensién y, por
tanto; sin energia potencial los dos lados) cuando Ia panicula se encuentre
, 024x(s
l)
024
\ 21
_ FIGURA 4.5
" ~
FIGURA 4.6
______""_ .._'f.__\;,:M
K 030er
/ , 031
s
' Fuente: Marion J., Dinémica Clésica de las particulas y Sistemas, 1981
- 22
La fuerza que acttia sobre la particula de masa ,M, es ,F, igual a,
F ' 0242z:(.-:l)$'crz(9),
024
como : = wxz + I3
de forma que,
_ 9 x I
Jen = 024
= 024 024 024 024 024
5 \/x3 + :2
F= 1 024.'+
\;A 035
(v 030x2
030
e :2 L) = zzcc 1 1
F (10)
H031031 V1 +G} 031
tendremos(Marion J.,1981),
23
Y cuando despreciamos todos los términos, salvo el primero, seré
k
F 5-! (11))?
' _ vcd N 3 3
» F(X)§ 0242 x030
024[rc(l
024l- 024d,/I
]x
expresién,
F _ dU(x)
_ dx
_ 55 _ xii-d) 3 _ __.s'u-{ac}
A 2(1)x [ F J) 030
024 .2; ' (12)
_E E 030
U(x)-z(2(l'))x i[j""F 034"
2 +4 x.
4 031] (13)
24
Luego considerando que el sistema se mueve solamente en forma horizontal
d ''at.\ (al.) 1 2
030 " T 024 024
= 030
- + ~ 15
d: (.6471) aa, 0 Q3
031 ( )
Donde:
horizontal.
\ 25
030 ' Reemplazando el Lagragiano en la ecuacién dinémica y las fuerzas
3 1 1 -,2 030
_:(_; .2 _1li((§-t}402j
.4) =_ 030g
¢:(a2(zM? 034 any 2(2 J1 4 :3 1 b(a't)+
(16)
+2 .1 +
030 = 024b¢_ +034§)
F0 Cos(wt)
(17)
A la ecuacién,
3 ii
}402 030
E . L-#3 .a_5 .
1:1 + M (dc) + 2 (in) X + [ 131:: J1 " M COS(wt) (18)
_ 26
Se hacen |as siguientes representaciones,
- i . F
.1: + 23:: + (A)(j.22. +
0316x3 = B 024:_ 031Cas(_wt) (19)
a}401rmar
que la ecuacién dinémica corresponde a un Sistema débilmente no
27
4.4. Marco conceptual
a) Caos
}401able
de su comportamiento a pesar de de tener un carécter determinista, aun
de }402uidos,
rayos léser, crecimiento de poblaclones, transmislén de
diagramas de bifurcacién.
La teoria del Caos esté dividida en dos amplias éreas de estudios, sistemas
\ 28
' c) Linealidad y no Iineaiidad
J., 2013).
\ 29
V. MATERIALES Y METODOS
5.1. Materiales
5.1.1. De ejecucién
especi}401camente
a la Fisica computacional, por ello se usé ordenadores
90 para Runge Kutta de 410 orden (Press W., Teukolsky 3., Vettering
investigacién.
investigacién.
con el programa.
5412! De la impresibn
materiales:
o Elementos de escritorio
X 30
5.2. La poblacién y muestra
poblacién y muestra.
5.3 Las
030 técnicas, procedimientos e instrumentos de recoleccién de datos
Teérica.
5.5. La Metodologia
FIGURA 4.7
K .: \\\s
M V L X \ \«
030
__.____ ________., 030M
I. 030 9 030I;
31
024La
recopilacién de datos se hace a través de un programa mencionado
11 ax 024.
024f
+ 2;? (- 024)
+ W 030):
+ 6.1.63 = F(.'ns(mt) , (20)
dz 030 dt
tiempo, I-'Cas(wt),
dz: _ '
dL 3
iii _ _9 _ 2 __ 3 ,
at ..,B at w 3: 6x + I'Co.s(cot)
(21)
K 32 .
l
i
resulta,
2 2 .
dia: 23 das _ 5 _3 Po
- 024;
+ + It + 024zx
= :2-C0S(QT) (23)
03 it ma dr mg mg, M
, . . . 2:3 031 2
-En el analisis se mantuvo constante |as cantidades , :3
024 ,
OJ.) 0
m = }402,}402
hasta
030lma
m = 3,4 mu con los datos que se muestra en el cuadro
. . . _ Z}402 s Fa
sugulente, mantenlendo constante las cantidades - 024,024 024 024g-, 030T
.
can an we M
\ 33
VI. RESULTADOS
TABLA N 6.1
034 .
1 9031 9 = I;
034O2 1:932!" - W0
034E111
Ejij
E
u
j
j
j
j
ET
\ 34
GRAFICA N 6.1
035
0.215798E+O1
o
\ -1.. 030 I, f
\ 034\1 030:,-
-10 V ,
in
now
030 2 4 5 s in 12 14 16 13
\ 35
GRAFICA N° 62
20
M an 024 024
5 030f 03
1: I r 030. )
.» 1 « 2 c, 030
030S % . . \ /I
\ \ V. . 030
.10 . V/4,4, *" 0307' M
\ 030/'
\ 36
GRAFICA N° 6.3
20 /\
AA\\.~v~'v*~~A
f W
° 1: I)))))! 031) ) 030l)
Va 1* 3"! 031 034031
034'
. ' A \
-212 H
_ _-1s__A_ -1o -5 0 5 1a 15 20
X 37
n I// ., h 034N§Cn:gRT|
5 0301 030T;33:::9P»:~§
1 034 .
Q sra}401h}401q}401r
035 5 %
agbtenidoks 030
2° }402
. "5 um»
Ill 030-I 031
\\V"77>
I K
1 030
\,/
V 030
4% 030:-I
"
X 39
g / %%
05$ 035?
g9F 031:Fof 030 /1
\;
4] 030? 20
A9
%
VII. DISCUSION
investigacion.
Amortiguado y Forzado.
\ 41
r
SUBROUTINE derivs(x,y.dydx)
REAL x,y(").dydx('),A,C,E,D
A=0.01
c=o.o2 030 .
D=3O.0
E=0.08
dydx(1)=y(2)
! dydx(2)=-((w1"2) 030y(1))
dydX(2)=-A"Y(2)-Y(1)-C*((Y(1))**3)+D 030C0S(E*X)
I dydx(3)=y(4)
I dydx(4)= 024((w2**2)*v(3))
return
END
! rkdumb.for
SUBROUTINE rkdumb(vstart,nvar,x1,x2,nstep,derivs)
INTEGER nstep,nvar,NMAX,NSTPMX
PARAMETER (NMAX=50,NSTPMX=100000)
EXTERNAL derivs
INTEGER i,k
REAL h,x,dv(NMAX),v(NMAX)
do 11 i=1,nvar
v(i)=vstar1(i)
V(i.1)=V(i)
11 continue
xx(1)=x1
44
031 x=x1 030
h=(x2-x1)/nstep
do 13 k=1,nstep
call derivs(x,v,dv)
call rk4(v,dv,nvar,x,h.v.derivs)
x=x+h
xx(k+1)=x
do 12 i=1,nvar
y(i,k+1)=v(i)
12 continue
13 continue
return
END
! rk4.for
SUBROUTINE rk4(y,dydx,n,x,h,yout,derivs)
INTEGER n,NMAX
EXTERNAL derivs
PARAMETER (NMAX=50)
INTEGER i
REAL h6,hh,xh,dym(NMAX).dyt(NMAX),yt(NMAX)
hh=h"0.5 -
V h6=h/6.
xh=x+hh
do 11 i=1,n
Yt(i)=Y(i)+hh*dYdX(i) '
11 continue
cail derivs(xh,yt,dyt) .
do 12 i=1,n -
Vt(i)=y(i)+hh"dyt(i)
45
12 continue
call derivs(xh,yt,dym)
do 13 i=1,n
Yt(i)=y(i)+h*dYm(i)
dym(i)=dyt(i)+dym(i)
13 continue
call derivs(x+h,yt,dyt)
do 14 i=1,n
14 continue
return
END
1985)
program ode
nn=n"(n+1)=
nn1=nn+1
external fcn
dimension y(nn1),yprime(nn1),v(nn1),A(nn1),B(nn1),C(nn1)
dimension D(nn1),cum(n),znorm(n),gsc(n)
open(1,}401|e="Lout.txt")
read(*,*) nstep,stpsze,irate,io
write(1,*) nstep,stpsze,irate.io
write(1,*)
v(1)=10.0
v(2)=1.0
v(3)=1.0
K 46
tme=0.0
do 10 i=n+1,nn
v(i)=0.0
10 continue
do 20 i=1,n
v((n+1)*i)=1.0
cum(i)=0.0
20 continue
do 100 m=1,nstep
if ((m/100)*100.EQ.m) wn'te(*,*)m
do 25 j=1,irate
do 26 i=1,nn
030/(i)=V(i)
26 continue
t=tme
do 27 i=1,nn
A(i)=yprime(i)
27 continue
do 28 i=1,nn
y(i)=v(i)+(stpsze*A(i))/2.0
28 continue
t=tme+stpszel2.0
do 29 i=1,nn
B(i)=yprime(i)
29 continue
do 30 i=1,nn
y(i)=v(i)+(stpsze'B(i))/2.0
30 continue
t=tme+stpsze/2.0
X 47
call fcn(t, y, yprime)
do 31 i=1,nn
C(i)=yprime(i)
31 continue
031 do 32 i=1,nn
y(i)=V(i)+(StPSZe*C(i))
32 continue
t=tme+stpsze
call fcn(t, y, yprime) 030
do 33 i=1,nn
D(i)=yprime(i)
33 continue
do 34 i=1,nn
v(i)=v(i)+stpsze*(A(i)+D(i)+2.0*(B(i)+C(i)))/6.0
34 continue
tme=tme+stpsze
25 continue
znorm(1)=0.0
do 38 j=1,n
znorm(1)=znorm(1)+v(n*j+1) 030*2
38 continue
znorm(1)=sqrt(znorm(1))
do 40 j=1,n
v(n*j+1)=v(n*j+1)/znorm(1)
40 continue
do 80 j=2,n '
do 50 k=1 ,1-1
gsc(k)=0.0
do so |=1,n _
gsc(k)=gsc(k)+v(n*I+j)*v(n*|+k)
50 continue
\ 48
do 60 k=1,n
do 60 l=1,j-1
v(n*k+j)=v(n*k+j)-gsc(|)"v(n*k+|)
60 continue
znorm(j)=0.0
do 70 k=1,n
znorm(j)=znorm(j)+v(n*k+j)**2
70 continue
znorm(j)=sqrt(znorm(j))
do 80 k=1,n
v(n*k+j)=v(n*k+j)/znorm(j)
80 continue
do 90 k=1,n
cum(k)=cum(k)+d|og(znorm(k))/a|og(2.0)
90 continue
write(1,334) tme,(cum(k)/tme,k=1_n)
334 format(1x,f12.6,2x,e12.6,2x,e12.6,2x,e12.6)
_ 100 continue
end
*m*m**..m*....Mmm...m..........m..............*........*
subroutine fcn(t,y,yprime)
b=4.0
sg=16.0
r=45.92
VPrime(1)=S9*(V(2)-Y(1))
yprime(2)=-y(1)*y(3)+r*y(1)-y(2)
VPrime(3)=Y(1)*Y(2)-b*Y(3)
do 10 i=0,2
vprime(4+i)=S9"(V(7+i)-Y(4+i))
vprime(7+i)=(r-y(3))*v(4+i)-v(7+i)-v(1) 030v(10+i)
\ 49
' 030 vr>rime(10+i)=y(2)*y(4+i)+y(1)*y(7+i)-b*v(10+i)
10 continue
. return
_ . end
031 V i X _ 50
/ X. ANEXOS
Matnz de Consistencia
51
UNIVERSIDAD NACIONAL DEL CALLAO
NATURALES Y MATEMATICA
SQ B 1 2013 f
{IE 6 030; 1 030)
:2 E 034»\§ )2
'l\> 030D:A\
55:5 (03
0 //$50
EM; 2018}
INFORME FINAL DEL PROYECTO DE INVESTIGACION
034CAOS
EN LA DINAMICA DE UN OSCILADOR NO LINEAL
(Periodo de Ejec}401ciénz
De 01 de enero del 2017 al 12 de Diciembre del 2017)
IJNDICE 1
1.1.|ndice de tablas 3
1.2.Indice de }401guras 4
1.3.1ndice de gra}401cas 5
lI.RESUMEN Y ABSTRACT
2.1. Resumen 6
III. INTRODUCCION
a) Seccién de Poincaré 15
b) Espectros de Potencia 15
c) Exponentes de Lyapunov 16
a) Caos 28
c)Linealidad y no linealidad 29
V. MATERIALES Y METODOS
5.1. Materiales 30
5.1.1. De ejecucién 30
5.1.2. De impresién 30
030n\
5.3. Las Técnicas de anélisis o Ios medios estadisticos 31
5.5. Metodologiaa 31
VI. RESULTADOS 34
VII. DISCUSION 41
VIII. REFERENCIALES 42
IX. APENDICES 43
X. ANEXOS 51
Matriz de Consistencia
2
TABLA DE CONTENIDO
3
TABLA DE FIGURAS
4 .
TABLA DE GRAFICAS
Gra}401ca
6.2 Espacio de fase sin Caos frecuencia 0.05 Pag.36
Gra}401ca
6.4 Espacio de fase con Caos frecuencia 0.1 Pag.38
Gra}401ca
6.5 Espacio de fase oon Caos frecuencia 0.3 Pag.39
6.6
Gra}401ca Espacio de fase oon Caos frecuencia 0.6 Pag.40
5
RESUMEN
fisico Oscilador no lineal amortiguado y forzado con una fuerza externa de tipo
' través de un programa desarrollado por Wolf et al. Se obtuvo valores del
trabajo de investigacion.
W
' 6
ABSTRACT
physical system with an external force of type sinusoidal through the indicator
oscillator. }401rst
determine the phase space where it is observed if it presents
Positive values of the Lyapunov exponents were obtained for disordered dynamics
-% 7
Ill. INTRODUCCION
Con el avance de la Ciencia se sabe que son mas sistemas no lineaies que los
\
031 Es importante determinar cuéndo un sistema presenta comportamiento de Caos
El problema identi}401cado
es la detenninacién de| estado de Coas en Ios
propuesto por el ruso matemético del mismo nombre que propone que el
a) Problema General
siguiente,
b) Problemas Especificos:
N 9
c) Objetivo General y Alcances de la investigacion
-Objetivo General
-Objetivos Especificos
dice que el Caos en ubicuo. Los trabajos de modelamiento que realizamos, |os
_ 030 que proponemos para estudiar otros fenémenos de la naturaleza debe ser
X 10
031 controlando Ios parémetros propios del sistema para saber |as condiciones
030 030 11
IV. MARCO TEORICO
' Zhang e Harne (2016) propone el estudio del Oscilador con efecto
un sistema. 030
12
i
Usualmente en }401sica
aprendimos que las leyes son descritas por ecuaciones
2013).
de la Mecanica por ejemplo las Leyes de Newton; por ello se dice también
Saletan E.,1998).
}401nales
del siglo XIX que eran contran'o al determinismo(Ra}401ada
A.,1990).
fuera ei caso de que las leyes naturales ya no tuvieran ningiin secreto para
. nosotros, pero no es asi en los resultados numéricos que nos describen que
& 13
solo .podn'amos conocer Ia situacién inicial aproximadamente. Si nos
E.,1998).
La teoria de| caos, es una teoria que trata sobre el fenomeno del
es por ello que en esta teoria se han desarrollado indicadores de| Caos para
\ 14
determinar si el sistema muestra el caos , entre los indicadores tenemos :
a) Seccién de Poincaré
2006).
FIGURA N° 4.1
3&1 ~
-}401?:§
mo» 2
53$
. _ _ _ $9; r
:<-on
Fuente: Ra}401ada.
Dinémica CIésica.1989
b) Espectros de Potencias
una transformada de Fourier continua, para ello se tiene que hacer variar
A.,1990)
FIGURA N°4.2
PERIODICO Y EN CAOS
030 ; I. q -_ .; ;:',-;2;": 030_ 031r::-_;'- .
o . .. - .- ,. 3 o _v ». '
Fuente: Ra}401ada,
Dinémica 030CIésica,
1989
c) Exponentes de Lyapunov
\ b 16
La evolucién temporal se rige por un espectro continuo de frecuencias
caos. Estudiemos el caso de sistemas con una sola variable que asumen dos
A.,1990)
los dos valores iniciales, Ia diferencia d,. entre los dos valores xn después de
aproximadamente.
«:1 = 034X0
+ s) ~ r(x.,) 2 5(: 024:)KO (2)
\ 17
Xn+1 = f(Xn)=-f(f(Xn-1)): f(f(f(Xn-2)))= = f(f(f(-»-(f(Xo)))))- (3)
es. 4
I d = f(x,1 + 2)
034 f(x,,) = f(x0 + 2)
024 f}401xo)
024 = sen" (4)
obtenemos I
In{[f n (X0 + 8) _
f S034
031 030_O
=
In[e 21}.] =
nix
-
(5)
Donde,
«u M, {T}{-- 024m(,,
}-..{ 024m
-nio_ 024dF1(X;1_1) dFl{Xn_7) dfutxo)
& 18
Reempiazando la ecuacién 7 en la ecuacién 6 para calcular el exponente
dzx dz 3
024,+B 024 024x+8x
=FCos(mt) (9)
dt dt
19
FiGURA N° 4.3
(3) (b)
6.0 6.0
030 = * if u 031
U
030(Swan
031 4035.0 u)6u.o
031 ""94.a 0.0 4,0
I x
FIGURA N 4.4
034 '
ESPACIO DE FASE
is .- J ' I xi}
4 _ I '7 O :.%
2004
20
La Dinémica esté compuesta por dos campos de estudio, Ia Dinémica Lineal
signi}401cativamente
el comportamiento dinémico del sistema convirtiendo los
sin rotar), suspendido entre dos muelles iguales, tal como se muestra en la
figura 5. Si ambos muelles estén sin de deformar(o sea, sin tensién y, por
tanto; sin energia potencial los dos lados) cuando Ia panicula se encuentre
, 024x(s
l)
024
\ 21
_ FIGURA 4.5
" ~
FIGURA 4.6
______""_ .._'f.__\;,:M
K 030er
/ , 031
s
' Fuente: Marion J., Dinémica Clésica de las particulas y Sistemas, 1981
- 22
La fuerza que acttia sobre la particula de masa ,M, es ,F, igual a,
F ' 0242z:(.-:l)$'crz(9),
024
como : = wxz + I3
de forma que,
_ 9 x I
Jen = 024
= 024 024 024 024 024
5 \/x3 + :2
F= 1 024.'+
\;A 035
(v 030x2
030
e :2 L) = zzcc 1 1
F (10)
H031031 V1 +G} 031
tendremos(Marion J.,1981),
23
Y cuando despreciamos todos los términos, salvo el primero, seré
k
F 5-! (11))?
' _ vcd N 3 3
» F(X)§ 0242 x030
024[rc(l
024l- 024d,/I
]x
expresién,
F _ dU(x)
_ dx
_ 55 _ xii-d) 3 _ __.s'u-{ac}
A 2(1)x [ F J) 030
024 .2; ' (12)
_E E 030
U(x)-z(2(l'))x i[j""F 034"
2 +4 x.
4 031] (13)
24
Luego considerando que el sistema se mueve solamente en forma horizontal
d ''at.\ (al.) 1 2
030 " T 024 024
= 030
- + ~ 15
d: (.6471) aa, 0 Q3
031 ( )
Donde:
horizontal.
\ 25
030 ' Reemplazando el Lagragiano en la ecuacién dinémica y las fuerzas
3 1 1 -,2 030
_:(_; .2 _1li((§-t}402j
.4) =_ 030g
¢:(a2(zM? 034 any 2(2 J1 4 :3 1 b(a't)+
(16)
+2 .1 +
030 = 024b¢_ +034§)
F0 Cos(wt)
(17)
A la ecuacién,
3 ii
}402 030
E . L-#3 .a_5 .
1:1 + M (dc) + 2 (in) X + [ 131:: J1 " M COS(wt) (18)
_ 26
Se hacen |as siguientes representaciones,
- i . F
.1: + 23:: + (A)(j.22. +
0316x3 = B 024:_ 031Cas(_wt) (19)
a}401rmar
que la ecuacién dinémica corresponde a un Sistema débilmente no
27
4.4. Marco conceptual
a) Caos
}401able
de su comportamiento a pesar de de tener un carécter determinista, aun
de }402uidos,
rayos léser, crecimiento de poblaclones, transmislén de
diagramas de bifurcacién.
La teoria del Caos esté dividida en dos amplias éreas de estudios, sistemas
\ 28
' c) Linealidad y no Iineaiidad
J., 2013).
\ 29
V. MATERIALES Y METODOS
5.1. Materiales
5.1.1. De ejecucién
especi}401camente
a la Fisica computacional, por ello se usé ordenadores
90 para Runge Kutta de 410 orden (Press W., Teukolsky 3., Vettering
investigacién.
investigacién.
con el programa.
5412! De la impresibn
materiales:
o Elementos de escritorio
X 30
5.2. La poblacién y muestra
poblacién y muestra.
5.3 Las
030 técnicas, procedimientos e instrumentos de recoleccién de datos
Teérica.
5.5. La Metodologia
FIGURA 4.7
K .: \\\s
M V L X \ \«
030
__.____ ________., 030M
I. 030 9 030I;
31
024La
recopilacién de datos se hace a través de un programa mencionado
11 ax 024.
024f
+ 2;? (- 024)
+ W 030):
+ 6.1.63 = F(.'ns(mt) , (20)
dz 030 dt
tiempo, I-'Cas(wt),
dz: _ '
dL 3
iii _ _9 _ 2 __ 3 ,
at ..,B at w 3: 6x + I'Co.s(cot)
(21)
K 32 .
l
i
resulta,
2 2 .
dia: 23 das _ 5 _3 Po
- 024;
+ + It + 024zx
= :2-C0S(QT) (23)
03 it ma dr mg mg, M
, . . . 2:3 031 2
-En el analisis se mantuvo constante |as cantidades , :3
024 ,
OJ.) 0
m = }402,}402
hasta
030lma
m = 3,4 mu con los datos que se muestra en el cuadro
. . . _ Z}402 s Fa
sugulente, mantenlendo constante las cantidades - 024,024 024 024g-, 030T
.
can an we M
\ 33
VI. RESULTADOS
TABLA N 6.1
034 .
1 9031 9 = I;
034O2 1:932!" - W0
034E111
Ejij
E
u
j
j
j
j
ET
\ 34
GRAFICA N 6.1
035
0.215798E+O1
o
\ -1.. 030 I, f
\ 034\1 030:,-
-10 V ,
in
now
030 2 4 5 s in 12 14 16 13
\ 35
GRAFICA N° 62
20
M an 024 024
5 030f 03
1: I r 030. )
.» 1 « 2 c, 030
030S % . . \ /I
\ \ V. . 030
.10 . V/4,4, *" 0307' M
\ 030/'
\ 36
GRAFICA N° 6.3
20 /\
AA\\.~v~'v*~~A
f W
° 1: I)))))! 031) ) 030l)
Va 1* 3"! 031 034031
034'
. ' A \
-212 H
_ _-1s__A_ -1o -5 0 5 1a 15 20
X 37
n I// ., h 034N§Cn:gRT|
5 0301 030T;33:::9P»:~§
1 034 .
Q sra}401h}401q}401r
035 5 %
agbtenidoks 030
2° }402
. "5 um»
Ill 030-I 031
\\V"77>
I K
1 030
\,/
V 030
4% 030:-I
"
X 39
g / %%
05$ 035?
g9F 031:Fof 030 /1
\;
4] 030? 20
A9
%
VII. DISCUSION
investigacion.
Amortiguado y Forzado.
\ 41
r
SUBROUTINE derivs(x,y.dydx)
REAL x,y(").dydx('),A,C,E,D
A=0.01
c=o.o2 030 .
D=3O.0
E=0.08
dydx(1)=y(2)
! dydx(2)=-((w1"2) 030y(1))
dydX(2)=-A"Y(2)-Y(1)-C*((Y(1))**3)+D 030C0S(E*X)
I dydx(3)=y(4)
I dydx(4)= 024((w2**2)*v(3))
return
END
! rkdumb.for
SUBROUTINE rkdumb(vstart,nvar,x1,x2,nstep,derivs)
INTEGER nstep,nvar,NMAX,NSTPMX
PARAMETER (NMAX=50,NSTPMX=100000)
EXTERNAL derivs
INTEGER i,k
REAL h,x,dv(NMAX),v(NMAX)
do 11 i=1,nvar
v(i)=vstar1(i)
V(i.1)=V(i)
11 continue
xx(1)=x1
44
031 x=x1 030
h=(x2-x1)/nstep
do 13 k=1,nstep
call derivs(x,v,dv)
call rk4(v,dv,nvar,x,h.v.derivs)
x=x+h
xx(k+1)=x
do 12 i=1,nvar
y(i,k+1)=v(i)
12 continue
13 continue
return
END
! rk4.for
SUBROUTINE rk4(y,dydx,n,x,h,yout,derivs)
INTEGER n,NMAX
EXTERNAL derivs
PARAMETER (NMAX=50)
INTEGER i
REAL h6,hh,xh,dym(NMAX).dyt(NMAX),yt(NMAX)
hh=h"0.5 -
V h6=h/6.
xh=x+hh
do 11 i=1,n
Yt(i)=Y(i)+hh*dYdX(i) '
11 continue
cail derivs(xh,yt,dyt) .
do 12 i=1,n -
Vt(i)=y(i)+hh"dyt(i)
45
12 continue
call derivs(xh,yt,dym)
do 13 i=1,n
Yt(i)=y(i)+h*dYm(i)
dym(i)=dyt(i)+dym(i)
13 continue
call derivs(x+h,yt,dyt)
do 14 i=1,n
14 continue
return
END
1985)
program ode
nn=n"(n+1)=
nn1=nn+1
external fcn
dimension y(nn1),yprime(nn1),v(nn1),A(nn1),B(nn1),C(nn1)
dimension D(nn1),cum(n),znorm(n),gsc(n)
open(1,}401|e="Lout.txt")
read(*,*) nstep,stpsze,irate,io
write(1,*) nstep,stpsze,irate.io
write(1,*)
v(1)=10.0
v(2)=1.0
v(3)=1.0
K 46
tme=0.0
do 10 i=n+1,nn
v(i)=0.0
10 continue
do 20 i=1,n
v((n+1)*i)=1.0
cum(i)=0.0
20 continue
do 100 m=1,nstep
if ((m/100)*100.EQ.m) wn'te(*,*)m
do 25 j=1,irate
do 26 i=1,nn
030/(i)=V(i)
26 continue
t=tme
do 27 i=1,nn
A(i)=yprime(i)
27 continue
do 28 i=1,nn
y(i)=v(i)+(stpsze*A(i))/2.0
28 continue
t=tme+stpszel2.0
do 29 i=1,nn
B(i)=yprime(i)
29 continue
do 30 i=1,nn
y(i)=v(i)+(stpsze'B(i))/2.0
30 continue
t=tme+stpsze/2.0
X 47
call fcn(t, y, yprime)
do 31 i=1,nn
C(i)=yprime(i)
31 continue
031 do 32 i=1,nn
y(i)=V(i)+(StPSZe*C(i))
32 continue
t=tme+stpsze
call fcn(t, y, yprime) 030
do 33 i=1,nn
D(i)=yprime(i)
33 continue
do 34 i=1,nn
v(i)=v(i)+stpsze*(A(i)+D(i)+2.0*(B(i)+C(i)))/6.0
34 continue
tme=tme+stpsze
25 continue
znorm(1)=0.0
do 38 j=1,n
znorm(1)=znorm(1)+v(n*j+1) 030*2
38 continue
znorm(1)=sqrt(znorm(1))
do 40 j=1,n
v(n*j+1)=v(n*j+1)/znorm(1)
40 continue
do 80 j=2,n '
do 50 k=1 ,1-1
gsc(k)=0.0
do so |=1,n _
gsc(k)=gsc(k)+v(n*I+j)*v(n*|+k)
50 continue
\ 48
do 60 k=1,n
do 60 l=1,j-1
v(n*k+j)=v(n*k+j)-gsc(|)"v(n*k+|)
60 continue
znorm(j)=0.0
do 70 k=1,n
znorm(j)=znorm(j)+v(n*k+j)**2
70 continue
znorm(j)=sqrt(znorm(j))
do 80 k=1,n
v(n*k+j)=v(n*k+j)/znorm(j)
80 continue
do 90 k=1,n
cum(k)=cum(k)+d|og(znorm(k))/a|og(2.0)
90 continue
write(1,334) tme,(cum(k)/tme,k=1_n)
334 format(1x,f12.6,2x,e12.6,2x,e12.6,2x,e12.6)
_ 100 continue
end
*m*m**..m*....Mmm...m..........m..............*........*
subroutine fcn(t,y,yprime)
b=4.0
sg=16.0
r=45.92
VPrime(1)=S9*(V(2)-Y(1))
yprime(2)=-y(1)*y(3)+r*y(1)-y(2)
VPrime(3)=Y(1)*Y(2)-b*Y(3)
do 10 i=0,2
vprime(4+i)=S9"(V(7+i)-Y(4+i))
vprime(7+i)=(r-y(3))*v(4+i)-v(7+i)-v(1) 030v(10+i)
\ 49
' 030 vr>rime(10+i)=y(2)*y(4+i)+y(1)*y(7+i)-b*v(10+i)
10 continue
. return
_ . end
031 V i X _ 50
/ X. ANEXOS
Matnz de Consistencia
51