发一个我做的OGY与大家共享
clearclc
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*; % 响应矢量矩阵
ju=xF+(xF.^2+B).^0.5; % 系统在xF点处不稳定方向上的特征根
eu=; % 系统在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的稳定方向上，即使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.');
关于ogy法的一些问题，希望得到大家的帮助
原帖由 ieluxinhua 于 2006-7-23 10:40 发表 http://www.chinavib.com/forum/images/common/back.gifclear
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]