MATLAB语言写一段代码计算由35°27'N139°38'E驶向37°48'N122°25'W的大圆航程和恒向线航程及二者的差值

MATLAB是一种强大的数值计算软件，对于大圆航线的计算通常会涉及到球面三角学。以下是一个简单的示例，假设地球半径为6371公里，使用`atan2`函数进行弧度计算：

% 定义起点和终点经纬度
lat1 = deg2rad(35.45); % 起点纬度 (35°27')
lon1 = deg2rad(139.6333); % 起点经度 (139°38' E)
lat2 = deg2rad(37.8); % 终点纬度 (37°48')
lon2 = deg2rad(-122.4167); % 终点经度 (122°25' W)

% 地球半径
R = 6371;

% 大圆航程 (Great Circle Distance)
dLat = lat2 - lat1;
dLon = lon2 - lon1;
a = sin(dLat/2)^2 + cos(lat1) * cos(lat2) * sin(dLon/2)^2;
c = 2 * atan2(sqrt(a), sqrt(1-a));
distance_GC = R * c; % 单位：公里

% 恒向线航程 (Rhumb Line Distance)
distance_RL = R * abs(lat2 - lat1); % 纬度差直接乘以半径（简化公式）

% 差值
distance_difference = distance_GC - distance_RL;

% 结果输出
fprintf('大圆航程: %.2f 公里\n', distance_GC);
fprintf('恒向线航程: %.2f 公里\n', distance_RL);
fprintf('两者差值: %.2f 公里\n', distance_difference);

这段代码首先将经纬度转换为弧度形式，然后分别计算大圆航程和恒向线航程。最后，它显示了这两种距离的结果以及它们之间的差异。