nn.m

来自「人工神经网络」· M 代码 · 共 14 行

%x1=-x1+pf(x1)-sf(x2)+asin(2piT/t);x2=-x2+sf(x1)+pf(x2)
clear;clc;
x1=0.14;x2=-0.1;
p=1.2;s=2;a=1.5;T=4;t=0.0001;
for i=1:12000
    y1=0.5*(abs(x1+1)-abs(x1-1));
    y2=0.5*(abs(x2+1)-abs(x2-1));
    x11=x1+(-x1+p*y1-s*y2+a*sin(2*pi*T/t))*0.01;
    x22=x2+(-x2+s*y1+p*y2)*0.01;
    X(i)=x11;Y(i)=x22;
    x1=x11;x2=x22;
    t=t+0.001;
end
plot(X,Y);