数据建模 - Floyd算法求最短路径

    Floyd算法又称为弗洛伊德算法、插点法,是一种利用动态规划的思想寻找给定的加权图中多源点之间最短路径的算法,与Dijkstra算法类似。该算法名称以创始人之一、1978年图灵奖获得者、斯坦福大学计算机科学系教授罗伯特·弗洛伊德命名。

课题名称:设备更新,使总的支付费用最少

    工厂的某台机器可连续工作四年,决策者每年年初都要决定机器是否需要更新。若更新,就要支付购置费用;若不更新,则要支付维修与运行费用,且随着机器使用年限的增加维修与运行费用逐年增多。计划期(4年)中每年年初的购置价格及各个年限内维修与运行费用由表1给出,

(1)试制订今后4年的机器更新计划,使总的支付费用最少。

表1

第 i 年初

第1年初

第2年初

第3年初

第4年初

购置费(万元)

扫描二维码关注公众号,回复: 2342760 查看本文章

2.5

2.6

2.8

3.1

使用年限

1

2

3

4

使用年限下对应的每年的维修与运行费(万元)

1

1.5

2

4

可将此问题看成一个最短路问题。设 v1 和 v5 分别表示计划期的始点和终点(v5 可理解为第4年年末)。下图中各边的值(vi,vj)表示在第 i 年初购进的机器使用到第 j 年初(即第 j-1 年底),(vi,vj)的值可由表1的数据计算得到。因此,把求最优设备更新问题转化为求从 i 到 j 的最短路问题。

(vi,vj)的值计算示例:

(v1,v2) = 2.5(第1年初的购置费)+1(使用第一年的维修费) = 3.5;

(v3,v5) = 2.8(第3年初的购置费)+1(使用第一年的维修费)+1.5(使用第二年的维修费) = 5.3;

             


Matlab求解

Matlab程序:

a=[0 3.5 5 7 11;inf 0 3.6 5.1 7.1;inf inf 0 3.8 5.3;inf inf inf 0 4.1;inf inf inf inf 0];
n=size(a,1);
D=a;
path=zeros(n,n);
for i=1:n
    for j=1:n
        if D(i,j)~=inf
            path(i,j)=j;
        end
    end
end
for k=1:n
    for i=1:n
        for j=1:n
            if D(i,k)+D(k,j)<D(i,j)
                D(i,j)=D(i,k)+D(k,j);
                path(i,j)=path(i,k);
            end
        end
    end
end
D;path;
Matlab运行结果:
D =
         0    3.5000    5.0000    7.0000   10.3000
       Inf         0    3.6000    5.1000    7.1000
       Inf       Inf         0    3.8000    5.3000
       Inf       Inf       Inf         0    4.1000
       Inf       Inf       Inf       Inf         0
>> path
path =
     1     2     3     4     3
     0     2     3     4     5
     0     0     3     4     5
     0     0     0     4     5
     0     0     0     0     5
解得方案:最短路(v1,v3,v5),即计划期内机器更新最优计划为第1年、第3年初各购进一台新机器,4年总的支付费用为10.3万元。


新问题来了!

按照上述方案,第3年初购置新机器,那么旧机器是否可以低价处理掉呢?这样也更加符合实际。如果已知不同役龄机器年末的处理价格如表2所示,那么在这计划期内机器的最优更新计划又会怎样?

表2

年度

第一年末

    第二年末

第三年末

第四年末

机器处理价(万元)

2.0

1.6

1.3

1.1

类似上面的处理方式,可转化为求下图中从 v1 到 v5 的最短路问题。

       

Matlab求解

Matlab程序:

a=[0 1.5 3.4 5.7 9.5;inf 0 1.6 3.5 5.8;inf inf 0 1.8 3.7;inf inf inf 0 3.1;inf inf inf inf 0];
n=size(a,1);
D=a;
path=zeros(n,n);
for i=1:n
    for j=1:n
        if D(i,j)~=inf
            path(i,j)=j;
        end
    end
end
for k=1:n
    for i=1:n
        for j=1:n
            if D(i,k)+D(k,j)<D(i,j)
                D(i,j)=D(i,k)+D(k,j);
                path(i,j)=path(i,k);
            end
        end
    end
end
D;path;
Matlab运行结果:
D =
         0    1.5000    3.1000    4.9000    6.8000
       Inf         0    1.6000    3.4000    5.3000
       Inf       Inf         0    1.8000    3.7000
       Inf       Inf       Inf         0    3.1000
       Inf       Inf       Inf       Inf         0
>> path
path =
     1     2     2     2     2
     0     2     3     3     3
     0     0     3     4     5
     0     0     0     4     5
     0     0     0     0     5
即最后解得:最短路 (v1,v2,v3,v5),即计划期内机器更新最优计划为第1年初购进一台新机器并在第2年初处理掉,第2年初购进一台新机器并在第3年初处理掉,第3年初购进一台新机器,4年总的支付费用为6.8万元。


内容很多,消化一下吧,希望可以为你所用!


猜你喜欢

转载自blog.csdn.net/weixin_40683253/article/details/80985599
今日推荐