function xdot=lorenzeq(t, x)
xdot=[-16*x(1)+16*x(2);
-x(1)*x(3)+45.92*x(1)-x(2);
-4*x(3)+x(1)*x(2);
-16*x(4)+16*x(5);
-x(1)*x(6)+45.92*x(1)-x(5);
-4*x(6)+x(1)*x(5)];
function ydot=lorenzeq(t, y)
ydot=[-16*y(1)+16*y(2);
-y(1)*y(3)+45.92*y(1)-y(2);
-4*y(3)+y(1)*y(2);
-16*y(4)+16*y(2);
-y(4)*y(6)+45.92*y(4)-y(5);
-4*y(6)+y(4)*y(2)];
主程序:
t_final=50;
%x 做驱动时的 PC 混沌同步
x0=[-10;-10;15;0;-25;-6];
[t, x]=ode45('lorenzeq',[0,t_final],x0);
plot(t, x(:,1),t, x(:,2),t, x(:,3));
figure;
plot3(x(:,1),x(:,2),x(:,3));
axis([-5 20 -5 20 -10 30]);
figure;
plot(t, x(:,2),'-',t, x(:,5),'--');
xlabel('t');
legend('—y','- -y’')
axis([0 3 -80 60]);
figure;
plot(t, x(:,3),'-',t, x(:,6),'--');
xlabel('t');
legend('—z','- -z’')
axis([0 3 -30 100]);
figure;
%y 做驱动时的 PC 混沌同步
y0=[-10;-8;16;0;12;-4];
[t, y]=ode45('lorenzeq1',[0,t_final],y0);
plot(t, y(:,1),t, y(:,2),t, y(:,3));
figure;
plot3(y(:,1),y(:,2),y(:,3));
axis([-5 20 -5 20 -10 30]);
figure;
plot(t, y (:,1),'-',t, y(:,4),'--');
xlabel('t');
评论1