j.m

杂波的一个MATLAB仿真

%何友-p93例子——杂波点和真实点的态势图的仿真 
%2009.4.20

%五个做匀速直线运动的目标的初始位置
x1(1)=55000;
y1(1)=55000;
x2(1)=45000;
y2(1)=45000;
x3(1)=35000;
y3(1)=35000;
x4(1)=45000;
y4(1)=25000;
x5(1)=15000;
y5(1)=55000;

%5个目标的速度
vx=500;
vy=0;

%采样周期
t=5;

%连续四个扫描周期
i=0:t:15;
n=size(i,2);

%5个目标的真实轨迹
for k=2:n
x1(k)=x1(k-1)+vx*t;
y1(k)=y1(k-1)+vy*t;
x2(k)=x2(k-1)+vx*t;
y2(k)=y2(k-1)+vy*t;
x3(k)=x3(k-1)+vx*t;
y3(k)=y3(k-1)+vy*t;
x4(k)=x4(k-1)+vx*t;
y4(k)=y4(k-1)+vy*t;
x5(k)=x5(k-1)+vx*t;
y5(k)=y5(k-1)+vy*t;
end
plot(x1,y1,'o');
axis([0 100000 0 100000]);
hold on;
plot(x2,y2,'o');
hold on;
plot(x3,y3,'o');
hold on;
plot(x4,y4,'o');
hold on;
plot(x5,y5,'o');
hold on;   
str={'r*','bs','k+','g.'};
%杂波的生成
for i=1:4
    r=rand(1);    %生成（0，1）分布的随机数
    J=1;
    while (poisscdf(J-1,50)>=r)||(poisscdf(J,50)<r)
        J=J+1;       
    end
    n=J;%杂波的个数,n是不定的，那存入矩阵应该注意什么？
    plot(unifrnd(0,100000,1,n),unifrnd(0,100000,1,n),str{i})
    hold on;
end