用于设计和分析具有恒定近心点半径的低推力螺旋轨迹研究(Matlab代码实现)

 💥💥💞💞欢迎来到本博客❤️❤️💥💥

🏆博主优势:🌞🌞🌞博客内容尽量做到思维缜密,逻辑清晰,为了方便读者。

⛳️座右铭:行百里者,半于九十。

📋📋📋本文目录如下:🎁🎁🎁

目录

💥1 概述

📚2 运行结果

🎉3 参考文献

🌈4 Matlab代码实现


💥1 概述

本文可用于设计和分析具有恒定近心点半径的低推力螺旋轨迹。低推力轨迹设计是一种航天器轨道设计的方法,通过使用较小的推力在较长的时间内进行轨道转移,以实现目标轨道的变化。这种方法适用于需要节省燃料或对轨道精度要求较高的任务。

用户需要提供初始地球停泊轨道的经典轨道要素,这些要素包括轨道倾角、升交点赤经、近地点幅角、轨道倾角的变化率以及近地点幅角的变化率。此外,用户还需要提供最终或“目标”轨道的离心率,即目标轨道的形状。

航天器的轨道运动是使用修改的赤道坐标轨道要素进行建模。这些修改的赤道坐标轨道要素包括半长轴、离心率、倾角、升交点赤经和近地点幅角。通过使用这些轨道要素,可以准确地描述航天器在轨道上的位置和速度。

在进行轨道转移计算时,本脚本假设整个轨道转移过程中推力加速度是恒定的。这意味着航天器在整个轨道转移过程中以相同的推力加速度进行运动。这种简化假设可以简化计算过程,并得到近似的轨道解决方案。

此外,本脚本还假设轨道运动是开普勒运动或未受扰动的。这意味着在计算过程中不考虑其他天体或外力对航天器轨道的影响。这种假设适用于大多数低地球轨道任务,其中其他天体或外力的影响可以忽略不计。

可以设计和分析具有恒定近地点半径的低推力螺旋轨道,以满足特定任务的要求。这个代码提供了一个简单而有效的工具,可以帮助工程师和科学家在航天器轨道设计和分析方面取得更好的结果。

📚2 运行结果

 

 

部分代码:

function [r, v, smovrp, cosl, sinl, wmee, ssqrd] = mee2eci_lts(mu, mee)

% convert modified equinoctial orbital elements to eci position
% and velocity vectors and auxillary variables

% input

%  mu     = gravitational constant (km**3/sec**2)
%  mee(1) = semilatus rectum of orbit (kilometers)
%  mee(2) = f equinoctial element
%  mee(3) = g equinoctial element
%  mee(4) = h equinoctial element
%  mee(5) = k equinoctial element
%  mee(6) = true longitude (radians)

% output

%  r = eci position vector (kilometers)
%  v = eci velocity vector (kilometers/second)

% auxilliary output quantities

%  smovrp = sqrt(mu / pmee)
%  cosl   = cos(xlmee)
%  sinl   = sin(xlmee)
%  wmee   = 1.0 + fmee * cosl + gmee * sinl
%  ssqrd  = 1.0 + tani2s

% Orbital Mechanics with Matlab

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

% unload equinoctial orbital elements

pmee = mee(1);
fmee = mee(2);
gmee = mee(3);
hmee = mee(4);
xkmee = mee(5);
xlmee = mee(6);

smovrp = sqrt(mu / pmee);

tani2s = hmee^2 + xkmee^2;

cosl = cos(xlmee);

sinl = sin(xlmee);

wmee = 1.0 + fmee * cosl + gmee * sinl;

radius = pmee / wmee;

hsmks = hmee^2 - xkmee^2;

ssqrd = 1.0 + tani2s;

% compute eci position vector (kilometers)

🎉3 参考文献

文章中一些内容引自网络,会注明出处或引用为参考文献,难免有未尽之处,如有不妥,请随时联系删除。

[1]刘玥.火星探测器近心点制动与轨道保持优化设计[D].哈尔滨工业大学,2011.DOI:10.7666/d.D262275.

[2]夏勇俊.具有中心导体轴的螺旋槽慢波系统及注—波互作用的研究[D].电子科技大学[2023-08-27].DOI:CNKI:CDMD:2.2002.102390.

🌈4 Matlab代码实现

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

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

相关文章

【LeetCode】剑指 Offer <二刷>(3)

目录 题目:剑指 Offer 06. 从尾到头打印链表 - 力扣(LeetCode) 题目的接口: 解题思路: 代码: 过啦!!! 题目:剑指 Offer 07. 重建二叉树 - 力扣&#xf…

A 股个股资金流 API 数据接口

A 股个股资金流 API 数据接口 全量股票资金流数据,全量A股数据,最长30日历史数据 1. 产品功能 支持所有A股资金流数据查询;每日定时更新数据;支持多达 30 日历史数据查询;超高的查询效率,数据秒级返回&am…

探索树堆Treap和红黑树的优势和劣势

探索树堆Treap和红黑树的优势和劣势 一、背景知识二、树堆(Treap)的介绍三、红黑树(RB-Tree)的介绍四、树堆(Treap)与红黑树(RB-Tree)的比较总结 博主简介 💡一个热爱分享…

Elasticsearch文档多个输入字段组成ID实现方法

1、场景描述: 使用Elasticsearch时,有时会需要指定文档id的场景,当文档id需要多个字段组成时,这种业务怎么处理呢? 2、问题描述: 现有一个ElasticSearch文档,假设文档id由userid、 eventTime…

【Mysql问题集锦】:Can‘t create table ‘#sql-58d7_431d‘ (errno: 28)

问题描述: 问题原因: OSError: [Errno 28] No space left on device,即:磁盘空间不足,无法创建文件。因此,导致Mysql无法执行SQL语句。 问题解法: Step 1,查看有哪些目录占用了大量…

已解决“SyntaxError: invalid character in identifier“报错问题

本文摘要:本文已解决 Python FileNotFoundError 的相关报错问题,并总结提出了几种可用解决方案。同时结合人工智能GPT排除可能得隐患及错误。 😎 作者介绍:我是程序员洲洲,一个热爱写作的非著名程序员。CSDN全栈优质领…

vmware虚拟机(ubuntu)远程开发golang、python环境安装

目录 1. 下载vmware2. 下载ubuntu镜像3. 安装4. 做一些设置4.1 分辨率设置4.2 语言下载4.3 输入法设置4.4 时区设置 5. 直接切换管理员权限6. 网络6.1 看ip6.2 ssh 7. 本地编译器连接远程服务器7.1 创建远程部署的配置7.2 文件同步7.3 远程启动项目 8. ubuntu安装golang环境8.1…

Springboot - 5.Bean的生命周期

✍1. Bean的生命周期: 当然,我会详细描述每一步的作用。 🎷1. 实例化Bean: 这是Bean生命周期的第一步。Spring容器通过反射机制创建Bean的实例。public class ExampleBean {// ... }🎷2. 设置Bean的属性: Spring容器将根据配置…

SQL查询本年每月的数据

--一、以一行数据的形式,显示本年的12月的数据,本示例以2017年为例,根据统计日期字段判断,计算总和,查询语句如下:selectsum(case when datepart(month,统计日期)1 then 支付金额 else 0 end) as 1月, sum…

15个关于AI的Github库

这里是我们精选的创新项目列表(排名不分先后),这些项目正在机器学习和人工智能领域蓬勃发展 1:privateGPT 作者:imartinezGithub 星数:16.7K描述:利用LLM的力量,在没有互联网连接的情…

初探---Qt

目录 一、介绍Qt 二、软件安装 三、QT工具介绍 四、Assistant帮助文档的使用 五、设计师界面的介绍 ​编辑 六、QT工程项目各文件初始程序的介绍 1> 配置文件:.pro文件 2> 头文件 3> 源文件 4> 主程序 5> 各文件之间调用方式 七、第一个…

【LeetCode】双指针妙解有效三角形的个数

Problem: 611. 有效三角形的个数 文章目录 题目分析讲解算法原理复杂度Code 题目分析 首先我们来分析一下本题的思路 看到题目中给出的示例 题目的意思很简单,就是将给到的数字去做一个组合,然后看看这三条边是否可以构成三角形。那判断的方法不用我说&a…

【LeetCode75】第四十四题 省份数量

目录 题目: 示例: 分析: 代码: 题目: 示例: 分析: 给我们一个二维数组,表示城市之间的连通情况,连在一起的城市为一个省份,问我们一共有多少个省份。 这…

CNN 01(CNN简介)

一、卷积神经网络的发展 convolutional neural network 在计算机视觉领域,通常要做的就是指用机器程序替代人眼对目标图像进行识别等。那么神经网络也好还是卷积神经网络其实都是上个世纪就有的算法,只是近些年来电脑的计算能力已非当年的那种计算水平…

【爬虫】7.2. JavaScript动态渲染界面爬取-Selenium实战

JavaScript动态渲染界面爬取-Selenium实战 爬取的网页为:https://spa2.scrape.center,里面的内容都是通过Ajax渲染出来的,在分析xhr时候发现url里面有token参数,所有我们使用selenium自动化工具来爬取JavaScript渲染的界面。 fr…

Swift 如何从图片数据(Data)检测原图片类型?

功能需求 如果我们之前把图片对应的数据(Data)保持在内存或数据库中,那么怎么从 Data 对象检测出原来图片的类型呢? 如上图所示:我们将 11 张不同类型的图片转换为 Data 数据,然后从 Data 对象正确检测出了原图片类型。 目前,我们的代码可以检测出 jpeg(jpg), tiff,…

手撕 视觉slam14讲 ch7 / pose_estimation_3d2d.cpp (1)

首先理清我们需要实现什么功能,怎么实现,提供一份整体逻辑:包括主函数和功能函数 主函数逻辑: 1. 读图,两张rgb(cv::imread) 2. 找到两张rgb图中的特征点匹配对 2.1定义所需要的参数:keypoints…

【2511. 最多可以摧毁的敌人城堡数目】

来源:力扣(LeetCode) 描述: 给你一个长度为 n ,下标从 0 开始的整数数组 forts ,表示一些城堡。forts[i] 可以是 -1 ,0 或者 1 ,其中: -1 表示第 i 个位置 没有 城堡。…

VMware 安装 Centos7 超详细过程

CentOS系统,安装教程可参考以下: 哪些模型需要在Linux下运行,需提前预装Linux系统呢,评论区讨论吧 比如Noah-MP 5.0模型 1.软硬件准备 软件:推荐使用 VMware,我用的是 VMware 12 镜像:CentO…

【python爬虫】10.指挥浏览器自动工作(selenium)

文章目录 前言selenium是什么怎么用设置浏览器引擎获取数据解析与提取数据自动操作浏览器 实操运用确认目标分析过程代码实现 本关总结 前言 上一关,我们认识了cookies和session。 分别学习了它们的用法,以及区别。 还做了一个项目:带着小…