基于pca数字识别，成功率达百分之九十五以上，内带有测试、训练集，可直接运行.zip

global im; imgdata=[]; for i=0:9 for j=51:200 a=imread(strcat('D:\程序\数字识别\num\训练\',num2str(i),'_',num2str(j),'.bmp')); b=a(1:28*28); b=double(b); imgdata=[imgdata;b]; end; end; imgdata=imgdata'; imgmean=mean(imgdata,2); %每列的平均值 for i=1:1500 minus(:,i)=imgdata(:,i)-imgmean; %去均值操作 end; covx=minus'*minus; %得到协方差矩阵 [COEFF,latent,explained]=pcacov(covx'); i=1; proportion=0; while(proportion<95) proportion=proportion+explained(i);%特征向量所占比例累加，达到95%停止 i=i+1; end; p=i-1; i=1; while(i<=p&&latent(i)>0) base(:,i)=latent(i)^(-1/2)*minus*COEFF(:,i); i=i+1; end; reference=base'*minus; %[f,p]=uigetfile('*.bmp'); %filename=fullfile(p,f); %a=imread(filename); n=0; for i=0:9 for j=201:240 a=imread(strcat('D:\程序\数字识别\num\测试\0\',num2str(i),'_',num2str(j),'.bmp')); b=a(1:28*28); b=double(b); b=b'; object=base'*(b-imgmean); distance=norm(object-reference(:,1)); for num1=1:10 for num2=1:150 temp=norm(object-reference(:,(num1-1)*150+num2)); if(distance>temp) k=num1; distance=temp; end; end; end; if((k-1)==i) n=n+1; end; end; end; fprintf('正确率为%f\n',n/400);