matlab 方位角 经纬度,经纬度转方位角matlab程序

以下是经纬度转方位角的 Matlab 程序示例：

function [azimuth, elevation] = ll2az(lat1, lon1, lat2, lon2)
% 将经纬度坐标转换为方位角和仰角
% 输入参数：
% lat1: 起点纬度，单位为度
% lon1: 起点经度，单位为度
% lat2: 终点纬度，单位为度
% lon2: 终点经度，单位为度
% 输出参数：
% azimuth: 方位角，单位为度
% elevation: 仰角，单位为度

% 将经纬度转化为弧度
lat1 = deg2rad(lat1);
lon1 = deg2rad(lon1);
lat2 = deg2rad(lat2);
lon2 = deg2rad(lon2);

% 计算方位角和仰角
delta_lon = lon2 - lon1;
y = sin(delta_lon) * cos(lat2);
x = cos(lat1) * sin(lat2) - sin(lat1) * cos(lat2) * cos(delta_lon);
azimuth = atan2(y, x);
elevation = atan2(sin(lat2 - lat1), cos(lat1) * cos(lat2) * cos(delta_lon));

% 将弧度转化为度数
azimuth = rad2deg(azimuth);
elevation = rad2deg(elevation);

示例用法：

lat1 = 30.5225;
lon1 = 114.3586;
lat2 = 30.5178;
lon2 = 114.3553;
[azimuth, elevation] = ll2az(lat1, lon1, lat2, lon2);
fprintf('起点(%f, %f)，终点(%f, %f)的方位角为%f度，仰角为%f度。\n', lat1, lon1, lat2, lon2, azimuth, elevation);