%指定存儲在網絡中的目標平衡點
T = [+1 +1;-1 +1; -1 -1]
axis([-1 1 -1 1 -1 1])';
set(gca,'box','on');
axis manual;
hold on;
plot3(T(1,:),T(2,:),T(3,:),'r*');
title('Hopfield神經網絡狀態空間');
xlabel('a(1)');
ylabel('a(2)');
zlabel('a(3)');
view([37.5 30]);

%建立網絡，并得到權重和閾值
net=newhop(T);
W=net.LW{1,1}
b=net.b{1,1}

%使用一個隨機點仿真網絡，并繪出其到達穩定點的軌跡
a = {rands(3,1)};
[y,Pf,Af] = sim(net,{1 10},{},a);
figure;
axis([-1 1 -1 1 -1 1])';
set(gca,'box','on');
axis manual;
hold on;
plot3(T(1,:),T(2,:),T(3,:),'r*');
record=[cell2mat(a) cell2mat(y)];
start=cell2mat(a);
plot3(start(1,1),start(2,1),start(3,1),'bx',record(1,:),record(2,:),record(3,:));
view([37.5 30]);
title('Hopfield神經網絡狀態空間');
xlabel('a(1)');
ylabel('a(2)');
zlabel('a(3)');

%用多個隨機點仿真網絡，并繪出相應的軌跡
figure;
axis([-1 1 -1 1 -1 1])';
set(gca,'box','on');
axis manual;
hold on;
view([37.5 30]);
plot3(T(1,:),T(2,:),T(3,:),'r*');
color = 'rgbmy';
for i=1:25
    a = {rands(3,1)};
    [y,Pf,Af] = sim(net,{1 10},{},a);
    record=[cell2mat(a) cell2mat(y)];
    start=cell2mat(a);
    plot3(start(1,1),start(2,1),start(3,1),'kx',record(1,:),record(2,:),record(3,:),color(rem(i,5)+1))
    drawnow
end
title('Hopfield神經網絡狀態空間');
xlabel('a(1)');
ylabel('a(2)');
zlabel('a(3)');

%指定一些特殊的初始點進行仿真，并繪出相應的軌跡
figure;
axis([-1 1 -1 1 -1 1])';
set(gca,'box','on');
axis manual;
hold on;
P = [ 1.0  -1.0  -0.5  1.00  1.00  0.0; 0.0   0.0   0.0 
 0.00  0.00 -0.0;-1.0   1.0   0.5 -1.01 -1.00  0.0];
view([37.5 30]);
color ='rgbmy';
for i=1:6
    a = {P(:,i)}; 
    [y,Pf,Af] = sim(net,{1 10},{},a);
    record=[cell2mat(a) cell2mat(y)];
    start=cell2mat(a);
    plot3(start(1,1),start(2,1),start(3,1),'kx',record(1,:),record(2,:),record(3,:),color(rem(i,5)+1))
    drawnow
end
title('Hopfield神經網絡狀態空間');
xlabel('a(1)');
ylabel('a(2)');
zlabel('a(3)');