2005年华东师大参加全国数模竞赛的候选队作品3课程.docVIP

上海到其他35个大城市的距离问题 模型的建立 建立以地球球心为坐标原点O,赤道平面为XOY平面,地球半径为r,并且使得上海在XOY平面上的投影在OX轴正向上的坐标系.设上海(点A)的纬度为A,经度为A;球面上任意点B的纬度为,经度为.再设点A的坐标为(xA,yA,zA),点B的坐标为(x,y,z).将B点投影到xoy平面上,设投影点为C,如图1. 则:,,. 下面考察过A,B两点的大圆,设圆心为O＇,如图2.A,B的弦长为: ;夹角为:,弧AB长s为:.(二)程序具体程序见附录1.程序中的jingwei为36*2的矩阵,第一列是各大城市的纬度,第二列是经度,并记上海为36号.为计算方便,把经(纬)度的分转换成小数. (三)结果(单位:千米)1阿姆斯特丹 2安卡拉 3雅典 4奥克兰 5曼谷 8880 7738 8545 9378 2887 6巴萨罗那 7北京 8柏林 9布鲁塞尔 10布达佩斯 9831 1074 8403 9013 8348 11开罗 12哥本哈根 13哈瓦那 14赫尔辛基 15香港 8349 8258 13464 7379 1220 16约翰内斯堡 17吉隆坡 18伦敦 19澳门 20墨西哥城 11767 3743 9198 1265 12914 21莫斯科 22新德里 23纽约 24奥斯陆 25巴黎 6823 4250 11858 8093 9273 26罗马 27斯德哥尔摩 28悉尼 29台北 30东京 9125 7774 7880 677 1766 31多伦多 32维也纳 33华沙 34惠灵顿 35苏黎世 11421 8482 7973 9735 9018 上海到5个大城市的路线所过最近城市问题 (一)模型的建立 仅以上海到伦敦为例建立模型.由于题目给了许多的候选城市,故考虑做初步筛选,得到较少且距离该路线较近的城市.我们借助了一款3D World Map 2.1(华军软件园)的软件,利用软件提供的环境,根据经度和纬度的基本范围,筛选出了30多个在路线附近,或是在上海,伦敦旁边的候选城市. 设上海为A点,坐标(),伦敦为B点,坐标(),经过A,B的大圆为.这两点球面距离为9198千米,还不到大圆圆弧的四分之一,弦AB远小于直径.如图3, 过AB的直线L的参数方程为: . 则L上任一点C到圆心O的距离. 直线OC交于D点,坐标(),则: . 设与L平行的直径为MN.当时,C在弦AB上,故D在圆弧AB上;当时,C趋向BA延长线的无穷远处,此时DM;同理,时,DN. 由此可知,D的取点范围为半圆弧MN.取一个候选城市,记为E点,坐标()E在上的投影点记为F点,如图4. 通过适当的筛选城市,首先,可以保证F不可能在直径MN上,因为如果这样,那么E点必在与垂直的大圆上,由于圆弧AB太短(还不到大圆圆弧的四分之一),无论E点怎么取,其到圆弧AB的球面距离总是大于我们之前所选择的候选城市.其次,还可以保证F点在的过直径MN和弦AB的半圆上.(如果F点在另一个半圆,那么E点到圆弧AB的球面距离还是太大). 设E点到D点球面距离为S,则:,故求的问题转化为求. 令,. 由于圆弧MN光滑,故在上连续可微.,,由数分知识,(,. 另外,)在上有唯一有限极小值点,即为最小值点(证明见注1). 解方程:对 求导数(,实际上有三个根:.)在Matlab6.5.1上运行解方程的程序,(只)得到唯一根,因此. ①,点在圆弧AB上,此时E点到圆弧AB的最短球面距离即为;②,点在圆弧AM上,此时E点到圆弧AB的最短球面距离即为EA的球面距离; ③,点在圆弧BN上,此时E点到圆弧AB的最短球面距离即为EB的球面距离;缺少证明这样,就求出了E点到圆弧AB的最短球面距离,即得到了某个候选城市的最短距离.可以通过在程序中用＂for＂循环语句,得到所有候选城市到弧AB的最短距离S(i),然后再找到S(i)中的最小值S(n)及其对应的城市n的坐标和点的坐标(②,为A或B的坐标),由坐标可进一步得经(纬)度(算法见注2).注1:验证上有唯一有限极小值点,即为最小值点. 证明:,在E点给定的情况下,不等式右边为常量,等号成立当且仅当OF?D三点共线,此时D对应着,由于筛点保证了F不在直径MN上,故D不取MN点,从而为的唯一最小值点.设OD和OF的夹角为(), .当时,随着的增大,单调递减至0,递增至1,递减到最小值;当时,随着的增大,单调递增至, 递减至,递增到(为任意小的正数).由于为常量,因此和的随t变化规律大致相同.如图5 故上有唯一有限极小值点L(),即为最小值点. 注2:设球面上一点的坐标(),纬度为,经度为,上海的经度为A.则根据问题1中的坐标转换公式,,.(二)