www.pudn.com > Target_Tracking_Simulation_Program.rar > realTrack.m, change:2005-06-29,size:1185b

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