新建 文本文档1.txt

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

load d:\exercise\2dvow\train.txt;
load d:\exercise\2dvow\test.txt;
a0=train(:,1:2);
b0=test(:,1:2);
a=(a0-ones(338,1)*mean(a0))./(ones(338,1)*std(a0));
b=(b0-ones(333,1)*mean(a0))./(ones(333,1)*std(a0));
a1=diag((a-ones(338,1)*mean(a(1:35,:)))*inv(cov(a(1:35,:)))*(a-ones(338,1)*mean(a(1:35,:)))');
a1(:,2:2)=diag((a-ones(338,1)*mean(a(36:72,:)))*inv(cov(a(36:72,:)))*(a-ones(338,1)*mean(a(36:72,:)))');
a1(:,3:3)=diag((a-ones(338,1)*mean(a(73:112,:)))*inv(cov(a(73:112,:)))*(a-ones(338,1)*mean(a(73:112,:)))');
a1(:,4:4)=diag((a-ones(338,1)*mean(a(113:148,:)))*inv(cov(a(113:148,:)))*(a-ones(338,1)*mean(a(13:148,:)))');
a1(:,5:5)=diag((a-ones(338,1)*mean(a(149:181,:)))*inv(cov(a(149:181,:)))*(a-ones(338,1)*mean(a(149:181,:)))');
a1(:,6:6)=diag((a-ones(338,1)*mean(a(182:204,:)))*inv(cov(a(182:204,:)))*(a-ones(338,1)*mean(a(182:204,:)))');
a1(:,7:7)=diag((a-ones(338,1)*mean(a(205:241,:)))*inv(cov(a(205:241,:)))*(a-ones(338,1)*mean(a(205:241,:)))');
a1(:,8:8)=diag((a-ones(338,1)*mean(a(242:276,:)))*inv(cov(a(242:276,:)))*(a-ones(338,1)*mean(a(242:276,:)))');
a1(:,9:9)=diag((a-ones(338,1)*mean(a(277:309,:)))*inv(cov(a(277:309,:)))*(a-ones(338,1)*mean(a(277:309,:)))');
a1(:,10:10)=diag((a-ones(338,1)*mean(a(310:338,:)))*inv(cov(a(310:338,:)))*(a-ones(338,1)*mean(a(310:338,:)))');
a2=exp(-a1/2);
a2(:,1:1)=a2(:,1:1)/((2*pi)*sqrt(det(cov(a(1:35,:)))));
a2(:,2:2)=a2(:,2:2)/((2*pi)*sqrt(det(cov(a(36:72,:)))));
a2(:,3:3)=a2(:,3:3)/((2*pi)*sqrt(det(cov(a(73:112,:)))));
a2(:,4:4)=a2(:,4:4)/((2*pi)*sqrt(det(cov(a(113:148,:)))));
a2(:,5:5)=a2(:,5:5)/((2*pi)*sqrt(det(cov(a(149:181,:)))));
a2(:,6:6)=a2(:,6:6)/((2*pi)*sqrt(det(cov(a(182:204,:)))));
a2(:,7:7)=a2(:,7:7)/((2*pi)*sqrt(det(cov(a(205:241,:)))));
a2(:,8:8)=a2(:,8:8)/((2*pi)*sqrt(det(cov(a(242:276,:)))));
a2(:,9:9)=a2(:,9:9)/((2*pi)*sqrt(det(cov(a(277:309,:)))));
a2(:,10:10)=a2(:,10:10)/((2*pi)*sqrt(det(cov(a(310:338,:)))));
a2(:,11:11)=train(:,13:13);
fid=fopen('d:\exercise\2dvow\resulttr.txt','w');
fprintf(fid,'%g %g %g %g %g %g %g %g %g %g %g\n',a2');
fid=fopen('d:\exercise\2dvow\resulttr.txt','w');
fprintf(fid,'%g %g %g %g %g %g %g %g %g %g %g\n',a2');
b1=diag((b-ones(333,1)*mean(a(1:35,:)))*inv(cov(a(1:35,:)))*(b-ones(333,1)*mean(a(1:35,:)))');
b1(:,2:2)=diag((b-ones(333,1)*mean(a(36:72,:)))*inv(cov(a(36:72,:)))*(b-ones(333,1)*mean(a(36:72,:)))');
b1(:,3:3)=diag((b-ones(333,1)*mean(a(73:112,:)))*inv(cov(a(73:112,:)))*(b-ones(333,1)*mean(a(73:112,:)))');
b1(:,4:4)=diag((b-ones(333,1)*mean(a(113:148,:)))*inv(cov(a(113:148,:)))*(b-ones(333,1)*mean(a(13:148,:)))');
b1(:,5:5)=diag((b-ones(333,1)*mean(a(149:181,:)))*inv(cov(a(149:181,:)))*(b-ones(333,1)*mean(a(149:181,:)))');
b1(:,6:6)=diag((b-ones(333,1)*mean(a(182:204,:)))*inv(cov(a(182:204,:)))*(b-ones(333,1)*mean(a(182:204,:)))');
b1(:,7:7)=diag((b-ones(333,1)*mean(a(205:241,:)))*inv(cov(a(205:241,:)))*(b-ones(333,1)*mean(a(205:241,:)))');
b1(:,8:8)=diag((b-ones(333,1)*mean(a(242:276,:)))*inv(cov(a(242:276,:)))*(b-ones(333,1)*mean(a(242:276,:)))');
b1(:,9:9)=diag((b-ones(333,1)*mean(a(277:309,:)))*inv(cov(a(277:309,:)))*(b-ones(333,1)*mean(a(277:309,:)))');
b1(:,10:10)=diag((b-ones(333,1)*mean(a(310:338,:)))*inv(cov(a(310:338,:)))*(b-ones(333,1)*mean(a(310:338,:)))');
b2=exp(-b1/2);
b2(:,1:1)=b2(:,1:1)/((2*pi)*sqrt(det(cov(a(1:35,:)))));
b2(:,2:2)=b2(:,2:2)/((2*pi)*sqrt(det(cov(a(36:72,:)))));
b2(:,3:3)=b2(:,3:3)/((2*pi)*sqrt(det(cov(a(73:112,:)))));
b2(:,4:4)=b2(:,4:4)/((2*pi)*sqrt(det(cov(a(113:148,:)))));
b2(:,5:5)=b2(:,5:5)/((2*pi)*sqrt(det(cov(a(149:181,:)))));
b2(:,6:6)=b2(:,6:6)/((2*pi)*sqrt(det(cov(a(182:204,:)))));
b2(:,7:7)=b2(:,7:7)/((2*pi)*sqrt(det(cov(a(205:241,:)))));
b2(:,8:8)=b2(:,8:8)/((2*pi)*sqrt(det(cov(a(242:276,:)))));
b2(:,9:9)=b2(:,9:9)/((2*pi)*sqrt(det(cov(a(277:309,:)))));
b2(:,10:10)=b2(:,10:10)/((2*pi)*sqrt(det(cov(a(310:338,:)))));
b2(:,11:11)=test(:,13:13);
fid=fopen('d:\exercise\2dvow\resultte.txt','w');
fprintf(fid,'%g %g %g %g %g %g %g %g %g %g %g\n',b2');
fid=fopen('d:\exercise\2dvow\resultte.txt','w');
fprintf(fid,'%g %g %g %g %g %g %g %g %g %g %g\n',b2');

load d:\exercise\2dvow\resulttr.txt;
load d:\exercise\2dvow\resultte.txt;
fid=fopen('d:\exercise\2dvow\resultr.txt','w');
fprintf(fid,'%7.3f %7.3f %7.3f %7.3f %7.3f %7.3f %7.3f %7.3f %7.3f %7.3f %g\n',resulttr');
fid=fopen('d:\exercise\2dvow\resulte.txt','w');
fprintf(fid,'%7.3f %7.3f %7.3f %7.3f %7.3f %7.3f %7.3f %7.3f %7.3f %7.3f %g\n',resultte');