korba 0.7.0 → 0.7.1

checksums.yaml

@@ -1,7 +1,7 @@
 ---
 SHA256:
-  metadata.gz: a8a9dc1d59d577592a4ef91601bb34a4921e2bb1eb731c0fe866698759717d09
-  data.tar.gz: bb572412946fd8f08a00ab844c4a4975d08953871984a142225ffe09b786417d
+  metadata.gz: d032139625c6722891cef771fc76ffb90e8b082f89146ebb860a5f8034004b15
+  data.tar.gz: 354e2f0f7c382059bf7a424ffd952c32b56b6ae616b5c7bd42f50a148d2f42c3
 SHA512:
-  metadata.gz: a6f13b11b06a7dc3159b5219204a3c83791a5708f7fc8d64c0d72734bd751c655b8438ef7d29fbb2599f94dfa6657e9c282235dbf86239fab49ad12a2eafff9d
-  data.tar.gz: bcc2a2d8b20d9accbe8b6079334ce6c4610a22e78363d7c9e6243f39a1608c46a90adc0b7d0e52c53adb9fb7c81b7a4dc60f1a291e523e3166bc54b68f9a6897
+  metadata.gz: 6f73a3a19efff8c89299173d67dbfa1aaf754ce1c0828148f7bd3933cd57d8de8ba74352cc12849fcd405882c964f345353e1f6d4204b4e451d00a8a6ebf5337
+  data.tar.gz: 293456b12a49007bb1f1c1acc6201fcf99678e5427b1c81942d9f1e53ab5c9cc8716ffdeaaa562846c45ef2b245400f7d73758a05494da91df1528da9e06cba0

data/CHANGELOG.md

@@ -1,3 +1,7 @@
+## [0.7.1] - 2026-04-05
+
+- FIX: Fixed calculation when orbital inclination and eccentricity are 0
+
 ## [0.7.0] - 2026-04-01
 
 - FIX: Use SGP4 for TLE to Cartesian conversion

data/lib/korba/kep.rb

@@ -19,13 +19,23 @@ module Korba
     end
 
     def to_car
-      vector_n = Vector[
-        Math.sin(deg_to_rad(inclination)) * Math.sin(deg_to_rad(ra_of_asc_node)),
-        -Math.sin(deg_to_rad(inclination)) * Math.cos(deg_to_rad(ra_of_asc_node)),
-        Math.cos(deg_to_rad(inclination))
-      ]
-      vector_omega = Vector[Math.cos(deg_to_rad(ra_of_asc_node)), Math.sin(deg_to_rad(ra_of_asc_node)), 0]
-      vector_m = vector_n.cross(vector_omega)
+      # 特異点対応: inclinationが0（または極めて小さい）場合の処理
+      is_equatorial = inclination.abs < 1e-10
+
+      if is_equatorial
+        # 赤道面内の場合、x軸を基準にする
+        vector_omega = Vector[1.0, 0.0, 0.0]
+        vector_m = Vector[0.0, 1.0, 0.0]
+        vector_n = Vector[0.0, 0.0, 1.0]
+      else
+        vector_n = Vector[
+          Math.sin(deg_to_rad(inclination)) * Math.sin(deg_to_rad(ra_of_asc_node)),
+          -Math.sin(deg_to_rad(inclination)) * Math.cos(deg_to_rad(ra_of_asc_node)),
+          Math.cos(deg_to_rad(inclination))
+        ]
+        vector_omega = Vector[Math.cos(deg_to_rad(ra_of_asc_node)), Math.sin(deg_to_rad(ra_of_asc_node)), 0.0]
+        vector_m = vector_n.cross(vector_omega)
+      end
 
       r_angle_factor = deg_to_rad(arg_of_pericenter + true_anomaly)
       vector_r = distance * (Math.cos(r_angle_factor) * vector_omega + Math.sin(r_angle_factor) * vector_m)

data/lib/korba/orbit_utils.rb

@@ -4,7 +4,7 @@ module Korba
   module OrbitUtils
     def semi_major_axis
       # a = (μ / n^2)^(1/3) m
-      (Constant::GME / (mean_motion * 2 * Math::PI / 86400.0) ** 2.0) ** (1.0 / 3.0)
+      (Constant::GME / (mean_motion * 2.0 * Math::PI / 86400.0) ** 2.0) ** (1.0 / 3.0)
     end
 
     def height_at_perigee
@@ -21,8 +21,12 @@ module Korba
     end
 
     def true_anomaly
-      factor = (Math.cos(deg_to_rad(eccentric_anomaly)) - eccentricity) / (1 - eccentricity * Math.cos(deg_to_rad(eccentric_anomaly)))
-      rad_to_deg(Math.acos(factor))
+      e_rad = deg_to_rad(eccentric_anomaly)
+      y = Math.sqrt(1 - eccentricity ** 2) * Math.sin(e_rad)
+      x = Math.cos(e_rad) - eccentricity
+
+      # atan2(y, x) で正確な位相（-pi から pi）を求める
+      rad_to_deg(Math.atan2(y, x))
     end
 
     def distance
@@ -30,10 +34,12 @@ module Korba
     end
 
     def velocity
-      Math.sqrt(Constant::GME * (2 / distance - 1 / semi_major_axis))
+      Math.sqrt(Constant::GME * (2.0 / distance - 1.0 / semi_major_axis))
     end
 
     def path_angle
+      return 0.0 if eccentricity < 1e-15
+
       factor = Math.sqrt(Constant::GME * semi_major_axis * (1 - eccentricity ** 2)) / (distance * velocity)
       rad_to_deg(Math.acos(factor))
     end
@@ -44,8 +50,8 @@ module Korba
     end
 
     def normalize_rad(rad)
-      rad = rad + 2 * Math::PI if rad < 0
-      normalize_rad = rad > 2 * Math::PI ? rad - 2 * Math::PI : 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
     end
@@ -56,8 +62,8 @@ module Korba
     end
 
     def normalize_deg(deg)
-      deg = deg + 360 if deg < 0
-      normalized_deg = deg > 360 ? deg - 360 : 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
     end

data/lib/korba/version.rb

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

metadata

@@ -1,7 +1,7 @@
 --- !ruby/object:Gem::Specification
 name: korba
 version: !ruby/object:Gem::Version
-  version: 0.7.0
+  version: 0.7.1
 platform: ruby
 authors:
 - kk0000-kk