myplot.m

雷达目标跟踪matlab源程序 希望对大家能够有所帮助

s1=size(all_x)
s1=s1(2);
s2=s1/2;%仿真总时间

plot(all_x(1,1:2:s1),all_x(3,1:2:s1),'r+',all_mes(1,1:12:s2*12),all_mes(2,1:12:s2*12),'k*',track_history(1).x(1,1:s2),track_history(1).x(3,1:s2),'g*')
legend('目标1真实','量测','immjpda')
u1=u_all(:,1:2:s1)%目标1跟踪门向量
u2=u_all(:,2:2:s1)


for i=1:53
  temp=[u1(1,i),u1(3,i),u1(2,i)-u1(1,i),u1(4,i)-u1(3,i)]
  rectangle('Position',temp);
end

pause;
plot(1:s2,all_x(1,1:2:s1)-all_mes(1,1:12:s2*12),'r-*',1:s2,all_x(1,1:2:s1)-track_history(1).x(1,1:s2),'g-*');
pause;
plot(1:s2,all_x(3,1:2:s1)-all_mes(2,1:12:s2*12),'k-*',1:s2,all_x(3,1:2:s1)-track_history(1).x(3,1:s2),'y-*');

legend('目标1真实-量测-x','真实-滤波-x','目标1真实-量测-y','真实-滤波-y')
pause;
plot(1:s2,all_x(1,1:2:s1)-track_history(1).x(1,1:s2),'g-*',1:s2,all_x(3,1:2:s1)-track_history(1).x(3,1:s2),'y-*',1:s2,all_x(2,1:2:s1)-track_history(1).x(2,1:s2),'b-*',1:s2,all_x(4,1:2:s1)-track_history(1).x(4,1:s2),'k-*');
legend('真实-滤波位置-x','真实-滤波位置y','真实-滤波速度-x','真实-滤波速度y')


%目标2
pause;
plot(all_x(1,2:2:s1),all_x(3,2:2:s1),'r+',all_mes(1,2:12:s2*12),all_mes(2,2:12:s2*12),'k*',track_history(2).x(1,1:s2),track_history(2).x(3,1:s2),'g*')
legend('目标2真实','量测','immjpda')
u1=u_all(:,1:2:s1)%目标1跟踪门向量
u2=u_all(:,2:2:s1)


for i=1:s2
  temp=[u2(1,i),u2(3,i),u2(2,i)-u2(1,i),u2(4,i)-u2(3,i)]
  rectangle('Position',temp);
end
plot(1:s2,all_x(1,2:2:s1)-all_mes(1,2:12:s2*12),'r-*',1:s2,all_x(1,2:2:s1)-track_history(2).x(1,1:s2),'g-*');

plot(1:s2,all_x(3,2:2:s1)-all_mes(2,2:12:s2*12),'k-*',1:s2,all_x(3,2:2:s1)-track_history(2).x(3,1:s2),'y-*');

legend('目标1真实-量测-x','真实-滤波-x','目标1真实-量测-y','真实-滤波-y')