%Examp1e52 Sim
clear all
clc;
p(1:256,1)=1;
p1=ones(16,16);%初始化16.16的二值圖像像素值（全白)
load E52net net;%加載訓練后的BP網絡
test=input('p1ease input s test image;','s');%提示輸入測試樣本圖像文件名
x=imread(test,'bmp');%讀入測試樣本圖像
bw=im2bw(x,0.5);%將讀入的訓練樣本圖像轉換為二值圖像
[i,j]=find(bw==0);%尋找二值圖像中像素值為。的行號和列號
imin=min(i);%尋找二值圖像中像素值為。的最小行號
imax=max(i);%尋找二值圖像中像素值為0的最大行號
jmin=min(j);%尋找二值圖像中像素值為。的最小列號
jmax=max(j);%尋找二值圖像中像素值為0的最大列號
bw1=bw(imin:imax,jmin:jmax);%截取圖像像素值為0的最大矩形區域
rate=16/max(size(bw1));%計算轉換為16.16的二值圖像的縮放比例
bw1=imresize(bw1,rate);%將截取圖像轉換為16 x 16的二值圖像
[i,j]=size(bw1);%轉換圖像的大小
i1=round((16-i)/2);%計算轉換圖像的寬度與16的差距
j1=round((16-j)/2);%計算轉換圖像的高度與16的差距
p1(i1+1:i1+i,j1+1:j1+j)=bw1;%將截取圖像轉換為標準的16.16的圖像
p1=-1.*p1+ones(16,16);%反色處理
 for m =0:15
        p(m*16+1:(m+1)*16,1)=p1(1:16,m+1);%形成鍘試樣本輸入向量
 end
    [a,Pf,Af]=sim(net,p);%網絡仿真
    imshow(p1);%顯示鍘試樣本圖像
    a=round(a)%輸出識別結果