将经纬度数组转换为xy数组

可以使用循环遍历每个经纬度点，并将其转换为对应的xy坐标点，最终得到一个xy坐标数组。以下是示例代码：

% 定义地球半径和转换参数
R = 6378137; % 地球半径（单位：米）
lat0 = 0; % 参考纬度
lon0 = 0; % 参考经度

% 输入经纬度数组
latlon = [31.2304, 121.4737; 39.9042, 116.4074; 23.1291, 113.2644]; % 每行为一个经纬度点

% 将经纬度转换为xy坐标
xy = zeros(size(latlon));
for i = 1:size(latlon, 1)
    % 将经纬度转换为弧度
    lat = latlon(i, 1) * pi / 180;
    lon = latlon(i, 2) * pi / 180;
    lat0 = lat0 * pi / 180;
    lon0 = lon0 * pi / 180;

    % 计算转换参数
    k = R * cos(lat0) / sqrt(1 - sin(lat0)^2 * sin(lon - lon0)^2);

    % 计算xy坐标
    x = k * cos(lat) * (lon - lon0);
    y = k * sin(lat);
    
    % 存储xy坐标
    xy(i, :) = [x, y];
end

其中，`latlon`为输入的经纬度数组，`xy`为转换后的xy坐标数组。遍历每个经纬度点，使用前面提到的转换公式计算出对应的xy坐标，并存储在`xy`数组中。最终得到的`xy`数组与`latlon`数组行数相同，每行对应一个xy坐标点。