%例5.5, hop5.m
%
T=[1 1;-1 1;-1 -1];
net=newhop(T);                        %創建Hopfield網絡
w=net.lw{1,1},b=net.b{1}                %輸出權值和偏差
P=[1 -1 -0.5 1 1 0;
   0 0 0 0 0.01 -0.2;
   -1 1 0.5 -1.01 -1 0];
Ptest={P};
[Z,Pf,Af] = sim(net,{6,60},[],Ptest);         %仿真輸出
z5=Z{5},z50=Z{50},                     %5和50次仿真后的輸出結果
plot3(T(1,:),T(2,:),T(3,:),'r*')               %作目標節點圖
hold on
plot3(P(1,:),P(2,:),P(3,:),'r+')               %作輸入節點圖 
axis square; 
grid on
xlabel('Z(1)');                           %坐標軸標題
ylabel('Z(2)');
zlabel('Z(3)');
for i=1:6
	D{i}=P(:,i) ;                         %初始化D{i}，共6個矩陣
end
for i=1:60
    plot3(Z{i}(1,:),Z{i}(2,:),Z{i}(3,:),'ro')   
                                       %作訓練過程中點的變化
    for j=1:6
        D{j}=[D{j} Z{i}(:,j)];             %把Z中60個矩陣中的同一位置的元素賦給同
%一矩陣   
    end
end
for j=1:6
    plot3(D{j}(1,:),D{j}(2,:),D{j}(3,:))  
                                       %畫每個點的變化直線圖
end
hold off