该楼层疑似违规已被系统折叠 隐藏此楼查看此楼
m文件
function dx=c7exdde(t,x,z)
xlag1=z(:,1);xlag2=z(:,2);
dx=[36*xlag1(2)-36*xlag1(1)+x(4);
-x(1)*x(3)+20*xlag1(2);
x(1)*x(2)-3*xlag2(3);
x(1)*x(3)+1.3*xlag2(4); 35*xlag1(6)-35*xlag1(5)-x(4)-x(9)*xlag1(2)-x(9)*xlag1(1)-x(13)*xlag1(6)-x(13)*xlag1(5)-x(1)-x(5);
7*xlag1(5)+12*xlag1(6)+x(1)*x(3)-xlag1(2)*x(10)-xlag1(6)*x(14)-xlag1(6)*x(14)-xlag1(5)*x(15)-x(2)-x(6);
-3*xlag2(7)-x(1)*x(2)+xlag2(3)*x(11)+xlag2(7)*x(16)-x(3)-x(7);
0.5*xlag2(8)-x(1)*x(3)-xlag2(4)*x(12)-xlag2(8)*x(17)-x(4)-x(8);
xlag1(2)*x(1)+xlag1(2)*x(5)-xlag1(1)*x(1)-xlag1(1)*x(5);
xlag1(2)*x(2)+xlag1(2)*x(6);
-xlag2(3)*x(3)-xlag2(3)*x(7);
xlag2(4)*x(4)-xlag2(4)*x(8);
xlag1(6)*x(1)+xlag1(6)*x(5)-xlag1(5)*x(1)-xlag1(5)*x(5);
xlag1(6)*x(2)+xlag1(6)*x(6);
xlag1(5)*x(2)+ xlag1(5)*x(6);
-xlag2(7)*x(3)- xlag2(7)*x(7);
xlag2(8)*x(4)+ xlag2(8)*x(8)];
下面为执行函数
lags=[1 1];
tx=dde23('c7exdde',lags,[5;8;-1;-3;-2;-3;8;10;5;9;7;11;2;6;10;4;8],[0,3]);
plot(tx.x,(tx.y(3,:)+tx.y(7,:)))