利用MATLAB写这样一段代码：在一个平面内有26个点，每一个点的位置已知，我到达每个点的时间已知，对每个点的时间进行排序，先从时间最小的点为圆心找到距离该点小于固定阈值距离的其他点，并且用小于固定阈值距离的其他点与圆心的距离差除以他们的到达时间差得到一速度矩阵

% 假设点的位置已知，到达时间存储在矩阵T中，阈值为r
% 假设点的位置存储在矩阵P中，每个点的位置存储在P的一行中

n = size(P,1); % n为点的数量
V = zeros(n,n); % 初始化速度矩阵为0

for i = 1:n % 对每个点进行循环
[~,index] = sort(T); % 对到达时间进行排序
neighbors = index(index~=i); % 找到除自己以外的其他点
for j = 1:length(neighbors) % 对每个邻居点进行循环
if norm(P(i,:)-P(neighbors(j),:)) < r % 如果距离小于阈值
V(i,neighbors(j)) = (norm(P(i,:)-P(neighbors(j),:))-r)/(T(neighbors(j))-T(i)); % 计算速度矩阵
end
end
end