飞越北极.pdf

第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千米的椭圆旋转而得的旋转椭球面上过给定两点的短程线即测地线.本文应用微分几何知识,给出了测地线满足的微分方程,并借助数学软件求得多数航线段的长度.该方法有一定的特点,是可取的.我们选取了论文的这一部分内容,予以发表.若使计算更加精确,应将地理纬度转化为归化纬度(详见本期 飞越北极的数学计算模型 一文).摘要:本文对“飞机从北京出发、飞越北极直达底特律的所需时间,可比原航线节省多少时间”的问题进行讨论,并将航线选择归结为寻求曲面上的最短弧.应用“曲面上最短弧为测地线”的事实进行了讨论.模型(一)假设地球是球体,我们可通过单位向量的点乘与夹角的关系,加以解决;对于模型(二)设地球是旋转椭球体,我们利用微分几何学中测地线方程加以解决,并且把球面的纬度转化为旋转椭球面纬度.对于4组较特殊的点,纬度几乎相等或相近,或者两者之间的经度差过大时,用测地线计算比较困难,我们用椭圆弧(长)代替测地线长,结合数学软件M athematica的数值积分功能,可求得测地线长.1问题的重述(略)2模型一(略)3模型二311假设(1)地球为旋转椭球体;(2)为简化模型,设赤道半径+飞行高度为单位长度1;(3)所给纬度为理想球面上的纬度;(4)从北京到底特律中途不需加油.312参数说明a:为子午线半径与飞行高度之和跟赤道半径与飞行高度之和的比值(a=6367?6388);:为经度数值,单位为度;:为纬度数值,单位为度;L:为弧长,单位为公里.313模型解答首先,飞机的航线一定是曲面上两点(站)之间的最短曲线,那么航线的问题可以归结为求测地线及其弧长的问题.根据1中定理2114,“在曲面上两点间的最短曲线C必为测地线”,用参数曲面方程来表示旋转椭球体可表示为详见2 P169:r=(Cos,Sin,f()其中=1?1+tg2?a2,(因为tg=f(u)?u=a1-u2?u,因为都是北纬,tg取正数值).因为椭圆上的标准式为:u21+f2(u)a2=1,f2(u)=(1-u2)a2,故:f(u)=au1-u2,于是,曲面的第一基本形式为:I=1+f2(u)du2+u2d2详见书2P169例题或该书第三章第三节其中E=1+f2(u);F=0(正交);G=u2.若在曲面上有一条曲线C,它的方程是关于参数t的函数:u=u(t);=(t);(t0tt1)则它的弧长公式就为:1111期何永强等:飞越北极L=t1t0E(u(t),(t)dudt2+G(u(t),(t)ddt212dt令u=t=(u),则ddt=ddu;L=uu01+f2(u)+u2ddu212du=uu01+(au)21-u2+u2ddu212du(1)根据2P170公式(10),测地线满足下面方程:ddu=c1+(f(u)2uu2-c2=c1+(au)21-u2uu2-c2(c:为常数)(2)把(2)代入(1)式得Lm in=uu01+(au)21-u2+c21+(au)21-u2u2-c212du(3)现在,我们先来确定常数c.由(2)式得=c1+uu0c1+(f(u)2uu2-c2du取u=u0,c1=0,所以=0+uu0c1+(f(u)2uu2-c2du=0+uu0c1+(au)21-u2uu2-c2du其中u=1?1+tg2?a2(是纬度).特别当u=u1时,=1,于是有方程1=0+u1u0c1+(au)21-u2uu2-c2du因为,各点的 、的值是确定.于是,可以通过M athematica数学软件中的FindRoot命令输入式,可得以下数据:常数C单位椭球时的弧长L弧长L(单位:公里)北京A10.371057L1=0.174480L1=1114.58A1A20.786575L2=0.275349L2=1758.93A2A30.600374L3=0.723465L3=4621.49A3A40.350192L4=0.209363L4=1337.41A4A50.417523L5=0.100230L5=640.27A5A60.304690L6=0.084212L6=537.946A6A70.314597L7=0.101895L7=650.903A7A80.489572L8=0.0778326L8=497.195因为其余四组:A8A9,A9A10,A10底特律,北京底特律之间的纬度相等或相近,或者两者之间的经度差过大,经过计算验证,上述方法对此四组数据无法处理,所以在计算此四组数据时,由于测地线计算的困难,我们考虑经过两点、以球心为中心的椭圆,并以两211数学的实践与认识31卷点间的椭圆弧(长)代替测地线(长).椭圆弧长设球面上两点(向量)A=(x,y,z),B=(x3,y3,z3),则连接它们的大圆圆弧的方程为r(t)=tA+(1-t)BtA+(1-t)B(0t1)即r(t)=tx+(1-t)x3tA+(1-t)B,ty+(1-t)y3tA+(1-t)B,tz+(1-t)z3tA+(1-t)B曲线C的弧长公式为L=10 x(t)2+y(t)2+z(t)2下面的映射将球面变成旋转椭球面,(x,y,z)(x,y,az).如图,球面上点的纬度 7 和变换后旋转椭球面上点的纬度 满足 7=arctg(tg?a),而经度不变.由球面曲线的弧长公式,可得旋转椭球的弧长公式为L=10 x(t)2+y(t)2+(az(t)2dt利用M athematica软件的数值积分功能,可得单位弧长L弧长L(单位:公里)A8A9L9=0.03564422474L9=227.695A9A10L10=0.4397685785L10=2809.24A10 底特律L11=0.05091055547L11=325.217北京 底特律L=1.669532606L=10665.0综上所述,可得从北京飞经十站到底特律的距离为:S=611i=1Li=1452019(公里),因此可得节约时间为:31935小时.314模型二的总结由计算可知,模型二所得结果比模型一更为合理,所节约时间更多,更符合实际情况.如考虑中途停靠加油,所节约时间将不止此数.参考文献:1方德值.微分几何基础,科学出版社,1984.2陈维桓.微分几何初步.北京大学出版社,1999.3陈省身,陈维桓.微分几何讲义.北京大学出版社,1999.4M artin M,杨正清、李世杰、黄锦能译.微分几何的理论和习题.上海科学技术出版社,1989.Fly over the Arctic CirellHE Yong2qiang,LU Xin2gen,SHEN Chong2huan(ZhejiangW anliU niversity,N ingbo315101)3111期何永强等:飞越北极第31卷第1期2001年1月数学的实践与认识MA THEMA T ICS I N PRACT ICE AND THEORYVol131No11Jan.2001Abstract:A s to the problem of the ti me needed for an airplane to start from Beijing,fly overarctic pole,and reach Detroit,this article discusses how much ti me can be saved in the modelsthat established in the article in comparison w ith the original flight route.A nd it summarizesselection of the flight route for searching the shortest arc in the surface.D iscussion is based onthe fact that the shortest arc on surface is geodesic.M odel 1 is on the assumption that theEarth is a sphere.It can be solved by the relation between inner2product and included angle oftwo unit vectors.M odel 2 is on the assumption that the Earth is a revolving ellipsoid.It can besolved by the geodesic equation in differential geometry,which turns latitude of the Earth intothat of ellipsoid.For the 4 pairs of special points,their latitudes or longitudes are too close tocalculate geodesic,so we replace geodesic w ith ellipse arc,and use software M athematica toobtain the length.航 程 计 算 的 数 学 模 型谭永基(复旦大学,上海400433)摘要:本文对飞机航线飞行距离计算的数学模型进行了概述,并对2000年全国大学生数学建模竞赛的C题答卷进行了评述.假定飞机保持飞行高度10千米作匀速飞行,忽略起飞、降落和地球自转和公转的影响,C题可以归结为求飞越通过指定各点的球面或旋转椭球面上的短程线(或测地线)的航线与飞越直接连结北京上空10千米至底特律上空10千米的经过北极圈的新航线的时间差.又由于假设飞机作时速为980千米?小时的匀速飞行,问题又可归结为求相应的航程差.1地球为球体的情形取直角坐标系如下:以球心为原点,z轴指向北极,x轴通过赤道上经度为0 和180 的两点,正向指向0,y轴垂直于x轴和z轴,构成右手坐标系.在半径r=6381(千米)的球面上建立球面坐标系(,),由于航线在北半球,我们取 和北纬度一致,和东经度一致.因此航线上某处的地理坐标为(f,l),可用以下方法得到对应的球面坐标(,):=f?180=l,l为东经360-l,l为西经=?180应有