A陶子-数学模型-食品安全的抽检问题论文2

1、最新A陶子-数学模型-食品平安的抽检问题论文2食品平安的抽检问题摘要 食品的质量和卫生问题是关系到民生的大问题，因此，对食品的检查显得非常重要。本文结合实际，应用AHP方法、分层抽样和线性目标规化方法，建立了集时间、费用和效果为一体的数学模型，具体如下。 对于问题一，我们首先将主要食品进行分类，然后将影响食品平安的因素主要分为生物性污染、化学性污染、物理性污染三大类，并将这三类污染所造成的主要危害归纳为七类，接着采用AHP法对问题进行定量分析，最后通过一致性检验并得出其危害性的大小，得到结果细菌危害最严重，食品添加剂导致的危害次之等。 对于问题二，针对局部主要产品，我们先采用了分层抽样的方法对

2、不同品牌不同批次的产品进行抽检，建立了样本分配率、样本方差、总体抽样率、分层抽样率等函数方程，然后对上一步所抽到的批次利用线性目标规划的方法，建立了集时间较短、本钱费用较低和抽样效果较好的抽检模型线性目标规划模型，并利用统计学原理对检测误差进行分析。最后，我们根据模型针对乳制品中的酸奶进行模拟检验，检验的结果误差百分比为4.24%5%，可靠性较高。对于问题三，我们利用问题二所建立的模型制订了一种较为合理的抽检方案根据假设总共抽检79个批次，每个批次抽检2个工程。然后，我们进行了可靠性分析，抽检的误差百分比为1.15%5%，可靠性较高。对于问题四，它实际是在问题三的根底上，对面粉进行屡次跟踪抽检

3、。我们对问题二所建立的模型进行了改良，引入新的变量建立函数关系，并运用MATLAB优化工具箱进行求解，得出了最正确的抽检策略和抽检数量结果为跟踪抽检3次，共抽检113个批次，所得结果可靠性较高、本钱较低，且工时比拟少，用计算机进行模拟检验时效果比拟乐观。最后，我们对模型的优缺点进行了评价，讨论了其推广应用的价值，并主管部门写了一份报告，提出了一些解决问题的可行性建议，可为主管部门和市民提供一些参考。关键词 ：AHP法，分层抽样，目标规划、统计分析、可靠性 一 问题的重述改革开放三十年来，我国人民生活水平在不断地提高，食品平安和卫生问题越来越受到人们的关注。近几年来，先后出现了苏丹红、瘦肉精、三

4、聚氰胺等事件，以及各种不利于健康的食品添加剂、强化剂问题的出现，食品平安和卫生的检测已成为全社会，乃至政府有关部门重点关注的问题之一。食品的质量和卫生问题涉及到原材料的使用、生产加工、运输与贮存、流通与销售等环节，在每一个环节上出现过失，都会导致食品出现平安和卫生问题，食品质量和卫生的检测工作在实际显得非常重要。但是，由于食品的种类、品牌和批次繁多，从生产加工到销售食用中间环节复杂，质检部门不可能对所有食品做到全面的质量检测，一般只能做一定的抽检。当然，对食品进行抽检也需要一定人力、物力和财力即本钱费用，抽检的越多检测效果就越好，但需要的时间就越长，其本钱费用也就越高。为此，应该如何抽检，既能

5、保证较好的检测效果，又能节省时间和本钱费用？请根据实际情况，建立数学模型分析研究以下问题：1根据主要食品的分类，试分析影响各类食品平安的因素，对其可能的危害性做出定量的比拟评估分析。2针对局部主要食品，结合实际建立合理的抽检模型，给出检测误差即检测的可靠性分析的方法，并对模型进行模拟检验。3面粉是我国中北部地区主要的主食原材料，不妨就已推广食用的“营养强化面粉抽检问题进行讨论，“营养强化面粉的配方标准如表1所示。假设某地区现有12个品牌的营养强化面粉产品，每个品牌每月将有不少于60个批次即同一企业、同一条生产线、同一批投料、同一班次生产的产品为1个批次的产品在市场上销售，质检部门要做一次全面的

6、质量检查，请你帮助制订一种合理的抽检方案，并分析其检测的可靠性。所需要检验工程、标准、本钱和工时如表2所示。4针对问题3，如果质检部门需要连续进行屡次跟踪抽检，请你给出相应的抽检策略和最正确的抽检数量，使其检测可靠性尽量高、本钱尽量低、工时尽量少，并用计算机进行模拟检验。5请根据你们的研究，深入分析食品平安存在的隐患和根源，并提出有效可行的解决问题方法和建议，可供主管部门和市民参考。二 模型的假设1、假设模型求解过程中所用的数据都是合理的；2、假设同一类食品每次抽检的工程数相同；3、假设不考虑抽检人数的限制；4、假设每天抽检的最大时间为法定工作时间8小时，一个月工作22天；5、假设抽检的最大费

7、用有限制；三 符号约定-总体均值的估计量st表示分层；-层权,i=1,2N;L-抽样系统中总的生产批次；-第i层生产批次数；-第i层的样本均值；-第i层抽取的样本量；-第i层的方差；-第i层的抽样比/；-为L在1-置信水平下的绝对误差限；-标准整天分布的双侧分位数；-第i层每个样本单位的平均费用；-第i层每个样本单位的平均检验工时；-固定费用；-每个抽检批次的检查工程；-抽检的最长时间；-抽检的最大费用；-抽检的总工程数；-每次每项抽检的平均时间；-每次每项抽检的平均费用；四 问题的分析与模型的建立4.1对问题一的分析与模型建立问题一的分析为了分析与比拟影响食品平安因素所产生的危害程度，首先我

8、们将主要食品分为：肉制品、乳制品、酿造食品和蛋白质含量较高食品等；接着分析确定影响食品平安的因素，我们从微生物污染、化学污染和物理污染三方面进行分析与确定；然后对于各因素影响食品产生的危害，其一是由微生物污染导致的细菌危害、病毒危害和寄生虫危害；其二是由化学污染引起的食品添加剂如非食用添加剂导致的危害、环境污染如农药，包括灭鼠药和兽医用药导致的化学危害和天然存在的化学危害；其三是由物理性污染引起的重金属中毒等危害。因此，问题一就可以归结为各因素影响食品导致的危害对食品质量的影响权重问题。解决这类问题首先要统计分析各类影响指标的数值特征。然后再对其进行归一化处理，并利用层次分析法对其权重进行赋值

9、。最终得到各因素对食品平安的影响排名。问题一的模型建立1根据影响食品平安的各因素与危害间的关系，我们建立的不完全层次分析算法的根本结构图如下： 2根据对问题一的分析，我们建立了如下成比照拟矩阵假设要比拟某一层的个因素对上一个因素O的影响，取两个因素和，用表示和对的影响之比，全部比拟的结果可用成比照拟矩阵：，表示，并称以上矩阵为正互反矩阵。3一致性检验及权重向量确实定成比照拟矩阵通常不是一致阵，为了能用它的最大特征根的特征向量作为被比拟因素的权向量，其不一致程度需要控制在一定的范围内。计算方法如下：定义：一致性指标；一致性比率为;本文中，对应的.当时，认为矩阵的不一致程度在容许的范围内。其中最大

10、特征根的特征向量即为相应的权重向量。4.2问题二的分析与模型建立问题二的分析从问题二的题设我们知道，为了建立合理的抽检模型，我们首先需要对现有的统计信息进行定量分析。由于问题二的求解几乎需要考虑文中提供的所有信息，所以为了简化数学模型的计算复杂度，我们可以在某些直观问题上进行简单的定性分析。例如，原文附录中给出了主要食品和主要抽检工程，为此我们可以主观上认为我们要考虑的局部食品就是附录1中的局部主要食品，要抽检的工程就是附录2中给的8个工程。在根本确定主要食品和抽检工程后，首先我们利用分层抽样模型确定这些主要食品中每种品牌所要抽检的批次数的最优值，然后，我们利用法定的工作时间确定最大抽检时间和

11、费用。最后计算出每种品牌中要抽的批次数及每个批次的药抽检的工程数，对其优劣进行综合评价。问题二的模型一分层抽样通过上述对问题二的分析，我们利用分层抽样法建立以下函数关系： 设整个抽样为一个系统L，系统中含有种品牌，每个品牌的生产批次数为，那么总体均值的估计量： 1第i层的方差： 2由误差知识可得： 3联立上述各式可得： 4调查费用函数： 5各层样本量的最优分配： 6由于每个样本的调查费用相同，故上式可以化为 7联立4和7可得，总体抽样率为 8也可得到分层抽样率为 9二模型的建立线性目标规化在以上函数关系以及抽检费用有限制的条件下，模型以抽检工程数最多为目标函数所建立的模型如下： S.T. 10

12、2检测误差即检测的可靠性分析的方法: 114.3问题三的分析面粉是一种十分重要的原材料，因此对面份的抽检具有代表性。根据题目的要求，质检部门要做一次全面的质量调查，其抽检方案是问题二模型的具体应用。所以，我们只需将问题二建立的模型具体化，得出具体的抽检方案，并对检测的可靠性进行分析即可。4.4问题四的分析与模型的建立问题四的分析对于问题四，它其实是在问题三的根底上，对面粉的质量进行屡次跟踪抽检，我们对问题二所建立的模型进行了改良，引入新的变量时间，建立新的函数关系，并运用MATLAB优化工具箱进行求解，得出了最正确的抽检策略和抽检数量，使检测可靠性较高、本钱较低，且工时比拟少，用计算机进行模拟

13、检验时效果比拟乐观。问题四的模型建立由屡次跟踪抽检可知：抽检批次数与时间是成反相关的，故我们可以作L与t的函数：,分别使抽检的效果最好、费用最低、时间最少。目标函数为：五 模型的求解5.1模型一的求解5.1.1确定准那么层对目标层的权重向量A=,=3.0385，权重向量为w=0.6370，0.2583，0.10475.1.2确定方案层对准那么层的权重向量=,= ，=1.求解结果如下：编程程序见附件一矩阵特征值一致性指标CI随机一致性指标RI一致性比率CR权重向量3.03850.019250.580.0330.10.6370，0.2583，0.10473.06490.032450.580.016

14、2250.10.6491，0.2790，0.071910000.11由上表中0.1可知矩阵、都通过了一致性检验。5.1.3确定影响因素的权重向量组合权重：(i)=(0.6370*0.6370,0.2583*0.6370,0.1047*0.6370,0.6491*0.2583,0.2790*0.2583,0.0719*0.2583,1*0.1047)=(0.405769,0.1645371,0.0666939, 0.16766253, 0.0720657, 0.01857177,0.1047)=0.4058，0.1645，0.0667，0.1677，0.0721，0.0186，0.1047CR=

15、(0.01925*0.6370+0.03245*0.2583+0*0.1047)/(0.58*0.6370+0.58*0.2583+0*0.1047)=0.039760.1,通过了组合一致性检验。5.1.4排名由上述组合权重可知，各因素对食品平安的影响程度由强到弱依次为0.40580.16770.16450.10470.07210.06670.0186细菌污染食品添加剂导致的危害病毒危害重金属中毒等环境污染导致的危害寄生虫危害天然存在的化学危害5.2问题二的求解层次抽样模型:主要食品包括肉制品、乳制品、酿造食品及高蛋白食品等，针对局部主要食品进行模拟抽检。首先针对局部主要食品乳制品中的酸奶讨论

16、，进行分层。第一层为：蒙牛、伊利、光明、百万庄园、三元、娃哈哈六个品牌的批次数，然后对每层随机抽样，得到样本。分层抽样图如下：局部参数的假设：各品牌生产批次序列为：置信度取95%，根据对其他参数的讨论，按照ISO质量管理体系和AQL抽样标准得到实际中抽检概率在8%10%之间能够到达抽检效果，标准差取值范围在0.0370.042之间，并在此区间内等距选取六个数进行计算，六个数分别为，d取0.01。 由标准正态分布表查表可知：=1.64。其他参数的计算：N = =440 批 所以由有： ； 分层抽样率的计算由9计算得到各层的抽样率如下表：酸奶品牌蒙牛伊利光明百万庄园三元哇哈哈抽样率0.08130.

17、08350.08570.08790.09010.0923假设t=1.5,p=60,Pmax=6000，Tmax=176, =11工程数为2.56413项。误差：=440/6=73.3=10.8013=1.6530误差百分比：1.6530/39=4.24%5%权重=0.2175，0.2175，0.0725，0.0725，0.0435，0.0435，0.0362，0.0311，0.2175，0.0242，0.0242排序：三聚氰胺,维生素A,维生素B11，维生素B1、维生素B2、铁、锌、维生素B5、钙、水分、蛋白精、硫氰酸钠随及抽检三个工程，前三项抽检概率最大，分别为三聚氰胺,维生素A,维生素B1

18、1.5.3 问题三的求解与结果分析 问题三的求解局部参数的假设：各品牌生产批次序列为：置信度取95%，d取0.01。根据对其他参数的讨论，按照ISO质量管理体系和AQL抽样标准得到实际中抽检概率在8%10%之间能够到达抽检效果，标准差取值范围在0.03250.0365之间，并在此区间内选取12个数进行计算，12个数分别为5其他参数的计算：N = =880批所以由有： ； 用MATLAB对函数进行优化可得：工程数为1.01272项。误差：=880/12=73.3=8.4782=0.9106误差百分比：0.9106/79=1.15% A=1 3 5;1/3 1 3;1/5 1/3 1A = 1.0

19、000 3.0000 5.0000 0.3333 1.0000 3.0000 0.2000 0.3333 1.0000 x,y=eig(A)x = 0.9161 0.9161 0.9161 0.3715 -0.1857 + 0.3217i -0.1857 - 0.3217i 0.1506 -0.0753 - 0.1304i -0.0753 + 0.1304iy = 3.0385 0 0 0 -0.0193 + 0.3415i 0 0 0 -0.0193 - 0.3415i w=x(:,1)/sum(x(:,1)w = 0.6370 0.25830.1047=3.0385，权重向量为w=0.63

20、70，0.2583，0.1047=(3.0385-3)/(3-1)=0.01925RI=0.58=0.01925/0.58=0.033 B1=1 3 5;1/3 1 3;1/5 1/3 1B1= 1.0000 3.0000 5.0000 0.3333 1.0000 3.0000 0.2000 0.3333 1.0000 x,y=eig(B1)x = 0.9161 0.9161 0.9161 0.3715 -0.1857 + 0.3217i -0.1857 - 0.3217i 0.1506 -0.0753 - 0.1304i -0.0753 + 0.1304iy = 3.0385 0 0 0 -

21、0.0193 + 0.3415i 0 0 0 -0.0193 - 0.3415i w=x(:,1)/sum(x(:,1)w = 0.6370 0.25830.10471=3.0385，权重向量w1=0.6370，0.2583，0.1047 B2=1 3 7;1/3 1 5;1/7 1/5 1B2 = 1.0000 3.0000 7.0000 0.3333 1.0000 5.0000 0.1429 0.2000 1.0000 x,y=eig(B2)x = 0.9140 0.9140 0.9140 0.3928 -0.1964 + 0.3402i -0.1964 - 0.3402i 0.1013

22、-0.0506 - 0.0877i -0.0506 + 0.0877iy = 3.0649 0 0 0 -0.0324 + 0.4448i 0 0 0 -0.0324 - 0.4448i w2=x(:,1)/sum(x(:,1)w2 = 0.6491 0.27900.07192=3.0649，权重向量w2=0.6491，0.2790，0.0719附件2： f=0.037*(9/44*0.037+2/11*0.038+3/22*0.039+15/88*0.04+13/88*0.041+7/44*0.042)/(9/44*0.0372+2/11*0.0382+3/22*0.0392+15/88*0

23、.042+13/88*0.0412+7/44*0.0422)+440*(0.01/1.64)2)f = 0.0813 f=0.038*(9/44*0.037+2/11*0.038+3/22*0.039+15/88*0.04+13/88*0.041+7/44*0.042)/(9/44*0.0372+2/11*0.0382+3/22*0.0392+15/88*0.042+13/88*0.0412+7/44*0.0422)+440*(0.01/1.64)2)f = 0.0835 f=0.039*(9/44*0.037+2/11*0.038+3/22*0.039+15/88*0.04+13/88*0.

24、041+7/44*0.042)/(9/44*0.0372+2/11*0.0382+3/22*0.0392+15/88*0.042+13/88*0.0412+7/44*0.0422)+440*(0.01/1.64)2)f = 0.0857 f=0.04*(9/44*0.037+2/11*0.038+3/22*0.039+15/88*0.04+13/88*0.041+7/44*0.042)/(9/44*0.0372+2/11*0.0382+3/22*0.0392+15/88*0.042+13/88*0.0412+7/44*0.0422)+440*(0.01/1.64)2)f = 0.0879 f=

25、0.041*(9/44*0.037+2/11*0.038+3/22*0.039+15/88*0.04+13/88*0.041+7/44*0.042)/(9/44*0.0372+2/11*0.0382+3/22*0.0392+15/88*0.042+13/88*0.0412+7/44*0.0422)+440*(0.01/1.64)2)f = 0.0901 f=0.042*(9/44*0.037+2/11*0.038+3/22*0.039+15/88*0.04+13/88*0.041+7/44*0.042)/(9/44*0.0372+2/11*0.0382+3/22*0.0392+15/88*0.

26、042+13/88*0.0412+7/44*0.0422)+440*(0.01/1.64)2)f = 0.092390，80，60，75，65，707.317，6.68，5.142，6.5925，5.8565，6.461=38.327n=8+7+5+7+6+6=39 假设t=1.5,p=60,Pmax=6000，Tmax=176, =11附件3：c=-39; A=39;39;1; b=117.3;100;11; Aeq=; beq=; vlb=0;0;0; vub=;x,fval=linprog(c,A,b,Aeq,beq,vlb,vub)x = 2.5641fval = -100.0000

27、w=1 1 3 3 5 5 6 7 1 9 9;1 1 3 3 5 5 6 7 1 9 9;1/3 1/3 1 1 5/3 5/3 2 7/3 1/3 3 3;1/3 1/3 1 1 5/3 5/3 2 7/3 1/3 3 3;1/5 1/5 3/5 3/5 1 1 6/5 7/5 1/5 9/5 9/5;1/5 1/5 3/5 3/5 1 1 6/5 7/5 1/5 9/5 9/5;1/6 1/6 1/2 1/2 5/6 5/6 1 7/6 1/6 3/2 3/2;1/7 1/7 3/7 3/7 5/7 5/7 6/7 1 1/7 9/7 9/7;1 1 3 3 5 5 6 7 1 9 9;

28、1/9 1/9 1/3 1/3 5/9 5/9 2/3 7/9 1/9 1 1;1/9 1/9 1/3 1/3 5/9 5/9 2/3 7/9 1/9 1 1;w = 1.0000 1.0000 3.0000 3.0000 5.0000 5.0000 6.0000 7.0000 1.0000 9.0000 9.0000 1.0000 1.0000 3.0000 3.0000 5.0000 5.0000 6.0000 7.0000 1.0000 9.0000 9.0000 0.3333 0.3333 1.0000 1.0000 1.6667 1.6667 2.0000 2.3333 0.3333

29、 3.0000 3.0000 0.3333 0.3333 1.0000 1.0000 1.6667 1.6667 2.0000 2.3333 0.3333 3.0000 3.0000 0.2000 0.2000 0.6000 0.6000 1.0000 1.0000 1.2000 1.4000 0.2000 1.8000 1.8000 0.2000 0.2000 0.6000 0.6000 1.0000 1.0000 1.2000 1.4000 0.2000 1.8000 1.8000 0.1667 0.1667 0.5000 0.5000 0.8333 0.8333 1.0000 1.166

30、7 0.1667 1.5000 1.5000 0.1429 0.1429 0.4286 0.4286 0.7143 0.7143 0.8571 1.0000 0.1429 1.2857 1.2857 1.0000 1.0000 3.0000 3.0000 5.0000 5.0000 6.0000 7.0000 1.0000 9.0000 9.0000 0.1111 0.1111 0.3333 0.3333 0.5556 0.5556 0.6667 0.7778 0.1111 1.0000 1.0000 0.1111 0.1111 0.3333 0.3333 0.5556 0.5556 0.66

31、67 0.7778 0.1111 1.0000 1.0000 x,y=eig(w)x = Columns 1 through 7 -0.9883 -0.5443 0.7667 0.7667 -0.9887 0.9780 0.9841 0.0988 -0.5443 -0.2898 - 0.4162i -0.2898 + 0.4162i 0.1232 -0.0768 0.0071 0.0329 -0.1814 0.1120 + 0.0340i 0.1120 - 0.0340i 0.0352 -0.0552 -0.0970 0.0329 -0.1814 0.1120 + 0.0340i 0.1120

32、 - 0.0340i 0.0352 -0.0552 -0.0970 0.0198 -0.1089 -0.1255 + 0.0472i -0.1255 - 0.0472i -0.0104 -0.0576 -0.0638 0.0198 -0.1089 -0.1255 + 0.0472i -0.1255 - 0.0472i -0.0104 -0.0576 -0.0638 0.0165 -0.0907 0.0560 + 0.0170i 0.0560 - 0.0170i 0.0176 -0.0276 -0.0485 0.0141 -0.0778 0.0471 + 0.0559i 0.0471 - 0.0

33、559i 0.0203 -0.0342 0.0111 0.0988 -0.5443 -0.2684 - 0.0868i -0.2684 + 0.0868i -0.0443 0.1501 0.0331 0.0110 -0.0605 -0.0162 - 0.0370i -0.0162 + 0.0370i 0.0309 0.0145 0.0227 0.0110 -0.0605 -0.0162 - 0.0370i -0.0162 + 0.0370i 0.0309 0.0145 0.0227 Columns 8 through 11 0.9883 0.9883 0.9883 0.9883 -0.0988

34、 + 0.0000i -0.0988 - 0.0000i -0.0988 -0.0988 -0.0329 + 0.0000i -0.0329 - 0.0000i -0.0329 -0.0329 -0.0329 - 0.0000i -0.0329 + 0.0000i -0.0329 -0.0329 -0.0198 + 0.0000i -0.0198 - 0.0000i -0.0198 -0.0198 -0.0198 - 0.0000i -0.0198 + 0.0000i -0.0198 -0.0198 -0.0165 + 0.0000i -0.0165 - 0.0000i -0.0165 -0.

35、0165 -0.0141 + 0.0000i -0.0141 - 0.0000i -0.0141 -0.0141 -0.0988 - 0.0000i -0.0988 + 0.0000i -0.0988 -0.0988 -0.0110 - 0.0000i -0.0110 + 0.0000i -0.0110 -0.0110 -0.0110 - 0.0000i -0.0110 + 0.0000i -0.0110 -0.0110 y = Columns 1 through 7 0 0 0 0 0 0 0 0 11.0000 0 0 0 0 0 0 0 0.0000 + 0.0000i 0 0 0 0

36、0 0 0 0.0000 - 0.0000i 0 0 0 0 0 0 0 0.0000 0 0 0 0 0 0 0 -0.0000 0 0 0 0 0 0 0 -0.0000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 Columns 8 through 11 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0.0000 + 0.0000i 0 0 0 0 0.0000 - 0.0000i 0 0 0 0 0.0000 0 0 0 0 0.0000 w1=x(:,

37、2)/sum(x(:,2)w1 = 0.2175 0.2175 0.0725 0.0725 0.0435 0.0435 0.0362 0.0311 0.2175 0.0242 0.0242=11，CI=0,CR=0,通过了一致性检验附件4：f=0.0325*(9/88*0.0325+1/11*0.033+3/44*0.0335+15/176*0.034+13/176*0.034+7/88*0.0345+1/1*0.0345+7/8*0.035+13/176*0.035+9/110*0.0355+41/440*0.036+71/880*0.0365)/(9/88*0.03252+1/11*0.0

38、332+3/44*0.03352+15/176*0.0342+13/176*0.0342+7/88*0.03452+1/1*0.03452+7/8*0.0352+13/176*0.0352+9/110*0.03552+41/440*0.0362+71/880*0.03652)+880*(0.01/1.64)2)f = 0.0847 f=0.033*(9/88*0.0325+1/11*0.033+3/44*0.0335+15/176*0.034+13/176*0.034+7/88*0.0345+1/1*0.0345+7/8*0.035+13/176*0.035+9/110*0.0355+41/4

39、40*0.036+71/880*0.0365)/(9/88*0.03252+1/11*0.0332+3/44*0.03352+15/176*0.0342+13/176*0.0342+7/88*0.03452+1/1*0.03452+7/8*0.0352+13/176*0.0352+9/110*0.03552+41/440*0.0362+71/880*0.03652)+880*(0.01/1.64)2)f = 0.0860 f=0.0335*(9/88*0.0325+1/11*0.033+3/44*0.0335+15/176*0.034+13/176*0.034+7/88*0.0345+1/1*

40、0.0345+7/8*0.035+13/176*0.035+9/110*0.0355+41/440*0.036+71/880*0.0365)/(9/88*0.03252+1/11*0.0332+3/44*0.03352+15/176*0.0342+13/176*0.0342+7/88*0.03452+1/1*0.03452+7/8*0.0352+13/176*0.0352+9/110*0.03552+41/440*0.0362+71/880*0.03652)+880*(0.01/1.64)2)f = 0.0873 f=0.034*(9/88*0.0325+1/11*0.033+3/44*0.0

41、335+15/176*0.034+13/176*0.034+7/88*0.0345+1/1*0.0345+7/8*0.035+13/176*0.035+9/110*0.0355+41/440*0.036+71/880*0.0365)/(9/88*0.03252+1/11*0.0332+3/44*0.03352+15/176*0.0342+13/176*0.0342+7/88*0.03452+1/1*0.03452+7/8*0.0352+13/176*0.0352+9/110*0.03552+41/440*0.0362+71/880*0.03652)+880*(0.01/1.64)2)f = 0

42、.0886 f=0.0345*(9/88*0.0325+1/11*0.033+3/44*0.0335+15/176*0.034+13/176*0.034+7/88*0.0345+1/1*0.0345+7/8*0.035+13/176*0.035+9/110*0.0355+41/440*0.036+71/880*0.0365)/(9/88*0.03252+1/11*0.0332+3/44*0.03352+15/176*0.0342+13/176*0.0342+7/88*0.03452+1/1*0.03452+7/8*0.0352+13/176*0.0352+9/110*0.03552+41/44

43、0*0.0362+71/880*0.03652)+880*(0.01/1.64)2)f = 0.0899 f=0.035*(9/88*0.0325+1/11*0.033+3/44*0.0335+15/176*0.034+13/176*0.034+7/88*0.0345+1/1*0.0345+7/8*0.035+13/176*0.035+9/110*0.0355+41/440*0.036+71/880*0.0365)/(9/88*0.03252+1/11*0.0332+3/44*0.03352+15/176*0.0342+13/176*0.0342+7/88*0.03452+1/1*0.0345

44、2+7/8*0.0352+13/176*0.0352+9/110*0.03552+41/440*0.0362+71/880*0.03652)+880*(0.01/1.64)2)f = 0.0912 f=0.0355*(9/88*0.0325+1/11*0.033+3/44*0.0335+15/176*0.034+13/176*0.034+7/88*0.0345+1/1*0.0345+7/8*0.035+13/176*0.035+9/110*0.0355+41/440*0.036+71/880*0.0365)/(9/88*0.03252+1/11*0.0332+3/44*0.03352+15/1

45、76*0.0342+13/176*0.0342+7/88*0.03452+1/1*0.03452+7/8*0.0352+13/176*0.0352+9/110*0.03552+41/440*0.0362+71/880*0.03652)+880*(0.01/1.64)2)f = 0.0925 f=0.036*(9/88*0.0325+1/11*0.033+3/44*0.0335+15/176*0.034+13/176*0.034+7/88*0.0345+1/1*0.0345+7/8*0.035+13/176*0.035+9/110*0.0355+41/440*0.036+71/880*0.036

46、5)/(9/88*0.03252+1/11*0.0332+3/44*0.03352+15/176*0.0342+13/176*0.0342+7/88*0.03452+1/1*0.03452+7/8*0.0352+13/176*0.0352+9/110*0.03552+41/440*0.0362+71/880*0.03652)+880*(0.01/1.64)2)f = 0.0938 f=0.0365*(9/88*0.0325+1/11*0.033+3/44*0.0335+15/176*0.034+13/176*0.034+7/88*0.0345+1/1*0.0345+7/8*0.035+13/1

47、76*0.035+9/110*0.0355+41/440*0.036+71/880*0.0365)/(9/88*0.03252+1/11*0.0332+3/44*0.03352+15/176*0.0342+13/176*0.0342+7/88*0.03452+1/1*0.03452+7/8*0.0352+13/176*0.0352+9/110*0.03552+41/440*0.0362+71/880*0.03652)+880*(0.01/1.64)2)f = 0.095190，80，60，75，65，70，80，70，65，72，82，717.623,6.88,5.238,6.645,5.75

48、9,6.293,7.192,6.384,5.928,6.66,7.6916,6.7521n=8+7+5+6+6+6+7+6+6+7+8+7=79附件5：S.T. 假设t=2,p=50,Pmax=4000，Tmax=176, =10c=-79; A=79;79;1; b=88;80;10; Aeq=; beq=; vlb=0;0;0; vub=;x,fval=linprog(c,A,b,Aeq,beq,vlb,vub)x = 1.0127fval = -80.0000附件6：w = 1.0000 1.0000 3.0000 3.0000 5.0000 5.0000 6.0000 7.0000 9

49、.0000 9.0000 1.0000 1.0000 3.0000 3.0000 5.0000 5.0000 6.0000 7.0000 9.0000 9.0000 0.3333 0.3333 1.0000 1.0000 1.6667 1.6667 2.0000 2.3333 3.0000 3.0000 0.3333 0.3333 1.0000 1.0000 1.6667 1.6667 2.0000 2.3333 3.0000 3.0000 0.2000 0.2000 0.6000 0.6000 1.0000 1.0000 1.2000 1.4000 1.8000 1.8000 0.2000 0.2000 0.6000 0.6000 1.0000 1.0000 1.2000 1.4000 1.8000 1.8000 0.1667 0.1667 0.5000 0.5000 0.8333 0.8333 1.0000 1.1667 1.5000 1.5000 0.1429 0.1429 0.4286 0.4286 0.7143 0.7143 0.8571 1.0000 1.2857 1.2857 0.1111 0.