clc;clear;
%功能在给定的输入输出以及矩阵m的情况下,求解出最短路径,并输出路径曲线以及最短路径
disp('请输入矩阵m')
m=input('m=');
[~,Array]=size(m);
dmin=zeros(Array,1);
P=zeros(Array,1);
path=zeros(Array,1);
%%起点初始化
start=1;
P(start)=1;
dmin(start)=0;
while sum(P)<Array
tb=find(P==0);%找到还未找到最短路径的点
fb=find(P);%找出已找到最短路径的点
min=inf;
for i=1:length(fb)
for j=1:length(tb)
longh=dmin(fb(i))+m(fb(i),tb(j));
if longh<min
min=longh;
dmin(tb(j))=min;
end
end
end
pb=find(dmin==min);
P(pb)=1;
end