2000年C题飞越北极_的数学模型.pdf
第31卷第1期2001年1月数学的实践与认识MA THEMA T ICS I N PRACT ICE AND THEORYVol131No11Jan.2001MathematicalM odel for Flying Across the North PoleZHON G Shao2jun,LUO Feng2yin,WAN G Guo2gang(Huanggang TeachersCollege,Huanggang438000)Abstract:By carefully analysing we give a definite answer to the question 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 shortestdistance 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小时”的问题,在飞行速度恒定的条件下,转化为计算飞机航程的问题.根据题目的要求建立两个模型.在球体模型?中,利用几何知识推出飞机航程和经纬度之间的直接关系,进而算得飞行节约的时间为410504小时.在旋转椭球体模型 中,解法?利用测量学中的贝赛尔方法,给出了飞机航程的近似计算公式,算得飞行节约的时间为41041小时.解法 则构造了一个简单的弧长作为两地间的近似航程,利用积分给出了弧长的精确计算公式和近似计算公式,算得飞行节约的时间分别为410535小时和410531小时.这些结果解释了原题中“节约4小时”的估计.1问题的提出(略)2模型的假设11 不考虑地球的自转.21 飞机每经相邻两地的航程,均以曲面上两点间最短距离进行计算.31 飞机飞行中途不需降落加油,同时忽略升降时间.http:/41 开辟新航线后,飞机由北京经过北极上空直飞底特律.3数据的说明在以下计算中,北京的坐标用A0(40,116),底特律的坐标用A11(43,83),飞机原航线途中站点经纬度用表1的数据.4模型的建立与求解模型?:地球是半径为6371千米的球体.(略)模型:地球是一旋转椭球体,赤道半径为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.由上式求导,得dydx=-b2xa2y,(1)曲线在点P处的一阶导数就是该处切线的斜率,即dydx=tan(90+B)=-cotB(2)由(1)、(2)两式得:yx=b2a2tanB,由椭圆的几何作图原理知,图2中,OP=a,OP=b.于是601数学的实践与认识31卷x=O Pcosu=acosuy=O Psinu=bsinu.代入yx=b2a2tanB得:tanu=aybx=batanB.所以归化纬度与大地纬度的关系为tanu=batanB(3)(2)计算公式的导出长距离大地主题解算的贝赛尔方法中,椭球面上两点P1(L1,u1)与P2(L2,u2)间的大地线S与其辅助球面上大圆弧 的微分关系式为dS=a1-e2cos2ud.由于地球的扁率=a-ba01003,则计算中略去 2项,于是dS=a1-2cos2ud=a(1-+sin2u)d(4)球面直角极三角形的球面三角公式为sinu=sinuncos,将其代入(4)得:dS=a1-+sin2un1+cos22d(5)对(5)式两端积分得:S=a1-+2sin2un+14asin2un(sin22-sin21)(6)由于un,1,2未知,将此用u1,u2和表示,其中=2-1.由球面三角公式有sinu1=sinuncos1sinu2=sinuncos2(7)令U=(sinu1+sinu2)2,V=(sinu1-sinu2)2,可得sin2un=U2(1+cos)+V2(1-cos)(8)sin2un(sin22-sin21)=sinU1+cos-V1-cos(9)是辅助球面上两点的大圆弧长,有如下近似公式:cos=sinu1sinu2+cosu1cosu2cos(L1-L2).根据对的近似处理,(6)式可简化为如下近似公式:S=a-a2sin2un+a4sin2un(sin22-sin21)将(8)(9)代入上式得:S=a-14a-sin1+cosU+sin1-cosV(10)若令H=a4=a-b4,M=-sin1+cos,N=+sin1-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-b4,M=-sin1+cos,N=-sin1+cos(12)其中ui,ui+1可由(3)式计算得出,且飞机在地球上空10千米的椭球面上运动,a=6388千米,b=6367千米.计算得,飞机原来的总航程S=14605136(千米).飞机由北京直飞底特律的新航程S=106451562(千米)飞机飞行节约的时间t=S-SV=14605136-10645136980=41041(小时)(13)(13)式的结果解释了“北京至底特律可节约4小时”的估计.解法:在控制测量学中有较多的计算大地线长度的方法,对于那些方法,已推导得比较严密,尽管精度较高,但计算方法过于繁琐.现在,我们运用数学的方法,将大地坐标转化为空间直角坐标后,用微分法求弧长的方法来近似计算旋转椭球面上两点间最短的距离.首先引进椭球面上点的空间直角坐标:以椭球体中心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=a2a2+b2tan2b,易得平面OAiAi+1的方程为n1x+n2y-n3z=0(15)n1=ba2(sinLitanBi+1-sinLi+1tanBi),n2=ba2(cosLi+1tanBi-cosLitanBi+1),n3=-sin(Li-Li+1).将(14)式代入(15)式得:n1a2cosL+n2a2sinL=n3b2tanB故tanB=n1a2n3b2cosL+n2a2n3b2sinL=a2b2nsin(L+0)(16)其中n=n21+n22n3,0=arctann1n2.将(16)式代回(14)式得:801数学的实践与认识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)=a2a2b2n2+12sin2(L+0)+(n2+1)cos2(L+0)ab2n2sin2(L+0)+13(18)所以椭球面上两点Ai,Ai+1的弧长AiAi+1=L2L1F(L)dL(19)(19)式可作为航程Si的近似公式,但由于计算较繁,先对它作出简化处理.因a与b的值差很小,我们以a2b2n2+1替换(18)式中出现的n2+1,(18)式化为近似式F(L)=a2a2b2n2+1ab2n2sin2(L+0)+12(20)将(20)代入(19)易得航程Si的一个近似计算公式:Si=aarctanba2n2+b2cot(Li+1+0)-arctanba2n2+b2cot(Li+0)(21)n1=ba2(sinLitanBi+1-sinLi+1tanBi),n2=ba2(cosLi+1tanBi-cosLitanBi+1),n3=-sin(L1-L2),n=n21+n22n3,0=arctann1n2由以上一组公式,用计算机编程可直接代入大地坐标得到飞机原来的总航程S=14610189(千米)飞机从北京直飞底特律的新航程S=10638179(千米)所以飞行节约的时间t=S-SV=14610.89-10638.79980=4.0531(小时)(22)在上面的计算中,我们对积分式作了近似处理.如不作近似,根据(18)和(19)并用M atthematica 410软件直接算得结果为:原航程S=146041448(千米),新航程S=106311986(千米).于是飞行节约的时间t=S-SV=14604.448-10631.986980=4.0535(小时)(23)式(22)和(23)的两个结果只相差010004小时,即相当于1144秒,其误差完全可以忽略.这说明了我们上面所作的近似计算是合理的,同时,这一结果解释了原题中“节约4小时”的估计.9011期仲银花等:“飞越北极”的数学模型第31卷第1期2001年1月数学的实践与认识MA THEMA T ICS I N PRACT ICE AND THEORYVol131No11Jan.20015模型的评价(略)参考文献:1姜启源 1 数学模型(第二版)1 北京:高等教育出版社,199312同济大学数学教研室 1 高等数学(第四版)1 北京:高等教育出版社,199613武汉测绘学院控制测量教研室等 1 控制测量学(下册)1 北京:测绘出版社,199514孔祥元,梅是义 1 控制测量学 1 北京:测绘出版社,199515杨珏,梅是义 1 控制测量学 1 北京:测绘出版社,19951The mathematicalM odels of Flying over the North PoleZHON G Yin2hua,L IL i2jun,ZHAN G Q in(L ianyungang College of Chem ical Technology,L ianyungang222001)Abstract:This paper intends to explain the problem of saving flight ti me by four hours,which is converted into the problem of the flight distance in a constant flying speed.U nderthese conditions,two mathematical models are established.(1)In the spherical model,thedirect relation between the flight distance and the latitude and longitude are clarified by applyingthe know ledge of geometry.For the model of revolving ellipsoid,the Bessel Theory is appliedto work out the approxi mate formula,by which the flight ti me is calculated for saving 41041hours.(2)A si mple curve is constructed,which is adopted as the approxi mate flight distance.By the application of integration,2 formulae,approxi mate and accurate,are established so theflight ti me save 410535 and 410531 hours,respectively.The calculation results explain theproblem of flight ti me.飞越北极何永强,陆新根,沈重欢指导老师:数模组(浙江万里学院,宁波315101)编者按:若将地球视作旋转椭球,飞机的航线应为长、短半轴分别为6388和6367千米的椭圆旋转而得的旋转椭球面上过给定两点的短程线即测地线.本文应用微分几何知识,给出了测地线满足的微分方程,并借助数学软件求得多数航线段的长度.该方法有一定的特点,是可取的.我们选取了论文的这一部分内容,予以发表.若使计算更加精确,应将地理纬度转化为归化纬度(详见本期 飞越北极的数学计算模型 一文).摘要:本文对“飞机从北京出发、飞越北极直达底特律的所需时间,可比原航线节省多少时间”的问题进行讨论,并将航线选择归结为寻求曲面上的最短弧.应用“曲面上最短弧为测地线”的事实进行了讨论.模型(一)假设地球是球体,我们可通过单位向量的
