aihara.m

这个是经典Aihara混沌神经网络模型的程序

clc
clear
k1=0.89;
alpha=0.017;
gama=0.007;
T=0.6;
u=250;
L=0.07;
A=1;B=1;
I=ones(10,10);
ee=2.71828;
xx(1)=0.4000; yy(1)=0.4493;
xx(2)=0.2493; yy(2)=0.1463;
xx(3)=0.1707; yy(3)=0.2293;
xx(4)=0.2293; yy(4)=0.7610;
xx(5)=0.5171; yy(5)=0.9414;
xx(6)=0.8732; yy(6)=0.6536;
xx(7)=0.6878; yy(7)=0.5219;
xx(8)=0.8488; yy(8)=0.3609;
xx(9)=0.6683; yy(9)=0.2536;
xx(10)=0.6195; yy(10)=0.2634;

n=10;

%conclude distance%
for i=1:n
for j=1:n
dista(i,j)=distance(xx(i),yy(i),xx(j),yy(j));
if i==j
dista(i,j)=0;
end
end
end
% chaos progress % 
e1=0;e11=0;e22=0;e2=0;e4=0;
y1=rand(10,10)*2-1;
f=0;
m=0;
s=0;
t=0;
for w=1:1000
v=I./(1+exp(-(u*y1)));

for x=1:n
for i=1:n
e1=e1+v(x,i); 
end
e11=e11+(e1-1)^2;
end

for i=1:n
for x=1:n
e2=e2+v(x,i); 
end
e22=e22+(e2-1)^2;
end

for x=1:n 
for y=x+1:n
for i=1:n
if i==1

e4=e4+v(x,i)*dista(x,y)*(v(y,2)+v(y,n));
elseif i==n
e4=e4+v(x,i)*dista(x,y)*(v(y,n-1)+v(y,1));
else
e4=e4+v(x,i)*dista(x,y)*(v(y,i-1)+v(y,i+1));
end
end
end
end
e=A/2*(e1+e2)+B/2*e4;


for i=1:n
for j=1:n 
v1(i)=sum(v(i,:))-1;
end
end
v1=sum(v1);
for j=1:n
for i=1:n
v2(j)=sum(v(:,j))-1;
end
end

v2=sum(v2);

for k=1:n 
for y=k+1:n 
for i=1:n
if i==1
v31(i)=v(y,2)+v(y,10);
elseif i==n
v31(i)=v(y,1)+v(y,9);
else
v31(i)=v(y,i+1)+v(y,i-1);

end
v32(y)=sum(v31);
end
v33(k)=dista(k,y)*sum(v32);
if k==10
v33(10)=dista(10,1)*sum(v32);
end
end 
end
v3=sum(v33);
DE=-A*(v1+v2-3)-B*v3;
y2=k1*y1+alpha*DE-L*(v-T);
L=L/(log2(2+gama*(1-L)));
y1=y2;
end