葡萄酒质量评价模型

1、葡萄酒质量评价模型摘 要葡萄酒质量的高低评估是通过评酒专家对葡萄酒的感官评分来体现。酿酒葡萄和葡萄酒的理化指标一定程度上反映了葡萄酒的质量。问题一，首先对附件1的数据进行预处理，分别求得评酒员关于样品酒的4组平均得分，在此基础上，利用f检验，发现不管对于红葡萄酒还是白葡萄酒，两组评酒专家的评分结果都存在显著的差异。此外，建立了评价可信度的层次分析模型，发现第二组评酒员的评分更加可信。问题二，运用主成分分析对酿酒葡萄的30个理化指标进行降维，主成分降维后减少了变量间的重叠部分，然后通过q型聚类对酿酒葡萄酒的样品进行归类，利用问题一中第二组评分数据，得到每一类样品的平均得分，通过得分的大小来分等级

2、。问题三，建立了酿酒葡萄与葡萄酒理化指标的典型相关分析模型，得出酿酒葡萄与葡萄酒理化指标之间有着密切的联系。如：红葡萄与红葡萄酒的理化指标的第一典型相关系数，第一典型变量可以解释29.9%红葡萄理化指标组内变差，并解释39%红葡萄酒理化指标的变差；其两者的相关系数相互解释每组内的变差。问题四，对于酿酒葡萄和葡萄酒的理化指标对葡萄酒质量的影响，本文先通过线性回归做初步的分析，然后运用topsis模型进行了进一步的分析，得到葡萄和葡萄酒的理化指标不一定能评价葡萄酒的质量，但有一定的联系。关键词：f检验；主成分分析；q型聚类；样品典型相关分析；topsis模型 1、问题提出葡萄酒是用新鲜的葡萄或者葡

3、萄汁经发酵酿成的酒精饮料。质量评价主要通过外观、香气、口味、典型性体现。所以确定葡萄酒的质量一般通过聘请一批有资深的评酒员对葡萄酒进行品尝后对其分类指标打分，然后求和得到总分，从而确定葡萄酒的质量。酿酒葡萄的好坏与所酿葡萄酒的质量有直接的关系，葡萄酒和酿酒葡萄检测的理化指标会在一定程度上反映葡萄酒和葡萄的质量。附件1给出了某一年份一些葡萄酒的评价结果，附件2和附件3分别给出了该年份这些葡萄酒的和酿酒葡萄的成分数据。建立数学模型讨论下列问题：1. 分析附件1中两组评酒员的评价结果有无显著性差异，哪一组结果更可信？2. 根据酿酒葡萄的理化指标和葡萄酒的质量对这些酿酒葡萄进行分级。3. 分析酿酒葡萄

4、与葡萄酒的理化指标之间的联系。4分析酿酒葡萄和葡萄酒的理化指标对葡萄酒质量的影响，并论证能否用葡萄和葡萄酒的理化指标来评价葡萄酒的质量？2、问题分析问题一：要求找出两组评酒专家的评分结果的差异性，可以选用方差分析当中的f检验体现评酒结果的显著性。若无差异性则都可信，若存在差异性是可以通过引入评价可信度的方法，找出到底哪一组更加的可信。问题二：葡萄酒质量高低与酿酒葡萄的优劣有直接的关系，通过主成分的方法进行变量的降维，后对样品进行q聚类，求出每类的均值后进行评级。问题三：酿酒葡萄跟葡萄酒的理化指标关联密切，个别的理化指标是其重要的成分，此过程通过了样品典型相关分析来分析得出。问题四：酿酒葡萄和葡

5、萄酒的理化指标对葡萄酒质量的影响，可以通过酿酒葡萄和葡萄酒的理化指标与葡萄质量之间的相关性得到，若存在相关性可以通过回归思想实现。3、符号说明：第组的第个样品的颜色葡萄酒总平均分 ：第组颜色葡萄酒总平均分的均值：葡萄颜色，n=1 表示红色，n=2 表示白色：第i个样品的n颜色酿酒葡萄的第j个一级理化指标的含量4、模型基本假设（1）每一种葡萄酒的生产工艺是大致相同。（2）葡萄的质量决定葡萄酒的质量。（3）相应的葡萄酒是由相应的酿酒葡萄酿制得到。（4）评酒员的个数足够多，评酒总平均分数能充分反映葡萄酒的质量。5、模型的建立与求解5.1数据的预处理和初步分析 通过对附件一中数据表的初步的观察，我们发

6、现，表中存在好几处的数据异常的情况，比如：第一组白葡萄酒中的样品三7号评酒员在浓度指标中的评分为77，明显的不正确。本文遇到类似的情况用本行的均值代替原来数据。为了评价一种葡萄酒的好坏，通常的做法是由感官评酒专家根据国际葡萄与葡萄组积（oiv）的评价方法对葡萄酒的澄清度、色度、纯净度等方面进行打分，最后将个方面得到总分的方法对葡萄酒的好坏进行评价。 通过对附件1的分析，对同一个样品分两组评酒专家每组10名，分别对27个红葡萄酒的样品和28个白葡萄酒的样品进行了评价。本模型采取10个评酒专家对同一个样品的总评价得分的平均值作为该样品的最后得分。数据部分如表1所示。表1：样品评分均值样品一组红葡萄

7、酒一组白葡萄酒二组红葡萄酒二组白葡萄酒162.7082.0068.1077.90280.3074.2074.0075.80380.4085.3074.6075.602673.8081.3072.0074.302773.0064.8071.5077.002881.3079.60图1和图2可以直观地反映上述数据：图1：红葡萄酒的平均得分图2：白葡萄酒的平均得分由于上述数据是由专家的感官评价得分获得的，在实践中，由于各种因素的共同影响下，专家组成员间和组间之间存在异质性。造成异质性的主要原因有：评价尺度的差异、评价位置的差异 、评价方向的差异这三个方面。5.2方差分析f检验 对同一个样品由两组不同的

8、评酒专家分开独立地进行评价。为了反映出两组评酒专家评价的结果是否存在较大的差异性，本文利用f检验对两组评酒专家的评价结果作显著性检验。5.2.1模型的建立step1:模型的假设两组评酒专家的评价得分可作为同一因素下的不同水平，对不同样品的评分可作为样本观察值。取原假设为： ； 样本的观察值可以分解为： 其中： 红葡萄酒的情况下：n=1时，白葡萄酒时： n=2时 step2:构造f统计量 是第k组数据的组平均值，是总平均值。考察全体数据对的偏差平方和： 对上述的式子分解可得： 记： 是各组均值对总方差的偏差平方和，反映两组品酒员间的差异是各组内的数据对均值偏差平方和的总和。则表示在同一组品酒员下

9、的随机误差的大小。由分布的可加性知： 对进一步分析可得： 当成立时，该比值服从自由度1，的f分布，即： 为检验，给定显著性水平，记f分布的分为数为 若：则接受，否则拒绝。5.2.2模型的求解本模型f检验的显著性水平取： ，由matlab求解可得： 拒绝对于红葡萄酒两组的评分存在显著性差异。 拒绝对于白葡萄酒两组的评分存在显著性差异。综上所述，对于两种葡萄酒，两组专家的给分都存在着显著性的差异。5.3判断可信度的层次分析模型从上述的模型可以看出第一组品酒员和第二组品酒员对于红葡萄酒和白葡萄酒的评价具有显著性差异，葡萄酒质量评价结果可信程度直接关系到消费者的利益和市场对葡萄酒的科学管理。那么对于上

10、述两组评酒专家的评价结果，我们关心的是到底哪一组的可信度比较高，从而选择接受哪一组的评价的结果，为了能更好地反映出两组数据的可信度，我们引入了判断可信度的层次分析模型，对、分配一个可信度权重从而得到一个排序，可以得到四种情况下的可信度排序。5.3.1模型的建立step1：根据评分的极差矩阵和可行度评估标度确定判断矩阵本模型在上述的数据的基础上选择每一种情况下的极差为：其中：， 为每一种情况下的最大和最小上述的极差能够很好地反映出每一组评酒专家的总评分的差异，这个差异可以再一定的程度上体现出这组评酒专家的可信度。根据每两种情况下的差异矩阵根据可行度评估标度，如表2所示。表2：可行度评估标度评估标

11、度值含义0判断完全无把握0.2判断非常无把握0.4判断比较无把握0.6判断比较有把握0.8判断非常有把握1判断完全有把握0.1,0.3,0.5,0.7,0.9表示相邻评估的中间值根据差异矩阵的每个数据的差异程度的大小反映在可信度评估标度表上，以差异越小相互可信度越大为标准可以得到判断矩阵step2:权重的分配 判断矩阵的最大特征值为： 相应的特征向量为： 所以得到权重向量：step3:一致性检验 计算一次性指标：相应的随机一次性指标可以通过查表获得。 计算一次性比率： 若，那么对一次性检验通过。5.3.2模型的求解通过matlab的计算可以得到极差矩阵为：通过上述的差异根据可行度评估标度得到了

12、判断矩阵为：可以得到：，权重向量为：上述的结果表明：无论是红葡萄酒还是白葡萄酒第二组的分配的权重都比对应的第一组的权重都要高，所以可以得到最后的结论为：第二组的评价结果比第一组的评价结果更为可信。5.4主成分模型附件二给出了和红葡萄27个样品和白葡萄28个样品的30项的理化指标数据，我们发现有部分的指标之间有一定的相关性，为了尽量排除指标间重叠部分对葡萄质量的重复的影响带来的误差，本文采用主成分的方法相对30个变量进行降维，提取出30个变量的主要成分。5.4.1模型的建立step1:数据标准化处理 表示第个样品第个指标的取值其中：分别表示红葡萄和白葡萄标准化后的数据为：其中： 为第理化指标的均

13、值，为第理化指标的方差step2:计算相关系数矩阵 根据相关系数的定义：可以得到每两个指标之间的相关系数：step3:寻找主成分 相关系数矩阵排序后的的特征值为：特征值对应的特征向量为：用新的变量作为主成分取代原来的数据： 得到个主成分矩阵形式为：，累加贡献率为：5.4.2模型的求解 （1）对红葡萄数据的求解通过matlab的求解可以得到红葡萄指标可以降到14个主成分部分数据如下所示： 对红葡萄的主成分解释如表3所示：表3：红葡萄主成分解析表主成分包含成分名称累计贡献率第一主成分花色苷、dpph自由基、总酚、23.23%第二主成分总糖、干物质含量39.70%第三主成分白藜芦醇、可滴定酸52.1

14、5%第四主成分苹果酸61.61%第五主成分固酸比、果穗质量68.28%第六主成分果皮颜色74.08%第七主成分黄酮醇、果皮质量、b*（+黄；-蓝）78.81%第八主成分柠檬酸83.04%第九主成分vc含量、葡萄总黄酮、出汁率86.25%第十主成分单宁、还原糖、可溶性固形物88.71%第十一主成分多酚氧化酶活力、百粒质量91.01%第十二主成分氨基酸总量、褐变度92.73%第十三主成分蛋白质、果梗比、l*94.37%第十四主成分酒石酸、ph值95.61%（2）对白葡萄酒的求解通过matlab的求解可以得到红葡萄指标可以降到14个主成分部分数据如下所示： 对红葡萄的主成分解释如表4所示：表4：白葡

15、萄主成分解析表主成分包含成分名称累计贡献率第一主成分可溶性固形物、干物质含量19.58%第二主成分总酚、葡萄总黄酮36.61%第三主成分可滴定酸、固酸比、l*、b*（+黄；-蓝）48.82%第四主成分氨基酸总量、褐变度55.81%第五主成分酒石酸、苹果酸、黄酮醇62.14%第六主成分白藜芦醇67.46%第七主成分ph值72.49%第八主成分单宁、果皮质量76.90%第九主成分果梗比、果皮颜色80.88%第十主成分出汁率84.35%第十一主成分花色苷、柠檬酸87.48%第十二主成分蛋白质、dpph自由基、还原糖、百粒质量90.16%第十三主成分vc含量、果穗质量92.31%第十四主成分多酚氧化酶

16、活力94.00%第十五主成分总糖95.30%（3）对主成分的解析 本文根据需要分别对红葡萄和白葡萄的30个理化指标经主成分分析后，得到14个主成分和15个主成分，累计贡献率可达到95.61%和95.30%。表中数据表明酿酒葡萄的理化性质对酿酒葡萄的优劣有很大的贡献作用。我们进一步对表3进行分析可以知道：对于红葡萄，第一主成分是花色苷、dpph自由基、总酚；第二主成分是总糖、干物质含量；第三主成分是白藜芦醇、可滴定酸说明这些物质比较具有代表性。对表4进行进一步分析可知：第一主成分是可溶性固形物、干物质含量；第二主成分是总酚、葡萄总黄酮；第三主成分是可滴定酸、固酸比、l*、b*（+黄；-蓝）。说明

17、这些物质依次对白葡萄具有代表性。5.5聚类模型我们的目的是为了对不同的样品进行评级，本模型先通过最远距离法分别对红葡萄和白葡萄的样品进行聚类，体现样品之间的差异后，然后对聚类的样品分在同一个等级。5.5.1q聚类模型的建立每个样品可以看做维空间的一个向量本模型使用欧氏距离定义每两个样品之间的距离每两个样品之间的距离定义为：两个类别之间的距离使用最长距离法：则： 5.5.2模型的求解运用spss软件求解可得聚类图，如图3与图4所示图3：红葡萄的样品聚类图图4：白葡萄的的样品聚类图5.6葡萄质量的等级模型酿酒葡萄的好坏与所酿葡萄酒的质量有直接的关系，从上述的聚类图可以根据葡萄酒样品的聚类结果和葡萄

18、酒的品分均值进行了等级的划分。表5：红葡萄酒样品聚类表类别样品平均得分14、6、7、11、12、15、18、19、20、22、23、369.1721、2、5、8、9、10、13、14、16、17、24、25、26、27、2171.16表6：白葡萄酒样品聚类表类别样品平均得分11、6、7、8、11、12、13、14、16、17、18、19、2274.982375.60321、23、2777.87415、24、2878.0352、4、5、9、10、20、25、2678.10上述的归类可以把红葡萄酒样品分为好、差两类，白葡萄酒归为五类为：很好、好、中等、差、很差。5.7典型相关系数模型5.7.1模型

19、的建立step1：计算相关系数矩阵将相应剖分为其中：分别为n颜色葡萄理化指标变量和n颜色葡萄酒理化指标变量的相关系数阵。为n颜色葡萄理化指标与n颜色葡萄酒变量的相关系数阵。step2:求典型相关系数及典型变量可得：的特征根和特征向量。 的特征根，特征向量，则有则典型变量为step3:典型相关系数的显著性检验。step4:典型结构与典型冗余分析。根据典型结构可以计算任一个典型变量或解释本组变量x (或y)总变差的百分比。同时可求得前t 个典型变量（或）解释本组变量x (或y)总变差的累计百分比典型冗余分析用来研究典型变量解释另一组变量总变差百分比的问题。第二组典型变量解释第一组变量x 总变差的百

20、分比 。5.7.2模型的求解 酿酒葡萄与葡萄酒的理化指标的联系分析 将附录2 经过处理的指标数据利用matlab软件的canoncorr函数进行处理，得出如下结果。典型相关系数及其检验如表 7所示：序号123456789典型相关系数红 0.990.920.890.860.620.490.440.280.11白 0.920.810.660.580.420.260.05表7：典型变量相关系数由上表可知，红葡萄和红葡萄酒理化指标的前两个典型相关系数均较高，表明相应典型变量之间密切相关。白葡萄与白葡萄酒的第一个典型相关系数较高，表明相应典型变量之间密切相关。但要确定典型变量相关性的显著程度，尚要进行相

21、关系数的统计量检验，具体做法是：比较统计量计算值与临界值的大小，据比较结果判定典型变量相关性的显著程度。其结果如表8所示 表8 相关系数检验表序号123456789概率红0.00020.0150.0870.3610.8780.8990.8330.8210.634白0.0410.4310.7930.8880.9590.9670.97从上表看着9、7对典型变量中葡萄与葡萄酒的理化指标的典型变量中红色的第一和第二典型变量、白色的第一典型变量均通过了统计量检验。红葡萄和红葡萄酒的理化指标的相关性很高，第一典型相关系数为0.99.它比红葡萄和红葡萄酒的理化的任一相关系数都高。检验总体中除了第一和第二典型

22、相关系数外。其他的都没有通过检验。第二典型系数为0.92。因此，两组变量的相关性的研究转化为研究第一对和第二对典型相关变量的相关性。白葡萄和白葡萄酒的理化指标的相关性很高，第一典型相关系数为0.92.它比白葡萄和白葡萄酒的理化指标的任一相关系数都要高。检验总体中除了第一典型相关系数外。其他的都没有通过检验。因此，两组变量的相关性的研究转化为研究第一对典型相关变量的相关性。 鉴于原始变量的计量单位不同，不宜直接比较，本文才用标准化的典型系数，给出典型相关模型，如表9所示：表9 ：典型相关模型1，2由表9 第一 组典型相关方程可知，红葡萄酒的理化指标的第一典型变量与呈高度相关，说明在红葡萄酒的理化

23、指标中，占的比重比较大。白葡萄酒的理化指标的第一典型变量与呈高度相关，说明在红葡萄根据第二组典型相关方程，是红葡萄理化指标的主要成分。所以总体上红葡萄理化指标的主要因素按重要程度依次是 ；反应红葡萄酒主要理化指标是典型冗余分析：通过典型变量解释另一组变量总变差百分比的关系，来解释本组变量的信息，还解释另一组变量的信息，典型冗余分析结果如表10所示： 表10 ：典型分析结果123456789mu红0.2990.1530.0760.1680.0770.0410.0630.0830.038白0.1280.2120.0980.1060.0810.1160.134mv红0.2920.1280.0650.

24、1240.0290.0100.0120.0060.001白0.1070.1380.0430.0360.0140.0080.000nu红0.3900.1440.0460.0720.0310.0210.0070.0020.001白0.0980.1370.1130.0230.0230.0040.000nv红0.4010.1710.0580.0980.0820.0880.0370.0280.037白0.1180.2110.2560.0690.1320.0530.069 注：mu：x组原始变量被u_i解释的方差比例，mv：x组原始变量被v_i解释的方差比例，nu：y组原始变量被u_i解释的方差比例，nv

25、：y组原始变量被v_i解释的方差比例） 由表10的典型冗余分析的结果，我们来分析标准化的方差，第一典型变量可以解释29.9%红葡萄理化指标组内变差，并解释39%红葡萄酒理化指标的变差；而典型变量可以解释40.1%红葡萄酒理化指标组内变差，并解释29.2%红葡萄理化指标的变差；第一典型变量可以解释29.2%白葡萄理化指标组内变差，并解释9.8%白葡萄酒理化指标的变差；变量可以解释11.8%白葡萄酒理化指标组内变差，并解释10.7%白葡萄理化指标的变差。5.8酿酒葡萄与葡萄酒的理化指标对葡萄酒质量的影响的分析模型5.8.1多元线性回归模型的初步分析（1）多元回归模型数据的建立 由附录一，评酒专家是

26、通过葡萄酒酒的外观分析、香气分析、口感分析这三个主要的方面进行评分来确定葡萄酒的好坏，通过数据的搜索和总结，经过分析不难发现：附录二中葡萄和葡萄酒的某些理化指标以及附录三中的数据可以提升为一些重要的理化指标。数据的提取如表11所示表11：数据的提取与提升表附录二中的提取附录三种的提升葡萄指标酒石酸、还原糖、ph值、单宁、总酚、干涉出物、色泽（）酒精度葡萄酒指标单宁、总酚、酒总黄酮、白藜芦醇、dpph、色泽（）酒精度酒精度的获得：对于葡萄酒来说酒精量是一个很重要的理化指标，而附录二对该指标的缺失，我们从附录三中用醇类的总量作为酒精度的量化。（2）多元线性回归模型的spss求解结果 利用统计软件s

27、pss本文对上述的葡萄指标和葡萄酒指标分别对专家的评分作线性回归。红葡萄的结果如下所示：上述的回归模型不能通过检验，大部分的变量没有通过检验，但是较大。对于其他的情况下的数据具有类似的结果。我们不能通过上述的模型说明酿酒葡萄和葡萄酒的理化指标对葡萄酒质量有影响，同时也不能说明影响不存在。5.8.2 topsis模型的进一步探讨 对于上述的结果我们没有明确的结论。本文通过topsis模型对该问题作进一步的探讨。（1）数据的获得 不失一般性，本模型联合附录二和附录三中的数据，把附录三中的每一个芳香物质提升为一个理化指标，联合附录二中的数据得到一个数据，并利用上述的主成分模型（累计贡献率为0.8）进

28、行降维，得到个主成分的新的数据。（2） topsis模型的建立step1：用向量规划的方法求得规范决策矩阵。设多属性决策问题的决策矩阵为：a=规范化决策矩阵为：其中 ，。step2：构成加权规范阵 其中：，。各个主成分的权重由主成分的贡献率得到： step3：确定正理想和负理想。设正理想解的第j个属性值为，负理想解第j个属性值为，则 正理想解 负理想解step4：计算各方案到正理想解与负理想解的距离。备选方案到正理想的距离为： ；备选方案到负理想解的距离为：。step5：计算各方案的排队指标值(即综合评价指数）： 。step6：按由大到小排列方案的优劣次序。（3） topsis模型的求解通过m

29、atlab的求解得到原来每个样品的排序和用topsis模型求解结果的排序。表12：样本的排序样品（红）1232627评分排序722241715红葡萄topsis排序2322182113红葡萄酒topsis排序16101527样品（白）1232728评分排序191091624白葡萄topsis排序22145283白葡萄酒topsis排序232218133图5：红葡萄与红葡萄酒topsis排序与评分的排序差图6：白葡萄与红葡萄酒topsis排序与评分的排序差 上述的图5和图6反映出排序的差异有很大的波动性，说明量化指标的变化有很大的不确定性，而评分的好坏是十分稳定的。 上述的两个模型都反映出葡萄和

30、葡萄酒的理化指标不一定能评价葡萄酒的质量，但有一定的联系。6、模型的评价与推广本文运用的各种模型包括主成分分析、q类聚类、典型相关分析、topsis模型等。主成分的维作用可以推广到很多方面。可以应用于影响天气好坏的多方面因素进行降维得到影响天气的几个主要因素。典型相关模型的应用例子：康复俱乐部里成员生理指标与训练指标的相关分析、城市竞争力与基础设施的典型相关分析、家庭特征与消费模式之间的关系。topsis模型模型可以用在样本的排序的各种问题当中，可以解决学生的评优问题。7、参考文献1 韩中庚.数学建模方法及其应用m.北京:高等教育出版社， 2005年6月.2 姜启源，谢金星，叶俊. 数学模型m

31、.北京:高等教育出版社， 2004年.3 蔡锁章. 数学建模原理与方法m. 北京：海洋出版社，2000.4 司守奎，孙玺菁.数学建模算法与应用m.北京:国防工业出版社，2011年8月.5 李华,刘曙东,王华,张予林.葡萄酒感官评价结果的统计分析方法研究j.中国食品学报,2006,6(2):126-130.6 李运,李记明,姜忠军.统计分析在葡萄酒质量评价中的应用j.酿酒科技,2009,178(4):79-82.附录：（1）数据的预处理样品一组红葡萄酒一组白葡萄酒二组红葡萄酒二组白葡萄酒162.7082.0068.1077.90280.3074.2074.0075.80380.4085.3074

32、.6075.60468.6079.4071.2076.90573.3071.0072.1081.50672.2068.4066.3075.50771.5077.5065.3074.20872.3071.4066.0072.30981.5072.9078.2080.401074.2074.3068.8079.801170.1072.3061.6071.401253.9063.3068.3072.401374.6065.9068.8073.901473.0072.0072.6077.101558.7072.4065.7078.401674.9074.0069.9067.301779.3078.80

33、74.5080.301859.9073.1065.4076.701978.6072.2072.6076.402078.6077.8075.8076.602177.1076.4072.2079.202277.2071.0071.6079.402385.6075.9077.1077.402478.0073.3071.5076.102569.2077.1068.2079.502673.8081.3072.0074.302773.0064.8071.5077.002881.3079.60（2）f检验的matlab代码x=;p1=anova1(x);p1（3）主成分降维的matlab代码%输入数据r=;

34、%数据的标准化clcr_var=var(r,1);r_mean=mean(r);a,b=size(r);r_s=ones(a,b);for i=1:a for j=1:b r_s(i,j)=(r(i,j)-r_mean(j)/r_var(j); endend%主成分分析z=r_s;rr=corrcoef(z);tx,tz=eig(rr);tz=diag(tz);so,id=sort(tz,descend);su=sum(so);ss=0;i=1;xl=;while ss0.8 ss=ss+so(i)/su; xl=xl;tx(:,id(i); i=i+1;end%新数据的获得zr=ones(size(r,1),size(xl,1);for i=1:size(r,1) for j=1:size(xl,1) zr(i,j)=sum(r(i,:).*xl(j,:); endend（4）样品典型相关分析代码x=;y=;n1=size(x,2);n2=size(y,2);x=zscore(x);y=zscore(y); %标准化数据n=size(x,1);%a,b返回的是典型变量的系数，r返回的是典型相关系数%u,v返回的是典型变量的值，stats返回的是假设检验的一些统计量的值a,b,r,u,v,stats=canoncorr(x,y)x_u