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看看能不誉漏旦能找到点感觉
欢迎分享,转载请注明来源:内存溢出
评论列表(0条)