加载"lab1data.txt".里面是50个节点的位置信息，用matlab写一个gossiping协议，每个节点随机向邻居4个节点发送数据包，计算sink节点到达所有节点的跳数和数据包

传输次数。

首先，读取节点位置信息和邻接矩阵：

data = load('lab1data.txt');
nodes = size(data, 1);
adj = zeros(nodes, nodes);
for i = 1:nodes
    for j = 1:nodes
        if i ~= j && norm(data(i,:) - data(j,:)) <= 0.2
            adj(i,j) = 1;
        end
    end
end

然后，实现gossiping协议：

sink = randi(nodes); % 选择一个随机的sink节点
hop_count = zeros(1, nodes); % 记录每个节点到达sink节点的跳数
data_sent = zeros(1, nodes); % 记录每个节点发送的数据包数量
for i = 1:nodes
    if i ~= sink
        neighbors = find(adj(i,:)); % 找到邻居节点
        for j = 1:4 % 随机选择4个邻居节点发送数据包
            if ~isempty(neighbors)
                k = randi(length(neighbors));
                neighbor = neighbors(k);
                neighbors(k) = [];
                hop_count(i) = max(hop_count(i), hop_count(neighbor) + 1);
                data_sent(i) = data_sent(i) + 1;
            end
        end
    end
end

最后，输出结果：

fprintf('Sink node: %d\n', sink);
fprintf('Max hop count: %d\n', max(hop_count));
fprintf('Total data sent: %d\n', sum(data_sent));

完整代码如下：

data = load('lab1data.txt');
nodes = size(data, 1);
adj = zeros(nodes, nodes);
for i = 1:nodes
    for j = 1:nodes
        if i ~= j && norm(data(i,:) - data(j,:)) <= 0.2
            adj(i,j) = 1;
        end
    end
end

sink = randi(nodes);
hop_count = zeros(1, nodes);
data_sent = zeros(1, nodes);
for i = 1:nodes
    if i ~= sink
        neighbors = find(adj(i,:));
        for j = 1:4
            if ~isempty(neighbors)
                k = randi(length(neighbors));
                neighbor = neighbors(k);
                neighbors(k) = [];
                hop_count(i) = max(hop_count(i), hop_count(neighbor) + 1);
                data_sent(i) = data_sent(i) + 1;
            end
        end
    end
end

fprintf('Sink node: %d\n', sink);
fprintf('Max hop count: %d\n', max(hop_count));
fprintf('Total data sent: %d\n', sum(data_sent));