请问为什么混沌控制成周期轨了

lll5008

function ydot=lorenz1(t,y)
ydot=zeros(3,1);
ydot(1)=10*(y(2)-y(1));
ydot(2)=28*y(1)-y(2)-y(1)*y(3)-0.1*(y(1)-72^(1/2))+72^(1/2)*(y(3)-27);
ydot(3)=y(1)*y(2)-8*y(3)/3;

>> syms t0 y0 t y
>> t0=[0 100];
>> y0=[1;2;3];
>> [t,y]=ode45(@lorenz1,t0,y0);
>> plot3(y(:,3),y(:,1),y(:,2))

是应为方程里的控制系数0.1取的太小了，所以误差引起的么？

superstones

你对它进行了控制，肯定就不再是原来的lorenz系统了
ydot(1)=10*(y(2)-y(1));
ydot(2)=28*y(1)-y(2)-y(1)*y(3)-0.1*(y(1)-72^(1/2))+72^(1/2)*(y(3)-27);
ydot(3)=y(1)*y(2)-8*y(3)/3
和
ydot(1)=10*(y(2)-y(1));
ydot(2)=28*y(1)-y(2)-y(1)*y(3)
ydot(3)=y(1)*y(2)-8*y(3)/3
的相图肯定是不同的

无水1324

原帖由 lll5008 于 2007-4-19 16:57 发表

                               
登录/注册后可看大图

function ydot=lorenz1(t,y)
ydot=zeros(3,1);
ydot(1)=10*(y(2)-y(1));
ydot(2)=28*y(1)-y(2)-y(1)*y(3)-0.1*(y(1)-72^(1/2))+72^(1/2)*(y(3)-27);
ydot(3)=y(1)*y(2)-8*y(3)/3;

>> syms t0 y0 t y
> ...

看你一下你的程序

syms 这句话是不是多余的？


还有你的程序施加了控制吗？因为你的标题“为什么混沌控制成周期...”
看不懂！

xjzuo

如楼上所言, syms t0 y0 t y  这句可以去掉.
我想LZ所说的＂控制＂大概是指调节控制参数，如调节为0.1.

＂混沌控制成周期轨了＂－－－这个要求LZ查查资料，选择合适的控制参数．

无水1324

那这个也不能说是控制了混沌呀
只能说在参数区域内存在周期轨道，此参数为0.1；

lll5008

哦,可是从理论上参数只要>0就可以实现控制到稳定点的,现在画的图确是周期轨,当参数到3是控制到稳定点