for dd=1:50
dt=0.01;
N=301.6/dt;
dxdt=zeros(N,2);
dydt=zeros(N,2);
x=zeros(N,2);
y=zeros(N,2);
a=[randn(1) randn(1)];
for i=1:160
x(i,1)=a(1);
x(i,2)=a(2);
y(i,1)=a(1);
y(i,2)=a(2);
end
tau=0.5;
n=normrnd(0,1,2,N-160);
for i=160:N-1
if abs(x(i,1))<1
a11=1.715;
a12=1.96;
b11=-1.768;
b12=-0.104;
d11=2;
else
a11=1.715;
a12=1.96;
b11=-1.768;
b12=-0.104;
d11=2;
end
if abs(x(i,2))<1
a21=1.8;
a22=-0.098;
b21=-1.8;
b22=-1.56;
d22=2;
else
a21=-1.568;
a22=-0.098;
b21=-1.8;
b22=-1.56;
d22=2;
end
if abs(y(i,1))<1
p11=1.8375;
p12=-0.105;
q11=-1.615;
q12=-0.095;
r11=2;
else
p11=1.8375;
p12=-0.105;
q11=-1.615;
q12=-0.095;
r11=2;
end
if abs(y(i,2))<1
p21=-1.68;
p22=2.1;
q21=0.095;
q22=-1.425;
r22=2;
else
p21=-1.68;
p22=2.1;
q21=0.095;
q22=-1.425;
r22=2;
end
ct=(i-160)*dt;
dxdt(i,1)=-d11*x(i,1)+a11*tanh(y(i,1))+a12*tanh(y(i,2))+b11*tanh(y(i-fix((tau/2+0.5*tau*sin(ct))/dt),1))+b12*tanh(y(i-fix((tau/2+0.5*tau*sin(ct))/dt),2));
dxdt(i,2)=-d22*x(i,2)+a21*tanh(y(i,1))+a22*tanh(y(i,2))+b21*tanh(y(i-fix((tau/2+0.5*tau*sin(ct))/dt),1))+b22*tanh(y(i-fix((tau/2+0.5*tau*sin(ct))/dt),2));
dydt(i,1)=-r11*y(i,1)+p11*tanh(x(i,1))+p12*tanh(x(i,2))+q11*tanh(x(i-fix((tau/2+0.5*tau*sin(ct))/dt),1))+q12*tanh(x(i-fix((tau/2+0.5*tau*sin(ct))/dt),2));
dydt(i,2)=-r22*y(i,2)+p21*tanh(x(i,1))+p22*tanh(x(i,2))+q21*tanh(x(i-fix((tau/2+0.5*tau*sin(ct))/dt),1))+q22*tanh(x(i-fix((tau/2+0.5*tau*sin(ct))/dt),2));
x(i+1,1)=x(i,1)+dxdt(i,1)*dt;
x(i+1,2)=x(i,2)+dxdt(i,2)*dt;
y(i+1,1)=y(i,1)+dydt(i,1)*dt;
y(i+1,2)=y(i,2)+dydt(i,2)*dt;
end
x1(:,dd)=x(:,1);
x2(:,dd)=x(:,2);
y1(:,dd)=y(:,1);
y2(:,dd)=y(:,2);
figure(1);
xlabel('x1','fontsize',12,'fontweight','b')
ylabel('x2','fontsize',12,'fontweight','b')
plot(x1(101:30160,dd),x2(101:30160,dd),'b');
hold on
figure(2);
xlabel('y1','fontsize',12,'fontweight','b')
ylabel('y2','fontsize',12,'fontweight','b')
plot(y1(101:6160,dd),y2(101:6160,dd),'b');
hold on
figure(3);
plot([0:0.01:60.59]',y1(101:6160,dd),'b');
hold on
end