matlab计算大圆航程

目录

1.问题描述 

2.问题分析

 3.公式简化与求解

 4.例题以及MATLAB程序

---

1.问题描述 

给定起始点和到达点的纬度和经度，分别记为（，），（，）计算大圆航程（弧度表示）

 

2.问题分析

已知条件：起始点纬度，到达点纬度，起始点经度，到达点经度

可以推出：起始点余纬90°-，到达点余纬90°-，起始点与到达点的经差 = -

因此，已获得两边及其夹角，可以利用球面三角形的边余弦公式进行求解。

 

 3.公式简化与求解

球面三角形的边余弦公式如下：

代入已知条件得：

化简得：

因此：

 4.例题以及MATLAB程序

某船从33°53.3′S、18°23.1'E，驶向40°27.1'N、73°49.'W,拟采用大圆航线航行，求大圆航程 

 程序运行结果：
大圆航程为112.712035°，即6762.722073′

注意：

本文采用的是单独编写matlab函数文件，然后再命令行窗口进行调用函数的方式。
输入经纬度时，东经北纬为（+），西经南纬为（-）