已经知道有八百多个经纬度的中心点和一千多个离散的已知经纬度的点，要求中心点附近五千米内有多少个离散的点并统计总数，请给出matlab代码

下面是一个使用 Matlab 实现的示例代码，用于计算中心点附近五千米范围内的离散点数量并统计总数：

% 中心点的经纬度
center_lat = 40.7128; % 维度
center_lon = -74.0060; % 经度

% 离散点的经纬度
discrete_data = load('discrete_data.txt');
discrete_lat = discrete_data(:, 1);
discrete_lon = discrete_data(:, 2);

% 计算距离
R = 6371; % 地球半径（单位：千米）
dlat = deg2rad(discrete_lat - center_lat);
dlon = deg2rad(discrete_lon - center_lon);
a = sin(dlat/2).^2 + cos(deg2rad(center_lat)) .* cos(deg2rad(discrete_lat)) .* sin(dlon/2).^2;
c = 2 * atan2(sqrt(a), sqrt(1-a));
dist = R * c; % 距离（单位：千米）

% 找出距离在五千米范围内的离散点
indices = find(dist <= 5);

% 统计总数
num_points = numel(indices);

% 输出结果
disp(['There are ', num2str(num_points), ' points within 5 kilometers of the center point.']);

请确保将上述代码中的`discrete_data.txt`替换为实际存储离散点数据的文件名，并且数据文件应该是一个包含经纬度的矩阵，每一行代表一个离散点的经纬度信息。这个示例代码使用 Haversine 公式来计算中心点与离散点之间的距离，并找出距离在五千米范围内的离散点的数量。最后，统计数量并输出结果。