1、模式识别程序及其运行结果第一题:二类协方差不等:clc;clear all;A(1).t=1 1 2;1 0 -1;A(2).t=-1 -1 -2;1 0 -1;pw=1/2,1/2;x=2;0;for i=1:2a(i).t=mean(A(i).t,2);b(i).t=A(i).t-a(i).t*ones(1,3);sigma(i).t=b(i).t*(b(i).t)/2;endsig_in1=inv(sigma(1).t);sig_in2=inv(sigma(2).t);f=0.5*(x-a(1).t)*sig_in1*(x-a(1).t)-0.5*(x-a(2).t)*sig_in2*(
2、x-a(2).t)+0.5*log(det(sigma(1).t)/det(sigma(2).t);i=1;if f0i=2enddisp(二类协方差相等时 x 属于的类别为:);disp(i);%界面方程y=0.5*(sym(x1;x2)-a(1).t)*sig_in1*(sym(x1;x2)-a(1).t)-0.5*(sym(x1;x2)-a(2).t)*sig_in2*(sym(x1;x2)-a(2).t)+0.5*log(det(sigma(1).t)/det(sigma(2).t);%绘制分界线L1=ezplot(y,-10,10);set(L1,Color,r);title(二类协
3、方差相等时的分界线);legend(第 1、2 类分类线,location,NorthEastOutside);hold on;plot(A(1).t(1,:),A(1).t(2,:),r+);hold on;plot(A(2).t(1,:),A(2).t(2,:),go);hold on;plot(x(1,:),x(2,:),k.,linewidth,2);第二题:三类协方差相等:clc;clear all;A(1).t=0 2 1;0 1 0;A(2).t=-1 -2 -2;1 0 -1;A(3).t=0 0 1;-2 -1 -2;pw=1/3,1/3,1/3;x=-2;2;for i=1
4、:3a(i).t=mean(A(i).t,2);b(i).t=A(i).t-a(i).t*ones(1,3);sigma(i).t=b(i).t*(b(i).t)/3;endsig=sigma(1).t+sigma(2).t+sigma(3).t;sig_in=inv(sig)f=zeros(1,3);for i=1:3f(i)=-0.5*(x-a(i).t)*sig_in*(x-a(i).t);endm=max(f);i=1;while(f(i)=m)i=i+1;enddisp(三类协方差相等时 x 属于的类别为:);disp(i);%界面方程y1=0.5*(sym(x1;x2)-a(1).
5、t)*sig_in*(sym(x1;x2)-a(1).t);y2=0.5*(sym(x1;x2)-a(2).t)*sig_in*(sym(x1;x2)-a(2).t);y3=0.5*(sym(x1;x2)-a(3).t)*sig_in*(sym(x1;x2)-a(3).t);y12=y1-y2;y13=y1-y3;y23=y2-y3;%绘制分界线 L1=ezplot(y12,-10,10);set(L1,Color,r);hold on;L2=ezplot(y13,-10,10);set(L2,Color,g);hold on;L3=ezplot(y23,-10,10);set(L3,Color,b);hold offtitle(三类协方差相等时的分界线);legend( 红色为第 1、2 类的分类线 ,绿色为第 1、3 类的分类线,蓝色为第 2、 3 类的分类线,location,NorthEastOutside);hold on;plot(A(1).t(1,:),A(1).t(2,:),r+);hold on;plot(A(2).t(1,:),A(2).t(2,:),go);hold on;plot(A(3).t(1,:),A(3).t(2,:),b*);hold on;plot(x(1,:),x(2,:),k.,linewidth,2);