korba 0.8.0 → 0.9.0

checksums.yaml

@@ -1,7 +1,7 @@
 ---
 SHA256:
-  metadata.gz: 5dc3a4abbdd189f4388673e3960a615940a4feb68c359cce37a29ff47fa7efcb
-  data.tar.gz: c312b2be706e321242d0f4c762134ad82ea681348e5631ec36a13ed91c36db11
+  metadata.gz: 82d326cac0ea342432bbb082f3abf220a6104073c070bc9c634bab9bb981c388
+  data.tar.gz: 0254dd92e20901966e424acf29a1c1cdc4c44128499e2c85b3fac3cb10fe89f9
 SHA512:
-  metadata.gz: f0cc15b9fa16cec125d3902e3b52e1336acf4246c72860049955c32456d1eaa0e56a2fd99953477d36141da23b3b8696ed0b8a611c80bc51e1249368c7f07f84
-  data.tar.gz: ffcbf5ed21c86f069d83276227c9db29e5b59d257d1630f5a1fb15eaff0550b881b1e7753fcc551ce6c647c49faf2bc84ee52efb015fbd9ad8efa9355b2da6ad
+  metadata.gz: a8c434b74a17a64f5ad4367a41391ec85345021f1b1de5ad741b9632a49ae92d123eb8aac9dd22e262a01953012862862a396fe9e5b6429dd5a126412d813bc4
+  data.tar.gz: f6e6c90bc931c98854823baa0089b6dff556b70def4a1a0d72f1544e633cc9745097506b00f422812ddbb4f598d640df53d18dfd6576145a2b19361cdc75270c

data/CHANGELOG.md

@@ -1,3 +1,7 @@
+## [0.9.0] - 2026-05-04
+
+- Enable orbit propagation with Keplerian elements
+
 ## [0.8.0] - 2026-04-19
 
 - Enable orbit propagation using 4th-order Runge-Kutta

data/lib/korba/orbit_utils.rb

@@ -7,6 +7,10 @@ module Korba
       (Constant::GME / (mean_motion * 2.0 * Math::PI / 86400.0) ** 2.0) ** (1.0 / 3.0)
     end
 
+    def period
+      2.0 * Math::PI * Math.sqrt(semi_major_axis ** 3 / Constant::GME)
+    end
+
     def height_at_perigee
       semi_major_axis * (1 - eccentricity) - Constant::EARTH_RADIUS
     end
@@ -50,10 +54,7 @@ module Korba
     end
 
     def normalize_rad(rad)
-      rad = rad + 2.0 * Math::PI if rad < 0
-      normalize_rad = rad > 2.0 * Math::PI ? rad - 2.0 * Math::PI : rad
-      normalize_rad(normalize_rad) if normalize_rad != rad
-      normalize_rad
+      rad % (2.0 * Math::PI)
     end
 
     def rad_to_deg(rad)
@@ -62,10 +63,7 @@ module Korba
     end
 
     def normalize_deg(deg)
-      deg = deg + 360.0 if deg < 0
-      normalized_deg = deg > 360.0 ? deg - 360.0 : deg
-      normalize_deg(normalized_deg) if normalized_deg != deg
-      normalized_deg
+      deg % 360.0
     end
   end
 end

data/lib/korba/propagator/kepler.rb

@@ -0,0 +1,54 @@
+module Korba
+  module Propagator
+    class Kepler
+      include OrbitUtils
+
+      def initialize(initial_kep, disable_j2: false)
+        @initial_kep = initial_kep
+        @disable_j2 = disable_j2
+      end
+
+      def propagate(seconds_after_epoch)
+        mean_motion = Math.sqrt(Constant::GME / @initial_kep.semi_major_axis ** 3)
+        mean_anomaly = rad_to_deg(deg_to_rad(@initial_kep.mean_anomaly) + mean_motion * seconds_after_epoch)
+
+        Kep.new(
+          object_name: @initial_kep.object_name,
+          epoch: @initial_kep.epoch + seconds_after_epoch,
+          semi_major_axis: @initial_kep.semi_major_axis,
+          eccentricity: @initial_kep.eccentricity,
+          inclination: @initial_kep.inclination,
+          ra_of_asc_node: ra_of_asc_node(seconds_after_epoch),
+          arg_of_pericenter: arg_of_pericenter(seconds_after_epoch),
+          mean_anomaly: mean_anomaly,
+        )
+      end
+
+      private
+
+      def ra_of_asc_node(seconds_after_epoch)
+        return @initial_kep.ra_of_asc_node if @disable_j2
+
+        rad_to_deg(deg_to_rad(@initial_kep.ra_of_asc_node) + ra_of_asc_node_dot_radian * seconds_after_epoch)
+      end
+
+      def ra_of_asc_node_dot_radian
+        -(3 * Math::PI * Constant::J2 *
+          (Constant::EARTH_RADIUS / @initial_kep.semi_major_axis * (1 - @initial_kep.eccentricity ** 2)) ** 2 *
+          Math.cos(deg_to_rad(@initial_kep.inclination))) / @initial_kep.period
+      end
+
+      def arg_of_pericenter(seconds_after_epoch)
+        return @initial_kep.arg_of_pericenter if @disable_j2
+
+        rad_to_deg(deg_to_rad(@initial_kep.arg_of_pericenter) + arg_of_pericenter_dot_radian * seconds_after_epoch)
+      end
+
+      def arg_of_pericenter_dot_radian
+        (3 * Math::PI / 2.0 * Constant::J2 *
+         (Constant::EARTH_RADIUS / @initial_kep.semi_major_axis * (1 - @initial_kep.eccentricity ** 2)) ** 2 *
+         (5 * (Math.cos(deg_to_rad(@initial_kep.inclination)) ** 2) - 1)) / @initial_kep.period
+      end
+    end
+  end
+end

data/lib/korba/version.rb

@@ -1,5 +1,5 @@
 # frozen_string_literal: true
 
 module Korba
-  VERSION = "0.8.0"
+  VERSION = "0.9.0"
 end

data/lib/korba.rb

@@ -10,6 +10,7 @@ require_relative "korba/keplers_equation"
 require_relative "korba/sgp4/elset_rec"
 require_relative "korba/propagator/sgp4"
 require_relative "korba/propagator/rk4"
+require_relative "korba/propagator/kepler"
 
 module Korba
 end

metadata

@@ -1,7 +1,7 @@
 --- !ruby/object:Gem::Specification
 name: korba
 version: !ruby/object:Gem::Version
-  version: 0.8.0
+  version: 0.9.0
 platform: ruby
 authors:
 - kk0000-kk
@@ -26,6 +26,7 @@ files:
 - lib/korba/kep.rb
 - lib/korba/keplers_equation.rb
 - lib/korba/orbit_utils.rb
+- lib/korba/propagator/kepler.rb
 - lib/korba/propagator/rk4.rb
 - lib/korba/propagator/sgp4.rb
 - lib/korba/sgp4/elset_rec.rb