新建 文本文档2.txt

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

load d:\exercise\iris\train.txt;
load d:\exercise\iris\test.txt;
a=train(:,1:4);
a0=(a-ones(75,1)*mean(a))./(ones(75,1)*std(a));
b=test(:,1:4);
a1=diag((a0-ones(75,1)*mean(a0(1:25,:)))*inv(cov(a0(1:25,:)))*(a0-ones(75,1)*mean(a0(1:25,:)))');
a1(:,2:2)=diag((a0-ones(75,1)*mean(a0(26:50,:)))*inv(cov(a0(26:50,:)))*(a0-ones(75,1)*mean(a0(26:50,:)))');
a1(:,3:3)=diag((a0-ones(75,1)*mean(a0(51:75,:)))*inv(cov(a0(51:75,:)))*(a0-ones(75,1)*mean(a0(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 %g\n',a2'); 
b0=(b-ones(75,1)*mean(a))./(ones(75,1)*std(a));
b1=diag((b0-ones(75,1)*mean(a0(1:25,:)))*inv(cov(a0(1:25,:)))*(b0-ones(75,1)*mean(a0(1:25,:)))');
b1(:,2:2)=diag((b0-ones(75,1)*mean(a0(26:50,:)))*inv(cov(a0(26:50,:)))*(b0-ones(75,1)*mean(a0(26:50,:)))');
b1(:,3:3)=diag((b0-ones(75,1)*mean(a0(51:75,:)))*inv(cov(a0(51:75,:)))*(b0-ones(75,1)*mean(a0(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 %g\n',b2');