首先需要创建一个零矩阵zeros(n),其中n为图论问题中节点的个数。再创建一个稀疏矩阵sparse(),
其中对sparse赋个变量,比如s=sparse([起点集合],[对应终点集合],[对应权重集合])。最后令d=sparse(s)输出结果
s=sparse([1 1 1 2 2],[2 3 4 3 4],[2 3 8 6 6]);
g=sparse(s)
再加个函数
view(biograph(G,[],'ShowWeights','on'))得到图形
再加入graphallshortestpaths(G)得出每个节点与对应各个节点的距离,G为稀疏矩阵的变量。
但这样我们还需要自己计算路径总和
可以考虑通过另一种函数来求解
% 构造邻接矩阵
G = zeros(6);
G = graph([6 1 2 2 3 4 4 5 5 6 1],[2 6 3 5 4 1 6 3 4 3 5],[41 99 51 32 15 45 38 32 36 29 21])
plot(G,'EdgeLabel',G.Edges.Weight)
[P,d] = shortestpath(G,1,6)
通过这组代码不仅能做出图形还能直接求出需要的距离,不过这组代码是来求解无向图的。
有向图代码如下:
function [min,path]=dijkstra(w,start,terminal)
n=size(w,1); label(start)=0; f(start)=start;
for i=1:n
if i~=start
label(i)=inf;
end, end
s(1)=start; u=start;
while length(s)<n
for i=1:n
ins=0;
for j=1:length(s)
if i==s(j)
ins=1;
end
end
if ins==0
v=i;
if label(v)>(label(u)+w(u,v))
label(v)=(label(u)+w(u,v));
f(v)=u;
end
end
end
v1=0;
k=inf;
for i=1:n
ins=0;
for j=1:length(s)
if i==s(j)
ins=1;
end
end
if ins==0
v=i;
if k>label(v)
k=label(v); v1=v;
end
end
end
s(length(s)+1)=v1;
u=v1;
end
min=label(terminal); path(1)=terminal;
i=1;
while path(i)~=start
path(i+1)=f(path(i));
i=i+1 ;
end
path(i)=start;
L=length(path);
path=path(L:-1:1);
这段代码为定义一个函数,保存在文件中即可。
% 构造邻接矩阵
a = zeros(6);
a = sparse([6 1 2 2 3 4 4 5 5 6 1],[2 6 3 5 4 1 6 3 4 3 5],[41 99 51 32 15 45 38 32 36 29 21])
a = a + a';
a(a==0) = inf; % 零元素换成inf
a(eye(6,6)==1)=0; % 对角线换成 0
view(biograph(a,[],'ShowWeights','on'))
[min,path]=dijkstra(a,1,6) % dijkstra模型求解节点一到节点六最短路径
在做题时只需要修改这段代码的参数即可。