calc_sat 0.1.0

checksums.yaml

@@ -0,0 +1,7 @@
+---
+SHA256:
+  metadata.gz: 1b4600e63a5f0471ccd31d140a5d3f9b9286851b428bede929d85463f1649b3c
+  data.tar.gz: 47741ed0974d8d2b70464086d938d7f22e326207411b5c9ce4432110d3cff02a
+SHA512:
+  metadata.gz: fce58a2cb12fbc75e3f4a7f195982f96db841128e0e0c5063dde9f654fa69a5fb68b4e639dd32995a71aac498f5ed939c961344155ba2b2d566806b6ed3facf8
+  data.tar.gz: d6fde12a94e37c0fabde245eb9d0d10aeeeb3d60e3a4eee9ea9e3528a36677da45d5a6c8f82ffc72242af2c2e47b2681bf727a9269108fb2165e950269476764

data/.circleci/config.yml

@@ -0,0 +1,13 @@
+version: 2.1
+jobs:
+  build:
+    docker:
+      - image: ruby:3.0.1
+    steps:
+      - checkout
+      - run:
+          name: Run the default task
+          command: |
+            gem install bundler -v 2.2.17
+            bundle install
+            bundle exec rake

data/.gitignore

@@ -0,0 +1,12 @@
+/.bundle/
+/.yardoc
+/_yardoc/
+/coverage/
+/doc/
+/pkg/
+/spec/reports/
+/tmp/
+Gemfile.lock
+
+# rspec failure tracking
+.rspec_status

data/.rspec

@@ -0,0 +1,3 @@
+--format documentation
+--color
+--require spec_helper

data/.rubocop.yml

@@ -0,0 +1,13 @@
+AllCops:
+  TargetRubyVersion: 2.4
+
+Style/StringLiterals:
+  Enabled: true
+  EnforcedStyle: double_quotes
+
+Style/StringLiteralsInInterpolation:
+  Enabled: true
+  EnforcedStyle: double_quotes
+
+Layout/LineLength:
+  Max: 120

data/CHANGELOG.md

@@ -0,0 +1,5 @@
+## [Unreleased]
+
+## [0.1.0] - 2021-05-23
+
+- Initial release

data/CODE_OF_CONDUCT.md

@@ -0,0 +1,84 @@
+# Contributor Covenant Code of Conduct
+
+## Our Pledge
+
+We as members, contributors, and leaders pledge to make participation in our community a harassment-free experience for everyone, regardless of age, body size, visible or invisible disability, ethnicity, sex characteristics, gender identity and expression, level of experience, education, socio-economic status, nationality, personal appearance, race, religion, or sexual identity and orientation.
+
+We pledge to act and interact in ways that contribute to an open, welcoming, diverse, inclusive, and healthy community.
+
+## Our Standards
+
+Examples of behavior that contributes to a positive environment for our community include:
+
+* Demonstrating empathy and kindness toward other people
+* Being respectful of differing opinions, viewpoints, and experiences
+* Giving and gracefully accepting constructive feedback
+* Accepting responsibility and apologizing to those affected by our mistakes, and learning from the experience
+* Focusing on what is best not just for us as individuals, but for the overall community
+
+Examples of unacceptable behavior include:
+
+* The use of sexualized language or imagery, and sexual attention or
+  advances of any kind
+* Trolling, insulting or derogatory comments, and personal or political attacks
+* Public or private harassment
+* Publishing others' private information, such as a physical or email
+  address, without their explicit permission
+* Other conduct which could reasonably be considered inappropriate in a
+  professional setting
+
+## Enforcement Responsibilities
+
+Community leaders are responsible for clarifying and enforcing our standards of acceptable behavior and will take appropriate and fair corrective action in response to any behavior that they deem inappropriate, threatening, offensive, or harmful.
+
+Community leaders have the right and responsibility to remove, edit, or reject comments, commits, code, wiki edits, issues, and other contributions that are not aligned to this Code of Conduct, and will communicate reasons for moderation decisions when appropriate.
+
+## Scope
+
+This Code of Conduct applies within all community spaces, and also applies when an individual is officially representing the community in public spaces. Examples of representing our community include using an official e-mail address, posting via an official social media account, or acting as an appointed representative at an online or offline event.
+
+## Enforcement
+
+Instances of abusive, harassing, or otherwise unacceptable behavior may be reported to the community leaders responsible for enforcement at OhmachiSeiya@gmail.com. All complaints will be reviewed and investigated promptly and fairly.
+
+All community leaders are obligated to respect the privacy and security of the reporter of any incident.
+
+## Enforcement Guidelines
+
+Community leaders will follow these Community Impact Guidelines in determining the consequences for any action they deem in violation of this Code of Conduct:
+
+### 1. Correction
+
+**Community Impact**: Use of inappropriate language or other behavior deemed unprofessional or unwelcome in the community.
+
+**Consequence**: A private, written warning from community leaders, providing clarity around the nature of the violation and an explanation of why the behavior was inappropriate. A public apology may be requested.
+
+### 2. Warning
+
+**Community Impact**: A violation through a single incident or series of actions.
+
+**Consequence**: A warning with consequences for continued behavior. No interaction with the people involved, including unsolicited interaction with those enforcing the Code of Conduct, for a specified period of time. This includes avoiding interactions in community spaces as well as external channels like social media. Violating these terms may lead to a temporary or permanent ban.
+
+### 3. Temporary Ban
+
+**Community Impact**: A serious violation of community standards, including sustained inappropriate behavior.
+
+**Consequence**: A temporary ban from any sort of interaction or public communication with the community for a specified period of time. No public or private interaction with the people involved, including unsolicited interaction with those enforcing the Code of Conduct, is allowed during this period. Violating these terms may lead to a permanent ban.
+
+### 4. Permanent Ban
+
+**Community Impact**: Demonstrating a pattern of violation of community standards, including sustained inappropriate behavior,  harassment of an individual, or aggression toward or disparagement of classes of individuals.
+
+**Consequence**: A permanent ban from any sort of public interaction within the community.
+
+## Attribution
+
+This Code of Conduct is adapted from the [Contributor Covenant][homepage], version 2.0,
+available at https://www.contributor-covenant.org/version/2/0/code_of_conduct.html.
+
+Community Impact Guidelines were inspired by [Mozilla's code of conduct enforcement ladder](https://github.com/mozilla/diversity).
+
+[homepage]: https://www.contributor-covenant.org
+
+For answers to common questions about this code of conduct, see the FAQ at
+https://www.contributor-covenant.org/faq. Translations are available at https://www.contributor-covenant.org/translations.

data/Gemfile

@@ -0,0 +1,12 @@
+# frozen_string_literal: true
+
+source "https://rubygems.org"
+
+# Specify your gem's dependencies in calc_sat.gemspec
+gemspec
+
+gem "rake", "~> 13.0"
+
+gem "rspec", "~> 3.0"
+
+gem "rubocop", "~> 1.7"

data/LICENSE.txt

@@ -0,0 +1,21 @@
+The MIT License (MIT)
+
+Copyright (c) 2021 OhmachiSeiya
+
+Permission is hereby granted, free of charge, to any person obtaining a copy
+of this software and associated documentation files (the "Software"), to deal
+in the Software without restriction, including without limitation the rights
+to use, copy, modify, merge, publish, distribute, sublicense, and/or sell
+copies of the Software, and to permit persons to whom the Software is
+furnished to do so, subject to the following conditions:
+
+The above copyright notice and this permission notice shall be included in
+all copies or substantial portions of the Software.
+
+THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
+IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
+FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE
+AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
+LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,
+OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN
+THE SOFTWARE.

data/README.md

@@ -0,0 +1,43 @@
+# CalcSat
+
+Welcome to your new gem! In this directory, you'll find the files you need to be able to package up your Ruby library into a gem. Put your Ruby code in the file `lib/calc_sat`. To experiment with that code, run `bin/console` for an interactive prompt.
+
+TODO: Delete this and the text above, and describe your gem
+
+## Installation
+
+Add this line to your application's Gemfile:
+
+```ruby
+gem 'calc_sat'
+```
+
+And then execute:
+
+    $ bundle install
+
+Or install it yourself as:
+
+    $ gem install calc_sat
+
+## Usage
+
+TODO: Write usage instructions here
+
+## 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.
+
+To install this gem onto your local machine, run `bundle exec rake install`. To release a new version, update the version number in `version.rb`, and then run `bundle exec rake release`, which will create a git tag for the version, push git commits and the created tag, and push the `.gem` file to [rubygems.org](https://rubygems.org).
+
+## Contributing
+
+Bug reports and pull requests are welcome on GitHub at https://github.com/[USERNAME]/calc_sat. This project is intended to be a safe, welcoming space for collaboration, and contributors are expected to adhere to the [code of conduct](https://github.com/[USERNAME]/calc_sat/blob/master/CODE_OF_CONDUCT.md).
+
+## License
+
+The gem is available as open source under the terms of the [MIT License](https://opensource.org/licenses/MIT).
+
+## Code of Conduct
+
+Everyone interacting in the CalcSat project's codebases, issue trackers, chat rooms and mailing lists is expected to follow the [code of conduct](https://github.com/[USERNAME]/calc_sat/blob/master/CODE_OF_CONDUCT.md).

data/Rakefile

@@ -0,0 +1,12 @@
+# frozen_string_literal: true
+
+require "bundler/gem_tasks"
+require "rspec/core/rake_task"
+
+RSpec::Core::RakeTask.new(:spec)
+
+require "rubocop/rake_task"
+
+RuboCop::RakeTask.new
+
+task default: %i[spec rubocop]

data/bin/console

@@ -0,0 +1,15 @@
+#!/usr/bin/env ruby
+# frozen_string_literal: true
+
+require "bundler/setup"
+require "calc_sat"
+
+# You can add fixtures and/or initialization code here to make experimenting
+# with your gem easier. You can also use a different console, if you like.
+
+# (If you use this, don't forget to add pry to your Gemfile!)
+# require "pry"
+# Pry.start
+
+require "irb"
+IRB.start(__FILE__)

data/bin/setup

@@ -0,0 +1,8 @@
+#!/usr/bin/env bash
+set -euo pipefail
+IFS=$'\n\t'
+set -vx
+
+bundle install
+
+# Do any other automated setup that you need to do here

data/calc_sat.gemspec

@@ -0,0 +1,37 @@
+# frozen_string_literal: true
+
+require_relative "lib/calc_sat/version"
+
+Gem::Specification.new do |spec|
+  spec.name          = "calc_sat"
+  spec.version       = CalcSat::VERSION
+  spec.authors       = ["OhmachiSeiya"]
+  spec.email         = ["OhmachiSeiya@gmail.com"]
+
+  spec.summary       = "satellite calculation"
+  spec.description   = "satellite orbit calculation, power calculation, etc."
+  spec.homepage      = "https://github.com/OhmachiSeiya/calc_sat"
+  spec.license       = "MIT"
+  spec.required_ruby_version = Gem::Requirement.new(">= 2.4.0")
+
+  spec.metadata["homepage_uri"] = spec.homepage
+  # spec.metadata["source_code_uri"] = "TODO: Put your gem's public repo URL here."
+  # spec.metadata["changelog_uri"] = "TODO: Put your gem's CHANGELOG.md URL here."
+
+  # Specify which files should be added to the gem when it is released.
+  # The `git ls-files -z` loads the files in the RubyGem that have been added into git.
+  spec.files = Dir.chdir(File.expand_path(__dir__)) do
+    `git ls-files -z`.split("\x0").reject { |f| f.match(%r{\A(?:test|spec|features)/}) }
+  end
+  spec.bindir        = "exe"
+  spec.executables   = spec.files.grep(%r{\Aexe/}) { |f| File.basename(f) }
+  spec.require_paths = ["lib"]
+
+  # Uncomment to register a new dependency of your gem
+  # spec.add_dependency "example-gem", "~> 1.0"
+  spec.add_dependency "numo-linalg", "~> 0.1.5"
+  spec.add_dependency "numo-narray", "~>0.9.1.9"
+
+  # For more information and examples about making a new gem, checkout our
+  # guide at: https://bundler.io/guides/creating_gem.html
+end

data/lib/calc_sat.rb

@@ -0,0 +1,177 @@
+# frozen_string_literal: true
+
+require_relative "calc_sat/version"
+
+require 'numo/narray'
+
+module CalcSat
+  # 参考文献：
+  # 「人工衛星をつくる（宮崎康行 著）」
+  # 「CanSat 超小型模擬人工衛星（大学宇宙工学コンソーシアム）」
+  class Orbita
+    # 軌道計算
+    #  a … 楕円軌道の長半径
+    #  e … 楕円の離心率（地球が楕円の中心からどれだけの割合ずれているかを表す）
+    #  i … 軌道傾斜角（楕円軌道面と赤道面がなす角度）
+    #  o … 昇交点赤経
+    #  w … 近地点引数
+    #  theta … 真近点角
+    def initialize()
+      @m_e = 3.986 * 1014  # 地球の重力定数（3.986*1014 [m3/s2]）
+      @r_e = 6378          # 地球の半径(6378 [km])
+      @j_2 = 0.001082628   # 地球の扁平性を表す定数(0.001082628)
+    end
+  
+    def velocity(a,e,theta)
+      # 速度
+      v = Math.sqrt(
+        @m_e * (1 + e**2 + 2*e*Math.cos(theta)) / 
+        a * (1 - e**2)
+      )
+      return v
+    end
+  
+    def periodic_time(a)
+      # 周期
+      t = 2*Math::PI*(a**1.5 / Math.sqrt(@m_e))
+      return t
+    end
+  
+    def sso_check(a, e, i, tolerance=10**-5)
+      # 太陽同期軌道（sun_synchronous_orbit）となる条件式の正誤を判定する
+      # 太陽同期軌道 … iとaを調整して同じ経度に戻るようにした時の軌道
+      # tolerance: 許容誤差
+  
+      comparison = 1.99 * 10**(-7)
+      res = -(
+        (3*@j_2*@r_e**2) / (2*a**2*(1-e**2)**2) *
+        Math.sqrt(@m_e / a**3) *
+        Math.cos(i)
+      )
+      diff = (res - comparison).abs
+      cond = diff < tolerance
+      
+      if cond
+        return "OK"
+      else
+        return "iとaを修正してください。"
+      end
+    end
+  end
+
+  class Communication
+    # 通信
+    # p_sat: 衛星からの電波の電力
+    # l_tr: 衛星からの電波が地上に届くまでのロス
+    # c_p: 搬送波電力
+    # n_0: 雑音電力密度
+    def initialize()
+      @c = 2.9979 * 10**6 # 電波の速さ 2.9979 * 10**6 [km/sec]
+    end
+
+    def line_establised?(p_sat, l_tr, n_0, cn_req)
+      # 通信回線が成立するか判定
+      c_p = p_sat - l_tr
+      cn = c_p / n_0
+      cond = cn - cn_req > 3
+      return cond
+    end
+
+    def fspl(d, f)
+      #自由空間伝搬損失 Free Space Path Loss
+      #電波が地上に届くまでに弱まってしまう割合
+
+      l_d = (@c / 4*Math::PI*f*d)**2
+      return l_d
+    end
+
+    def system_noise_temp(t_a,t_f,t_e,l_frk)
+      t_s = (t_a/l_frk) + (1 - (1/l_frk))*t_f + t_e
+      return t_s
+    end
+  end
+
+  class Posture
+    # 姿勢
+    def torque(m,b_s, y_air,f_air, y_sun,f_sun)
+      # 外乱トルク
+      n_mag = m * b_s
+      n_air = y_air * f_air
+      n_sun = y_sun * f_sun
+      dict = {"n_mag":n_mag, "n_air":n_air, "n_sun":n_sun}
+      return dict
+    end
+  end
+
+  class Kalman
+    # カルマンフィルタ:
+    # 状態方程式を計算して姿勢を予測しつつ、実測データと繰り返し比較して予測値を修正していく
+    include Numo
+  
+    attr_accessor :ft0, :pt1, :qt0, :gt0, :xht1, :ht0, :rt0
+  
+    def predict(ut0)
+      ut0 = 0 if not ut0
+      @xht0 = @ft0.dot(@xht1) + ut0
+      @pt0 = @ft0.dot(@pt1).dot(@ft0.transpose)+@gt0.dot(@qt0).dot(@gt0.transpose)
+      [@xht0, @pt0]
+    end
+  
+    def update(zt0)
+      @xht1 = @xht0 + kt0.dot(et0(zt0))
+      ktht = kt0.dot(@ht0)
+      @pt1 = (Int32.eye(*ktht.shape) - ktht).dot(@pt0)
+      [@xht1, @pt1]
+    end
+  
+    def observe(zt0, ut0=nil)
+      predict(ut0)
+      update(zt0)
+    end
+  
+    def et0(zt0)
+      zt0 - @ht0.dot(@xht0)
+    end
+  
+    def st0
+      @rt0 + @ht0.dot(@pt0).dot(@ht0.transpose)
+    end
+  
+    def kt0
+      @pt0.dot(@ht0.transpose).dot(Matrix[*st0.to_a].inv.to_a)
+    end
+  end
+
+  class Electricity
+    def simple_quote(s_eff, s_cell, p_cell, t_d)
+      # 簡単な見積もりとして、一回の日照時間中の太陽電池パネル全体の発電量u_spを計算する
+      # t_d: 平均日照時間
+      # s_eff: 平均的な面積
+      # s_cell: 太陽電池セル一枚あたりの面積
+      # p_cell: 太陽電池セル一枚あたりの発電電力
+
+      u_sp = (s_eff / s_cell) * p_cell * t_d
+    end
+  end
+
+  class Heat
+    def initialize
+      @sigma = 5.669 * 10**-8 # [W/(m**2/K**4)] ステファン・ボルツマン定数
+    end
+    def thermodynamic_equilibrium_part()
+      # パーツ一つで近似する場合の熱平衡方程式
+      # m: 衛星の質量
+      # c: 衛星材料全体の平均的な比熱
+      # t_capital: 衛星の温度
+      # t: 時間
+      # q_t: 衛星の外部から入ってくるトータルの熱量と衛星地震の発熱の総和
+      # epsilon: 太陽電池パネルの平均的な放射率
+
+      left = m * c * (d*t_capital/ d*t)
+      right = q_t - epsilon * @sigma * a_capital * f_s * t_capital**4
+      cond = left == right
+      return cond
+    end
+  end
+
+end

data/lib/calc_sat/version.rb

@@ -0,0 +1,5 @@
+# frozen_string_literal: true
+
+module CalcSat
+  VERSION = "0.1.0"
+end

metadata

@@ -0,0 +1,87 @@
+--- !ruby/object:Gem::Specification
+name: calc_sat
+version: !ruby/object:Gem::Version
+  version: 0.1.0
+platform: ruby
+authors:
+- OhmachiSeiya
+autorequire:
+bindir: exe
+cert_chain: []
+date: 2021-05-22 00:00:00.000000000 Z
+dependencies:
+- !ruby/object:Gem::Dependency
+  name: numo-linalg
+  requirement: !ruby/object:Gem::Requirement
+    requirements:
+    - - "~>"
+      - !ruby/object:Gem::Version
+        version: 0.1.5
+  type: :runtime
+  prerelease: false
+  version_requirements: !ruby/object:Gem::Requirement
+    requirements:
+    - - "~>"
+      - !ruby/object:Gem::Version
+        version: 0.1.5
+- !ruby/object:Gem::Dependency
+  name: numo-narray
+  requirement: !ruby/object:Gem::Requirement
+    requirements:
+    - - "~>"
+      - !ruby/object:Gem::Version
+        version: 0.9.1.9
+  type: :runtime
+  prerelease: false
+  version_requirements: !ruby/object:Gem::Requirement
+    requirements:
+    - - "~>"
+      - !ruby/object:Gem::Version
+        version: 0.9.1.9
+description: satellite orbit calculation, power calculation, etc.
+email:
+- OhmachiSeiya@gmail.com
+executables: []
+extensions: []
+extra_rdoc_files: []
+files:
+- ".circleci/config.yml"
+- ".gitignore"
+- ".rspec"
+- ".rubocop.yml"
+- CHANGELOG.md
+- CODE_OF_CONDUCT.md
+- Gemfile
+- LICENSE.txt
+- README.md
+- Rakefile
+- bin/console
+- bin/setup
+- calc_sat.gemspec
+- lib/calc_sat.rb
+- lib/calc_sat/version.rb
+homepage: https://github.com/OhmachiSeiya/calc_sat
+licenses:
+- MIT
+metadata:
+  homepage_uri: https://github.com/OhmachiSeiya/calc_sat
+post_install_message:
+rdoc_options: []
+require_paths:
+- lib
+required_ruby_version: !ruby/object:Gem::Requirement
+  requirements:
+  - - ">="
+    - !ruby/object:Gem::Version
+      version: 2.4.0
+required_rubygems_version: !ruby/object:Gem::Requirement
+  requirements:
+  - - ">="
+    - !ruby/object:Gem::Version
+      version: '0'
+requirements: []
+rubygems_version: 3.2.15
+signing_key:
+specification_version: 4
+summary: satellite calculation
+test_files: []