STK中的光照计算模型

本文简要阐述STK中光照计算的模型。

在航天任务中,通常需要分析地面站、飞行器在一定时间内的光照情况,具体包括:

  1. 地面站处在光照区和阴影区的具体时间范围;
  2. 考虑地形遮挡后,地面站的光照区和阴影区的变化情况;
  3. 飞行器绕地飞行过程中,处于光照区和阴影区的具体时间范围;
  4. 地面站当地水平系下,太阳方位角、高度角的变化;
  5. 飞行器轨道器或本体系下,太阳方位角、高度角的变化。

光照计算时,主要考虑的是地球对太阳的遮挡(其它天体也是类似的),见下图所示地面站和飞行器的光照示意图。

对于地面站(左图),考虑地球遮挡,当太阳在当地水平面之上时(对应的太阳的高度角大于0),即地面站为光照状态;反之则为阴影状态。

对于空间飞行器(上面右图),如卫星,当飞行器与太阳连线不被地球遮挡时,则飞行器为光照状态;反之则为阴影状态。

实际计算时,地球形状考虑为椭球体。

下面根据以上所涉及到的光照场景,详细讨论。

地面站和飞行器的光照示意图

坐标系及太阳方位角、高度角

涉及到太阳位置计算时,通常涉及到太阳的方位角和高度角概念,因此首先确定好常用的坐标系和太阳方位角和高度角的定义。

地面站地平坐标系

地面站观测太阳时,采用当地水平坐标系(也称北东地坐标系),简称LH坐标系(Local Horizontal),其定义如下:

  1. X轴指向当地北方向;
  2. Y轴指向当地东方向;
  3. Z轴指向当地天底方向。

XY平面为当地水平面,垂直于地球椭球体法线,见下图。
地面站当地水平坐标系定义
下图给出了地面站LH系下的太阳方位角和高度角。
地面站地平坐标系下的太阳方位角和高度角

飞行器VVLH系

飞行器观测太阳时,通常采用轨道坐标系,简称VVLH坐标系(Vehicle Velocity, Local Horizontal),其定义如下:

  1. X轴约束在惯性系速度方向(由Y叉乘Z得到);
  2. Y轴指向轨道面负法向;
  3. Z轴指向地心方向。

下图为飞行器的VVLH坐标系以及太阳方向矢量的高度角和方位角示意图。

飞行器VVLH系下的太阳方位角和高度角

太阳方位角和高度角

在地面站地平坐标系或者飞行器VVLH坐标系中,太阳方位角(Azimuth,简称Az)和高度角(Elevation,简称El)的定义如下图。
LH和VVLH坐标系下的太阳方位角和高度角定义

  1. 方位角定义为:X轴与太阳方向矢量在XY平面内的投影矢量的夹角,+X轴为零点,向+Y轴方向为正;
  2. 高度角定义为:太阳方向矢量与XY平面的夹角,-Z轴方向为正。

注意,对于飞行器本体坐标系(Body),则高度角以+Z轴方向为正。

根据以上定义,实际计算时,首先求得某时刻太阳(通常为视太阳)在地面站LH系或飞行器VVLH系的位置,设为 R s \textbf{R}_s Rs
R s = [ X s , Y s , Z s ] T \textbf{R}_s=\left[X_s,Y_s,Z_s\right]^T Rs=[Xs,Ys,Zs]T

则太阳方位角Az和高度角El计算如下:
{ A z = t a n 2 − 1 ( Y s , X s ) E l = s i n − 1 ( − Z s / R s ) \left\{ \begin{aligned} Az=tan2^{-1}(Y_s,X_s) \\ El=sin^{-1}(-Z_s/R_s) \end{aligned} \right. {Az=tan21(Ys,Xs)El=sin1(Zs/Rs)

光照计算模型

无论是地面站还是飞行器,在精确计算光照和阴影的时间时,必须考虑到以下因素:

  1. 太阳圆盘的大小,以及被遮挡的部分大小;
  2. 对于地面站,当地水平面附近地形遮挡的影响;
  3. 对于空间飞行器,考虑地球的遮挡。

同时,太阳的光照状态分为以下三种情形:

  1. 光照:太阳圆盘完全不被遮挡,地面站或飞行器处于完全光照状态,太阳光照强度因子为1;
  2. 半影:太阳圆盘部分被遮挡,地面站或飞行器处于半影状态,太阳光照强度因为0-1之间的小数;
  3. 全影:太阳圆盘完全被遮挡,地面站或飞行器处于完全阴影状态,太阳光照强度因子为0。

首先给出地面站的光照计算模型,见下图。

某时刻,以观测点为中心,视太阳方向的地形最大仰角为 α t \alpha_t αt,视太阳中心(即太阳位置)方向与地形最大仰角方向的夹角为 α g \alpha_g αg,太阳圆盘视半径为 α s \alpha_s αs

不考虑地形时( α t = 0 \alpha_t=0 αt=0), α g \alpha_g αg即为视太阳方向与当地水平面的夹角,即太阳仰角El。

太阳圆盘视半径 α s \alpha_s αs由下式给出:
α s = s i n − 1 R o R s \alpha_s=sin^{-1}{\frac{R_o}{R_s}} αs=sin1RsRo
上式中, R o R_o Ro为太阳圆盘半径,取值为695700km,对应的视半径约为0.27°,具体数值与太阳的距离变化而稍有不同。
地面目标的光照计算模型
下图为飞行器的光照计算模型示意图。与地面站不同的是,不需要考虑地形的遮挡,转而考虑地球的遮挡。

某时刻,以观测点为中心,地心方向与地球边缘方向的夹角(称为地球视半径)为 α t \alpha_t αt,视太阳中心(即太阳位置)方向与地球边缘方向的夹角为 α g \alpha_g αg,太阳圆盘视半径为 α s \alpha_s αs。观测点到地球的距离向量为 R E \textbf{R}_E RE
飞行器目标的光照计算模型
地球视半径 α t \alpha_t αt可由下式给出:
α t = s i n − 1 R e R E \alpha_t=sin^{-1}{\frac{R_e}{R_E}} αt=sin1RERe
上式中, R e R_e Re为地球赤道半径,常取6378.14km。

根据上述两种光照模型,太阳光照状态的判别依据如下:
{ 光照 : α g > α s 阴影 : α g < − α s 半影 : − α s ≤ α g ≤ α s \left\{ \begin{aligned} 光照: & \alpha_g>\alpha_s \\ 阴影: & \alpha_g<-\alpha_s \\ 半影: & -\alpha_s \leq \alpha_g\leq\alpha_s \end{aligned} \right. 光照:阴影:半影:αg>αsαg<αsαsαgαs

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.rhkb.cn/news/328638.html

如若内容造成侵权/违法违规/事实不符,请联系长河编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

MYSQL-9.问题排查

问题排查的思路与方向 问题排查思路 分析问题&#xff1a;根据理论知识经验分析问题&#xff0c;判断问题可能出现的位置或可能引起问题的原因&#xff0c;将目标缩小到一定范围&#xff1b;排查问题&#xff1a;基于上一步的结果&#xff0c;从引发问题的“可疑性”角度出发…

低空经济:无人机竞赛详解

无人机竞赛市场近年来呈现出蓬勃发展的态势&#xff0c;其市场价值不仅体现在竞赛本身&#xff0c;还体现在推动无人机技术创新、拓展应用场景以及促进产业链发展等多个方面。 一、比赛项目介绍 无人机竞赛通常分为多个项目&#xff0c;包括竞速赛、技巧赛、航拍赛等。每个项目…

新手也能看懂的前端单元测试框架:Vitest

单元测试的概念及作用 1.什么是单元测试&#xff1f; 单元测试是测试中的一个重要环节&#xff0c;它针对软件中的最小可测试单元进行验证&#xff0c;通常是指对代码中的单个函数、方法或模块进行测试。 单元测试旨在确定特定部分代码的行为是否符合预期&#xff0c;通过针…

【谷粒商城】01-环境准备

1.下载和安装VirtualBox 地址&#xff1a;https://www.virtualbox.org/wiki/Downloads 傻瓜式安装VirtualBox 2.下载和安装Vagrant官方镜像 地址&#xff1a;https://app.vagrantup.com/boxes/search 傻瓜式安装 验证是否安装成功 打开CMD,输入vagrant命令&#xff0c;是否…

Linux的常用指令 和 基础知识穿插巩固(巩固知识必看)

目录 前言 ls ls 扩展知识 ls -l ls -a ls -al cd cd 目录名 cd .. cd ~ cd - pwd 扩展知识 路径 / cp [选项] “源文件名” “目标文件名” mv [选项] “源文件名” “目标文件名” rm 作用 用法 ./"可执行程序名" mkdir rmdir touch m…

单位个人怎样向报社的报纸投稿?

作为一名单位的信息宣传员,我肩负着每月定期在媒体上投稿发表文章的重任。然而,在投稿的道路上,我经历了不少波折和挫折。 一开始,我天真地以为只要将稿件发送到报社的投稿邮箱,就能轻松完成任务。然而,现实却远比我想象的复杂。邮箱投稿的竞争异常激烈,编辑们会在众多稿件中挑…

什么是直接内存(NIO)

直接内存不受 JVM 内存回收管理&#xff0c;是虚拟机的系统内存&#xff0c;常见于 NIO 操作时&#xff0c;用于数据 缓冲区&#xff0c;分配回收成本较高&#xff0c;但读写性能高&#xff0c;不受 JVM 内存回收管理。 举例 当上传一个较大文件&#xff08;200M&#xff09;…

【传知代码】VRT: 关于视频修复的模型(论文复现)

前言&#xff1a;随着数字媒体技术的普及&#xff0c;制作和传播视频内容变得日益普遍。但是&#xff0c;视频中由于多种因素&#xff0c;例如传输、存储和录制设备等&#xff0c;经常出现质量上的问题&#xff0c;如图像模糊、噪声干扰和低清晰度等。这类问题对用户的体验和观…

【Ubuntu20.04安装java-8-openjdk】

1 下载 官网下载链接&#xff1a; https://www.oracle.com/java/technologies/downloads/#java8 下载 最后一行 jdk-8u411-linux-x64.tar.gz&#xff0c;并解压&#xff1a; tar -zxvf jdk-8u411-linux-x64.tar.gz2 环境配置 1、打开~/.bashrc文件 sudo gedit ~/.bashrc2、…

【解决】:git clone项目报错fatal: fetch-pack: invalid index-pack output

象&#xff1a;之前一直使用gitee将个人学习和工作相关记录上传到个人gitee仓库&#xff0c;一直没出现过问题。直到有一天换电脑重新拉取代码发现出了问题&#xff0c;具体如下图&#xff1a; 原因分析&#xff1a; 经过查询发现主要原因是因为git clone的远程仓库的项目过大…

NVM安装及VUE创建项目的N种方式

VUE 参考官网&#xff1a;https://cli.vuejs.org/zh/guide/ 目录 NVM安装 1.卸载node.js 2.安装nvm ​编辑​ 3.配置 4.使用nvm安装node.js 5.nvm常用命令 创建VUE项目 1.使用vue init 创建vue2&#xff08;不推荐&#xff09; 2.使用vue create创建vue2和3&#xff…

Docker安装Mosquitto

在物联网项目中&#xff0c;我们经常用到MQTT协议&#xff0c;用MQTT协议做交互就需要部署一个MQTT服务&#xff0c;而mosquitto是一个常用的MQTT应用服务&#xff0c; Mosquitto是一个实现了消息推送协议MQTT v3.1的开源消息代理软件。MQTT&#xff08;Message Queuing Teleme…

knife4j案例

1.导入 <dependency><groupId>com.github.xiaoymin</groupId><artifactId>knife4j-spring-boot-starter</artifactId> </dependency>2.在配置类中加入 knife4j 相关配置并设置静态资源映射&#xff08;否则接口文档页面无法访问&#xff…

大模型LLM 结合联网搜索增强isou

参考&#xff1a; https://github.com/yokingma/search_with_ai 在线使用网址&#xff1a; https://isou.chat/ 安装github下载&#xff0c;运行docker compose 如果一直报下面错误&#xff1a; 解决方法https://github.com/yokingma/search_with_ai/pull/7 默认打开&a…

【PG数据库】PostgreSQL 日志归档详细操作流程

1.1 日志归档的目的 pg数据库日志归档是将PostgreSQL数据库的日志文件进行归档的过程。 归档的主要目的是为了保留历史数据&#xff0c;确保数据的一致性和完整性&#xff0c;同时为数据恢复提供必要的支持。 pg数据库日志归档的目的包括&#xff1a; 1.数据恢复&#xff1…

[链表专题]力扣141, 142

1. 力扣141 : 环形链表 题 : 给你一个链表的头节点 head &#xff0c;判断链表中是否有环。 如果链表中有某个节点&#xff0c;可以通过连续跟踪 next 指针再次到达&#xff0c;则链表中存在环。 为了表示给定链表中的环&#xff0c;评测系统内部使用整数 pos 来表示链表尾…

洛谷P1364 医院设置

P1364 医院设置 题目描述 设有一棵二叉树&#xff0c;如图&#xff1a; 其中&#xff0c;圈中的数字表示结点中居民的人口。圈边上数字表示结点编号&#xff0c;现在要求在某个结点上建立一个医院&#xff0c;使所有居民所走的路程之和为最小&#xff0c;同时约定&#xff0c…

20240511每日运维----聊聊nignx改配置所有的nginx改完unknow

1、改配置所有的nginx改完unknow src/core/nginx.h src/http/ngx_http_header_filter_module.c src/http/ngx_http_special_response.c src/http/v2/ngx_http_v2_filter_module.c 2、make 3、去objs里面把nginx文件替换过去sbin/nginx

从JSON数据到Pandas DataFrame:如何解析出所需字段

目录 一、引言 二、JSON数据的基本结构 三、使用Pandas从JSON数据中读取数据 四、从DataFrame中解析出所需字段 解析对象字段 解析嵌套对象字段 解析数组字段 五、案例与代码示例 六、总结 一、引言 在数据分析和处理的日常工作中&#xff0c;我们经常需要从各种…

Python悬置动刚度模拟及复数绘制

Python悬置动刚度模拟及复数绘制 1、复数绘制极坐标图2、动刚度的计算公式3、悬置动刚度的影响因素4、 AVL Excite 悬置动刚度的模拟 1、复数绘制极坐标图 # _*_ coding:UTF-8 _*_import matplotlib.pyplot as plt import numpy as np# 定义复数数组 complexNums [1.5 1.2j,…