分享Dijkstra_算法_Matlab_实现.docVIP

分享 Dijkstra 算法的 Matlab 实现 function?[D I] = dijkstra(W,k) %??Dijkstra?算法??——??by Glorevo %??找出赋权图第?k?个顶点到其余各个顶点的最小权以及最短路 %??W?为带权邻接矩阵 %??返回值?D?为最小权，I?为最短路 ? n = size(W,1);??%???顶点个数 V = 1:n;????????%???初始化待作为终点的顶点向量 V(k) = 0;???????%???0?表示不再作为终点 temp = k;???????%???当前插入顶点 D(1:n) = inf;???%???初始化最小权向量 D(k) = 0;???????%???顶点?V(k)?到本身的权为?0 I(1:n) = k;?????%???初始化最短路向量 while?sum(V) 0??????????%???尚存在可作为终点的顶点，则继续迭代 ????v = find(V~=0);???????%???v?是可作为终点的顶点向量，用作循环变量的取值范围 ????for?i = v?????????????%???分别以?v?中所有的顶点作为终点，插入一个顶点?temp ????????if?D(i) D(temp)+W(temp,i)?????%???判断是否比之前更优 ????????????D(i) = D(temp)+W(temp,i);???%???保存更小权，更新最小权向量 ????????????I(i) = temp;????????????????%???修改当前路径的终点的前驱 ????????end ????end ? ??%???在可作为终点的顶点向量?v?中找到当前使?D?最小的顶点 ????curmin = find(D==min(D(v))); ????for?j = 1:length(curmin)?????%???因为以上找到的顶点可能有多个，这里只取一个 ????????if(any(curmin(j)==v)) %??不能直接用?temp=curmin(1)?或者?temp=curmin(2)?之类的表达 %??因为?curmin(1)等可能只是?D?却不是?v?的一个合法索引！ ????????????temp = curmin(j); ????????????break;???????????????%???将当前使?D?最小的顶点作为下次迭代的插入顶点 ????????end ????end ????V(temp) = 0;?????????????????%???将当前使?D?最小的顶点标记为不可再作为终点 end