现代葡萄酒业中统计方法的应用大学生数学建模竞赛论文.docx

1、现代葡萄酒业中统计方法的应用大学生数学建模竞赛论文现代葡萄酒业中统计方法的应用摘要：“现代葡萄酒业中统计方法的应用”数学模型运用了一些经典的统计原理和统计方法，探讨在已知葡萄及葡萄酒的理化指标数据的情况下，如何有效地利用来对葡萄酒进行一个简明初步的分级，并对葡萄与葡萄酒的理化指标之间的联系和两组不同品酒师对相同葡萄酒评价的差异性和可信度进行分析，再通过以葡萄及葡萄酒的理化指标对葡萄酒的分级与品酒师对葡萄酒的评定比较，最终得出对葡萄酒鉴定的合理方案。该数学模型借助SPSS软件运用T检验对两组品酒师对酒的评价的差异性进行了分析，又通过主成分分析对葡萄和葡萄酒的理化指标进行了归类和简化，进而对葡萄和

2、葡萄酒的理化指标对葡萄酒品质的贡献做了综合性的分析；肯德尔和谐系数、方差分析的运用和数据相关性系数的计算又为评价品酒师对酒品打分的客观性和数据的关联性提供了有力的依据。通过模型我们可以得出结论：两组品酒员评价结果无显著性差异，第二组品酒师的结果更可信；仅靠葡萄和葡萄酒的理化指标不能准确评价葡萄酒的质量。关键词：理化指标 统计方法 T检验 主成分分析 肯德尔和谐系数 相关系数 一、 问题重述与分析1.1问题的重述葡萄酒是酒中的精魂，飘渺于生命的内在，又令心灵交融；她是有生命的，来自于天地恩赐，日月光泽；她恍如披着柔幔的处女，娇憨清新，浓淡相宜。葡萄酒的质量是葡萄酒的一种特性，它是表示葡萄酒优秀的

3、程度。复杂性和协调性构成了葡萄酒质量的主要属性。复杂性意味着葡萄酒有浓郁的令人愉快的香气与滋味，及它的潜力与精巧性。协调性是指各种感官成分的结合是平衡的，每一种成分与其它成分相比，它的存在及其含量是非常适宜的，而且它的复杂性、颜色强度、香气、滋味、后味是稳定一致的。确定葡萄酒质量时一般是通过聘请一批有资质的评酒员进行品评。每个评酒员在对葡萄酒进行品尝后对其分类指标打分，然后求和得到其总分，从而确定葡萄酒的质量。酿酒葡萄的好坏与所酿葡萄酒的质量有直接的关系，葡萄酒和酿酒葡萄检测的理化指标会在一定程度上反映葡萄酒和葡萄的质量。在拥有附件1（某一年份一些葡萄酒的评价结果）、附件2（该年份这些葡萄酒的

4、成分数据）和附件3（酿酒葡萄的成分数据）的情况下，我们需要尝试建立数学模型讨论下列问题：1两组评酒员的评价结果有无显著性差异，哪一组结果更可信？2根据酿酒葡萄的理化指标和葡萄酒的质量对这些酿酒葡萄进行分级。3酿酒葡萄与葡萄酒的理化指标之间有什么样的理化的联系。4酿酒葡萄和葡萄酒的理化指标对葡萄酒质量的影响以及能否用葡萄和葡萄酒的理化指标来评价葡萄酒的质量？ 1.2问题的分析 由题意可知，目的就是为了建立数学模型，探讨在已知葡萄及葡萄酒的理化指标数据的情况下，如何有效地利用来对葡萄酒进行一个简明初步的分级，并对葡萄与葡萄酒的理化指标之间的联系和两组不同品酒师对相同葡萄酒评价的差异性和可信度进行分

5、析，再通过以葡萄及葡萄酒的理化指标对葡萄酒的分级与品酒师对葡萄酒的评定比较，最终得出对葡萄酒鉴定的合理方案。问题一中，首先可以用SPSS软件对两组品酒师评价进行T检验以确定两组品酒师评价的差异性，然后再利用SPSS求出两组品酒师评价结果的肯德尔和谐系数和方差以评定两组品酒师评价的可信度。而问题二，可以用主成分分析的方法对葡萄的理化指标和葡萄酒的质量做综合性的分析，以完成对酿酒葡萄的分级。由题意可知葡萄酒的质量由品酒师的打分确定，故需要对两组品酒师的打分结果进行综合评定。问题三，采用求相关系数的方法，求解出两种葡萄与对应的两种葡萄酒的理化指标之间的相关系数以确定酿酒葡萄与葡萄酒的理化指标之间有什

6、么样的理化的联系。问题四中，首先对葡萄和葡萄酒的理化指标进行综合的分级，然后求综合分级与品酒员对酒品评价高低的相关性，可知若相关性较高，则说明用葡萄和葡萄酒的理化指标来评价葡萄酒的质量具有一定的可行性，反之则说明不能用葡萄和葡萄酒的理化指标来评价葡萄酒的质量。二、 建模的过程1） 问题一1、模型假设：1、两组品酒员无信息交流，彼此之间打分无影响。2、样本足够大。 3、假设葡萄酒的酿造时期和被品尝时期相同。2、定义符号说明：R相关性系数3、模型的建立与求解差异性：如下图，对于红白葡萄酒质量评分，两组20名品酒员给出的结果如下：1.红葡萄样品123456789A62.780.380.468.673

7、.372.271.572.381.5B68.17474.671.272.166.365.36678.2样品101112131415161718A74.270.153.974.67358.774.979.359.9B68.861.668.368.872.665.769.974.565.4样品192021222324252627A78.678.677.177.285.67869.273.873B72.675.872.271.677.171.568.27271.52.白葡萄样品12345678910A8274.285.379.47168.477.571.472.974.3B77.975.875.67

8、6.981.575.574.272.380.479.8样品11121314151617181920A72.363.365.97272.47478.873.172.277.8B71.472.473.977.178.467.380.376.776.476.6样品2122232425262728A76.47175.973.377.181.364.881.3B79.279.477.476.179.574.37779.6对上述两表进行T检验，结果如下：1.白葡萄：组统计量AN均值标准差均值的标准误VAR000011.001074.26008.458742.674892.001076.53205.1661

9、21.63367独立样本检验方差方程的 Levene 检验均值方程的 t 检验FSig.tdfSig.(双侧)均值差值标准误差值差分的 95% 置信区间下限上限VAR00001假设方差相等2.808.111-.72518.478-2.272003.13431-8.856944.31294假设方差不相等-.72514.894.480-2.272003.13431-8.956764.41276对于白葡萄酒来说，F值为2.808，相伴概率为0.111，大于显著性水平0.05，不能拒绝方差相等的假设，可以认为第一组和第二组品酒员的评价结果无明显差异。然后看方差相等时的T检验结果，T统计量的相伴概率为0

10、.478，大于显著性水平0.05，不能拒绝T检验的零假设。也就是说，品酒员对白葡萄酒的评价结果不存在明显差异。另外从样本的均值差的95%置信区间看，区间跨0，这也说明两组品酒员对白葡萄酒的评价无显著性差异。2.红葡萄：组统计量AN均值标准差均值的标准误VAR000011.001073.66805.098071.612152.001073.57803.06119.96803独立样本检验方差方程的 Levene 检验均值方程的 t 检验FSig.tdfSig.(双侧)均值差值标准误差值差分的 95% 置信区间下限上限VAR00001假设方差相等2.426.137.04818.962.090001.

11、88046-3.860704.04070假设方差不相等.04814.743.962.090001.88046-3.924194.10419对于红葡萄酒来说，F值为2.426，相伴概率为0.137，大于显著性水平0.05，不能拒绝方差相等的假设，可以认为第一组和第二组品酒员的评价结果无明显差异。然后看方差相等时的T检验结果，T统计量的相伴概率为0.962，大于显著性水平0.05，不能拒绝T检验的零假设。也就是说，品酒员对红葡萄酒的评价结果不存在明显差异。另外从样本的均值差的95%置信区间看，区间跨0，这也说明两组品酒员对红葡萄酒的评价无显著性差异。相关性检验取两组的平均分：一组62.780.38

12、0.468.673.372.271.572.381.574.270.1二组68.17474.671.272.166.365.36678.268.861.6一组53.974.67358.774.979.359.978.678.677.177.2二组68.368.872.665.769.974.565.472.675.872.271.6一组85.67869.273.8738274.285.379.47168.4二组77.171.568.27271.577.975.875.676.981.575.5一组77.571.472.974.372.363.365.97272.47478.8二组74.272.

13、380.479.871.472.473.977.178.467.380.3一组73.172.277.876.47175.973.377.181.364.881.3二组76.776.476.679.279.477.476.179.574.37779.6相关性VAR00001VAR00002VAR00001Pearson 相关性1.454*显著性（双侧）.001N5555VAR00002Pearson 相关性.454*1显著性（双侧）.001N5555*. 在 .01 水平（双侧）上显著相关。R=0.4540.4 有一定的相关性总评：统计量AN均值标准差均值的标准误V211073.675.0981

14、.61221073.583.061.968独立样本检验方差方程的 Levene 检验均值方程的 t 检验FSig.tdfSig.(双侧)均值差值标准误差值差分的 95% 置信区间下限上限V2假设方差相等2.426.137.04818.962.0901.880-3.8614.041假设方差不相等.04814.743.962.0901.880-3.9244.104对于葡萄酒来说，F值为，相伴概率为2.426，大于显著性水平0.05，不能拒绝方差相等的假设，可以认为第一组和第二组品酒员的评价结果无明显差异。然后看方差相等时的T检验结果，T统计量的相伴概率为0.962，远大于显著性水平0.05，不能拒

15、绝T检验的零假设。也就是说，品酒员对葡萄酒的评价结果不存在明显差异。另外从样本的均值差的95%置信区间看，区间跨0，这也说明两组品酒员对葡萄酒的评价无显著性差异。可信度:1.对于同一种酒的同一指标的评价，在两组无显著性差异的情况下，一组10人的评价分数相差越小则可信度越高，评价分数的散聚程度由标准差的大小来评价。共27（28）组，标准差大小如下:x=1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28;y1=28.92 18.92 20.31 31.18 23.62 23.19 30.54 19.9

16、0 17.22 16.54 25.24 26.77 20.11 18.00 27.75 12.76 28.14 20.61 20.65 15.31 32.32 21.34 17.10 25.96 24.12 16.78 21.17 0;y2=28.81 42.54 57.32 20.06 33.73 38.27 18.77 40.65 28.89 43.75 39.93 32.28 39.20 32.06 34.42 40.02 36.02 37.54 20.44 24.07 39.43 35.33 19.82 31.62 17.46 25.61 36.05 26.91;y3=27.15 12

17、.08 16.63 19.28 11.09 13.79 23.75 24.21 15.22 18.04 18.50 15.04 11.73 14.44 19.29 13.45 9.08 21.27 22.28 18.75 17.88 14.78 14.93 9.82 19.84 19.34 13.58 0;y4=15.26 21.01 35.81 19.47 15.38 14.30 19.48 16.74 30.93 25.17 28.11 35.50 20.52 11.95 22.05 27.20 18.60 16.50 15.31 21.22 24.07 21.96 10.22 18.63

18、 30.96 30.43 17.89 15.11;比较直观的从下图来分析：上图中四条曲线分别表示A,B两组对红白葡萄酒的总体评价的标准差的大小，由蓝色和红色的两条曲线可以看出A组对红葡萄酒的评价差较B组来说更大，同样有棕色和浅蓝色的两条曲线得出白葡萄酒仍然是A族的标准差大。2.做两组品酒员所得结果的Alpha系数A组的Alpha系数：案例处理汇总N%案例有效55100.0已排除a0.0总计55100.0可靠性统计量Cronbachs Alpha项数0.82510B组的Alpha系数：案例处理汇总N%案例有效55100.0已排除a0.0总计55100.0可靠性统计量Cronbachs Alpha

19、项数0.83410由上图可以看出两组的Alpha系数均大于0.8，及说明数据的可靠性都较强，但第二组的Alpha系数略好于第一组，所以第二组的可靠性略好于第一组。所以综合以上两个方面，第二组的可靠性更高。2)问题二1、模型假设：1、假设品酒师对葡萄酒的打分作为酿酒葡萄质量的指标之一。2、假设酿酒葡萄理化指标中的数据不受单位的影响。4、假设样本的量足够大且具有代表性。2、定义符号说明：Xi-第i个理化指标ZXi-原始变量经过标准化处理之后的数据Hi-红葡萄样品iBi-白葡萄样品i3、模型建立与求解：对于酿酒葡萄的分级问题，先通过因子分析对葡萄的众多指标进行提炼，用少数几个变量来替代酿酒葡萄的变量

20、。表2.1将表2.1中的数据运用SPSS软件处理，对酿酒葡萄的26项指标进行因子分析。KMO值大于0.5，根据Kaiser给出的KMO度量标准可知原有的26个变量适合做因子分析。现得到“解释的总方差”结果：表2.2分析结果可得表2.3由表2.3可知因子变量1替代了总酚、葡萄酒黄酮、花色苷、单宁、蛋白质、DPPH自由基、黄酮醇，因子变量2替代了可溶性固形物、总糖、氨基酸总量，因子变量3替代了可滴定酸、固酸比、PH值三种指标，因子变量4替代了果皮质量、果穗质量、百粒质量三种指标，因子变量5替代了总分、酒石酸、果梗比三种指标，因子变量6替代了多酚氧化酶活力、褐变度两种指标，因子变量7替代了柠檬酸，因

21、子变量8出汁率、白藜芦醇两种指标。所以决定用这八个新变量来替代原有的26个变量。最后通过主成分分析模型来解决问题运用SPSS将红葡萄的各项指标进行标准化,得到如图2.4表2.4记表2.4中的数据为一矩阵HZ（27乘26）运用SPSS软件得到特征向量Ai（26乘8）将得到的特征向量与标准化后的数据相乘然后可以得出主成分表达式Fi=HZ*Ai其中Fi为一27乘8的矩阵，第i列表示27种酿酒葡萄（红）在第i个因子下的排名，其结果如表2.5所示表2.5之后以每个主成分所对应的特征值占所提取主成分总的特征值之和的比例作为权重计算主成分综合模型。结果如图2.5所示。同理，对白葡萄样品数据做同样的处理，所得

22、结果如图2.6所示。图2.62） 问题三1、 模型假设：1、 测理化性质所用的葡萄的其他条件相同。2、 所给数据真实可靠。2、定义符号说明： R-相关性系数3、模型的建立与求解：下列表格列举出了在葡萄酒与葡萄中都有的物质，白藜芦醇单宁花色苷葡萄总黄酮与酒总黄酮总酚取其中有相关性的作图1红葡萄酒单宁2红葡萄酒花色苷3红葡萄酒总酚4红葡萄酒黄酮5白葡萄酒总黄酮从以上图中，我们可以看出，红葡萄酒单宁，红葡萄酒花色苷，红葡萄酒总酚，红葡萄酒总黄酮，白葡萄酒总黄酮，均有很大的相关性。相同物质具有相关性必然具有一定的因果性质，所以以上指标在葡萄酒与酿酒葡萄之间有较强的联系，即酿酒葡萄中的量影响着葡萄酒中的

23、量。同一类的物质可以由作图得出但是不同的物质不能一一作图，所以由程序（代码见附录）得出相关性系数表（表1）。从表中提炼信息得出红葡萄花色苷单宁总酚酒总黄酮DPPH半抑制体积色泽l花色苷0.920.720.770.710.67-0.83褐变度0.770.450.460.440.38-0.56DPPH自由基0.570.750.810.760.78-0.71总酚0.610.820.880.880.87-0.75单宁0.660.720.740.70.7-0.68葡萄总黄酮0.440.680.820.820.81-0.61白葡萄酒总黄酮色泽l色泽b葡萄总黄酮0.70.1-0.2干物质含量0.11-0.7

24、20.68出汁率-0.050.69-0.76即上述涂有黄格的数据皆有相关性（R=0.7）。4）问题四1、模型假设： 1.葡萄酒的质量只由品酒师的打分确定。 2.葡萄的理化指标和葡萄酒的理化指标对葡萄酒的质量的贡献是一样大的。2、模型的建立与求解 利用SPSS求出葡萄和葡萄酒的理化指标与葡萄酒质量的相关系数，如下结果：红葡萄10.263290.065570.327850.2051130.0951410.12120.27830-0.195580.19286-0.072114-0.072142-0.00330.1931-0.16459-0.18287-0.2451150.0771430.064-0.

25、28320.19760-0.12388-0.0491160.017144-0.07650.39233-0.1761-0.079890.287117-0.0111450.16660.21434-0.135620.2690-0.2471180.0381460.01570.113350.17563-0.132910.151119-0.1671470.088-0.23536-0.19464-0.188920.0711200.2641480.0390.13837-0.00365-0.078930.002121-0.0961490.039100.29438-0.147660.03694-0.017122-0.0091500.0711-0.06939-0.205670.15395-0.1751230.121510.135120.17140-0.016680.09396-0.196124-0.031152-0.04513-0.226410.122690.36597-0.0281250.1041530.152140.03842-0.192700.065980.063126-0.061540.07150.0343-0.055710.279