%例5.3, hop3.m
%
clear
T=[1 -1;-1 1];
net=newhop(T);                          %創建Hopfield網絡
w=net.lw{1,1},b=net.b{1}                 %輸出權值和偏差
Ai ={T};                              
[Y,Pf,Af] = sim(net,{2,3},[],Ai);             %仿真輸出
A1=Y{1},A2=Y{2},A3=Y{3}               %給出三次循環的結果

P=[0.5621 0.3577 0.8694 0.0388 -0.9309 0.0594;
   -0.9095 0.3586 -0.2330 0.6619 0.8931 0.3423];
Ptest={P};                               %輸入檢測樣本
[Z,Pf,Af] = sim(net,{6,25},[],Ptest);           %仿真輸出
Ap=Z{25}                               %給出網絡循環的最后一次結果
plot(T(1,:),T(2,:),'r*')                       %作目標點
hold on
plot(P(1,:),P(2,:),'r+')                       %作起始點
for i=1:6
	D{i}=P(:,i) ;                          %初始化D{i}，共6個矩陣
end
for i=1:25
    plot(Z{i}(1,:),Z{i}(2,:),'ro')              %作訓練過程中點的變化
    for j=1:6
        D{j}=[D{j} Z{i}(:,j)];              %把Z中60個矩陣中的同一位置的元素賦給
%同一矩陣   
    end
end
for j=1:6
    plot(D{j}(1,:),D{j}(2,:))                 %畫每個點的變化直線圖
end
hold off