phasedraw.m

来自「陆地作战的模拟程序」· M 代码 · 共 39 行

lags = [17];
tx = dde23('dde', lags, 'f2', [-20, 400]);
%plot(tx.x,tx.y(1,:));
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














% 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);