matlab经纬度距离计算

在 MATLAB 中，计算两个经纬度之间的距离通常使用 Haversine 公式，也称为球面两点间的大圆距离公式。这个公式考虑了地球并非完美的球体，而是稍微扁平的椭球体。以下是一个简单的步骤来计算经纬度之间的距离：

1. 定义经纬度：你需要两个点的经度（longitude）和纬度（latitude），例如 `lat1`, `lon1` 和 `lat2`, `lon2`。

lat1 = 37.7749; % 假设第一个点的纬度
lon1 = -122.4194; % 假设第一个点的经度
lat2 = 40.7128; % 假设第二个点的纬度
lon2 = -74.0060; % 假设第二个点的经度

2. 使用 Haversine 函数：MATLAB 提供了一个名为 `haversine` 或 `geodist` 的函数，可以直接计算两个点之间的距离。

[distance, azimuth] = geodist(lat1, lon1, lat2, lon2, 'deg'); % 使用默认单位为度

`distance` 将返回两个点之间的距离（通常以公里为单位），而 `azimuth` 返回两个点之间的方位角，但在此场景下可能不需要。

如果你想要自己实现 Haversine 公式，你可以使用下面的代码：

R = 6371; % 地球平均半径，单位为公里
dLat = deg2rad(lat2-lat1);
dLon = deg2rad(lon2-lon1);
a = sin(dLat/2)^2 + cos(deg2rad(lat1)) * cos(deg2rad(lat2)) * sin(dLon/2)^2;
c = 2 * atan2(sqrt(a), sqrt(1-a));
distance = R * c;