卫星轨迹

一.问题提出

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

其中是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. echart饼状图没有数据的时候显示暂无数据_Python数据结构可视化 day 5
  2. MySQL数据库如何解决大数据量存储问题
  3. c++不允许使用不完整的类型
  4. python 多层嵌套 json_python解析多层json操作示例
  5. sql 注入神器sqlmap 源码分析之调试sqlmap
  6. div+css布局必了解的列表元素ul ol li dl dt dd详解
  7. 加菲猫的人生歪理~ 看完果然开心,哈哈~
  8. c语言 不显示dos窗口,运行DOS批处理不显示DOS窗口的方法
  9. 刀片服务器更换纽扣电池系统起不来,莆田四龄童将纽扣电池塞入鼻孔 医生建议手术遭拒...
  10. 彩虹易支付源码v8.1.1-支持多通道轮训/API退款/实名认证/完善手动提现t0t1
  11. 完美卸载SQL Server 2008的方案
  12. SQL注入语法类型和sql盲注
  13. 爱加密加固病毒分析-破解篇
  14. 软件测试结果分析和质量报告
  15. 文明重启战局服务器维护中,王牌战争文明重启8月23日更新公告
  16. vue run dev报错 缺少package.json文件、missing dev 命令 解决办法
  17. 海贼王83名刀:无上大快刀12工、大快刀21工、良快刀50工
  18. 1500个用户代理User agent,用于随机UA代理
  19. Django 重写authenticate实现输入账号、邮箱、手机号登录验证
  20. 三次函数求近似解(牛顿迭代法/二分法)

热门文章

  1. 手机鸿蒙系统APP,手机鸿蒙操作系统
  2. CSS中图标字体的引用及3种使用方式——阿里巴巴矢量库
  3. 分治法基本思想(汉诺塔问题 Tower of Hanoi)
  4. Java设计模式—适配器模式(adapter pattern)
  5. OpenOCD编译安装
  6. mysql 创建索引、删除索引、查看索引sql语句
  7. OpenSSL密码库算法笔记——第5.4.1章 椭圆曲线上点的仿射坐标表示
  8. 时代中国CIO张智:ITSM护航数字化,这是更好的时代!
  9. 目标检测制作数据集时画框的原则
  10. github使用AccessToken登录配置教程