linzx270 发表于 2008-4-22 08:43

关于如何画时滞系统的问题

用matlab的什么语句能画出时滞系统中的t-(t-T)和t-x的图象?
谢谢了可以加我qq45574874
程序如下:
a=35;
b=3;
c=28;
tau=2;
tspan = ;
opts = ddeset('RelTol',1e-5,'AbsTol',1e-8);
history = [-3;-4;2];
sol = dde23(@dde_chenf,tau,history,tspan,opts,a,b,c)
t=0:0.1:4000;
plot3(t,sol.y(1,4000:end))

function dydt = dde_chenf(t,y,Z,a,b,c)
ylag = Z(:,1);
dydt = [ -a*y(1)+a*y(2)-ylag(1)
      (c-a)*y(1)+c*y(2)-ylag(2)-y(1)*y(3)
      (-b)*y(3)+ylag(3)+y(1)*y(2)      ];
问题是plot3(t,sol.y(1,4000:end))出错

[ 本帖最后由 linzx270 于 2008-4-22 10:51 编辑 ]

earth123654 发表于 2013-3-10 21:04

请问楼主的问题解决了吗
页: [1]
查看完整版本: 关于如何画时滞系统的问题