卫星轨迹

一.问题提出

设卫星在空中运行的运动方程为:

其中是k 重力系数(k=401408km3/s)。卫星轨道采用极坐标表示,通过仿真,研究发射速度对卫星轨道的影响。实验将作出卫星在地球表面(r=6400KM ,θ=0)分别以v=8KM/s,v=10KM/s, v=12KM/s 发射时,卫星绕地球运行的轨迹。

二.问题分析

1.卫星运动方程一个二阶微分方程组,应用Matlab 的常微分方程求解命令ode45求解时,首先需要将二阶微分方程组转换成一阶微分方程组。若设 ,则有:

2.建立极坐标如上图所示,初值分别为:卫星径向初始位置,即地球半径:y(1,1)=6400;卫星初始角度位置:y(2,1)=0;卫星初始径向线速度:y(3,1)=0;卫星初始周向角速度:y(4,1)=v/6400。

3.将上述一阶微分方程及其初值带入常微分方程求解命令ode45求解,可得到一定时间间隔的卫星的径向坐标值y(1)向量;周向角度坐标值y(2)向量;径向线速度y(3)向量;周向角速度y(4)向量。

4.通过以上步骤所求得的是极坐标下的解,若需要在直角坐标系下绘制卫星的运动轨迹,还需要进行坐标变换,将径向坐标值y(1)向量;周向角度坐标值y(2)向量通过以下方程转换为直角坐标下的横纵坐标值X,Y 。

5.卫星发射速度速度的不同 将导致卫星的运动轨迹不同,实验将绘制卫星分别以v=8KM/s ,v=10KM/s ,v=12KM/s 的初速度发射的运动轨迹。

三.Matlab 程序及注释

1.主程序

v=input('请输入卫星发射速度单位Km/s :\nv='); %卫星发射速度输入。

axis([-26400 7000 -10000 42400 ]); %定制图形输出坐标范围。

%为了直观表达卫星轨迹,以下语句将绘制三维地球。

[x1,y1,z1]=sphere(15); %绘制单位球。

x1=x1*6400; y1=y1*6400;

???????-=+-=dt d dt dr r dt d dt d r r k dt r d θ

θθ2)(2

22222θ==)2(,)1(y r y ?????????????**-=**+*-===)1(/)4()3(2)4()4()4()1()1()1()3()4()2()

3()1(y y y dt dy y y y y y k dt dy y dt dy y dt dy ???*=*=)]2(sin[)1(Y )]2(cos[)1(X y y y

y

matlab给定坐标画卫星轨道,Matlab仿真实例-卫星轨迹相关推荐

  1. MATLAB STK构建卫星轨道姿态控制仿真演示系统

    MATLAB & STK构建卫星轨道姿态控制仿真演示系统 杜耀珂 郑科宇 摘要:通过在MATLAB中建立卫星的轨道运动和姿态动力学仿真模型就可以分析在不同的控制作用下,卫星的轨道和姿态的变化情 ...

  2. 卫星姿态的matlab编程,MATLABSTK构建卫星轨道姿态控制仿真演示系统.pdf

    MATLABkSTK构建卫星轨道姿态控制仿真演示系统 杜耀珂郑科宇 (上海航天控制[j翻研究所上海200233) 摘要:通过在MATLAB中建立卫星的轨道运动和姿态动力学仿真模型就可以分析在不同的控制 ...

  3. matlab 画卫星轨道,Matlab仿真实例-卫星轨迹

    卫星轨迹 一.问题提出 设卫星在空中运行的运动方程为: d2rkd 2 r() 22 dtrdt 2 d 2drd rdtdt dt2 其中是k重力系数(k=401408km3/s).卫星轨道采用极坐 ...

  4. 基于Matlab的汽车安全应用轨道融合仿真(附源码)

    目录 一.目的 二.轨道到轨道架构 三.定义方案 4.1 在仿真开始时分析跟踪 4.2 分析街道边行人的跟踪 ​ 4.3 避免谣言传播 五.总结 六.程序 此示例演示如何融合两辆车的履带,以提供比每辆 ...

  5. matlab坐标画正六边形,Matlab任意正六边形随机撒点实现

    任意正六边形随机撒点 是通信仿真里经常用到的,无奈网上资料很少,参照http://www.ilovematlab.cn/thread-136761-1-1.html,我又苦思敏想里一下午,这个问题才得 ...

  6. matlab画p轨道,matlab绘制原子轨道图.pdf

    matlab绘制原子轨道图 课程论文 结构化学 2012 年4 月 course work structural chemistry April ,2012 球谐函数,原子轨道,电子云空间分布等值线图 ...

  7. matlab画p轨道,matlab pcolor画图问题

    pcolor画图去掉网格 gca = pcolor(t_axis,h_axis,d3_4); set(gca, 'LineStyle','none'); % 去掉pcolor网格 eg: pcolor ...

  8. matlab中如何画圆锥,使用Matlab画出圆台圆锥圆柱

    用Matlab画圆柱圆台圆锥 自己建一个M文件,命名任意,如tu.m 在命令窗口输入:tu %以下为tu.m文件 clc clear all %圆柱 [X,Y,Z]=cylinder(15.7,20) ...

  9. matlab imshow加画网格,[转载]matlab中imshow函数

    函数简介 imshow是matlab中显示图像的函数.在matlab的命令窗口中输入doc imshow即可得到关于该函数的帮助信息. 调用方式 imshow(I,n) imshow(I,[low h ...

最新文章

  1. Apache 与Tomcat
  2. Tensorboard安装和访问(pytorch+MobaXterm)
  3. Draft-微软出品的云原生下的本地开发辅助工具
  4. 苹果11如何设置9宫格_4月9日三只可转债回落卖出条件单结果,10日永高转债上市,手把手教你如何设置华宝智投的涨跌幅条件单...
  5. linux vim基本操作,vim基本操作笔记
  6. 在GDI+中如何实现以左下角为原点的笛卡尔坐标系
  7. 2.6.24内核中对S3C2440的引导启动分析
  8. 联想Y9000P2021双SSD系统win10下ubantu20.04+ROS的安装(记录非教程)
  9. Collections.unmodifiableList方法的使用与场景
  10. hive 计算周几_HIVE 计算指定日期本周的第一天和最后一天
  11. mac安装win7流程 win7提示无法验证数字签名的解决
  12. gwt 同步和异步_使用GWT和RESTful Web服务构建动态的组织树
  13. leetcode 717. 1比特与2比特字符(python)
  14. 输入90输公因数java_关于Java中的几个数论中的知识
  15. ROM、PROM、EPROM、EEPROM、FLASH ROM
  16. 用webbroswer控件自动点击确定按钮
  17. 中国高校计算机大赛--网络技术挑战赛(C4-Network Technology Challenge)参加分享
  18. 【新书推荐】图神经网络导论,清华大学刘知远老师著作
  19. HADOOP集群大数据词频统计及设计比较(完整教程)
  20. 排序:quicksort,quicksort3way

热门文章

  1. Lazada知识局:lazada运营经验分享,基础知识学习
  2. Linux 从当前目录下切换到根目录下
  3. 老路用得上的商学课-21-40学习(读书)笔记
  4. springcloud H版+alibaba cloud4
  5. 卷积层、池化层等作用
  6. idea2018 2020_2020年8月21日全军武器装备采购信息网 采购信息筛选
  7. 整洁代码之道 12 迭进 Emergence
  8. ios html a标签,iOS解析html标签内容-OCGumbo
  9. 扬州大学微型计算机控制技术期末考试,计算机控制技术
  10. 手机APP-小米手机无障碍功能开启失败