参数说明:
m:带权邻接矩阵
start:起点设置
last:终点设置
m求解方法:
clc;clear;
%功能在给定的输入输出以及矩阵m的情况下,求解出最短路径,并输出路径曲线以及最短路径
disp('请输入矩阵m')
m=input('m=');
[Array,Array]=size(m);
dmin=zeros(Array,1); %用于存放到每个点的最短距离
P=zeros(Array,1); %确定那一点已经被标记了,主要用于主循环
path=zeros(Array,1); %存放路径
%%终点设置
last=9;
%%起点初始化
start=1;
P(start)=1;
dmin(start)=0;
path(1)=start;
n=2;
%主判断循环
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;
path(n)=pb;
n=n+1;
if(pb==last)
break;
end
end
disp('最短距离为:')
disp(min);
disp('最短路径规划为:')
disp(path);