怎样用matlab作出混沌

matlab命令如下：

global x

global y

global k

y=1

x=1

p=plot(x,y,'.','EraseMode','none','MarkerSize',3)

axis([0 2 -2 2])

hold on

for x=1:200

for k=1:500

y=1-x*y*y/100set(p,'Xdata',x/100,'Ydata',y)

drawnow

set(p,'Xdata',x/100,'Ydata',y)

drawnow

end

end

你找到某种混沌系统的表达式，作图就行了，比如经典的洛伦兹吸引子的程序如下：

% 第一个程序：lorenzfun.m

function dy=lorenzfun(t,y)

p=10

r=28

b=8/3

%p=16

%b=4.0

%r=45.92

dy=zeros(3,1)

dy(1,:)=p*(-y(1)+y(2))

dy(2,:)=r*y(1)-y(2)-y(1)*y(3)

dy(3,:)=y(1)*y(2)-b*y(3)

% 第一个程序：lorenztest.m

[t,yy]=ode45('lorenzfun',[0:0.01:100],[-1

0 1]) %模拟t=0~100步长0.01的10001个点，初值[-1

01]

x=yy(:,1)

y=yy(:,2)

z=yy(:,3)

figure(1)

plot3(x,y,z)

xlabel('x(t)'),ylabel('y(t)'),zlabel('z(t)')

title('Lorenz吸引子图')

x=x(1001:end)%抛弃前面的一些点

你在matlab下运行lorenztest看看能不能找到点感觉

---