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

siyaoming

这是一时滞方程的分岔图。
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=[pi/2,2*pi/3];T=2*pi/omega;
y1=[];G=[];y2=[];
for h=1:length(f)
sol=dde23(@yanchi,lags,[0.1 0],[0:T/100:40*T],f(h));
y1=[y1,sol.y(1,2000:100:end)];y2=[y2,sol.y(2,2000:100:end)];
m=length(y1);
for j=1:m/h
G=[G,f(h)];   
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=[y(2);-w0^2*y(1)-a*y(2)-a*b*y(1)^2*y(2)-w0^2*r*y(1)^3+f*cos(omega*t)+g1*Z(1,1)+g2*Z(2,2)]

siyaoming

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

Error in ==> yanchi at 10
dy=[y(2);-w0^2*y(1)-a*y(2)-a*b*y(1)^2*y(2)-w0^2*r*y(1)^3+f*cos(omega*t)+g1*Z(1,1)+g2*Z(2,2)]
Error in ==> dde23 at 220
f0 = feval(ddefun,t0,y0,Z0,varargin{:});

Error in ==> yanchi1 at 6
sol=dde23(@yanchi,lags,[0.1 0],[0:T/100:40*T],f(h));

无水1324

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

siyaoming

院长，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

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

fanxinga

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

569340337

这是时滞微分方程，解法是用dde23代替ode45执行，