function [X,Y]=realTrack(T,totalTime)
%產生真實軌跡[X,Y]
%雷達掃描周期T
%v2.0@李敏 2005.6.22%數據初始化,X,Y表示以T為間隔采樣前的數據
x=zeros(totalTime,1);
y=zeros(totalTime,1);
%x,y表示以T為間隔采樣后的數據
X=zeros(ceil(totalTime/T),1);
Y=zeros(ceil(totalTime/T),1);
%非機動段：t=0:400s;速度vx,vy(m/s)
x0=2000;%起始點坐標
y0=10000;
vx=0;
vy=-15;
%目標真實軌跡1
for t=0:400
    x(t+1)=x0+vx*t;
    y(t+1)=y0+vy*t;
end
%開始機動1:t=400:600s
ax=0.075;
ay=-0.075;
%目標真實軌跡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);
%目標真實軌跡3:t=600:610s勻速運動
for t=0:10
    x(t+601)=x(601)+vx*t;
    y(t+601)=y(601)+vy*t;
end
%開始機動2:t=610:660s
ax=-0.3;
ay=-0.3;
%目標真實軌跡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);
%目標真實軌跡5:660s以后勻速運動
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
%顯示真實軌跡
%plot(X,Y,'r');