Oscilador de duffing forzado - coficacin en fortran 90
?
0 likes?569 views
The document describes a Fortran 90 program that simulates an oscillator. It defines variables like arrays to store time, position, and velocity values. It initializes the oscillator at t=0.1, x=0.3, and y=0.1. In a loop, it calculates the next time step and updates position and velocity values based on the oscillator equations. The values are written to a file and screen at each step. Graphs are shown of position vs time for different values of the Z parameter, showing the oscillator behavior changes as Z increases above 1.
Oscilador de duffing forzado - coficacin en fortran 90
- 1. Jose Luis Leon Aguirre Cdigo fortran 90 PROGRAM OSCILADOR INTEGER::I REAL::X(10000),Y(10000),T(100000),H,E,F,Z OPEN(1,FILE='DATA4.TXT') X(0)=0.3 Y(0)=0.1 T(0)=0.1 H=0.01 E=1.0 F=1.3 Z=0.6 DO I=1,10000 T(I+1)=T(I)+H Y(I+1)=Y(I)+(-X(I)-E*(X(I)**3)+F*SIN(Z*T(I)))*H X(I+1)=X(I)+Y(I)*H WRITE(*,*)T(I),X(I),Y(I) WRITE(1,*)T(I),X(I),Y(I) END DO END
- 2. -2.50E+00 -2.00E+00 -1.50E+00 -1.00E+00 -5.00E-01 0.00E+00 5.00E-01 1.00E+00 1.50E+00 2.00E+00 2.50E+00 -2.00E+00-1.50E+00-1.00E+00-5.00E-010.00E+005.00E-011.00E+001.50E+002.00E+00 Z=0.7 -2.00E+00 -1.50E+00 -1.00E+00 -5.00E-01 0.00E+00 5.00E-01 1.00E+00 1.50E+00 2.00E+00 -2.00E+00-1.50E+00-1.00E+00 -5.00E-01 0.00E+00 5.00E-01 1.00E+00 1.50E+00 2.00E+00 Z=0.6
- 3. Anlisis tiempo vs distancia -2.50E+00 -2.00E+00 -1.50E+00 -1.00E+00 -5.00E-01 0.00E+00 5.00E-01 1.00E+00 1.50E+00 2.00E+00 2.50E+00 -2.00E+00-1.50E+00-1.00E+00 -5.00E-01 0.00E+00 5.00E-01 1.00E+00 1.50E+00 2.00E+00 Z=0.5 -2.00E+00 -1.50E+00 -1.00E+00 -5.00E-01 0.00E+00 5.00E-01 1.00E+00 1.50E+00 2.00E+00 0.00E+00 2.00E+01 4.00E+01 6.00E+01 8.00E+01 1.00E+02 1.20E+02 Z=0.5
- 4. Anlisis tiempo vs distancia Cuando se toma valores de Z=1.1 para ms arriba los valores toman grandes valores -8.00E+00 -6.00E+00 -4.00E+00 -2.00E+00 0.00E+00 2.00E+00 4.00E+00 6.00E+00 8.00E+00 -4.00E+00-3.00E+00-2.00E+00-1.00E+00 0.00E+00 1.00E+00 2.00E+00 3.00E+00 4.00E+00 Z=1.0 -4.00E+00 -3.00E+00 -2.00E+00 -1.00E+00 0.00E+00 1.00E+00 2.00E+00 3.00E+00 4.00E+00 0.00E+00 2.00E+01 4.00E+01 6.00E+01 8.00E+01 1.00E+02 1.20E+02 Z=1.0
- 5. Anlisis tiempo vs distancia -1.60E+22 -1.40E+22 -1.20E+22 -1.00E+22 -8.00E+21 -6.00E+21 -4.00E+21 -2.00E+21 0.00E+00 2.00E+21 -1.00E+13 0.00E+00 1.00E+13 2.00E+13 3.00E+13 Z=1.1 -5.00E+12 0.00E+00 5.00E+12 1.00E+13 1.50E+13 2.00E+13 2.50E+13 0.00E+00 2.00E+01 4.00E+01 6.00E+01 8.00E+01 1.00E+02 Z=1.1