数学建模-小行星的轨道模型.docVIP

小行星的轨道模型 问题: 天文学家要确定一颗小行星绕太阳运行的轨道。他在轨道平面内建立以太阳为原点的直角坐标系，在两坐标轴上取天文测量单位（1.4959787e11m）。在5个不同的时间对小行星作了5次观察，测的轨道上5个点的坐标数据如表： x1 x2 x3 x4 x5 X坐标 5.764 6.286 6.759 7.168 7.408 y1 y2 y3 y4 y5 Y坐标 0.648 1.202 1.823 2.526 3.360 一、问题分析与建立模型 求解椭圆一般方程的系数 由开普勒第一定律知，小行星的轨道为椭圆，现需要建立椭圆方程()以供研究。 天文学家确定小行星运动的轨迹时，他的依据是轨道上5个点的坐标数据。由椭圆轨道属于二次曲线，是一般方程。为了确定方程中的5个待定系数，需要将上述5个点的坐标代入上面的方程,得： ， 将这一包含5个未知数的线性方程组，写成矩阵的形式 = 。 求解这一线性方程组，即可得到曲线方程的系数。 （2）求解转化为标准方程后的系数 为了知道小行星轨道的一些参数，还必须将二次曲线方程化为椭圆的标准方程形式： . 由于太阳的位置是小行星的一个焦点，这时可以根据椭圆的长半轴a和短半轴b计算出小行星的近日点和远日点距离，以及椭圆周长L. 根据二次曲线理论，可得椭圆经过旋转和平移两种变换后的方程如下： . , , ,为C的特征值。 所以，椭圆长半轴： ；椭圆短半轴： ；椭圆半焦距： ；面积S=；离心率，则周长近似值为L= 。 二、计算求解 （1）利用matlab软件求解：首先由5个点的坐标数据形成线性方程组的系数矩阵 ： 求： X0=[5.764 6.286 6.759 7.168 7.408]; Y0=[0.648 1.202 1.823 2.526 3.360]; A=zeros(5);X0(1); for i=1:5 A(i,1)=X0(i)^2;A(i,2)=2*X0(i)*Y0(i);A(i,3)=Y0(i)^2;A(i,4)=2*X0(i);A(i,5)=2*Y0(i); end; A 得： A = 33.2237 7.4701 0.4199 11.5280 1.2960 39.5138 15.1115 1.4448 12.5720 2.4040 45.6841 24.6433 3.3233 13.5180 3.6460 51.3802 36.2127 6.3807 14.3360 5.0520 54.8785 49.7818 11.2896 14.8160 6.7200 又知A*X=B，求解B: B=[-1 -1 -1 -1 -1];X=A\B 得： X=[0.0508 -0.0351 0.0381 -0.2265 0.1321] 即（a1，a2，a3，a4，a5）= （0.0508，-0.0351，0.0381，-0.2265，0.1321） 然后求C、D 由 ；，用matlab计算: for i=1:5 a(i)=X(i); end C=[a(1) a(2);a(2) a(3)] det(C) 得 C = 0.0508 -0.0351 -0.0351 0.0381 = 7.0286e-004 求特征值: [v d]=eig(C,nobalance) 得 d = 0.0088 0 0 0.0801 所以特征值=0.0088 ，=0.0801 . 求D D=[a(1) a(2) a(3);a(2) a(3) a(5); a(4) a(5) 1] det(D) 得 D = 0.0508 -0.0351 0.0381 -0.0351 0.0381 0.1321 -0.2265 0.1321 1.0000 = 0.0010 于是，椭圆长半轴a=12.7158，短半轴b=4.2147，半焦距c=11.9970 。小行星近日点距离和远日点距离分别为 h=a-c=0.7188 和 H=24.7128 。 最后，计算离心率、面积、周长： c=11.9970; e=c/12.7158, S=pi*12.7158*4.2147, L=2*pi*12.7158*(1-e^2/4-3*e^4/64) 得 离心率 e = 0.9435 面积 S = 168.3683 周长 L = 59.1487 最终结果： 椭圆方程为： ； 离心率：e=0