已知一个BA无标度网络的节点编号和坐标x，y，z都保存在node_coordinates.txt文件中，已知该网络的邻接矩阵a1（a1不是txt文件）。如何通过连通中心性找出该网络的重要节点并输出它的坐标呢（前十个）？MATLAB实现，不要用现成的函数和库，不要用load函数读取

以下是MATLAB实现代码，其中假设节点编号和坐标分别保存在 nodeID 和 coordinates 两个变量中，邻接矩阵保存在 a1 变量中：
% 读取节点编号和坐标
fid = fopen('node_coordinates.txt', 'r');
nodeID = fscanf(fid, '%d');
coordinates = fscanf(fid, '%f');
fclose(fid);

% 将节点坐标转换为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 a1(i,j) == 1
C(i) = C(i) + 1/(n-1);
end
end
end

% 找出前10个重要节点
[~, idx] = sort(C, 'descend');
top10 = idx(1:10);

% 输出前10个重要节点的坐标
for i = 1:10
fprintf('Node %d: (%f, %f, %f)\n', nodeID(top10(i)), coordinates(top10(i),1), coordinates(top10(i),2), coordinates(top10(i),3));
end
登录后复制
该代码首先读取节点编号和坐标，并将坐标转换为3列矩阵。然后计算连通中心性，并找出前10个重要节点。最后输出前10个重要节点的坐标。注意，这里没有使用 load 函数读取文件，而是使用了 fopen 和 fscanf 函数手动读取文件。