main100b2.m

用9*9的矩阵显示字符0~9。 对BP网络进行训练后

clear
sampin=zeros(10,10,10);
sampin(:,:,1)=  [0,0,0,1,1,1,1,0,0,0
                        0,0,1,1,1,1,1,1,0,0
                        0,1,1,0,0,0,0,1,1,0
                        1,1,0,0,0,0,0,0,1,1
                        1,1,0,0,0,0,0,0,1,1
                        1,1,0,0,0,0,0,0,1,1
                        1,1,0,0,0,0,0,0,1,1
                        0,1,1,0,0,0,0,1,1,0
                        0,0,1,1,1,1,1,1,0,0
                        0,0,0,1,1,1,1,0,0,0];       
                    
sampin(:,:,2)= [ 0,0,0,0,1,1,0,0,0,0
                        0,0,0,1,1,1,0,0,0,0
                        0,0,0,1,1,1,0,0,0,0
                        0,0,0,0,1,1,0,0,0,0
                        0,0,0,0,1,1,0,0,0,0
                        0,0,0,0,1,1,0,0,0,0
                        0,0,0,0,1,1,0,0,0,0
                        0,0,0,0,1,1,0,0,0,0
                        0,0,0,1,1,1,1,0,0,0
                        0,0,1,1,1,1,1,1,0,0];      
                    
sampin(:,:,3)=[  0,0,0,1,1,1,1,0,0,0
                        0,0,1,1,1,1,1,1,0,0
                        0,1,1,0,0,0,0,1,1,0
                        1,1,0,0,0,0,0,0,1,1
                        0,0,0,0,0,0,0,1,1,0
                        0,0,0,0,0,0,1,1,0,0
                        0,0,0,0,0,1,1,0,0,0
                        0,0,0,0,1,1,0,0,0,0
                        0,0,1,1,1,1,1,1,1,0
                        1,1,1,1,1,1,1,1,1,1];      
                    
sampin(:,:,4)= [ 0,0,0,1,1,1,1,0,0,0
                        0,0,1,1,0,0,1,1,0,0
                        0,1,1,0,0,0,0,1,1,0
                        0,0,0,0,0,0,1,1,0,0
                        0,0,0,1,1,1,1,0,0,0
                        0,0,0,1,1,1,1,0,0,0
                        0,0,0,0,0,0,1,1,0,0
                        0,1,1,0,0,0,0,1,1,0
                        0,0,1,1,0,0,1,1,0,0
                        0,0,0,1,1,1,1,0,0,0];      
                    
sampin(:,:,5)= [ 0,0,0,0,0,1,1,0,0,0
                        0,0,0,0,1,1,1,0,0,0
                        0,0,0,1,1,1,1,0,0,0
                        0,0,1,1,0,1,1,0,0,0
                        0,1,1,0,0,1,1,0,0,0
                        1,1,1,1,1,1,1,1,1,1
                        1,1,1,1,1,1,1,1,1,1
                        0,0,0,0,0,1,1,0,0,0
                        0,0,0,0,0,1,1,0,0,0
                        0,0,0,0,0,1,1,0,0,0]; 
                    
sampin(:,:,6)=[  1,1,1,1,1,1,1,1,1,0
                        0,1,1,1,1,1,1,1,1,0
                        0,0,1,1,0,0,0,0,0,0
                        0,1,1,1,1,1,0,0,0,0
                        1,1,1,1,1,1,1,0,0,0
                        0,0,0,0,0,0,1,1,0,0
                        0,0,0,0,0,0,0,1,1,1
                        1,1,0,0,0,0,0,0,1,1
                        0,1,1,0,0,0,0,1,1,0
                        0,0,1,1,1,1,1,1,0,0];
                    
sampin(:,:,7)=[  0,0,0,1,1,1,1,0,0,0
                        0,0,1,1,0,0,1,1,0,0
                        0,1,1,0,0,0,0,1,1,0
                        1,1,0,0,0,0,0,0,0,0
                        1,1,0,1,1,1,1,0,0,0
                        1,1,1,1,0,0,1,1,0,0
                        1,1,0,0,0,0,0,0,1,1
                        0,1,1,0,0,0,0,1,1,0
                        0,0,1,1,0,0,1,1,0,0
                        0,0,0,1,1,1,1,0,0,0];
                    
sampin(:,:,8)= [ 1,1,1,1,1,1,1,1,1,1
                        1,1,1,1,1,1,1,1,1,1
                        0,0,0,0,0,0,0,0,1,1
                        0,0,0,0,0,0,1,1,0,0
                        0,0,0,0,0,1,1,0,0,0
                        0,0,0,0,1,1,0,0,0,0
                        0,0,0,1,1,0,0,0,0,0
                        0,0,1,1,0,0,0,0,0,0
                        0,1,1,0,0,0,0,0,0,0
                        0,1,1,0,0,0,0,0,0,0];
                    
sampin(:,:,9)= [ 0,0,0,1,1,1,1,0,0,0
                        0,0,1,1,1,1,1,1,0,0
                        0,1,1,0,0,0,0,1,1,0
                        0,0,1,1,0,0,1,1,0,0
                        0,0,0,1,1,1,1,0,0,0
                        0,0,0,1,1,1,1,0,0,0
                        0,0,1,1,0,0,1,1,0,0
                        0,1,1,0,0,0,0,1,1,0
                        0,0,1,1,0,0,1,1,0,0
                        0,0,0,1,1,1,1,0,0,0];
                    
sampin(:,:,10)= [0,0,0,1,1,1,1,0,0,0
                         0,0,1,1,0,0,1,1,0,0
                         0,1,1,0,0,0,0,1,1,0
                         0,0,1,1,0,0,1,1,0,0
                         0,0,0,1,1,1,1,0,0,0
                         0,0,0,0,1,1,0,0,0,0
                         0,0,0,1,1,0,0,0,0,0
                         0,0,1,1,0,0,0,0,0,0
                         0,1,1,0,0,0,0,0,0,0
                         1,1,1,0,0,0,0,0,0,0];      
         
sampout=eye(10);
w12=rand(9,100).*0.2-0.1;
w23=rand(10,9).*0.2-0.1;
ep=0.1;
e=ones(10,1);
qn=1;
N=5000;
para=zeros(19,109);
para(1:9,1:100)=w12;
para(10:19,101:109)=w23;
n=1;
e=GetErr100b2(sampin,sampout,para);
[emx(n),i]=max(e);
while(emx>ep&n<N)
    input=sampin(:,:,i);
    output=sampout(:,i);
    para=training100b2(para,output,input,qn,ep);
    e=GetErr100b2(sampin,sampout,para);
    n=n+1;
    [emx(n),i]=max(e);
   
end
input0  =  sampin(:,:,6)
input   =  sampin(:,:,6)+rand(10,10)*0.1

O2=zeros(9,1);
O3=zeros(10,1);
w12=para(1:9,1:100);
w23=para(10:19,101:109);
In=reshape(input,100,1);

net2=w12*In;
O2=sigmf(net2,[1 0]);
net3=w23*O2;
O3=sigmf(net3,[1,0]);

for i=1:10
    if O3(i)>0.5
        output(i)=1;
    else
        output(i)=0;
    end
end

output'
plot(emx)