发一个我做的OGY与大家共享

ieluxinhua

clear
clc
A=1.4;                              % 系统在A=1.4，B=0.3时呈现混沌运动
B=0.3;
d=0.08;
x=0.32;
y=0.32;% 系统初始值 x=y=0.32
xF=0.8839;   % 不动点（xF,xF)
f=[-2*xF B;1 0];                    % 系统在xF处的导数
g=((1-B).^2+4*A).^0.5*[1;1];        % 响应矢量矩阵
ju=xF+(xF.^2+B).^0.5;              % 系统在xF点处不稳定方向上的特征根
eu=[1;-((xF.^2+B).^0.5-xF)/B];      % 系统在xF点处不稳定方向上单位向量
hu=[((xF.^2+B).^0.5+xF)/(2*((xF.^2+B).^0.5)) B/(-2*((xF.^2+B).^0.5))];%系统在不稳定方向上正交对偶矢量
z=zeros(1,1000);
P=zeros(1,1000);
for n=1:1000
      xpre=x;
      ypre=y;
      diag=[x-xF;y-xF];         % 控制的目的是使x移到点xF的稳定方向上，即使diag在不稳定方向上的分量为零     
   if n<100
      P(n)=0;
   else
      P(n)=(ju*hu*diag)/((ju-1)*(hu*g)); % 控制率表达
  end
      x=A+B*xpre*ypre-xpre.^2+P(n);      % 系统动力学方程   
      y=xpre;
      z(n)=z(n)+x;
   end
  plot(z,'-k.');

findingnemo

原帖由 ieluxinhua 于 2006-7-23 10:40 发表

                               
登录/注册后可看大图

clear
clc
A=1.4;                              % 系统在A=1.4，B=0.3时呈现混沌运动
B=0.3;
d=0.08;
x=0.32;
y=0.32;% 系统初始值 x=y=0.32
xF=0.8839;   % 不动点（xF,xF)
f=[-2*xF B;1 0];              ...

首先感谢ieluxinhua 的ogy法帖子。
我也是初学混沌控制，这个帖子有一些看不懂的地方，希望大家指点迷津：
1 请问f=[-2*xF B;1 0];  是由
x=A+B*xpre*ypre-xpre.^2+P(n);      % 系统动力学方程   
y=xpre;
得到的吗，我推出的f=[-2*xF+B*xF B;1 0];，怎么解释呢？
2   文中的g=((1-B).^2+4*A).^0.5*[1;1];        % 响应矢量矩阵
怎么得到的？

sizhiyuan2006

谢谢分享