matlab数据分析方法 ppt,MATLAB数据分析方法-(3).ppt

称 称 为n维向量x，y之间的闵可夫斯基距离，其中 为常数。 即当 到 的马氏距离不超过到 的马氏距离时，判 来自 ；反之，判来自 . 解：假定两总体的协方差相等，源程序如下： apf=[1.14,1.78;1.18,1.96;1.20,1.86;1.26,2.;1.28,2;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]; % 输入原始数据 m1=mean(apf); m2=mean(af); s1=cov(apf); s2=cov(af); s=(5*s1+8*s2)/13; % 计算样本均值与协方差矩阵 for i=1:3 W(i)=(x(i,:)-1/2*(m1+m2))*inv(s)*(m1-m2)'; % 计算判别函数值 end; 直接调用MATLAB的判别分析命令classify。 apf=[1.14,1.78;1.18,1.96;1.20,1.86;1.26,2.;1.28,2;1.30,1.96]; %总体apf 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]; %总体af training=[apf;af]; %合并两个总体形成训练集 n1=size(apf,1); %总体apf中样本的行数 n2=size(af,1); %总体af中样本的行数 group=[ones(1,n1), 2*ones(1,n2)]; %apf中样本与af中样本类属 x=[1.24,1.8;1.28,1.84; 1.4,2.04]; % 输入原始待判数据即sample class = classify(x, training,group) %判别分析 输出结果为： W = 1.7611 3.8812 3.6468 由判别准则(4.1.17)可知，三个蠓虫均属于Af. 解：根据例2.2.3的结论，可以认为三类总体协方差 矩阵相等. A=[2607540183101223021320643917;…… 26013539292804037172501173616]; G1=A(:,1:4);G2=A(:,5:8);G3=A(:,9:12); %三类总体数据 x=[190 67 30 17;315 100 35 19;240 60 37 18]; %待判定的数据 m(1,:)=mean(G1);m(2,:)=mean(G2);m(3,:)=mean(G3); s1=cov(G1);s2=cov(G2);s3=cov(G3); s=19*(s1+s2+s3)/57; for i=1:3 for j=1:3 for k=1:3 w(j,k)=(x(i,:)-1/2*(m(j,:)+m(k,:)))*inv(s)*(m(j,:)-m(k,:))'; if w(j,k)<0 q=0;break; el