1、问题分三个步骤完成，第一步对整体业务员三月份订单进行分析，获得业务员订单的变化。第二步基于保证业务员订单数量的完成，达到工厂利益最大同时机器不空转。第三步对十名业务员三月份订单数量进行方差模拟，得到最优解。2.2 问题二的分析问题二在保证所有订单都能如期交货的情况下，采用随机数法，在大量随机数据下，寻得最优排产使得生产费用最低。问题分三个步骤完成，第一步，将订单进行排序，量化时间，设出每台机器对不同订单的贡献时段。第二步，用随机数法，产生足够数量且符合约束条件的每台机器对不同订单的贡献时段。第三步，综合第二步所生成的结果，以生产费用最低为标准，得出合理的排产计划。三 模型假设1. 工厂生产机器
2、可 24 小时运转且不出现故障，不产生废料。2. 客户下单日机器即开始工作，直至交货日期。3. 机器转化生产型号时无需耗时。4.上个月订单全部完成，不拖欠到本月。四符号说明第 i 项订单开始的时刻（即接单时间）第 i 项订单结束的时刻（即交货时间）第 j 台机器开始生产第 i 项订单的时刻第 j 台机器结束生产第 i 项订单的时刻判断第 j 台机器从第 i - 1 项订单转为第 i 项订单是否对第记录第 j 台机器从第 i - 1 项订单转为第 i 项订单的转化各型号漆包线订单总吨数第 i 项订单对应的铜漆包线类型i订单号五模型的建立与求解5.1问题 1 的模型建立与求解通过细化额度,我们现在
3、三月份将各个型号的漆包线（0.85mm，0.90mm，0.95mm，1.00mm，1.32mm，2.00mm，2.05mm，2.65mm，3.35mm，3.75mm）进行分析，把业务员的订单情况汇总为表格：（详见附录1）表5.1.1 ：每个业务员对不同型号的订单数量及总量业务员（人）数量（吨）（吨）型号（mm）0.850.90.951.001.322.02.52.653.353.75总计11080502023055256031004511012067708909380260445210265290195335在三月份，如果机器不空转的情况下，最多能生产2700吨，但是三月份业务员的订单总额却是
4、2900吨，根据附件可得超出订单量可在4月份生产，不影响本月的销售计划。我们仅对三月份的销售员的订单情况通过MATLAB计算方差运算，程序如下：（详见附录2）a=250 220 320 200 420 440 260 380 210 200;b=sum（a）/size（a,1）; c=0;for i=1:10 c=（a（i）-b）2+c;end d=sqrt（c/（size（a,1）-1）;disp（d）;最终我们得到业务员在每月275吨左右订量时，方差最小，即收益最高。同时也符合工厂的最优利益。5.2问题 2 的模型建立与求解我们已知该工厂车间拥有 20 台机器，由于机器不能在漆包线规格0.
5、85mm 1.32mm与规格2.0mm 3.75mm之间互相转化。因此，我们将规格种类0.85mm 1.32mm与2.0mm 3.75mm分别记为甲品种漆包线与乙品种漆包线。现在我们对两种类型的漆包线订单量进行讨论，通过对附录表格的数据进行计算，我们可以轻易地得出，甲类型与乙类型漆包线在3月份的总订单量分别为 1555t与1335t. 经计算，9台机器在 36 天（工厂订单共涵盖36天）的产量为 1458t，10台机器 36 天产量为1620t，11台机器 36 天产量为1782t. 由此，我们可以进行合理假设，工厂各指派 10 台机器分别加工甲类型漆包线与乙类型漆包线.并且可以知道每台机器每
6、小时的工作量为 0.1875t.。对甲类型订单的分析对附录表格进行分析整理，将所有订单按照接单日期的先后顺序，将天数小时化（3月 2 号到 4 月 6 号细分为 0 到 864 小时），并记 3 月 2 日 0 点为时间轴起点，我们可以得到如下时间图表 （图 5.2.1），并分析得到订单属性 （图 5.2.2） 。图5.2.1横坐标为时间轴，纵坐标为订单序号图 5.2.2我们用表示第j台机器生产第i项订单的时刻，用第i台机器结束生产第i 项订单的时刻。由订单属性可得不等式：其中 i = 1, 2, , 24 ，j = 1, 2, , 10 .由于机器只有结束第 i 项订单的生产后，才能转向生产
7、第 （i + 1） 项订单，可得方程：问题二要求每项订单均能按时完成，并最小化费用，而（-）表示第j台机器对第 i 项订单贡献的小时数，故得方程：其中 i = 1, 2, , 24.工厂生产费用考虑三个方面：生产类型转化的费用、总空转费用、总加工费.由于生产总量固定，故而总加工费为一常数为 15378 元，且在 36 天内的总空转天数数也为常数（1620-1555）/4.5=14.45，所以空转费用为 34667元.因此只需考虑生产类型转化的费用。因此可以定义用以判断第j 台机器从第 i - 1 项订单转为第 i 项订单后，是否对第 i 项订单有贡献. 定义用以判断第j 台机器从第 i - 1
8、 项订单转为第 i 项订单后，是否需要转化费. 故而，其中 i = 2,3, , 24 ，j = 1, 2, , 10 .于是可以得到第 j 台机器从第 i - 1 项订单转为第 i 项订单后所需要的转化费为：所以总转化费用为：此即为目标函数，我们希望得到它的最小值。采用计算机随机生成 10000 组数据进行模拟（详见附录 3 ），选取其中一组代表性数据置于（附录4）经随机模拟后得到的甲类型大致最小总转化费为 88300 元,将空转费以及加工费计算进去得到生产最小总费用为88300+15378+34667=138345（元）。对乙类型订单的分析思路。方法同甲类型，这里不再赘述。只将原始数据分析
9、图表图 5.2.2横坐标为时间轴，纵坐标为订单序号图 5.2.3经计算，B 类型生产最小总费用为 124600+14388+152000=290988（元）综上，按照附录中的排产方式，得到最小总费用为 138345+290988=429333（元）六模型评价与改进6.1 模型评价（1）本文基于对业务员工作量的规划以及工厂订单的合理安排，采取了方差模拟，MATLAB运算等，表现出了较强的数据处理能力，可以在现实情况使用对应求解方法使得工厂与业务员的利益最大化与平衡化。（2）充分结合现实生产经验与合理假设，提出了如随机模拟、多目标方程和最优问题求解等思想。6.2 模型改进（1）计算复杂度的优化由于
10、前面所采用的算法都基于方差模拟、随机生成数的思想，当工厂订单量递增时，计算复杂度将呈现指数增长。对此，我们可以采取更加优秀的搜索算法，例如使用分治算法进行二分查找等。（2）准确度的优化由于所采用方差模拟与随机数法均存在波动性，所得数据精确度可进一步提高。对此，我们可以采取深度搜索算法，使之得到更为确切的排产方式，并且可以为之为参考，在其附近寻求更为精确的排产方式。七参考文献1 韩中庚，数学建模方法及其应用，北京: 高等教育出版社, 20092 卓金武，MATLAB在数学建模中的应用，北京：北京航空航天出版社 2011附录1：附录2：附录3#include cstdiocmathcstdlibc
11、timeusing namespace std;int t2511, k2511;int u2511, v2511;int temp125 =0,0,0,24,48,48,72,72,96,144,168,192,192,192,312,336,336,360,360,408,432,504,528,552,624;int temp225 =0,336,216,336,360,312,432,360,432,504,408,552,456,456,672,600,576,696,672,720,672,720,768,864,840;double m25 =0,10,30,30,60,30,1
12、10,60,70,55,60,80,20,100,100,80,80,50,70,80,60,60,50,120,90;double s25 =0,0.85,0.9,1,0.95,1.32,0.85,0.95,0.9,0.95,0.85,0.9,1,1.32,0.95,0.85,1,0.9,1.32,1,0.9,1.32,0.95,0.95,0.85;/随机生成数据double random（double start,double end）return start + （end - start） * rand（） / （RAND_MAX + 1.0）;int main（） /随机生成t int
13、 x;srand（ unsigned（ time（0） ） ）;/srand（time（NULL）;for（int i = 1; i = 24; i+） ui1 = temp1i;vi1 = temp2i;for（int j = 1; j = mi / 0.1875） flag = 1;else flag = 0;continue;, kij）;/生成p , 生成q int p2511, q2511, min = 0;for（int i = 2;if（kij = tij） qij = 0;else qij = 1;if（si = si - 1） pij = 0;else pij = 400;min += qij + pij;%dn, min）;return 0;附录4：
