matlab dde23原理,求教高手！帮忙改个错！关于dde23的用法

该楼层疑似违规已被系统折叠 隐藏此楼查看此楼

仿真如下例子：

x1'(t)=a*(x2(t)-m1*x1(t))+f1(x1(t))-c*x1(t-τ)+f2(x1(t-τ))

x2'(t)=x1(t)-x2(t)+x3(t)-c*x1(t-τ)

x3'(t)=-b*x2(t)+c*(2*x1(t-τ)-x3(t-τ))

其中：

f1(x1(t))=-a/2*(m0-m1)*(abs(x1(t)+1)-abs(x1(t)-1))

f2(x1(t-τ))=-2*c*sin(0.1*x1(t-τ))

已知：

m0=-1/7 m1=2/7 a=9 b=14.28 c=0.01 τ=0.1

初值： x(0)=[-0.2 -0.3 0.2]

我的M文件的程序为：

---------------------------function F-------------------------------------

function F=CHua_circuit1(t,y,Z)

m0=-1/7;

m1=2/7;

a=9;

b=14.28;

c=0.01;

ylag=Z(:,1);

f1=-a/2*(m0-m1)*(abs(y(1)+1)-abs(y(1)-1));

f2=-2*c*sin(0.1*ylag(1));

F=[a*(y(2)-m1*y(1))+f1-c*ylag(1)+f2

y(1)-y(2)+y(3)-c*ylag(1)

-b*y(2)+c*(2*ylag(1)-ylag(3))]

----------------------------function compact---------------------------

function compact

tf=100;

tspan=[0:tf/10000:100]

sol=dde23(@CHua_circuit1,[0.1],[-0.2 -0.3 0.2],tspan);

figure(1);

clf;

plot3(sol.y(:,1),sol.y(:,2),sol.y(:,3))

运行没有错误，但是出来的图形不正确

各位高手帮帮忙啊！小弟实在是查不出错在哪儿！