《第二章-开普勒方程PPT课件》由会员分享,可在线阅读,更多相关《第二章-开普勒方程PPT课件(43页珍藏版)》请在金锄头文库上搜索。
1、1,第二章 开普勒方程,主讲教师:杏建军 2020年9月28日,2,授课内容,卫星轨迹的预测 开普勒方程的求解 开普勒轨道根数,3,1.1 轨道预报,已知t0时刻卫星的位置r0,速度v0,预报时刻t,卫星的位置r和速度v 在极坐标下,如何计算的?,4,1.1 轨道预报,5,1.1 轨道预报,带入位置r和速度v的表达式中,回到原始问题,已知t0,r0,v0,求时刻t的r和速度v,6,1.1 轨道预报,需要,开普勒第二定律,主要求解出上述方程,就可以得到真近点角f与时间t之间的关系,进而预报卫星在轨道的位置和速度。,7,授课内容,卫星轨迹的预测 开普勒方程的求解 开普勒轨道根数,8,2.1 开普勒
2、方程,当0e1时,引入一个新的角度E(偏近点角),在以C点(椭圆中心)为原点的坐标系中,9,2.1 开普勒方程,在以F点(椭圆一个焦点)为原点的坐标系时,椭圆参数方程为,坐标系原点平移到C点后,椭圆参数方程为,得到r与E 的数学关系,10,2.1 开普勒方程,进一步,可以得到,11,2.1 开普勒方程,上述两式相除,再开平方,得到,上式两边同时对时间求导数,12,2.1 开普勒方程,由开普勒第二定理,积分,得到,著名的开普勒方程,表示了时间与真近点角的函数关系,其中 是一个新的积分常数,13,2.1 开普勒方程,定义平近点角M:,通过开普勒方程,可以得到运动时间 t 与偏近点角E,根据E与真近
3、点角 f 的关系,得到 f ,进而进行卫星轨道的预报。,现在的问题是如何求解开普勒方程?,14,2.2 开普勒方程的求解,第一种方法:序列迭代法,迭代停止条件,15,2.2 开普勒方程的求解,迭代格式是否收敛?,课堂练习:编制开普勒迭代求解matlab程序,16,2.2 开普勒方程的求解,M = pi/2; E0 = 0; E1 = 1; tem =0; e = 0.2; n = 0; while(abs(E1-E0)1e-6) E0 = tem; E1 = M + e*sin(E0); tem = E1; n = n+1; end n,当n=6时,E收敛,17,2.2 开普勒方程的求解,第二
4、种方法:拉格朗日方法(1770年),考虑函数,为一个小参数,如椭圆偏心率,18,2.2 开普勒方程的求解,19,2.2 开普勒方程的求解,开普勒方程的拉格朗日级数解,E = M +e*sin(M)+e2/2*2*cos(M)*sin(M) + e3/6*(6*cos(M)2*sin(M)- 3*sin(M)3) + e4/24*(24*cos(M)3*sin(M) - 40*cos(M)*sin(M)3)+ e5/120*(65*sin(M)5 + 120*cos(M)4*sin(M) - 440*cos(M)2*sin(M)3),20,2.2 开普勒方程的求解,拉格朗日系数求解问题,n =5
5、; syms M,n; ff = diff(sin(M)n,n-1),拉格朗日级数收敛问题,21,2.2 开普勒方程的求解,第三种方法:牛顿迭代法,给定方程解的一个初始猜测值x0,迭代停止条件,22,2.2 开普勒方程的求解,课堂练习题:推导开普勒方程的牛顿迭代格式,并编制Matlab程序,比较与第一种迭代格式的区别,第一种迭代格式,23,2.2 开普勒方程的求解,M = pi/2; E0 = 0; E1 = 1; tem =0; e = 0.8; n = 0; while(abs(E1-E0)1e-6) E0 = tem; E1 = M + e*sin(E0); % E1 = E0 - (E
6、0-M-e*sin(E0)/(1-e*cos(E0); tem = E1; n = n+1; end n,e=0.8时,第一种迭代格式需要20次,牛顿迭代法需要6次,24,2.2 开普勒方程的求解,其他方法: 傅里叶级数法 高斯法,e=1,抛物线;e1,双曲线有类似的结果,25,2.2 开普勒方程的求解,课后作业 应用三种方法求解开普勒方程,计算M = 260,e = 0.8时的偏近点角E和真近点角f,并对结果进行对比分析。,26,2.3 卫星轨道的预报,已知t0,r0,v0,求时刻t的r和速度v,27,2.3 卫星轨道的预报,课后作业 以x0 = -5292392.072;-4862.201
7、380;3111662.355; -4136.781314; 3101114.660;-4147.028008人造地球轨道卫星的初始条件,应用本节的知识,预测2天后卫星的位置和速度,并与Matlab的结果进行比较。,28,授课内容,卫星轨迹的预测 开普勒方程的求解 开普勒轨道根数,29,3.1 开普勒轨道根数,已经得到的二体运动方程积分常数 轨道角动量矢量h,偏心率矢量e,卫星过近地点时刻,30,3.1 开普勒轨道根数,31,3.1 开普勒轨道根数,32,3.2 开普勒轨道根数与星下点轨迹,星下点轨迹:航天器质心与地心连线与地球表面的交点,航天器轨道周期,33,3.2 开普勒轨道根数与星下点轨
8、迹,星下点轨迹与轨道半长轴 地球同步轨道,34,3.2 开普勒轨道根数与星下点轨迹,星下点轨迹与轨道倾角,星下点轨迹的最高纬度就是航天器轨道倾角,35,3.2 开普勒轨道根数与星下点轨迹,星下点轨迹与轨道偏心率,36,3.3 开普勒轨道根数的计算,已知位置矢量和速度矢量,求6个轨道根数,37,3.3 开普勒轨道根数的计算,已知6个轨道根数,求卫星位置矢量和速度矢量,38,3.3 开普勒轨道根数的计算,39,3.3 开普勒轨道根数的计算,第一章作业 function robit_computer mu = 3.986004418e+14; %地球引力常数 tspan = 0:60:86400;
9、options = odeset(AbsTol,1e-15,RelTol,1e-12,NormControl,on); x0 = -5292392.072;-4862.201380;3111662.355; -4136.781314; 3101114.660;-4147.028008; T,Y = ode45(orbit,tspan,x0,options,mu);,课后作业:将第一章作业中计算的每一个时刻的位置速度转换6个轨道根数,然后再转换为位置速度,并于原来结果进行比较,画出二者的误差图。,40,本章小结,航天器轨道的预报,41,本章小结,开普勒方程,42,本章小结,开普勒轨道根数(6个),43,请批评指正!,