数学建模课程设计报告--飞越北极.docVIP

数学建模课程设计报告 设计题目： 飞越北极 小组编号： 07 组 班 级： 信息与计算科学 组 员： 摘 要 本文针对“中美航线飞越北极，底特律至北京的飞行时间可节省4个小时”的问题，在飞行速度恒定，时间已给出的条件下，将其转化为计算飞机航程的问题。并以飞机飞行航线所经地点的地理坐标为依据，建立球面数学模型，通过数学方程式将地理坐标转换为空间坐标以及MATLAB数据循环分析等方法获取模型参数，得到相应数据后可计算最优弧和位速度位移公式。 对于模型一，假设地球是球体，由数学球面模型绘出的飞机飞行航线，采用MATLAB 编程以及立体几何知识计算球面距离得出从北京直接到底特律的时间为10.18734h ，而按飞机的原航线则至少需14.17793h ，所以至少节省时间为3.9059h。对于模型二，假设地球是旋转椭球体，采用参数方程得出纬度与经度之间的函数关系，然后再用积分方法求得球面距离，最后得出节省时间为3.917h。 通过以上方法，对飞行航线和球体的分析可证明“飞越北极，可节省时间为4 小时”的命题。 关键词：球面数学模型；飞行时间；最优弧 一、问题重述、模型假设、 图2-大地纬度、归化纬度和地心纬度之间的关系 五、法线的斜率为tan f，它与切线斜率之间的关系为: 或者可以写成: 若已知大地纬度 f，则归化纬度为: α是地心纬度，则： 若已知大地纬度 f，则地心纬度： 旋转椭球面上两点，过这两点以及地心的平面与椭球面的交线是一条曲线，由于地球的长半轴6378千米与短半轴6367千米相比差别很小，可以近似地认为飞机沿着椭球面两点之间的“大圆弧”飞行时所需时间最短，于是可以用定积分的办法来近似计算。 旋转椭球面的直角坐标方程为： 它可以表示成以归化纬度θ和经度φ为参数的参数方程形式： 式中φ是经度（西经转换为东经），θ是归化纬度。 我们认为题目所给纬度是大地纬度 f，先把它转换成归化纬度θ，再计算旋转椭球面上两点间最短航程。 假设A1(θ1,φ1),A2(θ2,φ2)是北半球上的两点，过地心O和点A1,A2的平面与椭球面的交线是A1到A2的最短航线。 设A1,A2的直角坐标分别为(x1,y1,z1)和(x2,y2,z2)，则过地心O和点A1,A2 三个点的平面方程为： 它可以写成Ax+By+Cz=0的形式，其中:A=y1z2-y2z1，B=x2z1-x1z2，C=x1y2-x2y1， 将参数方程代入并化简得： 上式是平面方程的参数形式，从它可以求出θ与φ的关系式 以及导数： 参数方程对φ求导数得： 根据三维空间中弧长的定积分计算公式，两个点A1(θ1,φ1),A2(θ2,φ2)之间的弧长为： 两点之间的最短航程近似计算的表达式为： 用MATLAB软件计算得到的个点之间的最短航程和飞行时间如下表： i 起点 终点 最短航程（千米） 飞行时间（小时） 0 A0 A1 1117.4680 1.1403 1 A1 A2 1761.6744 1.7976 2 A2 A3 4634.2581 4.7288 3 A3 A4 1342.2460 1.3696 4 A4 A5 643.1888 0.6563 5 A5 A6 539.7118 0.5507 6 A6 A7 652.4363 0.6658 7 A7 A8 498.3833 0.5086 8 A8 A9 228.4989 0.2332 9 A9 A10 2818.1137 2.8756 10 A10 A11 331.8756 0.3386 合计 14567.8549 14.8652 飞越北极 A0 A11 10729.1774 10.9481 飞越北极可以节省飞行时间Δt=3.917小时。 六、 七、模型推广 上述模型的建立与求解依据立体几何、参数方程和定积分的原理和方法，并用MATLAB 计算，方法科学，符合实际，结果基本精确合理。模型1中，由于地球是一球体，所以优点显而易见，计算步骤少，方法相对简单。模型2 的计算过程就繁琐复杂得一些，但模型2 却比模型1 更切合实际。由于在两模型中，均不考虑地球自转等多种因素对飞机的影响，所以结果会产生一定误差，建模中引进了参数方程的定积分求解比较精确地求出了时间T。 参考文献 [1] 包晔等，“飞越北极”时间节省模型，第13卷第4期， 2001年。 [2] 姜启源等，数学模型，北京，高等教育出版社，2011年。 [