将经纬度转换为笛卡尔坐标通常涉及地理坐标到地球固定坐标（ECEF，地球中心地固坐标系）的转换。以下是MATLAB中用于此目的的代码示例。假设您使用的是WGS84椭球体参数，这是GPS系统使用的标准。 ```matlab function [x, y, z] = latlon_to_ecef(lat, lon, alt, a, f) % lat: 纬度 in degrees % lon: 经度in degrees % alt: 高度in meters (相对于椭球体) % a: 椭球体的长半轴 in meters % f: 扁率 % 将角度转换为弧度 lat_rad = deg2rad(lat); lon_rad = deg2rad(lon); % 计算椭球体的短半轴和偏心率平方 b = a * (1 - f); e2 = 2 * f - f^2; % 使用椭球体的参数计算N (卯酉圈曲率半径) N = a / sqrt(1 - e2 * sin(lat_rad).^2); % 计算ECEF坐标 x = (N + alt) .* cos(lat_rad) .* cos(lon_rad); y = (N + alt) .* cos(lat_rad) .* sin(lon_rad); z = ((b^2 / a) * sqrt(1 - e2 * sin(lat_rad).^2) + alt) .* sin(lat_rad); end ``` 使用此函数时，您需要提供纬度、经度、高度（相对于椭球体表面的高度，通常为0表示海平面），以及椭球体的长半轴和扁率。对于WGS84椭球体，长半轴`a`大约为6378137米，扁率`f`大约为1/298.257223563。 这是一个使用该函数的示例： ```matlab % WGS84椭球体参数 a = 6378137; % 长半轴 in meters f = 1/298.257223563; % 扁率 % 示例经纬度和高度 lat = 40.7128; % 纬度in degrees (例如，纽约的纬度) lon = -74.0060; % 经度in degrees (例如，纽约的经度) alt = 0; % 高度in meters (海平面) % 转换到ECEF坐标 [x, y, z] = latlon_to_ecef(lat, lon, alt, a, f); % 显示结果 fprintf('ECEF coordinates: x = %f, y = %f, z = %f\n', x, y, z); ```
