数学建模数据分析题(共25页)

1、精选优质文档-倾情为你奉上中国矿业大学数学建模常规赛竞赛承 诺 书我们仔细阅读了中国矿业大学数学建模常规赛论文格式规范和2016年中国矿业大学数学建模常规赛通知。我们完全明白，在竞赛开始后参赛队员不能以任何方式（包括电话、电子邮件、网上咨询等）与队外的任何人（包括指导教师）研究、讨论与赛题有关的问题。我们知道，抄袭别人的成果是违反竞赛章程和参赛规则的，如果引用别人的成果或资料（包括网上资料），必须按照规定的参考文献的表述方式列出，并在正文引用处予以标注。在网上交流和下载他人的论文是严重违规违纪行为。我们以中国矿业大学大学生名誉和诚信郑重承诺，严格遵守竞赛章程和参赛规则，以保证竞赛的公正、公平性

2、。如有违反竞赛章程和参赛规则的行为，我们将受到严肃处理。我们授权中国矿业大学数学建模协会，可将我们的论文以任何形式进行公开展示（包括进行网上公示，在书籍、期刊和其他媒体进行正式或非正式发表等）。我们的参赛队号： 25 参赛队员 (打印并签名)：1. 易阳俊 2. 令月霞 3. 刘景瑞 日期： 2016 年 10 月 日（请勿改动此页内容和格式。此承诺书打印签名后作为纸质论文的封面。以上内容请仔细核对，如填写错误，论文可能被取消评奖资格。）中国矿业大学数学建模常规赛竞赛编 号 专 用 页评阅记录（可供评阅时使用）：评阅人备注评阅统一编号（数学建模协会填写）：题目：数据的分析问题摘要本文需要解决的

3、问题是如何根据就诊人员体内7种元素含量来判别某人是否患有疾病G和确定哪些指标是影响人们患疾病G的主要因素。通过解读题目可知，此类问题为典型的分析判别问题。我们先对数据进行了预处理，剔除了有异常数据的样本，然后采用元素分布判别法、马氏距离判别法和Fisher判别法，应用Excel、SPSS和MATLAB等软件来对某人是否患病进行判别，并通过绘制7种元素含量的折线图等来确定患该疾病的主要因素，最后应用综合判别法对之前的结论进行了检验。对于问题一，在对数据预处理之后，我们删除了序号为10这个高度异常数据样本，然后我们分别采用元素分布判别法、马氏距离判别法和Fisher判别法对49个已知病例进行判别。

4、对于元素分布判别法，我们通过数据预处理知道7种元素含量分布均符合正态分布，然后我们确定了以均值为大致中心的元素正常含量范围，得出其判别准确度为96%；对于马氏距离判别法，通过编写MATLAB程序（见附录）来进行判别，得出其判别准确度为90%；对于Fisher判别法，通过SPSS软件来进行判别，得到线性判别函数，其判别准确度为96%；针对问题二：我们运用问题一中建立的三个判别模型对25名就诊人员（见附录）的化验结果进行检验，判别结果如下表1:表1判别方法元素分布判别法马氏距离判别法Fisher判别法健康人数151411患病人数101114对于问题三，我们运用Excel对健康人和患者的7种元素分别

5、绘制折线图进行对分析，我们初步判定元素4与元素5是影响人们患疾病G的主要因素，然后用方法一的三种判别方法进行检验，其准确度在85%以上；对于问题四，我们根据问题三得出的主要因素，分别用三种判别方法对25名就诊人员进行判别，再与问题二的判别结果进行对比，可知它们判断结果之间的差异性最高为24%。对于问题五，由于三种判别法都有不足，所以我们采用了综合判别法，将三种判别方法的结果进行综合判断，最终我们通过主要因素进行判别的差异性下降到了12%，与问题一的判断结果的一致性达到了88%。关键词：马氏距离判别，Fisher判别，综合判别，MATLAB，SPSS一 、问题重述.随着大数据时代的到来，人们在处

6、理问题时往往要借助一些实验或检测数据的分析；病人到医院就诊时，医生通常要通过化验和检测得到的数据分析来协助诊断。在医疗诊断中，要诊断一个人是否患上疾病G时,通常要检测人体内7种元素的含量。表1是50个确诊病例的检测结果，其中125号病例是已经确诊为疾病G的病人的检测结果；2650号病例是已经确定为健康人的检测结果。表2是25个就诊人员的检测结果。试解决下列问题：问题1：根据表1中的数据，提出一种或多种简便的判别方法，判别属于疾病G的病人或健康人的方法，并检验你提出方法的正确性。问题2：按照问题1提出的方法，对表2中的25名就诊人员的检测结果进行判别，判定他（她）们是疾病G的病人还是健康人。问题

7、3：能否根据表1的数据特征，确定哪些指标是影响人们患疾病G的关键或主要因素，以便减少化验的指标。问题4：根据问题3的结果，重复问题2的工作。问题5：对问题2和问题4的结果作进一步的分析。二、模型假设(1)假设题目中所给的数据绝大多数真实可靠，不排除出现少数异常数据的可能；(2)假设就诊人员在化验前不会采取增加体内这7种元素的措施；(3)题目中所给的样本只有因这7种元素含量而患G病的患者或者为健康人员，不会受到其他疾病和其他元素的干扰；(4)不考虑各种元素之间的相互作用对机体产生的影响； 三、问题的分析对于此题，我们需要通过对已知病例数据进行分析，从而找到疾病确诊方法。我们首先运用Excel对健

8、康人的7种元素含量绘制散点图进行分析和W检验，得知这7种元素的分布均符合正态分布，再通过拉依达准则法排除了高度异常数据10号样本。针对问题一，我们建立了元素分布判别、马氏距离判别和Fisher判别三种模型。首先，已知7种元素含量分布均符合正态分布，于是我们以健康人员的7种元素含量的均值为大致中心值，确立了一个元素正常含量范围，对于只要存在一个不在正常范围内的元素含量的待诊人员，我们就判定为患者，从而建立了元素分布判别模型；其次，我们利用马氏距离判别法，对原始数据直接进行分析归类，并排除了元素之间的相关性的干扰，运用MATLAB编程中得到的判别函数来对患者和健康人员这两组样本进行马氏判别分析，从

9、而建立了马氏距离判别模型；然后，我们运用SPSS软件对患者和健康人员这两组样本进行Fisher判别，得到区分人员是否患病的线性判别函数和判断准确率，从而建立了Fisher判别模型。最后，我们用已知的49个病例数据对这三个模型分别进行了检验，得到的判别准确率均在90%以上。 针对问题二，我们利用问题一中建立的三个判别模型分别对25位就诊人员是否患病进行判别，结果见表1: 表1判别方法元素分布判别法马氏距离判别法Fisher判别法健康人数151411患病人数101114针对问题三，为了确定影响人们患病G的关键或主要因素，我们运用Excel图表功能对7种元素分别绘制折线图进行对比。通过分析，可知绝大

10、多数健康人员的元素4和元素5含量均大于患者在这两种元素含量的最大值，健康人员与患者含量呈现明显数值差异；而健康人与患者在其他5种元素上的分布差距没有元素4与元素5如此明显，所以我们初步判定元素4与元素5是影响人们患疾病G的主要因素。然后我们用问题一中的三种判别模型对已知病例进行检验，得到的判别准确率均在85%以上。针对问题四，虽然利用问题三中得到的主要元素来进行分析判别的准确率在85%以上，但在与问题二的判别结果进行对比后，可知两种判别结果的差异性（见下表2）。 表2判别方法差异性元素分布判别法16%马氏距离判别法24%Fisher判别法12%由上表可知，问题二与问题四的判别结果的最大差异性有

11、24%，表明此时利用这两个主要因素进行判别结果无法较好达到问题二的判别效果。针对问题五，已知问题四中得到的两个主要因素的判断结果与问题二的判断结果差异性大，而且这三种判别方法均有各自的缺点，于是我们建立了一个综合判别模型，通过三种判别方法来得到最后的判别结果，最终使得两个判别结果的差异性降到了12%，一致性达到了88%。四、符号及变量说明符号意义显著性水平马氏距离判别函数Fisher判别线性函数Fisher判别临界值五、数据预处理5.1数据分析处理实验数据的时候，我们常常会遇到个别数据值偏离预期或大量统计数据值结果的情况，如果我们把这些数据值和正常数据值放在一起进行统计，可能会影响实验结果的正

12、确性。我们首先通过Excel做出了26-50号健康人的7种元素的散点图（如下图）：图1 图2 图3 图4图5 图6 图7由图可知：1、序号为10的健康人员的元素数据十分异常，我们暂且忽略这个样本。2、我们可以发现忽略掉少数几个异常数据外，每种元素数值大都集中在中间某个数值附近，则我们可以初步判断这7种元素的含量可能服从正态分布，于是我们认为可以通过统计学知识中的W检验来对这些数据分别进行正态性检验。5.2 W检验W检验全称Shapiro-Wilk检验，是一种基于相关性的算法，通过计算可得到一个相关系数，它越接近1就越表明数据和正态分布拟合得越好。W检验是检验样本容量8n 50时，样本是否符合正

13、态分布的一种方法。我们已知25个病人数据和24个健康人员数据，计算公式为:其检验步骤如下： 5.3 正态性检验通过编写MATLAB程序（见附录），我们可以得到健康人员的7种元素的W值（见下表3）： 表3元素j1234567W值0.94290.66740.62620.90510.95480.91250.9187 通过进一步对元素2和元素3的散点图（见表2和表3）的观察可知，其中的一些异常数据可能干扰了判断结果，在先忽略掉序号为19的健康人员的元素2数据和序号为11、21的健康人员的元素3的数据后，再对它们进行一次W检验（见下表3）： 表4元素序号i230.95850.8617认为它也应该满足正态

14、分布。5.4 数据检测拉依达准则法是最常用的异常值判定与剔除准则，优点是简单，无需查表，测量次数较多或要求不高时用。在这种情况下，异常值是指一组测定值中与的超过两倍的测定值。与平均值的偏差超过三倍标准差的测定值，称为高度异常的异常值。在处理时，应剔除高度异常的异常值。异常值是否剔除，视具体情况而定。在统计检验时，指定为检出异常值的显著性水平=0.05，称为检出水平；指定为检出高度异常的异常值的显著性水平=0.01，称为舍弃水平，又称剔除水平(reject level)。由于我们已知7种元素含量分布均符合正态分布，所以我们可以用以下计算公式来判断异常值是否存在，公式如下：当s>0时，为异常

15、值；当s<0时，为正常值；通过Excel的简单编程计算，我们得到以下异常数据表（见表4）: 表5序号i元素j含量类型序号i元素j含量类型11265.4高度异常102135高度异常193308高度异常103322高度异常91247异常1046747高度异常1245361异常1051090高度异常187899异常106228异常197874异常当c=3时，序号为10的健康人员在4个元素检测上均为高度异常值，故我们可以排除该样本。当c=2时，有6个健康人员也存在一些相对异常的数据，我们可以在对这些健康人员进行数据分析处理时，选择忽略掉这些异常数据，以免对最终分析结果产生较大的影响。 六、问题模

16、型建立与求解6.1问题一的模型建立6.1.1 元素分布判别模型的建立我们已知这7个元素分布符合正态分布，由正态分布关于平均值对称的性质，于是我们以平均值为大致中心值，同时尽量保证两端点值与均值的距离相差不大，以此来确定元素正常含量范围，建立元素分布判别模型，从而判断待诊人员是否患病；只要待诊人员的某一种元素含量在该范围之外时，我们就判定他为病人。根据元素的散点图分布和拉依达准则法得到的异常数据表(见上表5)，我们得到元素正常含量范围:表6元素序号i1234567均值（取整数）1851857226526793388最大值（不含异常值）235251553879432182739最小值MIN1336

17、.438.65102599.13167.3正常含量范围133-2356.43-258.65-1051025-387999-43231-18267.3-739然后我们用已知的49个病例来对模型进行检验，得到综合准确率为96%的判别模型，检验结果如下(见表7)： 表7已知病例检验成功人数检验成功率健康人员2395.8%患者2496%6.2马氏距离判别模型6.2.1基本思想首先根据已知分类的数据，分别计算各类中心即分组（类）的均值。判别准则是对任意给的一次观测，若它与第i类的中心距离最近，就认为它来自第i类。6.2.2 建立过程,按照距离最近准则判别归类，则可以写成： 6.2.3 马氏距离判别分析和

18、检验我们通过MATLAB编程（程序见附录）来对已知的49名数据样本进行检测，得到马氏距离检测法的综合检测准确率为90%，检测结果如下： 表8已知病例健康人员患者检测成功人数2420总人数2425检测成功率100%80%6.3 Fisher判别模型6.3.1基本思想从k个总体中抽取p个指标的样品观测数据，借助方差分析的构造一个线性断函数，其中系数确定的原则是使得总体之间区别最大，而是每个总体内部的离差最小。有了线性判断函数U后，对于一个新的样品，将它的p个指标值带入以上线性判断函数式子中求出U值，然后根据一定的判别规则，就可以判断新的样品属于哪个总体了6.3.2 算法步骤6.3.3、进行Fish

19、er判别分析和检验我们运用SPSS软件对已知病例数据进行Fisher判别分析，得到判别结果如下（详细表格见附录）：表9分類結果a 组别預測的群組成員資格總計1.02.0原始計數1.0241252.012324%1.096.04.0100.02.04.295.8100.0a. 95.9% 個原始分組觀察值已正確地分類。表10標準化典型區別函數係數函數1元素1-.183元素2-.023元素3.521元素4.546元素5.724元素6.348元素7-.841（1）由分类结果可知：通过对已知的49名病例样本进行检测，对健康人与患病者的检测准确率都达到了96%，故综合检测准确率为96%。（2）由标准化典

20、型区别函数系数表可知，线性判别函数D应为 :虽然我们得到了Fisher判断法的线性判断函数，但我们还需找到一个判定标准来对健康人员与患者进行区分，通过表9可知，分别有一名已知的健康人员与患者在运用Fisher进行判别时出现错误，所以我们可以根据这两名病例的判断得分来得到判定标准。我们运用Excel进行简单编程，求出已知49名病例的Fisher判断函数的得分，结果如下：图8由图可知：病例号为1-25的患者整体得分明显少于病例号为26-50的健康人员整体得分，显然病例号为23和38的得分情况是Fisher判别所出现的两个判断错误；所以病例号为23的病例在Fisher判别法认为是健康人，而得分略低的

21、病例号为18的病例被 Fisher判别法认为是患者，对于我们所设置的判别临界值S，他应该属于165.6958241.8355之间，不妨令6.4 问题二的求解 综述以上三种判别方法，我们得到了它们的判别准确率，如下表：表11准确度/种类元素分布判别马氏距离判别Fisher判别法病人95.83%100%95.83%健康人员96%80%96%综合96%90%96%我们采用问题一中的元素分布判别法、马氏距离判别法和Fisher判别法三种方法，运用Excel和MATLAB软件(代码详见附录)分别这25名就诊人员进行判别，再将结果整理如下：表12元素分布判别法马氏距离判别法Fisher判别法患病人员51-

22、59、61-63、66、73、7451-54、56、61-63、66、73、7551-56、58、59、61-63、66、73、75健康人员60、64、65、67-72、7555、57-60、62、64、65、67-72、7457、60、64、65、67-72、74患病与健康人员的人数之比15：1011：1414：116.5问题三的模型建立与问题四的模型求解6.5.1 主要元素的确定我们运用Excel分别对健康人与患者的7种元素含量进行了升序排列，然后对健康人和患者的7种元素分别绘制折线图进行对比，通过对元素4与元素5的折线图（见下图，其他5种元素图见附录）观察分析可知，绝大多数健康人的这两种

23、元素含量均大于患者在这两种元素含量的最大值，呈现出明显的数值大小差异，而健康人与患者在其他5种元素上的分布，他们之间的差距没有元素4与元素5明显，所以我们初步判定元素4与元素5是影响人们患疾病G的主要因素。图9 图10我们进一步通过对图9的观察可知，患者的元素4含量大多数在100-160之间，而健康人仅有两名在160以下，于是我们对问题一中元素分布判别法进行修改，对于元素4，我们以162为新的边界值来进行判别，即元素4的新范围为162-432。6.5.2 利用主要元素进行判别和检验我们利用问题三的主要元素，通过问题一中三种判别法来重新对就诊人员进行判别，经过整理后，判别结果如下： 表13准确度

24、/种类元素分布判别马氏距离判别Fisher判别法病人95.83%91.6%100%健康人91.6%79.2%79.2%综合93.7%85.4%89.6% 表14元素分布判别法马氏距离判别法Fisher判别法患病人员51、52、54-59、62、63、66、73、7551、52、54-56、59、63、66、69、73、7551-59、61-63、66、67、69、73、75健康人员53、61、64、65、67-72、7453、57、58、60-62、64、65、67、68、70-72、7460、64、65、68、70-72、74患病与人数之比14：1111：1417：8我们进一步将该判别结果与

25、问题2的判别结果进行对比，如下表：表14判别差异病例号差异性元素分布判别法53、61、74、7516%马氏距离判别法53、55、59、61、62、6924%Fisher判别法57、67、6912%由上表可知，虽然只利用主要因素的三种判别法的有效率仍然在85%以上，但是三种判别法的结果与问题二中的判别结果仍然有不小的差异，虽然Fisher判别法差异性最小只有12%，但马氏距离判别法的差异性却有24%，我们需要进一步降低差异率来确保只利用主要因素进行判别的判断率。 6.6问题五的分析虽然问题二与问题三的三种判别法的有效率都在85%以上，但是两个问题中的三种判别法的结果仍然有不小的差异，由于三种判别

26、法均存在各自的问题，于是我们建立了一个综合判别模型，通过三种判别方法来得到最后的判别结果。我们先将问题二与问题四的判别结果进行整理，为了简化表格，给出如下定义：当某个判别方法判定就诊人员为健康人员时，得分加1；当某个判别方法判定就诊人员为患者时，不得分。表15问题二问题四序号元素分布马氏距离Fisher判别综合得分元素分布马氏距离Fisher判别综合得分51000000005200000000530000110254000000005501010000560000000057011201015801010101590101000060111311126100001102620000010163

27、00000000641113111365111311136600000000671113110268111311136911131001701113111371111311137211131113730000000074011211137510010000由上表可知：有一些就诊人员在三种判别法中均被判定为患者，也有一些人在三种判别方法法中，均被判定为健康人员。所以我们根据得分的高低情况可以将上述就诊人员分为以下四类水平：完全健康（3分）、大概率健康（2分）、大概率患病（1分）、完全患病（0分），见下表16。表16范围完全健康大概率健康大概率患病安全患病问题二60、64、65、67-7257、7

28、455、58、59、7551-54、56、61-63、66、73问题四64、65、68、70-72、7453、60、61、6757、58、62、6951、52、54-56、59、63、66、73、75由于我们必须判断出就诊人员是否患病，所以我们根据大概率原则，可以判断完全健康、大概率健康这两种水平的就诊人员为健康，而完全患病、大概率患病这两种水平的就诊人员为患病，所以我们得到以下表格：表17问题二问题四差异病例号健康患病健康患病53、57、69（差异性为12%）57、60、64、65、67-72、7451-56、58、59、61-63、66、73、7553、60、61、64、65、67、68、

29、70-72、7451-52、54-59、62、63、66、69、73、75所以通过综合判别法，我们使问题二与问题四的差异性降到了12%，一致率达到了88%，所以此时元素4、元素5可以较好的作为主要元素来对就诊人员进行检测。 七、模型评价及推广7.1 模型评价7.1.1优点：针对问题一，为提出一种简便的判别方法，我们提出元素分布判别法马氏距离判别法和Fisher判别法，并进行了准确度的检验，得出元素分布判别法和Fisher判别法均有96%的准确度，马氏距离判别法也有90%的准确度，三种判别方法的准确度很高，其大致符合医院的实际情况，而且施行相对简便准确。针对问题二，在问题一的条件下可将就诊人员的

30、化验结果用元素分布判别法马氏距离判别法和Fisher判别法分别进行判别，从而得出结果，具有良好的通用性和实用性。针对问题三，我们运用Excel图表功能对7种元素分别绘制折线图进行对比，初步判定元素4与元素5是影响人们患疾病G的主要因素，然后三种判别模型进行检验，得到的判别准确率均在85%以上，方法简便且效率高。针对问题四和问题五， 已知问题四中得到的两个主要因素的判断结果与问题二的判断结果差异性大，而且这三种判别方法均有各自的缺点，于是我们建立了一个综合判别模型，通过三种判别方法来得到最后的较好判别结果，判别方法进行组合，优势互补，判断准确性高 7.1.2缺点：由于题中所给数据相对较少，且化验

31、元素种类也较少，对有一定影响，以至元素分布判别法马氏距离判别法和Fisher判别法所得结果与实际有差别。7.1.3模型推广多种判别方法相结合的综合判别法不仅可以用于对某种疾病的诊断，而且可以用于在经济学中根据人均国民收入、人均工农业产值、人均消费水平等指标判断一个国家的经济发展程度。多种判别方法相结合的综合判别法的结合可以处理有较大量的数据资料，且机理不甚清楚的问题，并能够进行有效的分析和提炼出关键因素，从而找到内在规律，对问题做出科学判断。 八、参考文献1 卓金武，MATLAB在数学建模中的应用，北京：北京航空航天大学出版社，20102 茆诗松程依明濮晓龙，概率论与数理统计，北京：高等教育出

32、版社，20043 姜启源 谢金星 数学模型，北京：高等教育出版社，20034 邰淑彩，应用数理统计，武汉：武汉大学出版社，20055韩中庚，数学建模方法及其应用，北京：高等教育出版社，2005九、附录9.1 MATLAB代码1：W检验程序%健康人员的元素1数据A=8.65 14.117 17.718.619.823.42728.229.8313535.5 3636.2 3761.762.864.28686.990.8;b=0;n=length(A);for p=1:n b=b+A(p);endb=b/n;%另一组正态性检验统计量W的系数数值表，a=0.4493 0.3098 0.2554 0.

33、2145 0.1807 0.1512 0.1245 0.0997 0.0764 0.0539 0.0321 0.0107;a=0.459 0.3156 0.2571 0.2131 0.1764 0.1443 0.115 0.0878 0.0618 0.0368 0.0122;s=0;for k=1:n/2d(k)=A(n-k)-A(k);s=s+a(k)*(A(k)-A(n-k+1);endm=0;for i=1:n m=m+(A(i)-b)*(A(i)-b);end%W检验的计算结果w=s*s/m9.2 MATLAB代码2：马氏距离判别法程序clcclear%矩阵A、B是已知病例的数据，X是

34、就诊人员的数据A= 16615.824.5700112179513;18515.731.5701125184427;1939.8025.9541163128642;15914.239.789699.2239726;22616.223.860615270.3218;1719.299.2930718745.5257;20113.326.655110149.4141;14714.530.0659102154680;1728.857.8655175.798.4318;15611.532.5639107103552;13215.917.757892.413141372;18211.311.3767111

35、264672;1869.2637.195823373.0347;1628.2327.162510862.4465;1506.6321.0627140179639;15910.711.761219098.5390;11716.17.0498895.5136572;18110.14.041437184101542;14620.723.812321281501092;42.310.39.7062993.7439888;28.212.453.137044.1454852;15413.853.3621105160723;17912.217.9113915045.2218;13.53.3616.81353

36、2.651.6182;1755.8424.980712355.6126;B= 21319.136.2222024940.0168;17013.929.8128522647.9330;16213.219.8152116636.2133;20313.090.8154416298.90394;16713.114.1227821246.3134;16412.918.6299319736.394.5;16715.027.0205626064.6237;15814.437.0102510144.672.5;13322.831.01633401180899;15613532267471090228810;1

37、698.00308106899.153.0289;24717.38.65255424177.9373;1668.1062.81233252134649;2096.4386.9215728874.0219;1826.4961.73870432143367;23515.623.4180616668.8188;17319.117.0249729565.8287;15119.764.22031403182874;19165.435.05361392137688;22324.486.0360335397.7479;22120.11553172368150739;21725.028.22343373110

38、494;16422.235.52212281153549;1738.9936.01624216103257;20218.617.7378522531.067.3;X= 16615.824.5700112179513;18515.731.5701125184427;1939.8025.9541163128642;15914.239.789699.2239726;22616.223.860615270.3218;1719.299.2930718745.5257;20113.326.655110149.4141;14714.530.0659102154680;1728.857.8655175.798

39、.4318;15611.532.5639107103552;13215.917.757892.413141372;18211.311.3767111264672;1869.2637.195823373.0347;1628.2327.162510862.4465;1506.6321.0627140179639;15910.711.761219098.5390;11716.17.0498895.5136572;18110.14.041437184101542;14620.723.812321281501092;42.310.39.7062993.7439888;28.212.453.137044.1454852;15413.853.3621105160723;17912.217.9113915045.2218;13.53.3616.813532.651.6182;1755.8424.980712355.6126;Ca=cov(A); Cb=cov(B); CA=inv(Ca); CB=inv(Cb); Ua=sum(A)./25; Ub=sum(B)./24; for i=1:25 Xa=X(i,:)-Ua; DA(i)=Xa*CA*Xa' Xb=X(i,:)-Ub; DB(i)=Xb*CB*Xb' w(i)=DA(i)-DB(i); end for i=1:25