实验三人造卫星轨道课件

15,/12,实验三,:,探月卫星的飞行速度计算,问题背景与数据处理,椭圆周长的近似计算,实验任务与操作,思考题与练习题,实验准备：Matlab函数,符号积分(精确)函数：int,int(fun,var,a,b),int(4*x*t,x,2,sin(t),2*sin(t)2*t-8*t,数值积分(近似)函数：quad,q = quad(fun,a,b),Examples:,F = inline(1./(x.3-2*x-5);,Q = quad(F,0,2);,长征三号甲,运载火箭提供给,卫星在近地点的速度大约为,10.3（km/s）,。这一速度比较低，不足以将卫星送往月球轨道,需要将速度提高到约,10.9（km/s）,才能使卫星奔向月球。为了达到奔月速度，中国航天工程师使用了卫星变轨调速技术,中国探月,2007,年,10,月,24,日,18,时,05,分,，,中国第一颗探月卫星嫦娥一号在西昌卫星发射中心发射升空。卫星进入的初始轨道是周期为16小时的地球同步轨道。,卫星轨道长度计算,近地点距离,h,=200 km,远地点距离,H,=,51000,km,地球半径,R,=6378 km.,卫星绕地球一周,16,小时,我国探月卫星的初始轨道,周长:,椭圆,:,x,=,a,cos,t,y,=,b,sin,t,0,t,2,问题,1:,椭圆长半轴,a = ?,椭圆短半轴,b = ?,问题,2:,椭圆周长的积分计算,L = ？,近地点距离,h =,200 km,远地点距离,H,=,51000,km,地球半径,R,=6378 km,椭圆长半轴,:,a=,(,h+H,+2,R,)/2,原始数据处理,椭圆半焦距,:,c=,(,H h,)/2,椭圆短半轴,:,椭圆偏心率,:,椭圆周长,由于,所以,椭圆轨道周长,级数展开,近似计算,运载火箭提供给,卫星的速度到底是多少？,h=200;H=51000;R=6378;,%,输入重要数据,a=(h+H+2*R)/2;c=(H-h)/2;,%,计算长半轴和半焦距,e1=c/a; b=sqrt(a*a-c*c);,%,计算离心率和短半轴,T=16*3600;S1=a*b*pi/T;,syms e2 t,%,定义两个符号变量,f=sqrt(1-e2*cos(t)2);,%,定义符号表达式,ft=subs(f,e2,e1*e1);,%,替换离心率数据,S=int(ft,0,pi/2);,%,计算积分,L=4*a*double(S),%,符号数据转换为数值,V=L/T,%,计算平均速度,Vmax=2*S1/(R+h),Vmax = 10.3023,开普勒天上的立法者,行星运动第一定律:,行星在通过太阳的平面内沿椭圆轨道运行，太阳位于椭圆的一个焦点上。又叫“轨道定律”。,行星运动第二定律:,行星在椭圆轨道上运行的行星速度不是常数，而是在相等时间内，行星与太阳的联线所扫过的面积相等。又叫“面积定律”。,行星运动第三定律:,太阳系内所有行星公转周期的平方同行星轨道半长径的立方之比为一常数，这一定律也叫“调和定律”。,根据,Kepler,第二定律,从地球到人造卫星的向径在相等的时间内扫过相等的面积,由此可知,人造卫星在近地点的速度达到最大,而在远地点的速度达到最小,.,向径在每秒钟扫过的面积,近地点,地球,远地点,Vmax:=2*s1/(R + h),Vmin:=2*s1/(R + H),轨道名称 近地点距离,h,远地点距离,H,初始轨道,200km 51000km,16,小时轨道,600km 51000km,24,小时轨道,600km,71400,km,48,小时轨道,600km,128000,km,地月转移轨道,600km,370000,km,实验任务,:,计算各轨道周长和卫星平均速度,练习题与思考题,1.,简述开普列第一定律和第二定律的内容,;,2.,卫星运行平均速度与第一宇宙速度相比如何,？,3.,卫星运行最大速度在变轨过程中的变化规律,？,4.,地月转移轨道的周期是否是116小时？,5. 说明面积计算公式,的根据，并证明椭圆面积公式,%,实验参考程序,R=6378;,h=200,600,600,600,600;,H=51000,51000,71000,128000,370000;,a=(h+H+2*R)/2;,c=(H-h)/2;E2=(c./a).2;,L=2*pi*a.*(1-E2/4-3*E2.2/64),Time=16,15.63,23.3,50.5,225*3600;,Times=16,16,24,48,232*3600;,V0=L./Times,b=sqrt(a.*a-c.*c);,S=a.*b.*pi./Times;,Vmax=2*S./(R+h),%近地距离,%远地距离,%长半轴,%椭圆周长,%运行周期,%平均速度,%最大速度,function mlab31(),R=6378;,h=200;H=51000;,a=(h+H+2*R)/2;,c=(H-h)/2;,b=sqrt(a*a-c*c);,t=(0:.05:2)*pi;,x=a*cos(t);y=b*sin(t);,plot(x,y,k,-c,0,ob),hold on,Mdistance=51000,71000,128000,370000;,h=600;c0=c;,for k=1:4,H=Mdistance(k);,a=(h+H+2*R)/2;,c=(H-h)/2; b=sqrt(a*a-c*c);,s=c-c0;,if k=1,s=-s;end,x=s+a*cos(t);y=b*sin(t);,plot(x,y,r,LineWidth,1),end,