卫星经纬度计算仰角的matlab代码

在 MATLAB 中，计算地球表面某点（基于给定的经纬度）到卫星位置的仰角（即视线与水平面的夹角）通常涉及使用球面三角学。以下是一个简单的示例代码，假设你已经有了卫星的纬度和经度以及观察点的经纬度：

% 定义地球半径（以公里为单位）
R = 6371;

% 卫星经纬度（假设为固定值）
satellite_lat = 40; % 卫星纬度（度）
satellite_lon = -120; % 卫星经度（度）

% 观察点经纬度（用户输入或变量）
observer_lat = input('请输入观察点的纬度（度）: ');
observer_lon = input('请输入观察点的经度（度）: ');

% 将角度转换为弧度
lat_sat = deg2rad(satellite_lat);
lon_sat = deg2rad(satellite_lon);
lat_obs = deg2rad(observer_lat);
lon_obs = deg2rad(observer_lon);

% 计算两点之间的初始距离差（忽略高度）
dlon = lon_sat - lon_obs;
dlat = lat_sat - lat_obs;

% 使用球面余弦公式计算距离
a = sin(dlat/2)^2 + cos(lat_sat) * cos(lat_obs) * sin(dlon/2)^2;
c = 2 * atan2(sqrt(a), sqrt(1-a));

% 仰角计算（注意方向，通常视线向上是正角度，向下是负角度）
elevation_angle = asin(sin(lat_sat) * sin(lat_obs) + cos(lat_sat) * cos(lat_obs) * cos(c)) * R;

% 显示仰角结果（单位：度）
elevation_angle_degrees = rad2deg(elevation_angle);
disp(['观察点到卫星的仰角约为 ' num2str(elevation_angle_degrees, '%.2f') ' 度']);

% 相关问题--