2012年电子科技大学中山学院优秀论文葡萄酒等级划分体系模型的探究 摘 要 针对目前葡萄酒评价体系不完善的现状，本文对葡萄酒评价体系作出探究。 对于问题一，运用单因素方差分析法，利用Matlab软件，以Anoval函数求解。求出p-value，显著性水平取0.05作为标准来判断那组有显著性，以及通过比较方差来判断那组数据更加可信。对于问题二，在问题一中得到第二组评分更可信，因此根据该组的评分进行分级，通过用Matlab软件的Corrcoef和Regress函数对该组成分进行相关性验证和用EXCEL画出图表进行分析，找出影响葡萄酒分级的成分，然后在酿酒葡萄数据中找出与影响葡萄酒分级相同的成分，再结合葡萄酒评分对葡萄样品进行分级，得出葡萄样品成分的排列，结合成分的量和葡萄酒分级得出影响酿酒葡萄分级成分的范围。 对于问题三，通过问题二的解答，可以知道葡萄酒和酿酒葡萄的划分级别，利用附件二的资料，对每一种理化指标的数据，根据对应的含量建立模型，运用matlab软件拟合数据，作出拟合线性图，并采用多元回归分析法进行回归分析，最后根据拟合线性图和回归系数来分析两类理化指标之间的关系。对于问题四，分析酿酒葡萄和葡萄酒的理化指标对葡萄酒质量的影响，并论证能否用葡萄和葡萄酒的理化指标来评价葡萄酒的质量。结合题目给出芬香物质的数据，对感官指标和理化指标进行综合分析，用MABTLE拟合感官指标和理化指标的数据，得出结论：需要结合葡萄酒的理化指标和感官指标对葡萄酒的质量进行综合评价。关键字：方差分析法 分级 理化指标 线性相关 回归分析一、问题的重述随着我国经济的快速发展，葡萄酒市场竞争也异常激烈和无序“三精一水”、假年份、假产地酒、假酒庄，影响消费者的健康，虽然我国的GB15037-2006葡萄酒国家标准对葡萄酒的质量作了规定，但由于相应规范的制定工作限制，我国关于葡萄酒质量等级分划的标准还未完善，国家迫切需要制定统一的质量等级制度。确定葡萄酒质量时一般是通过聘请一批有资质的评酒员进行品评。每个评酒员在对葡萄酒进行品尝后对其分类指标打分，然后求和得到其总分，从而确定葡萄酒的质量。酿酒葡萄的好坏与所酿葡萄酒的质量有直接的关系，葡萄酒和酿酒葡萄检测的理化指标会在一定程度上反映葡萄酒和葡萄的质量，文章给出了某一年份一些葡萄酒的评价结果及该年份这些葡萄酒的和酿酒葡萄的成分数据。本文尝试解决以下问题：问题一： 由于评酒师对葡萄酒的评分存在主观性，需对评酒师的分数进行客观分析，分析两组评酒员的评价结果有无显著性差异，哪一组结果更可信？问题二：葡萄酒的质量离不开原料酿酒葡萄的质量，所以酿酒葡萄的理化指标至关重要。需根据酿酒葡萄的理化指标和葡萄酒的质量对这些酿酒葡萄进行分级。问题三：酿酒葡萄与葡萄酒的理化指标之间的联系可能影响着葡萄酒质量，所以需建立模型，酿酒葡萄与葡萄酒的理化指标之间的联系。问题四：分析酿酒葡萄和葡萄酒的理化指标对葡萄酒质量的影响，并论证能否用葡萄和葡萄酒的理化指标来评价葡萄酒的质量，能否综合感官指标和理化指标，建立模型，来评价葡萄酒的质量是问题关键所在。 二、模型假设1. 品酒员打分相互之间没有影响；2.品酒员对样品的给的总分是他对该样品所有方面评分的总和，并且该样品的最终得分可认为是10位品酒员打分的平均值；3.题目所给的数据真实可靠；4.酿酒方式及酿酒过程对葡萄酒的质量没有影响；5.不同种类葡萄酒的成份数据值统一标准没有差异 ；6.所有样品的酿造过程相同。三、符号说明n测试数量r测试水平量A因素SS各类数据源的平方和Df各类数据相应的自由度MS各类的均方值F统计量P大于F的概率各组均值对总方差的偏差平方和各组数据对均值偏差平方和的总和 四、问题分析问题一的分析 我们要根据附件1的数据可知：评酒员对红酒27组样品，和白酒28组样品进行评分，每件样品都进行了两次评分，即是有两组评分数据，题目要求分析两组评酒员的评分结果有无显著性差异，以及那一组数据更加可信，对于显著性的判断，我们采用单因素方差分析法（Analysis Of Variance）。对于每件样品，评酒员对外观，香气，口感，及其整体评价进行打分，每一组的每件样品都有十名品酒员进行评分，故求每个品酒员对样品酒的总分，之后求出这十名品酒员给的总分的平均分，此平均分就是该样品的总分，葡萄酒分为白酒和红酒，我们对第一组的红酒和第二组的红酒进行方差分析法，运用matlab软件中的anova1函数可得出p-value，及F值，通过分析就可知道那组更加具有显著性。方差是考察数据的波动性的，方差小就说明数据比较稳定，方差大就是波动性比较大 ，故通过比较两组数据的方差大小，就知道那一组数据更加可信。问题二的分析根据问题一可知，第二组的评酒员的评酒分数更可靠，所以选择第二组葡萄酒的数据进行处理。从评酒员对葡萄酒评分的分数入手，用逆向思维反推葡萄的等级。首先将第一问中第二组的白葡萄酒和红葡萄酒的每一种样品的评分进行分等级，依次分为四个等级,然后用EXCEL将每个等级的样品酒的理化指标画成曲线图，忽略异常数据点，观察各等级间的理化指标有没有相关性，如果有相关性，找出影响葡萄酒质量的相关因素，跟酿酒葡萄的理化指标数据进行对照，得出酿酒葡萄的分级依据。 问题三的分析 结合葡萄酒和酿酒葡萄的理化指标，作出每两个理化指标间的直观趋势图，观察两者之间的大体关系，根据曲线拟合的方法得出两者间的函数关系。问题四的分析由第三问求解可得出酿酒葡萄与葡萄酒的理化指标之间是呈线性相关的，因此我们要证明酿酒葡萄和葡萄酒的理化指标对葡萄酒质量是有影响的，只需证明酿酒葡萄的理化指标对葡萄酒质量是有影响。在综合附录3给出的芬香物质，用MABTLE拟合出理化指标和感官指标的关系图呈相关性，所以要综合葡萄酒的理化指标和感官指标一起来评价葡萄酒的质量。 五、模型建立与求解5.1问题一的模型建立和求解 对于两组评酒员的评价结果有无显著性差异，我们采用单因素方差分析法去解决。单因素方差分析法: 只考虑一个因素A 对所关心的指标的影响，A 取几个水平，在每个水平上作若干个试验，试验过程中除A 外其它影响指标的因素都保持不变（只有随机因素存在),我们的任务是从试验结果推断，因素A 对指标有无显著影响，即当A 取不同水平时指标有无显著差别。A 取某个水平下的指标视为随机变量,判断A 取不同水平时指标有无显著差别，相当于检验若干总体的均值是否相等。设 A取n 个水平,在水平下总体 服从正态分步N(,),i=1,.,n,这里,未知，可以互不相同，但假定有相同的方差，又设在每个水平下作了次独立试验，即从中抽取容量为的样本，记作服从N(,)，i=1,n,j=1, 且且相互独立。将这些数据列成表1（单因素试验数据表）的形式。 表5.1 单因素试验数据表分值第一组红酒第二组红酒第一组白酒第二组白酒A1X12X21X12X21A2X21X22X21X22.A3X31X32X31X32根据上述理论，首先我们对数据进行处理，附件1里有四组数据：红葡萄酒和白葡萄酒各有两组数据，每种酒都有两组人进行对其进行评分，每件样品酒有十名品酒员号打分，采用单因素方差分析法，我们将样品酒的总分作为唯一考虑的因素A，运用matlab软件编程求出品酒员对每组样品打的总分的平均分，见下表： 表5.2组样品红酒和白酒的总分样品号第一组红葡萄酒品尝综合得分评分第一组白葡萄酒品尝综合得分评分第二组红葡萄酒品尝综合得分评分第二组白葡萄酒品尝综合得分评分162.78268.177.9280.374.27475.8380.479.774.675.6468.679.471.276.9573.37172.181.5672.268.466.375.5771.577.565.374.2872.371.46672.3981.572.978.280.41074.274.368.879.81170.172.361.671.41253.963.368.372.41374.665.968.873.914737272.677.11558.772.465.778.41674.97469.967.31779.378.874.580.31859.973.165.476.71978.672.272.676.42079.577.875.876.62177.176.472.279.22277.27171.679.42385.675.977.177.4247873.371.576.12569.277.168.279.52673.881.37274.3277364.871.5772881.379.6 对这四组数据，我们将白酒和红酒分开来判断其有无显著性，即第一组红酒与第二组红酒，第一组白酒和第二组白酒比较。运用matlab软件对数据处理编程得出以下结果，标准ANOVA表分析见下表：表5.3白葡萄酒ANOVA表图5.1 白葡萄酒盒型（box）图表5.4 红葡萄酒ANOVA表图5.2 红葡萄酒的盒型（box）图表5.5 方差分析表：方差来源平方和(SS)自由度 (df)均方（MS）1-P分数位F概率p因素A r-1误差n-r总和 通常情况下，实验结果p达到0.05水平或0.01水平，才可以说数据之间具备了差异显著或是极显著。在作结论时，应确实描述方向性（例如显著大于或显著小于）。sig值通常用 P0.05 表示差异性不显著。在此我们去0.05作为显著性水平标准，红酒中的ANOVA表中ProbF栏p值为0.0278 0.05，接受Ho，故没有显著性。 对于那组数据更加可信，我们知道方差是考察数据的波动性的，方差小就说明数据比较稳定，方差大就是波动性比较大。故我们将红酒，白酒每组样品酒一一对应，第一组的红酒中样品一与第二组红酒中的样品一进行方差分析，以此类推，我们将所求到的方差用matlab进行画图。图5.3 白葡萄酒的方差图图5.4 红葡萄酒的方差图从两副图中