2012年数学建模A题范文 葡萄酒的评价.doc

2013高教社杯全国大学生数学建模竞赛承 诺 书我们仔细阅读了全国大学生数学建模竞赛章程和全国大学生数学建模竞赛参赛规则（以下简称为“竞赛章程和参赛规则”，可从全国大学生数学建模竞赛网站下载）。我们完全明白，在竞赛开始后参赛队员不能以任何方式（包括电话、电子邮件、网上咨询等）与队外的任何人（包括指导教师）研究、讨论与赛题有关的问题。我们知道，抄袭别人的成果是违反竞赛章程和参赛规则的，如果引用别人的成果或其他公开的资料（包括网上查到的资料），必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。我们郑重承诺，严格遵守竞赛章程和参赛规则，以保证竞赛的公正、公平性。如有违反竞赛章程和参赛规则的行为，我们将受到严肃处理。我们授权全国大学生数学建模竞赛组委会，可将我们的论文以任何形式进行公开展示（包括进行网上公示，在书籍、期刊和其他媒体进行正式或非正式发表等）。我们参赛选择的题号是（从A/B/C/D中选择一项填写）： A 我们的参赛报名号为（如果赛区设置报名号的话）： 所属学校（请填写完整的全名）： 参赛队员 (打印并签名) ：1. 2. 3. 指导教师或指导教师组负责人 (打印并签名)： （论文纸质版与电子版中的以上信息必须一致，只是电子版中无需签名。以上内容请仔细核对，提交后将不再允许做任何修改。如填写错误，论文可能被取消评奖资格。） 日期： 年 月 日赛区评阅编号（由赛区组委会评阅前进行编号）：2013高教社杯全国大学生数学建模竞赛编 号 专 用 页赛区评阅编号（由赛区组委会评阅前进行编号）：赛区评阅记录（可供赛区评阅时使用）：评阅人评分备注全国统一编号（由赛区组委会送交全国前编号）：全国评阅编号（由全国组委会评阅前进行编号）：27葡萄酒的评价 摘 要本文是关于葡萄酒评价方面的问题，具体分析阐述了对评酒员评价结果的分析检验和酿酒葡萄的分级，应用matlab，excel，spss等软件，结合F检验、聚类分析、典型相关分析等各种统计学方法，研究了葡萄、葡萄酒的理化指标及葡萄酒质量之间的联系。针对问题一，本文利用excel软件对各组评酒员的评分进行了整理和分析。在此基础上，对各个样品计算两组评价员评分的样本方差和方差D，运用F检验公式：检验是否有显著性差异，结果为：两组评酒员对红葡萄酒样品5、17、24的评价结果有显著性差异，其它样品无显著性差异；对白葡萄酒样品2、5、6、8、14、18、27的评价结果有显著性差异，对其它样品无显著性差异；且第二组评酒员的结果更加可信。针对问题二，考虑到葡萄酒的质量受酿酒葡萄的直接影响，我们首先根据葡萄酒的得分通过聚类分析得到初步排名，之后集合主成分分析法和相关性分析得出影响葡萄酒质量的主要理化指标加权求和得到新的排名，结合两个方面结果是，红酒中好酒有3、9、23差酒有10、25；白葡萄酒中好酒有5、28差酒有16、13、8。针对问题三，我们先使用spss软件对各个理化指标进行典型相关性分析，筛选出有较大相关性的成分，然后运用matlab软件进行一元线性回归，得到一元回归模型：即葡萄的理化指标X与葡萄酒的理花指标Y之间的线性关系。 针对问题四，可以与三题相似的方法，即典型相关性分析，找出理化指标中与葡萄酒的质量有较大相关性的物质，然后用这些物质与葡萄酒的质量进行多元线性回归，得到如下回归函数：并且通过验证，我们认为可以用葡萄的蛋白质、DPPH自归基、葡萄总黄酮、PH值、果皮颜色a、果皮颜色b和葡萄酒的总黄酮、白藜芦醇、DPPH半抑制体来评价葡萄酒的质量。关键词：F检验；K-均值分类；主成分分析；典型相关性分析；一元线性回归；多元线性回归；一 问题重述 这次论文要解决的问题是葡萄酒的评价问题，葡萄酒的评价是个很抽象和模糊的概念，没有固定的指标可以测量，每个人对同一种葡萄酒的评价都会不一样，而且，这对一个人的味觉要求特别高，必须要经过一定的训练才能掌握这样技巧，于是，评酒员这一职业就诞生了，实际生活中，确定葡萄酒质量时一般是通过聘请一批有资质的评酒员进行品评。每个评酒员在对葡萄酒进行品尝后对其分类指标打分，然后求和得到其总分，从而确定葡萄酒的质量。酿酒葡萄的好坏与所酿葡萄酒的质量有直接的关系，葡萄酒和酿酒葡萄检测的理化指标会在一定程度上反映葡萄酒和葡萄的质量。附件1给出了某一年份一些葡萄酒的评价结果，附件2和附件3分别给出了该年份这些葡萄酒的和酿酒葡萄的成分数据。有了以上数据以后，需要我们完成以问题：1. 分析附件1中两组评酒员的评价结果有无显著性差异，哪一组结果更可信？2. 根据酿酒葡萄的理化指标和葡萄酒的质量对这些酿酒葡萄进行分级。3. 分析酿酒葡萄与葡萄酒的理化指标之间的联系。4分析酿酒葡萄和葡萄酒的理化指标对葡萄酒质量的影响，并论证能否用葡萄和葡萄酒的理化指标来评价葡萄酒的质量？二 基本假设1、 所有的葡萄均正常，没有腐烂变质的现象。2、 所有葡萄酒在由相同工艺酿制而成。3、 所有葡萄酒都储存在相同材料的容器中。4、 所有的理化指标均测量无误。三 符号说明葡萄中各理化指标葡萄酒中各理化指标检验方差 样本方差四 问题分析葡萄酒的评价是一个复杂的过程，需要综合考虑不同评价员的评分，而且葡萄酒和葡萄的组成成分非常复杂，它们也要影响葡萄酒的质量，对如此繁多的数据，我们就必须依靠计算机工具，运用数学统计学知识对它们进行处理，并找出各个含量之间的关系，联系生活实际，对葡萄酒作出有理有据的评价。对于问题一：要想得到两组评价员的评价结果有无显著差异，并对它们的可靠性作出判断，我们首先就应该将两组评价员的对27组红葡萄酒和28组白葡萄酒的评价结果整理出来，再运用统计学中的F检验进行假设与检验，再通过计算各组评价员的评价结果的方差和区分度，进而判断谁的结果更加可信。 对于问题二：需要对葡萄进行分级，由于葡萄酒的质量与酿酒葡萄的好坏有直接关系，所以我们可以根据葡萄酒的质量对酿酒葡萄做一个简单的分级，之后，我们用主成分分析法算出每一组样本葡萄的哪些指标该葡萄的主成分，然后通过数据分析判断出这些成分哪些对葡萄酒的质量作出了贡献，筛选出主要成分后，对不同葡萄的成分做加权求和，以此作为葡萄分级的另一个依据。对于问题三：要想得到葡萄与葡萄酒的指标间的联系，即得到它们之间的函数关系表达式，我们必须对各种指标进行插值和拟合处理。但是，由于它们各自的指标太多，也不知道互相之间谁与谁有关系，有什么样的关系，所以不能草率地进行拟合。为此，首先必须要对数据进行筛选和比较，为此可以采用spss软件进行典型相关性分析，找出酿酒葡萄内部和葡萄酒内部各指标之间的相关性，和酿酒葡萄的指标与葡萄酒的指标间的相关性，从而选出具有代表性的指标进行拟合，得到一元线性回归方程。对于问题四：题中想要求出理化指标对质量的影响，即各理化指标与质量的线性或非线性关系，但是，由于理化指标太多，并且并非没个理化指标都会对葡萄酒的质量造成影响，所以首先必须进行数据的筛选，这里我们使用spss软件进行典型相关性分析，找出哪些指标与质量有较大的关系，然后将这些指标设为自变量，将质量设为因变量，对它们进行多元线性拟合，最后得到一个多元表达式以后，我们就可以通过这个方程来对葡萄酒的质量进行验证，如果验证的结果与评价员打分的结果基本吻合的话，就说明可以用葡萄与葡萄酒的理化指标来对葡萄酒的质量进行评价。五 模型的建立与求解5.1 模型的准备为了方便下面模型的建立与求解，我们需要对附件中的数据进行预处理，计算出下面需要用到的数据，例如，各个数据的平均值，方差，样本方差等，并制成表格，对于葡萄的分级，我们可以上网参考一下国家标准或者国外葡萄酒工业发达国家的标准进行分级。5.2 问题1模型的建立与求解根据题一的要求，要求得两组评价员的评价结果有无显著性差异，就必须求出两组评价员的评价结果。对于红葡萄酒样品一，我们先求出第一组评价员的所打分数的平均值，如下表所示： 评酒员组数12345678910第一组51664954776172617462第二组68718052537671737067表一：对红葡萄酒样品一的评价表 对于评价结果我们采用F检验法，采用公式：其中 ， 我们假设:第一组数与第二组数无明显差异，计算F的值，得到F=1.134618，取显著性水平，查F分布表得和。,由于0.24811.1346180.85 break; end end %记下累积贡献率大85%的特征值的序号放入newi中fprintf(主成分数：%gnn,length(newi);fprintf(主成分载荷：n)for p=1:length(newi) for q=1:length(y) result(q,p)=sqrt(newval(newi(p)*vec(q,newi(p); endend %计算载荷disp(result)1.3 cwscore.m%cwscore.m,计算得分function score=cwscore(vector1,vector2);sco=vector1*vector2;csum=sum(sco,2);newcsum,i=sort(-1*csum);newi,j=sort(i);fprintf(计算得分：n)score=sco,csum,j %得分矩阵：sco为各主成分得分；csum为综合得分；j为排序结果 1.4 cwprint.m%cwprint.mfunction print=cwprint(filename,a,b); %filename为文本文件文件名，a为矩阵行数(样本数)，b为矩阵列数(变量指标数)fid=fopen(filename,r)vector=fscanf(fid,%g,a b);fprintf(标准化结果如下：n)v1=cwstd(vector)result=cwfac(v1);cwscore(v1,result);二、部分运行结果贡献率：newrate = Columns 1 through 4 0.2281 0.1839 0.1348 0.0839 Columns 5 through 8 0.0637 0.0590 0.0445 0.0385 Columns 9 through 12 0.0294 0.0226 0.0216 0.0169 Columns 13 through 16 0.0152 0.0113 0.0107 0.0089 Columns 17 through 20 0.0071 0.0067 0.0042 0.0027 Columns 21 through 24 0.0021 0.0014 0.0011 0.0007 Columns 25 through 28 0.0005 0.0002 0.0000 -0.0000 Columns 29 through 32 0.0000 + 0.0000i 0.0000 - 0.0000i 0.0000 -0.0000 + 0.0000i Columns 33 through 34 -0.0000 - 0.0000i 0.0000 主成分数：9主成分载荷： Columns 1 through 7 -0.1509 -0.5142 0.3486 -0.4608 0.3600 0.0206 -0.2335 -0.5840 -0.1415 -0.5419 -0.2622 -0.0966 0.2063 -0.1004 0.0094 0.2886 -0.3378 -0.0019 0.2767 -0.4931 0.0374 -0.7480 -0.3931 -0.0936 0.3005 -0.2293 -0.0848 0.1108 -0.1412 -0.4424 -0.2045 -0.3621 -0.0355 0.4147 -0.0812 -0.1026 -0.5522 0.0592 0.6642 -0.2070 -0.1377 0.2258 -0.0312 -0.4832 -0.1482 0.3896 -0.1059 0.3710 -0.1266 -0.3395 -0.0948 0.2010 0.5897 0.0016 0.2968 -0.1745 -0.5163 -0.2875 -0.1375 0.7070 0.0646 0.0079 0.0211 -0.7272 -0.1921 -0.4267 -0.2159 -0.1833 -0.1628 -0.2165 -0.8010 -0.3353 -0.1077 -0.2307 -0.1900 -0.1720 0.0102 -0.7342 -0.2501 -0.0213 0.0557 -0.1876 -0.3658 -0.1721 -0.7034 -0.2022 -0.1793 -0.2901 -0.3314 -0.2304 -0.0650 0.3291 -0.4653 -0.6619 -0.0621 0.1520 -0.2104 -0.2065 -0.4525 -0.3175 -0.0541 0.0656 0.3118 0.0770 -0.6542 -0.0033 -0.5517 0.6230 -0.2682 0.0092 0.1028 -0.0039 0.1375 -0.4450 0.6150 -0.1323 0.0574 0.1683 0.0004 -0.0664 -0.4548 0.7216 -0.1562 -0.0577 0.0913 -0.0596 -0.2614 -0.0476 -0.3433 -0.6873 0.0353 0.2228 0.1850 0.2035 0.1368 0.7254 -0.0199 -0.0643 -0.4823 -0.2252 -0.1749 -0.3985 -0.3508 0.0261 -0.1615 0.5453 0.2641 -0.0903 -0.6464 0.6906 -0.1012 -0.0120 0.0830 -0.0636 0.0806 0.5295 -0.1609 -0.0532 -0.5797 0.3133 -0.2108 0.1848 0.6798 0.0895 -0.0840 -0.4928 -0.0412 -0.2181 -0.5371 -0.2126 -0.2912 0.2106 0.5200 -0.0588 -0.2434 -0.5757 -0.1140 -0.0287 -0.1809 -0.3153 -0.2782 0.2244 -0.0804 .5507 0.2462 0.1063 -0.4428 0.0863 -0.5440 0.4508 0.4130 -0.3731 0.0525 0.1954 0.2538 -0.2277 0.7334 -0.4691 -0.3724 -0.0248 -0.0413 -0.1890 -0.1389 0.7725 -0.4815 -0.2956 -0.0565 -0.0280 -0.1377 -0.0609 0.7868 -0.3703 -0.3411 -0.0319 -0.0638 -0.0486 -0.1435 -0.6685 0.6007 0.2193 -0.0455 0.2612 0.1152 0.1294 -0.4820 0.6528 0.1095 0.0325 0.2992 0.2130 -0.0667 -0.6540 0.6146 -0.0059 -0.0832 0.3410 -0