function [X,Y]=realTrack(T,totalTime)
%產(chǎn)生真實(shí)軌跡[X,Y]
%雷達(dá)掃描周期T
%v2.0@李敏 2005.6.22%數(shù)據(jù)初始化,X,Y表示以T為間隔采樣前的數(shù)據(jù)
x=zeros(totalTime,1);
y=zeros(totalTime,1);
%x,y表示以T為間隔采樣后的數(shù)據(jù)
X=zeros(ceil(totalTime/T),1);
Y=zeros(ceil(totalTime/T),1);
%非機(jī)動段：t=0:400s;速度vx,vy(m/s)
x0=2000;%起始點(diǎn)坐標(biāo)
y0=10000;
vx=0;
vy=-15;
%目標(biāo)真實(shí)軌跡1
for t=0:400
    x(t+1)=x0+vx*t;
    y(t+1)=y0+vy*t;
end
%開始機(jī)動1:t=400:600s
ax=0.075;
ay=-0.075;
%目標(biāo)真實(shí)軌跡2
for t=0:200
    x(t+401)=x(401)+vx*t+ax*t*t/2;
    y(t+401)=y(401)+vy*t+ay*t*t/2;
end
vx=vx+ax*(600-400);%600s末的速度
vy=vy+ay*(600-400);
%目標(biāo)真實(shí)軌跡3:t=600:610s勻速運(yùn)動
for t=0:10
    x(t+601)=x(601)+vx*t;
    y(t+601)=y(601)+vy*t;
end
%開始機(jī)動2:t=610:660s
ax=-0.3;
ay=-0.3;
%目標(biāo)真實(shí)軌跡4
for t=0:50
    x(t+611)=x(611)+vx*t+ax*t*t/2;
    y(t+611)=y(611)+vy*t+ay*t*t/2;
end
vx=vx+ax*(660-610);%660s末速度
vy=vy+ay*(660-610);
%目標(biāo)真實(shí)軌跡5:660s以后勻速運(yùn)動
for t=0:(totalTime-660)
    x(t+661)=x(661)+vx*t;
    y(t+661)=y(661)+vy*t;
end
%以T為間隔采樣
for n=0:T:totalTime
    X(n/T+1)=x(n+1);
    Y(n/T+1)=y(n+1);
end
%顯示真實(shí)軌跡
%plot(X,Y,'r');