实战之利用穷举迭代法处理动态规划生产与储备问题

问题如下：

生产与存贮问题

    一个生产项目，在一定时期内,增大生产量可以降低成本费,但如果超过市场的需求量,就会因积压增加存贮费而造成损失。相反,如果减少生产量,虽然可以降低存贮费,但又会增加生产的成本费,同样会造成损失。因此,如何正确地制定生产计划,使得在一定时期内,生产的成本费与库存费之和最小,这是厂家最关心的优化指标，这就是生产与存贮问题。

    假设某车间每月底都要供应总装车间一定数量的部件。但由于生产条件的变化,该车间每月生产单位部件所耗费的工时不同,每月的生产量除供本月需要外,剩余部分可存入仓库备用。今已知半年内,各月份的需求量及生产该部件每单位数所需工时数如下所示:

    月份( k)：      1 2 3  4 5 6

    月需求量(bk)：   8  5  3 2 7 4

    单位工时(ak)：   11   18    13    17    20    10

    设库存容量H = 9,开始时库存量为2,期终库存量为0。要求制定一个半年逐月生产计划,使得既满足需求和库存容量的限制,又使得总耗费工时数最少。      

题目解读：进行生产与储备的规划，叫我们满足最大库存容量的条件下去求最小工时数。

穷举迭代法

方法原理：找出所有可能的生产情况，以及耗费的工时数，再去找出最小工时数下，对应的最优生产计划。

具体实现代码（matlab）：

% 需求量和单位工时数
b = [8, 5, 3, 2, 7, 4];
a = [11, 18, 13, 17, 20, 10];

% 库存容量和期初库存量
H = 9;     %最大存量
s = 2;     
s1 = 2;    %初始库存
s2 = 8-2;  %减去初始库存后的需求量
c1 = [];   %装取所有可能的生产部件数
c2 = [];   %装取所有可能的生产部件数下对应的工时数

for i = 8-2:8+9-2
    %剩余库存量
    l1 = 2+i-8;
    %第二层c2
    for h = 5-l1:5+9-l1
        l2 = l1+h-5;
        %第三层c3
        for u = 3-l2:3+9-l2
            l3 = l2+u-3;
            %第四层c4
            for k = 2-l3:2+9-l3
                l4 = l3+k-2;
                %第五层c5
                for o = 7-l4:7+9-l4
                    l5 = l4+o-7;
                    %第六层c6
                    if 4-l5 >= 0
                        g = 4-l5;
                        if i >= 0 && h >= 0 && u >= 0 && k >= 0 && o >= 0 && g >= 0
                            fprintf('生产计划为：%d\n',i,h,u,k,o,g);
                            c1 = [c1,i];
                            c1 = [c1,h];
                            c1 = [c1,u];
                            c1 = [c1,k];
                            c1 = [c1,o];
                            c1 = [c1,g];
                            wt = 11*i+18*h+13*u+17*k+20*o+10*g;
                            c2 = [c2,wt];
                            fprintf('总耗费工时数：%d\n',wt)
                            if wt == 309   %经过所有迭代后取得的c2最小值
                                fprintf('最优生产计划：%d\n',i,h,u,k,o,g);
                                fprintf('最小耗费工时数：%d\n',wt);
                            else
                                break 
                            end   
                        else
                            break
                        end
                    else
                        break

                    end
                end
            end
        end
    end
end



          
ssr = length(c1)/6;
fprintf('总的生产可能情况种数为：%d\n',ssr);

%生产计划中各个月的生产件数
x1 = [];
x2 = [];
x3 = [];
x4 = [];
x5 = [];
x6 = [];

x1 = [x1,c1(1:6:end)];
x2 = [x2,c1(2:6:end)];
x3 = [x3,c1(3:6:end)];
x4 = [x4,c1(4:6:end)];
x5 = [x5,c1(5:6:end)];
x6 = [x6,c1(6:6:end)];

%disp(length(x1));
%disp(length(c4));

%可视化第一个月与最小工时的关系图
scatter(x1,c4);
   

 

mathlab求解结果：

最优生产计划：15（第一个月）
最优生产计划：0（第二个月）
最优生产计划：8（第三个月）
最优生产计划：0（第四个月）
最优生产计划：0（第五个月）
最优生产计划：4（第六个月）
最小耗费工时数：309
总的生产可能情况种数为：3101

第一个月的所有生产情况与对应的总耗费工时的可视化关系图

因此，动态规划方法下，求得此题的解为，最优生产计划，第一个月生产15件，第三个月生产8件，第六个月生产4件，其余月份生产0件，最小工时为309。