korba 0.4.0 → 0.5.0

checksums.yaml

@@ -1,7 +1,7 @@
 ---
 SHA256:
-  metadata.gz: e02130eb6662759d8869d606d53c5d7ecc64c26ed59a5605487fdf5e8d785f01
-  data.tar.gz: 90d8fe4257842deffd1f45ca5fab34013f9decb275ff03f9d59c369f3182e11e
+  metadata.gz: 4d844574b98c6fe33241edac6be121c39c8a509b405ab1d5478aca0de1426465
+  data.tar.gz: 72ba006521490f7e0670ebb800457b34738df3ebbc6fd75cdcfa56259816610e
 SHA512:
-  metadata.gz: 95b4ade626d6622599c3d7e6b533db46a49b820c81cdbc3a8f57838901220950dfdb2643de622d8a414c63a4a24fa15831a8dc489c821ec15c1d10335a929d6c
-  data.tar.gz: 84b650670d1d10fc81b0cd9e9b5005a2dba965bb119df8913a7e38d0997b0c0f31f53ba105d9c3d13198ae4b768df5c7f23b5eb8e3d4694ea86bc7bf0b3e1ff6
+  metadata.gz: 3f169b773a67dd872046d08f605c7c6319930a069133c4b8732a8e7bf6953d0a64fd065b75404990cf69ac2d2892e014f6977385ae07b652e20267e2b1ebcc3e
+  data.tar.gz: 135c8c52c2cc1b25aa543c087ec4a6c9f3ac784c3042864398363f5258b08a9e1e0463dc2a6afe6dc150765e1183d6131e9bc32db39d5fb7b362c0f48055c6a4

data/CHANGELOG.md

@@ -1,3 +1,7 @@
+## [0.5.0] - 2024-12-15
+
+- Enable transformation to Cartesian elements
+
 ## [0.4.0] - 2024-12-09
 
 - Enable save tle_string

data/README.md

@@ -38,7 +38,7 @@ orbit.height_at_perigee
 orbit.height_at_apogee
 # => 420372.03623388056
 
-orbit.to_kep
+kep = orbit.to_kep
 # =>
 # #<Korba::Kep:xxxxxxxxxxxxxxxx
 #  @arg_of_pericenter=314.0303,
@@ -49,8 +49,28 @@ orbit.to_kep
 #  @object_name="ISS (ZARYA)",
 #  @ra_of_asc_node=174.9565,
 #  @semi_major_axis=6793877.649839985>
+
+kep.to_car
+# =>
+# #<Korba::Car:xxxxxxxxxxxxxxxx
+#  @epoch="2024-12-07T20:37:24.085055",
+#  @object_name="ISS (ZARYA)",
+#  @vx=6150.772410883998,
+#  @vy=2489.3298780751356,
+#  @vz=-3816.0301666253677,
+#  @x=4019753.8621700387,
+#  @y=-3623966.518673545,
+#  @z=4114361.6934797494>
+
 ```
 
+## References
+
+- 宇宙工学シリーズ 3 人工衛星と宇宙探査機
+  - https://www.coronasha.co.jp/np/isbn/9784339012316/
+- 人工衛星の軌道 概論
+  - https://www.coronasha.co.jp/np/isbn/9784339046403/
+
 ## Development
 
 After checking out the repo, run `bin/setup` to install dependencies. Then, run `rake spec` to run the tests. You can also run `bin/console` for an interactive prompt that will allow you to experiment.

data/lib/korba/car.rb

@@ -0,0 +1,18 @@
+# frozen_string_literal: true
+
+module Korba
+  class Car
+    attr_reader :object_name, :epoch, :x, :y, :z, :vx, :vy, :vz
+
+    def initialize(object_name:, epoch:, x:, y:, z:, vx:, vy:, vz:)
+      @object_name = object_name
+      @epoch = epoch
+      @x = x
+      @y = y
+      @z = z
+      @vx = vx
+      @vy = vy
+      @vz = vz
+    end
+  end
+end

data/lib/korba/kep.rb

@@ -1,9 +1,9 @@
 # frozen_string_literal: true
-require_relative "orbit_utils"
+require "matrix"
 
 module Korba
   class Kep
-    include Korba::OrbitUtils
+    include OrbitUtils
 
     attr_reader :object_name, :epoch, :semi_major_axis, :eccentricity, :inclination, :ra_of_asc_node, :arg_of_pericenter, :mean_anomaly
 
@@ -17,5 +17,22 @@ module Korba
       @arg_of_pericenter = arg_of_pericenter
       @mean_anomaly = mean_anomaly
     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)
+
+      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)
+      v_angle_factor = deg_to_rad(arg_of_pericenter + true_anomaly - path_angle)
+      vector_v = velocity * (-Math.sin(v_angle_factor) * vector_omega + (Math.cos(v_angle_factor)) * vector_m)
+
+      Car.new(object_name:, epoch:, x: vector_r[0], y: vector_r[1], z: vector_r[2], vx: vector_v[0], vy: vector_v[1], vz: vector_v[2])
+    end
   end
 end

data/lib/korba/kepler_equation_fucntion.rb

@@ -0,0 +1,20 @@
+# frozen_string_literal: true
+require_relative "newton_function"
+
+module Korba
+  class KeplerEquationFunction < NewtonFunction
+    include OrbitUtils
+
+    attr_reader :eccentricity, :mean_anomaly
+
+    def initialize(eccentricity:, mean_anomaly:)
+      super()
+      @eccentricity = eccentricity
+      @mean_anomaly = mean_anomaly
+    end
+
+    def values(x)
+      [x[0] - eccentricity * Math.sin(x[0]) - deg_to_rad(mean_anomaly)]
+    end
+  end
+end

data/lib/korba/newton_function.rb

@@ -0,0 +1,25 @@
+# frozen_string_literal: true
+require "bigdecimal/newton"
+include Newton
+
+module Korba
+  class NewtonFunction
+    def initialize()
+      @zero = BigDecimal("0.0")
+      @one = BigDecimal("1.0")
+      @two = BigDecimal("2.0")
+      @ten = BigDecimal("10.0")
+      @eps = BigDecimal("1.0e-14")
+    end
+
+    def zero; @zero; end
+    def one; @one; end
+    def two; @two; end
+    def ten; @ten; end
+    def eps; @eps; end
+
+    def values(x)
+      raise NotImplementedError
+    end
+  end
+end

data/lib/korba/orbit_utils.rb

@@ -4,15 +4,64 @@ module Korba
   module OrbitUtils
     def semi_major_axis
       # a = (μ / n^2)^(1/3) m
-      (Korba::Constant::GME / (mean_motion * 2 * Math::PI / 86400.0) ** 2.0) ** (1.0 / 3.0)
+      (Constant::GME / (mean_motion * 2 * Math::PI / 86400.0) ** 2.0) ** (1.0 / 3.0)
     end
 
     def height_at_perigee
-      semi_major_axis * (1 - eccentricity) - Korba::Constant::EARTH_RADIUS
+      semi_major_axis * (1 - eccentricity) - Constant::EARTH_RADIUS
     end
 
     def height_at_apogee
-      semi_major_axis * (1 + eccentricity) - Korba::Constant::EARTH_RADIUS
+      semi_major_axis * (1 + eccentricity) - Constant::EARTH_RADIUS
+    end
+
+    def eccentric_anomaly
+      f = KeplerEquationFunction.new(eccentricity:, mean_anomaly:)
+      x = [deg_to_rad(mean_anomaly)]
+      nlsolve(f, x)
+      rad_to_deg(x[0])
+    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))
+    end
+
+    def distance
+      semi_major_axis * (1 - eccentricity * Math.cos(deg_to_rad(eccentric_anomaly)))
+    end
+
+    def velocity
+      Math.sqrt(Constant::GME * (2 / distance - 1 / semi_major_axis))
+    end
+
+    def path_angle
+      factor = Math.sqrt(Constant::GME * semi_major_axis * (1 - eccentricity ** 2)) / (distance * velocity)
+      rad_to_deg(Math.acos(factor))
+    end
+
+    def deg_to_rad(deg)
+      rad = deg * Math::PI / 180.0
+      normalize_rad(rad)
+    end
+
+    def normalize_rad(rad)
+      rad = rad + 2 * Math::PI if rad < 0
+      normalize_rad = rad > 2 * Math::PI ? rad - 2 * Math::PI : rad
+      normalize_rad(normalize_rad) if normalize_rad != rad
+      normalize_rad
+    end
+
+    def rad_to_deg(rad)
+      deg = rad * 180.0 / Math::PI
+      normalize_deg(deg)
+    end
+
+    def normalize_deg(deg)
+      deg = deg + 360 if deg < 0
+      normalized_deg = deg > 360 ? deg - 360 : deg
+      normalize_deg(normalized_deg) if normalized_deg != deg
+      normalized_deg
     end
   end
 end

data/lib/korba/tle.rb

@@ -3,7 +3,7 @@ require_relative "orbit_utils"
 
 module Korba
   class Tle
-    include Korba::OrbitUtils
+    include OrbitUtils
 
     attr_reader :tle_json, :tle_string, :object_name, :object_id, :epoch, :mean_motion, :eccentricity, :inclination, :ra_of_asc_node, :arg_of_pericenter, :mean_anomaly
 
@@ -19,14 +19,19 @@ module Korba
     end
 
     def to_kep
-      Korba::Kep.new(object_name:,
-                     epoch:,
-                     semi_major_axis:,
-                     eccentricity:,
-                     inclination:,
-                     ra_of_asc_node:,
-                     arg_of_pericenter:,
-                     mean_anomaly:)
+      Kep.new(object_name:,
+              epoch:,
+              semi_major_axis:,
+              eccentricity:,
+              inclination:,
+              ra_of_asc_node:,
+              arg_of_pericenter:,
+              mean_anomaly:)
+    end
+
+    def to_car
+      kep = to_kep
+      kep.to_car
     end
 
     private

data/lib/korba/version.rb

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

data/lib/korba.rb

@@ -4,6 +4,9 @@ require_relative "korba/version"
 require_relative "korba/constant"
 require_relative "korba/tle"
 require_relative "korba/kep"
+require_relative "korba/car"
+require_relative "korba/orbit_utils"
+require_relative "korba/kepler_equation_fucntion"
 
 module Korba
 end

metadata

@@ -1,14 +1,14 @@
 --- !ruby/object:Gem::Specification
 name: korba
 version: !ruby/object:Gem::Version
-  version: 0.4.0
+  version: 0.5.0
 platform: ruby
 authors:
 - kk0000-kk
 autorequire:
 bindir: exe
 cert_chain: []
-date: 2024-12-09 00:00:00.000000000 Z
+date: 2024-12-15 00:00:00.000000000 Z
 dependencies: []
 description: Orbital elements calculation library
 email:
@@ -22,8 +22,11 @@ files:
 - README.md
 - Rakefile
 - lib/korba.rb
+- lib/korba/car.rb
 - lib/korba/constant.rb
 - lib/korba/kep.rb
+- lib/korba/kepler_equation_fucntion.rb
+- lib/korba/newton_function.rb
 - lib/korba/orbit_utils.rb
 - lib/korba/tle.rb
 - lib/korba/version.rb
@@ -50,7 +53,7 @@ required_rubygems_version: !ruby/object:Gem::Requirement
     - !ruby/object:Gem::Version
       version: '0'
 requirements: []
-rubygems_version: 3.5.9
+rubygems_version: 3.5.23
 signing_key:
 specification_version: 4
 summary: Orbital elements calculation library