《现代优化技术-靳志宏》作业1-最短路-DU.docxVIP

作业1：Dijkstra最短路径算法—Matlab实现1120140359 杜剑交通运输规划与管理问题描述Dijkstra最短路算法用于寻找路网中两点间权数最小的路径，最短路问题的实际应用包括：设备更新、工厂布局、运输路径选择、线路安排等。Matlab是处理矩阵运算的有效工具，此处采用Matlab进行求解。Matlab编码数值输入与参数设计：输入路网矩阵G，计算路网中节点数N=size(G,1)。设置v1为临时标号阵列，v0为永久标号阵列，当某节点转为永久编号时其数值由1变为0，节点加入v0。编码：v1=ones(1,N);v0=1; v1(v0)=0;输出为：v0 = 1v1 =0 1 1 1 1 1 1 1算法终止条件：算法循环使用while语句，考虑路网不存在最短路、及最多循环N次两个终止情况，设置终止条件，编码如下：if size(v0,2)==N,break;end跳出while循环if isinf(Dmin),error,end报错for j=1:N %设置两层循环，加两条判断语句。 if v1(j) for i=1:N if v1(i)==0 D(j)=min(D(j),D(i)+G(i,j)); end D(j)=min(D(j),G(v0(1),i)+G(i,j)); end endend外层中if v1(j)，表明j为临时编号点时进入第二层循环，为永久编号点时j=j+1。内层循环是在j为临时编号点条件下，确定通过其他各个点到达j点的最短路D(j)。if v1(i)==0 D(j)=min(D(j),D(i)+G(i,j)); end D(j)=min(D(j),G(v0(1),i)+G(i,j));其中，D(j)、D(i)指起点到达j\i点的最短路长度；G(i，j)为路网中第i行j列，即由i点到j点的直线距离；v0为永久编号阵列；v1为临时编号阵列，if v1(i)==0 用于判断i点是否为永久编号点。for k=1：N，通过求 min D(k)，确定下一个永久编号点。D=G(v0,:); DT=D; for i=1:N if v1(i)==0 %对i是否为永久编号做判断 DT(i)=inf; endend [Dmin id]=min(DT);v0=[v0 id];记录最短路径长度，修改永久编号阵列。算例输入参数G，v0(1)如下：G=[inf 2 8 1 inf inf inf inf; 2 inf 6 3 1 inf inf inf; 8 6 inf 7 5 1 2 inf; 1 3 7 inf inf inf 9 inf; inf 1 5 inf inf 3 inf 8; inf inf 1 inf 3 inf 4 6; inf inf 2 9 inf 4 inf 3; inf inf inf inf 8 6 3 inf]; v0=1; 运行程序：输出结果：最短路径:v0 = 1 4 2 5 6 3 7 8路径长度:D = Inf 2 7 1 3 6 9 11附录%% Dijskra问题的Matlab实现clearclcN=size(G,1);v1=ones(1,N); v1(v0)=0;D=G(v0,:);while 1 DT=D;fork=1:Nif v1(k)==0 DT(k)=inf; endend[Dmin id]=min(DT);if isinf(Dmin),error,endv0=[v0 id];v1(id)=0; if size(v0,2)==N,break;end id=0;for j=1:N if v1(j) for i=1:Nif v1(i)==0 D(j)=min(D(j),D(i)+G(i,j));end D(j)=min(D(j),G(v0(1),i)+G(i,j));endendendif isinf(Dmin),error,endenddisp(最短路径:)v0disp(路径长度:)D