帮我写一个经典DV-Hop算法的代码，用matlab

好的，以下是一份经典DV-Hop算法的Matlab代码：

function dv_hop(nodes, anchors, th)
% nodes: 所有节点的位置，Nx2矩阵，每行代表一个节点的坐标
% anchors: 锚点的位置，Mx2矩阵
% th: 阈值，用于计算每个节点的跳数

% 计算节点之间的距离
dist = pdist2(nodes, nodes);

% 计算每个节点到锚点的距离
dist_anchors = pdist2(nodes, anchors);

% 计算每个节点的跳数
hop = ceil(dist_anchors / th);

% 计算节点之间的平均跳数
avg_hop = mean(hop, 2);

% 绘制节点和锚点的位置
figure;
plot(nodes(:,1), nodes(:,2), 'ro'); hold on;
plot(anchors(:,1), anchors(:,2), 'bs');
legend('Nodes', 'Anchors');

% 根据跳数绘制节点之间的连线
for i = 1:size(nodes, 1)
    for j = i+1:size(nodes, 1)
        if hop(i,j) <= max(avg_hop)
            plot([nodes(i,1), nodes(j,1)], [nodes(i,2), nodes(j,2)], 'k--');
        end
    end
end

end

运行该函数需要输入三个参数：所有节点的位置、所有锚点的位置和阈值。函数会根据DV-Hop算法计算出每个节点之间的跳数，并根据跳数绘制节点之间的连线。代码中使用了Matlab自带的函数pdist2来计算距离矩阵。