siyaoming 发表于 2009-8-17 15:17

帮忙分析一下程序应该改哪里啊

这是一时滞方程的分岔图。
clf
w0=1;a=-1;b=-0.3;r=0.2;omega=1.1;g1=1;g2=1;f=0.1:0.1:5;
lags=;T=2*pi/omega;
y1=[];G=[];y2=[];
for h=1:length(f)
sol=dde23(@yanchi,lags,,,f(h));
y1=;y2=;
m=length(y1);
for j=1:m/h
G=;   
end;p=length(G);
N=length(f);
Over_to_step=N-i;
end
figure(1)
plot(G,y1,'k.','markersize',3)
hold on


function dy=yanchi(t,y,Z,f)
w0=1;
a=-1;
b=-0.3;
r=0.2;
omega=1.1;
g1=1;
g2=1;
f=0.1:0.1:5;
dy=

siyaoming 发表于 2009-8-17 15:17

这是运行是的matlab分析

??? Error using ==> vertcat
CAT arguments dimensions are not consistent.

Error in ==> yanchi at 10
dy=
Error in ==> dde23 at 220
f0 = feval(ddefun,t0,y0,Z0,varargin{:});

Error in ==> yanchi1 at 6
sol=dde23(@yanchi,lags,,,f(h));

无水1324 发表于 2009-8-17 16:16

回复 沙发 siyaoming 的帖子

你这里面Z与时间是什么关系，如果按照你这个样子，这个方程就不是简单单自由度系统了，刚清楚他们之间的关系再说吧

siyaoming 发表于 2009-8-17 16:41

院长，Z是表示时滞项的一种表示方法。
原方程是x''+a(1+b*x^2)*x'+w0^2(1+r*x^2)*x=f*cos(omega*t)+g1*x(t-tao1)+g2*x'(t-tao2)

[ 本帖最后由 无水1324 于 2009-8-18 22:15 编辑 ]

无水1324 发表于 2009-8-18 22:11

我怎么感觉你这里面的时滞量都是变化的呢？
难道它不是一个确定的值吗？

fanxinga 发表于 2011-6-6 21:59

这是时滞微分方程，解法是用dde23代替ode45执行，
却不料一个晚上弄不清原因。

569340337 发表于 2011-6-8 19:00

这是时滞微分方程，解法是用dde23代替ode45执行，
页: [1]
查看完整版本: 帮忙分析一下程序应该改哪里啊