马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
在网上看了一篇论文:孙海云,曹庆杰.相空间重构延迟时间的确定.山东工业大学学报.2000.30(2),有个例子方程是这样的dx/dt=A-(B+1)*x+x^2*y+a*cos(x*t);
dy/dt=B*x-x^2*y;
然后用matlab用最基本的方法和相空间重构的方法分布进行了相空间绘制,我把基本的相图作法程序贴出,但重构的那种方法我一直做不出来请高手指教
Brusselator.m文件为
function dy=Brusselator(t,y);
dy=zeros(2,1);
dy(1)=0.4-2.2*y(1)+y(1)^2*y(2)+0.08*cos(0.86*t);
dy(2)=1.2*y(1)-y(1)^2*y(2);
运行下面代码
% [T,Y]=ode45('Brusselator',[0 1000],[1 0]);
% plot(Y(:,1),Y(:,2))
请问重构的程序怎么做,我想画出如文章中的图
[ 本帖最后由 jyxhappy 于 2008-7-12 14:47 编辑 ] |