实验报告：混沌同步控制与图像加密(3)

tspan = tstart:0.5*tstep:(tstart + tstep*steps); [t,Y1] = ode45(@Lorenz_driven,tspan,y); y = Y1(3,:); y(1)=S(i+1,2);

tstart = tstart + tstep*steps; results(i,:)=Y1(3,:); end

subplot(3,1,2);

plot(results(:,1),results(:,3),'b-'); title('响应'); subplot(3,1,3);

plot(Y(:,1),Y(:,3),'k-',results(:,1),results(:,3),'b-'); title('驱动-响应叠加 '); for i=1:iteratetimes TT(1,i)=i*tstep; end figure(2); subplot(2,1,1);

plot(T,Y(:,3),'k-');title('驱动'); subplot(2,1,2);

plot(TT,results(:,3),'k-');title('响应');

3.驱动系统的输出： function s=output;

global Y; global T;

tstart = 0; tstep = 0.1; wholetimes = 200; tspan=tstart:tstep:wholetimes*tstep;

[T,Y] = ode45(@Lorenz_driven,tspan,[0.1 0.1 0.1]); s=[T,Y]; figure(1); subplot(3,1,1);

plot(Y(:,1),Y(:,3),'k-'); title('驱动');

其中，我们设置驱动系统的初值为[0.1,0.1,0.1],而响应系统的初值为[0.1,0.1,20],显然Z的初值相差很大。在MATLAB界面窗口输入指令：Lorenz_response ，可以看到以下X-Z相图图像：