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

如果需要将多个经纬度转换为 xy 坐标，可以使用以下 Matlab 代码：

% 输入经纬度，单位为度
lat = [39.9, 40.0, 40.1];
lon = [116.3, 116.4, 116.5];

% 将经纬度转换为弧度
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 坐标。需要注意的是，使用点乘 `.*` 和向量乘 `*` 运算符可以使程序正确处理多个点的情况。