实验三人造卫星轨道课件

《实验三人造卫星轨道课件》由会员分享，可在线阅读，更多相关《实验三人造卫星轨道课件（15页珍藏版）》请在装配图网上搜索。

1、1/12实验三实验三:探月卫星的飞行速度计算探月卫星的飞行速度计算 问题背景与数据处理问题背景与数据处理椭圆周长的近似计算椭圆周长的近似计算实验任务与操作实验任务与操作思考题与练习题思考题与练习题2/12实验准备：实验准备：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);3/12长征三号甲长征三号甲运载火箭

2、提供给运载火箭提供给卫星在近地点的速度大卫星在近地点的速度大约为约为10.3（km/s）。这一速度比较低，不足以将卫。这一速度比较低，不足以将卫星送往月球轨道星送往月球轨道,需要将速度提高到约需要将速度提高到约10.9（km/s）才能使卫星奔向月球。为了达到奔月速度，中国航才能使卫星奔向月球。为了达到奔月速度，中国航天工程师使用了卫星变轨调速技天工程师使用了卫星变轨调速技术术 中国探月中国探月2007年年10月月24日日18时时05分分，中国第中国第一颗探月卫星嫦娥一号在西昌卫星一颗探月卫星嫦娥一号在西昌卫星发射中心发射升空。卫星进入的初发射中心发射升空。卫星进入的初始轨道是周期为始轨道是周期

3、为1616小时的地球同步小时的地球同步轨道。轨道。4/12卫星轨道长度计算卫星轨道长度计算 近地点距离近地点距离 h=200 km,远地点距离远地点距离 H=51000 km,地球半径地球半径 R=6378 km.卫星绕地球一周卫星绕地球一周16小时小时我国探月卫星的初始轨道我国探月卫星的初始轨道周长周长:椭圆椭圆:x=a cos t y=b sin t0 t 2 2/02222cossin4 dttbtaL问题问题1:椭圆长半轴椭圆长半轴 a=?椭圆短半轴椭圆短半轴 b=?问题问题2:椭圆周长的积分计算椭圆周长的积分计算 L=？5/12近地点距离近地点距离 h=200 km,远地点距离远地点

4、距离 H=51000 km,地球半径地球半径 R=6378 km椭圆长半轴椭圆长半轴:a=(h+H+2R)/2原始数据处理原始数据处理椭圆半焦距椭圆半焦距:c=(H h)/2椭圆短半轴椭圆短半轴:22cab 椭圆偏心率椭圆偏心率:abaace22 6/12 2/02222cossin4 dttbtaLtbtatbta22222222cos)cos1(cossin tbaa2222cos)(tca222cos cos1 2222taca cos1 222tea 2/022cos14 dtteaL椭圆周长椭圆周长由于由于所以所以7/12tetete442222cos81cos211cos1 2/0

5、4422cos81cos211 dttete1283411 242ee 2/022cos14 dtteaL 2/02222cossin4 dttbtaL椭圆轨道周长椭圆轨道周长1283411 242eeaL 级数展开级数展开近似计算近似计算8/12运载火箭提供给运载火箭提供给卫星的速度到底是多少？卫星的速度到底是多少？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

6、;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.30239/12开普勒开普勒天上的立法者天上的立法者 行星运动第一定律行星运动第一定律:行星在通过太阳的平行星在通过太阳的平面内沿椭圆轨道运行，太阳位于椭圆的一面内沿椭圆轨道运行，太阳位于椭圆的一个焦点

7、上。又叫个焦点上。又叫“轨道定律轨道定律”。行星运动第二定律行星运动第二定律:行星在椭圆轨道上运行的行星速行星在椭圆轨道上运行的行星速度不是常数，而是在相等时间内，行星与太阳的联度不是常数，而是在相等时间内，行星与太阳的联线所扫过的面积相等。又叫线所扫过的面积相等。又叫“面积定律面积定律”。行星运动第三定律行星运动第三定律:太阳系内所有行星公转周期的平太阳系内所有行星公转周期的平方同行星轨道半长径的立方之比为一常数，这一定律方同行星轨道半长径的立方之比为一常数，这一定律也叫也叫“调和定律调和定律”。10/12根据根据Kepler第二定律第二定律,从地球到人造卫星的向径在相从地球到人造卫星的向径

8、在相等的时间内扫过相等的面积等的时间内扫过相等的面积,由此可知由此可知,人造卫星在人造卫星在近地点的速度达到最大近地点的速度达到最大,而在远地点的速度达到最小而在远地点的速度达到最小.Tabs/1 向径在每秒钟扫过的面积向径在每秒钟扫过的面积近地点近地点地球地球远地点远地点Vmax:=2*s1/(R+h)Vmin:=2*s1/(R+H)hlhRs)(5.01 HlHRs)(5.01 11/12轨道名称轨道名称 近地点距离近地点距离h 远地点距离远地点距离H初始轨道初始轨道 200km 51000km16小时轨道小时轨道 600km 51000km24小时轨道小时轨道 600km 71400km

9、48小时轨道小时轨道 600km 128000km地月转移轨道地月转移轨道 600km 370000km实验任务实验任务:计算各轨道周长和卫星平均速度计算各轨道周长和卫星平均速度 12/12练习题与思考题练习题与思考题1.简述开普列第一定律和第二定律的内容简述开普列第一定律和第二定律的内容;2.卫星运行平均速度与第一宇宙速度相比如何卫星运行平均速度与第一宇宙速度相比如何？3.卫星运行最大速度在变轨过程中的变化规律卫星运行最大速度在变轨过程中的变化规律？4.地月转移轨道的周期是否是地月转移轨道的周期是否是116小时？小时？5.说明面积计算公式说明面积计算公式 xdyydxS21的根据，并证明椭圆

10、面积公式的根据，并证明椭圆面积公式 abS 13/12%实验参考程序实验参考程序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)%近地距离近地距离%远

11、地距离远地距离%长半轴长半轴%椭圆周长椭圆周长%运行周期运行周期%平均速度平均速度%最大速度最大速度14/12function 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 onMdistance=51000,71000,128000,370000;h=600;c0=c;15/12for 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-0.500.511.522.533.54x 105-1-0.500.51x 105