数学建模葡萄酒的评价

葡萄酒的评价摘要葡萄拥有很高的营养价值，本文通过对葡萄酒的评价，以及酿酒葡萄和葡萄酒的理化指标之间的关系进行讨论分析，对不同的酿酒葡萄进行了分类，并更深入讨论两者的理化指标是否影响葡萄酒质量。针对问题一，我们首先分别计算每类葡萄酒样品在两组组评酒师评价下的综合得分，以此作为每组评酒师的最终评价结果。再运用统计学中的T检验进行假设与检验，得出两组评价结果具有显着性差异。最后通过计算各组评价员的评价结果的标准差，以此推算稳定性指标值P，P值较大的可信度较高，得出2pp红1红与2PP白1白，进而得出第二组的评价结果更加可信。针对问题二，我们分别对两组葡萄进行分类。在这里我们采用聚类分析法和主成分分析法，在matlab中实现对酿酒葡萄的分类。针对问题三，根据xZ对附件2中的数据进行标准化处理，排除单位不同的影响。以酿酒葡萄的30个一级理化指标作为自变量X，葡萄酒9个一级的理化指标作为因变量y,建立多元线性回归模型Xy,得出酿酒葡萄的理化指标与葡萄酒的理化指标之间的联系即回归系数矩阵。针对问题四，用灰色关联度分析对两者的关系进行度量，求得理化指标对样品酒的的关联系数。然后根据葡萄酒综合得分及指标的相关系数得出样品酒的综合指标，通过MATLAB软件对综合指标与第二问中葡萄酒的分数进行指数拟合，拟合效果不佳，因此不能定量的用葡萄和葡萄酒的理化指标来评价葡萄酒的质量，只能根据图像大致猜测综合指标与葡萄酒的质量负相关。关键词：T检验聚类分析法主成分分析法Z分数多元线性回归一、问题重述确定葡萄酒质量时一般是通过聘请一批有资质的评酒员进行品评。每个评酒员在对葡萄酒进行品尝后对其分类指标打分，然后求和得到其总分，从而确定葡萄酒的质量。酿酒葡萄的好坏与所酿葡萄酒的质量有直接的关系，葡萄酒和酿酒葡萄检测的理化指标会在一定程度上反映葡萄酒和葡萄的质量。附件1给出了某一年份一些葡萄酒的评价结果，附件2和附件3分别给出了该年份这些葡萄酒的和酿酒葡萄的成分数据。请尝试建立数学模型讨论下列问题：1.分析附件1中两组评酒员的评价结果有无显着性差异，哪一组结果更可信？2.根据酿酒葡萄的理化指标和葡萄酒的质量对这些酿酒葡萄进行分级。3.分析酿酒葡萄与葡萄酒的理化指标之间的联系。4.分析酿酒葡萄和葡萄酒的理化指标对葡萄酒质量的影响，并论证能否用葡萄和葡萄酒的理化指标来评价葡萄酒的质量？二、问题分析葡萄酒的评价是一个复杂的过程，需要综合考虑不同评价员的评分，而且葡萄酒和葡萄的组成成分非常复杂，它们也要影响葡萄酒的质量，对如此繁多的数据，我们就必须依靠计算机工具，运用数学统计学知识对它们进行处理，并找出各个含量之间的关系，联系生活实际，对葡萄酒作出有理有据的评价。对于问题一：要想得到两组评价员的评价结果有无显着差异，并对它们的可靠性作出判断，我们首先就应该将两组评价员的对27组红葡萄酒和28组白葡萄酒的评价结果整理出来，求得葡萄酒的综合得分，再运用统计学中的T检验进行假设与检验，判断两组是否存在显着性差异，再通过计算各组评价员的评价结果的标准差和稳定性指标，进而判断谁的结果更加可信。对于问题二：需要对葡萄进行分级，由于葡萄酒的质量与酿酒葡萄的好坏有直接关系，所以我们可以根据葡萄酒的质量对酿酒葡萄做一个简单的分级，之后，我们用主成分分析法算出每一组样本葡萄的哪些指标该葡萄的主成分，然后通过数据分析判断出这些成分哪些对葡萄酒的质量作出了贡献，筛选出主要成分后，对不同葡萄的成分做加权求和，以此作为葡萄分级的另一个依据。对于问题三：要想得到葡萄与葡萄酒的指标间的联系，即得到它们之间的函数关系表达式，必须求出两者指标之间的相关系数。但是，由于它们各自的指标太多，此处仅以一级指标作为相关因素进行分析。令酿酒葡萄的30个一级指标作为自变量，葡萄酒的9个一级指标作为因变量，建立线性回归模型，通过最小二乘法计算出回归系数，即酿酒葡萄的指标与葡萄酒的指标间的相关性。对于问题四：题中想要求出理化指标对质量的影响，即各理化指标与质量的线性或非线性关系，但是，由于理化指标太多，并且并非没个理化指标都会对葡萄酒的质量造成影响，所以首先必须进行数据的筛选，这里我们使用spss软件进行典型相关性分析，找出哪些指标与质量有较大的关系，然后将这些指标设为自变量，将质量设为因变量，对它们进行多元线性拟合，最后得到一个多元表达式以后，我们就可以通过这个方程来对葡萄酒的质量进行验证，如果验证的结果与评价员打分的结果基本吻合的话，就说明可以用葡萄与葡萄酒的理化指标来对葡萄酒的质量进行评价。三、基本假设1、假设评酒员对每种葡萄酒的评价结果是大致符合正态分布的；2、假设酿酒葡萄与葡萄酒中的芳香物质主要成分是：低醇、酯类、苯等，其余成份忽略；3、假设酿酒葡萄与葡萄酒的理化指标中一级指标为主要影响。4、假设酿酒葡萄中存在的而葡萄酒中不存在的理化指标也会影响葡萄酒的理化指标及质量；5、假设不考虑多种葡萄可制成一种酒，只考虑一种葡萄制成一种酒；6、假设只考虑红葡萄制成红葡萄酒，白葡萄制成白葡萄酒，忽略去皮红葡萄可酿制白葡萄酒；7、假设质量高的葡萄酒一定由质量好的酿酒葡萄制成，但是质量好的酿酒葡萄不一定能酿制成质量高的葡萄酒；8、ijA表示第i瓶酒的第j个指标无量纲化后的值9、ijB表示第i种酿酒葡萄的第j个指标无量纲化后的值10、iM表示第i瓶酒的综合指标四符号说明:T统计量T:khija第k组序号为h的样品第i个指标第j个品酒师的给分:khia序号为h的样品中第i个指标第k组10位品酒师给分的平均值:khiS第k组序号为h的样品第i个指标10位品酒师评分的标准差kib：第k组第i个指标所占权重:khx第k组序号为h的样品的稳定性指标k:p红第k组红葡萄酒的评分总平均稳定性指标k:P白第k组白葡萄酒的评分总平均稳定性指标ijX:为第i个样品的第j个指标is:第i个葡萄样品的总得分i:第i个样品葡萄理化指标得分为其中：第一个指标指澄清度，第二个指标指色调，第三个指标指香气纯正度，第四个指标指香气浓度，第五个指标指香气质量，第六个指标指口感纯正度，第七个指标指口感浓度，第八个指标指持久性，第九个指标指口感质量，第十个指标指平衡/整体评价。五模型建立与求解5.1问题一：葡萄酒评价结果的显着性差异及可信度分析5.1.1葡萄酒评价结果数据预处理对附件1中数据通过Excel筛选观察时可发现某些数据错误，如：第一组红葡萄酒品尝评分中酒样品20号下4号品酒员对于外观分析的色调评价数据缺失；第一组白葡萄酒品尝评分中酒样品3号下7号品酒员对于口感分析的持久性评价数据为77，明显超过该项上限8；第一组白葡萄酒品尝评分中酒样品8号下9号品酒员对于口感分析的持久性评价数据为16，明显超过该项上限8等。对这些异常数据为减少其对于总体评价结果的影响，采取预处理：取该酒样对应误差项目其余品酒员评价结果平均值替代该异常数据。经过数据预处理可得出每一种类葡萄酒的综合得分，建立表1与表2。表1红葡萄酒总得分平均值红酒n12345678910第一组第二组746611121314151617181920第一组73第二组21222324252627第一组7873第二组72根据表1，用excel作出两组评酒师对每一类葡萄酒的评分折线图。图1表2红葡萄酒总得分平均值白酒n12345678910第一组8271第二组11121314151617181920第一组7274第二组2122232425262728第一组71第二组77根据表2，用excel作出两组评酒师对每一类葡萄酒的评分折线图。图2根据图1、图2可初步简单看出两组评酒师的评价结果存在有显着性差异。5．1．2葡萄酒评价结果差异性分析与可信度分析模型建立与求解(1)t检验模型建立首先假定两个总体平均数间没有显着差异，即210:H查T值表，比较计算得到的T值与理论T值，推断发生概率（一般为95%）。两个正态总体的均值检验模型假设nXXX,...,,21是来自总体211,N的样本nYYY,...,,21是来自总体222,N的样本，且两样本独立。设1，2和2221,均未知，其检验问题为210:H.且2t~11)(2121321nnnnSYX.当0H为真时，统计量T的计算公式2~1121213nntnnSYXT.式中，211212222113nnSnSnS.查T值表，比较计算得到的T值与理论T值，推断发生概率（一般为95%），其中为显着性水平，05.010095-1因此当05.0T则认为0H不成立，两组评酒员对红葡萄酒的评价结果有显着性差异。（2）两组评酒员对红葡萄酒的评价结果比较：分别计算出7.3426S,73.0556,2711Xn05.00210.0T,说明该两组评酒员对红葡萄酒的评价结果有显着性差异。（3）两组评酒员对白葡萄酒的评价结果比较：分别计算出4.8266S,73.9786,2811Xn05.00129.0T,说明该两组评酒员对白葡萄酒的评价结果有显着性差异。5.1.3可信度分析模型建立与求解：第k组序号为h的样品第i个指标10位品酒师给分的平均值第k组序号为h的样品第i个指标10位品酒师的标准差算出第k组序号为h的样品的稳定性指标第k组红，白葡萄酒的评分总平均稳定性指标计算求得：比较红葡萄酒的两组总平均稳定性指标，因为2pp红1红,所以第二组品酒师的评价结果更可信。同样，比较白葡萄酒的总平均稳定性指标，因为2PP白1白，所以第二组品酒师的评价结果可信度更高。问题二：根据酿酒葡萄的理化指标和葡萄酒的质量对这些酿酒葡萄进行分级。问题二求根据酿酒葡萄的理化指标和葡萄酒的质量对酿酒葡萄进行分级，葡萄酒由酿酒葡萄酿制而成，则酿酒葡萄的质量与葡萄酒的质量有着直接的关系，则可以根据葡萄酒的质量对酿酒葡萄做一个简单的分级，在根据主成分分析从葡萄的理化指标中筛选出对葡萄质量产生影响的主要因素，根据所得各主要因素的贡献率给个因素加权作为系数，求出葡萄中主成分的含量，并进行排名，之后将此排名与之前根据葡萄酒质量所得出的排名综合，进而得出较准确的对酿酒葡萄的分级。K均值法聚类分析模型k均值法的基本步骤：(1)选择k个葡萄酒样品作为初始凝聚点，或者将所有葡萄酒样品分成k个初始类，然后将这k个类的重心(均值)作为初始凝聚点。(2)对除凝聚点之外的所有葡萄酒样品逐个归类，将每个葡萄酒样品归入凝聚点离它最近的那个类（通常采用欧氏距离），该类的凝聚点更新为这一类目前的均值，直至所有葡萄酒样品都归了类。(3)重复步骤(2)，直至所有的葡萄酒样品都不能再分配为止。最终的聚类结果在一定程度上依赖于初始凝聚点或初始分类的选择。经验表明，聚类过程中的绝大多数重要变化均发生在第一次再分配中。也就是：先算各类的均值再算各类中样本到本类及其他类的均值的绝对值距离（欧氏距离）将葡萄酒样本重新归类到欧氏距离较小的类中（重新归类就得算均值）首先，根据第一问得出的结果，我们采用第二组评酒员的结果作为判断葡萄酒质量的依据，根据各葡萄酒的分数，我们得出了红葡萄酒和白葡萄酒的排名，虽然是葡萄酒质量的排名，但由于葡萄酒的质量由酿酒葡萄的质量决定，所以上表可以看作是葡萄质量的排名，以上表中葡萄酒的分数作为酿酒葡萄质量的分数，可以对酿酒葡萄作出初步的分级，针对葡萄酒的成绩，我们用聚类分析的方法，得出了葡萄的初步分级，运行的得到的图样如下：图3图4根据上述结果，得出红、白葡萄酒的等级分类，建立表3，表4.表3红葡萄酒等级分类等级酒样品号A1，10,12,13,16，25B4,5,14,19,21,22,24,26,27C6,7,8,11,15,18D2,3,9,17,20,23表4白葡萄酒等级分类等级酒样品号A5，9,10,15,17,21,22,25,28B1,2,3,4,6,14,18,19,20,23,24,27C7,8,11,12,13,26D16主成分—权值分级模型虽然酿酒葡萄所对应葡萄酒的质量能在一定程度上反映