新建 文本文档1.txt

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

load d:\exercise\vow\train.txt;
load d:\exercise\vow\test.txt;
a0=train(:,1:10);
b0=test(:,1:10);
a=0.75*a0;
for i=1:11,
a1(:,i:i)=diag((a-ones(528,1)*mean(a((i-1)*48+1:i*48,:)))*inv(cov(a((i-1)*48+1:i*48,:)))*(a-ones(528,1)*mean(a((i-1)*48+1:i*48,:)))');
end
a2=exp(-a1/2);
for i=1:11;
a2(:,i:i)=a2(:,i:i)/(2*pi*2*pi*2*pi*2*pi*2*pi*sqrt(det(cov(a((i-1)*48+1:i*48,:)))));
end
a2(:,12:12)=train(:,22:22);
fid=fopen('d:\exercise\vow\resultr.txt','w');
fprintf(fid,'%-8.4f %-8.4f %-8.4f %-8.4f %-8.4f %-8.4f %-8.4f %-8.4f %-8.4f %-8.4f %-8.4f %g\n',a2');
b=0.75*b0;
for i=1:11,
b1(:,i:i)=diag((b-ones(462,1)*mean(a((i-1)*48+1:i*48,:)))*inv(cov(a((i-1)*48+1:i*48,:)))*(b-ones(462,1)*mean(a((i-1)*48+1:i*48,:)))');
end
b2=exp(-b1/2);
for i=1:11;
b2(:,i:i)=b2(:,i:i)/(2*pi*2*pi*2*pi*2*pi*2*pi*sqrt(det(cov(a((i-1)*48+1:i*48,:)))));
end
b2(:,12:12)=test(:,22:22);
fid=fopen('d:\exercise\vow\resulte.txt','w');
fprintf(fid,'%-8.4f %-8.4f %-8.4f %-8.4f %-8.4f %-8.4f %-8.4f %-8.4f %-8.4f %-8.4f %-8.4f %g\n',b2');