《数学建模--葡萄酒的分级(正式版)》由会员分享，可在线阅读，更多相关《数学建模--葡萄酒的分级(正式版)（23页珍藏版）》请在金锄头文库上搜索。
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 个红葡萄酒样品(1)iJ：第 i 个红葡萄酒样品的得分ijaT1：第一组评酒员全体T2：第二组评酒员全体五、模型的建立与求解5.1 模型一：基于 t 检验建立差异评估模型我们采用假设性检验验证是否能接受两组评酒员的评价结果存在显著性差异的假设。然后用方差
8、分析两组评酒员评价数据的波动，认为较平稳的一组比较可靠。5.1.1、数据预处理我们在分析数据是发现了几个显著性的异常数据：第一组红酒数据样品 20色调评酒员 4 号 数据缺失第一组白酒数据样品 3持久性评酒员 7 号 怀疑多了一个 7第一组白酒数据样品 8口感分析评酒员 2 号 数据明显异常因为随机样本在均值附近振荡，所以我们选用均值来代替异常数据以求误差最小。5.1.2 t 检验模型的建立分别代表第一，第二组整体，分别对红葡萄酒 （i=1,2，27）21,T iR和白葡萄酒 （i=1,2,27）进行感官评价， 的评价结果通过组内的每一iW1T评酒员的评分的均值来表示。同样的，T2 的评价结果
9、通过组内的每一评酒员的评分的均值来表示。从而得到两组评论员分别对红葡萄酒的评价结果见表一：红 葡 萄 酒 样 品 1 2 3 4 5 6 7 8 9 10 11 12 13 14第 一 组 评 分 62.7 80.3 80.4 68.6 73.3 72.2 71.5 72.3 81.5 74.2 70.1 53.9 74.6 73第 二 组 评 分 68.1 74 74.6 71.2 72.1 66.3 65.3 66 78.2 68.8 61.6 68.3 68.8 72.6D -5.4 6.3 5.8 -2.6 1.2 5.9 6.2 6.3 3.3 5.4 8.5 -14.4 5.8 0
10、.4红 葡 萄 酒 样 品 15 16 17 18 19 20 21 22 23 24 25 26 27第 一 组 评 分 58.7 74.9 79.3 59.9 78.6 79.15 77.1 77.2 85.6 78 69.2 73.8 73第 二 组 评 分 65.7 69.9 74.5 65.4 72.6 75.8 72.2 71.6 77.1 71.5 68.2 72 71.5D -7 5 4.8 -5.5 6 3.35 4.9 5.6 8.5 6.5 1 1.8 1.5表 1 红葡萄酒的评价结果表中对于同一酒样品的两个评价差异是由两个评酒员引起的，为鉴定他们的评价结果有无显著性差异
11、，可对两组数据对同一样品的差值进行分析，既表中的 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,且 ， 未知，基于这一样本检验假设：D2D2（1）0:,:10H分别记 的样本均值和样本方差的观测值为 ， 。对127,L d2Ds进行单个均值的 t 检验，检验问题的拒绝域为（显著水平为 ）：127, . （2） 2(1)Ddttns当 的值不落在拒绝域内时，接受 ，既认为
12、两组评价无显著性差异。否t 0H则，两组评价有显著性差异。对白葡萄酒的处理同红葡萄酒。2）模型的求解 现以红葡萄酒为例求解，首先，作出同一酒样品 分别由(1)iJ,27)L两组品酒员 、 得到的评价结果之差,列于表一的第三行，根据建立的模型检1T2验假设：. （3）01:,:0DDH 取 =0.05，运用 spss 软件求解得到表二： 表二 t 检验求解结果根据上表得到的 Sig=0.0200 及相应的正交化单位特征向量：则 X 的第 i 个主成分为 Fi = aiX i = 1，2， p。步骤三：选择主成分在已确定的全部 p 个主成分中合理选择 m 个来实现最终的评价分析。一般用方差贡献率（
13、10）解释主成分 Fi 所反映的信息量的大小， m 的确定以累计贡献率（11）达到足够大（一般在 85%以上）为原则。步骤四：计算主成分得分计算 n 个样品在 m 个主成分上的得分： i=1,2,3,，m （12）5.2.1.3 模型的求解利用 MATLAB 软件编程，对酿酒葡萄的理化指标进行主成分分析（以红葡萄的指标为例），根据累计贡献度大于 85%的原则筛选，得到的前 12 个特征值及其贡献度率如表所示：表三 酿酒红葡萄理化指标的主成分分析结果pppp aaa.,. 212121211/piik1()/pmikG12.iiipiFaXaX序号 特征值 贡献度1 9.79705 0.1661
14、2 8.02727 0.13613 6.79664 0.11524 5.31105 0.09005 5.15878 0.08746 3.58111 0.06077 2.7598 0.04688 2.42551 0.04119 2.25743 0.038310 1.9615 0.033211 1.58293 0.026812 1.47825 0.0251得分表达式为：(13)1=0.4320 10.38752+0.2533592=0.36191+0.00192+0.044759 12=0.01231+0.12292+0.172059 得分如下所示：品 种 1 2 3 4 5 6 7 8 9得 分 -2.1822 -0.3618 0.1567 1.2699 -0.0978 -0.3388 0.2743 -1.259 -0.4818品 种 10 11 12 13 14 15 16 17 18得 分 -0.5547 -1.2477 0.8408 -0.101 -0.2899 -0.7194 -0.0704 -0.1212 -1.0973品 种 19 20 21 22 23 24 25 26 27得 分 0.422 0.5899 -0.342 0.700
