椭圆周长所以近似计算-read.pptVIP

* */12 * 实验三 探月卫星速度计算 问题背景与数据处理 椭圆周长的近似计算 实验任务与操作 思考题与练习题 ? ? ? ? 应用数学学院：汪小平 wxiaoping325@163.com 长征三号甲运载火箭提供给卫星在近地点的速度大约为10.3（km/s）。这一速度比较低，不足以将卫星送往月球轨道,需要将速度提高到约10.9（km/s）才能使卫星奔向月球。为了达到奔月速度，中国航天工程师使用了卫星变轨调速技术············ 问题背景 --- 中国探月卫星的发射 2007年10月24日18时05分，中国第一颗探月卫星嫦娥一号在西昌卫星发射中心发射升空。卫星进入的初始轨道是周期为16小时的地球同步轨道。 卫星轨道长度计算 近地点距离 h=200 km,远地点距离 H=51000 km,地球半径 R=6378 km。卫星绕地球一周16小时。 我国探月卫星的初始轨道 周长: 椭圆: 问题1: 椭圆长半轴 a = ?椭圆短半轴 b = ? 问题2: 椭圆周长的积分计算 L = ？ 0≤ t ≤2? x =a cos t y =b sin t 近地点距离 h = 200 km, 远地点距离 H = 51000 km, 地球半径 R =6378 km 椭圆长半轴: a=(h+H+2R)/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) ? 利用三角形面积做近拟计算 近地点 地球 远地点 实验任务: 计算各轨道周长和卫星平均速度 370000 600 116 128000 600 48 71400 600 24 51000 600 16 51000 200 16 平均速度 (km/s) 轨道周长 (km) 远地距离 (km) 近地距离 (km) 轨道时长 (小时)