realtrack.m

雷达数据预测

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