matlab混沌系统-混沌 matlab

本文目录一览:

如何用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;

0; 1]

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看看能不能找到点感觉

混沌系统D2x+0.1*Dx+x^6=sin(t) 的matlab编程

楼主不是开玩笑吧,居然还会向楼上这样明显胡闹的人追问。

这个方程无法求出解析解(绝大多数混沌系统应该都是球不出解析解的吧),只能用数值方法求解。以下给出参考代码:

f=inline('[x(2); -0.1*x(2)-x(1)^5+6*sin(t)]','t','x')

[t,y]=ode45(f,100,[2; 3]);

comet(y(:,1),y(:,2))代码很短,看看相关函数的参考应该能搞明白,如果还有疑问再追问吧。

捎带说一下,你出的错是把“function”写成了“fuction”。

混沌系统的同步,你有用matlab编写相关的仿真程序吗?

只有 用matlab写的两次同步的误差系统的模拟图像。

function y=wuchaxitong(t,x);

a=35;b=4;c=25;d=5;e=35;f=100;

k1=1277;k2=1388;k3=55;k4=1;

y=zeros(8,1);

y(1)=a*((x(2)-x(6))-(x(1)-x(5)))+e*(x(2)-x(6))*(x(3)-x(7))+(x(3)-x(7))+e*(x(2)-x(6))*x(7)+e*x(6)*(x(3)-x(7))-k1*(x(1)-x(5));

y(2)=c*(x(1)-x(5))-d*(x(1)-x(5))*(x(3)-x(7))-d*x(5)*(x(3)-x(7))-d*x(7)*(x(1)-x(5))-(x(2)-x(6))+(x(4)-x(8))-k2*(x(2)-x(6));

y(3)=(x(1)-x(5))*(x(2)-x(6))+(x(1)-x(5))*x(6)+x(5)*(x(2)-x(6))-b*(x(3)-x(7))-k3*(x(3)-x(7));

y(4)=-f*(x(2)-x(6))-(k4+1)*(x(4)-x(8));

tspan=[0:0.001:4];%求解时间为2秒也可以 tspan=[0 2]表示自动时间步长

x0=[-1;30;2;10;-3;4;0.2;6];%初值

[T,Y]=ode45('wuchaxitong',tspan,x0);%用龙格库塔法求解

subplot(4,1,1),plot(T,Y(:,1)),xlabel('t/s'),ylabel('e1')

subplot(4,1,2),plot(T,Y(:,2)),xlabel('t/s'),ylabel('e2')

subplot(4,1,3),plot(T,Y(:,3)),xlabel('t/s'),ylabel('e3')

subplot(4,1,4),plot(T,Y(:,4)),xlabel('t/s'),ylabel('e4')

上面的是误差系统,下面的是画的误差系统的图像,不知道是不是你需要的,我做的是四维的系统

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

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

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

3、最后抓住细节填充颜色,即可绘制出奇异吸引子Matlab。


原文链接:http://527256.com/29438.html

相关文章

访客
访客
发布于 2022-08-26 02:01:24  回复
35;f=100;k1=1277;k2=1388;k3=55;k4=1;y=zeros(8,1);y(1)=a*((x(2)-x(6))-(x(1)-x(5)))+e*(x(2)-x(6))*(x(3)-x(7))+(x(3)-x(7))+e*(x(2)-x(6))*x(7)+e*x(6)*(
访客
访客
发布于 2022-08-26 07:26:42  回复
0:0.001:4];%求解时间为2秒也可以 tspan=[0 2]表示自动时间步长x0=[-1;30;2;10;-3;4;0.2;6];%初值[T,Y]=ode45('w
访客
访客
发布于 2022-08-26 06:39:22  回复
; 1]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吸引

发表评论:

◎欢迎参与讨论,请在这里发表您的看法、交流您的观点。

返回顶部