数学建模论文---葡萄酒的评价

数学建模论文---葡萄酒的评价
承诺书
我们仔细阅读了中国大学生数学建模竞赛的竞赛规则.
我们完全明白，在竞赛开始后参赛队员不能以任何方式（包括电话、电子邮件、网上咨询等）与队外的任何人（包括指导教师）研究、讨论与赛题有关的问题。

我们知道，抄袭别人的成果是违反竞赛规则的, 如果引用别人的成果或其他公开的资料（包括网上查到的资料），必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。

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

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

我们授权全国大学生数学建模竞赛组委会，可将我们的论文以任何形式进行公开展示（包括进行网上公示，在书籍、期刊和其他媒体进行正式或非正式发表等）。

我们参赛选择的题号是（从A/B/C/D中选择一项填写）： A
我们的参赛报名号为（如果赛区设置报名号的话）：4198
所属学校（请填写完整的全名）：广东医学院（东莞校区）
参赛队员(打印并签名) ：1. 黄洁
2. 顾家荣
3. 陈婉君
指导教师或指导教师组负责人(打印并签名)：唐国平
日期：2013年9月 9日
编号专用页
赛区评阅编号（由赛区组委会评阅前进行编号）：
全国统一编号（由赛区组委会送交全国前编号）：全国评阅编号（由全国组委会评阅前进行编号）：
葡萄酒质量的评价模型
摘要
本文主要讨论了关于葡萄酒与葡萄之间关系的研究，主要分析了附件1中两组
评酒员的评价结果有无显著性差异，并判断哪一组结果更可信；还根据酿酒葡萄的
理化指标和葡萄酒的质量把这些酿酒葡萄分为3个等级；分析了酿酒葡萄与葡萄酒
的理化指标之间的联系和酿酒葡萄和葡萄酒的理化指标对葡萄酒质量的影响，并论
证用葡萄和葡萄酒的理化指标来评价葡萄酒的质量。

通过这些分析有益于对葡萄酒
行业的发展有一定的贡献。

对于问题一，用10个品酒员对每种酒样品的总评分的来代表这种酒样品的质量，建立单因子数学模型，分别对两个水平进方差分析，由U检验，取置信区间为95%，最终得出两组品酒员对红葡萄酒的评分有显著性差异，对白葡萄酒的评分没有显著性差异。

并通过对两组品酒员对两种葡萄酒评分的方差分析可得知第二组品酒员的结果更可信。

对于问题二，，我们考虑将众多指标数据经过转换，统一成与感官排序一样的排序类型数据，并把转换后的指标即可直接用来对葡萄进行分级。

对所考虑的众多变量用数学统计方法，经过正交化处理，变成一些相互独立、为数较少的综合指标（即主导因子)后根据酿酒葡萄的理化指标和葡萄酒的质量用聚类分析- K-均值法（快速聚类法）对酿酒葡萄分为三类。

对于问题三，我们对附件二的数据用EXCEL处理并查阅资料可以得出红葡萄酒的理化指标主要影响因素有花色苷、单宁、总酚、酒总黄酮，白葡萄酒的理化指标的主要影响因素有单宁、总酚、酒总黄酮。

在此基础上利用matlab软件分析酿酒葡萄与葡萄酒的理化指标之间的偏相关系数从而得出酿酒葡萄与葡萄酒的理化指标之间的联系。

对于问题四，我们利用SPSS双变量相关性分析法分析酿酒葡萄和葡萄酒的理化指标对葡萄酒质量的相关关系，可知酿酒葡萄和葡萄酒的理化指标与葡萄酒质量相关显著，从而得出能用葡萄和萄酒的理化指标来评价葡萄酒的质量的结论。

关键词葡萄酒单因子方差分析 U检验法聚类分析- K-均值法（快速聚类法） SPSS双变量相关分析线性回归分析
一、问题重述
葡萄酒是用新鲜的葡萄或葡萄汁经发酵酿成的酒精饮料。

通常分红葡萄酒和白
葡萄酒两种。

前者是红葡萄带皮浸渍发酵而成；后者是葡萄汁发酵而成的。

医学研
究表明：葡萄的营养很高，而以葡萄为原料的葡萄酒也蕴藏了多种氨基酸、矿物质
和维生素，这些物质都是人体必须补充和吸收的营养品。

而确定葡萄酒质量时一般
是通过聘请一批有资质的评酒员进行品评。

每个评酒员在对葡萄酒进行品尝后对其
分类指标打分，然后求和得到其总分，从而确定葡萄酒的质量。

酿酒葡萄的好坏与
所酿葡萄酒的质量有直接的关系，葡萄酒和酿酒葡萄检测的理化指标会在一定程度
上反映葡萄酒和葡萄的质量。

附件1给出了某一年份一些葡萄酒的评价结果，附件
2和附件3分别给出了该年份这些葡萄酒的和酿酒葡萄的成分数据。

请尝试建立数学模型讨论下列问题：
1. 分析附件1中两组评酒员的评价结果有无显著性差异，哪一组结果更可信？
2. 根据酿酒葡萄的理化指标和葡萄酒的质量对这些酿酒葡萄进行分级。

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

4．分析酿酒葡萄和葡萄酒的理化指标对葡萄酒质量的影响，并论证能否用葡萄和萄酒的理化指标来评价葡萄酒的质量？
二、问题的分析
针对问题一，要讨论的是一、二两组品酒员对葡萄酒的评分结果有无显著性差
异，我们可以对附件一给出的数据进行分析处理，我们对每个品酒员分别对每个样
品酒的指标评分，对红、白酒单独处理，且处理方法相同从而得出可靠结论。

针对问题二，，我们考虑将众多指标数据经过转换，统一成与感官排序一样的排序类型数据，并把转换后的指标即可直接用来对葡萄进行分级。

对所考虑的众多变量用数学统计方法，经过正交化处理，变成一些相互独立、为数较少的综合指标（即主导因子)后根据酿酒葡萄的理化指标和葡萄酒的质量用聚类分析- K-均值法（快速聚类法）对酿酒葡萄分为三级。

针对问题三，我们对附件二的数据进行处理并经过查阅资料可以得出红葡萄酒的理化指标主要影响因素有花色苷、单宁、总酚、酒总黄酮，白葡萄酒的理化指标的主要影响因素有单宁、总酚、酒总黄酮。

在此基础上利用matlab软件分析酿酒葡萄与葡萄酒的理化指标之间的偏相关系数从而得出酿酒葡萄与葡萄酒的理化指
标之间的联系。

针对问题四，我们利用SPSS双变量相关性分析法分析酿酒葡萄和葡萄酒的理化指标对葡萄酒质量的相关关系，可知酿酒葡萄和葡萄酒的理化指标与葡萄酒质量相关显著，从而得出能用葡萄和萄酒的理化指标来评价葡萄酒的质量的结论。

三、基本假设
由于葡萄酒的主要质量指标大体可分为感官指标和理化指标两大类。

感官指标主要指色泽、香气、滋味和典型性方面的要求，理化指标主要指酒精含量（酒精度）、酸度和糖分指标。

根据各个指标的不同影响，应要求，对葡萄酒的质量进行评价，为了更好的对葡萄酒的质量做出评价，减少部分不必要的影响，可作出以下假设：1、假设酿酒葡萄在利用前，果肉和果汁等主要成分及各种有效成分都极小损失（即忽略不计）
2、假设在酿酒工具不影响酒的质量
3、酿酒方式以及酿酒过程对酒的品质无明显影响
4、不同种类酿酒葡萄的成分数据值统一标准无差异
5、不同种类葡萄酒的成分数据值统一标准无差异
6、品酒的先后打分没有影响
7、检测理化指标为标准值无差异、
8、假设两个水平下的数据都大致服从正态分布。

9、假设各品酒员的评分数据都相互独立。

10、假设酿造葡萄酒的环境是相同的。

11、假设不考虑多种葡萄可制成一种酒，只考虑一种葡萄制成一种酒；
12、假设本文所引用的数据、资料均真实可靠。

四、符号说明
问题一：
1234212223第一组品酒员对红葡萄酒的评分平均分；第二组品酒员对红葡萄酒的评分平均分；第一组品酒员对白葡萄酒的评分平均分；第二组品酒员对白葡萄酒的评分平均分；第一组品酒员对红葡萄酒的评分方差；
第二组品酒员对红葡萄酒的评分方差；第一组品酒------------------------------------------------------------------------m m m m v v v 2412员对白葡萄酒的评分方差；第二组品酒员对白葡萄酒的评分方差；
（、）表示两组品酒员评分呈正态分布的值；---------------------=i v u
问题三：
123456^
1^2^3^412345n n n n n n h h h h DPPH x x x x x ---------------------------------------------------------------------------白葡萄单宁白葡萄总酚白葡萄酒总黄酮白葡萄蛋白质白葡萄果梗比
白葡萄DPPH 半抑制体积白酒单宁白酒总酚白酒酒总黄酮白酒半抑制体积红葡萄花色苷红葡萄总酚红葡萄单宁红葡萄苹果酸红葡萄67^1^2^3^4x x y y y y ------------------------------总糖红葡萄多酚红葡萄酒总黄酮红酒花色苷红酒单宁红酒总酚红酒酒总黄酮
五、模型的建立与求解
（一）对于问题一
1、对于红葡萄酒或白葡萄酒，10个品酒员对每种酒样品的总平均分就是这个组的最终评分。

2 、求解如下：采用将这两组数据看作一个分别看作是两个单因子变量实验，水平有两个分别是第一组、第二组，以下分别是两组品酒员对红葡酒和白葡萄酒评分的平均分折线图。

从附件一的数据我们利用MATLAB软件编程
红酒
m =[62.7 80.3 80.4 68.6 73.3 72.2 71.5 72.3 71.8 74.2 70.1 53.9 74.6 73 58.7 74.9 79.3 59.9 78.6 79.4 77.1 77.2 85.6 78 69.2 73.8 73]';
y=[68.1 73.7 74.6 71.2 72.1 66.3 65.3 66 78.2 68.8 61.6 68.3 68.8 72.6 65.7 69.9 73.6 65.4 72.6 75.8 72.2 71.6 77.1 71.5 68.2 72 71.5]';
z=[x,y]; m=mean(z) v=var(z) cz=cov(z) cv=diag(cz)
corz=corrcoef(z) ex1307 平均数
1m = 72.7259 2m =70.4704 方差
21v = 51.4851 2
2v =15.5006
白酒
x=[82 74.2 78.3 79.4 71 68.4 77.5 71.4 72.9 74.3 72.3 63.3 65.9 72 72.4 74 78.8 73.1 72.2 77.8 76.4 71 75.9 73.3 77.1 81.3 64.8 81.3]';
y=[77.9 75.8 75.6 76.9 81.5 75.5 74.2 72.3 80.4 79.8 71.4 72.4 73.9 77.1 78.4 67.3 80.3 76.7 76.4 76.6 79.2 77.4 76.1 79.5 74.3 77 79.6]'; ex1307 平均数
3m = 74.0107 4m =76.5321
方差
23v = 23.0788 2
4v = 10.0549
我们利用U 检验法取置信水平为95%利用公式可计算得出对于红酒它的U 值为1.43对数据进行处理以及检验
012,112012:,1.432.3181
由题意，可设原假设H 备选假设：在下，同理得=≠=
==m m H m m H u u
111212
2
2
1.96,,u
u u u u α
αα-
--=<>明显可以看出而所以对于红酒两组的评分显著性差异
而对于白酒两个组的评分没有有显著性差异。

同时通过比较两组品酒员分别对两种酒的评分方差可知都是第二组的可信度高。

（二）对于问题二的分析、模型建立与求解
根据酿酒葡萄的理化指标和葡萄酒的质量对酿酒葡萄进行分级——聚类分析- K-均值法（快速聚类法）
根据题意要将葡萄的理化指标与葡萄酒的质量统一结合作为参考。

而葡萄酒质
量则是通过问题一中感官评价的得分反应的。

由于理化指标过多，因此在解决本问时，首先应该完成对指标的处理，尤其是怎样将附表三的芳香物质与附表二中的理化指标结合起来。

由于指标的繁杂，且难以确定指标是偏大型还是偏小型，因此，可考虑将众多指标数据经过转换，统一成与感官排序一样的排序类型数据，这样，转换后的指标即可直接用来对葡萄进行分级。

题目要求我们根据酿酒葡萄的理化指标和葡萄酒的质量对酿酒葡萄进行分级。

经验告诉我们，葡萄的理化指标越合理、葡萄酒的质量越好，该酿酒葡萄的质量也就越好。

这就要求我们分析葡萄的具体理化对葡萄的综合得分的贡献，并结合所酿葡萄酒的得分去评价葡萄的等级。

在葡萄品质的评价过程中，如果将葡萄所具备的每个理化指标不分主次进行评判不仅会增加工作量，也极有可能对评判结果产生比较大的影响。

因此，必须对所考虑的众多变量用数学统计方法，经过正交化处理，变成一些相互独立、为数较少的综合指标（即主导因子)。

2、对此提，我们根据酿酒葡萄的理化指标和葡萄酒的质量对酿酒葡萄进行分级——聚类分析- K-均值法（快速聚类法）
以下是利用酿酒红葡萄的理化指标对酿酒葡萄酒以快速聚类法进行分级：
表2.1
输出结果中，表2.2表示的是初始聚类的中心，也就是种子点
表2.2
表2.3表示的是迭代历史记录
表2.3
表2.4
表2.4是聚类表，表示的是每个个案的分类情况：第3列“聚类”表示的是该案例属于哪一类，第4列“距离”表示该案例与其所属类别重心之间的距离。

分析上表可知，若采用“K-均值聚类法（快速聚类法）”
第1类包括样品1、5、8、10、13、14、16、17、24、25、26、27
第2类包括样品2、4、6、7、9、11、12、15、18、19、20、20、22、23
第3类包括样品3 21。

表2.5
表2.5表示的是最终聚类中心，可以看出，第3类的平均指标最高，第2类居中，第3类最低
表2.6
表2.6表示的是
(1)最终聚类中心间的距离，可以看出，第1类与第2类之间的距离要比第1类与第3类之间的距离小。

(2)每个聚类中的案例数，可以看出，第1类有12个样本，第2类中有13个样本，第3类有2个样品
（2）对酿酒白葡萄进行分类利用酿酒白葡萄指标对白葡萄进行分级，：
在SPSS中运行程序（相关程序见附件一），得到如下图表
聚类成员如下：
第1类包括样品：3、7、8、11、12、13、16、26
第2类包括样品：1、2、4、5、6、9、10、14、15、17、18、19、20、21、22、23、
（三）对于问题三的模型建立、分析与求解
对附件二的数据进行主成分分析可以得出红葡萄酒的理化指标的主要影响因素有花色苷、单宁、总酚、酒总黄酮，白葡萄酒的理化指标的主要影响因素有单宁、总酚、酒总黄酮。

在此基础上利用matlab软件分析酿酒红葡萄与红葡萄酒的理化指标之间的偏相关系数如下
对于白酒单宁
n1=[2.947,2.239,2.990,3.148,2.626,4.502,4.729,1.672,4.434,6.781,3.312,3 .212,2.129,2.388,2.751,2.228,2.247,5.783,2.217,3.141,1.952,6.463,3.389, 8.506,2.757,5.517,6.251,4.583];
n3=[2.557,1.334,3.809,1.548,2.379,6.360,6.824,2.838,2.785,3.487,4.404,4 .314,7.044,2.393,5.190,3.135,2.323,6.261,2.152,2.341,0.744,2.710,2.508, 7.923,2.625,2.816,9.526,3.885];
h=[1.620,1.233,2.009,2.017,1.595,1.289,1.374,1.513,1.844,2.058,1.415,2. 307,1.515,1.320,2.530,1.279,1.549,1.330,1.963,2.676,1.204,1.897,1.330,4 .473,1.505,1.569,3.375,2.029];
save ex1330_data n1 n3 h
load ex1330_data
n=[ones(28,1),n1',n3'];
[b,bint,r,rint,stats]=regress(h',n);
b
bint
stats
ex1940
相关参数b =
0.8555
0.1783
0.0863
13223431234126150.85550.17830.08635.29580.37220.05570.0085;0.71700.10770.01070.1093;
0.02930.00350.00310.0600.0783;h h h n n n n n n n n DPP n n h H n n ∧
∧
∧
∧
=++=-+-+=+-+=-++++即同理得白酒总黄酚白酒回归方程为酿酒葡萄与葡萄酒单宁指标总酚白之间酒半的联系。

抑制体积红123572123411237;4123 3.65760.00670.24820.09950.02830.0716;28.5476 1.9550 3.8159 6.075211.5306;-0.5219+0.00590.2378x 0.0070.1456=1.8507+0.0058+0.0391+0x x x x x x x x x x x x y x y y y ∧
∧∧∧
=-++++-=-+-++=+++酒单宁红酒花色苷红酒总黄酚红酒总酚37.0801-0.1044;
从而得出酿酒红葡萄与红葡萄酒各指标间的联系：
①酿酒红葡萄花色苷、总酚、单宁、总糖与红葡萄酒的单宁呈正相关性，而酒总黄酚与红葡萄酒的单宁呈负相关性而且相关性显著；
②酿酒红葡萄的花色苷、总酚、单宁、苹果酸与红葡萄酒的花色苷的相关性极强； ③酿酒红葡萄的花色苷、总酚、单宁、酒总黄酮与红葡萄酒的酒总黄酮程正相关性，相关性显著；
④酿酒红葡萄的花色苷、总酚、单宁、酒总黄酮与红葡萄酒的总酚相关性显著；
酿酒白葡萄与白葡萄酒各指标间的联系：
①酿酒白葡萄的单宁、酒总黄酮与白葡萄酒的单宁呈正相关性,且相关性较显著；
②酿酒白葡萄的总酚、蛋白质与白葡萄酒的酒总黄酮呈正相关性，酿酒白葡萄酒总黄酮与白葡萄酒的酒总黄酮呈负相关性，但相关性都较显著；
③酿酒白葡萄的单宁、酒总黄酮与白葡萄酒的总酚呈正相关性，酿酒白葡萄总酚与白葡萄酒的总酚呈负相关性，但相关性都较显著；
④酿酒白葡萄的总酚、单宁与白葡萄酒的DPPH 半抑制体积呈正相关性，且相关性显著。

（四）对于问题四的模型建立、分析与求解 利用SPSS 对附件二、三的数据进行分析结果如下：
表4.1：SPSS数据分析结果
表4.2
再利用SPSS双变量相关性分析法分析酿酒葡萄和葡萄酒的理化指标对葡萄酒质量的相关关系如下：
由表可看出酿酒葡萄和葡萄酒的理化指标对葡萄酒质量相关显著，故可用葡和葡萄酒的理化指标来评价葡萄酒的质量。

六、结果分析、验证与模型检验及修正
本文主要应用数理统计、多元回归分析、相关分析、聚类分析等知识。

根据数理统计知识，我们利用SPSS进行数据处理研究，判断出两组评酒员的评酒结果有显著性差异，并选择出数据较为可靠的一组。

根据聚类分析对数据处理研究，基本得到需要的答案。

但依旧存在很多不足之处：
（1）假设评酒员对每种葡萄酒的评价结果大致符合正太分布；
（2）假设酿造葡萄酒的环境是相同的；
（3）只考虑低部分芳香物质影响酒的香气，忽略了其他成分的影响。

七、模型的推广
本文利用聚类分析对给定的数据进行处理，在SPSS中实现等级的划分。

该模型用于生活实践中，也可以解决很多实际问题，例如医学实践中根据各种化验结果、
疾病症状、体征判断患者患的是什么病；体育选材中根据运动员的体形、运动成绩、生理指标心理素质指标、遗传因素判断是否选入运动队继续培养，等等。

它在生活中有广泛的适用性。

八、模型的评价与改进
(一)模型优缺点分析
1、模型的优点：
简单快速地把葡萄酒进行分级，还有利于厂家再生产葡萄酒过程中要注意的问题。

2、模型的缺点:
在建模过程中，由于葡萄酒在生产过程中及装酒的各种要求，要注意一下问题：
（1）葡萄酒在包装运输过程要注意的是：瓶装酒须装入绿色、棕色或无色玻璃瓶中，瓶口封闭严密，不得有漏气、漏酒现象。

（2）酒瓶外部要贴有整齐干净的标签。

（3）输酒设备、灌装设备、瓶塞等应彻底清洗消毒
空瓶，尤其是回收的空瓶，应先经挑拣，剔除有异臭和不易洗净的空瓶，并进行严格的清洗和消毒。

九、参考文献
[1] 孙祝岭徐晓岭，数理统计，北京：高等教育出版社，2009；
[2] 茆诗松王静龙，数理统计，上海：华东师范大学出版社，1990；
[3] 王正林龚纯何倩，精通MATLAB科学计算（第2版）北京：电子工业出版社，2009.8；
[4] 卢纹岱，统计分析，北京：电子工业出版社，2000；
[5] 王颉，试验设计与SPSS应用，北京：化学工业出版社，2007；
[6] 袁志发，周静芋，多元统计分析，北京：科学出版社，2002。

十、附录
附件一：在SPSS中运行聚类分析的程序，对酿酒白葡萄进行分类
聚类成员：
附件二：
白葡萄酒酒总黄酚：
n2=[5.336,5.090,6.972,5.248,6.323,10.541,10.267,5.134,5.814,7.728,7.854 ,8.483,11.774,5.324,8.871,5.007,6.575,11.957,4.725,5.251,4.367,6.409,5. 127,10.755,7.666,5.816,16.965,6.567];
n3=[2.557,1.334,3.809,1.548,2.379,6.360,6.824,2.838,2.785,3.487,4.404,4
.314,7.044,2.393,5.190,3.135,2.323,6.261,2.152,2.341,0.744,2.710,2.508, 7.923,2.625,2.816,9.526,3.885];
n4=[496.457,538.451,467.239,496.201,467.203,499.080,560.343,493.739,482 .855,515.756,557.441,457.649,459.397,524.857,546.446
464.819,416.876,581.913,455.675,479.934,585.359,402.156,505.339,629.801 ,516.862,496.835,642.373,450.458];
h=[0.105,0.510,3.669,1.132,1.414,0.079,3.931,0.577,0.100,1.563,2.257,1. 492,2.036,2.544,0.942,1.923,0.500,2.878,0.408,0.901,0.541,0.089,0.100,3 .305,2.334,0.865,7.655,0.423];
save ex1330_data n2 n3 n4 h
load ex1330_data
n=[ones(28,1),n2',n3',n4'];
[b,bint,r,rint,stats]=regress(h',n);
b
bint
stats
ex1940
b =
-5.2958
0.3722
-0.0557
0.0085
bint =
-9.3458 -1.2458
-0.0763 0.8208
-0.6540 0.5425
-0.0002 0.0173
stats =
0.6183 12.9567 0.0000 1.1690
白葡萄酒总酚：
n1=[2.947,2.239,2.990,3.148,2.626,4.502,4.729,1.672,4.434,6.781,3.312,3 .212,2.129,2.388,2.751,2.228,2.247,5.783,2.217,3.141,1.952,6.463,3.389, 8.506,2.757,5.517,6.251,4.583];
n2=[5.336,5.090,6.972,5.248,6.323,10.541,10.267,5.134,5.814,7.728,7.854 ,8.483,11.774,5.324,8.871,5.007,6.575,11.957,4.725,5.251,4.367,6.409,5. 127,10.755,7.666,5.816,16.965,6.567];
n3=[2.557,1.334,3.809,1.548,2.379,6.360,6.824,2.838,2.785,3.487,4.404,4 .314,7.044,2.393,5.190,3.135,2.323,6.261,2.152,2.341,0.744,2.710,2.508, 7.923,2.625,2.816,9.526,3.885];
h=[1.264,1.104,1.820,1.485,1.537,1.176,1.202,0.472,1.287,1.325,1.276,1. 998,1.356,1.320,1.807,1.307, 1.269, 1.343, 1.343, 1.315, 1.029, 1.380, 1.114, 3.434, 1.459, 1.258, 2.539, 1.544];
save ex1330_data n1 n2 n3 h
load ex1330_data
n=[ones(28,1),n1',n2',n3'];
[b,bint,r,rint,stats]=regress(h',n);
b
bint
stats
ex1940
b =
0.7170
0.1077
-0.0107
0.1093
bint =
0.1098 1.3243
-0.0062 0.2217
-0.1823 0.1609
-0.1284 0.3471
stats =
0.4358 6.1791 0.0029 0.1757
白葡萄酒DPPH半抑制体积：
n1=[2.947,2.239,2.990,3.148,2.626,4.502,4.729,1.672,4.434,6.781,3.312,3 .212,2.129,2.388,2.751,2.228,2.247,5.783,2.217,3.141,1.952,6.463,3.389, 8.506,2.757,5.517,6.251,4.583];
n2=[5.336,5.090,6.972,5.248,6.323,10.541,10.267,5.134,5.814,7.728,7.854 ,8.483,11.774,5.324,8.871,5.007,6.575,11.957,4.725,5.251,4.367,6.409,5. 127,10.755,7.666,5.816,16.965,6.567];
n5=[4.14,3.96,3.74,3.64,4.15,4.99,2.86,3.75,3.85,3.66,3.26,2.41,3.40,3. 50,3.59,5.15,6.27,3.35,3.06,5.19,3.66,4.74,5.39,2.64,2.92,3.85,2.66,3.5 7];
n6=[0.3267,0.3061,0.2837,0.3685,0.2462,0.3369,0.3709,0.1045,0.3787,0.34 21,0.2320,0.4490,0.4094,0.4177,0.3540,0.0976,0.2599,0.2113,0.2004,0.335 0,0.2526,0.3660,0.3600,0.4066,0.2125,0.2998,0.3850,0.3654];
h=[0.0348,0.0331,0.0474,0.0526,0.0406,0.0420,0.0522,0.0392,0.0400,0.064 0,0.0243,0.0817,0.0470,0.0491,0.0738,0.0315,0.1321,0.0385,0.0373,0.0544 ,0.0464,0.0498,0.0382,0.1434,0.0306,0.0441,0.1031,0.0541];
save ex1330_data n1 n2 n5 n6 h
load ex1330_data
n=[ones(28,1),n1',n2',n5',n6'];
[b,bint,r,rint,stats]=regress(h',n);
b
bint
stats
ex1940
b =
-0.0293
0.0035
0.0031
0.0060
0.0783
bint =
-0.1045 0.0459
-0.0036 0.0107
-0.0013 0.0076
-0.0065 0.0186
-0.0491 0.2058
stats =
0.3092 2.5737 0.0647 0.0007
红葡萄酒单宁：
x1=[408.028 224.367 157.939 79.685 120.606 46.186 60.767 241.397 240.843 44.203 7.787 32.343 65.324 140.257 52.792 60.660 59.424 40.228 115.704 23.523 89.282 74.027 172.626 144.881 49.643 58.469 34.190];
x2=[23.604 26.875 21.685 10.698 17.618 10.671 9.214 15.241 30.114 9.476 6.075 12.059 14.385 14.657 11.901 11.214 15.336 7.381 17.426 12.677 16.192 16.442 29.704 8.751 11.502 7.348 8.897];
x3=[22.019,23.361,20.373,8.638,14.486,15.173,5.619,22.489,24.362,16.688 ,4.543,7.169,9.822,13.941,25.417,10.086,15.730,5.388,13.700,8.115,13.61 3,12.155,24.257,14.417,9.324,3.778,10.310];
x5=[208.175,205.000,256.190,189.722,209.663,244.385,209.861,198.849,193 .690,167.202,209.563,247.659,197.857,191.508,179.107,204.008,212.738,22 6.032,205.794,193.194,205.794,224.147,207.679,201.825,150.337,173.353,1 96.667];
x7=[9.480 13.806 10.794 4.482 10.275 6.838 3.468 8.483 20.490 4.631 2.517 3.897 7.330 7.809 5.511 9.157 8.701 5.245 9.454 8.155 7.515 7.846 24.295 8.206 5.373 3.383 4.711];
y=[11.030 11.078 13.259 6.477 5.849 7.354 4.014 12.028 12.933 5.567 4.588 6.458 6.385 6.073 3.985 4.832 9.170 4.447 5.981 5.864 10.090 7.105 10.888 5.747 5.406 3.615 5.961];
save ex1328_data x1 x2 x3 x5 x7 y
load ex1328_data
x=[ones(27,1),x1',x2',x3',x5',x7'];
[b,bint,r,rint,stats]=regress(y',x);
b
bint
stats
ex1329
b =
-3.6576
0.0067
0.2482
0.0995
0.0283
-0.0716
bint =
-9.5259 2.2106
-0.0048 0.0182
-0.0341 0.5305
-0.0582 0.2571
-0.0004 0.0569
-0.3825 0.2394
stats =
0.7645 13.6332 0.0000 2.4598
红葡萄酒总酚：
x1=[408.028 224.367 157.939 79.685 120.606 46.186 60.767 241.397 240.843 44.203 7.787 32.343 65.324 140.257 52.792 60.660 59.424 40.228 115.704 23.523 89.282 74.027 172.626 144.881 49.643 58.469 34.190];
x2=[23.604 26.875 21.685 10.698 17.618 10.671 9.214 15.241 30.114 9.476 6.075 12.059 14.385 14.657 11.901 11.214 15.336 7.381 17.426 12.677 16.192 16.442 29.704 8.751 11.502 7.348 8.897];
x3=[22.019,23.361,20.373,8.638,14.486,15.173,5.619,22.489,24.362,16.688
,4.543,7.169,9.822,13.941,25.417,10.086,15.730,5.388,13.700,8.115,13.61 3,12.155,24.257,14.417,9.324,3.778,10.310];
x7=[9.480 13.806 10.794 4.482 10.275 6.838 3.468 8.483 20.490 4.631 2.517 3.897 7.330 7.809 5.511 9.157 8.701 5.245 9.454 8.155 7.515 7.846 24.295 8.206 5.373 3.383 4.711];
y=[11.030 11.078 13.259 6.477 5.849 7.354 4.014 12.028 12.933 5.567 4.588 6.458 6.385 6.073 3.985 4.832 9.170 4.447 5.981 5.864 10.090 7.105 10.888 5.747 5.406 3.615 5.961];
save ex1330_data x1 x2 x3 x7 y
load ex1330_data
x=[ones(27,1),x1',x2',x3',x7'];
[b,bint,r,rint,stats]=regress(y',x);
b
bint
stats
ex1329
b =
1.8507
0.0058
0.3091
0.0801
-0.1044
bint =
-0.0973 3.7987
-0.0065 0.0181
0.0150 0.6033
-0.0869 0.2471
-0.4345 0.2257
stats =
0.7172 13.9467 0.0000 2.8196
红葡萄酒酒总黄酮：
x1=[408.028 224.367 157.939 79.685 120.606 46.186 60.767 241.397 240.843 44.203 7.787 32.343 65.324 140.257 52.792 60.660 59.424 40.228 115.704 23.523 89.282 74.027 172.626 144.881 49.643 58.469 34.190];
x2=[23.604 26.875 21.685 10.698 17.618 10.671 9.214 15.241 30.114 9.476 6.075 12.059 14.385 14.657 11.901 11.214 15.336 7.381 17.426 12.677 16.192 16.442 29.704 8.751 11.502 7.348 8.897];
x3=[22.019,23.361,20.373,8.638,14.486,15.173,5.619,22.489,24.362,16.688 ,4.543,7.169,9.822,13.941,25.417,10.086,15.730,5.388,13.700,8.115,13.61 3,12.155,24.257,14.417,9.324,3.778,10.310];
x7=[9.480 13.806 10.794 4.482 10.275 6.838 3.468 8.483 20.490 4.631 2.517 3.897 7.330 7.809 5.511 9.157 8.701 5.245 9.454 8.155 7.515 7.846 24.295 8.206 5.373 3.383 4.711];
y=[8.020 13.300 7.368 4.306 3.644 4.445 2.765 7.748 9.905 3.145 2.103 2.986 3.957 3.068 1.836 2.668 4.912 3.531 3.875 4.044 4.440 5.827 12.144 3.731 3.022 2.154 3.284];
save ex1330_data x1 x2 x3 x7 y
load ex1330_data
x=[ones(27,1),x1',x2',x3',x7'];
[b,bint,r,rint,stats]=regress(y',x);
b
bint
stats
ex1329
b =
-0.5219
0.0059
0.2378
0.0071
0.1456
bint =
-2.2063 1.1624
-0.0047 0.0166
-0.0165 0.4922
-0.1373 0.1515
-0.1398 0.4310
stats =
0.7998 21.9748 0.0000 2.1081
红葡萄酒花色苷：
x1=[408.028 224.367 157.939 79.685 120.606 46.186 60.767 241.397 240.843 44.203 7.787 32.343 65.324 140.257 52.792 60.660 59.424 40.228 115.704 23.523 89.282 74.027 172.626 144.881 49.643 58.469 34.190];
x2=[23.604 26.875 21.685 10.698 17.618 10.671 9.214 15.241 30.114 9.476 6.075 12.059 14.385 14.657 11.901 11.214 15.336 7.381 17.426 12.677 16.192 16.442 29.704 8.751 11.502 7.348 8.897];
x3=[22.019,23.361,20.373,8.638,14.486,15.173,5.619,22.489,24.362,16.688 ,4.543,7.169,9.822,13.941,25.417,10.086,15.730,5.388,13.700,8.115,13.61 3,12.155,24.257,14.417,9.324,3.778,10.310];
x4=[18.210 4.750 2.960 5.230 3.770 2.210 7.740 13.550 4.120 2.300 8.610 5.330
0.830 4.120 3.630 7.280 5.110 5.590 4.270 0.920 2.930 7.730 5.200 4.600 2.480
1.400 1.390];
y=[973.878 517.581 398.770 183.519 280.190 117.026 90.825 918.688 387.765 138.714 11.838 84.079 200.080 251.570 122.592 171.502 234.420 71.902 198.614 74.377 313.784 251.017 413.940 270.108 158.569 151.481 138.455];
save ex1328_data x1 x2 x3 x4 y
load ex1328_data
x=[ones(27,1),x1',x2',x3',x4'];
[b,bint,r,rint,stats]=regress(y',x);
b
bint
stats
ex1329
b =
-28.5476
1.9550
-3.8159
6.0752
11.5306
bint =
-147.8532 90.7581
1.0831
2.8269
-13.5231 5.8912
-2.4231 14.5735
-2.2468 25.3079
stats =
1.0e+003 *
0.0009 0.0415 0.0000 7.3186
13。