global a b k d
h=0.005;
a=20;
b=15;
c=32;
k=1;
d=1;
N=6000;
q=0.9;
x0=[0.1 0.2 2];
[t,y]=Fratalsim(h,N,x0,q,c);
figure(1)
plot(y(1,3000:end),y(3,3000:end));
%xlim([-12 12]);
xlabel('\itx');
ylabel('\itz');
figure(2)
plot(y(2,:),y(3,:));
%xlim([-12 12]);
xlabel('\ity');
ylabel('\itz');
figure(3)
plot(y(1,:),y(2,:));
%xlim([-12 12]);
xlabel('\itx');
ylabel('\ity');
figure(4)
plot3(y(3,:),y(2,:),y(1,:));
%xlim([-12 12]);
xlabel('\itz');
ylabel('\ity');
zlabel('\itx');
评论5