本文目录一览:
- 1、如何用matlab模拟混沌系统或者动力学系统
- 2、混沌系统D2x+0.1*Dx+x^6=sin(t) 的matlab编程
- 3、混沌系统的同步,你有用matlab编写相关的仿真程序吗?
- 4、如何绘制忆阻器混沌系统的奇异吸引子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。
发布于 2022-08-26 02:01:24 回复
发布于 2022-08-26 07:26:42 回复
发布于 2022-08-26 06:39:22 回复