基于理化指标的葡萄酒评价分析模型

1、 . . 38/41基于理化指标的葡萄酒评价分析模型 日期：2012高教社杯全国大学生数学建模竞赛编 号 专 用 页赛区评阅编号（由赛区组委会评阅前进行编号）：赛区评阅记录（可供赛区评阅时使用）：评阅人评分备注全国统一编号（由赛区组委会送交全国前编号）：全国评阅编号（由全国组委会评阅前进行编号）：基于理化指标的葡萄酒评价分析模型摘要本文在利用主成分分析法提取出主成分理化指标后，结合偏最小二乘回归分析法建立了对葡萄和葡萄酒与其理化指标的相关分析和评价模型。 对问题1，在对附件1表格中的酒样品序号按由小到大排序后，直接利用Excel中的函数命令，按每个评酒员对各项的打分进行求和，得到了每种酒样品的

2、总评分。最后又对10名评酒员的各自总评分进行均值处理，得到了每种酒样品的平均得分和平均得分方差。从平均得分和平均得分方差的比较，不同葡萄酒平均得分的排序和图表三方面得出了两组结果均存在明显差异。后又通过方差和总体分布检验，结合图表得出的第二组结果更可信，并给出了红白葡萄酒的质量排序。在问题2中，首先将芳香物质视为葡萄的理化指标。再通过网上查询对一些重要的理化指标进行了预先排序，后对所有理化指标实施了主成分分析法，分别提取出了9种和10种主成分红白葡萄理化指标，累积贡献率均达到了85%，并给出了对红白葡萄的综合评价值与其排序，详见表7和表10。在问题3中，首先通过主成分分析法提取出了红白葡萄酒中

3、的3种和4种主成分理化指标，与问题2中的红白葡萄主成分理化指标相互讨论，利用偏最小二乘回归分析法建立了两者之间联系关系式，并通过回归系数直方图和预测图进一步说明了建立的联系。下面是红葡萄与红葡萄酒理化指标的回归方程在问题4中，还是利用主成分分析法把得出的红白葡萄和葡萄酒的主成分理化指标，与可信度强的第二组红白葡萄酒的质量平均得分，视为自变量和因变量，通过偏最小二乘回归分析法，建立了两者之间的影响关系式，并结合回归系数直方图和预测图得出能够用葡萄和葡萄酒的理化指标来评价葡萄酒质量的结论。关键词：方差 总体分布检验 主成分分析法 偏最小二乘回归分析法1 问题重述确定葡萄酒质量时一般是通过聘请一批有

4、资质的评酒员进行品评。每个评酒员在对葡萄酒进行品尝后对其分类指标打分，然后求和得到其总分，从而确定葡萄酒的质量。酿酒葡萄的好坏与所酿葡萄酒的质量有直接的关系，葡萄酒和酿酒葡萄检测的理化指标会在一定程度上反映葡萄酒和葡萄的质量。附件1给出了某一年份一些葡萄酒的评价结果，附件2和附件3分别给出了该年份这些葡萄酒的和酿酒葡萄的成分数据。请尝试建立数学模型讨论下列问题：1. 分析附件1中两组评酒员的评价结果有无显著性差异，哪一组结果更可信？2. 根据酿酒葡萄的理化指标和葡萄酒的质量对这些酿酒葡萄进行分级。3. 分析酿酒葡萄与葡萄酒的理化指标之间的联系。4分析酿酒葡萄和葡萄酒的理化指标对葡萄酒质量的影响

5、，并论证能否用葡萄和葡萄酒的理化指标来评价葡萄酒的质量？2 问题的分析2.1 问题1对于第一问：要讨论两组评酒员评价结果的差异性，先对附件1中4个表格里的样酒按序号进行了由小到大的重新排序，因为它们排序都是打乱的，不方便实现差异性分析。其次，根据每个评酒员在打分后，求和得其总分来确定葡萄酒的质量得分。而葡萄酒在澄清度、色调等10个指标方面的总分又是一致的，在此直接利用Excel中函数命令按每个评酒员对各项的打分进行求和，得到了每种样酒的总评分。最后又对10名评酒员的各自总评分进行均值处理，得到了每种样酒的平均得分。由此还可求出10名评酒员各自总评分的方差，通过对平均得分和平均分方差的比较、排序

6、，结合图表可对差异性实现分析。对于第二问：两组评酒员评价结果的可信度分析。利用第一问得出的平均得分和平均得分方差，通过数据比较之后，会有一组平均得分的方差要大，说明它的偏离程度就大，可信度就降低。为进一步对偏离程度进行分析，可对27种红葡萄酒和28种白葡萄酒的每名评酒员得分进行总体分布检验。因为它们的酒质量是未知的，方差也是未知的，属于非参数检验。故可直接利用Matlab工具箱中的normfit命令进行非参数检验的总体分布检验，来确定哪一组评价结果偏离程度大，实现可信度分析。2.2 问题2首先我们也将芳香物质视为葡萄的理化理化指标。因为它是构成葡萄与葡萄酒质量的主要因素,决定着葡萄酒的风味和典

7、型性。题目要求根据酿酒葡萄的理化指标和葡萄酒的质量对这些酿酒葡萄进行分级。可考虑先从网上查询酿酒葡萄中的理化指标对葡萄酒的作用与影响，粗略了解一下各项指标在葡萄酒中的所占的比重，把比重明显靠前的放到前面，不重要靠后，主观的对所有理化指标进行排序。然后，根据所有理化指标的排序，通过编程，利用主成分分析法提取出贡献率大的n个理化指标，使其累积贡献率达到85%以上，之后就用这n个作为主成分进行分析。最终能得出红白葡萄综合评价值与其排序结果，来实现分级（在讨论白葡萄时，去掉花色苷理化指标，因白葡萄酒中不含有花色苷）。在运用主成分分析法前，得先对附件3中的葡萄酒样品进行重新排序，还可能得对一些理化指标的

8、单位进行统一换算处理。2.3 问题3由题意知，可用偏最小二乘回归分析法分析酿酒葡萄与葡萄酒的理化指标之间的联系。但由附件2知道，酿酒葡萄的理化指标含有28种，红白葡萄的理化指标各为7种和6种。可见要全部讨论两者理化指标之间的联系，数据处理较大，过程可能也会较复杂。又由问题2知，能通过主成分分析法提取出红白葡萄理化指标中的几个主成分，且其累积贡献率在85%以上，若用这几个主成分理化指标来代替红白葡萄的全部理化指标进行分析，在指标的个数上大大减少，并且主要信息却基本没丢，故方法是可行的。所以我们首先也要通过主成分分析法对红白葡萄酒的理化指标进行分析，提取出累积贡献率达到85%以上的几个理化指标，用

9、来代替全部指标。然后用这两部分主成分理化指标，利用偏最小二乘回归分析法可实现对酿酒葡萄与葡萄酒的理化指标之间的联系分析。2.4 问题4 第一问：从2.3可知，完全可用由主成分分析法提取出来的主成分理化指标来代替全部理化指标进行分析，故本问还用那些主成分理化指标来分析对葡萄酒质量的产生的影响。对红白葡萄酒质量的衡量，在此我们选取第一问中可信度较强的那一组结果。这样因变量确定了，自变量为葡萄和葡萄酒的主成分理化指标，故还可用偏最小二乘回归分析法实现对葡萄酒质量的影响分析，同时还给出了两都之间的一个关系式。第二问：利用第一问中给出的葡萄和葡萄酒的主成分理化指标对葡萄酒质量的影响关系式，结合回归系数直

10、方图和预测图，完全可判断出能否用葡萄和葡萄酒的理化指标来评价葡萄酒的质量。3 模型假设1.每个评酒员在对葡萄酒进行品尝后的分类打分都是客观的；2.在第一组白葡萄酒的品尝评分中，评酒员7和9分别对酒样3和8中口感分析的持久性（总分8）给出了77和16的错误打分，我们视为打印错误，根据均值改为7和6；3.在评酒员评分过程中，设葡萄酒在外观、香气、口感和平衡/整体评价四方面共10个指标的总分是一致的；4.芳香物质是构成葡萄与葡萄酒质量的主要因素,决定着葡萄酒的风味和典型性，故将葡萄和葡萄酒的芳香物质视为酿酒葡萄和葡萄酒的理化指标。5.设附件3葡萄和葡萄酒芳香物质的各项成分所给数据的单位均为mmol/

11、L；6.设附件2中理化指标的顺序是随机排列的，可以对其打乱重排；7.设附件2中的理化指标之间不存在可完全替代的两种理化指标，基本彼此线性无关。4 模型的建立与求解4.1两组评酒员评价结果的有无显著性差异分析和可信度分析首先，在附件1（葡萄酒品尝评分表）中，我们发现4个表格里的酒样排序都是打乱的，没有固定的次序，无法很好的实现差异性分析，故对酒样进行了序号由小到大的重新排序。其次，考虑到是根据每个评酒员在打分后，求和得其总分来确定葡萄酒的质量得分的。又葡萄酒在澄清度、色调等10个指标方面的总分是一致的，故没有必要再用各项的打分与总分的比值来进行计算。在此，我们直接利用Excel中函数命令按每个评

12、酒员对各项的打分进行求和，得到了每种样酒的总评分。后又对10名评酒员的各自总评分进行均值处理，得到了每种样酒的平均得分。最后，为了差异性的讨论又对10名评酒员的各自总评分求出了方差，同时对各项指标也进行了均值和方差处理，详见附件1（葡萄酒品尝评分表）中的4个表。4.1.1两组评酒员对红葡萄酒评价结果的有无显著性差异分析下面是两组评酒员对27种红葡萄酒样品的平均得分和平均分方差与它们的差（第一组减第二组）情况：红葡萄酒样品第一组第二组两组平均分之差两组平均分的方差之差平均得分平均分方差平均得分平均分方差162.792.9068.181.88-5.411.02280.339.797416.226.

13、323.57380.445.8274.630.715.815.11468.6108.0471.241.29-2.666.76573.362.0172.113.661.248.36672.259.7366.321.125.938.61771.5103.6165.362.686.240.93872.344.016665.116.3-21.10981.532.9478.225.733.37.211074.230.4068.836.185.4-5.781170.170.7761.638.048.532.721253.979.6668.325.12-14.454.531374.644.9368.815.

14、295.829.64147336.0072.623.160.412.841558.785.5765.741.34-744.221674.918.1069.920.105-2.001779.388.0174.59.174.878.841859.947.2165.450.27-5.5-3.061978.647.3872.655.166-7.782078.626.0475.839.072.8-13.022177.1116.1072.235.514.980.592277.250.6271.624.275.626.362385.632.4977.124.778.57.72247874.8971.510.

15、726.564.172569.264.6268.243.73120.892673.831.297241.561.8-10.27277349.7871.520.501.529.28表1 两组评酒员对红葡萄酒样品的平均得分和平均分方差与平均分之差和平均分方差之差表从表1的两组平均分之差列可以看出，只有5个是负值，并且还有13个数值都大于5，说明第一组评酒员的平均得分基本都要比第二组的平均得分要高，最多高出8.5分。第一组平均得分列的平均分73.06，比第二组的70.51也高出2.55分。最后的一列两组平均分方差之差中，也只有7个是负值，大部分方差之差都要高出20，显然第一组的平均得分偏离均值的程度

16、要大。它们的平均分方差的平均分差为58.62-33.79=24.83。为了更好说明两组评酒员在红葡萄酒样品中的差异，通过程序（附录1：ti1tu1.m），作出了两组评酒员的平均得分图，红圆代表第一组得分，蓝星代表第二组得分，见图1。图1 两组评酒员对红葡萄酒的平均得分图4.1.2两组评酒员对白葡萄酒评价结果的有无显著性差异分析下面是两组评酒员对28种白葡萄酒样品的平均得分和平均分方差与它们的差（第一组减第二组）情况：白葡萄酒样品第一组第二组两组的平均分之差两组平均分的方差之差平均得分平均分方差平均得分平均分方差18292.2277.925.884.166.34274.2201.0775.849

17、.07-1.6152.00378.366.4675.6142.492.7-76.03479.444.7176.942.102.52.61571126.4481.526.28-10.5100.17668.4162.7175.522.72-7.1139.99777.539.1774.242.183.3-3.01870.4163.3872.331.12-1.9132.26972.992.7780.4106.27-7.5-13.501074.3212.6879.870.40-5.5142.281172.3177.1271.487.820.989.301263.3115.7972.4140.04-9.1

18、-24.261365.9170.7773.946.77-8124.001472114.2277.115.88-5.198.341572.4131.6078.454.04-677.561674178.0067.382.236.795.771778.8144.1880.338.46-1.5105.721873.1156.5476.730.23-3.6126.311972.246.4076.426.04-4.220.362077.864.4076.650.041.214.362176.4172.7179.264.40-2.8108.312271138.6779.453.60-8.485.072375

19、.943.6677.411.60-1.532.062473.3111.1276.138.54-2.872.582577.133.8879.5106.50-2.4-72.622681.372.9074.3102.907-30.002764.8144.407735.56-12.2108.842881.380.4679.625.381.755.08表2 两组评酒员对白葡萄酒样品的平均得分和平均分方差与平均分之差和平均分方差之差表从表2的两组平均分之差列可以看出，只有9个是正值，负值中有10个数值小于-5，说明第一组评酒员的平均得分基本都要比第二组的平均得分低，最多低12.2分。第一组平均得分列的平均

20、分73.96，比第二组的76.53也低了2.57分。但最后的一列两组平均分方差之差中只有6个是负值，且大部分方差之差都要高出50，也是第一组的平均得分偏离均值的程度要大。它们的平均分方差的平均分差为117.80-56.02=61.78。为了更好说明两组评酒员在白葡萄酒样品中的差异，通过程序（附件2：ti1tu2.m），作出了两组评酒员的平均得分图，红圆代表第一组得分，蓝星代表第二组得分，见图2。图2 两组评酒员对白葡萄酒的平均得分图为了更进一步的说明两组评酒员对红白两种葡萄酒的评分差异，我们又根据它们的平均得分对它们进行了排序，结果见下面表3：27种红葡萄酒的排序28种白葡萄酒的排序第一组第二

21、组第一组第二组红葡萄酒2385.6红葡萄酒978.2白葡萄酒182白葡萄酒581.5红葡萄酒981.5红葡萄酒2377.1白葡萄酒2681.3白葡萄酒980.4红葡萄酒380.4红葡萄酒2075.8白葡萄酒2881.3白葡萄酒1780.3红葡萄酒280.3红葡萄酒374.6白葡萄酒479.4白葡萄酒1079.8红葡萄酒1779.3红葡萄酒1774.5白葡萄酒1778.8白葡萄酒2879.6红葡萄酒1978.6红葡萄酒274白葡萄酒378.3白葡萄酒2579.5红葡萄酒2078.6红葡萄酒1472.6白葡萄酒2077.8白葡萄酒2279.4红葡萄酒2478红葡萄酒1972.6白葡萄酒777.5白

22、葡萄酒2179.2红葡萄酒2277.2红葡萄酒2172.2白葡萄酒2577.1白葡萄酒1578.4红葡萄酒2177.1红葡萄酒572.1白葡萄酒2176.4白葡萄酒177.9红葡萄酒1674.9红葡萄酒2672白葡萄酒2375.9白葡萄酒2377.4红葡萄酒1374.6红葡萄酒2271.6白葡萄酒1074.3白葡萄酒1477.1红葡萄酒1074.2红葡萄酒2471.5白葡萄酒274.2白葡萄酒2777红葡萄酒2673.8红葡萄酒2771.5白葡萄酒1674白葡萄酒476.9红葡萄酒573.3红葡萄酒471.2白葡萄酒2473.3白葡萄酒1876.7红葡萄酒1473红葡萄酒1669.9白葡萄酒1

23、873.1白葡萄酒2076.6红葡萄酒2773红葡萄酒1068.8白葡萄酒972.9白葡萄酒1976.4红葡萄酒872.3红葡萄酒1368.8白葡萄酒1572.4白葡萄酒2476.1红葡萄酒672.2红葡萄酒1268.3白葡萄酒1172.3白葡萄酒275.8红葡萄酒771.5红葡萄酒2568.2白葡萄酒1972.2白葡萄酒375.6红葡萄酒1170.1红葡萄酒168.1白葡萄酒1472白葡萄酒675.5红葡萄酒2569.2红葡萄酒666.3白葡萄酒571白葡萄酒2674.3红葡萄酒468.6红葡萄酒866白葡萄酒2271白葡萄酒774.2红葡萄酒162.7红葡萄酒1565.7白葡萄酒870.4

24、白葡萄酒1373.9红葡萄酒1859.9红葡萄酒1865.4白葡萄酒668.4白葡萄酒1272.4红葡萄酒1558.7红葡萄酒765.3白葡萄酒1365.9白葡萄酒872.3红葡萄酒1253.9红葡萄酒1161.6白葡萄酒2764.8白葡萄酒1171.4白葡萄酒1263.3白葡萄酒1667.3表3 红白两种葡萄酒的排序4.1.3两组评酒员评价结果的可信度分析 由上面分析结果知道，两组评酒员的平均得分是存在明显差异的，并且不论红白葡萄酒第一组评酒员的评分结果方差都比第二组大，在这55种葡萄酒中竟多达42个，比例为76.36%，最大的平均分方差之差为152，说明第一组平均得分结果偏离均值的程度要大

25、。已可说明第二组的评价结果要比第一组的评价结果要可信。为了更好的说明问题，我们又用了非参数检验中的总体分布检验，对偏离程度进行了分析。如果只对每种葡萄酒样品的两组数据进行讨论，因数据太少结果不明显。在这里我们对27种红葡萄酒和28种白葡萄酒的每名评酒员得分进行讨论。因为它们的质量是未知的，肯定存在差异，所以方差也是未知的，属于非参数检验。为了简便，我们直接Matlab工具箱中的normfit命令进行总体分布检验，用程序（附录3：fenbu1.m，附录4：fenbu2.m）分别给出两组红、白葡萄酒样品的总体分布图，见下面图3，图4。图3 两组红葡萄酒样品的总体分布图图4 两组白葡萄酒样品的总体分

26、布图从上面两个图的下方4个线形图可以看出，数据基本都分布在一条直线上，说明两组评酒员的评价结果确实都是客观的。但上方的4个柱形图告诉我们，第二组的数据没有第一组集中，又知第一组红白葡萄酒评价结果在均值为0.95的置信区间为（71.83，74.28）、（72.64，75.31），第一组红白葡萄酒评价结果在均值为0.95的置信区间为（69.70，71.33）、（75.62，77.45），显然第一组的数据偏离程度大。进一步说明第二组的评价结果要比第一组的评价结果要可信。 4.2酿酒葡萄的分级方案4.2.1酿酒葡萄中的理化指标对葡萄酒的作用与影响芳香物质：是具挥发性的、能够产生一定气味的含香物质的总称

27、。葡萄与葡萄酒中芳香物质大部分都是结构简单的小分子有机物,含量少,种类多。它们是构成葡萄与葡萄酒质量的主要因素,决定着葡萄酒的风味和典型性，芳香物质的种类、含量、感觉阈值与其之间的相互作用决定着葡萄酒的感官质量，芳香物质是造就葡萄酒风味的物质之一，葡萄酒中的芳香物质越多，葡萄酒的风味就浓厚。氨基酸：在葡萄酒的酿制过程中，通过发酵菌的生成作用，由氨基酸的协助将葡萄中的糖分解转化成酒精。并在此过程中保留了原来的氨基酸并且有新的氨基酸，生成的这些氨基酸和多种微量元素都是生命必须的营养物质。蛋白质：蛋白质主要从葡萄籽、葡萄皮提取的。蛋白质在葡萄酒中起着多种作用：对酒石的稳定、对蛋白的稳定、芳香物质的支

28、撑体、稳定多酚物质。而多酚物质对红葡萄酒的感官质量起着决定性的作用。在酒精发酵过程中与当葡萄酒在酒泥上酿时，蛋白质可以影响葡萄酒的香气质量，改变芳香物质的挥发性，还可提高葡萄酒的酒石稳定性，同时可提高葡萄酒的“肥硕感”与香气的馥郁性。Vc：具有抗氧化的作用同时还具有促氧的作用，但是由于它会使葡萄酒产生沉淀物质，所以不宜使用。花色苷:主要是皮里的物质花色苷，在葡萄酒中起呈色作用，通过浸渍葡萄皮中的部分花色苷被转移到葡萄酒中，使得葡萄酒呈出紫红色。在酝酿过程中，生成稳定的色素，葡萄酒的颜色由紫色变为砖红，是红葡萄酒颜色的主要形式，白葡萄酒中不含有。单宁：主要是由葡萄籽、皮与梗浸泡发酵而来，或是因为

29、存于橡木桶而萃取橡木的单宁而来。它是葡萄酒中所含有的二种酚化合物其中的一种物质，尤其在红葡萄酒中含量较多，作用是使葡萄酒入口后口腔感觉干涩，口腔黏膜会有褶皱感。单宁的多少可以决定酒的风味、结构与质地。缺乏单宁的红酒质地轻薄，没有厚实的感觉。单宁丰富的红酒可以存放经年，并且逐渐酝酿出香醇细致的年风味。总酚：含有多种酚类成分，广泛分布在葡萄皮、梗和籽里面。这些成分在葡萄酒中起决定性的作用。首先，影响葡萄酒的口感，使葡萄酒具有苦味和收敛性，尤其是红酒；第二，影响葡萄酒的颜色；第三，起到防腐的作用，延长酒的寿命；第四，在葡萄酒中的抗氧化作用。白藜芦醇：主要存在于葡萄的果皮中，其次是种子。在葡萄中的含量

30、甚微，只有在酿造葡萄酒的发酵过程中与酶结合，才会使它的数量骤然增加。率能以游离态(顺式、反式) 和糖苷结合态(顺式、反式)2种形式在植物中分布与生物合成，且均具有抗氧化效能。葡萄酒中的白藜芦醇含量的高低已成为衡量优质酒和劣质酒的标准，但并不是所有的葡萄酒中都含有这种成分。总黄酮：主要是从葡萄皮和葡萄籽中提取的，在葡萄酒中的类黄酮类物质对葡萄酒的颜色、收敛性以与苦味等都具有重要的作用，还具有很强的抗氧化能力以与降血糖、抗病毒、抗肿瘤、免疫调节等作用,还可以作为食品添加剂来避免亚硝酸盐导致的致病隐患等作用。 酒石酸、苹果酸和柠檬酸：为葡萄果实中的主要有机酸，有近70的有机酸分布在中、果皮(即果肉)

31、中，其组成与含量是决定葡萄果实风味和品质的重要因素之一.有机酸的种类与其浓度调节着酸碱的平衡,必然要影响葡萄酒的酸味,进而也会影响葡萄酒的风味.此外,有机酸的防腐特性还可以增强葡萄酒的微生物和物化的稳定性,酸度较高的年限短的白葡萄酒更具有酿潜力. 因此它们对葡萄酒的组成成分、稳定性和感官品质都起着重要的作用.4.2.2主成分分析法的原理和步骤主成分分析（Principal Component Analysis）是利用降维的思想，将多个变量转化为少数几个综合变量（即主成分），其中每个主成分都是原始变量的线性组全，各主成分之间互不相关，从而这些主成分能够反映始变量的绝大部分信息，且所含的信息互不重

32、叠。采用这种方法可以从多达30多种酿酒葡萄的理化指标中，提取出贡献率大的决定性理化指标，然后对其分类排序给出分级，这样避免了都分析的复杂性，结果较为科学，信息也较为准确。下面为主成分分析法的评价步骤：（1）对原始数据进行标准化处理设所考察的n个量中含有m个指标变量（即成分），则由此可构造原始数据矩阵，（，），将各指标值转换成标准化指标，（，）其中，()，即为第个指标的样本均值和样本标准差。对应地，称，（）为标准化指标变量。（2）计算相关系数矩阵相关系数矩阵，（）式中，是第个指标与第个指标的相关系数。（3）计算特征值和特征向量计算相关系数矩阵的特征值，与对应的标准化特征向量，其中，由特征向量组成

33、m个新的指标变量式中是第1主成分，是第2主成分，是第m主成分。（4）选择（）个主成分，计算综合评价值 = 1 * roman i）计算特征值的信息贡献率和累积贡献率。称（）为主成分的信息贡献率；为主成分的累积贡献率，当接近于1（）时，则选择前个指标变量作为个主成分，代替原来m个指标变量，从而可对个主成分进行综合分析。 = 2 * roman ii）计算综合得分其中为第个主成分的信息贡献率，根据综合得分值就可进行评价。4.2.3利用主成分分析法对酿酒葡萄质量的分析与分级方案因为是从葡萄酒的质量角度出发对酿酒葡萄进行分级的，故根据5.2.1中理化指标对葡萄酒的作用与影响，我们先主观对一些理化指标进

34、行了排序，结果如下：12345678芳香物质氨基酸单宁总酚花色苷蛋白质酒石酸苹果酸910111213141516柠檬酸葡萄总黄酮白藜芦醇黄酮醇总糖还原糖可溶性固形物PH值表4 酿酒葡萄中理化指标的主观排序为了运用主成分分析法，先对数据进行初步处理。利用根据1摩尔=1分子量克，对葡萄和葡萄酒中的芳香物质成分进行单位统一化。对含有二级指标的，与一级指标求和再取平均值。有多次测试值的，也取平均值。下面通过主成分分析法分别对红白葡萄进行排序、分级。按表4中的理化指标排序，利用MATLAB软件编程序（附录5：fenjihong.m）,求得红葡萄相关系数矩阵的前10个特征根与其贡献率,见下表5。序号理化指

35、标特征值贡献率累计贡献率1芳香物质7.0322.6722.672氨基酸4.9816.0638.723单宁3.7712.1550.884总酚2.869.2460.125花色苷2.056.6166.736蛋白质1.745.6272.357酒石酸1.424.5876.948苹果酸1.284.1281.059柠檬酸1.113.5784.6310葡萄总黄酮0.892.8887.51表5 红葡萄分级的主成分分析结果可以看出，前9个特征根的累计贡献率就达到84.63%以上，已很接近85%，而增加一个成分贡献率也就增加2.88，故选前9个作为主成分进行分析。下面是前9个特征根对应的特征向量，见表6。第1特征向

36、量第2特征向量第3特征向量第4特征向量第5特征向量第6特征向量第7特征向量第8特征向量第9特征向量-0.1010.096-0.1020.1120.204-0.032-0.0530.1440.6350.1380.2460.0130.255-0.195-0.2250.146-0.021-0.0640.286-0.059-0.148-0.046-0.1170.1830.199-0.0300.113.-0.211-0.1540.151-0.0090.064-0.2360.0540.3210.029-0.1230.1110.3870.037-0.0200.2230.2130.186-0.077-0.05

37、30.2090.319-0.0080.1170.3500.2510.069-0.099表6 红葡萄主成分分析法的前9个特征根对就的特征向量由此可得9个主成分分别为. . . . . . . . . . . . . . . . . . . 分别以9个主成分的贡献率为权重，构建主成分综合评价模型：将红葡萄酒的9种主成分值代入上式，可以得到各种红葡萄的综合评价值与其排序结果，见下表7所示。名次123456789红葡萄样品3样品9样品11样品1样品2样品21样品23样品22样品8综合评价值1.581.171.061.060.990.950.340.280.28名次101112131415161718红

38、葡萄样品16样品12样品14样品6样品24样品17样品5样品13样品19综合评价值0.140.120.11-0.12-0.15-0.18-0.20-0.22-0.28名次192021222324252627红葡萄样品7样品18样品15样品4样品20样品27样品26样品10样品25综合评价值-0.30-0.33-0.40-0.65-0.69-0.80-1.01-1.17-1.58表7 红葡萄的综合评价值与其排序结果故按不含花色苷的表4中的理化指标排序（因白葡萄酒中不含有花色苷），利用MATLAB软件编程序（附录6：fenjibai.m）,求得白葡萄相关系数矩阵的前11个特征根与其贡献率,见下表8

39、。序号理化指标特征值贡献率累计贡献率1芳香物质5.72253619.7328819.73288312氨基酸4.76366716.4264436.15932013单宁3.57986412.3443648.50367894总酚2.1358427.36497355.86865185蛋白质1.9834646.83953262.70818376酒石酸1.6242075.60071568.30889887苹果酸1.4430474.97602473.28492258柠檬酸1.1832924.08031677.36523869葡萄总黄酮1.1655584.01916681.384404710白藜芦醇1.030

40、0233.55180484.936209111黄酮醇0.8156622.81262687.7488355表8 白葡萄分级的主成分分析结果可以看出，前10个特征根的累计贡献率就达到84.94%以上，约为85%，而增加一个成分贡献率也就增加2.81，故选前10个作为主成分进行分析。下面是前10个特征根对应的特征向量，见表9。第1特征向量第2特征向量第3特征向量第4特征向量第5特征向量第6特征向量第7特征向量第8特征向量第9特征向量第10特征向量0.02-0.050.02-0.51-0.04-0.25-0.130.37-0.060.040.240.110.070.170.340.02-0.09-0.

41、070.210.070.170.210.12-0.04-0.040.37-0.01-0.18-0.35-0.18.0.190.20-0.270.15-0.16-0.08-0.230.000.08-0.22-0.09-0.240.310.050.020.070.12-0.12-0.300.240.250.10-0.330.12-0.08-0.15-0.07-0.020.19-0.13表9 白葡萄主成分分析法的前10个特征根对就的特征向量由此可得10个主成分分别为. . . . . . . . . . . . . . . . . . . 分别以10个主成分的贡献率为权重，构建主成分综合评价模型：把

42、白葡萄酒的10个主成分值代入上式，可以得到各白葡萄的综合评价值与其排序结果，见下表10所示。白葡萄样品27样品24样品3样品5样品28样品9样品26名次1.0002.0003.0004.0005.0006.0007.000综合评价值1.8081.2930.7140.6590.6250.3450.334白葡萄样品20样品10样品25样品7样品15样品6样品23名次8.0009.00010.00011.00012.00013.00014.000综合评价值0.3130.2980.2700.1730.1450.0960.052白葡萄样品12样品18样品2样品4样品21样品22样品13名次15.00016.00017.00018.00019.00020.00021.000综合评价值