%第二层连通介数中心性 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换一种数据读取方式

这里可以使用函数`importdata`来读取数据文件，它可以自动根据数据格式读取数据。

首先读取节点坐标数据：

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

然后计算连通中心性：

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

相关问题

fid = fopen('node_coordinates.txt'); C = textscan(fid, 'Node %d: (%f,%f,%f)'); fclose(fid); nodes = [C{2}, C{3}, C{4}]; node_ids = C{1}; node_pos = nodes(:,1:2); [BC1_sorted, BC1_idx] = sort(BC1, 'descend'); % 将介数中心性从高到低排序并记录排序后的索引 top10_idx = BC1_idx(1:10); % 取前5个节点的索引 fid = fopen('resultjieshu1.txt', 'w'); % 打开文件 for i = 1:10 node_id = idx(top10_idx(i)); % 使用排序后的索引来获取节点编号 node_bc = BC1(node_id); node_x = node_pos(top10_idx(i), 1); % 使用未排序的索引来获取节点坐标 node_y = node_pos(top10_idx(i), 2); fprintf('节点 %d，介数中心性为 %f，坐标为 (%f,%f)\n', node_id, node_bc, node_x, node_y); end 改成输出xyz轴坐标

要将节点的xyz轴坐标输出，只需要将node_pos变量改为nodes变量即可。修改后的代码如下：

fid = fopen('node_coordinates.txt');
C = textscan(fid, 'Node %d: (%f,%f,%f)');
fclose(fid);

nodes = [C{2}, C{3}, C{4}];
node_ids = C{1};
node_pos = nodes(:,1:3);

[BC1_sorted, BC1_idx] = sort(BC1, 'descend'); % 将介数中心性从高到低排序并记录排序后的索引
top10_idx = BC1_idx(1:10); % 取前5个节点的索引

fid = fopen('resultjieshu1.txt', 'w'); % 打开文件
for i = 1:10
    node_id = idx(top10_idx(i)); % 使用排序后的索引来获取节点编号
    node_bc = BC1(node_id);
    node_x = node_pos(top10_idx(i), 1); % 使用未排序的索引来获取节点x坐标
    node_y = node_pos(top10_idx(i), 2); % 使用未排序的索引来获取节点y坐标
    node_z = node_pos(top10_idx(i), 3); % 使用未排序的索引来获取节点z坐标
    fprintf('节点 %d，介数中心性为 %f，坐标为 (%f,%f,%f)\n', node_id, node_bc, node_x, node_y, node_z);
end

这样修改后，程序会输出每个节点的编号、介数中心性和xyz轴坐标信息。

fid = fopen('node_coordinates555.txt', 'r');怎么查看内容

可以使用fscanf函数读取文件内容并打印到命令行窗口上。示例代码如下：

fid = fopen('node_coordinates555.txt', 'r');
if fid == -1
    error('Failed to open file!');
end

while ~feof(fid)  % 判断文件是否读取完毕
    line = fgetl(fid);  % 读取一行内容
    disp(line);  % 打印到命令行窗口上
end

fclose(fid);  % 关闭文件

这段代码会逐行读取文件内容并打印到命令行窗口上，直到文件读取完毕。需要注意的是，打印出来的内容可能比较多，可能需要在命令行窗口中手动滚动查看。