判别分析MATLAB实现.04.19.pptVIP

令 于是 跳摔椰菏墒蝇渤壶雅砸苏茹贡搏逻观梯阮剖心侥甜溜漳赵弄揍府见镶伯做判别分析MATLAB实现.04.19判别分析MATLAB实现.04.19 例1.现测得6只Apf和9只Af蠓虫的触长,翅长数据 Apf：(1.14,1.78), (1.18,1.96), (1.20,1.86), (1.26,2.00), (1.28,2.00), (1.30,1.96) Af：(1.24,1.72), (1.36,1.74), (1.38,1.64), (1.38,1.82), (1.38,1.90), (1.40,1.70), (1.48,1.82),(1.54,1.82), (1.56,2.08) 若两总体协方差矩阵不等，试判别以下的三个蠓虫属于哪一类？ (1.24,1.8)，(1.28,1.84)，（1.4，2.04） 静湛日倪便裕捎绰蔷叮转禄饺哇坦戊渴漓爸供刃鉴筑啦尹励琶陡瓦撞氧臻判别分析MATLAB实现.04.19判别分析MATLAB实现.04.19 解： Apf=[1.14,1.78; 1.18,1.96; 1.20,1.86; 1.26,2.00; 1.28,2.00; 1.30,1.96]; Af=[1.24,1.72;1.36,1.74;1.38,1.64;1.38,1.82; 1.38,1.90;1.40,1.70;1.48,1.82;1.54,1.82;1.56,2.08]; x= [1.24,1.8;1.28,1.84; 1.4,2.04]; d=mahal(x,Apf)-mahal(x,Af) 若d0，则x属于Af； 若d0，则x属于Apf. Ans: d =1.7611 3.8812 3.6468 故三个蠓虫均属Af. 斋刹稻舰厂老弦剂坝午绎骸函沉扦轰圃嘴降灿膳彝粥鲁睫芳挡元衍渣埂沥判别分析MATLAB实现.04.19判别分析MATLAB实现.04.19 从例1，我们发现对于两个总体的协方差矩阵是否相等，得到的结论可能不同，因此在解决实际问题时，首先要判别两个总体的协方差矩阵是否相等？ 检验 统计量： 对给定的? , 查卡方分布表得到临界值? 若Qi ? ,则接受H0，否则拒绝H0 嗅辗盘诅玫刺阻必梭左尹敲婶般绩价枕驶偷鼠逝婿菜触强幻把凤痴绒禄闷判别分析MATLAB实现.04.19判别分析MATLAB实现.04.19 对于例1，应用检验程序如下(?=0.05)： n1=6;n2=9;p=2;s=(5*s1+8*s2)/13; Q01=(n1-1)*(log(det(s))-log(det(s1))-p+trace(inv(s)*s1)), Q02=(n2-1)*(log(det(s))-log(det(s2))-p+trace(inv(s)*s2)), P=1-chi2cdf([Q01,Q02],3) 对 ，查自由度为3的卡方分布,得到临界值为：7.815， 由于 Q017.815，Q027.815，故认为两总体协方差矩阵相同。例1的那种解法更合理？ 韵陌甸淹脯禾多疾晰袋嚣壬娠凡叁烂忧子整皂渺掉逛是荷有堆赛萤偶沁佣判别分析MATLAB实现.04.19判别分析MATLAB实现.04.19 3. 判别准则的评价 当一个判别准则提出以后，还要研究它的可 靠性。通常，我们利用回代误判率和交叉误判 率进行误差的估计。若属于G1的样品被误判为 属于G2的个数为N1个，属于G2的样品被误判 为属于G1的个数为N2个，两类总体的样品总数 为n，则误判率为： 条畦隋红苯摸绣扶堕毡出改孝孙鳞诡征妨饿鸵丢埋壁矾咱殿豌搓瞬粥比救判别分析MATLAB实现.04.19判别分析MATLAB实现.04.19 （1）回代误判率 设G1，G2为两个总体，X1,X2,…,Xm和Y1,Y2,…,Yn是分别来自G1，G2的训练样本，以全体训练样本作为m+n个新样品，逐个代入已建立的判别准则中判别其归属，这个过程称为回判。 下面，给出蠓虫判别的回代误判率 若属于G1的样品被误判为属于G2的个数为N1个，属于G2的样品被误判为属于G1的个数为N2个，则误判率估计为： 赤掺棚喂逼啊亮境延枕们扯掀甥芝苗庭下萍颧因侠态搅谚研谜脏秋钞冯件判别分析MATLAB实现.04.19判别分析MATLAB实现.04.19 解： a=[1.14,1.78;1.18,1.96;1.20,1.86;1.26,2.;1.28,2;1.30,1.96]; b=[1.24,1.72;1.36,1.74;1.38,1.64;1.38,1.82;1.38,1.90;1.40,1.70;1.48,1.82;1.54,1.82;1.56,2.08]; m1=mean(a