使用MATLAB分别画出以下混沌系统的吸引盆和动力学地图并给出代码 function dy=ZY1(t,y) dy=zeros(5,1); global k1 k2 k3 a1 a2 b1 b2 c1 c2 w11 w22 W1 W2 k1=3;a1=2.8;b1=1.5;c1=5;k2=3.2;a2=-2.8;b2=1.5;c2=5;k3=0.01;w11=3;w22=1.4; W1=a1-b1*tanh(y(3))+c1*tanh(y(4));W2=a2-b2*tanh(y(4))+c2*tanh(y(3)); dy(1)=-y(1)+w11*tanh(y(1))+k2*W2*tanh(y(2))+k3*(y(1)-y(2))*tanh(y(5)); dy(2)=-y(2)+k1*W1*tanh(y(1))+w22*tanh(y(2))-k3*(y(1)-y(2))*tanh(y(5)); dy(3)=-y(3)+tanh(y(1)); dy(4)=-y(4)+tanh(y(2)); dy(5)=y(1)-y(2)-y(5);

好的，以下是MATLAB代码：

% 定义参数
global k1 k2 k3 a1 a2 b1 b2 c1 c2 w11 w22 W1 W2
k1=3;a1=2.8;b1=1.5;c1=5;
k2=3.2;a2=-2.8;b2=1.5;c2=5;
k3=0.01;w11=3;w22=1.4;

% 定义ODE
fun = @ZY1;
tspan = [0 100];
y0 = [0.1 0.1 0.1 0.1 0.1];
[t,y] = ode45(fun,tspan,y0);

% 画吸引盆
figure(1)
plot3(y(:,3),y(:,4),y(:,5))
xlabel('y_3')
ylabel('y_4')
zlabel('y_5')
title('吸引盆')

% 画动力学地图
figure(2)
plot3(y(:,1),y(:,2),y(:,5))
xlabel('y_1')
ylabel('y_2')
zlabel('y_5')
title('动力学地图')

% 定义ODE函数
function dy=ZY1(t,y)
dy=zeros(5,1);
global k1 k2 k3 a1 a2 b1 b2 c1 c2 w11 w22 W1 W2
W1=a1-b1*tanh(y(3))+c1*tanh(y(4));
W2=a2-b2*tanh(y(4))+c2*tanh(y(3));
dy(1)=-y(1)+w11*tanh(y(1))+k2*W2*tanh(y(2))+k3*(y(1)-y(2))*tanh(y(5));
dy(2)=-y(2)+k1*W1*tanh(y(1))+w22*tanh(y(2))-k3*(y(1)-y(2))*tanh(y(5));
dy(3)=-y(3)+tanh(y(1));
dy(4)=-y(4)+tanh(y(2));
dy(5)=y(1)-y(2)-y(5);
end