卫星轨迹

一.问题提出

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

d2rkd 2 r() 22 dtrdt 2 d 2drd

rdtdt dt2

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

二.问题分析

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

dy(1) dt y(3)

dy(2) y(4) dt dy(3)k y(1) y(4) y(4)dty(1) y(1) dy(4) 2 y(3) y(4)/y(1) dt

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。

X y(1) cos[y

(2)] Y y(1) sin[y(2)]

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;

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

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

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

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

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

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

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

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

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

  5. matlab给定坐标画卫星轨道,Matlab仿真实例-卫星轨迹

    卫星轨迹 一.问题提出 设卫星在空中运行的运动方程为: 其中是k 重力系数(k=401408km3/s).卫星轨道采用极坐标表示,通过仿真,研究发射速度对卫星轨道的影响.实验将作出卫星在地球表面(r= ...

  6. 怎样用MATLAB画二次函数曲线,matlab画二次函数图像

    [8 70 118 100 9 0 5]; 以上是每一个 X 和 Y 对应的坐标,请问如何编程能够绘制平滑曲线,这个图形就像二次函数一样的 如果要在图中绘制一条直线加上 y=...... MATLAB ...

  7. 用matlab画干涉条纹图,matlab干涉条纹处理

    基于matlab的等倾干涉实 验仿真黄强煌 2017.12.25 目录 Cont... [5] 符运良. MATLAB 在光学教学中的应用 . 华南热带农业大学学报, [J] 2004, (3 : 1 ...

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

    <使用Matlab画出圆台圆锥圆柱>由会员分享,可在线阅读,更多相关<使用Matlab画出圆台圆锥圆柱(5页珍藏版)>请在人人文库网上搜索. 1.用Matlab画圆柱圆台圆锥自 ...

  9. matlab 画非线性曲线,MATLAB实例:非线性曲线拟合

    MATLAB实例:非线性曲线拟合 用最小二乘法拟合非线性曲线,给出两种方法:(1)指定非线性函数,(2)用傅里叶函数拟合曲线 1. MATLAB程序 clear clc xdata=[0.1732;0 ...

最新文章

  1. ncbi查找目的基因序列_使用NCBI设计qPCR引物方法
  2. 补psp进度(11月4号-9号)
  3. iPhone13下周三发布,提前看剧透:刘海缩小、120Hz高刷屏、Mini又续一年…
  4. crosstool-ng配置中的Architecture level、Emit assembly for CPU 、Tune for CPU
  5. centos minimal Bind 主从服务器部署
  6. java上传文件以流方式判断类型
  7. 35.NO module named zlib
  8. HSRP冗余热备份协议
  9. 审稿较快的2区计算机sci期刊,哪些SCI期刊的审稿速度快
  10. 计算机缺少什么无法正常启动,电脑无法正常启动
  11. JS中attr和prop区别
  12. C#里的进制与ASCII转换
  13. [hadoop3.x]HDFS中的内存存储支持(七)概述
  14. 【Hive】总体标准差stddev()
  15. AM335x启动流程(BootRom-MLO-Uboot)
  16. 若依前后端分离版本集成Mybatis-plus
  17. Python 趋势:当今最热门语言的热门话题
  18. csbte路点机器人_《cs1.6》awp地图
  19. c++多态(已知测试main函数,设计符合要求的Money类。 类Money: 数据成员:yuan,jiao,fen 成员函数:operator )
  20. 模糊控制应用(一)模糊PID在一个三阶系统上的应用

热门文章

  1. DDR4、GDDR5、GDDR6内存的区别
  2. iOS内容含有隐藏功能、第三方支付等被拒4次的心路历程
  3. 东邪、西毒、南帝、北丐、中神通五大操作系统之华山论剑(独评)
  4. 采用两块半球体U235金属的原子弹
  5. 爬取QQ音乐中一首歌的相关信息及评论(破解反爬虫、多协程队列爬虫)
  6. H5 Vue 视频 video 支持预览图 poster
  7. linux 时钟分频,浅析AD9522时钟分频电路原理
  8. vrchat模型保存_轻松简单自己上传VRChat的Avatar
  9. 编程练习:既是完全平方数又有两位数字相同的三位数
  10. TN.STN液晶屏常见问题及解决办法