实验三探月卫星速度

《实验三探月卫星速度》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、实验三:探月卫星速度计算问题背景与数据处理椭圆周长的近似计算实验任务与操作思考题与练习题1/15中国探月工程2007年10月24日18时05分，中国第一颗探月卫星嫦娥一号在西昌卫星发射中心发射升空。卫星进入的初始轨道是周期为16小时的地球同步轨道。长征三号甲运载火箭提供给卫星在近地点的速度大约为10.3（km/s）。这一速度不足以将卫星送往月球,需要提速到10.9（km/s）才能使卫星奔向月球。为了达到奔月速度，中国航天工程师使用了卫星变轨调速技术············2/153/15轨道名称近地点距离h远地点距离H(16小时)初始轨道

2、200km51000km16小时轨道600km51000km24小时轨道600km71400km48小时轨道600km128000km116小时（地月转600km370000km移）轨道4/15卫星轨道长度计算我国探月卫星的初始轨道近地点距离h=200km,远地点距离H=51000km,地球半径R=6378km.卫星绕地球一周16小时x=acost椭圆:0≤t≤2πy=bsint/22222周长:L4asintbcostdt0问题1:椭圆长半轴a=?椭圆短半轴b=?问题2:椭圆周长L=?5/15原始数据处理近地点距离h=200km,远地点

3、距离H=51000km,地球半径R=6378km椭圆长半轴:椭圆半焦距:ahHR(2)/2cHh()/2椭圆短半轴:椭圆偏心率:22cabbac22eaa6/15π/22222椭圆周长L=4∫asint+bcostdt022222222asint+bcost=a(1−cost)+bcost2222222=a−(a−b)cost=a−ccost22c2222=a[1−cost]=a[1−ecost]2aπ/222所以L=4a∫1−ecostdt01234近似计算L≈2aπ[1−e−e]4647/15taylor(f,n,x)—n

4、-1次麦克劳林多项式展开例3.21求椭圆积分近似表达式/2221ecostdt0symste2xf=sqrt(1-e2*x)F=taylor(f,3,x)g=subs(F,x,cos(t)^2)int(g,0,pi/2)ans=1/2*pi-1/8*e2*pi-3/128*e2^2*pi/2131e22costdt1e2e4024648/15运载火箭提供给卫星的速度到底是多少？h=200;H=51000;R=6378;%输入重要数据a=(h+H+2*R)/2;c=(H-h)/2;%计算长半轴和半焦距e1=

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

6、又叫“轨道定律”。行星运动第二定律:行星在椭圆轨道上运行的行星速度不是常数，而是在相等时间内，行星与太阳的联线所扫过的面积相等。又叫“面积定律”。行星运动第三定律:太阳系内所有行星公转周期的平方同行星轨道半长径的立方之比为一常数，这一定律也叫“调和定律”。10/15根据开普列第二定律,从地球到卫星的向径在单位时间内扫过相等面积.所以,卫星在近地点速度达到最大,在远地点速度达到最小.设轨道周期为T向径每秒钟扫过的面积SabT11S1()RhlhVmax:=2*S1/(R+h)21S()RHlVmin:=2*S1/(R+H)1H2Hh远

7、地点近地点地球三角形底边做为速度近似值11/15实验任务:计算轨道周长和卫星平均速度最大速度轨道名称近地点距离h远地点距离H初始轨道200km51000km16小时轨道600km51000km24小时轨道600km71400km48小时轨道600km128000km116小时轨道600km370000km12/15%实验参考程序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.

8、/a).^2;L=2*pi*a.*(1-E2/4-3*E2.^2/64)%椭圆周长Times=[16,16,24,48,116]*360