MATLAB写出邻接矩阵A并计算出可达矩阵M

假设我们已经有了一个有向图的邻接矩阵 A，其中 A(i,j) = 1 表示从节点 i 可以到达节点 j，A(i,j) = 0 表示从节点 i 不能到达节点 j。则可以使用 MATLAB 计算可达矩阵 M 的算法如下：

1. 初始化可达矩阵 M 为邻接矩阵 A，即 M = A；

2. 重复以下步骤，直到 M 不再发生变化：

a. 计算 M 的新值，即 M = M + M * A；

b. 将 M 中的所有非零元素设为 1，即 M(M ~= 0) = 1。

下面是一个 MATLAB 代码示例：

% 假设有一个邻接矩阵 A，表示有向图的连通情况
A = [0 1 0 0 0;
     1 0 1 1 0;
     0 1 0 0 1;
     0 0 0 0 1;
     0 0 0 1 0];

% 初始化可达矩阵 M 为 A
M = A;

% 计算可达矩阵 M
while true
    M_new = M + M * A;
    if isequal(M_new, M)
        break;
    end
    M = M_new;
    M(M ~= 0) = 1;
end

% 输出可达矩阵 M
M

在这个示例中，我们假设有一个 5 个节点的有向图，并给出了其邻接矩阵 A。然后，我们使用上述算法计算可达矩阵 M，并输出结果。