经纬度转化为xy坐标系matlab具体代码

在 Matlab 中，可以使用以下代码将经纬度转换为 xy 坐标：

% 输入经纬度，单位为度
lat = 39.9;
lon = 116.3;

% 将经纬度转换为弧度
lat_rad = deg2rad(lat);
lon_rad = deg2rad(lon);

% 计算地球的半径
R = 6371e3;

% 定义参考经度和参考纬度
lat_ref = 39.9;
lon_ref = 116.3;

% 将参考经纬度转换为弧度
lat_ref_rad = deg2rad(lat_ref);
lon_ref_rad = deg2rad(lon_ref);

% 计算两点之间的距离
d = R * acos(sin(lat_rad)*sin(lat_ref_rad) + cos(lat_rad)*cos(lat_ref_rad)*cos(lon_rad-lon_ref_rad));

% 计算 xy 坐标
x = d * cos(lat_rad) * sin(lon_rad-lon_ref_rad);
y = d * (cos(lat_ref_rad)*sin(lat_rad) - sin(lat_ref_rad)*cos(lat_rad)*cos(lon_rad-lon_ref_rad));

其中，`lat` 和 `lon` 分别为输入的纬度和经度，`lat_ref` 和 `lon_ref` 是参考经纬度，可以选择与输入经纬度相同。`d` 表示两点之间的距离，`x` 和 `y` 分别为转换后的 xy 坐标。