%第二层连通介数中心性 fid = fopen('node_coordinates2.txt'); C = textscan(fid, 'Node %d: (%f,%f,%f)'); fclose(fid); coordinates = [C{2}, C{3}, C{4}]; nodeID = C{1}; % 将节点坐标转换为3列矩阵 coordinates = reshape(coordinates, [3, length(nodeID)]); coordinates = coordinates'; % 计算连通中心性 n = length(nodeID); C = zeros(n,1); for i = 1:n for j = 1:n if a2(i,j) == 1 C(i) = C(i) + 1/(n-1); end end end % 找出前10个重要节点 [~, idx] = sort(C, 'descend'); if length(idx) >= 10 top10 = idx(1:10); else top10 = idx; end % 输出前10个重要节点的坐标 for i = 1:length(top10) fprintf('Node %d: (%f, %f, %f)\n', nodeID(top10(i)), coordinates(top10(i),1), coordinates(top10(i),2), coordinates(top10(i),3)); end换一种数据读取方式，不要用datadata

假设你的节点坐标数据是存储在一个文本文件中，每一行对应一个节点的ID和三个坐标值，格式如下：

1 0.1 0.2 0.3
2 0.4 0.5 0.6
3 0.7 0.8 0.9
...

你可以使用`importdata`函数读取这个文本文件，然后提取出节点ID和坐标值：

data = importdata('node_coordinates.txt');
nodeID = data(:, 1);
coordinates = data(:, 2:4);

接下来，你可以使用这些节点坐标计算连通中心性。这里给出一个计算连通中心性的函数：

function C = compute_closeness_centrality(adjacency_matrix)
% adjacency_matrix: 节点的邻接矩阵
n = size(adjacency_matrix, 1);
C = zeros(n, 1);

for i = 1:n
    % 使用BFS算法计算节点i到其他节点的最短路径长度
    d = bfs(adjacency_matrix, i);
    % 计算节点i到其他节点的平均最短路径长度
    C(i) = (n - 1) / sum(d);
end

end

function d = bfs(adjacency_matrix, start_node)
% 使用BFS算法计算节点start_node到其他节点的最短路径长度
n = size(adjacency_matrix, 1);
d = Inf(n, 1);
d(start_node) = 0;
visited = false(n, 1);
queue = start_node;
while ~isempty(queue)
    u = queue(1);
    queue(1) = [];
    visited(u) = true;
    for v = 1:n
        if adjacency_matrix(u, v) == 1 && ~visited(v)
            d(v) = d(u) + 1;
            visited(v) = true;
            queue(end+1) = v;
        end
    end
end

end

这个函数接受一个节点的邻接矩阵作为输入，返回每个节点的连通中心性。你可以将这个函数应用到你的邻接矩阵上：

adjacency_matrix = your_adjacency_matrix;
C = compute_closeness_centrality(adjacency_matrix);

最后，你可以找到连通中心性排名前10的节点，并输出它们的ID和坐标：

[~, idx] = sort(C, 'descend');
if length(idx) >= 10
    top10 = idx(1:10);
else
    top10 = idx;
end

for i = 1:length(top10)
    fprintf('Node %d: (%f, %f, %f)\n', nodeID(top10(i)), coordinates(top10(i),1), coordinates(top10(i),2), coordinates(top10(i),3));
end