数学建模葡萄酒的分级正式版.docx

《数学建模葡萄酒的分级正式版.docx》由会员分享，可在线阅读，更多相关《数学建模葡萄酒的分级正式版.docx（17页珍藏版）》请在得力文库 - 分享文档赚钱的网站上搜索。

1、承 诺 书我们仔细阅读了中国大学生数学建模竞赛的竞赛规则.我们完全明白，在竞赛开始后参赛队员不能以任何方式（包括电话、电子邮件、网上咨询等）与队外的任何人（包括指导教师）研究、讨论与赛题有关的问题。我们知道，抄袭别人的成果是违反竞赛规则的, 如果引用别人的成果或其他公开的资料（包括网上查到的资料），必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。我们郑重承诺，严格遵守竞赛规则，以保证竞赛的公正、公平性。如有违反竞赛规则的行为，我们将受到严肃处理。我们参赛选择的题号是（从A/B/C/D中选择一项填写）： A 我们的参赛报名号为（如果赛区设置报名号的话）： 所属学校（请填写完整的

2、全名）： 西安理工大学 参赛队员 (打印并签名) ：1. 郑晓东 2. 罗璐 3. 宫维静 指导教师或指导教师组负责人 (打印并签名)： 日期： 2013 年 05月 10 日葡萄酒质量的综合评价分析摘要近年来，随着人们生活水平的提高，葡萄酒也随之受到人们的喜爱，加之食品科学技术的提高，人们对葡萄酒的品质也有了更高的要求，本文就针对葡萄酒品质的相关问题进行建模，求解和有关分析。对问题一，首先基于两组评酒员对同一批葡萄酒的评价分数数据，采用假设检验中的t检验法建立评估两组数据差异的模型，运用Spss软件求解，得到两组数据存在显著性差异的结论，其次，通过计算两组数据的方差，用以比较稳定性，得到第二

3、组更可信的结论。对问题二，首先对酿酒葡萄理化指标数据进行标准化处理，经过主成分分析法将葡萄分为四个等级，其次，按可信度高的一组（第二组）得分将葡萄酒分为五级，综合两种分级，将酿酒葡萄分为了级。对问题三，首先同问题二对酿酒葡萄与葡萄酒的理化指标进行主成分分析，用Matlab的曲线拟合得到葡萄酒的得分,分别与酿酒葡萄与葡萄酒的理化指标之间的函数关系，再进行反解即得到酿酒葡萄与葡萄酒的理化指标之间联系。对问题四，采用灰色关联度分析的方法进行求解，分别求出酿酒葡萄的理化指标与葡萄酒质量的关联度、葡萄酒理化指标与其质量的关联度，通过关联度值的大小，即可看出酿酒葡萄和葡萄酒的理化指标对葡萄酒质量的影响大小

4、，并以此为基准来论证酿酒葡萄和葡萄酒的理化指标能否用来评价葡萄酒的质量。关键词：t检验 主成分分析 曲线拟合 灰色关联度分析一、问题重述确定葡萄酒质量时一般是通过聘请一批有资质的评酒员进行品评。每个评酒员在对葡萄酒进行品尝后对其分类指标打分，然后求和得到其总分，从而确定葡萄酒的质量。酿酒葡萄的好坏与所酿葡萄酒的质量有直接的关系，葡萄酒和酿酒葡萄检测的理化指标会在一定程度上反映葡萄酒和葡萄的质量。附件1给出了某一年份一些葡萄酒的评价结果，附件2和附件3分别给出了该年份这些葡萄酒的和酿酒葡萄的成分数据。请尝试建立数学模型讨论下列问题：1. 分析附件1中两组评酒员的评价结果有无显著性差异，哪一组结果

5、更可信？2. 根据酿酒葡萄的理化指标和葡萄酒的质量对这些酿酒葡萄进行分级。3. 分析酿酒葡萄与葡萄酒的理化指标之间的联系。4分析酿酒葡萄和葡萄酒的理化指标对葡萄酒质量的影响，并论证能否用葡萄和葡萄酒的理化指标来评价葡萄酒的质量？二、问题分析 2.1 问题一本题给出了两组评酒员对同一批葡萄酒的评价分数，在本文，采用假设检验中的t检验建立评估两组数据差异性的模型，研究两组评论员的评价是否存在差异，判断能否接受它们存在显著性差异的假设。若接受，则继续第二步：可靠性分析，分别对两组数据求方差，方差小的说明波动小，既评酒员的评价较稳定，可靠性高。2.2 问题二首先，我们利用问题一得到的结果，对可靠性高的

6、一组数据进行处理，降低评论员之间的差异，提高葡萄酒样品最终得分的可靠度。按得分对葡萄酒进行分级。然后，用标准化处理后的酿酒葡萄的理化指标对葡萄进行主成分分析。最后，结合葡萄酒的分级对酿酒葡萄进行分级。2.3 问题三首先，用处理酿酒葡萄的理化指标的方法对葡萄酒的理化指标做同样的处理，得到葡萄酒理化指标的主成分。然后，分别根据主成分获得红葡萄和红葡萄酒的的得分。通过曲线拟合，分别建立红葡萄得分和专家的评分之间的关系；红葡萄酒得分和专家评分之间的关系。最后，根据两种理化指标和专家的评分之间的关系，建立两种理化指标之间的关系。2.4 问题四运用灰色关联度分析的方法，定量描述酿酒葡萄与葡萄酒的理化指标对

7、葡萄酒质量的影响，以此为基准来论证酿酒葡萄和葡萄酒的理化指标能否用来评价葡萄酒的质量。三、问题假设1同种葡萄酒在同一组评酒员的得分下成正态分布。2一种葡萄对应酿制一种葡萄酒。3葡萄的成分充分转化为葡萄酒里的成分，不存在意外的浪费和挥发。四、符号说明这里只列出主模型的全局参数，其他局部参数见具体模型。：第i个红葡萄酒样品：第i个红葡萄酒样品的得分T1：第一组评酒员全体T2：第二组评酒员全体五、模型的建立与求解5.1 模型一：基于t检验建立差异评估模型我们采用假设性检验验证是否能接受两组评酒员的评价结果存在显著性差异的假设。然后用方差分析两组评酒员评价数据的波动，认为较平稳的一组比较可靠。5.1.

8、1、数据预处理我们在分析数据是发现了几个显著性的异常数据：第一组红酒数据样品20色调评酒员4号 数据缺失第一组白酒数据样品3持久性评酒员7号 怀疑多了一个7第一组白酒数据样品8口感分析评酒员2号 数据明显异常因为随机样本在均值附近振荡，所以我们选用均值来代替异常数据以求误差最小。5.1.2 t检验模型的建立分别代表第一，第二组整体，分别对红葡萄酒（i=1,2，27）和白葡萄酒（i=1,2,27）进行感官评价，的评价结果通过组内的每一评酒员的评分的均值来表示。同样的，T2的评价结果通过组内的每一评酒员的评分的均值来表示。从而得到两组评论员分别对红葡萄酒的评价结果见表一：表1 红葡萄酒的评价结果表

9、中对于同一酒样品的两个评价差异是由两个评酒员引起的，为鉴定他们的评价结果有无显著性差异，可对两组数据对同一样品的差值进行分析，既表中的D。以红葡萄酒为例：有27对相互独立的评价结果（X1，Y1）（X2，Y2）（X27,Y27）,D1=X1-Y1,D2=X2-Y2,D27=X27-Y27，由于Di（i=1,2，,27）是由同一因素造成的，可认为它们服从同一分布。现假设DiN（，），i=1,2,27,且，未知，基于这一样本检验假设： （1）分别记的样本均值和样本方差的观测值为，。对进行单个均值的t检验，检验问题的拒绝域为（显著水平为）：. （2） 当的值不落在拒绝域内时，接受，既认为两组评价无显著

10、性差异。否则，两组评价有显著性差异。对白葡萄酒的处理同红葡萄酒。2）模型的求解 现以红葡萄酒为例求解，首先，作出同一酒样品分别由两组品酒员、得到的评价结果之差,列于表一的第三行，根据建立的模型检验假设： . （3） 取=0.05，运用spss软件求解得到表二： 表二 t检验求解结果根据上表得到的Sig=0.0200及相应的正交化单位特征向量：则X的第i个主成分为Fi = aiX i = 1，2，p。步骤三：选择主成分 在已确定的全部p个主成分中合理选择m个来实现最终的评价分析。一般用方差贡献率 （10） 解释主成分Fi所反映的信息量的大小，m的确定以累计贡献率 （11）达到足够大（一般在85%

11、以上）为原则。步骤四：计算主成分得分计算n个样品在m个主成分上的得分： i=1,2,3,，m （12）5.2.1.3 模型的求解 利用MATLAB软件编程，对酿酒葡萄的理化指标进行主成分分析（以红葡萄的指标为例），根据累计贡献度大于85%的原则筛选，得到的前12个特征值及其贡献度率如表所示： 表三 酿酒红葡萄理化指标的主成分分析结果序号特征值贡献度19.79705 0.166128.02727 0.136136.79664 0.115245.31105 0.090055.15878 0.087463.58111 0.060772.7598 0.046882.42551 0.041192.257

12、43 0.0383101.9615 0.0332111.58293 0.0268121.47825 0.0251得分表达式为： (13)得分如下所示：表四 酿酒葡萄的得分根据得分的大小对酿酒葡萄进行分级中上级酿酒葡萄：4,12,22,20 中级酿酒葡萄：19,7,27,3,25,26,23,16,5,13,17,24,14 中下级酿酒葡萄：6,21,2,9,10,15, 18 下级酿酒葡萄：11,8,15.2.2 根据评酒师的评分对葡萄酒的质量进行分级对于“根据评酒师的评分对葡萄酒的质量进行分级”这一问题，我们认为品质优良的葡萄酿出来的葡萄酒的品质也应是优良的。它们之间存在着一一对应关系。所以

13、可以通过专家评委们的打分对葡萄酒进行分级。根据我们对问题一模型的求解，得知第二组评委的评分更可信，所以针对于第二组评委对每种葡萄酒的综合得分的平均值进行排序，如下所示：品种92325262224得分78.277.175.874.674.57472.672.672.272.17271.671.52741611871171.571.269.968.868.868.368.268.166.36665.765.465.361.6表五 葡萄酒的评分葡萄酒样品的分级标准：8085分：高级葡萄酒7580分：中上级葡萄酒7075分：中级葡萄酒6570分：中下级葡萄酒6065分：下级葡萄酒所以采用上面的分级标准

14、，可将27中葡萄酒分为如下级别：中上级葡萄酒：9，23,20 中级葡萄酒：3,17，2,14,19,21,5,26，22,24,27,4中下级葡萄酒：16,10,13,12,25,1,6,8,15,18,7 下级葡萄酒：115.2.3 综合两种因素，对酿酒葡萄进行分级综合两种因素分级的结果，对酿酒葡萄进行综合评价，并得到酿酒葡萄的分级。具体步骤如下：1. 分别对两种分级结果进行编号：1) 中上级葡萄酒编号为02) 中级葡萄酒编号为13) 中下级葡萄酒编号为24) 下级葡萄酒编号为3对酿酒葡萄的分级结果编号同上。2.计算同一样品在两种不同情况下的编号与编号0的差值，分别记为，并计算总差值=+。3

15、.根据总差值的大小对酿酒葡萄进行分级。酿酒葡萄的分级标准：01：中上级酿酒葡萄 2：中级酿酒葡萄 35：中下级酿酒葡萄 6：下级酿酒葡萄 分级结果为：中上级酿酒葡萄：20,23,22,4 中级酿酒葡萄：12,9,3,19,17,27,24,26,5,14，中下级酿酒葡萄：7,25,16,13,6,21,2,10,15,1811,8 下级酿酒葡萄：15.3 模型三对于问题三，我们同样采用主成分分析法，得到葡萄酒的主成分记为Gi i=1,2，n，n为主成分的个数。应用Matlab软件编程，对葡萄酒的主成分进行求解（以红葡萄酒为例）根据累计贡献度大于85%的原则筛选，得到的前5个特征值及其贡献度率如

16、表所示：表六 红葡萄酒理化指标的主成分分析结果序号特征值贡献度15.99601 0.399723.18036 0.212031.76514 0.117741.42734 0.095250.881763 0.0588得分表达式：得分如下所示： 表七 葡萄酒的得分经分析直接求解葡萄酒理化指标与酿酒葡萄理化指标之间的联系难度较大，但我们发现它们都可与葡萄酒质量建立数学关系，将葡萄酒质量作为因变量，葡萄酒理化指标与酿酒葡萄理化指标的主成分分别作为自变量，采用曲线拟合，即可分别得到两理化指标与葡萄酒质量间的数学关系，再反解出两指标间的数学关系，这样不但简化了求解过程而且可减小误差。根据曲线拟合的方法对酿

17、酒葡萄理化指标的主成分和评委的得分之间建立函数关系为y=0.62598x8+1.7747x7-0.74728x6-4.457x5- 0.70899x4+2.8925x3+0.34017x2-0.36971x+3.5688图一 酿酒葡萄理化指标的主成分与评委得分的关系图根据曲线拟合的方法对葡萄酒理化指标的主成分和评委的得分之间建立函数关系为y=-4550.179x7-6212.6221x6-2244.2726x5+120.8559x4+160.1517x3+4.9978x2-3.2929x+3.5441图二 葡萄酒理化指标的主成分和评委的得分的关系图图三 酿酒葡萄理化指标的主成分和葡萄酒理化指标

18、的主成分的关系图5.4 模型四： 运用灰色关联度求解相互联系5.4.1 问题分析该问题需要分析酿酒葡萄和葡萄酒的理化指标对葡萄酒质量的影响，为方便运算，我们分别考虑酿酒葡萄的理化指标对葡萄酒质量的影响与葡萄酒的理化指标对酿酒葡萄的影响，为此我们引入灰色关联度的概念。灰色关联度能够定量描述事物或因素相互变化的情况，即变化的大小，方向和速度方面的关系。关联度越大，代表相互间联系越紧密，变化态势越一致，反之，若变化态势越不一致，则关联度越小。5.4.2 模型的建立与求解灰色关联度方法的计算介绍：1) 原始数据的标准化由于原始数据存在数量级和量纲的差异性，所以先进行标准化而方便进行计算。 =，= k=

19、1，2，3n2) 关联度的计算经数据处理后的参考数列为: 比较数列为; i=1,2,3m 从几何角度看，关联程度实质是参考数列与比较数列曲线形状的相似程度，参考数列与比较数列曲线形状接近，则两者关联度较大；反之参考数列与比较数列曲线形状相差较大，则两者间的关联度较小。因此，可用曲线间的差值大小作为关联度的衡量指标。则： k=1，2，3，n两极最大值和最小值： =，=关联系数：=，k=1，2，n，i=1，2，m式中为分辨系数，用来削弱过大而使关联系数失真的影响。人为引入这个系数是为了提高关联系数之间的差异显著性。3) 关联度的计算与比较由于每个比较数列与参数数列的关联程度是通过n个关联系数来反映

20、的，关联系数分散，不便于从整体上进行比较。因此，有必要对关联信息做集中处理，求平均值则是一种比较好的信息集中的方式。即用比较数列与参数数列各个时期的关联系数的平均值来定量反映两个数列的关联程度，其计算公式为：接下来用灰色关联度进行分析求解：（1） 对数据进行标准化处理 表八 标准化处理数据（2） 计算参考数列与比较数列在同一时期的绝对差，得到下表：花色苷 单宁 总酚 酒总黄酮 白藜芦醇 DPPH L*(D65) a*(D65) b*(D65) H(D65) C(D65)2.7144 0.5420 0.6175 0.6617 0.3043 0.6276 0.9155 0.6563 0.8002

21、0.6942 0.67770.8677 0.4310 0.4323 1.6222 0.0886 0.9653 0.7465 0.1850 0.0040 0.3861 0.16250.4529 0.7666 0.3064 0.4462 0.3867 0.7132 0.6593 0.1045 0.1899 0.4230 0.06100.2308 0.0347 0.0286 0.0469 0.1181 0.1351 0.1033 0.2556 0.2852 0.0751 0.24890.0957 0.1610 0.0029 0.2220 0.4104 0.0410 0.1273 0.2283 0.1

22、232 0.0055 0.20060.5481 0.0206 0.0564 0.0839 0.2291 0.0463 0.3946 0.0890 0.0759 0.1036 0.07650.7125 0.5043 0.4409 0.4920 0.5552 0.5559 0.3793 0.0875 0.3957 0.5244 0.00502.5535 0.7276 0.6904 0.6543 0.6479 0.6217 0.7184 0.1563 0.2629 0.1067 0.17970.3942 0.7047 0.7305 0.9473 0.0120 0.6501 0.7303 0.1602

23、 0.0203 0.3663 0.13800.4432 0.2026 0.2756 0.3265 0.0747 0.3596 0.4207 0.1838 0.6065 1.9912 0.08650.8289 0.2423 0.2317 0.4444 0.7687 0.4035 1.2872 0.6328 0.0111 0.8897 0.45960.6133 0.0431 0.1630 0.3221 0.3361 0.3014 0.3746 0.0696 0.4534 0.3400 0.13040.2560 0.1355 0.2272 0.2062 0.6454 0.2695 0.0019 0.

24、1517 0.1266 0.1697 0.10060.0680 0.1855 0.2212 0.3948 0.4247 0.2928 0.4318 0.0938 0.5774 0.4565 0.17320.5653 0.4814 0.3811 0.6550 0.7850 0.7246 0.2589 0.1190 0.1653 0.1690 0.06740.3799 0.3648 0.3844 0.4851 0.7097 0.5061 0.1986 0.1504 0.2056 0.2688 0.08970.0521 0.3216 0.0441 0.0627 0.4858 0.4471 0.062

25、6 0.1719 0.1972 0.0871 0.16520.6949 0.3554 0.2726 0.2465 0.4882 0.3489 0.4486 0.1189 0.0538 0.0284 0.09830.2573 0.1868 0.1867 0.2187 1.4766 0.2644 0.1511 0.2791 0.0738 0.2356 0.21700.7409 0.2157 0.2473 0.1969 0.7571 0.3177 0.8874 0.4988 0.3040 0.4245 0.46820.1394 0.3390 0.3776 0.1430 1.3727 0.5505 0

26、.5263 0.0014 0.5450 0.5282 0.09310.1580 0.1313 0.1198 0.0806 0.1211 0.0755 0.1223 0.0362 0.0773 0.3298 0.02770.5445 0.4744 0.9758 1.4557 2.4693 1.5097 0.6688 0.0938 0.1113 0.3315 0.06720.0093 0.2233 0.1532 0.2524 0.8779 0.2762 0.0285 0.1583 0.2135 0.3157 0.09590.4147 0.2716 0.3093 0.3985 0.3052 0.27

27、72 0.2072 0.2506 0.4934 1.2301 0.15140.3622 0.4387 0.3155 0.4963 0.1823 0.5962 0.1209 0.2956 0.3853 0.1287 0.2984 表九 参考数列与比较数列的绝对差（3） 通过上表可以计算出两极最大差与最小差，从而计算出关联系数，通过查阅文献得到得到关联系数表如下：葡萄酒样品花色苷 单宁 总酚 酒总黄酮 白藜芦醇 DPPH L*(D65) a*(D65) b*(D65) H(D65) C(D65)1 0.1671 0.5017 0.4691 0.4518 0.6425 0.4650 0.3732 0

28、.4539 0.4052 0.4400 0.4459葡萄酒样品 花色苷 单宁 总酚 酒总黄酮 白藜芦醇 DPPH L*(D65) a*(D65) b*(D65) H(D65) C(D65)1 0.1671 0.5017 0.4691 0.4518 0.6425 0.4650 0.3732 0.4539 0.4052 0.4400 0.44592 0.3859 0.5589 0.5581 0.2514 0.8619 0.3609 0.4221 0.7478 0.9952 0.5859 0.7716 2 0.3859 0.5589 0.5581 0.2514 0.8619 0.3609 0.422

29、1 0.7478 0.9952 0.5859 0.77163 0.5466 0.4156 0.6409 0.5503 0.5855 0.4333 0.4527 0.8407 0.7428 0.5635 0.90133 0.5466 0.4156 0.6409 0.5503 0.5855 0.4333 0.4527 0.8407 0.7428 0.5635 0.90134 0.7035 0.9423 0.9524 0.9229 0.8234 0.8028 0.8423 0.6816 0.6573 0.8807 0.68744 0.7035 0.9423 0.9524 0.9229 0.8234

30、0.8028 0.8423 0.6816 0.6573 0.8807 0.68745 0.8523 0.7733 0.9973 0.7116 0.5710 0.9322 0.8121 0.7058 0.8171 0.9925 0.73215 0.8523 0.7733 0.9973 0.7116 0.5710 0.9322 0.8121 0.7058 0.8171 0.9925 0.73216 0.4989 0.9659 0.9082 0.8684 0.7050 0.9238 0.5806 0.8614 0.8796 0.8419 0.8787 6 0.4989 0.9659 0.9082 0

31、.8684 0.7050 0.9238 0.5806 0.8614 0.8796 0.8419 0.87877 0.4336 0.5198 0.5533 0.5259 0.4957 0.4953 0.5902 0.8634 0.5799 0.5100 0.99347 0.4336 0.5198 0.5533 0.5259 0.4957 0.4953 0.5902 0.8634 0.5799 0.5100 0.99348 0.1758 0.4284 0.4413 0.4546 0.4571 0.4674 0.4315 0.7785 0.6755 0.8379 0.7532 8 0.1758 0.

32、4284 0.4413 0.4546 0.4571 0.4674 0.4315 0.7785 0.6755 0.8379 0.75329 0.5808 0.4363 0.4274 0.3652 0.9809 0.4562 0.4275 0.7741 0.9664 0.5986 0.79949 0.5808 0.4363 0.4274 0.3652 0.9809 0.4562 0.4275 0.7741 0.9664 0.5986 0.799410 0.5520 0.7301 0.6650 0.6261 0.8813 0.6031 0.5649 0.7490 0.4735 0.2148 0.86

33、48 10 0.5520 0.7301 0.6650 0.6261 0.8813 0.6031 0.5649 0.7490 0.4735 0.2148 0.864811 0.3968 0.6932 0.7027 0.5513 0.4150 0.5751 0.2974 0.4629 0.9825 0.3799 0.542911 0.3968 0.6932 0.7027 0.5513 0.4150 0.5751 0.2974 0.4629 0.9825 0.3799 0.542912 0.4708 0.9288 0.7711 0.6292 0.6192 0.6447 0.5932 0.8886 0

34、.5463 0.6165 0.808412 0.4708 0.9288 0.7711 0.6292 0.6192 0.6447 0.5932 0.8886 0.5463 0.6165 0.808413 0.6813 0.8023 0.7068 0.7266 0.4580 0.6700 0.9991 0.7836 0.8130 0.7638 0.845813 0.6813 0.8023 0.7068 0.7266 0.4580 0.6700 0.9991 0.7836 0.8130 0.7638 0.845814 0.8910 0.7472 0.7123 0.5805 0.5625 0.6513 0.5584 0.8549 0.4858 0.5446 0.760114 0.8910 0.7472 0.7123 0.5805 0.5625 0.6513 0.558