如何用matlab模拟混沌系统或者动力学系统

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

% 第一个程序：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看看能不誉漏旦能找到点感觉

clearclfhold on

axis([0,3,-3,3])grid

for a=0:0.005:3 x=[0.1]for i=2:150

x(i)=a*sin(pi*x(i-1))end pause(0.1) for i=101:150

plot(a,x(i),'k.')end end

---