First Commit in this repository.

Author: KernelPhilosopher

Simulación orbital/Controllers/Controlador.py

@@ -0,0 +1,9 @@
+class Controlador:
+    def __init__(self, simulador, vista):
+        self.simulador = simulador
+        self.vista = vista
+
+    def ejecutar_simulacion(self):
+        self.simulador.simular()
+        self.vista.configurar_grafico()
+        self.vista.animar()

Simulación orbital/Controllers/__pycache__/Controlador.cpython-313.pyc

@@ -0,0 +1,11 @@
+from Models.EulerSimulador import EulerSimulador
+from Models.VerletSimulador import VerletSimulador
+from Views.VistaOrbital import VistaOrbital
+from Controllers.Controlador import Controlador
+
+# Elegir método (EulerSimulador o VerletSimulador)
+simulador = VerletSimulador(dt=86400, velocidad_inicial=29783)
+vista = VistaOrbital(simulador)
+controlador = Controlador(simulador, vista)
+
+controlador.ejecutar_simulacion()

Simulación orbital/Main.py

@@ -0,0 +1,15 @@
+from .SimuladorBase import SimuladorBase
+
+
+class EulerSimulador(SimuladorBase):
+    def simular(self):
+        t = 0
+        while t < self.t_total:
+            ax, ay = self.calcular_aceleracion(self.x, self.y)
+            self.vx += ax * self.dt
+            self.vy += ay * self.dt
+            self.x += self.vx * self.dt
+            self.y += self.vy * self.dt
+            self.x_vals.append(self.x)
+            self.y_vals.append(self.y)
+            t += self.dt

Simulación orbital/Models/EulerSimulador.py

@@ -0,0 +1,27 @@
+import math
+
+
+class SimuladorBase:
+    G = 6.6748e-11  # Constante gravitacional
+    M_SOL = 1.989e30  # Masa del Sol (kg)
+    UA = 1.496e11  # 1 UA en metros
+
+    def __init__(self, dt=86400, velocidad_inicial=29783):
+        self.dt = dt  # Paso de tiempo (1 día)
+        self.x = self.UA  # Posición inicial en X
+        self.y = 0.0  # Posición inicial en Y
+        self.vx = 0.0  # Velocidad inicial en X
+        self.vy = velocidad_inicial  # Velocidad inicial en Y
+        self.t_total = dt * 365  # 1 año en segundos
+        self.x_vals = []
+        self.y_vals = []
+
+    def calcular_aceleracion(self, x, y):
+        r = math.sqrt(x**2 + y**2)
+        a_mag = -self.G * self.M_SOL / r**3
+        ax = a_mag * x
+        ay = a_mag * y
+        return ax, ay
+
+    def simular(self):
+        raise NotImplementedError("Método simular debe ser implementado en subclases")

Simulación orbital/Models/SimuladorBase.py

@@ -0,0 +1,25 @@
+from .SimuladorBase import SimuladorBase
+
+
+class VerletSimulador(SimuladorBase):
+    def simular(self):
+        t = 0
+        while t < self.t_total:
+            ax, ay = self.calcular_aceleracion(self.x, self.y)
+
+            # Actualizar posiciones
+            x_nuevo = self.x + self.vx * self.dt + 0.5 * ax * self.dt**2
+            y_nuevo = self.y + self.vy * self.dt + 0.5 * ay * self.dt**2
+
+            # Calcular nueva aceleración
+            ax_nuevo, ay_nuevo = self.calcular_aceleracion(x_nuevo, y_nuevo)
+
+            # Actualizar velocidades
+            self.vx += 0.5 * (ax + ax_nuevo) * self.dt
+            self.vy += 0.5 * (ay + ay_nuevo) * self.dt
+
+            self.x = x_nuevo
+            self.y = y_nuevo
+            self.x_vals.append(self.x)
+            self.y_vals.append(self.y)
+            t += self.dt

Simulación orbital/Models/VerletSimulador.py

@@ -0,0 +1,31 @@
+import matplotlib.pyplot as plt
+import matplotlib.animation as animation
+
+
+class VistaOrbital:
+    def __init__(self, simulador):
+        self.simulador = simulador
+        self.fig, self.ax = plt.subplots(figsize=(8, 8))
+
+    def configurar_grafico(self):
+        self.ax.set_xlim(-1.6 * self.simulador.UA, 1.6 * self.simulador.UA)
+        self.ax.set_ylim(-1.6 * self.simulador.UA, 1.6 * self.simulador.UA)
+        self.ax.scatter(0, 0, color="yellow", s=100, label="Sol")
+        self.ax.set_title("Simulación Orbital")
+        self.ax.legend()
+
+    def animar(self):
+        def update(frame):
+            self.ax.clear()
+            self.configurar_grafico()
+            self.ax.plot(
+                self.simulador.x_vals[:frame],
+                self.simulador.y_vals[:frame],
+                color="blue",
+                label="Órbita",
+            )
+
+        ani = animation.FuncAnimation(
+            self.fig, update, frames=len(self.simulador.x_vals), interval=20
+        )
+        plt.show()