模式识别课matlab数字识别程序

1、名 称：模式识别题 目 数字 3 '和 4的识别实验目的与要求：利用的数字样本 3 和 4，提取样本特征，并确定分类准那么， 在用测试样本对分类确 定准那么的错误率进行分析。 进一步加深对模式识别方法的理解， 强化利用电脑实现模式识别。 实验原理：1. 特征提取原理 :利用 MATLAN 软件把图片变为一个二维矩阵，然后对该矩阵进行二值化处理。由于“3的下半局部在横轴上的投影比“ 4的下半局部在横轴上的投影宽，所以可以统计3'4'在横轴上投影的 1'的个数作为一个特征。又由于 4'中间纵向比 3'的中间 1'的个 数多，所以可以统计 4&

2、#39;和 3'中间区域 1'的个数作为另外一个特征，又考虑4'的纵向可能会有点偏，所以在统计一的个数的时候，取的范围稍微大点，但不能太大。2. 分类准那么原理：利用最近邻对测试样本进行分类实验步骤1.利用 MATLAN 软件把前 30 个图片变为一个二维矩阵，然后对该矩阵进行二值化处理。 2.利用上述矩阵生成特征向量3. 利用 MATLAN 软件把后 5 个图片变为一个二维矩阵，然后对该矩阵进行二值化处理。4. 对测试样本进行分类，用 F 矩阵表示结果，如果是 1'表示分类正确， 0'表示分类错 误。5. 对分类错误率分析实验原始程序：f=zeros(

3、5,2)w=zeros(35,2)q=zeros(35,2)for i=1:35filename_1='D:MATLAB6p5toolboximagesimdemos3'filename_2='.bmp'a= num2str (i)b=strcat(filename_1,a)c=strcat(b,filename_2)d=imread(c)e=im2bw(d)n=0for u=1:20m=0for t=32:36if(e(t,u)=0)m=m+1endendif(m<5)n=n+1endendw(i,1)=nn=0for u=1:36for t=10:18

4、 n=n+e(u,t) end end w(i,2)=n filename_1='D:MATLAB6p5toolboximagesimdemos4' filename_2='.bmp'a= num2str(i) b=strcat(filename_1,a) c=strcat(b,filename_2) d=imread(c) e=im2bw(d)n=0for u=1:20m=0for t=32:36 if(e(t,u)=0) m=m+1 end end if(m<5) n=n+1 end end q(i,1)=n n=0for u=1:36for t=10

5、:18 n=n+e(u,t) end end q(i,2)=n endz=zeros(5,2) x=zeros(5,2)for i=1:5 filename_1='D:MATLAB6p5toolboximagesimdemos3' filename_2='.bmp'a= num2str (i+35)b=strcat(filename_1,a) c=strcat(b,filename_2)d=imread(c) e=im2bw(d) n=0 for u=1:20 m=0 for t=32:36 if(e(t,u)=0) m=m+1 end end if(m<

6、5) n=n+1 end end z(i,1)=n n=0 for u=1:36 for t=10:18 n=n+e(u,t) end end z(i,2)=n filename_1='D:MATLAB6p5toolboximagesimdemos4' filename_2='.bmp' a= num2str (i) b=strcat(filename_1,a) c=strcat(b,filename_2) d=imread(c) e=im2bw(d) n=0 for u=1:20 m=0 for t=32:36 if(e(t,u)=0) m=m+1 end e

7、nd if(m<5) n=n+1 end end x(i,1)=n n=0 for u=1:36for t=10:18n=n+e(u,t)endendx(i,2)=nendd仁 10000*o nes(5,2)d2=10000*o nes(5,2)for i=1:5for j=1:35a=(w(j,1)-z(i,1)*(w(j,1)-z(i,1)+(w(j,2)-z(i,2)*(w(j,2)-z(i,2) b=(q(j,1)-z(i,1)*(q(j,1)-z(i,1)+(q(j,2)-z(i,2)*(q(j,2)-z(i,2) if(a<d1(i,1)d1(i,1)=aendif(b<d1(i,2)d1(i,2)=bendendif(d1(i,1)<d1(i,2)f(i,1)=1endfor j=1:35c=(w(j,1)-x(i,1)*(w(j,1)-x(i,1)+(w(j,2)-x(i,2)*(w(j,2)-x(i,2) d=(q(j,1)-x(i,1)*(q(j,1)-x(i,1)+(q(j,2)-x(i,2)*(q(j,2)-x(i,2) if(d2(i,1)>c)d2(i,1)=cendif(d2(i,2)>d)d2(i,2)=dendendif(d2(