新建 文本文档3.txt

详细讲解了模式识别的算法理论

load d:\exercise\iris\train.txt;
load d:\exercise\iris\test.txt;
a=train(:,1:4);
b=test(:,1:4);
a1=diag((a-ones(75,1)*mean(a(1:25,:)))*inv(cov(a(1:25,:)))*(a-ones(75,1)*mean(a(1:25,:)))');
a1(:,2:2)=diag((a-ones(75,1)*mean(a(26:50,:)))*inv(cov(a(26:50,:)))*(a-ones(75,1)*mean(a(26:50,:)))');
a1(:,3:3)=diag((a-ones(75,1)*mean(a(51:75,:)))*inv(cov(a(51:75,:)))*(a-ones(75,1)*mean(a(51:75,:)))');
a2=exp(-a1/2);
a2(:,4:4)=train(:,8:8);
a2(:,1:1)=a2(:,1:1)/((2*pi)*(2*pi)*sqrt(det(cov(a(1:25,:)))));
a2(:,2:2)=a2(:,2:2)/((2*pi)*(2*pi)*sqrt(det(cov(a(26:50,:)))));
a2(:,3:3)=a2(:,3:3)/((2*pi)*(2*pi)*sqrt(det(cov(a(51:75,:)))));
fid=fopen('d:\exercise\iris\resulttr.txt','w');
fprintf(fid,'%g %g %g %g\n',a2'); 
b1=diag((b-ones(75,1)*mean(a(1:25,:)))*inv(cov(a(1:25,:)))*(b-ones(75,1)*mean(a(1:25,:)))');
b1(:,2:2)=diag((b-ones(75,1)*mean(a(26:50,:)))*inv(cov(a(26:50,:)))*(b-ones(75,1)*mean(a(26:50,:)))');
b1(:,3:3)=diag((b-ones(75,1)*mean(a(51:75,:)))*inv(cov(a(51:75,:)))*(b-ones(75,1)*mean(a(51:75,:)))');
b2=exp(-b1/2);
b2(:,1:1)=b2(:,1:1)/((2*pi)*(2*pi)*sqrt(det(cov(a(1:25,:)))));
b2(:,2:2)=b2(:,2:2)/((2*pi)*(2*pi)*sqrt(det(cov(a(26:50,:)))));
b2(:,3:3)=b2(:,3:3)/((2*pi)*(2*pi)*sqrt(det(cov(a(51:75,:)))));
b2(:,4:4)=test(:,8:8);
fid=fopen('d:\exercise\iris\resultte.txt','w');
fprintf(fid,'%g %g %g %g\n',b2');
b2

b2 =

    2.1389    0.0000    0.0000
    2.9367    0.0000    0.0000
   13.2024    0.0000    0.0000
    8.7521    0.0000    0.0000
    7.8305    0.0000    0.0000
    8.3579    0.0000    0.0000
    0.2340    0.0000    0.0000
    0.0036    0.0000    0.0000
    0.3230    0.0000    0.0000
    4.6323    0.0000    0.0000
    1.7610    0.0000    0.0000
    0.7824    0.0000    0.0000
    4.6323    0.0000    0.0000
    4.6392    0.0000    0.0000
   14.7704    0.0000    0.0000
    9.1711    0.0000    0.0000
    0.0000    0.0000    0.0000
    4.0869    0.0000    0.0000
    0.0037    0.0000    0.0000
    0.2313    0.0000    0.0000
    0.9785    0.0000    0.0000
    1.3524    0.0000    0.0000
   10.9982    0.0000    0.0000
    8.8849    0.0000    0.0000
   13.3270    0.0000    0.0000
    0.0000    1.4324    0.0000
    0.0000    1.0301    0.0001
    0.0000    0.5964    0.0168
    0.0000    3.0454    0.0077
    0.0000    0.2939    0.0000
    0.0000    1.9205    0.0000
    0.0000    1.0726    0.0000
    0.0000    2.4763    0.0000
    0.0000    0.0564    0.3569
    0.0000    0.1785    0.0049
    0.0000    0.4112    0.0003
    0.0000    1.8853    0.0001
    0.0000    0.4330    0.0001
    0.0000    1.2075    0.0001
    0.0000    2.4358    0.0012
    0.0000    0.3246    0.0015
    0.0000    2.3052    0.0013
    0.0000    3.3088    0.0000
    0.0000    0.3283    0.0000
    0.0000    2.8360    0.0015
    0.0000    0.5741    0.0000
    0.0000    2.2986    0.0002
    0.0000    3.4096    0.0000
    0.0000    0.0131    0.0000
    0.0000    3.3362    0.0002
    0.0000    0.0010    0.2338
    0.0000    0.0815    0.2095
    0.0000    0.1765    0.3969
    0.0000    0.0000    1.8708
    0.0000    0.0026    0.0583
    0.0000    0.0000    0.5099
    0.0000    0.0001    0.0052
    0.0000    0.0000    1.2239
    0.0000    0.2368    0.0890
    0.0000    0.0000    0.0029
    0.0000    0.0000    0.0177
    0.0000    0.0000    0.3596
    0.0000    0.0132    0.3886
    0.0000    0.1623    0.3213
    0.0000    0.0007    0.2150
    0.0000    0.0000    0.4005
    0.0000    0.0000    0.0009
    0.0000    0.0012    1.1992
    0.0000    0.0000    1.4643
    0.0000    0.0000    0.3437
    0.0000    0.0000    0.0259
    0.0000    0.0005    0.1867
    0.0000    0.0055    0.7302
    0.0000    0.0000    0.4384
    0.0000    0.0283    0.3892