在地理学中，大圆距离是地球上两点沿最短路径（大圆弧）的距离。WGS84是一个广泛使用的地球椭球模型，用于全球定位系统（GPS）。以下是使用MATLAB计算基于WGS84椭球模型的两点间大圆距离的示例代码： ```matlab function distance = calculate_great_circle_distance(lat1, lon1, lat2, lon2) % 将经纬度从度转换为弧度 lat1_rad = deg2rad(lat1); lon1_rad = deg2rad(lon1); lat2_rad = deg2rad(lat2); lon2_rad = deg2rad(lon2); % 使用Haversine公式计算大圆距离 % Haversine公式是一种计算球面上两点间大圆距离的近似方法 dlon = lon2_rad - lon1_rad; dlat = lat2_rad - lat1_rad; a = sin(dlat/2).^2 + cos(lat1_rad) .* cos(lat2_rad) .* sin(dlon/2).^2; c = 2 * atan2(sqrt(a), sqrt(1-a)); % 使用WGS84椭球模型的平均半径（单位：米） R = 6371000; % 平均半径，单位：米 % 计算距离 distance = R * c; end ``` 要使用此函数，只需提供两个点的纬度和经度（以度为单位）作为输入，然后调用该函数即可获得两点间的大圆距离（以米为单位）。例如： ```matlab lat1 = 40.7128; % 第一个点的纬度 lon1 = -74.0060; % 第一个点的经度 lat2 = 51.5074; % 第二个点的纬度 lon2 = -0.1278; % 第二个点的经度 distance = calculate_great_circle_distance(lat1, lon1, lat2, lon2); fprintf('两点间的大圆距离为：%.2f 米\n', distance); ```
