混沌分形----三翅鹰.pptVIP

混沌分形-----三翅鹰 一、迭代方程及其物理意义 1980年, 物理学家 Gumowski, I. 和 Mira, C. 尝试计算模拟基本粒子轨迹在加速器中的行为 ： X(n+1) = B*Y + F(X) Y(n+1) = -X + F(X(n+1)) 其中： F(X) = A*X + (1-A)*2*X*X/(1+X*X) 二、主要的matlab程序 function aaa clear b=0.93; x=1;y=1; for k=1:20000;%迭代次数 x(k+1)=b*y(k)+fff(x(k)); y(k+1)=-x(k)+fff(x(k+1)); end plot(x,y,‘.’,‘markersize’,2） 三、程序运行结果 （1）当a=-0.45;b=0.93;x=-15;y=-15; （2）当a=-0.45;b=0.93;x=1;y=1; （3）当a=0;b=0.93;x=1;y=1; 四、总结 混沌系统具有初值敏感性 当初值微小变化时，它的响应将发生很大的变化。 当初值发生较大的变化时，它的轨迹将向外膨胀，螺旋发散至无穷。 Logistic方程混沌同步 一、logistic方程 驱动系统： x(n+1)=a*x(n)-a*x(n)*x(n); 响应系统： y(n+1)=b*y(n)-b*y(n)^2+u(n); 控制率： u(n)=k*(x(n)-y(n))-p*(x(n)*x(n)-y(n)*y(n)); 当u=[3.57，4]时，系统处于混沌状态。 二、参数相同时混沌同步 （1）参数a=4;b=4;k=p=4; E为系统稳态误差 选择不同的初值：x=0.2,y=0; （2）改变控制参数， 当k=p=4.5时 ： 三、参数未知时的混沌同步 （1）系统混沌同步的条件 ： 系统一：x(n+1)=a*x(n)-a*x(n)*x(n); 系统二：y(n+1)=u(n)*x(n)-u(n)*x(n) )*x(n); 其中：u(n+1)=k*(x(n+1)-y(n+1))+u(n); 根据稳态误差 ?0 ，经过计算得出： 0k8 是系统同步的充分条件。 （2）当a=3.9; 初值x=0.1 ; k=5时 ： 当控制参数 k 减小到 k=1 时: E—n曲线： 四、参考程序 function abc y=0;u=0; x=0.1; E=0; k=input(请输入控制参数k=); for n=1:100 u(n+1)=k*(x(n)-y(n))+u(n); x(n+1)=3.9*x(n)-3.9*x(n)*x(n); y(n+1)=u(n)*x(n)-u(n)*x(n)*x(n); E(n)=abs(y(n)-x(n));%系统误差； end plot(E) 五、小结 当两个系统的参数已知且相等时： 系统能否达到同步的关键在于两个系统的初始值和控制器的系数大小； 当一个系统的参数未知时，采用自适应混沌控制方法简单，快速。只需满足一个条件0k8即可；当k越大时，系统同步的速度就越快。 * * -------采用参数自适应方法 可见，在初始误差相同的情况下，响应随着控制参数大小而改变。 误差（E）—n曲线： *