数学建模--葡萄酒的分级(正式版)

承诺书

我们仔细阅读了中国大学生数学建模竞赛的竞赛规则.

我们完全明白，在竞赛开始后参赛队员不能以任何方式（包括电话、电子

邮件、网上咨询等）与队外的任何人（包括指导教师）研究、讨论与赛题有关

的问题。

我们知道，抄袭别人的成果是违反竞赛规则的, 如果引用别人的成果或其

他公开的资料（包括网上查到的资料），必须按照规定的参考文献的表述方式

在正文引用处和参考文献中明确列出。

我们郑重承诺，严格遵守竞赛规则，以保证竞赛的公正、公平性。如有违

反竞赛规则的行为，我们将受到严肃处理。

我们参赛选择的题号是（从A/B/C/D中选择一项填写）： A 我们的参赛报名号为（如果赛区设置报名号的话）：

所属学校（请填写完整的全名）：西安理工大学

参赛队员 (打印并签名) ：1. 郑晓东

2. 罗璐

3. 宫维静

指导教师或指导教师组负责人 (打印并签名)：

日期： 2013 年 05月

10 日

葡萄酒质量的综合评价分析

摘要

近年来，随着人们生活水平的提高，葡萄酒也随之受到人们的喜爱，加之食品科学技术的提高，人们对葡萄酒的品质也有了更高的要求，本文就针对葡萄酒品质的相关问题进行建模，求解和有关分析。

对问题一，首先基于两组评酒员对同一批葡萄酒的评价分数数据，采用假设检验中的t检验法建立评估两组数据差异的模型，运用Spss软件求解，得到两组数据存在显著性差异的结论，其次，通过计算两组数据的方差，用以比较稳定性，得到第二组更可信的结论。

对问题二，首先对酿酒葡萄理化指标数据进行标准化处理，经过主成分分析法将葡萄分为四个等级，其次，按可信度高的一组（第二组）得分将葡萄酒分为五级，综合两种分级，将酿酒葡萄分为了——级。

对问题三，首先同问题二对酿酒葡萄与葡萄酒的理化指标进行主成分分析，用Matlab的曲线拟合得到葡萄酒的得分,分别与酿酒葡萄与葡萄酒的理化指标之间的函数关系，再进行反解即得到酿酒葡萄与葡萄酒的理化指标之间联系。

对问题四，采用灰色关联度分析的方法进行求解，分别求出酿酒葡萄的理化指标与葡萄酒质量的关联度、葡萄酒理化指标与其质量的关联度，通过关联度值的大小，即可看出酿酒葡萄和葡萄酒的理化指标对葡萄酒质量的影响大小，并以此为基准来论证酿酒葡萄和葡萄酒的理化指标能否用来评价葡萄酒的质量。

关键词：t检验主成分分析曲线拟合灰色关联度分析

一、问题重述

确定葡萄酒质量时一般是通过聘请一批有资质的评酒员进行品评。每个评酒员在对葡萄酒进行品尝后对其分类指标打分，然后求和得到其总分，从而确定葡萄酒的质量。酿酒葡萄的好坏与所酿葡萄酒的质量有直接的关系，葡萄酒和酿酒葡萄检测的理化指标会在一定程度上反映葡萄酒和葡萄的质量。附件1给出了某一年份一些葡萄酒的评价结果，附件2和附件3分别给出了该年份这些葡萄酒的和酿酒葡萄的成分数据。请尝试建立数学模型讨论下列问题：

1. 分析附件1中两组评酒员的评价结果有无显著性差异，哪一组结果更可信？

2. 根据酿酒葡萄的理化指标和葡萄酒的质量对这些酿酒葡萄进行分级。

3. 分析酿酒葡萄与葡萄酒的理化指标之间的联系。

4．分析酿酒葡萄和葡萄酒的理化指标对葡萄酒质量的影响，并论证能否用葡萄和葡萄酒的理化指标来评价葡萄酒的质量？

二、问题分析

2.1 问题一

本题给出了两组评酒员对同一批葡萄酒的评价分数，在本文，采用假设检验中的t检验建立评估两组数据差异性的模型，研究两组评论员的评价是否存在差异，判断能否接受它们存在显著性差异的假设。若接受，则继续第二步：可靠性分析，分别对两组数据求方差，方差小的说明波动小，既评酒员的评价较稳定，可靠性高。

2.2 问题二

首先，我们利用问题一得到的结果，对可靠性高的一组数据进行处理，降低评论员之间的差异，提高葡萄酒样品最终得分的可靠度。按得分对葡萄酒进行分级。

然后，用标准化处理后的酿酒葡萄的理化指标对葡萄进行主成分分析。

最后，结合葡萄酒的分级对酿酒葡萄进行分级。

2.3 问题三

首先，用处理酿酒葡萄的理化指标的方法对葡萄酒的理化指标做同样的处理，得到葡萄酒理化指标的主成分。

然后，分别根据主成分获得红葡萄和红葡萄酒的的得分。通过曲线拟合，分别建立红葡萄得分和专家的评分之间的关系；红葡萄酒得分和专家评分之间的关系。

最后，根据两种理化指标和专家的评分之间的关系，建立两种理化指标之间的关系。

2.4 问题四

运用灰色关联度分析的方法，定量描述酿酒葡萄与葡萄酒的理化指标对葡萄酒质量的影响，以此为基准来论证酿酒葡萄和葡萄酒的理化指标能否用来评价葡萄酒的质量。

三、问题假设

1．同种葡萄酒在同一组评酒员的得分下成正态分布。 2．一种葡萄对应酿制一种葡萄酒。

3．葡萄的成分充分转化为葡萄酒里的成分，不存在意外的浪费和挥发。

四、符号说明

这里只列出主模型的全局参数，其他局部参数见具体模型。

(1)i J ：第i 个红葡萄酒样品

(1)

ij a ：第i 个红葡萄酒样品的得分

T1：第一组评酒员全体 T2：第二组评酒员全体

五、模型的建立与求解

5.1 模型一：基于t 检验建立差异评估模型

我们采用假设性检验验证是否能接受两组评酒员的评价结果存在显著性差异的假设。然后用方差分析两组评酒员评价数据的波动，认为较平稳的一组比较可靠。

5.1.1、数据预处理

我们在分析数据是发现了几个显著性的异常数据：

第一组红酒数据——样品20——色调——评酒员4号 数据缺失 第一组白酒数据——样品3——持久性——评酒员7号 怀疑多了一个7 第一组白酒数据——样品8——口感分析——评酒员2号 数据明显异常 因为随机样本在均值附近振荡，所以我们选用均值来代替异常数据以求误差最小。

5.1.2 t 检验模型的建立

21,T T 分别代表第一，第二组整体，分别对红葡萄酒i R （i=1,2， (27)

和白葡萄酒i W

（i=1,2,…,27）进行感官评价，1T 的评价结果通过组内的每一评酒员的评分的均值来表示。同样的，T2的评价结果通过组内的每一评酒员的评分的均值来表示。

从而得到两组评论员分别对红葡萄酒的评价结果见表一：

表1 红葡萄酒的评价结果

表中对于同一酒样品的两个评价差异是由两个评酒员引起的，为鉴定他们的评价结果有无显著性差异，可对两组数据对同一样品的差值进行分析，既表中的D 。

以红葡萄酒为例：有27对相互独立的评价结果（X1，Y1）（X2，Y2）...（X27,Y27）,D1=X1-Y1,D2=X2-Y2,...,D27=X27-Y27，由于Di （i=1,2， (27)

是由同一因素造成的，可认为它们服从同一分布。现假设Di~N （D μ，2D σ），

i=1,2…,27,且D μ，2

D σ未知，基于这一样本检验假设：

0:,0:100≠=D H H μμ （1）

分别记1227

,,,D D D 的样本均值和样本方差的观测值为d ，2

D s 。对

1227

,,,D D D 进行单个均值的t 检验，检验问题的拒绝域为（显著水平为

α）：

(1)t t n α=

≥-. （2）

当t 的值不落在拒绝域内时，接受0H ，既认为两组评价无显著性差异。否则，两组评价有显著性差异。

对白葡萄酒的处理同红葡萄酒。 2）模型的求解

现以红葡萄酒为例求解，首先，作出同一酒样品(1)i J (1,2,,27)i =分别由两组品酒员1T 、2T 得到的评价结果之差,列于表一的第三行，根据建立的模型检验假设：

01:0,:0D D H H μμ=≠. （3）

取α=0.05，运用spss 软件求解得到表二：

表二 t 检验求解结果

根据上表得到的Sig=0.020<0.05,所以拒绝接受，即认为两组品酒员

的评价结果有显著性差异。 5.1.3 可信度定量分析

记第一组10位品酒员对红葡萄酒样品

(1)i J (1,2,,27)i =的评分为

(1)ij a (1,2,

,10)j =

10(1)

(1)1

110i

ij j a a ==∑，10(1)(1)2

(1)2111()10i ij i j s a a ==-∑ （4）

其中，(1)

i a 表示第一组品酒员对红葡萄酒样品(1)i J 的评分均值，(1)2

1i s 表示

(1)i J 的评分方差；同样，第二组对红葡萄酒样品(1)i J 的评分均值和方差分别为

10(1)

(1)1

110i

ij j c c ==∑，10(1)(1)2

(1)2211()10i ij i j s c c ==-∑ （5）

从而对每一组品酒员得到一个评分方差向量

(1)2(1)2(1)2(1)211112127(,,

,)S s s s = (1)2

(1)2(1)2(1)222122227(,,

,)S s s s =

同理可求得白葡萄酒的(2)21S ，(2)22S 。再对(1)21S 和(1)22S 中的元素分别求和得

到总方差，对于同一批红葡萄酒用总方差来代表两组不同的评价水平。总方差

小的稳定性好，评价结果是更可信的。

运用excel 软件可以求解得到(1)21S ，(1)22S ，(2)21S 和(2)22S 。得到(1)21S = 1415.513，(1)22S =821.11，(2)21S =2986.77，(2)2

2S =1411.69。不管是红葡萄酒还是

白葡萄酒，第一组的总方差总是远远大于第二组。说明第二组的评价结果更为可信。

5.2 模型二：

对于问题二，是要基于酿酒葡萄的理化性质和葡萄酒的质量对酿酒葡萄进行分级，因此，对于模型二可分为三步进行，即：

1) 根据酿酒葡萄的理化指标对酿酒葡萄进行分级； 2) 根据评酒师的评分对葡萄酒的质量进行分级； 3) 综合两种因素，对酿酒葡萄进行分级。

5.2.1 根据酿酒葡萄的理化指标对酿酒葡萄进行分级

根据附录给出的酿酒葡萄的理化指标，可以看出，有些理化指标含量很低，有些理化指标含量很高。所以对于此种情况，我们采用主成分分析法对附录中的理化指标进行处理，将理化指标分为几种主成分，然后根据主成分对酿酒葡萄进行打分，通过得分对酿酒葡萄进行分级。

5.2.1.1 对于不同的理化指标可能存在着不同的量纲，因此在进行主成分分析之前应对酿酒葡萄的理化指标进行标准化处理。处理方法如下： 将原始数据标准化，即做如下数据变换：

（6）

其中 ， ，j = 1，2，…，p 。标准化后的数 据阵记为X *，其中每个列向量（标准化变量）的均值为0，标准差为1，数据无量纲。

标准化后变量的协方差矩阵（Covariance Matrix ）Σ = (s ij )p ?p ，即原变量的相关系数矩阵（Correlation Matrix ）R= (r ij )p ?p ：

i ，j = 1，2，…，p （7） 此时n 个样品在m 个主成分上的得分应为：

F j = a 1j X 1* + a 2j X 2* +...+ a pj X p * j = 1，2，…，m （8）

*1,2,...,;1,2,...,ij j

ij j x x x i n j p

s -===1

1n

j ij i x x n ==∑221

1()1n j ij j i s x x n ==--∑**1

1

()()

1

1

11n

ki

i kj j n

n

ij ki kj ij

k k x

x x x x x x

s x x r n n ==---===

=--∑∑

5.2.1.2 主成分分析法的步骤如下： 步骤一：计算协方差矩阵

计算样品数据的协方差矩阵：Σ = (s ij )p ?p ，其中

i ，j = 1，2，…，p （9）

步骤二：求出Σ的特征值及相应的特征向量

求出协方差矩阵Σ的特征值λ1≥λ2≥…λp >0及相应的正交化单位特征向量：

则X 的第i 个主成分为F i = a i 'X i = 1，2，…，p 。 步骤三：选择主成分

在已确定的全部p 个主成分中合理选择m 个来实现最终的评价分析。一般用方差贡献率 （10） 解释主成分F i 所反映的信息量的大小，m 的确定以累计贡献率 （11）

达到足够大（一般在85%以上）为原则。 步骤四：计算主成分得分

计算n 个样品在m 个主成分上的得分：

i=1,2,3,……，m （12）

1

1()()

1

n

ij ki

i kj j k s x x x x n ==

---∑?

?

????

?

??=??????? ??=??????? ??=pp p p p p p a a a a a a a a a a a a ......,,...,...212221*********/p

i i k

k αλλ==∑1

1

()/p

m

i k i k G m λλ===∑∑1122...i i i pi p

F a X a X a X =+++

5.2.1.3 模型的求解

利用MATLAB软件编程，对酿酒葡萄的理化指标进行主成分分析（以红葡萄的指标为例），根据累计贡献度大于85%的原则筛选，得到的前12个特征值及其贡献度率如表所示：

表三酿酒红葡萄理化指标的主成分分析结果

序号特征值贡献度

1 9.79705 0.1661

2 8.02727 0.1361

3 6.7966

4 0.1152

4 5.3110

5 0.0900

5 5.15878 0.0874

6 3.58111 0.0607

7 2.7598 0.0468

8 2.42551 0.0411

9 2.25743 0.0383

10 1.9615 0.0332

11 1.58293 0.0268

12 1.47825 0.0251

得分表达式为：

(13)

得分如下所示：

品种123456789得分-2.1822-0.36180.1567 1.2699-0.0978-0.33880.2743-1.259-0.4818

品种101112131415161718得分-0.5547-1.24770.8408-0.101-0.2899-0.7194-0.0704-0.1212-1.0973

品种192021222324252627得分0.4220.5899-0.3420.70010.009-0.12480.15270.02040.2406

表四酿酒葡萄的得分

根据得分的大小对酿酒葡萄进行分级

中上级酿酒葡萄：4,12,22,20

中级酿酒葡萄：19,7,27,3,25,26,23,16,5,13,17,24,14 中下级酿酒葡萄：6,21,2,9,10,15, 18

下级酿酒葡萄：11,8,1

5.2.2 根据评酒师的评分对葡萄酒的质量进行分级

对于“根据评酒师的评分对葡萄酒的质量进行分级”这一问题，我们认为品质优良的葡萄酿出来的葡萄酒的品质也应是优良的。它们之间存在着一一对应关系。所以可以通过专家评委们的打分对葡萄酒进行分级。

根据我们对问题一模型的求解，得知第二组评委的评分更可信，所以针对于第二组评委对每种葡萄酒的综合得分的平均值进行排序，如下所示：

葡萄酒样品的分级标准：

80~85分：高级葡萄酒

75~80分：中上级葡萄酒

70~75分：中级葡萄酒

65~70分：中下级葡萄酒

60~65分：下级葡萄酒

所以采用上面的分级标准，可将27中葡萄酒分为如下级别：

中上级葡萄酒：9，23,20

中级葡萄酒：3,17，2,14,19,21,5,26，22,24,27,4

中下级葡萄酒：16,10,13,12,25,1,6,8,15,18,7

下级葡萄酒：11

5.2.3 综合两种因素，对酿酒葡萄进行分级

综合两种因素分级的结果，对酿酒葡萄进行综合评价，并得到酿酒葡萄的分级。具体步骤如下：

1.分别对两种分级结果进行编号：

1)中上级葡萄酒编号为0

2)中级葡萄酒编号为1

3)中下级葡萄酒编号为2

4)下级葡萄酒编号为3

对酿酒葡萄的分级结果编号同上。

2.计算同一样品在两种不同情况下的编号与编号0的差值，分别记为

，，并计算总差值=+。

3.根据总差值的大小对酿酒葡萄进行分级。

酿酒葡萄的分级标准：

0~1：中上级酿酒葡萄

2：中级酿酒葡萄

3~5：中下级酿酒葡萄

6：下级酿酒葡萄

分级结果为：

中上级酿酒葡萄：20,23,22,4

中级酿酒葡萄：12,9,3,19,17,27,24,26,5,14，

中下级酿酒葡萄：7,25,16,13,6,21,2,10,15,1811,8

下级酿酒葡萄：1

5.3 模型三

对于问题三，我们同样采用主成分分析法，得到葡萄酒的主成分记为G

i i=1,2，n，n为主成分的个数。

应用Matlab软件编程，对葡萄酒的主成分进行求解（以红葡萄酒为例）根据累计贡献度大于85%的原则筛选，得到的前5个特征值及其贡献度率如表所示：

表六红葡萄酒理化指标的主成分分析结果

序号特征值贡献度

1 5.99601 0.3997

2 3.18036 0.2120

3 1.7651

4 0.1177

4 1.42734 0.0952

5 0.881763 0.0588

得分表达式：

得分如下所示：

品种123456789得分-0.613-1.2807-2.0259-0.6906-1.007-0.8185-0.4586-1.6665-1.989品种101112131415161718得分-1.79260.0572-0.2931-1.1967-1.1978-0.6577-0.5249-0.9953-0.5421品种192021222324252627得分-1.6224-0.8721-1.6659-0.9015-2.777-0.9881-1.0773-0.2149-0.8774

经分析直接求解葡萄酒理化指标与酿酒葡萄理化指标之间的联系难度较大，但我们发现它们都可与葡萄酒质量建立数学关系，将葡萄酒质量作为因变量，葡萄酒理化指标与酿酒葡萄理化指标的主成分分别作为自变量，采用曲线拟合，即可分别得到两理化指标与葡萄酒质量间的数学关系，再反解出两指标间的数学关系，这样不但简化了求解过程而且可减小误差。

根据曲线拟合的方法对酿酒葡萄理化指标的主成分和评委的得分之间建立函数关系为

y=0.62598x^8+1.7747x^7-0.74728x^6-4.457x^5- 0.70899x^4+2.8925x^3+0.34017x^2-0.36971x+3.5688

-2.5-2-1.5-1-0.500.51 1.5

2.6

2.83

3.23.43.63.84

4.2

图一 酿酒葡萄理化指标的主成分与评委得分的关系图

根据曲线拟合的方法对葡萄酒理化指标的主成分和评委的得分之间建立函数关系为

y=-4550.179x^7-6212.6221x^6-2244.2726x^5+120.8559x^4+160.1517x^3+4.9978x^2-3.2929x+3.5441

-0.8

-0.6

-0.4

-0.2

0.2

0.4

0.6

33.13.23.33.43.53.63.73.83.9

4

图二 葡萄酒理化指标的主成分和评委的得分的关系图

-1.5

-1-0.500.51 1.5

01

2

3

4

5

6

24

图三 酿酒葡萄理化指标的主成分和葡萄酒理化指标的主成分的关系图

5.4 模型四：运用灰色关联度求解相互联系

5.4.1问题分析

该问题需要分析酿酒葡萄和葡萄酒的理化指标对葡萄酒质量的影响，为方便运算，我们分别考虑酿酒葡萄的理化指标对葡萄酒质量的影响与葡萄酒的理化指标对酿酒葡萄的影响，为此我们引入灰色关联度的概念。灰色关联度能够定量描述事物或因素相互变化的情况，即变化的大小，方向和速度方面的关系。关联度越大，代表相互间联系越紧密，变化态势越一致，反之，若变化态势越不一致，则关联度越小。

5.4.2模型的建立与求解

灰色关联度方法的计算介绍：

1)原始数据的标准化

由于原始数据存在数量级和量纲的差异性，所以先进行标准化而方便进行计算。

=，= k=1，2，3……n

2)关联度的计算

经数据处理后的参考数列为:

比较数列为;

i=1,2,3……m

从几何角度看，关联程度实质是参考数列与比较数列曲线形状的相似程度，参考数列与比较数列曲线形状接近，则两者关联度较大；反之参考数列与比较数列曲线形状相差较大，则两者间的关联度较小。因此，可用曲线间的差值大小作为关联度的衡量指标。

则：

k=1，2，3，……，n

两极最大值和最小值：

=，=

关联系数：

=，，k=1，2，…，n，i=1，2，…，m

式中为分辨系数，用来削弱过大而使关联系数失真的影响。人为引入这个系数是为了提高关联系数之间的差异显著性。

3)关联度的计算与比较

由于每个比较数列与参数数列的关联程度是通过n个关联系数来反映的，关联系数分散，不便于从整体上进行比较。因此，有必要对关联信息做集中处理，求平均值则是一种比较好的信息集中的方式。即用比较数列与参数数列各个时期的关联系数的平均值来定量反映两个数列的关联程度，其计算公式为：

接下来用灰色关联度进行分析求解：

（1）对数据进行标准化处理

葡萄酒质量花色苷单宁总酚酒总黄酮白藜芦醇DPPH L*(D65)a*（D65）b*（D65H(D65)C(D65)

0.97593.6903 1.518 1.5934 1.6380.6721.60350.06040.31960.1757 1.67010.2982

1.09361.9613 1.525 1.5259

2.716 1.0052.05890.34710.9086 1.08960.70750.9311 1.05821.5111 1.825 1.3646 1.504 1.4451.77140.39890.9537 1.24810.63520.9972 0.92620.69540.8920.95480.8790.8080.79111.0295 1.1818 1.21140.8511 1.1751 0.966 1.06170.8050.96310.744 1.3760.9250.8387 1.1943 1.08920.9715 1.1666

0.99150.4434 1.0120.93510.908 1.2210.94521.3861 1.0805 1.06740.8879 1.068

1.05670.34420.5520.61580.5650.5020.5008 1.4360.9692 1.45240.5323 1.0517

0.92773.4812 1.655 1.6181 1.5820.28 1.54940.20930.77140.6648 1.03440.748

1.07521.4694 1.78 1.8057

2.023 1.0631.72530.34490.915 1.09550.70890.9372 0.96880.52560.7660.69320.6420.8940.60921.3895 1.15260.3623 2.96 1.0553 0.87380.04490.6320.64210.4290.1050.4703 2.1610.2410.8849-0.01590.4142

0.93190.31860.8890.76890.610.5960.63051.3065 1.0015 1.38530.5919 1.0623

1.01420.75820.8790.7870.8080.3690.74471.0123 1.16590.8876 1.1839 1.1148 1.02130.95330.8360.80010.6270.5970.72850.5895 1.1151 1.59870.5648 1.1945 1.02980.46450.5480.64870.3750.2450.30521.2887 1.14880.8645 1.1988 1.0972 1.02980.64990.6650.64540.5450.320.52371.2284 1.18020.8242 1.2986 1.1195 0.94040.8883 1.2620.9845 1.0030.4551.3875 1.003 1.1123 1.13760.8533 1.1056

0.96740.27250.6120.69480.7210.4790.6185 1.416 1.0863 1.02120.939 1.0657

1.00990.75260.8230.82320.791

2.4870.7455 1.161 1.2890.9361 1.2455 1.2269 1.02270.28180.8070.77540.8260.2660.705 1.91010.52390.71870.59820.5545 1.0496 1.189 1.389 1.42720.907 2.4221.60010.5233 1.0482 1.59460.5214 1.1427 1.10920.95120.9780.9894 1.19 1.23 1.03370.9869 1.073 1.18650.7794 1.0815 1.02411.5686 1.499 1.9999 2.48

3.4932.53380.35530.9303 1.13540.69260.9569 1.01421.02350.7910.8610.762 1.8920.738 1.0427 1.17250.8007 1.3299 1.1101 1.01560.60090.7440.70630.6170.710.73841.2228 1.26620.5222 2.2457 1.167

（2）计算参考数列与比较数列在同一时期的绝对差，得到下表：

花色苷单宁总酚酒总黄酮白藜芦醇 DPPH L*(D65) a*(D65) b*(D65) H(D65) C(D65) 2.7144 0.5420 0.6175 0.6617 0.3043 0.6276 0.9155 0.6563 0.8002 0.6942 0.6777 0.8677 0.4310 0.4323 1.6222 0.0886 0.9653 0.7465 0.1850 0.0040 0.3861 0.1625 0.4529 0.7666 0.3064 0.4462 0.3867 0.7132 0.6593 0.1045 0.1899 0.4230 0.0610 0.2308 0.0347 0.0286 0.0469 0.1181 0.1351 0.1033 0.2556 0.2852 0.0751 0.2489 0.0957 0.1610 0.0029 0.2220 0.4104 0.0410 0.1273 0.2283 0.1232 0.0055 0.2006 0.5481 0.0206 0.0564 0.0839 0.2291 0.0463 0.3946 0.0890 0.0759 0.1036 0.0765 0.7125 0.5043 0.4409 0.4920 0.5552 0.5559 0.3793 0.0875 0.3957 0.5244 0.0050 2.5535 0.7276 0.6904 0.6543 0.6479 0.6217 0.7184 0.1563 0.2629 0.1067 0.1797 0.3942 0.7047 0.7305 0.9473 0.0120 0.6501 0.7303 0.1602 0.0203 0.3663 0.1380 0.4432 0.2026 0.2756 0.3265 0.0747 0.3596 0.4207 0.1838 0.6065 1.9912 0.0865 0.8289 0.2423 0.2317 0.4444 0.7687 0.4035 1.2872 0.6328 0.0111 0.8897 0.4596 0.6133 0.0431 0.1630 0.3221 0.3361 0.3014 0.3746 0.0696 0.4534 0.3400 0.1304 0.2560 0.1355 0.2272 0.2062 0.6454 0.2695 0.0019 0.1517 0.1266 0.1697 0.1006 0.0680 0.1855 0.2212 0.3948 0.4247 0.2928 0.4318 0.0938 0.5774 0.4565 0.1732 0.5653 0.4814 0.3811 0.6550 0.7850 0.7246 0.2589 0.1190 0.1653 0.1690 0.0674 0.3799 0.3648 0.3844 0.4851 0.7097 0.5061 0.1986 0.1504 0.2056 0.2688 0.0897 0.0521 0.3216 0.0441 0.0627 0.4858 0.4471 0.0626 0.1719 0.1972 0.0871 0.1652 0.6949 0.3554 0.2726 0.2465 0.4882 0.3489 0.4486 0.1189 0.0538 0.0284 0.0983 0.2573 0.1868 0.1867 0.2187 1.4766 0.2644 0.1511 0.2791 0.0738 0.2356 0.2170 0.7409 0.2157 0.2473 0.1969 0.7571 0.3177 0.8874 0.4988 0.3040 0.4245 0.4682 0.1394 0.3390 0.3776 0.1430 1.3727 0.5505 0.5263 0.0014 0.5450 0.5282 0.0931 0.1580 0.1313 0.1198 0.0806 0.1211 0.0755 0.1223 0.0362 0.0773 0.3298 0.0277 0.5445 0.4744 0.9758 1.4557 2.4693 1.5097 0.6688 0.0938 0.1113 0.3315 0.0672 0.0093 0.2233 0.1532 0.2524 0.8779 0.2762 0.0285 0.1583 0.2135 0.3157 0.0959 0.4147 0.2716 0.3093 0.3985 0.3052 0.2772 0.2072 0.2506 0.4934 1.2301 0.1514 0.3622 0.4387 0.3155 0.4963 0.1823 0.5962 0.1209 0.2956 0.3853 0.1287 0.2984

表九参考数列与比较数列的绝对差

（3）通过上表可以计算出两极最大差与最小差，从而计算出关联系数，通过查阅文献得到得到关联系数表如下：

葡萄酒样品花色苷单宁总酚酒总黄酮白藜芦醇 DPPH L*(D65) a*(D65) b*(D65) H(D65) C(D65)

1 0.1671 0.5017 0.4691 0.4518 0.6425 0.4650 0.373

2 0.4539 0.4052 0.4400 0.4459

2 0.3859 0.5589 0.5581 0.2514 0.8619 0.3609 0.4221 0.7478 0.9952 0.5859 0.7716

3 0.5466 0.4156 0.6409 0.5503 0.5855 0.4333 0.4527 0.8407 0.7428 0.5635 0.9013

4 0.703

5 0.9423 0.9524 0.9229 0.8234 0.8028 0.8423 0.681

6 0.6573 0.880

7 0.6874

5 0.8523 0.7733 0.9973 0.711

6 0.5710 0.9322 0.8121 0.7058 0.8171 0.9925 0.7321

6 0.4989 0.9659 0.9082 0.8684 0.7050 0.9238 0.5806 0.8614 0.8796 0.8419 0.8787

7 0.4336 0.5198 0.5533 0.5259 0.4957 0.4953 0.5902 0.8634 0.5799 0.5100 0.9934

8 0.1758 0.4284 0.4413 0.4546 0.4571 0.4674 0.4315 0.7785 0.6755 0.8379 0.7532

9 0.5808 0.4363 0.4274 0.3652 0.9809 0.4562 0.4275 0.7741 0.9664 0.5986 0.7994

10 0.5520 0.7301 0.6650 0.6261 0.8813 0.6031 0.5649 0.7490 0.4735 0.2148 0.8648

11 0.3968 0.6932 0.7027 0.5513 0.4150 0.5751 0.2974 0.4629 0.9825 0.3799 0.5429

12 0.4708 0.9288 0.7711 0.6292 0.6192 0.6447 0.5932 0.8886 0.5463 0.6165 0.8084

13 0.6813 0.8023 0.7068 0.7266 0.4580 0.6700 0.9991 0.7836 0.8130 0.7638 0.8458

14 0.8910 0.7472 0.7123 0.5805 0.5625 0.6513 0.5584 0.8549 0.4858 0.5446 0.7601

15 0.4911 0.5314 0.5891 0.4544 0.4099 0.4294 0.6788 0.8223 0.7686 0.7646 0.8919

16 0.5898 0.5996 0.5870 0.5295 0.4345 0.5189 0.7340 0.7851 0.7272 0.6706 0.8604

17 0.9148 0.6296 0.9273 0.8988 0.5291 0.5498 0.8989 0.7615 0.7354 0.8640 0.7687

18 0.4397 0.6059 0.6674 0.6895 0.5279 0.6103 0.5490 0.8224 0.9122 0.9527 0.8489

19 0.6802 0.7459 0.7460 0.7147 0.2695 0.6742 0.7843 0.6622 0.8826 0.6992 0.7163

20 0.4240 0.7175 0.6888 0.7357 0.4187 0.6325 0.3805 0.5225 0.6427 0.5626 0.5383

21 0.7977 0.6172 0.5913 0.7935 0.2841 0.4978 0.5091 1.0000 0.5003 0.5082 0.8558

22 0.7766 0.8073 0.8213 0.8730 0.8197 0.8802 0.8182 0.9399 0.8776 0.6237 0.9539

23 0.5005 0.5350 0.3584 0.2723 0.1807 0.2652 0.4492 0.8549 0.8320 0.6225 0.8921

24 0.9857 0.7104 0.7819 0.6844 0.3831 0.6645 0.9526 0.7762 0.7196 0.6339 0.8521

25 0.5684 0.6683 0.6387 0.5782 0.6418 0.6637 0.7256 0.6859 0.5252 0.3070 0.7839

26 0.6014 0.5545 0.6341 0.5238 0.7505 0.4778 0.8200 0.6491 0.5864 0.8104 0.6470

27 0.5475 0.7793 0.7133 0.6417 0.6168 0.6464 0.4931 0.9871 0.6846 0.7037 0.9142

表十关联系数计算表

（4）计算关联度并对关联度进行分析

由表十得到了关联系数表，对各个时期的关联系数求平均值即得到了各个因素与葡萄酒的关联度如下表：

表十一关联度表

由上表可以得到，单宁，总酚，a*(D65)，b*(D65)，C(D65)对葡萄酒的质量有较大的影响，而H(D65)，酒总黄酮等对葡萄酒质量的影响则次之。

用同样的方法可以得到葡萄的理化指标与葡萄酒质量的关联度：

由上表可以看出乙醇，乙酸乙酯，单宁，总酚，a*(D65),b*(D65),C(D65)对葡萄酒质量的影响较大，而苯乙醇和酒总黄酮对葡萄酒质量的影响次之。

综上所述，虽然酿酒葡萄和葡萄酒的理化指标对葡萄酒的质量有影响，但是关联度都在0.5—0.7之间，所以不能用单独用理化指标对葡萄酒的质量进行评价，还要结合其他指标对葡萄酒的质量进行分析。

六、模型的评价

6.1模型的优点

1)主成分分析法可以很好的解决决定因素多但是又无法全部考虑的问题，很

大程度降低了计算量与计算难度。

2)灰色关联度分析法能够把两个变量之间的关联程度量化，可以很好的解决

决定因素和事物之间相互联系的问题。

6.2模型的缺点

1)本模型没有把葡萄酒具体的质量指标如：香气，气味等进行对应理化指标

的分析。

七、参考文献

（1）田民，刘思峰，朴志坤，灰色关联度算法模型的关联综述，《统计与决策》，卷期号：25—27，2008

（2）杨小平，刘刚，赵晋，《统计分析方法与SPSS应用程序》，北京：清华大学出版社，2008

（3）刘伟江，杨友社，寇光兴，《概率论与数理统计》，北京：清华大学出版社，2004

（4）曾建军，李世行，王永国，叶仁玉，夏惠异，《MATAB语言与数学建模》，安徽：安徽大学出版社