资源描述
实验三:探月卫星的飞行速度计算,问题背景与数据处理椭圆周长的近似计算实验任务与操作思考题与练习题,实验准备：Matlab函数,符号积分(精确)函数：intint(fun,var,a,b),int(4*x*t,x,2,sin(t)2*sin(t)2*t-8*t,数值积分(近似)函数：quadq=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=200km,远地点距离H=51000km,地球半径R=6378km.卫星绕地球一周16小时,我国探月卫星的初始轨道,周长:,椭圆:,x=acosty=bsint,0t2,问题1:椭圆长半轴a=?椭圆短半轴b=?问题2:椭圆周长的积分计算L=？,近地点距离h=200km,远地点距离H=51000km,地球半径R=6378km椭圆长半轴:a=(h+H+2R)/2,原始数据处理,椭圆半焦距:c=(Hh)/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;symse2t%定义两个符号变量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初始轨道200km51000km16小时轨道600km51000km24小时轨道600km71400km48小时轨道600km128000km地月转移轨道600km370000km,实验任务:计算各轨道周长和卫星平均速度,练习题与思考题,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./Timesb=sqrt(a.*a-c.*c);S=a.*b.*pi./Times;Vmax=2*S./(R+h),%近地距离%远地距离%长半轴%椭圆周长%运行周期%平均速度%最大速度,functionmlab31()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),holdonMdistance=51000,71000,128000,370000;h=600;c0=c;,fork=1:4H=Mdistance(k);a=(h+H+2*R)/2;c=(H-h)/2;b=sqrt(a*a-c*c);s=c-c0;ifk=1,s=-s;endx=s+a*cos(t);y=b*sin(t);plot(x,y,r,LineWidth,1)end,
展开阅读全文