delaye.m

陆地作战的模拟程序

InlayLength = 4; % 4 6 8

lags = [17];
tx = dde23('dde', lags, 'f2', [0, 400]);
%plot(tx.x,tx.y(1,:));

% producing time series data ----------- begin --------------
tint = 0:0.05:400;
yint = deval(tx, tint);
fp = fopen('ddets.txt','w');
for i = 1:8000 - InlayLength
    for j = 1:InlayLength
        fprintf(fp,'%f\t',yint(i-1+j));
    end
    fprintf(fp,'\n');
end
fclose(fp);

% producing time series data ----------- end --------------


% Phase chart----------------------------begin---------------------
% tint1=17:0.1:400;
% yint1 = deval(tx, tint1);
% tint2=-0:0.1:383;
% yint2 = deval(tx, tint2);
% plot(yint1,yint2);
% hold on
% svt=load('sv.txt');
% svx=deval(tx, svt(:,1));
% svt1=svt(:,1)-17;
% svx1=deval(tx, svt1(:));
% plot(svx(:),svx1(:),'r.');
% hold off
% Phase chart---------------------------- end ----------------------

% tint = 0:400;
% yint = deval(tx, tint);
% plot(tint,yint);
% fp = fopen('dde.txt','w');
% fprintf(fp, '%f\n',tint);
% fprintf(fp, '%f\n',yint);
% fclose(fp);