如何绘制忆阻器混沌系统的奇异吸引子Matlab?

1、首先忆阻器混沌系统具有非常丰富的动态特性，通过绘制分析，不仅能观察到四翼混沌差敏吸引子。

2、其次观察到三翼混沌吸引子以及不同周期伏则和吸引子形状的周期轨，虚厅枝进行简单描绘。

3、最后抓住细节填充颜色，即可绘制出奇异吸引子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看看能不誉漏旦能找到点感觉

伊侬映射是答谨比较简单的混沌吸引子，离散映射方程为：xn+1=1+b*yn-a*xn^2yn+1=xn

我弄了个程序，不知是销州不是你想要的

clear all

clc

a=1.4b=0.3u=1.0e-4

x(1)=0.2y(1)=0.5

x1(1)=x(1)+uy1(1)=y(1)%初值敏感性亏举蔽

p=5000q=40

for m=2:p

x(m)=1+b*y(m-1)-a*(x(m-1))^2

x1(m)=1+b*y1(m-1)-a*(x1(m-1))^2

y(m)=x(m-1)

y1(m)=x1(m)-1

end

h=1:p

plot(x(h),y(h))；

---