xueyongzhou 发表于 2007-8-19 21:58

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

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

谢谢大侠帮忙!!!

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

shenyongjun 发表于 2007-8-19 22:02

回复 #1 xueyongzhou 的帖子

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

xueyongzhou 发表于 2007-8-19 22:04

回复 #2 shenyongjun 的帖子

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

shenyongjun 发表于 2007-8-19 22:05

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

sssssxxxxx921 发表于 2007-8-19 22:11

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

sssssxxxxx921 发表于 2007-8-19 22:12

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

无水1324 发表于 2007-8-19 22:26

回复 #5 sssssxxxxx921 的帖子

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

咕噜噜 发表于 2007-8-20 08:22

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

xueyongzhou 发表于 2007-8-20 09:01

我的程序如下: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,,);
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 发表于 2007-8-20 10:35

回复 #9 xueyongzhou 的帖子

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

无水1324 发表于 2007-8-20 10:36

回复 #9 xueyongzhou 的帖子

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

xueyongzhou 发表于 2007-8-20 11:01

回复 #11 无水1324 的帖子

是的

xueyongzhou 发表于 2007-8-20 11:02

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

无水1324 发表于 2007-8-20 11:17

回复 #13 xueyongzhou 的帖子

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

octopussheng 发表于 2007-8-20 15:40

回复 #13 xueyongzhou 的帖子

仅仅从相图结构去分析系统是远远不够的，还是要考虑其他方法的综合，如时程曲线、功率谱、Poincare截面、Lyapunov指数等等，广从相轨迹上是看不出太多东西的！
页: [1] 2
查看完整版本: 看看我的图形是什么意思?