《现代优化技术-靳志宏》DU-现代优化技术.docxVIP

现代优化技术 现代优化技术 期末作业 交通运输规划与管理 杜剑 1120140359  作业1：Dijkstra最短路径算法—Matlab实现 问题描述 Dijkstra最短路算法用于寻找路网中两点间权数最小的路径，最短路问题的实际应用包括：设备更新、工厂布局、运输路径选择、线路安排等。 Matlab是处理矩阵运算的有效工具，此处采用Matlab进行求解。 Matlab编码 数值输入与参数设计： 输入路网矩阵 G，计算路网中节点数 N=size(G,1)。 设置v1为临时标号阵列，v0为永久标号阵列，当某节点转为永久编号时其数值由1变为0，节点加入v0。编码：v1=ones(1,N); v0=1; v1(v0)=0; 输出为： v0 = 1 v1 = 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 end end 外层中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; end end [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实现 clear clc N=size(G,1); v1=ones(1,N); v1(v0)=0; D=G(v0,:); while 1 DT=D; for k=1:N if v1(k)==0 DT(k)=inf; end end [Dmin id]=min(DT); if isinf(Dmin),error,end v0=[v0 id]; v1(id)=0; if size(v0,2)==N,break;end id=0; for j=1:N if v1(j) for i=1:N if v1