%例5.4, hop4.m
%
clear
T=[1 -1;-1 1];
net=newhop(T);                      %創建Hopfield網絡
w=net.lw{1,1},b=net.b{1}              %輸出權值和偏差
plot(T(1,:),T(2,:),'r*')                  %作目標節點圖
P=[-1 -0.5 0 0.5 1.0;
    -1 -0.5 0 0.5 1.0];
hold on
plot(P(1,:),P(2,:),'r+')                  %作輸入節點圖 
Ptest={P};
[Z,Pf,Af] = sim(net,{5,60},[],Ptest);      %仿真輸出
for i=1:5
	D{i}=P(:,i) ;                     %初始化D{i}，共6個矩陣
end
for i=1:60
    plot(Z{i}(1,:),Z{i}(2,:),'ro')         %作訓練過程中點的變化
    for j=1:5
        D{j}=[D{j} Z{i}(:,j)];         %把Z中60個矩陣中的同一位置的元素賦給同一矩陣   
    end
end
for j=1:5
    plot(D{j}(1,:),D{j}(2,:))            %畫每個點的變化直線圖
end
hold off