用matlab写一段程序，两个地理AB坐标LBH，已知AB两点的高度之间受地球曲率的影响，计算在地球曲率影响下从A点看到B点最远能看多少公里

在MATLAB中，我们可以使用地球曲率和两点之间的地理坐标来计算从A点看到B点的最远距离。以下是一个简单的MATLAB程序示例，用于计算这个距离：

function max_distance = calculate_max_distance(lat1, lon1, alt1, lat2, lon2, alt2)
    % 输入参数:
    % lat1, lon1, alt1 - A点的纬度、经度和高度（单位：度和米）
    % lat2, lon2, alt2 - B点的纬度、经度和高度（单位：度和米）
    
    % 地球半径（单位：米）
    R = 6371000;
    
    % 将角度转换为弧度
    lat1 = deg2rad(lat1);
    lon1 = deg2rad(lon1);
    lat2 = deg2rad(lat2);
    lon2 = deg2rad(lon2);
    
    % 计算两点间的大圆距离（单位：米）
    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 = R * c;
    
    % 计算高度差对距离的影响
    height_diff = abs(alt2 - alt1);
    adjusted_distance = sqrt(distance^2 + height_diff^2);
    
    % 输出最远可视距离（单位：公里）
    max_distance = adjusted_distance / 1000;
end

在这个函数中，我们首先定义了地球的半径，然后将输入的经纬度从度转换为弧度。接着，我们使用球面三角学公式计算两点之间的大圆距离。然后，我们考虑了两点之间的高度差，并调整了距离以考虑地球曲率的影响。最后，我们将结果从米转换为公里并返回。

要使用这个函数，只需提供A点和B点的纬度、经度和高度即可。例如：

lat1 = 34.0522; % A点的纬度
lon1 = -118.2437; % A点的经度
alt1 = 100; % A点的高度（米）
lat2 = 40.7128; % B点的纬度
lon2 = -74.0060; % B点的经度
alt2 = 30; % B点的高度（米）

max_dist = calculate_max_distance(lat1, lon1, alt1, lat2, lon2, alt2);
fprintf('从A点看到B点的最远距离是 %.2f 公里\n', max_dist);

这将计算并打印从洛杉矶到纽约的最远距离，考虑到地球曲率和高度差的影响。