看看我的图形是什么意思?

xueyongzhou

1.doc (38.5 KB, 下载次数: 26)

2007-8-19 22:01 上传

点击文件名下载附件

我用MATLAB得到相图(见附件).我想知道这是多周期的吗?会产生混沌吗?见笑了!!!

谢谢大侠帮忙!!!

[ 本帖最后由 xueyongzhou 于 2007-8-19 22:02 编辑 ]

shenyongjun

好像是混沌的，但是最好用Poincare截面做一下！

xueyongzhou

用Poincare截面怎么做?我的是时滞方程!

shenyongjun

论坛中关于Poincare截面的讨论和程序很多，你找一找吧

sssssxxxxx921

我把他的图形发上来 看着更直接点  
因为我看不出来几周期    只能打个下手   让大牛看看  也省大家的时间

sssssxxxxx921

这还是个三维的   你能不能把你的程序发上来让大家看看啊

无水1324

从这个图确实看不出什么东西来

咕噜噜

还是画成2维的看的清楚，这样还真是看不明白

xueyongzhou

我的程序如下:clc;clear all
global r K beta mu c m delta h tau
r=1; K=500; beta=0.001; mu=0.03; c=8; m=150; delta=0.2; h=0.5; tau=40;
sol=dde23('Hopfbifurcation',tau,[1.50;1.00;0.20],[0,3000]);
figure(1);
plot(sol.x,sol.y(1,:),'r');xlabel('t');ylabel('T');
figure(2);
plot(sol.x,sol.y(2,:),'r');xlabel('t');ylabel('I');
figure(3);
plot(sol.x,sol.y(3,:),'r');xlabel('t');ylabel('V');
figure(4);
plot3(sol.y(1,:),sol.y(2,:),sol.y(3,:));
xlabel('T');ylabel('I'),zlabel('V');




function y = Hopfbifurcation(t,x,z)
global r K beta mu c m delta h tau
xlag=z(:,1);
y=zeros(3,1);    % a column vector
y(1)=r.*x(1).*(1-x(1)/K)-beta.*x(1).*xlag(2);  
y(2)=beta.*x(1).*xlag(2)-mu.*x(2)-c.*x(2).*x(3)/(m.*x(3)+x(2));
y(3)=delta.*x(3).*(1-h.*x(3)/x(2));

无水1324

这不是周期的，你可以调整一下时间延迟量，应该就可以出现混沌、周期状态

无水1324

对于时滞系统一般是研究时滞量对系统响应的影响吗？

xueyongzhou

是的

xueyongzhou

对于时间延迟方程怎么去研究混沌、周期状态?

无水1324

1、定性的分析，你可以参考胡海岩教授的一些文献；
2、对于时滞系统，研究其混沌周期，也是从时间位移曲线、相图、Poincare相图、分岔图等去分析；

octopussheng

仅仅从相图结构去分析系统是远远不够的，还是要考虑其他方法的综合，如时程曲线、功率谱、Poincare截面、Lyapunov指数等等，广从相轨迹上是看不出太多东西的！