TH方程学习(5)

2024-06-04 00:52
文章标签 学习 方程 th

本文主要是介绍TH方程学习(5),希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

一、内容介绍

在本节中,将会利用TH方程的状态转移矩阵,求解双脉冲轨道转移问题:给定转移时间T,通过状态转移矩阵求解初始速度增量\Delta \boldsymbol{v}_{1}和末端速度\Delta \boldsymbol{v}_{2},使得追踪星在初始速度增量\Delta \boldsymbol{v}_{1}的作用下转移到目标位置\boldsymbol{r_{f}},追踪星到达\boldsymbol{r_{f}}后,在\Delta \boldsymbol{v}_{2}的作用下达到期望的相对速度\boldsymbol{v_{f}}

在第四节中,在给定目标星初始偏心率,近地点高度,真近地点角以及转移时间的情况下,可以求出其状态转移矩阵\Phi_{}。根据状态转移矩阵的性质,可以写为如下形式

\begin{bmatrix} \boldsymbol{r}_f\\\boldsymbol{v}_f \end{bmatrix}=\begin{bmatrix} \boldsymbol{\Phi}_{rr} & \boldsymbol{\Phi}_{rv}\\\boldsymbol{\Phi}_{vr} & \boldsymbol{\Phi}_{vv} \end{bmatrix}\begin{bmatrix} \boldsymbol{r}_0\\\boldsymbol{v}_0 \end{bmatrix}

因此,对于第一次需要的脉冲,可以通过

\boldsymbol{r}_{f}=\boldsymbol{\Phi}_{rr}\boldsymbol{r}_0+\boldsymbol{\Phi}_{rv}\boldsymbol{v}_0^{+}\\ \boldsymbol{v}_0^{+}=\boldsymbol{\Phi}_{rv}^{-1}(\boldsymbol{r}_f-\boldsymbol{\Phi}_{rr}\boldsymbol{r}_0)\\ \Delta \boldsymbol{v}_{1}=\boldsymbol{v}_0^{+}-\boldsymbol{v}_0

对于第二次需要施加的脉冲,可以得到

\boldsymbol{v}_f^{-}=\boldsymbol{\Phi}_{vr}\boldsymbol{r}_0+\boldsymbol{\Phi}_{vv}\boldsymbol{v}_0^{+}\\\Delta \boldsymbol{v}_2=\boldsymbol{v}_f-\boldsymbol{v}_f^{-}

二、案例仿真

本文使用STK进行了一次案例仿真,代码如下:TH_solver函数在第四节给出

% 使用STK验证VVLH坐标系
clc;clear
uiApplication = actxGetRunningServer('STK12.application');
root = uiApplication.Personality2;
checkempty = root.Children.Count;
if checkempty ~= 0root.CurrentScenario.Unloadroot.CloseScenario;
end
root.NewScenario('VVLH');
StartTime = '26 Jan 2024 04:00:00.000';    % 场景开始时间
StopTime = '10 Feb 2024 04:00:00.000';     % 场景结束时间
root.ExecuteCommand(['SetAnalysisTimePeriod * "',StartTime,'" "',StopTime,'"']);
root.ExecuteCommand(' Animate * Reset');
SatName = 'Target';       %  SAR_   GX_   Sat_  GX_1_  SAR_1_
satellite = root.CurrentScenario.Children.New('eSatellite', SatName);
satellite.SetPropagatorType('ePropagatorAstrogator');   %  不设置的时候默认为二体模型  ePropagatorJ4Perturbation
satellite.Propagator;
% 目标星初始状态
Perigee = 500;
T       = 60;
% 追踪星在VVLH坐下的相对位置
delta_r = [0.1;0.01;0.01];
delta_v = [0.0001;0.0001;0.0001];
Perige  = 6378.137+Perigee;
ecc     = 0.1;
sma     = Perige/(1-ecc);
Inc     = 30;
w       = 0;
RAAN    = 0;
TA      = 45;
root.ExecuteCommand(['Astrogator */Satellite/',SatName,' SetValue MainSequence.SegmentList Initial_State Propagate']);
InitialState=satellite.Propagator.MainSequence.Item(0);
%% 初始化卫星参数
root.ExecuteCommand(['Astrogator */Satellite/',SatName,' SetValue MainSequence.SegmentList.Initial_State.CoordinateType Modified Keplerian']);
root.ExecuteCommand(['Astrogator */Satellite/',SatName,' SetValue MainSequence.SegmentList.Initial_State.InitialState.Epoch ',StartTime,' UTCG']);
root.ExecuteCommand(['Astrogator */Satellite/',SatName,' SetValue MainSequence.SegmentList.Initial_State.InitialState.Keplerian.sma ',num2str(sma),' km']);
root.ExecuteCommand(['Astrogator */Satellite/',SatName,' SetValue MainSequence.SegmentList.Initial_State.InitialState.Keplerian.ecc ',num2str(ecc)]);
root.ExecuteCommand(['Astrogator */Satellite/',SatName,' SetValue MainSequence.SegmentList.Initial_State.InitialState.Keplerian.inc ',num2str(Inc),' deg']);
root.ExecuteCommand(['Astrogator */Satellite/',SatName,' SetValue MainSequence.SegmentList.Initial_State.InitialState.Keplerian.w ',num2str(w),' deg']);
root.ExecuteCommand(['Astrogator */Satellite/',SatName,' SetValue MainSequence.SegmentList.Initial_State.InitialState.Keplerian.RAAN ',num2str(RAAN),' deg']);
root.ExecuteCommand(['Astrogator */Satellite/',SatName,' SetValue MainSequence.SegmentList.Initial_State.InitialState.Keplerian.TA ',num2str(TA),' deg']);
%% 二体传播
Propagate=satellite.Propagator.MainSequence.Item(1);
Propagate.PropagatorName='Earth Point Mass';
root.ExecuteCommand(['Astrogator */Satellite/',SatName,' RunMCS']);% 插入目标星
SatName2 = 'Chaser';       
satellite2 = root.CurrentScenario.Children.New('eSatellite', SatName2);
satellite2.SetPropagatorType('ePropagatorAstrogator');   %  不设置的时候默认为二体模型  ePropagatorJ4Perturbation
satellite2.Propagator;
InitialState2=satellite2.Propagator.MainSequence.Item(0);
InitialState2.CoordSystemName='Satellite/Target VVLH';
InitialState2.Element.X=delta_r(1);
InitialState2.Element.Y=delta_r(2);
InitialState2.Element.Z=delta_r(3);
InitialState2.Element.Vx=delta_v(1);
InitialState2.Element.Vy=delta_v(2);
InitialState2.Element.Vz=delta_v(3);
Propagate2=satellite2.Propagator.MainSequence.Item(1);
Propagate2.PropagatorName='Earth Point Mass';
satellite2.Propagator.MainSequence.Insert('eVASegmentTypeManeuver','Maneuver1','Propagate');
Maneuver=satellite2.Propagator.MainSequence.Item(1);
root.ExecuteCommand(['Astrogator */Satellite/',SatName2,' SetValue MainSequence.SegmentList.Maneuver1.ImpulsiveMnvr.AttitudeControl Thrust Vector']);
Maneuver.Maneuver.AttitudeControl.ThrustAxesName='Satellite VVLH.Axes';satellite2.Propagator.MainSequence.Insert('eVASegmentTypeManeuver','Maneuver2','Propagate');
satellite2.Propagator.MainSequence.Insert('eVASegmentTypePropagate','Propagate1','Maneuver2');
Propagate3=satellite2.Propagator.MainSequence.Item(2);
Propagate3.PropagatorName='Earth Point Mass';
Propagate3.Properties.Color=255;
Propagate3.StoppingConditions.Item(0).Properties.Trip = T;
root.ExecuteCommand(['Astrogator */Satellite/',SatName2,' SetValue MainSequence.SegmentList.Maneuver2.ImpulsiveMnvr.AttitudeControl Thrust Vector']);
Maneuver2=satellite2.Propagator.MainSequence.Item(3);
Maneuver2.Maneuver.AttitudeControl.ThrustAxesName='Satellite VVLH.Axes';
%% 本代码旨在使用TH方程求解双脉冲转移问题% 期望到达的位置
r_f     =    zeros(3,1);
v_f     =    zeros(3,1);
% TH求出状态转移矩阵
[v,Phi,vv]=TH_solver(ecc,Perigee,TA,delta_r,delta_v,T);
% 求出状态转移矩阵的每个部分
Phi_rr    =  Phi(1:3,1:3);
Phi_rv    =  Phi(1:3,4:6);
Phi_vr    =  Phi(4:6,1:3);
Phi_vv    =  Phi(4:6,4:6);
% 求出第一次脉冲施加的大小
vv_0      =  inv(Phi_rv)*(r_f-Phi_rr*delta_r);
dv1       =  vv_0-delta_v;
% 求出第二次施加脉冲的大小
vv_f      =  Phi_vr*delta_r+Phi_vv*vv_0;
dv2       =  v_f-vv_f;
Maneuver.Maneuver.AttitudeControl.X=dv1(1)*1000;
Maneuver.Maneuver.AttitudeControl.Y=dv1(2)*1000;
Maneuver.Maneuver.AttitudeControl.Z=dv1(3)*1000;
Maneuver2.Maneuver.AttitudeControl.X=dv2(1)*1000;
Maneuver2.Maneuver.AttitudeControl.Y=dv2(2)*1000;
Maneuver2.Maneuver.AttitudeControl.Z=dv2(3)*1000;
root.ExecuteCommand(['Astrogator */Satellite/',SatName2,' RunMCS']);
% 报告二颗卫星的三维关系
satellite.VO.OrbitSystems.InertialByWindow.IsVisible=0;
satellite2.VO.OrbitSystems.InertialByWindow.IsVisible=0;
satellite2.VO.OrbitSystems.Add('Satellite/Target VVLH System')
satellite.VO.Vector.RefCrdns.Item(2).Visible=1;targetdata=root.ExecuteCommand(['Report_RM */Satellite/Target  Style "VVLH" TimePeriod "26 Jan 2024 04:00:00.000" "26 Jan 2024 16:00:00.000" TimeStep 10']);
Num=targetdata.Count;
root.ExecuteCommand('Astrogator */Satellite/Target ClearDWCGraphics');
root.ExecuteCommand('Astrogator */Satellite/Chaser ClearDWCGraphics');
for j=1:Num-2struct=regexp(targetdata.Item(j),',','split');Tar_x(j)=str2double(struct{2});Tar_y(j)=str2double(struct{3});Tar_z(j)=str2double(struct{4});
endfigure(1)
plot3(Tar_x(1:12),Tar_y(1:12),Tar_z(1:12),'LineWidth',1);
axis([-0.1 0.1 -0.1 0.1 -0.1 0.1])
set(gca,'XDir','reverse');
set(gca,'YDir','reverse');
set(gca,'ZDir','reverse');
xlabel('X axis(km)','FontName','Times New Roman')
ylabel('Y axis(km)','FontName','Times New Roman')
zlabel('Z axis(km)','FontName','Times New Roman')
title('e=0.1,Perigee=500km','FontName','Times New Roman')
grid on
hold on
plot3(0,0,0,'g.')
plot3(delta_r(1),delta_r(2),delta_r(3),'r.')
legend('transfer trajectory(TH)','Original','Final','Location','Northeast')

这篇关于TH方程学习(5)的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/1028598

相关文章

HarmonyOS学习(七)——UI(五)常用布局总结

自适应布局 1.1、线性布局(LinearLayout) 通过线性容器Row和Column实现线性布局。Column容器内的子组件按照垂直方向排列,Row组件中的子组件按照水平方向排列。 属性说明space通过space参数设置主轴上子组件的间距,达到各子组件在排列上的等间距效果alignItems设置子组件在交叉轴上的对齐方式,且在各类尺寸屏幕上表现一致,其中交叉轴为垂直时,取值为Vert

Ilya-AI分享的他在OpenAI学习到的15个提示工程技巧

Ilya(不是本人,claude AI)在社交媒体上分享了他在OpenAI学习到的15个Prompt撰写技巧。 以下是详细的内容: 提示精确化:在编写提示时,力求表达清晰准确。清楚地阐述任务需求和概念定义至关重要。例:不用"分析文本",而用"判断这段话的情感倾向:积极、消极还是中性"。 快速迭代:善于快速连续调整提示。熟练的提示工程师能够灵活地进行多轮优化。例:从"总结文章"到"用

【前端学习】AntV G6-08 深入图形与图形分组、自定义节点、节点动画(下)

【课程链接】 AntV G6:深入图形与图形分组、自定义节点、节点动画(下)_哔哩哔哩_bilibili 本章十吾老师讲解了一个复杂的自定义节点中,应该怎样去计算和绘制图形,如何给一个图形制作不间断的动画,以及在鼠标事件之后产生动画。(有点难,需要好好理解) <!DOCTYPE html><html><head><meta charset="UTF-8"><title>06

学习hash总结

2014/1/29/   最近刚开始学hash,名字很陌生,但是hash的思想却很熟悉,以前早就做过此类的题,但是不知道这就是hash思想而已,说白了hash就是一个映射,往往灵活利用数组的下标来实现算法,hash的作用:1、判重;2、统计次数;

零基础学习Redis(10) -- zset类型命令使用

zset是有序集合,内部除了存储元素外,还会存储一个score,存储在zset中的元素会按照score的大小升序排列,不同元素的score可以重复,score相同的元素会按照元素的字典序排列。 1. zset常用命令 1.1 zadd  zadd key [NX | XX] [GT | LT]   [CH] [INCR] score member [score member ...]

【机器学习】高斯过程的基本概念和应用领域以及在python中的实例

引言 高斯过程(Gaussian Process,简称GP)是一种概率模型,用于描述一组随机变量的联合概率分布,其中任何一个有限维度的子集都具有高斯分布 文章目录 引言一、高斯过程1.1 基本定义1.1.1 随机过程1.1.2 高斯分布 1.2 高斯过程的特性1.2.1 联合高斯性1.2.2 均值函数1.2.3 协方差函数(或核函数) 1.3 核函数1.4 高斯过程回归(Gauss

【生成模型系列(初级)】嵌入(Embedding)方程——自然语言处理的数学灵魂【通俗理解】

【通俗理解】嵌入(Embedding)方程——自然语言处理的数学灵魂 关键词提炼 #嵌入方程 #自然语言处理 #词向量 #机器学习 #神经网络 #向量空间模型 #Siri #Google翻译 #AlexNet 第一节:嵌入方程的类比与核心概念【尽可能通俗】 嵌入方程可以被看作是自然语言处理中的“翻译机”,它将文本中的单词或短语转换成计算机能够理解的数学形式,即向量。 正如翻译机将一种语言

【学习笔记】 陈强-机器学习-Python-Ch15 人工神经网络(1)sklearn

系列文章目录 监督学习:参数方法 【学习笔记】 陈强-机器学习-Python-Ch4 线性回归 【学习笔记】 陈强-机器学习-Python-Ch5 逻辑回归 【课后题练习】 陈强-机器学习-Python-Ch5 逻辑回归(SAheart.csv) 【学习笔记】 陈强-机器学习-Python-Ch6 多项逻辑回归 【学习笔记 及 课后题练习】 陈强-机器学习-Python-Ch7 判别分析 【学

系统架构师考试学习笔记第三篇——架构设计高级知识(20)通信系统架构设计理论与实践

本章知识考点:         第20课时主要学习通信系统架构设计的理论和工作中的实践。根据新版考试大纲,本课时知识点会涉及案例分析题(25分),而在历年考试中,案例题对该部分内容的考查并不多,虽在综合知识选择题目中经常考查,但分值也不高。本课时内容侧重于对知识点的记忆和理解,按照以往的出题规律,通信系统架构设计基础知识点多来源于教材内的基础网络设备、网络架构和教材外最新时事热点技术。本课时知识

线性代数|机器学习-P36在图中找聚类

文章目录 1. 常见图结构2. 谱聚类 感觉后面几节课的内容跨越太大,需要补充太多的知识点,教授讲得内容跨越较大,一般一节课的内容是书本上的一章节内容,所以看视频比较吃力,需要先预习课本内容后才能够很好的理解教授讲解的知识点。 1. 常见图结构 假设我们有如下图结构: Adjacency Matrix:行和列表示的是节点的位置,A[i,j]表示的第 i 个节点和第 j 个