《_飞越北极_的数学模型》由会员分享，可在线阅读，更多相关《_飞越北极_的数学模型（6页珍藏版）》请在金锄头文库上搜索。
1、第31卷第1期 2001年1月数学的实践与认识 MA THEMA T ICS I N PRACT ICE AND THEORYVol131 No11 Jan.2001 MathematicalM odel for Flying Across the North PoleZHON G Shao2jun, LUO Feng2yin, WAN G Guo2gang(Huanggang TeachersCollege, Huanggang 438000)Abstract: By carefully analysing we give a definite answer to the question
2、posed in Yang ZiN ight N ew s that the ti me for flying from Beijing to Detroit w ill be four hours shorter if theplane can fly across the North Pole. A reasonably mathematicalmadel is framed in both casesthat the earth is a sphere or an ellipsoid. W e also give an algorithm for finding the shortest
3、distance between two points in a surface(especially in an ellipsoid).Two approxi matealgorithm s, constriction ratio method and i m itating search method, are added.“飞越北极” 的数学模型仲银花, 李利军, 张 琴 指导老师: 数模教练组(连云港化工高等专科学校,连云港 222001)编者按: 本文对视地球为椭球的情况采用了测量学中的贝塞尔方法和用积分计算弧长的方法,针对地球扁率很小的特点对积分求弧长给出了一个方便的近似公式,使计
4、算大为简化,这是本文的一个特点.在处理实际问题时,在满足一定精度的前提下,合理地作简化近似可以用较少的计算达到预期的目的.本文的另一个特点是对大地纬度和归化纬度的关系交待很清楚,反映了同学查阅背景材料是出色的.摘要: 本文将 “飞行时间节约4小时” 的问题,在飞行速度恒定的条件下,转化为计算飞机航程的问题.根据题目的要求建立两个模型.在球体模型? 中,利用几何知识推出飞机航程和经纬度之间的直接关系,进而算得飞行节约的时间为410504小时.在旋转椭球体模型 中,解法? 利用测量学中的贝赛尔方法,给出了飞机航程的近似计算公式,算得飞行节约的时间为41041小时.解法 则构造了一个简单的弧长作为两
5、地间的近似航程,利用积分给出了弧长的精确计算公式和近似计算公式,算得飞行节约的时间分别为410535小时和410531小时.这些结果解释了原题中 “节约4小时” 的估计.1 问题的提出(略)2 模型的假设11 不考虑地球的自转.21 飞机每经相邻两地的航程,均以曲面上两点间最短距离进行计算.31 飞机飞行中途不需降落加油,同时忽略升降时间.41 开辟新航线后,飞机由北京经过北极上空直飞底特律.3 数据的说明在以下计算中,北京的坐标用A0(40, 116),底特律的坐标用A11(43, 83),飞机原航 线途中站点经纬度用表1的数据.4 模型的建立与求解模型?:地球是半径为6371千米的球体.(
6、略) 模型:地球是一旋转椭球体,赤道半径为6378千米,子午线短半轴为6357千米.几个基本概念: (1)大地经度:如图1,P点的子午面N PS与N GS起始子午面所构成的二面角L.(2)大地纬度:如图1,P点的法线Pn与赤道面的夹角B.(3)归化纬度:如图2,把点P上的y线向上延伸,与以a为半径的大圆弧A N 相交于P,则OP 与X轴夹角叫归化纬度. (4)大地线:椭球面上两点间的最短曲线.图1图2解法?: 在控制测量学中,常使用归化纬度计算两已知点之间的大地线长度.(1)归化纬度与大地纬度的转化如图2,P点所在子午椭圆的方程式为x2a2+y2b2= 1.由上式求导,得dy dx= -b2x
7、 a2y,(1)曲线在点P处的一阶导数就是该处切线的斜率,即dy dx= tan (90+B) = -cotB(2)由(1)、(2)两式得:y x=b2 a2tanB,由椭圆的几何作图原理知,图2中,OP=a,OP=b.于是601数 学 的 实 践 与 认 识31卷x=O Pcosu=acosuy=O Psinu=bsinu.代入y x=b2 a2tanB得: tanu=ay bx=b atanB.所以归化纬度与大地纬度的关系为tanu=b atanB(3)(2)计算公式的导出长距离大地主题解算的贝赛尔方法中,椭球面上两点P1(L1,u1)与P2(L2,u2)间的大地线S与其辅助球面上大圆弧
8、的微分关系式为dS=a1-e2cos2ud.由于地球的扁率 =a-b a01003,则计算中略去 2项,于是dS=a1 -2cos2ud=a(1 -+sin2u)d(4)球面直角极三角形的球面三角公式为sinu= sinuncos,将其代入(4)得:dS=a1 -+sin2un1 + cos2 2d(5)对(5)式两端积分得:S=a1 -+ 2sin2un+1 4asin2un(sin22-sin21)(6)由于un,1,2未知,将此用u1,u2和表示,其中=2-1.由球面三角公式有sinu1= sinuncos1sinu2= sinuncos2(7)令U= (sinu1+ sinu2)2,V
9、= (sinu1- sinu2)2,可得sin2un=U 2(1 + cos)+V 2(1 -cos)(8)sin2un(sin22-sin21) = sinU 1 + cos-V 1 -cos(9)是辅助球面上两点的大圆弧长,有如下近似公式:cos= sinu1sinu2+ cosu1cosu2cos(L1-L2).根据对的近似处理, (6)式可简化为如下近似公式:S=a-a 2sin2un+a 4sin2un(sin22-sin21)将(8) (9)代入上式得:S=a-1 4a-sin 1 + cosU+ sin 1 -cosV(10)若令H=a 4=a-b 4,M=- sin 1+ co
10、s,N=+ sin 1- cos,则S=a-H(M U-N V)(11)对于飞机航线相邻两站点Ai(Bi,Li),Ai+ 1(Bi+ 1,Li+ 1)之间的航程有如下近似计算公式:Si= -a-H(M U+N V),cos= sinuisinui+ 1+ cosuicosui+ 1cos(Li+ 1-cosLi),U=(sinui+ sinui+ 1)2,V=(sinui-sinui+ 1)2,7011期仲银花等:“飞越北极” 的数学模型H=a-b 4,M=-sin 1 + cos,N=-sin 1 + cos(12)其中ui,ui+ 1可由(3)式计算得出,且飞机在地球上空10千米的椭球面上
11、运动,a= 6388千 米,b= 6367千米.计算得,飞机原来的总航程S= 14605136(千米).飞机由北京直飞底特律的新航程S= 106451562(千米) 飞机飞行节约的时间t=S-S V=14605136 -10645136 980= 41041(小时)(13)(13)式的结果解释了 “北京至底特律可节约4小时” 的估计.解法: 在控制测量学中有较多的计算大地线长度的方法,对于那些方法,已推导得比较严密, 尽管精度较高,但计算方法过于繁琐.现在,我们运用数学的方法,将大地坐标转化为空间 直角坐标后,用微分法求弧长的方法来近似计算旋转椭球面上两点间最短的距离. 首先引进椭球面上点的空
12、间直角坐标:以椭球体中心O为原点,起始子午面与赤道交 线为X轴,在赤道面上与X轴正交的方向为Y轴,椭球体的旋转轴为Z轴,构成空间直角 坐标系,则(x,y,z)为点的空间直角坐标. 根据大地坐标及空间直角坐标的定义,可以得x=a2a2+b2tan2BcosLy=a2a2+b2tan2BsinLz=b2tanBa2+b2tan2B(14)设飞机航线上相邻两站点为Ai(Bi,Li),Ai+ 1(Bi+ 1,Li+ 1),显然其处于同一椭球体表 面.利用(14)式计算出空间坐标.我们用过点Ai,Ai+ 1和原点的平面截椭球表面所得弧长AiAi+ 1近似作为Ai,Ai+ 1两地间航程Si.方便起见,令D
13、=a2a2+b2tan2b,易得平面OAiAi+ 1的方程为n1x+n2y-n3z= 0(15)n1=b a2(sinLitanBi+ 1-sinLi+ 1tanBi),n2=b a2(cosLi+ 1tanBi-cosLitanBi+ 1),n3= -sin (Li-Li+ 1).将(14)式代入(15)式得:n1a2cosL+n2a2sinL=n3b2tanB 故tanB=n1a2 n3b2cosL+n2a2 n3b2sinL=a2 b2nsin (L+0)(16)其中n=n21+n22 n3,0= arctann1 n2.将(16)式代回(14)式得:801数 学 的 实 践 与 认 识
14、31卷x=aba2n2sin2(L+0) +b2cosLy=aba2n2sin2(L+0) +b2sinLz=abnsin (L+0)a2n2sin2(L+0) +b2(17)于是,F(L)=x2(L) +y2(L) +z2(L)=a2a2 b2n2+ 12 sin2(L+0) +(n2+ 1)cos2(L+0)a b2 n2sin2(L+0) + 13(18)所以椭球面上两点Ai,Ai+ 1的弧长AiAi+ 1=L2L1F(L) dL(19)(19)式可作为航程Si的近似公式,但由于计算较繁,先对它作出简化处理.因a与b的值差很小,我们以a2b2n2+ 1替换(18)式中出现的n2+ 1,
15、(18)式化为近似式F(L) =a2a2 b2n2+ 1a b2 n2sin2(L+0) + 12(20)将(20)代入(19)易得航程Si的一个近似计算公式:Si=aarctanba2n2+b2cot(Li+ 1+0)-arctanba2n2+b2cot(Li+0)(21)n1=b a2(sinLitanBi+ 1-sinLi+ 1tanBi),n2=b a2(cosLi+ 1tanBi-cosLitanBi+ 1),n3= -sin (L1-L2),n=n21+n22 n3,0= arctann1 n2由以上一组公式,用计算机编程可直接代入大地坐标得到 飞机原来的总航程S= 1461018
16、9(千米) 飞机从北京直飞底特律的新航程S= 10638179(千米) 所以飞行节约的时间t=S-S V=14610. 89 -10638. 79 980= 4. 0531(小时)(22)在上面的计算中,我们对积分式作了近似处理.如不作近似,根据(18)和(19)并用M atthematica 410软件直接算得结果为:原航程S= 146041448 (千米),新航程S=106311986(千米).于是飞行节约的时间t=S-S V=14604. 448 -10631. 986 980= 4. 0535(小时)(23)式(22)和(23)的两个结果只相差010004小时,即相当于1144秒,其误差完全可以忽略.这说 明了我们上面所作的近似计算是合理的,同时,这一结果解释了原题中 “节约4小时” 的估计.
