%% 清空環(huán)境變量

clc

clear

%% 數(shù)據(jù)導入

load data1 array_one

load data2 array_two

%% 訓練樣本（目標向量）

 T=[array_one;array_two]';

%% 創(chuàng)建網(wǎng)絡

 net=newhop(T);

%% 數(shù)字1和2的帶噪聲數(shù)字點陣（隨機法）

noisy_array_one=array_one;

noisy_array_two=array_two;

for i=1:100

    a=rand;

    if a<0.1

       noisy_array_one(i)=-array_one(i);

       noisy_array_two(i)=-array_two(i);

    end

end

%% 數(shù)字識別

% identify_one=sim(net,10,[],noisy_array_one');

noisy_one={(noisy_array_one)'};

identify_one=sim(net,{10,10},{},noisy_one);

identify_one{10}';

noisy_two={(noisy_array_two)'};

identify_two=sim(net,{10,10},{},noisy_two);

identify_two{10}';

%% 結果顯示

Array_one=imresize(array_one,20);

subplot(3,2,1)

imshow(Array_one)

title('標準(數(shù)字1)') 

Array_two=imresize(array_two,20);

subplot(3,2,2)

imshow(Array_two)

title('標準(數(shù)字2)') 

subplot(3,2,3)

Noisy_array_one=imresize(noisy_array_one,20);

imshow(Noisy_array_one)

title('噪聲(數(shù)字1)') 

subplot(3,2,4)

Noisy_array_two=imresize(noisy_array_two,20);

imshow(Noisy_array_two)

title('噪聲(數(shù)字2)')

subplot(3,2,5)

imshow(imresize(identify_one{10}',20))

title('識別(數(shù)字1)')

subplot(3,2,6)

imshow(imresize(identify_two{10}',20))

title('識別(數(shù)字2)')