已知两个点的经纬度，求第一个点到第二个点与0度经纬度交点的连线之间的距离的matlab代码

在Matlab中，计算两点间的大圆航线（Great Circle Distance）通常涉及到球面三角函数。假设你有两个经度`lat1`, `lon1`和纬度`lat2`, `lon2`的点，可以使用Haversine公式来估算这个距离。以下是一个简单的示例代码：

% 定义地球半径 (单位通常是公里)
R = 6371;

% 给定的经纬度坐标
lat1 = deg2rad(double(lat1)); % 将角度转换为弧度
lon1 = deg2rad(double(lon1));
lat2 = deg2rad(double(lat2));
lon2 = deg2rad(double(lon2));

% Haversine公式
deltaLat = lat2 - lat1;
deltaLon = lon2 - lon1;

a = sin(deltaLat/2)^2 + cos(lat1) * cos(lat2) * sin(deltaLon/2)^2;
c = 2 * atan2(sqrt(a), sqrt(1-a));

% 计算大圆距离
distance = R * c; % 单位：公里

% 如果需要到0度经线的距离，这通常是指从赤道到另一个点的距离，可以简单地加上或减去另一半地球周长
if abs(lon1) == 180 || abs(lon2) == 180
    distance = distance + R * (pi / 2); % 赤道到极点的一半圆周
else
    if lon1 > 0 && lon2 > 0 || lon1 < 0 && lon2 < 0 % 同半球
        distance = R * pi / 2; % 半圆周
    end
end

% 返回距离结果
disp(['The great circle distance is approximately ', num2str(distance), ' kilometers.']);