matlab计算地图上两个点的球面距离

可以使用MATLAB内置的 `distance` 函数来计算地球表面上两个点的球面距离。这个函数需要输入两个经纬度坐标点，可以使用标准的[WGS 84](https://en.wikipedia.org/wiki/World_Geodetic_System)椭球模型进行计算。

例如，假设我们要计算两个点 `(lat1, lon1)` 和 `(lat2, lon2)` 之间的球面距离，可以使用以下代码：

% 设置 WGS 84 的椭球参数
a = 6378137;  % 地球长半轴
f = 1/298.257223563;  % 扁率

% 将经纬度转换为弧度
lat1_rad = deg2rad(lat1);
lon1_rad = deg2rad(lon1);
lat2_rad = deg2rad(lat2);
lon2_rad = deg2rad(lon2);

% 使用 distance 函数计算球面距离
d = distance(lat1_rad, lon1_rad, lat2_rad, lon2_rad, a, f);

其中，`lat1` 和 `lon1` 表示第一个点的纬度和经度，`lat2` 和 `lon2` 表示第二个点的纬度和经度。`a` 和 `f` 是椭球参数，可以根据需要进行调整。最后的球面距离 `d` 的单位是米。