使用MATLAB计算北京到纽约的新旧航线里程

数学实验的PPT课件参见百度文库
http://wenku.baidu.com/view/2eda008884868762caaed5a6.html
本文给出了实验的示例程序以及输出结果。

以下是mlab21.m的代码，用于计算北京到纽约的航线例程：

执行mlab21后得到以下结果

Dmatrix =
             0       1144.90       2155.64       9608.53      10993.77
       1144.90             0       1766.78       9936.44      11870.27
       2155.64       1766.78             0       8283.06      10764.64
       9608.53       9936.44       8283.06             0       4061.47
      10993.77      11870.27      10764.64       4061.47             0

可以清楚的看到，北京直飞纽约的航线例程为10993.77km，
而旧航线里程为1144.90+1766.78+8283.06+4061.47= 15256.21 (km).
假设飞机按照统计的平均时速980km/h匀速前进，那么新航线比旧航线至少节约了4.35小时，这尚未包含在各个中转站花费的候机时间。

以下是mlab22.m的代码，用于在球面上绘制新旧航线：


这一段程序调用了skyway(p1,p2,color)函数，代码如下：


执行mlab.22之后，可以在Figure(1)窗口中看到绘制的新旧航线图：

视角一(示太平洋)

视角二(示北极)