声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 3509|回复: 2

[分形与混沌] 发一个我做的OGY与大家共享

[复制链接]
发表于 2006-7-23 10:40 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

x
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.');
回复
分享到:

使用道具 举报

发表于 2008-12-14 11:17 | 显示全部楼层

关于ogy法的一些问题,希望得到大家的帮助

原帖由 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];        % 响应矢量矩阵
怎么得到的?
发表于 2017-10-9 14:20 | 显示全部楼层
谢谢分享         
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-11-13 09:48 , Processed in 0.063518 second(s), 18 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表