设计一段matlab代码，将经纬度坐标转换为xy平面坐标

要将经纬度坐标转换为 xy 平面坐标，需要先将经纬度坐标转换为投影坐标系下的坐标，然后再进行平面坐标的转换。以下是一段 Matlab 代码，可以实现将经纬度坐标转换为 xy 平面坐标的功能：

% 定义经纬度坐标
lat = 39.91; % 纬度
lon = 116.40; % 经度

% 定义投影坐标系
utmzone = utmzone(lon, lat); % 获取 UTM 投影带
proj = utmproj(utmzone); % 获取 UTM 投影坐标系

% 将经纬度坐标转换为投影坐标系下的坐标
[x, y] = projfwd(proj, lat, lon);

% 将投影坐标系下的坐标转换为 xy 平面坐标
x0 = 0; % 定义原点 x 坐标
y0 = 0; % 定义原点 y 坐标
xy = [x-x0, y-y0]; % 计算 xy 平面坐标

% 显示转换结果
disp(['经纬度坐标：(', num2str(lat), ', ', num2str(lon), ')']);
disp(['xy 平面坐标：(', num2str(xy(1)), ', ', num2str(xy(2)), ')']);

这段代码首先定义了经纬度坐标，然后根据经纬度坐标获取 UTM 投影带和投影坐标系。接下来，使用 `projfwd` 函数将经纬度坐标转换为投影坐标系下的坐标。最后，将投影坐标系下的坐标转换为 xy 平面坐标，并显示转换结果。您可以将经纬度坐标替换为您的数据，然后运行这段代码，即可得到相应的 xy 平面坐标。