除了显示的代码中没有23行以外:)，问题实际上就在这一行中
距离= (deg2km(距离(lat1 (z), lon1 (z), lat2, lon2)));%计算距离
因为你只保存每次循环中计算的当前值，所以你只保存最后一个值。然后，同样的问题出现在下一行。
类似的
．..
距离= 0(元素个数(stname), 1);%预分配输出数组
为z = 1:元素个数(stname)
距离(z) = deg2km(距离(lat1 (z), lon1 (z), lat2, lon2));
结束
Withinkm = distance <= 2;找到距离小于2公里的行索引。
%从所有点集合中只提取那些关闭点。
closePointskm = [lat2(withinkm)， lon2(withinkm)， c(withinkm)， distance (withinkm)， find(withinkm)];
．..
这将写入数组中找到/保存了相应值的行。
然后，您还需要修复调用的语法
writematrix
