实验三MATLAB探月卫星速度ppt课件

1、1/18实验三实验三: :探月卫星速度计算探月卫星速度计算 问题背景与数据处理问题背景与数据处理椭圆周长的近似计算椭圆周长的近似计算实验任务与操作实验任务与操作思考题与练习题思考题与练习题2/18 长征三号甲运载火箭提供给卫星在近地点的速长征三号甲运载火箭提供给卫星在近地点的速度大约为度大约为10.310.3km/skm/s）。这一速度不足以将卫星送）。这一速度不足以将卫星送往月球轨道往月球轨道, ,需要将速度提高到约需要将速度提高到约10.910.9km/skm/s才才能使卫星奔向月球。为了达到奔月速度，中国航天能使卫星奔向月球。为了达到奔月速度，中国航天工程师使用了卫星变轨调速技工程师使用

2、了卫星变轨调速技术术 中国探月工程中国探月工程 2019年年10月月24日日18时时05分，中分，中国第一颗探月卫星嫦娥一号在西昌国第一颗探月卫星嫦娥一号在西昌卫星发射中心发射升空。卫星进入卫星发射中心发射升空。卫星进入的初始轨道是周期为的初始轨道是周期为16小时的地球小时的地球同步轨道。同步轨道。3/18卫星轨道长度计算卫星轨道长度计算 近地点距离近地点距离 h=200 km,远地点距离远地点距离 H=51000 km,地球地球半径半径 R=6378 km.卫星绕地球一周卫星绕地球一周16小时小时我国探月卫星的初始轨道我国探月卫星的初始轨道周长周长:椭圆椭圆:x =a cos t y =b

3、sin t0 t 2 2/02222cossin4 dttbtaL问题问题1: 1: 椭圆长半轴椭圆长半轴 a = ? a = ?椭圆短半轴椭圆短半轴 b = ? b = ?问题问题2: 2: 椭圆周长椭圆周长 L = L = ？4/18近地点距离近地点距离 h = 200 km,远地点距离远地点距离 H = 51000 km,地球半径地球半径 R =6378 km椭圆长半轴椭圆长半轴: a=(h+H+2R)/2原始数据处理原始数据处理椭圆半焦距椭圆半焦距: c=( H h )/2椭圆短半轴椭圆短半轴:22cab 椭圆偏心率椭圆偏心率:abaace22 5/18 2/02222cossin4

4、dttbtaLtbtatbta22222222cos)cos1(cossin tbaa2222cos)( tca222cos cos1 2222taca cos1 222tea 2/022cos14 dtteaL椭圆周长椭圆周长所以所以643411 242eeaL 近似计算近似计算6/18运载火箭提供给卫星的速度到底是多少？运载火箭提供给卫星的速度到底是多少？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); %计算离心率和短半轴计算离心率

5、和短半轴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.30237/18开普列开普列天上的立法者天上的立法者 行星运动第一定律行星运动第一定律: : 行星在通过太阳的平行星在通过太阳的平面内沿

6、椭圆轨道运行，太阳位于椭圆的一面内沿椭圆轨道运行，太阳位于椭圆的一个焦点上。又叫个焦点上。又叫“轨道定律轨道定律”。行星运动第二定律行星运动第二定律: :行星在椭圆轨道上运行的行星速行星在椭圆轨道上运行的行星速度不是常数，而是在相等时间内，行星与太阳的联度不是常数，而是在相等时间内，行星与太阳的联线所扫过的面积相等。又叫线所扫过的面积相等。又叫“面积定律面积定律”。行星运动第三定律行星运动第三定律 : :太阳系内所有行星公转周期的平太阳系内所有行星公转周期的平方同行星轨道半长径的立方之比为一常数，这一定律方同行星轨道半长径的立方之比为一常数，这一定律也叫也叫“调和定律调和定律”。8/18根据根

7、据Kepler第二定律第二定律,从地球到卫星的向径在单位时从地球到卫星的向径在单位时间内扫过相等面积间内扫过相等面积.所以所以,卫星在近地点速度达到最卫星在近地点速度达到最大大.设轨道周期为设轨道周期为T,每秒钟扫过的面积每秒钟扫过的面积TabS/1 hlhRS)(211 近地点近地点地球地球远地点远地点Hh)(2maxhRTabV hlV max9/18轨道名称轨道名称 近地点距离近地点距离h h 远地点距远地点距离离H H初始轨道初始轨道 200km 200km 51000km51000km1616小时轨道小时轨道 600km 600km 51000km51000km2424小时轨道小时轨

8、道 600km 600km 71400km71400km4848小时轨道小时轨道 600km 600km 128000km128000km116116小时轨道小时轨道 600km 600km 370000km370000km实验任务实验任务: : 计算轨道周长和卫星平均速度最大速度计算轨道周长和卫星平均速度最大速度 10/18%实验参考程序实验参考程序IR=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.

9、2/64)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)Bar(Vmax)% %近地距离近地距离% %远地距离远地距离% %长半轴长半轴% %椭圆周长椭圆周长% %运行周期运行周期% %平均速度平均速度% %最大速度最大速度11/18123459101112Vmax=10.30, 10.06, 10.21, 11.27, 10.59 )(2maxhRTabV )(2maxRhVabT 123459101112设计各轨道最大速度为设计各轨道最大速度为:10.30 1

10、0.45 10.60 10.75 10.90 各轨道周期各轨道周期 = ? = ? ？ ？ ？ ？设各轨道周期设各轨道周期: 16,16,24,48,23212/18R=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;b=sqrt(a.*a-c.*c);Vmax=linspace(10.3,10.9,5);bar(Vmax,c)S=a.*b.*pi;Times=2*S./Vmax./(R+h)/3600Times = 16.00 15.41 23.1

11、3 50.36 225.4113/18练习题与思考题练习题与思考题1.简述开普列第一定律和第二定律的内容简述开普列第一定律和第二定律的内容;2.变轨过程中几个轨道的周期是否是准确的？变轨过程中几个轨道的周期是否是准确的？ 3.卫星运行速度与第一、第二宇宙速度相比如何？卫星运行速度与第一、第二宇宙速度相比如何？4.卫星运行最大速度在变轨过程中的变化规律？卫星运行最大速度在变轨过程中的变化规律？5.地月转移轨道的周期是否是地月转移轨道的周期是否是116小时？小时？6. 说明面积计算公式说明面积计算公式 xdyydxS21的根据，并证明椭圆面积公式的根据，并证明椭圆面积公式 abS 14/18附注

12、NoteIn what follows, is the area the satellite swipe, is the element of displacement, ( ) is the distance betweenthe satellite and the surface of earth.1, ( )1( )21( )2sdlh tdsdlSv tdtdtdsRh t dldsRh tdt)2 1( )( )2 1max2 1mindldtSv tRh tSVRhSVRH15/18附注 Note22201()2Stokes formula says1()( )2cos( ),si

13、n( )11()(sin ( )cos ( )221211()22ydxxdyddxdydydxxdydxdySxatybtydxxdyabtt dtabdtydxxdyabdtab 16/18function satel1() R=6378;h=200,600,600,600,600;H=51000,51000,71000,128000,370000;a=(h+H+2*R)/2;c=(H-h)/2;b=sqrt(a.*a-c.*c);t=linspace(-pi,pi,500);x=a*cos(t);y=b*sin(t);s=c-c(1);s(2)=-s(2); %从近地轨道到从近地轨道到1

14、6小时轨道为小时轨道为远地点变焦，而其他均为近地点变焦远地点变焦，而其他均为近地点变焦S=s*ones(size(t);x=S+x;X=x;Y=y;figure(1),plot(X,Y,k,-c(1),0,ob)X=reshape(X,2500,1);Y=reshape(Y,2500,1);figure(2),comet(X,Y)17/18function orbit()R=6378;h=200,600,600,600,600;H=51000,51000,71000,128000,370000;a=(h+H+2*R)/2;c=(H-h)/2;b=sqrt(a.*a-c.*c);t=linspace(-pi,pi,5000);t1=linspace(-pi,0,5000);t2=linspace(0,pi,5000);x=a*cos(t);y=b*sin(t);x(1,:)=a(1)*cos(t1);y(1,:)=b(1)*sin(t1);x(2,:)=a(2)*cos(t2);y(2,:)=b(2)*sin(t2);x(5,:)=a(5)*cos(t1);y(5,:)=b(5)*sin(t1);s=c